Theoretical studies of enzyme mechanisms involving high-valent iron intermediates

The Apparently Unreactive Substrate Facilitates the Electron Transfer for Dioxygen Activation in Rieske Dioxygenases

Rieske dioxygenases belong to the non‐heme iron family of oxygenases and catalyze important cis‐dihydroxylation as well as O‐/N‐dealkylation and oxidative cyclization reactions for a wide range of substrates. The lack of substrate coordination at the non‐heme ferrous iron center, however, makes it particularly challenging to delineate the role of the substrate for productive O2 activation. Here, we studied the role of the substrate in the key elementary reaction leading to O2 activation from a theoretical perspective by systematically considering (i) the 6‐coordinate to 5‐coordinate conversion of the non‐heme Fe^II upon abstraction of a water ligand, (ii) binding of O2 , and (iii) transfer of an electron from the Rieske cluster. We systematically evaluated the spin‐state‐dependent reaction energies and structural effects at the active site for all combinations of the three elementary processes in the presence and absence of substrate using naphthalene dioxygenase as a prototypical Rieske dioxygenase. We find that reaction energies for the generation of a coordination vacancy at the non‐heme FeII center through thermoneutral H₂O reorientation and exothermic O2 binding prior to Rieske cluster oxidation are largely insensitive to the presence of naphthalene and do not lead to formation of any of the known reactive Fe‐oxygen species. By contrast, the role of the substrate becomes evident after Rieske cluster oxidation and exclusively for the 6‐coordinate non‐heme FeII sites in that the additional electron is found at the substrate instead of at the iron and oxygen atoms. Our results imply an allosteric control of the substrate on Rieske dioxygenase reactivity to happen prior to changes at the non‐heme FeII in agreement with a strategy that avoids unproductive O2 activation.

Which is the real oxidant in the competitive ligand self-hydroxylation and substrate oxidation, a biomimetic iron(II)-hydroperoxo species or an oxo-iron(IV)-hydroxy one?

Nonheme iron(II)-hydroperoxo species (FeII-(η2-OOH)) 1 and the concomitant oxo-iron(IV)-hydroxyl one 2 are proposed as the key intermediates of a large class of 2-oxoglutarate dependent dioxygenases (e.g., isopenicillin N synthase). Extensive...

Epigenetic Metalloenzymes

Epigenetics controls the expression of genes and is responsible for cellular phenotypes. The fundamental basis of these mechanisms involves in part the post-translational modifications (PTMs) of DNA and proteins, in particular, the nuclear histones. DNA can be methylated or demethylated on cytosine. Histones are marked by several modifications including acetylation and/or methylation, and of particular importance are the covalent modifications of lysine. There exists a balance between addition and removal of these PTMs, leading to three groups of enzymes involved in these processes: the writers adding marks, the erasers removing them, and the readers able to detect these marks and participating in the recruitment of transcription factors. The stimulation or the repression in the expression of genes is thus the result of a subtle equilibrium between all the possibilities coming from the combinations of these PTMs. Indeed, these mechanisms can be deregulated and then participate in the appearance, development and maintenance of various human diseases, including cancers, neurological and metabolic disorders. Some of the key players in epigenetics are metalloenzymes, belonging mostly to the group of erasers: the zinc-dependent histone deacetylases (HDACs), the iron-dependent lysine demethylases of the Jumonji family (JMJ or KDM) and for DNA the iron-dependent ten-eleven-translocation enzymes (TET) responsible for the oxidation of methylcytosine prior to the demethylation of DNA. This review presents these metalloenzymes, their importance in human disease and their inhibitors.

Fe2 O+ Cation Mediated Propane Oxidation by Dioxygen in the Gas Phase

The mass-selected Fe₂ O⁺ cation mediated propane oxidation by O₂ was investigated by mass spectrometry and density functional theory calculations. In the reaction of Fe₂ O⁺ with C₃ H₈ , H₂ was liberated by C-H bond activation to give Fe₂ OC₃ H₆⁺ . Interestingly, when a mixture of C₃ H₈ /O₂ was introduced into the reactor, an intense signal that corresponded to the Fe₂ O₂⁺ cation was present; the experiments indicated that O₂ was activated in its reaction with Fe₂ O(C₃ H₆ )⁺ to give Fe₂ O₂⁺ and C₃ H₆ O (acetone or propanal). A Langmuir-Hinshelwood-like mechanism was adopted in the propane oxidation reaction by O₂ on gas-phase Fe₂ O⁺ cations. In comparison with the absence of Fe₂ O₂⁺ in the reaction of Fe₂ O⁺ with O₂ , the ligand effect of C₃ H₆ on Fe₂ OC₃ H₆⁺ is important in the oxygen activation reaction. The theoretical results are consistent with the experimental observations. The propane oxidation by O₂ in the presence of Fe₂ O⁺ might be applied as a model for alkane and O₂ activations over iron oxide catalysts, and the mechanisms and kinetic data are useful for understanding corresponding heterogeneous reactions.

Mechanistic insights into intramolecular ortho-amination/hydroxylation by nonheme FeIV[double bond, length as m-dash]NTs/FeIV[double bond, length as m-dash]O species: the σ vs. the π channels

Comparative oxidative abilities of nonheme Fe^IV[double bond, length as m-dash]NTs and Fe^IV[double bond, length as m-dash]O species using DFT has been explored. Our calculations reveal that the Fe^IV[double bond, length as m-dash]NTs is found to be a stronger oxidant in two electron transfer reactions and react exclusively via π channels while the Fe^IV[double bond, length as m-dash]O species is found to be a stronger oxidant when the σ-pathway is activated such as in HAT reactions.

Asymmetric abstraction of two chemically-equivalent methylene hydrogens: significant enantioselectivity of endoperoxide presented by fumitremorgin B endoperoxidase

The unique enantioselectivity for an R-chiral endoperoxy ring in verruculogen biosynthesis originates from asymmetric abstraction of two chemically-equivalent methylene hydrogens.

Interaction with the Redox Cofactor MYW and Functional Role of a Mobile Arginine in Eukaryotic Catalase-Peroxidase

Catalase-peroxidases (KatGs) are unique bifunctional heme peroxidases with an additional posttranslationally formed redox-active Met-Tyr-Trp cofactor that is essential for catalase activity. On the basis of studies of bacterial KatGs, controversial mechanisms of hydrogen peroxide oxidation were proposed. The recent discovery of eukaryotic KatGs with differing pH optima of catalase activity now allows us to scrutinize those postulated reaction mechanisms. In our study, secreted KatG from the fungus Magnaporthe grisea (MagKatG2) was used to analyze the role of a remote KatG-typical mobile arginine that was shown to interact with the Met-Tyr-Trp adduct in a pH-dependent manner in bacterial KatGs. Here we present crystal structures of MagKatG2 at pH 3.0, 5.5, and 7.0 and investigate the mobility of Arg461 by molecular dynamics simulation. Data suggest that at pH ≥4.5 Arg461 mostly interacts with the deprotonated adduct Tyr. Elimination of Arg461 by mutation to Ala slightly increases the thermal stability but does not alter the active site architecture or the kinetics of cyanide binding. However, the variant Arg461Ala lost the wild-type-typical optimum of catalase activity at pH 5.25 (k_cat = 6450 s^–1) but exhibits a broad plateau between pH 4.5 and 7.5 (k_cat = 270 s^–1 at pH 5.5). Moreover, significant differences in the kinetics of interconversion of redox intermediates of wild-type and mutant protein mixed with either peroxyacetic acid or hydrogen peroxide are observed. These findings together with published data from bacterial KatGs allow us to propose a role of Arg461 in the H₂O₂ oxidation reaction of KatG.

Benchmark study on methanol C-H and O-H bond activation by bare [Fe(IV)O](2+)

We present a high-level computational study on methanol C-H and O-H bond cleavages by bare [Fe(IV)O](2+), as well as benchmarks of various density functional theory (DFT) methods. We considered direct and concerted hydrogen transfer (DHT and CHT) pathways, respectively. The potential energy surfaces were constructed at the CCSD(T)/def2-TZVPP//B3LYP/def2-TZVP level of theory. Mechanistically, (1) the C-H bond cleavage is dominant and the O-H activation only plays minor role on the PESs; (2) the DHT from methyl should be the most practical channel; and (3) electronic structure analysis demonstrates the proton and electron transfer coupling behavior along the reaction coordinates. The solvent effect is evident and plays distinct roles in regulating the two bond activations in different mechanisms during the catalysis. The effect of optimizing the geometries using different density functionals was also studied, showing that it is not meaningful to discuss which DFT method could give the accurate prediction of the geometries, especially for transition structures. Furthermore, the gold-standard CCSD(T) method was used to benchmark 19 different density functionals with different Hartree-Fock exchange fractions. The results revealed that (i) the structural factor plays a minor role in the single point energy (SPE) calculations; (ii) reaction energy prediction is quite challenging for DFT methods; (iii) the mean absolute deviations (MADs) reflect the problematic description of the DFs when dealing with metal oxidation state change, giving a strong correlation on the HF exchange in the DFs. Knowledge from this study should be of great value for computational chemistry, especially for the de novo design of transition metal catalysts.

Mechanism of benzylic hydroxylation by 4-hydroxymandelate synthase. A computational study

Hydroxymandelate synthase (HMS) and 4-hydroxyphenylpyruvate dioxygenase (HPPD) are highly related enzymes using the same substrates but catalyzing hydroxylation reactions yielding different products. The first steps of the HMS and HPPD catalytic reactions are believed to proceed in the same way and lead to an Fe(IV)═O-hydroxyphenylacetate (HPA) intermediate. Further down the catalytic cycles, HMS uses Fe(IV)═O to perform hydroxylation of the benzylic carbon, whereas in HPPD, the reactive oxoferryl intermediate attacks the aromatic ring of HPA. This study focuses on this part of the HMS catalytic cycle that starts from the oxoferryl intermediate and aims to identify interactions within the active site that are responsible for enzyme specificity. To this end, a HMS-Fe(IV)═O-HPA complex was modeled with molecular dynamics simulations. On the basis of the molecular dynamics-equilibrated structure, an active site model suitable for quantum chemical investigations was constructed and used for density functional theory (B3LYP) calculations of the mechanism of the native reaction of HMS, i.e., benzylic hydroxylation, and the alternative electrophilic attack on the ring, which is a step of the HPPD catalytic cycle. The most important result of this study is the finding that the conformation of the Ser201 side chain in the second coordination shell has a key role in directing the reaction of Fe(IV)═O into either the HMS or the HPPD channel.

Hydroxylation catalysis by mononuclear and dinuclear iron oxo catalysts: a methane monooxygenase model system versus the Fenton reagent Fe(IV)O(H2O)5(2+)

Hydroxylation of aliphatic C-H bonds is a chemically and biologically important reaction, which is catalyzed by the oxidoiron group FeO(2+) in both mononuclear (heme and nonheme) and dinuclear complexes. We investigate the similarities and dissimilarities of the action of the FeO(2+) group in these two configurations, using the Fenton-type reagent [FeO(2+) in a water solution, FeO(H(2)O)(5)(2+)] and a model system for the methane monooxygenase (MMO) enzyme as representatives. The high-valent iron oxo intermediate MMOH(Q) (compound Q) is regarded as the active species in methane oxidation. We show that the electronic structure of compound Q can be understood as a dimer of two Fe(IV)O(2+) units. This implies that the insights from the past years in the oxidative action of this ubiquitous moiety in oxidation catalysis can be applied immediately to MMOH(Q). Electronically the dinuclear system is not fundamentally different from the mononuclear system. However, there is an important difference of MMOH(Q) from FeO(H(2)O)(5)(2+): the largest contribution to the transition state (TS) barrier in the case of MMOH(Q) is not the activation strain (which is in this case the energy for the C-H bond lengthening to the TS value), but it is the steric hindrance of the incoming CH(4) with the ligands representing glutamate residues. The importance of the steric factor in the dinuclear system suggests that it may be exploited, through variation in the ligand framework, to build a synthetic oxidation catalyst with the desired selectivity for the methane substrate.

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.

Reaction mechanism of apocarotenoid oxygenase (ACO): a DFT study

The mechanism of the oxidative cleavage catalyzed by apocarotenoid oxygenase (ACO) was studied by using a quantum chemical (DFT: B3 LYP) method. Based on the available crystal structure, relatively large models of the unusual active-site region, in which a ferrous ion is coordinated by four histidines and no negatively charged ligand, were selected and used in the computational investigation of the reaction mechanism. The results suggest that binding of dioxygen to the ferrous ion in the active site promotes one-electron oxidation of carotenoid leading to a substrate radical cation and a Fe-bound superoxide radical. Recombination of the two radicals, which can be realized in at least two different ways, yields a reactive peroxo species that subsequently evolves into either a dioxetane or an epoxide intermediate. The former easily decays into the final aldehyde products, whereas the oxidation of the epoxide to the proper products of the reaction requires involvement of a water molecule. The calculated activation barriers favor the dioxetane mechanism, yet the mechanism involving the epoxide intermediate cannot be ruled out.

Protein Environment Facilitates O2 Binding in Non-Heme Iron Enzyme. An Insight from ONIOM Calculations on Isopenicillin N Synthase (IPNS)

Binding of dioxygen to a non-heme enzyme has been modeled using the ONIOM combined quantum mechanical/molecular mechanical (QM/MM) method. For the present system, isopenicillin N synthase (IPNS), binding of dioxygen is stabilized by 8-10 kcal/mol for a QM:MM (B3LYP:Amber) protein model compared to a quantum mechanical model of the active site only. In the protein system, the free energy change of O2 binding is close to zero. Two major factors consistently stabilize O2 binding. The first effect, evaluated at the QM level, originates from a change in coordination geometry of the iron center. The active-site model artificially favors the deoxy state (O2 not bound) because it allows too-large rearrangements of the five-coordinate iron site. This error is corrected when the protein is included. The corresponding effect on binding energies is 3-6 kcal/mol, depending on the coordination mode of O2 (side-on or end-on). The second major factor that stabilizes O2 binding is van der Waals interactions between dioxygen and the surrounding enzyme. These interactions, 3-4 kcal/mol at the MM level, are neglected in models that include only the active site. Polarization of the active site by surrounding amino acids does not have a significant effect on the binding energy in the present system.

The structures and electronic configuration of compound I intermediates of Helicobacter pylori and Penicillium vitale catalases determined by X-ray crystallography and QM/MM density functional theory calculations

The structures of Helicobacter pylori (HPC) and Penicillium vitale (PVC) catalases, each with two subunits in the crystal asymmetric unit, oxidized with peroxoacetic acid are reported at 1.8 and 1.7 A resolution, respectively. Despite the similar oxidation conditions employed, the iron-oxygen coordination length is 1.72 A for PVC, close to what is expected for a Fe=O double bond, and 1.80 and 1.85 A for HPC, suggestive of a Fe-O single bond. The structure and electronic configuration of the oxoferryl heme and immediate protein environment is investigated further by QM/MM density functional theory calculations. Four different active site electronic configurations are considered, Por*+-FeIV=O, Por*+-FeIV=O...HisH+, Por*+-FeIV-OH+ and Por-FeIV-OH (a protein radical is assumed in the latter configuration). The electronic structure of the primary oxidized species, Por*+-FeIV=O, differs qualitatively between HPC and PVC with an A2u-like porphyrin radical delocalized on the porphyrin in HPC and a mixed A1u-like "fluctuating" radical partially delocalized over the essential distal histidine, the porphyrin, and, to a lesser extent, the proximal tyrosine residue. This difference is rationalized in terms of HPC containing heme b and PVC containing heme d. It is concluded that compound I of PVC contains an oxoferryl Por*+-FeIV=O species with partial protonation of the distal histidine and compound I of HPC contains a hydroxoferryl Por-FeIV-OH with the second oxidation equivalent delocalized as a protein radical. The findings support the idea that there is a relation between radical migration to the protein and protonation of the oxoferryl bond in catalase.

Oxidative properties of FeO2+: electronic structure and solvation effects

An electronic structure analysis is provided of the action of solvated FeO(2+), [FeO(H(2)O)(5)](2+), as a hydroxylation catalyst. It is emphasized that the oxo end of FeO(2+) does not form hydrogen bonds (as electron donor and H-bond acceptor) with H-bond donors nor with aliphatic C-H bonds, but it activates C-H bonds as an electron acceptor. It is extremely electrophilic, to the extent that it can activate even such poor electron donors as aliphatic C-H bonds, the C-H bond orbital acting as electron donor in a charge transfer type of interaction. Lower lying O-H bonding orbitals are less easily activated. The primary electron accepting orbital in a water environment is the 3sigma*alpha orbital, an antibonding combination of Fe-3d(z(2)) and O-2p(z), which is very low-lying relative to the pi*alpha compared with, for example, the sigma* orbital in O(2) relative to its pi*. This is ascribed to relatively small Fe-3d(z(2)) with O-2p(z) overlap, due to the nodal structure of the 3d(z(2)).The H-abstraction barrier is very low in the gas phase, but it is considerably enhanced in water solvent. This is shown to be due to strong screening effects of the dielectric medium, leading to relative destabilization of the levels of the charged [FeO(H(2)O)(5)](2+) species compared to those of the neutral substrate molecules, making it a less effective electron acceptor. The solvent directly affects the orbital interactions responsible for the catalytic reaction.

Transition metal spin state energetics and noninnocent systems: challenges for DFT in the bioinorganic arena

Although density functional theory (DFT) provides a generally good description of transition metal systems, we have identified several cases, involving Fe(III) porphyrins and related systems, where common functionals fail to correctly describe the energetics of the different low-lying spin states. The question of metal- versus ligand-centered oxidation in high-valent transition metal complexes is also a challenging one for DFT calculations, as I have tried to illustrate with examples from among porphyrin, corrole, biliverdine, and NO complexes. In a number of cases, I have compared results obtained with different exchange-correlation functionals; in addition, I have added a discussion on the relative performance of pure versus hybrid functionals. Finally, I have offered some thoughts on the role that traditional wavefunction-based ab initio methods, now essentially absent from the bioinorganic arena, might play in the future.