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

Effects of protein structure on iron-polypeptide vibrational dynamic coupling in cytochrome c

Cytochrome c (Cyt c) has a heme covalently bound to the polypeptide via a Cys-X-X-Cys-His (CXXCH) linker that is located in the interface region for protein-protein interactions. To determine whether the polypeptide matrix influences iron vibrational dynamics, nuclear resonance vibrational spectroscopy (NRVS) measurements were performed on (57)Fe-labeled ferric Hydrogenobacter thermophilus cytochrome c-552, and variants M13V, M13V/K22M, and A7F, which have structural modifications that alter the composition or environment of the CXXCH pentapeptide loop. Simulations of the NRVS data indicate that the 150-325 cm(-1) region is dominated by NHis-Fe-SMet axial ligand and polypeptide motions, while the 325-400 cm(-1) region shows dominant contributions from ν(Fe-NPyr) (Pyr = pyrrole) and other heme-based modes. Diagnostic spectral signatures that directly relate to structural features of the heme active site are identified using a quantum chemistry-centered normal coordinate analysis (QCC-NCA). In particular, spectral features that directly correlate with CXXCH loop stiffness, the strength of the Fe-His interaction, and the degree of heme distortion are identified. Cumulative results from our investigation suggest that compared to the wild type (wt), variants M13V and M13V/K22M have a more rigid CXXCH pentapeptide segment, a stronger Fe-NHis interaction, and a more ruffled heme. Conversely, the A7F variant has a more planar heme and a weaker Fe-NHis bond. These results are correlated to the observed changes in reduction potential between wt protein and the variants studied here. Implications of these results for Cyt c biogenesis and electron transfer are also discussed.

How to make a porphyrin flip: dynamics of asymmetric porphyrin oligomers

The pathways and energy barriers from computational and experimental methods reveal the dynamics of asymmetric porphyrin oligomers.

Cytochrome c: A Multifunctional Protein Combining Conformational Rigidity with Flexibility

Cytochrome has served as a model system for studying redox reactions, protein folding, and more recently peroxidase activity induced by partial unfolding on membranes. This review illuminates some important aspects of the research on this biomolecule. The first part summarizes the results of structural analyses of its active site. Owing to heme-protein interactions the heme group is subject to both in-plane and out-of-plane deformations. The unfolding of the protein as discussed in detail in the second part of this review can be induced by changes of pH and temperature and most prominently by the addition of denaturing agents. Both the kinetic and thermodynamic folding and unfolding involve intermediate states with regard to all unfolding conditions. If allowed to sit at alkaline pH (11.5) for a week, the protein does not return to its folding state when the solvent is switched back to neutral pH. It rather adopts a misfolded state that is prone to aggregation via domain swapping. On the surface of cardiolipin containing liposomes, the protein can adopt a variety of partially unfolded states. Apparently, ferricytochrome c can perform biological functions even if it is only partially folded.

Crystallographic and spectroscopic insights into heme degradation by Mycobacterium tuberculosis MhuD

Mycobacterium heme utilization degrader (MhuD) is a heme-degrading protein from Mycobacterium tuberculosis responsible for extracting the essential nutrient iron from host-derived heme. MhuD has been previously shown to produce unique organic products compared to those of canonical heme oxygenases (HOs) as well as those of the IsdG/I heme-degrading enzymes from Staphylococcus aureus. Here, we report the X-ray crystal structure of cyanide-inhibited MhuD (MhuD-heme-CN) as well as detailed (1)H nuclear magnetic resonance (NMR), UV/vis absorption, and magnetic circular dichroism (MCD) spectroscopic characterization of this species. There is no evidence for an ordered network of water molecules on the distal side of the heme substrate in the X-ray crystal structure, as was previously reported for canonical HOs. The degree of heme ruffling in the crystal structure of MhuD is greater than that observed for HO and less than that observed for IsdI. As a consequence, the Fe 3dxz-, 3dyz-, and 3dxy-based MOs are very close in energy, and the room-temperature (1)H NMR spectrum of MhuD-heme-CN is consistent with population of both a (2)Eg electronic state with a (dxy)(2)(dxz,dyz)(3) electron configuration, similar to the ground state of canonical HOs, and a (2)B2g state with a (dxz,dyz)(4)(dxy)(1) electron configuration, similar to the ground state of cyanide-inhibited IsdI. Variable temperature, variable field MCD saturation magnetization data establishes that MhuD-heme-CN has a (2)B2g electronic ground state with a low-lying (2)Eg excited state. Our crystallographic and spectroscopic data suggest that there are both structural and electronic contributions to the α-meso regioselectivity of MhuD-catalyzed heme cleavage. The structural distortion of the heme substrate observed in the X-ray crystal structure of MhuD-heme-CN is likely to favor cleavage at the α- and γ-meso carbons, whereas the spin density distribution may favor selective oxygenation of the α-meso carbon.

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.

Redox state dependence of axial ligand dynamics in Nitrosomonas europaea cytochrome c552

Analysis of NMR spectra reveals that the heme axial Met ligand orientation and dynamics in Nitrosomonas europaea cytochrome c552 (Ne cyt c) are dependent on the heme redox state. In the oxidized state, the heme axial Met is fluxional, interconverting between two conformers related to each other by inversion through the Met δS atom. In the reduced state, there is no evidence of fluxionality, with the Met occupying one conformation similar to that seen in the homologous Pseudomonas aeruginosa cytochrome c551. Comparison of the observed and calculated pseudocontact shifts for oxidized Ne cyt c using the reduced protein structure as a reference structure reveals a redox-dependent change in the structure of the loop bearing the axial Met (loop 3). Analysis of nuclear Overhauser effects (NOEs) and existing structural data provides further support for the redox state dependence of the loop 3 structure. Implications for electron transfer function are discussed.

Structural characterization of nitrosomonas europaea cytochrome c-552 variants with marked differences in electronic structure

Nitrosomonas europaea cytochrome c-552 (Ne c-552) variants with the same His/Met axial ligand set but with different EPR spectra have been characterized structurally, to aid understanding of how molecular structure determines heme electronic structure. Visible light absorption, Raman, and resonance Raman spectroscopy of the protein crystals was performed along with structure determination. The structures solved are those of Ne c-552, which displays a "HALS" (or highly anisotropic low-spin) EPR spectrum, and of the deletion mutant Ne N64Δ, which has a rhombic EPR spectrum. Two X-ray crystal structures of wild-type Ne c-552 are reported; one is of the protein isolated from N. europaea cells (Ne c-552n, 2.35 Å resolution), and the other is of recombinant protein expressed in Escherichia coli (Ne c-552r, 1.63 Å resolution). Ne N64Δ crystallized in two different space groups, and two structures are reported [monoclinic (2.1 Å resolution) and hexagonal (2.3 Å resolution)]. Comparison of the structures of the wild-type and mutant proteins reveals that heme ruffling is increased in the mutant; increased ruffling is predicted to yield a more rhombic EPR spectrum. The 2.35 Å Ne c-552n structure shows 18 molecules in the asymmetric unit; analysis of the structure is consistent with population of more than one axial Met configuration, as seen previously by NMR. Finally, the mutation was shown to yield a more hydrophobic heme pocket and to expel water molecules from near the axial Met. These structures reveal that heme pocket residue 64 plays multiple roles in regulating the axial ligand orientation and the interaction of water with the heme. These results support the hypothesis that more ruffled hemes lead to more rhombic EPR signals in cytochromes c with His/Met axial ligation.

Comparative molecular field analysis and comparative molecular similarity index analysis studies on 1H NMR chemical shift of NH group of diaryl triazene derivatives

Comparative molecular field analysis (CoMFA), comparative molecular field analysis region focusing (CoMFA-RF) for optimizing the region for the final partial least square analysis, and comparative molecular similarity indices analysis (CoMSIA) methods were employed to develop three-dimensional quantitative structure-activity relationship (3D-QSAR) models of (1)H NMR chemical shift of NH proton of diaryl triazene derivatives. The best orientation was searched by all-orientation search (AOS) strategy to minimize the effect of the initial orientation of the structures. The predictive abilities of CoMFA-RF and CoMSIA models were determined using a test set of ten compounds affording predictive correlation coefficients of 0.721 and 0.754, respectively, indicating good predictive power. For further model validation, cross validation (leave one out), progressive scrambling, and bootstrapping were also applied. The accuracy and speed of obtained 3D-QSAR models for the prediction of (1)H NMR chemical shifts of NH group of diaryl triazene derivatives were greater compared to some computational well-known procedures.

Exploring the active site structure of a photoreceptor protein by Raman optical activity

We have developed a near-infrared excited Raman optical activity (ROA) spectrometer and report the first measurement of near-infrared ROA spectra of a light-driven proton pump, bacteriorhodopsin. Our results demonstrate that a near-infrared excitation enables us to measure the ROA spectra of the chromophore within a protein environment. Furthermore, the ROA spectra of the all-trans, 15-anti and 13-cis, 15-syn isomers differ significantly, indicating a high structural sensitivity of the ROA spectra. We therefore expect that future applications of the near-infrared ROA will allow the experimental elucidation of the active site structures in other proteins as well as reaction intermediates.

Conformational change and human cytochrome c function: mutation of residue 41 modulates caspase activation and destabilizes Met-80 coordination

Cytochrome c is a highly conserved protein, with 20 residues identical in all eukaryotic cytochromes c. Gly-41 is one of these invariant residues, and is the position of the only reported naturally occurring mutation in cytochrome c (human G41S). The basis, if any, for the conservation of Gly-41 is unknown. The mutation of Gly-41 to Ser enhances the apoptotic activity of cytochrome c without altering its role in mitochondrial electron transport. Here we have studied additional residue 41 variants and determined their effects on cytochrome c functions and conformation. A G41T mutation decreased the ability of cytochrome c to induce caspase activation and decreased the redox potential, whereas a G41A mutation had no impact on caspase induction but the redox potential increased. All residue 41 variants decreased the pK (a) of a structural transition of oxidized cytochrome c to the alkaline conformation, and this correlated with a destabilization of the interaction of Met-80 with the heme iron(III) at physiological pH. In reduced cytochrome c the G41T and G41S mutations had distinct effects on a network of hydrogen bonds involving Met-80, and in G41T the conformational mobility of two Ω-loops was altered. These results suggest the impact of residue 41 on the conformation of cytochrome c influences its ability to act in both of its physiological roles, electron transport and caspase activation.

Marked difference in the electronic structure of cyanide-ligated ferric protoglobins and myoglobin due to heme ruffling

Electron paramagnetic resonance experiments reveal a significant difference between the principal g values (and hence ligand-field parameters) of the ferric cyanide-ligated form of different variants of the protoglobin of Methanosarcina acetivorans (MaPgb) and of horse heart myoglobin (hhMb). The largest principal g value of the ferric cyanide-ligated MaPgb variants is found to be significantly lower than for any of the other globins reported so far. This is at least partially caused by the strong heme distortions as proven by the determination of the hyperfine interaction of the heme nitrogens and mesoprotons. Furthermore, the experiments confirm recent theoretical predictions [Forti, F.; Boechi, L., Bikiel, D., Martí, M.A.; Nardini, M.; Bolognesi, M.; Viappiani, C.; Estrin, D.; Luque, F. J. J. Phys. Chem. B 2011, 115, 13771-13780] that Phe(G8)145 plays a crucial role in the ligand modulation in MaPgb. Finally, the influence of the N-terminal 20 amino-acid chain on the heme pocket in these protoglobins is also proven.

Steric bulkiness of pyrrole substituents and the out-of-plane deformations of porphyrins: nickel(II) octaisopropylporphyrin and its meso-nitro derivative

Sterically bulky substituents at the β-carbons of the pyrrole rings of porphyrins are sufficient to cause large out-of-plane porphyrin distortions even in the absence of substituent groups at the meso carbons. It is well established that substituents at the meso-positions only or at both the β-pyrrole and the meso-positions are sufficiently bulky to result in large non-planar distortions of the macrocycle. However, no systematic studies of the effects of bulky β-pyrrole substituents alone have been reported. Herein, molecular simulations and X-ray crystallography of nickel(II) 2,3,7,8,12,13,17,18-octa(isopropyl)porphyrin reveal that large out-of-plane distortions (>1.5 Å) are induced by the steric repulsion of the β-isopropyl groups but fail to lead to a single strongly energetically favored conformer. The molecular simulations indicate that multiple conformers differing in the orientation of the isopropyl groups and the macrocycle conformation coexist in solution and this is confirmed by resonance Raman spectroscopy. Large downshifts in the structure-sensitive lines result from the non-planar distortion, and line broadenings result from structural heterogeneity. The heterogeneity originates from tradeoffs between energy contributions from steric repulsion and macrocycle distortion. Simulations for 5-nitro-2,3,7,8,12,13,17,18-octa(isopropyl)porphyrin suggest two possible orientations of the nitro group with respect to the macrocycle mean plane — one nearly vertical (as in the crystal structure) and another that is nearly parallel. INDO/S semiempirical calculations indicate an orbital of the NO2 group resides between the porphyrin frontier orbitals with significant mixing of the nitro and porphyrin orbitals.KEYWORDS: porphyrin, non-planar, resonance Raman, X-ray crystallography, crystal structure, isopropyl, nitro, conformer, molecular mechanics, molecular simulations, density functional theory, steric crowding, conformational heterogeneity.

Nuclear inelastic scattering and Mössbauer spectroscopy as local probes for ligand binding modes and electronic properties in proteins: vibrational behavior of a ferriheme center inside a β-barrel protein

In this work, we present a study of the influence of the protein matrix on its ability to tune the binding of small ligands such as NO, cyanide (CN(-)), and histamine to the ferric heme iron center in the NO-storage and -transport protein Nitrophorin 2 (NP2) from the salivary glands of the blood-sucking insect Rhodnius prolixus. Conventional Mössbauer spectroscopy shows a diamagnetic ground state of the NP2-NO complex and Type I and II electronic ground states of the NP2-CN(-) and NP2-histamine complex, respectively. The change in the vibrational signature of the protein upon ligand binding has been monitored by Nuclear Inelastic Scattering (NIS), also called Nuclear Resonant Vibrational Spectroscopy (NRVS). The NIS data thus obtained have also been calculated by quantum mechanical (QM) density functional theory (DFT) coupled with molecular mechanics (MM) methods. The calculations presented here show that the heme ruffling in NP2 is a consequence of the interaction with the protein matrix. Structure optimizations of the heme and its ligands with DFT retain the characteristic saddling and ruffling only if the protein matrix is taken into account. Furthermore, simulations of the NIS data by QM/MM calculations suggest that the pH dependence of the binding of NO, but not of CN(-) and histamine, might be a consequence of the protonation state of the heme carboxyls.

Revealing substituent effects on the electronic structure and planarity of Ni-porphyrins

Using density functional theory, we have studied the effects on structural and electronic consequences (including HOMO-LUMO energy gaps, vertical ionization potentials (IP_v), and vertical electron affinities (EA_v)) of the following two factors: (a) meso- and β-substituents acting as inductive donors (CH₃), inductive acceptors that are electron-donating through resonance (Br), inductive electron acceptors (CF₃), and resonance enabled acceptors (NO₂); and (b) complete replacement of pyrrole nitrogens with P-atoms. The principal results of the study are: (1) For the bare Ni-porphyrin, the solvents were found not to affect the HOMO-LUMO gaps but to change the IP_v and EA_v noticeably. (2) In the series CH₃ → Br → CF₃ → NO₂ the HOMO-LUMO energy gaps, IP_v, and EA_v increase for both meso- and β-substituents. The ruffling distortion of the porphyrin core is retained, and becomes stronger for the two acceptor groups. In general, effects of meso-substituents on the ruffling distortion of the porphyrin core is more pronounced. (3) Most significantly, complete replacement of pyrrole nitrogens in the NiP with phosphorus atoms produces the species, NiP(P)₄, with the structural and electronic features drastically different from the original NiP. This implies that NiP(P)₄ can possess interesting and unusual novel properties, including aromaticity and reactivity, leading to its various beneficial potential applications. Furthermore, NiP(P)₄ high stability both in the gas phase and different solvents was shown, implying the feasibility of its synthesis.

Modulation of ligand-field parameters by heme ruffling in cytochromes c revealed by EPR spectroscopy

Electron paramagnetic resonance (EPR) spectra of variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht c-552) and Pseudomonas aeruginosa cytochrome c(551) (Pa c-551) are analyzed to determine the effect of heme ruffling on ligand-field parameters. Mutations introduced at positions 13 and 22 in Ht c-552 were previously demonstrated to influence hydrogen bonding in the proximal heme pocket and to tune reduction potential (E(m)) over a range of 80 mV [Michel, L. V.; Ye, T.; Bowman, S. E. J.; Levin, B. D.; Hahn, M. A.; Russell, B. S.; Elliott, S. J.; Bren, K. L. Biochemistry 2007, 46, 11753-11760]. These mutations are shown here to also increase heme ruffling as E(m) decreases. The primary effect on electronic structure of increasing heme ruffling is found to be a decrease in the axial ligand-field term Δ/λ, which is proposed to arise from an increase in the energy of the d(xy) orbital. Mutations at position 7, previously demonstrated to influence heme ruffling in Pa c-551 and Ht c-552, are utilized to test this correlation between molecular and electronic structure. In conclusion, the structure of the proximal heme pocket of cytochromes c is shown to play a role in determining heme conformation and electronic structure.

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.

Comparing substrate specificity between cytochrome c maturation and cytochrome c heme lyase systems for cytochrome c biogenesis

Hemes c are characterized by their covalent attachment to a polypeptide via a widely conserved CXXCH motif. There are multiple biological systems that facilitate heme c biogenesis. System I, the cytochrome c maturation (CCM) system, is found in many bacteria and is commonly employed in the maturation of bacterial cytochromes c in Escherichia coli-based expression systems. System III, cytochrome c heme lyase (CCHL), is an enzyme found in the mitochondria of many eukaryotes and is used for heterologous expression of mitochondrial holocytochromes c. To test CCM specificity, a series of Hydrogenobacter thermophilus cytochrome c(552) variants was successfully expressed and matured by the CCM system with CX(n)CH motifs where n = 1-4, further extending the known substrate flexibility of the CCM system by successful maturation of a bacterial cytochrome c with a novel CXCH motif. Horse cytochrome c variants with both expanded and contracted attachment motifs (n = 1-3) were also tested for expression and maturation by both CCM and CCHL, allowing direct comparison of CCM and CCHL substrate specificities. Successful maturation of horse cytochrome c by CCHL with an extended CXXXCH motif was observed, demonstrating that CCHL shares the ability of CCM to mature hemes c with extended heme attachment motifs. In contrast, two single amino acid mutants were found in horse cytochrome c that severely limit maturation by CCHL, yet were efficiently matured with CCM. These results identify potentially important residues for the substrate recognition of CCHL.

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.

Variation and analysis of second-sphere interactions and axial histidinate character in c-type cytochromes

The electron-donating properties of the axial His ligand to heme iron in cytochromes c (cyts c) are found to be correlated with the midpoint reduction potential (E(m)) in variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht cyt c(552)) in which mutations have been made in and near the Cys-X-X-Cys-His (CXXCH) heme-binding motif. To probe the strength of the His-Fe(III) interaction, we have measured (13)C nuclear magnetic resonance (NMR) chemical shifts for (13)CN(-) bound to heme iron trans to the axial His in Ht Fe(III) cyt c(552) variants. We observe a linear relationship between these (13)C chemical shifts and E(m), indicating that the His-Fe(III) bond strength correlates with E(m). To probe a conserved hydrogen bonding interaction between the axial His Hdelta1 and the backbone carbonyl of a Pro residue, we measured the pK(a) of the axial His Hdelta1 proton (pK(a(2))), which we propose to relate to the His-Fe(III) interaction, reduction potential, and local electrostatic effects. The observed linear relationship between the axial His (13)Cbeta chemical shift and E(m) is proposed to reflect histidinate (anionic) character of the ligand. A linear relationship also is seen between the average heme methyl (1)H chemical shift and E(m) which may reflect variation in axial His electron-donating properties or in the ruffling distortion of the heme plane. In summary, chemical shifts of the axial His and exogenous CN(-) bound trans to His are shown to be sensitive probes of the His-Fe(III) interaction in variants of Ht cyt c(552) and display trends that correlate with E(m).