Hydration water and bulk water in proteins have distinct properties in radial distributions calculated from 105 atomic resolution crystal structures

Water plays a critical role in the structure and function of proteins, although the experimental properties of water around protein structures are not well understood. The water can be classified by the separation from the protein surface into bulk water and hydration water. Hydration water interacts closely with the protein and contributes to protein folding, stability, and dynamics, as well as interacting with the bulk water. Water potential functions are often parametrized to fit bulk water properties because of the limited experimental data for hydration water. Therefore, the structural and energetic properties of the hydration water were assessed for 105 atomic resolution (相似文献   

Hydrogen bond breaking mechanism and water reorientational dynamics in the hydration layer of lysozyme

The mechanism and the rate of hydrogen bond-breaking in the hydration layer surrounding an aqueous protein are important ingredients required to understand the various aspects of protein dynamics, its function, and stability. Here, we use computer simulation and a time correlation function technique to understand these aspects in the hydration layer of lysozyme. Water molecules in the layer are found to exhibit three distinct bond-breaking mechanisms. A large angle orientational jump of the donor water molecule is common among all of them. In the most common ( approximately 80%) bond-breaking event in the layer, the new acceptor water molecule comes from the first coordination shell (initially within 3.5 A of the donor), and the old acceptor water molecule remains within the first coordination shell, even after the bond-breaking. This is in contrast to that in bulk water, in which both of the acceptor molecules involve the second coordination shell. Additionally, the motion of the incoming and the outgoing acceptor molecules involved is not diffusive in the hydration layer, in contrast to their observed diffusive motion in the bulk. The difference in rotational dynamics between the bulk and the hydration layer water molecules is clearly manifested in the calculated time-dependent angular van Hove self-correlation function ( G(theta, t)) which has a pronounced two-peak structure in the layer, and this can be traced to the constrained translational motion in the layer. The longevity of the surrounding hydrogen bond network is found to be significantly enhanced near a hydrophilic residue.     

Water structure about the dimer and hexamer repeat units of amylose from molecular dynamics computer simulations

We have analyzed a set of molecular dynamics (MD) trajectories of maltose in vacuum and water for solute imposed structuring on the solvent. To do this, we used a novel technique to calculate water probability densities to locate the areas in which the solvent is most populated in the maltose solution. We found that only the layer of water within the first maltose hydration shell has a probability density 50% and greater than that of bulk water. On investigating this water layer using Voronoi polyhedra (VP) analysis it was seen that only the waters adjacent to the hydrophobic (CH and CH₂) groups are more structured than bulk water. We found that in a maltose solution of approximately 1.0 g/cm³ the solute does not disrupt the structure of the surrounding water beyond the first hydration shell. Next we performed a 700‐ps MD simulation of a maltohexaose strand in a box of 4096 SPC/E waters. The water probability density calculations and the VP analysis of the maltohexaose solution show that the larger amylose repeat unit decreases the solvent configurational entropy of the water beyond the first hydration shell. Analysis of this trajectory reveals that the helical conformation of the maltohexaose strand is preserved via bridging intermolecular water hydrogen bonds, indicating that a single amylose helical turn in water is preserved by hydrophilic and not hydrophobic interactions. Using VP analysis we present a method to accurately determine the number of water molecules in the first hydration shell of dissolved solutes. In the case of maltose, there are 40 water molecules in this shell, while for maltohexaose the number is 98. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 445–456, 2001     

On the structure and dynamics of the hydrated sulfite ion in aqueous solution--an ab initio QMCF MD simulation and large angle X-ray scattering study

Theoretical ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism has been applied in conjunction to experimental large angle X-ray scattering to study the structure and dynamics of the hydrated sulfite ion in aqueous solution. The results show that there is a considerable effect of the lone electron-pair on sulfur concerning structure and dynamics in comparison with the sulfate ion with higher oxidation number and symmetry of the hydration shell. The S-O bond distance in the hydrated sulfite ion has been determined to 1.53(1) ? by both methods. The hydrogen bonds between the three water molecules bound to each sulfite oxygen are only slightly stronger than those in bulk water. The sulfite ion can therefore be regarded as a weak structure maker. The water exchange rate is somewhat slower for the sulfite ion than for the sulfate ion, τ(0.5) = 3.2 and 2.6 ps, respectively. An even more striking observation in the angular radial distribution (ARD) functions is that the for sulfite ion the water exchange takes place in close vicinity of the lone electron-pair directed at its sides, while in principle no water exchange did take place of the water molecules hydrogen bound to sulfite oxygens during the simulation time. This is also confirmed when detailed pathway analysis is conducted. The simulation showed that the water molecules hydrogen bound to the sulfite oxygens can move inside the hydration shell to the area outside the lone electron-pair and there be exchanged. On the other hand, for the hydrated sulfate ion in aqueous solution one can clearly see from the ARD that the distribution of exchange events is symmetrical around the entire hydration sphere.     

Long-range influence of carbohydrates on the solvation dynamics of water--answers from terahertz absorption measurements and molecular modeling simulations

We present new terahertz (THz) spectroscopic measurements of solvated sugars and compare the effect of two disaccharides (trehalose and lactose) and one monosaccharide (glucose) with respect to the solute-induced changes in the sub-picosecond network dynamics of the hydration water. We found that the solute affects the fast collective network motions of the solvent, even beyond the first solvation layer. For all three carbohydrates, we find an increase of 2-4% in the THz absorption coefficient of the hydration water in comparison to bulk water. Concentration-dependent changes in the THz absorption between 2.1 and 2.8 THz of the solute-water mixture were measured with a precision better than 1% and were used to deduce a dynamical hydration shell, which extends from the surface up to 5.7 +/- 0.4 and 6.5 +/- 0.9 A for the disaccharides lactose and trehalose, respectively, and 3.7 +/- 0.9 A for the glucose. This exceeds the values for the static hydration shell as determined, for example, by scattering, where the long-range structure was found to be not significantly affected by the solute beyond the first hydration shell. When comparing all three carbohydrates, we found that the solute-induced change in the THz absorption depends on the product of molar concentration of the solute and the number of hydrogen bonds between the carbohydrate and water molecules. We can conclude that the long-range influence on the sub-picosecond collective water network motions of the hydration water is directly correlated with the average number of hydrogen bonds between the molecule and adjacent water molecules for carbohydrates. This implies that monosaccharides have a smaller influence on the surrounding water molecules than disaccharides. This could explain the bioprotection mechanism of sugar-water mixtures, which has been found to be more effective for disaccharides than for monosaccharides.     

Structure and dynamics of the hydration shells of the Al3+ ion

First principles simulations of the hydration shells surrounding Al3+ ions are reported for temperatures near 300 degrees C. The predicted six water molecules in the octahedral first hydration shell were found to be trigonally coordinated via hydrogen bonds to 12 s shell water molecules in agreement with the putative structure used to analyze the x-ray data, but in disagreement with the results reported from conventional molecular dynamics using two-and three-body potentials. Bond lengths and angles of the water molecules in the first and second hydration shells and the average radii of these shells also agreed very well with the results of the x-ray analysis. Water transfers into and out of the second solvation shell were observed to occur on a picosecond time scale via a dissociative mechanism. Beyond the second shell the bonding pattern substantially returned to the tetrahedral structure of bulk water. Most of the simulations were done with 64 solvating water molecules (20 ps). Limited simulations with 128 water molecules (7 ps) were also carried out. Results agreed as to the general structure of the solvation region and were essentially the same for the first and second shell. However, there were differences in hydrogen bonding and Al-O radial distribution function in the region just beyond the second shell. At the end of the second shell a nearly zero minimum in the Al-O radial distribution was found for the 128 water system. This minimum is less pronounced minimum found for the 64 water system, which may indicate that sizes larger than 64 may be required to reliably predict behavior in this region.     

The structure of water in the hydration shell of cations from x-ray Raman and small angle x-ray scattering measurements

X-ray Raman scattering (XRS) spectroscopy and small angle x-ray scattering (SAXS) are used to study water in aqueous solutions of NaCl, MgCl(2), and AlCl(3) with the particular aim to provide information about the structure of the hydration shells of the cations. The XRS spectra show that Na(+) weakens the hydrogen bonds of water molecules in its vicinity, similar to the effect of increased temperature and pressure. Mg(2+) and Al(3+), on the other hand, cause the formation of short and strong hydrogen bonds between the surrounding water molecules. The SAXS data show that Mg(2+) and Al(3+) form tightly bound hydration shells that give a large density contrast in the scattering data. From the form factors extracted from the SAXS data, we found that Mg(2+) and Al(3+) have, respectively, an equivalent of one and one and a half stable hydration shells that appear as a density contrast. In addition, we estimated that the density of water in the hydration shells of Mg(2+) and Al(3+) is, respectively, ～61% and ～71% higher than in bulk water.     

Hydrogen bond dynamics and water structure in glucose-water solutions by depolarized Rayleigh scattering and low-frequency Raman spectroscopy

The effect of glucose on the relaxation process of water at picosecond time scales has been investigated by depolarized Rayleigh scattering (DRS) experiments. The process is assigned to the fast hydrogen bonding dynamics of the water network. In DRS spectra this contribution can be safely separated from the slower relaxation process due to the sugar. The detected relaxation time is studied at different glucose concentrations and modeled considering bulk and hydrating water contributions. As a result, it is found that in diluted conditions the hydrogen bond lifetime of proximal water molecules becomes about three times slower than that of the bulk. The effect of the sugar on the hydrogen bond water structure is investigated by analyzing the low-frequency Raman (LFR) spectrum sensitive to intermolecular modes. The addition of glucose strongly reduces the intensity of the band at 170 cm(-1) assigned to a collective stretching mode of water molecules arranged in cooperative tetrahedral domains. These findings indicate that proximal water molecules partially lose the tetrahedral ordering typical of the bulk leading to the formation of high density environments around the sugar. Thus the glucose imposes a new local order among water molecules localized in its hydration shell in which the hydrogen bond breaking dynamics is sensitively retarded. This work provides new experimental evidences that support recent molecular dynamics simulation and thermodynamics results.     

Orientational dynamics for an amphiphilic-solvent solution

In this work, we performed Monte Carlo simulations on a lattice model for spontaneous amphiphilic aggregation, in order to study the orientational and hydrogen-bonding dynamics of water on different regions inside the micellar solution. We employed an associating lattice gas model that mimics the aqueous solvent, which presents a rich phase diagram with first- and second-order transition lines. Even though this is a simplified model, it makes possible to investigate the orientational dynamics of water in an equilibrium solution of amphiphiles, as well as the influence of the different phases of the solvent in the interfacial and bulk water dynamics. By means of extensive simulations, we showed that, at high temperatures, the behavior of the orientational relaxation and hydrogen bonding of water molecules in the bulk, first, and second hydration shells are considerable different. We observe the appearance of a very slow component for water molecules in the first hydration shell of micelles when the system reaches a high-density phase, consistent with previous theoretical and experimental studies concerning biological water. Also, at high temperatures, we find that water molecules in the second hydration shell of micelles have an orientational decay similar to that of bulk water, but with a generally slower dynamics. Otherwise, at low temperatures, we have two components for the orientational relaxation of bulk water in the low density liquid phase, and only a single component in the high density liquid (HDL) phase, which reflect the symmetry properties of the different phases of the solvent model. In the very dense region of water molecules in the first hydration shell of micelles at low temperatures, we find two components for the orientational relaxation on both liquid phases, one of them much slower than that in the single component of bulk water in the HDL phase. This happens even though our model does not present any hindrance to the water rotational freedom caused by the presence of the amphiphiles.     

Potential of mean force of hydrophobic association: dependence on solute size

The potentials of mean force (PMFs) were determined for systems involving formation of nonpolar dimers composed of methane, ethane, propane, isobutane, and neopentane, respectively, in water, using the TIP3P water model, and in vacuo. A series of umbrella-sampling molecular dynamics simulations with the AMBER force field was carried out for each pair in either water or in vacuo. The PMFs were calculated by using the weighted histogram analysis method (WHAM). The shape of the PMFs for dimers of all five nonpolar molecules is characteristic of hydrophobic interactions with contact and solvent-separated minima and desolvation maxima. The positions of all these minima and maxima change with the size of the nonpolar molecule, that is, for larger molecules they shift toward larger distances. The PMF of the neopentane dimer is similar to those of other small nonpolar molecules studied in this work, and hence the neopentane dimer is too small to be treated as a nanoscale hydrophobic object. The solvent contribution to the PMF was also computed by subtracting the PMF determined in vacuo from the PMF in explicit solvent. The molecular surface area model correctly describes the solvent contribution to the PMF together with the changes of the height and positions of the desolvation barrier for all dimers investigated. The water molecules in the first solvation sphere of the dimer are more ordered compared to bulk water, with their dipole moments pointing away from the surface of the dimer. The average number of hydrogen bonds per water molecule in this first hydration shell is smaller compared to that in bulk water, which can be explained by coordination of water molecules to the hydrocarbon surface. In the second hydration shell, the average number of hydrogen bonds is greater compared to bulk water, which can be explained by increased ordering of water from the first hydration shell; the net effect is more efficient hydrogen bonding between the water molecules in the first and second hydration shells.     

Equatorial and apical solvent shells of the UO2 2+ ion

First principles molecular dynamics simulations of the hydration shells surrounding UO(2)(2+) ions are reported for temperatures near 300 K. Most of the simulations were done with 64 solvating water molecules (22 ps). Simulations with 122 water molecules (9 ps) were also carried out. The hydration structure predicted from the simulations was found to agree with very well-known results from x-ray data. The average U=O bond length was found to be 1.77 A. The first hydration shell contained five trigonally coordinated water molecules that were equatorially oriented about the O-U-O axis with the hydrogen atoms oriented away from the uranium atom. The five waters in the first shell were located at an average distance of 2.44 A (2.46 A, 122 water simulation). The second hydration shell was composed of distinct equatorial and apical regions resulting in a peak in the U-O radial distribution function at 4.59 A. The equatorial second shell contained ten water molecules hydrogen bonded to the five first shell molecules. Above and below the UO(2)(2+) ion, the water molecules were found to be significantly less structured. In these apical regions, water molecules were found to sporadically hydrogen bond to the oxygen atoms of the UO(2)(2+), oriented in such a way as to have their protons pointed toward the cation. While the number of apical waters varied greatly, an average of five to six waters was found in this region. Many water transfers into and out of the equatorial and apical second solvation shells were observed to occur on a picosecond time scale via dissociative mechanisms. Beyond these shells, the bonding pattern substantially returned to the tetrahedral structure of bulk water.     

Thermal signature of hydrophobic hydration dynamics

Hydrophobic hydration, the perturbation of the aqueous solvent near an apolar solute or interface, is a fundamental ingredient in many chemical and biological processes. Both bulk water and aqueous solutions of apolar solutes behave anomalously at low temperatures for reasons that are not fully understood. Here, we use (2)H NMR relaxation to characterize the rotational dynamics in hydrophobic hydration shells over a wide temperature range, extending down to 243 K. We examine four partly hydrophobic solutes: the peptides N-acetyl-glycine-N'-methylamide and N-acetyl-leucine-N'-methylamide, and the osmolytes trimethylamine N-oxide and tetramethylurea. For all four solutes, we find that water rotates with lower activation energy in the hydration shell than in bulk water below 255 +/- 2 K. At still lower temperatures, water rotation is predicted to be faster in the shell than in bulk. We rationalize this behavior in terms of the geometric constraints imposed by the solute. These findings reverse the classical "iceberg" view of hydrophobic hydration by indicating that hydrophobic hydration water is less ice-like than bulk water. Our results also challenge the "structural temperature" concept. The two investigated osmolytes have opposite effects on protein stability but have virtually the same effect on water dynamics, suggesting that they do not act indirectly via solvent perturbations. The NMR-derived picture of hydrophobic hydration dynamics differs substantially from views emerging from recent quasielastic neutron scattering and pump-probe infrared spectroscopy studies of the same solutes. We discuss the possible reasons for these discrepancies.     

Anomalously Enhanced Hydration of Aqueous Electrolyte Solution in Hydrophobic Carbon Nanotubes to Maintain Stability

An understanding of the structure and behavior of electrolyte solutions in nanoenvironements is crucial not only for a wide variety of applications, but also for the development of physical, chemical, and biological processes. We demonstrate the structure and stability of electrolyte in carbon nanotubes using hybrid reverse Monte Carlo simulations of X‐ray diffraction patterns. Hydrogen bonds between water are adequately formed in carbon nanotubes, although some hydrogen bonds are restricted by the interfaces of carbon nanotubes. The hydrogen bonding network of water in electrolyte in the carbon nanotubes is further weakened. On the other hand, formation of the ion hydration shell is significantly enhanced in the electrolyte in the carbon nanotubes in comparison to ion hydration in bulk electrolyte. The significant hydrogen bond and hydration shell formation are a result of gaining stability in the hydrophobic nanoenvironment.     

Hydration structure of water confined between mica surfaces

We report further molecular dynamics simulations on the structure of bound hydration layers under extreme confinement between mica surfaces. We find that the liquid phase of water is maintained down to 2 monolayer (ML) thick, whereas the structure of the K(+) ion hydration shell is close to the bulk structure even under D = 0.92 nm confinement. Unexpectedly, the density of confined water remains approximately the bulk value or less, whereas the diffusion of water molecules decreases dramatically. Further increase in confinement leads to a transition to a bilayer ice, whose density is much less than that of ice Ih due to the formation of a specific hydrogen-bonding network.     

Low-frequency vibrational spectrum of water in the hydration layer of a protein: a molecular dynamics simulation study

An atomistic molecular dynamics simulation has been carried out to understand the low-frequency intermolecular vibrational spectrum of water present in the hydration layer of the protein villin headpiece subdomain or HP-36. An attempt is made to explore how the heterogeneous rigidity of the hydration layers of different segments (three alpha helices) of the protein, strength of the protein-water hydrogen bonds, and their differential relaxation behavior influence the distribution of the intermolecular vibrational density of states of water in the hydration layers. The calculations revealed that compared to bulk water these bands are nonuniformly blue-shifted for water near the helices, the extent of shifts being more pronounced for water molecules hydrogen bonded to the protein residues. It is further noticed that the larger blue shift observed for the water molecules hydrogen bonded to helix 2 residues correlates excellently with the slowest structural relaxation of these hydrogen bonds. These results can be verified by suitable experimental measurements.     

Tritium/Protium Fractionation Near and Inside DNA

The fractionation factor of tritium between water and DNA as well as between water and the first hydration shell of DNA is determined. For this purpose the sublimation kinetics of water from DNA dissolved in water were determined at -200 °C and tritium was measured in the remaining water free DNA. The last sublimating water fractions showed a tritium level of about 1.4, the residual water free DNA about 1.9 units above the bulk water. The tritium accumulation inside and near DNA is attributed to the thermodynamic triton-proton exchange isotope effect existing between the strong hydrogen bridges of water and weaker hydrogen bridges as well as inside DNA as between the first hydration shell and DNA.     

水环境对GC和AT碱基对质子转移的影响

采用B3LYP/DZP++的方法研究了第一水化层作用和连续化处理的水溶剂作用对鸟嘌呤-胞嘧啶(GC)碱基对和腺嘌呤-胸腺嘧啶(AT)碱基对质子转移反应的影响. GC和AT碱基对在连续化水溶剂作用下,均发生单质子转移(SPT1)和分步的双质子转移(DPT),而在第一水化层5 个水分子的作用下(GC·5H₂O,AT·5H₂O)或同时考虑第一水化层作用和连续化水溶剂作用(GC·5H₂O+PCM,AT·5H₂O+PCM)时,GC和AT碱基对的质子转移均只得到单质子转移反应(SPT1). 单质子转移过程中的活化能变化情况表明：第一水化层对GC和AT碱基对结构和质子转移影响较大,水环境对碱基对的作用主要发生在第一水化层.     

How protein surfaces induce anomalous dynamics of hydration water

Water around biomolecules slows down with respect to pure water, and both rotation and translation exhibit anomalous time dependence in the hydration shell. The origin of such behavior remains elusive. We use molecular dynamics simulations of water dynamics around several designed protein models to establish the connection between the appearance of the anomalous dynamics and water-protein interactions. For the first time we quantify the separate effect of protein topological and energetic disorder on the hydration water dynamics. When a static protein structure is simulated, we show that both types of disorder contribute to slow down water diffusion, and that allowing for protein motion, increasing the spatial dimensionality of the interface, reduces the anomalous character of hydration water. The rotation of water is, instead, altered by the energetic disorder only; indeed, when electrostatic interactions between the protein and water are switched off, water reorients even faster than in the bulk. The dynamics of water is also related to the collective structure--à voir the hydrogen bond (H-bond) network--formed by the solvent enclosing the protein surface. We show that, as expected for a full hydrated protein, when the protein surface offers pinning sites (charged or polar sites), the superficial water-water H-bond network percolates throughout the whole surface, hindering the water diffusion, whereas it does not when the protein surface lacks electrostatic interactions with water and the water diffusion is enhanced.     

Tl(I)‐the strongest structure‐breaking metal ion in water? A quantum mechanical/molecular mechanical simulation study

Structural and dynamical properties of the Tl(I) ion in dilute aqueous solution have been investigated by ab initio quantum mechanics in combination with molecular mechanics. The first shell plus a part of the second shell were treated by quantum mechanics at Hartree-Fock level, the rest of the system was described by an ab initio constructed potential. The radial distribution functions indicate two different bond lengths (2.79 and 3.16 A) in the first hydration shell, in good agreement with large-angle X-ray scattering and extended X-ray absorption fine structure spectroscopy results. The average first shell coordination number was found as 5.9, and several other structural parameters such as coordination number distributions, angular distribution functions, and tilt- and theta-angle distributions were evaluated. The ion-ligand vibration spectrum and reorientational times were obtained via velocity auto correlation functions. The Tl-O stretching force constant is very weak with 5.0 N m(-1). During the simulation, numerous water exchange processes took place between first and second hydration shell and between second shell and bulk. The mean ligand residence times for the first and second shell were determined as 1.3 and 1.5 ps, respectively, indicating Tl(I) to be a typical "structure-breaker". The calculated hydration energy of -84 +/- 16 kcal mol(-1) agrees well with the experimental value of -81 kcal mol(-1). All data obtained for structure and dynamics of hydrated Tl(I) characterize this ion as a very special case among all monovalent metal ions, being the most potent "structure-breaker", but at the same time forming a distinct second hydration shell and thus having a far-reaching influence on the solvent structure.