Crystal Structures and Functional Studies of T4moD, the Toluene 4-Monooxygenase Catalytic Effector Protein,

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/.

Purification and Characterization of the Isoprene Monooxygenase from Rhodococcus sp. Strain AD45

Isoprene (2-methyl-1,3-butadiene) is a climate-active gas released to the atmosphere in large quantities, comparable to methane in magnitude. Several bacteria have been isolated which can grow on isoprene as a sole carbon and energy source, but very little information is available about the degradation of isoprene by these bacteria at the biochemical level. Isoprene utilization is dependent on a multistep pathway, with the first step being the oxidation of isoprene to epoxy-isoprene. This is catalyzed by a four-component soluble diiron monooxygenase, isoprene monooxygenase (IsoMO). IsoMO is a six-protein complex comprising an oxygenase (IsoABE), containing the di-iron active site, a Rieske-type ferredoxin (IsoC), a NADH reductase (IsoF), and a coupling/effector protein (IsoD), homologous to the soluble methane monooxygenase and alkene/aromatic monooxygenases. Here, we describe the purification of the IsoMO components from Rhodococcus sp. AD45 and reconstitution of isoprene-oxidation activity in vitro. Some IsoMO components were expressed and purified from the homologous host Rhodococcus sp. AD45-ID, a Rhodococcus sp. AD45 strain lacking the megaplasmid which contains the isoprene metabolic gene cluster. Others were expressed in Escherichia coli and purified as fusion proteins. We describe the characterization of these purified components and demonstrate their activity when combined with Rhodococcus sp. AD45 cell lysate. Demonstration of IsoMO activity in vitro provides a platform for further biochemical and biophysical characterization of this novel soluble diiron center monooxygenase, facilitating new insights into the enzymatic basis for the bacterial degradation of isoprene. IMPORTANCE Isoprene is a highly abundant climate-active gas and a carbon source for some bacteria. Analyses of the genes encoding isoprene monooxygenase (IsoMO) indicate this enzyme is a soluble diiron center monooxygenase in the same family of oxygenases as soluble methane monooxygenase, alkene monooxygenase, and toluene monooxygenase. We report the initial biochemical characterization of IsoMO from Rhodococcus, the first from any bacterium, describing the challenging purification and reconstitution of in vitro activity of its four components. This study lays the foundation for future detailed mechanistic studies of IsoMO, a key enzyme in the global isoprene cycle.

An Overview of the Electron-Transfer Proteins That Activate Alkane Monooxygenase (AlkB)

Alkane-oxidizing enzymes play an important role in the global carbon cycle. Alkane monooxygenase (AlkB) oxidizes most of the medium-chain length alkanes in the environment. The first AlkB identified was from P. putida GPo1 (initially known as P. oleovorans) in the early 1970s, and it continues to be the family member about which the most is known. This AlkB is found as part of the OCT operon, in which all of the key proteins required for growth on alkanes are present. The AlkB catalytic cycle requires that the diiron active site be reduced. In P. putida GPo1, electrons originate from NADH and arrive at AlkB via the intermediacy of a flavin reductase and an iron–sulfur protein (a rubredoxin). In this Mini Review, we will review what is known about the canonical arrangement of electron-transfer proteins that activate AlkB and, more importantly, point to several other arrangements that are possible. These other arrangements include the presence of a simpler rubredoxin than what is found in the canonical arrangement, as well as two other classes of AlkBs with fused electron-transfer partners. In one class, a rubredoxin is fused to the hydroxylase and in another less well-explored class, a ferredoxin reductase and a ferredoxin are fused to the hydroxylase. We review what is known about the biochemistry of these electron-transfer proteins, speculate on the biological significance of this diversity, and point to key questions for future research.

BesC Initiates C-C Cleavage through a Substrate-Triggered and Reactive Diferric-Peroxo Intermediate

BesC catalyzes the iron- and O₂-dependent cleavage of 4-chloro-l-lysine to form 4-chloro-l-allylglycine, formaldehyde, and ammonia. This process is a critical step for a biosynthetic pathway that generates a terminal alkyne amino acid which can be leveraged as a useful bio-orthogonal handle for protein labeling. As a member of an emerging family of diiron enzymes that are typified by their heme oxygenase-like fold and a very similar set of coordinating ligands, recently termed HDOs, BesC performs an unusual type of carbon-carbon cleavage reaction that is a significant departure from reactions catalyzed by canonical dinuclear-iron enzymes. Here, we show that BesC activates O₂ in a substrate-gated manner to generate a diferric-peroxo intermediate. Examination of the reactivity of the peroxo intermediate with a series of lysine derivatives demonstrates that BesC initiates this unique reaction trajectory via cleavage of the C4-H bond; this process represents the rate-limiting step in a single turnover reaction. The observed reactivity of BesC represents the first example of a dinuclear-iron enzyme that utilizes a diferric-peroxo intermediate to capably cleave a C-H bond as part of its native function, thus circumventing the formation of a high-valent intermediate more commonly associated with substrate monooxygenations.

MhpA Is a Hydroxylase Catalyzing the Initial Reaction of 3-(3-Hydroxyphenyl)Propionate Catabolism in Escherichia coli K-12

Escherichia coli K-12 and some other strains have been reported to be capable of utilizing 3-(3-hydroxyphenyl)propionate (3HPP), one of the phenylpropanoids from lignin. Although other enzymes involved in 3HPP catabolism and their corresponding genes from its degraders have been identified, 3HPP 2-hydroxylase, catalyzing the first step of its catabolism, has yet to be functionally identified at biochemical and genetic levels. In this study, we investigated the function and characteristics of MhpA from E. coli strain K-12 (MhpA_K-12). Gene deletion and complementation showed that mhpA was vital for its growth on 3HPP, but the mhpA deletion strain was still able to grow on 3-(2,3-dihydroxyphenyl)propionate (DHPP), the hydroxylation product transformed from 3HPP by MhpA_K-12 MhpA_K-12 was overexpressed and purified, and it was likely a polymer and tightly bound with an approximately equal number of moles of FAD. Using NADH or NADPH as a cofactor, purified MhpA_K-12 catalyzed the conversion of 3HPP to DHPP at a similar efficiency. The conversion from 3HPP to DHPP by purified MhpA_K-12 was confirmed using high-performance liquid chromatography and liquid chromatography-mass spectrometry. Bioinformatics analysis indicated that MhpA_K-12 and its putative homologues belonged to taxa that were phylogenetically distant from functionally identified FAD-containing monooxygenases (hydroxylases). Interestingly, MhpA_K-12 has approximately an extra 150 residues at its C terminus in comparison to its close homologues, but its truncated versions MhpA_K-12 ⁴⁰⁰ and MhpA_K-12 ⁴⁸⁰ (with 154 and 74 residues deleted from the C terminus, respectively) both lost their activities. Thus, MhpA_K-12 has been confirmed to be a 3HPP 2-hydroxylase catalyzing the conversion of 3HPP to DHPP, the initial reaction of 3HPP degradation.IMPORTANCE Phenylpropionate and its hydroxylated derivatives resulted from lignin degradation ubiquitously exist on the Earth. A number of bacterial strains have the ability to grow on 3HPP, one of the above derivatives. The hydroxylation was thought to be the initial and vital step for its aerobic catabolism via the meta pathway. The significance of our research is the functional identification and characterization of the purified 3HPP 2-hydroxylase MhpA from Escherichia coli K-12 at biochemical and genetic levels, since this enzyme has not previously been expressed from its encoding gene, purified, and characterized in any bacteria. It will not only fill a gap in our understanding of 3HPP 2-hydroxylase and its corresponding gene for the critical step in microbial 3HPP catabolism but also provide another example of the diversity of microbial degradation of plant-derived phenylpropionate and its hydroxylated derivatives.

Deciphering Cryptic Binding Sites on Proteins by Mixed-Solvent Molecular Dynamics

In recent years, molecular dynamics simulations of proteins in explicit mixed solvents have been applied to various problems in protein biophysics and drug discovery, including protein folding, protein surface characterization, fragment screening, allostery, and druggability assessment. In this study, we perform a systematic study on how mixtures of organic solvent probes in water can reveal cryptic ligand binding pockets that are not evident in crystal structures of apo proteins. We examine a diverse set of eight PDB proteins that show pocket opening induced by ligand binding and investigate whether solvent MD simulations on the apo structures can induce the binding site observed in the holo structures. The cosolvent simulations were found to induce conformational changes on the protein surface, which were characterized and compared with the holo structures. Analyses of the biological systems, choice of probes and concentrations, druggability of the resulting induced pockets, and application to drug discovery are discussed here.

Structure/function correlations over binuclear non-heme iron active sites

Binuclear non-heme iron enzymes activate O2 to perform diverse chemistries. Three different structural mechanisms of O2 binding to a coupled binuclear iron site have been identified utilizing variable-temperature, variable-field magnetic circular dichroism spectroscopy (VTVH MCD). For the μ-OH-bridged Fe(II)2 site in hemerythrin, O2 binds terminally to a five-coordinate Fe(II) center as hydroperoxide with the proton deriving from the μ-OH bridge and the second electron transferring through the resulting μ-oxo superexchange pathway from the second coordinatively saturated Fe(II) center in a proton-coupled electron transfer process. For carboxylate-only-bridged Fe(II)2 sites, O2 binding as a bridged peroxide requires both Fe(II) centers to be coordinatively unsaturated and has good frontier orbital overlap with the two orthogonal O2 π* orbitals to form peroxo-bridged Fe(III)2 intermediates. Alternatively, carboxylate-only-bridged Fe(II)2 sites with only a single open coordination position on an Fe(II) enable the one-electron formation of Fe(III)-O2 (-) or Fe(III)-NO(-) species. Finally, for the peroxo-bridged Fe(III)2 intermediates, further activation is necessary for their reactivities in one-electron reduction and electrophilic aromatic substitution, and a strategy consistent with existing spectral data is discussed.

Structure of T4moF, the Toluene 4-Monooxygenase Ferredoxin Oxidoreductase

The 1.6 Å crystal structure of toluene 4-monooxygenase reductase T4moF is reported. The structure includes ferredoxin, flavin, and NADH binding domains. The position of the ferredoxin domain relative to the other two domains represents a new configuration for the iron-sulfur flavoprotein family. Close contacts between the C8 methyl group of FAD and [2Fe-2S] ligand Cys36-O represent a plausible pathway for electron transfer between the redox cofactors. Energy-minimized docking of NADH and calculation of hingelike motions between domains suggest how simple coordinated shifts of residues at the C-terminus of the enzyme could expose the N5 position of FAD for productive interaction with the nicotinamide ring. The domain configuration revealed by the T4moF structure provides an excellent steric and electrostatic match to the obligate electron acceptor, Rieske-type [2Fe-2S] ferredoxin T4moC. Protein-protein docking and energy minimization of the T4moFC complex indicate that T4moF [2Fe-2S] ligand Cys41 and T4moC [2Fe-2S] ligand His67, along with other electrostatic interactions between the protein partners, form the functional electron transfer interface.

Atomic picture of ligand migration in toluene 4-monooxygenase

Computational modeling combined with mutational and activity assays was used to underline the substrate migration pathways in toluene 4-monooxygenase, a member of the important family of bacterial multicomponent monooxygenases (BMMs). In all structurally defined BMM hydroxylases, several hydrophobic cavities in the α-subunit map a preserved path from the protein surface to the diiron active site. Our results confirm the presence of two pathways by which different aromatic molecules can enter/escape the active site. While the substrate is observed to enter from both channels, the more hydrophilic product is withdrawn mainly from the shorter channel ending at residues D285 and E214. The long channel ends in the vicinity of S395, whose variants have been seen to affect activity and specificity. These mutational effects are clearly reproduced and rationalized by the in silico studies. Furthermore, the combined computational and experimental results highlight the importance of residue F269, which is located at the intersection of the two channels.

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.

Structural consequences of effector protein complex formation in a diiron hydroxylase

Carboxylate-bridged diiron hydroxylases are multicomponent enzyme complexes responsible for the catabolism of a wide range of hydrocarbons and as such have drawn attention for their mechanism of action and potential uses in bioremediation and enzymatic synthesis. These enzyme complexes use a small molecular weight effector protein to modulate the function of the hydroxylase. However, the origin of these functional changes is poorly understood. Here, we report the structures of the biologically relevant effector protein-hydroxylase complex of toluene 4-monooxygenase in 2 redox states. The structures reveal a number of coordinated changes that occur up to 25 A from the active site and poise the diiron center for catalysis. The results provide a structural basis for the changes observed in a number of the measurable properties associated with effector protein binding. This description provides insight into the functional role of effector protein binding in all carboxylate-bridged diiron hydroxylases.

Soluble expression and purification of the oxidoreductase component of toluene 4-monooxygenase

Toluene 4-monooxygenase (T4MO) is a member of the bacterial multicomponent monooxygenases, an enzyme family that utilizes a soluble diiron hydroxylase to oxidize a variety of hydrocarbons as the initial step in their metabolism. The hydroxylases obtain reducing equivalents from NAD(P)H via an electron transfer chain that is initiated by an oxidoreductase containing an N-terminal ferredoxin domain and C-terminal flavin- and NAD-binding domains. T4moF, the NADH oxidoreductase of T4MO, was expressed as a soluble protein in Escherichia coli BL21(DE3) from the pUC-derived expression vector pRS205. This vector contains a lac promoter instead of a T7 promoter. A three step purification from the soluble cell lysate yielded approximately 1 mg of T4moF per gram of wet cell paste with greater than 90% purity. The purified protein contained 1 mol of FAD and 2 mol of Fe per mol of T4moF; quantitative EPR spectroscopy showed approximately 1 mol of the S=1/2 signal from the reduced [2Fe-2S] cluster per mol of T4moF. Steady state kinetic analysis of p-cresol formation activity treating T4moF as the variable substrate while all other proteins and substrates were held constant gave apparent K(M-) and apparent k(cat)-values of 0.15 microM and 3.0 s(-1), respectively. This expression system and purification allows for the recovery of the soluble oxidoreductase in yields that facilitate further biochemical and structural characterizations.

X-ray structure of a hydroxylase-regulatory protein complex from a hydrocarbon-oxidizing multicomponent monooxygenase, Pseudomonas sp. OX1 phenol hydroxylase

Phenol hydroxylase (PH) belongs to a family of bacterial multicomponent monooxygenases (BMMs) with carboxylate-bridged diiron active sites. Included are toluene/o-xylene (ToMO) and soluble methane (sMMO) monooxygenase. PH hydroxylates aromatic compounds, but unlike sMMO, it cannot oxidize alkanes despite having a similar dinuclear iron active site. Important for activity is formation of a complex between the hydroxylase and a regulatory protein component. To address how structural features of BMM hydroxylases and their component complexes may facilitate the catalytic mechanism and choice of substrate, we determined X-ray structures of native and SeMet forms of the PH hydroxylase (PHH) in complex with its regulatory protein (PHM) to 2.3 A resolution. PHM binds in a canyon on one side of the (alphabetagamma)2 PHH dimer, contacting alpha-subunit helices A, E, and F approximately 12 A above the diiron core. The structure of the dinuclear iron center in PHH resembles that of mixed-valent MMOH, suggesting an Fe(II)Fe(III) oxidation state. Helix E, which comprises part of the iron-coordinating four-helix bundle, has more pi-helical character than analogous E helices in MMOH and ToMOH lacking a bound regulatory protein. Consequently, conserved active site Thr and Asn residues translocate to the protein surface, and an approximately 6 A pore opens through the four-helix bundle. Of likely functional significance is a specific hydrogen bond formed between this Asn residue and a conserved Ser side chain on PHM. The PHM protein covers a putative docking site on PHH for the PH reductase, which transfers electrons to the PHH diiron center prior to O2 activation, suggesting that the regulatory component may function to block undesired reduction of oxygenated intermediates during the catalytic cycle. A series of hydrophobic cavities through the PHH alpha-subunit, analogous to those in MMOH, may facilitate movement of the substrate to and/or product from the active site pocket. Comparisons between the ToMOH and PHH structures provide insights into their substrate regiospecificities.

X-ray crystal structures of manganese(II)-reconstituted and native toluene/o-xylene monooxygenase hydroxylase reveal rotamer shifts in conserved residues and an enhanced view of the protein interior

We report the X-ray crystal structures of native and manganese(II)-reconstituted toluene/o-xylene monooxygenase hydroxylase (ToMOH) from Pseudomonas stutzeri OX1 to 1.85 and 2.20 A resolution, respectively. The structures reveal that reduction of the dimetallic active site is accompanied by a carboxylate shift and alteration of the coordination environment for dioxygen binding and activation. A rotamer shift in a strategically placed asparagine 202 accompanies dimetallic center reduction and is proposed to influence protein component interactions. This rotamer shift is conserved between ToMOH and the corresponding residue in methane monooxygenase hydroxylase (MMOH). Previously unidentified hydrophobic pockets similar to those present in MMOH are assigned.