Hydration-coupled dynamics in proteins studied by neutron scattering and NMR: the case of the typical EF-hand calcium-binding parvalbumin

The influence of hydration on the internal dynamics of a typical EF-hand calciprotein, parvalbumin, was investigated by incoherent quasi-elastic neutron scattering (IQNS) and solid-state 13C-NMR spectroscopy using the powdered protein at different hydration levels. Both approaches establish an increase in protein dynamics upon progressive hydration above a threshold that only corresponds to partial coverage of the protein surface by the water molecules. Selective motions are apparent by NMR in the 10-ns time scale at the level of the polar lysyl side chains (externally located), as well as of more internally located side chains (from Ala and Ile), whereas IQNS monitors diffusive motions of hydrogen atoms in the protein at time scales up to 20 ps. Hydration-induced dynamics at the level of the abundant lysyl residues mainly involve the ammonium extremity of the side chain, as shown by NMR. The combined results suggest that peripheral water-protein interactions influence the protein dynamics in a global manner. There is a progressive induction of mobility at increasing hydration from the periphery toward the protein interior. This study gives a microscopic view of the structural and dynamic events following the hydration of a globular protein.     

Dynamic transition associated with the thermal denaturation of a small Beta protein

We studied the temperature dependence of the picosecond internal dynamics of an all-beta protein, neocarzinostatin, by incoherent quasielastic neutron scattering. Measurements were made between 20 degrees C and 71 degrees C in heavy water solution. At 20 degrees C, only 33% of the nonexchanged hydrogen atoms show detectable dynamics, a number very close to the fraction of protons involved in the side chains of random coil structures, therefore suggesting a rigid structure in which the only detectable diffusive movements are those involving the side chains of random coil structures. At 61.8 degrees C, although the protein structure is still native, slight dynamic changes are detected that could reflect enhanced backbone and beta-sheet side-chain motions at this higher temperature. Conversely, all internal dynamics parameters (amplitude of diffusive motions, fraction of immobile scatterers, mean-squared vibration amplitude) rapidly change during heat-induced unfolding, indicating a major loss of rigidity of the beta-sandwich structure. The number of protons with diffusive motion increases markedly, whereas the volume occupied by the diffusive motion of protons is reduced. At the half-transition temperature (T = 71 degrees C) most of backbone and beta-sheet side-chain hydrogen atoms are involved in picosecond dynamics.     

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 immobilized and native Escherichia coli dihydrofolate reductase by quasielastic neutron scattering

The internal dynamics of native and immobilized Escherichia coli dihydrofolate reductase (DHFR) have been examined using incoherent quasielastic neutron scattering. These results reveal no difference between the high frequency vibration mean-square displacement of the native and the immobilized E. coli DHFR. However, length-scale-dependent, picosecond dynamical changes are found. On longer length scales, the dynamics are comparable for both DHFR samples. On shorter length scales, the dynamics is dominated by local jump motions over potential barriers. The residence time for the protons to stay in a potential well is tau = 7.95 +/- 1.02 ps for the native DHFR and tau = 20.36 +/- 1.80 ps for the immobilized DHFR. The average height of the potential barrier to the local motions is increased in the immobilized DHFR, and may increase the activation energy for the activity reaction, decreasing the rate as observed experimentally. These results suggest that the local motions on the picosecond timescale may act as a lubricant for those associated with DHFR activity occurring on a slower millisecond timescale. Experiments indicate a significantly slower catalytic reaction rate for the immobilized E. coli DHFR. However, the immobilization of the DHFR is on the exterior of the enzyme and essentially distal to the active site, thus this phenomenon has broad implications for the action of drugs distal to the active site.     

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.     

Picosecond molecular motions in bacteriorhodopsin from neutron scattering.

The characteristics of internal molecular motions of bacteriorhodopsin in the purple membrane have been studied by quasielastic incoherent neutron scattering. Because of the quasihomogeneous distribution of hydrogen atoms in biological molecules, this technique enables one to study a wide variety of intramolecular motions, especially those occurring in the picosecond to nanosecond time scale. We performed measurements at different energy resolutions with samples at various hydration levels within a temperature range of 10-300 K. The analysis of the data revealed a dynamical transition at temperatures Td between 180 K and 220 K for all motions resolved at time scales ranging from 0.1 to a few hundred picoseconds. Whereas below Td the motions are purely vibrational, they are predominantly diffusive above Td, characterized by an enormously broad distribution of correlation times. The variation of the hydration level, on the other hand, mainly affects motions slower than a few picoseconds.     

Probing the flexibility of the bacterial reaction center: the wild-type protein is more rigid than two site-specific mutants

Experimental and theoretical studies have stressed the importance of flexibility for protein function. However, more local studies of protein dynamics, using temperature factors from crystallographic data or elastic models of protein mechanics, suggest that active sites are among the most rigid parts of proteins. We have used quasielastic neutron scattering to study the native reaction center protein from the purple bacterium Rhodobacter sphaeroides, over a temperature range of 4-260 K, in parallel with two nonfunctional mutants both carrying the mutations L212Glu/L213Asp --> Ala/Ala (one mutant carrying, in addition, the M249Ala --> Tyr mutation). The so-called dynamical transition temperature, Td, remains the same for the three proteins around 230 K. Below Td the mean square displacement, u2, and the dynamical structure factor, S(Q,omega), as measured respectively by backscattering and time-of-flight techniques are identical. However, we report that above Td, where anharmonicity and diffusive motions take place, the native protein is more rigid than the two nonfunctional mutants. The higher flexibility of both mutant proteins is demonstrated by either their higher u2 values or the notable quasielastic broadening of S(Q,omega) that reveals the diffusive nature of the motions involved. Remarkably, we demonstrate here that in proteins, point genetic mutations may notably affect the overall protein dynamics, and this effect can be quantified by neutron scattering. Our results suggest a new direction of investigation for further understanding of the relationship between fast dynamics and activity in proteins. Brownian dynamics simulations we have carried out are consistent with the neutron experiments, suggesting that a rigid core within the native protein is specifically softened by distant point mutations. L212Glu, which is systematically conserved in all photosynthetic bacteria, seems to be one of the key residues that exerts a distant control over the rigidity of the core of the protein.     

The effect of hydration on protein flexibility in photosystem II of green plants studied by quasielastic neutron scattering

The effect of hydration on protein dynamics in photosystem II (PS II) membrane fragments from spinach has been investigated by using the method of quasielastic neutron scattering (QENS) at room temperature. The QENS data obtained indicate that the protein dynamics is strongly dependent on the extent of hydration. In particular, the hydration-induced activation of localized diffusive protein motions and Q_A⁻ reoxidation by Q_B in PS II appear to be correlated in their onset at a hydration value of about 45% relative humidity (r.h.). These findings underline the crucial functional relevance of localized diffusive protein motions on the picosecond-timescale for the reactions of light-induced photosynthetic water splitting under formation of plastoquinol and molecular oxygen in PS II of green plants. 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.     

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.     

From Powder to Solution: Hydration Dependence of Human Hemoglobin Dynamics Correlated to Body Temperature

A transition in hemoglobin (Hb), involving partial unfolding and aggregation, has been shown previously by various biophysical methods. The correlation between the transition temperature and body temperature for Hb from different species, suggested that it might be significant for biological function. To focus on such biologically relevant human Hb dynamics, we studied the protein internal picosecond motions as a response to hydration, by elastic and quasielastic neutron scattering. Rates of fast diffusive motions were found to be significantly enhanced with increasing hydration from fully hydrated powder to concentrated Hb solution. In concentrated protein solution, the data showed that amino acid side chains can explore larger volumes above body temperature than expected from normal temperature dependence. The body temperature transition in protein dynamics was absent in fully hydrated powder, indicating that picosecond protein dynamics responsible for the transition is activated only at a sufficient level of hydration. A collateral result from the study is that fully hydrated protein powder samples do not accurately describe all aspects of protein picosecond dynamics that might be necessary for biological function.     

Temperature- and hydration-dependent protein dynamics in photosystem II of green plants studied by quasielastic neutron scattering

Protein dynamics in hydrated and vacuum-dried photosystem II (PS II) membrane fragments from spinach has been investigated by quasielastic neutron scattering (QENS) in the temperature range between 5 and 300 K. Three distinct temperature ranges can be clearly distinguished by active type(s) of protein dynamics: (A) At low temperatures (T < 120 K), the protein dynamics of both dry and hydrated PS II is characterized by harmonic vibrational motions. (B) In the intermediate temperature range (120 < T < 240 K), the total mean square displacement total slightly deviates from the predicted linear behavior. The QENS data indicate that this deviation, which is virtually independent of the extent of hydration, is due to a partial onset of diffusive protein motions. (C) At temperatures above 240 K, the protein flexibility drastically changes because of the onset of diffusive (large-amplitude) protein motions. This dynamical transition is clearly hydration-dependent since it is strongly suppressed in dry PS II. The thermally activated onset of protein flexibility as monitored by QENS is found to be strictly correlated with the temperature-dependent increase of the electron transport efficiency from Q(A)(-) to QB (Garbers et al. (1998) Biochemistry 37, 11399-11404). Analogously, the freezing of protein mobility by dehydration in dry PS II appears to be responsible for the blockage of Q(A)(-) reoxidation by Q(B) at hydration values lower than 45% r.h. (Kaminskaya et al. (2003) Biochemistry 42, 8119-8132). Similar effects were observed for reactions of the water-oxidizing complex as outlined in the Discussion section.     

Picosecond dynamical changes on denaturation of yeast phosphoglycerate kinase revealed by quasielastic neutron scattering

Quasielastic neutron scattering experiments performed on yeast phosphoglycerate kinase in the native form and denatured in 1.5 M guanidinium chloride reveal a change in the fast (picosecond time scale) diffusive internal dynamics of the protein. The momentum and energy transfer dependences of the scattering for both states are fitted by an analytical model in which, on the experimentally accessible picosecond time scale and angstrom length scale, the dynamics of a fraction of the nonexchangeable hydrogens in the protein is described as a superposition of vibrations with uniform diffusion in a sphere, the rest of the hydrogens undergoing only vibrational motion. The fraction diffusing changes, from ≈60% in the native protein to ≈82% in the denatured protein. The radius of the sphere also changes slightly, from ≈1.8 Å in the native protein to ≈2.2 Å in the denatured protein. Possible implications of these results for the general protein folding problem are discussed. Proteins 28:380–387, 1997 © 1997 Wiley-Liss, Inc.     

In vivo measurement of internal and global macromolecular motions in Escherichia coli

We present direct quasielastic neutron scattering measurements, in vivo, of macromolecular dynamics in Escherichia coli. The experiments were performed on a wide range of timescales to cover the large panel of internal and self-diffusion motions. Three major internal processes were extracted at physiological temperature: a fast picosecond process that corresponded to restricted jump diffusion motions and two slower processes that resulted from reorientational motions occurring in ∼40 ps and 90 ps, respectively. The analysis of the fast process revealed that the cellular environment leads to an appreciable increase in internal molecular flexibility and diffusive motion rates compared with those evaluated in fully hydrated powders. The result showed that the amount of cell water plays a decisive role in internal molecular dynamics. Macromolecular interactions and confinement, however, attenuate slightly the lubricating effect of water, as revealed by the decrease of the in vivo parameters compared with those measured in solution. The study demonstrated that standard sample preparations do not mimic accurately the physiological environment and suggested that intracellular complexity participates in functional dynamics necessary for biological activity. Furthermore, the method allowed the extraction of the self-diffusion of E. coli macromolecules, which presented similar parameters as those extracted for hemoglobin in red blood cells.     

Temperature dependence of dynamics of hydrated myoglobin. Comparison of force field calculations with neutron scattering data

Molecular dynamics is used to probe the atomic motions of the carboxy-myoglobin protein as a function of temperature. Simulations of 150 picoseconds in length are carried out on the protein at 20, 60, 100, 180, 220, 240, 260, 280, 300, 320 and 340 K. The simulations attempt to mimic neutron scattering experiments very closely by including a partial hydration shell around the protein. Theoretical elastic, quasielastic and inelastic neutron scattering data are derived from the trajectories and directly compared with experiment. Compared to experiment, the simulation-derived elastic scattering curves show a decrease in intensity as a function of the scattering wavevector, q2. The inelastic and quasielastic spectra show that the inelastic peak is shifted to lower frequency than the experimental value, while quasielastic behavior is in good agreement with experiment. This suggests that the theoretical model is too flexible in the harmonic limit (low temperature), but accurately reproduces high-temperature behavior. Time correlation functions of the intermediate scattering function are determined. At low temperature there is one fast decay process, and at high temperatures there is an additional slow relaxation process that is due to quasielastic scattering. The average atomic fluctuations show that the protein behaves harmonically at low temperatures. At approximately 210 K, a glass-like transition in atomic fluctuations is seen. Above the transition temperature, the atomic fluctuations exhibit both harmonic and anharmonic behavior. Comparison of protein mobility behavior with experiment indicate the fluctuations derived from simulations are larger in the harmonic region. However, the anharmonic region agrees very well with experiment. The anharmonicity is large at all temperatures, with a gradual monotonic increase from 0.5 at 20 K to greater than 0.7 at 340 K without a noticeable change at the glass transition temperature. Heavy-atom dihedral transitions are monitored as a function of temperature. Trends in the type of dihedral transitions that occur with temperature are clearly visible. Dihedral transitions involving backbone atoms occur only above the glass transition temperature. The overall protein behavior results suggest that at low temperatures there is purely vibrational motion with one fast decay process, and above the glass transition temperature there is more anharmonic motion with a fast and a slower relaxation process occurring simultaneously.(ABSTRACT TRUNCATED AT 400 WORDS)     

Self-similar dynamics of proteins under hydrostatic pressure—Computer simulations and experiments

Different experimental techniques, such as kinetic studies of ligand binding and fluorescence correlation spectroscopy, have revealed that the diffusive, internal dynamics of proteins exhibits autosimilarity on the time scale from microseconds to hours. Computer simulations have demonstrated that this type of dynamics is already established on the much shorter nanosecond time scale, which is also covered by quasielastic neutron scattering experiments. The autosimilarity of protein dynamics is reflected in long-time memory effects in the underlying diffusion processes, which lead to a non-exponential decay of the observed time correlation functions. Fractional Brownian dynamics is an empirical model which is able to capture the essential aspects of internal protein dynamics. Here we give a brief introduction into the theory and show how the model can be used to interpret neutron scattering experiments and molecular dynamics simulation of proteins in solution under hydrostatic pressure.     

Neutron spin echo studies on ferritin: free-particle diffusion and interacting solutions

The dynamics of proteins are often studied by means of quasielastic neutron scattering (QENS), for example by time-of-flight methods. The spatial dimensions (10-20 nm) present in protein solutions are accessible by neutron scattering. In this article, a systematic study of diffusive dynamics of ferritin and apoferritin (=ferritin without iron core) is presented. Apoferritin consists of a spherical shell built of 24 protein units and carries net negative charge at pH 5. We have studied diffusive dynamics of ferritin solutions by neutron spin echo (NSE). We pay attention to an important feature of this technique compared to other QENS methods, which being the usage of a broad wavelength band. Using a more sophisticated fit function than usually used in NSE, we find as expected in low concentrated systems that the diffusion coefficient approaches the free-particle value of apoferritin and coincides with the diameter of the apoferritin shell (12.2 nm). In interacting solutions, the NSE results reveal that the dynamic picture of this complex liquid is dominated by slowing down of the dynamics. In low-salt solutions, a structure factor peak appears due to ordering of the ferritin molecules on the length scale of several intermolecular distances. We discuss the usage of different NSE fit functions for interacting solutions near the structure factor peak. Comparison of the dependence of elastic and dynamic data on the scattering vector value shows the influence of indirect interactions on the dynamic picture, irrespective of the way of data analysis, which being necessary due to the broad wavelength spectrum.     

Illuminating the Off-Pathway Nature of the Molten Globule Folding Intermediate of an α-β Parallel Protein

Partially folded protein species transiently form during folding of most proteins. Often, these species are molten globules, which may be on- or off-pathway to the native state. Molten globules are ensembles of interconverting protein conformers that have a substantial amount of secondary structure, but lack virtually all tertiary side-chain packing characteristics of natively folded proteins. Due to solvent-exposed hydrophobic groups, molten globules are prone to aggregation, which can have detrimental effects on organisms. The molten globule observed during folding of the 179-residue apoflavodoxin from Azotobacter vinelandii is off-pathway, as it has to unfold before native protein can form. Here, we study folding of apoflavodoxin and characterize its molten globule using fluorescence spectroscopy and Förster Resonance Energy Transfer (FRET). Apoflavodoxin is site-specifically labeled with fluorescent donor and acceptor dyes, utilizing dye-inaccessibility of Cys69 in cofactor-bound protein. Donor (i.e., Alexa Fluor 488) is covalently attached to Cys69 in all apoflavodoxin variants used. Acceptor (i.e., Alexa Fluor 568) is coupled to Cys1, Cys131 and Cys178, respectively. Our FRET data show that apoflavodoxin’s molten globule forms in a non-cooperative manner and that its N-terminal 69 residues fold last. In addition, striking conformational differences between molten globule and native protein are revealed, because the inter-label distances sampled in the 111-residue C-terminal segment of the molten globule are shorter than observed for native apoflavodoxin. Thus, FRET sheds light on the off-pathway nature of the molten globule during folding of an α-β parallel protein.     

Structural equilibrium fluctuations in mesophilic and thermophilic alpha-amylase

By comparing a mesophilic alpha-amylase with its thermophilic homolog, we investigated the relationship between thermal stability and internal equilibrium fluctuations. Fourier transform infrared spectroscopy monitoring hydrogen/deuterium (H/D) exchange kinetics and incoherent neutron scattering measuring picosecond dynamics were used to study dynamic features of the folded state at room temperature. Fairly similar rates of slowly exchanging amide protons indicate about the same free energy of stabilization DeltaG(stab) for both enzymes at room temperature. With respect to motions on shorter time scales, the thermophilic enzyme is characterized by an unexpected higher structural flexibility as compared to the mesophilic counterpart. In particular, the picosecond dynamics revealed a higher degree of conformational freedom for the thermophilic alpha-amylase. The mechanism proposed for increasing thermal stability in the present case is characterized by entropic stabilization and by flattening of the curvature of DeltaG(stab) as a function of temperature.     

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.