Dioxygen activation in methane monooxygenase: a theoretical study

Superoxide to Peroxide Interconversion in Ni-TMC Complexes: The Significance of Structure and Spin States

A deeper comprehension of the characteristics of metal-superoxide and metal-peroxide chemical species is imperative, considering their pivotal functions in oxygen transport, enzymatic activation, and catalytic oxygenations. O2 activation is mediated by the interconversion of superoxide and peroxide species. Even though there are multiple studies on metal-superoxide and -peroxide intermediates, robust examples of their interconversion processes are scarce synthetically. For example, Ni-superoxide/peroxide complexes have been characterized with N-Tetramethylated Cyclam (TMC) ligands with different ring sizes, i.e., Nickel(II)-superoxide complex is characterized with 14-TMC while Nickel(III)-peroxide complex with 12-TMC. Later, both complexes were obtained with 13-TMC ligand by employing different bases; interestingly, no evidence of interconversion between them was identified. What are the factors influencing these processes and why is this preference? We attempt a computational analysis of this issue and provide arguments based on our conclusions. 2-dimensional potential energy scan is performed on the 12-TMC, 13-TMC, and 14-TMC systems to identify the reaction path connecting superoxide and peroxide species. Analyses indicate that structure and spin states play a significant role in determining the probability of interconversion. The superoxide-peroxide interconversion process appears to be bound by their propensity for distinct structural features and spin states.

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.

Substrate-Triggered μ-Peroxodiiron(III) Intermediate in the 4-Chloro-l-Lysine-Fragmenting Heme-Oxygenase-like Diiron Oxidase (HDO) BesC: Substrate Dissociation from, and C4 Targeting by, the Intermediate

The enzyme BesC from the β-ethynyl-l-serine biosynthetic pathway in Streptomyces cattleya fragments 4-chloro-l-lysine (produced from l-Lysine by BesD) to ammonia, formaldehyde, and 4-chloro-l-allylglycine and can analogously fragment l-Lys itself. BesC belongs to the emerging family of O2-activating non-heme-diiron enzymes with the "heme-oxygenase-like" protein fold (HDOs). Here, we show that the binding of l-Lys or an analogue triggers capture of O2 by the protein's diiron(II) cofactor to form a blue μ-peroxodiiron(III) intermediate analogous to those previously characterized in two other HDOs, the olefin-installing fatty acid decarboxylase, UndA, and the guanidino-N-oxygenase domain of SznF. The ∼5- and ∼30-fold faster decay of the intermediate in reactions with 4-thia-l-Lys and (4RS)-chloro-dl-lysine than in the reaction with l-Lys itself and the primary deuterium kinetic isotope effects (D-KIEs) on decay of the intermediate and production of l-allylglycine in the reaction with 4,4,5,5-[2H4]-l-Lys suggest that the peroxide intermediate or a reversibly connected successor complex abstracts a hydrogen atom from C4 to enable olefin formation. Surprisingly, the sluggish substrate l-Lys can dissociate after triggering intermediate formation, thereby allowing one of the better substrates to bind and react. The structure of apo BesC and the demonstrated linkage between Fe(II) and substrate binding suggest that the triggering event involves an induced ordering of ligand-providing helix 3 (α3) of the conditionally stable HDO core. As previously suggested for SznF, the dynamic α3 also likely initiates the spontaneous degradation of the diiron(III) product cluster after decay of the peroxide intermediate, a trait emerging as characteristic of the nascent HDO family.

Plant monounsaturated fatty acids: Diversity, biosynthesis, functions and uses

Monounsaturated fatty acids are straight-chain aliphatic monocarboxylic acids comprising a unique carbon‑carbon double bond, also termed unsaturation. More than 50 distinct molecular structures have been described in the plant kingdom, and more remain to be discovered. The evolution of land plants has apparently resulted in the convergent evolution of non-homologous enzymes catalyzing the dehydrogenation of saturated acyl chain substrates in a chemo-, regio- and stereoselective manner. Contrasted enzymatic characteristics and different subcellular localizations of these desaturases account for the diversity of existing fatty acid structures. Interestingly, the location and geometrical configuration of the unsaturation confer specific characteristics to these molecules found in a variety of membrane, storage, and surface lipids. An ongoing research effort aimed at exploring the links existing between fatty acid structures and their biological functions has already unraveled the importance of several monounsaturated fatty acids in various physiological and developmental contexts. What is more, the monounsaturated acyl chains found in the oils of seeds and fruits are widely and increasingly used in the food and chemical industries due to the physicochemical properties inherent in their structures. Breeders and plant biotechnologists therefore develop new crops with high monounsaturated contents for various agro-industrial purposes.

Computational study on epoxidation of propylene by dioxygen using the silanol-functionalized polyoxometalate-supported osmium oxide catalyst

The silanol-functionalized POM-supported single-site Os oxide catalyst has been theoretically considered for epoxidation of propylene in the presence of dioxygen based on density functional theory calculations.

Quantitative Assessment of rPM6 for Fluorine- and Chlorine-Containing Metal Complexes: Comparison with Experimental, First-Principles, and Other Semiempirical Results

We report a reparameterization of PM6 parameters for fluorine and chlorine using our training set containing transition metal complexes. Spin unrestricted calculations with the resulting rPM6 (UrPM6) were examined quantitatively using two test sets: (i) the description of magnetic interactions in 25 dinuclear metal complexes and (ii) the prediction of barrier heights and reaction energies for epoxidation and fluorination reactions catalyzed by high-valent manganese-oxo species. The conventional UPM6 and UPM7 methods were also evaluated for comparison on the basis of either experimental or computational (the UB3LYP/SVP level) outcomes. The merits of UrPM6 are highlighted by both the test sets. As regards magnetic exchange coupling constants, the UrPM6 method had the smallest mean absolute errors from the experimental data (19 cm-1), followed by UPM7 (119 cm-1) and UPM6 (373 cm-1). For the epoxidation and fluorination reactions, all of the transition state searches were successful using UrPM6, while the success rates obtained by UPM6 and UPM7 were only 50%. The UrPM6-optimized stationary points also agreed well with the reference UB3LYP-optimized geometries. The accuracy for estimating reaction energies was also greatly remedied.

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.

In-crystal reaction cycle of a toluene-bound diiron hydroxylase

Electrophilic aromatic substitution is one of the most important and recognizable classes of organic chemical transformation. Enzymes create the strong electrophiles that are needed for these highly energetic reactions by using O₂, electrons, and metals or other cofactors. Although the nature of the oxidants that carry out electrophilic aromatic substitution has been deduced from many approaches, it has been difficult to determine their structures. Here we show the structure of a diiron hydroxylase intermediate formed during a reaction with toluene. Density functional theory geometry optimizations of an active site model reveal that the intermediate is an arylperoxo Fe²⁺/Fe³⁺ species with delocalized aryl radical character. The structure suggests that a carboxylate ligand of the diiron centre may trigger homolytic cleavage of the O-O bond by transferring a proton from a metal-bound water. Our work provides the spatial and electronic constraints needed to propose a comprehensive mechanism for diiron enzyme arene hydroxylation that accounts for many prior experimental results.

Geometrical properties of the manganese(iv)/iron(iii) cofactor of Chlamydia trachomatis ribonucleotide reductase unveiled by simulations of XAS spectra

A combination of EXAFS simulations and DFT calculations, including a novel protocol to evaluate Debye–Waller factors, provide insights into the structure of the Mn(iv)/Fe(iii) cofactor ofCtR2.

A radical rebound mechanism for the methane oxidation reaction promoted by the dicopper center of a pMMO enzyme: a computational perspective

Hydrogen Atom Transfer (HAT) promoted by a triplet state of the bis-oxoCu₂(iii) core generates a new radical rebound mechanism for the hydroxylation of methane catalyzed by the binuclear copper site of a pMMO enzyme.

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.

Computational Examination on the Active Site Structure of a (Peroxo)diiron(III) Intermediate in the Amine Oxygenase AurF

In this work, we report the first computational investigation on the structure and properties of the (peroxo)diiron(III) intermediate of the AurF enzyme. Our calculations predict that, in the oxidized state of the AurF enzyme, the peroxo ligand is depicted in a μ-1,1-coordination mode with a protonated bridging ligand and is not in a μ-η(2):η(2) or μ-1,2 mode. Computed spectral data for the μ-1,1-coordination mode correlate well with experimental observations and unravel the potential of the energetics-spectroscopic approach adapted here.

Mechanistic insights into the reductive dehydroxylation pathway for the biosynthesis of isoprenoids promoted by the IspH enzyme

Here, we report an integrated quantum mechanics/molecular mechanics (QM/MM) study of the bio-organometallic reaction pathway of the 2H+/2e- reduction of (E)-4-hydroxy-3-methylbut-2-enyl pyrophosphate (HMBPP) into the so called universal terpenoid precursors isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP), promoted by the IspH enzyme. Our results support the viability of the bio-organometallic pathway through rotation of the OH group of HMBPP away from the [Fe4S4] cluster at the core of the catalytic site, to become engaged in a H-bond with Glu126. This rotation is synchronous with π-coordination of the C2[double bond, length as m-dash]C3 double bond of HMBPP to the apical Fe atom of the [Fe4S4] cluster. Dehydroxylation of HMBPP is triggered by a proton transfer from Glu126 to the OH group of HMBPP. The reaction pathway is completed by competitive proton transfer from the terminal phosphate group to the C2 or C4 atom of HMBPP.

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.

Compound Q is finally deciphered

The elusive compound Q in methane monooxygenase has been characterised and firmly established as a diamond-core bis-μ-oxo Fe^IVFe^IVcluster.

A 2.8 Å Fe-Fe separation in the Fe2(III/IV) intermediate, X, from Escherichia coli ribonucleotide reductase

A class Ia ribonucleotide reductase (RNR) employs a μ-oxo-Fe2(III/III)/tyrosyl radical cofactor in its β subunit to oxidize a cysteine residue ~35 Å away in its α subunit; the resultant cysteine radical initiates substrate reduction. During self-assembly of the Escherichia coli RNR-β cofactor, reaction of the protein's Fe2(II/II) complex with O2 results in accumulation of an Fe2(III/IV) cluster, termed X, which oxidizes the adjacent tyrosine (Y122) to the radical (Y122(•)) as the cluster is converted to the μ-oxo-Fe2(III/III) product. As the first high-valent non-heme-iron enzyme complex to be identified and the key activating intermediate of class Ia RNRs, X has been the focus of intensive efforts to determine its structure. Initial characterization by extended X-ray absorption fine structure (EXAFS) spectroscopy yielded a Fe-Fe separation (d(Fe-Fe)) of 2.5 Å, which was interpreted to imply the presence of three single-atom bridges (O(2-), HO(-), and/or μ-1,1-carboxylates). This short distance has been irreconcilable with computational and synthetic models, which all have d(Fe-Fe) ≥ 2.7 Å. To resolve this conundrum, we revisited the EXAFS characterization of X. Assuming that samples containing increased concentrations of the intermediate would yield EXAFS data of improved quality, we applied our recently developed method of generating O2 in situ from chlorite using the enzyme chlorite dismutase to prepare X at ~2.0 mM, more than 2.5 times the concentration realized in the previous EXAFS study. The measured d(Fe-Fe) = 2.78 Å is fully consistent with computational models containing a (μ-oxo)2-Fe2(III/IV) core. Correction of the d(Fe-Fe) brings the experimental data and computational models into full conformity and informs analysis of the mechanism by which X generates Y122(•).

The sulfur shift: an activation mechanism for periplasmic nitrate reductase and formate dehydrogenase

A structural rearrangement known as sulfur shift occurs in some Mo-containing enzymes of the DMSO reductase family. This mechanism is characterized by the displacement of a coordinating cysteine thiol (or SeCys in Fdh) from the first to the second shell of the Mo-coordination sphere metal. The hexa-coordinated Mo ion found in the as-isolated state cannot bind directly any exogenous ligand (substrate or inhibitors), while the penta-coordinated ion, attained upon sulfur shift, has a free binding site for direct coordination of the substrate. This rearrangement provides an efficient mechanism to keep a constant coordination number throughout an entire catalytic pathway. This mechanism is very similar to the carboxylate shift observed in Zn-dependent enzymes, and it has been recently detected by experimental means. In the present paper, we calculated the geometries and energies involved in the sulfur-shift mechanism using QM-methods (M06/(6-311++G(3df,2pd),SDD)//B3LYP/(6-31G(d),SDD)). The results indicated that the sulfur-shift mechanism provides an efficient way to enable the metal ion for substrate coordination.

Intermediate P* from soluble methane monooxygenase contains a diferrous cluster

During a single turnover of the hydroxylase component (MMOH) of soluble methane monooxygenase from Methylosinus trichosporium OB3b, several discrete intermediates are formed. The diiron cluster of MMOH is first reduced to the Fe(II)Fe(II) state (H(red)). O₂ binds rapidly at a site away from the cluster to form the Fe(II)Fe(II) intermediate O, which converts to an Fe(III)Fe(III)-peroxo intermediate P and finally to the Fe(IV)Fe(IV) intermediate Q. Q binds and reacts with methane to yield methanol and water. The rate constants for these steps are increased by a regulatory protein, MMOB. Previously reported transient kinetic studies have suggested that an intermediate P* forms between O and P in which the g = 16 EPR signal characteristic of the reduced diiron cluster of H(red) and O is lost. This was interpreted as signaling oxidation of the cluster, but a low level of accumulation of P* prevented further characterization. In this study, three methods for directly detecting and trapping P* are applied together to allow its spectroscopic and kinetic characterization. First, the MMOB mutant His33Ala is used to specifically slow the decay of P* without affecting its formation rate, leading to its nearly quantitative accumulation. Second, spectra-kinetic data collection is used to provide a sensitive measure of the formation and decay rate constants of intermediates as well as their optical spectra. Finally, the substrate furan is included to react with Q and quench its strong chromophore. The optical spectrum of P* closely mimics those of H(red) and O, but it is distinctly different from that of P. The reaction cycle rate constants allowed prediction of the times for maximal accumulation of the intermediates. Mössbauer spectra of rapid freeze-quench samples at these times show that the intermediates are formed at almost exactly the predicted levels. The Mössbauer spectra show that the diiron cluster of P*, quite unexpectedly, is in the Fe(II)Fe(II) state. Thus, the loss of the g = 16 EPR signal results from a change in the electronic structure of the Fe(II)Fe(II) center rather than oxidation. The similarity of the optical and Mössbauer spectra of H(red), O, and P* suggests that only subtle changes occur in the electronic and physical structure of the diiron cluster as P* forms. Nevertheless, the changes that do occur are necessary for O₂ to be activated for hydrocarbon oxidation.

Reactions of acids with naphthyridine-functionalized ferrocenes: protonation and metal extrusion

Reaction of 1,8-naphthyrid-2-yl-ferrocene (FcNP) with a variety of acids affords protonated salts at first, whereas longer reaction time leads to partial demetalation of FcNP resulting in a series of Fe complexes. The corresponding salts [FcNP·H][X] (X = BF(4) or CF(3)SO(3) (1)) are isolated for HBF(4) and CF(3)SO(3)H. Reaction of FcNP with equimolar amount of CF(3)CO(2)H for 12 h affords a neutral complex [Fe(FcNP)(2)(O(2)CCF(3))(2)(OH(2))(2)] (2). Use of excess acid gave a trinuclear Fe(II) complex [Fe(3)(H(2)O)(2)(O(2)CCF(3))(8)(FcNP·H)(2)] (3). Three linear iron atoms are held together by four bridging trifluoroacetates and two aqua ligands in a symmetric fashion. Reaction with ethereal solution of HCl afforded [(FcNP·H)(3)(Cl)][FeCl(4)](2) (4) irrespective of the amount of the acid used. Even the picric acid (HPic) led to metal extrusion giving rise to [Fe(2)(Cl)(2)(FcNP)(2)(Pic)(2)] (5) when crystallized from dichloromethane. Metal extrusion was also observed for CF(3)SO(3)H, but an analytically pure compound could not be isolated. The demetalation reaction proceeds with an initial proton attack to the distal nitrogen of the NP unit. Subsequently, coordination of the conjugate base to the electrophilic Fe facilitates the release of Cp rings from metal. The conjugate base plays an important role in the demetalation process and favors the isolation of the Fe complex as well. The 1,1'-bis(1,8-naphthyrid-2-yl)ferrocene (FcNP(2)) does not undergo demetalation under identical conditions. Two NP units share one positive charge causing the Fe-Cp bonds weakened to an extent that is not sufficient for demetalation. X-ray structure of the monoprotonated FcNP(2) reveals a discrete dimer [(FcNP(2)·H)](2)[OTf](2) (6) supported by two N-H···N hydrogen bonds. Crystal packing and dispersive forces associated with intra- and intermolecular π-π stacking interactions (NP···NP and Cp···NP) allow the formation of the dimer in the solid-state. The protonation and demetalation reactions of FcNP and FcNP(2) with a variety of acids are reported.

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.

O(2)-evolving chlorite dismutase as a tool for studying O(2)-utilizing enzymes

The direct interrogation of fleeting intermediates by rapid-mixing kinetic methods has significantly advanced our understanding of enzymes that utilize dioxygen. The gas's modest aqueous solubility (<2 mM at 1 atm) presents a technical challenge to this approach, because it limits the rate of formation and extent of accumulation of intermediates. This challenge can be overcome by use of the heme enzyme chlorite dismutase (Cld) for the rapid, in situ generation of O(2) at concentrations far exceeding 2 mM. This method was used to define the O(2) concentration dependence of the reaction of the class Ic ribonucleotide reductase (RNR) from Chlamydia trachomatis, in which the enzyme's Mn(IV)/Fe(III) cofactor forms from a Mn(II)/Fe(II) complex and O(2) via a Mn(IV)/Fe(IV) intermediate, at effective O(2) concentrations as high as ~10 mM. With a more soluble receptor, myoglobin, an O(2) adduct accumulated to a concentration of >6 mM in <15 ms. Finally, the C-H-bond-cleaving Fe(IV)-oxo complex, J, in taurine:α-ketoglutarate dioxygenase and superoxo-Fe(2)(III/III) complex, G, in myo-inositol oxygenase, and the tyrosyl-radical-generating Fe(2)(III/IV) intermediate, X, in Escherichia coli RNR, were all accumulated to yields more than twice those previously attained. This means of in situ O(2) evolution permits a >5 mM "pulse" of O(2) to be generated in <1 ms at the easily accessible Cld concentration of 50 μM. It should therefore significantly extend the range of kinetic and spectroscopic experiments that can routinely be undertaken in the study of these enzymes and could also facilitate resolution of mechanistic pathways in cases of either sluggish or thermodynamically unfavorable O(2) addition steps.

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.

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.

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.

A transition path ensemble study reveals a linchpin role for Mg(2+) during rate-limiting ADP release from protein kinase A

Protein kinases are key regulators of diverse signaling networks critical for growth and development. Protein kinase A (PKA) is an important kinase prototype that phosphorylates protein targets at Ser and Thr residues by converting ATP to ADP. Mg²⁺ ions play a crucial role in regulating phosphoryl transfer and can limit overall enzyme turnover by affecting ADP release. However, the mechanism by which Mg²⁺ participates in ADP release is poorly understood. Here we use a novel transition path ensemble technique, the harmonic Fourier beads method, to explore the atomic and energetic details of the Mg²⁺-dependent ADP binding and release. Our studies demonstrate that adenine-driven ADP binding to PKA creates three ion-binding sites at the ADP/PKA interface that are absent otherwise. Two of these sites bind the previously characterized Mg²⁺ ions, whereas the third site binds a monovalent cation with high affinity. This third site can bind the P-3 residue of substrate proteins and may serve as a reporter of the active site occupation. Binding of Mg²⁺ ions restricts mobility of the Gly-rich loop that closes over the active site. We find that simultaneous release of ADP with Mg²⁺ ions from the active site is unfeasible. Thus, we conclude that Mg²⁺ ions act as a linchpin and that at least one ion must be removed prior to pyrophosphate-driven ADP release. The results of the present study enhance understanding of Mg²⁺-dependent association of nucleotides with protein kinases.