Dynamics of large elongated RNA by NMR carbon relaxation

We present an NMR strategy for characterizing picosecond-to-nanosecond internal motions in uniformly 13C/15N-labeled RNAs that combines measurements of R1, R1rho, and heteronuclear 13C{1H} NOEs for protonated base (C2, C5, C6, and C8) and sugar (C1') carbons with a domain elongation strategy for decoupling internal from overall motions and residual dipolar coupling (RDC) measurements for determining the average RNA global conformation and orientation of the principal axis of the axially symmetric rotational diffusion. TROSY-detected pulse sequences are presented for the accurate measurement of nucleobase carbon R1 and R1rho rates in large RNAs. The relaxation data is analyzed using a model free formalism which takes into account the very high anisotropy of overall rotational diffusion (Dratio approximately 4.7), asymmetry of the nucleobase CSAs and noncollinearity of C-C, C-H dipolar and CSA interactions under the assumption that all interaction tensors for a given carbon experience identical isotropic internal motions. The approach is demonstrated and validated on an elongated HIV-1 TAR RNA (taum approximately 18 ns) both in free form and bound to the ligand argininamide (ARG). Results show that, while ARG binding reduces the amplitude of collective helix motions and local mobility at the binding pocket, it leads to a drastic increase in the local mobility of "spacer" bulge residues linking the two helices which undergo virtually unrestricted internal motions (S2 approximately 0.2) in the ARG bound state. Our results establish the ability to quantitatively study the dynamics of RNAs which are significantly larger and more anisotropic than customarily studied by NMR carbon relaxation.     

Implications of fast-time scale dynamics of human DNA/RNA cytosine methyltransferases (DNMTs) for protein function

The role of protein dynamics in the control of substrate recognition, catalysis, and protein–protein interactions is often underestimated. Recently, a number of studies have examined the contribution of protein dynamics to the thermodynamics of ligand binding in detail, mostly using NMR relaxation measurements and molecular dynamics (MD) simulations. The results unequivocally demonstrate that conformational dynamics play a pivotal role in the properties and functions of proteins, and ignoring this contribution is likely to lead to substantial errors when explaining the biological function of proteins and in predictions of the binding affinities of their cognate ligands. However, the details of the interplay between structure and dynamics and the way it affects the biological function of the target protein remain poorly understood. In this study, the changes in fast (picosecond-to-nanosecond time scale) dynamics of catalytic domains of four human cytosine DNA methyltransferases (DNMTs) were studied using molecular dynamics (MD) simulations. The results provide insight into the protein dynamics changes that occur upon binding of the cofactor, S-adenosylmethionine (SAM). Contrary to expectations, increased amplitude of motions of backbone amide (N–H) and terminal heavy atom (C–C) bond vectors was observed in all studied DNMTs upon binding of SAM. These results imply that the cofactor binding causes a global increase in the extent of protein dynamics in the short time scale. This global dynamic change constitutes a favourable entropic contribution to the free energy of SAM binding. These results suggest that cytosine DNA methyltransferases may exploit changes in their fast scale dynamics to reduce the entropic cost of the substrate binding.     

Backbone motions of free and pheromone-bound major urinary protein I studied by molecular dynamics simulation

Molecular motions of free and pheromone-bound mouse major urinary protein I, previously investigated by NMR relaxation, were simulated in 30 ns molecular dynamics (MD) runs. The backbone flexibility was described in terms of order parameters and correlation times, commonly used in the NMR relaxation analysis. Special attention was paid to the effect of conformational changes on the nanosecond time scale. Time-dependent order parameters were determined in order to separate motions occurring on different time scales. As an alternative approach, slow conformational changes were identified from the backbone torsion angle variances, and "conformationally filtered" order parameters were calculated for well-defined conformation states. A comparison of the data obtained for the free and pheromone-bound protein showed that some residues are more rigid in the bound form, but a larger portion of the protein becomes more flexible upon the pheromone binding. This finding is in general agreement with the NMR results. The higher flexibility observed on the fast (fs-ps) time scale was typically observed for the residues exhibiting higher conformational freedom on the ns time scale. An inspection of the hydrogen bond network provided a structural explanation for the flexibility differences between the free and pheromone-bound proteins in the simulations.     

Characterization of the dynamics of an essential helix in the U1A protein by time-resolved fluorescence measurements

The RNA recognition motif (RRM), one of the most common RNA-binding domains, recognizes single-stranded RNA. A C-terminal helix that undergoes conformational changes upon binding is often an important contributor to RNA recognition. The N-terminal RRM of the U1A protein contains a C-terminal helix (helix C) that interacts with the RNA-binding surface of a beta-sheet in the free protein (closed conformation), but is directed away from this beta-sheet in the complex with RNA (open conformation). The dynamics of helix C in the free protein have been proposed to contribute to binding affinity and specificity. We report here a direct investigation of the dynamics of helix C in the free U1A protein on the nanosecond time scale using time-resolved fluorescence anisotropy. The results indicate that helix C is dynamic on a 2-3 ns time scale within a 20 degrees range of motion. Steady-state fluorescence experiments and molecular dynamics simulations suggest that the dynamical motion of helix C occurs within the closed conformation. Mutation of a residue on the beta-sheet that contacts helix C in the closed conformation dramatically destabilizes the complex (Phe56Ala) and alters the steady-state fluorescence, but not the time-resolved fluorescence anisotropy, of a Trp in helix C. Mutation of Asp90 in the hinge region between helix C and the remainder of the protein to Ala or Gly subtly alters the dynamics of the U1A protein and destabilizes the complex. Together these results show that helix C maintains a dynamic closed conformation that is stable to these targeted protein modifications and does not equilibrate with the open conformation on the nanosecond time scale.     

13C NMR relaxation rates in the ionic liquid 1-ethyl-3-methylimidazolium butanesulfonate

A new method of obtaining molecular reorientational dynamics from 13C spin-lattice relaxation data of aromatic carbons in viscous solutions is applied to 13C relaxation data of both the cation and anion in the ionic liquid, 1-ethyl-3-methylimidazolium butanesulfonate ([EMIM]BSO3). 13C pseudorotational correlation times are used to calculate corrected maximum NOE factors from a combined isotropic dipolar and nuclear Overhauser effect (NOE) equation. These corrected maximum NOE factors are then used to determine the dipolar relaxation rate part of the total relaxation rate for each aromatic 13C nucleus in the imidazolium ring. Rotational correlation times are compared with viscosity data and indicate several [EMIM]BSO3 phase changes over the temperature range from 278 to 328 K. Modifications of the Stokes-Einstein-Debye (SED) model are used to determine molecular radii for the 1-ethyl-3-methylimidazolium cation. The Hu-Zwanzig correction yields a cationic radius that compares favorably with a DFT gas-phase calculation, B3LYP/(6-311+G(2d,p)). Chemical shift anisotropy values, Deltasigma, are obtained for the ring and immediately adjacent methylene and methyl carbons in the imidazolium cation and for the three carbon atoms nearest to the sulfonate group in the anion.     

Slow internal dynamics in proteins: application of NMR relaxation dispersion spectroscopy to methyl groups in a cavity mutant of T4 lysozyme

Recently developed carbon transverse relaxation dispersion experiments (Skrynnikov, N. R.; et al. J. Am. Chem. Soc. 2001, 123, 4556-4566) were applied to the study of millisecond to microsecond time scale motions in a cavity mutant of T4 lysozyme (L99A) using methyl groups as probes of dynamics. Protein expressed in E. coli cells with (13)CH(3)-pyruvate as the sole carbon source contained high levels of (13)C enrichment at a total of 80 Val gamma, Leu delta, Ile gamma (2), Ala beta, and Met epsilon methyl positions with little extraneous incorporation. Data for 72 methyl groups were available for analysis. Dispersion profiles with large amplitudes were measured for many of these residues and were well fit to a two-state exchange model. The interconversion rates and populations of the states, obtained from fitting relaxation dispersion profiles of each individual probe, were remarkably homogeneous and data for nearly all methyl groups in the protein could be collectively fit to a single cooperative conformational transition. The present study demonstrates the general applicability of methyl relaxation dispersion measurements for the investigation of millisecond time scale protein motions at a large number of side-chain positions. Potential artifacts associated with the experiments are described and methods to minimize their effects presented. These experiments should be particularly well suited for probing dynamics in high molecular weight systems due to the favorable NMR spectroscopic properties of methyl groups.     

Structural characterization of drawn (ethylene–vinyl alcohol) copolymer by high-resolution solid-state 13C NMR

High-resolution solid-state ¹³C NMR spectra of drawn (ethylene–vinyl alcohol) copolymer (EVOH) and their ¹³C spin-lattice relaxation times (T₁) were measured in order to elucidate the structure and dynamics of EVOH in the drawn state. It is found that the ¹³C chemical shift of the CH carbon for drawn EVOH moved upfield as compared with that of undrawn EVOH. Further, the mobility for the CH carbons of low molecular motion components are increased by drawing the EVOH sample. From these experimental results, the changes in the structure and dynamics of the drawn EVOH sample can be elucidated.     

13C NMR relaxation rates in the ionic liquid 1-methyl-3-nonylimidazolium hexafluorophosphate

A new method of obtaining molecular reorientational dynamics from 13C spin-lattice relaxation data of aromatic carbons in viscous solutions is applied to 13C relaxation data of the ionic liquid, 1-methyl-3-nonylimidazolium hexafluorophosphate ([MNIM]PF6). Spin-lattice relaxation times (13C) are used to determine pseudorotational correlation times for the [MNIM]PF6 ionic liquid. Pseudorotational correlation times are used to calculate corrected maximum NOE factors from a combined isotropic dipolar and nuclear Overhauser effect (NOE) equation. These corrected maximum NOE factors are then used to determine the dipolar relaxation rate part of the total relaxation rate for each aromatic 13C nucleus in the imidazolium ring. Rotational correlation times are compared with viscosity data and indicate several [MNIM]PF6 phase changes over the temperature range from 282 to 362 K. Modifications of the Stokes-Einstein-Debye (SED) model are used to determine molecular radii for the 1-methyl-3-nonylimidazolium cation. The Hu-Zwanzig correction yields a cationic radius that compares favorably with a DFT gas-phase calculation, B3LYP/(6-311+G(2d,p)). Chemical shift anisotropy values, delta sigma, are obtained for the ring and immediately adjacent methylene and methyl carbons in the imidazolium cation. The average delta sigma values for the imidazolium ring carbons are similar to those of pyrimidine in liquid crystal solutions.     

New probes of ligand flexibility in drug design: transferred (13)C CSA-dipolar cross-correlated relaxation at natural abundance

Understanding the impact of molecular flexibility remains an important outstanding problem in rational drug design. Toward this end, we present new NMR relaxation methods that describe ligand flexibility at the atomic level. Specifically, we measure natural abundance (13)C cross-correlated relaxation parameters for ligands in rapid exchange between the free and receptor-bound states. The rapid exchange transfers the bound state relaxation parameters to the free state, such that a comparison of relaxation rates in the absence and presence of protein receptor yields site-specific information concerning the bound ligand flexibility. We perform these measurements for aromatic carbons, which are highly prevalent in drug-like molecules and demonstrate significant cross-correlated relaxation between the (13)C-(1)H dipole-dipole (DD) and (13)C chemical shift anisotropy (CSA) relaxation mechanisms. Our use of natural abundance measurements addresses the practical difficulties of obtaining isotope-labeled ligands in pharmaceutical research settings. We demonstrate our methods on a small ligand of the 42 kDa kinase domain of the p38 MAP kinase. We show that exchange-transferred cross-correlated relaxation measurements are not only sensitive probes of bound ligand flexibility but also offer complementary advantages over standard R(1) = 1/T(1) and R(2) = 1/T(2) measurements. The ligand flexibility profiles obtained from the relaxation data can help assess the influence of dynamics on ligand potency or pharmacokinetic properties or both, and thereby include inherent molecular flexibility in drug design.     

A solid-state NMR study of the fast and slow dynamics of collagen fibrils at varying hydration levels

We report solid-state NMR investigations of the effect of temperature and hydration on the molecular mobility of collagen isolated from bovine achilles tendon. (13)C cross-polarization magic angle spinning (MAS) experiments were performed on samples at natural abundance, using NMR methods that detect motionally averaged dipolar interactions and chemical shift anisotropies and also slow reorientational processes. Fast motions with correlation times much shorter than 40 micro s scale dipolar couplings and chemical shift anisotropies of the carbon sites in collagen. These motionally averaged anisotropic interactions provide a measure of the amplitudes of the segmental motions expressed by a molecular order parameter. The data reveal that increasing hydration has a much stronger effect on the amplitude of the molecular processes than increasing temperature. In particular, the Cgamma carbons of the hydroxyproline residues exhibit a strong dependence of the amplitude of motion on the hydration level. This could be correlated with the effect of hydration on the hydrogen bonding structure in collagen, for which this residue is known to play a crucial role. The applicability of 1D MAS exchange experiments to investigate motions on the millisecond time-scale is discussed and first results are presented. Slow motions with correlation times of the order of milliseconds have also been detected for hydrated collagen.     

Modelling the active SARS-CoV-2 helicase complex as a basis for structure-based inhibitor design

The RNA helicase (non-structural protein 13, NSP13) of SARS-CoV-2 is essential for viral replication, and it is highly conserved among the coronaviridae family, thus a prominent drug target to treat COVID-19. We present here structural models and dynamics of the helicase in complex with its native substrates based on thorough analysis of homologous sequences and existing experimental structures. We performed and analysed microseconds of molecular dynamics (MD) simulations, and our model provides valuable insights to the binding of the ATP and ssRNA at the atomic level. We identify the principal motions characterising the enzyme and highlight the effect of the natural substrates on this dynamics. Furthermore, allosteric binding sites are suggested by our pocket analysis. Our obtained structural and dynamical insights are important for subsequent studies of the catalytic function and for the development of specific inhibitors at our characterised binding pockets for this promising COVID-19 drug target.

The RNA helicase (non-structural protein 13, NSP13) of SARS-CoV-2 is essential for viral replication, and it is highly conserved among the coronaviridae family, thus a prominent drug target to treat COVID-19.     

Chemical shift tensors of protonated base carbons in helical RNA and DNA from NMR relaxation and liquid crystal measurements

Knowledge of (13)C chemical shift anisotropy (CSA) tensors in nucleotide bases is important for interpretation of NMR relaxation data in terms of local dynamic properties of nucleic acids and for analysis of residual chemical shift anisotropy (RCSA) resulting from weak alignment. CSA tensors for protonated nucleic acid base carbons have been derived from measurements on a uniformly (13)C-enriched helical A-form RNA segment and a helical B-form DNA dodecamer at natural (13)C abundance. The magnitudes of the derived CSA principal values are tightly restricted by the magnetic field dependencies of the (13)C transverse relaxation rates, whereas the tensor orientation and asymmetry follow from quantitative measurements of interference between (13)C-{(1)H} dipolar and (13)C CSA relaxation mechanisms. Changes in the chemical shift between the isotropic and aligned states, Deltadelta, complement these measurements and permit cross-validation. The CSA tensors are determined from the experimental Deltadelta values and relaxation rates, under the assumption that the CSA tensor of any specific carbon in a given type of base is independent of the base position in either the RNA or DNA helix. However, the experimental data indicate that for pyrimidine C(6) carbons in A-form RNA the CSA magnitude is considerably larger than in B-form DNA. This result is supported by quantum chemical calculations and is attributed in part to the close proximity between intranucleotide C(6)H and O(5)' atoms in RNA. The magnitudes of the measured CSA tensors, on average, agree better with previous solid-state NMR results obtained on powdered nucleosides than with prior results from quantum chemical calculations on isolated bases, which depend rather strongly on the level of theory at which the calculations are carried out. In contrast, previously computed orientations of the chemical shift tensors agree well with the present experimental results and exhibit less dependence on the level of theory at which the computations are performed.     

13C spin relaxation studies of the basic pancreatic trypsin inhibitor

The longitudinal ¹³C spin relaxation times T₁ and the ¹³C{¹H} nuclear Overhauser enhancement were measured in a concentrated aqueous solution of the basic pancreatic trypsin inhibitor. The correlation time for overall rotational motions of the basic pancreatic trypsin inhibitor molecules was found to be τ_R ≈ 2 × 10^?8 s. In connection with previous ¹H n.m.r. studies of intramolecular motions of the aromatics, we were particularly interested in the correlation times τ_G for intramolecular segmental motions of the aromatic rings. The present experiments revealed no manifestation of intramolecular motions of the aromatics, indicating that τ_G ? 2 × 10^?8 s for the aromatic ring carbon atoms. On the other hand, rapid segmental motions were evidenced for the peripheral carbon atoms of aliphatic amino acid sidechains. Comparison of the ¹H and ¹³C n.m.r. data on the basic pancreatic trypsin inhibitor indicates that the time scale of high resolution ¹H n.m.r. at high fields may in many instances be more appropriate for studies of the molecular dynamics in globular proteins than the time scale of spin relaxation measurements.     

Direct detection of structurally resolved dynamics in a multiconformation receptor-ligand complex

Structure-based drug design relies on static protein structures despite significant evidence for the need to include protein dynamics as a serious consideration. In practice, dynamic motions are neglected because they are not understood well enough to model, a situation resulting from a lack of explicit experimental examples of dynamic receptor-ligand complexes. Here, we report high-resolution details of pronounced ~1 ms time scale motions of a receptor-small molecule complex using a combination of NMR and X-ray crystallography. Large conformational dynamics in Escherichia coli dihydrofolate reductase are driven by internal switching motions of the drug-like, nanomolar-affinity inhibitor. Carr-Purcell-Meiboom-Gill relaxation dispersion experiments and NOEs revealed the crystal structure to contain critical elements of the high energy protein-ligand conformation. The availability of accurate, structurally resolved dynamics in a protein-ligand complex should serve as a valuable benchmark for modeling dynamics in other receptor-ligand complexes and prediction of binding affinities.     

Probing methyl dynamics from 13C autocorrelated and cross-correlated relaxation

An understanding of side-chain motions in protein is of great interest since side chains often play an important role in protein folding and intermolecular interactions. A novel method for measuring dipole-dipole cross-correlated relaxation in methyl groups of uniformly 13C-labeled proteins without deuteration has been developed by our group. The excellent agreement between dynamic parameters of methyl groups in ubiquitin obtained from the cross-correlated relaxation and 13C spin-lattice relaxation and those derived previously from 2H relaxation data demonstrates the reliability of the method. This method was applied to the study of side-chain dynamics of human intestinal fatty acid binding protein (IFABP) with and without its ligand. Binding of oleic acid to the protein results in decreased mobility of most of the methyl groups in the binding region, whereas no significant change in mobility was observed for methyl groups in the nonbinding region.     

Temperature dependence of anisotropic protein backbone dynamics

The measurement of (15)N NMR spin relaxation, which reports the (15)N-(1)H vector reorientational dynamics, is a widely used experimental method to assess the motion of the protein backbone. Here, we investigate whether the (15)N-(1)H vector motions are representative of the overall backbone motions, by analyzing the temperature dependence of the (15)N-(1)H and (13)CO-(13)C(alpha) reorientational dynamics for the small proteins binase and ubiquitin. The latter dynamics were measured using NMR cross-correlated relaxation experiments. The data show that, on average, the (15)N-(1)H order parameters decrease only by 2.5% between 5 and 30 degrees C. In contrast, the (13)CO-(13)C(alpha) order parameters decrease by 10% over the same temperature trajectory. This strongly indicates that there are polypeptide-backbone motions activated at room temperature that are not sensed by the (15)N-(1)H vector. Our findings are at variance with the common crank-shaft model for protein backbone dynamics, which predicts the opposite behavior. This study suggests that investigation of the (15)N relaxation alone would lead to underestimation of the dynamics of the protein backbone and the entropy contained therein.     

Interpretation of spin relaxation in polyisobutylene in terms of specific motions

Porton and carbon spin-lattice relaxation times T₁ and nuclear Overhauser enhancements are interpreted in terms of motions likely in linear polyisobutylene. Most of the interpretation is based on relaxation data in the literature, but some additional ¹H and ¹³C pulse Fourier transform experiments were conducted to resolve a disagreement in the literature concerning cross relaxation between the two types of protons present in polyisobutylene. Spin relaxation in solution and the bulk is accounted for by three specific motions considered as independent sources of motional modulation of the dipole–dipole interaction. The first motion is overall isotropic rotatory diffusion which has a known dependence on molecular weight, intrinsic viscosity, and solvent viscosity for polymers in solution, and a known dependence on molecular weight and viscosity for bulk polymers. The effects of overall tumbling account for a decrease of T₁ for the methylene and methyl carbons with increasing molecular weight in solution and increase of T₁ of methylene carbons with molecular weight in bulk. The second motion considered is backbone rearrangements caused by the three-bond jump. This motion dominates relaxation of the methylene carbons either in solution or in the bulk allowing for the determination of the associated correlation time. The correlation time characterizing the occurrence of the three-bond jump in a 5% (wt/vol) solution in CCI₄ at 45°C is 58 psec, and in the bulk at 45°C it is 11 nsec. The last motion included in the model is methyl-group rotation about the threefold symmetry axis. The methyl-group rotational correlation time is 0.20 nsec in a 5% (wt/vol) solution in CCI₄ at 45°C and 0.33 nsec in the bulk at 45°C. The concentration dependence of the backbone motion contrasts strongly with the corresponding dependence of methyl-group rotation.     

The lipidated membrane anchor of full length N-Ras protein shows an extensive dynamics as revealed by solid-state NMR spectroscopy

Many proteins involved in signal transduction are equipped with covalently attached lipid chains providing a hydrophobic anchor targeting these molecules to membranes. Despite the considerable biological significance of this membrane binding mechanism for 5-10% of all cellular proteins, to date very little is known about structural and dynamical features of lipidated membrane binding domains. Here we report the first comprehensive study of the molecular dynamics of the C-terminus of membrane-associated full-length lipidated Ras protein determined by solid-state NMR. Fully functional lipid-modified N-Ras protein was obtained by chemical-biological synthesis ligating the expressed water soluble N-terminus with a chemically synthesized (2)H or (13)C labeled lipidated heptapeptide. Dynamical parameters for the lipid chain modification at Cys 181 were determined from static (2)H NMR order parameter and relaxation measurements. Order parameters describing the amplitude of motion in the protein backbone and the side chain were determined from site-specific measurements of (1)H-(13)C dipolar couplings for all seven amino acids in the membrane anchor of Ras. Finally, the correlation times of motion were determined from temperature dependent relaxation time measurements and analyzed using a modified Lipari Szabo approach. Overall, the C-terminus of Ras shows a versatile dynamics with segmental fluctuations and axially symmetric overall motions on the membrane surface. In particular, the lipid chain modifications are highly flexible in the membrane.     

Local chain motions of poly(β-hydroxyoctanoate) in the bulk. 13C NMR relaxation study

Variable-temperature ¹³C NMR spin-lattice relaxation times (T₁) and Nuclear Overhauser Enhancements (NOE) at two magnetic fields have been used to study the dynamics of the amorphous part of a semicrystalline sample (33% of crystallinity) of poly(β-hydroxyoctanoate) (PHO). The interpretation of the relaxation data of the backbone carbons was made by employing a number of motional models. Among these, the DLM model offered the best interpretation of the relaxation data in terms of conformational transitions and librational motions of the backbone C? H vectors, and proved to be superior to unimodal distribution functions. The interpretation of temperature- and frequency-dependent T₁ and NOE data of the carbon nuclei in the n-pentyl side chain was made by employing a newly developed composite spectral density function for multiple internal C? C bond rotations of restricted amplitude and chain segmental motion. The temperature dependence of the linewidths of the various protonated carbon resonances of PHO has been discussed in terms of the semicrystalline character of this polymer. © 1995 John Wiley & Sons, Inc.     

Mapping the hydration dynamics of ubiquitin

The nature of water's interaction with biomolecules such as proteins has been difficult to examine in detail at atomic resolution. Solution NMR spectroscopy is potentially a powerful method for characterizing both the structural and temporal aspects of protein hydration but has been plagued by artifacts. Encapsulation of the protein of interest within the aqueous core of a reverse micelle particle results in a general slowing of water dynamics, significant reduction in hydrogen exchange chemistry and elimination of contributions from bulk water thereby enabling the use of nuclear Overhauser effects to quantify interactions between the protein surface and hydration water. Here we extend this approach to allow use of dipolar interactions between hydration water and hydrogens bonded to protein carbon atoms. By manipulating the molecular reorientation time of the reverse micelle particle through use of low viscosity liquid propane, the T(1ρ) relaxation time constants of (1)H bonded to (13)C were sufficiently lengthened to allow high quality rotating frame nuclear Overhauser effects to be obtained. These data supplement previous results obtained from dipolar interactions between the protein and hydrogens bonded to nitrogen and in aggregate cover the majority of the molecular surface of the protein. A wide range of hydration dynamics is observed. Clustering of hydration dynamics on the molecular surface is also seen. Regions of long-lived hydration water correspond with regions of the protein that participate in molecular recognition of binding partners suggesting that the contribution of the solvent entropy to the entropy of binding has been maximized through evolution.