Heme Electron Transfer in Peroxidases: The Propionate e-Pathway

Tyrosyl Radical-Mediated Sequential Oxidative Decarboxylation of Coproporphyrinogen III through PCET: Theoretical Insights into the Mechanism of Coproheme Decarboxylase ChdC

The peroxide-dependent coproheme decarboxylase ChdC from Geobacillus stearothermophilus catalyzes two key steps in the synthesis of heme b, i.e., two sequential oxidative decarboxylations of coproporphyrinogen III (coproheme III) at propionate groups P2 and P4. In the binding site of coproheme III, P2 and P4 are anchored by different residues (Tyr144, Arg217, and Ser222 for P2 and Tyr113, Lys148, and Trp156 for P4); however, strong experimental evidence supports that the generated Tyr144 radical acts as an unique intermediary for hydrogen atom transfer (HAT) from both reactive propionates. So far, the reaction details are still unclear. Herein, we carried out quantum mechanics/molecular mechanics calculations to explore the decarboxylation mechanism of coproheme III. In our calculations, the coproheme Cpd I, Fe(IV) = O coupled to a porphyrin radical cation (por^•+) with four propionate groups, was used as a reactant model. Our calculations reveal that Tyr144 is directly involved in the decarboxylation of propionate group P2. First, the proton-coupled electron transfer (PCET) occurs from Tyr144 to P2, generating a Tyr144 radical, which then abstracts a hydrogen atom from the Cβ of P2. The β-H extraction was calculated to be the rate-limiting step of decarboxylation. It is the porphyrin radical cation (por^•+) that makes the PCET from Tyr144 to P2 to be quite easy to initiate the decarboxylation. Finally, the electron transfers from the Cβ^• through the porphyrin to the iron center, leading to the decarboxylation of P2. Importantly, the decarboxylation of P4 mediated by Lys148 was calculated to be very difficult, which suggests that after the P2 decarboxylation, the generated harderoheme III intermediate should rebind or rotate in the active site so that the propionate P4 occupies the binding site of P2, and Tyr144 again mediates the decarboxylation of P4. Thus, our calculations support the fact that Tyr144 is responsible for the decarboxylation of both P2 and P4.

Visualizing the protons in a metalloenzyme electron proton transfer pathway

In redox metalloenzymes, the process of electron transfer often involves the concerted movement of a proton. These processes are referred to as proton-coupled electron transfer, and they underpin a wide variety of biological processes, including respiration, energy conversion, photosynthesis, and metalloenzyme catalysis. The mechanisms of proton delivery are incompletely understood, in part due to an absence of information on exact proton locations and hydrogen bonding structures in a bona fide metalloenzyme proton pathway. Here, we present a 2.1-Å neutron crystal structure of the complex formed between a redox metalloenzyme (ascorbate peroxidase) and its reducing substrate (ascorbate). In the neutron structure of the complex, the protonation states of the electron/proton donor (ascorbate) and all of the residues involved in the electron/proton transfer pathway are directly observed. This information sheds light on possible proton movements during heme-catalyzed oxygen activation, as well as on ascorbate oxidation.

A Heme Propionate Staples the Structure of Cytochrome c for Methionine Ligation to the Heme Iron

Ligand-switch reactions at the heme iron are common in biological systems, but their mechanisms and the features of the polypeptide fold that support dual ligation are not well understood. In cytochrome c (cyt c), two low-stability loops (Ω-loop C and Ω-loop D) are connected by the heme propionate HP6. At alkaline pH, the native Met80 ligand from Ω-loop D switches to a Lys residue from the same loop. Deprotonation of an as yet unknown group triggers the alkaline transition. We have created the two cyt c variants T49V/K79G and T78V/K79G with altered connections of these two loops to HP6. Electronic absorption, NMR, and EPR studies demonstrate that at pH 7.4 ferric forms of these variants are Lys-ligated, whereas ferrous forms maintain the native Met80 ligation. Measurements of protein stability, cyclic voltammetry, pH-jump and gated electron-transfer kinetics have revealed that these Thr to Val substitutions greatly affect the alkaline transition in both ferric and ferrous proteins. The substitutions modify the stability of the Met-ligated species and reduction potentials of the heme iron. The kinetics of ligand-switch processes are also altered, and analyses of these effects implicate redox-dependent differences in metal-ligand interactions and the role of the protein dynamics, including cross-talk between the two Ω-loops. With the two destabilized variants, it is possible to map energy levels for the Met- and Lys-ligated species in both ferric and ferrous proteins and assess the role of the protein scaffold in redox-dependent preferences for these two ligands. The estimated shift in the heme iron reduction potential upon deprotonation of the "trigger" group is consistent with those associated with deprotonation of an HP, suggesting that HP6, on its own or as a part of a hydrogen-bonded cluster, is a likely "trigger" for the Met to Lys ligand switch.

Oxygen Activation and Radical Transformations in Heme Proteins and Metalloporphyrins

As a result of the adaptation of life to an aerobic environment, nature has evolved a panoply of metalloproteins for oxidative metabolism and protection against reactive oxygen species. Despite the diverse structures and functions of these proteins, they share common mechanistic grounds. An open-shell transition metal like iron or copper is employed to interact with O2 and its derived intermediates such as hydrogen peroxide to afford a variety of metal-oxygen intermediates. These reactive intermediates, including metal-superoxo, -(hydro)peroxo, and high-valent metal-oxo species, are the basis for the various biological functions of O2-utilizing metalloproteins. Collectively, these processes are called oxygen activation. Much of our understanding of the reactivity of these reactive intermediates has come from the study of heme-containing proteins and related metalloporphyrin compounds. These studies not only have deepened our understanding of various functions of heme proteins, such as O2 storage and transport, degradation of reactive oxygen species, redox signaling, and biological oxygenation, etc., but also have driven the development of bioinorganic chemistry and biomimetic catalysis. In this review, we survey the range of O2 activation processes mediated by heme proteins and model compounds with a focus on recent progress in the characterization and reactivity of important iron-oxygen intermediates. Representative reactions initiated by these reactive intermediates as well as some context from prior decades will also be presented. We will discuss the fundamental mechanistic features of these transformations and delineate the underlying structural and electronic factors that contribute to the spectrum of reactivities that has been observed in nature as well as those that have been invented using these paradigms. Given the recent developments in biocatalysis for non-natural chemistries and the renaissance of radical chemistry in organic synthesis, we envision that new enzymatic and synthetic transformations will emerge based on the radical processes mediated by metalloproteins and their synthetic analogs.

Mapping the Long-Range Electron Transfer Route in Ligninolytic Peroxidases

Combining a computational analysis with site-directed mutagenesis, we have studied the long-range electron transfer pathway in versatile and lignin peroxidases, two enzymes of biotechnological interest that play a key role for fungal degradation of the bulky lignin molecule in plant biomass. The in silico study established two possible electron transfer routes starting at the surface tryptophan residue previously identified as responsible for oxidation of the bulky lignin polymer. Moreover, in both enzymes, a second buried tryptophan residue appears as a top electron transfer carrier, indicating the prevalence of one pathway. Site-directed mutagenesis of versatile peroxidase (from Pleurotus eryngii) allowed us to corroborate the computational analysis and the role played by the buried tryptophan (Trp244) and a neighbor phenylalanine residue (Phe198), together with the surface tryptophan, in the electron transfer. These three aromatic residues are highly conserved in all the sequences analyzed (up to a total of 169). The importance of the surface (Trp171) and buried (Trp251) tryptophan residues in lignin peroxidase has been also confirmed by directed mutagenesis of the Phanerochaete chrysosporium enzyme. Overall, the combined procedure identifies analogous electron transfer pathways in the long-range oxidation mechanism for both ligninolytic peroxidases, constituting a good example of how computational analysis avoids making extensive trial-error mutagenic experiments.

Structure-Based Mechanism for Oxidative Decarboxylation Reactions Mediated by Amino Acids and Heme Propionates in Coproheme Decarboxylase (HemQ)

Coproheme decarboxylase catalyzes two sequential oxidative decarboxylations with H2O2 as the oxidant, coproheme III as substrate and cofactor, and heme b as the product. Each reaction breaks a C-C bond and results in net loss of hydride, via steps that are not clear. Solution and solid-state structural characterization of the protein in complex with a substrate analog revealed a highly unconventional H2O2-activating distal environment with the reactive propionic acids (2 and 4) on the opposite side of the porphyrin plane. This suggested that, in contrast to direct C-H bond cleavage catalyzed by a high-valent iron intermediate, the coproheme oxidations must occur through mediating amino acid residues. A tyrosine that hydrogen bonds to propionate 2 in a position analogous to the substrate in ascorbate peroxidase is essential for both decarboxylations, while a lysine that salt bridges to propionate 4 is required solely for the second. A mechanism is proposed in which propionate 2 relays an oxidizing equivalent from a coproheme compound I intermediate to the reactive deprotonated tyrosine, forming Tyr•. This residue then abstracts a net hydrogen atom (H•) from propionate 2, followed by migration of the unpaired propionyl electron to the coproheme iron to yield the ferric harderoheme and CO2 products. A similar pathway is proposed for decarboxylation of propionate 4, but with a lysine residue as an essential proton shuttle. The proposed reaction suggests an extended relay of heme-mediated e-/H+ transfers and a novel route for the conversion of carboxylic acids to alkenes.

Structural basis for cellobiose dehydrogenase action during oxidative cellulose degradation

A new paradigm for cellulose depolymerization by fungi focuses on an oxidative mechanism involving cellobiose dehydrogenases (CDH) and copper-dependent lytic polysaccharide monooxygenases (LPMO); however, mechanistic studies have been hampered by the lack of structural information regarding CDH. CDH contains a haem-binding cytochrome (CYT) connected via a flexible linker to a flavin-dependent dehydrogenase (DH). Electrons are generated from cellobiose oxidation catalysed by DH and shuttled via CYT to LPMO. Here we present structural analyses that provide a comprehensive picture of CDH conformers, which govern the electron transfer between redox centres. Using structure-based site-directed mutagenesis, rapid kinetics analysis and molecular docking, we demonstrate that flavin-to-haem interdomain electron transfer (IET) is enabled by a haem propionate group and that rapid IET requires a closed CDH state in which the propionate is tightly enfolded by DH. Following haem reduction, CYT reduces LPMO to initiate oxygen activation at the copper centre and subsequent cellulose depolymerization.

Insights into Laccase Engineering from Molecular Simulations: Toward a Binding-Focused Strategy

Understanding the molecular determinants of enzyme performance is of primary importance for the rational design of ad hoc mutants. A novel approach, which combines efficient conformational sampling and quick reactivity scoring, is used here to shed light on how substrate oxidation was improved during the directed evolution experiment of a fungal laccase (from Pycnoporus cinnabarinus), an industrially relevant class of oxidoreductases. It is found that the enhanced activity of the evolved enzyme is mainly the result of substrate arrangement in the active site, with no important change in the redox potential of the T1 copper. Mutations at the active site shift the binding mode into a more buried substrate position and provide a more favorable electrostatic environment for substrate oxidation. As a consequence, engineering the binding event seems to be a viable way to in silico evolution of oxidoreductases.

Formation of a tyrosine adduct involved in lignin degradation by Trametopsis cervina lignin peroxidase: a novel peroxidase activation mechanism

LiP (lignin peroxidase) from Trametopsis cervina has an exposed catalytic tyrosine residue (Tyr181) instead of the tryptophan conserved in other lignin-degrading peroxidases. Pristine LiP showed a lag period in VA (veratryl alcohol) oxidation. However, VA-LiP (LiP after treatment with H2O2 and VA) lacked this lag, and H2O2-LiP (H2O2-treated LiP) was inactive. MS analyses revealed that VA-LiP includes one VA molecule covalently bound to the side chain of Tyr181, whereas H2O2-LiP contains a hydroxylated Tyr181. No adduct is formed in the Y171N variant. Molecular docking showed that VA binding is favoured by sandwich π stacking with Tyr181 and Phe89. EPR spectroscopy after peroxide activation of the pre-treated LiPs showed protein radicals other than the tyrosine radical found in pristine LiP, which were assigned to a tyrosine–VA adduct radical in VA-LiP and a dihydroxyphenyalanine radical in H2O2-LiP. Both radicals are able to oxidize large low-redox-potential substrates, but H2O2-LiP is unable to oxidize high-redox-potential substrates. Transient-state kinetics showed that the tyrosine–VA adduct strongly promotes (>100-fold) substrate oxidation by compound II, the rate-limiting step in catalysis. The novel activation mechanism is involved in ligninolysis, as demonstrated using lignin model substrates. The present paper is the first report on autocatalytic modification, resulting in functional alteration, among class II peroxidases.

Two tyrosyl radicals stabilize high oxidation states in cytochrome C oxidase for efficient energy conservation and proton translocation

The reaction of oxidized bovine cytochrome c oxidase (bCcO) with hydrogen peroxide (H(2)O(2)) was studied by electron paramagnetic resonance (EPR) to determine the properties of radical intermediates. Two distinct radicals with widths of 12 and 46 G are directly observed by X-band EPR in the reaction of bCcO with H(2)O(2) at pH 6 and pH 8. High-frequency EPR (D-band) provides assignments to tyrosine for both radicals based on well-resolved g-tensors. The wide radical (46 G) exhibits g-values similar to a radical generated on L-Tyr by UV-irradiation and to tyrosyl radicals identified in many other enzyme systems. In contrast, the g-values of the narrow radical (12 G) deviate from L-Tyr in a trend akin to the radicals on tyrosines with substitutions at the ortho position. X-band EPR demonstrates that the two tyrosyl radicals differ in the orientation of their β-methylene protons. The 12 G wide radical has minimal hyperfine structure and can be fit using parameters unique to the post-translationally modified Y244 in bCcO. The 46 G wide radical has extensive hyperfine structure and can be fit with parameters consistent with Y129. The results are supported by mixed quantum mechanics and molecular mechanics calculations. In addition to providing spectroscopic evidence of a radical formed on the post-translationally modified tyrosine in CcO, this study resolves the much debated controversy of whether the wide radical seen at low pH in the bovine enzyme is a tyrosine or tryptophan. The possible role of radical formation and migration in proton translocation is discussed.

Multireference and multiconfiguration ab initio methods in heme-related systems: what have we learned so far?

This work reviews the recent applications of ab initio multireference/multiconfiguration (MR/MC) electronic structure methods to heme-related systems, involving tetra-, penta-, and hexa-coordinate species, as well as the high-valent iron-oxo species. The current accuracy of these methods in the various systems is discussed, with special attention to potential sources of systematic errors. Thus, the review summarizes and tries to rationalize the key elements of MR/MC calculations, namely, the choice of the employed active space, especially the so-called double-shell effect that has already been recognized to be important in transition-metal-containing systems, and the impact of these elements on the spin-state energetics of heme species, as well as on the bonding mechanism of small molecules to the heme. It is shown that expansion of the MC wave function into one based on localized orbitals provides a compact and insightful view on some otherwise complex electronic structures. The effects of protein environment on the MR/MC results are summarized for the few available quantum mechanical/molecular mechanical (QM/MM) studies. Comparisons with corresponding DFT results are also made wherever available. Potential future directions are proposed.

Proton-coupled electron transfer reactions at a heme-propionate in an iron-protoporphyrin-IX model compound

A heme model system has been developed in which the heme-propionate is the only proton donating/accepting site, using protoporphyrin IX-monomethyl esters (PPIX(MME)) and N-methylimidazole (MeIm). Proton-coupled electron transfer (PCET) reactions of these model compounds have been examined in acetonitrile solvent. (PPIX(MME))Fe(III)(MeIm)(2)-propionate (Fe(III)~CO(2)) is readily reduced by the ascorbate derivative 5,6-isopropylidine ascorbate to give (PPIX(MME))Fe(II)(MeIm)(2)-propionic acid (Fe(II)~CO(2)H). An excess of the hydroxylamine TEMPOH or of hydroquinone similarly reduces Fe(III)~CO(2), and TEMPO and benzoquinone oxidize Fe(II)~CO(2)H to return to Fe(III)~CO(2). The measured equilibrium constants, and the determined pK(a) and E(1/2) values, indicate that Fe(II)~CO(2)H has an effective bond dissociation free energy (BDFE) of 67.8 ± 0.6 kcal mol(-1). In these PPIX models, electron transfer occurs at the iron center and proton transfer occurs at the remote heme propionate. According to thermochemical and other arguments, the TEMPOH reaction occurs by concerted proton-electron transfer (CPET), and a similar pathway is indicated for the ascorbate derivative. Based on these results, heme propionates should be considered as potential key components of PCET/CPET active sites in heme proteins.

Crystallographic, kinetic, and spectroscopic study of the first ligninolytic peroxidase presenting a catalytic tyrosine

Trametes cervina lignin peroxidase (LiP) is a unique enzyme lacking the catalytic tryptophan strictly conserved in all other LiPs and versatile peroxidases (more than 30 sequences available). Recombinant T. cervina LiP and site-directed variants were investigated by crystallographic, kinetic, and spectroscopic techniques. The crystal structure shows three substrate oxidation site candidates involving His-170, Asp-146, and Tyr-181. Steady-state kinetics for oxidation of veratryl alcohol (the typical LiP substrate) by variants at the above three residues reveals a crucial role of Tyr-181 in LiP activity. Moreover, assays with ferrocytochrome c show that its ability to oxidize large molecules (a requisite property for oxidation of the lignin polymer) originates in Tyr-181. This residue is also involved in the oxidation of 1,4-dimethoxybenzene, a reaction initiated by the one-electron abstraction with formation of substrate cation radical, as described for the well known Phanerochaete chrysosporium LiP. Detailed spectroscopic and kinetic investigations, including low temperature EPR, show that the porphyrin radical in the two-electron activated T. cervina LiP is unstable and rapidly receives one electron from Tyr-181, forming a catalytic protein radical, which is identified as an H-bonded neutral tyrosyl radical. The crystal structure reveals a partially exposed location of Tyr-181, compatible with its catalytic role, and several neighbor residues probably contributing to catalysis: (i) by enabling substrate recognition by aromatic interactions; (ii) by acting as proton acceptor/donor from Tyr-181 or H-bonding the radical form; and (iii) by providing the acidic environment that would facilitate oxidation. This is the first structure-function study of the only ligninolytic peroxidase described to date that has a catalytic tyrosine.

Nitric oxide synthase stabilizes the tetrahydrobiopterin cofactor radical by controlling its protonation state

Nitric oxide synthase (NOS), a homodimeric enzyme with a flavin reductase domain and a P450-type heme-containing oxygenase domain, catalyzes the formation of NO from L-arginine, NADPH, and O(2) in a two-step reaction sequence. In the first step, a tetrahydrobiopterin (H(4)B) cofactor bound near one of the heme propionate groups acts as an electron donor to the P450-type heme active site, yielding a one-electron oxidized radical that is subsequently re-reduced. In solution, H(4)B undergoes two-electron oxidation, showing that the enzyme significantly alters the proton- and electron-transfer properties of the cofactor. Multifrequency EPR and ENDOR spectroscopy were used to determine magnetic parameters, and from them the (de)protonation state of the H(4)B radical in the oxygenase domain dimer of inducible NO synthase that was trapped by rapid freeze quench. From 9.5 and 330-416 GHz EPR and from 34 GHz (1)H ENDOR spectroscopy, the g tensor of the radical and the hyperfine tensors of several N and H nuclei in the radical were obtained. Density functional theory calculations at the PBE0/EPR-II level for H(4)B radical models predict different spin density distributions and g and hyperfine tensors for different protonation states. Comparison of the predicted and experimental values leads to the conclusion that the radical is cationic H(4)B(*+), suggesting that NOS stabilizes this protonated form to utilize the cofactor in a unique dual one-electron redox role, where it can deliver an electron to the active site for reductive oxygen activation and also remove an electron from the active site to generate NO and not NO(-). The protein environment also prevents further oxidation and subsequent loss of function of the cofactor, thus enabling the enzyme to perform the unusual catalytic one-electron chemistry.

Electronic structures of heme a of cytochrome c oxidase in the redox states--charge density migration to the propionate groups of heme a

The electronic structures of heme a of cytochrome c oxidase in the redox states were studied, using hybrid density functional theory with a polarizable continuum model and a point charge model. We found that the most stable electronic configurations of the d electrons of the Fe ion are determined by the orbital interactions of the d orbitals of the Fe ion with the pi orbitals of the porphyrin ring and the His residues. The redox reaction of the Fe ion influences the charge density on the formyl group through the pi conjugation of the porphyrin ring. In addition, we found the charge transfer from the Fe ion to the propionate group of heme a in the redox change despite the lack of the pi-conjugation. We elucidated that the charge propagation originates from the heme a structure itself and that the origin of the charge delocalization to the heme propionate is the orbital interactions between the d orbital of the Fe ion and the p orbitals of the carboxylate part of the heme propionate via the pi conjugation of the porphyrin ring and the sigma* orbital of the C-C bond of the propionate group. The electrostatic effect by surrounding proteins enhances the charge transfer from the Fe ion to the propionate group. These results indicate that heme propionate groups serve electron mediators in electron transfer as well as electrostatic anchors, and that proteins surrounding the active site reinforce the congenital abilities of the cofactors.

QM/MM methods: looking inside heme proteins biochemistry

Mixed quantum mechanics/molecular mechanics methods offer a valuable computational tool for understanding biochemical events. When combined with conformational sampling techniques, they allow for an exhaustive exploration of the enzymatic mechanism. Heme proteins are ubiquitous and essential for every organism. In this review we summarize our efforts towards the understanding of heme biochemistry. We present: 1) results on ligand migration on globins coupled to the ligand binding event, 2) results on the localization of the spin density in compound I of cytochromes and peroxidases, 3) novel methodologies for mapping the electron transfer pathways and 4) novel data on Tryptophan 2,3-dioxygenase. For this enzyme our results strongly indicate that the distal oxygen will end up on the C3 indole carbon, whereas the proximal oxygen will end up in the C2 position. Interestingly, the process involves the formation of an epoxide and a heme ferryl intermediate. The overall energy profile indicates an energy barrier of approximately 18 kcal/mol and an exothermic driving force of almost 80 kcal/mol.