Molecular Insights into the Unusual Structure of an Antifreeze Protein with a Hydrated Core

Preferential Ordering and Organization of Hydration Water Favor Nucleation of Ice by Ice-Nucleating Proteins over Antifreeze Proteins

Bacteria containing ice-nucleating proteins (INPs) evolved in nature to nucleate ice at the high sub-zero ambiance. The ability of the INPs to induce order in the hydration layer and their aggregation propensity appear to be key factors of their ice nucleation abilities. However, the mechanism of the process of ice nucleation by INPs is yet to be understood clearly. Here, we have performed all-atom molecular dynamics simulations and analyzed the structure and dynamics of the hydration layer around the proposed ice-nucleating surface of a model INP. Results are compared with the hydration of a topologically similar non-ice-binding protein (non-IBP) and another ice-growth inhibitory antifreeze protein (sbwAFP). We observed that the hydration structure around the ice-nucleating surface of INP is highly ordered and the dynamics of the hydration water are slower, compared to the non-IBP. Even the ordering of the hydration layer is more evident around the ice-binding surface of INP, compared to the antifreeze protein sbwAFP. Particularly with increasing repeat units of INP, we observe an increased population of ice-like water. Interestingly, the distances between the hydroxyl groups of the threonine ladder and its associated channel water of the ice-binding surface (IBS) of INP in the X and Y direction mimic the oxygen atom distances of the basal plane of hexagonal ice. However, the structural synergies between the hydroxyl group distances of the threonine ladder and its associated channel water of the IBS of sbwAFP and oxygen atom distances of the basal plane are less evident. This difference makes the IBS of the INP a better template for ice nucleation than AFP, although both of them bind to the ice surface efficiently.

Molecular Factors of Ice Growth Inhibition for Hyperactive and Globular Antifreeze Proteins: Insights from Molecular Dynamics Simulation

The molecular mechanism behind the ice growth inhibition by antifreeze proteins (AFPs) is yet to be understood completely. Also, what physical parameters differentiate between the AFP and non-AFP are largely unknown. Thus, to get an atomistic overview of the differential antifreeze activities of different classes of AFPs, we have studied ice growth from different ice surfaces in the presence of a moderately active globular type III AFP and a hyperactive spruce budworm (sbw) AFP. Results are compared with the observations of ice growth simulations in the presence of topologically similar non-AFPs using all-atom molecular dynamics simulations. Simulation data suggest that the ice surface coverage is a critical factor in ice growth inhibition. Due to the presence of an ice binding surface (IBS), AFPs form a high affinity complex with ice, accompanied by a transition of hydration water around the IBS from clathrate-like to ice-like. Several residues around the periphery of the IBS anchor the AFP to the curved ice surface mediated by multiple strong hydrogen bonds, stabilizing the complex immensely. In the high surface coverage regime, the slow unbinding kinetics dominates over the ice growth kinetics and thus facilitates the ice growth inhibition. Due to the non-availability of a proper IBS, non-AFPs form a low-affinity complex with the growing ice surface. As a result, the non-AFPs are continuously repelled by the surface. If the concentration of AFPs is low, then the effective surface coverage is reduced significantly. In this low surface coverage regime, AFPs can also behave like impurities and are engulfed by the growing ice crystal.

Structural characteristics of low-density environments in liquid water

The existence of two structural forms in liquid water has been a point of discussion for a long time. A phase transition between these two forms of liquid water has been proposed based on evidence from molecular simulations, and experiments have also been very recently able to track the proposed transition of the low-density liquid form to the high-density liquid form. We propose to use the average angle an oxygen atom makes with its neighbors to describe the structural environment of a water molecule. The distribution of this order parameter is observed to have two peaks with one peak at ∼109.5^{∘}, corresponding to the internal angle of a regular tetrahedron, indicating tetrahedral arrangement. The other peak corresponds to an environment with a tighter arrangement of neighboring molecules. The distribution of O-O-O angles is decomposed into two skewed distributions to estimate the fractions of the two liquid forms in water. A good similarity is observed between the temperature and pressure trends of fractions of locally favored tetrahedral structure (LFTS) form estimated using the new order parameter and the reports in the literature, over a range of temperatures and pressures. We also compare the structural environments indicated by different order parameters and find that the order parameter proposed in this paper captures the structure of first solvation shell of the LFTS accurately.

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.

Deciphering the Role of the Non-ice-binding Surface in the Antifreeze Activity of Hyperactive Antifreeze Proteins

Antifreeze proteins (AFPs) show thermal hysteresis through specific interaction with the ice crystal. Hyperactive AFPs interact with the ice surface through a threonine-rich motif present at their ice-binding surface (IBS). Ordering of water around the IBS was extensively investigated. However, the role of non-IBS in ice growth inhibition is yet to be understood completely. The present study explores the nature of hydration and its length-scale evaluation around the non-IBS for hyperactive AFPs. We observed that the hydration layer of non-IBS is liquid-like, even in highly supercooled conditions, and the nature of hydration is drastically different from the hydration pattern of non-AFP surfaces. In similar conditions, the hydration layer around the IBS is ice-like ordered. Non-IBS of the hyperactive AFP exposes toward the bulk and is able to maintain the liquid-like character of its hydration water up to 15 Å. We also find that the amino acid compositions and their spatial distribution on the non-IBS are markedly different from those of the IBS and non-AFP surfaces. These results elucidate the combined role of IBS and non-IBS in ice-growth inhibition. While IBS is required to adsorb on ice efficiently, the exposed non-IBS may prevent ice nucleation/growth on top of the bound AFPs.

Factors Promoting the Formation of Clathrate-Like Ordering of Water in Biomolecular Structure at Ambient Temperature and Pressure

Clathrate hydrate forms when a hydrophobic molecule is entrapped inside a water cage or cavity. Although biomolecular structures also have hydrophobic patches, clathrate-like water is found in only a limited number of biomolecules. Also, while clathrate hydrates form at low temperature and moderately higher pressure, clathrate-like water is observed in biomolecular structure at ambient temperature and pressure. These indicate presence of other factors along with hydrophobic environment behind the formation of clathrate-like water in biomolecules. In the current study, we presented a systematic approach to explore the factors behind the formation of clathrate-like water in biomolecules by means of molecular dynamics simulation of a model protein, maxi, which is a naturally occurring nanopore and has clathrate-like water inside the pore. Removal of either confinement or hydrophobic environment results in the disappearance of clathrate-like water ordering, indicating a coupled role of these two factors. Apart from these two factors, clathrate-like water ordering also requires anchoring groups that can stabilize the clathrate-like water through hydrogen bonding. Our results uncover crucial factors for the stabilization of clathrate-like ordering in biomolecular structure which can be used for the development of new biomolecular structure promoting clathrate formation.

Ordered hydration layer mediated ice adsorption of a globular antifreeze protein: mechanistic insight

The ice binding surface of a type III AFP induces water ordering at lower temperature, which mediates its adsorption on the ice surface.

Influence of glycerol on the cooling effect of pair hydrophobicity in water: relevance to proteins’ stabilization at low temperature

Glycerol reduces the cooling effect of pair hydrophobicity (reduction of hydrophobicity with decreasing temperature) in water.