Role for Threonine 201 in the Catalytic Cycle of the Soluble Diiron Hydroxylase Toluene 4-Monooxygenase,

HMDB: A curated database of genes involved in hydrocarbon monooxygenation reaction with homologous genes as background

The investigation of hydrocarbon degradation potential of environmental microorganisms is an important research topic, whether for the global carbon cycle or oil pollution remediation. Under aerobic conditions, the microorganisms employ a range of monooxygenases to use hydrocarbons substrates as a source of carbon and energy. With the explosion of sequencing data, mining genes in genomes or metagenomes has become computationally expensive and time-consuming. We proposed the HMDB, which is a professional gene database of hydrocarbon monooxygenases. HMDB contains 38 genes, which encode 11 monooxygenases responsible for the hydroxylation of 8 hydrocarbons. To reduce false positives, the strategy of using homologous genes as background noise was applied for HMDB. We added 10,095 gene sequences of homologous enzymes which took non-hydrocarbons as substrates to HMDB. The classic BLAST method and best-hit strategy were recommended for HMDB usage, but not limited. The performance of HMDB was validated using 264,402 prokaryote genomes from RefSeq and 51 metagenomes from SRA. The results showed that HMDB database had high sensitivity and low false positive rate. We release the HMDB database here, hoping to speed up the process for investigation of hydrocarbon monooxygenases in massive metagenomic data. HMDB is freely available at http://www.orgene.net/HMDB/.

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

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

In Vitro Reconstitution Reveals a Central Role for the N-Oxygenase PvfB in (Dihydro)pyrazine-N-oxide and Valdiazen Biosynthesis

The Pseudomonas virulence factor (pvf) operon is essential for the biosynthesis of two very different natural product scaffolds: the (dihydro)pyrazine-N-oxides and the diazeniumdiolate, valdiazen. PvfB is a member of the non-heme diiron N-oxygenase enzyme family that commonly convert anilines to their nitroaromatic counterparts. In contrast, we show that PvfB catalyzes N-oxygenation of the α-amine of valine, first to the hydroxylamine and then the nitroso, while linked to the carrier protein of PvfC. PvfB modification of PvfC-tethered valine was observed directly by protein NMR spectroscopy, establishing the intermediacy of the hydroxylamine. This work reveals a central role for PvfB in the biosynthesis of (dihydro)pyrazine-N-oxides and valdiazen.

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.

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.

Coupling Oxygen Consumption with Hydrocarbon Oxidation in Bacterial Multicomponent Monooxygenases

A fundamental goal in catalysis is the coupling of multiple reactions to yield a desired product. Enzymes have evolved elegant approaches to address this grand challenge. A salient example is the biological conversion of methane to methanol catalyzed by soluble methane monooxygenase (sMMO), a member of the bacterial multicomponent monooxygenase (BMM) superfamily. sMMO is a dynamic protein complex of three components: a hydroxylase, a reductase, and a regulatory protein. The active site, a carboxylate-rich non-heme diiron center, is buried inside the 251 kDa hydroxylase component. The enzyme processes four substrates: O2, protons, electrons, and methane. To couple O2 activation to methane oxidation, timely control of substrate access to the active site is critical. Recent studies of sMMO, as well as its homologues in the BMM superfamily, have begun to unravel the mechanism. The emerging and unifying picture reveals that each substrate gains access to the active site along a specific pathway through the hydroxylase. Electrons and protons are delivered via a three-amino-acid pore located adjacent to the diiron center; O2 migrates via a series of hydrophobic cavities; and hydrocarbon substrates reach the active site through a channel or linked set of cavities. The gating of these pathways mediates entry of each substrate to the diiron active site in a timed sequence and is coordinated by dynamic interactions with the other component proteins. The result is coupling of dioxygen consumption with hydrocarbon oxidation, avoiding unproductive oxidation of the reductant rather than the desired hydrocarbon. To initiate catalysis, the reductase delivers two electrons to the diiron(III) center by binding over the pore of the hydroxylase. The regulatory component then displaces the reductase, docking onto the same surface of the hydroxylase. Formation of the hydroxylase-regulatory component complex (i) induces conformational changes of pore residues that may bring protons to the active site; (ii) connects hydrophobic cavities in the hydroxylase leading from the exterior to the diiron active site, providing a pathway for O2 and methane, in the case of sMMO, to the reduced diiron center for O2 activation and substrate hydroxylation; (iii) closes the pore, as well as a channel in the case of four-component BMM enzymes, restricting proton access to the diiron center during formation of "Fe2O2" intermediates required for hydrocarbon oxidation; and (iv) inhibits undesired electron transfer to the Fe2O2 intermediates by blocking reductase binding during O2 activation. This mechanism is quite different from that adopted by cytochromes P450, a large class of heme-containing monooxygenases that catalyze reactions very similar to those catalyzed by the BMM enzymes. Understanding the timed enzyme control of substrate access has implications for designing artificial catalysts. To achieve multiple turnovers and tight coupling, synthetic models must also control substrate access, a major challenge considering that nature requires large, multimeric, dynamic protein complexes to accomplish this feat.

Structural basis for biomolecular recognition in overlapping binding sites in a diiron enzyme system

Productive biomolecular recognition requires exquisite control of affinity and specificity. Accordingly, nature has devised many strategies to achieve proper binding interactions. Bacterial multicomponent monooxygenases provide a fascinating example, where a diiron hydroxylase must reversibly interact with both ferredoxin and catalytic effector in order to achieve electron transfer and O2 activation during catalysis. Because these two accessory proteins have distinct structures, and because the hydroxylase-effector complex covers the entire surface closest to the hydroxylase diiron centre, how ferredoxin binds to the hydroxylase has been unclear. Here we present high-resolution structures of toluene 4-monooxygenase hydroxylase complexed with its electron transfer ferredoxin and compare them with the hydroxylase-effector structure. These structures reveal that ferredoxin or effector protein binding produce different arrangements of conserved residues and customized interfaces on the hydroxylase in order to achieve different aspects of catalysis.

A flexible glutamine regulates the catalytic activity of toluene o-xylene monooxygenase

Toluene/o-xylene monooxygenase (ToMO) is a bacterial multicomponent monooxygenase capable of oxidizing aromatic substrates. The carboxylate-rich diiron active site is located in the hydroxylase component of ToMO (ToMOH), buried 12 Å from the surface of the protein. A small, hydrophilic pore is the shortest pathway between the diiron active site and the protein exterior. In this study of ToMOH from Pseudomonas sp. OX1, the functions of two residues lining this pore, N202 and Q228, were investigated using site-directed mutagenesis. Steady-state characterization of WT and the three mutant enzymes demonstrates that residues N202 and Q228 are critical for turnover. Kinetic isotope effects and pH profiles reveal that these residues govern the kinetics of water egress and prevent quenching of activated oxygen intermediates formed at the diiron active site. We propose that this activity arises from movement of these residues, opening and closing the pore during catalysis, as seen in previous X-ray crystallographic studies. In addition, N202 and Q228 are important for the interactions of the reductase and regulatory components to ToMOH, suggesting that they bind competitively to the hydroxylase. The role of the pore in the hydroxylase components of other bacterial multicomponent monooxygenases within the superfamily is discussed in light of these conclusions.

Crystallographic analysis of active site contributions to regiospecificity in the diiron enzyme toluene 4-monooxygenase

Crystal structures of toluene 4-monooxygenase hydroxylase in complex with reaction products and effector protein reveal active site interactions leading to regiospecificity. Complexes with phenolic products yield an asymmetric μ-phenoxo-bridged diiron center and a shift of diiron ligand E231 into a hydrogen bonding position with conserved T201. In contrast, complexes with inhibitors p-NH(2)-benzoate and p-Br-benzoate showed a μ-1,1 coordination of carboxylate oxygen between the iron atoms and only a partial shift in the position of E231. Among active site residues, F176 trapped the aromatic ring of products against a surface of the active site cavity formed by G103, E104 and A107, while F196 positioned the aromatic ring against this surface via a π-stacking interaction. The proximity of G103 and F176 to the para substituent of the substrate aromatic ring and the structure of G103L T4moHD suggest how changes in regiospecificity arise from mutations at G103. Although effector protein binding produced significant shifts in the positions of residues along the outer portion of the active site (T201, N202, and Q228) and in some iron ligands (E231 and E197), surprisingly minor shifts (<1 Å) were produced in F176, F196, and other interior residues of the active site. Likewise, products bound to the diiron center in either the presence or absence of effector protein did not significantly shift the position of the interior residues, suggesting that positioning of the cognate substrates will not be strongly influenced by effector protein binding. Thus, changes in product distributions in the absence of the effector protein are proposed to arise from differences in rates of chemical steps of the reaction relative to motion of substrates within the active site channel of the uncomplexed, less efficient enzyme, while structural changes in diiron ligand geometry associated with cycling between diferrous and diferric states are discussed for their potential contribution to product release.

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

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

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.

Multiple roles of component proteins in bacterial multicomponent monooxygenases: phenol hydroxylase and toluene/o-xylene monooxygenase from Pseudomonas sp. OX1

Phenol hydroxylase (PH) and toluene/o-xylene monooxygenase (ToMO) from Pseudomonas sp. OX1 require three or four protein components to activate dioxygen for the oxidation of aromatic substrates at a carboxylate-bridged diiron center. In this study, we investigated the influence of the hydroxylases, regulatory proteins, and electron-transfer components of these systems on substrate (phenol; NADH) consumption and product (catechol; H(2)O(2)) generation. Single-turnover experiments revealed that only complete systems containing all three or four protein components are capable of oxidizing phenol, a major substrate for both enzymes. Under ideal conditions, the hydroxylated product yield was ∼50% of the diiron centers for both systems, suggesting that these enzymes operate by half-sites reactivity mechanisms. Single-turnover studies indicated that the PH and ToMO electron-transfer components exert regulatory effects on substrate oxidation processes taking place at the hydroxylase actives sites, most likely through allostery. Steady state NADH consumption assays showed that the regulatory proteins facilitate the electron-transfer step in the hydrocarbon oxidation cycle in the absence of phenol. Under these conditions, electron consumption is coupled to H(2)O(2) formation in a hydroxylase-dependent manner. Mechanistic implications of these results are discussed.

Active site threonine facilitates proton transfer during dioxygen activation at the diiron center of toluene/o-xylene monooxygenase hydroxylase

Toluene/o-xylene monooxygenase hydroxylase (ToMOH), a diiron-containing enzyme, can activate dioxygen to oxidize aromatic substrates. To elucidate the role of a strictly conserved T201 residue during dioxygen activation of the enzyme, T201S, T201G, T201C, and T201V variants of ToMOH were prepared by site-directed mutagenesis. X-ray crystal structures of all the variants were obtained. Steady-state activity, regiospecificity, and single-turnover yields were also determined for the T201 mutants. Dioxygen activation by the reduced T201 variants was explored by stopped-flow UV-vis and Mössbauer spectroscopy. These studies demonstrate that the dioxygen activation mechanism is preserved in all T201 variants; however, both the formation and decay kinetics of a peroxodiiron(III) intermediate, T201(peroxo), were greatly altered, revealing that T201 is critically involved in dioxygen activation. A comparison of the kinetics of O(2) activation in the T201S, T201C, and T201G variants under various reaction conditions revealed that T201 plays a major role in proton transfer, which is required to generate the peroxodiiron(III) intermediate. A mechanism is postulated for dioxygen activation, and possible structures of oxygenated intermediates are discussed.

Recombinant expression, purification, and characterization of ThmD, the oxidoreductase component of tetrahydrofuran monooxygenase

Tetrahydrofuran monooxygenase (Thm) catalyzes the NADH-and oxygen-dependent hydroxylation of tetrahydrofuran to 2-hydroxytetrahydrofuran. Thm is composed of a hydroxylase enzyme, a regulatory subunit, and an oxidoreductase named ThmD. ThmD was expressed in Escherichia coli as a fusion to maltose-binding protein (MBP) and isolated to homogeneity after removal of the MBP. Purified ThmD contains covalently bound FAD, [2Fe-2S] center, and was shown to use ferricyanide, cytochrome c, 2,6-dichloroindophenol, and to a lesser extent, oxygen as surrogate electron acceptors. ThmD displays 160-fold preference for NADH over NADPH and functions as a monomer. The flavin-binding domain of ThmD (ThmD-FD) was purified and characterized. ThmD-FD displayed similar activity as the full-length ThmD and showed a unique flavin spectrum with a major peak at 463nm and a small peak at 396 nm. Computational modeling and mutagenesis analyses suggest a novel three-dimensional fold or covalent flavin attachment in ThmD.