A Density Functional Study of Possible Intermediates of the Reaction of Dioxygen Molecule with Non-Heme Iron Complexes. 1. N-Side versus O-Side Mechanism with Water-Free Model

Structure-Spectroscopy Correlations for Intermediate Q of Soluble Methane Monooxygenase: Insights from QM/MM Calculations

The determination of the diiron core intermediate structures involved in the catalytic cycle of soluble methane monooxygenase (sMMO), the enzyme that selectively catalyzes the conversion of methane to methanol, has been a subject of intense interest within the bioinorganic scientific community. Particularly, the specific geometry and electronic structure of the intermediate that precedes methane binding, known as intermediate Q (or MMOHQ), has been debated for over 30 years. Some reported studies support a bis-μ-oxo-bridged Fe(IV)2O2 closed-core conformation Fe(IV)2O2 core, whereas others favor an open-core geometry, with a longer Fe-Fe distance. The lack of consensus calls for a thorough re-examination and reinterpretation of the spectroscopic data available on the MMOHQ intermediate. Herein, we report extensive simulations based on a hybrid quantum mechanics/molecular mechanics approach (QM/MM) approach that takes into account the complete enzyme to explore possible conformations for intermediates MMOHox and MMOHQ of the sMMOH catalytic cycle. High-level quantum chemical approaches are used to correlate specific structural motifs with geometric parameters for comparison with crystallographic and EXAFS data, as well as with spectroscopic data from Mössbauer spectroscopy, Fe K-edge high-energy resolution X-ray absorption spectroscopy (HERFD XAS), and resonance Raman 16O-18O difference spectroscopy. The results provide strong support for an open-core-type configuration in MMOHQ, with the most likely topology involving mono-oxo-bridged Fe ions and alternate terminal Fe-oxo and Fe-hydroxo groups that interact via intramolecular hydrogen bonding. The implications of an open-core intermediate Q on the reaction mechanism of sMMO are discussed.

Intermediates in dioxygen activation by methane monooxygenase: a QM/MM study

Protein effects in the activation of dioxygen by methane monooxygenase (MMO) were investigated by using combined QM/MM and broken-symmetry Density Functional Theory (DFT) methods. The effects of a novel empirical scheme recently developed by our group on the relative DFT energies of the various intermediates in the catalytic cycle are investigated. Inclusion of the protein leads to much better agreement between the experimental and computed geometric structures for the reduced form (MMOH(red)). Analysis of the electronic structure of MMOH(red) reveals that the two iron atoms have distinct environments. Different coordination geometries tested for the MMOH(peroxo) intermediate reveal that, in the protein environment, the mu-eta2,eta2 structure is more stable than the others. Our analysis also shows that the protein helps to drive reactants toward products along the reaction path. Furthermore, these results demonstrate the importance of including the protein environment in our models and the usefulness of the QM/MM approach for accurate modeling of enzymatic reactions. A discrepancy remains in our calculation of the Fe-Fe distance in our model of HQ as compared to EXAFS data obtained several years ago, for which we currently do not have an explanation.

Hydroxylation of methane by non-heme diiron enzymes: molecular orbital analysis of C-H bond activation by reactive intermediate Q

The electronic structures of key species involved in methane hydroxylation performed by the hydroxylase component of soluble methane monooxygenase (sMMO), as proposed previously on the basis of high-level density functional theory, were investigated. The reaction starts with initial approach of methane at one of the bridging oxo atoms in intermediate Q, a di(mu-oxo)diiron(IV) unit. This step is accompanied by a proton-coupled outer-sphere transfer of the first electron from a C-H sigma-bond in methane to one of the metal centers. The second electron transfer, also an outer-sphere electron transfer process, occurs along a two-component reaction pathway. Both redox reactions are strongly coupled to structural distortions of the diiron core. The electronic consequence and driving force of these distortions are intuitively explained by using the computed Kohn-Sham orbitals in the broken-symmetry framework to incorporate the experimentally observed antiferromagnetic coupling of the unpaired electrons at the metal centers. The broken-symmetry orbital scheme is essential for describing the C-H bond activation process in a consistent and complete manner, enabling derivation of both an intuitive and quantitative understanding of the most salient electronic features that govern the details of the hydroxylation.

Theoretical study of the mechanism of alkane hydroxylation and ethylene epoxidation reactions catalyzed by diiron bis-oxo complexes. The effect of substrate molecules

The hybrid density functional method B3LYP was used to study the mechanism of the hydrocarbon (methane, ethane, methyl fluoride, and ethylene) oxidation reaction catalyzed by the complexes cis-(H(2)O)(NH(2))Fe(mu-O)(2)(eta(2)-HCOO)(2)Fe(NH(2))(H(2)O), I, and cis-(HCOO)(Imd)Fe(mu-O)(2)(eta(2)-HCOO)(2)Fe(Imd)(HCOO) (Imd = Imidazole), I_m, the "small" and "medium" model of compound Q of the methane monooxygenase (MMO). The improvement of the model from "small" to "medium" did not change the qualitative conclusions but significantly changed the calculated energetics. As in the case of methane oxidation reported by the authors previously, the reaction of all the substrates studied here is shown to start by coordination of the substrate molecule to the bridging oxygen atom, O(1) of I, an Fe(IV)-Fe(IV) complex, followed by the H-atom abstraction at the transition state III leading to the bound hydroxy alkyl intermediate IV of Fe(III)-Fe(IV) core. IV undergoes a very exothermic coupling of alkyl and hydroxy groups to give the alcohol complex VI of Fe(III)-Fe(III) core, from which alcohol dissociates. The H(b)-atom abstraction (or C-H bond activation) barrier at transition state III is found to be a few kcal/mol lower for C(2)H(6) and CH(3)F than for CH(4). The calculated trend in the H(b)-abstraction barrier, CH(4) (21.8 kcal/mol) > CH(3)F (18.8 kcal/mol) > or = C(2)H(6) (18.5 kcal/mol), is consistent with the C-H(b) bond strength in these substrates. Thus, the weaker the C-H(b) bond, the lower is the H(b)-abstraction barrier. It was shown that the replacement of a H-atom in a methane molecule with a more electronegative group tends to make the H(b)-abstraction transition state less "reactant-like". In contrast, the replacement of the H-atom in CH(4) with a less electronegative group makes the H(b)-abstraction transition state more "reactant-like". The epoxidation of ethylene by complex I is found to proceed without barrier and is a highly exothermic process. Thus, in the reaction of ethylene with complex I the only product is expected to be ethylene oxide, which is consistent with the experiment.

Computational studies of reaction mechanisms of methane monooxygenase and ribonucleotide reductase

An overview of the computational efforts made by our group during the last few years in the field of nonheme diiron proteins is presented. Through application of ab initio methodology to a reasonable set of molecular models, significant progress is made in understanding how the soluble Methane Monooxygenase system achieves the hydroxylation of methane and how the catalytic cycle of Ribonucleotide Reductase is initiated. In particular, the current studies reveal in more detail (1) the nature of key intermediates in the reaction cycles of these two metalloenzymes, (2) details of how the iron centers regulate the systems, and (3) important aspects of how the carboxylate ligands in the active sites may tailor the enzymatic needs of the metalloprotein. This knowledge also leads to novel connections between the two enzymes. The coordinative unsaturation and carboxylate shifts investigated herein are two properties that are likely to be of more general impact in nonheme proteins. The control of the redox chemistry of the enzyme by the binuclear metal center, also analyzed here, should find common ground among other bimetallic systems as well.