The physical nature underlying intermolecular interactions between two rod-like winter flounder antifreeze protein (AFP) molecules and their implication for the mechanism of antifreeze function are examined in this work using molecular dynamics simulations, augmented with free energy calculations employing a continuum solvation model. The energetics for different modes of interactions of two AFP molecules is examined in both vacuum and aqueous phases along with the water distribution in the region encapsulated by two antiparallel AFP backbones. The results show that in a vacuum two AFP molecules intrinsically attract each other in the antiparallel fashion, where their complementary charge side chains face each other directly. In the aqueous environment, this attraction is counteracted by both screening and entropic effects. Therefore, two nearly energetically degenerate states, an aggregated state and a dissociated state, result as a new aspect of intermolecular interaction in the paradigm for the mechanism of action of AFP. The relevance of these findings to the mechanism of function of freezing inhibition in the context of our work on Antarctic cod antifreeze glycoprotein (Nguyen et al., Biophysical Journal, 2002, Vol. 82, pp. 2892-2905) is discussed. (C) 2004 Wiley Periodicals, Inc.

Antifreeze proteins (AFPs) and antifreeze glycoproteins (AFGPs), collectively abbreviated as AF(G)Ps, are synthesized by various organisms to enable their cells to survive in subzero environments. Although the AF(G)Ps are markedly diverse in structure, they all function by adsorbing to the surface of embryonic ice crystals to inhibit their growth. This adsorption results in a freezing temperature depression without an appreciable change in the melting temperature. The difference between the melting and freezing temperatures, termed thermal hysteresis (TH), is used to detect and quantify the antifreeze activity. Insights from crystallographic structures of a number of AFPs have led to a good understanding of the ice-protein interaction features. Computational studies have focused either on verifying a specific model of AFP-ice interaction or on understanding the protein-induced changes in the ice crystal morphology. In order to explain the origin of TH, we propose a novel two-dimensional adsorption kinetic model between AFPs and ice crystal surfaces. The validity of the model has been demonstrated by reproducing the TH curve on two different beta-helical AFPs upon increasing the protein concentration. In particular, this model is able to accommodate the change in the TH behavior observed experimentally when the size of the AFPs is increased systematically. Our results suggest that in addition to the specificity of the AFPs for the ice, the coverage of the AFPs on the ice surface is an equally necessary condition for their TH activity. (c) 2006 American Institute of Physics.