A unified picture of protein hydration: prediction of hydrodynamic properties from known structures.

Hydration is essential for the structural and functional integrity of globular proteins. How much hydration water is required for that integrity? A number of techniques such as X-ray diffraction, nuclear magnetic resonance (NMR) spectroscopy, calorimetry, infrared spectroscopy, and molecular dynamics (MD) simulations indicate that the hydration level is 0.3-0.5 g of water per gram of protein for medium sized proteins. Hydrodynamic properties, when accounted for by modeling proteins as ellipsoids, appear to give a wide range of hydration levels. In this paper we describe an alternative numerical technique for hydrodynamic calculations that takes account of the detailed protein structures. This is made possible by relating hydrodynamic properties (translational and rotational diffusion constants and intrinsic viscosity) to electrostatic properties (capacitance and polarizability). We show that the use of detailed protein structures in predicting hydrodynamic properties leads to hydration levels in agreement with other techniques. A unified picture of protein hydration emerges. There are preferred hydration sites around a protein surface. These sites are occupied nearly all the time, but by different water molecules at different times. Thus, though a given water molecule may have a very short residence time (approximately 100-500 ps from NMR spectroscopy and MD simulations) in a particular site, the site appears fully occupied in experiments in which time-averaged properties are measured.     

Dynamics of hydration in hen egg white lysozyme

We investigate the hydration dynamics of a small globular protein, hen egg-white lysozyme. Extensive simulations (two trajectories of 9 ns each) were carried out to identify the time-scales and mechanism of water attachment to this protein. The location of the surface and integral water molecules in lysozyme was also investigated. Three peculiar temporal scales of the hydration dynamics can be discerned: two among these, with sub-nanosecond mean residence time, tau(w), are characteristic of surface hydration water; the slower time-scale (tau(w) approximately 2/3 ns) is associated with buried water molecules in hydrophilic pores and in superficial clefts. The computed tau(w) values in the two independent runs fall in a similar range and are consistent with each other, thus adding extra weight to our result. The tau(w) of surface water obtained from the two independent trajectories is 20 and 24 ps. In both simulations only three water molecules are bound to lysozyme for the entire length of the trajectories, in agreement with nuclear magnetic relaxation dispersion estimates. Locations other than those identified in the protein crystal are found to be possible for these long-residing water molecules. The dynamics of the hydration water molecules observed in our simulations implies that each water molecule visits a multitude of residues during the lifetime of its bound with the protein. The number of residues seen by a single water molecule increases with the time-scale of its residence time and, on average, is equal to one only for the water molecules with shorter residence time. Thus, tau(w) values obtained from inelastic neutron scattering and based on jump-diffusion models are likely not to account for the contribution of water molecules with longer residence time.     

Thermodynamic analysis of water molecules at the surface of proteins and applications to binding site prediction and characterization

Water plays an essential role in determining the structure and function of all biological systems. Recent methodological advances allow for an accurate and efficient estimation of the thermodynamic properties of water molecules at the surface of proteins. In this work, we characterize these thermodynamic properties and relate them to various structural and functional characteristics of the protein. We find that high-energy hydration sites often exist near protein motifs typically characterized as hydrophilic, such as backbone amide groups. We also find that waters around alpha helices and beta sheets tend to be less stable than waters around loops. Furthermore, we find no significant correlation between the hydration site-free energy and the solvent accessible surface area of the site. In addition, we find that the distribution of high-energy hydration sites on the protein surface can be used to identify the location of binding sites and that binding sites of druggable targets tend to have a greater density of thermodynamically unstable hydration sites. Using this information, we characterize the FKBP12 protein and show good agreement between fragment screening hit rates from NMR spectroscopy and hydration site energetics. Finally, we show that water molecules observed in crystal structures are less stable on average than bulk water as a consequence of the high degree of spatial localization, thereby resulting in a significant loss in entropy. These findings should help to better understand the characteristics of waters at the surface of proteins and are expected to lead to insights that can guide structure-based drug design efforts.     

Hydration dynamics of the collagen triple helix by NMR

The hydration of the collagen-like Ac-(Gly-Pro-Hyp)(6)-NH(2) triple-helical peptide in solution was investigated using an integrated set of high-resolution NMR hydration experiments, including different recently developed exchange-network editing methods. This approach was designed to explore the hydration dynamics in the proximity of labile groups, such as the hydroxyproline hydroxyl group, and revealed that the first shell of hydration in collagen-like triple helices is kinetically labile with upper limits for water molecule residence times in the nanosecond to sub-nanosecond range. This result is consistent with a "hopping" hydration model in which solvent molecules are exchanged in and out of solvation sites at a rate that is not directly correlated to the degree of site localization. The hopping model thus reconciles the dynamic view of hydration revealed by NMR with the previously suggested partially ordered semi-clathrate-like cylinder of hydration. In addition, the nanosecond to sub-nanosecond upper limits for water molecule residence times imply that hydration-dehydration events are not likely to be the rate-limiting step for triple helix self-recognition, complementing previous investigations on water dynamics in collagen fibers. This study has also revealed labile proton features expected to facilitate the characterization of the structure and folding of triple helices in collagen peptides.     

Residence times of water molecules in the hydration sites of myoglobin

Hydration sites are high-density regions in the three-dimensional time-averaged solvent structure in molecular dynamics simulations and diffraction experiments. In a simulation of sperm whale myoglobin, we found 294 such high-density regions. Their positions appear to agree reasonably well with the distributions of waters of hydration found in 38 x-ray and 1 neutron high-resolution structures of this protein. The hydration sites are characterized by an average occupancy and a combination of residence time parameters designed to approximate a distribution of residence times. It appears that although the occupancy and residence times of the majority of sites are rather bulk-like, the residence time distribution is shifted toward the longer components, relative to bulk. The sites with particularly long residence times are located only in the cavities and clefts of the protein. This indicates that other factors, such as hydrogen bonds and hydrophobicity of underlying protein residues, play a lesser role in determining the residence times of the longest-lived sites.     

NMR observation of individual molecules of hydration water bound to DNA duplexes: direct evidence for a spine of hydration water present in aqueous solution.

The residence times of individual hydration water molecules in the major and minor grooves of DNA were measured by nuclear magnetic resonance (NMR) spectroscopy in aqueous solutions of d-(CGCGAATTCGCG)2 and d-(AAAAATTTTT)2. The experimental observations were nuclear Overhauser effects (NOE) between water protons and the protons of the DNA. The positive sign of NOEs with the thymine methyl groups shows that the residence times of the hydration water molecules near these protons in the major groove of the DNA must be shorter than about 500 ps, which coincides with the behavior of surface hydration water in peptides and proteins. Negative NOEs were observed with the hydrogen atoms in position 2 of adenine in both duplexes studied. This indicates that a 'spine of hydration' in the minor groove, as observed by X-ray diffraction in DNA crystals, is present also in solution, with residence times significantly longer than 1 ns. Such residence times are reminiscent of 'interior' hydration water molecules in globular proteins, which are an integral part of the molecular architecture both in solution and in crystals.     

Protein-water interactions in ribonuclease A and angiogenin: a molecular dynamics study

It is known that water molecules play an important role in the biological functioning of proteins. The members of the ribonuclease A (RNase A) family of proteins, which are sequentially and structurally similar, are known to carry out the obligatory function of cleaving RNA and individually perform other diverse biological functions. Our focus is on elucidating whether the sequence and structural similarity lead to common hydration patterns, what the common hydration sites are and what the differences are. Extensive molecular dynamics simulations followed by a detailed analysis of protein-water interactions have been carried out on two members of the ribonuclease A superfamily-RNase A and angiogenin. The water residence times are analyzed and their relationship with the characteristic properties of the protein polar atoms, such as their accessible surface area and mean hydration, is studied. The capacity of the polar atoms to form hydrogen bonds with water molecules and participate in protein-water networks are investigated. The locations of such networks are identified for both proteins.     

Molecular dynamics simulation of solvated azurin: correlation between surface solvent accessibility and water residence times

A system containing the globular protein azurin and 3,658 water molecules has been simulated to investigate the influence on water dynamics exerted by a protein surface. Evaluation of water mean residence time for elements having different secondary structure did not show any correlation. Identically, comparison of solvent residence time for atoms having different charge and polarity did not show any clear trend. The main factor influencing water residence time in proximity to a specific site was found to be its solvent accessibility. In detail for atoms belonging to lateral chains and having solvent-accessible surface lower than approximately 16 A(2)a relation is found for which charged and polar atoms are surrounded by water molecules characterized by residence times longer than the non polar ones. The involvement of the low accessible protein atom in an intraprotein hydrogen bond further modulates the length of the water residence time. On the other hand for surfaces having high solvent accessibility, all atoms, independently of their character, are surrounded by water molecules which rapidly exchange with the bulk solvent. Proteins 2000;39:56-67.     

Cluster analysis of hydration waters around the active sites of bacterial alanine racemase using a 2-ns MD simulation

Structural data produced by a 2-ns molecular dynamics (MD) simulation on Geobacillus alanine racemase (AlaR; PDB: 1SFT) was used to study hydration around the two AlaR active sites. AlaR is a crucial enzyme for bacterial cell wall biosynthesis. It has been shown previously that the potency of an inhibitor can be increased by incorporating a functional group or atom that displaces hydration sites close to the substrate binding pocket of its target enzyme. The complete linkage algorithm was used for cluster analysis of the active site water positions from 126 solvent configurations sampled at regular intervals from the 2-ns MD simulation. Crystal waters in the 1SFT X-ray structure occupy most of the tightly bound water sites that were discovered. We show here that tightly bound water sites can be identified by cluster analysis of MD-generated coordinates starting with data supplied by a single X-ray structure, and we predict a highly conserved hydration site close to the carboxyl oxygen of L-Ala substrate. This approach holds promise for accelerating the drug design process. We also discuss an analysis of the well-known notion of residence time and introduce a new measure called retention time.     

Dissection of the structural and functional role of a conserved hydration site in RNase T1

The reoccurrence of water molecules in crystal structures of RNase T1 was investigated. Five waters were found to be invariant in RNase T1 as well as in six other related fungal RNases. The structural, dynamical, and functional characteristics of one of these conserved hydration sites (WAT1) were analyzed by protein engineering, X-ray crystallography, and (17)O and 2H nuclear magnetic relaxation dispersion (NMRD). The position of WAT1 and its surrounding hydrogen bond network are unaffected by deletions of two neighboring side chains. In the mutant Thr93Gln, the Gln93N epsilon2 nitrogen replaces WAT1 and participates in a similar hydrogen bond network involving Cys6, Asn9, Asp76, and Thr91. The ability of WAT1 to form four hydrogen bonds may explain why evolution has preserved a water molecule, rather than a side-chain atom, at the center of this intricate hydrogen bond network. Comparison of the (17)O NMRD profiles from wild-type and Thr93Gln RNase T1 yield a mean residence time of 7 ns at 27 degrees C and an orientational order parameter of 0.45. The effects of mutations around WAT1 on the kinetic parameters of RNase T1 are small but significant and probably relate to the dynamics of the active site.     

Continuum electrostatic analysis of preferred solvation sites around proteins in solution

To understand water-protein interactions in solution, the electrostatic field is calculated by solving the Poisson-Boltzmann equation, and the free energy surface of water is mapped by translating and rotating an explicit water molecule around the protein. The calculation is applied to T4 lysozyme with data available on the conservation of solvent binding sites in 18 crystallographically independent molecules. The free energy maps around the ordered water sites provide information on the relationship between water positions in crystal structure and in solution. Results show that almost all conserved sites and the majority of nonconserved sites are within 1.3 A of local free energy minima. This finding is in sharp contrast to the behavior of randomly placed water molecules in the boundary layer, which, on the average, must travel more than 3 A to the nearest free energy minimum. Thus, the solvation sites are at least partially determined by protein-water interactions rather than by crystal packing alone. The characteristic water residence times, obtained from the free energies at the local minima, are in good agreement with nuclear magnetic resonance experiments. Only about half of the potential sites show up as ordered water in the 1.7 A resolution X-ray structure. Crystal packing interactions can stabilize weak or mobile potential sites (in fact, some ordered water positions are not close to free energy minima) or can prevent water from occupying certain sites. Apart from a few buried water molecules that are strong binders, the free energies are not very different for conserved and nonconserved sites. We show that conservation of a water site between two crystals occurs if the positions of protein atoms, primarily contributing to the free energy at the local minimum, do not substantially change from one structure to the other. This requirement can be correlated with the nature of the side chain contacting the water molecule in the site.     

Orientational order and dynamics of hydration water in a single crystal of bovine pancreatic trypsin inhibitor.

The orientational order and dynamics of the water molecules in form II crystals of bovine pancreatic trypsin inhibitor (BPTI) are studied by (2)H NMR in the temperature range 6-50 degrees C. From the orientation dependence of the single crystal quadrupole splitting and linewidth, the principal components of the motionally averaged quadrupole interaction tensor and the irreducible linewidth components for the orthorhombic crystal are determined. With the aid of water orientations derived from neutron and x-ray diffraction, it is shown that the NMR data can be accounted for by a small number of highly ordered crystal waters, some of which have residence times in the microsecond range. Most of these specific hydration sites must be located at intermolecular contacts. The surface hydration layer that is also present in dilute solution is likely to be only weakly ordered and would then not contribute significantly to the splitting and linewidth from the protein crystal. To probe water dynamics on shorter time scales, the (2)H longitudinal relaxation dispersion is measured for a polycrystalline BPTI sample. The observed dispersion is dominated by rapidly exchanging deuterons in protein side chains, undergoing restricted rotational motions on a time scale of 10 ns.     

Hydration of a hydrophobic cavity and its functional role: a simulation study of human interleukin-1beta

Molecular dynamics simulations reveal that the hydrophobic cavity in human cytokine Interleukin-1beta is hydrated and can dynamically accommodate between one and four water molecules. These waters have residence times > 500 ps and can give rise to detectable NOEs, in agreement with NMR observations of Ernst et al. (Science 1995; 267:1813-1817). The waters also display high positional disorder within the cavity, which explains why they have not been resolved crystallographically. The average distribution of water molecules over time within the cavity matches well the low resolution electron density extracted by Yu et al. (Proc Natl Acad Sci 1999; 96:103-108). The water molecules hydrate the hydrophobic cavity preferentially as complex clusters. These clusters result from a combination of hydrogen bonds between the waters and stabilizing interactions between the waters and aromatic rings forming the cavity. Free energy estimates suggest that it takes 4-waters to hydrate the cavity in a thermodynamically stable manner leading to a gain in free energy of transfer from bulk of approximately approximately 3.6 kcal/mol. This arises from the existence of the water clusters in multiple hydrogen bonded states. In addition, the waters are found to migrate either individually or as clusters out of the cavity through several pathways. The upper limit for one-dimensional diffusion of the waters within the protein matrix is 4 A/ps (relative to 6 A/ps for bulk). Simulations reveal pathways in addition to those identified crystallographically, with motions controlled by the rotations of sidechains. We find that only when the hydrophobic cavity is hydrated, do correlated motions couple distant sites with the sites that make contact with the receptor and this data partly offers an explanation of experimental mutagenesis data. Simulations, together with recent observations based on mutagenesis by Heidary et al. (J Mol Biol 2005; 353:1187-1198) that hydrogen bond networks couple motions across long distances in interleukin-1beta, lead us to hypothesize that the hydration of the cavity (conserved across mammals) can thermodynamically enhance hydrogen bond networks to enable coupling across long distances by acting as a plug and this in turn enables a kinetic control of the rate of transmission of signals.     

Structure and hydration of BamHI DNA recognition site: a molecular dynamics investigation

The results of a 3-ns molecular dynamics simulation of the dodecamer duplex d(TATGGATCCATA)(2) recognized by the BamHI endonuclease are presented here. The DNA has been simulated as a flexible molecule using an AMBER force field and the Ewald summation method, which eliminates the undesired effects of truncation and permits evaluation of the full effects of electrostatic forces. The starting B conformation evolves toward a configuration quite close to that observed through x-ray diffraction in its complex with BamHI. This configuration is fairly stable and the Watson-Crick hydrogen bonds are well maintained over the simulation trajectory. Hydration analysis indicates a preferential hydration for the phosphate rather than for the ester oxygens. Hydration shells in both the major and minor groove were observed. In both grooves the C-G pairs were found to be more hydrated than A-T pairs. The "spine of hydration" in the minor groove was clear. Water residence times are longer in the minor groove than in the major groove, although relatively short in both cases. No special long values are observed for sites where water molecules were observed by x-ray diffraction, indicating that water molecules having a high probability of being located in a specific site are also fast-exchanging.     

Structural characteristics of hydration sites in lysozyme

A new method is presented for determining the hydration site of proteins, where the effect of structural fluctuations in both protein and hydration water is explicitly considered by using molecular dynamics simulation (MDS). The whole hydration sites (HS) of lysozyme are composed of 195 single HSs and 38 clustered ones (CHS), and divided into 231 external HSs (EHS) and 2 internal ones (IHS). The largest CHSs, ‘Hg’ and ‘Lβ’, are the IHSs having 2.54 and 1.35 mean internal hydration waters respectively. The largest EHS, ‘Clft’, is located in the cleft region. The real hydration structure of a CHS is an ensemble of multiple structures. The transition between two structures occurs through recombinations of some H-bonds. The number of the experimental X-ray crystal waters is nearly the same as that of the estimated MDS hydration waters for 70% of the HSs, but significantly different for the rest of HSs.     

Role of hydration in collagen recognition by bacterial adhesins

Protein-protein recognition regulates the vast majority of physiological or pathological processes. We investigated the role of hydration in collagen recognition by bacterial adhesin CNA by means of first principle molecular-dynamics samplings. Our characterization of the hydration properties of the isolated partners highlights dewetting-prone areas on the surface of CNA that closely match the key regions involved in hydrophobic intermolecular interactions upon complex formation, suggesting that the hydration state of the ligand-free CNA predisposes the protein to the collagen recognition. Moreover, hydration maps of the CNA-collagen complex reveal the presence of a number of structured water molecules that mediate intermolecular interactions at the interface between the two proteins. These hydration sites feature long residence times, significant binding free energies, and a geometrical distribution that closely resembles the hydration pattern of the isolated collagen triple helix. These findings are striking evidence that CNA recognizes the collagen triple helix as a hydrated molecule. For this structural motif, the exposure of several unsatisfied backbone carbonyl groups results in a strong interplay with the solvent, which is shown to also play a role in collagen recognition.