A DFT study of the cis-dihydroxylation of nitroaromatic compounds catalyzed by nitrobenzene dioxygenase

Custom tuning of Rieske oxygenase reactivity

Rieske oxygenases use a Rieske-type [2Fe-2S] cluster and a mononuclear iron center to initiate a range of chemical transformations. However, few details exist regarding how this catalytic scaffold can be predictively tuned to catalyze divergent reactions. Therefore, in this work, using a combination of structural analyses, as well as substrate and rational protein-based engineering campaigns, we elucidate the architectural trends that govern catalytic outcome in the Rieske monooxygenase TsaM. We identify structural features that permit a substrate to be functionalized by TsaM and pinpoint active-site residues that can be targeted to manipulate reactivity. Exploiting these findings allowed for custom tuning of TsaM reactivity: substrates are identified that support divergent TsaM-catalyzed reactions and variants are created that exclusively catalyze dioxygenation or sequential monooxygenation chemistry. Importantly, we further leverage these trends to tune the reactivity of additional monooxygenase and dioxygenase enzymes, and thereby provide strategies to custom tune Rieske oxygenase reaction outcomes.

Rieske Oxygenases and Other Ferredoxin-Dependent Enzymes: Electron Transfer Principles and Catalytic Capabilities

Enzymes that depend on sophisticated electron transfer via ferredoxins (Fds) exhibit outstanding catalytic capabilities, but despite decades of research, many of them are still not well understood or exploited for synthetic applications. This review aims to provide a general overview of the most important Fd-dependent enzymes and the electron transfer processes involved. While several examples are discussed, we focus in particular on the family of Rieske non-heme iron-dependent oxygenases (ROs). In addition to illustrating their electron transfer principles and catalytic potential, the current state of knowledge on structure-function relationships and the mode of interaction between the redox partner proteins is reviewed. Moreover, we highlight several key catalyzed transformations, but also take a deeper dive into their engineerability for biocatalytic applications. The overall findings from these case studies highlight the catalytic capabilities of these biocatalysts and could stimulate future interest in developing additional Fd-dependent enzyme classes for synthetic applications.

Theoretical Kinetic Isotope Effects in Establishing the Precise Biodegradation Mechanisms of Organic Pollutants

Compound-specific isotope analysis (CSIA) for natural isotope ratios has been recognized as a promising tool to elucidate biodegradation pathways of organic pollutants by microbial enzymes by relating reported kinetic isotope effects (KIEs) to apparent KIEs (AKIEs) derived from bulk isotope fractionations (ε_bulk). However, for many environmental reactions, neither are the reference KIE ranges sufficiently narrow nor are the mechanisms elucidated to the point that rate-determining steps have been identified unequivocally. In this work, besides providing reference KIEs and rationalizing AKIEs, good relationships have been explained by DFT computations for diverse biodegradation pathways with known enzymatic models between the theoretical isotope fractionations (ε_bulk^') from intrinsic KIEs on the rate-determining steps and the observed ε_bulk. (1) To confirm the mechanistic details of previously reported pathway-dependent CSIA, it includes isotope changes in MTBE biodegradation between hydroxylation by CYP450 and S_N2 reaction by cobalamin-dependent methyltransferase, the regioselectivity of toluene biodegradation by CYP450, and the rate-determining step in toluene biodegradation by benzylsuccinate synthase. (2) To yield new fundamental insights into some unclear biodegradation pathways, it consists of the oxidative function of toluene dioxygenase in biodegradation of TCE, the epoxidation mode in biodegradation of TCE by toluene 4-monooxygenase, and the weighted average mechanism in biodegradation of cDCE by CYP450.

Comparative oxidative ability of mononuclear and dinuclear high-valent iron-oxo species towards the activation of methane: does the axial/bridge atom modulate the reactivity?

Over the years, mononuclear Fe^IVO species have been extensively studied, but the presence of dinuclear Fe^IVO species in soluble methane monooxygenase (sMMO) has inspired the development of biomimic models that could activate inert substrates such as methane. There are some successful attempts; particularly the [(Por)(m-CBA) Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- species has been reported to activate methane and yield decent catalytic turnover numbers and therefore regarded as the closest to the sMMO enzyme functional model, as no mononuclear Fe^IVO analogues could achieve this feat. In this work, we have studied a series of mono and dinuclear models using DFT and ab initio DLPNO-CCSD(T) calculations to probe the importance of nuclearity in enhancing the reactivity. We have probed the catalytic activities of four complexes: [(HO)Fe^IV(O)(Por)]^- (1), [(HO)Fe^IV(O)(Por˙⁺)] (2), μ-oxo dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-O)Fe^IV(O) (Por˙⁺)]^- (3) and N-bridged dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- (4) towards the activation of methane. Additionally, calculations were performed on the mononuclear models [(X)Fe^IV(O)(Por˙⁺)]ⁿ {X = N 4a (n = -2), NH 4b (n = -1) and NH₂4c (n = 0)} to understand the role of nuclearity in the reactivity. DFT calculations performed on species 1-4 suggest an interesting variation among them, with species 1-3 possessing an intermediate spin (S = 1) as a ground state and species 4 possessing a high-spin (S = 2) as a ground state. Furthermore, the two Fe^IV centres in species 3 and 4 are antiferromagnetically coupled, yielding a singlet state with a distinct difference in their electronic structure. On the other hand, species 2 exhibits a ferromagnetic coupling between the Fe^IV and the Por˙⁺ moiety. Our calculations suggest that the higher barriers for the C-H bond activation of methane and the rebound step for species 1 and 3 are very high in energy, rendering them unreactive towards methane, while species 2 and 4 have lower barriers, suggesting their reactivity towards methane. Studies on the system reveal that model 4a has multiple FeN bonds facilitating greater reactivity, whereas the other two models have longer Fe-N bonds and less radical character with steeper barriers. Strong electronic cooperativity is found to be facilitated by the bridging nitride atom, and this cooperativity is suppressed by substituents such as oxygen, rendering them inactive. Thus, our study unravels that apart from enhancing the nuclearity, bridging atoms that facilitate strong cooperation between the metals are required to activate very inert substrates such as methane, and our results are broadly in agreement with earlier experimental findings.

Elucidating the Role of O2 Uncoupling in the Oxidative Biodegradation of Organic Contaminants by Rieske Non-heme Iron Dioxygenases

Oxygenations of aromatic soil and water contaminants with molecular O₂ catalyzed by Rieske dioxygenases are frequent initial steps of biodegradation in natural and engineered environments. Many of these non-heme ferrous iron enzymes are known to be involved in contaminant metabolism, but the understanding of enzyme–substrate interactions that lead to successful biodegradation is still elusive. Here, we studied the mechanisms of O₂ activation and substrate hydroxylation of two nitroarene dioxygenases to evaluate enzyme- and substrate-specific factors that determine the efficiency of oxygenated product formation. Experiments in enzyme assays of 2-nitrotoluene dioxygenase (2NTDO) and nitrobenzene dioxygenase (NBDO) with methyl-, fluoro-, chloro-, and hydroxy-substituted nitroaromatic substrates reveal that typically 20–100% of the enzyme’s activity involves unproductive paths of O₂ activation with generation of reactive oxygen species through so-called O₂ uncoupling. The ¹⁸O and ¹³C kinetic isotope effects of O₂ activation and nitroaromatic substrate hydroxylation, respectively, suggest that O₂ uncoupling occurs after generation of Fe^III-(hydro)peroxo species in the catalytic cycle. While 2NTDO hydroxylates ortho-substituted nitroaromatic substrates more efficiently, NBDO favors meta-substituted, presumably due to distinct active site residues of the two enzymes. Our data implies, however, that the O₂ uncoupling and hydroxylation activity cannot be assessed from simple structure–reactivity relationships. By quantifying O₂ uncoupling by Rieske dioxygenases, our work provides a mechanistic link between contaminant biodegradation, the generation of reactive oxygen species, and possible adaptation strategies of microorganisms to the exposure of new contaminants.

Substrate-Specific Coupling of O2 Activation to Hydroxylations of Aromatic Compounds by Rieske Non-heme Iron Dioxygenases

Rieske dioxygenases catalyze the initial steps in the hydroxylation of aromatic compounds and are critical for the metabolism of xenobiotic substances. Because substrates do not bind to the mononuclear non-heme Fe^II center, elementary steps leading to O₂ activation and substrate hydroxylation are difficult to delineate, thus making it challenging to rationalize divergent observations on enzyme mechanisms, reactivity, and substrate specificity. Here, we show for nitrobenzene dioxygenase, a Rieske dioxygenase capable of transforming nitroarenes to nitrite and substituted catechols, that unproductive O₂ activation with the release of the unreacted substrate and reactive oxygen species represents an important path in the catalytic cycle. Through correlation of O₂ uncoupling for a series of substituted nitroaromatic compounds with ¹⁸O and ¹³C kinetic isotope effects of dissolved O₂ and aromatic substrates, respectively, we show that O₂ uncoupling occurs after the rate-limiting formation of Fe^III-(hydro)peroxo species from which substrates are hydroxylated. Substituent effects on the extent of O₂ uncoupling suggest that the positioning of the substrate in the active site rather than the susceptibility of the substrate for attack by electrophilic oxygen species is responsible for unproductive O₂ uncoupling. The proposed catalytic cycle provides a mechanistic basis for assessing the very different efficiencies of substrate hydroxylation vs unproductive O₂ activation and generation of reactive oxygen species in reactions catalyzed by Rieske dioxygenases.

Investigating the molecular mechanism of hydroxylated bromdiphenyl ethers to inhibit the thyroid hormone sulfotransferase SULT1A1

Hydroxylated bromodiphenyl ethers (OH-BDEs) have raised great concern due to their potential endocrine disrupting effects on humans. In vitro experiments have indicated OH-BDEs can inhibit the activity of thyroid hormone (TH) sulfotransferases (SULTs); however, the molecular mechanism has not been investigated in depth. In this work, we employed 17 OH-BDEs with five or fewer Br atoms, and performed integrated computational simulations to unravel the possible inhibition mechanism of OH-BDEs on human SULT1A1. The molecular docking results demonstrate that OH-BDEs form hydrogen bonds with residues Lys106 and His108, and the neutral OH-BDEs show comparable binding energies with their anionic counterparts. The further hybrid quantum mechanical/molecular mechanical (QM/MM) calculations unravel a metabolic mechanism of OH-BDEs comprised by proton abstraction and sulfation steps. This mechanism is involved in the SULT1A1 inhibition by some OH-BDEs comprised of three or fewer Br atoms, while other OH-BDEs likely only form ternary complexes to competitively inhibit SULT1A1 activity. Moreover, the effect of the hydroxyl group of OH-BDEs on SULT1A1 inhibition potential follows the order of ortho-OH BDE > meta-OH BDE > para-OH BDE. These results provide an insight into the inhibition mechanism of OH-BDEs to SULT1A1 at the molecular level, which are beneficial in illuminating the molecular initiating events involved in the TH disruption of OH-BDEs.

Molecular Basis for Metabolic Regioselectivity and Mechanism of Cytochrome P450s toward Carcinogenic 4-(Methylnitrosamino)-(3-pyridyl)-1-butanone

As an abundantly present tobacco component, carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) has also been detected in atmospheric particulate matter, suggesting the ineluctable exposure risk of this contaminant. NNK metabolic activation by cytochrome P450 enzymes (CYPs) is a prerequisite to exerting its genotoxicity, but the metabolic regioselectivity and mechanism are still unknown. Here the binding feature and regioselectivity of CYPs 1A1, 1A2, 2A6, 2A13, 2B6, and 3A4 toward NNK are unraveled through molecular docking and molecular dynamics (MD) simulations. Binding mode analyses reveal that 1A2 and 2B6 have definite preferences for NNK α-methyl hydroxylation, while the other four CYPs preferentially catalyze α-methylene hydroxylation. The binding affinities between NNK and CYPs evaluated by the binding free energies follow the order 2A13 > 2B6 > 1A2 > 2A6 > 1A1 > 3A4. Density functional theory (DFT) calculations are further performed to characterize the mechanism of NNK biotransformation. Results show that the α-hydroxyNNK generated from α-hydroxylation may undergo nonenzymatic decomposition to form genotoxic diazohydroxide and aldehyde, and further oxidation by P450 to yield nitrosamide, which mainly contributes to NNK toxification capacity. Meanwhile the pyridine N-oxidation and denitrosation of C_α-radical intermediate play an important role in detoxifying NNK. Overall, the present study provides the molecular basis for CYP-catalyzed regioselectivity and mechanism of NNK biotransformation, which can enable the identification of metabolites for assessing the health risk of individual NNK exposure.

Bacterial steroid hydroxylases: enzyme classes, their functions and comparison of their catalytic mechanisms

The steroid superfamily includes a wide range of compounds that are essential for living organisms of the animal and plant kingdoms. Structural modifications of steroids highly affect their biological activity. In this review, we focus on hydroxylation of steroids by bacterial hydroxylases, which take part in steroid catabolic pathways and play an important role in steroid degradation. We compare three distinct classes of metalloenzymes responsible for aerobic or anaerobic hydroxylation of steroids, namely: cytochrome P450, Rieske-type monooxygenase 3-ketosteroid 9α-hydroxylase, and molybdenum-containing steroid C25 dehydrogenases. We analyze the available literature data on reactivity, regioselectivity, and potential application of these enzymes in organic synthesis of hydroxysteroids. Moreover, we describe mechanistic hypotheses proposed for all three classes of enzymes along with experimental and theoretical evidences, which have provided grounds for their formulation. In case of the 3-ketosteroid 9α-hydroxylase, such a mechanistic hypothesis is formulated for the first time in the literature based on studies conducted for other Rieske monooxygenases. Finally, we provide comparative analysis of similarities and differences in the reaction mechanisms utilized by bacterial steroid hydroxylases.

NRVS Studies of the Peroxide Shunt Intermediate in a Rieske Dioxygenase and Its Relation to the Native FeII O2 Reaction

The Rieske dioxygenases are a major subclass of mononuclear nonheme iron enzymes that play an important role in bioremediation. Recently, a high-spin FeIII-(hydro)peroxy intermediate (BZDOp) has been trapped in the peroxide shunt reaction of benzoate 1,2-dioxygenase. Defining the structure of this intermediate is essential to understanding the reactivity of these enzymes. Nuclear resonance vibrational spectroscopy (NRVS) is a recently developed synchrotron technique that is ideal for obtaining vibrational, and thus structural, information on Fe sites, as it gives complete information on all vibrational normal modes containing Fe displacement. In this study, we present NRVS data on BZDOp and assign its structure using these data coupled to experimentally calibrated density functional theory calculations. From this NRVS structure, we define the mechanism for the peroxide shunt reaction. The relevance of the peroxide shunt to the native FeII/O2 reaction is evaluated. For the native FeII/O2 reaction, an FeIII-superoxo intermediate is found to react directly with substrate. This process, while uphill thermodynamically, is found to be driven by the highly favorable thermodynamics of proton-coupled electron transfer with an electron provided by the Rieske [2Fe-2S] center at a later step in the reaction. These results offer important insight into the relative reactivities of FeIII-superoxo and FeIII-hydroperoxo species in nonheme Fe biochemistry.

Resolving Discrepancy between Theory and Experiment in 4-Nitrotoluene Oxidation

We have performed calculations of possible oxidation pathways of 4-nitrotoluene (4NT) by permanganate anion and evaluated relative contributions of oxidation of the methyl group and aromatic ring. Although a few theory levels matched the experimental results obtained by compound specific isotope analysis (CSIA) for 4NT, they failed in reproducing results for other nitrotoluene derivatives studied previously [Wijker, R.S.; Adamczyk, P.; Bolotin, J.; Paneth, P.; Hofstetter, T.B. Environ. Sci. Technol., 2013, 47, 13459-13468]. This discrepancy prompted us to reevaluate the experimental isotopic fractionation of carbon and hydrogen for 4NT on which the relative contributions of the oxidation channels has been based. Using position specific isotope analysis (PSIA) for hydrogen isotopic fractionation we have found that the previously determined value was incorrect. Reexamination of theory levels that are in agreement with these new findings indicated that while a better agreement for this particular case can be reached, overall, the previously used B3LYP functional expressed in the 6-31+G(d,p) basis set with inclusion of the polarized continuum model (PCM) of aqueous solution remains the theoretical level of choice in modeling oxidation of nitrotoluene derivatives by permanganate.

Theoretical Study of the Mechanism of the Nonheme Iron Enzyme EgtB

EgtB is a nonheme iron enzyme catalyzing the C-S bond formation between γ-glutamyl cysteine (γGC) and N-α-trimethyl histidine (TMH) in the ergothioneine biosynthesis. Density functional calculations were performed to elucidate and delineate the reaction mechanism of this enzyme. Two different mechanisms were considered, depending on whether the sulfoxidation or the S-C bond formation takes place first. The calculations suggest that the S-O bond formation occurs first between the thiolate and the ferric superoxide, followed by homolytic O-O bond cleavage, very similar to the case of cysteine dioxygenase. Subsequently, proton transfer from a second-shell residue Tyr377 to the newly generated iron-oxo moiety takes place, which is followed by proton transfer from the TMH imidazole to Tyr377, facilitated by two crystallographically observed water molecules. Next, the S-C bond is formed between γGC and TMH, followed by proton transfer from the imidazole CH moiety to Tyr377, which was calculated to be the rate-limiting step for the whole reaction, with a barrier of 17.9 kcal/mol in the quintet state. The calculated barrier for the rate-limiting step agrees quite well with experimental kinetic data. Finally, this proton is transferred back to the imidazole nitrogen to form the product. The alternative thiyl radical attack mechanism has a very high barrier, being 25.8 kcal/mol, ruling out this possibility.

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.

Substrate and Enzyme Specificity of the Kinetic Isotope Effects Associated with the Dioxygenation of Nitroaromatic Contaminants

Compound-specific isotope analysis (CSIA) is a promising approach for tracking biotransformation of organic pollutants, but isotope fractionation associated with aromatic oxygenations is only poorly understood. We investigated the dioxygenation of a series of nitroaromatic compounds to the corresponding catechols by two enzymes, namely, nitrobenzene and 2-nitrotoluene dioxygenase (NBDO and 2NTDO) to elucidate the enzyme- and substrate-specificity of C and H isotope fractionation. While the apparent (13)C- and (2)H-kinetic isotope effects of nitrobenzene, nitrotoluene isomers, 2,6-dinitrotoluene, and naphthalene dioxygenation by NBDO varied considerably, the correlation of C and H isotope fractionation revealed a common mechanism for nitrobenzene and nitrotoluenes. Similar observations were made for the dioxygenation of these substrates by 2NTDO. Evaluation of reaction kinetics, isotope effects, and commitment-to-catalysis based on experiment and theory showed that rates of dioxygenation are determined by the enzymatic O2 activation and aromatic C oxygenation. The contribution of enzymatic O2 activation to the reaction rate varies for different nitroaromatic substrates of NBDO and 2NTDO. Because aromatic dioxygenation by nonheme iron dioxygenases is frequently the initial step of biodegradation, O2 activation kinetics may also have been responsible for the minor isotope fractionation reported for the oxygenation of other aromatic contaminants.

Insights into the Reaction Mechanism of Aromatic Ring Cleavage by Homogentisate Dioxygenase: A Quantum Mechanical/Molecular Mechanical Study

To elucidate the reaction mechanism of the ring cleavage of homogentisate by homogentisate dioxygenase, quantum mechanical/molecular mechanical (QM/MM) calculations were carried out by using two systems in different protonation states of the substrate C2 hydroxyl group. When the substrate C2 hydroxyl group is ionized (the ionized pathway), the superoxo attack on the substrate is the rate-limiting step in the catalytic cycle, with a barrier of 15.9 kcal/mol. Glu396 was found to play an important role in stabilizing the bridge species and its O-O cleavage product by donating a proton via a hydrogen-bonded water molecule. When the substrate C2 hydroxyl group is not ionized (the nonionized pathway), the O-O bond cleavage of the bridge species is the rate-limiting step, with a barrier of 15.3 kcal/mol. The QM/MM-optimized geometries for the dioxygen and alkylperoxo complexes using the nonionized model (for the C2 hydroxyl group) are in agreement with the experimental crystal structures, suggesting that the C2 hydroxyl group is more likely to be nonionized.

Measurement of oxygen isotope ratios ((18)O/(16)O) of aqueous O2 in small samples by gas chromatography/isotope ratio mass spectrometry

RATIONALE

Oxygen isotope fractionation of molecular O2 is an important process for the study of aerobic metabolism, photosynthesis, and formation of reactive oxygen species. The latter is of particular interest for investigating the mechanism of enzyme-catalyzed reactions, such as the oxygenation of organic pollutants, which is an important detoxification mechanism.

METHODS

We developed a simple method to measure the δ(18) O values of dissolved O2 in small samples using automated split injection for gas chromatography coupled to isotope ratio mass spectrometry (GC/IRMS). After creating a N2 headspace, the dissolved O2 partitions from aqueous solution to the headspace, from which it can be injected into the gas chromatograph.

RESULTS

In aqueous samples of 10 mL and in diluted air samples, we quantified the δ(18) O values at O2 concentrations of 16 μM and 86 μM, respectively. The chromatographic separation of O2 and N2 with a molecular sieve column made it possible to use N2 as the headspace gas for the extraction of dissolved O2 from water. We were therefore able to apply a rigorous δ(18) O blank correction for the quantification of (18) O/(16) O ratios in 20 nmol of injected O2 .

CONCLUSIONS

The successful quantification of (18) O-kinetic isotope effects associated with enzymatic and chemical reduction of dissolved O2 illustrates how the proposed method can be applied for studying enzymatic O2 activation mechanisms in a variety of (bio)chemical processes.

Rate-Determining Attack on Substrate Precedes Rieske Cluster Oxidation during Cis-Dihydroxylation by Benzoate Dioxygenase

Rieske dearomatizing dioxygenases utilize a Rieske iron-sulfur cluster and a mononuclear Fe(II) located 15 Å across a subunit boundary to catalyze O2-dependent formation of cis-dihydrodiol products from aromatic substrates. During catalysis, O2 binds to the Fe(II) while the substrate binds nearby. Single-turnover reactions have shown that one electron from each metal center is required for catalysis. This finding suggested that the reactive intermediate is Fe(III)-(H)peroxo or HO-Fe(V)═O formed by O-O bond scission. Surprisingly, several kinetic phases were observed during the single-turnover Rieske cluster oxidation. Here, the Rieske cluster oxidation and product formation steps of a single turnover of benzoate 1,2-dioxygenase are investigated using benzoate and three fluorinated analogues. It is shown that the rate constant for product formation correlates with the reciprocal relaxation time of only the fastest kinetic phase (RRT-1) for each substrate, suggesting that the slower phases are not mechanistically relevant. RRT-1 is strongly dependent on substrate type, suggesting a role for substrate in electron transfer from the Rieske cluster to the mononuclear iron site. This insight, together with the substrate and O2 concentration dependencies of RRT-1, indicates that a reactive species is formed after substrate and O2 binding but before electron transfer from the Rieske cluster. Computational studies show that RRT-1 is correlated with the electron density at the substrate carbon closest to the Fe(II), consistent with initial electrophilic attack by an Fe(III)-superoxo intermediate. The resulting Fe(III)-peroxo-aryl radical species would then readily accept an electron from the Rieske cluster to complete the cis-dihydroxylation reaction.

Mechanism and selectivity of the dinuclear iron benzoyl-coenzyme A epoxidase BoxB

DFT calculations are used to elucidate the reaction mechanism and selectivity of BoxB catalyzed benzoyl-CoA epoxidation.

Benzoyl-CoA epoxidase is a dinuclear iron enzyme that catalyzes the epoxidation reaction of the aromatic ring of benzoyl-CoA with chemo-, regio- and stereo-selectivity. It has been suggested that this enzyme may also catalyze the deoxygenation reaction of epoxide, suggesting a unique bifunctionality among the diiron enzymes. We report a density functional theory study of this enzyme aimed at elucidating its mechanism and the various selectivities. The epoxidation is suggested to start with the binding of the O₂ molecule to the diferrous center to generate a diferric peroxide complex, followed by concerted O–O bond cleavage and epoxide formation. Two different pathways have been located, leading to (2S,3R)-epoxy and (2R,3S)-epoxy products, with barriers of 17.6 and 20.4 kcal mol^–1, respectively. The barrier difference is 2.8 kcal mol^–1, corresponding to a diastereomeric excess of about 99 : 1. Further isomerization from epoxide to phenol is found to have quite a high barrier, which cannot compete with the product release step. After product release into solution, fast epoxide–oxepin isomerization and racemization can take place easily, leading to a racemic mixture of (2S,3R) and (2R,3S) products. The deoxygenation of epoxide to regenerate benzoyl-CoA by a diferrous form of the enzyme proceeds via a stepwise mechanism. The C2–O bond cleavage happens first, coupled with one electron transfer from one iron center to the substrate, to form a radical intermediate, which is followed by the second C3–O bond cleavage. The first step is rate-limiting with a barrier of only 10.8 kcal mol^–1. Further experimental studies are encouraged to verify our results.

Oxidation of methane by an N-bridged high-valent diiron–oxo species: electronic structure implications on the reactivity

Methane activation by dinuclear high-valent iron–oxo species: do we need two metals to activate such inert bonds? Our theoretical study using DFT methods where electronic structure details and mechanistic aspects are established answers this intriguing question.

Isotope effects of enzymatic dioxygenation of nitrobenzene and 2-nitrotoluene by nitrobenzene dioxygenase

Oxygenation of aromatic rings is a frequent initial step in the biodegradation of persistent contaminants, and the accompanying isotope fractionation is increasingly used to assess the extent of transformation in the environment. Here, we systematically investigated the dioxygenation of two nitroaromatic compounds (nitrobenzene and 2-nitrotoluene) by nitrobenzene dioxygenase (NBDO) to obtain insights into the factors governing its C, H, and N isotope fractionation. Experiments were carried out at different levels of biological complexity from whole bacterial cells to pure enzyme. C, H, and N isotope enrichment factors and kinetic isotope effects (KIEs) were derived from the compound-specific isotope analysis of nitroarenes, whereas C isotope fractionation was also quantified in the oxygenated reaction products. Dioxygenation of nitrobenzene to catechol and 2-nitrotoluene to 3-methylcatechol showed large C isotope enrichment factors, ϵC, of -4.1 ± 0.2‰ and -2.5 ± 0.2‰, respectively, and was observed consistently in the substrates and dioxygenation products. ϵH- and ϵN-values were smaller, that is -5.7 ± 1.3‰ and -1.0 ± 0.3‰, respectively. C isotope fractionation was also identical in experiments with whole bacterial cells and pure enzymes. The corresponding (13)C-KIEs for the dioxygenation of nitrobenzene and 2-nitrotoluene were 1.025 ± 0.001 and 1.018 ± 0.001 and suggest a moderate substrate specificity. Our study illustrates that dioxygenation of nitroaromatic contaminants exhibits a large C isotope fractionation, which is not masked by substrate transport and uptake processes and larger than dioxygenation of other aromatic hydrocarbons.

A DFT and ONIOM study of C–H hydroxylation catalyzed by nitrobenzene 1,2-dioxygenase

The DFT and ONIOM calculations show that C–H hydroxylation by nitrobenzene 1,2-dioxygenase proceeds through a HO–Fe^VO intermediate.