Hydration dynamics near a model protein surface

The evolution of water dynamics from dilute to very high concentration solutions of a prototypical hydrophobic amino acid with its polar backbone, N-acetyl-leucine-methylamide (NALMA), is studied by quasi-elastic neutron scattering (QENS) and molecular dynamics (MD) simulation for both the completely deuterated and completely hydrogenated leucine monomer. The NALMA-water system and the QENS data together provide a unique study for characterizing the dynamics of different hydration layers near a prototypical hydrophobic side chain and the backbone of which it is attached. We observe several unexpected features in the dynamics of these biological solutions under ambient conditions. The NALMA dynamics shows evidence of de Gennes narrowing, an indication of coherent long timescale structural relaxation dynamics. The translational and rotational water dynamics at the highest solute concentrations are found to be highly suppressed as characterized by long residential time and slow diffusion coefficients. The analysis of the more dilute concentration solutions models the first hydration shell with the 2.0 M spectra. We find that for outer layer hydration dynamics that the translational diffusion dynamics is still suppressed, although the rotational relaxation time and residential time are converged to bulk-water values. Molecular dynamics analysis of the first hydration shell water dynamics shows spatially heterogeneous water dynamics, with fast water motions near the hydrophobic side chain, and much slower water motions near the hydrophilic backbone. We discuss the hydration dynamics results of this model protein system in the context of protein function and protein-protein recognition.     

Influence of water clustering on the dynamics of hydration water at the surface of a lysozyme

Dynamics of hydration water at the surface of a lysozyme molecule is studied by computer simulations at various hydration levels in relation with water clustering and percolation transition. Increase of the translational mobility of water molecules at the surface of a rigid lysozyme molecule upon hydration is governed by the water-water interactions. Lysozyme dynamics strongly affect translational motions of water and this dynamic coupling is maximal at hydration levels, corresponding to the formation of a spanning water network. Anomalous diffusion of hydration water does not depend on hydration level up to monolayer coverage and reflects spatial disorder. Rotational dynamics of water molecules show stretched exponential decay at low hydrations. With increasing hydration, we observe appearance of weakly bound water molecules with bulklike rotational dynamics, whose fraction achieves 20-25% at the percolation threshold.     

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 structure and dynamics of B- and Z-DNA in the presence of counterions via molecular dynamics simulations

Following our previous attempts at understanding the structural and dynamical properties of water and counterions hydrating nucleic acids, we have performed molecular dynamics simulations for B- and Z-DNA. In these simulations, the nucleic acids were held rigid. In the case of B-DNA, one turn of B-DNA double helix was considered in the presence of 1500 water molecules and 20 counterions (K⁺). The simulations were performed for 4.0 ps after equilibrating the system. For Z-DNA, we considered one turn of the double helix in the presence of 1851 water molecules and 24 counterions (K⁺). The simulations were carried out for 3.5 ps after equilibration. The average temperature of these simulations was ～ 360 K for Z-DNA and ～ 345 K for B-DNA. In these simulations the hydrogen atoms were explicitly taken into account. For both simulations, a fifth-order predictor-corrector was used for solving the translational equations of motion. The rotational motion of the water molecules was represented in terms of quaternion algebra and the rotational equations of motion were solved with a second-order quaternion method using a sixth-order predictor-corrector method. A time step of 0.5 · 10^?15 s was used in these simulations. The structural and the dynamical properties of water solvating the counterions, and the phosphate groups of the DNA, were computed to understand the hydration structure. Diffusion coefficients and velocity correlation functions were calculated for both ions and the water molecules. The velocity correlation functions for the ions exhibit a caged behavior. The dipole correlation functions for the water molecules indicate that the water molecules close to the helix retain the memory of their initial orientations for longer periods of time than those away from the helix. During the time period of our simulation (3–4 ps) the ion probability distributions show a well-defined pattern and suggest limited mobility for the ions, being close to the helix.     

A model for water motion in crystals of lysozyme based on an incoherent quasielastic neutron-scattering study

This paper reports an incoherent quasielastic neutron scattering study of the single particle, diffusive motions of water molecules surrounding a globular protein, the hen egg-white lysozyme. For the first time such an analysis has been done on protein crystals. It can thus be directly related and compared with a recent structural study of the same sample. The measurement temperature ranged from 100 to 300 K, but focus was on the room temperature analysis. The very good agreement between the structural and dynamical studies suggested a model for the dynamics of water in triclinic crystals of lysozyme in the time range approximately 330 ps and at 300 K. Herein, the dynamics of all water molecules is affected by the presence of the protein, and the water molecules can be divided into two populations. The first mainly corresponds to the first hydration shell, in which water molecules reorient themselves fivefold to 10-fold slower than in bulk solvent, and diffuse by jumps from hydration site to hydration site. The long-range diffusion coefficient is five to sixfold less than for bulk solvent. The second group corresponds to water molecules further away from the surface of the protein, in a second incomplete hydration layer, confined between hydrated macromolecules. Within the time scale probed they undergo a translational diffusion with a self-diffusion coefficient reduced approximately 50-fold compared with bulk solvent. As protein crystals have a highly crowded arrangement close to the packing of macromolecules in cells, our conclusion can be discussed with respect to solvent behavior in intracellular media: as the mobility is highest next to the surface, it suggests that under some crowding conditions, a two-dimensional motion for the transport of metabolites can be dominant.     

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.     

Hydration dependent dynamics in RNA

The essential role played by local and collective motions in RNA function has led to a growing interest in the characterization of RNA dynamics. Recent investigations have revealed that even relatively simple RNAs experience complex motions over multiple time scales covering the entire ms–ps motional range. In this work, we use deuterium solid-state NMR to systematically investigate motions in HIV-1 TAR RNA as a function of hydration. We probe dynamics at three uridine residues in different structural environments ranging from helical to completely unrestrained. We observe distinct and substantial changes in ²H solid-state relaxation times and lineshapes at each site as hydration levels increase. By comparing solid-state and solution state ¹³C relaxation measurements, we establish that ns–μs motions that may be indicative of collective dynamics suddenly arise in the RNA as hydration reaches a critical point coincident with the onset of bulk hydration. Beyond that point, we observe smaller changes in relaxation rates and lineshapes in these highly hydrated solid samples, compared to the dramatic activation of motion occurring at moderate hydration.     

Evolution of the internal dynamics of two globular proteins from dry powder to solution.

Myoglobin and lysozyme picosecond internal dynamics in solution is compared to that in hydrated powders by quasielastic incoherent neutron scattering. This technique is sensitive to the motions of the nonexchangeable hydrogen atoms in a sample. Because these are homogeneously distributed throughout the protein structure, the average dynamics of the protein is described. We first propose an original data treatment to deal with the protein global motions in the case of solution samples. The validity of this treatment is checked by comparison with classical measurements of the diffusion constants. The evolution with the scattering vector of the width and relative contribution of the quasielastic component was then used to derive information on the amount of local diffusive motions and their characteristic average relaxation time. From dry powder to coverage by one water layer, the surface side chains progressively acquire the possibility to diffuse locally. On subsequent hydration, the main effect of water is to improve the rate of these diffusive motions. Motions with higher average amplitude occur in solution, about three times more than for a hydrated powder at complete coverage, with a shorter average relaxation time, approximately 4.5 ps compared to 9.4 ps for one water monolayer.     

Protein-bound water molecule counting by resolution of (1)H spin-lattice relaxation mechanisms

Water proton spin-lattice relaxation is studied in dilute solutions of bovine serum albumin as a function of magnetic field strength, oxygen concentration, and solvent deuteration. In contrast to previous studies conducted at high protein concentrations, the observed relaxation dispersion is accurately Lorentzian with an effective correlation time of 41 +/- 3 ns when measured at low proton and low protein concentrations to minimize protein aggregation. Elimination of oxygen flattens the relaxation dispersion profile above the rotational inflection frequency, nearly eliminating the high field tail previously attributed to a distribution of exchange times for either whole water molecules or individual protons at the protein-water interface. The small high-field dispersion that remains is attributed to motion of the bound water molecules on the protein or to internal protein motions on a time scale of order one ns. Measurements as a function of isotope composition permit separation of intramolecular and intermolecular relaxation contributions. The magnitude of the intramolecular proton-proton relaxation rate constant is interpreted in terms of 25 +/- 4 water molecules that are bound rigidly to the protein for a time long compared with the rotational correlation time of 42 ns. This number of bound water molecules neglects the possibility of local motions of the water in the binding site; inclusion of these effects may increase the number of bound water molecules by 50%.     

H and C NMR studies on the temperature-dependent water and protein dynamics in hydrated elastin,myoglobin and collagen

²H NMR spin-lattice relaxation and line-shape analyses are performed to study the temperature-dependent dynamics of water in the hydration shells of myoglobin, elastin, and collagen. The results show that the dynamical behaviors of the hydration waters are similar for these proteins when using comparable hydration levels of h = 0.25–0.43. Since water dynamics is characterized by strongly nonexponential correlation functions, we use a Cole–Cole spectral density for spin-lattice relaxation analysis, leading to correlation times, which are in nice agreement with results for the main dielectric relaxation process observed for various proteins in the literature. The temperature dependence can roughly be described by an Arrhenius law, with the possibility of a weak crossover in the vicinity of 220 K. Near ambient temperatures, the results substantially depend on the exact shape of the spectral density so that deviations from an Arrhenius behavior cannot be excluded in the high-temperature regime. However, for the studied proteins, the data give no evidence for the existence of a sharp fragile-to-strong transition reported for lysozyme at about 220 K. Line-shape analysis reveals that the mechanism for the rotational motion of hydration waters changes in the vicinity of 220 K. For myoglobin, we observe an isotropic motion at high temperatures and an anisotropic large-amplitude motion at low temperatures. Both mechanisms coexist in the vicinity of 220 K. ¹³C CP MAS spectra show that hydration results in enhanced elastin dynamics at ambient temperatures, where the enhancement varies among different amino acids. Upon cooling, the enhanced mobility decreases. Comparison of ²H and ¹³C NMR data reveals that the observed protein dynamics is slower than the water dynamics.     

Hydration effects on dynamics of polyglycine and sodium poly(L-glutamate).

Solid state nmr methods were applied to the study of the motions and structural heterogeneity in polyglycine, sodium poly(L-glutamate), and poly(L-alanine). The response of both the main-chain and side-chain resonances to the addition of water was studied using static and magic-angle sample-spinning line shapes as well as the carbon spin-lattice relaxation times, the proton spin-lattice relaxation time in the rotating frame, and the proton-carbon cross-polarization time. The polyglycine motions are not drastically affected by the addition of water when the polymer is in the 3(1)-helix or the beta-sheet structure. The sodium poly(L-glutamate), however, responds to increased hydration with little motion in the main-chain carbon atoms, but considerable flexibility of the side-chain atoms. The greatest motions are reported for the C5 carbon with rotational amplitudes about the C4-C5 bond of about of 50 degrees. In addition, motions somewhat less than half this size are required closer to the main chain.     

A Direct Coupling between Global and Internal Motions in a Single Domain Protein? MD Investigation of Extreme Scenarios

Proteins are not rigid molecules, but exhibit internal motions on timescales ranging from femto- to milliseconds and beyond. In solution, proteins also experience global translational and rotational motions, sometimes on timescales comparable to those of the internal fluctuations. The possibility that internal and global motions may be directly coupled has intriguing implications, given that enzymes and cell signaling proteins typically associate with binding partners and cellular scaffolds. Such processes alter their global motion and may affect protein function. Here, we present molecular dynamics simulations of extreme case scenarios to examine whether a possible relationship exists. In our model protein, a ubiquitin-like RhoGTPase binding domain of plexin-B1, we removed either internal or global motions. Comparisons with unrestrained simulations show that internal and global motions are not appreciably coupled in this single-domain protein. This lack of coupling is consistent with the observation that the dynamics of water around the protein, which is thought to permit, if not stimulate, internal dynamics, is also largely independent of global motion. We discuss implications of these results for the structure and function of proteins.     

Dynamics of tRNA at Different Levels of Hydration

The influence of hydration on the nanosecond timescale dynamics of tRNA is investigated using neutron scattering spectroscopy. Unlike protein dynamics, the dynamics of tRNA is not affected by methyl group rotation. This allows for a simpler analysis of the influence of hydration on the conformational motions in RNA. We find that hydration affects the dynamics of tRNA significantly more than that of lysozyme. Both the characteristic length scale and the timescale of the conformational motions in tRNA depend strongly on hydration. Even the characteristic temperature of the so-called “dynamical transition” appears to be hydration-dependent in tRNA. The amplitude of the conformational motions in fully hydrated tRNA is almost twice as large as in hydrated lysozyme. We ascribe these differences to a more open and flexible structure of hydrated RNA, and to a larger fraction and different nature of hydrophilic sites. The latter leads to a higher density of water that makes the biomolecule more flexible. All-atom molecular-dynamics simulations are used to show that the extent of hydration is greater in tRNA than in lysozyme. We propose that water acts as a “lubricant” in facilitating enhanced motion in solvated RNA molecules.     

Temperature dependence of protein-hydration hydrodynamics by molecular dynamics simulations

The dynamics of water molecules near the protein surface are different from those of bulk water and influence the structure and dynamics of the protein itself. To elucidate the temperature dependence hydration dynamics of water molecules, we present results from the molecular dynamic simulation of the water molecules surrounding two proteins (Carboxypeptidase inhibitor and Ovomucoid) at seven different temperatures (T=273 to 303 K, in increments of 5 K). Translational diffusion coefficients of the surface water and bulk water molecules were estimated from 2 ns molecular dynamics simulation trajectories. Temperature dependence of the estimated bulk water diffusion closely reflects the experimental values, while hydration water diffusion is retarded significantly due to the protein. Protein surface induced scaling of translational dynamics of the hydration waters is uniform over the temperature range studied, suggesting the importance protein-water interactions.     

Dynamics of water molecules buried in cavities of apolipoprotein E studied by molecular dynamics simulations and continuum electrostatic calculations

Molecular dynamics (MD) simulations of several nanoseconds each were used to monitor the dynamic behavior of the five crystal water molecules buried in the interior of the N-terminal domain of apolipoprotein E. These crystal water molecules are fairly well conserved in several apolipoprotein E structures, suggesting that they are not an artifact of the crystal and that they may have a structural and/or functional role for the protein. All five buried crystal water molecules leave the protein interior in the course of the longest simulations and exchange with water molecules from the bulk. The free energies of binding evaluated from the electrostatic binding free energy computed using a continuum model and estimates of the binding entropy changes represent shallow minima. The corresponding calculated residence times of the buried water molecules range from tens of picoseconds to hundreds of nanoseconds, which denote rather short times as for buried water molecules. Several water exchanges monitored in the simulations show that water molecules along the exit/entrance pathway use a relay of H bonds primarily formed with charged residues which helps either the exit or the entrance from or into the buried site. The exit/entrance of water molecules from/into the sites is permitted essentially by local motions of, at most, two side chains, indicating that, in these cases, complex correlated atomic motions are not needed to open the buried site toward the surface of the protein. This provides a possible explanation for the short residence times.     

Protein dynamics viewed by hydrogen exchange

To examine the relationship between protein structural dynamics and measurable hydrogen exchange (HX) data, the detailed exchange behavior of most of the backbone amide hydrogens of Staphylococcal nuclease was compared with that of their neighbors, with their structural environment, and with other information. Results show that H-bonded hydrogens are protected from exchange, with HX rate effectively zero, even when they are directly adjacent to solvent. The transition to exchange competence requires a dynamic structural excursion that removes H-bond protection and allows exposure to solvent HX catalyst. The detailed data often make clear the nature of the dynamic excursion required. These range from whole molecule unfolding, through smaller cooperative unfolding reactions of secondary structural elements, and down to local fluctuations that involve as little as a single peptide group or side chain or water molecule. The particular motion that dominates the exchange of any hydrogen is the one that allows the fastest HX rate. The motion and the rate it produces are determined by surrounding structure and not by nearness to solvent or the strength of the protecting H-bond itself or its acceptor type (main chain, side chain, structurally bound water). Many of these motions occur over time scales that are appropriate for biochemical function.     

Hydration dependent dynamics in sol–gel encapsulated myoglobin

In this work we study the effect of hydration on the dynamics of a protein in confined geometry, i.e. encapsulated in a porous silica matrix. Using elastic neutron scattering we investigate the temperature dependence of the mean square displacements of non-exchangeable hydrogen atoms of sol-gel encapsulated met-myoglobin. The study is extended to samples at 0.2, 0.3 and 0.5 g water/g protein fractions and comparison is made with met-myoglobin powders at the same average hydration and with a dry powder sample. Elastic data are analysed using a model of dynamical heterogeneity to take into account deviations of elastic intensity from gaussian behaviour in a large momentum transfer range and reveal a specific, model independent, effect of sol-gel confinement on protein dynamics, consisting mainly in a reduction of large-scale motions that are activated at temperatures larger than approximately 230 K. Surprisingly, the effect of confinement depends markedly on hydration and has a maximum at about 35% water/protein fraction corresponding to full first shell hydration. The presence of hydration-dependent MSD also in encapsulated met-Mb strongly supports the idea that the effect of sol-gel confinement on protein dynamics involves a modification of the structural/dynamical properties of the co-encapsulated solvent more than direct protein-matrix interactions.     

A theoretical study of rotational diffusion models for rod-shaped viruses. The influence of motion on 31P nuclear magnetic resonance lineshapes and transversal relaxation.

Information about the interaction between nucleic acids and coat proteins in intact virus particles may be obtained by studying the restricted backbone dynamics of the incapsulated nucleic acids using 31P nuclear magnetic resonance (NMR) spectroscopy. In this article, simulations are carried out to investigate how reorientation of a rod-shaped virus particle as a whole and isolated nucleic acid motions within the virion influence the 31P NMR lineshape and transversal relaxation dominated by the phosphorus chemical shift anisotropy. Two opposite cases are considered on a theoretical level. First, isotropic rotational diffusion is used as a model for mobile nucleic acids that are loosely or partially bound to the protein coat. The effect of this type of diffusion on lineshape and transversal relaxation is calculated by solving the stochastic Liouville equation by an expansion in spherical functions. Next, uniaxial rotational diffusion is assumed to represent the mobility of phosphorus in a virion that rotates as a rigid rod about its length axis. This type of diffusion is approximated by an exchange process among discrete sites. As turns out from these simulations, the amplitude and the frequency of the motion can only be unequivocally determined from experimental data by a combined analysis of the lineshape and the transversal relaxation. In the fast motional region both the isotropic and the uniaxial diffusion model predict the same transversal relaxation as the Redfield theory. For very slow motion, transversal relaxation resembles the nonexponential relaxation as observed for water molecules undergoing translational diffusion in a magnetic field gradient. In this frequency region T2e is inversely proportional to the cube root of the diffusion coefficient. In addition to the isotropic and uniaxial diffusion models, a third model is presented, in which fast restricted nucleic acid backbone motions dominating the lineshape are superimposed on a slow rotation of the virion about its length axis, dominating transversal relaxation. In an accompanying article the models are applied to the 31P NMR results obtained for bacteriophage M13 and tobacco mosaic virus.     

Femtosecond dynamics of intracellular water probed with nonlinear optical Kerr effect microspectroscopy

A nonlinear optical Kerr effect (OKE) microscope was developed and used to elucidate the ultra-fast diffusive motions of intracellular water molecules. In the OKE microscope, a pump-induced birefringence is sensed by a delayed probe pulse within a spatially confined volume that measures 0.5 microm in the lateral direction and 4.0 microm along the axial coordinate. This microscope allows the recording of time-resolved Kerr signals, which reflect the ultra-fast structural relaxation of the liquid, exclusively from intracellular aqueous domains. Because relaxation occurs on a picosecond time scale, only local diffusive motions are probed. The microscopic OKE signal is therefore insensitive to long-time-scale hindered translational motions enforced by intracellular mechanical barriers but probes the intrinsic orientational mobility of water molecules in cells instead. The Kerr response as determined from single intact mammalian cells under physiological conditions shows a structural relaxation time of 1.35 ps, which is 1.7 times slower than the Kerr decay observed in pure water. The data indicate that the mobility of water molecules in cellular domains is moderately restricted due to the high intracellular content of proteins and solutes.