Genome-scale metabolic models (GEMs) represent extensive knowledgebases that provide a platform for model simulations and integrative analysis of omics data. This study introduces Yeast8 and an associated ecosystem of models that represent a comprehensive computational resource for performing simulations of the metabolism of Saccharomyces cerevisiae––an important model organism and widely used cell-factory. Yeast8 tracks community development with version control, setting a standard for how GEMs can be continuously updated in a simple and reproducible way. We use Yeast8 to develop the derived models panYeast8 and coreYeast8, which in turn enable the reconstruction of GEMs for 1,011 different yeast strains. Through integration with enzyme constraints (ecYeast8) and protein 3D structures (proYeast8DB), Yeast8 further facilitates the exploration of yeast metabolism at a multi-scale level, enabling prediction of how single nucleotide variations translate to phenotypic traits.

Enzyme turnover numbers (kcat) are key to understanding cellular metabolism, proteome allocation and physiological diversity, but experimentally measured kcat data are sparse and noisy. Here we provide a deep learning approach (DLKcat) for high-throughput kcat prediction for metabolic enzymes from any organism merely from substrate structures and protein sequences. DLKcat can capture kcat changes for mutated enzymes and identify amino acid residues with a strong impact on kcat values. We applied this approach to predict genome-scale kcat values for more than 300 yeast species. Additionally, we designed a Bayesian pipeline to parameterize enzyme-constrained genome-scale metabolic models from predicted kcat values. The resulting models outperformed the corresponding original enzyme-constrained genome-scale metabolic models from previous pipelines in predicting phenotypes and proteomes, and enabled us to explain phenotypic differences. DLKcat and the enzyme-constrained genome-scale metabolic model construction pipeline are valuable tools to uncover global trends of enzyme kinetics and physiological diversity, and to further elucidate cellular metabolism on a large scale.

No abstract

BackgroundAspergillus niger is widely used for enzyme production and achievement of high enzyme production depends on the comprehensive understanding of cell’s metabolic regulation mechanisms.ResultsIn this paper, we investigate the metabolic differences and regulation mechanisms between a high glucoamylase-producing strain A. niger DS03043 and its wild-type parent strain A. niger CBS513.88 via an integrated isotope-assisted metabolomics and 13C metabolic flux analysis approach. We found that A. niger DS03043 had higher cell growth, glucose uptake, and glucoamylase production rates but lower oxalic acid and citric acid secretion rates. In response to above phenotype changes, A. niger DS03043 was characterized by an increased carbon flux directed to the oxidative pentose phosphate pathway in contrast to reduced flux through TCA cycle, which were confirmed by consistent changes in pool sizes of metabolites. A higher ratio of ATP over AMP in the high producing strain might contribute to the increase in the PP pathway flux as glucosephosphate isomerase was inhibited at higher ATP concentrations. A. niger CBS513.88, however, was in a higher redox state due to the imbalance of NADH regeneration and consumption, resulting in the secretion of oxalic acid and citric acid, as well as the accumulation of intracellular OAA and PEP, which may in turn result in the decrease in the glucose uptake rate.ConclusionsThe application of integrated metabolomics and 13C metabolic flux analysis highlights the regulation mechanisms of energy and redox metabolism on flux redistribution in A. niger. Graphical abstractAn integrated isotope-assisted metabolomics and 13C metabolic flux analysis was was firstly systematically performed in A. niger. In response to enzyme production, the metabolic flux in A. niger DS03043 (high-producing) was redistributed, characterized by an increased carbon flux directed to the oxidative pentose phosphate pathway as well as an increased pool size of pentose. The consistency in 13C metabolic flux analysis and metabolites quantification indicated that an imbalance of NADH formation and consumption led to the accumulation and secretion of organic acids in A. niger CBS513.88 (wild-type)Electronic supplementary materialThe online version of this article (doi:10.1186/s12934-015-0329-y) contains supplementary material, which is available to authorized users.

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome-scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence-level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.