A Dynamical Transition in the Protein Myoglobin Observed by Infrared Vibrational Echo Experiments

Role of hydration water in the onset of protein structural dynamics

Proteins are the molecular workhorses in a living organism. Their 3D structures are animated by a multitude of equilibrium fluctuations and specific out-of-equilibrium motions that are required for proteins to be biologically active. When studied as a function of temperature, functionally relevant dynamics are observed at and above the so-called protein dynamical transition (~240 K) in hydrated, but not in dry proteins. In this review we present and discuss the main experimental and computational results that provided evidence for the dynamical transition, with a focus on the role of hydration water dynamics in sustaining functional protein dynamics. The coupling and mutual influence of hydration water dynamics and protein dynamics are discussed and the hypotheses illustrated that have been put forward to explain the physical origin of their onsets.

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.

Dynamical transition in a small helical peptide and its implication for vibrational energy transport

The two-dimensional infrared spectrum of an octameric helical peptide in chloroform was measured as a function of temperature. Isotope labeling of the carbonyl group of one of the amino acids was used to obtain information for an isolated vibration. The antidiagonal width of the 2D-IR signal, which is a measure of the homogeneous dephasing time T(2), is constant from 220 to 260 K (within experimental error), and increases steeply above. The homogeneous dephasing time of the carbonyl vibration is attributed to the flexibility of the system and/or its immediate surrounding. The system undergoes a dynamical transition at about 270 K, with similarities to the protein dynamical transition. Furthermore, the temperature dependence of the antidiagonal width strongly resembles that of the efficiency of vibrational energy transport along the helix, which has been studied in a recent paper (J. Phys. Chem. B 2008, 112, 15487). The connection between the two processes, structural flexibility and energy transport mechanism, is discussed.

Applications of 2D IR spectroscopy to peptides, proteins, and hydrogen-bond dynamics

Following a survey of 2D IR principles, this article describes recent experiments on the hydrogen-bond dynamics of small ions, amide-I modes, nitrile probes, peptides, reverse transcriptase inhibitors, and amyloid fibrils.

Native and unfolded cytochrome c--comparison of dynamics using 2D-IR vibrational echo spectroscopy

Unfolded vs native CO-coordinated horse heart cytochrome c (h-cyt c) and a heme axial methionine mutant cyt c552 from Hydrogenobacter thermophilus ( Ht-M61A) are studied by IR absorption spectroscopy and ultrafast 2D-IR vibrational echo spectroscopy of the CO stretching mode. The unfolding is induced by guanidinium hydrochloride (GuHCl). The CO IR absorption spectra for both h-cyt c and Ht-M61A shift to the red as the GuHCl concentration is increased through the concentration region over which unfolding occurs. The spectra for the unfolded state are substantially broader than the spectra for the native proteins. A plot of the CO peak position vs GuHCl concentration produces a sigmoidal curve that overlays the concentration-dependent circular dichroism (CD) data of the CO-coordinated forms of both Ht-M61A and h-cyt c within experimental error. The coincidence of the CO peak shift curve with the CD curves demonstrates that the CO vibrational frequency is sensitive to the structural changes induced by the denaturant. 2D-IR vibrational echo experiments are performed on native Ht-M61A and on the protein in low- and high-concentration GuHCl solutions. The 2D-IR vibrational echo is sensitive to the global protein structural dynamics on time scales from subpicosecond to greater than 100 ps through the change in the shape of the 2D spectrum with time (spectral diffusion). At the high GuHCl concentration (5.1 M), at which Ht-M61A is essentially fully denatured as judged by CD, a very large reduction in dynamics is observed compared to the native protein within the approximately 100 ps time window of the experiment. The results suggest the denatured protein may be in a glassy-like state involving hydrophobic collapse around the heme.

Viscosity-dependent protein dynamics

Spectrally resolved stimulated vibrational echo spectroscopy is used to investigate the dependence of fast protein dynamics on bulk solution viscosity at room temperature in four heme proteins: hemoglobin, myoglobin, a myoglobin mutant with the distal histidine replaced by a valine (H64V), and a cytochrome c552 mutant with the distal methionine replaced by an alanine (M61A). Fructose is added to increase the viscosity of the aqueous protein solutions over many orders of magnitude. The fast dynamics of the four globular proteins were found to be sensitive to solution viscosity and asymptotically approached the dynamical behavior that was previously observed in room temperature sugar glasses. The viscosity-dependent protein dynamics are analyzed in the context of a viscoelastic relaxation model that treats the protein as a deformable breathing sphere. The viscoelastic model is in qualitative agreement with the experimental data but does not capture sufficient system detail to offer a quantitative description of the underlying fluctuation amplitudes and relaxation rates. A calibration method based on the near-infrared spectrum of water overtones was constructed to accurately determine the viscosity of small volumes of protein solutions.

Ultrafast dynamics of myoglobin without the distal histidine: stimulated vibrational echo experiments and molecular dynamics simulations

Ultrafast protein dynamics of the CO adduct of a myoglobin mutant with the polar distal histidine replaced by a nonpolar valine (H64V) have been investigated by spectrally resolved infrared stimulated vibrational echo experiments and molecular dynamics (MD) simulations. In aqueous solution at room temperature, the vibrational dephasing rate of CO in the mutant is reduced by approximately 50% relative to the native protein. This finding confirms that the dephasing of the CO vibration in the native protein is sensitive to the interaction between the ligand and the distal histidine. The stimulated vibrational echo observable is calculated from MD simulations of H64V within a model in which vibrational dephasing is driven by electrostatic forces. In agreement with experiment, calculated vibrational echoes show slower dephasing for the mutant than for the native protein. However, vibrational echoes calculated for H64V do not show the quantitative agreement with measurements demonstrated previously for the native protein.

Probing dynamics of complex molecular systems with ultrafast 2D IR vibrational echo spectroscopy

Ultrafast 2D IR vibrational echo spectroscopy is described and a number of experimental examples are given. Details of the experimental method including the pulse sequence, heterodyne detection, and determination of the absorptive component of the 2D spectrum are outlined. As an initial example, the 2D spectrum of the stretching mode of CO bound to the protein myoglobin (MbCO) is presented. The time dependence of the 2D spectrum of MbCO, which is caused by protein structural evolution, is presented and its relationship to the frequency-frequency correlation function is described and used to make protein structural assignments based on comparisons to molecular dynamics simulations. The 2D vibrational echo experiments on the protein horseradish peroxidase are presented. The time dependence of the 2D spectra of the enzyme in the free form and with a substrate bound at the active site are compared and used to examine the influence of substrate binding on the protein's structural dynamics. The application of 2D vibrational echo spectroscopy to the study of chemical exchange under thermal equilibrium conditions is described. 2D vibrational echo chemical exchange spectroscopy is applied to the study of formation and dissociation of organic solute-solvent complexes and to the isomerization around a carbon-carbon single bond of an ethane derivative.

The Influence of Aqueous versus Glassy Solvents on Protein Dynamics: Vibrational Echo Experiments and Molecular Dynamics Simulations

Spectrally resolved infrared stimulated vibrational echo measurements are used to measure the vibrational dephasing of the CO stretching mode of carbonmonoxy-hemoglobin (HbCO), a myoglobin mutant (H64V), and a bacterial cytochrome c(552) mutant (Ht-M61A) in aqueous solution and trehalose glasses. The vibrational dephasing of the heme-bound CO is significantly slower for all three proteins embedded in trehalose glasses compared to that of aqueous protein solutions. All three proteins exhibit persistent but notably slower spectral diffusion when the protein surface is fixed by the glassy solvent. Frequency-frequency correlation functions (FFCFs) of the CO are extracted from the vibrational echo data to reveal that the structural dynamics, as sensed by the CO, of the three proteins in trehalose and aqueous solution are dominated by fast (tens of femtoseconds), motionally narrowed fluctuations. MD simulations of H64V in dynamic and "static" water are presented as models of the aqueous and glassy environments. FFCFs are calculated from the H64V simulations and qualitatively reproduce the important features of the experimentally extracted FFCFs. The suppression of long time scale (picoseconds to tens of picoseconds) frequency fluctuations (spectral diffusion) in the glassy solvent is the result of a damping of atomic displacements throughout the protein structure and is not limited to structural dynamics that occur only at the protein surface. The analysis provides evidence that some dynamics are coupled to the hydration shell of water, supporting the idea that the bioprotection offered by trehalose is due to its ability to immobilize the protein surface through a thin, constrained layer of water.

Enzyme activity and flexibility at very low hydration

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

Bulk-solvent and hydration-shell fluctuations, similar to alpha- and beta-fluctuations in glasses, control protein motions and functions

The concept that proteins exist in numerous different conformations or conformational substates, described by an energy landscape, is now accepted, but the dynamics is incompletely explored. We have previously shown that large-scale protein motions, such as the exit of a ligand from the protein interior, follow the dielectric fluctuations in the bulk solvent. Here, we demonstrate, by using mean-square displacements (msd) from Mossbauer and neutron-scattering experiments, that fluctuations in the hydration shell control fast fluctuations in the protein. We call the first type solvent-slaved or alpha-fluctuations and the second type hydration-shell-coupled or beta-fluctuations. Solvent-slaved motions are similar to the alpha-fluctuations in glasses. Their temperature dependence can be approximated by a Vogel-Tammann-Fulcher relation and they are absent in a solid environment. Hydration-shell-coupled fluctuations are similar to the beta-relaxation in glasses. They can be approximated by a Ferry or an Arrhenius relation, are much reduced or absent in dehydrated proteins, and occur in hydrated proteins even if embedded in a solid. They can be responsible for internal processes such as the migration of ligands within myoglobin. The existence of two functionally important fluctuations in proteins, one slaved to bulk motions and the other coupled to hydration-shell fluctuations, implies that the environment can control protein functions through different avenues and that no real protein transition occurs at approximately 200 K. The large number of conformational substates is essential; proteins cannot function without this reservoir of entropy, which resides mainly in the hydration shell.

Vibrational echo correlation spectroscopy probes of hydrogen bond dynamics in water and methanol

Multidimensional vibrational echo correlation spectroscopy with full phase resolution is used to measure hydrogen bond dynamics in water and methanol. The OD hydroxyl stretches of methanol-OD oligomers in CCl(4) and HOD in H(2)O are studied using the shortest mid-IR pulses (<50 fs, <4 cycles of light) produced to date. The pulses have sufficient spectral bandwidth to span the very broad (>400 cm(-1)) spectrum of the 0-1 and 1-2 transitions. Hydrogen bond population dynamics are extricated with exceptional detail in MeOD oligomers because the different hydrogen bonded species are spectrally distinct. The experimental results along with detailed calculations indicate the strongest hydrogen bonds are selectively broken through a non-equilibrium relaxation pathway following vibrational relaxation of the hydroxyl stretch. The correlation spectra are also a sensitive probe of the fluctuations in water and provide a stringent test of water models that are widely used in simulations of aqueous systems. The analysis of the 2D band shapes demonstrates that different hydrogen bonded species are subject to distinct (wavelength dependent) ultrafast (~100 fs) local fluctuations and essentially identical slow (0.4 and ~2 ps) structural rearrangements. Observation of wavelength dependent dynamics demonstrates that standard theoretical approaches assuming Gaussian fluctuations cannot adequately describe water dynamics.

Myoglobin-CO substate structures and dynamics: multidimensional vibrational echoes and molecular dynamics simulations

Spectrally resolved infrared stimulated vibrational echo data were obtained for sperm whale carbonmonoxymyoglobin (MbCO) at 300 K. The measured dephasing dynamics of the CO ligand are in agreement with dephasing dynamics calculated with molecular dynamics (MD) simulations for MbCO with the residue histidine-64 (His64) having its imidazole epsilon nitrogen protonated (N(epsilon)-H). The two conformational substate structures B(epsilon) and R(epsilon) observed in the MD simulations are assigned to the spectroscopic A(1) and A(3) conformational substates of MbCO, respectively, based on the agreement between the experimentally measured and calculated dephasing dynamics for these substates. In the A(1) substate, the N(epsilon)-H proton and N(delta) of His64 are approximately equidistant from the CO ligand, while in the A(3) substate, the N(epsilon)-H of His64 is oriented toward the CO, and the N(delta) is on the surface of the protein. The MD simulations show that dynamics of His64 represent the major source of vibrational dephasing of the CO ligand in the A(3) state on both femtosecond and picosecond time scales. Dephasing in the A(1) state is controlled by His64 on femtosecond time scales, and by the rest of the protein and the water solvent on longer time scales.

Myoglobin-CO conformational substate dynamics: 2D vibrational echoes and MD simulations

Two-dimensional (2D) infrared vibrational echoes were performed on horse heart carbonmonoxymyoglobin (MbCO) in water over a range of temperatures. The A(1) and A(3) conformational substates of MbCO are found to have different dephasing rates with different temperature dependences. A frequency-frequency correlation function derived from molecular dynamics simulations on MbCO at 298 K is used to calculate the vibrational echo decay. The calculated decay shows substantial agreement with the experimentally measured decays. The 2D vibrational echo probes protein dynamics and provides an observable that can be used to test structural assignments for the MbCO conformational substates.

Fast protein dynamics probed with infrared vibrational echo experiments

IR vibrational echo experiments are used to study dynamics in myoglobin (Mb) by investigating the dephasing of the CO-stretching mode of CO bound at the active site of the protein (Mb-CO). The temperature dependence and the viscosity dependence of Mb-CO pure dephasing have been measured in several solvents. In low-temperature, glassy solvents, the pure dephasing has a power law temperature dependence, T(1.3), that reflects glasslike protein dynamics. In liquids, the temperature dependence is much steeper and arises from a combination of pure temperature dependence and the influence of decreasing solvent viscosity with increasing temperature. As the solvent viscosity decreases, the ability of the protein's surface to undergo topological fluctuations increases, which in turn increases the internal protein-structural fluctuations. The protein-structural motions are coupled to the CO bound at the active site by electric field fluctuations that accompany movements of polar residues. The dynamic electric field-coupling mechanism is tested by observing differences in the temperature dependence of the pure dephasing of Mb-CO mutations.

Multidimensional femtosecond correlation spectroscopies of electronic and vibrational excitations

Femtosecond visible and infrared analogues of multiple-pulse nuclear magnetic resonance techniques provide novel snapshot probes into the structure and electronic and vibrational dynamics of complex molecular assemblies such as photosynthetic antennae, proteins, and hydrogen-bonded liquids. A classical-oscillator description of these spectroscopies in terms of interacting quasiparticles (rather than transitions among global eigenstates) is developed and sets the stage for designing new pulse sequences and inverting the multidimensional signals to yield molecular structures. Considerable computational advantages and a clear physical insight into the origin of the response and the relevant coherence sizes are provided by a real-space analysis of the underlying coherence-transfer pathways in Liouville space.