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

Feeling the Strain: Quantifying Ligand Deformation in Photosynthesis

Structural distortion of protein-bound ligands can play a critical role in enzyme function by tuning the electronic and chemical properties of the ligand molecule. However, quantifying these effects is difficult due to the limited resolution of protein structures and the difficulty of generating accurate structural restraints for nonprotein ligands. Here, we seek to quantify these effects through a statistical analysis of ligand distortion in chlorophyll proteins (CP), where ring deformation is thought to play a role in energy and electron transfer. To assess the accuracy of ring-deformation estimates from available structural data, we take advantage of the C2 symmetry of photosystem II (PSII), comparing ring-deformation estimates for equivalent sites both within and between 113 distinct X-ray and cryogenic electron microscopy PSII structures. Significantly, we find that several deformation modes exhibit considerable variability in predictions, even for equivalent monomers, down to a 2 Å resolution, to an extent that probably prevents their utilization in optical calculations. We further find that refinement restraints play a critical role in determining deformation values to resolution as low as 2 Å. However, for those modes that are well-resolved in the structural data, ring deformation in PSII is strongly conserved across all species tested from cyanobacteria to algae. These results highlight both the opportunities and limitations inherent in structure-based analyses of the bioenergetic and optical properties of CPs and other protein-ligand complexes.

Long-Range Charge Delocalization Mediates the Ultrafast Ligand-to-Metal Charge Transfer Dynamics at the Cu2+-Active Site in Azurin

The blue color metalloprotein in azurin has traditionally been attributed to the intense cysteine-to-Cu²⁺ ligand-to-metal charge transfer transition centered at 628 nm. Although resonance Raman measurements of the Cu²⁺ active site have implied that the LMCT transition electronically couples to the protein scaffold well beyond its primary metal-ligand coordination shell, the structural extent of this electronic coupling and visualization of the protein-mediated charge transfer dynamics have remained elusive. Here, using femtosecond broadband transient absorption and impulsive Raman spectroscopy, we provide direct evidence for a rapid relaxation between two distinct charge transfer states, having different spatial delocalization, within ∼300 fs followed by recombination of charges in subpicosecond time scales. We invoke the formation of a protein-centered radical cation, possibly Trp48 or a Phe residue, within 100 fs substantiating the long-range electronic coupling for the first time beyond the traditional copper active site. The Raman spectra of the excited CT state show the presence of protein-centric vibrations along with the vibrational modes assigned to the copper active site. Our results demonstrate a large delocalization length scale of the initially populated CT state, thereby highlighting the possibility of exploiting azurin photochemistry for energy conversion techniques.

Ultra-broadband terahertz fingerprint spectrum of melatonin with vibrational mode analysis

Melatonin (MLT), as a neurotransmitter and an endogenous neurohormone, plays an important role in physiological functions through interactions with specific receptors. The conformations of MLT are closely related to its biological activities and functions. However, the internal relationship between the structure and interaction of MLT and its allosteric transition remains unclear. In this work, we obtain the broadband fingerprint terahertz (THz) spectrum of MLT in the range of 0.5-18 THz using the air-plasma terahertz time-domain spectroscopy (THz-TDS) system. DFT calculations are employed to analyze the vibration characteristics of MLT. The result shows that the low-frequency vibrations mainly come from the strong coupling between inter- and intramolecular vibrations, and the contribution of intramolecular vibrations gradually dominates with increasing frequency. Meanwhile, the local vibrations of the different functional groups distribute widely in the THz low-frequency band, relating to the diversity of conformational changes in the molecule. The intermolecular hydrogen bonds (HBs) have distinct resonant responses and play critical roles in the THz low-frequency vibrations. The study reveals the complex characteristics of the resonant coupling of MLT with THz electromagnetic waves. The results will help to understand the conformational preferences of MLT in neural signal transmission processes.

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.

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.

Ultrafast excited-state charge-transfer dynamics in laccase type I copper site

Femtosecond pump-probe spectroscopy was used to investigate the excited state dynamics of the T1 copper site of laccase from Pleurotus ostreatus, by exciting its 600 nm charge transfer band with a 15-fs pulse and probing over a broad range in the visible region. The decay of the pump-induced ground-state bleaching occurs in a single step and is modulated by clearly visible oscillations. Global analysis of the two-dimensional differential transmission map shows that the excited state exponentially decays with a time constant of 375 fs, thus featuring a decay rate slower than those occurring in quite all the investigated T1 copper site proteins. The ultrashort pump pulse induces a vibrational coherence in the protein, which is mainly assigned to ground state activity, as expected in a system with fast excited state decay. Vibrational features are discussed also in comparison with the traditional resonance Raman spectrum of the enzyme. The results indicate that both excited state dynamics and vibrational modes associated with the T1 Cu laccase charge transfer have main characteristics similar to those of all the T1 copper site-containing proteins. On the other hand, the differences observed for laccase from P. ostreatus further confirm the peculiar hypothesized trigonal T1 Cu site geometry.

Concepts in bio-molecular spectroscopy: vibrational case studies on metalloenzymes

Challenges and chances in bio-molecular spectroscopy are exemplified by vibrational case studies on metalloenzymes.

Ultrafast infrared spectroscopy reveals water-mediated coherent dynamics in an enzyme active site

Ultrafast infrared spectroscopy provides insights into the dynamic nature of water in the active sites of catalase and peroxidase enzymes.

Understanding the impact of fast dynamics upon the chemical processes occurring within the active sites of proteins and enzymes is a key challenge that continues to attract significant interest, though direct experimental insight in the solution phase remains sparse. Similar gaps in our knowledge exist in understanding the role played by water, either as a solvent or as a structural/dynamic component of the active site. In order to investigate further the potential biological roles of water, we have employed ultrafast multidimensional infrared spectroscopy experiments that directly probe the structural and vibrational dynamics of NO bound to the ferric haem of the catalase enzyme from Corynebacterium glutamicum in both H₂O and D₂O. Despite catalases having what is believed to be a solvent-inaccessible active site, an isotopic dependence of the spectral diffusion and vibrational lifetime parameters of the NO stretching vibration are observed, indicating that water molecules interact directly with the haem ligand. Furthermore, IR pump–probe data feature oscillations originating from the preparation of a coherent superposition of low-frequency vibrational modes in the active site of catalase that are coupled to the haem ligand stretching vibration. Comparisons with an exemplar of the closely-related peroxidase enzyme family shows that they too exhibit solvent-dependent active-site dynamics, supporting the presence of interactions between the haem ligand and water molecules in the active sites of both catalases and peroxidases that may be linked to proton transfer events leading to the formation of the ferryl intermediate Compound I. In addition, a strong, water-mediated, hydrogen bonding structure is suggested to occur in catalase that is not replicated in peroxidase; an observation that may shed light on the origins of the different functions of the two enzymes.

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 heme distortion, low-frequency vibrational excitations, and electron transfer in cytochrome c

Cytochrome (cyt) c is an important electron transfer protein. The ruffling deformation of its heme cofactor has been suggested to relate to its electron transfer rate. However, there is no direct experimental evidence demonstrating this correlation. In this work, we studied Pseudomonas aeruginosa cytochrome c551 and its F7A mutant. These two proteins, although similar in their X-ray crystal structure, display a significant difference in their heme out-of-plane deformations, mainly along the ruffling coordinate. Resonance Raman and vibrational coherence measurements also indicate significant differences in ruffling-sensitive modes, particularly the low-frequency γa mode found between ∼50-60 cm(-1). This supports previous assignments of γa as having a large ruffling content. Measurement of the photoreduction kinetics finds an order of magnitude decrease of the photoreduction cross-section in the F7A mutant, which has nearly twice the ruffling deformation as the WT. Additional measurements on cytochrome c demonstrate that heme ruffling is correlated exponentially with the electron transfer rates and suggest that ruffling could play an important role in redox control. A major relaxation of heme ruffling in cytochrome c, upon binding to the mitochondrial membrane, is discussed in this context.

The influence of heme ruffling on spin densities in ferricytochromes c probed by heme core 13C NMR

The heme in cytochromes c undergoes a conserved out-of-plane distortion known as ruffling. For cytochromes c from the bacteria Hydrogenobacter thermophilus and Pseudomonas aeruginosa , NMR and EPR spectra have been shown to be sensitive to the extent of heme ruffling and to provide insights into the effect of ruffling on the electronic structure. Through the use of mutants of each of these cytochromes that differ in the amount of heme ruffling, NMR characterization of the low-spin (S = ½) ferric proteins has confirmed and refined the developing understanding of how ruffling influences the spin distribution on heme. The chemical shifts of the core heme carbons were obtained through site-specific labeling of the heme via biosynthetic incorporation of (13)C-labeled 5-aminolevulinic acid derivatives. Analysis of the contact shifts of these core heme carbons allowed Fermi contact spin densities to be estimated and changes upon ruffling to be evaluated. The results allow a deconvolution of the contributions to heme hyperfine shifts and a test of the influence of heme ruffling on the electronic structure and hyperfine shifts. The data indicate that as heme ruffling increases, the spin densities on the β-pyrrole carbons decrease while the spin densities on the α-pyrrole carbons and meso carbons increase. Furthermore, increased ruffling is associated with stronger bonding to the heme axial His ligand.

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.

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.

NMR and DFT investigation of heme ruffling: functional implications for cytochrome c

Out-of-plane (OOP) deformations of the heme cofactor are found in numerous heme-containing proteins and the type of deformation tends to be conserved within functionally related classes of heme proteins. We demonstrate correlations between the heme ruffling OOP deformation and the (13)C and (1)H nuclear magnetic resonance (NMR) hyperfine shifts of heme aided by density functional theory (DFT) calculations. The degree of ruffling in the heme cofactor of Hydrogenobacter thermophilus cytochrome c(552) has been modified by a single amino acid mutation in the second coordination sphere of the cofactor. The (13)C and (1)H resonances of the cofactor have been assigned using one- and two-dimensional NMR spectroscopy aided by selective (13)C-enrichment of the heme. DFT has been used to predict the NMR hyperfine shifts and electron paramagnetic resonance (EPR) g-tensor at several points along the ruffling deformation coordinate. The DFT-predicted NMR and EPR parameters agree with the experimental observations, confirming that an accurate theoretical model of the electronic structure and its response to ruffling has been established. As the degree of ruffling increases, the heme methyl (1)H resonances move upfield while the heme methyl and meso (13)C resonances move downfield. These changes are a consequence of altered overlap of the Fe 3d and porphyrin pi orbitals, which destabilizes all three occupied Fe 3d-based molecular orbitals and decreases the positive and negative spin density on the beta-pyrrole and meso carbons, respectively. Consequently, the heme ruffling deformation decreases the electronic coupling of the cofactor with external redox partners and lowers the reduction potential of heme.

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.

Vibrational dynamics of iron in cytochrome C

Nuclear resonance vibrational spectroscopy (NRVS) and Raman spectroscopy on (54)Fe- and (57)Fe-enriched cytochrome c (cyt c) identify multiple bands involving vibrations of the heme Fe. Comparison with predictions from Fe isotope shifts reveals that 70% of the NRVS signal in the 300-450 cm(-1) frequency range corresponds to vibrations resolved in Soret-enhanced Raman spectra. This frequency range dominates the "stiffness", an effective force constant determined by the Fe vibrational density of states (VDOS), which measures the strength of nearest-neighbor interactions with Fe. The stiffness of the low-spin Fe environment in both oxidation states of cyt c significantly exceeds that for the high-spin Fe in deoxymyoglobin, where the 200-300 cm(-1) frequency range dominates the VDOS. This situation is reflected in the shorter Fe-ligand bond lengths in the former with respect to the latter. The longer Fe-S(Met80) in oxidized cyt c with respect to reduced cyt c leads to a decrease in the stiffness of the iron environment upon oxidation. Comparison with NRVS measurements allows us to assess assignments for vibrational modes resolved in this region of the heme Raman spectrum. We consider the possibility that the 372 cm(-1) band in reduced cyt c involves the Fe-S(Met80) bond.

Three-pulse photon echo peak shift spectroscopy as a probe of flexibility and conformational heterogeneity in protein folding

We investigate the equilibrium unfolding of Zn-cytochrome c in guanidine hydrochloride by three-pulse photon echo peak shift (3PEPS) spectroscopy. Unexpectedly, the measurements reveal that inhomogeneous broadening of the sample at the midpoint of the denaturation is larger than that of either native or unfolded states. To interpret this finding, we present simulations of the peak shift for both two-state and three-state unfolding models. Both the denaturant concentration dependence of the asymptotic peak shift (APS) and the wavelength dependence of the APS at the midpoint of the denaturation are different for the two models. Our data are consistent with two-state unfolding.

The chemistry and biochemistry of heme c: functional bases for covalent attachment

A discussion of the literature concerning the synthesis, function, and activity of heme c-containing proteins is presented. Comparison of the properties of heme c, which is covalently bound to protein, is made to heme b, which is bound noncovalently. A question of interest is why nature uses biochemically expensive heme c in many proteins when its properties are expected to be similar to heme b. Considering the effects of covalent heme attachment on heme conformation and on the proximal histidine interaction with iron, it is proposed that heme attachment influences both heme reduction potential and ligand-iron interactions.

Low-frequency mode activity of heme: femtosecond coherence spectroscopy of iron porphine halides and nitrophorin

The low-frequency mode activity of metalloporphyrins has been studied for iron porphine-halides (Fe(P)(X), X = Cl, Br) and nitrophorin 4 (NP4) using femtosecond coherence spectroscopy (FCS) in combination with polarized resonance Raman spectroscopy and density functional theory (DFT). It is confirmed that the mode symmetry selection rules for FCS are the same as for Raman scattering and that both Franck-Condon and Jahn-Teller mode activities are observed for Fe(P)(X) under Soret resonance conditions. The DFT-calculated low-frequency (20-400 cm (-1)) modes, and their frequency shifts upon halide substitution, are in good agreement with experimental Raman and coherence data, so that mode assignments can be made. The doming mode is located at approximately 80 cm (-1) for Fe(P)(Cl) and at approximately 60 cm (-1) for Fe(P)(Br). NP4 is also studied with coherence techniques, and the NO-bound species of ferric and ferrous NP4 display a mode at approximately 30-40 cm (-1) that is associated with transient heme doming motion following NO photolysis. The coherence spectra of three ferric derivatives of NP4 with different degrees of heme ruffling distortion are also investigated. We find a mode at approximately 60 cm (-1) whose relative intensity in the coherence spectra depends quadratically on the magnitude of the ruffling distortion. To quantitatively account for this correlation, a new "distortion-induced" Raman enhancement mechanism is presented. This mechanism is unique to low-frequency "soft modes" of the molecular framework that can be distorted by environmental forces. These results demonstrate the potential of FCS as a sensitive probe of dynamic and functionally important nonplanar heme vibrational excitations that are induced by the protein environmental forces or by the chemical reactions in the aqueous phase.

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.