Heme-solvent coupling: a Mössbauer study of myoglobin in sucrose

Reactivity Tracking of an Enzyme Progress Coordinate

The reactivity of individual solvent-coupled protein configurations is used to track and resolve the progress coordinate for the core reaction sequence of substrate radical rearrangement and hydrogen atom transfer in the ethanolamine ammonia-lyase (EAL) enzyme from Salmonella enterica. The first-order decay of the substrate radical intermediate is the monitored reaction. Heterogeneous confinement from sucrose hydrates in the mesophase solvent surrounding the cryotrapped protein introduces distributed kinetics in the non-native decay of the substrate radical pair capture substate, which arise from an ensemble of configurational microstates. Reaction rates increase by >103-fold across the distribution to approach that for the native enabled substate for radical rearrangement, which reacts with monotonic kinetics. The native progress coordinate thus involves a collapse of the configuration space to generate optimized reactivity. Reactivity tracking reveals fundamental features of solvent-protein-reaction configurational coupling and leads to a model that refines the ensemble paradigm of enzyme catalysis for strongly adiabatic chemical steps.

Temperature Dependence of Structural Dynamics at the Catalytic Cofactor of [FeFe]-hydrogenase

[FeFe]-hydrogenases are nature's blueprint for efficient hydrogen turnover. Understanding their enzymatic mechanism may improve technological H₂ fuel generation. The active-site cofactor (H-cluster) consists of a [4Fe-4S] cluster ([4Fe]_H), cysteine-linked to a diiron site ([2Fe]_H) carrying an azadithiolate (adt) group, terminal cyanide and carbon monoxide ligands, and a bridging carbon monoxide (μCO) in the oxidized protein (Hox). Recently, the debate on the structure of reduced H-cluster states was intensified by the assignment of new species under cryogenic conditions. We investigated temperature effects (4-280 K) in infrared (IR) and X-ray absorption spectroscopy (XAS) data of [FeFe]-hydrogenases using fit analyses and quantum-chemical calculations. IR data from our laboratory and literature sources were evaluated. At ambient temperatures, reduced H-cluster states with a bridging hydride (μH^-, in Hred and Hsred) or with an additional proton at [4Fe]_H (Hred') or at the distal iron of [2Fe]_H (Hhyd) prevail. At cryogenic temperatures, these species are largely replaced by states that hold a μCO, lack [4Fe]_H protonation, and bind an additional proton at the adt nitrogen (HredH⁺ and HsredH⁺). XAS revealed the atomic coordinate dispersion (i.e., the Debye-Waller parameter, 2σ²) of the iron-ligand bonds and Fe-Fe distances in the oxidized and reduced H-cluster. 2σ² showed a temperature dependence typical for the so-called protein-glass transition, with small changes below ∼200 K and a pronounced increase above this "breakpoint". This behavior is attributed to the freezing-out of larger-scale anharmonic motions of amino acid side chains and water species. We propose that protonation at [4Fe]_H as well as ligand rearrangement and μH^- binding at [2Fe]_H are impaired because of restricted molecular mobility at cryogenic temperatures so that protonation can be biased toward adt. We conclude that a H-cluster with a μCO, selective [4Fe]_H or [2Fe]_H protonation, and catalytic proton transfer via adt facilitates efficient H₂ conversion in [FeFe]-hydrogenase.

Temperature Induced Dynamical Transition of Biomolecules in Polarizable and Nonpolarizable TIP3P Water

Temperature induced dynamical transition (DT), associated with a sharp rise in molecular flexibility, is well-known to be exhibited between 270 and 280 K in glycerol to 200-230 K in hydrated biomolecules and is controlled by diffusivity (viscosity) of the solvation layer. In the molecular dynamics (MD) community, especially for water as a solvent, this has been an intense area of research despite decades of investigations. However, in general, water in these studies is described by empirical nonpolarizable force fields in which electronic polarizability is treated implicitly with effective charges and related parameters. This might have led to the present trait of discovery that DTs of biomolecules, irrespective of the potential functions for water models used, occur within a narrow band of temperature variation (30-40 K). Whereas a water molecule in a biomolecular surface and one in bulk are polarized differently, therefore explicit treatment of water polarizability would be a powerful approach toward the treatment of hydration water, believed to cause the DT manifestation. Using MD simulations, we investigated the effects of polarizable water on the DT of biomolecules and the dynamic properties of hydration water. We chose two types of solutes: globular protein (lysozyme) and more open and flexible RNAs (a hairpin and a riboswitch) with different natures of hydrophilic sites than proteins in general. We found that the characteristic temperature of DT ( T_DT) for the solutes in polarizable water is always higher than that in its nonpolarizable counterpart. In particular, for RNAs, the variations are found to be ∼45 K between the two water models, whereas for the more compact lysozyme, it is only ∼4 K. The results are discussed in light of the enormous increase in relaxation times of a liquid upon cooling in the paradigm of dynamic switchover in hydration water with liquid-liquid phase transition, derived from the existence of the second critical point. Our result supports the idea that structures of biomolecules and their interactions with the hydration water determines T_DT and provides evidence for the decisive role of polarizable water on the onset of DT, which has been hitherto ignored.

Role of Hydration Layer in Dynamical Transition in Proteins: Insights from Translational Self-Diffusivity

Elucidation of the role of hydration water underpinning dynamical crossover in proteins has proven challenging. Indeed, many contradictory findings in the literature seek to establish either causal or correlative links between water and protein behavior. Here, via molecular dynamics, we compute the temperature dependence of mean-square displacement and translational self-diffusivities for both hen egg white lysozyme and its hydration layer from 190 to 300 K. We find that the protein's mobility increases sharply at ∼230 K, indicating dynamical onset; concerted motion with hydration-water molecules is evident up to ∼285 K, confirming dynamical correlation between them. Exploring underlying mechanisms of such concerted motion, we scrutinize the water-protein hydrogen-bonding network as a function of temperature, noting sharp deviation from linearity of the hydrogen bond number's profile with temperature originating near the protein dynamical transition. Our studies reveal a common temperature profile/dependence of self-diffusivity values of the protein, hydration water, and the bulk solvent, originating from a common dependence on the bulk solvent viscosity, η_S. The key mechanistic role adopted by the protein-water hydrogen bond network in relation to the onset of proteins' dynamical transition is also discussed.

Dynamics of hydrated proteins and bio-protectants: Caged dynamics, β-relaxation, and α-relaxation

BACKGROUND

The properties of the three dynamic processes, α-relaxation, ν-relaxation, and caged dynamics in aqueous mixtures and hydrated proteins are analogous to corresponding processes found in van der Waals and polymeric glass-formers apart from minor differences.

METHODS

Collection of various experimental data enables us to characterize the structural α-relaxation of the protein coupled to hydration water (HW), the secondary or ν-relaxation of HW, and the caged HW process.

RESULTS

From the T-dependence of the ν-relaxation time of hydrated myoglobin, lysozyme, and bovine serum albumin, we obtain Ton at which it enters the experimental time windows of Mössbauer and neutron scattering spectroscopies, coinciding with protein dynamical transition (PDT) temperature Td. However, for all systems considered, the α-relaxation time at Ton or Td is many orders of magnitude longer. The other step change of the mean-square-displacement (MSD) at Tg_alpha originates from the coupling of the nearly constant loss (NCL) of caged HW to density. The coupling of the NCL to density is further demonstrated by another step change at the secondary glass temperature Tg_beta in two bio-protectants, trehalose and sucrose.

CONCLUSIONS

The structural α-relaxation plays no role in PDT. Since PDT is simply due to the ν-relaxation of HW, the term PDT is a misnomer. NCL of caged dynamics is coupled to density and show transitions at lower temperature, Tg_beta and Tg_alpha.

GENERAL SIGNIFICANCE

The so-called protein dynamical transition (PDT) of hydrated proteins is not caused by the structural α-relaxation of the protein but by the secondary ν-relaxation of hydration water. "This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

Atomistic details of protein dynamics and the role of hydration water

BACKGROUND

The importance of protein dynamics for their biological activity is now well recognized. Different experimental and computational techniques have been employed to study protein dynamics, hierarchy of different processes and the coupling between protein and hydration water dynamics. Yet, understanding the atomistic details of protein dynamics and the role of hydration water remains rather limited.

SCOOP OF REVIEW

Based on overview of neutron scattering, molecular dynamic simulations, NMR and dielectric spectroscopy results we present a general picture of protein dynamics covering time scales from faster than ps to microseconds and the influence of hydration water on different relaxation processes.

MAJOR CONCLUSIONS

Internal protein dynamics spread over a wide time range from faster than picosecond to longer than microseconds. We suggest that the structural relaxation in hydrated proteins appears on the microsecond time scale, while faster processes present mostly motion of side groups and some domains. Hydration water plays a crucial role in protein dynamics on all time scales. It controls the coupled protein-hydration water relaxation on 10-100ps time scale. This process defines the friction for slower protein dynamics. Analysis suggests that changes in amount of hydration water affect not only general friction, but also influence significantly the protein's energy landscape.

GENERAL SIGNIFICANCE

The proposed atomistic picture of protein dynamics provides deeper understanding of various relaxation processes and their hierarchy, similarity and differences between various biological macromolecules, including proteins, DNA and RNA. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

Protein dynamics: from rattling in a cage to structural relaxation

We present an overview of protein dynamics based mostly on results of neutron scattering, dielectric relaxation spectroscopy and molecular dynamics simulations. We identify several major classes of protein motions on the time scale from faster than picoseconds to several microseconds, and discuss the coupling of these processes to solvent dynamics. Our analysis suggests that the microsecond backbone relaxation process might be the main structural relaxation of the protein that defines its glass transition temperature, while faster processes present some localized secondary relaxations. Based on the overview, we formulate a general picture of protein dynamics and discuss the challenges in this field.

Local minimum in fragility for trehalose/glycerol mixtures: implications for biopharmaceutical stabilization

Approximately a decade ago it was observed that adding a small amount (5 wt %) of glycerol to trehalose could substantially improve the stability of enzymes stored in these glasses even though the final glass transition temperature (Tg) was reduced by ∼20 K. This finding inspired great interest in the fast dynamics of dehydrated trehalose/glycerol mixtures, leading to the observation that suppression of fast dynamics was optimal in the presence of ∼5 wt % of glycerol. It was also recognized that the fast dynamics should, in theory, be related to the fragility of these glass formers, but experimental confirmation of this hypothesis has been lacking for trehalose/glycerol mixtures or any other mixtures of this nature. In the present study a dynamic mechanical analyzer (DMA) was used to determine both the Tg and the kinetic fragility index (m) of trehalose/glycerol mixtures within the mass fraction range of 80-100 wt % of trehalose. It was found that the fragility index correlated with the mass fraction of trehalose in a nonmonotonic manner, with a local minimum between 87.5 and 95 wt % of trehalose, whereas the composition dependence of Tg was found to follow a Gordon-Taylor-like relationship, with no local minimum. The composition of 5-12.5 wt % glycerol in trehalose thus yielded a matrix that maximized the strong glass-forming contribution of glycerol, while minimizing its Tg lowering effect. This quantitative evidence supports speculation about the fragility characteristics of these mixtures that has been ongoing for the past decade. The DMA-based Tg and fragility determination method developed in this study represents a new approach for identifying optimal compositions for preservation of biologics.

One role of hydration water in proteins: key to the "softening" of short time intraprotein collective vibrations of a specific length scale

High resolution inelastic X-ray scattering (IXS) experiments show that the "phonon energy softening" and "phonon population enhancement" observed in a hydrated native protein when increasing the temperature from 200 K to physiological temperature are not directly related to the protein structure. Such phenomena were also observed in a denatured sample without a defined tertiary structure and with a limited residual secondary structure. However, in a dry sample, such "softening" is strongly suppressed. These facts suggest that the above-mentioned protein "softening" phenomenon is water-induced. In addition, increasing the hydration level can also induce "phonon energy softening" at room temperature, but not at 200 K. This change may be due to a qualitative difference in the dynamics of hydration water at 200 K and at room temperature.

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.

Protein surface and core dynamics show concerted hydration-dependent activation

By specifically labeling leucine/valine methyl groups and lysine side chains "inside" and "outside" dynamics of proteins on the nanosecond timescale are compared using neutron scattering. Surprisingly, both groups display similar dynamics as a function of temperature, and the buried hydrophobic core is sensitive to hydration and undergoes a dynamical transition.

Temperature-dependent macromolecular X-ray crystallography

X-ray crystallography provides structural details of biological macromolecules. Whereas routine data are collected close to 100 K in order to mitigate radiation damage, more exotic temperature-controlled experiments in a broader temperature range from 15 K to room temperature can provide both dynamical and structural insights. Here, the dynamical behaviour of crystalline macromolecules and their surrounding solvent as a function of cryo-temperature is reviewed. Experimental strategies of kinetic crystallography are discussed that have allowed the generation and trapping of macromolecular intermediate states by combining reaction initiation in the crystalline state with appropriate temperature profiles. A particular focus is on recruiting X-ray-induced changes for reaction initiation, thus unveiling useful aspects of radiation damage, which otherwise has to be minimized in macromolecular crystallography.

Dynamical transition of protein-hydration water

Thin layers of water on biomolecular and other nanostructured surfaces can be supercooled to temperatures not accessible with bulk water. Chen et al. [Proc. Natl. Acad. Sci. U.S.A. 103, 9012 (2006)]10.1073/pnas.0602474103 suggested that anomalies near 220 K observed by quasielastic neutron scattering can be explained by a hidden critical point of bulk water. Based on more sensitive measurements of water on perdeuterated phycocyanin, using the new neutron backscattering spectrometer SPHERES, and an improved data analysis, we present results that show no sign of such a fragile-to-strong transition. The inflection of the elastic intensity at 220 K has a dynamic origin that is compatible with a calorimetric glass transition at 170 K. The temperature dependence of the relaxation times is highly sensitive to data evaluation; it can be brought into perfect agreement with the results of other techniques, without any anomaly.

Thermal denaturation of myoglobin in water--disaccharide matrixes: relation with the glass transition of the system

Proteins embedded in glassy saccharide systems are protected against adverse environmental conditions [Crowe et al. Annu. Rev. Physiol. 1998, 60, 73-103]. To further characterize this process, we studied the relationship between the glass transition temperature of the protein-containing saccharide system (T(g)) and the temperature of thermal denaturation of the embedded protein (T(den)). To this end, we studied by differential scanning calorimetry the thermal denaturation of ferric myoglobin in water/disaccharide mixtures containing nonreducing (trehalose, sucrose) or reducing (maltose, lactose) disaccharides. All the samples studied are, at room temperature, liquid systems whose viscosity varies from very low to very large values, depending on the water content. At a high water/saccharide mole ratio, homogeneous glass formation does not occur; regions of glass form, whose T(g) does not vary by varying the saccharide content, and the disaccharide barely affects the myoglobin denaturation temperature. At a suitably low water/saccharide mole ratio, by lowering the temperature, the systems undergo transition to the glassy state whose T(g) is determined by the water content; the Gordon-Taylor relationship between T(g) and the water/disaccharide mole ratio is obeyed; and T(den) increases by decreasing the hydration regardless of the disaccharide, such effect being entropy-driven. The presence of the protein was found to lower the T(g). Furthermore, for nonreducing disaccharides, plots of T(den) vs T(g) give linear correlations, whereas for reducing disaccharides, data exhibit an erratic behavior below a critical water/disaccharide ratio. We ascribe this behavior to the likelihood that in the latter samples, proteins have undergone Maillard reaction before thermal denaturation.

Frontiers in biotransport: water transport and hydration

Biotransport, by its nature, is concerned with the motions of molecules in biological systems while water remains as the most important and the most commonly studied molecule across all disciplines. In this review, we focus on biopreservation and thermal therapies from the perspective of water, exploring how its molecular motions, properties, kinetic, and thermodynamic transitions govern biotransport phenomena and enable preservation or controlled destruction of biological systems.

Frontiers in Biotransport: Water Transport and Hydration

Biotransport, by its nature, is concerned with the motions of molecules in biological systems while water remains as the most important and the most commonly studied molecule across all disciplines. In this review, we focus on biopreservation and thermal therapies from the perspective of water, exploring how its molecular motions, properties, kinetic, and thermodynamic transitions govern biotransport phenomena and enable preservation or controlled destruction of biological systems.

Observation of a low-temperature, dynamically driven structural transition in a polypeptide by solid-state NMR spectroscopy

At reduced temperatures, proteins and other biomolecules are generally found to exhibit dynamic as well as structural transitions. This includes a so-called protein glass transition that is universally observed in systems cooled between 200 and 230 K, and which is generally attributed to interactions between hydrating solvent molecules and protein side chains. However, there is also experimental and theoretical evidence for a low-temperature transition in the intrinsic dynamics of the protein itself, absent any solvent. Here, we use low-temperature solid-state NMR to examine site-specific fluctuations in atomic structure and dynamics in the absence of solvents. In particular, we employ magic angle spinning NMR to examine a structural phase transition associated with dynamic processes in a solvent-free polypeptide, N-f-MLF-OH, lattice at temperatures as low as 90 K. This transition is characterized by the appearance of an extra set of lines in 1D (15)N spectra as well as additional cross peaks in 2D (13)C-(13)C and (13)C-(15)N spectra. Interestingly, the gradual, temperature-dependent appearance of the new spectral component is not accompanied by the line broadening typical of dynamic transitions. A direct comparison between the spectra of N-f-MLF-OH and the analog N-f-MLF-OMe, which does not display this transition, indicates a correlation of the structural transition to the temperature dependent motion of the aromatic phenylalanine side chain. Several quantitative solid state NMR experiments were employed to provide site-specific measurements of structural and motional features of the observed transition.

Resilience of the iron environment in heme proteins

Conformational flexibility is essential to the functional behavior of proteins. We use an effective force constant introduced by Zaccai, the resilience, to quantify this flexibility. Site-selective experimental and computational methods allow us to determine the resilience of heme protein active sites. The vibrational density of states of the heme Fe determined using nuclear resonance vibrational spectroscopy provides a direct experimental measure of the resilience of the Fe environment, which we compare quantitatively with values derived from the temperature dependence of atomic mean-squared displacements in molecular dynamics simulations. Vibrational normal modes in the THz frequency range dominate the resilience. Both experimental and computational methods find a higher resilience for cytochrome c than for myoglobin, which we attribute to the increased number of covalent links to the peptide in the former protein. For myoglobin, the resilience of the iron environment is larger than the average resilience previously determined for hydrogen sites using neutron scattering. Experimental results suggest a slightly reduced resilience for cytochrome c upon oxidation, although the change is smaller than reported in previous Mössbauer investigations on a bacterial cytochrome c, and is not reproduced by the simulations. Oxidation state also has no significant influence on the compressibility calculated for cyt c, although a slightly larger compressibility is predicted for myoglobin.

Protein dynamical transition does not require protein structure

Terahertz time domain spectroscopy shows that the protein dynamical transition, the rapid increase in protein dynamics occurring at approximately 200 K, needs neither tertiary nor secondary structure. Further, short chain alanine studies find a dynamical transition down to penta-alanine, with no transition observed for di-alanine or tri-alanine. These results reveal the temperature dependence arises strictly from the side-chain interaction with the solvent. The lack of a transition for shorter chain peptides may indicate a qualitative change in this interaction occurs at a specific peptide chain length.

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.

Coupling of protein and hydration-water dynamics in biological membranes

The dynamical coupling between proteins and their hydration water is important for the understanding of macromolecular function in a cellular context. In the case of membrane proteins, the environment is heterogeneous, composed of lipids and hydration water, and the dynamical coupling might be more complex than in the case of the extensively studied soluble proteins. Here, we examine the dynamical coupling between a biological membrane, the purple membrane (PM), and its hydration water by a combination of elastic incoherent neutron scattering, specific deuteration, and molecular dynamics simulations. Examining completely deuterated PM, hydrated in H(2)O, allowed the direct experimental exploration of water dynamics. The study of natural abundance PM in D(2)O focused on membrane dynamics. The temperature-dependence of atomic mean-square displacements shows inflections at 120 K and 260 K for the membrane and at 200 K and 260 K for the hydration water. Because transition temperatures are different for PM and hydration water, we conclude that ps-ns hydration water dynamics are not directly coupled to membrane motions on the same time scale at temperatures <260 K. Molecular-dynamics simulations of hydrated PM in the temperature range from 100 to 296 K revealed an onset of hydration-water translational diffusion at approximately 200 K, but no transition in the PM at the same temperature. Our results suggest that, in contrast to soluble proteins, the dynamics of the membrane protein is not controlled by that of hydration water at temperatures <260 K. Lipid dynamics may have a stronger impact on membrane protein dynamics than hydration water.

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

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

Coupling between lysozyme and trehalose dynamics: microscopic insights from molecular-dynamics simulations

We have carried out molecular-dynamics simulations on fully flexible all-atom models of the protein lysozyme immersed in trehalose, an effective biopreservative, with the purpose of exploring the nature and extent of the dynamical coupling between them. Our study shows a strong coupling over a wide range of temperatures. We found that the onset of anharmonic behavior was dictated by changes in the dynamics and relaxation processes in the trehalose glass. The physical origin of protein-trehalose coupling was traced to the hydrogen bonds formed at the interface between the protein and the solvent. Moreover, protein-solvent hydrogen bonding was found to control the structural relaxation of the protein. The dynamics of the protein was found to be heterogeneous; the motions of surface and core atoms had different dependencies on temperature and, in addition, the surface atoms were more sensitive to the dynamics of the solvent than the core atoms. From the solvent perspective we found that the dynamics near the protein surface showed an unexpected enhanced mobility compared to the bulk. These results shed some light on the microscopic origins of the dynamical coupling in protein-solvent systems.

A comparative study of carboxy myoglobin in saccharide-water systems by molecular dynamics simulation

Results from room-temperature molecular dynamics simulation on a system containing carboxy-myoglobin, water, and maltose molecules are reported. Protein atomic fluctuations, protein-solvent and solvent-solvent hydrogen bonding have been analyzed and compared to the ones in trehalose-water and sucrose-water systems (Proteins 2005, 59, 291-302). Results help in rationalizing, at a molecular level, the effects of homologues disaccharides on protein structure/dynamics experimentally observed. Furthermore, the effectiveness of disaccharides in bioprotection in terms of peculiar protein-matrix coupling is also discussed.

Molecular level probing of preferential hydration and its modulation by osmolytes through the use of pyranine complexed to hemoglobin

Two spectroscopic probes are used to expose molecular level changes in hydration shell water interactions that directly relate to such issues as preferential hydration and protein stability. The major focus of the present study is on the use of pyranine (HPT) fluorescence to probe as a function of added osmolytes (PEG, urea, trehalose, and magnesium), the extent to which glycerol is preferentially excluded from the hydration shell of free HPT and HPT localized in the diphosphoglycerate (DPG) binding site of hemoglobin in both solution and in Sol-Gel matrices. The pyranine study is complemented by the use of vibronic side band luminescence from the gadolinium cation that directly exposes the changes in hydrogen bonding between first and second shell waters as a function of added osmolytes. Together the results form the basis for a water partitioning model that can account for both preferential hydration and water/osmolyte-mediated conformational changes in protein structure.

Picosecond dynamics of T and R forms of aspartate transcarbamylase: A neutron scattering study

E. coli aspartate transcarbamylase (ATCase) is a 310 kDa allosteric enzyme which catalyses the first committed step in pyrimidine biosynthesis. The binding of its substrates, carbamylphosphate and aspartate, induces significant conformational changes. This enzyme shows homotropic cooperative interactions between the catalytic sites for the binding of aspartate. This property is explained by a quaternary structure transition from T state (aspartate low affinity) to R state (aspartate high affinity) accompanied by a 5% increase of radius of gyration of ATCase. The same quaternary structure change is observed upon binding of the bisubstrate analogue PALA (N-(phosphonacetyl)-L-aspartate. Owing to the large incoherent neutron scattering cross-section of the hydrogen atom and the abundance of this element in proteins, inelastic neutron scattering gives a global view of protein dynamics as sensed via the individual motions of its hydrogen atoms. We present neutron scattering results of the local dynamics (few angstroms), at short time (few tens of picoseconds), of ATCase in T and R forms. Compared to the T form, we observe an increased mobility of the protein in the R form that we associate to an increase of accessible surface area to the solvent. Beyond this specific result, this highlights the key role of the accessible surface area (ASA) in dynamic contribution to inelastic neutron data in the picosecond time scale. In particular, we want to stress out (i) that a difference at the picosecond time scale does not allow to conclude to a difference in the dynamics at a longer time scale and to address whether the T state is looser than the R state (ii) how challenging is, any comparison in terms of general dynamics (tense or relaxed) between dynamic values deduced from experimental neutron data on proteins with different sequences and therefore ASA. This caveat holds particularly when comparing dynamics of a mesophile with the corresponding extremophile.

Observation of fragile-to-strong dynamic crossover in protein hydration water

At low temperatures, proteins exist in a glassy state, a state that has no conformational flexibility and shows no biological functions. In a hydrated protein, at temperatures greater-- similar 220 K, this flexibility is restored, and the protein is able to sample more conformational substates, thus becoming biologically functional. This "dynamical" transition of protein is believed to be triggered by its strong coupling with the hydration water, which also shows a similar dynamic transition. Here we demonstrate experimentally that this sudden switch in dynamic behavior of the hydration water on lysozyme occurs precisely at 220 K and can be described as a fragile-to-strong dynamic crossover. At the fragile-to-strong dynamic crossover, the structure of hydration water makes a transition from predominantly high-density (more fluid state) to low-density (less fluid state) forms derived from the existence of the second critical point at an elevated pressure.

Role of solvent on protein-matrix coupling in MbCO embedded in water-saccharide systems: a Fourier transform infrared spectroscopy study

Embedding protein in sugar systems of low water content enables one to investigate the protein dynamic-structure function in matrixes whose rigidity is modulated by varying the content of residual water. Accordingly, studying the dynamics and structure thermal evolution of a protein in sugar systems of different hydration constitutes a tool for disentangling solvent rigidity from temperature effects. Furthermore, studies performed using different sugars may give information on how the detailed composition of the surrounding solvent affects the internal protein dynamics and structural evolution. In this work, we compare Fourier transform infrared spectroscopy measurements (300-20 K) on MbCO embedded in trehalose, sucrose, maltose, raffinose, and glucose matrixes of different water content. At all the water contents investigated, the protein-solvent coupling was tighter in trehalose than in the other sugars, thus suggesting a molecular basis for the trehalose peculiarity. These results are in line with the observation that protein-matrix phase separation takes place in lysozyme-lactose, whereas it is absent in lysozyme-trehalose systems; indeed, these behaviors may respectively be due to the lack or presence of suitable water-mediated hydrogen-bond networks, which match the protein surface to the surroundings. The above processes might be at the basis of pattern recognition in crowded living systems; indeed, hydration shells structural and dynamic matching is first needed for successful come together of interacting biomolecules.

Low-temperature glass transitions of quenched and annealed bovine serum albumin aqueous solutions

To investigate the glass transition behaviors of a 20% (w/w) aqueous solution of bovine serum albumin, heat capacities and enthalpy relaxation rates were measured by adiabatic calorimetry at temperatures ranging from 80 to 300 K. One series of measurements was carried out after quenching from 300 down to 80 K and another after annealing in 200-240 K. The quenched sample showed a heat capacity jump indicating a glass transition temperature T(g) = 170 K, and the annealed sample showed a smaller jump with the T(g) shifted toward the higher temperature side. The temperature dependence of the enthalpy relaxation rates for the quenched sample indicated the presence of two enthalpy relaxation effects: one at around 110 K and the other over a wide temperature range (120-190 K). The annealed sample showed three separate relaxation effects giving 1) T(g) = 110 K, 2) 135 K, and 3) temperature higher than 180 K, whereas nothing around 170 K. These effects were thought to originate, respectively, from the rearrangement motions of 1) primary hydrate water forming a direct hydrogen bond with the protein, 2) part of the internal water localized in the opening of a protein structure, and 3) the disordered region in the protein.

Dynamic Transition in tRNA is Solvent Induced

Dynamics of tRNA was studied using neutron scattering spectroscopy. Despite vast differences in the architecture and backbone structure of proteins and RNA, hydrated tRNA undergoes the dynamic transition at the same temperature as hydrated lysozyme. The similarity of the dynamic transition in RNA and proteins supports the idea that it is solvent induced. Because tRNA essentially has no methyl groups, the results also suggest that methyl groups are not the main contributor of the dynamic transition in biological macromolecules. However, they may explain strong differences in the dynamics of tRNA and lysozyme observed at low temperatures.

Picosecond-time-scale fluctuations of proteins in glassy matrices: the role of viscosity

Through elastic neutron scattering we investigated the fast dynamics of lysozyme in hydrated powder form or embedded in glycerol-water and glucose-water matrices. We calculated the relaxational contribution to the mean square displacements of protein hydrogen atoms. We found that the inverse of this quantity is linearly proportional to the logarithm of the viscosity of the solvent glassy matrix. This relationship suggests a close connection between the picosecond-time-scale dynamics of protein side chains and the solvent structural relaxation.

Liquid-like water confined in stacks of biological membranes at 200 k and its relation to protein dynamics

Confined water is of considerable current interest owing to its biophysical importance and relevance to cryopreservation. It can be studied in its amorphous or supercooled state in the "no-man's land", i.e., in the temperature range between 150 and 235 K, in which bulk water is always crystalline. Amorphous deuterium oxide (D(2)O) was obtained in the intermembrane spaces of a stack of purple membranes from Halobacterium salinarum by flash cooling to 77 K. Neutron diffraction showed that upon heating to 200 K the intermembrane water space decreased sharply with an associated strengthening of ice diffraction, indicating that water beyond the first membrane hydration layer flowed out of the intermembrane space to form crystalline ice. It was concluded that the confined water undergoes a glass transition at or below 200 K to adopt an ultraviscous liquid state from which it crystallizes to form ice as soon as it finds itself in an unconfined, bulk-water environment. Our results provide model-free evidence for translational diffusion of confined water in the no-man's land. Potential effects of the confined-water glass transition on nanosecond membrane dynamics were investigated by incoherent elastic neutron scattering experiments. These revealed no differences between flash-cooled and slow-cooled samples (in the latter, the intermembrane space at temperatures <250 K is occupied only by the first membrane hydration layers), with dynamical transitions at 150 and 260 K, but not at 200 K, suggesting that nanosecond membrane dynamics are not sensitive to the state of the water beyond the first hydration shell at cryotemperatures.

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.

“Glass Transition” near 200 K in the Bacterial Photosynthetic Reaction Center Protein Detected by Studying the Distances in the Transient P+QA- Radical Pair

The transient radical pair P(+)Q(A)(-) in the photosynthetic reaction center from Rhodobacter sphaeroides R26 was studied over a wide temperature range using out-of-phase electron spin-echo envelope modulation (ESEEM) spectroscopy. This method is sensitive to the magnetic dipole-dipole interaction between the two electron spins of the pair and allows precise determination of the distance in the pair P(+)Q(A)(-). The out-of-phase data were complemented by normal in-phase ESEEM spectra from the two stable radicals of P(+) and Q(A)(-). The results seem to indicate that the radical pair undergoes a noticeable molecular motion around 200 K that may be characterized by a change in the distance in the pair by approximately 0.3 nm. As the two cofactors, P(+) and Q(A)(-), are held in a well-defined relative position by the reaction center protein, this means that the protein becomes flexible at 200 K. This effect may be ascribed to a dynamic glass transition around 200 K. The relation with the temperature dependence of the back reaction of P(+)Q(A)(-) is discussed.

Molecular dynamics simulation of sucrose- and trehalose-coated carboxy-myoglobin

We performed a room temperature molecular dynamics (MD) simulation on a system containing 1 carboxy-myoglobin (MbCO) molecule in a sucrose-water matrix of identical composition (89% [sucrose/(sucrose + water)] w/w) as for a previous trehalose-water-MbCO simulation (Cottone et al., Biophys J 2001;80:931-938). Results show that, as for trehalose, the amplitude of protein atomic mean-square fluctuations, on the nanosecond timescale, is reduced with respect to aqueous solutions also in sucrose. A detailed comparison as a function of residue number evidences mobility differences along the protein backbone, which can be related to a different efficacy in bioprotection. Different heme pocket structures are observed in the 2 systems. The joint distribution of the magnitude of the electric field at the CO oxygen atom and of the angle between the field and the CO unit vector shows a secondary maximum in sucrose, absent in trehalose. This can explain the CO stretching band profile (A substates distribution) differences evidenced by infrared spectroscopy in sucrose- and trehalose-coated MbCO (Giuffrida et al., J Phys Chem B 2004;108:15415-15421), and in particular the appearance of a further substate in sucrose. Analysis of hydrogen bonds at the protein-solvent interface shows that the fraction of water molecules shared between the protein and the sugar is lower in sucrose than in trehalose, in spite of a larger number of water molecules bound to the protein in the former system, thus indicating a lower protein-matrix coupling, as recently observed by Fourier transform infrared (FTIR) experiments (Giuffrida et al., J Phys Chem B 2004;108:15415-15421).

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.

Fast dynamics and stabilization of proteins: binary glasses of trehalose and glycerol

We present elastic and inelastic incoherent neutron scattering data from a series of trehalose glasses diluted with glycerol. A strong correlation with recently published protein stability data in the same series of glasses illustrates that the dynamics at Q >or= 0.71 A(-1) and omega > 200 MHz are important to stabilization of horseradish peroxidase and yeast alcohol dehydrogenase in these glasses. To the best of our knowledge, this is the first direct evidence that enzyme stability in a room temperature glass depends upon suppressing these short-length scale, high-frequency dynamics within the glass. We briefly discuss the coupling of protein motions to the local dynamics of the glass. Also, we show that T(g) alone is not a good indicator for the protein stability in this series of glasses; the glass that confers the maximum room-temperature stability does not have the highest T(g).