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

Direct Methane Oxidation by Copper- and Iron-Dependent Methane Monooxygenases

Methane is a potent greenhouse gas that contributes significantly to climate change and is primarily regulated in Nature by methanotrophic bacteria, which consume methane gas as their source of energy and carbon, first by oxidizing it to methanol. The direct oxidation of methane to methanol is a chemically difficult transformation, accomplished in methanotrophs by complex methane monooxygenase (MMO) enzyme systems. These enzymes use iron or copper metallocofactors and have been the subject of detailed investigation. While the structure, function, and active site architecture of the copper-dependent particulate methane monooxygenase (pMMO) have been investigated extensively, its putative quaternary interactions, regulation, requisite cofactors, and mechanism remain enigmatic. The iron-dependent soluble methane monooxygenase (sMMO) has been characterized biochemically, structurally, spectroscopically, and, for the most part, mechanistically. Here, we review the history of MMO research, focusing on recent developments and providing an outlook for future directions of the field. Engineered biological catalysis systems and bioinspired synthetic catalysts may continue to emerge along with a deeper understanding of the molecular mechanisms of biological methane oxidation. Harnessing the power of these enzymes will necessitate combined efforts in biochemistry, structural biology, inorganic chemistry, microbiology, computational biology, and engineering.

DFT mechanistic study of biomimetic diiron complex catalyzed dehydrogenation: Unexpected Fe(III)Fe(III)-1,1-μ-hydroperoxy active species for hydride abstraction

The diiron active site is pivotal in catalyzing transformations in both biological and chemical systems. Recently, a range of biomimetic diiron catalysts have been synthesized, drawing inspiration from the active architecture of soluble methane monooxygenase (sMMO). These catalysts have been successfully deployed for the dehydrogenation of indolines, marking a significant advancement in the field. Using density functional theory (DFT) calculations, we have identified a novel mechanistic pathway that governs the dehydrogenation of indolines catalyzed by a biomimetic diiron complex. Specifically, this reaction is facilitated by the transfer of a hybrid atom from the C1 position of the substrate to the distal oxygen atom of the Fe(III)Fe(III)-1,1-μ-hydroperoxy active species. This transfer serves as the rate-limiting step for the heterolytic cleavage of the OO bond, ultimately generating the substrate cation. The mechanism we propose aligns well with mechanistic investigations incorporating both kinetic isotope effect (KIE) measurements and evaluations of stereochemical selectivity. This research contributes to the broader scientific understanding of catalysis involving biomimetic diiron complexes and offers valuable insights into the catalytic behaviors of non-heme diiron metalloenzymes.

Rebound or Cage Escape? The Role of the Rebound Barrier for the Reactivity of Non-Heme High-Valent FeIV =O Species

Owing to their high reactivity and selectivity, variations in the spin ground state and a range of possible pathways, high-valent FeIV =O species are popular models with potential bioinspired applications. An interesting example of a structure-reactivity pattern is the detailed study with five nonheme amine-pyridine pentadentate ligand FeIV =O species, including N4py: [(L1 )FeIV =O]2+ (1), bntpen: [(L2 )FeIV =O]2+ (2), py2 tacn: [(L3 )FeIV =O]2+ (3), and two isomeric bispidine derivatives: [(L4 )FeIV =O]2+ (4) and [(L5 )FeIV =O]2+ (5). In this set, the order of increasing reactivity in the hydroxylation of cyclohexane differs from that with cyclohexadiene as substrate. A comprehensive DFT, ab initio CASSCF/NEVPT2 and DLPNO-CCSD(T) study is presented to untangle the observed patterns. These are well reproduced when both activation barriers for the C-H abstraction and the OH rebound are taken into account. An MO, NBO and deformation energy analysis reveals the importance of π(pyr) → π*xz (FeIII -OH) electron donation for weakening the FeIII -OH bond and thus reducing the rebound barrier. This requires that pyridine rings are oriented perpendicularly to the FeIII -OH bond and this is a subtle but crucial point in ligand design for non-heme iron alkane hydroxylation.

Unravelling the Effect of Acid-Driven Electron Transfer in High-Valent FeIV =O/MnIV =O Species and Its Implications for Reactivity

The electron transfer (ET) step is one of the crucial processes in biochemical redox reactions that occur in nature and has been established as a key step in dictating the reactivity of high-valent metal-oxo species. Although metalloenzymes possessing metal-oxo units at their active site are typically associated with outer-sphere electron transfer (OSET) processes, biomimetic models, in contrast, have been found to manifest either an inner-sphere electron transfer (ISET) or OSET mechanism. This distinction is clearly illustrated through the behaviour of [(N4Py)MnIV (O)]2+ (1) and [(N4Py)FeIV (O)]2+ (2) complexes, where complex 1 showcases an OSET mechanism, while complex 2 exhibits an ISET mechanism, especially evident in their reactions involving C-H bond activation and oxygen atom transfer reactions in the presence of a Lewis/Bronsted acid. However, the precise reason for this puzzling difference remains elusive. This work unveils the origin of the perplexing inner-sphere vs outer-sphere electron transfer process (ISET vs OSET) in [(N4Py)MnIV (O)]2+ (1) and [(N4Py)FeIV (O)]2+ (2) species in the presence of Bronsted acid. The calculations indicate that when the substrate (toluene) approaches both 1 and 2 that is hydrogen bonded with two HOTf molecules (denoted as 1-HOTf and 2-HOTf, respectively), proton transfer from one of the HOTf molecules to the metal-oxo unit is triggered and a simultaneous electron transfer occurs from toluene to the metal centre. Interestingly, the preference for OSET by 1-HOTf is found to originate from the choice of MnIV =O centre to abstract spin-down (β) electron from toluene to its δ(dxy ) orbital. On the other hand, in 2-HOTf, a spin state inversion from triplet to quintet state takes place during the proton (from HOTf) coupled electron transfer (from toluene) preferring a spin-up (α) electron abstraction to its σ* (dz 2 ) orbital mediated by HOTf giving rise to ISET. In addition, 2-HOTf was calculated to possess a larger reorganisation energy, which facilitates the ISET process via the acid. The absence of spin-inversion and smaller reorganisation energy switch the mechanism to OSET for 1-HOTf. Therefore, for the first time, the significance of spin-state and spin-inversion in the electron transfer process has been identified and demonstrated within the realm of high-valent metal-oxo chemistry. This discovery holds implications for the potential involvement of high-valent Mn-oxo species in performing similar transformative processes within Photosystem II.

Nuclear Resonance Vibrational Spectroscopic Definition of the Fe(IV)2 Intermediate Q in Methane Monooxygenase and Its Reactivity

Methanotrophic bacteria utilize the nonheme diiron enzyme soluble methane monooxygenase (sMMO) to convert methane to methanol in the first step of their metabolic cycle under copper-limiting conditions. The structure of the sMMO Fe(IV)₂ intermediate Q responsible for activating the inert C-H bond of methane (BDE = 104 kcal/mol) remains controversial, with recent studies suggesting both "open" and "closed" core geometries for its active site. In this study, we employ nuclear resonance vibrational spectroscopy (NRVS) to probe the geometric and electronic structure of intermediate Q at cryogenic temperatures. These data demonstrate that Q decays rapidly during the NRVS experiment. Combining data from several years of measurements, we derive the NRVS vibrational features of intermediate Q as well as its cryoreduced decay product. A library of 90 open and closed core models of intermediate Q is generated using density functional theory to analyze the NRVS data of Q and its cryoreduced product as well as prior spectroscopic data on Q. Our analysis reveals that a subset of closed core models reproduce these newly acquired NRVS data as well as prior data. The reaction coordinate with methane is also evaluated using both closed and open core models of Q. These studies show that the potent reactivity of Q toward methane resides in the "spectator oxo" of its Fe(IV)₂O₂ core, in contrast to nonheme mononuclear Fe(IV)═O enzyme intermediates that H atoms abstract from weaker C-H bonds.

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.

DFT Study on the Catalytic Activity of ALD-Grown Diiron Oxide Nanoclusters for Partial Oxidation of Methane to Methanol

Using density functional theory (DFT), we studied the catalytic activity of iron oxide nanoclusters that mimic the structure of the active site in the soluble form of methane monooxygenase (sMMO) for the partial oxidation of methane to methanol. Using N₂O as the oxidant, we consider a radical-rebound mechanism and a concerted mechanism for the oxidation of methane on either a bridging oxygen (O_b) or a terminal oxygen (O_t) active site. We find that the radical-rebound pathway is preferred over the concerted pathway by 40-50 kJ/mol, but the desorption of methanol and the regeneration of the oxygen site are found to be the highest barriers for the direct conversion of methane to methanol with these catalysts. As demonstrated by a population analysis, the O_x (x = b or t) site behaves as an oxygen radical during the H abstraction, and the [Fe⁺-O_x^-] site behaves as a Lewis acid-base pair during the concerted C-H cleavage. Molecular orbital decomposition analysis further demonstrates electron transfer during the oxidation and reduction steps of the reaction. High-level multireference calculations were also carried out to further assess the DFT results. Understanding how these systems behave during the proposed reaction pathways provides new insights into how they can be tuned for methane partial oxidation.

Soluble Methane Monooxygenase

Aerobic life is possible because the molecular structure of oxygen (O₂) makes direct reaction with most organic materials at ambient temperatures an exceptionally slow process. Of course, these reactions are inherently very favorable, and they occur rapidly with the release of a great deal of energy at high temperature. Nature has been able to tap this sequestered reservoir of energy with great spatial and temporal selectivity at ambient temperatures through the evolution of oxidase and oxygenase enzymes. One mechanism used by these enzymes for O₂ activation has been studied in detail for the soluble form of the enzyme methane monooxygenase. These studies have revealed the step-by-step process of O₂ activation and insertion into the ultimately stable C-H bond of methane. Additionally, an elegant regulatory mechanism has been defined that enlists size selection and quantum tunneling to allow methane oxidation to occur specifically in the presence of more easily oxidized substrates.

Alkane Oxidation: Methane Monooxygenases, Related Enzymes, and Their Biomimetics

Methane monooxygenases (MMOs) mediate the facile conversion of methane into methanol in methanotrophic bacteria with high efficiency under ambient conditions. Because the selective oxidation of methane is extremely challenging, there is considerable interest in understanding how these enzymes carry out this difficult chemistry. The impetus of these efforts is to learn from the microbes to develop a biomimetic catalyst to accomplish the same chemical transformation. Here, we review the progress made over the past two to three decades toward delineating the structures and functions of the catalytic sites in two MMOs: soluble methane monooxygenase (sMMO) and particulate methane monooxygenase (pMMO). sMMO is a water-soluble three-component protein complex consisting of a hydroxylase with a nonheme diiron catalytic site; pMMO is a membrane-bound metalloenzyme with a unique tricopper cluster as the site of hydroxylation. The metal cluster in each of these MMOs harnesses O₂ to functionalize the C-H bond using different chemistry. We highlight some of the common basic principles that they share. Finally, the development of functional models of the catalytic sites of MMOs is described. These efforts have culminated in the first successful biomimetic catalyst capable of efficient methane oxidation without overoxidation at room temperature.

VTST/MT studies of the catalytic mechanism of C-H activation by transition metal complexes with [Cu2(μ-O2)], [Fe2(μ-O2)] and Fe(IV)-O cores based on DFT potential energy surfaces

High-valent Cu and Fe species, which are generated from dioxygen activation in metalloenzymes, carry out the functionalization of strong C-H bonds. Understanding the atomic details of the catalytic mechanism has long been one of the main objectives of bioinorganic chemistry. Large H/D kinetic isotope effects (KIEs) were observed in the C-H activation by high-valent non-heme Cu or Fe complexes in enzymes and their synthetic models. The H/D KIE depends significantly on the transition state properties, such as structure, energies, frequencies, and shape of the potential energy surface, when the tunneling effect is large. Therefore, theoretical predictions of kinetic parameters such as rate constants and KIEs can provide a reliable link between atomic-level quantum mechanical mechanisms and experiments. The accurate prediction of the tunneling effect is essential to reproduce the kinetic parameters. The rate constants and HD/KIE have been calculated using the variational transition-state theory including multidimensional tunneling based on DFT potential energy surfaces along the reaction coordinate. Excellent agreement was observed between the predicted and experimental results, which assures the validity of the DFT potential energy surfaces and, therefore, the proposed atomic-level mechanisms. The [Cu₂(μ-O)₂], [Fe₂(μ-O)₂], and Fe(IV)-oxo species were employed for C-H activation, and their role as catalysts was discussed at an atomic level.

Pathways for Arene Oxidation in Non-Heme Diiron Enzymes: Lessons from Computational Studies on Benzoyl Coenzyme A Epoxidase

Oxygenation of aromatic rings using O₂ is catalyzed by several non-heme carboxylate-bridged diiron enzymes. In order to provide a general mechanistic description for these reactions, computational studies were carried out at the ONIOM(B3LYP/BP86/Amber) level on the non-heme diiron enzyme benzoyl coenzyme A epoxidase, BoxB. The calculations revealed four possible pathways for attacking the aromatic ring: (a) electrophilic (2e^-) attack by a bis(μ-oxo)-diiron(IV) species (Q pathway); (b) electrophilic (2e^-) attack via the σ* orbital of a μ-η²:η²-peroxo-diiron(III) intermediate (Pσ* pathway); (c) radical (1e^-) attack via the π*-orbital of a superoxo-diiron(II,III) species (Pπ* pathway); (d) radical (1e^-) attack of a partially quenched bis(μ-oxo)-diiron(IV) intermediate (Q' pathway). The results allowed earlier work of de Visser on olefin epoxidation by diiron complexes and QM-cluster studies of Liao and Siegbahn on BoxB to be put into a broader perspective. Parallels with epoxidation using organic peracids were also examined. Specifically for the BoxB enzyme, the Q pathway was found to be the most preferred, but the corresponding bis(μ-oxo)-diiron(IV) species is significantly destabilized and not expected to be directly observable. Epoxidation via the Pσ* pathway represents an energetically somewhat higher lying alternative; possible strategies for experimental discrimination are discussed. The selectivity toward epoxidation is shown to stem from a combination of inherent electronic properties of the thioacyl substituent and enzymatic constraints. Possible implications of the results for toluene monooxygenases are considered as well.

A heterometallic (Fe6Na8) cage-like silsesquioxane: synthesis, structure, spin glass behavior and high catalytic activity

The exotic “Asian Lantern” heterometallic cage silsesquioxane [(PhSiO_1.5)₂₀(FeO_1.5)₆(NaO_0.5)₈(n-BuOH)_9.6(C₇H₈)] (I) was obtained and characterized by X-ray diffraction, EXAFS, topological analyses and DFT calculation.

Successful application of the DBLOC method to the hydroxylation of camphor by cytochrome p450

The activation barrier for the hydroxylation of camphor by cytochrome P450 was computed using a mixed quantum mechanics/molecular mechanics (QM/MM) model of the full protein-ligand system and a fully QM calculation using a cluster model of the active site at the B3LYP/LACVP*/LACV3P** level of theory, which consisted of B3LYP/LACV3P** single point energies computed at B3LYP/LACVP* optimized geometries. From the QM/MM calculation, a barrier height of 17.5 kcal/mol was obtained, while the experimental value was known to be less than or equal to 10 kcal/mol. This process was repeated using the D3 correction for hybrid DFT in order to investigate whether the inadequate treatment of dispersion interaction was responsible for the overestimation of the barrier. While the D3 correction does reduce the computed barrier to 13.3 kcal/mol, it was still in disagreement with experiment. After application of a series of transition metal optimized localized orbital corrections (DBLOC) and without any refitting of parameters, the barrier was further reduced to 10.0 kcal/mol, which was consistent with the experimental results. The DBLOC method to CH bond activation in methane monooxygenase (MMO) was also applied, as a second, independent test. The barrier in MMO was known, by experiment, to be 15.4 kcal/mol. After application of the DBLOC corrections to the MMO barrier compute by B3LYP, in a previous study, and accounting for dispersion with Grimme's D3 method, the unsigned deviation from experiment was improved from 3.2 to 2.3 kcal/mol. These results suggested that the combination of dispersion plus localized orbital corrections could yield significant quantitative improvements in modeling the catalytic chemistry of transition-metal containing enzymes, within the limitations of the statistical errors of the model, which appear to be on the order of approximately 2 kcal/mole.

Substrate-dependent H/D kinetic isotope effects and the role of the di(μ-oxo)diiron(IV) core in soluble methane monooxygenase: a theoretical study

Soluble methane monooxygenase (sMMO) is an enzyme that converts alkanes to alcohols using a di(μ-oxo)diiron(IV) intermediate Q at the active site. Very large kinetic isotope effects (KIEs) indicative of significant tunneling are observed for the hydrogen transfer (H-transfer) of CH4 and CH3 CN; however, a relatively small KIE is observed for CH3NO2. The detailed mechanism of the enzymatic H-transfer responsible for the diverse range of KIEs is not yet fully understood. In this study, variational transition-state theory including the multidimensional tunneling approximation is used to calculate rate constants to predict KIEs based on the quantum-mechanically generated intrinsic reaction coordinates of the H-transfer by the di(μ-oxo)diiron(IV) complex. The results of our study reveal that the role of the di(μ-oxo)diiron(IV) core and the H-transfer mechanism are dependent on the substrate. For CH4 , substrate binding induces an electron transfer from the oxygen to one Fe(IV) center, which in turn makes the μ-O ligand more electrophilic and assists the H-transfer by abstracting an electron from the C-H σ orbital. For CH3CN, the reduction of Fe(IV) to Fe(III) occurs gradually with substrate binding and H-transfer. The charge density and electrophilicity of the μ-O ligand hardly change upon substrate binding; however, for CH3NO2, there seems to be no electron movement from μ-O to Fe(IV) during the H-transfer. Thus, the μ-O ligand appears to abstract a proton without an electron from the C-H σ orbital. The calculated KIEs for CH4, CH3CN, and CH3NO2 are 24.4, 49.0, and 8.27, respectively, at 293 K, in remarkably good agreement with the experimental values. This study reveals that diverse KIE values originate mainly from tunneling to the same di(μ-oxo)diiron(IV) core for all substrates, and demonstrate that the reaction dynamics are essential for reproducing experimental results and understanding the role of the diiron core for methane oxidation in sMMO.

DFT study of the mechanism for methane hydroxylation by soluble methane monooxygenase (sMMO): effects of oxidation state, spin state, and coordination number

The exact structure of the active site of intermediate Q, the methane-oxidizing species of soluble methane monooxygenase (sMMO), and the reaction mechanism of Q with methane molecule are still not fully clear. To gain further insights into the structure and reaction mechanism, five diiron models of Q that differ in shape, oxidation state, spin state, and coordination number of the two iron centers are studied. Different mechanisms in different spin states were explored. Density functional theory (DFT) calculations show that Fe(III)Fe(IV)(μ-O)(μ-OH) is more reactive than Fe(IV)(2)(μ-O)(2) in the oxygen-rich environment and that the reactivity of the active core of sMMO-Q is not enhanced by converting its oxo bridge into a terminal ligand. A four-coordinated diiron model is the most effective for methane hydroxylation. Both radical and non-radical intermediates are involved in the reactions for the four-coordinated diiron model.

Triptycene-based Bis(benzimidazole) Carboxylate-Bridged Biomimetic Diiron(II) Complexes

A triptycene-based bis(benzimidazole) ester ligand, L3, was designed to enhance the electron donating ability of the heterocyclic nitrogen atoms relative to those of the first generation bis(benzoxazole) analogs, L1 and L2. A convergent synthesis of L3 was designed and executed. Three-component titration experiments using UV-visible spectroscopy revealed that the desired diiron(II) complex could be obtained with a 1:2:1 ratio of L3:Fe(OTf)2(MeCN)2:external carboxylate reactants. X-ray crystallographic studies of two diiron complexes derived in this manner from L3 revealed their formulas to be [Fe2L3(μ-OH)(μ-O2CR)(OTf)2], where R = 2,6-bis(p-tolyl)benzoate (7) or triphenylacetate (8). The structures are similar to that of a diiron complex derived from L1, [Fe2L1(μ-OH)(μ-O2CArTol)(OTf)2] (9) with a notable difference being that, in 7 and 8, the geometry at iron more closely resembles square-pyramidal than trigonal-bipyramidal. Mössbauer spectroscopic analyses of 7 and 8 indicate the presence of high-spin diiron(II) cores. These results demonstrate the importance of substituting benzimidazole for benzoxazole for assembling biomimetic diiron complexes with syn disposition of two N-donor ligands, as found in O2-activating carboxylate-bridged diiron centers in biology.

Nuclear resonance vibrational spectroscopic and computational study of high-valent diiron complexes relevant to enzyme intermediates

High-valent intermediates of binuclear nonheme iron enzymes are structurally unknown despite their importance for understanding enzyme reactivity. Nuclear resonance vibrational spectroscopy combined with density functional theory calculations has been applied to structurally well-characterized high-valent mono- and di-oxo bridged binuclear Fe model complexes. Low-frequency vibrational modes of these high-valent diiron complexes involving Fe motion have been observed and assigned. These are independent of Fe oxidation state and show a strong dependence on spin state. It is important to note that they are sensitive to the nature of the Fe2 core bridges and provide the basis for interpreting parallel nuclear resonance vibrational spectroscopy data on the high-valent oxo intermediates in the binuclear nonheme iron enzymes.

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.

Evolution of strategies to prepare synthetic mimics of carboxylate-bridged diiron protein active sites

We present a comprehensive review of research conducted in our laboratory in pursuit of the long-term goal of reproducing the structures and reactivity of carboxylate-bridged diiron centers used in biology to activate dioxygen for the conversion of hydrocarbons to alcohols and related products. This article describes the evolution of strategies devised to achieve these goals and illustrates the challenges in getting there. Particular emphasis is placed on controlling the geometry and coordination environment of the diiron core, preventing formation of polynuclear iron clusters, maintaining the structural integrity of model complexes during reactions with dioxygen, and tuning the ligand framework to stabilize desired oxygenated diiron species. Studies of the various model systems have improved our understanding of the electronic and physical characteristics of carboxylate-bridged diiron units and their reactivity toward molecular oxygen and organic moieties. The principles and lessons that have emerged from these investigations will guide future efforts to develop more sophisticated diiron protein model complexes.

Characterization of a synthetic peroxodiiron(III) protein model complex by nuclear resonance vibrational spectroscopy

The vibrational spectrum of an η(1),η(1)-1,2-peroxodiiron(III) complex was measured by nuclear resonance vibrational spectroscopy and fit using an empirical force field analysis. Isotopic (18)O(2) labelling studies revealed a feature involving motion of the {Fe(2)(O(2))}(4+) core that was not previously observed by resonance Raman spectroscopy.

Mechanistic studies of reactions of peroxodiiron(III) intermediates in T201 variants of toluene/o-xylene monooxygenase hydroxylase

Site-directed mutagenesis studies of a strictly conserved T201 residue in the active site of toluene/o-xylene monooxygenase hydroxylase (ToMOH) revealed that a single mutation can facilitate kinetic isolation of two distinctive peroxodiiron(III) species, designated T201(peroxo) and ToMOH(peroxo), during dioxygen activation. Previously, we characterized both oxygenated intermediates by UV-vis and Mössbauer spectroscopy, proposed structures from DFT and QM/MM computational studies, and elucidated chemical steps involved in dioxygen activation through the kinetic studies of T201(peroxo) formation. In this study, we investigate the kinetics of T201(peroxo) decay to explore the reaction mechanism of the oxygenated intermediates following O(2) activation. The decay rates of T201(peroxo) were monitored in the absence and presence of external (phenol) or internal (tryptophan residue in an I100W variant) substrates under pre-steady-state conditions. Three possible reaction models for the formation and decay of T201(peroxo) were evaluated, and the results demonstrate that this species is on the pathway of arene oxidation and appears to be in equilibrium with ToMOH(peroxo).

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.

Current challenges of modeling diiron enzyme active sites for dioxygen activation by biomimetic synthetic complexes

This tutorial review describes recent progress in modeling the active sites of carboxylate-rich non-heme diiron enzymes that activate dioxygen to carry out several key reactions in Nature. The chemistry of soluble methane monooxygenase, which catalyzes the selective oxidation of methane to methanol, is of particular interest for (bio)technological applications. Novel synthetic diiron complexes that mimic structural, and, to a lesser extent, functional features of these diiron enzymes are discussed. The chemistry of the enzymes is also briefly summarized. A particular focus of this review is on models that mimic characteristics of the diiron systems that were previously not emphasized, including systems that contain (i) aqua ligands, (ii) different substrates tethered to the ligand framework, (iii) dendrimers attached to carboxylates to mimic the protein environment, (iv) two N-donors in a syn-orientation with respect to the iron-iron vector, and (v) a N-rich ligand environment capable of accessing oxygenated high-valent diiron intermediates.

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.

Modeling the syn disposition of nitrogen donors in non-heme diiron enzymes. Synthesis, characterization, and hydrogen peroxide reactivity of diiron(III) complexes with the syn N-donor ligand H2BPG2DEV

In order to model the syn disposition of histidine residues in carboxylate-bridged non-heme diiron enzymes, we prepared a new dinucleating ligand, H(2)BPG(2)DEV, that provides this geometric feature. The ligand incorporates biologically relevant carboxylate functionalities, which have not been explored as extensively as nitrogen-only analogues. Three novel oxo-bridged diiron(III) complexes, [Fe(2)(mu-O)(H(2)O)(2)(BPG(2)DEV)](ClO(4))(2) (6), [Fe(2)(mu-O)(mu-O(2)CAr(iPrO))(BPG(2)DEV)](ClO(4)) (7), and [Fe(2)(mu-O)(mu-CO(3))(BPG(2)DEV)] (8), were prepared. Single-crystal X-ray structural characterization confirms that two pyridyl groups are bound syn with respect to the Fe-Fe vector in these compounds. The carbonato-bridged complex 8 forms quantitatively from 6 in a rapid reaction with gaseous CO(2) in organic solvents. A common maroon-colored intermediate (lambda(max) = 490 nm; epsilon = 1500 M(-1) cm(-1)) forms in reactions of 6, 7, or 8 with H(2)O(2) and NEt(3) in CH(3)CN/H(2)O solutions. Mass spectrometric analyses of this species, formed using (18)O-labeled H(2)O(2), indicate the presence of a peroxide ligand bound to the oxo-bridged diiron(III) center. The Mossbauer spectrum at 90 K of the EPR-silent intermediate exhibits a quadrupole doublet with delta = 0.58 mm/s and DeltaE(Q) = 0.58 mm/s. The isomer shift is typical for a peroxodiiron(III) species, but the quadrupole splitting parameter is unusually small compared to those of related complexes. These Mossbauer parameters are comparable to those observed for a peroxo intermediate formed in the reaction of reduced toluene/o-xylene monooxygenase hydroxylase with dioxygen. Resonance Raman studies reveal an unusually low-energy O-O stretching mode in the peroxo intermediate that is consistent with a short diiron distance. Although peroxodiiron(III) intermediates generated from 6, 7, and 8 are poor O-atom-transfer catalysts, they display highly efficient catalase activity, with turnover numbers up to 10,000. In contrast to hydrogen peroxide reactions of diiron(III) complexes that lack a dinucleating ligand, the intermediates generated here could be re-formed in significant quantities after a second addition of H(2)O(2), as observed spectroscopically and by mass spectrometry.

Mechanisms of water oxidation catalyzed by ruthenium diimine complexes

(18)O-isotope-labeling studies have led to the conclusion that there exist two major pathways for water oxidation catalyzed by dimeric ruthenium ions of the general type cis, cis-[L2Ru(III)(OH2)]2O(4+). We have proposed that both pathways involve concerted addition of H and OH fragments derived from H 2O to the complexes in their four-electron-oxidized states, i.e., [L2Ru(V)(O)]2O(4+), ultimately generating bound peroxy intermediates that decay with the evolution of O2. The pathways differ primarily in the site of addition of the OH fragment, which is either a ruthenyl O atom or a bipyridine ligand. In the former case, water addition is thought to give rise to a critical intermediate whose structure is L2Ru(IV)(OH)ORu(IV)(OOH)L2(4+); the structures of intermediates involved in the other pathway are less well defined but may involve bipyridine OH adducts of the type L2Ru(V)(O)ORu(IV)(OH)(L(*)OH)L(4+), which could react further to generate unstable dioxetanes or similar endoperoxides. Published experimental and theoretical support for these pathways is reviewed within the broader context of water oxidation catalysis and related reactions reported for other diruthenium and group 8 monomeric diimine-based catalysts. New experiments that are designed to probe the issue of bipyridine ligand "noninnocence" in catalysis are described. Specifically, the relative contributions of the two pathways have been shown to correlate with substituent effects in 4,4'- and 5,5'-substituted bipyridine complexes in a manner consistent with the formation of a reactive OH-adduct intermediate in one of the pathways, and the formation of OH-bipyridine adducts during catalytic turnover has been directly confirmed by optical spectroscopy. Finally, a photosensitized system for catalyzed water oxidation has been developed that allows assessment of the catalytic efficiencies of the complex ions under neutral and alkaline conditions; these studies show that the ions are far better catalysts than had previously been assumed based upon reported catalytic parameters obtained with strong oxidants in acidic media.

Structural Model Studies for the High-Valent Intermediate Q of Methane Monooxygenase from Broken-Symmetry Density Functional Calculations

Mössbauer isomer shift parameters have been obtained for both density functional theory (DFT) OPBE and OLYP functionals by linear regressions between the measured isomer shifts and calculated electron densities at Fe nuclei for a number of Fe(2+,2.5+) and Fe(2.5+,3+,3.5+,4+) complexes grouped separately. The calculated isomer shifts and quadrupole splittings on the sample Fe complexes from OPBE and OLYP functionals are similar to those of PW91 calculations (J. Comput. Chem. 27 (2006) 1292), however the fit parameters from the linear regressions differ between PW91 and OPBE, OLYP. Four models for the active site structure of the hydroxylase component of soluble methane monooxygenase (MMOH) have been studied, using three DFT functionals OPBE, OLYP, and PW91, incorporated with broken-symmetry methodology and the conductor-like screening (COSMO) solvation model. The calculated properties, including optimized geometries, electronic energies, pK(a)'s, Fe net spin populations, and Mössbauer isomer shifts and quadrupole splittings, have been reported and compared with available experimental values. The high-spin antiferromagnetically (AF) coupled Fe(4+) sites are correctly predicted by OPBE and OLYP methods for all active site models. PW91 potential overestimates the Fe-ligand covalencies for some of the models because of spin crossover. Our calculations and data analysis support the structure (our current model II shown in Figure 8) proposed by Friesner and Lippard's group (J. Am. Chem. Soc. 123 (2001) 3836-3837), which contains an Fe(4+)(μ-O)(2)Fe(4+) center, one axial water which also H-bonds to both side chains of Glu243 and Glu114, and one bidentate carboxylate group from the side chain of Glu144, which is likely to represent the active site of MMOH-Q. A new model structure (model IV shown in Figure 9), which has a terminal hydroxo and a protonated His147 which is dissociated from a nearby Fe, is more asymmetric in its Fe(μ-O)(2)Fe diamond core, and is another very good candidate for intermediate Q.

A synthetic precedent for the [FeIV2(mu-O)2] diamond core proposed for methane monooxygenase intermediate Q

Intermediate Q, the methane-oxidizing species of soluble methane monooxygenase, is proposed to have an [Fe(IV)(2)(mu-O)(2)] diamond core. In an effort to obtain a synthetic precedent for such a core, bulk electrolysis at 900 mV (versus Fc(+/0)) has been performed in MeCN at -40 degrees C on a valence-delocalized [Fe(III)Fe(IV)(mu-O)(2)(L(b))(2)](3+) complex (1b) (E(1/2) = 760 mV versus Fc(+/0)). Oxidation of 1b results in the near-quantitative formation of a deep red complex, designated 2b, that exhibits a visible spectrum with lambda(max) at 485 nm (9,800 M(-1).cm(-1)) and 875 nm (2,200 M(-1).cm(-1)). The 4.2 K Mössbauer spectrum of 2b exhibits a quadrupole doublet with delta = -0.04(1) mm.s(-1) and DeltaE(Q) = 2.09(2) mm.s(-1), parameters typical of an iron(IV) center. The Mössbauer patterns observed in strong applied fields show that 2b is an antiferromagnetically coupled diiron(IV) center. Resonance Raman studies reveal the diagnostic vibration mode of the [Fe(2)(mu-O)(2)] core at 674 cm(-1), downshifting 30 cm(-1) upon (18)O labeling. Extended x-ray absorption fine structure (EXAFS) analysis shows two O/N scatterers at 1.78 A and an Fe scatterer at 2.73 A. Based on the accumulated spectroscopic evidence, 2b thus can be formulated as [Fe(IV)(2)(mu-O)(2)(L(b))(2)](4+), the first synthetic complex with an [Fe(IV)(2)(mu-O)(2)] core. A comparison of 2b and its mononuclear analog [Fe(IV)(O)(L(b))(NCMe)](2+) (4b) reveals that 4b is 100-fold more reactive than 2b in oxidizing weak C H bonds. This surprising observation may shed further light on how intermediate Q carries out the hydroxylation of methane.

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.

Resonance Raman spectroscopy of chloroperoxidase compound II provides direct evidence for the existence of an iron(IV)-hydroxide

We report direct evidence for the existence of an iron(IV)-hydroxide. Resonance Raman measurements on chloroperoxidase compound II (CPO-II) reveal an isotope ((18)O and (2)H)-sensitive band at nu(Fe-O) = 565 cm(-1). Preparation of CPO-II in H(2)O using H(2)(18)O(2) results in a red-shift of 22 cm(-1), while preparation of CPO-II in (2)H(2)O using H(2)O(2) results in a red-shift of 13 cm(-1). These values are in good agreement with the isotopic shifts predicted (23 and 12 cm(-1), respectively) for an Fe-OH harmonic oscillator. The measured Fe-O stretching frequency is also in good agreement with the 1.82-A Fe-O bond reported for CPO-II. A Badger's rule analysis of this distance provides an Fe-O stretching frequency of nu(Badger) = 563 cm(-1). We also present X-band electron nuclear double resonance (ENDOR) data for cryoreduced CPO-II. Cryogenic reduction (77 K) of the EPR-silent Fe(IV)OH center in CPO-II results in an EPR-active Fe(III)OH species with a strongly coupled (13.4 MHz) exchangeable proton. Based on comparisons with alkaline myoglobin, we assign this resonance to the hydroxide proton of cryoreduced CPO-II.

Theoretical modeling of the hydroxylation of methane as mediated by the particulate methane monooxygenase

We present here the results of density functional theory (DFT) calculations directed toward elucidation of the CH bond activation mechanism that might be adopted by the particulate methane monooxygenase (pMMO) in the hydroxylation of methane and related small alkanes. In these calculations, we considered three of the most probable models for the transition metal active site mediating the "oxo-transfer": (i) the trinuclear copper cluster bis(mu(3)-oxo)trinuclear copper(II, II, III) complex 1, recently proposed by Chan et al. [S.I. Chan, K.H.-C. Chen, S.S.-F. Yu, C.-L. Chen, S.S.-J. Kuo, Biochemistry 43 (2004) 4421-4430.]; (ii) the most frequently used model complex, bis(mu-oxo)Cu(III)(2) complex 2; and (iii) the mixed-valence bis(mu-oxo)Cu(II)Cu(III) complex 3. The results obtained indicate that the methane hydroxylation chemistry mediated by the trinuclear copper cluster bis(mu(3)-oxo)trinuclear copper(II, II, III) complex 1 offers the most facile pathway for methane hydroxylation, and this model yields KIE values that are in good agreement with experiment. In this mechanism, the reaction proceeds along a "singlet" potential surface and a "singlet oxene" is directly inserted across a CH bond in a concerted manner. Kinetic isotope effects (k(H)/k(D) or KIE) associated with the concerted oxene insertion process mediated by complex 1 are calculated to be 5.2 at 300K when tunneling effects are included. Overall rate constants for the methane hydroxylation by the three models have been calculated as a function of temperature, and the rates are at least 5-6 orders of magnitude more facile when the chemistry is mediated by complex 1 compared to complex 2 or complex 3.

Theoretical Studies of Manganese and Iron Superoxide Dismutases: Superoxide Binding and Superoxide Oxidation

Density-functional calculations indicate that the second sphere of coordination around the metal centers of manganese and iron superoxide dismutases (MnSODs and FeSODs) plays an important role in the binding of O2(-). In these systems, O2(-) prefers to bind to Mn or Fe in end-on configurations. For human and E. coli MnSODs, the bound O2(-) forms hydrogen bonds with the tyrosine and glutamine amino acid residues in the second sphere of coordination. In the cases of E. coli and T. elongates FeSODs, hydrogen bonding occurs between the bound O2(-) and the tyrosine amino acid only because the glutamine is too far away for an effective bonding interaction. The manner in which the O2(-) binds to the metal center in MnSODs and FeSODs can affect the rate of subsequent protonation and determine the mechanism for the formation of H2O2. Both Mn- and Fe-containing superoxide dismutases contain a metal-bound solvent molecule that has been suggested to be involved in the uptake of a H+ upon reduction of the metal center [Bull, C.; Fee, J. A. J. Am. Chem. Soc. 1985, 107, 3295; Miller, A.-F.; Padmakumar, K.; Sorkin, D. L.; Karapetian, A.; Vance, C. K. J. Inorg. Biochem. 2003, 93, 71]. Using density-functional theory, we confirm this suggestion and show the involvement of the second sphere of coordination in the process. We show that the oxidation of superoxide by Mn- or Fe-containing superoxide dismutases is facilitated by a cooperative effect between superoxide binding, protonation of the OH- bound to the metal, and electron transfer from the superoxide molecule to the oxidized metal. In particular, proton transfer through tyrosine-34 on the absence of a bound superoxide is uphill while, once superoxide is bound, the energetic barrier is lowered. It is this barrier that likely keeps the resting state (Mn(III)SOD) of the enzyme with a bound hydroxide, instead of a water. This work provides a model for the mechanism of reaction of superoxide with the oxidized form of the metal within Mn- and FeSODs.

Synthesis of diethynyltriptycene-linked dipyridyl ligands

[reaction: see text] An efficient route to a new family of dinucleating ligands has been developed. A convergent strategy to these ligands involved dual Sonogashira cross-coupling of 2,3-diethynyltriptycene with a variety of functionally diverse 5-bromopyridines. The resultant ligands were accessed in four steps and 40-50% overall yield from 1,2,4,5-tetrabromobenzene. Synthesis of an imidazole and a quinoline derivative by this method is also described.

AB INITIO QUANTUM CHEMICAL AND MIXED QUANTUM MECHANICS/MOLECULAR MECHANICS (QM/MM) METHODS FOR STUDYING ENZYMATIC CATALYSIS

We describe large scale ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic reactions. First, technical aspects of the methodology are reviewed, including the hybrid density functional theory (DFT) methods that are typically employed for the QM aspect of the calculations, and various approaches to defining the interface between the QM and MM regions in QM/MM approaches. The modeling of the enzymatic catalytic cycle for three examples--methane monooxygenase, cytochrome P450, and triose phosphate isomerase--are discussed in some depth, followed by a brief summary of other systems that have been investigated by ab initio methods over the past several years. Finally, a discussion of the qualitative and quantitative conclusions concerning enzymatic catalysis that are available from modern ab initio approaches is presented, followed by a conclusion briefly summarizing future prospects.

Hydroxylation of C-H bonds at carboxylate-bridged diiron centres

Nature uses carboxylate-bridged diiron centres at the active sites of enzymes that catalyse the selective hydroxylation of hydrocarbons to alcohols. The resting diiron(III) state of the hydroxylase component of soluble methane monooxygenase enzyme is converted by two-electron transfer from an NADH-requiring reductase into the active diiron(II) form, which subsequently reacts with O2 to generate a high-valent diiron(IV) oxo species (Q) that converts CH4 into CH3OH. In this step, C-H bond activation is achieved through a transition state having a linear C...H...O unit involving a bound methyl radical. Kinetic studies of the reaction of Q with substrates CH3X, where X=H, D, CH3, NO2, CN or OH, reveal two classes of reactivity depending upon whether binding to the enzyme or C-H bond activation is rate-limiting. Access of substrates to the carboxylate-bridged diiron active site in the hydroxylase (MMOH) occurs through a series of hydrophobic pockets. In the hydroxylase component of the closely related enzyme toluene/o-xylene monooxygenase (ToMOH), substrates enter through a wide channel in the alpha-subunit of the protein that tracks a course identical to that found in the structurally homologous MMOH. Synthetic models for the carboxylate-bridged diiron centres in MMOH and ToMOH have been prepared that reproduce the stoichiometry and key geometric and physical properties of the reduced and oxidized forms of the proteins. Reactions of the diiron(II) model complexes with dioxygen similarly generate reactive intermediates, including high-valent species capable not only of hydroxylating pendant C-H bonds but also of oxidizing phosphine and sulphide groups.