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

Machine Learning-Enhanced Quantum Chemistry-Assisted Refinement of the Active Site Structure of Metalloproteins

Understanding the fine structural details of inhibitor binding at the active site of metalloenzymes can have a profound impact on the rational drug design targeted to this broad class of biomolecules. Structural techniques such as NMR, cryo-EM, and X-ray crystallography can provide bond lengths and angles, but the uncertainties in these measurements can be as large as the range of values that have been observed for these quantities in all the published structures. This uncertainty is far too large to allow for reliable calculations at the quantum chemical (QC) levels for developing precise structure-activity relationships or for improving the energetic considerations in protein-inhibitor studies. Therefore, the need arises to rely upon computational methods to refine the active site structures well beyond the resolution obtained with routine application of structural methods. In a recent paper, we have shown that it is possible to refine the active site of cobalt(II)-substituted MMP12, a metalloprotein that is a relevant drug target, by matching to the experimental pseudocontact shifts (PCS) those calculated using multireference ab initio QC methods. The computational cost of this methodology becomes a significant bottleneck when the starting structure is not sufficiently close to the final one, which is often the case with biomolecular structures. To tackle this problem, we have developed an approach based on a neural network (NN) and a support vector regression (SVR) and applied it to the refinement of the active site structure of oxalate-inhibited human carbonic anhydrase 2 (hCAII), another prototypical metalloprotein target. The refined structure gives a remarkably good agreement between the QC-calculated and the experimental PCS. This study not only contributes to the knowledge of CAII but also demonstrates the utility of combining machine learning (ML) algorithms with QC calculations, offering a promising avenue for investigating other drug targets and complex biological systems in general.

Cytochrome P460 Cofactor Maturation Proceeds via Peroxide-Dependent Post-translational Modification

Cytochrome P460s are heme enzymes that oxidize hydroxylamine to nitrous oxide. They bear specialized "heme P460" cofactors that are cross-linked to their host polypeptides by a post-translationally modified lysine residue. Wild-type N. europaea cytochrome P460 may be isolated as a cross-link-deficient proenzyme following anaerobic overexpression in E. coli. When treated with peroxide, this proenzyme undergoes maturation to active enzyme with spectroscopic and catalytic properties that match wild-type cyt P460. This maturation reactivity requires no chaperones─it is intrinsic to the protein. This behavior extends to the broader cytochrome c'β superfamily. Accumulated data reveal key contributions from the secondary coordination sphere that enable selective, complete maturation. Spectroscopic data support the intermediacy of a ferryl species along the maturation pathway.

Heme protonation affects iron-NO binding in the NO transport protein nitrophorin

The nitrophorins (NPs) comprise an unusual group of heme proteins with stable ferric heme iron nitric oxide (Fe-NO) complexes. They are found in the salivary glands of the blood-sucking kissing bug Rhodnius prolixus, which uses the NPs to transport the highly reactive signaling molecule NO. Nuclear resonance vibrational spectroscopy (NRVS) of both isoform NP2 and a mutant NP2(Leu132Val) show, after addition of NO, a strong structured vibrational band at around 600 cm^-1, which is due to modes with significant Fe-NO bending and stretching contribution. Based on a hybrid calculation method, which uses density functional theory and molecular mechanics, it is demonstrated that protonation of the heme carboxyl groups does influence both the vibrational properties of the Fe-NO entity and its electronic ground state. Moreover, heme protonation causes a significant increase of the gap between the highest occupied and lowest unoccupied molecular orbital by almost one order of magnitude leading to a stabilization of the Fe-NO bond.

Ligand-Based Regulation of Dynamics and Reactivity of Hemoproteins

Hemoproteins include several heme-binding proteins with distinct structure and function. The presence of the heme group confers specific reactivity and spectroscopic properties to hemoproteins. In this review, we provide an overview of five families of hemoproteins in terms of dynamics and reactivity. First, we describe how ligands modulate cooperativity and reactivity in globins, such as myoglobin and hemoglobin. Second, we move on to another family of hemoproteins devoted to electron transport, such as cytochromes. Later, we consider heme-based reactivity in hemopexin, the main heme-scavenging protein. Then, we focus on heme-albumin, a chronosteric hemoprotein with peculiar spectroscopic and enzymatic properties. Eventually, we analyze the reactivity and dynamics of the most recently discovered family of hemoproteins, i.e., nitrobindins.

Heme-Protein Interactions and Functional Relevant Heme Deformations: The Cytochrome c Case

Heme proteins are known to perform a plethora of biologically important functions. This article reviews work that has been conducted on various class I cytochrome c proteins over a period of nearly 50 years. The article focuses on the relevance of symmetry-lowering heme-protein interactions that affect the function of the electron transfer protein cytochrome c. The article provides an overview of various, mostly spectroscopic studies that explored the electronic structure of the heme group in these proteins and how it is affected by symmetry-lowering deformations. In addition to discussing a large variety of spectroscopic studies, the article provides a theoretical framework that should enable a comprehensive understanding of the physical chemistry that underlies the function not only of cytochrome c but of all heme proteins.

Electric field-induced functional changes in electrode-immobilized mutant species of human cytochrome c

Post-translational modifications and naturally occurring mutations of cytochrome c have been recognized as a regulatory mechanism to control its biology. In this work, we investigate the effect of such in vivo chemical modifications of human cytochrome c on its redox properties in the adsorbed state onto an electrode. In particular, tyrosines 48 and 97 have been replaced by the non-canonical amino acid p-carboxymethyl-L-phenylalanine (pCMF), thus mimicking tyrosine phosphorylation. Additionally, tyrosine 48 has been replaced by a histidine producing the natural Y48H pathogenic mutant. Thermodynamics and kinetics of the interfacial electron transfer of wild-type cytochrome c and herein produced variants, adsorbed electrostatically under different local interfacial electric fields, were determined by means of variable temperature cyclic film voltammetry. It is shown that non-native cytochrome c variants immobilized under a low interfacial electric field display redox thermodynamics and kinetics similar to those of wild-type cytochrome c. However, upon increasing the strength of the electric field, the redox thermodynamics and kinetics of the modified proteins markedly differ from those of the wild-type species. The mutations promote stabilization of the oxidized form and a significant increase in the activation enthalpy values that can be ascribed to a subtle distortion of the heme cofactor and/or difference of the amino acid rearrangements rather than to a coarse protein structural change. Overall, these results point to a combined effect of the single point mutations at positions 48 and 97 and the strength of electrostatic binding on the regulatory mechanism of mitochondrial membrane activity, when acting as a redox shuttle protein.

Bioinorganic Chemistry on Electrodes: Methods to Functional Modeling

One of the major goals of bioinorganic chemistry has been to mimic the function of elegant metalloenzymes. Such functional modeling has been difficult to attain in solution, in particular, for reactions that require multiple protons and multiple electrons (nH⁺/ne^-). Using a combination of heterogeneous electrochemistry, electrode and molecule design one may control both electron transfer (ET) and proton transfer (PT) of these nH⁺/ne^- reactions. Such control can allow functional modeling of hydrogenases (H⁺ + e^- → 1/2 H₂), cytochrome c oxidase (O₂ + 4 e^- + 4 H⁺ → 2 H₂O), monooxygenases (RR'CH₂ + O₂ + 2 e^- + 2 H⁺ → RR'CHOH + H₂O) and dioxygenases (S + O₂ → SO₂; S = organic substrate) in aqueous medium and at room temperatures. In addition, these heterogeneous constructs allow probing unnatural bioinspired reactions and estimation of the inner- and outer-sphere reorganization energy of small molecules and proteins.

Molecular Structure of Nickel Octamethylporphyrin—Rare Experimental Evidence of a Ruffling Effect in Gas Phase

The structure of a free nickel (II) octamethylporphyrin (NiOMP) molecule was determined for the first time through a combined gas-phase electron diffraction (GED) and mass spectrometry (MS) experiment, as well as through quantum chemical (QC) calculations. Density functional theory (DFT) calculations do not provide an unambiguous answer about the planarity or non-planar distortion of the NiOMP skeleton. The GED refinement in such cases is non-trivial. Several approaches to the inverse problem solution were used. The obtained results allow us to argue that the ruffling effect is manifested in the NiOMP molecule. The minimal critical distance between the central atom of the metal and nitrogen atoms of the coordination cavity that provokes ruffling distortion in metal porphyrins is about 1.96 Å.

The Biologically Relevant Coordination Chemistry of Iron and Nitric Oxide: Electronic Structure and Reactivity

Nitric oxide (NO) is an important signaling molecule that is involved in a wide range of physiological and pathological events in biology. Metal coordination chemistry, especially with iron, is at the heart of many biological transformations involving NO. A series of heme proteins, nitric oxide synthases (NOS), soluble guanylate cyclase (sGC), and nitrophorins, are responsible for the biosynthesis, sensing, and transport of NO. Alternatively, NO can be generated from nitrite by heme- and copper-containing nitrite reductases (NIRs). The NO-bearing small molecules such as nitrosothiols and dinitrosyl iron complexes (DNICs) can serve as an alternative vehicle for NO storage and transport. Once NO is formed, the rich reaction chemistry of NO leads to a wide variety of biological activities including reduction of NO by heme or non-heme iron-containing NO reductases and protein post-translational modifications by DNICs. Much of our understanding of the reactivity of metal sites in biology with NO and the mechanisms of these transformations has come from the elucidation of the geometric and electronic structures and chemical reactivity of synthetic model systems, in synergy with biochemical and biophysical studies on the relevant proteins themselves. This review focuses on recent advancements from studies on proteins and model complexes that not only have improved our understanding of the biological roles of NO but also have provided foundations for biomedical research and for bio-inspired catalyst design in energy science.

Comparing Properties of Common Bioinorganic Ligands with Switchable Variants of Cytochrome c

Ligand substitution at the metal center is common in catalysis and signal transduction of metalloproteins. Understanding the effects of particular ligands, as well as the polypeptide surrounding, is critical for uncovering mechanisms of these biological processes and exploiting them in the design of bioinspired catalysts and molecular devices. A series of switchable K79G/M80X/F82C (X = Met, His, or Lys) variants of cytochrome (cyt) c was employed to directly compare the stability of differently ligated proteins and activation barriers for Met, His, and Lys replacement at the ferric heme iron. Studies of these variants and their nonswitchable counterparts K79G/M80X have revealed stability trends Met < Lys < His and Lys < His < Met for the protein Fe^III-X and Fe^II-X species, respectively. The differences in the hydrogen-bonding interactions in folded proteins and in solvation of unbound X in the unfolded proteins explain these trends. Calculations of free energy of ligand dissociation in small heme model complexes reveal that the ease of the Fe^III-X bond breaking increases in the series amine < imidazole < thioether, mirroring trends in hardness of these ligands. Experimental rate constants for X dissociation in differently ligated cyt c variants are consistent with this sequence, but the differences between Met and His dissociation rates are attenuated because the former process is limited by the heme crevice opening. Analyses of activation parameters and comparisons to those for the Lys-to-Met ligand switch in the alkaline transition suggest that ligand dissociation is entropically driven in all the variants and accompanied by Lys protonation at neutral pH. The described thiolate redox-linked switches have offered a wealth of new information about interactions of different protein-derived ligands with the heme iron in cyt c model proteins, and we anticipate that the strategy of employing these switches could benefit studies of other redox metalloproteins and model complexes.

Contributions to cytochrome c inner- and outer-sphere reorganization energy

Cytochromes c are small water-soluble proteins that catalyze electron transfer in metabolism and energy conversion processes. Hydrogenobacter thermophilus cytochrome c₅₅₂ presents a curious case in displaying fluxionality of its heme axial methionine ligand; this behavior is altered by single point mutation of the Q64 residue to N64 or V64, which fixes the ligand in a single configuration. The reorganization energy (λ) of these cytochrome c₅₅₂ variants is experimentally determined using a combination of rotating disc electrochemistry, chronoamperometry and cyclic voltammetry. The differences between the λ determined from these complementary techniques helps to deconvolute the contribution of the active site and its immediate environment to the overall λ (λ_Total). The experimentally determined λ values in conjunction with DFT calculations indicate that the differences in λ among the protein variants are mainly due to the differences in contributions from the protein environment and not just inner-sphere λ. DFT calculations indicate that the position of residue 64, responsible for the orientation of the axial methionine, determines the geometric relaxation of the redox active molecular orbital (RAMO). The orientation of the RAMO with respect to the heme is key to determining electron transfer coupling (H_AB) which results in higher ET rates in the wild-type protein relative to the Q64V mutant despite a 150 mV higher λ_Total in the former.

Efficient delocalization of the redox-active molecular orbital (RAMO) in HtWT results in an increase in H_AB value which in turn accelerates the electron transfer (ET) rate in spite of the higher reorganization energy (λ) than the HtQ64V mutant.

A High-Resolution View of the Coordination Environment in a Paramagnetic Metalloprotein from its Magnetic Properties

Metalloproteins constitute a significant fraction of the proteome of all organisms and their characterization is critical for both basic sciences and biomedical applications. A large portion of metalloproteins bind paramagnetic metal ions, and paramagnetic NMR spectroscopy has been widely used in their structural characterization. However, the signals of nuclei in the immediate vicinity of the metal center are often broadened beyond detection. In this work, we show that it is possible to determine the coordination environment of the paramagnetic metal in the protein at a resolution inaccessible to other techniques. Taking the structure of a diamagnetic analogue as a starting point, a geometry optimization is carried out by fitting the pseudocontact shifts obtained from first principles quantum chemical calculations to the experimental ones.

Ultrafast dynamics of heme distortion in the O2-sensor of a thermophilic anaerobe bacterium

Heme-Nitric oxide and Oxygen binding protein domains (H-NOX) are found in signaling pathways of both prokaryotes and eukaryotes and share sequence homology with soluble guanylate cyclase, the mammalian NO receptor. In bacteria, H-NOX is associated with kinase or methyl accepting chemotaxis domains. In the O₂-sensor of the strict anaerobe Caldanaerobacter tengcongensis (Ct H-NOX) the heme appears highly distorted after O₂ binding, but the role of heme distortion in allosteric transitions was not yet evidenced. Here, we measure the dynamics of the heme distortion triggered by the dissociation of diatomics from Ct H-NOX using transient electronic absorption spectroscopy in the picosecond to millisecond time range. We obtained a spectroscopic signature of the heme flattening upon O₂ dissociation. The heme distortion is immediately (<1 ps) released after O₂ dissociation to produce a relaxed state. This heme conformational change occurs with different proportions depending on diatomics as follows: CO < NO < O₂. Our time-resolved data demonstrate that the primary structural event of allostery is the heme distortion in the Ct H-NOX sensor, contrastingly with hemoglobin and the human NO receptor, in which the primary structural events are respectively the motion of the proximal histidine and the rupture of the iron-histidine bond.

Characterization of a nitrite-reducing octaheme hydroxylamine oxidoreductase that lacks the tyrosine cross-link

The hydroxylamine oxidoreductase (HAO) family consists of octaheme proteins that harbor seven bis-His ligated electron-transferring hemes and one 5-coordinate catalytic heme with His axial ligation. Oxidative HAOs have a homotrimeric configuration with the monomers covalently attached to each other via a unique double cross-link between a Tyr residue and the catalytic heme moiety of an adjacent subunit. This cross-linked active site heme, termed the P460 cofactor, has been hypothesized to modulate enzyme reactivity toward oxidative catalysis. Conversely, the absence of this cross-link is predicted to favor reductive catalysis. However, this prediction has not been directly tested. In this study, an HAO homolog that lacks the heme-Tyr cross-link (HAOr) was purified to homogeneity from the nitrite-dependent anaerobic ammonium-oxidizing (anammox) bacterium Kuenenia stuttgartiensis, and its catalytic and spectroscopic properties were assessed. We show that HAOr reduced nitrite to nitric oxide and also reduced nitric oxide and hydroxylamine as nonphysiological substrates. In contrast, HAOr was not able to oxidize hydroxylamine or hydrazine supporting the notion that cross-link-deficient HAO enzymes are reductases. Compared with oxidative HAOs, we found that HAOr harbors an active site heme with a higher (at least 80 mV) midpoint potential and a much lower degree of porphyrin ruffling. Based on the physiology of anammox bacteria and our results, we propose that HAOr reduces nitrite to nitric oxide in vivo, providing anammox bacteria with NO, which they use to activate ammonium in the absence of oxygen.

Ab Initio Evaluation of the Redox Potential of Cytochrome c

Various biochemical activities of metabolism and biosynthesis are fulfilled by redox processes with explicit electron exchange, which furnish redox enzymes with high chemical reactivity. However, theoretical investigation of a redox process, which simultaneously involves a complex electronic change at a redox metal center and conformational reorganization of the surrounding protein environment coupled to the electronic change, requires computationally conflicting approaches, highly accurate quantum chemical calculations, and long-time molecular dynamics (MD) simulations, limiting the physicochemical understanding of biological redox processes. Here, we theoretically examined a redox process of cytochrome c by means of a hybrid molecular simulation technique, which enables one to consistently treat the redox center at the ab initio quantum chemistry level of theory and the protein reorganization with long-time MD simulations on the microsecond timescale. The calculations successfully evaluated a large absolute redox potential, 4.34 eV, with errors of only 0.03 to 0.34 eV to the experimental ones without any problem-specific empirical parameters. Through the long-time MD sampling, large and nonlinear reorganization of the protein environment was unveiled and the molecular determinants for the redox potential were identified. The present ab initio approach significantly expands the applicability of theoretical investigation to biological redox systems with more electronically complicated redox centers such as polynuclear transition metal complexes.

A Quantum Chemistry View on Two Archetypical Paramagnetic Pentacoordinate Nickel(II) Complexes Offers a Fresh Look on Their NMR Spectra

Quantum chemical methods for calculating paramagnetic NMR observables are becoming increasingly accessible and are being included in the inorganic chemistry practice. Here, we test the performance of these methods in the prediction of proton hyperfine shifts of two archetypical high-spin pentacoordinate nickel(II) complexes (NiSAL-MeDPT and NiSAL-HDPT), which, for a variety of reasons, turned out to be perfectly suited to challenge the predictions to the finest level of detail. For NiSAL-MeDPT, new NMR experiments yield an assignment that perfectly matches the calculations. The slightly different hyperfine shifts from the two “halves” of the molecules related by a pseudo-C₂ axis, which are experimentally divided into two well-defined spin systems, are also straightforwardly distinguished by the calculations. In the case of NiSAL-HDPT, for which no X-ray structure is available, the quality of the calculations allowed us to refine its structure using as a starting template the structure of NiSAL-MeDPT.

State-of-the-art quantum chemical methods and paramagnetism-tailored NMR experiments provide a deep insight on the relation between the spectra and the electronic structure for two paramagnetic pentacoordinate nickel(II) complexes.

Histidine versus Cysteine-Bearing Heme-Dependent Halogen Peroxidases: Parallels and Differences for Cl- Oxidation

The homodimeric myeloperoxidase (MPO) features a histidine as a proximal ligand and a sulfonium linkage covalently attaching the heme porphyrin ring to the protein. MPO is able to catalyze Cl^- oxidation with about the same efficiency as chloroperoxidase at pH 7.0. In this study, we seek to explore the parallels and differences between the histidine and cysteine heme-dependent halogen peroxidases. Transition states, reaction barriers, and relevant thermodynamic properties are calculated on protein models. Together with electronic structure calculations, it gives an overview of the reaction mechanisms and of the factors that determine the selectivity between one- and two-electron paths. Conclusions point to the innate oxidizing nature of MPO with the ester and sulfonium linkages hiking up the reactivity to enable chloride oxidation. The installation of a deprotonated imidazolate as a proximal ligand does not shift the equilibrium from one- to two-electron events without influencing the chemistry of the oxidation reaction.

Spin cascade and doming in ferric hemes: Femtosecond X-ray absorption and X-ray emission studies

The structure-function relationship is at the heart of biology, and major protein deformations are correlated to specific functions. For ferrous heme proteins, doming is associated with the respiratory function in hemoglobin and myoglobins. Cytochrome c (Cyt c) has evolved to become an important electron-transfer protein in humans. In its ferrous form, it undergoes ligand release and doming upon photoexcitation, but its ferric form does not release the distal ligand, while the return to the ground state has been attributed to thermal relaxation. Here, by combining femtosecond Fe Kα and Kβ X-ray emission spectroscopy (XES) with Fe K-edge X-ray absorption near-edge structure (XANES), we demonstrate that the photocycle of ferric Cyt c is entirely due to a cascade among excited spin states of the iron ion, causing the ferric heme to undergo doming, which we identify. We also argue that this pattern is common to a wide diversity of ferric heme proteins, raising the question of the biological relevance of doming in such proteins.

Paramagnetic Shifts and Guest Exchange Kinetics in ConFe4-n Metal-Organic Capsules

We investigate the magnetic resonance properties and exchange kinetics of guest molecules in a series of hetero-bimetallic capsules, [Co_nFe_4-nL₆]^4- (n = 1-3), where L^2- = 4,4'-bis[(2-pyridinylmethylene)amino]-[1,1'-biphenyl]-2,2'-disulfonate. H bond networks between capsule sulfonates and guanidinium cations promote the crystallization of [Co_nFe_4-nL₆]^4-. The following four isostructural crystals are reported: two guest-free forms, (C(NH₂)₃)₄[Co_1.8Fe_2.2L₆]·69H₂O (1) and (C(NH₂)₃)₄[Co_2.7Fe_1.3L₆]·73H₂O (2), and two Xe- and CFCl₃-encapsulated forms, (C(NH₂)₃)₄[(Xe)_0.8Co_1.8Fe_2.2L₆]·69H₂O (3) and (C(NH₂)₃)₄[(CFCl₃)Co_2.0Fe_2.0L₆]·73H₂O (4), respectively. Structural analyses reveal that Xe induces negligible structural changes in 3, while the angles between neighboring phenyl groups expand by ca. 3° to accommodate the much larger guest, CFCl₃, in 4. These guest-encapsulated [Co_nFe_4-nL₆]^4- molecules reveal ¹²⁹Xe and ¹⁹F chemical shift changes of ca. -22 and -10 ppm at 298 K, respectively, per substitution of low-spin Fe^II by high-spin Co^II. Likewise, the temperature dependence of the ¹²⁹Xe and ¹⁹F NMR resonances increases by 0.1 and 0.06 ppm/K, respectively, with each additional paramagnetic Co^II center. The optimal temperature for hyperpolarized (hp) ¹²⁹Xe chemical exchange saturation transfer (hyper-CEST) with [Co_nFe_4-nL₆]^4- capsules was found to be inversely proportional to the number of Co^II centers, n, which is consistent with the Xe chemical exchange accelerating as the portals expand. The systematic study was facilitated by the tunability of the [M₄L₆]^4- capsules, further highlighting these metal-organic systems for developing responsive sensors with highly shifted ¹²⁹Xe resonances.

Biological and Bioinspired Inorganic N-N Bond-Forming Reactions

The metallobiochemistry underlying the formation of the inorganic N-N-bond-containing molecules nitrous oxide (N2O), dinitrogen (N2), and hydrazine (N2H4) is essential to the lifestyles of diverse organisms. Similar reactions hold promise as means to use N-based fuels as alternative carbon-free energy sources. This review discusses research efforts to understand the mechanisms underlying biological N-N bond formation in primary metabolism and how the associated reactions are tied to energy transduction and organismal survival. These efforts comprise studies of both natural and engineered metalloenzymes as well as synthetic model complexes.

The Heme-Lys Cross-Link in Cytochrome P460 Promotes Catalysis by Enforcing Secondary Coordination Sphere Architecture

Cytochrome (cyt) P460 is a c-type monoheme enzyme found in ammonia-oxidizing bacteria (AOB) and methanotrophs; additionally, genes encoding it have been found in some pathogenic bacteria. Cyt P460 is defined by a unique post-translational modification to the heme macrocycle, where a lysine (Lys) residue covalently attaches to the 13' meso carbon of the porphyrin, modifying this heme macrocycle into the enzyme's eponymous P460 cofactor, similar to the cofactor found in the enzyme hydroxylamine oxidoreductase. This cross-link imbues the protein with unique spectroscopic properties, the most obvious of which is the enzyme's green color in solution. Cyt P460 from the AOB Nitrosomonas europaea is a homodimeric redox enzyme that produces nitrous oxide (N₂O) from 2 equiv of hydroxylamine. Mutation of the Lys cross-link results in spectroscopic features that are more similar to those of standard cyt c' proteins and renders the enzyme catalytically incompetent for NH₂OH oxidation. Recently, the necessity of a second-sphere glutamate (Glu) residue for redox catalysis was established; it plausibly serves as proton relay during the first oxidative half of the catalytic cycle. Herein, we report the first crystal structure of a cross-link deficient cyt P460. This structure shows that the positioning of the catalytically essential Glu changes by approximately 0.8 Å when compared to a cross-linked, catalytically competent cyt P460. It appears that the heme-Lys cross-link affects the relative position of the P460 cofactor with respect to the second-sphere Glu residue, therefore dictating the catalytic competency of the enzyme.

Iron transitions during activation of allosteric heme proteins in cell signaling

Allosteric heme proteins can fulfill a very large number of different functions thanks to the remarkable chemical versatility of heme through the entire living kingdom. Their efficacy resides in the ability of heme to transmit both iron coordination changes and iron redox state changes to the protein structure. Besides the properties of iron, proteins may impose a particular heme geometry leading to distortion, which allows selection or modulation of the electronic properties of heme. This review focusses on the mechanisms of allosteric protein activation triggered by heme coordination changes following diatomic binding to proteins as diverse as the human NO-receptor, cytochromes, NO-transporters and sensors, and a heme-activated potassium channel. It describes at the molecular level the chemical capabilities of heme to achieve very different tasks and emphasizes how the properties of heme are determined by the protein structure. Particularly, this reviews aims at giving an overview of the exquisite adaptability of heme, from bacteria to mammals.

ONIOM investigations of the heme degradation mechanism by MhuD: the critical function of heme ruffling

A unique ruffling conformation of hydroxyheme in MhuD inhibits its “on-site” monooxygenation but induces “remote-site” dioxygenation.

Ruffling drives coproheme decarboxylation by facilitating PCET: a theoretical investigation of ChdC

Coproheme decarboxylase (ChdC) is an essential enzyme in the coproporphyrin-dependent heme synthesis pathway, which catalyzes oxidative decarboxylation of coproheme at the positions p2 and p4 to generate heme b under the action of hydrogen peroxide.

Spectroscopic Evidence for Electronic Control of Heme Hydroxylation by IsdG

Staphylococcus aureus IsdG catalyzes a unique trioxygenation of heme to staphylobilin, and the data presented in this article elucidate the mechanism of the novel chemical transformation. More specifically, the roles of the second-sphere Asn and Trp residues in the monooxygenation of ferric-peroxoheme have been clarified via spectroscopic characterization of the ferric-azidoheme analogue. Analysis of UV/vis absorption data quantified the strength of the hydrogen bond that exists between the Asn7 side chain and the azide moiety of ferric-azidoheme. X-band electron paramagnetic resonance data were acquired and analyzed, which revealed that this hydrogen bond weakens the π-donor strength of the azide, resulting in perturbations of the Fe 3d based orbitals. Finally, nuclear magnetic resonance characterization of ¹³C-enriched samples demonstrated that the Asn7···N₃ hydrogen bond triggers partial porphyrin to iron electron transfer, resulting in spin density delocalization onto the heme meso carbons. These spectroscopic experiments were complemented by combined quantum mechanics/molecular mechanics computational modeling, which strongly suggested that the electronic structure changes observed for the N7A variant arose from loss of the Asn7···N₃ hydrogen bond as opposed to a decrease in porphyrin ruffling. From these data a fascinating picture emerges where an Asn7···N₃ hydrogen bond is communicated through four bonds, resulting in meso carbons with partial cationic radical character that are poised for hydroxylation. This chemistry is not observed in other heme proteins because Asn7 and Trp67 must work in concert to trigger the requisite electronic structure change.

A 192-heme electron transfer network in the hydrazine dehydrogenase complex

Anaerobic ammonium oxidation (anammox) is a major process in the biogeochemical nitrogen cycle in which nitrite and ammonium are converted to dinitrogen gas and water through the highly reactive intermediate hydrazine. So far, it is unknown how anammox organisms convert the toxic hydrazine into nitrogen and harvest the extremely low potential electrons (-750 mV) released in this process. We report the crystal structure and cryo electron microscopy structures of the responsible enzyme, hydrazine dehydrogenase, which is a 1.7 MDa multiprotein complex containing an extended electron transfer network of 192 heme groups spanning the entire complex. This unique molecular arrangement suggests a way in which the protein stores and releases the electrons obtained from hydrazine conversion, the final step in the globally important anammox process.

A 60-heme reductase complex from an anammox bacterium shows an extended electron transfer pathway

The hydroxylamine oxidoreductase/hydrazine dehydrogenase (HAO/HDH) protein family constitutes an important group of octaheme cytochromes c (OCCs). The majority of these proteins form homotrimers, with their subunits being covalently attached to each other via a rare cross-link between the catalytic heme moiety and a conserved tyrosine residue in an adjacent subunit. This covalent cross-link has been proposed to modulate the active-site heme towards oxidative catalysis by distorting the heme plane. In this study, the crystal structure of a stable complex of an HAO homologue (KsHAOr) with its diheme cytochrome c redox partner (KsDH) from the anammox bacterium Kuenenia stuttgartiensis was determined. KsHAOr lacks the tyrosine cross-link and is therefore tuned to reductive catalysis. The molecular model of the KsHAOr–KsDH complex at 2.6 Å resolution shows a heterododecameric (α6β6) assembly, which was also shown to be the oligomeric state in solution by analytical ultracentrifugation and multi-angle static light scattering. The 60-heme-containing protein complex reveals a unique extended electron transfer pathway and provides deeper insights into catalysis and electron transfer in reductive OCCs.

One fold, two functions: cytochrome P460 and cytochrome c'-β from the methanotroph Methylococcus capsulatus (Bath)

Nature is adept at utilising highly similar protein folds to carry out very different functions, yet the mechanisms by which this functional divergence occurs remain poorly characterised. In certain methanotrophic bacteria, two homologous pentacoordinate c-type heme proteins have been identified: a cytochrome P460 (cyt P460) and a cytochrome c'-β (cyt cp-β). Cytochromes P460 are able to convert hydroxylamine to nitrous oxide (N₂O), a potent greenhouse gas. This reactivity is similar to that of hydroxylamine oxidoreductase (HAO), which is a key enzyme in nitrifying and methanotrophic bacteria. Cyt P460 and HAO both have unusual protein-heme cross-links, formed by a Tyr residue in HAO and a Lys in cyt P460. In contrast, cyts cp-β (the only known cytochromes c' with a β-sheet fold) lack this crosslink and appears to be optimized for binding non-polar molecules (including NO and CO) without enzymatic conversion. Our bioinformatics analysis supports the proposal that cyt cp-β may have evolved from cyt P460 via a gene duplication event. Using high-resolution X-ray crystallography, UV-visible absorption, electron paramagnetic resonance (EPR) and resonance Raman spectroscopy, we have characterized the overall protein folding and active site structures of cyt cp-β and cyt P460 from the obligate methanotroph, Methylococcus capsulatus (Bath). These proteins display a similar β-sheet protein fold, together with a pattern of changes to the heme pocket regions and localised tertiary structure that have converted a hydroxylamine oxidizing enzyme into a gas-binding protein. Structural comparisons provide insights relevant to enzyme redesign for synthetic enzymology and engineering of gas sensor proteins. We also show the widespread occurrence of cyts cp-β and characterise their phylogeny.

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.

Catalytic iron-carbene intermediate revealed in a cytochrome c carbene transferase

Recently, heme proteins have been discovered and engineered by directed evolution to catalyze chemical transformations that are biochemically unprecedented. Many of these nonnatural enzyme-catalyzed reactions are assumed to proceed through a catalytic iron porphyrin carbene (IPC) intermediate, although this intermediate has never been observed in a protein. Using crystallographic, spectroscopic, and computational methods, we have captured and studied a catalytic IPC intermediate in the active site of an enzyme derived from thermostable Rhodothermus marinus (Rma) cytochrome c High-resolution crystal structures and computational methods reveal how directed evolution created an active site for carbene transfer in an electron transfer protein and how the laboratory-evolved enzyme achieves perfect carbene transfer stereoselectivity by holding the catalytic IPC in a single orientation. We also discovered that the IPC in Rma cytochrome c has a singlet ground electronic state and that the protein environment uses geometrical constraints and noncovalent interactions to influence different IPC electronic states. This information helps us to understand the impressive reactivity and selectivity of carbene transfer enzymes and offers insights that will guide and inspire future engineering efforts.

The Eponymous Cofactors in Cytochrome P460s from Ammonia-Oxidizing Bacteria Are Iron Porphyrinoids Whose Macrocycles Are Dibasic

The enzymes hydroxylamine oxidoreductase and cytochrome (cyt) P460 contain related unconventional "heme P460" cofactors. These cofactors are unusual in their inclusion of nonstandard cross-links between amino acid side chains and the heme macrocycle. Mutagenesis studies performed on the Nitrosomonas europaea cyt P460 that remove its lysine-heme cross-link show that the cross-link is key to defining the spectroscopic properties and kinetic competence of the enzyme. However, exactly how this cross-link confers these features remains unclear. Here we report the 1.45 Å crystal structure of cyt P460 from Nitrosomonas sp. AL212 and conclude that the cross-link does not lead to a change in hybridization of the heme carbon participating in the cross-link but rather enforces structural distortions to the macrocycle away from planarity. Time-dependent density functional theory coupled to experimental structural and spectroscopic analysis suggest that this geometric distortion is sufficient to define the spectroscopic properties of the heme P460 cofactor and provide clues toward establishing a relationship between heme P460 electronic structure and function.

Influences of the heme-lysine crosslink in cytochrome P460 over redox catalysis and nitric oxide sensitivity

Ammonia (NH₃)-oxidizing bacteria (AOB) derive total energy for life from the multi-electron oxidation of NH₃ to nitrite (NO₂^-). One obligate intermediate of this metabolism is hydroxylamine (NH₂OH), which can be oxidized to the potent greenhouse agent nitrous oxide (N₂O) by the AOB enzyme cytochrome (cyt) P460. We have now spectroscopically characterized a 6-coordinate (6c) {FeNO}⁷ intermediate on the NH₂OH oxidation pathway of cyt P460. This species has two fates: it can either be oxidized to the {FeNO}⁶ that then undergoes attack by NH₂OH to ultimately generate N₂O, or it can lose its axial His ligand, thus generating a stable, off-pathway 5-coordinate (5c) {FeNO}⁷ species. We show that the wild type (WT) cyt P460 exhibits a slow nitric oxide (NO)-independent conversion (k_His-off = 2.90 × 10^-3 s^-1), whereas a cross-link-deficient Lys70Tyr cyt P460 mutant protein underwent His dissociation via both a NO-independent (k_His-off = 3.8 × 10^-4 s^-1) and a NO-dependent pathway [k_His-off(NO) = 790 M^-1 s^-1]. Eyring analyses of the NO-independent pathways for these two proteins revealed a significantly larger (ca. 27 cal mol^-1 K^-1) activation entropy (ΔS^‡) in the cross-link-deficient mutant. Our results suggest that the Lys-heme cross-link confers rigidity to the positioning of the heme P460 cofactor to avoid the fast NO-dependent His dissociation pathway and subsequent formation of the off-pathway 5c {FeNO}⁷ species. The relevance of these findings to NO signaling proteins such as heme-nitric oxide/oxygen binding (H-NOX) is also 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.

NADH reduction of nitroaromatics as a probe for residual ferric form high-spin in a cytochrome P450

The existence of a substrate-sensitive equilibrium between high spin (S=5/2) and low spin (S=1/2) ferric iron is a well-established phenomenon in the cytochrome P450 (CYP) superfamily, although its origins are still a subject of discussion. A series of mutations that strongly perturb the spin state equilibrium in the camphor hydroxylase CYP101A1 were recently described (Colthart et al., Sci. Rep. 6, 22035 (2016)). Wild type CYP101A1 as well as some CYP101A1 mutants are herein shown to be capable of catalyzing the reduction of nitroacetophenones by NADH to the corresponding anilino compounds (nitroreductase or NRase activity). The distinguishing characteristic between those mutants that catalyze the reduction and those that cannot appears to be the extent to which residual high spin form exists in the absence of the native substrate d-camphor, with those showing the largest spin state shifts upon camphor binding also exhibiting NRase activity. Optical and EPR spectroscopy was used to further examine these phenomena. These results suggest that reduction of nitroaromatics may provide a useful probe of residual high spin states in the CYP superfamily. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Elucidation of the heme active site electronic structure affecting the unprecedented nitrite dismutase activity of the ferriheme b proteins, the nitrophorins

Nitrophorins (NPs) catalyze the nitrite dismutation reaction that is unprecedented in ferriheme proteins. Despite progress in studying the reaction mechanism, fundamental issues regarding the correlation of the structural features with the nitrite dismutase activity of NPs remain elusive. On the other hand, it has been shown that the nitrite complexes of NPs are unique among those of the ferriheme proteins since some of their electron paramagnetic resonance (EPR) spectra show significant highly anisotropic low spin (HALS) signals with large gmax values over 3.2. The origin of HALS signals in ferriheme proteins or models is not well understood, especially in cases where axial ligands other than histidine are present. In this study several mutations were introduced in NP4. The related nitrite coordination and dismutation reaction were investigated. As a result, the EPR spectra of the NP-nitrite complexes were found to be tightly correlated with the extent of heme ruffling and protonation state of the proximal His ligand-dictated by an extended H-bonding network at the heme active site. Furthermore, it is established that the two factors are essential in determining the nitrite dismutase activity of NPs. These results may provide a valuable guide for identifying or designing novel heme proteins with similar activity.

Conformational control of cofactors in nature - the influence of protein-induced macrocycle distortion on the biological function of tetrapyrroles

Tetrapyrrole-containing proteins are one of the most fundamental classes of enzymes in nature and it remains an open question to give a chemical rationale for the multitude of biological reactions that can be catalyzed by these pigment-protein complexes. There are many fundamental processes where the same (i.e., chemically identical) porphyrin cofactor is involved in chemically quite distinct reactions. For example, heme is the active cofactor for oxygen transport and storage (hemoglobin, myoglobin) and for the incorporation of molecular oxygen in organic substrates (cytochrome P450). It is involved in the terminal oxidation (cytochrome c oxidase) and the metabolism of H2O2 (catalases and peroxidases) and catalyzes various electron transfer reactions in cytochromes. Likewise, in photosynthesis the same chlorophyll cofactor may function as a reaction center pigment (charge separation) or as an accessory pigment (exciton transfer) in light harvesting complexes (e.g., chlorophyll a). Whilst differences in the apoprotein sequences alone cannot explain the often drastic differences in physicochemical properties encountered for the same cofactor in diverse protein complexes, a critical factor for all biological functions must be the close structural interplay between bound cofactors and the respective apoprotein in addition to factors such as hydrogen bonding or electronic effects. Here, we explore how nature can use the same chemical molecule as a cofactor for chemically distinct reactions using the concept of conformational flexibility of tetrapyrroles. The multifaceted roles of tetrapyrroles are discussed in the context of the current knowledge on distorted porphyrins. Contemporary analytical methods now allow a more quantitative look at cofactors in protein complexes and the development of the field is illustrated by case studies on hemeproteins and photosynthetic complexes. Specific tetrapyrrole conformations are now used to prepare bioengineered designer proteins with specific catalytic or photochemical properties.

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.

Dynamic ruffling distortion of the heme substrate in non-canonical heme oxygenase enzymes

A double-well exists along the ruffling coordinate of cyanide-inhibited ferric heme, which explains the observation of “nested” VTVH MCD saturation magnetization curves.

Biological significance and applications of heme c proteins and peptides

Hemes are ubiquitous in biology and carry out a wide range of functions. The heme group is largely invariant across proteins with different functions, although there are a few variations seen in nature. The most common variant is heme c, which is formed by a post-translational modification in which heme is covalently linked to two Cys residues on the polypeptide via thioether bonds. In this Account, the influence of this covalent attachment on heme c properties and function is discussed, and examples of how covalent attachment has been used in selected applications are presented. Proteins that bind heme c are among the most well-characterized proteins in biochemistry. Most of these proteins are cytochromes c (cyts c) that serve as electron carriers in photosynthesis and respiration. Despite the intense study of cyts c, the functional significance of heme covalent attachment has remained elusive. One observation is that heme c reaches a lower reduction potential in nature than its noncovalently linked counterpart, heme b, when comparing proteins with the same axial ligands. Furthermore, covalent attachment is known to enhance protein stability and allow the heme to be relatively solvent exposed. However, an inorganic chemistry perspective on the effects of covalent attachment has been lacking. Spectroscopic measurements and computations on cyts c and model systems reveal a number of effects of covalent attachment on heme electronic structure and reactivity. One is that the predominant nonplanar ruffling distortion seen in heme c lowers heme reduction potential. Another is that covalent attachment influences the interaction of the heme iron with the proximal His ligand. Heme ruffling also has been shown to influence electronic coupling to redox partners and, therefore, electron transfer rates by altering the distribution of the orbital hole on the porphyrin in oxidized cyt c. Another consequence of heme covalent attachment is the strong vibrational coupling seen between the iron and the protein surface as revealed by nuclear resonance vibrational spectroscopy studies. Finally, heme covalent attachment is proposed to be an important feature supporting multiple roles of cyt c in programmed cell death (apoptosis). Heme covalent attachment is not only vital for the biological functions of cyt c but also provides a useful handle in a number of applications. For one, the engineering of heme c onto an exposed portion of a protein of interest has been shown to provide a visible affinity purification tag. In addition, peptides with covalently attached heme, known as microperoxidases, have been studied as model compounds and oxidation catalysts and, more recently, in applications for energy conversion and storage. The wealth of insight gained about heme c through fundamental studies of cyts c forms a basis for future efforts toward engineering natural and artificial cytochromes for a variety of applications.

The low-affinity complex of cytochrome c and its peroxidase

The complex of yeast cytochrome c peroxidase and cytochrome c is a paradigm of the biological electron transfer (ET). Building on seven decades of research, two different models have been proposed to explain its functional redox activity. One postulates that the intermolecular ET occurs only in the dominant, high-affinity protein–protein orientation, while the other posits formation of an additional, low-affinity complex, which is much more active than the dominant one. Unlike the high-affinity interaction—extensively studied by X-ray crystallography and NMR spectroscopy—until now the binding of cytochrome c to the low-affinity site has not been observed directly, but inferred mainly from kinetics experiments. Here we report the structure of this elusive, weak protein complex and show that it consists of a dominant, inactive bound species and an ensemble of minor, ET-competent protein–protein orientations, which summarily account for the experimentally determined value of the ET rate constant.

The redox activity of cytochrome c in complex with its peroxidase has been rationalized by two possible models; a readily observable high-affinity complex and a more elusive but potentially more active low-affinity complex. Here, the authors provide an NMR-based structural mapping of this low-affinity complex.