Fast protein dynamics probed with infrared vibrational echo experiments

Collective oscillations in the classical nonlinear response of a chaotic system

We establish a general semiquantitative phase-space picture of the classical nonlinear response in a strongly chaotic system. As opposed to the case of stable dynamics, the response functions decay exponentially at long times. Damped oscillations in response functions are attributed to collective resonances which do not correspond to any periodic classical motions. We calculate analytically the second-order response in a simple chaotic system and demonstrate the relevance of the concept for the interpretation of spectroscopic data.

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.

Investigations of vibrational coherence in the low-frequency region of ferric heme proteins

Femtosecond coherence spectroscopy is applied to a series of ferric heme protein samples. The low-frequency vibrational spectra that are revealed show dominant oscillations near 40 cm(-1). MbCN is taken as a typical example of a histidine-ligated, six-coordinate, ferric heme and a comprehensive spectroscopic analysis is carried out. The results of this analysis reveal a new heme photoproduct species, absorbing near 418 nm, which is consistent with the photolysis of the His(93) axial ligand. The photoproduct undergoes subsequent rebinding/recovery with a time constant of approximately 4 ps. The photoproduct lineshapes are consistent with a photolysis quantum yield of 75-100%, although the observation of a relatively strong six-coordinate heme coherence near 252 cm(-1) (assigned to nu(9) in the MbCN Raman spectrum) suggests that the 75% lower limit is much more likely. The phase and amplitude excitation profiles of the low-frequency mode at 40 cm(-1) suggest that this mode is strongly coupled to the MbCN photoproduct species and it is assigned to the doming mode of the transient penta-coordinated material. The absolute phase of the 40 cm(-1) mode is found to be pi/2 on the red side of 418 nm and it jumps to 3pi/2 as excitation is tuned to the blue side of 418 nm. The absolute phase of the 40 cm(-1) signal is not explained by the standard theory for resonant impulsive stimulated Raman scattering. New mechanisms that give a dominant momentum impulse to the resonant wavepacket, rather than a coordinate displacement, are discussed. The possibilities of heme iron atom recoil after photolysis, as well as ultrafast nonradiative decay, are explored as potential ways to generate the strong momentum impulse needed to understand the phase properties of the 40 cm(-1) mode.

Two-dimensional infrared spectra reveal relaxation of the nonnucleoside inhibitor TMC278 complexed with HIV-1 reverse transcriptase

The two nitrile groups at the wings of the nonnucleoside HIV-1 reverse transcriptase (RT) inhibitor TMC278 are both identified in high-sensitivity 2D IR spectroscopy experiments of the HIV-1 RT/TMC278 complex. The vibrational spectra indicate that the two arms of the inhibitor sense quite different environments within the hydrophobic pocket. The vibrational relaxation of the two arms are almost equal at 3 ps from model studies. The 2D IR spectra expose a significant distribution of nitrile frequencies that diffuse at equilibrium on ultrafast time scales ranging from hundreds of femtoseconds to tens of picoseconds. The slow spectral diffusion of the cyanovinyl arm of the inhibitor is attributed to its interaction with the backbone and side chains in the hydrophobic tunnel. The results show that the inhibitor cyano modes lose memory of their structural configurations relative to the hydrophobic pocket within tens of picoseconds. The cross-peaks between the two arms of the drug are tentatively attributed to relaxation of the nitrile state with both arms excited.

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.

Quantitative vibrational dynamics of iron in carbonyl porphyrins

We use nuclear resonance vibrational spectroscopy and computational predictions based on density functional theory (DFT) to explore the vibrational dynamics of (57)Fe in porphyrins that mimic the active sites of histidine-ligated heme proteins complexed with carbon monoxide. Nuclear resonance vibrational spectroscopy yields the complete vibrational spectrum of a Mössbauer isotope, and provides a valuable probe that is not only selective for protein active sites but quantifies the mean-squared amplitude and direction of the motion of the probe nucleus, in addition to vibrational frequencies. Quantitative comparison of the experimental results with DFT calculations provides a detailed, rigorous test of the vibrational predictions, which in turn provide a reliable description of the observed vibrational features. In addition to the well-studied stretching vibration of the Fe-CO bond, vibrations involving the Fe-imidazole bond, and the Fe-N(pyr) bonds to the pyrrole nitrogens of the porphyrin contribute prominently to the observed experimental signal. All of these frequencies show structural sensitivity to the corresponding bond lengths, but previous studies have failed to identify the latter vibrations, presumably because the coupling to the electronic excitation is too small in resonance Raman measurements. We also observe the FeCO bending vibrations, which are not Raman active for these unhindered model compounds. The observed Fe amplitude is strongly inconsistent with three-body oscillator descriptions of the FeCO fragment, but agrees quantitatively with DFT predictions. Over the past decade, quantum chemical calculations have suggested revised estimates of the importance of steric distortion of the bound CO in preventing poisoning of heme proteins by carbon monoxide. Quantitative agreement with the predicted frequency, amplitude, and direction of Fe motion for the FeCO bending vibrations provides direct experimental support for the quantum chemical description of the energetics of the FeCO unit.

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.

Classical and quantum mechanical infrared echoes from resonantly coupled molecular vibrations

The nonlinear response function associated with the infrared vibrational echo is calculated for a quantum mechanical model of resonantly coupled, anharmonic oscillators at zero temperature. The classical mechanical response function is determined from the quantum response function by setting variant Planck's over 2pi-->0, permitting the comparison of the effects of resonant vibrational coupling among an arbitrary number of anharmonic oscillators on quantum and classical vibrational echoes. The quantum response function displays a time dependence that reflects both anharmonicity and resonant coupling, while the classical response function depends on anharmonicity only through a time-independent amplitude, and shows a time dependence controlled only by the resonant coupling. In addition, the classical response function grows without bound in time, a phenomenon associated with the nonlinearity of classical mechanics, and absent in quantum mechanics. This unbounded growth was previously identified in the response function for a system without resonant vibrational energy transfer, and is observed to persist in the presence of resonant coupling among vibrations. Quantitative agreement between classical and quantum response functions is limited to a time scale of duration inversely proportional to the anharmonicity.

Fifth-order contributions to ultrafast spectrally resolved vibrational echoes: heme-CO proteins

The fifth order contributions to the signals of ultrafast infrared spectrally resolved stimulated vibrational echoes at high intensities have been investigated in carbonmonoxy heme proteins. High intensities are often required to obtain good data. Intensity dependent measurements are presented on hemoglobin-CO (Hb-CO) and a mutant of myoglobin, H64V-CO. The spectrally resolved vibrational echoes demonstrate that fifth order effects arise at both the 1-0 and the 2-1 emission frequencies of the stretching mode of the CO chromophore bound at the active site of heme proteins. Unlike one-dimensional experiments, in which the signal is integrated over all emission frequencies, spectrally resolving the signal shows that the fifth order contributions have a much more pronounced influence on the 2-1 transition than on the 1-0 transition. By spectrally isolating the 1-0 transition, the influence of fifth order contributions to vibrational echo data can be substantially reduced. Analysis of fifth order Feynman diagrams that contribute in the vibrational echo phase-matched direction demonstrates the reason for the greater influence of fifth order processes on the 1-2 transition, and that the fifth order contributions are heterodyne amplified by the third order signal. Finally, it is shown that the anharmonic oscillations in vibrational echo data of Hb-CO that previous work had attributed strictly to fifth order effects arise even without fifth order contributions.

A structure-based simulation approach for electron paramagnetic resonance spectra using molecular and stochastic dynamics simulations

Electron paramagnetic resonance (EPR) spectroscopy using site-directed spin-labeling is an appropriate technique to analyze the structure and dynamics of flexible protein regions as well as protein-protein interactions under native conditions. The analysis of a set of protein mutants with consecutive spin-label positions leads to the identification of secondary and tertiary structure elements. In the first place, continuous-wave EPR spectra reflect the motional freedom of the spin-label specifically linked to a desired site within the protein. EPR spectra calculations based on molecular dynamics (MD) and stochastic dynamics simulations facilitate verification or refinement of predicted computer-aided models of local protein conformations. The presented spectra simulation algorithm implies a specialized in vacuo MD simulation at 600 K with additional restrictions to sample the entire accessible space of the bound spin-label without large temporal effort. It is shown that the distribution of spin-label orientations obtained from such MD simulations at 600 K agrees well with the extrapolated motion behavior during a long timescale MD at 300 K with explicit water. The following potential-dependent stochastic dynamics simulation combines the MD data about the site-specific orientation probabilities of the spin-label with a realistic rotational diffusion coefficient yielding a set of trajectories, each more than 700 ns long, essential to calculate the EPR spectrum. Analyses of a structural model of the loop between helices E and F of bacteriorhodopsin are illustrated to demonstrate the applicability and potentials of the reported simulation approach. Furthermore, effects on the motional freedom of bound spin-labels induced by solubilization of bacteriorhodopsin with Triton X-100 are examined.

Infrared absorption study of the heme pocket dynamics of carbonmonoxyheme proteins

The temperature dependencies of the infrared absorption CO bands of carboxy complexes of horseradish peroxidase (HRP(CO)) in glycerol/water mixture at pH 6.0 and 9.3 are interpreted using the theory of optical absorption bandshape. The bands' anharmonic behavior is explained assuming that there is a higher-energy set of conformational substates (CSS(h)), which are populated upon heating and correspond to the protein substates with disordered water molecules in the heme pocket. Analysis of the second moments of the CO bands of the carboxy complexes of myoglobin (Mb(CO)) and hemoglobin (Hb(CO)), and of HRP(CO) with benzohydroxamic acid (HRP(CO)+BHA), shows that the low energy CSS(h) exists also in the open conformation of Mb(CO), where the heme pocket is spacious enough to accommodate a water molecule. In the HRP(CO)+BHA and closed conformations of Mb(CO) and Hb(CO), the heme pocket is packed with BHA and different amino acids, the CSS(h) has much higher energy and is hardly populated even at the highest temperatures. Therefore only motions of these amino acids contribute to the band broadening. These motions are linked to the protein surface and frozen in the glassy matrix, whereas in the liquid solvent they are harmonic. Thus the second moment of the CO band is temperature-independent in glass and is proportional to the temperature in liquid. The temperature dependence of the second moment of the CO peak of HRP(CO) in the trehalose glass exhibits linear coupling to an oscillator. This oscillator can be a moving water molecule locked in the heme pocket in the whole interval of temperatures or a trehalose molecule located in the heme pocket.

Modeling vibrational dephasing and energy relaxation of intramolecular anharmonic modes for multidimensional infrared spectroscopies

Starting from a system-bath Hamiltonian in a molecular coordinate representation, we examine an applicability of a stochastic multilevel model for vibrational dephasing and energy relaxation in multidimensional infrared spectroscopy. We consider an intramolecular anharmonic mode nonlinearly coupled to a colored noise bath at finite temperature. The system-bath interaction is assumed linear plus square in the system coordinate, but linear in the bath coordinates. The square-linear system-bath interaction leads to dephasing due to the frequency fluctuation of system vibration, while the linear-linear interaction contributes to energy relaxation and a part of dephasing arises from anharmonicity. To clarify the role and origin of vibrational dephasing and energy relaxation in the stochastic model, the system part is then transformed into an energy eigenstate representation without using the rotating wave approximation. Two-dimensional (2D) infrared spectra are then calculated by solving a low-temperature corrected quantum Fokker-Planck (LTC-QFP) equation for a colored noise bath and by the stochastic theory. In motional narrowing regime, the spectra from the stochastic model are quite different from those from the LTC-QFP. In spectral diffusion regime, however, the 2D line shapes from the stochastic model resemble those from the LTC-QFP besides the blueshifts caused by the dissipation from the colored noise bath. The preconditions for validity of the stochastic theory for molecular vibrational motion are also discussed.

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 proteins encapsulated in silica sol-gel glasses studied with IR vibrational echo spectroscopy

Spectrally resolved infrared stimulated vibrational echo spectroscopy is used to measure the fast dynamics of heme-bound CO in carbonmonoxy-myoglobin (MbCO) and -hemoglobin (HbCO) embedded in silica sol-gel glasses. On the time scale of approximately 100 fs to several picoseconds, the vibrational dephasing of the heme-bound CO is measurably slower for both MbCO and HbCO relative to that of aqueous protein solutions. The fast structural dynamics of MbCO, as sensed by the heme-bound CO, are influenced more by the sol-gel environment than those of HbCO. Longer time scale structural dynamics (tens of picoseconds), as measured by the extent of spectral diffusion, are the same for both proteins encapsulated in sol-gel glasses compared to that in aqueous solutions. A comparison of the sol-gel experimental results to viscosity-dependent vibrational echo data taken on various mixtures of water and fructose shows that the sol-gel-encapsulated MbCO exhibits dynamics that are the equivalent of the protein in a solution that is nearly 20 times more viscous than bulk water. In contrast, the HbCO dephasing in the sol-gel reflects only a 2-fold increase in viscosity. Attempts to alter the encapsulating pore size by varying the molar ratio of silane precursor to water (R value) used to prepare the sol-gel glasses were found to have no effect on the fast or steady-state spectroscopic results. The vibrational echo data are discussed in the context of solvent confinement and protein-pore wall interactions to provide insights into the influence of a confined environment on the fast structural dynamics experienced by a biomolecule.

Electronic and Vibrational Second-Order Nonlinear Optical Properties of Protein Secondary Structural Motifs

A perturbation theory approach was developed for predicting the vibrational and electronic second-order nonlinear optical (NLO) polarizabilities of materials and macromolecules comprised of many coupled chromophores, with an emphasis on common protein secondary structural motifs. The polarization-dependent NLO properties of electronic and vibrational transitions in assemblies of amide chromophores comprising the polypeptide backbones of proteins were found to be accurately recovered in quantum chemical calculations by treating the coupling between adjacent oscillators perturbatively. A novel diagrammatic approach was developed to provide an intuitive visual means of interpreting the results of the perturbation theory calculations. Using this approach, the chiral and achiral polarization-dependent electronic SHG, isotropic SFG, and vibrational SFG nonlinear optical activities of protein structures were predicted and interpreted within the context of simple orientational models.

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.

Effect of protein binding on ultrafast DNA dynamics: characterization of a DNA:APE1 complex

Synthetic oligonucleotides with a fluorescent coumarin group replacing a basepair have been used in recent time-resolved Stokes-shift experiments to measure DNA dynamics on the femtosecond to nanosecond timescales. Here, we show that the APE1 endonuclease cleaves such a modified oligonucleotide at the abasic site opposite the coumarin with only a fourfold reduction in rate. In addition, a noncatalytic mutant (D210N) binds tightly to the same oligonucleotide, albeit with an 85-fold reduction in binding constant relative to a native oligonucleotide containing a guanine opposite the abasic site. Thus, the modified oligonucleotide retains substantial biological activity and serves as a useful model of native DNA. In the complex of the coumarin-containing oligonucleotide and the noncatalytic APE1, the dye's absorption spectrum is shifted relative to its spectrum in either water or within the unbound oligonucleotide. Thus the dye occupies a site within the DNA:protein complex. This result is consistent with modeling, which shows that the complex accommodates coumarin at the site of the orphaned base with little distortion of the native structure. Stokes-shift measurements of the complex show surprisingly little change in the dynamics within the 40 ps-40 ns time range.

Exciton Analysis in 2D Electronic Spectroscopy

A theoretical description of femtosecond two-dimensional electronic spectroscopy of multichromophoric systems is presented. Applying the stationary phase approximation to the calculation of photon echo spectra and taking into account exciton relaxation processes, we obtain an analytic expression for numerical simulations of time- and frequency-resolved 2D photon echo signals. The delocalization of one-exciton states, spatial overlaps between the probability densities of different excitonic states, and their influences on both one- and two-dimensional electronic spectra are studied. The nature of the off-diagonal cross-peaks and the time evolution of both diagonal and off-diagonal peak amplitudes are discussed in detail by comparing experimentally measured and theoretically simulated 2D spectra of the natural Fenna-Matthews-Olson (FMO) photosynthetic light-harvesting complex. We find that there are two noncascading exciton energy relaxation pathways.

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.

Vibrational energy relaxation in proteins

An overview of theories related to vibrational energy relaxation (VER) in proteins is presented. VER of a selected mode in cytochrome c is studied by using two theoretical approaches. One approach is the equilibrium simulation approach with quantum correction factors, and the other is the reduced model approach, which describes the protein as an ensemble of normal modes interacting through nonlinear coupling elements. Both methods result in similar estimates of the VER time (subpicoseconds) for a CD stretching mode in the protein at room temperature. The theoretical predictions are in accord with previous experimental data. A perspective on directions for the detailed study of time scales and mechanisms of VER in proteins is presented.

Ultrafast molecular dynamics of liquid aromatic molecules and the mixtures with CCl4

The ultrafast molecular dynamics of liquid aromatic molecules, benzene, toluene, ethylbenzene, cumene, and 1,3-diphenylpropane, and the mixtures with CCl(4) have been investigated by means of femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy. The picosecond Kerr transients of benzene, toluene, ethylbenzene, and cumene and the mixtures with CCl(4) show a biexponential feature. 1,3-Diphenylpropane and the mixtures with CCl(4) show triexponential picosecond Kerr transients. The slow relaxation time constants of the aromatic molecules and the mixtures with CCl(4) are qualitatively described by the Stoke-Einstein-Debye hydrodynamic model. The ultrafast dynamics have been discussed based on the Kerr spectra in the frequency range of 0-800 cm(-1) obtained by the Fourier transform analysis of the Kerr transients. The line shapes of the low-frequency intermolecular spectra located at 0-180 cm(-1) frequency range have been analyzed by two Brownian oscillators ( approximately 11 cm(-1) and approximately 45 cm(-1) peaks) and an antisymmetric Gaussian function ( approximately 65 cm(-1) peak). The spectrum shape of 1,3-diphenylpropane is quite different from the spectrum shapes of the other aromatic molecules for the low magnitude of the low-frequency mode of 1,3-diphenylpropane and/or an intramolecular vibration. Although the concentration dependences of the low- and intermediate-frequency intermolecular modes (Brownian oscillators) do not show a significant trend, the width of high-frequency intermolecular mode (antisymmetric Gaussian) becomes narrower with the higher CCl(4) concentration for all the aromatics mixtures with CCl(4). The result indicates that the inhomogeneity of the intermolecular vibrational mode in aromatics/CCl(4) mixtures is decreasing with the lower concentration of aromatics. The intramolecular vibrational modes of the aromatic molecules observed in the Kerr spectra are also shown with the calculation results based on the density functional theory.

Probing electric fields in protein cavities by using the vibrational stark effect of carbon monoxide

To determine the magnitude and direction of the internal electric field in the Xe4 cavity of myoglobin mutant L29W-S108L, we have studied the vibrational Stark effect of carbon monoxide (CO) using infrared spectroscopy at cryogenic temperatures. CO was photodissociated from the heme iron and deposited selectively in Xe4. Its infrared spectrum exhibits Stark splitting into two bands associated with CO in opposite orientations. Two different photoproduct states can be distinguished, C' and C'', with markedly different properties. For C', characteristic temperature-dependent changes of the area, shift, and width were analyzed, based on a dynamic model in which the CO performs fast librations within a double-well model potential. For the barrier between the wells, a height of approximately 1.8 kJ/mol was obtained, in which the CO performs oscillations at an angular frequency of approximately 25 cm(-1). The magnitude of the electric field in the C' conformation was determined as 11.1 MV/cm; it is tilted by an angle of 29 degrees to the symmetry axis of the potential. Above 140 K, a protein relaxation leads to a significantly altered photoproduct, C'', with a smaller Stark splitting and a more confining potential (barrier >4 kJ/mol) governing the CO librations.

Interpreting nonlinear vibrational spectroscopy with the classical mechanical analogs of double-sided Feynman diagrams

Observables in coherent, multiple-pulse infrared spectroscopy may be computed from a vibrational nonlinear response function. This response function is conventionally calculated quantum-mechanically, but the challenges in applying quantum mechanics to large, anharmonic systems motivate the examination of classical mechanical vibrational nonlinear response functions. We present an approximate formulation of the classical mechanical third-order vibrational response function for an anharmonic solute oscillator interacting with a harmonic solvent, which establishes a clear connection between classical and quantum mechanical treatments. This formalism permits the identification of the classical mechanical analog of the pure dephasing of a quantum mechanical degree of freedom, and suggests the construction of classical mechanical analogs of the double-sided Feynman diagrams of quantum mechanics, which are widely applied to nonlinear spectroscopy. Application of a rotating wave approximation permits the analytic extraction of signals obeying particular spatial phase matching conditions from a classical-mechanical response function. Calculations of the third-order response function for an anharmonic oscillator coupled to a harmonic solvent are compared to numerically correct classical mechanical results.

Semiclassical calculation of the vibrational echo

The infrared echo measurement probes the time scales of the molecular motions that couple to a vibrational transition. Computation of the echo observable within rigorous quantum mechanics is problematic for systems with many degrees of freedom, motivating the development of semiclassical approximations to the nonlinear optical response. We present a semiclassical approximation to the echo observable, based on the Herman-Kluk propagator. This calculation requires averaging over a quantity generated by two pairs of classical trajectories and associated stability matrices, connected by a pair of phase-space jumps. Quantum, classical, and semiclassical echo calculations are compared for a thermal ensemble of noninteracting anharmonic oscillators. The semiclassical approach uses input from classical mechanics to reproduce the significant features of a complete, quantum mechanical calculation of the nonlinear response.

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.

Are the conformational dynamics and the ligand binding properties of myoglobin affected by exposure to microwave radiation?

The global uptake of mobile communication emphasizes the question about possible adverse consequences of the exposure to low-level radiofrequency radiation from mobile phones on human health as result of so-called "non-thermal effects". In order to state safety guidelines it seems appropriate to start by excluding, if possible, non-specific effects on structural and dynamic properties of fundamental biomolecules such as proteins. Proteins are flexible polyelectrolytes; thus, they are susceptible, in principle, to the action of electromagnetic fields. In this article, we investigated the effects of microwaves on structural and functional properties of Tunnus tynnus myoglobin at 1.95 GHz, a frequency used by new wireless microwave communication systems. The protein solution was exposed for 2.5 h to 51 mW/g SAR (specific absorption rate) level. Measurements of absorption spectroscopy, circular dichroism and fluorescence emission decay in the frequency domain do not exhibit any influence of the radiation on the native structural state of protein macromolecules.

Charge-transfer mechanism for cytochrome c adsorbed on nanometer thick films. Distinguishing frictional control from conformational gating

Using nanometer thick tunneling barriers with specifically attached cytochrome c, the electron-transfer rate constant was studied as a function of the SAM composition (alkane versus terthiophene), the omega-terminating group type (pyridine, imidazole, nitrile), and the solution viscosity. At large electrode-reactant separations, the pyridine terminated alkanethiols exhibit an exponential decline of the rate constant with increasing electron-transfer distance. At short separations, a plateau behavior, analogous to systems involving -COOH terminal groups to which cytochrome c can be attached electrostatically, is observed. The dependence of the rate constant in the plateau region on system properties is investigated. The rate constant is insensitive to the mode of attachment to the surface but displays a significant viscosity dependence, change with spacer composition (alkane versus terthiophene), and nature of the solvent (H(2)O versus D(2)O). Based on these findings and others, the conclusion is drawn that the charge-transfer rate constant at short distance is determined by polarization relaxation processes in the structure, rather than the electron tunneling probability or large-amplitude conformational rearrangement (gating). The transition in reaction mechanism with distance reflects a gradual transition between the tunneling and frictional mechanisms. This conclusion is consistent with data from a number of other sources as well.

Flexibility and molecular recognition in the immune system

Photon echo spectroscopy has been used to measure the response of three antibody-binding sites to perturbation from electronic excitation of a bound antigen, fluorescein. The three antibodies show motions that range in time scale from tens of femtoseconds to nanoseconds. Relative to the others, one antibody, 4-4-20, possesses a rigid binding site that likely results from a short and inflexible heavy chain complementarity-determining region 3 (HCDR3) loop and a critical Tyr that acts as a "molecular splint," rigidifying the antigen across its most flexible internal degree of freedom. The remaining two antibodies, 34F10 and 40G4, despite being generated against the same antigen, possess binding sites that are considerably more flexible. The more flexible combining sites likely result from longer HCDR3 loops and a deletion in the light chain complementarity-determining region 1 (LCDR1) that removes the critical Tyr residue. The binding site flexibilities may result in varying mechanisms of antigen recognition including lock-and-key, induced-fit, and conformational selection.

Effect of solution viscosity on intramolecular electron transfer in sulfite oxidase

Our previous studies have shown that the rate constant for intramolecular electron transfer (IET) between the heme and molybdenum centers of chicken liver sulfite oxidase varies from approximately 20 to 1400 s(-1) depending upon reaction conditions [Pacheco, A., Hazzard, J. T., Tollin, G., and Enemark, J. H. (1999) J. Biol. Inorg. Chem. 4, 390-401]. These two centers are linked by a flexible polypeptide loop, suggesting that conformational changes, which alter the Mo-Fe distance, may play an important role in the observed IET rates. In this study, we have investigated IET in sulfite oxidase using laser flash photolysis as a function of solution viscosity. The solution viscosity was varied over the range of 1.0-2.0 cP by addition of either polyethylene glycol 400 or sucrose. In the presence of either viscosogen, an appreciable decrease in the IET rate constant value is observed with an increase in the solvent viscosity. The IET rate constant exhibits a linear dependence on the negative 0.7th power of the viscosity. Steady-state kinetics and EPR experiments are consistent with the interpretation that viscosity, and not other properties of the added viscosogens, is responsible for the dependence of IET rates on the solvent composition. The results are consistent with the role of conformational changes on IET in sulfite oxidase, which helps to clarify the inconsistency between the large rate constant for IET between the Mo and Fe centers and the long distance (approximately 32 A) between these two metal centers observed in the crystal structure [Kisker, C., Schindelin, H., Pacheco, A., Wehbi, W., Garnett, R. M., Rajagopalan, K. V., Enemark, J. H., and Rees, D. C. (1997) Cell 91, 973-983].