Will It Be Beneficial To Simulate the Antifreeze Proteins at Ice Freezing Condition or at Lower Temperature?

Characteristic Structural Knowledge for Morphological Identification and Classification in Meso-Scale Simulations Using Principal Component Analysis

Meso-scale simulations have been widely used to probe aggregation caused by structural formation in macromolecular systems. However, the limitations of the long-length scale, resulting from its simulation box, cause difficulties in terms of morphological identification and insufficient classification. In this study, structural knowledge derived from meso-scale simulations based on parameters from atomistic simulations were analyzed in dissipative particle dynamic (DPD) simulations of PS-b-PI diblock copolymers. The radial distribution function and its Fourier-space counterpart or structure factor were proposed using principal component analysis (PCA) as key characteristics for morphological identification and classification. Disorder, discrete clusters, hexagonally packed cylinders, connected clusters, defected lamellae, lamellae and connected cylinders were effectively grouped.

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.

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.

Molecular dynamics study on the role of solvation water in the adsorption of hyperactive AFP to the ice surface

Using computer simulations, the early stages of the adsorption of the CfAFP molecule to the ice surface were analyzed. We found that the ice and the protein interact at least as early as when the protein is about 1 nm away from the ice surface. These interactions are mediated by interfacial solvation water and are possible thanks to the structural ordering of the solvent. This ordering leads to positional preference of the protein relative to the ice crystal before the final attachment to the ice surface takes place, accompanied by the solidification of the interfacial water. It is possible because the solvation water of the ice-binding plane of CfAFP is susceptible to the overlapping with the solvation water of ice and is mostly changeable into ice itself. These remote interactions significantly increase efficacy of the adsorption process by facilitating the geometric adjustment of the active region of the CfAFP molecule to the ice surface. Because of the ordered nature of the water molecules at the ice-binding plane, the energy of their interactions with the ice-binding surface of the protein does not change upon the ongoing solidification of solvation water. However, the structure of the solvation water is not strictly ice-like and the growth of ice in the interfacial water is not initiated at the side of the protein. On the contrary, we find that solvation water of CfAFP solidifies slower than solvation water of ice - the solidification of interfacial water starts at the surface of ice. Moreover, in the presence of the CfAFP molecule, also solvation water of ice solidifies slower compared to the situation when the protein is not present next to the ice surface. Additionally, the presence of the protein molecule shifts the ratio of cubic to hexagonal ice that spontaneously forms at the ice surface, by introducing another layer of ordered water molecules - opposite to the ice lattice, at the other side of the crystallizing layer of water.

Molecular dynamics studies show solvation structure of type III antifreeze protein is disrupted at low pH

Antifreeze proteins are a class of biological molecules of interest in many research and industrial applications due to their highly specialized function, but there is little information of their stability and properties under varied pH derived from computational studies. To gain novel insights in this area, we conducted molecular dynamics (MD) simulations with the antifreeze protein 1KDF at varied temperatures and pH. Water solvation and H-bond formation around specific residues - ASN14, THR18 and GLN44 - involved in its antifreeze activity were extensively studied. We found that at pH1 there was a disruption in water solvation around the basal and the ice binding surfaces of the molecule. This was induced by a small change in the secondary structure propensities of some titrable residues, particularly GLU35. This change explains the experimentally observed reduction in antifreeze activity previously reported for this protein at pH1. We also found that THR18 showed extremely low H-bond formation, and that the three antifreeze residues all had very low average H-bond lifetimes. Our results confirm long-standing assumptions that these small, compact molecules can maintain their antifreeze activity in a wide range of pH, while demonstrating the mechanism that may reduce antifreeze activity at low pH. This aspect is useful when considering industrial and commercial use of antifreeze proteins subject to extreme pH environments, in particular in food industrial applications.

Structure of solvation water around the active and inactive regions of a type III antifreeze protein and its mutants of lowered activity

Water molecules from the solvation shell of the ice-binding surface are considered important for the antifreeze proteins to perform their function properly. Herein, we discuss the problem whether the extent of changes of the mean properties of solvation water can be connected with the antifreeze activity of the protein. To this aim, the structure of solvation water of a type III antifreeze protein from Macrozoarces americanus (eel pout) is investigated. A wild type of the protein is used, along with its three mutants, with antifreeze activities equal to 54% or 10% of the activity of the native form. The solvation water of the ice-binding surface and the rest of the protein are analyzed separately. To characterize the structure of solvation shell, parameters describing radial and angular characteristics of the mutual arrangement of the molecules were employed. They take into account short-distance (first hydration shell) or long-distance (two solvation shells) effects. The obtained results and the comparison with the results obtained previously for a hyperactive antifreeze protein from Choristoneura fumiferana lead to the conclusion that the structure and amino acid composition of the active region of the protein evolved to achieve two goals. The first one is the modification of the properties of the solvation water. The second one is the geometrical adjustment of the protein surface to the specific crystallographic plane of ice. Both of these goals have to be achieved simultaneously in order for the protein to perform its function properly. However, they seem to be independent from one another in a sense that very small antifreeze activity does not imply that properties of water become different from the ones observed for the wild type. The proteins with significantly lower activity still modify the mean properties of solvation water in a right direction, in spite of the fact that the accuracy of the geometrical match with the ice lattice is lost because of the mutations. Therefore, we do not observe any correlation between the antifreeze activity and the extent of modification of the properties of solvation water.

Investigation of changes in structure and thermodynamic of spruce budworm antifreeze protein under subfreezing temperature

The aim of this theoretical work is to investigate of the changes in structure and thermodynamics of spruce budworm antifreeze protein (sbAFP) at low temperatures by using molecular dynamics simulation. The aqueous solution will form ice crystal network under the vaguely hexagonal shape at low temperature and fully represented the characteristics of hydrophobic interaction. Like ice crystal network, the cyclohexane region (including cyclohexane molecules) have enough of the characteristics of hydrophobic interaction. Therefore, in this research the cyclohexane region will be used as a representation of ice crystal network to investigate the interactions of sbAFP and ice crystal network at low temperature. The activity of sbAFP in subfreezing environment, therefore, can be clearly observed via the changes of the hydrophobic (cyclohexane region) and hydrophilic (water region) interactions. The obtained results from total energies, hydrogen bond lifetime correlation C(t), radial distribution function, mean square deviation and snapshots of sbAFP complexes indicated that sbAFP has some special changes in structure and interaction with water and cyclohexane regions at 278 K, as being transition temperature point of water molecules in sbAFP complex at low temperatures, which is more structured and support the experimental observation that the sbAFP complex becomes more rigid as the temperature is lowered.

Exploring the Effects of Subfreezing Temperature and Salt Concentration on Ice Growth Inhibition of Antarctic Gram-Negative Bacterium Marinomonas Primoryensis Using Coarse-Grained Simulation

The aim of this work is to study the freezing process of water molecules surrounding Antarctic Gram-negative bacterium Marinomonas primoryensis antifreeze protein (MpAFP) and the MpAFP interactions to the surface of ice crystals under various marine environments (at different NaCl concentrations of 0.3, 0.6, and 0.8 mol/l). Our result indicates that activating temperature region of MpAFPs reduced as NaCl concentration increased. Specifically, MpAFP was activated and functioned at 0.6 mol/l with temperatures equal or larger 278 K, and at 0.8 mol/l with temperatures equal or larger 270 K. Additionally, MpAFP was inhibited by ice crystal network from 268 to 274 K and solid-liquid hybrid from 276 to 282 K at 0.3 mol/l concentration. Our results shed lights on structural dynamics of MpAFP among different marine environments.