Detailed Analysis of the Ice Surface after Binding of an Insect Antifreeze Protein and Correlation with the Gibbs-Thomson Equation

Antifreeze proteins (AFPs) are able to influence the ice crystal growth and the recrystallization process due to the Gibbs-Thomson effect. The binding of the AFP leads to the formation of a curved ice surface and it is generally assumed that there is a critical radius between the proteins on the ice surface that determines the maximal thermal hysteresis. Up to now, this critical radius has not yet been proven beyond doubt or only in poor agreement with the Gibbs-Thomson equation. Using molecular dynamics (MD) simulations, the resulting three-dimensional surface structure is analyzed and the location of the critical radius is identified. Our results demonstrate that the correct analysis of the geometry of the ice surface is extremely important and cannot be guessed upfront a simulation. In contrary to earlier expectations from the literature, we could show that the critical radius is not located directly between the adsorbed proteins. In addition, we showed that the minimum temperature at which the system does not freeze is in very good agreement with the value calculated with the Gibbs-Thomson equation at the critical radius, as long as dynamic system conditions are taken into account. This proves on the one hand that the Gibbs-Thomson effect is the basis of thermal hysteresis and that MD simulations are suitable for the prediction of the melting point depression.

The influence of a type III antifreeze protein and its mutants on methane hydrate adsorption-inhibition: a molecular dynamics simulation study

Antifreeze proteins (AFPs) inhibit ice growth in various organisms at subzero temperature. Recently, AFPs as a hydrate inhibitor have been a topic of intense discussion, while the detailed mechanism remains obscure. The present work aims to explore molecular insight into the adsorption and inhibition of an AFP III on methane hydrate. Three polar, hydrophilic, and neutral amino acids (Asn14, Thr18, and Gln44) are mutated to elucidate the molecular mechanism of AFP III antifreeze activity. Another triple mutation is also designed to investigate the effect of the side chain. Atomistic molecular dynamics simulations provide detailed structural and dynamical aspects of protein residues and water molecules at the hydrate/water interface. Initially, it was proposed that the AFP III operates by the adsorption-inhibition mechanism on hydrates, almost similar to that of ice. The exchange of amide and hydroxyl groups by mutagenesis alters the shape of the side chain and the capability of hydrogen bonding and demonstrates that hydrogen bonds are not directly responsible for the AFP III antifreeze activity. Moreover, we deciphered that the length of the pendant group is an important factor in the entrapment of the AFP III on the hydrate cages, which is compatible with van der Waals interactions between the side chains and hydrate surface. The results suggest that this interaction is sensitive to the geometry and shape of the hydrate-binding surface (HBS) of the AFP, which implies that the interface between hydrates and the AFP is relatively rigid.

Anisotropy in Stable Conformations of Hydroxylate Ions between the {001} and {110} Planes of TiO2 Rutile Crystals for Glycolate, Lactate, and 2-Hydroxybutyrate Ions Studied by Metadynamics Method

Control over TiO2 rutile crystal growth and morphology using additives is essential for the development of functional materials. Computer simulation studies on the thermodynamically stable conformations of additives at the surfaces of rutile crystals contribute to understanding the mechanisms underlying this control. In this study, a metadynamics method was combined with molecular dynamics simulations to investigate the thermodynamically stable conformations of glycolate, lactate, and 2-hydroxybutyrate ions at the {001} and {110} planes of rutile crystals. Two simple atom-atom distances were selected as collective variables for the metadynamics method. At the {001} plane, a conformation in which the COO- group was oriented toward the surface was found to be the most stable for the lactate and 2-hydroxybutyrate ions, whereas a conformation in which the COO- group was oriented toward water was the most stable for the glycolate ion. At the {110} plane, a conformation in which the COO- group was oriented toward the surface was the most stable for all three hydroxylate ions, and a second most stable conformation was also observed for the lactate ion at positions close to the {110} plane. For all three hydroxylate ions (α-hydroxycarboxylate ions), the stability of the most stable conformation was higher for the {110} plane than for the {001} plane. At both planes, the stability of the most stable conformation was highest for the 2-hydroxybutyrate ion and lowest for the glycolate ion. Supposing that all three hydroxylate ions serve to decrease the surface free energy at the rutile surface and that a more stable conformation at the rutile surface leads to a greater decrease in the surface free energy, the present results partially explain experimentally observed differences in the changes in growth rate and morphology of rutile crystals in the presence of glycolic, lactic, and 2-hydroxybutyric acids.

Crystal-plane-dependent effects of antifreeze glycoprotein impurity for ice growth dynamics

An impurity effect on ice crystal growth in supercooled water is an important subject in relation to ice crystal formation in various conditions in the Earth's cryosphere regions. In this review, we consider antifreeze glycoprotein molecules as an impurity. These molecules are well known as functional molecules for controlling ice crystal growth by their adsorption on growing ice/water interfaces. Experiments on free growth of ice crystals in supercooled water containing an antifreeze protein were conducted on the ground and in the International Space Station, and the normal growth rates for the main crystallographic faces of ice, namely, basal and prismatic faces, were precisely measured as functions of growth conditions and time. The crystal-plane-dependent functions of AFGP molecules for ice crystal growth were clearly shown. Based on the magnitude relationship for normal growth rates among basal, prismatic and pyramidal faces, we explain the formation of a dodecahedral external shape of an ice crystal in relation to the key principle governing the growth of polyhedral crystals. Finally, we emphasize that the crystal-plane dependence of the function of antifreeze proteins on ice crystal growth relates to the freezing prevention of living organisms in sub-zero temperature conditions. This article is part of the theme issue 'The physics and chemistry of ice: scaffolding across scales, from the viability of life to the formation of planets'.

Molecular Dynamics Simulation of Ice Crystal Growth Inhibition by Hexadecyl-trimethyl-ammonium Bromide

Recent experiments have found hexadecyl-trimethyl-ammonium bromide (CTAB) to have superior ice nucleation inhibition properties [ J. Phys. Chem. B 121, 6580]. The mechanism of how the inhibition takes place remains unclear. Therefore, molecular dynamics was used to simulate ice crystallization of a water/CTAB/ice system. The ice crystallization rate for a pure water system was compared for the basal [0001], first prism [101̅0], and secondary prism plane [112̅0], where the basal plane grew the slowest followed by the first prism plane. When CTAB was added to the ice-liquid water system, crystallization was clearly impeded. Even when ice starts growing away from the CTAB molecule, the hydrophilic head would at some point protrude and get caught in the water/ice interface. Once the head of the CTAB was encapsulated in the advancing interface, the hydrophobic body would wriggle around and disrupt the formation of hydrogen bond networks that are essential for ice growth. When the interface clears the length of the body of the CTAB molecule, ice crystallization resumes at its normal pace. In summary, the inhibition of ice growth is a combination of the hydrophilic head acting as an anchor and the dynamic motion of the hydrophobic tail hindering stable hydrogen bonding for ice growth.

Operation of Kelvin Effect in the Activities of an Antifreeze Protein: A Molecular Dynamics Simulation Study

Ice growth and melting inhibition activities of antifreeze proteins (AFPs) are better explained by the adsorption-inhibition mechanism. Inhibition occurs as a result of the Kelvin effect induced by adsorbed protein molecules onto the surface of seed ice crystal. However, the Kelvin effect has not been explored by the state-of-the-art experimental techniques. In this work, atomistic molecular dynamics simulations have been carried out with Tenebrio molitor antifreeze protein ( TmAFP) placed at ice-water interface to probe the Kelvin effect in the mechanism of AFPs. Simulations show that, below equilibrium melting temperature, ice growth is inhibited through the convex ice-water interface formation toward the water phase and, above equilibrium melting temperature, ice melting is inhibited through the concave ice-water interface formation inward to ice phase. Simulations further reveal that the radius of curvature of the interface formed to stop the ice growth increases with decrease in the degree of supercooling. Our results are in qualitative agreement with the theoretical prediction of the Kelvin effect and thus reveal its operation in the activities of AFPs.

Adsorption of three-domain antifreeze proteins on ice: a study using LGMMAS theory and Monte Carlo simulations

In the present work, the adsorption of three-domain antifreeze proteins on ice is studied by combining a statistical thermodynamics based theory and Monte Carlo simulations. The three-domain protein is modeled by a trimer, and the ice surface is represented by a lattice of adsorption sites. The statistical theory, obtained from the exact partition function of non-interacting trimers adsorbed in one dimension and its extension to two dimensions, includes the configuration of the molecule in the adsorbed state, and allows the existence of multiple adsorption states for the protein. We called this theory "lattice-gas model of molecules with multiple adsorption states" (LGMMAS). The main thermodynamics functions (partial and total adsorption isotherms, Helmholtz free energy and configurational entropy) are obtained by solving a non-linear system of j equations, where j is the total number of possible adsorption states of the protein. The theoretical results are contrasted with Monte Carlo simulations, and a modified Langmuir model (MLM) where the arrangement of the adsorption sites in space is immaterial. The formalism introduced here provides exact results in one-dimensional lattices, and offers a very accurate description in two dimensions (2D). In addition, the scheme is capable of predicting the proportion between coverage degrees corresponding to different conformations in the same energetic state. In contrast, the MLM does not distinguish between different adsorption states, and shows severe discrepancies with the 2D simulation results. These findings indicate that the adsorbate structure and the lattice geometry play fundamental roles in determining the statistics of multistate adsorbed molecules, and consequently, must be included in the theory.

Effect of antifreeze protein on heterogeneous ice nucleation based on a two-dimensional random-field Ising model

Antifreeze proteins (AFPs) are the key biomolecules that protect many species from suffering the extreme conditions. Their unique properties of antifreezing provide the potential of a wide range of applications. Inspired by the present experimental approaches of creating an antifreeze surface by coating AFPs, here we present a two-dimensional random-field lattice Ising model to study the effect of AFPs on heterogeneous ice nucleation. The model shows that both the size and the free-energy effect of individual AFPs and their surface coverage dominate the antifreeze capacity of an AFP-coated surface. The simulation results are consistent with the recent experiments qualitatively, revealing the origin of the surprisingly low antifreeze capacity of an AFP-coated surface when the coverage is not particularly high as shown in experiment. These results will hopefully deepen our understanding of the antifreeze effects and thus be potentially useful for designing novel antifreeze coating materials based on biomolecules.

Interfacial Water Arrangement in the Ice-Bound State of an Antifreeze Protein: A Molecular Dynamics Simulation Study

Molecular dynamics (MD) simulations have been carried out to study the heterogeneous ice nucleation on modeled peptide surfaces. Simulations show that large peptide surfaces made by TxT (threonine-x-threonine) motifs with the arrangements of threonine (Thr) residues identical to the periodic arrangements of waters on either the basal or prism plane of ice are capable of ice nucleation. Nucleated ice plane is the (0001) basal plane of hexagonal ice (Ih) or (111) plane of cubic ice (Ic). However, due to predefined simulation cell dimensions, the ice growth is only observed on the surface where the Thr residues are arranged like the water arrangement on the basal plane of ice Ih. The γ-methyl and γ-hydroxyl groups of Thr residue are necessary for such ice formation. From this ice nucleation and growth simulation, the interfacial water arrangement in the ice-bound state of Tenebrio molitor antifreeze protein (TmAFP) has been determined. The interfacial water arrangement in the ice-bound state of TmAFP is characterized by five-membered hydrogen bonded rings, where each of the hydroxyl groups of the Thr residues on the ice-binding surface (IBS) of the protein is a ring member. It is found that the water arrangement at the protein-ice interface is distorted from that in bulk ice. Our analysis further reveals that the hydroxyl groups of Thr residues on the IBS of TmAFP form maximum three hydrogen bonds each with the waters in the bound state and methyl groups of Thr residues occupy wider spaces than the normal grooves on the (111) plane of ice Ic. Methyl groups are also located above and along the 3-fold rotational axes of the chair-formed hexagonal hydrogen bonded water rings on the (111) plane.

Anisotropy in geometrically rough structure of ice prismatic plane interface during growth: Development of a modified six-site model of H2O and a molecular dynamics simulation

This paper presents a modified version of the six-site model of H₂O [H. Nada and J. P. J. M. van der Eerden, J. Chem. Phys. 118, 7401 (2003)]. Although the original six-site model was optimized by assuming the cut-off of the Coulomb interaction at an intermolecular distance of 10 Å, the modified model is optimized by using the Ewald method for estimating the Coulomb interaction. Molecular dynamics (MD) simulations of an ice-water interface suggest that the melting point of ice at 1 atm in the modified model is approximately 274.5 K, in good agreement with the real melting point of 273.15 K. MD simulations of bulk ice and water suggest that the modified model reproduces not only the structures and density curves of ice and water, but also the diffusion coefficient of water molecules in water near the melting point at 1 atm. Using the modified model, a large-scale MD simulation of the growth at an ice-water interface of the prismatic plane is performed to elucidate the anisotropy in the interface structure during growth. Simulation results indicate that the geometrical roughness of the ice growth front at the interface is greater in the c-axis direction than in the direction normal to the c-axis when it is analyzed along the axes parallel to the prismatic plane. In addition, during the growth at the interface, the transient appearance of specific crystallographic planes, such as a {202¯1} pyramidal plane, occurs preferentially at the ice growth front. The effect of different ensembles with different simulation systems on the anisotropy in the interface structure is also investigated.

Oscillations and accelerations of ice crystal growth rates in microgravity in presence of antifreeze glycoprotein impurity in supercooled water

The free growth of ice crystals in supercooled bulk water containing an impurity of glycoprotein, a bio-macromolecule that functions as ‘antifreeze’ in living organisms in a subzero environment, was observed under microgravity conditions on the International Space Station. We observed the acceleration and oscillation of the normal growth rates as a result of the interfacial adsorption of these protein molecules, which is a newly discovered impurity effect for crystal growth. As the convection caused by gravity may mitigate or modify this effect, secure observations of this effect were first made possible by continuous measurements of normal growth rates under long-term microgravity condition realized only in the spacecraft. Our findings will lead to a better understanding of a novel kinetic process for growth oscillation in relation to growth promotion due to the adsorption of protein molecules and will shed light on the role that crystal growth kinetics has in the onset of the mysterious antifreeze effect in living organisms, namely, how this protein may prevent fish freezing.

Unusual structural properties of water within the hydration shell of hyperactive antifreeze protein

Many hypotheses can be encountered explaining the mechanism of action of antifreeze proteins. One widespread theory postulates that the similarity of structural properties of solvation water of antifreeze proteins to ice is crucial to the antifreeze activity of these agents. In order to investigate this problem, the structural properties of solvation water of the hyperactive antifreeze protein from Choristoneura fumiferana were analyzed and compared with the properties of solvation water present at the surface of ice. The most striking observations concerned the temperature dependence of changes in water structure. In the case of solvation water of the ice-binding plane, the difference between the overall structural ordering of solvation water and bulk water diminished with increasing temperature; in the case of solvation water of the rest of the protein, the trend was opposite. In this respect, the solvation water of the ice-binding plane roughly resembled the hydration layer of ice. Simultaneously, the whole solvation shell of the protein displayed some features that are typical for solvation shells of many other proteins and are not encountered in the solvation water of ice. In the first place, this is an increase in density of water around the protein. The opposite is true for the solvation water of ice - it is less dense than bulk water. Therefore, even though the structure of solvation water of ice-binding plane and the structure of solvation water of ice seem to share some similarities, densitywise they differ.

The biological function of an insect antifreeze protein simulated by molecular dynamics

Antifreeze proteins (AFPs) protect certain cold-adapted organisms from freezing to death by selectively adsorbing to internal ice crystals and inhibiting ice propagation. The molecular details of AFP adsorption-inhibition is uncertain but is proposed to involve the Gibbs–Thomson effect. Here we show by using unbiased molecular dynamics simulations a protein structure-function mechanism for the spruce budworm Choristoneura fumiferana AFP, including stereo-specific binding and consequential melting and freezing inhibition. The protein binds indirectly to the prism ice face through a linear array of ordered water molecules that are structurally distinct from the ice. Mutation of the ice binding surface disrupts water-ordering and abolishes activity. The adsorption is virtually irreversible, and we confirm the ice growth inhibition is consistent with the Gibbs–Thomson law.

DOI:http://dx.doi.org/10.7554/eLife.05142.001

Water expands as it freezes. If this happens to the water inside plants and animals, the resulting ice crystals can rupture cells. To prevent this, many plants and animals that live in cold climates have evolved ‘antifreeze proteins’. When a small particle of ice first starts to form, the antifreeze proteins bind to it and prevent the water around it freezing, hence preventing the growth of an ice crystal.

There are many different types of antifreeze protein, and some are more active than others. For example, some insects including the spruce budworm are exposed to extremely cold temperatures—sometimes below −30°C—and these insects have antifreeze proteins that are highly active.

It is not fully understood how different antifreeze proteins interact with ice and prevent the growth of ice crystals. This is largely because, as yet, there are no experimental techniques that make it possible to see how antifreeze proteins and water molecules arrange themselves at the surface of a growing particle of ice. Instead, scientists have developed computer simulations to investigate this process. While many of these studies have provided valuable information, the computational methods used have only recently become powerful enough to analyze how the antifreeze proteins approach the surface of the ice particle.

Kuiper et al. carried out simulations involving a highly active antifreeze protein from the spruce budworm. The results of these simulations revealed that this antifreeze protein does not bind directly to ice; instead, water molecules at the surface of the protein act as a bridge between the protein and the ice. These water molecules are highly ordered and though they have similarities with how water is structured in the ice, they are distinct from the ice lattice itself. Furthermore, this arrangement appears to be important for allowing the spruce budworm antifreeze protein to interact with the ice.

This study provides detailed insights as to how a highly active antifreeze protein helps to prevent ice crystals forming. In the future, the computational simulations used here may be extended to study the dynamics of other antifreeze proteins, and also how crystals of other materials form.

DOI:http://dx.doi.org/10.7554/eLife.05142.002

Crystallization of confined water pools with radii greater than 1 nm in AOT reverse micelles

Freezing of water pools inside aerosol sodium bis(2-ethylhexyl) sulfosuccinate (AOT) reverse micelles has been investigated. Previous freezing experiments suffer from collision and fusion of AOT micelles and resultant loss of water from the water pool by shedding out during the cooling process. These phenomena have restricted the formation of ice to only when the radius of the water pool (Rw) is below 1 nm, and only amorphous ice has been observed. To overcome the size limitation, a combination of rapid cooling and a custom-made cell allowing thin sample loading is applied for instantaneous and homogeneous freezing. The freezing process is monitored with attenuated total reflection infrared spectroscopy (ATR-IR) measurements. A cooling rate of ca. -100 K/min and a sample thickness of ca. 50 μm overcomes the limitations mentioned above and allows the crystallization of water pools with larger radii (Rw > 1 nm). The corresponding ATR-IR spectra of the frozen water pools with Rw < 2.0 nm show similar features to the spectrum of metastable cubic ice (Ic). Further increase of the radius of the water pool (Rw > 2.0 nm), unfortunately, drastically decreased the integrated area of the ν(OH) band observed just after freezing, indicating the breakup of the micellar structure and shedding out of the water pool. In addition, it was revealed that Ic ice can also be formed in flexible organic self-assembled AOT reverse micelles for at least Rw ≤ ca. 2 nm, as well as in inorganic and solid materials with a pore radius of ca. 2 nm. The dependence of the phase transition temperature on the curvature of the reverse micelles is discussed from the viewpoint of the Gibbs-Thomson effect.

Crystal structure of an insect antifreeze protein and its implications for ice binding

Antifreeze proteins (AFPs) help some organisms resist freezing by binding to ice crystals and inhibiting their growth. The molecular basis for how these proteins recognize and bind ice is not well understood. The longhorn beetle Rhagium inquisitor can supercool to below -25 °C, in part by synthesizing the most potent antifreeze protein studied thus far (RiAFP). We report the crystal structure of the 13-kDa RiAFP, determined at 1.21 Å resolution using direct methods. The structure, which contains 1,914 nonhydrogen protein atoms in the asymmetric unit, is the largest determined ab initio without heavy atoms. It reveals a compressed β-solenoid fold in which the top and bottom sheets are held together by a silk-like interdigitation of short side chains. RiAFP is perhaps the most regular structure yet observed. It is a second independently evolved AFP type in beetles. The two beetle AFPs have in common an extremely flat ice-binding surface comprising regular outward-projecting parallel arrays of threonine residues. The more active, wider RiAFP has four (rather than two) of these arrays between which the crystal structure shows the presence of ice-like waters. Molecular dynamics simulations independently reproduce the locations of these ordered crystallographic waters and predict additional waters that together provide an extensive view of the AFP interaction with ice. By matching several planes of hexagonal ice, these waters may help freeze the AFP to the ice surface, thus providing the molecular basis of ice binding.