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

The importance of thiamine (vitamin B1) in plant health: From crop yield to biofortification

Ensuring that people have access to sufficient and nutritious food is necessary for a healthy life and the core tenet of food security. With the global population set to reach 9.8 billion by 2050, and the compounding effects of climate change, the planet is facing challenges that necessitate significant and rapid changes in agricultural practices. In the effort to provide food in terms of calories, the essential contribution of micronutrients (vitamins and minerals) to nutrition is often overlooked. Here, we focus on the importance of thiamine (vitamin B1) in plant health and discuss its impact on human health. Vitamin B1 is an essential dietary component, and deficiencies in this micronutrient underlie several diseases, notably nervous system disorders. The predominant source of dietary vitamin B1 is plant-based foods. Moreover, vitamin B1 is also vital for plants themselves, and its benefits in plant health have received less attention than in the human health sphere. In general, vitamin B1 is well-characterized for its role as a coenzyme in metabolic pathways, particularly those involved in energy production and central metabolism, including carbon assimilation and respiration. Vitamin B1 is also emerging as an important component of plant stress responses, and several noncoenzyme roles of this vitamin are being characterized. We summarize the importance of vitamin B1 in plants from the perspective of food security, including its roles in plant disease resistance, stress tolerance, and crop yield, and review the potential benefits of biofortification of crops with increased vitamin B1 content to improve human health.

Comparative transcriptome analysis provides key insights into seedling development in switchgrass (Panicum virgatum L.)

BACKGROUND

Switchgrass (Panicum virgatum L.), a warm-season perennial C4 plant, can be used as a forage plant, a soil and water conservation plant, a windbreak plant, and as a good source of biofuels and alternative energy with low planting costs. However, switchgrass exhibits low rates of seedling development compared to other crops, which means it is typically out-competed by weeds. There is a large variation in seedling development rates among different plantlets in switchgrass, which limits its usefulness for large-scale cultivation. Little is currently known about the molecular reasons for slow seedling growth.

RESULTS

Characterization of the seedling development process via growth indices indicated a relatively stagnant growth stage in switchgrass. A total of 678 differentially expressed genes (DEGs) were identified from the comparison of transcriptomes from slowly developed (sd) and rapidly developed (rd) switchgrass seedlings. Gene ontology and pathway enrichment analysis showed that DEGs were enriched in diterpenoid biosynthesis, thiamine metabolism, and circadian rhythm. Transcription factor enrichment and expression analyses showed MYB-related, bHLH and NAC family genes were essential for seedling growth. The transcriptome results were consistent with those of quantitative real-time polymerase chain reaction. Then, the expression profiles of maize and switchgrass were compared during seedling leaf development. A total of 128 DEGs that play key roles in seedling growth were aligned to maize genes. Transcriptional information and physiological indices suggested that several genes involved in the circadian rhythm, thiamine metabolism, energy metabolism, gibberellic acid biosynthesis, and signal transduction played important roles in seedling development.

CONCLUSIONS

The seedling development process of switchgrass was characterized, and the molecular differences between slowly developed and rapidly developed seedlings were discussed. This study provides new insights into the reasons for slow seedling development in switchgrass and will be useful for the genetic improvement of switchgrass and other crops.

A Very Oil Yellow1 Modifier of the Oil Yellow1-N1989 Allele Uncovers a Cryptic Phenotypic Impact of Cis-regulatory Variation in Maize

Forward genetics determines the function of genes underlying trait variation by identifying the change in DNA responsible for changes in phenotype. Detecting phenotypically-relevant variation outside protein coding sequences and distinguishing this from neutral variants is not trivial; partly because the mechanisms by which DNA polymorphisms in the intergenic regions affect gene regulation are poorly understood. Here we utilized a dominant genetic reporter to investigate the effect of cis and trans-acting regulatory variation. We performed a forward genetic screen for natural variation that suppressed or enhanced the semi-dominant mutant allele Oy1-N1989, encoding the magnesium chelatase subunit I of maize. This mutant permits rapid phenotyping of leaf color as a reporter for chlorophyll accumulation, and mapping of natural variation in maize affecting chlorophyll metabolism. We identified a single modifier locus segregating between B73 and Mo17 that was linked to the reporter gene itself, which we call very oil yellow1 (vey1). Based on the variation in OY1 transcript abundance and genome-wide association data, vey1 is predicted to consist of multiple cis-acting regulatory sequence polymorphisms encoded at the wild-type oy1 alleles. The vey1 locus appears to be a common polymorphism in the maize germplasm that alters the expression level of a key gene in chlorophyll biosynthesis. These vey1 alleles have no discernable impact on leaf chlorophyll in the absence of the Oy1-N1989 reporter. Thus, the use of a mutant as a reporter for magnesium chelatase activity resulted in the detection of expression-level polymorphisms not readily visible in the laboratory.

Chronic Physical Disturbance Substantially Alters the Response of Biological Soil Crusts to a Wetting Pulse, as Characterized by Metatranscriptomic Sequencing

Biological soil crusts (biocrusts) are microbial communities that are a feature of arid surface soils worldwide. In drylands where precipitation is pulsed and ephemeral, the ability of biocrust microbiota to rapidly initiate metabolic activity is critical to their survival. Community gene expression was compared after a short duration (1 h) wetting pulse in both intact and soils disturbed by chronic foot trampling. Across the metatranscriptomes the majority of transcripts were cyanobacterial in origin, suggesting that cyanobacteria accounted for the bulk of the transcriptionally active cells. Chronic trampling substantially altered the functional profile of the metatranscriptomes, specifically resulting in a significant decrease in transcripts for nitrogen fixation. Soil depth (biocrust and below crust) was a relatively small factor in differentiating the metatranscriptomes, suggesting that the metabolically active bacteria were similar between shallow soil horizons. The dry samples were consistently enriched for hydrogenase genes, indicating that molecular hydrogen may serve as an energy source for the desiccated soil communities. The water pulse was associated with a restructuring of the metatranscriptome, particularly for the biocrusts. Biocrusts increased transcripts for photosynthesis and carbon fixation, suggesting a rapid resuscitation upon wetting. In contrast, the trampled surface soils showed a much smaller response to wetting, indicating that trampling altered the metabolic response of the community. Finally, several biogeochemical cycling genes in carbon and nitrogen cycling were assessed for their change in abundance due to wetting in the biocrusts. Different transcripts encoding the same gene product did not show a consensus response, with some more abundant in dry or wet biocrusts, highlighting the challenges in relating transcript abundance to biogeochemical cycling rates. These observations demonstrate that metatranscriptome sequencing was able to distinguish alterations in the function of arid soil microbial communities at two varying temporal scales, a long-term ecosystems disturbance through foot trampling, and a short term wetting pulse. Thus, community metatranscriptomes have the potential to inform studies on the response and resilience of biocrusts to various environmental perturbations.

Thiaminase I Provides a Growth Advantage by Salvaging Precursors from Environmental Thiamine and Its Analogs in Burkholderia thailandensis

Thiamine is essential to life, as it serves as a cofactor for enzymes involved in critical carbon transformations. Many bacteria can synthesize thiamine, while thiamine auxotrophs must obtain it or its precursors from the environment. Thiaminases degrade thiamine by catalyzing the base-exchange substitution of thiazole with a nucleophile, and thiaminase I specifically has been implicated in thiamine deficiency syndromes in animals. The biological role of this secreted enzyme has been a long-standing mystery. We used the thiaminase I-producing soil bacterium Burkholderia thailandensis as a model to ascertain its function. First, we generated thiamine auxotrophs, which are still able to use exogenous precursors (thiazole and hydroxymethyl pyrimidine), to synthesize thiamine. We found that thiaminase I extended the survival of these strains, when grown in defined media where thiamine was serially diluted out, compared to isogenic strains that could not produce thiaminase I. Thiamine auxotrophs grew better on thiamine precursors than thiamine itself, suggesting thiaminase I functions to convert thiamine to useful precursors. Furthermore, our findings demonstrate that thiaminase I cleaves phosphorylated thiamine and toxic analogs, which releases precursors that can then be used for thiamine synthesis. This study establishes a biological role for this perplexing enzyme and provides additional insight into the complicated nature of thiamine metabolism and how individual bacteria may manipulate the availability of a vital nutrient in the environment.IMPORTANCE The function of thiaminase I has remained a long-standing, unsolved mystery. The enzyme is only known to be produced by a small subset of microorganisms, although thiaminase I activity has been associated with numerous plants and animals, and is implicated in thiamine deficiencies brought on by consumption of organisms containing this enzyme. Genomic and biochemical analyses have shed light on potential roles for the enzyme. Using the genetically amenable thiaminase I-producing soil bacterium Burkholderia thailandensis, we were able to demonstrate that thiaminase I helps salvage precursors from thiamine derivatives in the environment and degrades thiamine to its precursors, which are preferentially used by B. thailandensis auxotrophs. Our study establishes a biological role for this perplexing enzyme and provides insight into the complicated nature of thiamine metabolism. It also establishes B. thailandensis as a robust model system for studying thiamine metabolism.

Thiaminase I Provides a Growth Advantage by Salvaging Precursors from Environmental Thiamine and Its Analogs in Burkholderia thailandensis

Thiamine is essential to life, as it serves as a cofactor for enzymes involved in critical carbon transformations. Many bacteria can synthesize thiamine, while thiamine auxotrophs must obtain it or its precursors from the environment. Thiaminases degrade thiamine by catalyzing the base-exchange substitution of thiazole with a nucleophile, and thiaminase I specifically has been implicated in thiamine deficiency syndromes in animals. The biological role of this secreted enzyme has been a long-standing mystery. We used the thiaminase I-producing soil bacterium Burkholderia thailandensis as a model to ascertain its function. First, we generated thiamine auxotrophs, which are still able to use exogenous precursors (thiazole and hydroxymethyl pyrimidine), to synthesize thiamine. We found that thiaminase I extended the survival of these strains, when grown in defined media where thiamine was serially diluted out, compared to isogenic strains that could not produce thiaminase I. Thiamine auxotrophs grew better on thiamine precursors than thiamine itself, suggesting thiaminase I functions to convert thiamine to useful precursors. Furthermore, our findings demonstrate that thiaminase I cleaves phosphorylated thiamine and toxic analogs, which releases precursors that can then be used for thiamine synthesis. This study establishes a biological role for this perplexing enzyme and provides additional insight into the complicated nature of thiamine metabolism and how individual bacteria may manipulate the availability of a vital nutrient in the environment.IMPORTANCE The function of thiaminase I has remained a long-standing, unsolved mystery. The enzyme is only known to be produced by a small subset of microorganisms, although thiaminase I activity has been associated with numerous plants and animals, and is implicated in thiamine deficiencies brought on by consumption of organisms containing this enzyme. Genomic and biochemical analyses have shed light on potential roles for the enzyme. Using the genetically amenable thiaminase I-producing soil bacterium Burkholderia thailandensis, we were able to demonstrate that thiaminase I helps salvage precursors from thiamine derivatives in the environment and degrades thiamine to its precursors, which are preferentially used by B. thailandensis auxotrophs. Our study establishes a biological role for this perplexing enzyme and provides insight into the complicated nature of thiamine metabolism. It also establishes B. thailandensis as a robust model system for studying thiamine metabolism.

A strictly monofunctional bacterial hydroxymethylpyrimidine phosphate kinase precludes damaging errors in thiamin biosynthesis

The canonical kinase (ThiD) that converts the thiamin biosynthesis intermediate hydroxymethylpyrimidine (HMP) monophosphate into the diphosphate can also very efficiently convert free HMP into the monophosphate in prokaryotes, plants, and fungi. This HMP kinase activity enables salvage of HMP, but it is not substrate-specific and so allows toxic HMP analogs and damage products to infiltrate the thiamin biosynthesis pathway. Comparative analysis of bacterial genomes uncovered a gene, thiD2, that is often fused to the thiamin synthesis gene thiE and could potentially encode a replacement for ThiD. Standalone ThiD2 proteins and ThiD2 fusion domains are small (∼130 residues) and do not belong to any previously known protein family. Genetic and biochemical analyses showed that representative standalone and fused ThiD2 proteins catalyze phosphorylation of HMP monophosphate, but not of HMP or its toxic analogs and damage products such as bacimethrin and 5-(hydroxymethyl)-2-methylpyrimidin-4-ol. As strictly monofunctional HMP monophosphate kinases, ThiD2 proteins eliminate a potentially fatal vulnerability of canonical ThiD, at the cost of the ability to reclaim HMP formed by thiamin turnover.

Genome-scale metabolic modelling common cofactors metabolism in microorganisms

The common cofactors ATP/ADP, NAD(P)(H), and acetyl-CoA/CoA are indispensable participants in biochemical reactions in industrial microbes. To systematically explore the effects of these cofactors on cell growth and metabolic phenotypes, the first genome-scale cofactor metabolic model, icmNX6434, including 6434 genes, 1782 metabolites, and 6877 reactions, was constructed from 14 genome-scale metabolic models of 14 industrial strains. The origin, consumption, and interactions of these common cofactors in microbial cells were elucidated by the icmNX6434 model, and they played important roles in cell growth. The essential cofactor modules contained 2480 genes and 2948 reactions; therefore, improving cofactor biosynthesis, directing these cofactors into essential metabolic pathways, as well as avoiding cofactor utilization during byproduct biosynthesis and futile cycles, are three ways to increase cell growth. The effects of these common cofactors on the distribution and rate of the carbon flux in four universal modes, as well as an optimized metabolic flux, could be obtained by manipulating cofactor availability and balance. Significant changes in the ATP, NAD(H), NADP(H), or acetyl-CoA concentrations triggered relevant metabolic responses to acidic, oxidative, heat, and osmotic stress. Globally, the model icmNX6434 provides a comprehensive platform to elucidate the physiological effects of these cofactors on cell growth, metabolic flux, and industrial robustness. Moreover, the results of this study are a further example of using a consensus genome-scale metabolic model to increase our understanding of key biological processes.

Functional Annotations of Paralogs: A Blessing and a Curse

Gene duplication followed by mutation is a classic mechanism of neofunctionalization, producing gene families with functional diversity. In some cases, a single point mutation is sufficient to change the substrate specificity and/or the chemistry performed by an enzyme, making it difficult to accurately separate enzymes with identical functions from homologs with different functions. Because sequence similarity is often used as a basis for assigning functional annotations to genes, non-isofunctional gene families pose a great challenge for genome annotation pipelines. Here we describe how integrating evolutionary and functional information such as genome context, phylogeny, metabolic reconstruction and signature motifs may be required to correctly annotate multifunctional families. These integrative analyses can also lead to the discovery of novel gene functions, as hints from specific subgroups can guide the functional characterization of other members of the family. We demonstrate how careful manual curation processes using comparative genomics can disambiguate subgroups within large multifunctional families and discover their functions. We present the COG0720 protein family as a case study. We also discuss strategies to automate this process to improve the accuracy of genome functional annotation pipelines.

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

Conserved active site cysteine residue of archaeal THI4 homolog is essential for thiamine biosynthesis in Haloferax volcanii

Background

Thiamine (vitamin B1) is synthesized de novo by certain yeast, fungi, plants, protozoans, bacteria and archaea. The pathway of thiamine biosynthesis by archaea is poorly understood, particularly the route of sulfur relay to form the thiazole ring. Archaea harbor structural homologs of both the bacterial (ThiS-ThiF) and eukaryotic (THI4) proteins that mobilize sulfur to thiazole ring precursors by distinct mechanisms.

Results

Based on comparative genome analysis, halophilic archaea are predicted to synthesize the pyrimidine moiety of thiamine by the bacterial pathway, initially suggesting that also a bacterial ThiS-ThiF type mechanism for synthesis of the thiazole ring is used in which the sulfur carrier ThiS is first activated by ThiF-catalyzed adenylation. The only ThiF homolog of Haloferax volcanii (UbaA) was deleted but this had no effect on growth in the absence of thiamine. Usage of the eukaryotic THI4-type sulfur relay was initially considered less likely for thiamine biosynthesis in archaea, since the active-site cysteine residue of yeast THI4p that donates the sulfur to the thiazole ring by a suicide mechanism is replaced by a histidine residue in many archaeal THI4 homologs and these are described as D-ribose-1,5-bisphosphate isomerases. The THI4 homolog of the halophilic archaea, including Hfx. volcanii (HVO_0665, HvThi4) was found to differ from that of methanogens and thermococci by having a cysteine residue (Cys165) corresponding to the conserved active site cysteine of yeast THI4p (Cys205). Deletion of HVO_0665 generated a thiamine auxotroph that was trans-complemented by a wild-type copy of HVO_0665, but not the modified gene encoding an HvThi4 C165A variant.

Conclusions

Based on our results, we conclude that the archaeon Hfx. volcanii uses a yeast THI4-type mechanism for sulfur relay to form the thiazole ring of thiamine. We extend this finding to a relatively large group of archaea, including haloarchaea, ammonium oxidizing archaea, and some methanogen and Pyrococcus species, by observing that these organisms code for THI4 homologs that have a conserved active site cysteine residue which is likely used in thiamine biosynthesis. Thus, archaeal members of IPR002922 THI4 family that have a conserved cysteine active site should be reexamined for a function in thiamine biosynthesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-014-0260-0) contains supplementary material, which is available to authorized users.

Salvage of the thiamin pyrimidine moiety by plant TenA proteins lacking an active-site cysteine

The TenA protein family occurs in prokaryotes, plants and fungi; it has two subfamilies, one (TenA_C) having an active-site cysteine, the other (TenA_E) not. TenA_C proteins participate in thiamin salvage by hydrolysing the thiamin breakdown product amino-HMP (4-amino-5-aminomethyl-2-methylpyrimidine) to HMP (4-amino-5-hydroxymethyl-2-methylpyrimidine); the function of TenA_E proteins is unknown. Comparative analysis of prokaryote and plant genomes predicted that (i) TenA_E has a salvage role similar to, but not identical with, that of TenA_C and (ii) that TenA_E and TenA_C also have non-salvage roles since they occur in organisms that cannot make thiamin. Recombinant Arabidopsis and maize TenA_E proteins (At3g16990, GRMZM2G080501) hydrolysed amino-HMP to HMP and, far more actively, hydrolysed the N-formyl derivative of amino-HMP to amino-HMP. Ablating the At3g16990 gene in a line with a null mutation in the HMP biosynthesis gene ThiC prevented its rescue by amino-HMP. Ablating At3g16990 in the wild-type increased sensitivity to paraquat-induced oxidative stress; HMP overcame this increased sensitivity. Furthermore, the expression of TenA_E and ThiC genes in Arabidopsis and maize was inversely correlated. These results indicate that TenA_E proteins mediate amidohydrolase and aminohydrolase steps in the salvage of thiamin breakdown products. As such products can be toxic, TenA_E proteins may also pre-empt toxicity.

Divisions of labor in the thiamin biosynthetic pathway among organs of maize

The B vitamin thiamin is essential for central metabolism in all cellular organisms including plants. While plants synthesize thiamin de novo, organs vary widely in their capacities for thiamin synthesis. We use a transcriptomics approach to appraise the distribution of de novo synthesis and thiamin salvage pathways among organs of maize. We identify at least six developmental contexts in which metabolically active, non-photosynthetic organs exhibit low expression of one or both branches of the de novo thiamin biosynthetic pathway indicating a dependence on inter-cellular transport of thiamin and/or thiamin precursors. Neither the thiazole (THI4) nor pyrimidine (THIC) branches of the pathway are expressed in developing pollen implying a dependence on import of thiamin from surrounding floral and inflorescence organs. Consistent with that hypothesis, organs of the male inflorescence and flowers are shown to have high relative expression of the thiamin biosynthetic pathway and comparatively high thiamin contents. By contrast, divergent patterns of THIC and THI4 expression occur in the shoot apical meristem, embyro sac, embryo, endosperm, and root-tips suggesting that these sink organs acquire significant amounts of thiamin via salvage pathways. In the root and shoot meristems, expression of THIC in the absence of THI4 indicates a capacity for thiamin synthesis via salvage of thiazole, whereas the opposite pattern obtains in embryo and endosperm implying that seed storage organs are poised for pyrimidine salvage. Finally, stable isotope labeling experiments set an upper limit on the rate of de novo thiamin biosynthesis in maize leaf explants. Overall, the observed patterns of thiamin biosynthetic gene expression mirror the strategies for thiamin acquisition that have evolved in bacteria.

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.

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.