Investigations of Coherent Vibrational Oscillations in Myoglobin

Tracking Energy Transfer across a Platinum Center

Rigid, conjugated alkyne bridges serve as important components in various transition-metal complexes used for energy conversion, charge separation, sensing, and molecular electronics. Alkyne stretching modes have potential for modulating charge separation in donor–bridge–acceptor compounds. Understanding the rules of energy relaxation and energy transfer across the metal center in such compounds can help optimize their electron transfer switching properties. We used relaxation-assisted two-dimensional infrared spectroscopy to track energy transfer across metal centers in platinum complexes featuring a triazole-terminated alkyne ligand of two or six carbons, a perfluorophenyl ligand, and two tri(p-tolyl)phosphine ligands. Comprehensive analyses of waiting-time dynamics for numerous cross and diagonal peaks were performed, focusing on coherent oscillation, energy transfer, and cooling parameters. These observables augmented with density functional theory computations of vibrational frequencies and anharmonic force constants enabled identification of different functional groups of the compounds. Computations of vibrational relaxation pathways and mode couplings were performed, and two regimes of intramolecular energy redistribution are described. One involves energy transfer between ligands via high-frequency modes; the transfer is efficient only if the modes involved are delocalized over both ligands. The energy transport pathways between the ligands are identified. Another regime involves redistribution via low-frequency delocalized modes, which does not lead to interligand energy transport.

Excited-State Energy Surfaces in Molecules Revealed by Impulsive Stimulated Raman Excitation Profiles

Photophysical and photochemical processes are ruled by the interplay between transient vibrational and electronic degrees of freedom, which are ultimately determined by the multidimensional potential energy surfaces (PESs). Differences between ground and excited PESs are encoded in the relative intensities of resonant Raman bands, but they are experimentally challenging to access, requiring measurements at multiple wavelengths under identical conditions. Herein, we perform a two-color impulsive vibrational scattering experiment to launch nuclear wavepacket motions by an impulsive pump and record their coupling with a targeted excited-state potential by resonant Raman processes with a delayed probe, generating in a single measurement background-free vibrational spectra across the entire sample absorption. Building on the interference between the multiple pathways resonant with the excited-state manifold that generate the Raman signal, we show how to experimentally tune their relative phase by varying the probe chirp, decoding nuclear displacements along different normal modes and revealing the multidimensional PESs. Our results are validated against time-dependent density functional theory.

Low-Frequency Protein Motions Coupled to Catalytic Sites

This review examines low-frequency vibrational modes of proteins and their coupling to enzyme catalytic sites. That protein motions are critical to enzyme function is clear, but the kinds of motions present in proteins and how they are involved in function remain unclear. Several models of enzyme-catalyzed reaction suggest that protein dynamics may be involved in the chemical step of the catalyzed reaction, but the evidence in support of such models is indirect. Spectroscopic studies of low-frequency protein vibrations consistently show that there are underdamped modes of the protein with frequencies in the tens of wavenumbers where overdamped behavior would be expected. Recent studies even show that such underdamped vibrations modulate enzyme active sites. These observations suggest that increasingly sophisticated spectroscopic methods will be able to unravel the link between low-frequency protein vibrations and enzyme function.

Advances and opportunities in ultrafast X-ray crystallography and ultrafast structural optical crystallography of nuclear and electronic protein dynamics

Both nuclear and electronic dynamics contribute to protein function and need multiple and complementary techniques to reveal their ultrafast structural dynamics response. Real-space information obtained from the measurement of electron density dynamics by X-ray crystallography provides aspects of both, while the molecular physics of coherence parameters and frequency-frequency correlation needs spectroscopy methods. Ultrafast pump-probe applications of protein dynamics in crystals provide real-space information through direct X-ray crystallographic structure analysis or through structural optical crystallographic analysis. A discussion of methods of analysis using ultrafast macromolecular X-ray crystallography and ultrafast nonlinear structural optical crystallography is presented. The current and future high repetition rate capabilities provided by X-ray free electron lasers for ultrafast diffraction studies provide opportunities for optical control and optical selection of nuclear coherence which may develop to access higher frequency dynamics through improvements of sensitivity and time resolution to reveal coherence directly. Specific selection of electronic coherence requires optical probes, which can provide real-space structural information through photoselection of oriented samples and specifically in birefringent crystals. Ultrafast structural optical crystallography of photosynthetic energy transfer has been demonstrated, and the theory of two-dimensional structural optical crystallography has shown a method for accessing the structural selection of electronic coherence.

Observation of Vibronic-Coupling-Mediated Energy Transfer in Light-Harvesting Nanotubes Stabilized in a Solid-State Matrix

Ultrafast vibrational spectroscopy is employed to obtain real-time structural information on energy transport in double-walled light-harvesting nanotubes at room temperature, stabilized in a host matrix to mimic the rigid scaffolds of natural light-harvesting systems. We observe evidence of a low-frequency vibrational mode at 315 cm^-1, which transfers excitons from the outer wall of the nanotubes to a crossing point through which energy transfer to the inner wall can occur. This mode is furthermore absent in solution phase. Importantly, the coherence of this mode is not transferred to the inner wall upon energy transfer and is only present on the outer wall's excited-state energy surface, highlighting that complete energy transfer between the outer and inner walls does not take place. Isolation of the individual walls of the nanotubes provides evidence that this mode corresponds to a supramolecular motion of the nanotubes. Our results emphasize the importance of the solid-state environment in modulating vibronic coupling and directing energy transfer in molecular light-harvesting systems.

Vibrational coherence transfer in the ultrafast intersystem crossing of a diplatinum complex in solution

We investigate the ultrafast transient absorption response of tetrakis(μ-pyrophosphito)diplatinate(II), [Pt₂(μ-P₂O₅H₂)₄]^4- [hereafter abbreviated Pt(pop)], in acetonitrile upon excitation of its lowest singlet ¹A_2u state. Compared with previously reported solvents [van der Veen RM, Cannizzo A, van Mourik F, Vlček A, Jr, Chergui M (2011) J Am Chem Soc 133:305-315], a significant shortening of the intersystem crossing (ISC) time (<1 ps) from the lowest singlet to the lowest triplet state is found, allowing for a transfer of vibrational coherence, observed in the course of an ISC in a polyatomic molecule in solution. Density functional theory (DFT) quantum mechanical/molecular mechanical (QM/MM) simulations of Pt(pop) in acetonitrile and ethanol show that high-lying, mostly triplet, states are strongly mixed and shifted to lower energies due to interactions with the solvent, providing an intermediate state (or manifold of states) for the ISC. This suggests that the larger the solvation energies of the intermediate state(s), the shorter the ISC time. Because the latter is smaller than the pure dephasing time of the vibrational wave packet, coherence is conserved during the spin transition. These results underscore the crucial role of the solvent in directing pathways of intramolecular energy flow.

Observing ground state vibrational coherence and excited state relaxation dynamics of a cyanine dye in pure solvents

Using a degenerate pump probe technique at 800 nm, Ground State Vibrational Coherence (GSVC) of a cyanine dye (IR780) is explored in various solvents.

Ultrafast CO Kinetics in Heme Proteins: Adiabatic Ligand Binding and Heavy Atom Tunneling

The ultrafast kinetics of CO rebinding to carbon monoxide oxidation activator protein (ChCooA) are measured over a wide temperature range and compared with the kinetics of CO binding in other heme systems such as myoglobin (Mb) and hemoglobin (Hb). The Arrhenius prefactor for CO binding to ChCooA and protoheme (∼10¹¹ s^-1) is similar to what is found for spin-allowed NO binding to heme proteins and is several orders of magnitude larger than the prefactor of Mb and Hb (∼10⁹ s^-1). This indicates that the CO binding reaction is adiabatic, in contrast to the commonly held view that it is nonadiabatic due to spin-forbidden (ΔS = 2) selection rules. Under the adiabatic condition, entropic factors, rather than spin-selection rules, are the source of the reduced Arrhenius prefactors associated with CO binding in Mb and Hb. The kinetic response of ChCooA-CO is nonexponential at all temperatures, including 298 K, and is described quantitatively using a distribution of enthalpic rebinding barriers associated with heterogeneity in the heme doming conformation. Above the solvent glass transition (T_g ∼ 180 K), the rebinding progress slows as temperature increases, and this is ascribed to an evolution of the distribution toward increased heme doming and larger enthalpic barriers. Between T_g and ∼60 K, the nonexponential rebinding slows down as the temperature is lowered and the survival fraction follows the predictions expected for a quenched barrier distribution. Below ∼60 K the rebinding kinetics do not follow these predictions unless quantum mechanical tunneling along the heme doming coordinate is also included as an active channel for CO binding.

Ultrafast Time-Resolved Emission and Absorption Spectra of meso-Pyridyl Porphyrins upon Soret Band Excitation Studied by Fluorescence Up-Conversion and Transient Absorption Spectroscopy

A comprehensive study of ultrafast molecular relaxation processes of isomeric meso-(pyridyl) porphyrins (TpyPs) has been carried out by using femtosecond time-resolved emission and absorption spectroscopic techniques upon pumping at 400 nm, Soret band (B band or S2), in 4:1 dichloromethane (DCM) and tetrahydrofuran (THF) solvent mixture. By combined studies of fluorescence up-conversion, time-correlated single photon counting, and transient absorption spectroscopic techniques, a complete model with different microscopic rate constants associated with elementary processes involved in electronic manifolds has been reported. Besides, a distinct coherent nuclear wave packet motion in Qy state is observed at low-frequency mode, ca. 26 cm(-1) region. Fluorescence up-conversion studies constitute ultrafast time-resolved emission spectra (TRES) over the whole emission range (430-710 nm) starting from S2 state to Qx state via Qy state. Careful analysis of time profiles of up-converted signals at different emission wavelengths helps to reveal detail molecular dynamics. The observed lifetimes are as indicated: A very fast decay component with 80 ± 20 fs observed at ∼435 nm is assigned to the lifetime of S2 (B) state, whereas being a rise component in the region of between 550 and 710 nm emission wavelength pertaining to Qy and Qx states, it is attributed to very fast internal conversion (IC) occurring from B → Qy and B → Qx as well. Two distinct components of Qy emission decay with ∼200-300 fs and ∼1-1.5 ps time constants are due to intramolecular vibrational redistribution (IVR) induced by solute-solvent inelastic collisions and vibrational redistribution induced by solute-solvent elastic collision, respectively. The weighted average of these two decay components is assigned as the characteristic lifetime of Qy, and it ranges between 0.3 and 0.5 ps. An additional ∼20 ± 2 ps rise component is observed in Qx emission, and it is assigned to the formation time of thermally equilibrated Qx state by vibrational cooling/relaxations of excess energy within solvent. This relaxed Qx state decays to ground as well as triplet state by 7-8 ns time scale. The femtosecond transient absorption studies of TpyPs in three different excitations at S2 (400 nm), Qy (515 nm), and Qx (590 nm) along with extensive global and target model analysis of TA data exclusively generate the true spectra of each excited species/state with their respective lifetimes along with microscopic rate constants associated with each state. The following five exponential components with lifetime values of 65-70 fs, ∼0.3-0.5 ps, ∼20 ± 2 ps, ∼7 ± 1 ns, and 1-2 μs are observed which are associated with S2, Qy, hot Qx, thermally relaxed Qx, and lowest triplet (T1) states, respectively, when excited at S2, and four (Qy, hot Qx, thermally relaxed Qx, and lowest triplet (T1) states) and three (hot Qx, thermally relaxed Qx, and lowest triplet (T1) states) states are obtained when excited at 515 nm (Qy) and 590 nm (Qx), respectively, as expected. The TA results parallel the fluorescence up-conversion studies, and both the results not only compliment each other but also unveil the ultrafast internal conversion from S2 to Qy, S2 to Qx, and Qy to Qx for all three isomers in a similar fashion with nearly equal characteristic decay times.

Kinetic Control of O2 Reactivity in H-NOX Domains

Transient absorption, resonance Raman, and vibrational coherence spectroscopies are used to investigate the mechanisms of NO and O2 binding to WT Tt H-NOX and its P115A mutant. Vibrational coherence spectra of the oxy complexes provide clear evidence for the enhancement of an iron-histidine mode near 217 cm(-1) following photoexcitation, which indicates that O2 can be dissociated in these proteins. However, the quantum yield of O2 photolysis is low, particularly in the wild type (≲3%). Geminate recombination of O2 and NO in both of these proteins is very fast (∼1.4 × 10(11) s(-1)) and highly efficient. We show that the distal heme pocket of the H-NOX system forms an efficient trap that limits the O2 off-rate and determines the overall affinity. The distal pocket hydrogen bond, which appears to be stronger in the P115A mutant, may help retard the O2 ligand from escaping into the solvent following either photoinduced or thermal dissociation. This, along with a strengthening of the Fe-O2 bond that is correlated with the significant heme ruffing and saddling distortions, explains the unusually high O2 affinity of WT Tt H-NOX and the even higher affinity found in the P115A mutant.

Direct observation of ultrafast collective motions in CO myoglobin upon ligand dissociation

The hemoprotein myoglobin is a model system for the study of protein dynamics. We used time-resolved serial femtosecond crystallography at an x-ray free-electron laser to resolve the ultrafast structural changes in the carbonmonoxy myoglobin complex upon photolysis of the Fe-CO bond. Structural changes appear throughout the protein within 500 femtoseconds, with the C, F, and H helices moving away from the heme cofactor and the E and A helices moving toward it. These collective movements are predicted by hybrid quantum mechanics/molecular mechanics simulations. Together with the observed oscillations of residues contacting the heme, our calculations support the prediction that an immediate collective response of the protein occurs upon ligand dissociation, as a result of heme vibrational modes coupling to global modes of the protein.

Investigations of the low frequency modes of ferric cytochrome c using vibrational coherence spectroscopy

Femtosecond vibrational coherence spectroscopy is used to investigate the low frequency vibrational dynamics of the electron transfer heme protein, cytochrome c (cyt c). The vibrational coherence spectra of ferric cyt c have been measured as a function of excitation wavelength within the Soret band. Vibrational coherence spectra obtained with excitation between 412 and 421 nm display a strong mode at ∼44 cm^–1 that has been assigned to have a significant contribution from heme ruffling motion in the electronic ground state. This assignment is based partially on the presence of a large heme ruffling distortion in the normal coordinate structural decomposition (NSD) analysis of the X-ray crystal structures. When the excitation wavelength is moved into the ∼421–435 nm region, the transient absorption increases along with the relative intensity of two modes near ∼55 and 30 cm^–1. The intensity of the mode near 44 cm^–1 appears to minimize in this region and then recover (but with an opposite phase compared to the blue excitation) when the laser is tuned to 443 nm. These observations are consistent with the superposition of both ground and excited state coherence in the 421–435 nm region due to the excitation of a weak porphyrin-to-iron charge transfer (CT) state, which has a lifetime long enough to observe vibrational coherence. The mode near 55 cm^–1 is suggested to arise from ruffling in a transient CT state that has a less ruffled heme due to its iron d⁶ configuration.

Investigations of the low-frequency spectral density of cytochrome c upon equilibrium unfolding

The equilibrium unfolding process of ferric horse heart cytochrome c (cyt c), induced by guanidinium hydrochloride (GdHCl), was studied using UV-vis absorption spectroscopy, resonance Raman spectroscopy, and vibrational coherence spectroscopy (VCS). The unfolding process was successfully fit using a three-state model which included the fully folded (N) and unfolded (U) states, along with an intermediate (I) assigned to a Lys bound heme. The VCS spectra revealed for the first time several low-frequency heme modes that are sensitive to cyt c unfolding: γ(a) (~50 cm(-1)), γ(b) (~80 cm(-1)), γ(c) (~100 cm(-1)), and ν(s)(His-Fe-His) at 205 cm(-1). These out-of-plane modes have potential functional relevance and are activated by protein-induced heme distortions. The free energies for the N-I and the I-U transitions at pH 7.0 and 20 °C were found to be 4.6 kcal/M and 11.6 kcal/M, respectively. Imidazole was also introduced to replace the methionine ligand so the unfolding can be modeled as a two-state system. The intensity of the mode γ(b)~80 cm(-1) remains nearly constant during the unfolding process, while the amplitudes of the other low frequency modes track with spectral changes observed at higher frequency. This confirms that the heme deformation changes are coupled to the protein tertiary structural changes that take place upon unfolding. These studies also reveal that damping of the coherent oscillations depends sensitively on the coupling between heme and the surrounding water solvent.

Fluorescence of tryptophan in designed hairpin and Trp-cage miniproteins: measurements of fluorescence yields and calculations by quantum mechanical molecular dynamics simulations

The quantum yield of tryptophan (Trp) fluorescence was measured in 30 designed miniproteins (17 β-hairpins and 13 Trp-cage peptides), each containing a single Trp residue. Measurements were made in D(2)O and H(2)O to distinguish between fluorescence quenching mechanisms involving electron and proton transfer in the hairpin peptides, and at two temperatures to check for effects of partial unfolding of the Trp-cage peptides. The extent of folding of all the peptides also was measured by NMR. The fluorescence yields ranged from 0.01 in some of the Trp-cage peptides to 0.27 in some hairpins. Fluorescence quenching was found to occur by electron transfer from the excited indole ring of the Trp to a backbone amide group or the protonated side chain of a nearby histidine, glutamate, aspartate, tyrosine, or cysteine residue. Ionized tyrosine side chains quenched strongly by resonance energy transfer or electron transfer to the excited indole ring. Hybrid classical/quantum mechanical molecular dynamics simulations were performed by a method that optimized induced electric dipoles separately for the ground and excited states in multiple π-π* and charge-transfer (CT) excitations. Twenty 0.5 ns trajectories in the tryptophan's lowest excited singlet π-π* state were run for each peptide, beginning by projections from trajectories in the ground state. Fluorescence quenching was correlated with the availability of a CT or exciton state that was strongly coupled to the π-π* state and that matched or fell below the π-π* state in energy. The fluorescence yields predicted by summing the calculated rates of charge and energy transfer are in good accord with the measured yields.

An investigation into the feasibility of myoglobin-based single-electron transistors

Myoglobin single-electron transistors were investigated using nanometer-gap platinum electrodes fabricated by electromigration at cryogenic temperatures. Apomyoglobin (myoglobin without the heme group) was used as a reference. The results suggest single-electron transport is mediated by resonant tunneling with the electronic and vibrational levels of the heme group in a single protein. They also represent a proof-of-principle that proteins with redox centers across nanometer-gap electrodes can be utilized to fabricate single-electron transistors. The protein orientation and conformation may significantly affect the conductance of these devices. Future improvements in device reproducibility and yield will require control of these factors.

Vibrational coherence spectroscopy of the heme domain in the CO-sensing transcriptional activator CooA

Femtosecond vibrational coherence spectroscopy was used to investigate the low-frequency vibrational dynamics of the heme in the carbon monoxide oxidation activator protein (CooA) from the thermophilic anaerobic bacterium Carboxydothermus hydrogenoformans (Ch-CooA). Low frequency vibrational modes are important because they are excited by the ambient thermal bath (k(B)T = 200 cm(-1)) and participate in thermally activated barrier crossing events. However, such modes are nearly impossible to detect in the aqueous phase using traditional spectroscopic methods. Here, we present the low frequency coherence spectra of the ferric, ferrous, and CO-bound forms of Ch-CooA in order to compare the protein-induced heme distortions in its active and inactive states. Distortions take place predominantly along the coordinates of low-frequency modes because of their weak force constants, and such distortions are reflected in the intensity of the vibrational coherence signals. A strong mode near ~90 cm(-1) in the ferrous form of Ch-CooA is suggested to contain a large component of heme ruffling, consistent with the imidazole-bound ferrous heme crystal structure, which shows a significant protein-induced heme distortion along this coordinate. A mode observed at ~228 cm(-1) in the six-coordinate ferrous state is proposed to be the ν(Fe-His) stretching vibration. The observation of the Fe-His mode indicates that photolysis of the N-terminal α-amino axial ligand takes place. This is followed by a rapid (~8.5 ps) transient absorption recovery, analogous to methionine rebinding in photolyzed ferrous cytochrome c. We have also studied CO photolysis in CooA, which revealed very strong photoproduct state coherent oscillations. The observation of heme-CO photoproduct oscillations is unusual because most other heme systems have CO rebinding kinetics that are too slow to make the measurement possible. The low frequency coherence spectrum of the CO-bound form of Ch-CooA shows a strong vibration at ~230 cm(-1) that is broadened and up-shifted compared to the ν(Fe-His) of Rr-CooA (216 cm(-1)). We propose that the stronger Fe-His bond is related to the enhanced thermal stability of Ch-CooA and that there is a smaller (time dependent) tilt of the histidine ring with respect to the heme plane in Ch-CooA. The appearance of strong modes at ~48 cm(-1) in both the ferrous and CO-bound forms of Ch-CooA is consistent with coupling of the heme doming distortion to the photolysis reaction in both samples. Upon CO binding and protein activation, a heme mode near 112 ± 5 cm(-1) disappears, probably indicating a decreased heme saddling distortion. This reflects changes in the heme environment and geometry that must be associated with the conformational transition activating the DNA-binding domain. Protein-specific DNA binding to the CO-bound form of Ch-CooA was also investigated, and although the CO rebinding kinetics are significantly perturbed, there are negligible changes in the low-frequency vibrational spectrum of the heme.

Spectroscopic identification of reactive porphyrin motions

Nuclear resonance vibrational spectroscopy (NRVS) reveals the vibrational dynamics of a Mössbauer probe nucleus. Here, (57)Fe NRVS measurements yield the complete spectrum of Fe vibrations in halide complexes of iron porphyrins. Iron porphine serves as a useful symmetric model for the more complex spectrum of asymmetric heme molecules that contribute to numerous essential biological processes. Quantitative comparison with the vibrational density of states (VDOS) predicted for the Fe atom by density functional theory calculations unambiguously identifies the correct sextet ground state in each case. These experimentally authenticated calculations then provide detailed normal mode descriptions for each observed vibration. All Fe-ligand vibrations are clearly identified despite the high symmetry of the Fe environment. Low frequency molecular distortions and acoustic lattice modes also contribute to the experimental signal. Correlation matrices compare vibrations between different molecules and yield a detailed picture of how heme vibrations evolve in response to (a) halide binding and (b) asymmetric placement of porphyrin side chains. The side chains strongly influence the energetics of heme doming motions that control Fe reactivity, which are easily observed in the experimental signal.

Femtosecond Coherence Spectroscopic Study of the Onset of Chemical Denaturation of Myoglobin upon Addition of Minor Amounts of Urea

The interaction of urea with myoglobin, as a benchmark system for heme-containing proteins, is studiedviafemtosecond coherence spectroscopy. The work focuses on the effect of urea on the appearance of low-wavenumber oscillations, which are a measure of the geometrical structure of the heme group and its interaction with the polypeptide chain. Pursuing this approach, structural alterations (i.e.changes in the vibrational dynamics of the heme group) are detected at denaturant concentrations below the full denaturation limit of 6 M urea for myoglobin. In particular, the low-wavenumber oscillation associated with the heme-doming (i.e.the out-off-plane vibration of the propyrin macrocycle) is found to appear spectrally shifted with a concentration of only 3 M urea. These results suggest that the local environment around the heme is already altered despite the fact that macroscopic unfolding as manifested in the thermodynamic properties of the polypeptide chain is not complete at these urea concentrations.

Investigation of the low frequency dynamics of heme proteins: native and mutant cytochrome P450(cam) and redox partner complexes

Vibrational coherence spectroscopy (VCS) is used to investigate the low-frequency dynamics of camphor-free and camphor-bound cytochrome P450(cam) (CYP 101) and its L358P mutant. The low-frequency heme vibrations are found to be perturbed upon binding to the electron transfer partner putidaredoxin (Pdx). A strong correlation between the "detuned" vibrational coherence spectrum, which monitors frequencies between 100 and 400 cm(-1), and the lower frequency part of the Raman spectrum is also demonstrated. The very low frequency region ≤200 cm(-1), uniquely accessed by open-band VCS measurements, reveals a mode near 103 cm(-1) in P450(cam) when camphor is not present in the distal pocket. This reflects the presence of a specific heme distortion, such as saddling or ruffling, in the substrate-free state where water is coordinated to the low-spin iron atom. Such distortions are likely to retard the rate of electron transfer to the substrate-free protein. The presence of strong mode near ∼33 cm(-1) in the camphor-bound form suggests a significant heme-doming distortion, which is supported by analysis using normal coordinate structural decomposition. Pdx also displays a strong coherent vibration near 30 cm(-1) that in principle could be involved in vibrational resonance with its electron transfer target. A splitting of the 33 cm(-1) feature and intensification of a mode near 78 cm(-1) appear when the P450(cam)/Pdx complex is formed. These observations are consistent with vibrational mixing and heme geometric distortions upon Pdx binding that are coincident with the increased thiolate electron donation to the heme. The appearance of a mode near 65 cm(-1) in the coherence spectra of the L358P mutant is comparable to the mode at 78 cm(-1) seen in the P450(cam)/Pdx complex and is consistent with the view that the heme and its environment in the L358P mutant are similar to the Pdx-bound native protein. Resonance Raman spectra are presented for both P450(cam) and the L358P mutant and the changes are correlated with an increased amount of thiolate electron donation to the heme in the mutant sample.

Investigations of low-frequency vibrational dynamics and ligand binding kinetics of cystathionine beta-synthase

Vibrational coherence spectroscopy is used to study the low frequency dynamics of the truncated dimer of human cystathionine beta-synthase (CBS). CBS is a pyridoxal-5'-phosphate-dependent heme enzyme with cysteine and histidine axial ligands that catalyzes the condensation of serine and homocysteine to form cystathionine. A strong correlation between the "detuned" coherence spectrum (which probes higher frequencies) and the Raman spectrum is demonstrated, and a rich pattern of modes below 200 cm(-1) is revealed. Normal coordinate structural decomposition (NSD) of the ferric CBS crystal structure predicts the enhancement of normal modes with significant heme "doming", "ruffling", and "saddling" content, and they are observed in the coherence spectra near approximately 40, approximately 60, and approximately 90 cm(-1). When pH is varied, the relative intensities and frequencies of the low frequency heme modes indicate the presence of a unique protein-induced heme structural perturbation near pH 7 that differs from what is observed at higher or lower pH. For ferric CBS, we observe a new mode near approximately 25 cm(-1), possibly involving the response of the protein, which exhibits a phase jump of approximately pi for excitation on the blue and red side of the Soret band maximum. The low frequency vibrational coherence spectrum of ferrous CBS is also presented, along with our efforts to probe its NO-bound complex. The CO geminate rebinding kinetics of CBS are similar to the CO-bound form of the gene activator protein CooA, but with the appearance of a significant additional kinetic inhomogeneity. Analysis of this inhomogeneity suggests that it arises from the two subunits of CBS and leads to a factor of approximately 20 for the ratio of the average CO geminate rebinding rates of the two subunits.

Measurements of heme relaxation and ligand recombination in strong magnetic fields

Heme cooling signals and diatomic ligand recombination kinetics are measured in strong magnetic fields (up to 10 T). We examined diatomic ligand recombination to heme model compounds (NO and CO), myoglobin (NO and O(2)), and horseradish peroxidase (NO). No magnetic field induced rate changes in any of the samples were observed within the experimental detection limit. However, in the case of CO binding to heme in glycerol and O(2) binding to myoglobin, we observe a small magnetic field dependent change in the early time amplitude of the optical response that is assigned to heme cooling. One possibility, consistent with this observation, is that there is a weak magnetic field dependence of the nonradiative branching ratio into the vibrationally hot electronic ground state during CO photolysis. Ancillary studies of the "spin-forbidden" CO binding reaction in a variety of heme compounds in the absence of magnetic field demonstrate a surprisingly wide range for the Arrhenius prefactor. We conclude that CO binding to heme is not always retarded by unfavorable spin selection rules involving a double spin-flip superexchange mechanism. In fact, it appears that the small prefactor ( approximately 10(9) s(-1)) found for CO rebinding to Mb may be anomalous, rather than the general rule for heme-CO rebinding. These results point to unresolved fundamental issues that underlie the theory of heme-ligand photolysis and rebinding.

Direct evidence for mode-specific vibrational energy relaxation from quantum time-dependent perturbation theory. III. The nu(4) and nu(7) modes of nonplanar nickel porphyrin models

The time scales and pathways of vibrational energy relaxation (VER) of the nu(4) and nu(7) modes of three nickel porphyrin models, nickel porphine (NiP), nickel protoporphyrin IX (Ni-heme), and nickel octaethylporphyrin (NiOEP), were studied using a non-Markovian time-dependent perturbation theory at the B3LYP/6-31G(d) level. When NiP is calculated with D(4h) symmetry, it has the planar structure and the same VER properties as ferrous iron porphine (FeP). The porphine cores of both Ni-heme and NiOEP were distorted from a planar geometry, assuming a nonplanar structure, similar to that of the heme structure in cytochrome c. The VER time scales of Ni-heme are found to be similar to those predicted for a planar iron heme, but the derived pathways have distinctly different features. In particular, the strong coupling between the nu(7) mode and the overtone of the approximately 350 cm(-1) gamma(7) mode, observed for planar porphyrins, is absent in both nonplanar nickel porphyrins. Direct energy exchange between the nu(4) and nu(7) modes is not observed in NiOEP, but is found to play an essential role in the VER of the nu(4) mode in Ni-heme. The Ni-heme isopropionate groups are involved in the dominant VER pathways of both the nu(4) and nu(7) modes of Ni-heme. However, in contrast with VER pathways derived in planar iron heme, the isopropionate groups are not observed to play an essential role relative to other side chains in spatially directing the vibrational energy flow.

Low-frequency dynamics of Caldariomyces fumago chloroperoxidase probed by femtosecond coherence spectroscopy

Ultrafast laser spectroscopy techniques are used to measure the low-frequency vibrational coherence spectra and nitric oxide rebinding kinetics of Caldariomyces fumago chloroperoxidase (CPO). Comparisons of the CPO coherence spectra with those of other heme species are made to gauge the protein-specific nature of the low-frequency spectra. The coherence spectrum of native CPO is dominated by a mode that appears near 32-33 cm(-1) at all excitation wavelengths, with a phase that is consistent with a ground-state Raman-excited vibrational wavepacket. On the basis of a normal coordinate structural decomposition (NSD) analysis, we assign this feature to the thiolate-bound heme doming mode. Spectral resolution of the probe pulse ("detuned" detection) reveals a mode at 349 cm(-1), which has been previously assigned using Raman spectroscopy to the Fe-S stretching mode of native CPO. The ferrous species displays a larger degree of spectral inhomogeneity than the ferric species, as reflected by multiple shoulders in the optical absorption spectra. The inhomogeneities are revealed by changes in the coherence spectra at different excitation wavelengths. The appearance of a mode close to 220 cm(-1) in the coherence spectrum of reduced CPO excited at 440 nm suggests that a subpopulation of five coordinated histidine-ligated hemes is present in the ferrous state at a physiologically relevant pH. A significant increase in the amplitude of the coherence signal is observed for the resonance with the 440 nm subpopulation. Kinetics measurements reveal that nitric oxide binding to ferric and ferrous CPO can be described as a single-exponential process, with rebinding time constants of 29.4 +/- 1 and 9.3 +/- 1 ps, respectively. This is very similar to results previously reported for nitric oxide binding to horseradish peroxidase.

Coherence spectroscopy investigations of the low-frequency vibrations of heme: effects of protein-specific perturbations

Femtosecond coherence spectroscopy is used to probe the low-frequency (20-200 cm(-1)) vibrational modes of heme proteins in solution. Horseradish peroxidase (HRP), myoglobin (Mb), and Campylobacter jejuni globin (Cgb) are compared and significant differences in the coherence spectra are revealed. It is concluded that hydrogen bonding and ligand charge do not strongly affect the low-frequency coherence spectra and that protein-specific deformations of the heme group lower its symmetry and control the relative spectral intensities. Such deformations potentially provide a means for proteins to tune heme reaction coordinates, so that they can perform a broad array of specific functions. Native HRP displays complex spectral behavior above approximately 50 cm(-1) and very weak activity below approximately 50 cm(-1). Binding of the substrate analog, benzhydroxamic acid, leads to distinct changes in the coherence and Raman spectra of HRP that are consistent with the stabilization of a heme water ligand. The CN derivatives of the three proteins are studied to make comparisons under conditions of uniform heme coordination and spin-state. MbCN is dominated by a doming mode near 40 cm(-1), while HRPCN displays a strong oscillation at higher frequency (96 cm(-1)) that can be correlated with the saddling distortion observed in the X-ray structure. In contrast, CgbCN displays low-frequency coherence spectra that contain strong modes near 30 and 80 cm(-1), probably associated with a combination of heme doming and ruffling. HRPNO displays a strong doming mode near 40 cm(-1) that is activated by photolysis. The damping of the coherent motions is significantly reduced when the heme is shielded from solvent fluctuations by the protein material and reduced still further when T approximately < 50 K, as pure dephasing processes due to the protein-solvent phonon bath are frozen out.

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.

Ultrafast dynamics of ligands within heme proteins

Physiological bond formation and bond breaking events between proteins and ligands and their immediate consequences are difficult to synchronize and study in general. However, diatomic ligands can be photodissociated from heme, and thus in heme proteins ligand release and rebinding dynamics and trajectories have been studied on timescales of the internal vibrations of the protein that drive many biochemical reactions, and longer. The rapidly expanding number of characterized heme proteins involved in a large variety of functions allows comparative dynamics-structure-function studies. In this review, an overview is given of recent progress in this field, and in particular on initial sensing processes in signaling proteins, and on ligand and electron transfer dynamics in oxidases and cytochromes.

Temperature-dependent heme kinetics with nonexponential binding and barrier relaxation in the absence of protein conformational substates

We present temperature-dependent kinetic measurements of ultrafast diatomic ligand binding to the "bare" protoheme (L(1)-FePPIX-L(2), where L(1) = H(2)O or 2-methyl imidazole and L(2) = CO or NO). We found that the binding of CO is temperature-dependent and nonexponential over many decades in time, whereas the binding of NO is exponential and temperature-independent. The nonexponential nature of CO binding to protoheme, as well as its relaxation above the solvent glass transition, mimics the kinetics of CO binding to myoglobin (Mb) but on faster time scales. This demonstrates that the nonexponential kinetic response observed for Mb is not necessarily due to the presence of protein conformational substates but rather is an inherent property of the solvated heme. The nonexponential kinetic data were analyzed by using a linear coupling model with a distribution of enthalpic barriers that fluctuate on slower time scales than the heme-CO recombination time. Below the solvent glass transition (T(g) approximately 180 K), the average enthalpic rebinding barrier for H(2)O-PPIX-CO was found to be approximately 1 kJ/mol. Above T(g), the barrier relaxes and is approximately 6 kJ/mol at 290 K. Values for the first two moments of the heme doming coordinate distribution extracted from the kinetic data suggest significant anharmonicity above T(g). In contrast to Mb, the protoheme shows no indication of the presence of "distal" enthalpic barriers. Moreover, the wide range of Arrhenius prefactors (10(9) to 10(11) s(-1)) observed for CO binding to heme under differing conditions suggests that entropic barriers may be an important source of control in this class of biochemical reactions.

Low frequency spectral density of ferrous heme: perturbations induced by axial ligation and protein insertion

Femtosecond coherence spectroscopy is used to probe low frequency (20-400 cm(-1)) modes of the ferrous heme group in solution, with and without 2-methyl imidazole (2MeIm) as an axial ligand. The results are compared to heme proteins (CPO, P450(cam), HRP, Mb) where insertion of the heme into the protein results in redistribution of the low frequency spectral density and in (approximately 60%) longer damping times for the coherent signals. The major effect of imidazole ligation to the ferrous heme is the "softening" of the low frequency force constants by a factor of approximately 0.6 +/- 0.1. The functional consequences of imidazole ligation are assessed and it is found that the enthalpic CO rebinding barrier is increased significantly when imidazole is bound. The force constant softening analysis, combined with the kinetics results, indicates that the iron is displaced by only approximately 0.2 A from the heme plane in the absence of the imidazole ligand, whereas it is displaced by approximately 0.4 A when imidazole (histidine) is present. This suggests that binding of imidazole (histidine) as an axial ligand, and the concomitant softening of the force constants, leads to an anharmonic distortion of the heme group that has significant functional consequences.

Hydrophobic Distal Pocket Affects NO−Heme Geminate Recombination Dynamics in Dehaloperoxidase and H64V Myoglobin

The recombination dynamics of NO with dehaloperoxidase (DHP) from Amphitrite ornata following photolysis were measured by femtosecond time-resolved absorption spectroscopy. Singular value decomposition (SVD) analysis reveals two important basis spectra. The first SVD basis spectrum reports on the population of photolyzed NO molecules and has the appearance of the equilibrium difference spectrum between the deoxy and NO forms of DHP. The first basis time course has two kinetic components with time constants of tau(11) approximately 9 ps and tau(12) approximately 50 ps that correspond to geminate recombination. The fast geminate process tau(11) arises from a contact pair with the heme iron in a bound state with S = 3/2 spin. The slow geminate process tau(12) corresponds to the recombination from a more remote docking site >3 A from the heme iron with the greater barrier corresponding to a S = 5/2 spin state. The second SVD basis spectrum represents a time-dependent Soret band shift indicative of heme photophysical processes and protein relaxation with time constants of tau(21) approximately 3 ps and tau(22) approximately 17 ps, respectively. A comparison between the more rapid rate constant of the slow geminate phase in DHP-NO and horse heart myoglobin (HHMbNO) or sperm whale myoglobin (SWMbNO) suggests that protein interactions with photolyzed NO are weaker in DHP than in the wild-type MbNOs, consistent with the hydrophobic distal pocket of DHP. The slower protein relaxation rate tau(22) in DHP-NO relative to HHMbNO implies less effective trapping in the docking site of the distal pocket and is consistent with a greater yield for the fast geminate process. The trends observed for DHP-NO also hold for the H64V mutant of SWMb (H64V MbNO), consistent with a more hydrophobic distal pocket for that protein as well. We examine the influence of solution viscosity on NO recombination by varying the glycerol content in the range from 0% to 90% (v/v). The dominant effect of increasing viscosity is the increase of the rate of the slow geminate process, tau(12), coupled with a population decrease of the slow geminate component. Both phenomena are similar to the effect of viscosity on wild-type Mb due to slowing of protein relaxation resulting from an increased solution viscosity and protein surface dehydration.

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.

Spectrally resolved femtosecond 2-colour 3-pulse photon echoes: a new spectroscopic tool to study molecular dynamics

We present a new multidimensional femtosecond spectroscopy technique based on spectrally resolved 2-colour 3-pulse photon echoes for investigating molecular dynamics in a variety of systems including proteins. In this technique the sample is illuminated by two femtosecond 'pump' pulses with wave vectors k1, k2 and wavelength lambda(pump) and a femtosecond 'probe' pulse with wave vector k3 and wavelength lambda(probe). Nonlinear signals are generated in the phase-matching directions k4 = - k1 + k2 + k3 and k6 = - k3 + k1 + k2. These signals are analysed in spectrometers equipped with CCD detectors and the spectra of the signals are recorded for various values of (i) the delay t12 between pulses 1 and 2, (ii) the delay t23 between pulses 2 and 3, and (iii) the wavelengths lambda(pump), lambda(probe). The technique has been used for studying vibrational and electronic dynamics of dye molecules, such as cresyl violet in methanol, and ultra-fast transient processes that occur during the photo-dissociation of carbonmonoxy myoglobin (MbCO) into myoglobin (Mb) and CO.

Time-resolved visible and infrared study of the cyano complexes of myoglobin and of hemoglobin I from Lucina pectinata

The dynamics of the ferric CN complexes of the heme proteins Myoglobin and Hemoglobin I from the clam Lucina pectinata upon Soret band excitation is monitored using infrared and broad band visible pump-probe spectroscopy. The transient response in the UV-vis spectral region does not depend on the heme pocket environment and is very similar to that known for ferrous proteins. The main feature is an instantaneous, broad, short-lived absorption signal that develops into a narrower red-shifted Soret band. Significant transient absorption is also observed in the 360-390 nm range. At all probe wavelengths the signal decays to zero with a longest time constant of 3.6 ps. The infrared data on MbCN reveal a bleaching of the C triple bond N stretch vibration of the heme-bound ligand, and the formation of a five-times weaker transient absorption band, 28 cm(-1) lower in energy, within the time resolution of the experiment. The MbC triple bond N stretch vibration provides a direct measure for the return of population to the ligated electronic (and vibrational) ground state with a 3-4 ps time constant. In addition, the CN-stretch frequency is sensitive to the excitation of low frequency heme modes, and yields independent information about vibrational cooling, which occurs on the same timescale.