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.     

Cytochrome c552 mutants: structure and dynamics at the active site probed by multidimensional NMR and vibration echo spectroscopy

Spectrally resolved infrared stimulated vibrational echo experiments are used to measure the vibrational dephasing of a CO ligand bound to the heme cofactor in two mutated forms of the cytochrome c552 from Hydrogenobacter thermophilus. The first mutant (Ht-M61A) is characterized by a single mutation of Met61 to an Ala (Ht-M61A), while the second variant is doubly modified to have Gln64 replaced by an Asn in addition to the M61A mutation (Ht-M61A/Q64N). Multidimensional NMR experiments determined that the geometry of residue 64 in the two mutants is consistent with a non-hydrogen-bonding and hydrogen-bonding interaction with the CO ligand for Ht-M61A and Ht-M61A/Q64N, respectively. The vibrational echo experiments reveal that the shortest time scale vibrational dephasing of the CO is faster in the Ht-M61A/Q64N mutant than that in Ht-M61A. Longer time scale dynamics, measured as spectral diffusion, are unchanged by the Q64N modification. Frequency-frequency correlation functions (FFCFs) of the CO are extracted from the vibrational echo data to confirm that the dynamical difference induced by the Q64N mutation is primarily an increase in the fast (hundreds of femtoseconds) frequency fluctuations, while the slower (tens of picoseconds) dynamics are nearly unaffected. We conclude that the faster dynamics in Ht-M61A/Q64N are due to the location of Asn64, which is a hydrogen bond donor, above the heme-bound CO. A similar difference in CO ligand dynamics has been observed in the comparison of the CO derivative of myoglobin (MbCO) and its H64V variant, which is caused by the difference in axial residue interactions with the CO ligand. The results suggest a general trend for rapid ligand vibrational dynamics in the presence of a hydrogen bond donor.     

Dynamics of hemoglobin in human erythrocytes and in solution: influence of viscosity studied by ultrafast vibrational echo experiments

Ultrafast spectrally resolved stimulated vibrational echo experiments are used to measure the vibrational dephasing of the CO stretching mode of hemoglobin-CO (HbCO) inside living human erythrocytes (red blood cells), in liquid solutions, and in a glassy matrix. A method is presented to overcome the adverse impact on the vibrational echo signal from the strong light scattering caused by the cells. The results from the cytoplasmic HbCO are compared to experiments on aqueous HbCO samples prepared in different buffers, solutions containing low and high concentrations of glycerol, and in a solid trehalose matrix. Measurements are also presented that provide an accurate determination of the viscosity at the very high Hb concentration that is found inside the cells. It is demonstrated that the dynamics of the protein, as sensed by the CO ligand, are the same inside the erythrocytes and in aqueous solution and are independent of the viscosity. In solutions that are predominantly glycerol, the dynamics are modified somewhat but are still independent of viscosity. The experiments in trehalose give the dynamics at infinite viscosity and are used to separate the viscosity-dependent dynamics from the viscosity-independent dynamics. Although the HbCO dynamics are the same in the red blood cell and in the equivalent aqueous solutions, differences in the absorption spectra show that the distribution of a protein's equilibrium substates is sensitive to small pH differences.     

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.     

Vibrational echo experiments on red blood cells: Comparison of the dynamics of cytoplasmic and aqueous hemoglobin

Ultrafast spectrally resolved stimulated vibrational echo experiments measure the dephasing of the CO stretching mode of hemoglobin-CO (Hb-CO) inside living human erythrocytes (red blood cells). A method is presented to overcome the adverse impact on the vibrational echo signal from the strong light scattering caused by the cells. The results are compared to experiments on Hb-CO aqueous solutions. It is demonstrated that the dynamics of the protein as sensed by the CO ligand are the same inside the erythrocytes and in aqueous solution, but differences in the absorption spectra show that the cell affects the protein's potential energy surface.     

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.     

Temperature dependent equilibrium native to unfolded protein dynamics and properties observed with IR absorption and 2D IR vibrational echo experiments

Dynamic and structural properties of carbonmonoxy (CO)-coordinated cytochrome c(552) from Hydrogenobacter thermophilus (Ht-M61A) at different temperatures under thermal equilibrium conditions were studied with infrared absorption spectroscopy and ultrafast two-dimensional infrared (2D IR) vibrational echo experiments using the heme-bound CO as the vibrational probe. Depending on the temperature, the stretching mode of CO shows two distinct bands corresponding to the native and unfolded proteins. As the temperature is increased from low temperature, a new absorption band for the unfolded protein grows in and the native band decreases in amplitude. Both the temperature-dependent circular dichroism and the IR absorption area ratio R(A)(T), defined as the ratio of the area under the unfolded band to the sum of the areas of the native and unfolded bands, suggest a two-state transition from the native to the unfolded protein. However, it is found that the absorption spectrum of the unfolded protein increases its inhomogeneous line width and the center frequency shifts as the temperature is increased. The changes in line width and center frequency demonstrate that the unfolding does not follow simple two-state behavior. The temperature-dependent 2D IR vibrational echo experiments show that the fast dynamics of the native protein are virtually temperature independent. In contrast, the fast dynamics of the unfolded protein are slower than those of the native protein, and the unfolded protein fast dynamics and at least a portion of the slower dynamics of the unfolded protein change significantly, becoming faster as the temperature is raised. The temperature dependence of the absorption spectrum and the changes in dynamics measured with the 2D IR experiments confirm that the unfolded ensemble of conformers continuously changes its nature as unfolding proceeds, in contrast to the native state, which displays a temperature-independent distribution of structures.     

Vibrational relaxation of CH3I in the gas phase and in solution

Transient electronic absorption measurements reveal the vibrational relaxation dynamics of CH(3)I following excitation of the C-H stretch overtone in the gas phase and in liquid solutions. The isolated molecule relaxes through two stages of intramolecular vibrational relaxation (IVR), a fast component that occurs in a few picoseconds and a slow component that takes place in about 400 ps. In contrast, a single 5-7 ps component of IVR precedes intermolecular energy transfer (IET) to the solvent, which dissipates energy from the molecule in 50 ps, 44 ps, and 16 ps for 1 M solutions of CH(3)I in CCl(4), CDCl(3), and (CD(3))(2)CO, respectively. The vibrational state structure suggests a model for the relaxation dynamics in which a fast component of IVR populates the states that are most strongly coupled to the initially excited C-H stretch overtone, regardless of the environment, and the remaining, weakly coupled states result in a secondary relaxation only in the absence of IET.     

Fast enzyme dynamics at the active site of formate dehydrogenase

The role of femtosecond-picosecond structural dynamics of proteins in enzyme-catalyzed reactions is a hotly debated topic. We report infrared photon echo measurement of the formate dehydrogenase-NAD+-azide ternary complex. In contrast to earlier studies of protein dynamics, the data show complete spectral diffusion on the femtosecond-picosecond time scale with no static heterogeneity. This result indicates that this transition-state analogue complex completely samples the distribution of structures that determine the distribution of azide vibrational frequencies within a few picoseconds and that there are no slower motions that perturb the H-bond network at the active site.     

Anisotropic structural relaxation and its correlation with the excess energy diffusion in the incipient process of photodissociated MbCO: high-resolution analysis via ensemble perturbation method

The effectiveness of the ensemble perturbation method, in which many pairs of perturbed and unperturbed molecular dynamics simulations are executed for the ensemble average, has been demonstrated by calculating the subtle anisotropic structural change of carbonmonoxy myoglobin (MbCO) triggered by ligand photolysis. The results show that Mb largely expands in the direction perpendicular to the heme plane and slightly contracts in the horizontal one. This agrees well with the report in the transient grating experiment. In addition, it is suggested that the expansion contributes strongly to the fast energy-transfer process to the water solvent because it is undergone almost within several picoseconds. The mechanical work done on the solvent by the expansion within 1 ps was thermodynamically estimated to be 4.8 kcal/mol.     

Water dynamics--the effects of ions and nanoconfinement

Hydrogen bond dynamics of water in highly concentrated NaBr salt solutions and reverse micelles are studied using ultrafast 2D-IR vibrational echo spectroscopy and polarization-selective IR pump-probe experiments performed on the OD hydroxyl stretch of dilute HOD in H(2)O. The vibrational echo experiments measure spectral diffusion, and the pump-probe experiments measure orientational relaxation. Both experimental observables are directly related to the structural dynamics of water's hydrogen bond network. The measurements performed on NaBr solutions as a function of concentration show that the hydrogen bond dynamics slow as the NaBr concentration increases. The most pronounced change is in the longest time scale dynamics which are related to the global rearrangement of the hydrogen bond structure. Complete hydrogen bond network randomization slows by a factor of approximately 3 in approximately 6 M NaBr solution compared to that in bulk water. The hydrogen bond dynamics of water in nanoscopically confined environments are studied by encapsulating water molecules in ionic head group (AOT) and nonionic head group (Igepal CO 520) reverse micelles. Water dynamics in the nanopools of AOT reverse micelles are studied as a function of size by observing orientational relaxation. Orientational relaxation dynamics deviate significantly from bulk water when the size of the reverse micelles is smaller than several nm and become nonexponential and slower as the size of the reverse micelles decreases. In the smallest reverse micelles, orientational relaxation (hydrogen bond structural randomization) is almost 20 times slower than that in bulk water. To determine if the changes in dynamics from bulk water are caused by the influence of the ionic head groups of AOT or the nanoconfinement, the water dynamics in 4 nm nanopools in AOT reverse micelles (ionic) and Igepal reverse micelles (nonionic) are compared. It is found that the water orientational relaxation in the 4 nm diameter nanopools of the two types of reverse micelles is almost identical, which indicates that confinement by an interface to form a nanoscopic water pool is a primary factor governing the dynamics of nanoscopic water rather than the presence of charged groups at the interface.     

Dynamics of nanoscopic water: vibrational echo and infrared pump-probe studies of reverse micelles

The dynamics of water in nanoscopic pools 1.7-4.0 nm in diameter in AOT reverse micelles were studied with ultrafast infrared spectrally resolved stimulated vibrational echo and pump-probe spectroscopies. The experiments were conducted on the OD hydroxyl stretch of low-concentration HOD in the H2O, providing a direct examination of the hydrogen-bond network dynamics. Pump-probe experiments show that the vibrational lifetime of the OD stretch mode increases as the size of the reverse micelle decreases. These experiments are also sensitive to hydrogen-bond dissociation and reformation dynamics, which are observed to change with reverse micelle size. Spectrally resolved vibrational echo data were obtained at several frequencies. The vibrational echo data are compared to data taken on bulk water and on a 6 M NaCl solution, which is used to examine the role of ionic strength on the water dynamics in reverse micelles. Two types of vibrational echo measurements are presented: the vibrational echo decays and the vibrational echo peak shifts. As the water nanopool size decreases, the vibrational echo decays become slower. Even the largest nanopool (4 nm, approximately 1000 water molecules) has dynamics that are substantially slower than bulk water. It is demonstrated that the slow dynamics in the reverse micelle water nanopools are a result of confinement rather than ionic strength. The data are fit using time-dependent diagrammatic perturbation theory to obtain the frequency-frequency correlation function (FFCF) for each reverse micelle. The results are compared to the FFCF of water and show that the largest differences are in the slowest time scale dynamics. In bulk water, the slowest time scale dynamics are caused by hydrogen-bond network equilibration, i.e., the making and breaking of hydrogen bonds. For the smallest nanopools, the longest time scale component of the water dynamics is approximately 10 times longer than the dynamics in bulk water. The vibrational echo data for the smallest reverse micelle displays a dependence on the detection wavelength, which may indicate that multiple ensembles of water molecules are being observed.     

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.     

二羰基茂铁二聚体[CpFe(CO)2]2的中红外泵浦探测光谱

利用一维稳态红外光谱和5-μm泵浦探测红外光谱手段,结合量子化学计算,以非桥连三价羰基为探针,研究了二羰基茂铁二聚体[CpFe(CO)2]2在二氯甲烷中的结构和振动动力学.结果表明,[CpFe(CO)2]2两个主要结构(顺式cis和反式trans摩尔比为1.7)的振动态寿命和转动动力学都有一定不同.两种结构的两个羰基振动激发态的指数衰减过程都有一个<1ps的快组分和一个~20ps的慢组分.我们认为前者与宽带激发所产生的振动相干态的快速失相过程有关,而后者属于典型的C≡O伸缩振动态寿命.此外,cis结构与溶剂的较强作用使得其转动衰减较慢.结果表明,非桥连羰基的红外吸收频率和振转动力学对分子结构和溶剂环境都非常敏感.     

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.     

Protein dynamics in cytochrome P450 molecular recognition and substrate specificity using 2D IR vibrational echo spectroscopy

Cytochrome (cyt) P450s hydroxylate a variety of substrates that can differ widely in their chemical structure. The importance of these enzymes in drug metabolism and other biological processes has motivated the study of the factors that enable their activity on diverse classes of molecules. Protein dynamics have been implicated in cyt P450 substrate specificity. Here, 2D IR vibrational echo spectroscopy is employed to measure the dynamics of cyt P450(cam) from Pseudomonas putida on fast time scales using CO bound at the active site as a vibrational probe. The substrate-free enzyme and the enzyme bound to both its natural substrate, camphor, and a series of related substrates are investigated to explicate the role of dynamics in molecular recognition in cyt P450(cam) and to delineate how the motions may contribute to hydroxylation specificity. In substrate-free cyt P450(cam), three conformational states are populated, and the structural fluctuations within a conformational state are relatively slow. Substrate binding selectively stabilizes one conformational state, and the dynamics become faster. Correlations in the observed dynamics with the specificity of hydroxylation of the substrates, the binding affinity, and the substrates' molecular volume suggest that motions on the hundreds of picosecond time scale contribute to the variation in activity of cyt P450(cam) toward different substrates.     

Ultrafast vibrational relaxation in photogenerated S1 α-terthiophene in solution by femtosecond time-resolved absorption/emission and picosecond time-resolved Raman spectroscopy

Ultrafast vibrational population redistribution dynamics and subsequent cooling processes have been identified and quantified for photogenerated S₁ α-terthiophene in solution by using femtosecond time-resolved absorption/emission spectroscopy combined with picosecond time-resolved Raman spectroscopy. The population redistribution among the S₁ vibrational levels dominates the dynamics for the first few picoseconds after the photoexcitation. The cooling of the solvent-solute system then takes place in the time scale of a few tens of picoseconds.     

Binary theory of dephasing in liquid solutions. I. The non-markovian theory of encounters

A non-markovian version of encounter theory is proposed which consistently takes into account fast dissipative processes. Kinetic equations of vibrational molecular states dephasing at binary encounters in solutions are derived in terms of this theory. A particular case of fast molecular rotation is considered. The characteristic impurity dephasing time is calculated for the case of weak molecular interactions.     

Photoinduced bimolecular electron transfer investigated by femtosecond time-resolved infrared spectroscopy

Ultrafast infrared transient absorption spectroscopy is used to study the photoinduced bimolecular electron transfer reaction between perylene in the first singlet excited state and 1,4-dicyanobenzene in acetonitrile and dichloromethane. Following vibrational marker modes on both donor and acceptor sides in real time provides direct insight into the structural dynamics during the reaction. A band narrowing on a time scale of a few tens of picoseconds observed on the antisymmetric CN stretching vibration of the dicyanobenzene radical anion indicates that a substantial part of the excess energy is channeled into vibrational modes of the product, despite the fact that the reaction is weakly exergonic. An additional narrowing of the same band on a time scale of several hundreds of picoseconds observed in acetonitrile only is interpreted as a signature of the dissociation of the geminate ion pairs into free ions.