The structural origin of nonplanar heme distortions in tetraheme ferricytochromes c3

Probing the versatility of cytochrome c by spectroscopic means: A Laudatio on resonance Raman spectroscopy

Over the last 50 years resonance Raman spectroscopy has become an invaluable tool for the exploration of chromophores in biological macromolecules. Among them, heme proteins and metal complexes have attracted considerable attention. This interest results from the fact that resonance Raman spectroscopy probes the vibrational dynamics of these chromophores without direct interference from the surrounding. However, the indirect influence via through-bond and through-space chromophore-protein interactions can be conveniently probed and analyzed. This review article illustrates this point by focusing on class 1 cytochrome c, a comparatively simple heme protein generally known as electron carrier in mitochondria. The article demonstrates how through selective excitation of resonance Raman active modes information about the ligation, the redox state and the spin state of the heme iron can be obtained from band positions in the Raman spectra. The investigation of intensities and depolarization ratios emerged as tools for the analysis of in-plane and out-of-plane deformations of the heme macrocycle. The article further shows how resonance Raman spectroscopy was used to characterize partially unfolded states of oxidized cytochrome c. Finally, it describes its use for exploring structural changes due to the protein's binding to anionic surfaces like cardiolipin containing membranes.

Modern optical approaches in redox biology: Genetically encoded sensors and Raman spectroscopy

The objective of the current review is to summarize the current state of optical methods in redox biology. It consists of two parts, the first is dedicated to genetically encoded fluorescent indicators and the second to Raman spectroscopy. In the first part, we provide a detailed classification of the currently available redox biosensors based on their target analytes. We thoroughly discuss the main architecture types of these proteins, the underlying engineering strategies for their development, the biochemical properties of existing tools and their advantages and disadvantages from a practical point of view. Particular attention is paid to fluorescence lifetime imaging microscopy as a possible readout technique, since it is less prone to certain artifacts than traditional intensiometric measurements. In the second part, the characteristic Raman peaks of the most important redox intermediates are listed, and examples of how this knowledge can be implemented in biological studies are given. This part covers such fields as estimation of the redox states and concentrations of Fe-S clusters, cytochromes, other heme-containing proteins, oxidative derivatives of thiols, lipids, and nucleotides. Finally, we touch on the issue of multiparameter imaging, in which biosensors are combined with other visualization methods for simultaneous assessment of several cellular parameters.

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.

SERS uncovers the link between conformation of cytochrome c heme and mitochondrial membrane potential

The balance between the mitochondrial respiratory chain activity and the cell's needs in ATP ensures optimal cellular function. Cytochrome c is an essential component of the electron transport chain (ETC), which regulates ETC activity, oxygen consumption, ATP synthesis and can initiate apoptosis. The impact of conformational changes in cytochrome c on its function is not understood for the lack of access to these changes in intact mitochondria. We have developed a novel sensor that uses unique properties of label-free surface-enhanced Raman spectroscopy (SERS) to identify conformational changes in heme of cytochrome c and to elucidate their role in functioning mitochondria. We have verified that molecule bond vibrations assessed by SERS are a reliable indicator of the heme conformation during changes in the inner mitochondrial membrane potential and ETC activity. We have demonstrated that cytochrome c heme reversibly switches between planar and ruffled conformations in response to the inner mitochondrial membrane potential (ΔΨ) and H⁺ concentration in the intermembrane space. This regulates the efficiency of the mitochondrial respiratory chain, thus, adjusting the mitochondrial respiration to the cell's consumption of ATP and the overall activity. We have found that under hypertensive conditions cytochrome c heme loses its sensitivity to ΔΨ that can affect the regulation of ETC activity. The ability of the proposed SERS-based sensor to track mitochondrial function opens broad perspectives in cell bioenergetics.

A 300-fold conductivity increase in microbial cytochrome nanowires due to temperature-induced restructuring of hydrogen bonding networks

Although proteins are considered as nonconductors that transfer electrons only up to 1 to 2 nanometers via tunneling, Geobacter sulfurreducens transports respiratory electrons over micrometers, to insoluble acceptors or syntrophic partner cells, via nanowires composed of polymerized cytochrome OmcS. However, the mechanism enabling this long-range conduction is unclear. Here, we demonstrate that individual nanowires exhibit theoretically predicted hopping conductance, at rate (>10¹⁰ s^-1) comparable to synthetic molecular wires, with negligible carrier loss over micrometers. Unexpectedly, nanowires show a 300-fold increase in their intrinsic conductance upon cooling, which vanishes upon deuteration. Computations show that cooling causes a massive rearrangement of hydrogen bonding networks in nanowires. Cooling makes hemes more planar, as revealed by Raman spectroscopy and simulations, and lowers their reduction potential. We find that the protein surrounding the hemes acts as a temperature-sensitive switch that controls charge transport by sensing environmental perturbations. Rational engineering of heme environments could enable systematic tuning of extracellular respiration.

Nature's nitrite-to-ammonia expressway, with no stop at dinitrogen

Since the characterization of cytochrome c₅₅₂ as a multiheme nitrite reductase, research on this enzyme has gained major interest. Today, it is known as pentaheme cytochrome c nitrite reductase (NrfA). Part of the NH₄⁺ produced from NO₂^- is released as NH₃ leading to nitrogen loss, similar to denitrification which generates NO, N₂O, and N₂. NH₄⁺ can also be used for assimilatory purposes, thus NrfA contributes to nitrogen retention. It catalyses the six-electron reduction of NO₂^- to NH₄⁺, hosting four His/His ligated c-type hemes for electron transfer and one structurally differentiated active site heme. Catalysis occurs at the distal side of a Fe(III) heme c proximally coordinated by lysine of a unique CXXCK motif (Sulfurospirillum deleyianum, Wolinella succinogenes) or, presumably, by the canonical histidine in Campylobacter jejeuni. Replacement of Lys by His in NrfA of W. succinogenes led to a significant loss of enzyme activity. NrfA forms homodimers as shown by high resolution X-ray crystallography, and there exist at least two distinct electron transfer systems to the enzyme. In γ-proteobacteria (Escherichia coli) NrfA is linked to the menaquinol pool in the cytoplasmic membrane through a pentaheme electron carrier (NrfB), in δ- and ε-proteobacteria (S. deleyianum, W. succinogenes), the NrfA dimer interacts with a tetraheme cytochrome c (NrfH). Both form a membrane-associated respiratory complex on the extracellular side of the cytoplasmic membrane to optimize electron transfer efficiency. This minireview traces important steps in understanding the nature of pentaheme cytochrome c nitrite reductases, and discusses their structural and functional features.

Protein Engineering of Electron Transfer Components from Electroactive Geobacter Bacteria

Electrogenic microorganisms possess unique redox biological features, being capable of transferring electrons to the cell exterior and converting highly toxic compounds into nonhazardous forms. These microorganisms have led to the development of Microbial Electrochemical Technologies (METs), which include applications in the fields of bioremediation and bioenergy production. The optimization of these technologies involves efforts from several different disciplines, ranging from microbiology to materials science. Geobacter bacteria have served as a model for understanding the mechanisms underlying the phenomenon of extracellular electron transfer, which is highly dependent on a multitude of multiheme cytochromes (MCs). MCs are, therefore, logical targets for rational protein engineering to improve the extracellular electron transfer rates of these bacteria. However, the presence of several heme groups complicates the detailed redox characterization of MCs. In this Review, the main characteristics of electroactive Geobacter bacteria, their potential to develop microbial electrochemical technologies and the main features of MCs are initially highlighted. This is followed by a detailed description of the current methodologies that assist the characterization of the functional redox networks in MCs. Finally, it is discussed how this information can be explored to design optimal Geobacter-mutated strains with improved capabilities in METs.

Structural Insights into Oxygen-Dependent Signal Transduction within Globin Coupled Sensors

In order to respond to external stimuli, bacteria have evolved sensor proteins linking external signals to intracellular outputs that can then regulate downstream pathways and phenotypes. Globin coupled sensor proteins (GCSs) serve to link environmental O₂ levels to cellular processes by coupling a heme-containing sensor globin domain to a catalytic output domain. However, the mechanism by which O₂ binding activates these proteins is currently unknown. To provide insights into the signaling mechanism, two distinct dimeric complexes of the isolated globin domain of the GCS from Bordetella pertussis ( BpeGlobin) were solved via X-ray crystallography in which differences in ligand-bound states were observed. Both monomers of one dimer contain Fe(II)-O₂ states, while the other dimer consists of the Fe(III)-H₂O and Fe(II)-O₂ states. These data provide the first molecular insights into the heme pocket conformation of the active Fe(II)-O₂ form of these enzymes. In addition, heme distortion modes and heme-protein interactions were found to correlate with the ligation state of the globin, suggesting that these conformational changes play a role in O₂-dependent signaling. Fourier transform infrared spectroscopy (FTIR) of the full-length GCS from B. pertussis ( BpeGReg) and the closely related GCS from Pectobacterium carotovorum ssp. carotovorum ( PccGCS) confirmed the importance of an ordered water within the heme pocket and two distal residues (Tyr43 and Ser68) as hydrogen-bond donors. Taken together, this work provides mechanistic insights into BpeGReg O₂ sensing and the signaling mechanisms of diguanylate cyclase-containing GCS proteins.

Influence of heme c attachment on heme conformation and potential

Heme c is characterized by its covalent attachment to a polypeptide. The attachment is typically to a CXXCH motif in which the two Cys form thioether bonds with the heme, "X" can be any amino acid other than Cys, and the His serves as a heme axial ligand. Some cytochromes c, however, contain heme attachment motifs with three or four intervening residues in a CX₃CH or CX₄CH motif. Here, the impacts of these variations in the heme attachment motif on heme ruffling and electronic structure are investigated by spectroscopically characterizing CX₃CH and CX₄CH variants of Hydrogenobacter thermophilus cytochrome c₅₅₂. In addition, a novel CXCH variant is studied. ¹H and ¹³C NMR, EPR, and resonance Raman spectra of the protein variants are analyzed to deduce the extent of ruffling using previously reported relationships between these spectral data and heme ruffling. In addition, the reduction potentials of these protein variants are measured using protein film voltammetry. The CXCH and CX₄CH variants are found to have enhanced heme ruffling and lower reduction potentials. Implications of these results for the use of these noncanonical motifs in nature, and for the engineering of novel heme peptide structures, are discussed.

New insight into the mechanism of mitochondrial cytochrome c function

We investigate functional role of the P⁷⁶GTKMIFA⁸³ fragment of the primary structure of cytochrome c. Based on the data obtained by the analysis of informational structure (ANIS), we propose a model of functioning of cytochrome c. According to this model, conformational rearrangements of the P⁷⁶GTKMIFA⁸³ loop fragment have a significant effect on conformational mobility of the heme. It is suggested that the conformational mobility of cytochrome c heme is responsible for its optimal orientation with respect to electron donor and acceptor within ubiquinol–cytochrome c oxidoreductase (complex III) and cytochrome c oxidase (complex IV), respectively, thus, ensuring electron transfer from complex III to complex IV. To validate the model, we design several mutant variants of horse cytochrome c with multiple substitutions of amino acid residues in the P⁷⁶GTKMIFA⁸³ sequence that reduce its ability to undergo conformational rearrangements. With this, we study the succinate–cytochrome c reductase and cytochrome c oxidase activities of rat liver mitoplasts in the presence of mutant variants of cytochrome c. The electron transport activity of the mutant variants decreases to different extent. Resonance Raman spectroscopy (RRS) and surface-enhanced Raman spectroscopy (SERS) data demonstrate, that all mutant cytochromes possess heme with the higher degree of ruffling deformation, than that of the wild-type (WT) cytochrome c. The increase in the ruffled deformation of the heme of oxidized cytochromes correlated with the decrease in the electron transport rate of ubiquinol–cytochrome c reductase (complex III). Besides, all mutant cytochromes have lower mobility of the pyrrol rings and methine bridges, than WT cytochrome c. We show that a decrease in electron transport activity in the mutant variants correlates with conformational changes and reduced mobility of heme porphyrin. This points to a significant role of the P⁷⁶GTKMIFA⁸³ fragment in the electron transport function of cytochrome c.

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.

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.

Resonance Raman spectroscopy reveals that substrate structure selectively impacts the heme-bound diatomic ligands of CYP17

An important function of steroidogenic cytochromes P450 is the transformation of cholesterol to produce androgens, estrogens, and the corticosteroids. The activities of cytochrome P450c17 (CYP17) are essential in sex hormone biosynthesis, with severe developmental defects being a consequence of deficiency or mutations. The first reaction catalyzed by this multifunctional P450 is the 17α-hydroxylation of pregnenolone (PREG) to 17α-hydroxypregnenolone (17-OH PREG) and progesterone (PROG) to 17α-hydroxyprogesterone (17-OH PROG). The hydroxylated products then either are used for production of corticoids or undergo a second CYP17 catalyzed transformation, representing the first committed step of androgen formation. While the hydroxylation reactions are catalyzed by the well-known Compound I intermediate, the lyase reaction is believed to involve nucleophilic attack of the earlier peroxo- intermediate on the C20-carbonyl. Herein, resonance Raman (rR) spectroscopy reveals that substrate structure does not impact heme structure for this set of physiologically important substrates. On the other hand, rR spectra obtained here for the ferrous CO adducts with these four substrates show that substrates do interact differently with the Fe-C-O fragment, with large differences between the spectra obtained for the samples containing 17-OH PROG and 17-OH PREG, the latter providing evidence for the presence of two Fe-C-O conformers. Collectively, these results demonstrate that individual substrates can differentially impact the disposition of a heme-bound ligand, including dioxygen, altering the reactivity patterns in such a way as to promote preferred chemical conversions, thereby avoiding the profound functional consequences of unwanted side reactions.

Charge parametrization of the DvH-c3 heme group: validation using constant-(pH,E) molecular dynamics simulations

We studied the effect of using different heme group charge parametrization methods and schemes (Merz-Kollman, CHELPG, and single- and multiconformational RESP) on the quality of the results produced by the constant-(pH,E) MD method, applied to the redox titration of Desulfovibrio vulgaris Hildenborough cytochrome c(3). These new and more accurate charge sets enabled us to overcome the previously reported dependence of the method's performance on the dielectric constant, ε, assigned to the protein region. In particular, we found a systematic, clear shift of the E(mod) toward more negative values than those previously reported, in agreement with an electrostatics based reasoning. The simulations showed strong coupling between protonating/redox sites. We were also able to capture significant direct and, especially, indirect interactions between hemes, such as those mediated by histidine 67. Our results highlight the importance of having a good quantum description of the system before deriving atomic partial charges for classic force fields.

Optimized end-stacking provides specificity of N-methyl mesoporphyrin IX for human telomeric G-quadruplex DNA

N-methyl mesoporphyrin IX (NMM) is exceptionally selective for G-quadruplexes (GQ) relative to duplex DNA and, as such, has found a wide range of applications in biology and chemistry. In addition, NMM is selective for parallel versus antiparallel GQ folds, as was recently demonstrated in our laboratory. Here, we present the X-ray crystal structure of a complex between NMM and human telomeric DNA dAGGG(TTAGGG)(3), Tel22, determined in two space groups, P2(1)2(1)2 and P6, at 1.65 and 2.15 Å resolution, respectively. The former is the highest resolution structure of the human telomeric GQ DNA reported to date. The biological unit contains a Tel22 dimer of 5'-5' stacked parallel-stranded quadruplexes capped on both ends with NMM, supporting the spectroscopically determined 1:1 stoichiometry. NMM is capable of adjusting its macrocycle geometry to closely match that of the terminal G-tetrad required for efficient π-π stacking. The out-of-plane N-methyl group of NMM fits perfectly into the center of the parallel GQ core where it aligns with potassium ions. In contrast, the interaction of the N-methyl group with duplex DNA or antiparallel GQ would lead to steric clashes that prevent NMM from binding to these structures, thus explaining its unique selectivity. On the basis of the biochemical data, binding of NMM to Tel22 does not rely on relatively nonspecific electrostatic interactions, which characterize most canonical GQ ligands, but rather it is hydrophobic in nature. The structural features observed in the NMM-Tel22 complex described here will serve as guidelines for developing new quadruplex ligands that have excellent affinity and precisely defined selectivity.

Inactivation of the heme degrading enzyme IsdI by an active site substitution that diminishes heme ruffling

IsdG and IsdI are paralogous heme degrading enzymes from the bacterium Staphylococcus aureus. Heme bound by these enzymes is extensively ruffled such that the meso-carbons at the sites of oxidation are distorted toward bound oxygen. In contrast, the canonical heme oxygenase family degrades heme that is bound with minimal distortion. Trp-66 is a conserved heme pocket residue in IsdI implicated in heme ruffling. IsdI variants with Trp-66 replaced with residues having less bulky aromatic and alkyl side chains were characterized with respect to catalytic activity, heme ruffling, and electrochemical properties. The heme degradation activity of the W66Y and W66F variants was approximately half that of the wild-type enzyme, whereas the W66L and W66A variants were inactive. A crystal structure and NMR spectroscopic analysis of the W66Y variant reveals that heme binds to this enzyme with less heme ruffling than observed for wild-type IsdI. The reduction potential of this variant (-96 ± 7 mV versus standard hydrogen electrode) is similar to that of wild-type IsdI (-89 ± 7 mV), so we attribute the diminished activity of this variant to the diminished heme ruffling observed for heme bound to this enzyme and conclude that Trp-66 is required for optimal catalytic activity.

Molecular simulations of porphyrins and heme proteins

An overview of the use of classical mechanical molecular simulations of porphyrins, hydroporphyrins and heme proteins is given. The topics cover molecular mechanics calculations of structures and conformer energies of porphyrins, energies of barriers for interconversion between stable conformers, molecular dynamics of porphyrins and heme proteins, and normal-coordinate structural analysis of experimental and calculated porphyrin structures. Molecular mechanics and dynamics are currently a fertile area of research on porphyrins. In the future, other computational methods such as Monte Carlo simulations, which have yet to be applied to porphyrins, will come into use and open new avenues of research into molecular simulations of porphyrins.

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.

Using spectroscopic tools to probe porphyrin deformation and porphyrin-protein interactions

The reactivity and functionality of heme proteins are to a significant extent determined by the conformation of their functional heme groups and by the interaction of axial ligands with their protein environment. This review focuses on experimental methods and theoretical concepts for elucidating symmetry lowering perturbations of the heme induced by the protein environment of the heme pocket. First, we discuss a variety of methods which can be used to probe the electric field at the heme, including spectral hole burning as well as low temperature absorption and room temperature circular dichroism spectroscopy. Second, we show how heme deformations can be described as superposition of deformations along normal coordinates, thereby using the irreducible representations of the D4h point group as a classification tool. Finally, resonance Raman spectroscopy is introduced as a tool to probe the deformations of metalloprophyrins in solution and in protein matrices by measuring and comparing intensities and depolarization properties rather than wavenumber positions.

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.

Electronic properties of the highly ruffled heme bound to the heme degrading enzyme IsdI

IsdI, a heme-degrading protein from Staphylococcus aureus, binds heme in a manner that distorts the normally planar heme prosthetic group to an extent greater than that observed so far for any other heme-binding protein. To understand better the relationship between this distinct structural characteristic and the functional properties of IsdI, spectroscopic, electrochemical, and crystallographic results are reported that provide evidence that this heme ruffling is essential to the catalytic activity of the protein and eliminates the need for the water cluster in the distal heme pocket that is essential for the activity of classical heme oxygenases. The lack of heme orientational disorder in (1)H-NMR spectra of the protein argues that the catalytic formation of β- and δ-biliverdin in nearly equal yield results from the ability of the protein to attack opposite sides of the heme ring rather than from binding of the heme substrate in two alternative orientations.

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.

Soluble guanylate cyclase is activated differently by excess NO and by YC-1: resonance Raman spectroscopic evidence

Modulation of soluble guanylate cyclase (sGC) activity by nitric oxide (NO) involves two distinct steps. Low-level activation of sGC is achieved by the stoichiometric binding of NO (1-NO) to the heme cofactor, while much higher activation is achieved by the binding of additional NO (xsNO) at a non-heme site. Addition of the allosteric activator YC-1 to the 1-NO form leads to activity comparable to that of the xsNO state. In this study, the mechanisms of sGC activation were investigated using electronic absorption and resonance Raman (RR) spectroscopic methods. RR spectroscopy confirmed that the 1-NO form contains five-coordinate NO-heme and showed that the addition of NO to the 1-NO form has no significant effect on the spectrum. In contrast, addition of YC-1 to either the 1-NO or xsNO forms alters the RR spectrum significantly, indicating a protein-induced change in the heme geometry. This change in the heme geometry was also observed when BAY 41-2272 was added to the xsNO form. Bands assigned to bending and stretching motions of the vinyl and propionate substituents undergo changes in intensity in a pattern suggesting altered tilting of the pyrrole rings to which they are attached. In addition, the N-O stretching frequency increases, with no change in the Fe-NO stretching frequency, an effect modeled via DFT calculations as resulting from a small opening of the Fe-N-O angle. These spectral differences demonstrate different mechanisms of activation by synthetic activators, such as YC-1 and BAY 41-2272, and excess NO.

Probing the function of heme distortion in the H-NOX family

Hemoproteins carry out diverse functions utilizing a wide range of chemical reactivity while employing the same heme prosthetic group. It is clear from high-resolution crystal structures and biochemical studies that protein-bound hemes are not planar and adopt diverse conformations. The crystal structure of an H-NOX domain from Thermoanaerobacter tengcongensis (Tt H-NOX) contains the most distorted heme reported to date. In this study, Tt H-NOX was engineered to adopt a flatter heme by mutating proline 115, a conserved residue in the H-NOX family, to alanine. Decreasing heme distortion in Tt H-NOX increases affinity for oxygen and decreases the reduction potential of the heme iron. Additionally, flattening the heme is associated with significant shifts in the N-terminus of the protein. These results show a clear link between the heme conformation and Tt H-NOX structure and demonstrate that heme distortion is an important determinant for maintaining biochemical properties in H-NOX proteins.

Bis-histidine-coordinated hemes in four-helix bundles: how the geometry of the bundle controls the axial imidazole plane orientations in transmembrane cytochromes of mitochondrial complexes II and III and related proteins

Early investigation of the electron paramagnetic resonance spectra of bis-histidine-coordinated membrane-bound ferriheme proteins led to the description of a spectral signal that had only one resolved feature. These became known as "highly anisotropic low-spin" or "large g(max)" ferriheme centers. Extensive work with small-molecule model heme complexes showed that this spectroscopic signature occurs in bis-imidazole ferrihemes in which the planes of the imidazole ligands are nearly perpendicular, deltaphi = 57-90 degrees. In the last decade protein crystallographic studies have revealed the atomic structures of a number of examples of bis-histidine heme proteins. A frequent characteristic of these large g(max) ferrihemes in membrane-bound proteins is the occurrence of the heme within a four-helix bundle with a left-handed twist. The histidine ligands occur at the same level on two diametrically opposed helices of the bundle. These ligands have the same side-chain conformation and ligate heme iron on the bundle axis, resulting in a quasi-twofold symmetric structure. The two non-ligand-bearing helices also obey this symmetry, and have a conserved small residue, usually glycine, where the edge of the heme ring makes contact with the helix backbones. In many cases this small residue is preceded by a threonine or serine residue whose side-chain hydroxyl oxygen acts as a hydrogen-bond acceptor from the N(delta1) atom of the heme-ligating histidine. The deltaphi angle is thus determined by the common histidine side-chain conformation and the crossing angle of the ligand-bearing helices, in some cases constrained by hydrogen bonds to the serine/threonine residues on the non-ligand-bearing helices.

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.

Ruffling of metalloporphyrins bound to IsdG and IsdI, two heme-degrading enzymes in Staphylococcus aureus

IsdG and IsdI are paralogous proteins that are intracellular components of a complex heme uptake system in Staphylococcus aureus. IsdG and IsdI were shown previously to reductively degrade hemin. Crystal structures of the apoproteins show that these proteins belong to a newly identified heme degradation family distinct from canonical eukaryotic and prokaryotic heme oxygenases. Here we report the crystal structures of an inactive N7A variant of IsdG in complex with Fe(3+)-protoporphyrin IX (IsdG-hemin) and of IsdI in complex with cobalt protoporphyrin IX (IsdI-CoPPIX) to 1.8 A or better resolution. These structures show that the metalloporphyrins are buried into similar deep clefts such that the propionic acids form salt bridges to two Arg residues. His(77) (IsdG) or His(76) (IsdI), a critical residue required for activity, is coordinated to the Fe(3+) or Co(3+) atoms, respectively. The bound porphyrin rings form extensive steric interactions in the binding cleft such that the rings are highly distorted from the plane. This distortion is best described as ruffled and places the beta- and delta-meso carbons proximal to the distal oxygen-binding site. In the IsdG-hemin structure, Fe(3+) is pentacoordinate, and the distal side is occluded by the side chain of Ile(55). However, in the structure of IsdI-CoPPIX, the distal side of the CoPPIX accommodates a chloride ion in a cavity formed through a conformational change in Ile(55). The chloride ion participates in a hydrogen bond to the side chain amide of Asn(6). Together the structures suggest a reaction mechanism in which a reactive peroxide intermediate proceeds with nucleophilic oxidation at the beta- or delta-meso carbon of the hemin.

Modulation of redox potential and alteration in reactivity via the peroxide shunt pathway by mutation of cytochrome P450 around the proximal heme ligand

In the thermophilic cytochrome P450 from the thermoacidophilic crenarchaeon Sulfolobus tokodaii strain 7 (P450st), a phenylalanine residue at position 310 and an alanine residue at position 320 are located close to the heme thiolate ligand, Cys317. Single site-directed mutants F310A and A320Q and double mutant F310A/A320Q have been constructed. All mutant enzymes as well as wild-type (WT) P450st were expressed at high levels. The substitution of F310 with Ala and of A320 with Gln induced shifts in redox potential and blue shifts in Soret absorption of ferrous-CO forms, while spectral characterization showed that in the resting state, the mutants almost retained the structural integrity of the active site. The redox potential of the heme varied as follows: -481 mV (WT), -477 mV (A320Q), -453 mV (F310A), and -450 mV (F310A/A320Q). The trend in the Soret band of the ferrous-CO form was as follows: 450 nm (WT) < 449 nm (A320Q) < 446 nm (F310A) < 444 nm (F310A/A320Q). These results established that the reduction potential and electron density on the heme iron are modulated by the Phe310 and Ala320 residues in P450st. The electron density on the heme decreases in the following order: WT > A320Q > F310A > F310A/A320Q. The electron density on the heme iron infers an essential role in P450 activity. The decrease in electron density interferes with the formation of a high-valent oxo-ferryl species called Compound I. However, steady-state turnover rates of styrene epoxidation with H2O2 show the following trend: WT approximately equal to A320Q < F310A approximately equal to F310A/A320Q. The shunt pathway which can provide the two electrons and oxygen required for a P450 reaction instead of NAD(P)H and dioxygen can rule out the first and second heme reduction in the catalytic process. Because the electron density on the heme iron might be deeply involved in the k cat values in this system, the intermediate Compound 0 which is the precursor species of Compound I mainly appears to participate dominantly in epoxidation with H2O2.

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

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

Heme attachment motif mobility tunes cytochrome c redox potential

Hydrogen exchange (HX) rates and midpoint potentials (Em) of variants of cytochrome c from Pseudomonas aeruginosa (Pa cyt c551) and Hydrogenobacter thermophilus (Ht cyt c552) have been characterized in an effort to develop an understanding of the impact of properties of the Cys-X-X-Cys-His pentapeptide c-heme attachment (CXXCH) motif on heme redox potential. Despite structural conservation of the CXXCH motif, Ht cyt c552 exhibits a low level of protection from HX for amide protons within this motif relative to Pa cyt c551. Site-directed mutants have been prepared to determine the structural basis for and functional implications of these variations on HX behavior. The double mutant Ht-M13V/K22M displays suppressed HX within the CXXCH motif as well as a decreased Em (by 81 mV), whereas the corresponding double mutant of Pa cyt c551 (V13M/M22K) exhibits enhanced HX within the CXXCH pentapeptide and a modest increase in Em (by 30 mV). The changes in Em correlate with changes in axial His chemical shifts in the ferric proteins reflecting the extent of histidinate character. Thus, the mobility of the CXXCH pentapeptide is found to impact the His-Fe(III) interaction and therefore the heme redox potential.

Resonance Raman fingerprinting of multiheme cytochromes from the cytochrome c 3 family

Resonance Raman (RR) spectroscopy was used to investigate conformational characteristics of the hemes of several ferricytochromes of the cytochrome c3 family, electron transfer proteins isolated from the periplasm and membranes of sulfate-reducing bacteria. Our analysis concentrated on the low-frequency region of the RR spectra, a fingerprint region that includes vibrations for heme-protein C-S bonds [nu(C(a)S)]. It has been proposed that these bonds are directly involved in the electron transfer process. The three groups of tetraheme cytochrome c3 analyzed, namely Type I cytochrome c (3) (TpIc (3)s), Type II cytochrome c (3) (TpIIc (3)s) and Desulfomicrobium cytochromes c3, display different frequency separations for the two nu(C(a)S) lines that are similar among members of each group. These spectral differences correlate with differences in protein structure observed among the three groups of cytochromes c3. Two larger cytochromes of the cytochrome c3 family display RR spectral characteristics for the nu(C(a)S) lines that are closer to TpIIc3 than to TpIc3. Two other multiheme cytochromes from Desulfovibrio that do not belong to the cytochrome c3 family display nu(C(a)S) lines with reverse relative areas in comparison with the latter family. This RR study shows that the small differences in protein structure observed among these cytochrome c3 correlate to differences on the heme-protein bonds, which are likely to have an impact upon the protein function, making RR spectroscopy a sensitive and useful tool for characterizing these cytochromes.

Heme carbonyls: environmental effects on nu(C-O) and Fe-C/C-O bond length correlations

The synthesis and characterization of four low-spin (carbonyl)iron(II) tetraphenylporphyrinates, [Fe(TPP)(CO)(L)], where L = 1-methylimidazole, 2-methylimidazole, 1,2-dimethylimidazole (unsolvated), and 1,2-dimethylimidazole (toluene solvate) are reported. The complexes show nearly the same value of nu(C-O) in toluene solution (1969-72 cm(-1)) but a large range of CO stretching frequencies in the solid-state (1926-1968 cm(-1)). The large solid-state variation results from CO interactions in the solid state, as shown by an examination of the crystal structures of the four complexes. The high precision of the four structures obtained allows us to make a number of structural and spectroscopic correlations that describe the Fe-C-O and N(Im)-Fe-CO units. The values of nu(C-O) and the Fe-C and C-O bond distances are strongly correlated and provide a structural, as well as a spectroscopic, correlation of the pi back-bonding model. The interactions of CO described are closely related to the large range of CO stretching frequencies observed in heme proteins and specific interactions observed in carbonylmyoglobin (MbCO).

Absorption and resonance Raman investigations of ligand rotation and nonplanar heme distortion in bis-base low-spin iron(II)-tetrakis(o-pivalamidophenyl)porphyrin complexes

The absorption and resonance Raman (RR) spectra of the bis-N-methylimidazole, bis-1,5-dicyclohexylimidazole, and bis-pyridine complexes of the meso-alphaalphabetabeta and meso-alphabetaalphabeta atropisomers of Fe(II)-tetrakis(o-pivalamidophenyl)porphyrins (Fe(II)TpivPP) were obtained in methylene chloride. The different spatial arrangements of the o-pivalamide pickets in these two Fe(II)TpivPP compounds are expected to control the absolute and relative positions of the axial ligand rings with respect to the Fe-N(pyrrole) bonds. In particular, the spectroscopic data obtained for the bis-N-methylimidazole and bis-dicyclohexylimidazole complexes of the Fe(II)[alphabetaalphabeta-TpivPP] derivative showed the most important differences. Redshifts of the B and Q absorption bands (+ 4-5 nm) as well as an upshift of the low frequency nu(8) RR mode (+ 5 cm(-)(1)) were observed. No shift of the skeletal high frequency modes was detected. These spectral effects were associated with a change in relative position of the axial imidazole rings from nearly parallel in the bis-N-methylimidazole complex to nearly perpendicular in the bis-dicyclohexylimidazole complex. On the basis of stereochemical considerations as well as previous spectroscopic investigations, the data were interpreted in terms of change in porphyrin structure from planar to saddled. Complementing to a parallel study on bis-base Fe(II) "basket handle" porphyrin complexes, this spectroscopic investigation provides an additional means to distinguish planar, ruffled, and saddled conformations for ferrous hemes included in proteins.

Origin of the red shifts in the optical absorption bands of nonplanar tetraalkylporphyrins

The view that the large red shifts seen in the UV-visible absorption bands of peripherally crowded nonplanar porphyrins are the result of nonplanar deformations of the macrocycle has recently been challenged by the suggestion that the red shifts arise from substituent-induced changes in the macrocycle bond lengths and bond angles, termed in-plane nuclear reorganization (IPNR). We have analyzed the contributions to the UV-visible band shifts in a series of nickel or zinc meso-tetraalkylporphyrins to establish the origins of the red shifts in these ruffled porphyrins. Structures were obtained using a molecular mechanics force field optimized for porphyrins, and the nonplanar deformations were quantified by using normal-coordinate structural decomposition (NSD). Transition energies were calculated by the INDO/S semiempirical method. These computational studies demonstrate conclusively that the large Soret band red shifts ( approximately 40 nm) seen for very nonplanar meso-tetra(tert-butyl)porphyrin compared to meso-tetra(methyl)porphyrin are primarily the result of nonplanar deformations and not IPNR. Strikingly, nonplanar deformations along the high-frequency 2B(1u) and 3B(1u) normal coordinates of the macrocycle are shown to contribute significantly to the observed red shifts, even though these deformations are an order of magnitude smaller than the observed ruffling (1B(1u)) deformation. Other structural and electronic influences on the UV-visible band shifts are discussed and problems with the recent studies are examined (e.g., the systematic underestimation of the 2B(1u) and 3B(1u) modes in artificially constrained porphyrin structures that leads to a mistaken attribution of the red shift to IPNR). The effect of nonplanar deformations on the UV-visible absorption bands is then probed experimentally with a series of novel bridled nickel chiroporphyrins. In these compounds, the substituent effect is essentially invariant and the amount of nonplanar deformation decreases as the length of the straps connecting adjacent meso-cyclopropyl substituents decreases (the opposite of the effect observed for conventional strapped porphyrins). Several spectroscopic markers for nonplanarity (UV-visible bands, resonance Raman lines, and (1)H NMR resonances) are found to correlate with time-averaged deformations obtained from an NSD analysis of molecular dynamics snapshot structures. These results suggest that UV-visible band shifts of tetrapyrroles in proteins are potentially useful indicators of changes in nonplanarity provided other structural and electronic factors can be eliminated.