Picosecond resonance Raman evidence for unrelaxed heme in the (carbonmonoxy)myoglobin photoproduct

Ultrafast myoglobin structural dynamics observed with an X-ray free-electron laser

Light absorption can trigger biologically relevant protein conformational changes. The light-induced structural rearrangement at the level of a photoexcited chromophore is known to occur in the femtosecond timescale and is expected to propagate through the protein as a quake-like intramolecular motion. Here we report direct experimental evidence of such 'proteinquake' observed in myoglobin through femtosecond X-ray solution scattering measurements performed at the Linac Coherent Light Source X-ray free-electron laser. An ultrafast increase of myoglobin radius of gyration occurs within 1 picosecond and is followed by a delayed protein expansion. As the system approaches equilibrium it undergoes damped oscillations with a ~3.6-picosecond time period. Our results unambiguously show how initially localized chemical changes can propagate at the level of the global protein conformation in the picosecond timescale.

Comparing the temperature dependence of photosynthetic electron transfer in Chloroflexus aurantiacus and Rhodobactor sphaeroides reaction centers

The process of electron transfer from the special pair, P, to the primary electron donor, H(A), in quinone-depleted reaction centers (RCs) of Chloroflexus (Cf.) aurantiacus has been investigated over the temperature range from 10 to 295 K using time-resolved pump-probe spectroscopic techniques. The kinetics of the electron transfer reaction, P* → P(+)H(A)(-), was found to be nonexponential, and the degree of nonexponentiality increased strongly as temperature decreased. The temperature-dependent behavior of electron transfer in Cf. aurantiacus RCs was compared with that of the purple bacterium Rhodobacter (Rb.) sphaeroides . Distinct transitions were found in the temperature-dependent kinetics of both Cf. aurantiacus and Rb. sphaeroides RCs, at around 220 and 160 K, respectively. Structural differences between these two RCs, which may be associated with those differences, are discussed. It is suggested that weaker protein-cofactor hydrogen bonding, stronger electrostatic interactions at the protein surface, and larger solvent interactions likely contribute to the higher transition temperature in Cf. aurantiacus RCs temperature-dependent kinetics compared with that of Rb. sphaeroides RCs. The reaction-diffusion model provides an accurate description for the room-temperature electron transfer kinetics in Cf. aurantiacus RCs with no free parameters, using coupling and reorganization energy values previously determined for Rb. sphaeroides , along with an experimental measure of protein conformational diffusion dynamics and an experimental literature value of the free energy gap between P* and P(+)H(A)(-).

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.

On the origin of heme absorption band shifts and associated protein structural relaxation in myoglobin following flash photolysis

The role of the protein structural change monitored by absorption band shifts following flash photolysis of CO from myoglobin is discussed in terms of structure-function relationships. Evidence is presented that the Soret band shift does not depend primarily on the covalent linkage of the heme iron to the protein by using the mutation H93G(L) in which the proximal histidine 93 is replaced by glycine and an exogenous ligand L, which coordinates the heme iron but is not covalently bound to the globin. While CO rebinding kinetics depend strongly on the nature of the exogenous ligand L in H93G(L), the magnitude and time evolution of the Soret band shift in a viscous buffer on the nanosecond time scale are hardly perturbed in all cases studied. Comparison of the Soret band and charge transfer Band III shifts demonstrates that both have a similar time dependence on the nanosecond to microsecond time scale following flash photolysis in viscous solvents. We conclude that the nonexponential kinetics of protein relaxation probed by transient absorption band position shifts involves primarily distal coordinates prior to ligand escape. This result agrees with earlier measurements of Soret band shifts in distal pocket mutants of myoglobin (1). We suggest that the band shifts are primarily a response to changes in the electrostatic field around the heme (a transient Stark shift) associated with changes in protein structure that occur following ligand photodissociation.

Transient and time-resolved resonance Raman investigation of photoinitiated electron transfer in ruthenated cytochromes c

Ruthenation of exterior amino acid residues of heme proteins provides an effective means by which biological ET reactions can be studied within the context of highly complex protein environments. Resonance Raman spectroscopy can probe both ET kinetics and structural dynamics at the molecular level. Here we present the first comprehensive use of time-resolved and transient resonance Raman spectroscopies to examine photoinduced ET in cytochromes. Two ruthenated cytochromes c, Ru(lys72)-cyt.c and Ru(cyt102)cyt.c, were studied with TRRS using 10 ns laser pulses and with TRRRS on a 10 ns to 10 ms time scale. It was found that resonance Raman protocols can effectively trace ET kinetics and associated heme--protein structural dynamics. Care must be exercised, however, when drawing comparisons to measurements made by other methods (i.e., transient absorbance). The TRRS studies directly probe the heme and its local environment and reveal that the heme dynamics accompanying ET are very rapid relative to phenomenological ET kinetics. The heme and its local environment evolve to their equilibrium (ferrous) structure in less than 10 ns subsequent to ET, with no evidence for the existence of metastable heme pocket geometries analogous to those observed in the dynamic response of hemoglobins and oxidases. Finally, species-specific differences are observed in the photoinduced ET kinetics and heme structural dynamics. However, these differences are confined to nanosecond or faster time scales.

Photolysis-induced structural changes in single crystals of carbonmonoxy myoglobin at 40 K

Myoglobin's reversible binding of oxygen is a model for studies of protein control of ligand binding and discrimination. Protein relaxation and geminate ligand rebinding subsequent to ligand photodissociation have been studied extensively by a variety of techniques. The ps to ns time scales for these processes are still much shorter than the ms time resolution of X-ray diffraction experiments, but it may be possible to trap these intermediates at low temperatures. We report here an X-ray diffraction investigation of structural changes induced by photolysis of carbonmonoxy myoglobin crystals at 40 K. Our results provide a structural basis for the interpretation of ambient and low temperature spectroscopic observations and molecular dynamics simulations of the ligand photodissociation and binding processes in haem proteins.

Conformational relaxation and ligand binding in myoglobin

Absorption spectroscopy with nanosecond time resolution shows that myoglobin undergoes conformational relaxation on the same time scale as geminate rebinding of carbon monoxide. Ligand rebinding following photodissociation of the heme-CO complex was measured from the amplitude of the average difference spectrum, while conformational changes were measured from changes in the detailed shape of the Soret spectra of the deoxyhemes. Experiments in which the solvent viscosity was varied between 1 and 300 cP and the temperature between 268 and 308 K were analyzed by fitting the multiwavelength kinetic data with both empirical and molecular models. Novel numerical techniques were employed in fitting the data, including the use of singular value decomposition to remove the effects of temperature and solvent on the spectra and of a Monte Carlo method to overcome the multiple minimum problem in searching parameter space. The molecular model is the minimal model that incorporates all of the major features of myoglobin kinetics at ambient temperatures, including a fast and slow rebinding conformation and two geminate states for each conformation. The results of fitting the kinetic data with this model indicate that the geminate-rebinding rates for the two conformations differ by at least a factor of 100. The differences between the spectra of the two conformations generated from the fits are similar to the differences between those of the R and T conformations of hemoglobin. In modeling the data, the dependence of the rates on temperature and viscosity was parametrized using a modification of Kramers theory which includes the contributions of both protein and solvent to the friction. The rate of the transition from the fast to the slow rebinding conformation is found to be inversely proportional to the viscosity when the viscosity exceeds about 30 cP and nearly viscosity independent at low viscosity. The viscosity dependence at high viscosities suggests that the two conformations differ by the global displacement of protein atoms on the proximal side of the heme observed by X-ray crystallography. We suggest that the conformational change observed in our experiments corresponds to the final portion of the nonexponential conformational relaxation recently observed by Anfinrud and co-workers, which begins on a picosecond time scale. Furthermore, extrapolation of our data to temperatures near that of the solvent glass transition suggests that this conformational relaxation may very well be the one postulated by Frauenfelder and co-workers to explain the decrease in the rate of geminate rebinding with increasing temperature above 180 K.

Picosecond study of the near infrared absorption band of hemoglobin after photolysis of carbonmonoxyhemoglobin

Picosecond absorption spectroscopy is used to examine the position and band shape of the near infrared absorption band of hemoglobin as a function of time after the photodissociation of CO from carbonmonoxyhemoglobin. For the earliest delay time probed, 35 ps, the peak of the transient spectrum is at 765 nm, red shifted by 6 nm from that characteristic of equilibrium deoxyhemoglobin. No evolution in either the peak position or band shape is observed for time delays up to 60 ns. In addition, the position and shape of the spectrum are independent of photolysis energies ranging from 15 microJ/pulse to 150 microJ/pulse, spanning conditions under which the photon/heme ratio is varied from 0.01 to 2.0. This indicates that the geometry in the heme group is unrelaxed and that equilibration of the surrounding protein structure occurs on a time scale longer than 60 ns.

Applications of ultrafast laser spectroscopy for the study of biological systems

The discovery of mode-locked laser operation now nearly two decades ago has started a development which enables researchers to probe the dynamics of ultrafast physical and chemical processes at the molecular level on shorter and shorter time scales. Naturally the first applications were in the fields of photophysics and photochemistry where it was then possible for the first time to probe electronic and vibrational relaxation processes on a sub-nanosecond timescale. The development went from lasers producing pulses of many picoseconds to the shortest pulses which are at present just a few femtoseconds long. Soon after their discovery ultrashort pulses were applied also to biological systems which has revealed a wealth of information contributing to our understanding of a broadrange of biological processes on the molecular level.It is the aim of this review to discuss the recent advances and point out some future trends in the study of ultrafast processes in biological systems using laser techniques. The emphasis will be mainly on new results obtained during the last 5 or 6 years. The term ultrafast means that I shall restrict myself to sub-nanosecond processes with a few exceptions.

Contribution of protein conformation to heme stereochemistry and reactivity. Low-temperature magnetic circular dichroism data

Visible and near infrared magnetic circular dichroism (MCD) spectra of heme proteins and enzymes as well as those of a protein-free heme bound to 2-methylimidazole were recorded and compared at 4.2 K in unrelaxed metastable and relaxed equilibrium heme stereochemistry. The relaxed and unrelaxed stereochemistries of a 5-coordinate ferrous heme were generated by chemical reduction of iron at room temperature before freezing the sample and by photolysis of CO or O2 complexes at 4.2 K, respectively. The results are discussed in terms of a protein contribution into energies of the Fe-N epsilon(His) and Fe-N(pyrrols) bonds and their change on a ligand binding. We observed and analyzed cases of weak (myoglobin, hemoglobin) and strong (leghemoglobin, peroxidases) constraints imposed by the protein conformation on the proximal heme stereochemistry by comparing the bond energies in proteins with those in the protoheme-(2-methylimidazole) model compound. The role of a protein moiety in modulating the ligand binding properties of leghemoglobin and the heme reactivity of horseradish peroxidase is discussed.

Effect of solvent viscosity on the heme-pocket dynamics of photolyzed (carbonmonoxy)hemoglobin

The heme-pocket dynamics subsequent to carbon monoxide photolysis from human hemoglobin have been monitored as a function of glycerol-water solvent composition with time-resolved resonance Raman spectroscopy. Prompt (geminate) ligand recombination rates and the transient heme-pocket geometry established within 10 ns after photolysis appear to be largely independent of solvent composition. The rate of relaxation of the transient geometry to an equilibrium deoxy configuration is, however, quite sensitive to solvent composition. These observations suggest that the former processes result from local, internal motions of the protein, while the relaxation dynamics of the proximal heme pocket are predicated upon more global protein motions that are dependent upon solvent viscosity.

Picosecond transient absorption study of photodissociated carboxy hemoglobin and myoglobin

The optical transient absorption spectra at 30 ps and 6.5 ns after photolysis are compared for both carboxy hemoglobin (HbCO) and carboxy myoglobin (MbCO). Both 355- and 532-nm excitation pulses were used. In all cases the shapes of the optical difference spectra thus generated are stationary over the complete time-scale studied. The photolysis spectra for MbCO are not significantly different from the equilibrium difference spectra generated on the same picosecond spectrometer when measured to an accuracy of +/- 0.5 nm. In addition, spectral parameters for delegated HbCO generated on the same spectrometer but detected by two different techniques, either by a Vidicon detector or point by point with photomultiplier tubes, are reported; the results are different from some of the previously reported picosecond experiments.

Heme-linked ionization and ligand binding produce identical changes of proximal heme stereochemistry in reduced horseradish peroxidase. Evidence for existence of two protein conformations

The visible and near infrared magnetic circular dichroism spectra of chemically reduced horseradish peroxidase at neutral and alkaline pH values and 5-coordinate protoheme-(2-methylimidazole) at pH 9.1 were compared at 4.2 K with those of photolysis products of their carbon monoxide complexes. From the results obtained we concluded that: (i) there are two protein conformations of HRP which determine the geometry of the Fe-N(His) bond; (ii) the transition from one conformation (heme stereochemistry) to another can be induced by either heme-linked ionization or ligand binding; (iii) a trigger mechanism for switching between two conformations has to exist.

Molecular dynamics simulations of cooling in laser-excited heme proteins

In transient optical experiments the absorbed photon raises the vibrational temperature of the chromophore. In heme proteins at room temperature conversion of a 530-nm photon into vibrational energy is estimated to raise the temperature of the heme by 500-700 K. Cooling of the heme is expected to occur mainly by interacting with the surrounding protein. We report molecular dynamics simulations for myoglobin and cytochrome c in vacuo that predict that this cooling occurs on the ps time scale. The decay of the vibrational temperature is nonexponential with about 50% loss occurring in 1-4 ps and with the remainder in 20-40 ps. These results predict the presence of nonequilibrium vibrational populations that would introduce ambiguity into the interpretation of transient ps absorption and Raman spectra and influence the kinetics of sub-ns geminate recombination.

Resonance Raman characterization of the 7-ns photoproduct of (carbonmonoxy)hemoglobin: implications for hemoglobin dynamics

Resonance Raman spectra are reported for deoxyhemoglobin (deoxyHb) and the (carbonmonoxy)hemoglobin (HbCO) photoproduct Hb by use of 7-ns YAG laser pulses at wavelengths of 416 and 532 nm, where enhancement is observed for totally symmetric and nontotally symmetric modes, respectively. The frequencies of the porphyrin skeletal modes v10, v2, v19, v11, and v3 have been determined to be 1602, 1559, 1553, 1542, and 1466 cm-1 in Hb. These frequencies are 2-3 cm-1 lower than the corresponding frequencies for deoxyHb. The v19 and v11 frequencies are at the expected values for a Ct-N distance of 2.057 A, the known core size for a 6-coordinate high-spin FeII-porphyrin complex. The remaining frequencies, however, deviate from the core size correlations for these modes in the same direction as do those of deoxyHb, suggesting that the porphyrin ring is domed in both species. Thus, the heme structure is similar for deoxyHb and Hb but is slightly expanded in the latter. The expanded heme in Hb implies a restraint on the full out-of-plane displacement of the Fe atom, by an estimated approximately 0.1 A relative to deoxyHb. This could result from a residual interaction with the CO molecule if the latter remains held by the protein against the Fe atom, in a high-spin 6-coordinate complex. The available spectroscopic evidence suggests that such a complex may be stabilized at 4 K but is unlikely to persist at room temperature beyond the electronic relaxation (0.35 ps) of the electronically excited heme.(ABSTRACT TRUNCATED AT 250 WORDS)

Transient resonance Raman spectroscopy shows unrelaxed heme following CO photodissociation from cytochrome-c peroxidase

The 7 ns 436 nm pulses of an H2-shifted YAG laser have been used to photolyze the CO adduct of cytochrome-c peroxidase and produce the resonance Raman spectrum of the photoproduct. A 3 cm-1 downshift, relative to the spectrum of reduced enzyme, was observed for the porphyrin C-N breathing mode, v4. The downshift diminishes with decreasing CO /protein ratio, implying, in conjunction with a recent study of CO binding, that the unrelaxed heme is associated with adduct having a tilted, H-bonded FeCO unit. The downshift is eliminated when the phosphate buffer concentration is increased from 0.01 to 0.1 M. It is proposed that the heme relaxation under study involves a transition between two conformations, B and A, differing in the disposition of the distal residues, and having different v4 frequencies for unligated Fe(II) heme. Conformation B allows H-bonding to bound CO, and is favored at high CO and phosphate concentrations, while conformation A, which is unfavorable to CO H-bonding, is favored at low CO and phosphate concentrations. The recently reported absence of unrelaxed frequencies in the 7 ns photo-product of the CO adduct of horseradish peroxidase has been confirmed, and is attributed to lower stability for conformation B and a smaller A - B v4 difference.

Metastable photoproducts from carbon monoxide myoglobin

The photoproduct of carbon monoxide myoglobin generated at 4 K and lower has a resonance Raman spectrum characteristic of a high-spin heme but in which the high-frequency core size-sensitive lines are at lower frequency than those in the deoxy preparation. Such differences are not detected in the photoproduct generated at higher temperatures (50 K) or in that generated at room temperature with 10-nsec pulses. The data indicate that at the low temperature (4 K), the heme in the photoproduct is not fully relaxed, and from the data we conclude that the photoproduct has an expanded porphyrin core. We infer that the core size exceeds that in deoxymyoglobin because the rigid protein prevents the highspin iron atom from moving to its full out-of-plane displacement at the very low temperatures.