'Nothing of chemistry disappears in biology': the Top 30 damage-prone endogenous metabolites.

Many common metabolites are intrinsically unstable and reactive, and hence prone to chemical (i.e. non-enzymatic) damage in vivo Although this fact is widely recognized, the purely chemical side-reactions of metabolic intermediates can be surprisingly hard to track down in the literature and are often treated in an unprioritized case-by-case way. Moreover, spontaneous chemical side-reactions tend to be overshadowed today by side-reactions mediated by promiscuous ('sloppy') enzymes even though chemical damage to metabolites may be even more prevalent than damage from enzyme sloppiness, has similar outcomes, and is held in check by similar biochemical repair or pre-emption mechanisms. To address these limitations and imbalances, here we draw together and systematically integrate information from the (bio)chemical literature, from cheminformatics, and from genome-scale metabolic models to objectively define a 'Top 30' list of damage-prone metabolites. A foundational part of this process was to derive general reaction rules for the damage chemistries involved. The criteria for a 'Top 30' metabolite included predicted chemical reactivity, essentiality, and occurrence in diverse organisms. We also explain how the damage chemistry reaction rules ('operators') are implemented in the Chemical-Damage-MINE (CD-MINE) database (minedatabase.mcs.anl.gov/#/top30) to provide a predictive tool for many additional potential metabolite damage products. Lastly, we illustrate how defining a 'Top 30' list can drive genomics-enabled discovery of the enzymes of previously unrecognized damage-control systems, and how applying chemical damage reaction rules can help identify previously unknown peaks in metabolomics profiles.

Systematic identification and analysis of frequent gene fusion events in metabolic pathways.

BACKGROUND: Gene fusions are the most powerful type of in silico-derived functional associations. However, many fusion compilations were made when <100 genomes were available, and algorithms for identifying fusions need updating to handle the current avalanche of sequenced genomes. The availability of a large fusion dataset would help probe functional associations and enable systematic analysis of where and why fusion events occur. RESULTS: Here we present a systematic analysis of fusions in prokaryotes. We manually generated two training sets: (i) 121 fusions in the model organism Escherichia coli; (ii) 131 fusions found in B vitamin metabolism. These sets were used to develop a fusion prediction algorithm that captured the training set fusions with only 7 % false negatives and 50 % false positives, a substantial improvement over existing approaches. This algorithm was then applied to identify 3.8 million potential fusions across 11,473 genomes. The results of the analysis are available in a searchable database at http://modelseed.org/projects/fusions/ . A functional analysis identified 3,000 reactions associated with frequent fusion events and revealed areas of metabolism where fusions are particularly prevalent. CONCLUSIONS: Customary definitions of fusions were shown to be ambiguous, and a stricter one was proposed. Exploring the genes participating in fusion events showed that they most commonly encode transporters, regulators, and metabolic enzymes. The major rationales for fusions between metabolic genes appear to be overcoming pathway bottlenecks, avoiding toxicity, controlling competing pathways, and facilitating expression and assembly of protein complexes. Finally, our fusion dataset provides powerful clues to decipher the biological activities of domains of unknown function.

Experimental and Metabolic Modeling Evidence for a Folate-Cleaving Side-Activity of Ketopantoate Hydroxymethyltransferase (PanB).

Tetrahydrofolate (THF) and its one-carbon derivatives, collectively termed folates, are essential cofactors, but are inherently unstable. While it is clear that chemical oxidation can cleave folates or damage their pterin precursors, very little is known about enzymatic damage to these molecules or about whether the folate biosynthesis pathway responds adaptively to damage to its end-products. The presence of a duplication of the gene encoding the folate biosynthesis enzyme 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (FolK) in many sequenced bacterial genomes combined with a strong chromosomal clustering of the folK gene with panB, encoding the 5,10-methylene-THF-dependent enzyme ketopantoate hydroxymethyltransferase, led us to infer that PanB has a side activity that cleaves 5,10-methylene-THF, yielding a pterin product that is recycled by FolK. Genetic and metabolic analyses of Escherichia coli strains showed that overexpression of PanB leads to accumulation of the likely folate cleavage product 6-hydroxymethylpterin and other pterins in cells and medium, and-unexpectedly-to a 46% increase in total folate content. In silico modeling of the folate biosynthesis pathway showed that these observations are consistent with the in vivo cleavage of 5,10-methylene-THF by a side-activity of PanB, with FolK-mediated recycling of the pterin cleavage product, and with regulation of folate biosynthesis by folates or their damage products.

Methionine biosynthesis is essential for infection in the rice blast fungus Magnaporthe oryzae.

Methionine is a sulfur amino acid standing at the crossroads of several biosynthetic pathways. In fungi, the last step of methionine biosynthesis is catalyzed by a cobalamine-independent methionine synthase (Met6, EC 2.1.1.14). In the present work, we studied the role of Met6 in the infection process of the rice blast fungus, Magnaporthe oryzae. To this end MET6 null mutants were obtained by targeted gene replacement. On minimum medium, MET6 null mutants were auxotrophic for methionine. Even when grown in presence of excess methionine, these mutants displayed developmental defects, such as reduced mycelium pigmentation, aerial hypha formation and sporulation. They also displayed characteristic metabolic signatures such as increased levels of cysteine, cystathionine, homocysteine, S-adenosylmethionine, S-adenosylhomocysteine while methionine and glutathione levels remained unchanged. These metabolic perturbations were associated with the over-expression of MgCBS1 involved in the reversed transsulfuration pathway that metabolizes homocysteine into cysteine and MgSAM1 and MgSAHH1 involved in the methyl cycle. This suggests a physiological adaptation of M. oryzae to metabolic defects induced by the loss of Met6, in particular an increase in homocysteine levels. Pathogenicity assays showed that MET6 null mutants were non-pathogenic on both barley and rice leaves. These mutants were defective in appressorium-mediated penetration and invasive infectious growth. These pathogenicity defects were rescued by addition of exogenous methionine and S-methylmethionine. These results show that M. oryzae cannot assimilate sufficient methionine from plant tissues and must synthesize this amino acid de novo to fulfill its sulfur amino acid requirement during infection.

Improved evidence-based genome-scale metabolic models for maize leaf, embryo, and endosperm.

There is a growing demand for genome-scale metabolic reconstructions for plants, fueled by the need to understand the metabolic basis of crop yield and by progress in genome and transcriptome sequencing. Methods are also required to enable the interpretation of plant transcriptome data to study how cellular metabolic activity varies under different growth conditions or even within different organs, tissues, and developmental stages. Such methods depend extensively on the accuracy with which genes have been mapped to the biochemical reactions in the plant metabolic pathways. Errors in these mappings lead to metabolic reconstructions with an inflated number of reactions and possible generation of unreliable metabolic phenotype predictions. Here we introduce a new evidence-based genome-scale metabolic reconstruction of maize, with significant improvements in the quality of the gene-reaction associations included within our model. We also present a new approach for applying our model to predict active metabolic genes based on transcriptome data. This method includes a minimal set of reactions associated with low expression genes to enable activity of a maximum number of reactions associated with high expression genes. We apply this method to construct an organ-specific model for the maize leaf, and tissue specific models for maize embryo and endosperm cells. We validate our models using fluxomics data for the endosperm and embryo, demonstrating an improved capacity of our models to fit the available fluxomics data. All models are publicly available via the DOE Systems Biology Knowledgebase and PlantSEED, and our new method is generally applicable for analysis transcript profiles from any plant, paving the way for further in silico studies with a wide variety of plant genomes.

A directed-overflow and damage-control N-glycosidase in riboflavin biosynthesis.

Plants and bacteria synthesize the essential human micronutrient riboflavin (vitamin B2) via the same multi-step pathway. The early intermediates of this pathway are notoriously reactive and may be overproduced in vivo because riboflavin biosynthesis enzymes lack feedback controls. In the present paper, we demonstrate disposal of riboflavin intermediates by COG3236 (DUF1768), a protein of previously unknown function that is fused to two different riboflavin pathway enzymes in plants and bacteria (RIBR and RibA respectively). We present cheminformatic, biochemical, genetic and genomic evidence to show that: (i) plant and bacterial COG3236 proteins cleave the N-glycosidic bond of the first two intermediates of riboflavin biosynthesis, yielding relatively innocuous products; (ii) certain COG3236 proteins are in a multi-enzyme riboflavin biosynthesis complex that gives them privileged access to riboflavin intermediates; and (iii) COG3236 action in Arabidopsis thaliana and Escherichia coli helps maintain flavin levels. COG3236 proteins thus illustrate two emerging principles in chemical biology: directed overflow metabolism, in which excess flux is diverted out of a pathway, and the pre-emption of damage from reactive metabolites.

Evidence that glutamine transaminase and omega-amidase potentially act in tandem to close the methionine salvage cycle in bacteria and plants.

S-Adenosylmethionine is converted enzymatically and non-enzymatically to methylthioadenosine, which is recycled to methionine (Met) via a salvage pathway. In plants and bacteria, enzymes for all steps in this pathway are known except the last: transamination of α-ketomethylthiobutyrate to give Met. In mammals, glutamine transaminase K (GTK) and ω-amidase (ω-Am) are thought to act in tandem to execute this step, with GTK forming α-ketoglutaramate, which ω-Am hydrolyzes. Comparative genomics indicated that GTK and ω-Am could function likewise in plants and bacteria because genes encoding GTK and ω-Am homologs (i) co-express with the Met salvage gene 5-methylthioribose kinase in Arabidopsis, and (ii) cluster on the chromosome with each other and with Met salvage genes in diverse bacteria. Consistent with this possibility, tomato, maize, and Bacillus subtilis GTK and ω-Am homologs had the predicted activities: GTK was specific for glutamine as amino donor and strongly preferred α-ketomethylthiobutyrate as amino acceptor, and ω-Am strongly preferred α-ketoglutaramate. Also consistent with this possibility, plant GTK and ω-Am were localized to the cytosol, where the Met salvage pathway resides, as well as to organelles. This multiple targeting was shown to result from use of alternative start codons. In B. subtilis, ablating GTK or ω-Am had a modest but significant inhibitory effect on growth on 5-methylthioribose as sole sulfur source. Collectively, these data indicate that while GTK, coupled with ω-Am, is positioned to support significant Met salvage flux in plants and bacteria, it can probably be replaced by other aminotransferases.

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

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.

A cross-kingdom Nudix enzyme that pre-empts damage in thiamin metabolism.

Genes specifying the thiamin monophosphate phosphatase and adenylated thiazole diphosphatase steps in fungal and plant thiamin biosynthesis remain unknown, as do genes for ThDP (thiamin diphosphate) hydrolysis in thiamin metabolism. A distinctive Nudix domain fused to Tnr3 (thiamin diphosphokinase) in Schizosaccharomyces pombe was evaluated as a candidate for these functions. Comparative genomic analysis predicted a role in thiamin metabolism, not biosynthesis, because free-standing homologues of this Nudix domain occur not only in fungi and plants, but also in proteobacteria (whose thiamin biosynthesis pathway has no adenylated thiazole or thiamin monophosphate hydrolysis steps) and animals (which do not make thiamin). Supporting this prediction, recombinant Tnr3 and its Saccharomyces cerevisiae, Arabidopsis and maize Nudix homologues lacked thiamin monophosphate phosphatase activity, but were active against ThDP, and up to 60-fold more active against diphosphates of the toxic thiamin degradation products oxy- and oxo-thiamin. Deleting the S. cerevisiae Nudix gene (YJR142W) lowered oxythiamin resistance, overexpressing it raised resistance, and expressing its plant or bacterial counterparts restored resistance to the YJR142W deletant. By converting the diphosphates of damaged forms of thiamin into monophosphates, the Tnr3 Nudix domain and its homologues can pre-empt the misincorporation of damaged diphosphates into ThDP-dependent enzymes, and the resulting toxicity.

Identification of mitochondrial coenzyme a transporters from maize and Arabidopsis.

Plants make coenzyme A (CoA) in the cytoplasm but use it for reactions in mitochondria, chloroplasts, and peroxisomes, implying that these organelles have CoA transporters. A plant peroxisomal CoA transporter is already known, but plant mitochondrial or chloroplastic CoA transporters are not. Mitochondrial CoA transporters belonging to the mitochondrial carrier family, however, have been identified in yeast (Saccharomyces cerevisiae; Leu-5p) and mammals (SLC25A42). Comparative genomic analysis indicated that angiosperms have two distinct homologs of these mitochondrial CoA transporters, whereas nonflowering plants have only one. The homologs from maize (Zea mays; GRMZM2G161299 and GRMZM2G420119) and Arabidopsis (Arabidopsis thaliana; At1g14560 and At4g26180) all complemented the growth defect of the yeast leu5Δ mitochondrial CoA carrier mutant and substantially restored its mitochondrial CoA level, confirming that these proteins have CoA transport activity. Dual-import assays with purified pea (Pisum sativum) mitochondria and chloroplasts, and subcellular localization of green fluorescent protein fusions in transiently transformed tobacco (Nicotiana tabacum) Bright Yellow-2 cells, showed that the maize and Arabidopsis proteins are targeted to mitochondria. Consistent with the ubiquitous importance of CoA, the maize and Arabidopsis mitochondrial CoA transporter genes are expressed at similar levels throughout the plant. These data show that representatives of both monocotyledons and eudicotyledons have twin, mitochondrially located mitochondrial carrier family carriers for CoA. The highly conserved nature of these carriers makes possible their reliable annotation in other angiosperm genomes.

Identification and characterization of the missing pyrimidine reductase in the plant riboflavin biosynthesis pathway.

Riboflavin (vitamin B2) is the precursor of the flavin coenzymes flavin mononucleotide and flavin adenine dinucleotide. In Escherichia coli and other bacteria, sequential deamination and reduction steps in riboflavin biosynthesis are catalyzed by RibD, a bifunctional protein with distinct pyrimidine deaminase and reductase domains. Plants have two diverged RibD homologs, PyrD and PyrR; PyrR proteins have an extra carboxyl-terminal domain (COG3236) of unknown function. Arabidopsis (Arabidopsis thaliana) PyrD (encoded by At4g20960) is known to be a monofunctional pyrimidine deaminase, but no pyrimidine reductase has been identified. Bioinformatic analyses indicated that plant PyrR proteins have a catalytically competent reductase domain but lack essential zinc-binding residues in the deaminase domain, and that the Arabidopsis PyrR gene (At3g47390) is coexpressed with riboflavin synthesis genes. These observations imply that PyrR is a pyrimidine reductase without deaminase activity. Consistent with this inference, Arabidopsis or maize (Zea mays) PyrR (At3g47390 or GRMZM2G090068) restored riboflavin prototrophy to an E. coli ribD deletant strain when coexpressed with the corresponding PyrD protein (At4g20960 or GRMZM2G320099) but not when expressed alone; the COG3236 domain was unnecessary for complementing activity. Furthermore, recombinant maize PyrR mediated NAD(P)H-dependent pyrimidine reduction in vitro. Import assays with pea (Pisum sativum) chloroplasts showed that PyrR and PyrD are taken up and proteolytically processed. Ablation of the maize PyrR gene caused early seed lethality. These data argue that PyrR is the missing plant pyrimidine reductase, that it is plastid localized, and that it is essential. The role of the COG3236 domain remains mysterious; no evidence was obtained for the possibility that it catalyzes the dephosphorylation that follows pyrimidine reduction.

Plant B vitamin pathways and their compartmentation: a guide for the perplexed.

The B vitamins and the cofactors derived from them are essential for life. B vitamin synthesis in plants is consequently as crucial to plants themselves as it is to humans and animals, whose B vitamin nutrition depends largely on plants. The synthesis and salvage pathways for the seven plant B vitamins are now broadly known, but certain enzymes and many transporters have yet to be identified, and the subcellular locations of various reactions are unclear. Although very substantial, what is not known about plant B vitamin pathways is regrettably difficult to discern from the literature or from biochemical pathway databases. Nor do databases accurately represent all that is known about B vitamin pathways-above all their compartmentation-because the facts are scattered throughout the literature, and thus hard to piece together. These problems (i) deter discoveries because newcomers to B vitamins cannot see which mysteries still need solving; and (ii) impede metabolic reconstruction and modelling of B vitamin pathways because genes for reactions or transport steps are missing. This review therefore takes a fresh approach to capture current knowledge of B vitamin pathways in plants. The synthesis pathways, key salvage routes, and their subcellular compartmentation are surveyed in depth, and encoded in the SEED database (http://pubseed.theseed.org/seedviewer.cgi?page=PlantGateway) for Arabidopsis and maize. The review itself and the encoded pathways specifically identify enigmatic or missing reactions, enzymes, and transporters. The SEED-encoded B vitamin pathway collection is a publicly available, expertly curated, one-stop resource for metabolic reconstruction and modeling.

Identification of mitochondrial thiamin diphosphate carriers from Arabidopsis and maize.

It is currently held that thiamin is made in chloroplasts and converted in the cytosol to the active cofactor thiamin diphosphate (ThDP), and that mitochondria and plastids import ThDP. The organellar transporters that mediate ThDP import in plants have not been identified. Comparative genomic analysis indicated that two members of the mitochondrial carrier family (MCF) in Arabidopsis (At5g48970 and At3g21390) and two in maize (GRMZM2G118515 and GRMZM2G124911) are related to the ThDP carriers of animals and Saccharomyces cerevisiae. Expression of each of these plant proteins in a S. cerevisiae ThDP carrier (TPC1) null mutant complemented the growth defect on fermentable carbon sources and restored the level of mitochondrial ThDP and the activity of the mitochondrial ThDP-dependent enzyme acetolactate synthase. The plant proteins were targeted to mitochondria as judged by dual import assays with purified pea mitochondria and chloroplasts, and by microscopic analysis of the subcellular localization of green fluorescent protein fusions in transiently transformed tobacco suspension cells. Both maize genes were shown to be expressed throughout the plant, which is consistent with the known ubiquity of mitochondrial ThDP-dependent enzymes. Collectively, these data establish that plants have mitochondrially located MCF carriers for ThDP, and indicate that these carriers are highly evolutionarily conserved. Our data provide a firm basis to propagate the functional annotation of mitochondrial ThDP carriers to other angiosperm genomes.

A 5-formyltetrahydrofolate cycloligase paralog from all domains of life: comparative genomic and experimental evidence for a cryptic role in thiamin metabolism.

A paralog (here termed COG0212) of the ATP-dependent folate salvage enzyme 5-formyltetrahydrofolate cycloligase (5-FCL) occurs in all domains of life and, although typically annotated as 5-FCL in pro- and eukaryotic genomes, is of unknown function. COG0212 is similar in overall structure to 5-FCL, particularly in the substrate binding region, and has distant similarity to other kinases. The Arabidopsis thaliana COG0212 protein was shown to be targeted to chloroplasts and to be required for embryo viability. Comparative genomic analysis revealed that a high proportion (19%) of archaeal and bacterial COG0212 genes are clustered on the chromosome with various genes implicated in thiamin metabolism or transport but showed no such association between COG0212 and folate metabolism. Consistent with the bioinformatic evidence for a role in thiamin metabolism, ablating COG0212 in the archaeon Haloferax volcanii caused accumulation of thiamin monophosphate. Biochemical and functional complementation tests of several known and hypothetical thiamin-related activities (involving thiamin, its breakdown products, and their phosphates) were, however, negative. Also consistent with the bioinformatic evidence, the COG0212 proteins from A. thaliana and prokaryote sources lacked 5-FCL activity in vitro and did not complement the growth defect or the characteristic 5-formyltetrahydrofolate accumulation of a 5-FCL-deficient (ΔygfA) Escherichia coli strain. We therefore propose (a) that COG0212 has an unrecognized yet sometimes crucial role in thiamin metabolism, most probably in salvage or detoxification, and (b) that is not a 5-FCL and should no longer be so annotated.