Structural Model Studies for the Peroxo Intermediate P and the Reaction Pathway from P → Q of Methane Monooxygenase Using Broken-Symmetry Density Functional Calculations

Mechanistic Insights into the N-Hydroxylations Catalyzed by the Binuclear Iron Domain of SznF Enzyme: Key Piece in the Synthesis of Streptozotocin

SznF, a member of the emerging family of heme-oxygenase-like (HO-like) di-iron oxidases and oxygenases, employs two distinct domains to catalyze the conversion of Nω-methyl-L-arginine (L-NMA) into N-nitroso-containing product, which can subsequently be transformed into streptozotocin. Using unrestricted density functional theory (UDFT) with the hybrid functional B3LYP, we have mechanistically investigated the two sequential hydroxylations of L-NMA catalyzed by SznF's binuclear iron central domain. Mechanism B primarily involves the O-O bond dissociation, forming Fe(IV)=O, induced by the H+/e- introduction to the FeA side of μ-1,2-peroxo-Fe2(III/III), the substrate hydrogen abstraction by Fe(IV)=O, and the hydroxyl rebound to the substrate N radical. The stochastic addition of H+/e- to the FeB side (mechanism C) can transition to mechanism B, thereby preventing enzyme deactivation. Two other competing mechanisms, involving the direct O-O bond dissociation (mechanism A) and the addition of H2O as a co-substrate (mechanism D), have been ruled out.

Second coordination sphere effects on the mechanistic pathways for dioxygen activation by a ferritin: involvement of a Tyr radical and the identification of a cation binding site

Ferritins are ubiquitous diiron enzymes involved in iron(II) detoxification and oxidative stress responses and can act as metabolic iron stores. The overall reaction mechanisms of ferritin enzymes are still unclear, particularly concerning the role of the conserved, near catalytic center Tyr residue. Thus, we carried out a computational study of a ferritin using a large cluster model of well over 300 atoms including its first- and second-coordination sphere. The calculations reveal important insight into the structure and reactivity of ferritins. Specifically, the active site Tyr residue delivers a proton and electron in the catalytic cycle prior to iron(II) oxidation. In addition, the calculations highlight a likely cation binding site at Asp65, which through long-range electrostatic interactions, influences the electronic configuration and charge distributions of the metal center. The results are consistent with experimental observations but reveal novel detail of early mechanistic steps that lead to an unusual mixed-valent iron(III)-iron(II) center.

Mössbauer Property Calculations on Fea33+∙∙∙H2O∙∙∙CuB2+ Dinuclear Center Models of the Resting Oxidized as-Isolated Cytochrome c Oxidase

Mössbauer isomer shift and quadrupole splitting properties have been calculated using the OLYP‐D3(BJ) density functional method on previously obtained (W.‐G. Han Du, et al., Inorg Chem. 2020, 59, 8906–8915) geometry optimized Fe_a3³⁺−H₂O−Cu_B²⁺ dinuclear center (DNC) clusters of the resting oxidized (O state) “as‐isolated” cytochrome c oxidase (CcO). The calculated results are highly consistent with the available experimental observations. The calculations have also shown that the structural heterogeneities of the O state DNCs implicated by the Mössbauer experiments are likely consequences of various factors, particularly the variable positions of the central H₂O molecule between the Fe_a3³⁺ and Cu_B²⁺ sites in different DNCs, whether or not this central H₂O molecule has H‐bonding interaction with another H₂O molecule, the different spin states having similar energies for the Fe_a3³⁺ sites, and whether the Fe_a3³⁺ and Cu_B²⁺ sites are ferromagnetically or antiferromagnetically spin‐coupled.

Understanding desaturation/hydroxylation activity of castor stearoyl Δ9-Desaturase through rational mutagenesis

•

Rationally designed mutations in the Δ⁹ desaturase promoted hydroxylation activity.

•

Proton and electron transfer to the active site is crucial for the Δ⁹D to desaturate

•

Detailed analysis of all enzymatic products of the Δ⁹D was carried out

•

Insight into the chemo-, and stereoselectivity of non-heme diiron enzymes was obtained

A recently proposed reaction mechanism of soluble Δ⁹ desaturase (Δ⁹D) allowed us to identify auxiliary residues His203, Asp101, Thr206 and Cys222 localized near the di-iron active site that are supposedly involved in the proton transfer (PT) to and from the active site. The PT, along with the electron transfer (ET), seems to be crucial for efficient desaturation. Thus, perturbing the major PT chains is expected to impair the native reaction and (potentially) amplify minor reaction channels, such as the substrate hydroxylation. To verify this hypothesis, we mutated the four residues mentioned above into their counterparts present in a soluble methane monooxygenase (sMMO), and determined the reaction products of mutants. We found that the mutations significantly promote residual monohydroxylation activities on stearoyl-CoA, often at the expense of native desaturation activity. The favored hydroxylation positions are C₉, followed by C₁₀ and C₁₁. Reactions with unsaturated substrate, oleoyl-CoA, yield erythro-9,10-diol, cis-9,10-epoxide and a mixture of allylic alcohols. Additionally, using 9- and 11-hydroxystearoyl-CoA, we showed that the desaturation reaction can proceed only with the hydroxyl group at position C₁₁, whereas the hydroxylation reaction is possible in both cases, i.e. with hydroxyl at position C₉ or C₁₁. Despite the fact that the overall outcome of hydroxylation is rather modest and that it is mostly the desaturation/hydroxylation ratio that is affected, our results broaden understanding of the origin of chemo- and stereoselectivity of the Δ⁹D and provide further insight into the catalytic action of the NHFe₂ enzymes.

Proton-Electron Transfer to the Active Site Is Essential for the Reaction Mechanism of Soluble Δ9-Desaturase

A full understanding of the catalytic action of non-heme iron (NHFe) and non-heme diiron (NHFe₂) enzymes is still beyond the grasp of contemporary computational and experimental techniques. Many of these enzymes exhibit fascinating chemo-, regio-, and stereoselectivity, in spite of employing highly reactive intermediates which are necessary for activations of most stable chemical bonds. Herein, we study in detail one intriguing representative of the NHFe₂ family of enzymes: soluble Δ⁹ desaturase (Δ⁹D), which desaturates rather than performing the thermodynamically favorable hydroxylation of substrate. Its catalytic mechanism has been explored in great detail by using QM(DFT)/MM and multireference wave function methods. Starting from the spectroscopically observed 1,2-μ-peroxo diferric P intermediate, the proton-electron uptake by P is the favored mechanism for catalytic activation, since it allows a significant reduction of the barrier of the initial (and rate-determining) H-atom abstraction from the stearoyl substrate as compared to the "proton-only activated" pathway. Also, we ruled out that a Q-like intermediate (high-valent diamond-core bis-μ-oxo-[Fe^IV]₂ unit) is involved in the reaction mechanism. Our mechanistic picture is consistent with the experimental data available for Δ⁹D and satisfies fairly stringent conditions required by Nature: the chemo-, stereo-, and regioselectivity of the desaturation of stearic acid. Finally, the mechanisms evaluated are placed into a broader context of NHFe₂ chemistry, provided by an amino acid sequence analysis through the families of the NHFe₂ enzymes. Our study thus represents an important contribution toward understanding the catalytic action of the NHFe₂ enzymes and may inspire further work in NHFe₍₂₎ biomimetic chemistry.

Approaches for Selective Oxidation of Methane to Methanol

Methane activation chemistry, despite being widely reported in literature, remains to date a subject of debate. The challenges in this reaction are not limited to methane activation but extend to stabilization of the intermediate species. The low C-H dissociation energy of intermediates vs. reactants leads to CO2 formation. For selective oxidation, nature presents methane monooxygenase as a benchmark. This enzyme selectively consumes methane by breaking it down into methanol. To assemble an active site similar to monooxygenase, the literature reports Cu-ZSM-5, Fe-ZSM-5, and Cu-MOR, using zeolites and systems like CeO2/Cu2O/Cu. However, the trade-off between methane activation and methanol selectivity remains a challenge. Density functional theory (DFT) calculations and spectroscopic studies indicate catalyst reducibility, oxygen mobility, and water as co-feed as primary factors that can assist in enabling higher selectivity. The use of chemical looping can further improve selectivity. However, in all systems, improvements in productivity per cycle are required in order to meet the economical/industrial standards.

Beneath the surface: Evolution of methane activity in the bacterial multicomponent monooxygenases

The bacterial multicomponent monooxygenase (BMM) family has evolved to oxidise a wide array of hydrocarbon substrates of importance to environmental emissions and biotechnology: foremost amongst these is methane, which requires among the most powerful oxidant in biology to activate. To understand how the BMM evolved methane oxidation activity, we investigated the changes in the enzyme family at different levels: operonic, phylogenetic analysis of the catalytic hydroxylase, subunit or folding factor presence, and sequence-function analysis across the entirety of the BMM phylogeny. Our results show that the BMM evolution of new activities was enabled by incremental increases in oxidative power of the active site, and these occur in multiple branches of the hydroxylase phylogenetic tree. While the hydroxylase primary sequence changes that resulted in increased oxidative power of the enzyme appear to be minor, the principle evolutionary advances enabling methane activity occurred in the other components of the BMM complex and in the recruitment of stability proteins. We propose that enzyme assembly and stabilization factors have independently-evolved multiple times in the BMM family to support enzymes that oxidise increasingly difficult substrates. Herein, we show an important example of evolution of catalytic function where modifications to the active site and substrate accessibility, which are the usual focus of enzyme evolution, are overshadowed by broader scale changes to structural stabilization and non-catalytic unit development. Retracing macroscale changes during enzyme evolution, as demonstrated here, should find ready application to other enzyme systems and in protein design.

High-Resolution Extended X-ray Absorption Fine Structure Analysis Provides Evidence for a Longer Fe···Fe Distance in the Q Intermediate of Methane Monooxygenase

Despite decades of intense research, the core structure of the methane C-H bond breaking diiron(IV) intermediate, Q, of soluble methane monooxygenase remains controversial, with conflicting reports supporting either a "diamond" diiron core structure or an open core structure. Early extended X-ray absorption fine structure (EXAFS) data assigned a short 2.46 Å Fe-Fe distance to Q (Shu et al. Science 1997, 275, 515 ) that is inconsistent with several theoretical studies and in conflict with our recent high-resolution Fe K-edge X-ray absorption spectroscopy (XAS) studies (Castillo et al. J. Am. Chem. Soc. 2017, 139, 18024 ). Herein, we revisit the EXAFS of Q using high-energy resolution fluorescence-detected extended X-ray absorption fine structure (HERFD-EXAFS) studies. The present data show no evidence for a short Fe-Fe distance, but rather a long 3.4 Å diiron distance, as observed in open core synthetic model complexes. The previously reported 2.46 Å feature plausibly arises from a background metallic iron contribution from the experimental setup, which is eliminated in HERFD-EXAFS due to the increased selectivity. Herein, we explore the origin of the short diiron feature in partial-fluorescent yield EXAFS measurements and discuss the diagnostic features of background metallic scattering contribution to the EXAFS of dilute biological samples. Lastly, differences in sample preparation and resultant sample inhomogeneity in rapid-freeze quenched samples for EXAFS analysis are discussed. The presented approaches have broad implications for EXAFS studies of all dilute iron-containing samples.

Dioxygen Activation by Nonheme Diiron Enzymes: Diverse Dioxygen Adducts, High-Valent Intermediates, and Related Model Complexes

A growing subset of metalloenzymes activates dioxygen with nonheme diiron active sites to effect substrate oxidations that range from the hydroxylation of methane and the desaturation of fatty acids to the deformylation of fatty aldehydes to produce alkanes and the six-electron oxidation of aminoarenes to nitroarenes in the biosynthesis of antibiotics. A common feature of their reaction mechanisms is the formation of O2 adducts that evolve into more reactive derivatives such as diiron(II,III)-superoxo, diiron(III)-peroxo, diiron(III,IV)-oxo, and diiron(IV)-oxo species, which carry out particular substrate oxidation tasks. In this review, we survey the various enzymes belonging to this unique subset and the mechanisms by which substrate oxidation is carried out. We examine the nature of the reactive intermediates, as revealed by X-ray crystallography and the application of various spectroscopic methods and their associated reactivity. We also discuss the structural and electronic properties of the model complexes that have been found to mimic salient aspects of these enzyme active sites. Much has been learned in the past 25 years, but key questions remain to be answered.

Unprecedented (μ-1,1-Peroxo)diferric Structure for the Ambiphilic Orange Peroxo Intermediate of the Nonheme N-Oxygenase CmlI

The final step in the biosynthesis of the antibiotic chloramphenicol is the oxidation of an aryl-amine substrate to an aryl-nitro product catalyzed by the N-oxygenase CmlI in three two-electron steps. The CmlI active site contains a diiron cluster ligated by three histidine and four glutamate residues and activates dioxygen to perform its role in the biosynthetic pathway. It was previously shown that the active oxidant used by CmlI to facilitate this chemistry is a peroxo-diferric intermediate (CmlIP). Spectroscopic characterization demonstrated that the peroxo binding geometry of CmlIP is not consistent with the μ-1,2 mode commonly observed in nonheme diiron systems. Its geometry was tentatively assigned as μ-η2:η1 based on comparison with resonance Raman (rR) features of mixed-metal model complexes in the absence of appropriate diiron models. Here, X-ray absorption spectroscopy (XAS) and rR studies have been used to establish a refined structure for the diferric cluster of CmlIP. The rR experiments carried out with isotopically labeled water identified the symmetric and asymmetric vibrations of an Fe-O-Fe unit in the active site at 485 and 780 cm-1, respectively, which was confirmed by the 1.83 Å Fe-O bond observed by XAS. In addition, a unique Fe···O scatterer at 2.82 Å observed from XAS analysis is assigned as arising from the distal O atom of a μ-1,1-peroxo ligand that is bound symmetrically between the irons. The (μ-oxo)(μ-1,1-peroxo)diferric core structure associated with CmlIP is unprecedented among diiron cluster-containing enzymes and corresponding biomimetic complexes. Importantly, it allows the peroxo-diferric intermediate to be ambiphilic, acting as an electrophilic oxidant in the initial N-hydroxylation of an arylamine and then becoming a nucleophilic oxidant in the final oxidation of an aryl-nitroso intermediate to the aryl-nitro product.

A tale of two methane monooxygenases

Methane monooxygenase (MMO) enzymes activate O2 for oxidation of methane. Two distinct MMOs exist in nature, a soluble form that uses a diiron active site (sMMO) and a membrane-bound form with a catalytic copper center (pMMO). Understanding the reaction mechanisms of these enzymes is of fundamental importance to biologists and chemists, and is also relevant to the development of new biocatalysts. The sMMO catalytic cycle has been elucidated in detail, including O2 activation intermediates and the nature of the methane-oxidizing species. By contrast, many aspects of pMMO catalysis remain unclear, most notably the nuclearity and molecular details of the copper active site. Here, we review the current state of knowledge for both enzymes, and consider pMMO O2 activation intermediates suggested by computational and synthetic studies in the context of existing biochemical data. Further work is needed on all fronts, with the ultimate goal of understanding how these two remarkable enzymes catalyze a reaction not readily achieved by any other metalloenzyme or biomimetic compound.

Structural differences between the active sites of the Ni-A and Ni-B states of the [NiFe] hydrogenase: an approach by quantum chemistry and single crystal ENDOR spectroscopy

The two resting forms of the active site of [NiFe] hydrogenase, Ni-A and Ni-B, have significantly different activation kinetics, but reveal nearly identical spectroscopic features which suggest the two states exhibit subtle structural differences. Previous studies have indicated that the states differ by the identity of the bridging ligand between Ni and Fe; proposals include OH(-), OOH(-), H2O, O(2-), accompanied by modified cysteine residues. In this study, we use single crystal ENDOR spectroscopy and quantum chemical calculations within the framework of density functional theory (DFT) to calculate vibrational frequencies, (1)H and (17)O hyperfine coupling constants and g values with the aim to compare these data to experimental results obtained by crystallography, FTIR and EPR/ENDOR spectroscopy. We find that the Ni-A and Ni-B states are constitutional isomers that differ in their fine structural details. Calculated vibrational frequencies for the cyano and carbonyl ligands and (1)H and (17)O hyperfine coupling constants indicate that the bridging ligand in both Ni-A and Ni-B is indeed an OH(-) ligand. The difference in the isotropic hyperfine coupling constants of the β-CH2 protons of Cys-549 is sensitive to the orientation of Cys-549; a difference of 0.5 MHz is observed experimentally for Ni-A and 1.9 MHz for Ni-B, which results from a rotation of 7 degrees about the Cα-Cβ-Sγ-Ni dihedral angle. Likewise, the difference of the intermediate g value is correlated with a rotation of Cys-546 of about 10 degrees.

Enzymatic oxidation of methane

Methane monooxygenases (MMOs) are enzymes that catalyze the oxidation of methane to methanol in methanotrophic bacteria. As potential targets for new gas-to-liquid methane bioconversion processes, MMOs have attracted intense attention in recent years. There are two distinct types of MMO, a soluble, cytoplasmic MMO (sMMO) and a membrane-bound, particulate MMO (pMMO). Both oxidize methane at metal centers within a complex, multisubunit scaffold, but the structures, active sites, and chemical mechanisms are completely different. This Current Topic review article focuses on the overall architectures, active site structures, substrate reactivities, protein-protein interactions, and chemical mechanisms of both MMOs, with an emphasis on fundamental aspects. In addition, recent advances, including new details of interactions between the sMMO components, characterization of sMMO intermediates, and progress toward understanding the pMMO metal centers are highlighted. The work summarized here provides a guide for those interested in exploiting MMOs for biotechnological applications.

Structure of the key species in the enzymatic oxidation of methane to methanol

Methane monooxygenase (MMO) catalyses the O2-dependent conversion of methane to methanol in methanotrophic bacteria, thereby preventing the atmospheric egress of approximately one billion tons of this potent greenhouse gas annually. The key reaction cycle intermediate of the soluble form of MMO (sMMO) is termed compound Q (Q). Q contains a unique dinuclear Fe(IV) cluster that reacts with methane to break an exceptionally strong 105 kcal mol(-1) C-H bond and insert one oxygen atom. No other biological oxidant, except that found in the particulate form of MMO, is capable of such catalysis. The structure of Q remains controversial despite numerous spectroscopic, computational and synthetic model studies. A definitive structural assignment can be made from resonance Raman vibrational spectroscopy but, despite efforts over the past two decades, no vibrational spectrum of Q has yet been obtained. Here we report the core structures of Q and the following product complex, compound T, using time-resolved resonance Raman spectroscopy (TR(3)). TR(3) permits fingerprinting of intermediates by their unique vibrational signatures through extended signal averaging for short-lived species. We report unambiguous evidence that Q possesses a bis-μ-oxo diamond core structure and show that both bridging oxygens originate from O2. This observation strongly supports a homolytic mechanism for O-O bond cleavage. We also show that T retains a single oxygen atom from O2 as a bridging ligand, while the other oxygen atom is incorporated into the product. Capture of the extreme oxidizing potential of Q is of great contemporary interest for bioremediation and the development of synthetic approaches to methane-based alternative fuels and chemical industry feedstocks. Insight into the formation and reactivity of Q from the structure reported here is an important step towards harnessing this potential.

Reactivity of the binuclear non-heme iron active site of Δ⁹ desaturase studied by large-scale multireference ab initio calculations

The results of density matrix renormalization group complete active space self-consistent field (DMRG-CASSCF) and second-order perturbation theory (DMRG-CASPT2) calculations are presented on various structural alternatives for the O-O and first C-H activating step of the catalytic cycle of the binuclear nonheme iron enzyme Δ(9) desaturase. This enzyme is capable of inserting a double bond into an alkyl chain by double hydrogen (H) atom abstraction using molecular O2. The reaction step studied here is presumably associated with the highest activation barrier along the full pathway; therefore, its quantitative assessment is of key importance to the understanding of the catalysis. The DMRG approach allows unprecedentedly large active spaces for the explicit correlation of electrons in the large part of the chemically important valence space, which is apparently conditio sine qua non for obtaining well-converged reaction energetics. The derived reaction mechanism involves protonation of the previously characterized 1,2-μ peroxy Fe(III)Fe(III) (P) intermediate to a 1,1-μ hydroperoxy species, which abstracts an H atom from the C10 site of the substrate. An Fe(IV)-oxo unit is generated concomitantly, supposedly capable of the second H atom abstraction from C9. In addition, several popular DFT functionals were compared to the computed DMRG-CASPT2 data. Notably, many of these show a preference for heterolytic C-H cleavage, erroneously predicting substrate hydroxylation. This study shows that, despite its limitations, DMRG-CASPT2 is a significant methodological advancement toward the accurate computational treatment of complex bioinorganic systems, such as those with the highly open-shell diiron active sites.

On Predicting Mössbauer Parameters of Iron-Containing Molecules with Density-Functional Theory

The performance of six frequently used density functional theory (DFT) methods (RPBE, OLYP, TPSS, B3LYP, B3LYP*, and TPSSh) in the prediction of Mössbauer isomer shifts(δ) and quadrupole splittings (ΔEQ) is studied for an extended and diverse set of Fe complexes. In addition to the influence of the applied density functional and the type of the basis set, the effect of the environment of the molecule, approximated with the conducting-like screening solvation model (COSMO) on the computed Mössbauer parameters, is also investigated. For the isomer shifts the COSMO-B3LYP method is found to provide accurate δ values for all 66 investigated complexes, with a mean absolute error (MAE) of 0.05 mm s-1 and a maximum deviation of 0.12 mm s-1. Obtaining accurate ΔEQ values presents a bigger challenge; however, with the selection of an appropriate DFT method, a reasonable agreement can be achieved between experiment and theory. Identifying the various chemical classes of compounds that need different treatment allowed us to construct a recipe for ΔEQ calculations; the application of this approach yields a MAE of 0.12 mm s-1 (7% error) and a maximum deviation of 0.55 mm s-1 (17% error). This accuracy should be sufficient for most chemical problems that concern Fe complexes. Furthermore, the reliability of the DFT approach is verified by extending the investigation to chemically relevant case studies which include geometric isomerism, phase transitions induced by variations of the electronic structure (e.g., spin crossover and inversion of the orbital ground state), and the description of electronically degenerate triplet and quintet states. Finally, the immense and often unexploited potential of utilizing the sign of the ΔEQ in characterizing distortions or in identifying the appropriate electronic state at the assignment of the spectral lines is also shown.

Protonation of a peroxodiiron(III) complex and conversion to a diiron(III/IV) intermediate: implications for proton-assisted O-O bond cleavage in nonheme diiron enzymes

Oxygenation of a diiron(II) complex, [Fe(II)(2)(μ-OH)(2)(BnBQA)(2)(NCMe)(2)](2+) [2, where BnBQA is N-benzyl-N,N-bis(2-quinolinylmethyl)amine], results in the formation of a metastable peroxodiferric intermediate, 3. The treatment of 3 with strong acid affords its conjugate acid, 4, in which the (μ-oxo)(μ-1,2-peroxo)diiron(III) core of 3 is protonated at the oxo bridge. The core structures of 3 and 4 are characterized in detail by UV-vis, Mössbauer, resonance Raman, and X-ray absorption spectroscopies. Complex 4 is shorter-lived than 3 and decays to generate in ~20% yield of a diiron(III/IV) species 5, which can be identified by electron paramagnetic resonance and Mössbauer spectroscopies. This reaction sequence demonstrates for the first time that protonation of the oxo bridge of a (μ-oxo)(μ-1,2-peroxo)diiron(III) complex leads to cleavage of the peroxo O-O bond and formation of a high-valent diiron complex, thereby mimicking the steps involved in the formation of intermediate X in the activation cycle of ribonucleotide reductase.

New insights into the mechanism of electron transfer within flavohemoglobins: tunnelling pathways, packing density, thermodynamic and kinetic analyses

Flavohemoglobins (FlavoHb) are metalloenzymes catalyzing the reaction of nitric oxide dioxygenation. The iron cation of the heme group needs to be preliminarily reduced to the ferrous state to be catalytically competent. This reduction is triggered by a flavin adenine dinucleotide (FAD) prosthetic group which is localized in a distinct domain of the protein. In this paper we obtain new insights into the internal long range electron transfer (over ca. 12 Å) using a combination of experimental and computational approaches. Employing a time-resolved pulse radiolysis technique we report the first direct measurement of the FADH˙→ HemeFe(III) electron transfer rate. A rate constant of (6.8 ± 0.5) × 10(3) s(-1) is found. A large panel of computational approaches are used to provide the first estimation of the thermodynamic characteristics of the internal electron transfer step within flavoHb: both the driving force and the reorganization energy are estimated as a function of the protonated state of the flavin semi-quinone. We also report an analysis of the electron pathways involved in the tunnelling of the electron through the aqueous interface between the globin and the flavin domains.

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.

THEORETICAL STUDY ON METHANE HYDROXYLATION BY MIMIC METHANE MONOOXYGENASE WITH bis(μ-OXO)DIMANGANESE CORE

The reaction mechanism for methane hydroxylation catalyzed by mimic methane monooxygenase (MMO) with bis(μ-oxo)dimanganese core has been investigated on the septet and nonet potential energy surfaces by hybrid density functional method B3LYP. The key reactive compound Q of MMO was modeled by trans- (H2CNH)(COOH) Mn(μ-O)2(μ-HCOO)2Mn(H2CNH)(COOH) . The ground state of Q is located on the septet state, which has a diamond-core structure with two Mn(IV) atoms. It is shown that the reaction proceeds via a radical-rebound mechanism, in which the step of C–H cleavage is the rate-determining step both in the gas phase and solution. Furthermore, the reaction may proceed more easily as the polarity of solution is larger. On the other hand, the kinetic isotope effects (KIEs) for H atom abstraction from methane are taken into account on the basis of transition state theory with Wigner tunneling corrections. The mimic MMO with bis(μ-oxo)dimanganese core might be an effective mimic catalyst for methane hydroxylation.

Structure and mechanism of the diiron benzoyl-coenzyme A epoxidase BoxB

The coenzyme A (CoA)-dependent aerobic benzoate metabolic pathway uses an unprecedented chemical strategy to overcome the high aromatic resonance energy by forming the non-aromatic 2,3-epoxybenzoyl-CoA. The crucial dearomatizing reaction is catalyzed by three enzymes, BoxABC, where BoxA is an NADPH-dependent reductase, BoxB is a benzoyl-CoA 2,3-epoxidase, and BoxC is an epoxide ring hydrolase. We characterized the key enzyme BoxB from Azoarcus evansii by structural and Mössbauer spectroscopic methods as a new member of class I diiron enzymes. Several family members were structurally studied with respect to the diiron center architecture, but no structure of an intact diiron enzyme with its natural substrate has been reported. X-ray structures between 1.9 and 2.5 Å resolution were determined for BoxB in the diferric state and with bound substrate benzoyl-CoA in the reduced state. The substrate-bound reduced state is distinguished from the diferric state by increased iron-ligand distances and the absence of directly bridging groups between them. The position of benzoyl-CoA inside a 20 Å long channel and the position of the phenyl ring relative to the diiron center are accurately defined. The C2 and C3 atoms of the phenyl ring are closer to one of the irons. Therefore, one oxygen of activated O(2) must be ligated predominantly to this proximate iron to be in a geometrically suitable position to attack the phenyl ring. Consistent with the observed iron/phenyl geometry, BoxB stereoselectively should form the 2S,3R-epoxide. We postulate a reaction cycle that allows a charge delocalization because of the phenyl ring and the electron-withdrawing CoA thioester.

Insights into the different dioxygen activation pathways of methane and toluene monooxygenase hydroxylases

The methane and toluene monooxygenase hydroxylases (MMOH and TMOH, respectively) have almost identical active sites, yet the physical and chemical properties of their oxygenated intermediates, designated P*, H(peroxo), Q, and Q* in MMOH and ToMOH(peroxo) in a subclass of TMOH, ToMOH, are substantially different. We review and compare the structural differences in the vicinity of the active sites of these enzymes and discuss which changes could give rise to the different behavior of H(peroxo) and Q. In particular, analysis of multiple crystal structures reveals that T213 in MMOH and the analogous T201 in TMOH, located in the immediate vicinity of the active site, have different rotatory configurations. We study the rotational energy profiles of these threonine residues with the use of molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) computational methods and put forward a hypothesis according to which T213 and T201 play an important role in the formation of different types of peroxodiiron(III) species in MMOH and ToMOH. The hypothesis is indirectly supported by the QM/MM calculations of the peroxodiiron(III) models of ToMOH and the theoretically computed Mössbauer spectra. It also helps explain the formation of two distinct peroxodiiron(III) species in the T201S mutant of ToMOH. Additionally, a role for the ToMOD regulatory protein, which is essential for intermediate formation and protein functioning in the ToMO system, is advanced. We find that the low quadrupole splitting parameter in the Mössbauer spectrum observed for a ToMOH(peroxo) intermediate can be explained by protonation of the peroxo moiety, possibly stabilized by the T201 residue. Finally, similarities between the oxygen activation mechanisms of the monooxygenases and cytochrome P450 are discussed.

Dioxygen activation in soluble methane monooxygenase

The controlled oxidation of methane to methanol is a chemical transformation of great value, particularly in the pursuit of alternative fuels, but the reaction remains underutilized industrially because of inefficient and costly synthetic procedures. In contrast, methane monooxygenase enzymes (MMOs) from methanotrophic bacteria achieve this chemistry efficiently under ambient conditions. In this Account, we discuss the first observable step in the oxidation of methane at the carboxylate-bridged diiron active site of the soluble MMO (sMMO), namely, the reductive activation of atmospheric O(2). The results provide benchmarks against which the dioxygen activation mechanisms of other bacterial multicomponent monooxygenases can be measured. Molecular oxygen reacts rapidly with the reduced diiron(II) cen-ter of the hydroxylase component of sMMO (MMOH). The first spectroscopically characterized intermediate that results from this process is a peroxodiiron(III) species, P*, in which the iron atoms have identical environments. P* converts to a second peroxodiiron(III) unit, H(peroxo), in a process accompanied by the transfer of a proton, probably with the assistance of a residue near the active site. Proton-promoted O-O bond scission and rearrangement of the diiron core then leads to a diiron(IV) unit, termed Q, that is directly responsible for the oxidation of methane to methanol. In one section of this Account, we provide a detailed discussion of these processes, with particular emphasis on possible structures of the intermediates. The geometries of P* and H(peroxo) are currently unknown, and recent synthetic modeling chemistry has highlighted the need for further structural characterization of Q, currently assigned as a di(μ-oxo)diiron(IV) "diamond core." In another section of the Account, we discuss in detail proton transfer during the O(2) activation events. The role of protons in promoting O-O bond cleavage, thereby initiating the conversion of H(peroxo) to Q, was previously a controversial topic. Recent studies of the mechanism, covering a range of pH values and in D(2)O instead of H(2)O, confirmed conclusively that the transfer of protons, possibly at or near the active site, is necessary for both P*-to-H(peroxo) and H(peroxo)-to-Q conversions. Specific mechanistic insights into these processes are provided. In the final section of the Account, we present our view of experiments that need to be done to further define crucial aspects of sMMO chemistry. Here our goal is to detail the challenges that we and others face in this research, particularly with respect to some long-standing questions about the system, as well as approaches that might be used to solve them.

DFT calculations for intermediate and active states of the diiron center with a tryptophan or tyrosine radical in Escherichia coli ribonucleotide reductase

Class Ia ribonucleotide reductase subunit R2 contains a diiron active site. In this paper, active-site models for the intermediate X-Trp48(•+) and X-Tyr122(•), the active Fe(III)Fe(III)-Tyr122(•), and the met Fe(III)Fe(III) states of Escherichia coli R2 are studied, using broken-symmetry density functional theory incorporated with the conductor-like screening solvation model. Different structural isomers and different protonation states have been explored. Calculated geometric, energetic, Mössbauer, hyperfine, and redox properties are compared with available experimental data. Feasible detailed structures of these intermediate and active states are proposed. Asp84 and Trp48 are most likely the main contributing residues to the result that the transient Fe(IV)Fe(IV) state is not observed in wild-type class Ia E. coli R2. Asp84 is proposed to serve as a proton-transfer conduit between the diiron cluster and Tyr122 in both the tyrosine radical activation pathway and the first steps of the catalytic proton-coupled electron-transfer pathway. Proton-coupled and simple redox potential calculations show that the kinetic control of proton transfer to Tyr122(•) plays a critical role in preventing reduction from the active Fe(III)Fe(III)-Tyr122(•) state to the met state, which is potentially the reason why Tyr122(•) in the active state can be stable over a very long period.

Modeling the cis-oxo-labile binding site motif of non-heme iron oxygenases: water exchange and oxidation reactivity of a non-heme iron(IV)-oxo compound bearing a tripodal tetradentate ligand

The spectroscopic and chemical characterization of a new synthetic non-heme iron(IV)-oxo species [Fe(IV)(O)((Me,H) Pytacn)(S)](2+) (2, (Me,H)Pytacn=1-(2'-pyridylmethyl)-4,7-dimethyl-1,4,7-triazacyclononane, S=CH(3)CN or H(2)O) is described. Complex 2 was prepared by reaction of [Fe(II)(CF(3)SO(3))(2)((Me,H) Pytacn)] (1) with peracetic acid. Complex 2 bears a tetradentate N(4) ligand that leaves two cis sites available for binding an oxo group and a second external ligand but, unlike the related iron(IV)-oxo species with tetradentate ligands, it is remarkably stable at room temperature (t(1/2)>2 h at 288 K). Its ability to exchange the oxygen atom of the oxo ligand with water has been analyzed in detail by means of kinetic studies, and a mechanism is proposed on the basis of DFT calculations. Hydrogen-atom abstraction from C-H bonds and oxygen-atom transfer to sulfides by 2 have also been studied. Despite its thermal stability, 2 proves to be a very powerful oxidant that is capable of breaking the strong C-H bond of cyclohexane (bond dissociation energy=99.3 kcal mol(-1)).

The prediction of Fe Mössbauer parameters by the density functional theory: a benchmark study

We report the performance of eight density functionals (B3LYP, BPW91, OLYP, O3LYP, M06, M06-2X, PBE, and SVWN5) in two Gaussian basis sets (Wachters and Partridge-1 on iron atoms; cc-pVDZ on the rest of atoms) for the prediction of the isomer shift (IS) and the quadrupole splitting (QS) parameters of Mössbauer spectroscopy. Two sources of geometry (density functional theory-optimized and X-ray) are used. Our data set consists of 31 iron-containing compounds (35 signals), the Mössbauer spectra of which were determined at liquid helium temperature and where the X-ray geometries are known. Our results indicate that the larger and uncontracted Partridge-1 basis set produces slightly more accurate linear correlations of electronic density used for the prediction of IS and noticeably more accurate results for the QS parameter. We confirm and discuss the earlier observation of Noodleman and co-workers that different oxidation states of iron produce different IS calibration lines. The B3LYP and O3LYP functionals have the lowest errors for either IS or QS. BPW91, OLYP, PBE, and M06 have a mixed success whereas SVWN5 and M06-2X demonstrate the worst performance. Finally, our calibrations and conclusions regarding the best functional to compute the Mössbauer characteristics are applied to candidate structures for the peroxo and Q intermediates of the enzyme methane monooxygenase hydroxylase (MMOH), and compared to experimental data in the literature.

Oxidation reactions performed by soluble methane monooxygenase hydroxylase intermediates H(peroxo) and Q proceed by distinct mechanisms

Soluble methane monooxygenase is a bacterial enzyme that converts methane to methanol at a carboxylate-bridged diiron center with exquisite control. Because the oxidizing power required for this transformation is demanding, it is not surprising that the enzyme is also capable of hydroxylating and epoxidizing a broad range of hydrocarbon substrates in addition to methane. In this work we took advantage of this promiscuity of the enzyme to gain insight into the mechanisms of action of H(peroxo) and Q, two oxidants that are generated sequentially during the reaction of reduced protein with O(2). Using double-mixing stopped-flow spectroscopy, we investigated the reactions of the two intermediate species with a panel of substrates of varying C-H bond strength. Three classes of substrates were identified according to the rate-determining step in the reaction. We show for the first time that an inverse trend exists between the rate constant of reaction with H(peroxo) and the C-H bond strength of the hydrocarbon examined for those substrates in which C-H bond activation is rate-determining. Deuterium kinetic isotope effects revealed that reactions performed by Q, but probably not H(peroxo), involve extensive quantum mechanical tunneling. This difference sheds light on the observation that H(peroxo) is not a sufficiently potent oxidant to hydroxylate methane, whereas Q can perform this reaction in a facile manner. In addition, the reaction of H(peroxo) with acetonitrile appears to proceed by a distinct mechanism in which a cyanomethide anionic intermediate is generated, bolstering the argument that H(peroxo) is an electrophilic oxidant that operates via two-electron transfer chemistry.

Density functional theory analysis of structure, energetics, and spectroscopy for the Mn-Fe active site of Chlamydia trachomatis ribonucleotide reductase in four oxidation states

Models for the Mn-Fe active site structure of ribonucleotide reductase (RNR) from pathogenic bacteria Chlamydia trachomatis (Ct) in different oxidation states have been studied in this paper, using broken-symmetry density functional theory (DFT) incorporated with the conductor like screening (COSMO) solvation model and also with finite-difference Poisson-Boltzmann self-consistent reaction field (PB-SCRF) calculations. The detailed structures for the reduced Mn(II)-Fe(II), the met Mn(III)-Fe(III), the oxidized Mn(IV)-Fe(III) and the superoxidized Mn(IV)-Fe(IV) states are predicted. The calculated properties, including geometries, (57)Fe Mossbauer isomer shifts and quadrupole splittings, and (57)Fe and (55)Mn electron nuclear double resonance (ENDOR) hyperfine coupling constants, are compared with the available experimental data. The Mössbauer and energetic calculations show that the (mu-oxo, mu-hydroxo) models better represent the structure of the Mn(IV)-Fe(III) state than the di-mu-oxo models. The predicted Mn(IV)-Fe(III) distances (2.95 and 2.98 A) in the (mu-oxo, mu-hydroxo) models are in agreement with the extended X-ray absorption fine structure (EXAFS) experimental value of 2.92 A (Younker et al. J. Am. Chem. Soc. 2008, 130, 15022-15027). The effect of the protein and solvent environment on the assignment of the Mn metal position is examined by comparing the relative energies of alternative mono-Mn(II) active site structures. It is proposed that if the Mn(II)-Fe(II) protein is prepared with prior addition of Mn(II) or with Mn(II) richer than Fe(II), Mn is likely positioned at metal site 2, which is further from Phe127.

Revisiting the mechanism of dioxygen activation in soluble methane monooxygenase from M. capsulatus (Bath): evidence for a multi-step, proton-dependent reaction pathway

Stopped-flow kinetic investigations of soluble methane monooxygenase (sMMO) from M. capsulatus (Bath) have clarified discrepancies that exist in the literature regarding several aspects of catalysis by this enzyme. The development of thorough kinetic analytical techniques has led to the discovery of two novel oxygenated iron species that accumulate in addition to the well-established intermediates H(peroxo) and Q. The first intermediate, P*, is a precursor to H(peroxo) and was identified when the reaction of reduced MMOH and MMOB with O(2) was carried out in the presence of >or=540 microM methane to suppress the dominating absorbance signal due to Q. The optical properties of P* are similar to those of H(peroxo), with epsilon(420) = 3500 M(-1) cm(-1) and epsilon(720) = 1250 M(-1) cm(-1). These values are suggestive of a peroxo-to-iron(III) charge-transfer transition and resemble those of peroxodiiron(III) intermediates characterized in other carboxylate-bridged diiron proteins and synthetic model complexes. The second identified intermediate, Q*, forms on the pathway of Q decay when reactions are performed in the absence of hydrocarbon substrate. Q* does not react with methane, forms independently of buffer composition, and displays a unique shoulder at 455 nm in its optical spectrum. Studies conducted at different pH values reveal that rate constants corresponding to P* decay/H(peroxo) formation and H(peroxo) decay/Q formation are both significantly retarded at high pH and indicate that both events require proton transfer. The processes exhibit normal kinetic solvent isotope effects (KSIEs) of 2.0 and 1.8, respectively, when the reactions are performed in D(2)O. Mechanisms are proposed to account for the observations of these novel intermediates and the proton dependencies of P* to H(peroxo) and H(peroxo) to Q conversion.

Density functional theory for transition metals and transition metal chemistry

We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids.

DFT calculations of comparative energetics and ENDOR/Mössbauer properties for two protonation states of the iron dimer cluster of ribonucleotide reductase intermediate X

Two models (I and II) for the active site structure of class-I ribonucleotide reductase (RNR) intermediate X in subunit R2 have been studied in this paper, using broken-symmetry density functional theory (DFT) incorporated with the conductor like screening (COSMO) solvation model and with the finite-difference Poisson-Boltzmann self-consistent reaction field (PB-SCRF) calculations. Only one of the bridging groups between the two iron centers is different between model-I and model-II. Model-I contains two mu-oxo bridges, while model-II has one bridging oxo and one bridging hydroxo. These are large active site models including up to the fourth coordination shell H-bonding residues. Mössbauer and ENDOR hyperfine property calculations show that model-I is more likely to represent the active site structure of RNR-X. However, energetically our pK(a) calculations at first highly favored the bridging oxo and hydroxo (in model-II) structure of the diiron center rather than having the di-oxo bridge (in model-I). Since the Arg236 and the nearby Lys42, which are very close to the diiron center, are on the protein surface of RNR-R2, it is highly feasible that one or two anion groups in solution would interact with the positively charged side chains of Arg236 and Lys42. The anion group(s) can be a reductant, phosphate, sulfate, nitrate, and other negatively charged groups existing in biological environments or in the buffer of the experiment. Since sulfate ions certainly exist in the buffer of the ENDOR experiment, we have examined the effect of the sulfate (SO(4)(2-), surrounded by explicit water molecules) H-bonding to the side chain of Arg236. We find that when sulfate interacts with Arg236, the carboxylate group of Asp237 tends to be protonated, and once Asp237 is protonated, the Fe(iii)Fe(iv) center in X favors the di-oxo bridge (model-I). This would explain that the ENDOR observed RNR-X active site structure is likely to be represented by model-I rather than model-II.

Quantum cluster size and solvent polarity effects on the geometries and Mössbauer properties of the active site model for ribonucleotide reductase intermediate X: a density functional theory study

In studying the properties of metalloproteins using ab initio quantum mechanical methods, one has to focus on the calculations on the active site. The bulk protein and solvent environment is often neglected, or is treated as a continuum dielectric medium with a certain dielectric constant. The size of the quantum cluster of the active site chosen for calculations can vary by including only the first-shell ligands which are directly bound to the metal centers, or including also the second-shell residues which are adjacent to and normally have H-bonding interactions with the first-shell ligands, or by including also further hydrogen bonding residues. It is not well understood how the size of the quantum cluster and the value of the dielectric constant chosen for the calculations will influence the calculated properties. In this paper, we have studied three models (A, B, and C) of different sizes for the active site of the ribonucleotide reductase intermediate X, using density functional theory (DFT) OPBE functional with broken-symmetry methodology. Each model is studied in gas-phase and in the conductor-like screening (COSMO) solvation model with different dielectric constants ε = 4, 10, 20, and 80, respectively. All the calculated Fe-ligand geometries, Heisenberg J coupling constants, and the Mössbauer isomer shifts, quadrupole splittings, and the (57)Fe, (1)H, and (17)O hyperfine tensors are compared. We find that the calculated isomer shifts are very stable. They are virtually unchanged with respect to the size of the cluster and the dielectric constant of the environment. On the other hand, certain Fe-ligand distances are sensitive to both the size of the cluster and the value of ε. ε = 4, which is normally used for the protein environment, appears too small when studying the diiron active site geometry with only the first-shell ligands as seen by comparisons with larger models.

Insights into the P-to-Q conversion in the catalytic cycle of methane monooxygenase from a synthetic model system

For the catalytic cycle of soluble methane monooxygenase (sMMO), it has been proposed that cleavage of the O–O bond in the (μ-peroxo)diiron(III) intermediate P gives rise to the diiron(IV) intermediate Q with an Fe2(μ–O)2 diamond core, which oxidizes methane to methanol. As a model for this conversion, (μ–oxo) diiron(III) complex 1 ([FeIII2(μ–O)(μ–O2H3)(L)2]3+, L = tris(3,5-dimethyl-4-methoxypyridyl-2-methyl)amine) has been treated consecutively with one eq of H2O2 and one eq of HClO4 to form 3 ([FeIV2(μ–O)2(L)2]4+). In the course of this reaction a new species, 2, can be observed before the protonation step; 2 gives rise to a cationic peak cluster by ESI-MS at m /z 1,399, corresponding to the {[Fe2O3L2H](OTf)2}+ ion in which 1 oxygen atom derives from 1 and the other two originate from H2O2. Mössbauer studies of 2 reveal the presence of two distinct, exchange coupled iron(IV) centers, and EXAFS fits indicate a short Fe–O bond at 1.66 Å and an Fe–Fe distance of 3.32 Å. Taken together, the spectroscopic data point to an HO-FeIV-O-FeIV = O core for 2. Protonation of 2 results in the loss of H2O and the formation of 3. Isotope labeling experiments show that the [FeIV2(μ–O)2] core of 3 can incorporate both oxygen atoms from H2O2. The reactions described here serve as the only biomimetic precedent for the conversion of intermediates P to Q in the sMMO reaction cycle and shed light on how a peroxodiiron(III) unit can transform into an [FeIV2(μ–O)2] core.

CD and MCD studies of the effects of component B variant binding on the biferrous active site of methane monooxygenase

The multicomponent soluble form of methane monooxygenase (sMMO) catalyzes the oxidation of methane through the activation of O 2 at a nonheme biferrous center in the hydroxylase component, MMOH. Reactivity is limited without binding of the sMMO effector protein, MMOB. Past studies show that mutations of specific MMOB surface residues cause large changes in the rates of individual steps in the MMOH reaction cycle. To define the structural and mechanistic bases for these observations, CD, MCD, and VTVH MCD spectroscopies coupled with ligand-field (LF) calculations are used to elucidate changes occurring near and at the MMOH biferrous cluster upon binding of MMOB and the MMOB variants. Perturbations to both the CD and MCD are observed upon binding wild-type MMOB and the MMOB variant that similarly increases O 2 reactivity. MMOB variants that do not greatly increase O 2 reactivity fail to cause one or both of these changes. LF calculations indicate that reorientation of the terminal glutamate on Fe2 reproduces the spectral perturbations in MCD. Although this structural change allows O 2 to bridge the diiron site and shifts the redox active orbitals for good overlap, it is not sufficient for enhanced O 2 reactivity of the enzyme. Binding of the T111Y-MMOB variant to MMOH induces the MCD, but not CD changes, and causes only a small increase in reactivity. Thus, both the geometric rearrangement at Fe2 (observed in MCD) coupled with a more global conformational change that may control O 2 access (probed by CD), induced by MMOB binding, are critical factors in the reactivity of sMMO.