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.     

X-ray absorption spectroscopic study of the reduced hydroxylases of methane monooxygenase and toluene/o-xylene monooxygenase: differences in active site structure and effects of the coupling proteins MMOB and ToMOD

The diiron active sites of the reduced hydroxylases from methane monooxygenase (MMOH(red)) and toluene/o-xylene monooxygenase (ToMOH(red)) have been investigated by X-ray absorption spectroscopy (XAS). Results of Fe K-edge and extended X-ray absorption fine structure analysis reveal subtle differences between the hydroxylases that may be correlated to access of the active site. XAS data were also recorded for each hydroxylase in the presence of its respective coupling protein. MMOB affects the outer-shell scattering contributions in the diiron site of MMOH(red), whereas ToMOD exerts its main effect on the first-shell ligation of ToMOH(red); it also causes a slight decrease in the Fe-Fe separation. These results provide an initial step toward delineating the differences in structure and reactivity in bacterial multicomponent monooxygenase proteins.     

Dioxygen activation in methane monooxygenase: a theoretical study

Using broken-symmetry unrestricted Density Functional Theory, the mechanism of enzymatic dioxygen activation by the hydroxylase component of soluble methane monooxygenase (MMOH) is determined to atomic detail. After a thorough examination of mechanistic alternatives, an optimal pathway was identified. The diiron(II) state H(red) reacts with dioxygen to give a ferromagnetically coupled diiron(II,III) H(superoxo) structure, which undergoes intersystem crossing to the antiferromagnetic surface and affords H(peroxo), a symmetric diiron(III) unit with a nonplanar mu-eta(2):eta(2)-O(2)(2)(-) binding mode. Homolytic cleavage of the O-O bond yields the catalytically competent intermediate Q, which has a di (mu-oxo)diiron(IV) core. A carboxylate shift involving Glu243 is essential to the formation of the symmetric H(peroxo) and Q structures. Both thermodynamic and kinetic features agree well with experimental data, and computed spin-exchange coupling constants are in accord with spectroscopic values. Evidence is presented for pH-independent decay of H(red) and H(peroxo). Key electron-transfer steps that occur in the course of generating Q from H(red) are also detailed and interpreted. In contrast to prior theoretical studies, a requisite large model has been employed, electron spins and couplings have been treated in a quantitative manner, potential energy surfaces have been extensively explored, and quantitative total energies have been determined along the reaction pathway.     

Mechanism of dioxygen cleavage in tetrahydrobiopterin-dependent amino acid hydroxylases

The reaction mechanism for the formation of the hydroxylating intermediate in aromatic amino acid hydroxylases (i.e., phenylalanine hydroxylase, tyrosine hydroxylase, tryptophan hydroxylase) was investigated by means of hybrid density functional theory. These enzymes use molecular oxygen to hydroxylate both the tetrahydrobiopterin cofactor and the aromatic amino acid. A mechanism is proposed in which dioxygen forms a bridging bond between the cofactor and iron. The product is an iron(II)-peroxy-pterin intermediate, and iron was found to be essential for the catalysis of this step. No stable intermediates involving a pterin radical cation and a superoxide ion O(2)(-) were found on the reaction pathway. Heterolysis of the O-O bond in the iron(II)-peroxy-pterin intermediate is promoted by one of the water molecules coordinated to iron and releases hydroxypterin and the high-valent iron oxo species Fe(IV)=O, which can carry out subsequent hydroxylation of aromatic rings. In the proposed mechanism, the formation of the bridging C-O bond is rate-limiting in the formation of Fe(IV)=O.     

甲烷单加氧酶化学模拟研究进展

本文综述了近年来对甲烷单加氧酶进行结构模拟和催化功能模拟研究工作的进展情况。就结构模拟而言，已合成和表征了很多模型化合物，并与天然酶进行了比较，获得了一些能较好地再现天然酶活性中心结构特征和诺学特征的模型化合物。在催化功能模拟方面，少数体系能将环己烷等经类转化为醇类或酮类。而真正能把甲烷转化成甲醇的体系仅有一例报道，这方面的工作比较薄弱，尚待加强。     

Capture of activated dioxygen intermediates at the copper-active site of a lytic polysaccharide monooxygenase

Metalloproteins perform a diverse array of redox-related reactions facilitated by the increased chemical functionality afforded by their metallocofactors. Lytic polysaccharide monooxygenases (LPMOs) are a class of copper-dependent enzymes that are responsible for the breakdown of recalcitrant polysaccharides via oxidative cleavage at the glycosidic bond. The activated copper-oxygen intermediates and their mechanism of formation remains to be established. Neutron protein crystallography which permits direct visualization of protonation states was used to investigate the initial steps of oxygen activation directly following active site copper reduction in Neurospora crassa LPMO9D. Herein, we cryo-trap an activated dioxygen intermediate in a mixture of superoxo and hydroperoxo states, and we identify the conserved second coordination shell residue His157 as the proton donor. Density functional theory calculations indicate that both superoxo and hydroperoxo active site states are stable. The hydroperoxo formed is potentially an early LPMO catalytic reaction intermediate or the first step in the mechanism of hydrogen peroxide formation in the absence of substrate. We observe that the N-terminal amino group of the copper coordinating His1 remains doubly protonated directly following molecular oxygen reduction by copper. Aided by molecular dynamics and mining minima free energy calculations we establish that the conserved second-shell His161 in MtPMO3* displays conformational flexibility in solution and that this flexibility is also observed, though to a lesser extent, in His157 of NcLPMO9D. The imidazolate form of His157 observed in our structure following oxygen intermediate protonation can be attributed to abolished His157 flexibility due steric hindrance in the crystal as well as the solvent-occluded active site environment due to crystal packing. A neutron crystal structure of NcLPMO9D at low pH further supports occlusion of the active site since His157 remains singly protonated even at acidic conditions.

Superoxo and hydroperoxo intermediates were cryotrapped at the copper active site of lytic polysaccharide monooxygenase using neutron protein crystallography.     

The quest for the particulate methane monooxygenase active site

Particulate methane monooxygenase is a copper-containing, membrane-bound metalloenzyme that converts methane to methanol in Nature. How pMMO accomplishes this difficult reaction under ambient conditions is one of the major unsolved problems in bioinorganic chemistry. Despite considerable research efforts in the past 20 years, the active site of the enzyme remains unknown. We recently solved the first crystal structure of pMMO to 2.8 è resolution, revealing the overall structure, oligomerization state, subunit ratio, and composition and location of the metal centers. Almost none of the key structural features were predicted. In this Perspective, we review the state of knowledge before and after the structure determination, emphasizing elucidation of the pMMO active site.     

A non-radical mechanism for methane hydroxylation at the diiron active site of soluble methane monooxygenase

We propose a non‐radical mechanism for the conversion of methane into methanol by soluble methane monooxygenase (sMMO), the active site of which involves a diiron active center. We assume the active site of the MMOH_Q intermediate, exhibiting direct reactivity with the methane substrate, to be a bis(μ‐oxo)diiron(IV ) complex in which one of the iron atoms is coordinatively unsaturated (five‐coordinate). Is it reasonable for such a diiron complex to be formed in the catalytic reaction of sMMO? The answer to this important question is positive from the viewpoint of energetics in density functional theory (DFT) calculations. Our model thus has a vacant coordination site for substrate methane. If MMOH_Q involves a coordinatively unsaturated iron atom at the active center, methane is effectively converted into methanol in the broken‐symmetry singlet state by a non‐radical mechanism; in the first step a methane C? H bond is dissociated via a four‐centered transition state (TS1) resulting in an important intermediate involving a hydroxo ligand and a methyl ligand, and in the second step the binding of the methyl ligand and the hydroxo ligand through a three‐centered transition state (TS2) results in the formation of a methanol complex. This mechanism is essentially identical to that of the methane–methanol conversion by the bare FeO⁺ complex and relevant transition metal–oxo complexes in the gas phase. Neither radical species nor ionic species are involved in this mechanism. We look in detail at kinetic isotope effects (KIEs) for H atom abstraction from methane on the basis of transition state theory with Wigner tunneling corrections.     

The FeIV 2(μ-O)2 cluster and bridge O2 activation at the active center of methane monooxygenase

The newest experimental data concerning soluble forms of methane monooxygenase (MMO) were analyzed taking into account the bridge mechanism of O₂ activation by this enzyme proposed earlier. The results confirm that the scheme suggested and the structures of the key intermediates are valid and show a basic difference between the mechanisms of activation for heme (cytochrome P-450) and non-heme (MMO) monooxygenases. The X-ray diffraction analysis of MMO allowed us to develop a more detailed scheme that reflects the dynamics of O₂ activation and the role of ligands in this process. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1676–1682, September, 1997.     

Reactions of methane monooxygenase intermediate Q with derivatized methanes

The reactivity of intermediate Q of soluble methane monooxygenase from Methylococcus capsulatus (Bath) with a series of derivatized methanes of general formula CH3-R, where R = CN, NO2, and OH, and their deuterated analogues is described. Despite relatively slow reactions, significant kinetic isotope effects, KIEs, are observed in reactions with acetonitrile and nitromethane; however, none is obtained with methanol. In addition, evidence of a substrate-binding step that occurs prior to substrate oxidation is observed.     

Dioxygen activation at non-heme diiron centers: characterization of intermediates in a mutant form of toluene/o-xylene monooxygenase hydroxylase

We report the generation and characterization of an intermediate in a mutant form of the toluene/o-xylene monooxygenase hydroxylase component from Pseudomonas stutzeri OX1. The reaction of chemically reduced I100W variant in the presence of the coupling protein, ToMOD, with dioxygen was monitored by stopped-flow UV/visible spectroscopy. Rapid-freeze quench (RFQ) samples were also generated for EPR and M?ssbauer spectroscopy. A transient species is observed in the UV/visible spectrum with an absorption maximum at 500 nm. EPR and M?ssbauer spectra of RFQ samples identified this species as a diiron(III,IV) cluster spin-coupled to a neutral W radical. A diamagnetic precursor to the mixed-valent diiron(III,IV) was also observed at an earlier time point, with M?ssbauer parameters typical of high-spin FeIII. We have tentatively assigned this antiferromagnetically coupled diiron(III) intermediate as a peroxo-bridged cluster, and this complex has also been observed in preliminary studies of the wild-type hydroxylase.     

Computational studies of reaction mechanisms of methane monooxygenase and ribonucleotide reductase

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

Hydroxylation mechanism of methane and its derivatives over designed methane monooxygenase model with peroxo dizinc core

The peroxo dizinc Zn(2)O(2) complex Q coordinated by imidazole and carboxylate groups for each Zn center has been designed to model the hydroxylase component of methane monooxygenase (MMO) enzyme, on the basis of the experimentally available structure information of enzyme with divalent zinc ion and the MMO with Fe(2)O(2) core. The reaction mechanism for the hydroxylation of methane and its derivatives catalyzed by Q has been investigated at the B3LYP*/cc-pVTZ, Lanl2tz level in protein solution environment. These hydroxylation reactions proceed via a radical-rebound mechanism, with the rate-determining step of the C-H bond cleavage. This radical-rebound reaction mechanism is analogous to the experimentally available MMOs with diamond Fe(2)O(2) core accompanied by a coordinate number of six for the hydroxylation of methane. The rate constants for the hydroxylation of substrates catalyzed by Q increase along CH(4) < CH(3)F < CH(3)CN ≈ CH(3)NO(2) < CH(3)CH(3). Both the activation strain ΔE(≠)(strain) and the stabilizing interaction ΔE(≠)(int) jointly affect the activation energy ΔE(≠). For the C-H cleavage of substrate CH(3)X, with the decrease of steric shielding for the substituted CH(3)X (X = F > H > CH(3) > NO(2) > CN) attacking the O center in Q, the activation strain ΔE(≠)(strain) decreases, whereas the stabilizing interaction ΔE(≠)(int) increases. It is predicted that the MMO with peroxo dizinc Zn(2)O(2) core should be a promising catalyst for the hydroxylation of methane and its derivatives.     

Reactions of the peroxo intermediate of soluble methane monooxygenase hydroxylase with ethers

Soluble methane monooxygenase (sMMO) isolated from Methylococcus capsulatus (Bath) utilizes a carboxylate-bridged diiron center and dioxygen to catalyze the conversion of methane to methanol. Previous studies revealed that a di(mu-oxo)diiron(IV) intermediate termed Q is responsible for the catalytic activity with hydrocarbons. In addition, the peroxodiiron(III) intermediate (H(peroxo)) that precedes Q formation in the catalytic cycle has been demonstrated to react with propylene, but its reactivity has not been extensively investigated. Given the burgeoning interest in the existence of multiple oxidants in metalloenzymes, a more exhaustive study of the reactivity of H(peroxo) was undertaken. The kinetics of single turnover reactions of the two intermediates with ethyl vinyl ether and diethyl ether were monitored by single- and double-mixing stopped-flow optical spectroscopy. For both substrates, the rate constants for reaction with H(peroxo) are greater than those for Q. An analytical model for explaining the transient kinetics is described and used successfully to fit the observed data. Activation parameters were determined through temperature-dependent studies, and the kinetic isotope effects for the reactions with diethyl ether were measured. The rate constants indicate that H(peroxo) is a more electrophilic oxidant than Q. We propose that H(peroxo) reacts via two-electron transfer mechanisms, and that Q reacts by single-electron transfer steps.     

Structure and catalytic mechanism of methane monooxygenase and approaches to its modelling

The mechanisms of activation of O₂ and CH₄ by methane monooxygenase (MMO) are discussed. A new concept for the catalytic cycle of MMO is suggested, and approaches to its chemical modelling are considered.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 6, pp. 1011–1020, June, 1995.The author is grateful to Academician A. E. Shilov for interest in this work and helpful discussions. The work was financially supported by the Russian Foundation for Basic Research (Project No. 94-03-08529), International Scientific Foundation (Grant REU000), and INTAS (Project No. 93-315).     

Structural and functional model of methane hydroxylase of membrane-bound methane monooxygenase from Methylococcus capsulatus (Bath)

Computer analysis of a wide range of primary sequences showed that -, -, and -peptides of membrane-bound methane hydroxylase contained 2, 7, and 6 transmembrane helices respectively. Conservative amino acid residues participating in complex formation were revealed. The - and -peptides are suggested to contain mononuclear copper ions with the ligand environment mainly consisting of His residues. The Cu sites are located in the hydrophilic region and are responsible for ESR signals. The active site of -peptide in which the activation of O₂ and oxidation of CH₄ occur is localized in the hydrophobic region close to the membrane surface. This site is formed by the amino acid residues of four transmembrane helices and one loop between them and is suggested to be a binuclear Cu—Fe or Fe—Fe center. The Cu site of -peptide transfers electrons to the active site of -peptide, and the Cu site of -peptide is either involved in this process or only stabilizes the protein structure.     

Evidence for oxygen binding at the active site of particulate methane monooxygenase

Particulate methane monooxygenase (pMMO) is an integral membrane metalloenzyme that converts methane to methanol in methanotrophic bacteria. The enzyme consists of three subunits, pmoB, pmoA, and pmoC, organized in an α(3)β(3)γ(3) trimer. Studies of intact pMMO and a recombinant soluble fragment of the pmoB subunit (denoted as spmoB) indicate that the active site is located within the soluble region of pmoB at the site of a crystallographically modeled dicopper center. In this work, we have investigated the reactivity of pMMO and spmoB with oxidants. Upon reduction and treatment of spmoB with O(2) or H(2)O(2) or pMMO with H(2)O(2), an absorbance feature at 345 nm is generated. The energy and intensity of this band are similar to those of the μ-η(2):η(2)-peroxo-Cu(II)(2) species formed in several dicopper enzymes and model compounds. The feature is not observed in inactive spmoB variants in which the dicopper center is disrupted, consistent with O(2) binding to the proposed active site. Reaction of the 345 nm species with CH(4) results in the disappearance of the spectroscopic feature, suggesting that this O(2) intermediate is mechanistically relevant. Taken together, these observations provide strong new support for the identity and location of the pMMO active site.     

Dynamics of alkane hydroxylation at the non-heme diiron center in methane monooxygenase

Semiclassical molecular dynamics simulations have been combined with quantum chemistry calculations to provide detailed modeling of the methane and ethane hydroxylation reactions catalyzed by the hydroxylase enzymes of the soluble methane monooxygenase system. The experimental distribution of enantiomeric alcohols in the reaction of ethanes made chiral by the use of hydrogen isotopes is quantitatively reproduced and explained. The reaction dynamics involve a mixture of concerted and bound radical trajectories, and we characterize each of these reactive channels in detail. Diffusion of the bound radical intermediate at the active site core determines the global rate constant. The results also provide a qualitative rationale for the lack of ring-opened products derived from certain radical clock substrate probes and for the relative rate constants and kinetic isotope effects exhibited by a variety of substrates.     

MnIIILx/t-BuOOH-induced activation of dioxygen for the oxygenation of cyclohexene

Several manganese (III) complexes (Mn^IIIL_x) in combination with tert-butyl hydroperoxide (t-BuOOH) activate dioxygen (O₂) to oxygenate cyclohexene (c-C₆H₁₀) to its ketone, alcohol, and epoxide. The product profiles depend on the ligand and solvent matrix. With picolinate (PA), bipyridine (bpy), or triphenylphosphine oxide (OPPh₃) as the ligand in py/HOAc (2:1 molar ratio) dominant product is the ketone [c-C₆H₈(O)] whereas Schiff–base complexes produce c-C₆H₈(O), c-C₆H₉(OH) and the epoxide in almost equal yields. However, in MeCN c-C₆H₈(O) is the dominant product for all of the complexes.     

Geometric and electronic structure studies of the binuclear nonheme ferrous active site of toluene-4-monooxygenase: parallels with methane monooxygenase and insight into the role of the effector proteins in O2 activation

Multicomponent monooxygenases, which carry out a variety of highly specific hydroxylation reactions, are of great interest as potential biocatalysts in a number of applications. These proteins share many similarities in structure and show a marked increase in O2 reactivity upon addition of an effector component. In this study, circular dichroism (CD), magnetic circular dichroism (MCD), and variable-temperature, variable-field (VTVH) MCD have been used to gain spectroscopic insight into the Fe(II)Fe(II) active site in the hydroxylase component of Toluene-4 monoxygenase (T4moH) and the complex of T4moH bound by its effector protein, T4moD. These results have been correlated to spectroscopic data and density functional theory (DFT) calculations on MmoH and its interaction with MmoB. Together, these data provide further insight into the geometric and electronic structure of these biferrous active sites and, in particular, the perturbation associated with component B/D binding. It is found that binding of the effector protein changes the geometry of one iron center and orientation of its redox active orbital to accommodate the binding of O2 in a bridged structure for efficient 2-electron transfer that can form a peroxo intermediate.