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

Wet Conformation of Prion-Like Domain and Intimate Correlation of Hydration and Conformational Fluctuations

Proteins with prion-like domains (PLDs) are involved in neurodegeneration-associated aggregation and are prevalent in liquid-like membrane-less organelles. These PLDs contain amyloidogenic stretches but can maintain dynamic disordered conformations, even in the condensed phase. However, the molecular mechanism underlying such intricate conformational properties of PLDs remains elusive. Here we employed molecular dynamics simulations to investigate the conformational properties of a prototypical PLD system (i.e., FUS PLD). According to our simulation results, PLD adopts a wet collapsed conformation, wherein most residues maintain sufficient hydration with the abundance of internal water. These internal water molecules can rapidly exchange between the protein interior and the bulk, enabling intensive coupling of the entire protein with its hydration environment. The dynamic exchange of water molecules is intimately correlated to the overall conformational fluctuations of PLD. Furthermore, the abundance of dynamic internal water suppresses the formation of aggregation-prone ordered structures. These results collectively elucidate the crucial role of internal water in sustaining the dynamic disordered conformation of the PLD and inhibiting its aggregation propensity.

Interfacial water molecules contribute to antibody binding to the receptor-binding domain of SARS-CoV-2 spike protein

Antibodies that recognize the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), especially the neutralizing antibodies, carry great hope in the treatment and final elimination of COVID-19. Driven by a synchronized global effort, thousands of antibodies against the spike protein have been identified during the past two years, with the structural information available at atomistic detail for hundreds of these antibodies. We developed an improved molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method including explicitly treated interfacial water to calculate the binding free energy between representative antibodies and the receptor binding domain (RBD) domain of SARS-COV-2 spike proteins. We discovered that explicit treatment of water molecules located at the interface between RBD and antibody effectively improves the results for the WT and variants of concern (VOC) systems. Interfacial water molecules, together with surface and internal water molecules, behave drastically from bulk water and exert peculiar impacts on protein dynamics and energy, and thus warrant explicit treatment to complement implicit solvent models. Our results illustrate the importance of including interfacial water molecules to approach efficient and reliable prediction of binding free energy.Communicated by Ramaswamy H. Sarma.

Antibody CDR amino acids underlying the functionality of antibody repertoires in recognizing diverse protein antigens

Antibodies recognize protein antigens with exquisite specificity in a complex aqueous environment, where interfacial waters are an integral part of the antibody–protein complex interfaces. In this work, we elucidate, with computational analyses, the principles governing the antibodies’ specificity and affinity towards their cognate protein antigens in the presence of explicit interfacial waters. Experimentally, in four model antibody–protein complexes, we compared the contributions of the interaction types in antibody–protein antigen complex interfaces with the antibody variants selected from phage-displayed synthetic antibody libraries. Evidently, the specific interactions involving a subset of aromatic CDR (complementarity determining region) residues largely form the predominant determinant underlying the specificity of the antibody–protein complexes in nature. The interfacial direct/water-mediated hydrogen bonds accompanying the CDR aromatic interactions are optimized locally but contribute little in determining the epitope location. The results provide insights into the phenomenon that natural antibodies with limited sequence and structural variations in an antibody repertoire can recognize seemingly unlimited protein antigens. Our work suggests guidelines in designing functional artificial antibody repertoires with practical applications in developing novel antibody-based therapeutics and diagnostics for treating and preventing human diseases.

Solvent effects on the NMR shieldings of stacked DNA base pairs

Stacking effects are among the most important effects in DNA. We have recently studied their influence in fragments of DNA through the analysis of NMR magnetic shieldings, firstly in vacuo. As a continuation of this line of research we show here the influence of solvent effects on the shieldings through the application of both explicit and implicit models. We found that the explicit solvent model is more appropriate for consideration due to the results matching better in general with experiments, as well as providing clear knowledge of the electronic origin of the value of the shieldings. Our study is grounded on a recently developed theoretical model of our own, by which we are able to learn about the magnetic effects of given fragments of DNA molecules on selected base pairs. We use the shieldings of the atoms of a central base pair (guanine-cytosine) of a selected fragment of DNA molecules as descriptors of physical effects, like π-stacking and solvent effects. They can be taken separately and altogether. The effect of π-stacking is introduced through the addition of some pairs above and below of the central base pair, and now, the solvent effect is considered including a network of water molecules that consist of two solvation layers, which were fixed in the calculations performed in all fragments. We show that the solvent effects enhance the stacking effects on the magnetic shieldings of atoms that belong to the external N-H bonds. The net effect is of deshielding on both atoms. There is also a deshielding effect on the carbon atoms that belong to CO bonds, for which the oxygen atom has an explicit hydrogen bond (HB) with a solvent water molecule. Solvent effects are found to be no higher than a few percent of the total value of the shieldings (between 1% and 5%) for most atoms, although there are few for which such an effect can be higher. There is one nitrogen atom, the acceptor of the HB between guanine and cytosine, that is more highly shielded (around 15 ppm or 10%) when the explicit solvent is considered. In a similar manner, the most external nitrogen atom of cytosine and the hydrogen atom that is bonded to it are highly deshielded (around 10 ppm for nitrogen and around 3 ppm for hydrogen).

The catalytic mechanism of the mitochondrial methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2)

Methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2) is a new drug target that is expressed in cancer cells but not in normal adult cells, which provides an Achilles heel to selectively kill cancer cells. Despite the availability of crystal structures of MTHFD2 in the inhibitor- and cofactor-bound forms, key information is missing due to technical limitations, including (a) the location of absolutely required Mg2+ ion, and (b) the substrate-bound form of MTHFD2. Using computational modeling and simulations, we propose that two magnesium ions are present at the active site whereby (i) Arg233, Asp225, and two water molecules coordinate MgA2+, while MgA2+ together with Arg233 stabilize the inorganic phosphate (Pi); (ii) Asp168 and three water molecules coordinate MgB2+, and MgB2+ further stabilizes Pi by forming a hydrogen bond with two oxygens of Pi; (iii) Arg201 directly coordinates the Pi; and (iv) through three water-mediated interactions, Asp168 contributes to the positioning and stabilization of MgA2+, MgB2+ and Pi. Our computational study at the empirical valence bond level allowed us also to elucidate the detailed reaction mechanisms. We found that the dehydrogenase activity features a proton-coupled electron transfer with charge redistribution connected to the reorganization of the surrounding water molecules which further facilitates the subsequent cyclohydrolase activity. The cyclohydrolase activity then drives the hydration of the imidazoline ring and the ring opening in a concerted way. Furthermore, we have uncovered that two key residues, Ser197/Arg233, are important factors in determining the cofactor (NADP+/NAD+) preference of the dehydrogenase activity. Our work sheds new light on the structural and kinetic framework of MTHFD2, which will be helpful to design small molecule inhibitors that can be used for cancer treatment.

How anisotropic and isotropic atomic displacement parameters monitor protein covalent bonds rigidity: isotropic B-factors underestimate bond rigidity

Under the assumption that covalent bonds are rigid, it is possible to compare the estimations of rigidity based on anisotropic and isotropic B-factors. This is done by computing the difference of the mean-square displacements (Delta-u) of atoms A and Z along the covalent bond A–Z, which must be close to zero for a rigid bond. The analysis of a high-quality set of protein structures, refined at a resolution better than (or equal to) 0.8 Angstroms, showed that Delta-u is significantly close to zero when anisotropic B-factors are used, with an average 60% Delta-u reduction. This reduction is larger for larger B-factors and this suggests that care should be taken in data-mining procedures that involve isotropic B-factors, especially at lower resolution, when anisotropic B-factors cannot be determined and when the average B-factor increases.

Split the Charge Difference in Two! A Rule of Thumb for Adding Proper Amounts of Ions in MD Simulations

Despite the relevance of properly setting ionic concentrations in Molecular Dynamics (MD) simulations, methods or practical rules to set ionic strength are scarce and rarely documented. Based on a recently proposed thermodynamics method we provide an accurate rule of thumb to define the electrolytic content in simulation boxes. Extending the use of good practices in setting up MD systems is promptly needed to ensure reproducibility and consistency in molecular simulations.

Dynamics of TMAO and urea in the hydration shell of the protein SNase

Using all-atom molecular dynamics simulations of aqueous solutions of the globular protein SNase, the dynamic behavior of water molecules and cosolvents (trimethylamine-N-oxide (TMAO) and urea) in the hydration shell of the protein was studied for different solvent compositions. TMAO is a potent protein-stabilizing osmolyte, whereas urea is known to destabilize proteins. For molecules that are initially located in successive narrow layers at a given distance from the protein, the mean displacements and the distribution of displacements for short time intervals are calculated. For molecules that are initially located in solvation shells of a given thickness around the protein, the characteristic residence times in these shells are determined to characterize the dynamic behavior of the solvent molecules as a function of the distance to the protein. A combined consideration of these characteristics allows to reveal additional features of the dynamics of the cosolvents. It is shown that TMAO molecules leave the nearest vicinity of the protein faster than urea molecules, despite the fact that the mobility of TMAO molecules, measured by their mean displacements, is lower than that of urea. Moreover, we show that the rate of release of TMAO molecules from the hydration shell is lower in ternary (TMAO + urea + H₂O) solvent mixtures than in the binary ones. This is consistent with a recent observation that the fraction of TMAO near the protein decreases in the presence of urea. From the analysis of the decay of the number of particles initially located in the region of the first peak of the distribution function of solvent molecules around the protein, we estimated that about 20 water molecules and 6-7 urea molecules stay near the protein for more than 1000 ps.

The impact of cryosolution thermal contraction on proteins and protein crystals: volumes, conformation and order

Cryocooling of macromolecular crystals is commonly employed to limit radiation damage during X-ray diffraction data collection. However, cooling itself affects macromolecular conformation and often damages crystals via poorly understood processes. Here, the effects of cryosolution thermal contraction on macromolecular conformation and crystal order in crystals ranging from 32 to 67% solvent content are systematically investigated. It is found that the solution thermal contraction affects macromolecule configurations and volumes, unit-cell volumes, crystal packing and crystal order. The effects occur through not only thermal contraction, but also pressure caused by the mismatched contraction of cryosolvent and pores. Higher solvent-content crystals are more affected. In some cases the solvent contraction can be adjusted to reduce mosaicity and increase the strength of diffraction. Ice formation in some crystals is found to cause damage via a reduction in unit-cell volume, which is interpreted through solvent transport out of unit cells during cooling. The results point to more deductive approaches to cryoprotection optimization by adjusting the cryosolution composition to reduce thermal contraction-induced stresses in the crystal with cooling.

The spatial range of protein hydration

Proteins interact with their aqueous surroundings, thereby modifying the physical properties of the solvent. The extent of this perturbation has been investigated by numerous methods in the past half-century, but a consensus has still not emerged regarding the spatial range of the perturbation. To a large extent, the disparate views found in the current literature can be traced to the lack of a rigorous definition of the perturbation range. Stating that a particular solvent property differs from its bulk value at a certain distance from the protein is not particularly helpful since such findings depend on the sensitivity and precision of the technique used to probe the system. What is needed is a well-defined decay length, an intrinsic property of the protein in a dilute aqueous solution, that specifies the length scale on which a given physical property approaches its bulk-water value. Based on molecular dynamics simulations of four small globular proteins, we present such an analysis of the structural and dynamic properties of the hydrogen-bonded solvent network. The results demonstrate unequivocally that the solvent perturbation is short-ranged, with all investigated properties having exponential decay lengths of less than one hydration shell. The short range of the perturbation is a consequence of the high energy density of bulk water, rendering this solvent highly resistant to structural perturbations. The electric field from the protein, which under certain conditions can be long-ranged, induces a weak alignment of water dipoles, which, however, is merely the linear dielectric response of bulk water and, therefore, should not be thought of as a structural perturbation. By decomposing the first hydration shell into polarity-based subsets, we find that the hydration structure of the nonpolar parts of the protein surface is similar to that of small nonpolar solutes. For all four examined proteins, the mean number of water-water hydrogen bonds in the nonpolar subset is within 1% of the value in bulk water, suggesting that the fragmentation and topography of the nonpolar protein-water interface has evolved to minimize the propensity for protein aggregation by reducing the unfavorable free energy of hydrophobic hydration.

Paradoxical Effect of Trehalose on the Aggregation of α-Synuclein: Expedites Onset of Aggregation yet Reduces Fibril Load

Aggregation of α-synuclein is closely connected to the pathology of Parkinson's disease. The phenomenon involves multiple steps, commenced by partial misfolding and eventually leading to mature amyloid fibril formation. Trehalose, a widely accepted osmolyte, has been shown previously to inhibit aggregation of various globular proteins owing to its ability to prevent the initial unfolding of protein. In this study, we have examined if it behaves in a similar fashion with intrinsically disordered protein α-synuclein and possesses the potential to act as therapeutic agent against Parkinson's disease. It was observed experimentally that samples coincubated with trehalose fibrillate faster compared to the case in its absence. Molecular dynamics simulations suggested that this initial acceleration is manifestation of trehalose's tendency to perturb the conformational transitions between different conformers of monomeric protein. It stabilizes the aggregation prone "extended" conformer of α-synuclein, by binding to its exposed acidic residues of the C terminus. It also favors the β-rich oligomers once formed. Interestingly, the total fibrils formed are still promisingly less since it accelerates the competing pathway toward formation of amorphous aggregates.

Recent advances in the UltraScan SOlution MOdeller (US-SOMO) hydrodynamic and small-angle scattering data analysis and simulation suite

The UltraScan SOlution MOdeller (US-SOMO) is a comprehensive, public domain, open-source suite of computer programs centred on hydrodynamic modelling and small-angle scattering (SAS) data analysis and simulation. We describe here the advances that have been implemented since its last official release (#3087, 2017), which are available from release #3141 for Windows, Linux and Mac operating systems. A major effort has been the transition from the legacy Qt3 cross platform software development and user interface library to the modern Qt5 release. Apart from improved graphical support, this has allowed the direct implementation of the newest, almost two-orders of magnitude faster version of the ZENO hydrodynamic computation algorithm for all operating systems. Coupled with the SoMo-generated bead models with overlaps, ZENO provides the most accurate translational friction computations from atomic-level structures available (Rocco and Byron Eur Biophys J 44:417-431, 2015a), with computational times comparable with or faster than those of other methods. In addition, it has allowed us to introduce the direct representation of each atom in a structure as a (hydrated) bead, opening interesting new modelling possibilities. In the small-angle scattering (SAS) part of the suite, an indirect Fourier transform Bayesian algorithm has been implemented for the computation of the pairwise distance distribution function from SAS data. Finally, the SAS HPLC module, recently upgraded with improved baseline correction and Gaussian decomposition of not baseline-resolved peaks and with advanced statistical evaluation tools (Brookes et al. J Appl Cryst 49:1827-1841, 2016), now allows automatic top-peak frame selection and averaging.

Effects of temperature, concentration, and isomer on the hydration structure in monosaccharide solutions

Water-monosaccharide coupled interactions are essential for the function, stability, and dynamics of all glycans. Using molecular dynamics simulations, we investigated the effects of temperature, concentration, and monosaccharide isomer on the hydration structure and water dynamics in the hydration shell of monosaccharides in solution. We found that perturbations of the hydrogen-bond (H-bond) network in the first hydration shell around each monosaccharide molecule can be separated into two regions: one rich in water molecules with donor H-bonds (in the 2.4-2.8 Å region) and the other rich in water molecules with abundant acceptor H-bonds (in the 2.8-3.3 Å region). Moreover, we investigated the dependencies of clustering and conversion of the conformers of the monosaccharides on temperature and concentration. Increasing the concentration enhances monosaccharide clustering in all the monosaccharide solutions, while cluster formation does not depend on temperature. In the clusters, some water molecules in the hydration shell are replaced with monosaccharide oxygen atoms, which contributes to the shrinkage of the hydration shell with increasing monosaccharide concentration. The monosaccharides basically adopt one of two conformers, the stable chair or the unstable boat conformer. We revealed that the hydration structures of the boat and chair conformers were dramatically different. As the temperature increases, the content of the chair conformer decreases. Thus, the conversion of conformers strongly affects the hydration structure around the monosaccharide. These results are critical to understand the important roles of the hydration structure in glycan solutions.

Hydrogen Bond Network of Water around Protein Investigated with Terahertz and Infrared Spectroscopy

The dynamical and structural properties of water at protein interfaces were characterized on the basis of the broadband complex dielectric constant (0.25 to 400 THz) of albumin aqueous solutions. Our analysis of the dielectric responses between 0.25 and 12 THz first revealed hydration water with retarded reorientational dynamics extending ∼8.5 Å (corresponding to three to four layers) out from the albumin surface. Second, the number of nonhydrogen-bonded water was decreased in the presence of the albumin solute, indicating protein inhibits the fragmentation of the water hydrogen-bond network. Finally, water molecules at the albumin interface were found to form a distorted hydrogen-bond structure due to topological and energetic disorder of the protein surface. In addition, the intramolecular O-H stretching vibration of water (∼100 THz), which is sensitive to hydrogen-bond environment, pointed to a trend that hydration water has a larger population of strongly hydrogen-bonded water molecules compared with that of bulk water. From these experimental results, we concluded that the "strengthened" water hydrogen bonds at the protein interface dynamically slow down the reorientational motion of water and form the less-defective hydrogen-bond network by inhibiting the fragmentation of water-water hydrogen bonds. Nevertheless, such a strengthened water hydrogen-bond network is composed of heterogeneous hydrogen-bond distances and angles, and thus characterized as structurally "distorted."

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

Are Waters around RNA More than Just a Solvent? - An Insight from Molecular Dynamics Simulations

Hydrating water molecules are believed to be an inherent part of the RNA structure and have a considerable impact on RNA conformation. However, the magnitude and mechanism of the interplay between water molecules and the RNA structure are still poorly understood. In principle, such hydration effects can be studied by molecular dynamics (MD) simulations. In our recent MD studies, we observed that the choice of water model has a visible impact on the predicted structure and structural dynamics of RNA and, in particular, has a larger effect than type, parametrization, and concentration of the ions. Furthermore, the water model effect is sequence dependent and modulates the sequence dependence of A-RNA helical parameters. Clearly, the sensitivity of A-RNA structural dynamics to the water model parametrization is a rather spurious effect that complicates MD studies of RNA molecules. These results nevertheless suggest that the sequence dependence of the A-RNA structure, usually attributed to base stacking, might be driven by the structural dynamics of specific hydration. Here, we present a systematic MD study that aimed to (i) clarify the atomistic mechanism of the water model sensitivity and (ii) discover whether and to what extent specific hydration modulates the A-RNA structural variability. We carried out an extended set of MD simulations of canonical A-RNA duplexes with TIP3P, TIP4P/2005, TIP5P, and SPC/E explicit water models and found that different water models provided a different extent of water bridging between 2'-OH groups across the minor groove, which in turn influences their distance and consequently also inclination, roll, and slide parameters. Minor groove hydration is also responsible for the sequence dependence of these helical parameters. Our simulations suggest that TIP5P is not optimal for RNA simulations.

Revisiting the conundrum of trehalose stabilization

Protein aggregation and loss of protein's biological functionality are manifestations of protein instability. Cosolvents, in particular trehalose, are widely accepted antidotes against such destabilization. Although numerous theories have been promulgated in the literature with regard to its mechanism of stabilization, the present scenario is still elusive in view of the discrepancies existing in them. To this end, we have revisited the conundrum and attempted to rationalize the mechanism by conducting thorough investigation of the effect of trehalose on the native, partially unfolded and denatured states of protein "Lysozyme" by means of molecular dynamic (MD) simulations under different temperature and concentration regimes. Two-dimensional contour plots along with principal component analysis suggest that trehalose molecules offer on-pathway stabilization unaltering the principal direction of protein's motion, although it slows down protein dynamics so that the protein gets trapped in the homogeneous ensemble of conformations closer to the native state. Free energy landscape reveals higher population of native compared to intermediate and denatured states. Delphi results and calculation of the preferential interaction parameter demonstrate that this relative stabilization of the native state can be ascribed to be the consequence of favourable interactions of trehalose with side chains of certain loci on the protein surface encompassing polar flexible residues. Stability of protein results from the observed difference in binding affinity of trehalose for native and denatured states of protein. Our findings are at variance with the common conception of relative destabilization of the denatured state. Rather, we provide evidence for relative stabilization of the native state. This stabilization is due to interplay of protein-trehalose, water-trehalose, water-water, protein-water and trehalose-trehalose interactions.

Water response to intense electric fields: A molecular dynamics study

This paper investigated polarization properties of water molecules in close proximity to an ionic charge in the presence of external electric fields by using an approach based on simulations at the atomic level. We chose sodium and chloride ions in water as examples of dilute ionic solutions and used molecular dynamics simulations to systematically investigate the influence of an external static electric field on structural, dipolar, and polarization properties of water near charged ions. Results showed that a threshold electric field higher than 10(8) V/m is needed to affect water polarization and increase mean dipole moment of water molecules close to the ion. A similar threshold holds for water permittivity profiles, although a field 10× higher is needed to ensure that water permittivity is almost constant independently of the position close to the ion. Electric fields of such intensities can greatly enhance polarizability of water in hydration shells around ions.

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.

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-protein interactions: the secret of protein dynamics

Water-protein interactions help to maintain flexible conformation conditions which are required for multifunctional protein recognition processes. The intimate relationship between the protein surface and hydration water can be analyzed by studying experimental water properties measured in protein systems in solution. In particular, proteins in solution modify the structure and the dynamics of the bulk water at the solute-solvent interface. The ordering effects of proteins on hydration water are extended for several angstroms. In this paper we propose a method for analyzing the dynamical properties of the water molecules present in the hydration shells of proteins. The approach is based on the analysis of the effects of protein-solvent interactions on water protons NMR relaxation parameters. NMR relaxation parameters, especially the nonselective (R₁^NS ) and selective (R₁^SE ) spin-lattice relaxation rates of water protons, are useful for investigating the solvent dynamics at the macromolecule-solvent interfaces as well as the perturbation effects caused by the water-macromolecule interactions on the solvent dynamical properties. In this paper we demonstrate that Nuclear Magnetic Resonance Spectroscopy can be used to determine the dynamical contributions of proteins to the water molecules belonging to their hydration shells.

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.