Local order, energy, and mobility of water molecules in the hydration shell of small peptides

Structure, energetics and dynamics in crowded amino acid solutions: a molecular dynamics study

A comprehensive understanding of crowding effects on biomolecular processes necessitates investigating the bulk thermodynamic and kinetic properties of the solutions with an accurate molecular representation of the crowded milieu. Recent studies have reparameterized the non-bonded dispersion interaction of solutes to precisely model intermolecular interactions, which would circumvent artificial aggregation as shown by the original force-fields. However, the performance of this reparameterization is yet to be assessed for concentrated crowded solutions in terms of investigating the hydration shell structure, energetics and dynamics. In this study, we perform molecular dynamics simulations of crowded aqueous solutions of five zwitterionic neutral amino acids (Gly, Ala, Thr, Pro, and Ser), mimicking the molecular crowding environment, using a modified AMBER ff99SB-ILDN force-field. We systematically examine and show that the reproducibility of the osmotic coefficients, density, viscosity and self-diffusivity of amino acids improves using the modified force-field in crowded concentrations. The modified force-field also improves the structuring of the solute solvation shells, solute interaction energy and convergence of tails of radial distribution functions, indicating reduction in the artificial aggregation. Our results also indicate that the hydrogen bonding network of water weakens and water molecules anomalously diffuse at small time scales in the crowded solutions. These results underscore the significance of examining the solution properties and anomalous hydration behaviour of water in crowded solutions, which have implications in shaping the structure and dynamics of biomolecules. The findings also illustrate the improvement in predicting bulk solution properties using the modified force-field, thereby providing an approach towards accurate modeling of crowded molecular solutions.

Flexibility of the Binding Regions of a Protein-DNA Complex and the Structure and Ordering of Interfacial Water

Noncovalent interactions between protein and DNA are important to comprehend different biological activities in living organisms. One important issue is how the protein identifies the target DNA and the influence of the resulting protein-DNA complex on the hydration environment around it. In this study, we have carried out atomistic molecular dynamics simulations of the protein-DNA complex formed by the dimeric form of the α-helical N-terminal domain of the λ-repressor protein with its operator DNA. Local heterogeneous flexibilities of the residues of the protein and the DNA components that are involved in binding and the microscopic structure and ordering of water around those have been investigated in detail. The calculations revealed concurrent existence of highly ordered as well as disordered water molecules at the interface. It is found that a fraction of doubly coordinated water molecules exhibit high degree of ordering at the interface, while the randomly oriented ones are coordinated with three water molecules. The effect has been found to be more around the protein and DNA residues that are in contact in the complexed state. We believe that such highly ordered two-coordinated water molecules are likely to act as an adhesive to facilitate the formation of a protein-DNA complex and maintain its structural stability.

Decoding signatures of structure, bulk thermodynamics, and solvation in three-body angle distributions of rigid water models

A tetrahedral structure resulting from hydrogen bonding is a hallmark of liquid water and plays a significant role in determining its unique thermophysical properties. This water feature has helped understand anomalous properties and physically interpret and model hydrophobic solvation thermodynamics. Tetrahedrality is well described by the geometric relationship of any central water molecule with two of its nearest neighbors in the first coordination shell, as defined by the corresponding "three-body" angle. While order parameters and even full water models have been developed using specific or average features of the three-body angle distribution, here we examine the distribution holistically, tracking its response to changes in temperature, density, and the presence of model solutes. Surprisingly, we find that the three-body distribution responds by varying primarily along a single degree of freedom, suggesting a remarkably simplified view of water structure. We characterize three-body angle distributions across temperature and density space and identify principal components of the variations with state conditions. We show that these principal components embed physical significance and trace out transitions between tetrahedral and simple-fluid-like behavior. Moreover, we find that the ways three-body angles vary within the hydration shells of model colloids of different types and sizes are nearly identical to the variations seen in bulk water across density and temperature. Importantly, through the principal directions of these variations, we find that perturbations to the hydration-water distributions well predict the thermodynamics associated with colloid solvation, in particular, the relative entropy of this process that captures indirect, solvent-mediated contributions to the hydration free energy.

Molecular Order Affects Interfacial Water Structure and Temperature-Dependent Hydrophobic Interactions between Nonpolar Self-Assembled Monolayers

Understanding how material properties affect hydrophobic interactions-the water-mediated interactions that drive the association of nonpolar materials-is vital to the design of materials in contact with water. Conventionally, the magnitude of the hydrophobic interactions between extended interfaces is attributed to interfacial chemical properties, such as the amount of nonpolar solvent-exposed surface area. However, recent experiments have demonstrated that the hydrophobic interactions between uniformly nonpolar self-assembled monolayers (SAMs) also depend on molecular-level SAM order. In this work, we use atomistic molecular dynamics simulations to investigate the relationship between SAM order, water structure, and hydrophobic interactions to explain these experimental observations. The SAM-SAM hydrophobic interactions calculated from the simulations increase in magnitude as SAM order increases, matching experimental observations. We explain this trend by showing that the molecular-level order of the SAM impacts the nanoscale structure of interfacial water molecules, leading to an increase in water structure near disordered SAMs. These findings are consistent with a decrease in the solvation entropy of disordered SAMs, which is confirmed by measuring the temperature dependence of hydrophobic interactions using both simulations and experiments. This study elucidates how hydrophobic interactions can be influenced by an interfacial physical property, which may guide the design of synthetic materials with fine-tuned interfacial hydrophobicity.

Understanding the microscopic origin behind heterogeneous properties of water confined in and around Aβ17-42 protofilaments

Aggregation of amyloid beta (Aβ) peptides in the brain is responsible for one of the most devastating neurodegenerative diseases, namely, Alzheimer's disease. In this study, we have carried out atomistic molecular dynamics simulations to explore the effects of non-uniform structural distortions of Aβ_17-42 pre-fibrillar aggregates of different sizes on the microscopic structure and ordering of water molecules confined within their amphiphilic nanocores. The calculations revealed non-uniform peptide-water interactions resulting in simultaneous existence of both highly ordered and disordered water molecules within the spatially heterogeneous confined environment of the protofilament cores. It is found that the high degree of ordering originates from a sizable fraction of doubly coordinated core water molecules, while the randomly oriented ones are those that are coordinated with three neighbors in their first coordination shells. Furthermore, it is quantitatively demonstrated that relative fractions of these two types of water molecules are correlated with the protofilament core topology and the degree of confinement within that. It is proposed that the ordered core waters are likely to stabilize the Aβ protofilaments by screening the residue charges and favoring water-mediated salt bridge formations, while the randomly oriented ones can drive further growth of the protofilaments by being displaced easily during the docking of additional peptides. In that way, both types of core water molecules can play equally important roles in controlling the growth and stability of the Aβ-aggregates.

Entropy connects water structure and dynamics in protein hydration layer

The enzyme Candida Antarctica lipase B (CALB) serves here as a model for understanding connections among hydration layer dynamics, solvation shell structure, and protein surface structure. The structure and dynamics of water molecules in the hydration layer were characterized for regions of the CALB surface, divided around each α-helix, β-sheet, and loop structure. Heterogeneous hydration dynamics were observed around the surface of the enzyme, in line with spectroscopic observations of other proteins. Regional differences in the structure of the biomolecular hydration layer were found to be concomitant with variations in dynamics. In particular, it was seen that regions of higher density exhibit faster water dynamics. This is analogous to the behavior of bulk water, where dynamics (diffusion coefficients) are connected to water structure (density and tetrahedrality) by excess (or pair) entropy, detailed in the Rosenfeld scaling relationship. Additionally, effects of protein surface topology and hydrophobicity on water structure and dynamics were evaluated using multiregression analysis, showing that topology has a somewhat larger effect on hydration layer structure-dynamics. Concave and hydrophobic protein surfaces favor a less dense and more tetrahedral solvation layer, akin to a more ice-like structure, with slower dynamics. Results show that pairwise entropies of local hydration layers, calculated from regional radial distribution functions, scale logarithmically with local hydration dynamics. Thus, the Rosenfeld relationship describes the heterogeneous structure-dynamics of the hydration layer around the enzyme CALB. These findings raise the question of whether this may be a general principle for understanding the structure-dynamics of biomolecular solvation.

Origin of the blueshift of water molecules at interfaces of hydrophilic cyclic compounds

Water molecules at interfaces of materials exhibit enigmatic properties. A variety of spectroscopic studies have observed a high-frequency motion in these water molecules, represented by a blueshift, at both hydrophobic and hydrophilic interfaces. However, the molecular mechanism behind this blueshift has remained unclear. Using Raman spectroscopy and ab initio molecular dynamics simulations, we reveal the molecular mechanism of the blueshift of water molecules around six monosaccharide isomers. In the first hydration shell, we found weak hydrogen-bonded water molecules that cannot have a stable tetrahedral water network. In the water molecules, the vibrational state of the OH bond oriented toward the bulk solvent strongly contributes to the observed blueshift. Our work suggests that the blueshift in various solutions originates from the vibrational motions of these observed water molecules.

The sensitivity of folding free energy landscapes of trpzips to mutations in the hydrophobic core

The sensitivity of the stability of folded states and free energy landscapes to the differences in the hydrophobic content of the core residues has been studied for the set of 16-residue trpzips, namely, Trpzip4, Trpzip5 and Trpzip6.

On the coupling between the dynamics of protein and water

The solvation entropy of flexible protein regions is higher than that of rigid regions and contributes differently to the overall thermodynamic stability.

Conformational features of the Aβ42 peptide monomer and its interaction with the surrounding solvent

Heterogeneous conformational flexibility of the Aβ monomers has been found to be correlated with the corresponding non-uniform entropy gains.

Comparison of hydration behavior and conformational preferences of the Trp-cage mini-protein in different rigid-body water models

Trp-cage unfolds at different temperatures in different water models revealing the sensitivity of conformational order metrics to the choice of water models.

Estimating Translational and Orientational Entropies Using the k-Nearest Neighbors Algorithm

Inhomogeneous fluid solvation theory (IFST) and free energy perturbation (FEP) calculations were performed for a set of 20 solutes to compute the hydration free energies. We identify the weakness of histogram methods in computing the IFST hydration entropy by showing that previously employed histogram methods overestimate the translational and orientational entropies and thus underestimate their contribution to the free energy by a significant amount. Conversely, we demonstrate the accuracy of the k-nearest neighbors (KNN) algorithm in computing these translational and orientational entropies. Implementing the KNN algorithm within the IFST framework produces a powerful method that can be used to calculate free-energy changes for large perturbations. We introduce a new KNN approach to compute the total solute-water entropy with six degrees of freedom, as well as the translational and orientational contributions. However, results suggest that both the solute-water and water-water entropy terms are significant and must be included. When they are combined, the IFST and FEP hydration free energies are highly correlated, with an R(2) of 0.999 and a mean unsigned difference of 0.9 kcal/mol. IFST predictions are also highly correlated with experimental hydration free energies, with an R(2) of 0.997 and a mean unsigned error of 1.2 kcal/mol. In summary, the KNN algorithm is shown to yield accurate estimates of the combined translational-orientational entropy and the novel approach of combining distance metrics that is developed here could be extended to provide a powerful method for entropy estimation in numerous contexts.

Free Energy Landscapes of Alanine Oligopeptides in Rigid-Body and Hybrid Water Models

Replica exchange molecular dynamics is used to study the effect of different rigid-body (mTIP3P, TIP4P, SPC/E) and hybrid (H1.56, H3.00) water models on the conformational free energy landscape of the alanine oligopeptides (acAnme and acA5nme), in conjunction with the CHARMM22 force field. The free energy landscape is mapped out as a function of the Ramachandran angles. In addition, various secondary structure metrics, solvation shell properties, and the number of peptide-solvent hydrogen bonds are monitored. Alanine dipeptide is found to have similar free energy landscapes in different solvent models, an insensitivity which may be due to the absence of possibilities for forming i-(i + 4) or i-(i + 3) intrapeptide hydrogen bonds. The pentapeptide, acA5nme, where there are three intrapeptide backbone hydrogen bonds, shows a conformational free energy landscape with a much greater degree of sensitivity to the choice of solvent model, though the three rigid-body water models differ only quantitatively. The pentapeptide prefers nonhelical, non-native PPII and β-sheet populations as the solvent is changed from SPC/E to the less tetrahedral liquid (H1.56) to an LJ-like liquid (H3.00). The pentapeptide conformational order metrics indicate a preference for open, solvent-exposed, non-native structures in hybrid solvent models at all temperatures of study. The possible correlations between the properties of solvent models and secondary structure preferences of alanine oligopeptides are discussed, and the competition between intrapeptide, peptide-solvent, and solvent-solvent hydrogen bonding is shown to be crucial in the relative free energies of different conformers.

Characterization of the Local Structure in Liquid Water by Various Order Parameters

A wide range of geometric order parameters have been suggested to characterize the local structure of liquid water and its tetrahedral arrangement, but their respective merits have remained elusive. Here, we consider a series of popular order parameters and analyze molecular dynamics simulations of water, in the bulk and in the hydration shell of a hydrophobic solute, at 298 and 260 K. We show that these parameters are weakly correlated and probe different distortions, for example the angular versus radial disorders. We first combine these complementary descriptions to analyze the structural rearrangements leading to the density maximum in liquid water. Our results reveal no sign of a heterogeneous mixture and show that the density maximum arises from the depletion in interstitial water molecules upon cooling. In the hydration shell of the hydrophobic moiety of propanol, the order parameters suggest that the water local structure is similar to that in the bulk, with only a very weak depletion in ordered configurations, thus confirming the absence of any iceberg-type structure. Finally, we show that the main structural fluctuations that affect water reorientation dynamics in the bulk are angular distortions, which we explain by the jump hydrogen-bond exchange mechanism.

The non-uniform early structural response of globular proteins to cold denaturing conditions: a case study with Yfh1

The mechanism of cold denaturation in proteins is often incompletely understood due to limitations in accessing the denatured states at extremely low temperatures. Using atomistic molecular dynamics simulations, we have compared early (nanosecond timescale) structural and solvation properties of yeast frataxin (Yfh1) at its temperature of maximum stability, 292 K (Ts), and the experimentally observed temperature of complete unfolding, 268 K (Tc). Within the simulated timescales, discernible "global" level structural loss at Tc is correlated with a distinct increase in surface hydration. However, the hydration and the unfolding events do not occur uniformly over the entire protein surface, but are sensitive to local structural propensity and hydrophobicity. Calculated infrared absorption spectra in the amide-I region of the whole protein show a distinct red shift at Tc in comparison to Ts. Domain specific calculations of IR spectra indicate that the red shift primarily arises from the beta strands. This is commensurate with a marked increase in solvent accessible surface area per residue for the beta-sheets at Tc. Detailed analyses of structure and dynamics of hydration water around the hydrophobic residues of the beta-sheets show a more bulk water like behavior at Tc due to preferential disruption of the hydrophobic effects around these domains. Our results indicate that in this protein, the surface exposed beta-sheet domains are more susceptible to cold denaturing conditions, in qualitative agreement with solution NMR experimental results.

Molecular dynamics simulation of a single-stranded DNA with heterogeneous distribution of nucleobases in aqueous medium

The DNA metabolic processes often involve single-stranded DNA (ss-DNA) molecules as important intermediates. In the absence of base complementarity, ss-DNAs are more flexible and interact strongly with water in aqueous media. Ss-DNA-water interactions are expected to control the conformational flexibility of the DNA strand, which in turn should influence the properties of the surrounding water molecules. We have performed room temperature molecular dynamics simulation of an aqueous solution containing the ss-DNA dodecamer, 5'-CGCGAATTCGCG-3'. The conformational flexibility of the DNA strand and the microscopic structure and ordering of water molecules around it have been explored. The simulation reveals transformation of the initial base-stacked form of the ss-DNA to a fluctuating collapsed coil-like conformation with the formation of a few non-sequentially stacked base pairs. A preliminary analysis shows further collapse of the DNA conformation in presence of additional salt (NaCl) due to screening of negative charges along the backbone by excess cations. Additionally, higher packing of water molecules within a short distance from the DNA strand is found to be associated with realignment of water molecules by breaking their regular tetrahedral ordering.

Comparing distance metrics for rotation using the k-nearest neighbors algorithm for entropy estimation

Distance metrics facilitate a number of methods for statistical analysis. For statistical mechanical applications, it is useful to be able to compute the distance between two different orientations of a molecule. However, a number of distance metrics for rotation have been employed, and in this study, we consider different distance metrics and their utility in entropy estimation using the k-nearest neighbors (KNN) algorithm. This approach shows a number of advantages over entropy estimation using a histogram method, and the different approaches are assessed using uniform randomly generated data, biased randomly generated data, and data from a molecular dynamics (MD) simulation of bulk water. The results identify quaternion metrics as superior to a metric based on the Euler angles. However, it is demonstrated that samples from MD simulation must be independent for effective use of the KNN algorithm and this finding impacts any application to time series data.

Sensitivity of local hydration behaviour and conformational preferences of peptides to choice of water model

Secondary structural preferences of the beta-hairpin of the 2GB1 protein in the folded and unfolded ensembles are shown to be sensitive to the choice of water model.

Water and water-like liquids: relationships between structure, entropy and mobility

Liquids with very diverse underlying interactions share the thermodynamic and transport anomalies of water, including metalloids, ionic melts and mesoscopic fluids. The generic feature that characterises such water-like liquids is a density-driven shift in the nature of local order in the condensed phases. The key semiquantitative relationships between structural order, thermodynamics and transport that are necessary in order to map out the consequences of this common qualitative feature for liquid-state properties and phase transformations of such systems are reviewed here. The application of these ideas to understand and model tetrahedral liquids, especially water, is discussed and possible extensions to other complex fluids are considered.

Assessing the accuracy of inhomogeneous fluid solvation theory in predicting hydration free energies of simple solutes

Accurate prediction of hydration free energies is a key objective of any free energy method that is applied to modeling and understanding interactions in the aqueous phase. Inhomogeneous fluid solvation theory (IFST) is a statistical mechanical method for calculating solvation free energies by quantifying the effect of a solute acting as a perturbation to bulk water. IFST has found wide application in understanding hydration phenomena in biological systems, but quantitative applications have not been comprehensively assessed. In this study, we report the hydration free energies of six simple solutes calculated using IFST and independently using free energy perturbation (FEP). This facilitates a validation of IFST that is independent of the accuracy of the force field. The results demonstrate that IFST shows good agreement with FEP, with an R² coefficient of determination of 0.99 and a mean unsigned difference of 0.7 kcal/mol. However, sampling is a major issue that plagues IFST calculations and the results suggest that a histogram method may require prohibitively long simulations to achieve convergence of the entropies, for bin sizes which effectively capture the underlying probability distributions. Results also highlight the sensitivity of IFST to the reference interaction energy of a water molecule in bulk, with a difference of 0.01 kcal/mol changing the predicted hydration free energies by approximately 2.4 kcal/mol for the systems studied here. One of the major advantages of IFST over perturbation methods such as FEP is that the systems are spatially decomposed to consider the contribution of specific regions to the total solvation free energies. Visualizing these contributions can yield detailed insights into solvation thermodynamics. An insight from this work is the identification and explanation of regions with unfavorable free energy density relative to bulk water. These regions contribute unfavorably to the hydration free energy. Further work is necessary before IFST can be extended to yield accurate predictions of binding free energies, but the work presented here demonstrates its potential.

Importance of protein conformational motions and electrostatic anchoring sites on the dynamics and hydrogen bond properties of hydration water

The microscopic dynamic properties of water molecules present in the vicinity of a protein are expected to be sensitive to its local conformational motions and the presence of polar and charged groups at the surface capable of anchoring water molecules through hydrogen bonds. In this work, we attempt to understand such sensitivity by performing detailed molecular dynamics simulations of the globular protein barstar solvated in aqueous medium. Our calculations demonstrate that enhanced confinement at the protein surface on freezing its local motions leads to increasingly restricted water mobility with long residence times around the secondary structures. It is found that the inability of the surface water molecules to bind with the protein residues by hydrogen bonds in the absence of protein-water (PW) electrostatic interactions is compensated by enhanced water-water hydrogen bonds around the protein with uniform bulklike behaviors. Importantly, it is further noticed that in contrast to the PW hydrogen bond relaxation time scale, the kinetics of the breaking and formation of such bonds are not affected on freezing the protein's conformational motions.

Water and other tetrahedral liquids: order, anomalies and solvation

In order to understand the common features of tetrahedral liquids with water-like anomalies, the relationship between local order and anomalies has been studied using molecular dynamics simulations for three categories of such liquids: (a) atomistic rigid-body models for water (TIP4P, TIP4P/2005, mTIP3P, SPC/E), (b) ionic melts, BeF(2) (TRIM model) and SiO(2) (BKS potential) and (c) Stillinger-Weber liquids parametrized to model water (mW) and silicon. Rigid-body, atomistic models for water and the Stillinger-Weber liquids show a strong correlation between tetrahedral and pair correlation order and the temperature for the onset of the density anomaly is close to the melting temperature. In contrast, the ionic melts show weaker and more variable degrees of correlation between tetrahedral and pair correlation metrics, and the onset temperature for the density anomaly is more than twice the melting temperature. In the case of water, the relationship between water-like anomalies and solvation is studied by examining the hydration of spherical solutes (Na(+), Cl(-), Ar) in water models with different temperature regimes of anomalies (SPC/E, TIP4P and mTIP3P). For both ionic and nonpolar solutes, the local structure and energy of water molecules is essentially the same as in bulk water beyond the second-neighbour shell. The local order and binding energy of water molecules are not perturbed by the presence of a hydrophobic solute. In the case of ionic solutes, the perturbation is largely localized within the first hydration shell. The binding energies for the ions are strongly dependent on the water models and clearly indicate that the geometry of the partial charge distributions, and the associated multipole moments, play an important role. However the anomalous behaviour of the water network has been found to be unimportant for polar solvation.

Conformational fluctuations of a protein-DNA complex and the structure and ordering of water around it

Protein-DNA binding is an important process responsible for the regulation of genetic activities in living organisms. The most crucial issue in this problem is how the protein recognizes the DNA and identifies its target base sequences. Water molecules present around the protein and DNA are also expected to play an important role in mediating the recognition process and controlling the structure of the complex. We have performed atomistic molecular dynamics simulations of an aqueous solution of the protein-DNA complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. The conformational fluctuations of the protein and DNA and the microscopic structure and ordering of water around them in the complex have been explored. In agreement with experimental studies, the calculations reveal conformational immobilization of the terminal segments of the protein on complexation. Importantly, it is discovered that both structural adaptations of the protein and DNA, and the subsequent correlation between them to bind, contribute to the net entropy loss associated with the complex formation. Further, it is found that water molecules around the DNA are more structured with significantly higher density and ordering than that around the protein in the complex.

Conformational flexibility of a protein-carbohydrate complex and the structure and ordering of surrounding water

Protein-carbohydrate non-covalent interactions are important to understand various biological processes in living organisms. One of the important issues in protein-carbohydrate binding is how the protein identifies the target carbohydrate and recognizes its conformational features. Surrounding water molecules are expected to play a critical role not only in mediating the recognition process but also in maintaining the structure of the complex. We carried out atomistic molecular dynamics (MD) simulations of an aqueous solution of the protein-carbohydrate complex formed between the hyaluronan binding domain (HABD) of the murine Cd44 protein and the octasaccharide hyaluronan (HA(8)). The conformational flexibilities of the protein and the carbohydrate, and the microscopic structure and ordering of water molecules around them in the complexed form have been explored. It is revealed that the formation of the complex is associated with significant immobilization of the monosaccharide units of the carbohydrate moiety that are involved in binding. Further, reduction in water densities around the binding residues of the two molecules in the complex with respect to their free forms clearly demonstrated that the recognition between the protein and the carbohydrate is facilitated by removal of a fraction of water molecules from regions around the binding domains.

Local order and mobility of water molecules around ambivalent helices

Water on a protein surface plays a key role in determining the structure and dynamics of proteins. Compared to the properties of bulk water, many aspects of the structure and dynamics of the water surrounding the proteins are less understood. It is interesting therefore to explore how the properties of the water within the solvation shell around the peptide molecule depend on its specific secondary structure. In this work we investigate the orientational order and residence times of the water molecules to characterize the structure, energetics, and dynamics of the hydration shell water around ambivalent peptides. Ambivalent sequences are identical sequences which display multiple secondary structures in different proteins. Molecular dynamics simulations of representative proteins containing variable helix, variable nonhelix, and conserved helix are also used to explore the local structure and mobility of water molecules in their vicinity. The results, for the first time, depict a different water distribution pattern around the conserved and variable helices. The water molecules surrounding the helical segments in variable helices are found to possess a less locally ordered structure compared to those around their corresponding nonhelical counterparts and conserved helices. The long conserved helices exhibit extremely high local residence times compared to the helical conformations of the variable helices, whereas the residence times of the nonhelical conformations of the variable helices are comparable to those of the short conserved helices. This differential pattern of the structure and dynamics of water molecules in the vicinity of conserved/variable helices may lend valuable insights for understanding the role of solvent effects in determining sequence ambivalency and help in improving the accuracy of water models used in the simulations of proteins.

Comparison of Tetrahedral Order, Liquid State Anomalies, and Hydration Behavior of mTIP3P and TIP4P Water Models

The relationship between local tetrahedral order, tagged particle potential energy, and coordination number is studied for mTIP3P and TIP4P models of water in the bulk as well as in the neighborhood of a small peptide. The tendency of water molecules with different binding or tagged particle potential energies to occupy environments with different degrees of disorder can be effectively illustrated by constructing tetrahedral order distributions and corresponding entropy metrics conditional on restricted ranges of local binding energy. At the state point corresponding to the onset of the density anomaly, the correlation between tetrahedral entropy versus tagged potential energy is strong and virtually identical for mTIP3P and TIP4P. In TIP4P, this correlation is retained up to temperatures as high as 300 K, while it is lost by 250 K in mTIP3P. In the 250-300 K regime that is important for biomolecular simulations, mTIP3P behaves essentially as a simple liquid while TIP4P shows the density and related anomalies characteristic of water. We also study the number of water molecules, the tetrahedral order, and the tagged molecule potential energies for water molecules as a function of the distance from the peptide for the 16-residue β-hairpin fragment of 2GB1 in mTIP3P and TIP4P solvents. The hydration shell coordination profiles (n(r)) of the number of water molecules are almost identical in the two solvents, but the radial variation in the local energies and local order show significant differences. The residue-wise variation in the tagged potential energy of water molecules within the first hydration shell is qualitatively similar in the two models. A comparison of the tetrahedral order distributions of water molecules lying at different distances from the biomolecular solute shows that the perturbation in the local tetrahedral order distributions of the bulk solvent due to the presence of the solute is marginal. Thus, in the 250-300 K regime, the mTIP3P and TIP4P water models show qualitatively different behavior in terms of the relationship between tetrahedral order and local energy, but as solvents in the neighborhood of a biomolecular solute, the differences between the two models are only quantitative and not qualitative.

Structure of hydration water in proteins: a comparison of molecular dynamics simulations and database analysis

Hydration layer water molecules play important structural and functional roles in proteins. Despite being a critical component in biomolecular systems, characterizing the properties of hydration water poses a challenge for both experiments and simulations. In this context we investigate the local structure of hydration water molecules as a function of the distance from the protein and water molecules respectively in 188 high resolution protein structures and compare it with those obtained from molecular dynamics simulations. Tetrahedral order parameter of water in proteins calculated from previous and present simulation studies show that the potential of bulk water overestimates the average tetrahedral order parameter compared to those calculated from crystal structures. Hydration waters are found to be more ordered at a distance between the first and second solvation shell from the protein surface. The values of the order parameter decrease sharply when the water molecules are located very near or far away from the protein surface. At small water-water distance, the values of order parameter of water are very low. The average order parameter records a maximum value at a distance equivalent to the first solvation layer with respect to the water-water radial distribution and asymptotically approaches a constant value at large distances. Results from present analysis will help to get a better insight into structure of hydration water around proteins. The analysis will also help to improve the accuracy of water models on the protein surface.

Thermodynamic, diffusional, and structural anomalies in rigid-body water models

Structural, density, entropy, and diffusivity anomalies of the TIP4P/2005 model of water are mapped out over a wide range of densities and temperatures. The locus of temperatures of maximum density (TMD) for this model is very close to the experimental TMD locus for temperatures between 250 and 275 K. Four different water models (mTIP3P, TIP4P, TIP5P, and SPC/E) are compared with the TIP4P/2005 model in terms of their anomalous behavior. For all the water models, the density regimes for anomalous behavior are bounded by a low-density limit at around 0.85-0.90 g cm(-3) and a high-density limit at about 1.10-1.15 g cm(-3). The onset temperatures of the density anomaly in the various models show a much greater variation, ranging from 202 K for mTIP3P to 289 K for TIP5P. The order maps for the various water models are qualitatively very similar with the structurally anomalous regions almost superimposable in the q(tet)-τ plane. Comparison of the phase diagrams of water models with the region of liquid-state anomalies shows that the crystalline phases are much more sensitive to the choice of water models than the liquid state anomalies; for example, SPC/E and TIP4P/2005 show qualitatively similar liquid state anomalies but very different phase diagrams. The anomalies in the liquid in all the models occur at much lower pressures than those at which the melting line changes from negative to positive slope. The results in this study demonstrate several aspects of structure-entropy-diffusivity relationships of water models that can be compared with experiment and used to develop better atomistic and coarse-grained models for water.

Hydration properties of α-, β-, and γ-cyclodextrins from molecular dynamics simulations

Atomistic molecular dynamics (MD) simulations of α-, β-, and γ-cyclodextrins (ACD, BCD, and GCD) in aqueous solutions have been performed. Detailed analyses were carried out to compare the microscopic properties of water confined within the cavities of these macromolecules and in the hydration layers around them. It is noticed that reduced tetrahedral ordering of water in and around the CD molecules are associated with their restricted motions. Interestingly, unlike the translational motions, the rotational motions of cavity water molecules are found to be highly dependent on cavity dimensions. Additionally, it is found that severely hindered rotational motion of cavity water molecules is the origin of drastically restricted structural relaxation of hydrogen bonds involving those water molecules. It is demonstrated that the geometrical constraints within the cavities of the CD molecules enhance the rate of reformation of broken hydrogen bonds, thereby resulting in rapid establishment of the breaking and reformation equilibria for hydrogen bonds involving cavity water molecules.

Entropy, diffusivity and the energy landscape of a waterlike fluid

Molecular dynamics simulations and instantaneous normal mode (INM) analysis of a fluid with core-softened pair interactions and waterlike liquid-state anomalies are performed to obtain an understanding of the relationship between thermodynamics, transport properties, and the potential energy landscape. Rosenfeld scaling of diffusivities with the thermodynamic excess and pair correlation entropy is demonstrated for this model. The INM spectra are shown to carry information about the dynamical consequences of the interplay between length scales characteristic of anomalous fluids, such as bimodality of the real and imaginary branches of the frequency distribution. The INM spectral information is used to partition the liquid entropy into two contributions associated with the real and imaginary frequency modes; only the entropy contribution from the imaginary branch captures the nonmonotonic behavior of the excess entropy and diffusivity in the anomalous regime of the fluid.