Functional annotation and characterization of 3-hydroxybenzoate 6-hydroxylase from Rhodococcus jostii RHA1

The genome of Rhodococcus jostii RHA1 contains an unusually large number of oxygenase encoding genes. Many of these genes have yet an unknown function, implying that a notable part of the biochemical and catabolic biodiversity of this Gram-positive soil actinomycete is still elusive. Here we present a multiple sequence alignment and phylogenetic analysis of putative R. jostii RHA1 flavoprotein hydroxylases. Out of 18 candidate sequences, three hydroxylases are absent in other available Rhodococcus genomes. In addition, we report the biochemical characterization of 3-hydroxybenzoate 6-hydroxylase (3HB6H), a gentisate-producing enzyme originally mis-annotated as salicylate hydroxylase. R. jostii RHA1 3HB6H expressed in Escherichia coli is a homodimer with each 47 kDa subunit containing a non-covalently bound FAD cofactor. The enzyme has a pH optimum around pH 8.3 and prefers NADH as external electron donor. 3HB6H is active with a series of 3-hydroxybenzoate analogues, bearing substituents in ortho- or meta-position of the aromatic ring. Gentisate, the physiological product, is a non-substrate effector of 3HB6H. This compound is not hydroxylated but strongly stimulates the NADH oxidase activity of the enzyme.     

Reduction kinetics of 3-hydroxybenzoate 6-hydroxylase from Rhodococcus jostii RHA1

3-Hydroxybenzoate 6-hydroxylase (3HB6H) from Rhodococcus jostii RHA1 is a nicotinamide adenine dinucleotide (NADH)-specific flavoprotein monooxygenase involved in microbial aromatic degradation. The enzyme catalyzes the para hydroxylation of 3-hydroxybenzoate (3-HB) to 2,5-dihydroxybenzoate (2,5-DHB), the ring-fission fuel of the gentisate pathway. In this study, the kinetics of reduction of the enzyme-bound flavin by NADH was investigated at pH 8.0 using a stopped-flow spectrophotometer, and the data were analyzed comprehensively according to kinetic derivations and simulations. Observed rate constants for reduction of the free enzyme by NADH under anaerobic conditions were linearly dependent on NADH concentrations, consistent with a one-step irreversible reduction model with a bimolecular rate constant of 43 ± 2 M(-1) s(-1). In the presence of 3-HB, observed rate constants for flavin reduction were hyperbolically dependent on NADH concentrations and approached a limiting value of 48 ± 2 s(-1). At saturating concentrations of NADH (10 mM) and 3-HB (10 mM), the reduction rate constant is ~51 s(-1), whereas without 3-HB, the rate constant is 0.43 s(-1) at a similar NADH concentration. A similar stimulation of flavin reduction was found for the enzyme-product (2,5-DHB) complex, with a rate constant of 45 ± 2 s(-1). The rate enhancement induced by aromatic ligands is not due to a thermodynamic driving force because Em 0 for the enzyme-substrate complex is -179 ± 1 mV compared to an E(m)(0) of -175 ± 2 mV for the free enzyme. It is proposed that the reduction mechanism of 3HB6H involves an isomerization of the initial enzyme-ligand complex to a fully activated form before flavin reduction takes place.     

Crystal Structure of 3-Hydroxybenzoate 6-Hydroxylase Uncovers Lipid-assisted Flavoprotein Strategy for Regioselective Aromatic Hydroxylation

3-Hydroxybenzoate 6-hydroxylase (3HB6H) from Rhodococcus jostii RHA1 is a dimeric flavoprotein that catalyzes the NADH- and oxygen-dependent para-hydroxylation of 3-hydroxybenzoate to 2,5-dihydroxybenzoate. In this study, we report the crystal structure of 3HB6H as expressed in Escherichia coli. The overall fold of 3HB6H is similar to that of p-hydroxybenzoate hydroxylase and other flavoprotein aromatic hydroxylases. Unexpectedly, a lipid ligand is bound to each 3HB6H monomer. Mass spectral analysis identified the ligand as a mixture of phosphatidylglycerol and phosphatidylethanolamine. The fatty acid chains occupy hydrophobic channels that deeply penetrate into the interior of the substrate-binding domain of each subunit, whereas the hydrophilic part is exposed on the protein surface, connecting the dimerization domains via a few interactions. Most remarkably, the terminal part of a phospholipid acyl chain is directly involved in the substrate-binding site. Co-crystallized chloride ion and the crystal structure of the H213S variant with bound 3-hydroxybenzoate provide hints about oxygen activation and substrate hydroxylation. Essential roles are played by His-213 in catalysis and Tyr-105 in substrate binding. This phospholipid-assisted strategy to control regioselective aromatic hydroxylation is of relevance for optimization of flavin-dependent biocatalysts.     

Comparative and functional genomics of Rhodococcus opacus PD630 for biofuels development

The Actinomycetales bacteria Rhodococcus opacus PD630 and Rhodococcus jostii RHA1 bioconvert a diverse range of organic substrates through lipid biosynthesis into large quantities of energy-rich triacylglycerols (TAGs). To describe the genetic basis of the Rhodococcus oleaginous metabolism, we sequenced and performed comparative analysis of the 9.27 Mb R. opacus PD630 genome. Metabolic-reconstruction assigned 2017 enzymatic reactions to the 8632 R. opacus PD630 genes we identified. Of these, 261 genes were implicated in the R. opacus PD630 TAGs cycle by metabolic reconstruction and gene family analysis. Rhodococcus synthesizes uncommon straight-chain odd-carbon fatty acids in high abundance and stores them as TAGs. We have identified these to be pentadecanoic, heptadecanoic, and cis-heptadecenoic acids. To identify bioconversion pathways, we screened R. opacus PD630, R. jostii RHA1, Ralstonia eutropha H16, and C. glutamicum 13032 for growth on 190 compounds. The results of the catabolic screen, phylogenetic analysis of the TAGs cycle enzymes, and metabolic product characterizations were integrated into a working model of prokaryotic oleaginy.     

Adventures in Rhodococcus - from steroids to explosives

Rhodococcus is a genus of mycolic-acid-containing actinomycetes that utilize a remarkable variety of organic compounds as growth substrates. This degradation helps maintain the global carbon cycle and has increasing applications ranging from the biodegradation of pollutants to the biocatalytic production of drugs and hormones. We have been using Rhodococcus jostii RHA1 as a model organism to understand the catabolic versatility of Rhodococcus and related bacteria. Our approach is exemplified by the discovery of a cluster of genes specifying the catabolism of cholesterol. This degradation proceeds via β-oxidative degradation of the side chain and O2-dependent cleavage of steroid ring A in a process similar to bacterial degradation of aromatic compounds. The pathway is widespread in Actinobacteria and is critical to the pathogenesis of Mycobacterium tuberculosis, arguably the world's most successful pathogen. The close similarity of some of these enzymes with biphenyl- and polychlorinated-biphenyl-degrading enzymes that we have characterized is facilitating inhibitor design. Our studies in RHA1 have also provided important insights into a number of novel metalloenzymes and their biosynthesis, such as acetonitrile hydratase (ANHase), a cobalt-containing enzyme with no significant sequence identity with characterized nitrile hydratases. Molecular genetic and biochemical studies have identified AnhE as a dimeric metallochaperone that delivers cobalt to ANHase, enabling its maturation in vivo. Other metalloenzymes we are characterizing include N-acetylmuramic acid hydroxylase, which catalyzes an unusual hydroxylation of the rhodococcal and mycobacterial peptidoglycan, and 2 RHA1 dye-decolorizing peroxidases. Using molecular genetic and biochemical approaches, we have demonstrated that one of these enzymes is involved in the degradation of lignin. Overall, our studies are providing fundamental insights into a range of catabolic processes that have a wide variety of applications.     

Expression in Escherichia coli of biphenyl 2,3-dioxygenase genes from a Gram-positive polychlorinated biphenyl degrader, Rhodococcus jostii RHA1

Rhodococcus jostii RHA1 is a polychlorinated biphenyl degrader. Multi-component biphenyl 2,3-dioxygenase (BphA) genes of RHA1 encode large and small subunits of oxygenase component and ferredoxin and reductase components. They did not express enzyme activity in Escherichia coli. To obtain BphA activity in E. coli, hybrid BphA gene derivatives were constructed by replacing ferredoxin and/or reductase component genes of RHA1 with those of Pseudomonas pseudoalcaligenes KF707. The results obtained indicate a lack of catalytic activity of the RHA1 ferredoxin component gene, bphAc in E. coli. To determine the cause of inability of RHA1 bphAc to express in E. coli, the bphAc gene was introduced into Rosetta (DE3) pLacI, which has extra tRNA genes for rare codons in E. coli. The resulting strain abundantly produced the bphAc product, and showed activity. These results suggest that codon usage bias is involved in inability of RHA1 bphAc to express its catalytic activity in E. coli.     

Vanillin catabolism in Rhodococcus jostii RHA1

Genes encoding vanillin dehydrogenase (vdh) and vanillate O-demethylase (vanAB) were identified in Rhodococcus jostii RHA1 using gene disruption and enzyme activities. During growth on vanillin or vanillate, vanA was highly upregulated while vdh was not. This study contributes to our understanding of lignin degradation by RHA1 and other actinomycetes.     

Engineering of a xylose metabolic pathway in Rhodococcus strains

The two metabolically versatile actinobacteria Rhodococcus opacus PD630 and R. jostii RHA1 can efficiently convert diverse organic substrates into neutral lipids mainly consisting of triacylglycerol (TAG), the precursor of energy-rich hydrocarbon. Neither, however, is able to utilize xylose, the important component present in lignocellulosic biomass, as the carbon source for growth and lipid accumulation. In order to broaden their substrate utilization range, the metabolic pathway of d-xylose utilization was introduced into these two strains. This was accomplished by heterogenous expression of two well-selected genes, xylA, encoding xylose isomerase, and xylB, encoding xylulokinase from Streptomyces lividans TK23, under the control of the tac promoter with an Escherichia coli-Rhodococcus shuttle vector. The recombinant R. jostii RHA1 bearing xylA could grow on xylose as the sole carbon source, and additional expression of xylB further improved the biomass yield. The recombinant could consume both glucose and xylose in the sugar mixture, although xylose metabolism was still affected by the presence of glucose. The xylose metabolic pathway was also introduced into the high-lipid-producing strain R. opacus PD630 by expression of xylA and xylB. Under nitrogen-limited conditions, the fatty acid composition was determined, and lipid produced from xylose by recombinants of R. jostii RHA1 and R. opacus PD630 carrying xylA and xylB represented up to 52.5% and 68.3% of the cell dry weight (CDW), respectively. This work demonstrates that it is feasible to produce lipid from the sugars, including xylose, derived from renewable feedstock by genetic modification of rhodococcus strains.     

The Reaction Kinetics of 3-Hydroxybenzoate 6-Hydroxylase from Rhodococcus jostii RHA1 Provide an Understanding of the para-Hydroxylation Enzyme Catalytic Cycle

3-Hydroxybenzoate 6-hydroxylase (3HB6H) from Rhodococcus jostii RHA1 is an NADH-specific flavoprotein monooxygenase that catalyzes the para-hydroxylation of 3-hydroxybenzoate (3HB) to form 2,5-dihydroxybenzoate (2,5-DHB). Based on results from stopped-flow spectrophotometry, the reduced enzyme-3HB complex reacts with oxygen to form a C4a-peroxy flavin with a rate constant of 1.13 ± 0.01 × 10⁶ m⁻¹ s⁻¹ (pH 8.0, 4 °C). This intermediate is subsequently protonated to form a C4a-hydroperoxyflavin with a rate constant of 96 ± 3 s⁻¹. This step shows a solvent kinetic isotope effect of 1.7. Based on rapid-quench measurements, the hydroxylation occurs with a rate constant of 36 ± 2 s⁻¹. 3HB6H does not exhibit substrate inhibition on the flavin oxidation step, a common characteristic found in most ortho-hydroxylation enzymes. The apparent k_cat at saturating concentrations of 3HB, NADH, and oxygen is 6.49 ± 0.02 s⁻¹. Pre-steady state and steady-state kinetic data were used to construct the catalytic cycle of the reaction. The data indicate that the steps of product release (11.7 s⁻¹) and hydroxylation (36 ± 2 s⁻¹) partially control the overall turnover.     

Identification of a toluene-degrading bacterium from a soil sample through H(2)(18)O DNA stable isotope probing

DNA stable isotope probing (DNA-SIP) with H(2)(18)O was used to identify a toluene-degrading bacterium in soil amended with 48 ppm toluene. After quantification of toluene degradation rates in soil, DNA was extracted from soil incubated with H(2)(18)O, H(2)(16)O, H(2)(16)O and 48 ppm toluene, or H(2)(18)O and 48 ppm toluene. A single DNA band formed along a cesium chloride gradient after isopycnic centrifugation of extracts from soils incubated with H(2)(16)O. With extracts from soils to which only H(2)(18)O was added, two distinct DNA bands formed, while three bands formed when DNA extracted from soil incubated with both H(2)(18)O and toluene was analyzed. We suggest that this third band formed because toluene does not contain any oxygen atoms and toluene-degrading organisms had to transfer oxygen atoms from H(2)(18)O into metabolic intermediates to form nucleic acids de novo. We extracted the third DNA band and amplified a large fraction of the bacterial 16S rRNA gene. Direct sequencing of the PCR product obtained from the labeled DNA, as well as cloned 16S rRNA amplicons, identified a known toluene degrader, Rhodococcus jostii RHA1. A toluene-degrading bacterial strain was subsequently isolated from soil and shown to be Rhodococcus jostii RHA1. Finally, quantitative real-time PCR analysis showed that the abundance of the 16S rRNA gene of Rhodococcus jostii RHA1 increased in soil after toluene exposure but not in soils from which toluene was withheld. This study indicates that H(2)(18)O DNA-SIP can be a useful method for identifying pollutant-degrading bacteria in soil.     

Inactivation of three monohydroxybenzoate mono-oxygenases from Rhodococcus erythropolis: The role of an arginine residue in the substrate-binding domain of p-hydroxybenzoate 3-hydroxylase

Abstract p-Hydroxybenzoate 3-hydroxylase from Rhodococcus erythropolis was inactivated by 2,3-butanedione (BD), phenylglyoxal (PGO), and other chemical reagents. p -Hydroxybenzoate and NADH protected the enzyme from inactivation by BD. Judging from the amino acid composition of BD-treated enzyme in the presence and absence of p -hydroxybenzoate, one essential arginine residue in substrate-binding domain of the enzyme was shown to be essential to the binding of p -hydrozybenzoate to the enzyme. Salicylate 5-hydroxylase and m -hydroxybenzoate 6-hydroxylase from R. erythropolis were hardly inactivated. Neither of these two enzymes was considered to have a functional arginine residue required for interaction with the substrate.     

3-Hydroxybenzoate 6-hydroxylase from Pseudomonas aeruginosa

An inducible 3-hydroxybenzoate 6-hydroxylase has been purified to homogeneity from $\text{Pseudomonas aeruginosa}$. It contains FAD as a prosthetic group. 3-Hydroxybenzoate is quantitatively hydroxylated to give gentisate with equimolar consumptions of NADH and O₂. NADPH will substitute as an electron donor, and several aromatic analogues of 3-hydroxybenzoate stimulate reduced nucleotide oxidation by the enzyme with formation of both hydrogen peroxide and hydroxylated products. Of various analogues of 3-hydroxybenzoate, those substituted in 2,4,5 and 6-positions are competent substrates; partial uncoupling of electron flow from hydroxylation with concomitant formation of hydrogen peroxide and “gentisates” occurs. The “natural” product of the reaction, gentisate, is an effector in that it stimulates NADH oxidation with the formation of hydrogen peroxide. 3-hydroxybenzoate 6-hydroxylase thus resembles other flavoprotein hydroxylases in the general regulatory properties dictated by their aromatic substrates, pseudosubstrates or effectors.     

Identification of DypB from Rhodococcus jostii RHA1 as a lignin peroxidase

Rhodococcus jostii RHA1, a polychlorinated biphenyl-degrading soil bacterium whose genome has been sequenced, shows lignin degrading activity in two recently developed spectrophotometric assays. Bioinformatic analysis reveals two unannotated peroxidase genes present in the genome of R. jostii RHA1 with sequence similarity to open reading frames in other lignin-degrading microbes. They are members of the Dyp peroxidase family and were annotated as DypA and DypB, on the basis of bioinformatic analysis. Assay of gene deletion mutants using a colorimetric lignin degradation assay reveals that a ΔdypB mutant shows greatly reduced lignin degradation activity, consistent with a role in lignin breakdown. Recombinant DypB protein shows activity in the colorimetric assay and shows Michaelis-Menten kinetic behavior using Kraft lignin as a substrate. DypB is activated by Mn(2+) by 5-23-fold using a range of assay substrates, and breakdown of wheat straw lignocellulose by recombinant DypB is observed over 24-48 h in the presence of 1 mM MnCl(2). Incubation of recombinant DypB with a β-aryl ether lignin model compound shows time-dependent turnover, giving vanillin as a product, indicating that C(α)-C(β) bond cleavage has taken place. This reaction is inhibited by addition of diaphorase, consistent with a radical mechanism for C-C bond cleavage. Stopped-flow kinetic analysis of the DypB-catalyzed reaction shows reaction between the intermediate compound I (397 nm) and either Mn(II) (k(obs) = 2.35 s(-1)) or the β-aryl ether (k(obs) = 3.10 s(-1)), in the latter case also showing a transient at 417 nm, consistent with a compound II intermediate. These results indicate that DypB has a significant role in lignin degradation in R. jostii RHA1, is able to oxidize both polymeric lignin and a lignin model compound, and appears to have both Mn(II) and lignin oxidation sites. This is the first detailed characterization of a recombinant bacterial lignin peroxidase.     

Polychlorinated biphenyl/biphenyl degrading gene clusters in Rhodococcus sp. K37, HA99, and TA431 are different from well-known bph gene clusters of Rhodococci

Four kinds of polychlorinated biphenyl (PCB)-degrading Rhodococcus sp. (TA421, TA431, HA99, and K37) have been isolated from termite ecosystem and under alkaline condition. The bph gene cluster involved in the degradation of PCB/biphenyl has been analyzed in strain TA421. This gene cluster was highly homologous to bph gene clusters in R. globerulus P6 and Rhodococcus sp. RHA1. In this study, we cloned and analyzed the bph gene cluster essential to PCB/biphenyl degradation from R. rhodochrous K37. The order of the genes and the sequence were different in K37 than in P6, RHA1, and TA421. The bphC8(K37) gene was more homologous to the meta-cleavage enzyme involved in phenanthrene metabolism than bphC genes involved in biphenyl metabolism. Two other Rhodococcus strains (HA99 and TA431) had PCB/biphenyl degradation gene clusters similar to that in K37. These findings suggest that these bph gene clusters evolved separately from the well-known bph gene clusters of PCB/biphenyl degraders.     

Molecular and biochemical characterization of 3-hydroxybenzoate 6-hydroxylase from Polaromonas naphthalenivorans CJ2

Prior research revealed that Polaromonas naphthalenivorans CJ2 carries and expresses genes encoding the gentisate metabolic pathway for naphthalene. These metabolic genes are split into two clusters, comprising nagRAaGHAbAcAdBFCQEDJI'-orf1-tnpA and nagR2-orf2I'KL (C. O. Jeon, M. Park, H. Ro, W. Park, and E. L. Madsen, Appl. Environ. Microbiol. 72:1086-1095, 2006). BLAST homology searches of sequences in GenBank indicated that the orf2 gene from the small cluster likely encoded a salicylate 5-hydroxylase, presumed to catalyze the conversion of salicylate into gentisate. Here, we report physiological and genetic evidence that orf2 does not encode salicylate 5-hydroxylase. Instead, we have found that orf2 encodes 3-hydroxybenzoate 6-hydroxylase, the enzyme which catalyzes the NADH-dependent conversion of 3-hydroxybenzoate into gentisate. Accordingly, we have renamed orf2 nagX. After expression in Escherichia coli, the NagX enzyme had an approximate molecular mass of 43 kDa, as estimated by gel filtration, and was probably a monomeric protein. The enzyme was able to convert 3-hydroxybenzoate into gentisate without salicylate 5-hydroxylase activity. Like other 3-hydroxybenzoate 6-hydroxylases, NagX utilized both NADH and NADPH as electron donors and exhibited a yellowish color, indicative of a bound flavin adenine dinucleotide. An engineered mutant of P. naphthalenivorans CJ2 defective in nagX failed to grow on 3-hydroxybenzoate but grew normally on naphthalene. These results indicate that the previously described small catabolic cluster in strain CJ2 may be multifunctional and is essential for the degradation of 3-hydroxybenzoate. Because nagX and an adjacent MarR-type regulatory gene are both closely related to homologues in Azoarcus species, this study raises questions about horizontal gene transfer events that contribute to operon evolution.     

Pseudomonas cepacia 3-hydroxybenzoate 6-hydroxylase: induction, purification, and characterization

A single strain of Pseudomonas cepacia cells was differentially induced to synthesize salicylate hydroxylase, 3-hydroxybenzoate 6-hydroxylase, or 4-hydroxybenzoate 3-hydroxylase. A procedure was developed for the purification of 3-hydroxybenzoate 6-hydroxylase to apparent homogeneity. The purified hydroxylase appears to be a monomer with a molecular weight of about 44,000 and exhibits optimal activity near pH 8. The hydroxylase contains one FAD per enzyme molecule and utilizes NADH and NADPH with similar efficiencies. The reaction stoichiometry for this enzyme has been determined. In comparison with other aromatic flavohydroxylases, this enzyme is unique in inserting a new hydroxyl group to the substrate at a position para to an existing one.     

Physicochemical and immunochemical characterization of salicylate 5-hydroxylase,m-hydroxybenzoate 6-hydroxylase and p-hydroxybenzoate 3-hydroxylase from Rhodococcus erythropolis

Salicylate 5-hydroxylase (SAL5H), m-hydroxybenzoate 6-hydroxylase (MHB6H), and p-hydroxybenzoate 3-hydroxylase (PHB3H) from Gram-positive Rhodococcus erythropolis strain S1 were characterized physicochemically and immunochemically. The subunit size and amino acid composition of SAL5H, MHB6H, and PHB3H from strain S1 showed properties similar to those of other flavin-containing aromatic compound monooxygenases such as p-hydroxybenzoate hydroxylase and salicylate 1-hydroxylase (SAL1H), belonging to p-hydroxybenzoate hydroxylase-class, except for homotetrameric structure and cofactor specficity. The N-terminal amino acid sequence of MHB6H from strain S1 indicated significant similarity of ADP-binding region in the N-terminal portion of the enzyme with that known for SAL1H from Pseudomonas putida. Immunochemical properties, determined while conducting serological experiments, showed SAL5H and MHB6H from strain S1 to be immunologically different from PHB3H from strain S1, while SAL5H and MHB6H to apparently share partial antigenic determinants.