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

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

Specific protein dynamics near the solvent glass transition assayed by radiation-induced structural changes

The nature of the dynamical coupling between a protein and its surrounding solvent is an important, yet open issue. Here we used temperature-dependent protein crystallography to study structural alterations that arise in the enzyme acetylcholinesterase upon X-ray irradiation at two temperatures: below and above the glass transition of the crystal solvent. A buried disulfide bond, a buried cysteine, and solvent exposed methionine residues show drastically increased radiation damage at 155 K, in comparison to 100 K. Additionally, the irradiation-induced unit cell volume increase is linear at 100 K, but not at 155 K, which is attributed to the increased solvent mobility at 155 K. Most importantly, we observed conformational changes in the catalytic triad at the active site at 155 K but not at 100 K. These changes lead to an inactive catalytic triad conformation and represent, therefore, the observation of radiation-inactivation of an enzyme at the atomic level. Our results show that at 155 K, the protein has acquired--at least locally--sufficient conformational flexibility to adapt to irradiation-induced alterations in the conformational energy landscape. The increased protein flexibility may be a direct consequence of the solvent glass transition, which expresses as dynamical changes in the enzyme's environment. Our results reveal the importance of protein and solvent dynamics in specific radiation damage to biological macromolecules, which in turn can serve as a tool to study protein flexibility and its relation to changes in a protein's environment.     

The 'glass transition' in protein dynamics: what it is, why it occurs, and how to exploit it

All proteins undergo a dramatic change in their dynamical properties at approximately 200 K. Above this temperature, their dynamic behavior is dominated by large-scale collective motions of bonded and nonbonded groups of atoms. At lower temperatures, simple harmonic vibrations predominate. The transition has been described as a 'glass transition' to emphasize certain similarities between the change in dynamic behavior of individual protein molecules and the changes in viscosity and other properties of liquids when they form a glass. The glass transition may reflect the intrinsic temperature dependence of the motions of atoms in the protein itself, in the bound solvent on the surface of the protein, or it may reflect contributions from both. Protein function is significantly altered below this transition temperature; a fact that can be exploited to trap normally unstable intermediates in enzyme-catalyzed reactions and stabilize them for periods long enough to permit their characterization by high-resolution protein crystallography.     

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

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

How a vicinal layer of solvent modulates the dynamics of proteins

The dynamics of a folded protein is studied in water and glycerol at a series of temperatures below and above their respective dynamical transition. The system is modeled in two distinct states whereby the protein is decoupled from the bulk solvent at low temperatures, and communicates with it through a vicinal layer at physiological temperatures. A linear viscoelastic model elucidates the less-than-expected increase in the relaxation times observed in the backbone dynamics of the protein. The model further explains the increase in the flexibility of the protein once the transition takes place and the differences in the flexibility under the different solvent environments. Coupling between the vicinal layer and the protein fluctuations is necessary to interpret these observations. The vicinal layer is postulated to form once a threshold for the volumetric fluctuations in the protein to accommodate solvents of different sizes is reached. Compensation of entropic-energetic contributions from the protein-coupled vicinal layer quantifies the scaling of the dynamical transition temperatures in various solvents. The protein adapts different conformational routes for organizing the required coupling to a specific solvent, which is achieved by adjusting the amount of conformational jumps in the surface-group dihedrals.     

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.     

Glycerol effects on protein flexibility: a tryptophan phosphorescence study.

In exploring the dynamic properties of protein structure, numerous studies have focussed on the dependence of structural fluctuations on solvent viscosity, but the emerging picture is still not well defined. Exploiting the sensitivity of the phosphorescence lifetime of tryptophan to the viscosity of its environment we have used the delayed emission as an intrinsic probe of protein flexibility and investigated the effects of glycerol as a viscogenic cosolvent. The phosphorescence lifetime of alcohol dehydrogenase, alkaline phosphatase, apoazurin and RNase T1, as a function of glycerol concentration was studied at various temperatures. Flexibility data, which refer to rather rigid sites of the globular structures, point out that, for some concentration ranges glycerol, effects on the rate of structural fluctuations of alcohol dehydrogenase and RNase T1 do not obey Kramers' a power law on solvent viscosity and emphasize that cosolvent-induced structural changes can be important, even for inner cores of the macromolecule. When the data is analyzed in terms of Kramers' model, for the temperature range 0-30 degrees C one derives frictional coefficients that are relatively large (0.6-0.7) for RNase T1, where the probe is in a flexible region near the surface of the macromolecule and much smaller, less than 0.2, for the rigid sites of the other proteins. For the latter sites the frictional coefficient rises sharply between 40 and 60 degrees C, and its value correlates weakly with molecular parameters such as the depth of burial or the rigidity of a particular site. For RNase T1, coupling to solvent viscosity increases at subzero temperatures, with the coefficient becoming as large as 1 at -20 degrees C.(ABSTRACT TRUNCATED AT 250 WORDS)     

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

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

Conformational change in the activation of lipase: an analysis in terms of low-frequency normal modes.

The interfacial activation of Rhizomucor miehei lipase (RmL) involves the motion of an alpha-helical region (residues 82-96) which acts as a "lid" over the active site of the enzyme, undergoing a displacement from a "closed" to an "open" conformation upon binding of substrate. Normal mode analyses performed in both low and high dielectric media reveal that low-frequency vibrational modes contribute significantly to the conformational transition between the closed and open conformations. In these modes, the lid displacement is coupled to local motions of active site loops as well as global breathing motions. Atomic fluctuations of the first hinge of the lid (residues 83-84) are substantially larger in the low dielectric medium than in the high dielectric medium. Our results also suggest that electrostatic interactions of Arg86 play an important role in terms of both the intrinsic stability of the lid and its displacement, through enhancement of hinge mobility in a high dielectric medium. Additional calculations demonstrate that the observed patterns of atomic fluctuations are an intrinsic feature of the protein structure and not dependent on the nature of specific energy minima.     

Switch from conventional to distributed kinetics in the bacteriorhodopsin photocycle

Below 195 K, the bacteriorhodopsin photocycle could not be adequately described with exponential kinetics [Dioumaev, A. K., and Lanyi, J. K. (2007) Proc. Natl. Acad. Sci. U.S.A. 104, 9621-9626] but required distributed kinetics, previously found in hemoglobin and myoglobin at temperatures below the vitrification point of the surrounding solvent. The aim of this study is to determine which factors cause the switch from this low-temperature regime to the conventional kinetics observed at ambient temperature. The photocycle was monitored by time-resolved FTIR between 180 and 280 K, using the D96N mutant. Depending on the temperature, decay and temporal redistribution of two or three intermediates (L, M, and N) were observed. Above approximately 245 K, an abrupt change in the kinetic behavior of the photocycle takes place. It does not affect the intermediates present but greatly accelerates their decay. Below approximately 240 K, a kinetic pattern with partial decay that cannot be explained by conventional kinetics, but suggesting distributed kinetics, was dominant, while above approximately 250 K, there were no significant deviations from exponential behavior. The approximately 245 K critical point is >/=10 K below the freezing point of interbilayer water, and we were unable to correlate it with any FTIR-detectable transition of the lipids. Therefore, we attribute the change from distributed to conventional kinetics to a thermodynamic phase transition in the protein. Most probably, it is related to the freezing and thawing of internal fluctuations of the protein, known as the dynamic phase transition, although in bacteriorhodopsin the latter is usually believed to take place at least 15 K below the observed critical temperature of approximately 245 K.     

The adaptive nature of protein residue networks

Protein residue networks PRNs are used to describe proteins. These networks are usually based on an average structure for the protein. However, proteins are dynamic entities that are affected by their surroundings. In this work, we study the effect of temperatures above and below the protein dynamical transition temperature(≈200 K), on three important network parameters gleaned from weighted PRNs for the solvated β‐lactamase inhibitory protein BLIP: the betweenness centrality B, the closeness centrality C, and the clustering coefficient CC. The B and C values will be extracted for each node from PRNs at six different temperatures: 150 K, 180 K, 200 K, 220 K, 250 K, and 310 K respectively. The average value for the CC for each PRN will also be calculated at each temperature, respectively. We find that at temperatures ≤200 K, the network nodes with the most significant B and C values tend to have lower relative solvent accessibility RSA values, and to fall within the protein secondary structure elements (α helices and β sheets). At temperatures >200 K, the significant nodes in terms of B and C tend to have larger RSA values, and to fall on the connecting loops in the protein. The average CC decreases in value for the PRNs up to 200 K, and then remains basically constant above 200 K. This clearly shows that any conclusions based on static PRNs should be handled with care. The dynamic nature of proteins and its coupling to the surrounding environment should be taken into consideration when using the PRN paradigm. Proteins 2017; 85:917–923. © 2016 Wiley Periodicals, Inc.     

Effects of temperature on protein structure and dynamics: X-ray crystallographic studies of the protein ribonuclease-A at nine different temperatures from 98 to 320 K.

Structures using X-ray diffraction data collected to 1.5-A resolution have been determined for the protein ribonuclease-A at nine different temperatures ranging from 98 to 320 K. It is determined that the protein molecule expands slightly (0.4% per 100 K) with increasing temperature and that this expansion is linear. The expansion is due primarily to subtle repacking of the molecule, with exposed and mobile loop regions exhibiting the largest movements. Individual atomic Debye-Waller factors exhibit predominantly biphasic behavior, with a small positive slope at low temperatures and a larger positive slope at higher temperatures. The break in this curve occurs at a characteristic temperature of 180-200 K, perhaps indicative of fundamental changes in the dynamical structure of the surrounding protein solvent. The distribution of protein Debye-Waller factors is observed to broaden as well as shift to higher values as the temperature is increased.     

On the temperature and pressure dependence of a range of properties of a type of water model commonly used in high-temperature protein unfolding simulations

Molecular dynamics simulations of protein folding and unfolding are often carried out at temperatures (400-600 K) that are much higher than physiological or room temperature to speed up the (un)folding process. Use of such high temperatures changes both the protein and solvent properties considerably, compared to physiological or room temperature. Water models designed for use in conjunction with biomolecules, such as the simple point charge (SPC) model, have generally been calibrated at room temperature and pressure. To determine the distortive effect of high simulation temperatures on the behavior of such "room temperature" water models, the structural, dynamic, and thermodynamic properties of the much-used SPC water model are investigated in the temperature range from 300 to 500 K. Both constant pressure and constant volume conditions, as used in protein simulations, were analyzed. We found that all properties analyzed change markedly with increasing temperature, but no phase transition in this temperature range was observed.     

Temperature dependence of the structure and dynamics of myoglobin. A simulation approach

The results of simulations of the structure and internal motions of carbonomonoxymyoglobin (MbCO) at two different temperatures (325 and 80 K) are presented and compared with experimental data. Properties calculated from the 120 ps trajectory at 325 K are used as a reference in the analysis of the motion of the protein at 80 K. Three separate 80 K molecular dynamics trajectories were calculated; they were started with different coordinate sets from the 325 K simulation and the lower temperature was achieved by scaling the velocities. The simulations yield results for the structural changes between 325 and 80 K that are in general accord with those from X-ray data. Both the experimental and calculated radii of gyration, distances from the center of mass and main-chain difference distance matrices show that there is a significant but inhomogeneous shrinkage with decreasing temperature. For the atomic fluctuations, by contrast, the calculated temperature dependence is very different from the X-ray results; i.e. the calculated root-mean-square backbone fluctuations decrease to 0.11 A at 80 K from 0.51 A at 325 K, while the fluctuations obtained from the X-ray B factors go from 0.56 A at 260 K to 0.47 A at 80 K. The smaller temperature dependence of the B factors suggests that there is significant conformational disorder in MbCO crystals at lower temperatures. This is in accord with the simulation results, which show that the protein is trapped in restricted regions of conformational space at 80 K, while at 325 K a much larger region is accessible to the protein. Analysis of the fluctuations at 325 K and 80 K shows that the room temperature flexibility of the protein is determined by the mobility of the loop regions and by side-chain torsional motions (in accord with earlier simulation results), while the low temperature fluctuations involve motion within a single well. Examination of the calculated iron atom fluctuations and comparison with Mossbauer data show good agreement. It is found that the dominant contribution to the iron motion arises from heme sliding; motion of the iron relative to the heme are much smaller.     

Temperature dependence of the free energy landscape of the src-SH3 protein domain

The temperature dependence of the free energy landscape of the src-SH3 protein domain is investigated through fully atomic simulations in explicit solvent. Simulations are performed above and below the folding transition temperature, enabling an analysis of both protein folding and unfolding. The transition state for folding and unfolding, identified from the free energy surfaces, is found to be very similar, with structure in the central hydrophobic sheet and little structure throughout the rest of the protein. This is a result of a polarized folding (unfolding) mechanism involving early formation (late loss) of the central hydrophobic sheet at the transition state. Unfolding simulations map qualitatively well onto low-temperature free energy surfaces but appear, however, to miss important features observed in folding simulations. In particular, details of the folding mechanism involving the opening and closing of the hydrophobic core are not captured by unfolding simulations performed under strongly denaturing conditions. In addition, free energy surfaces at high temperatures do not display a desolvation barrier found at lower temperatures, involving the expulsion of water molecules from the hydrophobic core.     

Temperature-dependent vibrational and conformational dynamics of photosystem II membrane fragments from spinach investigated by elastic and inelastic neutron scattering

Vibrational and conformational protein dynamics of photosystem II (PS II) membrane fragments from spinach were investigated by elastic and inelastic incoherent neutron scattering (EINS and IINS). As to the EINS experiments, the average atomic mean square displacement values of PS II membrane fragments hydrated at a relative humidity of 57% exhibit a dynamical transition at ~230K. In contrast, the dynamical transition was absent at a relative humidity of 44%. These findings are in agreement with previous studies which reported a "freezing" of protein mobility due to dehydration (Pieper et al. (2008) Eur. Biophys. J. 37: 657-663) and its correlation with an inhibition of electron transfer from Q(A)(-) to Q(B) (Kaminskaya et al. (2003) Biochemistry 42, 8119-8132). IINS spectra of a sample hydrated at a relative humidity of 57% show a distinct Boson peak at ~7.5meV at 20K, which shifts towards lower energy values upon temperature increase to 250K. This unexpected effect is interpreted in terms of a "softening" of the protein matrix along with the onset of conformational protein dynamics as revealed by the EINS experiments. Information on the density of vibrational states of pigment-protein complexes is important for a realistic calculation of excitation energy transfer kinetics and spectral lineshapes and is often routinely obtained by optical line-narrowing spectroscopy at liquid helium temperature. The data presented here demonstrate that IINS is a valuable experimental tool in determining the density of vibrational states not only at cryogenic, but also at nearly physiological temperatures up to 250K. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.     

Slow relaxation process in DNA at different levels of hydration

There are many speculations about the dynamic transition observed in hydrated bio-polymers at temperatures T 200 – 230 K being an important factor for enabling of their functions. The transition shows up as a sharp increase of atomic mean-squared displacements above this temperature. The nature of the dynamic transition is not yet clear. Using inelastic neutron scattering we show in this Note that the transition in DNA is related to the appearance of a slow relaxation process. Decrease in the hydration level suppresses the process and the dynamic transition. It is found that, in terms of dynamics, the decrease in water content is similar in effect to a decrease in temperature. The obtained results support the idea that the dynamic transition is mediated by the water of hydration since bulk water has a dynamic transition around the same temperature.     

Enzyme activity below the dynamical transition at 220 K.

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

The influence of solvent composition on global dynamics of human butyrylcholinesterase powders: a neutron-scattering study

A major result of incoherent elastic neutron-scattering experiments on protein powders is the strong dependence of the intramolecular dynamics on the sample environment. We performed a series of incoherent elastic neutron-scattering experiments on lyophilized human butyrylcholinesterase (HuBChE) powders under different conditions (solvent composition and hydration degree) in the temperature range from 20 to 285 K to elucidate the effect of the environment on the enzyme atomic mean-square displacements. Comparing D(2)O- with H(2)O-hydrated samples, we were able to investigate protein as well as hydration water molecular dynamics. HuBChE lyophilized from three distinct buffers showed completely different atomic mean-square displacements at temperatures above approximately 200 K: a salt-free sample and a sample containing Tris-HCl showed identical small-amplitude motions. A third sample, containing sodium phosphate, displayed highly reduced mean-square displacements at ambient temperature with respect to the other two samples. Below 200 K, all samples displayed similar mean-square displacements. We draw the conclusion that the reduction of intramolecular protein mean-square displacements on an Angstrom-nanosecond scale by the solvent depends not only on the presence of salt ions but also on their type.