Genetic and structural organization of the aminophenol catabolic operon and its implication for evolutionary process

Modulators of a robust and efficient metabolism: Perspective and insights from the Rid superfamily of proteins

Metabolism is an integrated network of biochemical pathways that assemble to generate the robust, responsive physiologies of microorganisms. Despite decades of fundamental studies on metabolic processes and pathways, our understanding of the nuance and complexity of metabolism remains incomplete. The ability to predict and model metabolic network structure, and its influence on cellular fitness, is complicated by the persistence of genes of unknown function, even in the best-studied model organisms. This review describes the definition and continuing study of the Rid superfamily of proteins. These studies are presented with a perspective that illustrates how metabolic complexity can complicate the assignment of function to uncharacterized genes. The Rid superfamily of proteins has been divided into eight subfamilies, including the well-studied RidA subfamily. Aside from the RidA proteins, which are present in all domains of life and prevent metabolic stress, most members of the Rid superfamily have no demonstrated physiological role. Recent progress on functional assignment supports the hypothesis that, overall, proteins in the Rid superfamily modulate metabolic processes to ensure optimal organismal fitness.

Construction of complete degradation pathway for nitrobenzene in Escherichia coli

Nitrobenzene is widely present in industrial wastewater and soil. Biodegradation has become an ideal method to remediate organic pollutants due to its low cost, high efficiency, and absence of secondary pollution. In the present study, 10 exogenous genes that can completely degrade nitrobenzene were introduced into Escherichia coli, and their successful expression in the strain was verified by fluorescence quantitative polymerase chain reaction and proteomic analysis. The results of the degradation experiment showed that the engineered strain could completely degrade 4 mM nitrobenzene within 8 h. The formation of intermediate metabolites was detected, and the final metabolites entered the E. coli tricarboxylic acid cycle smoothly. This process was discovered by isotope tracing method. Results indicated the integrality of the degradation pathway and the complete degradation of nitrobenzene. Finally, further experiments were conducted in soil to verify its degradation ability and showed that the engineered strain could also degrade 1 mM nitrobenzene within 10 h. In this study, engineered bacteria that can completely degrade nitrobenzene have been constructed successfully. The construction of remediation-engineered bacteria by synthetic biology laid the foundation for the industrial application of biological degradation of organic pollutants.

RidA Proteins Protect against Metabolic Damage by Reactive Intermediates

The Rid (YjgF/YER057c/UK114) protein superfamily was first defined by sequence homology with available protein sequences from bacteria, archaea, and eukaryotes (L. Parsons, N. Bonander, E. Eisenstein, M. Gilson, et al., Biochemistry 42:80-89, 2003, https://doi.org/10.1021/bi020541w). The archetypal subfamily, RidA (reactive intermediate deaminase A), is found in all domains of life, with the vast majority of free-living organisms carrying at least one RidA homolog. In over 2 decades, close to 100 reports have implicated Rid family members in cellular processes in prokaryotes, yeast, plants, and mammals. Functional roles have been proposed for Rid enzymes in amino acid biosynthesis, plant root development and nutrient acquisition, cellular respiration, and carcinogenesis. Despite the wealth of literature and over a dozen high-resolution structures of different RidA enzymes, their biochemical function remained elusive for decades. The function of the RidA protein was elucidated in a bacterial model system despite (i) a minimal phenotype of ridA mutants, (ii) the enzyme catalyzing a reaction believed to occur spontaneously, and (iii) confusing literature on the pleiotropic effects of RidA homologs in prokaryotes and eukaryotes. Subsequent work provided the physiological framework to support the RidA paradigm in Salmonella enterica by linking the phenotypes of mutants lacking ridA to the accumulation of the reactive metabolite 2-aminoacrylate (2AA), which damaged metabolic enzymes. Conservation of enamine/imine deaminase activity of RidA enzymes from all domains raises the likelihood that, despite the diverse phenotypes, the consequences when RidA is absent are due to accumulated 2AA (or a similar reactive enamine) and the diversity of metabolic phenotypes can be attributed to differences in metabolic network architecture. The discovery of the RidA paradigm in S. enterica laid a foundation for assessing the role of Rid enzymes in diverse organisms and contributed fundamental lessons on metabolic network evolution and diversity in microbes. This review describes the studies that defined the conserved function of RidA, the paradigm of enamine stress in S. enterica, and emerging studies that explore how this paradigm differs in other organisms. We focus primarily on the RidA subfamily, while remarking on our current understanding of the other Rid subfamilies. Finally, we describe the current status of the field and pose questions that will drive future studies on this widely conserved protein family to provide fundamental new metabolic information.

HosA, a MarR Family Transcriptional Regulator, Represses Nonoxidative Hydroxyarylic Acid Decarboxylase Operon and Is Modulated by 4-Hydroxybenzoic Acid

Members of the Multiple antibiotic resistance Regulator (MarR) family of DNA binding proteins regulate transcription of a wide array of genes required for virulence and pathogenicity of bacteria. The present study reports the molecular characterization of HosA (Homologue of SlyA), a MarR protein, with respect to its target gene, DNA recognition motif, and nature of its ligand. Through a comparative genomics approach, we demonstrate that hosA is in synteny with nonoxidative hydroxyarylic acid decarboxylase (HAD) operon and is present exclusively within the mutS-rpoS polymorphic region in nine different genera of Enterobacteriaceae family. Using molecular biology and biochemical approach, we demonstrate that HosA binds to a palindromic sequence downstream to the transcription start site of divergently transcribed nonoxidative HAD operon and represses its expression. Furthermore, in silico analysis showed that the recognition motif for HosA is highly conserved in the upstream region of divergently transcribed operon in different genera of Enterobacteriaceae family. A systematic chemical search for the physiological ligand revealed that 4-hydroxybenzoic acid (4-HBA) interacts with HosA and derepresses HosA mediated repression of the nonoxidative HAD operon. Based on our study, we propose a model for molecular mechanism underlying the regulation of nonoxidative HAD operon by HosA in Enterobacteriaceae family.

Indole-3-acetic acid in plant-microbe interactions

Indole-3-acetic acid (IAA) is an important phytohormone with the capacity to control plant development in both beneficial and deleterious ways. The ability to synthesize IAA is an attribute that many bacteria including both plant growth-promoters and phytopathogens possess. There are three main pathways through which IAA is synthesized; the indole-3-pyruvic acid, indole-3-acetamide and indole-3-acetonitrile pathways. This chapter reviews the factors that effect the production of this phytohormone, the role of IAA in bacterial physiology and in plant-microbe interactions including phytostimulation and phytopathogenesis.

Genomic and functional analyses of the 2-aminophenol catabolic pathway and partial conversion of its substrate into picolinic acid in Burkholderia xenovorans LB400

2-aminophenol (2-AP) is a toxic nitrogen-containing aromatic pollutant. Burkholderia xenovorans LB400 possess an amn gene cluster that encodes the 2-AP catabolic pathway. In this report, the functionality of the 2-aminophenol pathway of B. xenovorans strain LB400 was analyzed. The amnRJBACDFEHG cluster located at chromosome 1 encodes the enzymes for the degradation of 2-aminophenol. The absence of habA and habB genes in LB400 genome correlates with its no growth on nitrobenzene. RT-PCR analyses in strain LB400 showed the co-expression of amnJB, amnBAC, amnACD, amnDFE and amnEHG genes, suggesting that the amn cluster is an operon. RT-qPCR showed that the amnB gene expression was highly induced by 2-AP, whereas a basal constitutive expression was observed in glucose, indicating that these amn genes are regulated. We propose that the predicted MarR-type transcriptional regulator encoded by the amnR gene acts as repressor of the amn gene cluster using a MarR-type regulatory binding sequence. This report showed that LB400 resting cells degrade completely 2-AP. The amn gene cluster from strain LB400 is highly identical to the amn gene cluster from P. knackmussi strain B13, which could not grow on 2-AP. However, we demonstrate that B. xenovorans LB400 is able to grow using 2-AP as sole nitrogen source and glucose as sole carbon source. An amnBA (-) mutant of strain LB400 was unable to grow with 2-AP as nitrogen source and glucose as carbon source and to degrade 2-AP. This study showed that during LB400 growth on 2-AP this substrate was partially converted into picolinic acid (PA), a well-known antibiotic. The addition of PA at lag or mid-exponential phase inhibited LB400 growth. The MIC of PA for strain LB400 is 2 mM. Overall, these results demonstrate that B. xenovorans strain LB400 posses a functional 2-AP catabolic central pathway, which could lead to the production of picolinic acid.

Role of EctR as transcriptional regulator of ectoine biosynthesis genes in Methylophaga thalassica

In the halophilic aerobic methylotrophic bacterium Methylophaga thalassica, the genes encoding the enzymes for biosynthesis of the osmoprotectant ectoine were shown to be located in operon ectABC-ask. Transcription of the ect-operon was started from the two promoters homologous to the σ(70)-dependent promoter of Escherichia coli and regulated by protein EctR, whose encoding gene, ectR, is transcribed from three promoters. Genes homologous to ectR of methylotrophs were found in clusters of ectoine biosynthesis genes in some non-methylotrophic halophilic bacteria. EctR proteins of methylotrophic and heterotrophic halophiles belong to the MarR-family of transcriptional regulators but form a separate branch on the phylogenetic tree of the MarR proteins.

Responses of a novel salt-tolerant Streptomyces albidoflavus DUT_AHX capable of degrading nitrobenzene to salinity stress

A novel salt-tolerant strain DUT_AHX, which was capable of utilizing nitrobenzene (NB) as the sole carbon source, was isolated from NB-contaminated soil. Furthermore, it was identified as Streptomyces albidoflavus on the basis of physiological and biochemical tests and 16S ribosomal DNA (rDNA) sequence analysis. It can grow in the presence of NaCl up to 12% (w/v) or NB up to 900 mg/l in mineral salts basal (MSB) medium. The exogenously added osmoprotectants such as glycin, glutamic acid, proline, betaine and ectoine can improve growth of strain DUT_AHX in the presence of 10% (w/v) NaCl. NB-grown cells of strain DUT_AHX in modified MSB medium can degrade NB with the concomitant release of ammonia. Moreover, crude extracts of NB-grown strain DUT_AHX mainly contained 2-aminophenol 1,6-dioxygenase activity. These indicate that NB degradation by strain DUT_AHX might involve a partial reductive pathway. The proteins induced by salinity stress or NB were analyzed by native-gradient polyacrylamide gel electrophoresis (PAGE) and sodium dodecyl sulfate (SDS)-PAGE. In NB-induced proteins de novo, 141 kDa protein on the native-gradient PAGE gel was excised and electroeluted. Furthermore, enzyme tests exhibit the 2-aminophenol 1,6-dioxygenase activity of purified 141 kDa protein is 11-fold that of the cell-free extracts. The exploitation of strain DUT_AHX in salinity stress will be a remarkable improvement in NB bioremediation and wastewater treatment in high salinity.

Metabolism of 4-amino-3-hydroxybenzoic acid by Bordetella sp. strain 10d: A different modified meta-cleavage pathway for 2-aminophenols

Bordetella sp. strain 10d metabolizes 4-amino-3-hydroxybenzoic acid via 2-hydroxymuconic 6-semialdehyde. Cell extracts from 4-amino-3-hydroxybenzoate-grown cells showed high NAD(+)-dependent 2-hydroxymuconic 6-semialdehyde dehydrogenase, 4-oxalocrotonate tautomerase, 4-oxalocrotonate decarboxylase, and 2-oxopent-4-enoate hydratase activities, but no 2-hydroxymuconic 6-semialdehyde hydrolase activity. These enzymes involved in 4-amino-3-hydroxybenzoate metabolism were purified and characterized. When 2-hydroxymuconic 6-semialdehyde was used as substrate in a reaction mixture containing NAD(+) and cell extracts from 4-amino-3-hydroxybenzoate-grown cells, 4-oxalocrotonic acid, 2-oxopent-4-enoic acid, and 4-hydroxy-2-oxovaleric acid were identified as intermediates, and pyruvic acid was identified as the final product. A complete pathway for the metabolism of 4-amino-3-hydroxybenzoic acid in strain 10d is proposed. Strain 10d metabolized 2-hydroxymuconic 6-semialdehyde derived from 4-amino-3-hydroxybenzoic acid via a dehydrogenative route, not via a hydrolytic route. This proposed metabolic pathway differs considerably from the modified meta-cleavage pathway of 2-aminophenol and those previously reported for methyl- and chloro-derivatives.

Novel partial reductive pathway for 4-chloronitrobenzene and nitrobenzene degradation in Comamonas sp. strain CNB-1

Comamonas sp. strain CNB-1 grows on 4-chloronitrobenzene (4-CNB) and nitrobenzene as sole carbon and nitrogen sources. In this study, two genetic segments, cnbB-orf2-cnbA and cnbR-orf1-cnbCaCbDEFGHI, located on a newly isolated plasmid, pCNB1 (ca. 89 kb), and involved in 4-CNB/nitrobenzene degradation, were characterized. Seven genes (cnbA, cnbB, cnbCa, cnbCb, cnbD, cnbG, and cnbH) were cloned and functionally expressed in recombinant Escherichia coli, and they were identified as encoding 4-CNB nitroreductase (CnbA), 1-hydroxylaminobenzene mutase (CnbB), 2-aminophenol 1,6-dioxygenase (CnbCab), 2-amino-5-chloromuconic semialdehyde dehydrogenase (CnbD), 2-hydroxy-5-chloromuconic acid (2H5CM) tautomerase, and 2-amino-5-chloromuconic acid (2A5CM) deaminase (CnbH). In particular, the 2A5CM deaminase showed significant identities (31 to 38%) to subunit A of Asp-tRNAAsn/Glu-tRNAGln amidotransferase and not to the previously identified deaminases for nitroaromatic compound degradation. Genetic cloning and expression of cnbH in Escherichia coli revealed that CnbH catalyzed the conversion of 2A5CM into 2H5CM and ammonium. Four other genes (cnbR, cnbE, cnbF, and cnbI) were tentatively identified according to their high sequence identities to other functionally identified genes. It was proposed that CnbH might represent a novel type of deaminase and be involved in a novel partial reductive pathway for chloronitrobenzene or nitrobenzene degradation.

The clc element of Pseudomonas sp. strain B13, a genomic island with various catabolic properties

Pseudomonas sp. strain B13 is a bacterium known to degrade chloroaromatic compounds. The properties to use 3- and 4-chlorocatechol are determined by a self-transferable DNA element, the clc element, which normally resides at two locations in the cell's chromosome. Here we report the complete nucleotide sequence of the clc element, demonstrating the unique catabolic properties while showing its relatedness to genomic islands and integrative and conjugative elements rather than to other known catabolic plasmids. As far as catabolic functions, the clc element harbored, in addition to the genes for chlorocatechol degradation, a complete functional operon for 2-aminophenol degradation and genes for a putative aromatic compound transport protein and for a multicomponent aromatic ring dioxygenase similar to anthranilate hydroxylase. The genes for catabolic functions were inducible under various conditions, suggesting a network of catabolic pathway induction. For about half of the open reading frames (ORFs) on the clc element, no clear functional prediction could be given, although some indications were found for functions that were similar to plasmid conjugation. The region in which these ORFs were situated displayed a high overall conservation of nucleotide sequence and gene order to genomic regions in other recently completed bacterial genomes or to other genomic islands. Most notably, except for two discrete regions, the clc element was almost 100% identical over the whole length to a chromosomal region in Burkholderia xenovorans LB400. This indicates the dynamic evolution of this type of element and the continued transition between elements with a more pathogenic character and those with catabolic properties.

The naphthalene catabolic (nag) genes of Polaromonas naphthalenivorans CJ2: evolutionary implications for two gene clusters and novel regulatory control

Polaromonas naphthalenivorans CJ2, found to be responsible for the degradation of naphthalene in situ at a coal tar waste-contaminated site (C.-O. Jeon et al., Proc. Natl. Acad. Sci. USA 100:13591-13596, 2003), is able to grow on mineral salts agar media with naphthalene as the sole carbon source. Beginning from a 484-bp nagAc-like region, we used a genome walking strategy to sequence genes encoding the entire naphthalene degradation pathway andadditional flanking regions. We found that the naphthalene catabolic genes in P. naphthalenivorans CJ2 were divided into one large and one small gene cluster, separated by an unknown distance. The large gene cluster (nagRAaGHAbAcAdBFCQEDJI'ORF1tnpA) is bounded by a LysR-type regulator (nagR). The small cluster (nagR2ORF2I"KL) is bounded by a MarR-type regulator (nagR2). The catabolic genes of P. naphthalenivorans CJ2 were homologous to many of those of Ralstonia U2, which uses the gentisate pathway to convert naphthalene to central metabolites. However, three open reading frames (nagY, nagM, and nagN), present in Ralstonia U2, were absent. Also, P. naphthalenivorans carries two copies of gentisate dioxygenase (nagI) with 77.4% DNA sequence identity to one another and 82% amino acid identity to their homologue in Ralstonia sp. strain U2. Investigation of the operons using reverse transcription PCR showed that each cluster was controlled independently by its respective promoter. Insertional inactivation and lacZ reporter assays showed that nagR2 is a negative regulator and that expression of the small cluster is not induced by naphthalene, salicylate, or gentisate. Association of two putative Azoarcus-related transposases with the large cluster and one Azoarcus-related putative salicylate 5-hydroxylase gene (ORF2) in the small cluster suggests that mobile genetic elements were likely involved in creating the novel arrangement of catabolic and regulatory genes in P. naphthalenivorans.

Bacterial transcriptional regulators for degradation pathways of aromatic compounds

Human activities have resulted in the release and introduction into the environment of a plethora of aromatic chemicals. The interest in discovering how bacteria are dealing with hazardous environmental pollutants has driven a large research community and has resulted in important biochemical, genetic, and physiological knowledge about the degradation capacities of microorganisms and their application in bioremediation, green chemistry, or production of pharmacy synthons. In addition, regulation of catabolic pathway expression has attracted the interest of numerous different groups, and several catabolic pathway regulators have been exemplary for understanding transcription control mechanisms. More recently, information about regulatory systems has been used to construct whole-cell living bioreporters that are used to measure the quality of the aqueous, soil, and air environment. The topic of biodegradation is relatively coherent, and this review presents a coherent overview of the regulatory systems involved in the transcriptional control of catabolic pathways. This review summarizes the different regulatory systems involved in biodegradation pathways of aromatic compounds linking them to other known protein families. Specific attention has been paid to describing the genetic organization of the regulatory genes, promoters, and target operon(s) and to discussing present knowledge about signaling molecules, DNA binding properties, and operator characteristics, and evidence from regulatory mutants. For each regulator family, this information is combined with recently obtained protein structural information to arrive at a possible mechanism of transcription activation. This demonstrates the diversity of control mechanisms existing in catabolic pathways.

A novel coupled enzyme assay reveals an enzyme responsible for the deamination of a chemically unstable intermediate in the metabolic pathway of 4-amino-3-hydroxybenzoic acid in Bordetella sp. strain 10d

2-amino-5-carboxymuconic 6-semialdehyde is an unstable intermediate in the meta-cleavage pathway of 4-amino-3-hydroxybenzoic acid in Bordetella sp. strain 10d. In vitro, this compound is nonenzymatically converted to 2,5-pyridinedicarboxylic acid. Crude extracts of strain 10d grown on 4-amino-3-hydroxybenzoic acid converted 2-amino-5-carboxymuconic 6-semialdehyde formed from 4-amino-3-hydroxybenzoic acid by the first enzyme in the pathway, 4-amino-3-hydroxybenzoate 2,3-dioxygenase, to a yellow compound (epsilonmax = 375 nm). The enzyme in the crude extract carrying out the next step was purified to homogeneity. The yellow compound formed from 4-amino-3-hydroxybenzoic acid by this purified enzyme and purified 4-amino-3-hydroxybenzoate 2,3-dioxygenase in a coupled assay was identified as 2-hydroxymuconic 6-semialdehyde by GC-MS analysis. A mechanism for the formation of 2-hydroxymuconic 6-semialdehyde via enzymatic deamination and nonenzymatic decarboxylation is proposed based on results of spectrophotometric analyses. The purified enzyme, designated 2-amino-5-carboxymuconic 6-semialdehyde deaminase, is a new type of deaminase that differs from the 2-aminomuconate deaminases reported previously in that it primarily and specifically attacks 2-amino-5-carboxymuconic 6-semialdehyde. The deamination step in the proposed pathway differs from that in the pathways for 2-aminophenol and its derivatives.

Expression of the nitroarene dioxygenase genes in Comamonas sp. strain JS765 and Acidovorax sp. strain JS42 is induced by multiple aromatic compounds

This work reports a genetic analysis of the expression of nitrobenzene dioxygenase (NBDO) in Comamonas sp. strain JS765 and 2-nitrotoluene dioxygenase (2NTDO) in Acidovorax sp. strain JS42. Strains JS765 and JS42 possess identical LysR-type regulatory proteins, NbzR and NtdR, respectively. NbzR/NtdR is homologous to NahR, the positive salicylate-responsive transcriptional activator of the naphthalene degradation genes in Pseudomonas putida G7. The genes encoding NBDO and 2NTDO in each strain are cotranscribed, and transcription starts at the same site within identical promoter regions for each operon. Results from a lacZ reporter gene fusion demonstrated that expression of NBDO and 2NTDO is induced by multiple aromatic compounds, including an array of nitroaromatic compounds (nitrobenzene, 2-, 3-, and 4-nitrotoluene, 2,4- and 2,6-dinitrotoluene, and aminodinitrotoluenes), as well as salicylate and anthranilate. The nitroaromatic compounds appear to be the actual effector molecules. Analysis of beta-galactosidase and 2NTDO activities with strain JS42 demonstrated that NtdR was required for induction by all of the inducing compounds, high basal-level expression of 2NTDO, and complementation of a JS42 ntdR null mutant. Complementation with the closely related regulators NagR (from Ralstonia sp. strain U2) and NahR restored only induction by the archetype inducers, salicylate or salicylate and anthranilate, respectively, and did not restore the high basal level of expression of 2NTDO. The mechanism of 2NTDO gene regulation in JS42, and presumably that of NBDO gene regulation in JS765, appear similar to that of NahR-regulated genes in Pseudomonas putida G7. However, NbzR and NtdR appear to have evolved a broader specificity in JS42 and JS765, allowing for recognition of nitroaromatic compounds while retaining the ability to respond to salicylate and anthranilate. NtdR is also the first example of a nitroarene-responsive LysR-type transcriptional activator.

Prokaryotic homologs of the eukaryotic 3-hydroxyanthranilate 3,4-dioxygenase and 2-amino-3-carboxymuconate-6-semialdehyde decarboxylase in the 2-nitrobenzoate degradation pathway of Pseudomonas fluorescens strain KU-7

The 2-nitrobenzoic acid degradation pathway of Pseudomonas fluorescens strain KU-7 proceeds via a novel 3-hydroxyanthranilate intermediate. In this study, we cloned and sequenced a 19-kb DNA locus of strain KU-7 that encompasses the 3-hydroxyanthranilate meta-cleavage pathway genes. The gene cluster, designated nbaEXHJIGFCDR, is organized tightly and in the same direction. The nbaC and nbaD gene products were found to be novel homologs of the eukaryotic 3-hydroxyanthranilate 3,4-dioxygenase and 2-amino-3-carboxymuconate-6-semialdehyde decarboxylase, respectively. The NbaC enzyme carries out the oxidation of 3-hydroxyanthranilate to 2-amino-3-carboxymuconate-6-semialdehyde, while the NbaD enzyme catalyzes the decarboxylation of the latter compound to 2-aminomuconate-6-semialdehyde. The NbaC and NbaD proteins were overexpressed in Escherichia coli and characterized. The substrate specificity of the 23.8-kDa NbaC protein was found to be restricted to 3-hydroxyanthranilate. In E. coli, this enzyme oxidizes 3-hydroxyanthranilate with a specific activity of 8 U/mg of protein. Site-directed mutagenesis experiments revealed the essential role of two conserved histidine residues (His52 and His96) in the NbaC sequence. The NbaC activity is also dependent on the presence of Fe(2+) but is inhibited by other metal ions, such as Zn(2+), Cu(2+), and Cd(2+). The NbaD protein was overproduced as a 38.7-kDa protein, and its specific activity towards 2-amino-3-carboxymuconate-6-semialdehyde was 195 U/mg of protein. Further processing of 2-aminomuconate-6-semialdehyde to pyruvic acid and acetyl coenzyme A was predicted to proceed via the activities of NbaE, NbaF, NbaG, NbaH, NbaI, and NbaJ. The predicted amino acid sequences of these proteins are highly homologous to those of the corresponding proteins involved in the metabolism of 2-aminophenol (e.g., AmnCDEFGH in Pseudomonas sp. strain AP-3). The NbaR-encoding gene is predicted to have a regulatory function of the LysR family type. The function of the product of the small open reading frame, NbaX, like the homologous sequences in the nitrobenzene or 2-aminophenol metabolic pathway, remains elusive.

Computational identification of operons in microbial genomes

By applying graph representations to biochemical pathways, a new computational pipeline is proposed to find potential operons in microbial genomes. The algorithm relies on the fact that enzyme genes in operons tend to catalyze successive reactions in metabolic pathways. We applied this algorithm to 42 microbial genomes to identify putative operon structures. The predicted operons from Escherichia coli were compared with a selected metabolism-related operon dataset from the RegulonDB database, yielding a prediction sensitivity (89%) and specificity (87%) relative to this dataset. Several examples of detected operons are given and analyzed. Modular gene cluster transfer and operon fusion are observed. A further use of predicted operon data to assign function to putative genes was suggested and, as an example, a previous putative gene (MJ1604) from Methanococcus jannaschii is now annotated as a phosphofructokinase, which was regarded previously as a missing enzyme in this organism. GC content changes in the operon region and nonoperon region were examined. The results reveal a clear GC content transition at the boundaries of putative operons. We looked further into the conservation of operons across genomes. A trp operon alignment is analyzed in depth to show gene loss and rearrangement in different organisms during operon evolution.