Thomas D. Niehaus

MINEs: open access databases of computationally predicted enzyme promiscuity products for untargeted metabolomics

BackgroundIn spite of its great promise, metabolomics has proven difficult to execute in an untargeted and generalizable manner. Liquid chromatography–mass spectrometry (LC–MS) has made it possible to gather data on thousands of cellular metabolites. However, matching metabolites to their spectral features continues to be a bottleneck, meaning that much of the collected information remains uninterpreted and that new metabolites are seldom discovered in untargeted studies. These challenges require new approaches that consider compounds beyond those available in curated biochemistry databases.DescriptionHere we present Metabolic In silico Network Expansions (MINEs), an extension of known metabolite databases to include molecules that have not been observed, but are likely to occur based on known metabolites and common biochemical reactions. We utilize an algorithm called the Biochemical Network Integrated Computational Explorer (BNICE) and expert-curated reaction rules based on the Enzyme Commission classification system to propose the novel chemical structures and reactions that comprise MINE databases. Starting from the Kyoto Encyclopedia of Genes and Genomes (KEGG) COMPOUND database, the MINE contains over 571,000 compounds, of which 93% are not present in the PubChem database. However, these MINE compounds have on average higher structural similarity to natural products than compounds from KEGG or PubChem. MINE databases were able to propose annotations for 98.6% of a set of 667 MassBank spectra, 14% more than KEGG alone and equivalent to PubChem while returning far fewer candidates per spectra than PubChem (46 vs. 1715 median candidates). Application of MINEs to LC–MS accurate mass data enabled the identity of an unknown peak to be confidently predicted.ConclusionsMINE databases are freely accessible for non-commercial use via user-friendly web-tools at http://minedatabase.mcs.anl.gov and developer-friendly APIs. MINEs improve metabolomics peak identification as compared to general chemical databases whose results include irrelevant synthetic compounds. Furthermore, MINEs complement and expand on previous in silico generated compound databases that focus on human metabolism. We are actively developing the database; future versions of this resource will incorporate transformation rules for spontaneous chemical reactions and more advanced filtering and prioritization of candidate structures.

High-throughput comparison, functional annotation, and metabolic modeling of plant genomes using the PlantSEED resource.

Significance Genes must be annotated with their correct functions if genome data are to support hypothesis building and metabolic engineering. PlantSEED was developed to streamline the process of annotating plant genome sequences, to construct metabolic models based on genome annotations automatically, and to use models to test the annotation of these sequences, allowing the detection of gaps and errors in gene annotations and the prediction of new functions. PlantSEED is designed to grow in an iterative manner by including new plant genome sequences, new annotations harvested from the literature, and improved biochemical data, all of which are integrated in a consistent manner into the PlantSEED genomes and metabolic models. The increasing number of sequenced plant genomes is placing new demands on the methods applied to analyze, annotate, and model these genomes. Today’s annotation pipelines result in inconsistent gene assignments that complicate comparative analyses and prevent efficient construction of metabolic models. To overcome these problems, we have developed the PlantSEED, an integrated, metabolism-centric database to support subsystems-based annotation and metabolic model reconstruction for plant genomes. PlantSEED combines SEED subsystems technology, first developed for microbial genomes, with refined protein families and biochemical data to assign fully consistent functional annotations to orthologous genes, particularly those encoding primary metabolic pathways. Seamless integration with its parent, the prokaryotic SEED database, makes PlantSEED a unique environment for cross-kingdom comparative analysis of plant and bacterial genomes. The consistent annotations imposed by PlantSEED permit rapid reconstruction and modeling of primary metabolism for all plant genomes in the database. This feature opens the unique possibility of model-based assessment of the completeness and accuracy of gene annotation and thus allows computational identification of genes and pathways that are restricted to certain genomes or need better curation. We demonstrate the PlantSEED system by producing consistent annotations for 10 reference genomes. We also produce a functioning metabolic model for each genome, gapfilling to identify missing annotations and proposing gene candidates for missing annotations. Models are built around an extended biomass composition representing the most comprehensive published to date. To our knowledge, our models are the first to be published for seven of the genomes analyzed.

Arabidopsis and Maize RidA Proteins Preempt Reactive Enamine/Imine Damage to Branched-Chain Amino Acid Biosynthesis in Plastids

Plant RidA proteins protect an enzyme of branched-chain amino acid biosynthesis from inactivation by hydrolyzing reactive pathway intermediates before they can damage the enzyme. RidA proteins are thus crucial for the efficient biosynthesis of branched-chain amino acids in plants and provide an iconic example of the preemption of metabolite damage. RidA (for Reactive Intermediate Deaminase A) proteins are ubiquitous, yet their function in eukaryotes is unclear. It is known that deleting Salmonella enterica ridA causes Ser sensitivity and that S. enterica RidA and its homologs from other organisms hydrolyze the enamine/imine intermediates that Thr dehydratase forms from Ser or Thr. In S. enterica, the Ser-derived enamine/imine inactivates a branched-chain aminotransferase; RidA prevents this damage. Arabidopsis thaliana and maize (Zea mays) have a RidA homolog that is predicted to be plastidial. Expression of either homolog complemented the Ser sensitivity of the S. enterica ridA mutant. The purified proteins hydrolyzed the enamines/imines formed by Thr dehydratase from Ser or Thr and protected the Arabidopsis plastidial branched-chain aminotransferase BCAT3 from inactivation by the Ser-derived enamine/imine. In vitro chloroplast import assays and in vivo localization of green fluorescent protein fusions showed that Arabidopsis RidA and Thr dehydratase are chloroplast targeted. Disrupting Arabidopsis RidA reduced root growth and raised the root and shoot levels of the branched-chain amino acid biosynthesis intermediate 2-oxobutanoate; Ser treatment exacerbated these effects in roots. Supplying Ile reversed the root growth defect. These results indicate that plastidial RidA proteins can preempt damage to BCAT3 and Ile biosynthesis by hydrolyzing the Ser-derived enamine/imine product of Thr dehydratase.

Proteins of Unknown Biochemical Function: A Persistent Problem and a Roadmap to Help Overcome It

The functions of over 60% of plant enzymes are unknown; nearly one-half of these proteins have microbial homologs, enabling the use of cross-kingdom comparative genomics to predict their functions. The number of sequenced genomes is rapidly increasing, but functional annotation of the genes in these genomes lags far behind. Even in Arabidopsis (Arabidopsis thaliana), only approximately 40% of enzyme- and transporter-encoding genes have credible functional annotations, and this number is even lower in nonmodel plants. Functional characterization of unknown genes is a challenge, but various databases (e.g. for protein localization and coexpression) can be mined to provide clues. If homologous microbial genes exist—and about one-half the genes encoding unknown enzymes and transporters in Arabidopsis have microbial homologs—cross-kingdom comparative genomics can powerfully complement plant-based data. Multiple lines of evidence can strengthen predictions and warrant experimental characterization. In some cases, relatively quick tests in genetically tractable microbes can determine whether a prediction merits biochemical validation, which is costly and demands specialized skills.

Arabidopsis TH2 Encodes the Orphan Enzyme Thiamin Monophosphate Phosphatase

Genomic, transcriptomic, genetic, and biochemical evidence implicates Arabidopsis TH2 in dephosphorylation of thiamin monophosphate, a key step in biosynthesis of thiamin diphosphate, an essential cofactor for central metabolic enzymes. To synthesize the cofactor thiamin diphosphate (ThDP), plants must first hydrolyze thiamin monophosphate (ThMP) to thiamin, but dedicated enzymes for this hydrolysis step were unknown and widely doubted to exist. The classical thiamin-requiring th2-1 mutation in Arabidopsis thaliana was shown to reduce ThDP levels by half and to increase ThMP levels 5-fold, implying that the THIAMIN REQUIRING2 (TH2) gene product could be a dedicated ThMP phosphatase. Genomic and transcriptomic data indicated that TH2 corresponds to At5g32470, encoding a HAD (haloacid dehalogenase) family phosphatase fused to a TenA (thiamin salvage) family protein. Like the th2-1 mutant, an insertional mutant of At5g32470 accumulated ThMP, and the thiamin requirement of the th2-1 mutant was complemented by wild-type At5g32470. Complementation tests in Escherichia coli and enzyme assays with recombinant proteins confirmed that At5g32470 and its maize (Zea mays) orthologs GRMZM2G148896 and GRMZM2G078283 are ThMP-selective phosphatases whose activity resides in the HAD domain and that the At5g32470 TenA domain has the expected thiamin salvage activity. In vitro and in vivo experiments showed that alternative translation start sites direct the At5g32470 protein to the cytosol and potentially also to mitochondria. Our findings establish that plants have a dedicated ThMP phosphatase and indicate that modest (50%) ThDP depletion can produce severe deficiency symptoms.

Discovery of a widespread prokaryotic 5-oxoprolinase that was hiding in plain sight

5-Oxoproline (OP) is well-known as an enzymatic intermediate in the eukaryotic γ-glutamyl cycle, but it is also an unavoidable damage product formed spontaneously from glutamine and other sources. Eukaryotes metabolize OP via an ATP-dependent 5-oxoprolinase; most prokaryotes lack homologs of this enzyme (and the γ-glutamyl cycle) but are predicted to have some way to dispose of OP if its spontaneous formation in vivo is significant. Comparative analysis of prokaryotic genomes showed that the gene encoding pyroglutamyl peptidase, which removes N-terminal OP residues, clusters in diverse genomes with genes specifying homologs of a fungal lactamase (renamed prokaryotic 5-oxoprolinase A, pxpA) and homologs of allophanate hydrolase subunits (renamed pxpB and pxpC). Inactivation of Bacillus subtilis pxpA, pxpB, or pxpC genes slowed growth, caused OP accumulation in cells and medium, and prevented use of OP as a nitrogen source. Assays of cell lysates showed that ATP-dependent 5-oxoprolinase activity disappeared when pxpA, pxpB, or pxpC was inactivated. 5-Oxoprolinase activity could be reconstituted in vitro by mixing recombinant B. subtilis PxpA, PxpB, and PxpC proteins. In addition, overexpressing Escherichia coli pxpABC genes in E. coli increased 5-oxoprolinase activity in lysates ≥1700-fold. This work shows that OP is a major universal metabolite damage product and that OP disposal systems are common in all domains of life. Furthermore, it illustrates how easily metabolite damage and damage-control systems can be overlooked, even for central metabolites in model organisms.

Draft Nuclear Genome Sequence of the Liquid Hydrocarbon–Accumulating Green Microalga Botryococcus braunii Race B (Showa)

ABSTRACT Botryococcus braunii has long been known as a prodigious producer of liquid hydrocarbon oils that can be converted into combustion engine fuels. This draft genome for the B race of B. braunii will allow researchers to unravel important hydrocarbon biosynthetic pathways and identify possible regulatory networks controlling this unusual metabolism.

Mapping a kingdom-specific functional domain of squalene synthase.

Squalene synthase catalyzes the first committed step in sterol biosynthesis and consists of both an amino-terminal catalytic domain and a carboxy-terminal domain tethering the enzyme to the ER membrane. While the overall architecture of this enzyme is identical in eukaryotes, it was previously shown that plant and animal genes cannot complement a squalene synthase knockout mutation in yeast unless the carboxy-terminal domain is swapped for one of fungal origin. This implied a unique component of the fungal carboxy-terminal domain was responsible for the complementation phenotype. To identify this motif, we used Saccharomyces cerevisiae with a squalene synthase knockout mutation, and expressed intact and chimeric squalene synthases originating from fungi, plants, and animals. In contrast to previous observations, all enzymes tested could partially complement the knockout mutation when the genes were weakly expressed. However, when highly expressed, non-fungal squalene synthases could not complement the yeast mutation and instead led to the accumulation of a toxic intermediate(s) as defined by mutations of genes downstream in the ergosterol pathway. Restoration of the complete complementation phenotype was mapped to a 26-amino acid hinge region linking the catalytic and membrane-spanning domains specific to fungal squalene synthases. Over-expression of the C-terminal domain containing a hinge domain from fungi, not from animals or plants, led to growth inhibition of wild-type yeast. Because this hinge region is unique to and highly conserved within each kingdom of life, the data suggests that the hinge domain plays an essential functional role, such as assembly of ergosterol multi-enzyme complexes in fungi.

Carboxythiazole is a key microbial nutrient currency and critical component of thiamin biosynthesis

Almost all cells require thiamin, vitamin B1 (B1), which is synthesized via the coupling of thiazole and pyrimidine precursors. Here we demonstrate that 5-(2-hydroxyethyl)-4-methyl-1,3-thiazole-2-carboxylic acid (cHET) is a useful in vivo B1 precursor for representatives of ubiquitous marine picoeukaryotic phytoplankton and Escherichia coli – drawing attention to cHET as a valuable exogenous micronutrient for microorganisms with ecological, industrial, and biomedical value. Comparative utilization experiments with the terrestrial plant Arabidopsis thaliana revealed that it can also use exogenous cHET, but notably, picoeukaryotic marine phytoplankton and E. coli were adapted to grow on low (picomolar) concentrations of exogenous cHET. Our results call for the modification of the conventional B1 biosynthesis model to incorporate cHET as a key precursor for B1 biosynthesis in two domains of life, and for consideration of cHET as a microbial micronutrient currency modulating marine primary productivity and community interactions in human gut-hosted microbiomes.

Identification of a metabolic disposal route for the oncometabolite S-(2-succino)cysteine in Bacillus subtilis

Cellular thiols such as cysteine spontaneously and readily react with the respiratory intermediate fumarate, resulting in the formation of stable S-(2-succino)-adducts. Fumarate-mediated succination of thiols increases in certain tumors and in response to glucotoxicity associated with diabetes. Therefore, S-(2-succino)-adducts such as S-(2-succino)cysteine (2SC) are considered oncometabolites and biomarkers for human disease. No disposal routes for S-(2-succino)-compounds have been reported prior to this study. Here, we show that Bacillus subtilis metabolizes 2SC to cysteine using a pathway encoded by the yxe operon. The first step is N-acetylation of 2SC followed by an oxygenation that we propose results in the release of oxaloacetate and N-acetylcysteine, which is deacetylated to give cysteine. Knockouts of the genes predicted to mediate each step in the pathway lose the ability to grow on 2SC as the sulfur source and accumulate the expected upstream metabolite(s). We further show that N-acetylation of 2SC relieves toxicity. This is the first demonstration of a metabolic disposal route for any S-(2-succino)-compound, paving the way toward the identification of corresponding pathways in other species.