What Controls the Limit of Supercooling and Superheating of Pinned Ice Surfaces?

On the engulfment of antifreeze proteins by ice

Antifreeze proteins (AFPs) are remarkable biomolecules that suppress ice formation at trace concentrations. To inhibit ice growth, AFPs must not only bind to ice crystals, but also resist engulfment by ice. The highest supercooling, [Formula: see text], for which AFPs are able to resist engulfment is widely believed to scale as the inverse of the separation, [Formula: see text], between bound AFPs, whereas its dependence on the molecular characteristics of the AFP remains poorly understood. By using specialized molecular simulations and interfacial thermodynamics, here, we show that in contrast with conventional wisdom, [Formula: see text] scales as [Formula: see text] and not as [Formula: see text]. We further show that [Formula: see text] is proportional to AFP size and that diverse naturally occurring AFPs are optimal at resisting engulfment by ice. By facilitating the development of AFP structure-function relationships, we hope that our findings will pave the way for the rational design of AFPs.

Ice Recrystallization Unveils the Binding Mechanism Operating at a Diffused Interface

Recrystallization of ice is a natural phenomenon that causes adverse effects in cryopreservation, agriculture, and in frozen food industry. It has long been recognized that ice recrystallization occurs through the Ostwald ripening and accretion processes. However, neither of these processes has been explored in microscopic detail by state-of-the-art experimental techniques. We carried out atomistic molecular dynamics (MD) simulations to explore ice recrystallization through the accretion process. Attempts have been made to elucidate the binding mechanism that is operating at the diffused ice-water interface. It is demonstrated that two ice crystals spontaneously recognize each other and bind together to form a large crystal in liquid water, resulting in ice recrystallization by accretion. Interestingly, the study reveals that the binding occurs due to the freezing of the interfacial water layer present between the two ice planes, even at a temperature above the melting point of the ice crystal. The synergistically enhanced ordering effect of two ice surfaces on the interfacial water leads to such freezing occurring during the binding process. However, proper crystallographic alignment is not necessarily required for the binding of the two crystals. Simulations have also been carried out to study the binding between an ice crystal and the model ice-binding surface (IBS) of an antifreeze protein above the melting point of the ice crystal. It is found that such binding at the IBS is accompanied by freezing of the interfacial water. This establishes that the synergetic ordering-driven freezing of interfacial water is a common binding mechanism at the diffused surfaces of ice crystals. We believe that this mechanism will provide a microscopic understanding of the process of recrystallization inhibition and thus help in designing suitable materials for potent applications in recrystallization inhibition.

Unidirectional Freezing of Polymer Solution Droplets

Ice templating provides a means of generating textures with a well-defined topography. Recent applications involve the freezing of water droplets, with or without colloids, on flat or textured surfaces. An interesting feature of water droplets freezing on a substrate is the formation of a pointy tip at a constant angle, regardless of the substrate temperature, surface energy, or droplet volume. Here, by adding the polymer to water, we demonstrate how to manipulate and even prevent the formation of such an icy tip. We find that the sharpness of the tip decreases with increasing polymer concentration until completely disappearing above the overlap concentration, while the total freezing time increases concomitantly. Building on these observations, we combined simple geometrical arguments with heat flux measurements to model and connect the spatial and temporal evolution of polymer droplets under unidirectional freezing. Together our results provide new ways to control the shape of frozen droplets for ice templating or microstructure fabrication, with applications in tissue engineering, separation membranes, and soft robotics.

Does freezing induce self-assembly of polymers? A molecular dynamics study

This work investigates the freezing-induced self-assembly (FISA) of polyvinyl alcohol (PVA) and PVA-like polymers using molecular dynamics simulations. In particular, the effect of the degree of supercooling, degree of polymerization, polymer type, and initial local concentration on the FISA was studied. It was found that the preeminent factor responsible for FISA is not the diffusion of the polymers away from the nucleating ice front, but the increase in the polymer's local concentration upon freezing of the solvent (water). At a higher degree of supercooling, the polymers are engulfed by the growing ice front, impeding their diffusion into the supercooled solution and finally inhibiting their self-assembly. Conversely, at a relatively lower degree of supercooling, the rate of diffusion of the polymers into the supercooled solution is higher, which increases their local concentration and results in FISA. FISA was also observed to depend on the polymer-solvent interactions. Strongly favorable solute-solvent interactions hinder the self-assembly, whereas unfavorable solute-solvent interactions promote the self-assembly. The polymer and aggregate morphology were investigated using the radius of gyration, end-to-end distance, and asphericity analysis. This study brings molecular insights into the quintessential factors governing self-assembly via freezing of the solvent, which is a novel self-assembly technique especially suitable for biomedical applications.

Engulfment Avalanches and Thermal Hysteresis for Antifreeze Proteins on Supercooled Ice

Antifreeze proteins (AFPs) bind to the ice-water surface and prevent ice growth at temperatures below 0 °C through a Gibbs-Thomson effect. Each adsorbed AFP creates a metastable depression on the surface that locally resists ice growth, until ice engulfs the AFP. We recently predicted the susceptibility to engulfment as a function of AFP size, distance between AFPs, and supercooling [ J. Chem. Phys. 2023, 158, 094501]. For an ensemble of AFPs adsorbed on the ice surface, the most isolated AFPs are the most susceptible, and when an isolated AFP gets engulfed, its former neighbors become more isolated and more susceptible to engulfment. Thus, an initial engulfment event can trigger an avalanche of subsequent engulfment events, leading to a sudden surge of unrestrained ice growth. This work develops a model to predict the supercooling at which the first engulfment event will occur for an ensemble of randomly distributed AFP pinning sites on an ice surface. Specifically, we formulate an inhomogeneous survival probability that accounts for the AFP coverage, the distribution of AFP neighbor distances, the resulting ensemble of engulfment rates, the ice surface area, and the cooling rate. We use the model to predict thermal hysteresis trends and compare with experimental data.

Free energy barriers for anti-freeze protein engulfment in ice: Effects of supercooling, footprint size, and spatial separation

Anti-freeze proteins (AFPs) protect organisms at freezing conditions by attaching to the ice surface and arresting its growth. Each adsorbed AFP locally pins the ice surface, resulting in a metastable dimple for which the interfacial forces counteract the driving force for growth. As supercooling increases, these metastable dimples become deeper, until metastability is lost in an engulfment event where the ice irreversibly swallows the AFP. Engulfment resembles nucleation in some respects, and this paper develops a model for the "critical profile" and free energy barrier for the engulfment process. Specifically, we variationally optimize the ice-water interface and estimate the free energy barrier as a function of the supercooling, the AFP footprint size, and the distance to neighboring AFPs on the ice surface. Finally, we use symbolic regression to derive a simple closed-form expression for the free energy barrier as a function of two physically interpretable, dimensionless parameters.

Oriented attachment kinetics for rod-like particles at a flat surface: Buffon's needle at the nanoscale

The adsorption of large rod-like molecules or crystallites on a flat crystal face, similar to Buffon's needle, requires the rods to "land," with their binding sites in precise orientational alignment with matching sites on the surface. An example is provided by long, helical antifreeze proteins (AFPs), which bind at specific facets and orientations on the ice surface. The alignment constraint for adsorption, in combination with the loss in orientational freedom as the molecule diffuses toward the surface, results in an entropic barrier that hinders the adsorption. Prior kinetic models do not factor in the complete geometry of the molecule, nor explicitly enforce orientational constraints for adsorption. Here, we develop a diffusion-controlled adsorption theory for AFP molecules binding at specific orientations to flat ice surfaces. We formulate the diffusion equation with relevant boundary conditions and present analytical solutions to the attachment rate constant. The resulting rate constant is a function of the length and aspect ratio of the AFP, the distance threshold associated with binding, and solvent conditions such as temperature and viscosity. These results and methods of calculation may also be useful for predicting the kinetics of crystal growth through oriented attachment.

Structural and Mechanical Properties of Konjac Glucomannan Gels and Influence of Freezing-Thawing Treatments on Them

Freezing has been widely used for long-term food preservation. However, freezing-thawing (FT) treatment usually influences the texture and structure of food gels such as konjac. For their texture control after FT treatment, it is important to clarify the structural change of food gels during the FT process. In this study, we investigated the aggregated structures of konjac glucomannan (GM) gels during the FT process using simultaneous synchrotron small-angle X-ray/wide-angle X-ray scattering (SAXS/WAXS) techniques. The FT treatment resulted in more crystallization of GM, and consequently, a large increase in compressive stress. In-situ SAXS/WAXS measurements revealed the following findings: on freezing, water molecules came out of the aggregated phase of GM and after the thawing, they came back into the aggregated phase, but the aggregated structure did not return to the one before the freezing; the gel network enhanced the inhomogeneity due to the growth of ice crystals during freezing. Furthermore, we examined the influence of additives such as polyvinyl (alcohol) (PVA) and antifreeze glycoprotein (AFGP) on the mechanical and structural properties of freeze-thawed GM gels. Although the addition of PVA and AFGP suppressed the crystallization of GM, it could not prevent the growth of ice crystals and the increase in the inhomogeneity of the gel network. As a result, the compressive stresses for freeze-thawed GM gels containing PVA or AFGP were significantly higher compared with those of GM gels without FT treatments, although they were lower than those of freeze-thawed GM gels. The findings of this study may be useful for not only the texture control of freeze-thawed foods but also the improvement of the mechanical performance of the biomaterials.

Signatures of sluggish dynamics and local structural ordering during ice nucleation

We investigate the microscopic pathway of spontaneous crystallization in the ST2 model of water under deeply supercooled conditions via unbiased classical molecular dynamics simulations. After quenching below the liquid–liquid critical point, the ST2 model spontaneously separates into low-density liquid (LDL) and high-density liquid phases, respectively. The LDL phase, which is characterized by lower molecular mobility and enhanced structural order, fosters the formation of a sub-critical ice nucleus that, after a stabilization time, develops into the critical nucleus and grows. Polymorphic selection coincides with the development of the sub-critical nucleus and favors the formation of cubic (Ic) over hexagonal (Ih) ice. We rationalize polymorphic selection in terms of geometric arguments based on differences in the symmetry of second neighbor shells of ice Ic and Ih, which are posited to favor formation of the former. The rapidly growing critical nucleus absorbs both Ic and Ih crystallites dispersed in the liquid phase, a crystal with stacking faults. Our results are consistent with, and expand upon, recent observations of non-classical nucleation pathways in several systems.

Ice Recrystallization Inhibition by Amino Acids: The Curious Case of Alpha- and Beta-Alanine

Extremophiles produce macromolecules which inhibit ice recrystallization, but there is increasing interest in discovering and developing small molecules that can modulate ice growth. Realizing their potential requires an understanding of how these molecules function at the atomistic level. Here, we report the discovery that the amino acid l-α-alanine demonstrates ice recrystallization inhibition (IRI) activity, functioning at 100 mM (∼10 mg/mL). We combined experimental assays with molecular simulations to investigate this IRI agent, drawing comparison to β-alanine, an isomer of l-α-alanine which displays no IRI activity. We found that the difference in the IRI activity of these molecules does not originate from their ice binding affinity, but from their capacity to (not) become overgrown, dictated by the degree of structural (in)compatibility within the growing ice lattice. These findings shed new light on the microscopic mechanisms of small molecule cryoprotectants, particularly in terms of their molecular structure and overgrowth by ice.

Diffusion Attachment Model for Long Helical Antifreeze Proteins to Ice

Some of the most potent antifreeze proteins (AFPs) are approximately rigid helical structures that bind with one side in contact with the ice surface at specific orientations. These AFPs take random orientations in solution; however, most orientations become sterically inaccessible as the AFP approaches the ice surface. The effect of these inaccessible orientations on the rate of adsorption of AFP to ice has never been explored. Here, we present a diffusion-controlled theory of adsorption kinetics that accounts for these orientational restrictions to predict a rate constant for adsorption (k_on, in m/s) as a function of the length and width of the AFP molecules. We find that k_on decreases with length and diameter of the AFP and is almost proportional to the inverse of the area of the binding surface. We demonstrate that the restricted orientations create an entropic barrier to AFP adsorption, which we compute to be approximately 7 k_BT for most AFPs and up to 9 k_BT for Maxi, the largest known AFP. We compare the entropic resistance 1/k_on to resistances for diffusion through boundary layers and across typical distances in the extracellular matrix and find that these entropic and diffusion resistances could become comparable in the small confined spaces of biological environments.

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.

Tilting Behavior of Lamellar Ice Tip during Unidirectional Freezing of Aqueous Solutions

Freezing of ice has been largely reported from many aspects, especially its complex pattern formation. Ice grown from liquid phase is usually characteristic of lamellar morphology that plays a significant role in various domains. However, tilted growth of ice via transition from coplanar to noncoplanar in directional solidification has been paid little attention in previous studies and there was a misleading explanation of the formation of tilted lamellar ice. Here, we in situ investigated the variations of tilting behavior of lamellar ice tips under different conditions within a single ice crystal of manipulated orientation via unidirectional freezing of aqueous solutions. It is found that tilted growth of ice tips is sensitive to pulling velocity and solute type. These experimental results reveal intrinsic tilted growth behavior of lamellar ice, which is suggested to enrich our understanding of pattern formation of ice in relevant physical processes.

Molecular Properties and Chemical Transformations Near Interfaces

The properties of bulk water and aqueous solutions are known to change in the vicinity of an interface and/or in a confined environment, including the thermodynamics of ion selectivity at interfaces, transition states and pathways of chemical reactions, and nucleation events and phase growth. Here we describe joint progress in identifying unifying concepts about how air, liquid, and solid interfaces can alter molecular properties and chemical reactivity compared to bulk water and multicomponent solutions. We also discuss progress made in interfacial chemistry through advancements in new theory, molecular simulation, and experiments.

Is Ice Nucleation by Organic Crystals Nonclassical? An Assessment of the Monolayer Hypothesis of Ice Nucleation

Potent ice nucleating organic crystals display an increase in nucleation efficiency with pressure and memory effect after pressurization that set them apart from inorganic nucleants. These characteristics were proposed to arise from an ordered water monolayer at the organic-water interface. It was interpreted that ordering of the monolayer is the limiting step for ice nucleation on organic crystals, rendering their mechanism of nucleation nonclassical. Despite the importance of organics in atmospheric ice nucleation, that explanation has never been investigated. Here we elucidate the structure of interfacial water and its role in ice nucleation at ambient pressure on phloroglucinol dihydrate, the paradigmatic example of outstanding ice nucleating organic crystal, using molecular simulations. The simulations confirm the existence of an interfacial monolayer that orders on cooling and becomes fully ordered upon ice formation. The monolayer does not resemble any ice face but seamlessly connects the distinct hydrogen-bonding orders of ice and the organic surface. Although large ordered patches develop in the monolayer before ice nucleates, we find that the critical step is the formation of the ice crystallite, indicating that the mechanism is classical. We predict that the fully ordered, crystalline monolayer nucleates ice above -2 °C and could be responsible for the exceptional ice nucleation by the organic crystal at high pressures. The lifetime of the fully ordered monolayer around 0 °C, however, is too short to account for the memory effect reported in the experiments. The latter could arise from an increase in the melting temperature of ice confined by strongly ice-binding surfaces.

Single Ice Crystal Growth with Controlled Orientation during Directional Freezing

Ice growth has attracted great attention for its capability of fabricating hierarchically porous microstructure. However, the formation of tilted lamellar microstructure during freezing needs to be reconsidered due to the limited control of ice orientation with respect to the thermal gradient during in situ observations, which can greatly enrich our insight into architectural control of porous biomaterials. This paper provides an in situ study of the solid/liquid interface morphology evolution of directionally solidified single crystal ice with its C-axis (optical axis) perpendicular to directions of both the thermal gradient and the incident light in poly(vinyl alcohol, PVA) solutions. Multifaceted morphology and V-shaped lamellar morphology were clearly observed in situ for the first time. Quantitative characterizations on lamellar spacing, tilt angle, and tip undercooling of lamellar ice platelets provide a clearer insight into the inherent ice growth habit in polymeric aqueous systems and are suggested to exert significant impact on future design and optimization in porous biomaterials.

Genetic Algorithm Approach for the Optimization of Protein Antifreeze Activity Using Molecular Simulations

Antifreeze proteins (AFPs) are of much interest for their ability to inhibit ice growth at low concentrations. In this work, we present a genetic algorithm for the in silico design of AFP mutants with improved antifreeze activity, measured as the predicted thermal hysteresis at a fixed concentration, ΔT_C. Central to the algorithm is our recently developed neural network method for predicting ΔT_C from molecular simulations [Kozuch et al., PNAS, 115, 13252 (2018)]. Applying the algorithm to three structurally diverse AFPs, wfAFP, rQAE, and RiAFP, we find that significantly improved mutants are discovered for rQAE and RiAFP. Testing of the optimized mutants shows an increase in ΔT_C of 0.572 ± 0.11 K (262 ± 50.6%) and 1.33 ± 0.14 K (39.9 ± 4.19%) over the native structures for rQAE and RiAFP, respectively. Structural analysis of the optimized mutants reveals that the algorithm is able to exploit two pathways for enhancing the predicted antifreeze activity of the mutants: (1) increasing the local order of surface waters by encouraging the formation of internal water channels in the protein and (2) increasing the total ice-binding area by improving the planar structure of the ice-binding surface. Additionally, analysis of all mutants explored by the algorithm reveals that a subset of residues, mainly nonpolar, are particularly helpful in improving antifreeze activity at the ice-binding surface.

Water Freezes at Near-Zero Temperatures Using Carbon Nanotube-Based Electrodes under Static Electric Fields

Although static electric fields have been effective in controlling ice nucleation, the highest freezing temperature (T_f) of water that can be achieved in an electric field (E) is still uncertain. We performed a systematic study of the effect of an electric field on water freezing by varying the thickness of a dielectric layer and the voltage across it in an electrowetting system. Results show that T_f first increases sharply with E and then reaches saturation at -3.5 °C after a critical value E of 6 × 10⁶ V/m. Using classical heterogeneous nucleation theory, it is revealed that this behavior is due to saturation in the contact angle of the ice embryo with the underlying substrate. Finally, we show that it is possible to overcome this freezing saturation by controlling the uniformity of the electric field using carbon nanotubes. We achieve a T_f of -0.6 °C using carbon nanotube-based electrodes with an E of 3 × 10⁷ V/m. This work sheds new light on the control of ice nucleation and has the potential to impact many applications ranging from food freezing to ice production.

Antifreeze proteins and homogeneous nucleation: On the physical determinants impeding ice crystal growth

Antifreeze proteins (AFPs) are biopolymers capable of interfering with ice growth. Their antifreeze action is commonly understood considering that the AFPs, by pinning the ice surface, force the crystal-liquid interface to bend forming an ice meniscus, causing an increase in the surface free energy and resulting in a decrease in the freezing point ΔT^max. Here, we present an extensive computational study for a model protein adsorbed on a TIP4P/Ice crystal, computing ΔT^max as a function of the average distance d between AFPs, with simulations spanning over 1 µs. First, we show that the lower the d, the larger the ΔT^max. Then, we find that the water-ice-protein contact angle along the line ΔT^max(d) is always larger than 0°, and we provide a theoretical interpretation. We compute the curvature radius of the stable solid-liquid interface at a given supercooling ΔT ≤ ΔT^max, connecting it with the critical ice nucleus at ΔT. Finally, we discuss the antifreeze capability of AFPs in terms of the protein-water and protein-ice interactions. Our findings establish a unified description of the AFPs in the contest of homogeneous ice nucleation, elucidating key aspects of the antifreeze mechanisms and paving the way for the design of novel ice-controlling materials.

Modelling of short synthetic antifreeze peptides: Insights into ice-pinning mechanism

Organisms living in icy environments produce antifreeze proteins to control ice growth and recrystallization. It has been proposed that these molecules pin the surface of ice crystals, thus inducing the formation of a curved surface that arrests crystal growth. Such proteins are very appealing for many potential applications in food industry, material science and cryoconservation of organs and tissues. Unfortunately, their structural complexity has seriously hampered their practical use, while efficient and accessible synthetic analogues are highly desirable. In this paper, we used molecular dynamics based techniques to model the interaction of three short antifreeze synthetic peptides with an ice surface. The employed protocols succeeded in reproducing the ice pinning action of antifreeze peptides and the consequent ice growth arrest, as well as in distinguishing between antifreeze and control peptides, for which no such effect was observed. Principal components analysis of peptides trajectories in different simulation settings permitted to highlight the main structural features associated to antifreeze activity. Modeling results are highly correlated with experimentally measured properties, and insights on ice-peptide interactions and on conformational patterns favoring antifreeze activity will prompt the design of new and improved antifreeze peptides.

Gas hydrates in sustainable chemistry

This review includes the current state of the art understanding and advances in technical developments about various fields of gas hydrates, which are combined with expert perspectives and analyses.

Peptidic Antifreeze Materials: Prospects and Challenges

Necessitated by the subzero temperatures and seasonal exposure to ice, various organisms have developed a remarkably effective means to survive the harsh climate of their natural habitats. Their ice-binding (glyco)proteins keep the nucleation and growth of ice crystals in check by recognizing and binding to specific ice crystal faces, which arrests further ice growth and inhibits ice recrystallization (IRI). Inspired by the success of this adaptive strategy, various approaches have been proposed over the past decades to engineer materials that harness these cryoprotective features. In this review we discuss the prospects and challenges associated with these advances focusing in particular on peptidic antifreeze materials both identical and akin to natural ice-binding proteins (IBPs). We address the latest advances in their design, synthesis, characterization and application in preservation of biologics and foods. Particular attention is devoted to insights in structure-activity relations culminating in the synthesis of de novo peptide analogues. These are sequences that resemble but are not identical to naturally occurring IBPs. We also draw attention to impactful developments in solid-phase peptide synthesis and 'greener' synthesis routes, which may aid to overcome one of the major bottlenecks in the translation of this technology: unavailability of large quantities of low-cost antifreeze materials with excellent IRI activity at (sub)micromolar concentrations.

Role of the Solvation Water in Remote Interactions of Hyperactive Antifreeze Proteins with the Surface of Ice

Most protein molecules do not adsorb onto ice, one of the exceptions being so-called antifreeze proteins. In this paper, we describe that there is a force pushing an antifreeze protein molecule away from the ice surface when it is not oriented with its ice-binding plane toward the ice and that this pushing force may be also present even when the protein is oriented with its ice-binding plane toward the ice. This force is absent only when certain specific distance criteria are met, regarding the surface of ice and the protein. It acts at early stages of adsorption, prior to the solidification of water between the ice and the protein molecule nearby. We propose the water-originating mechanism of the generation of this force and also the mechanism of remote attachment of an antifreeze molecule to the ice surface. In liquid water, there exist locally favored structures, ordered and of high specific volume. The presence of a protein molecule usually shifts the equilibrium that exists in liquid water toward increasing the number of high-density, disordered structures and diminishing the number of low-density structures. Creation of the locally favored structures may be hampered not only near the non-ice-binding surfaces but also between the ice surface and the protein surface, if the distance between these surfaces does not allow these structures to develop because the available space is not sufficient for their proper formation. This conclusion is supported by the analysis of the mean geometry of a single hydrogen bond, as well as of the hydrogen bond network in the solvation layer and a structural order parameter that characterizes the separation between the first and second solvation shells of a water molecule.

Mimicking the Ice Recrystallization Activity of Biological Antifreezes. When is a New Polymer "Active"?

Antifreeze proteins and ice-binding proteins have been discovered in a diverse range of extremophiles and have the ability to modulate the growth and formation of ice crystals. Considering the importance of cryoscience across transport, biomedicine, and climate science, there is significant interest in developing synthetic macromolecular mimics of antifreeze proteins, in particular to reproduce their property of ice recrystallization inhibition (IRI). This activity is a continuum rather than an "on/off" property and there may be multiple molecular mechanisms which give rise to differences in this observable property; the limiting concentrations for ice growth vary by more than a thousand between an antifreeze glycoprotein and poly(vinyl alcohol), for example. The aim of this article is to provide a concise comparison of a range of natural and synthetic materials that are known to have IRI, thus providing a guide to see if a new synthetic mimic is active or not, including emerging materials which are comparatively weak compared to antifreeze proteins, but may have technological importance. The link between activity and the mechanisms involving either ice binding or amphiphilicity is discussed and known materials assigned into classes based on this.

Ice-binding proteins and the applicability and limitations of the kinetic pinning model

Ice-binding proteins (IBPs) are unique molecules that bind to and are active on the interface between two phases of water: ice and liquid water. This property allows them to affect ice growth in multiple ways: shaping ice crystals, suppressing the freezing point, inhibiting recrystallization and promoting nucleation. Advances in the protein's production technologies make these proteins promising agents for medical applications among others. Here, we focus on a special class of IBPs that suppress freezing by causing thermal hysteresis (TH): antifreeze proteins (AFPs). The kinetic pinning model describes the dynamics of a growing ice face with proteins binding to it, which eventually slow it down to a halt. We use the kinetic pinning model, with some adjustments made, to study the TH dependence on the solution's concentration of AFPs by fitting the model to published experimental data. We find this model describes the activity of (moderate) type III AFPs well, but is inadequate for the (hyperactive) Tenebrio molitor AFPs. We also find the engulfment resistance to be a key parameter, which depends on the protein's size. Finally, we explain intuitively how TH depends on the seeding time of the ice crystal in the protein solution. Using this insight, we explain the discrepancy in TH measurements between different assays. 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'.

How Do Surfactants Control the Agglomeration of Clathrate Hydrates?

Clathrate hydrates can spontaneously form under typical conditions found in oil and gas pipelines. The agglomeration of clathrates into large solid masses plugs the pipelines, posing adverse safety, economic, and environmental threats. Surfactants are customarily used to prevent the aggregation of clathrate particles and their coalescence with water droplets. It is generally assumed that a large contact angle between the surfactant-covered clathrate and water is a key predictor of the antiagglomerant performance of the surfactant. Here we use molecular dynamic simulations to investigate the structure and dynamics of surfactant films at the clathrate-oil interface, and their impact on the contact angle and coalescence between water droplets and hydrate particles. In agreement with the experiments, the simulations predict that surfactant-covered clathrate-oil interfaces are oil wet but super-hydrophobic to water. Although the water contact angle determines the driving force for coalescence, we find that a large contact angle is not sufficient to predict good antiagglomerant performance of a surfactant. We conclude that the length of the surfactant molecules, the density of the interfacial film, and the strength of binding of its molecules to the clathrate surface are the main factors in preventing the coalescence and agglomeration of clathrate particles with water droplets in oil. Our analysis provides a molecular foundation to guide the molecular design of effective clathrate antiagglomerants.

Combined molecular dynamics and neural network method for predicting protein antifreeze activity

Antifreeze proteins (AFPs) are a diverse class of proteins that depress the kinetically observable freezing point of water. AFPs have been of scientific interest for decades, but the lack of an accurate model for predicting AFP activity has hindered the logical design of novel antifreeze systems. To address this, we perform molecular dynamics simulation for a collection of well-studied AFPs. By analyzing both the dynamic behavior of water near the protein surface and the geometric structure of the protein, we introduce a method that automatically detects the ice binding face of AFPs. From these data, we construct a simple neural network that is capable of quantitatively predicting experimentally observed thermal hysteresis from a trio of relevant physical variables. The model's accuracy is tested against data for 17 known AFPs and 5 non-AFP controls.

The Clathrate-Water Interface Is Oleophilic

The slow nucleation of clathrate hydrates is a central challenge for their use in the storage and transportation of natural gas. Molecules that strongly adsorb to the clathrate-water interface decrease the crystal-water surface tension, lowering the barrier for clathrate nucleation. Surfactants are widely used to promote the nucleation and growth of clathrate hydrates. It has been proposed that these amphiphilic molecules bind to the clathrate surface via hydrogen bonding. However, recent studies reveal that PVCap, an amphiphilic polymer, binds to clathrates through hydrophobic moieties. Here we use molecular dynamic simulations and theory to investigate the mode and strength of binding of surfactants to the clathrate-water interface and their effect on the nucleation rate. We find that the surfactants bind to the clathrate-water interface exclusively through their hydrophobic tails. The binding is strong, driven by the entropy of dehydration of the alkyl chain, as it penetrates empty cavities at the hydrate surface. The hydrophobic attraction of alkyl groups to the clathrate surface also results in strong adsorption of alkanes. We identify two regimes for the binding of surfactants as a function of their density at the hydrate surface, which we interpret to correspond to the two steps of the Langmuir adsorption isotherm observed in experiments. Our results indicate that hydrophobic attraction to the clathrate-water interface is key for the design of soluble additives that promote the nucleation of hydrates. We use the calculated adsorption coefficients to estimate the concentration of sodium dodecyl sulfate (SDS) required to reach nucleation rates for methane hydrate consistent with those measured in experiments. To our knowledge, this study is the first to quantify the effect of surfactant concentration in the nucleation rate of clathrate hydrates.

Ice-Nucleating and Antifreeze Proteins Recognize Ice through a Diversity of Anchored Clathrate and Ice-like Motifs

Cold-adapted organisms produce antifreeze and ice-nucleating proteins to prevent and promote ice formation. The crystal structure of hyperactive bacterial antifreeze protein (AFP) MpAFP suggests that this protein binds ice through an anchored clathrate motif. It is not known whether other hyperactive AFPs and ice-nucleating proteins (INPs) use the same motif to recognize or nucleate ice. Here we use molecular simulations to elucidate the ice-binding motifs of hyperactive insect AFPs and a model INP of Pseudomonas syringae. We find that insect AFPs recognize ice through anchored clathrate motifs distinct from that of MpAFP. By performing simulations of ice nucleation by PsINP, we identify two distinct ice-binding sites on opposite sides of the β-helix. The ice-nucleating sequences identified in the simulations agree with those previously proposed for the closely related INP of Pseudomonas borealis based on the structure of the protein. The simulations indicate that these sites have comparable ice-nucleating efficiency, but distinct binding motifs, controlled by the amino acid sequence: one is an anchored clathrate and the other ice-like. We conclude that anchored clathrate and ice-like motifs can be equally effective for binding proteins to ice and promoting ice nucleation.