Hydration affects both harmonic and anharmonic nature of protein dynamics

To understand the effect of hydration on protein dynamics, inelastic neutron-scattering experiments were performed on staphylococcal nuclease samples at differing hydration levels: dehydrated, partially hydrated, and hydrated. At cryogenic temperatures, hydration affected the collective motions with energies lower than 5 meV, whereas the high-energy localized motions were independent of hydration. The prominent change was a shift of boson peak toward higher energy by hydration, suggesting a hardening of harmonic potential at local minima on the energy landscape. The 240 K transition was observed only for the hydrated protein. Significant quasielastic scattering at 300 K was observed only for the hydrated sample, indicating that the origin of the transition is the motion activated by hydration water. The neutron-scattering profile of the partially hydrated sample was quite similar to that of the hydrated sample at 100 K and 200 K, whereas it was close to the dehydrated sample at 300 K, indicating that partial hydration is sufficient to affect the harmonic nature of protein dynamics, and that there is a threshold hydration level to activate anharmonic motions. Thus, hydration water controls both harmonic and anharmonic protein dynamics by differing means.     

Temperature dependence of the low frequency dynamics of myoglobin. Measurement of the vibrational frequency distribution by inelastic neutron scattering.

Inelastic neutron scattering spectra of myoglobin hydrated to 0.33 g water (D2O)/g protein have been measured in the low frequency range (1-150 cm-1) at various temperatures between 100 and 350 K. The spectra at low temperatures show a well-resolved maximum in the incoherent dynamic structure factor Sinc(q, omega) at approximately 25 cm-1 and no elastic broadening. This maximum becomes gradually less distinct above 180 K due to the increasing amplitude of quasielastic scattering which extends out to 30 cm-1. The vibrational frequency distribution derived independently at 100 and 180 K are very similar, suggesting harmonic behavior at these temperatures. This result has been used to separate the vibrational motion from the quasielastic motion at temperatures above 180 K. The form of the density of states of myoglobin is discussed in relation to that of other amorphous systems, to theoretical calculations of low frequency modes in proteins, and to previous observations by electron-spin relaxation of fractal-like spectral properties of proteins. The onset of quasielastic scattering above 180 K is indicative of a dynamic transition of the system and correlates with an anomalous increase in the atomic mean-squared displacements observed by M?ssbauer spectroscopy (Parak, F., E. W. Knapp, and D. Kucheida. 1982. J. Mol. Biol. 161: 177-194.) and inelastic neutron scattering (Doster, W., S. Cusack, and W. Petry, 1989. Nature [Lond.]. 337: 754-756.) Similar behavior is observed for a hydrated powder of lysozyme suggesting that the low frequency dynamics of globular proteins have common features.     

Enzyme activity below the dynamical transition at 220 K.

Enzyme activity requires the activation of anharmonic motions, such as jumps between potential energy wells. However, in general, the forms and time scales of the functionally important anharmonic dynamics coupled to motion along the reaction coordinate remain to be determined. In particular, the question arises whether the temperature-dependent dynamical transition from harmonic to anharmonic motion in proteins, which has been observed experimentally and using molecular dynamics simulation, involves the activation of motions required for enzyme function. Here we present parallel measurements of the activity and dynamics of a cryosolution of glutamate dehydrogenase as a function of temperature. The dynamical atomic fluctuations faster than approximately 100 ps were determined using neutron scattering. The results show that the enzyme remains active below the dynamical transition observed at approximately 220 K, i.e., at temperatures where no anharmonic motion is detected. Furthermore, the activity shows no significant deviation from Arrhenius behavior down to 190 K. The results indicate that the observed transition in the enzyme's dynamics is decoupled from the rate-limiting step along the reaction coordinate.     

Temperature-dependent structural and functional features of a hyperthermostable enzyme using elastic neutron scattering

The dynamic behavior of an endoglucanase from the hyperthermophilic microorganism Pyrococcus furiosus was investigated using elastic neutron scattering. The temperature dependence of the atomic motions was correlated with conformational and functional characteristics of the enzyme. The onset of biological function at temperatures higher than approximately 25 degrees C (the hyperthermostable enzyme is essentially inactive at room temperature) was associated with a dynamical transition in the anharmonic motions domain. This transition from the nonactive to the enzymatically active conformation involved structurally similar conformational substates in the energy landscape. From the mean-square displacement of the protein atoms, the molecular flexibility and the effective force constants were calculated at different temperature zones. The results showed that the activity increases at higher temperatures where the intramolecular bonds are weakened and the overall rigidity of the protein is decreased. Further temperature increase resulted in significantly increased atomic fluctuations featuring heat denaturation of the protein.     

Radially softening diffusive motions in a globular protein

Molecular dynamics simulation, quasielastic neutron scattering and analytical theory are combined to characterize diffusive motions in a hydrated protein, C-phycocyanin. The simulation-derived scattering function is in approximate agreement with experiment and is decomposed to determine the essential contributions. It is found that the geometry of the atomic motions can be modeled as diffusion in spheres with a distribution of radii. The time dependence of the dynamics follows stretched exponential behavior, reflecting a distribution of relaxation times. The average side chain and backbone dynamics are quantified and compared. The dynamical parameters are shown to present a smooth variation with distance from the core of the protein. Moving outward from the center of the protein there is a progressive increase of the mean sphere size, accompanied by a narrowing and shifting to shorter times of the relaxation time distribution. This smooth, "radially softening" dynamics may have important consequences for protein function. It also raises the possibility that the dynamical or "glass" transition with temperature observed experimentally in proteins might be depth dependent, involving, as the temperature decreases, progressive freezing out of the anharmonic dynamics with increasing distance from the center of the protein.     

The low-temperature heat capacity of solid proteins.

Several harmonic models of protein fluctuations are used to calculate the heat capacity. They get the spectral density of conformational modes from inelastic neutron scattering, normal mode calculations, or macroscopic elasticity (Debye model). It is assumed that the low-frequency spectral density depends only weakly on temperature and protein species. The Debye model predicts temperatures below which modes are primarily in their ground states: 10 and 80 K for the lattice and conformational modes, respectively. The models differ most below 100 K. The mode calculations yield the most accurate predictions, though all three models are within twofold of the data. The heat capacity has the power law form aTb for T less than 30 K. The experimental b's of proteins are 1.6-1.8, and the theoretical, 1.1-1.3. One possible explanation for the discrepancy is the occurrence of transitions between discrete conformations. All of the models approach the measured data in the range 100-200 K. They are very similar above 200 K, where the heat capacity includes significant contributions from bond stretching and bending. This masks the possible anharmonic behavior of the conformational modes. Hydration substantially increases the heat capacity above 200 K. This effect seems to be a consequence of conformational transitions that have higher energy than the ones seen with low hydration. The analysis also predicts that denaturation with constant hydration produces a negligible increase of heat capacity. The larger increment in solution arises from the different hydration of the folded and unfolded states, and is responsible for the existence of cold denaturation. This phenomenon is thus predicted not to occur when the hydration is constant.     

Mean square fluctuations of hydrogen atoms and water-biopolymer interactions in hydrated saccharides.

We have used the elastic neutron scattering technique to investigate the dynamics of the two main saccharidic components of starch: amylose and amylopectin. The measurements were carried out in the temperature range of 20 to 320 K and at different hydration levels from the dry state up to 0.47 g saccharide/g D(2)O. In the dry samples, the atomic dynamics is harmonic up to approximately 300 K. In the hydrated samples a "glass-like" transition leading to an anharmonic dynamics is observed. The onset of the anharmonicity occurs at temperatures that increase from approximately 180 K to 260 K upon decreasing hydration from 0.5 to 0.1 g saccharide/g D(2)O. This behavior is qualitatively similar to that observed in hydrated globular proteins, but quantitative differences are present. Assuming a simple asymmetric double-well potential model, the temperature and hydration dependence of the transition have been described in terms of few physical parameters.     

Incoherent neutron scattering of copper azurin: a comparison with molecular dynamics simulation results

The low-frequency dynamics of copper azurin has been studied at different temperatures for a dry and deuterium hydrated sample by incoherent neutron scattering and the experimental results have been compared with molecular dynamics (MD) simulations carried out in the same temperature range. Experimental Debye-Waller factors are consistent with a dynamical transition at approximately 200 K which appears partially suppressed in the dry sample. Inelastic and quasielastic scattering indicate that hydration water modulates both vibrational and diffusive motions. The low-temperature experimental dynamical structure factor of the hydrated protein shows an excess of inelastic scattering peaking at about 3 meV and whose position is slightly shifted downwards in the dry sample. Such an excess is reminiscent of the “boson peak” observed in glass-like materials. This vibrational peak is quite well reproduced by MD simulations, although at a lower energy. The experimental quasielastic scattering of the two samples at 300 K shows a two-step relaxation behaviour with similar characteristic times, while the corresponding intensities differ only by a scale factor. Also, MD simulations confirm the two-step diffusive trend, but the slow process seems to be characterized by a decay faster than the experimental one. Comparison with incoherent neutron scattering studies carried out on proteins having different structure indicates that globular proteins display common elastic, quasielastic and inelastic features, with an almost similar hydration dependence, irrespective of their secondary and tertiary structure. Received: 12 October 1998 / Revised version: 19 February 1999 / Accepted: 1 March 1999     

Dynamics of hydrogen atoms in superoxide dismutase by quasielastic neutron scattering.

The low energy dynamic of the enzyme Cu,Zn superoxide dismutase have been investigated by means of quasielastic neutron scattering in the temperature range 4-320 K. Below 200 K the scattering is purely elastic, while above this temperature a pronounced decrease in the elastic intensity is observed, together with the onset of a small quasielastic component. This behavior is similar to that previously observed in other more flexible globular proteins, and can be attributed to transitions between slightly different conformational substates of the protein tertiary structure. The presence of only a small quasielastic component, whose intensity is < or = 25% of the total spectrum, is related to the high structural rigidity of this protein.     

Enzyme activity and flexibility at very low hydration

Recent measurements have demonstrated enzyme activity at hydrations as low as 3%. This raises the question of whether hydration-induced enzyme flexibility is important for activity. Here, to address this, picosecond dynamic neutron scattering experiments are performed on pig liver esterase powders at 0%, 3%, 12%, and 50% hydration by weight and at temperatures ranging from 120 to 300 K. At all temperatures and hydrations, significant quasielastic scattering intensity is found in the protein, indicating the presence of anharmonic, diffusive motion. As the hydration increases, a temperature-dependent dynamical transition appears and strengthens involving additional diffusive motion. The implication of these results is that, although the additional hydration-induced diffusive motion in the protein detected here may be related to increased activity, it is not required for the enzyme to function.     

Effect of the environment on the protein dynamical transition: a neutron scattering study

We performed an elastic neutron scattering investigation of the molecular dynamics of lysozyme solvated in glycerol, at different water contents h (grams of water/grams of lysozyme). The marked non-Gaussian behavior of the elastic intensity was studied in a wide experimental momentum transfer range, as a function of the temperature. The internal dynamics is well described in terms of the double-well jump model. At low temperature, the protein total mean square displacements exhibit an almost linear harmonic trend irrespective of the hydration level, whereas at the temperature T(d) a clear changeover toward an anharmonic regime marks a protein dynamical transition. The decrease of T(d) from approximately 238 K to approximately 195 K as a function of h is reminiscent of that found in the glass transition temperature of aqueous solutions of glycerol, thus suggesting that the protein internal dynamics as a whole is slave to the environment properties. Both T(d) and the total mean square displacements indicate that the protein flexibility strongly rises between 0.1 and 0.2h. This hydration-dependent dynamical activation, which is similar to that of hydrated lysozyme powders, is related to the specific interplay of the protein with the surrounding water and glycerol molecules.     

Hydration effect on low-frequency protein dynamics observed in simulated neutron scattering spectra

Hydration effects on protein dynamics were investigated by comparing the frequency dependence of the calculated neutron scattering spectra between full and minimal hydration states at temperatures between 100 and 300 K. The protein boson peak is observed in the frequency range 1-4 meV at 100 K in both states. The peak frequency in the minimal hydration state shifts to lower than that in the full hydration state. Protein motions with a frequency higher than 4 meV were shown to undergo almost harmonic motion in both states at all temperatures simulated, whereas those with a frequency lower than 1 meV dominate the total fluctuations above 220 K and contribute to the origin of the glass-like transition. At 300 K, the boson peak becomes buried in the quasielastic contributions in the full hydration state but is still observed in the minimal hydration state. The boson peak is observed when protein dynamics are trapped within a local minimum of its energy surface. Protein motions, which contribute to the boson peak, are distributed throughout the whole protein. The fine structure of the dynamics structure factor is expected to be detected by the experiment if a high resolution instrument (<∼20 μeV) is developed in the near future.     

Large-scale oscillatory motions in myoglobin

Elastic and inelastic components of neutron scattering from lyophilized myoglobin were measured at wavevector transfers, Q of 1.0, 1.4, 1.85, and 3.0 Å^-1 and at temperatures of 77 and 298 K. The elastic scattering determines an overall effective mean square atomic displacement which was employed in a simple harmonic Langevin formalism to obtain effective force constants describing molecular deformations. These force constants and the inelastic data are compatible with alpha helices as a vibrating species at room temperature, and with amino acid residues as the vibrating entities at 77 K. The present measurements are consistent with recent molecular dynamics calculations.     

Dynamics of Protein and its Hydration Water: Neutron Scattering Studies on Fully Deuterated GFP

We present a detailed analysis of the picosecond-to-nanosecond motions of green fluorescent protein (GFP) and its hydration water using neutron scattering spectroscopy and hydrogen/deuterium contrast. The analysis reveals that hydration water suppresses protein motions at lower temperatures (<∼200 K), and facilitates protein dynamics at high temperatures. Experimental data demonstrate that the hydration water is harmonic at temperatures <∼180–190 K and is not affected by the proteins’ methyl group rotations. The dynamics of the hydration water exhibits changes at ∼180–190 K that we ascribe to the glass transition in the hydrated protein. Our results confirm significant differences in the dynamics of protein and its hydration water at high temperatures: on the picosecond-to-nanosecond timescale, the hydration water exhibits diffusive dynamics, while the protein motions are localized to <∼3 Å. The diffusion of the GFP hydration water is similar to the behavior of hydration water previously observed for other proteins. Comparison with other globular proteins (e.g., lysozyme) reveals that on the timescale of 1 ns and at equivalent hydration level, GFP dynamics (mean-square displacements and quasielastic intensity) are of much smaller amplitude. Moreover, the suppression of the protein dynamics by the hydration water at low temperatures appears to be stronger in GFP than in other globular proteins. We ascribe this observation to the barrellike structure of GFP.     

Temperature dependence of protein dynamics: computer simulation analysis of neutron scattering properties

The temperature dependence of the internal dynamics of an isolated protein, bovine pancreatic trypsin inhibitor, is examined using normal mode analysis and molecular dynamics (MD) simulation. It is found that the protein exhibits marked anharmonic dynamics at temperatures of approximately 100-120 K, as evidenced by departure of the MD-derived average mean square displacement from that of the harmonic model. This activation of anharmonic dynamics is at lower temperatures than previously detected in proteins and is found in the absence of solvent molecules. The simulation data are also used to investigate neutron scattering properties. The effects are determined of instrumental energy resolution and of approximations commonly used to extract mean square displacement data from elastic scattering experiments. Both the presence of a distribution of mean square displacements in the protein and the use of the Gaussian approximation to the dynamic structure factor lead to quantified underestimation of the mean square displacement obtained.     

The dynamical transition of proteins,concepts and misconceptions

The dynamics of hydrated proteins and of protein crystals can be studied within a wide temperature range, since the water of hydration does not crystallize at low temperature. Instead it turns into an amorphous glassy state below 200 K. Extending the temperature range facilitates the spectral separation of different molecular processes. The conformational motions of proteins show an abrupt enhancement near 180 K, which has been called a "dynamical transition". In this contribution various aspects of the transition are critically reviewed: the role of the instrumental resolution function in extracting displacements from neutron elastic scattering data and the question of the appropriate dynamic model, discrete transitions between states of different energy versus continuous diffusion inside a harmonic well, are discussed. A decomposition of the transition involving two motional components is performed: rotational transitions of methyl groups and small scale librations of side-chains, induced by water at the protein surface. Both processes create an enhancement of the observed amplitude. The onset occurs, when their time scale becomes compatible with the resolution of the spectrometer. The reorientational rate of hydration water follows a super-Arrhenius temperature dependence, a characteristic feature of a dynamical transition. It occurs only with hydrated proteins, while the torsional motion of methyl groups takes place also in the dehydrated or solvent-vitrified system. Finally, the role of fast hydrogen bond fluctuations contributing to the amplitude enhancement is discussed.     

The temperature dependence of internal molecular motions in hydrated and dry alpha-amylase: the role of hydration water in the dynamical transition of proteins.

Internal molecular motions of proteins are strongly affected by environmental conditions, like temperature and hydration. As known from numerous studies, the dynamical behavior of hydrated proteins on the picosecond time scale is characterized by vibrational motions in the low-temperature regime and by an onset of stochastic large-amplitude fluctuations at a transition temperature of 180-230 K. The present study reports on the temperature dependence of internal molecular motions as measured with incoherent neutron scattering from the globular water-soluble protein alpha-amylase and from a protein-lipid complex of rhodopsin in disk membranes. Samples of alpha-amylase have been measured in a hydrated and dehydrated state. In contrast to the hydrated sample, which exhibits a pronounced dynamical transition near 200 K, the dehydrated alpha-amylase does not show an appreciable proportion of stochastic large-amplitude fluctuations and no dynamical transition in the measured temperature range of 140-300 K. The obtained results, which are compared to the dynamical behavior of protein-lipid complexes, are discussed with respect to the influence of hydration on the dynamical transition and in the framework of the glass transition.     

Dynamics of C-phycocyanin in various deuterated trehalose/water environments measured by quasielastic and elastic neutron scattering

The molecular understanding of protein stabilization by the disaccharide trehalose in extreme temperature or hydration conditions is still debated. In the present study, we investigated the role of trehalose on the dynamics of the protein C-phycocyanin (C-PC) by neutron scattering. To single out the motions of C-PC hydrogen (H) atoms in various trehalose/water environments, measurements were performed in deuterated trehalose and heavy water (D(2)O). We report that trehalose decreases the internal C-PC dynamics, as shown by a reduced diffusion coefficient of protein H atoms. By fitting the Elastic Incoherent Structure Factor-which gives access to the "geometry" of the internal proton motions-with the model of diffusion inside a sphere, we found that the presence of trehalose induces a significantly higher proportion of immobile C-PC hydrogens. We investigated, by elastic neutron scattering, the mean square displacements (MSDs) of deuterated trehalose/D(2)O-embedded C-PC as a function of temperature in the range of 40-318 K. Between 40 and approximately 225 K, harmonic MSDs of C-PC are slightly smaller in samples containing trehalose. Above a transition temperature of approximately 225 K, we observed anharmonic motions in all trehalose/water-coated C-PC samples. In the hydrated samples, MSDs are not significantly changed by addition of 15% trehalose but are slightly reduced by 30% trehalose. In opposition, no dynamical transition was detected in dry trehalose-embedded C-PC, whose hydrogen motions remain harmonic up to 318 K. These results suggest that a role of trehalose would be to stabilize proteins by inhibiting some fluctuations at the origin of protein unfolding and denaturation.     

Water-coupled low-frequency modes of myoglobin and lysozyme observed by inelastic neutron scattering.

Conformational changes of proteins often involve the relative motion of rigid structural domains. Normal mode analysis and molecular dynamics simulations of small globular proteins predict delocalized vibrations with frequencies below 20 cm(-1), which may be overdamped in solution due to solvent friction. In search of these modes, we have studied deuterium-exchanged myoglobin and lysozyme using inelastic neutron scattering in the low-frequency range at full and low hydration to modify the degree of damping. At room temperature, the hydrated samples exhibit a more pronounced quasielastic spectrum due to diffusive motions than the dehydrated samples. The analysis of the corresponding lineshapes suggests that water modifies mainly the amplitude, but not the characteristic time of fast protein motions. At low temperatures, in contrast, the dehydrated samples exhibit larger motional amplitudes than the hydrated ones. The excess scattering, culminating at 16 cm(-1), is suggested to reflect water-coupled librations of polar side chains that are depressed in the hydrated system by strong intermolecular hydrogen bonding. Both myoglobin and lysozyme exhibit ultra-low-frequency modes below 10 cm(-1) in the dry state, possibly related to the breathing modes predicted by harmonic analysis.     

Thermal fluctuations of DNA enclosed by glycerol–water glassy matrices: an elastic neutron scattering investigation

Through elastic neutron scattering measurements, we investigated the thermal fluctuations of DNA enclosed by glycerol–water glassy matrices, at different levels of hydration, over the wide temperature range from 20 to 300 K. For all the samples, the extracted hydrogen mean square displacements (MSD) show a purely vibrational harmonic trend at very low temperatures, and a first onset of anharmonic dynamics above ∼100 K. Such onset is consistent with the activation of DNA methyl group rotational motions. Then, at a certain temperature T _d, the MSD show a second onset of anharmonicity, which corresponds to the DNA dynamical transition. The T _d values vary as a function of the hydration degree of the environment. The crucial role of the solvent mobility to activate the DNA thermal fluctuations is proposed, together with a preferential hydration effect of the DNA phosphate groups. Finally, a comparison between the average mobility of homologous samples of DNA and the lysozyme protein is considered. Advanced neutron scattering and complementary techniques to study biological systems. Contributions from the meetings, “Neutrons in Biology”, STFC Rutherford Appleton Laboratory, Didcot, UK, 11–13 July and “Proteins At Work 2007”, Perugia, Italy, 28–30 May 2007.