Picosecond molecular motions in bacteriorhodopsin from neutron scattering

Probing the Hydrophobic Region of a Lipid Bilayer at Specific Depths Using Vibrational Spectroscopy

A novel spectroscopic approach for studying the flexibility and mobility in the hydrophobic interior of lipid bilayers at specific depths is proposed. A set of test compounds featuring an azido moiety and a cyano or carboxylic acid moiety, connected by an alkyl chain of different lengths, was synthesized. FTIR data and molecular dynamics calculations indicated that the test compounds in a bilayer are oriented so that the cyano or carboxylic acid moiety is located in the lipid head-group region, while the azido group stays inside the bilayer at the depth determined by its alkyl chain length. We found that the asymmetric stretching mode of the azido group (νN3) can serve as a reporter of the membrane interior dynamics. FTIR and two-dimensional infrared (2DIR) studies were performed at different temperatures, ranging from 22 to 45 °C, covering the Lβ-Lα phase transition temperature of dipalmitoylphosphatidylcholine (∼41 °C). The width of the νN3 peak was found to be very sensitive to the phase transition and to the temperature in general. We introduced an order parameter, SN3, which characterizes restrictions to motion inside the bilayer. 2DIR spectra of νN3 showed different extents of inhomogeneity at different depths in the bilayer, with the smallest inhomogeneity in the middle of the leaflet. The spectral diffusion dynamics of the N3 peak was found to be dependent on the depth of the N3 group location in the bilayer. The obtained results enhance our understanding of the bilayer dynamics and can be extended to investigate membranes with more complex compositions.

Quantifying Bacteriorhodopsin Activity as a Function of its Local Environment with a Raman-Based Assay

Bacteriorhodopsin (bR) is a transmembrane protein that functions as a light-driven proton pump in halophilic archaea. The bR photocycle has been well-characterized; however, these measurements almost exclusively measured purified bR, outside of its native membrane. To investigate what effect the cellular environment has on the bR photocycle, we have developed a Raman-based assay that can monitor the activity of the bR in a variety of conditions, including in its native membrane. The assay uses two continuous-wave lasers, one to initiate photochemistry and one to monitor bR activity. The excitation leads to the steady-state depletion of ground-state bR, which directly relates to the population of photocycle intermediate states. We have used this assay to monitor bR activity both in vitro and in vivo. Our in vitro measurements confirm that our assay is sensitive to bulk environmental changes reported in the literature. Our in vivo measurements show a decrease in bR activity with increasing extracellular pH for bR in its native membrane. The difference in activity with increasing pH indicates that the native membrane environment affects the function of bR. This assay opens the door to future measurements into understanding how the local environment of this transmembrane protein affects function.

Raman Marker Bands for Secondary Structure Changes of Frozen Therapeutic Monoclonal Antibody Formulations During Thawing

In this work we use Raman spectroscopy for protein characterization in the frozen state. We investigate the behavior of frozen therapeutic monoclonal antibody IgG1 formulation upon thawing by Raman spectroscopy. Secondary and tertiary structure of the protein in three different mab formulations in the frozen state are followed through observation of marker bands for α-helix, β-sheet and random coil. We identify the tyrosine intensity ratio I₈₅₆/I₈₃₀ as a marker for mab aggregation. Upon fast cooling (40 °C/min) to -80 °C we observe a significant increase of random coil and α -helical structures, while this is not the case for slower cooling (20 °C/min) to -80 °C. Most changes in the protein's secondary structure are observed in the course of thawing in the range up to -20 °C, when passing through the glass transitions and cold-crystallization of the two types of freeze-concentrated solutions formed through macro- and microcryoconcentration. An increase of protein concentration and the addition of mannitol suppress secondary structural changes but do no impact on aggregation.

Effect of water activity on the mechanical glass transition and dynamical transition of bacteria

The purpose of this study was to clarify the glass-transition behavior of bacteria (Cronobacter sakazakii) as a function of water activity (aw). From the water sorption isotherm (298 K) for C. sakazakii, monolayer water content and monolayer aw were determined to be 0.0724 g/g-dry matter and 0.252, respectively. Mechanical relaxation was investigated at 298 K. In a higher aw range of over 0.529, the degree of mechanical relaxation increased with an increase in aw. From the effect of aw on the degree of mechanical relaxation, the mechanical awc (aw at which mechanical glass transition occurs at 298 K) was determined to be 0.667. Mean-square displacement of atoms in the bacteria was investigated by incoherent elastic neutron scattering. The mean-square displacement increased gradually with an increase in temperature depending on the aw of samples. From the linear fitting, two or three dynamical transition temperatures (low, middle, and high Tds) were determined at each aw. The low-Td values (142-158 K) were almost independent from aw. There was a minor effect of aw on the middle Td (214-234 K) except for the anhydrous sample (261 K). The high Td (252-322 K) largely increased with the decrease in aw. From the aw dependence of the high Td, the dynamical awc was determined to be 0.675, which was almost equivalent to the mechanical awc. The high Td was assumed to be the glass-transition temperature (Tg), and anhydrous Tg was estimated to be 409 K. In addition, molecular relaxation time (τ) of the bacteria was calculated as a function of aw. From the result, it is suggested that the progress of metabolism in the bacterial system requires a lower τ than approximately 6 × 10-5 s.

Dynamics of Apolipoprotein B-100 in Interaction with Detergent Probed by Incoherent Neutron Scattering

Apolipoprotein B-100 (apo B-100) is the protein moiety of both low- and very-low-density lipoproteins, whose role is crucial to cholesterol and triglyceride transport. Aiming at the molecular dynamics' details of apo B-100, scarcely studied, we performed elastic and quasi-elastic incoherent neutron scattering (EINS, QENS) experiments combining different instruments and time scales. Similar to classical membrane proteins, the solubilization results in remaining detergent, here Nonidet P-40 (NP40). Therefore, we propose a framework for QENS studies of protein-detergent complexes, with the introduction of a combined model, including the experimental apo B-100/NP40 ratio. Relying on the simultaneous analysis of all QENS amplitudes, this approach is sensitive enough to separate both contributions. Its application identified two points: (i) apo B-100 slow dynamics and (ii) the acceleration of NP40 dynamics in the presence of apo B-100. Direct translation of the exposed methodology now makes the investigation of more membrane proteins by neutron spectroscopy achievable.

Dynamic Coupling of Tyrosine 185 with the Bacteriorhodopsin Photocycle, as Revealed by Chemical Shifts, Assisted AF-QM/MM Calculations and Molecular Dynamic Simulations

Aromatic residues are highly conserved in microbial photoreceptors and play crucial roles in the dynamic regulation of receptor functions. However, little is known about the dynamic mechanism of the functional role of those highly conserved aromatic residues during the receptor photocycle. Tyrosine 185 (Y185) is one of the highly conserved aromatic residues within the retinal binding pocket of bacteriorhodopsin (bR). In this study, we explored the molecular mechanism of its dynamic coupling with the bR photocycle by automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) calculations and molecular dynamic (MD) simulations based on chemical shifts obtained by 2D solid-state NMR correlation experiments. We observed that Y185 plays a significant role in regulating the retinal cis–trans thermal equilibrium, stabilizing the pentagonal H-bond network, participating in the orientation switch of Schiff Base (SB) nitrogen, and opening the F42 gate by interacting with the retinal and several key residues along the proton translocation channel. Our findings provide a detailed molecular mechanism of the dynamic couplings of Y185 and the bR photocycle from a structural perspective. The method used in this paper may be applied to the study of other microbial photoreceptors.

Probing Intermolecular Interactions in Phospholipid Bilayers by Far-Infrared Spectroscopy

Fast thermal fluctuations and low frequency phonon modes are thought to play a part in the dynamic mechanisms of many important biological functions in cell membranes. Here we report a detailed far-infrared study of the molecular subpicosecond motions of phospholipid bilayers at various hydrations. We show that these systems sustain several low frequency collective modes and deduce that they arise from vibrations of different lipids interacting through intermolecular van der Waals forces. Furthermore, we observe that the low frequency vibrations of lipid membrane have strong similarities with the subpicosecond motions of liquid water and suggest that resonance mechanisms are an important element to the dynamics coupling between membranes and their hydration water.

Picosecond to nanosecond dynamics provide a source of conformational entropy for protein folding

Myoglobin can be trapped in fully folded structures, partially folded molten globules, and unfolded states under stable equilibrium conditions.

Thermal fluctuations in chemically cross-linked polymers of cyclodextrins

The extent and nature of thermal fluctuations in the innovative class of cross-linked polymers called cyclodextrin nanosponges (CDNS) are investigated, on the picosecond time scale, through elastic and quasielastic neutron scattering experiments. Nanosponges are complex 3D polymer networks where covalent bonds connecting different cyclodextrin (CD) units and intra- and inter-molecular hydrogen-bond interactions cooperate to define the molecular architecture and fast dynamics of the polymer. The study presented here aims to clarify the nature of the conformational rearrangements activated by increasing temperature in the nanosponge polymer, and the constraints imposed by intra- and inter-molecular hydrogen-bond patterns on the internal dynamics of the macromolecule. The results suggest a picture, in which conformational rearrangements involving the torsion of the OH groups around the C-O bonds dominate the internal dynamics of the polymer over the picosecond time scale. Moreover, the estimated values of mean square displacements reveal that the motions of the hydrogen atoms in the nanosponges are progressively hampered as the cross-linking degree of the polymer is increased. Finally, the study of the molecular relaxations suggests a dynamical rearrangement of the hydrogen-bond networks, which is characterized by a jump diffusion motion of the more mobile hydrogen atoms belonging to the OH groups of the CD units. All these findings add further contribution to the rational comprehensive view of the dynamics of these macromolecules, which may be particularly beneficial in designing new drug-delivery systems with tuneable inclusion/release properties.

Does a dry protein undergo a glass transition?

Bovine serum albumin (BSA) with extremely low hydration level 0.04, which is usually defined as dry, has been investigated in the temperature range between 200 and 340 K by incoherent inelastic neutron scattering using the neutron time-of-flight spectrometer FOCUS (PSI, Switzerland). Anomalous temperature behavior has been revealed for relaxational and low-frequency vibrational dynamics of BSA in the vicinity of 250 K. The mean-square atomic displacement has been shown to exhibit a change in the slope of temperature dependence near the same temperature. The presented results point out that the glass-like transition occurs in the dry protein.

Scaling analysis of bio-molecular dynamics derived from elastic incoherent neutron scattering experiments

Numerous neutron scattering studies of bio-molecular dynamics employ a qualitative analysis of elastic scattering data and atomic mean square displacements. We provide a new quantitative approach showing that the intensity at zero energy exchange can be a rich source of information of bio-structural fluctuations on a pico- to nano-second time scale. Elastic intensity scans performed either as a function of the temperature (back-scattering) and∕or by varying the instrumental resolution (time of flight spectroscopy) yield the activation parameters of molecular motions and the approximate structural correlation function in the time domain. The two methods are unified by a scaling function, which depends on the ratio of correlation time and instrumental resolution time. The elastic scattering concept is illustrated with a dynamic characterization of alanine-dipeptide, protein hydration water, and water-coupled protein motions of lysozyme, per-deuterated c-phycocyanin (CPC) and hydrated myoglobin. The complete elastic scattering function versus temperature, momentum exchange, and instrumental resolution is analyzed instead of focusing on a single cross-over temperature of mean square displacements at the apparent onset temperature of an-harmonic motions. Our method predicts the protein dynamical transition (PDT) at Td from the collective (α) structural relaxation rates of the solvation shell as input. By contrast, the secondary (β) relaxation enhances the amplitude of fast local motions in the vicinity of the glass temperature Tg. The PDT is specified by step function in the elastic intensity leading from elastic to viscoelastic dynamic behavior at a transition temperature Td.

Derivation of mean-square displacements for protein dynamics from elastic incoherent neutron scattering

The derivation of mean-square displacements from elastic incoherent neutron scattering (EINS) of proteins is examined, with the aid of experiments on camphor-bound cytochrome P450cam and complementary molecular dynamics simulations. It is shown that a q(4) correction to the elastic incoherent structure factor (where q is the scattering vector) can be simply used to reliably estimate from the experiment both the average mean-square atomic displacement, <Δr(2)> of the nonexchanged hydrogen atoms in the protein and its variance, σ(2). The molecular dynamics simulation results are in broad agreement with the experimentally derived <Δr(2)> and σ(2) derived from EINS on instruments at two different energy resolutions, corresponding to dynamics on the ∼100 ps and ∼1 ns time scales. Significant dynamical heterogeneity is found to arise from methyl-group rotations. The easy-to-apply q(4) correction extends the information extracted from elastic incoherent neutron scattering experiments and should be of wide applicability.

Temperature- and hydration-dependent internal dynamics of stripped human erythrocyte vesicles studied by incoherent neutron scattering

BACKGROUND

We focus on temperature- and hydration-dependence of internal molecular motions in stripped human red blood cell (RBC) vesicles, widely used as a model system for more complex biomembranes.

METHODS

We singled out picosecond local motions of the non-exchangeable hydrogen atoms of RBC vesicles by performing elastic and quasielastic incoherent neutron scattering measurements in dry and heavy water (D₂O)-hydrated RBC powders.

RESULTS

In dry stripped RBCs, hydrogen motions remained harmonic all along the measured temperature range (100-310K) and mean-square displacements (MSDs) exhibited no temperature transition up to 310K. In contrast, MSDs of hydrated stripped RBCs (h ≈ 0.38g D₂O/g dry powder) exhibited a pronounced transition near 260K, with the sharp rise of anharmonic diffusive motions of hydrogen atoms. This transition at ~260K was correlated with both the onset of nonvibrational (harmonic and nonharmonic) motions and the melting of crystallized hydration water.

GENERAL SIGNIFICANCE

In conclusion, we have shown that MSDs in human RBC vesicles are temperature-and hydration-dependent. These results provide insight into biomembrane internal dynamics at picosecond timescale and nanometer length scale. Such motions have been shown to act as the "lubricant" of larger conformational changes on a slower, millisecond timescale that are necessary for important biological processes.

Macromolecular dynamics in red blood cells investigated using neutron spectroscopy

We present neutron scattering measurements on the dynamics of haemoglobin (Hb) in human red blood cells (RBCs) in vivo. Global and internal Hb dynamics were measured in the ps to ns time and Å length scales using quasi-elastic neutron backscattering spectroscopy. We observed the cross over from global Hb short-time to long-time self-diffusion. Both short- and long-time diffusion coefficients agree quantitatively with predicted values from the hydrodynamic theory of non-charged hard-sphere suspensions when a bound water fraction of around 0.23 gram H(2)O per gram Hb is taken into account. The higher amount of water in the cells facilitates internal protein fluctuations in the ps time scale when compared with fully hydrated Hb powder. Slower internal dynamics of Hb in RBCs in the ns time range were found to be rather similar to results obtained with fully hydrated protein powders, solutions and Escherichia coli cells.

Solvent-induced backbone fluctuations and the collective librational dynamics of lysozyme studied by terahertz spectroscopy

THz spectroscopy is used to investigate the dynamics of the globular protein hen egg white lysozyme under varying hydration and temperature conditions. An analysis of the experimental spectra has revealed that the amount of solvent in the hydration shell has a strong influence on the low-frequency protein conformational dynamics and also the arrangement of hydrogen bonds in the protein secondary structure. Furthermore at a hydration level >0.2 we identify collective backbone fluctuations in the protein secondary structure that are not present at low hydration. It is possible that these solvent induced modes are important for the biological function of the protein.

Atomic-scale dynamics inside living cells explored by neutron scattering

Single-particle neutron spectroscopy has contributed important experimental data on molecular dynamics in biological systems. The technique provides information on atomic and molecular motions in macromolecules on the picosecond to the nanosecond time scale, which are essential to biological function. Here, we report on recent neutron measurements performed directly in living cells by using isotope labelling to explore the dynamics of specific cellular components. The paper proposes an integrated view of results on atomic-scale cell water dynamics, internal and global macromolecular motions and solvent isotope effect on macromolecular dynamics. The work established the specific usefulness of the neutron scattering technique to get insight into biologically relevant dynamical features, in particular through comparative measurements. The method developed can now be applied to look for dynamical signatures related to cell characteristics in many different cell types and organelles.

From protons to OXPHOS supercomplexes and Alzheimer's disease: structure-dynamics-function relationships of energy-transducing membranes

By the elucidation of high-resolution structures the view of the bioenergetic processes has become more precise. But in the face of these fundamental advances, many problems are still unresolved. We have examined a variety of aspects of energy-transducing membranes from large protein complexes down to the level of protons and functional relevant picosecond protein dynamics. Based on the central role of the ATP synthase for supplying the biological fuel ATP, one main emphasis was put on this protein complex from both chloroplast and mitochondria. In particular the stoichiometry of protons required for the synthesis of one ATP molecule and the supramolecular organisation of ATP synthases were examined. Since formation of supercomplexes also concerns other complexes of the respiratory chain, our work was directed to unravel this kind of organisation, e.g. of the OXPHOS supercomplex I(1)III(2)IV(1), in terms of structure and function. Not only the large protein complexes or supercomplexes work as key players for biological energy conversion, but also small components as quinones which facilitate the transfer of electrons and protons. Therefore, their location in the membrane profile was determined by neutron diffraction. Physico-chemical features of the path of protons from the generators of the electrochemical gradient to the ATP synthase, as well as of their interaction with the membrane surface, could be elucidated by time-resolved absorption spectroscopy in combination with optical pH indicators. Diseases such as Alzheimer's dementia (AD) are triggered by perturbation of membranes and bioenergetics as demonstrated by our neutron scattering studies.

Infrared monitoring of interlayer water in stacks of purple membranes

The thermodynamic behavior of films of hydrated purple membranes from Halobacterium salinarum and the water confined in it was studied by Fourier transform infrared spectroscopy in the 180-280 K range. Unlike bulk water, water in the thin layers sandwiched between the biological membranes does not freeze at 273 K but will be supercooled to approximately 256 K. The melting point is unaffected, leading to hysteresis between 250 and 273 K. In its heating branch, a gradually increasing light-scattering by ice is observed with rate-limiting kinetics of tens of minutes. Infrared (IR) spectra decomposition provided extinction coefficients for the confined water vibrational bands and their changes upon freezing. Because of the hysteresis, at any given temperature in the 255-270 K range, the interbilayer water could be either liquid or frozen, depending on thermal history. We find that this difference affects the dynamics of the bacteriorhodopsin photocycle in the hysteresis range: the decay of the M and N states and the redistribution between them are different depending on whether or not the water was initially precooled to below the freezing point. However, freezing of interbilayer water does block the M to N transition. Unlike the water, the purple membrane lipids do not undergo any IR-detectable phase transition in the 180-280 K range.

Dynamical heterogeneity of specific amino acids in bacteriorhodopsin

Components of biological macromolecules, complexes and membranes are animated by motions occurring over a wide range of time and length scales, the synergy of which is at the basis of biological activity. Understanding biological function thus requires a detailed analysis of the underlying dynamical heterogeneity. Neutron scattering, using specific isotope labeling, and molecular dynamics simulations were combined in order to study the dynamics of specific amino acid types in bacteriorhodopsin within the purple membrane (PM) of Halobacterium salinarum. Motions of leucine, isoleucine and tyrosine residues on the pico- to nanosecond time scale were examined separately as a function of temperature from 20 to 300 K. The dynamics of the three residue types displayed different temperature dependence: isoleucine residues have larger displacements compared to the global PM above 120 K; leucine residues have displacements similar to that of PM in the entire temperature range studied; and tyrosine residues have displacements smaller than that of the average membrane in an intermediate temperature range. Experimental features were mostly well reproduced by molecular dynamics simulations performed at five temperatures, which allowed the dynamical characterisation of the amino acids under study as a function of local environment. The resulting dynamical map of bacteriorhodopsin revealed that movements of a specific residue are determined by both its environment and its residue type.

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 approximately 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.

Hydration dependence of active core fluctuations in bacteriorhodopsin

We used neutron scattering and specific hydrogen-deuterium labeling to investigate the thermal dynamics of isotope-labeled amino acids and retinal, predominantly in the active core and extracellular moiety of bacteriorhodopsin (BR) in the purple membrane and the dynamical response to hydration. Measurements on two neutron spectrometers allowed two populations of motions to be characterized. The lower amplitude motions were found to be the same for both the labeled amino acids and retinal of BR and the global membrane. The larger amplitude dynamics of the labeled part, however, were found to be more resilient than the average membrane, suggesting their functional importance. The response to hydration was characterized, showing that the labeled part of BR is not shielded from hydration effects. The results suggest that the inhibition of high-amplitude motions by lowering hydration may play a key role in the slowing down of the photocycle and the proton pumping activity of BR.

Dynamics of hydration water in deuterated purple membranes explored by neutron scattering

The function and dynamics of proteins depend on their direct environment, and much evidence has pointed to a strong coupling between water and protein motions. Recently however, neutron scattering measurements on deuterated and natural-abundance purple membrane (PM), hydrated in H2O and D2O, respectively, revealed that membrane and water motions on the ns-ps time scale are not directly coupled below 260 K (Wood et al. in Proc Natl Acad Sci USA 104:18049-18054, 2007). In the initial study, samples with a high level of hydration were measured. Here, we have measured the dynamics of PM and water separately, at a low-hydration level corresponding to the first layer of hydration water only. As in the case of the higher hydration samples previously studied, the dynamics of PM and water display different temperature dependencies, with a transition in the hydration water at 200 K not triggering a transition in the membrane at the same temperature. Furthermore, neutron diffraction experiments were carried out to monitor the lamellar spacing of a flash-cooled deuterated PM stack hydrated in H2O as a function of temperature. At 200 K, a sudden decrease in lamellar spacing indicated the onset of long-range translational water diffusion in the second hydration layer as has already been observed on flash-cooled natural-abundance PM stacks hydrated in D2O (Weik et al. in J Mol Biol 275:632-634, 2005), excluding thus a notable isotope effect. Our results reinforce the notion that membrane-protein dynamics may be less strongly coupled to hydration water motions than the dynamics of soluble proteins.

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.

Influence of hydration on the dynamics of lysozyme

Quasielastic neutron and light-scattering techniques along with molecular dynamics simulations were employed to study the influence of hydration on the internal dynamics of lysozyme. We identified three major relaxation processes that contribute to the observed dynamics in the picosecond to nanosecond time range: 1), fluctuations of methyl groups; 2), fast picosecond relaxation; and 3), a slow relaxation process. A low-temperature onset of anharmonicity at T approximately 100 K is ascribed to methyl-group dynamics that is not sensitive to hydration level. The increase of hydration level seems to first increase the fast relaxation process and then activate the slow relaxation process at h approximately 0.2. The quasielastic scattering intensity associated with the slow process increases sharply with an increase of hydration to above h approximately 0.2. Activation of the slow process is responsible for the dynamical transition at T approximately 200 K. The dependence of the slow process on hydration correlates with the hydration dependence of the enzymatic activity of lysozyme, whereas the dependence of the fast process seems to correlate with the hydration dependence of hydrogen exchange of lysozyme.

Thermodynamic and dynamic factors involved in the stability of native protein structure in amorphous solids in relation to levels of hydration

The internal, dynamical fluctuations of protein molecules exhibit many of the features typical of polymeric and bulk small molecule glass forming systems. The response of a protein's internal molecular mobility to temperature changes is similar to that of other amorphous systems, in that different types of motions freeze out at different temperatures, suggesting they exhibit the alpha-beta-modes of motion typical of polymeric glass formers. These modes of motion are attributed to the dynamic regimes that afford proteins the flexibility for function but that also develop into the large-scale collective motions that lead to unfolding. The protein dynamical transition, T(d), which has the same meaning as the T(g) value of other amorphous systems, is attributed to the temperature where protein activity is lost and the unfolding process is inhibited. This review describes how modulation of T(d) by hydration and lyoprotectants can determine the stability of protein molecules that have been processed as bulk, amorphous materials. It also examines the thermodynamic, dynamic, and molecular factors involved in stabilizing folded proteins, and the effects typical pharmaceutical processes can have on native protein structure in going from the solution state to the solid state.

Protein–water displacement distributions

The statistical properties of fast protein-water motions are analyzed by dynamic neutron scattering experiments. Using isotopic exchange, one probes either protein or water hydrogen displacements. A moment analysis of the scattering function in the time domain yields model-independent information such as time-resolved mean square displacements and the Gauss-deviation. From the moments, one can reconstruct the displacement distribution. Hydration water displays two dynamical components, related to librational motions and anomalous diffusion along the protein surface. Rotational transitions of side chains, in particular of methyl groups, persist in the dehydrated and in the solvent-vitrified protein structure. The interaction with water induces further continuous protein motions on a small scale. Water acts as a plasticizer of displacements, which couple to functional processes such as open-closed transitions and ligand exchange.

Methyl group dynamics as a probe of the protein dynamical transition

Hydrated proteins undergo a dynamical transition around 200 K from glasslike to liquidlike motion. Molecular dynamics simulations have been used to study the temperature dependence of the dynamics of ribonuclease A in the hydrated crystal, a model dehydrated powder, and aqueous solution. Changes in the dynamics accompanying the transition throughout the protein have been quantified in terms of the mean-squared fluctuations (MSFs) of methyl hydrogen atoms on the 100 ps time scale. In solution at 300 K the MSFs span a broad distribution, consistent with NMR relaxation measurements. The MSF distribution in the hydrated crystal at 300 K is qualitatively similar to the solution result, except for a slight shift to lower values, and dehydration results in a dramatic shift of the MSFs to lower values. As the temperature is lowered, the whole distribution of methyl group fluctuations in the hydrated crystal shifts to lower values. Most of the methyl groups in the hydrated protein display a nonlinear temperature dependence with a dynamical transition at approximately 200 K, but most methyl groups do not undergo a transition in the dehydrated protein. We conclude that the dynamical transition occurs throughout most of the protein and that solvent is required for the transition.

Structure, dynamics and reactions of protein hydration water

The apparent simplicity of the water molecule belies the wide range of fascinating protein phenomena in which it participates. We review recent computer simulation work on buried, internal water molecules, discussing the thermodynamics of water molecule binding and the participation of water in proton transfer reactions. Surface water molecules are also considered, with emphasis on the modification of average solvent structure on a protein surface, the role of water in the protein dynamical 'glass' transition and a simplified description of the protein motions thereby activated.

Neutron frequency windows and the protein dynamical transition

Proteins undergo an apparent dynamical transition on temperature variation that has been correlated with the onset of function. The transition in the mean-square displacement, , that is observed using a spectrometer or computer simulation, depends on the relationship between the timescales of the relaxation processes activated and the timescale accessible to the instrument or simulation. Models are described of two extreme situations---an "equilibrium" model, in which the long-time dynamics changes with temperature and all motions are resolved by the instrument used; and a "frequency window" model, in which there is no change in the long-time dynamics but as the temperature increases, the relaxation frequencies move into the instrumental range. Here we demonstrate that the latter, frequency-window model can describe the temperature and timescale dependences of both the intermediate neutron scattering function and derived from molecular dynamics simulations of a small protein in a cryosolution. The frequency-window model also describes the energy-resolution and temperature-dependences of obtained from experimental neutron scattering on glutamate dehydrogenase in the same solvent. Although equilibrium effects should also contribute to dynamical transitions in proteins, the present results suggests that frequency-window effects can play a role in the simulations and experiments examined. Finally, misquotations of previous findings are discussed in the context of solvent activation of protein dynamics and the possible relationship of this to activity.

Picosecond internal dynamics of lysozyme as affected by thermal unfolding in nonaqueous environment

A neutron-scattering investigation of the internal picosecond dynamics of lysozyme solvated in glycerol as a function of temperature in the range 200-410 K has been undertaken. The inelastic contribution to the measured intensity is characterized by the presence of a bump generally known as "boson peak", clearly distinguishable at low temperature. When the temperature is increased the quasielastic component of the spectrum becomes more and more intrusive and progressively overwhelms the vibrational bump. This happens especially for T > 345 K when the protein goes through an unfolding process, which leads to the complete denaturation. The quasielastic term is the superposition of two components whose intensities and linewidths have been studied as a function of temperature. The slower component describes motions with characteristic times of approximately 4 ps corresponding to reorientations of polypeptide side chains. Both the intensity and linewidth of this kind of relaxations show two distinct regimes with a crossover in the temperature range where the melting process occurs, thus suggesting the presence of a dynamical transition correlated to the protein unfolding. Conversely the faster component might be ascribed to the local dynamics of hydrogen atoms caged by the nearest neighbors with characteristic time of approximately 0.3 ps.

The role of dynamics in enzyme activity

Although protein function is thought to depend on flexibility, precisely how the dynamics of the molecule and its environment contribute to catalytic mechanisms is unclear. We review experimental and computational work relating to enzyme dynamics and function, including the role of solvent. The evidence suggests that fast motions on the 100 ps timescale, and any motions coupled to these, are not required for enzyme function. Proteins where the function is electron transfer, proton tunneling, or ligand binding may have different dynamical dependencies from those for enzymes, and enzymes with large turnover numbers may have different dynamical dependencies from those that turn over more slowly. The timescale differences between the fastest anharmonic fluctuations and the barrier-crossing rate point to the need to develop methods to resolve the range of motions present in enzymes on different time- and lengthscales.

Translational hydration water dynamics drives the protein glass transition

Experimental and computer simulation studies have revealed the presence of a glass-like transition in the internal dynamics of hydrated proteins at approximately 200 K involving an increase of the amplitude of anharmonic dynamics. This increase in flexibility has been correlated with the onset of protein activity. Here, we determine the driving force behind the protein transition by performing molecular dynamics simulations of myoglobin surrounded by a shell of water. A dual heat bath method is used with which, in any given simulation, the protein and solvent are held at different temperatures, and sets of simulations are performed varying the temperature of the components. The results show that the protein transition is driven by a dynamical transition in the hydration water that induces increased fluctuations primarily in side chains in the external regions of the protein. The water transition involves activation of translational diffusion and occurs even in simulations where the protein atoms are held fixed.

Molecular dynamics decomposition of temperature-dependent elastic neutron scattering by a protein solution

Molecular dynamics simulations are performed of bovine pancreatic trypsin inhibitor in a cryosolution over a range of temperatures from 80 to 300 K and the origins identified of elastic dynamic neutron scattering from the solution. The elastic scattering and mean-square displacement calculated from the molecular dynamics trajectories are in reasonable agreement with experiments on a larger protein in the same solvent. The solvent and protein contributions to the scattering from the simulation model are determined. At lower temperatures (< approximately 200 K) or on shorter timescales ( approximately 10 ps) the scattering contributions are proportional to the isotopic nuclear scattering cross-sections of each component. However, for T > 200 K marked deviations from these cross-sections are seen due to differences in the dynamics of the components of the solution. Rapid activation of solvent diffusion leads to the variation with temperature of the total elastic intensity being determined largely by that of the solvent. At higher temperatures (>240 K) and longer times ( approximately 100 ps) the protein makes the only significant contribution to the scattering, the solvent scattering having moved out of the accessible time-space window. Decomposition of the protein mean-square displacement shows that the observed dynamical transition in the solution at 200-220 K involves activation of both internal motions and external whole-molecule rotational and translational diffusion. The proportion that the external dynamics contributes to the protein mean-square displacement increases to approximately 30 and 60% at 300 K on the 10- and 100-ps timescales, respectively.

Energy resolution and dynamical heterogeneity effects on elastic incoherent neutron scattering from molecular systems

Incoherent neutron scattering is widely used to probe picosecond-nanosecond time scale dynamics of molecular systems. In systems of spatially confined atoms the relatively high intensity of elastic incoherent neutron scattering is often used to obtain a first estimate of the dynamics present. For many complex systems, however, experimental elastic scattering is difficult to interpret unambiguously using analytical dynamical models that go beyond the determination of an average mean-square displacement. To circumvent this problem a description of the scattering is derived here that encompasses a variety of analytical models in a common framework. The framework describes the time-converged part of the dynamic structure factor [the elastic incoherent scattering function (EISF)] and lends itself to practical use by explicitly incorporating effects due to the finite energy resolution of the instrument used. The dependence of the elastic scattering on wave vector is examined, and it is shown how heterogeneity in the distribution of mean-square displacements can be related to deviations of the scattering from Gaussian behavior. In this case, a correction to fourth order in the scattering vector can be used to extract the variance of the distribution of mean-square displacements. The formalism is used in a discussion of measurements on dynamics accompanying the glass transition in molecular systems. By fitting to experimental data obtained on a protein solution the present methodology is used to show how the existence of a temperature-dependent relaxation frequency can lead to a transition in the measured mean-square displacement in the absence of an EISF change.