Reconstruction and validation of Saccharomyces cerevisiae iND750, a fully compartmentalized genome-scale metabolic model

Optimizing Metabolite Production with Neighborhood-Based Binary Quantum-Behaved Particle Swarm Optimization and Flux Balance Analysis

Metabolic engineering is a rapidly evolving field that involves optimizing microbial cell factories to overproduce various industrial products. To achieve this, several tools, leveraging constraint-based stoichiometric models and metaheuristic algorithms like particle swarm optimization (PSO), have been developed. However, PSO can potentially get trapped in local optima. Quantum-behaved PSO (QPSO) overcomes this limitation, and our study further enhances its binary version (BQPSO) with a neighborhood topology, leading to the advanced neighborhood-based BQPSO (NBQPSO). Combined with flux balance analysis (FBA), this forms an innovative approach, NBQPSO-FBA, for identifying optimal knockout strategies to maximize the desired metabolite production. Additionally, we introduced a novel encoding strategy suitable for large-scale genome-scale metabolic models (GSMMs). Evaluated on four E. coli GSMMs (iJR904, iAF1260, iJO1366, and iML1515), NBQPSO-FBA matches or surpasses established bi-level linear programming (LP) and heuristic methods in metabolite production optimization. Notably, it achieved 90.69% realization of the theoretical maximum in acetate production and demonstrated comparable performance with leading algorithms in lactate production. The efficiency of NBQPSO-FBA, which requires fewer knockouts, makes it a practical and effective tool for optimizing microbial cell factories. This addresses the rising demand for microbial products across various industries.

Auxotrophy-based curation improves the consensus genome-scale metabolic model of yeast

Saccharomyces cerevisiae, a widely utilized model organism, has seen continuous updates to its genome-scale metabolic model (GEM) to enhance the prediction performance for metabolic engineering and systems biology. This study presents an auxotrophy-based curation of the yeast GEM, enabling facile upgrades to yeast GEMs in future endeavors. We illustrated that the curation bolstered the predictive capability of the yeast GEM particularly in predicting auxotrophs without compromising accuracy in other simulations, and thus could be an effective manner for GEM refinement. Last, we leveraged the curated yeast GEM to systematically predict auxotrophs, thereby furnishing a valuable reference for the design of nutrient-dependent cell factories and synthetic yeast consortia.

Genome-Scale Metabolic Models in Fungal Pathogens: Past, Present, and Future

Fungi are diverse organisms with various characteristics and functions. Some play a role in recycling essential elements, such as nitrogen and carbon, while others are utilized in the food and drink production industry. Some others are known to cause diseases in various organisms, including humans. Fungal pathogens cause superficial, subcutaneous, and systemic infections. Consequently, many scientists have focused on studying the factors contributing to the development of human diseases. Therefore, multiple approaches have been assessed to examine the biology of these intriguing organisms. The genome-scale metabolic models (GEMs) have demonstrated many advantages to microbial metabolism studies and the ability to propose novel therapeutic alternatives. Despite significant advancements, much remains to be elucidated regarding the use of this tool for investigating fungal metabolism. This review aims to compile the data provided by the published GEMs of human fungal pathogens. It gives specific examples of the most significant contributions made by these models, examines the advantages and difficulties associated with using such models, and explores the novel approaches suggested to enhance and refine their development.

Yeast9: a consensus genome-scale metabolic model for S. cerevisiae curated by the community

Genome-scale metabolic models (GEMs) can facilitate metabolism-focused multi-omics integrative analysis. Since Yeast8, the yeast-GEM of Saccharomyces cerevisiae, published in 2019, has been continuously updated by the community. This has increased the quality and scope of the model, culminating now in Yeast9. To evaluate its predictive performance, we generated 163 condition-specific GEMs constrained by single-cell transcriptomics from osmotic pressure or reference conditions. Comparative flux analysis showed that yeast adapting to high osmotic pressure benefits from upregulating fluxes through central carbon metabolism. Furthermore, combining Yeast9 with proteomics revealed metabolic rewiring underlying its preference for nitrogen sources. Lastly, we created strain-specific GEMs (ssGEMs) constrained by transcriptomics for 1229 mutant strains. Well able to predict the strains' growth rates, fluxomics from those large-scale ssGEMs outperformed transcriptomics in predicting functional categories for all studied genes in machine learning models. Based on those findings we anticipate that Yeast9 will continue to empower systems biology studies of yeast metabolism.

Genome-scale metabolic models in translational medicine: the current status and potential of machine learning in improving the effectiveness of the models

The genome-scale metabolic model (GEM) has emerged as one of the leading modeling approaches for systems-level metabolic studies and has been widely explored for a broad range of organisms and applications. Owing to the development of genome sequencing technologies and available biochemical data, it is possible to reconstruct GEMs for model and non-model microorganisms as well as for multicellular organisms such as humans and animal models. GEMs will evolve in parallel with the availability of biological data, new mathematical modeling techniques and the development of automated GEM reconstruction tools. The use of high-quality, context-specific GEMs, a subset of the original GEM in which inactive reactions are removed while maintaining metabolic functions in the extracted model, for model organisms along with machine learning (ML) techniques could increase their applications and effectiveness in translational research in the near future. Here, we briefly review the current state of GEMs, discuss the potential contributions of ML approaches for more efficient and frequent application of these models in translational research, and explore the extension of GEMs to integrative cellular models.

Genome-scale metabolic network models for industrial microorganisms metabolic engineering: Current advances and future prospects

The construction of high-performance microbial cell factories (MCFs) is the centerpiece of biomanufacturing. However, the complex metabolic regulatory network of microorganisms poses great challenges for the efficient design and construction of MCFs. The genome-scale metabolic network models (GSMs) can systematically simulate the metabolic regulation process of microorganisms in silico, providing effective guidance for the rapid design and construction of MCFs. In this review, we summarized the development status of 16 important industrial microbial GSMs, and further outline the technologies or methods that continuously promote high-quality GSMs construction from five aspects: I) Databases and modeling tools facilitate GSMs reconstruction; II) evolving gap-filling technologies; III) constraint-based model reconstruction; IV) advances in algorithms; and V) developed visualization tools. In addition, we also summarized the applications of GSMs in guiding metabolic engineering from four aspects: I) exploring and explaining metabolic features; II) predicting the effects of genetic perturbations on metabolism; III) predicting the optimal phenotype; IV) guiding cell factories construction in practical experiment. Finally, we discussed the development of GSMs, aiming to provide a reference for efficiently reconstructing GSMs and guiding metabolic engineering.

Microbial Pathway Thermodynamics: Stoichiometric Models Unveil Anabolic and Catabolic Processes

The biotechnological exploitation of microorganisms enables the use of metabolism for the production of economically valuable substances, such as drugs or food. It is, thus, unsurprising that the investigation of microbial metabolism and its regulation has been an active research field for many decades. As a result, several theories and techniques were developed that allow for the prediction of metabolic fluxes and yields as biotechnologically relevant output parameters. One important approach is to derive macrochemical equations that describe the overall metabolic conversion of an organism and basically treat microbial metabolism as a black box. The opposite approach is to include all known metabolic reactions of an organism to assemble a genome-scale metabolic model. Interestingly, both approaches are rather successful at characterizing and predicting the expected product yield. Over the years, macrochemical equations especially have been extensively characterized in terms of their thermodynamic properties. However, a common challenge when characterizing microbial metabolism by a single equation is to split this equation into two, describing the two modes of metabolism, anabolism and catabolism. Here, we present strategies to systematically identify separate equations for anabolism and catabolism. Based on metabolic models, we systematically identify all theoretically possible catabolic routes and determine their thermodynamic efficiency. We then show how anabolic routes can be derived, and we use these to approximate biomass yield. Finally, we challenge the view of metabolism as a linear energy converter, in which the free energy gradient of catabolism drives the anabolic reactions.

FastKnock: an efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed by making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies (with a predefined maximum number of reaction deletions) for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using various Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple- and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more beneficial and important practical solutions. The availability of all the solutions provides the opportunity to further characterize, rank and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock .

ICON-GEMs: integration of co-expression network in genome-scale metabolic models, shedding light through systems biology

BACKGROUND

Flux Balance Analysis (FBA) is a key metabolic modeling method used to simulate cellular metabolism under steady-state conditions. Its simplicity and versatility have led to various strategies incorporating transcriptomic and proteomic data into FBA, successfully predicting flux distribution and phenotypic results. However, despite these advances, the untapped potential lies in leveraging gene-related connections like co-expression patterns for valuable insights.

RESULTS

To fill this gap, we introduce ICON-GEMs, an innovative constraint-based model to incorporate gene co-expression network into the FBA model, facilitating more precise determination of flux distributions and functional pathways. In this study, transcriptomic data from both Escherichia coli and Saccharomyces cerevisiae were integrated into their respective genome-scale metabolic models. A comprehensive gene co-expression network was constructed as a global view of metabolic mechanism of the cell. By leveraging quadratic programming, we maximized the alignment between pairs of reaction fluxes and the correlation of their corresponding genes in the co-expression network. The outcomes notably demonstrated that ICON-GEMs outperformed existing methodologies in predictive accuracy. Flux variabilities over subsystems and functional modules also demonstrate promising results. Furthermore, a comparison involving different types of biological networks, including protein-protein interactions and random networks, reveals insights into the utilization of the co-expression network in genome-scale metabolic engineering.

CONCLUSION

ICON-GEMs introduce an innovative constrained model capable of simultaneous integration of gene co-expression networks, ready for board application across diverse transcriptomic data sets and multiple organisms. It is freely available as open-source at https://github.com/ThummaratPaklao/ICOM-GEMs.git .

Applications of genome-scale metabolic models to investigate microbial metabolic adaptations in response to genetic or environmental perturbations

Environmental perturbations are encountered by microorganisms regularly and will require metabolic adaptations to ensure an organism can survive in the newly presenting conditions. In order to study the mechanisms of metabolic adaptation in such conditions, various experimental and computational approaches have been used. Genome-scale metabolic models (GEMs) are one of the most powerful approaches to study metabolism, providing a platform to study the systems level adaptations of an organism to different environments which could otherwise be infeasible experimentally. In this review, we are describing the application of GEMs in understanding how microbes reprogram their metabolic system as a result of environmental variation. In particular, we provide the details of metabolic model reconstruction approaches, various algorithms and tools for model simulation, consequences of genetic perturbations, integration of '-omics' datasets for creating context-specific models and their application in studying metabolic adaptation due to the change in environmental conditions.

Genome-scale metabolic models reveal determinants of phenotypic differences in non-Saccharomyces yeasts

BACKGROUND

Use of alternative non-Saccharomyces yeasts in wine and beer brewing has gained more attention the recent years. This is both due to the desire to obtain a wider variety of flavours in the product and to reduce the final alcohol content. Given the metabolic differences between the yeast species, we wanted to account for some of the differences by using in silico models.

RESULTS

We created and studied genome-scale metabolic models of five different non-Saccharomyces species using an automated processes. These were: Metschnikowia pulcherrima, Lachancea thermotolerans, Hanseniaspora osmophila, Torulaspora delbrueckii and Kluyveromyces lactis. Using the models, we predicted that M. pulcherrima, when compared to the other species, conducts more respiration and thus produces less fermentation products, a finding which agrees with experimental data. Complex I of the electron transport chain was to be present in M. pulcherrima, but absent in the others. The predicted importance of Complex I was diminished when we incorporated constraints on the amount of enzymatic protein, as this shifts the metabolism towards fermentation.

CONCLUSIONS

Our results suggest that Complex I in the electron transport chain is a key differentiator between Metschnikowia pulcherrima and the other yeasts considered. Yet, more annotations and experimental data have the potential to improve model quality in order to increase fidelity and confidence in these results. Further experiments should be conducted to confirm the in vivo effect of Complex I in M. pulcherrima and its respiratory metabolism.

The limitations of phenotype prediction in metabolism

Phenotype prediction is at the center of many questions in biology. Prediction is often achieved by determining statistical associations between genetic and phenotypic variation, ignoring the exact processes that cause the phenotype. Here, we present a framework based on genome-scale metabolic reconstructions to reveal the mechanisms behind the associations. We calculated a polygenic score (PGS) that identifies a set of enzymes as predictors of growth, the phenotype. This set arises from the synergy of the functional mode of metabolism in a particular setting and its evolutionary history, and is suitable to infer the phenotype across a variety of conditions. We also find that there is optimal genetic variation for predictability and demonstrate how the linear PGS can still explain phenotypes generated by the underlying nonlinear biochemistry. Therefore, the explicit model interprets the black box statistical associations of the genotype-to-phenotype map and helps to discover what limits the prediction in metabolism.

A structured evaluation of genome-scale constraint-based modeling tools for microbial consortia

Harnessing the power of microbial consortia is integral to a diverse range of sectors, from healthcare to biotechnology to environmental remediation. To fully realize this potential, it is critical to understand the mechanisms behind the interactions that structure microbial consortia and determine their functions. Constraint-based reconstruction and analysis (COBRA) approaches, employing genome-scale metabolic models (GEMs), have emerged as the state-of-the-art tool to simulate the behavior of microbial communities from their constituent genomes. In the last decade, many tools have been developed that use COBRA approaches to simulate multi-species consortia, under either steady-state, dynamic, or spatiotemporally varying scenarios. Yet, these tools have not been systematically evaluated regarding their software quality, most suitable application, and predictive power. Hence, it is uncertain which tools users should apply to their system and what are the most urgent directions that developers should take in the future to improve existing capacities. This study conducted a systematic evaluation of COBRA-based tools for microbial communities using datasets from two-member communities as test cases. First, we performed a qualitative assessment in which we evaluated 24 published tools based on a list of FAIR (Findability, Accessibility, Interoperability, and Reusability) features essential for software quality. Next, we quantitatively tested the predictions in a subset of 14 of these tools against experimental data from three different case studies: a) syngas fermentation by C. autoethanogenum and C. kluyveri for the static tools, b) glucose/xylose fermentation with engineered E. coli and S. cerevisiae for the dynamic tools, and c) a Petri dish of E. coli and S. enterica for tools incorporating spatiotemporal variation. Our results show varying performance levels of the best qualitatively assessed tools when examining the different categories of tools. The differences in the mathematical formulation of the approaches and their relation to the results were also discussed. Ultimately, we provide recommendations for refining future GEM microbial modeling tools.

FastKnock: An efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using three Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more useful and important practical solutions. The availability of all the solutions provides the opportunity to further characterize and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock.

Moving beyond DNA: towards functional analysis of the vaginal microbiome by non-sequencing-based methods

Over the last two decades, sequencing-based methods have revolutionised our understanding of niche-specific microbial complexity. In the lower female reproductive tract, these approaches have enabled identification of bacterial compositional structures associated with health and disease. Application of metagenomics and metatranscriptomics strategies have provided insight into the putative function of these communities but it is increasingly clear that direct measures of microbial and host cell function are required to understand the contribution of microbe-host interactions to pathophysiology. Here we explore and discuss current methods and approaches, many of which rely upon mass-spectrometry, being used to capture functional insight into the vaginal mucosal interface. In addition to improving mechanistic understanding, these methods offer innovative solutions for the development of diagnostic and therapeutic strategies designed to improve women's health.

In silico-guided metabolic engineering of Bacillus subtilis for efficient biosynthesis of purine nucleosides by blocking the key backflow nodes

Background

Purine nucleosides play essential roles in cellular physiological processes and have a wide range of applications in the fields of antitumor/antiviral drugs and food. However, microbial overproduction of purine nucleosides by de novo metabolic engineering remains a great challenge due to their strict and complex regulatory machinery involved in biosynthetic pathways.

Results

In this study, we designed an in silico-guided strategy for overproducing purine nucleosides based on a genome-scale metabolic network model in Bacillus subtilis. The metabolic flux was analyzed to predict two key backflow nodes, Drm (purine nucleotides toward PPP) and YwjH (PPP–EMP), to resolve the competitive relationship between biomass and purine nucleotide synthesis. In terms of the purine synthesis pathway, the first backflow node Drm was inactivated to block the degradation of purine nucleotides, which greatly increased the inosine production to 13.98–14.47 g/L without affecting cell growth. Furthermore, releasing feedback inhibition of the purine operon by promoter replacement enhanced the accumulation of purine nucleotides. In terms of the central carbon metabolic pathways, the deletion of the second backflow node YwjH and overexpression of Zwf were combined to increase inosine production to 22.01 ± 1.18 g/L by enhancing the metabolic flow of PPP. By switching on the flux node of the glucose-6-phosphate to PPP or EMP, the final inosine engineered strain produced up to 25.81 ± 1.23 g/L inosine by a pgi-based metabolic switch with a yield of 0.126 mol/mol glucose, a productivity of 0.358 g/L/h and a synthesis rate of 0.088 mmol/gDW/h, representing the highest yield in de novo engineered inosine bacteria. Under the guidance of this in silico-designed strategy, a general chassis bacterium was generated, for the first time, to efficiently synthesize inosine, adenosine, guanosine, IMP and GMP, which provides sufficient precursors for the synthesis of various purine intermediates.

Conclusions

Our study reveals that in silico-guided metabolic engineering successfully optimized the purine synthesis pathway by exploring efficient targets, which could be applied as a superior strategy for efficient biosynthesis of biotechnological products.

A pipeline for the reconstruction and evaluation of context-specific human metabolic models at a large-scale

Constraint-based (CB) metabolic models provide a mathematical framework and scaffold for in silico cell metabolism analysis and manipulation. In the past decade, significant efforts have been done to model human metabolism, enabled by the increased availability of multi-omics datasets and curated genome-scale reconstructions, as well as the development of several algorithms for context-specific model (CSM) reconstruction. Although CSM reconstruction has revealed insights on the deregulated metabolism of several pathologies, the process of reconstructing representative models of human tissues still lacks benchmarks and appropriate integrated software frameworks, since many tools required for this process are still disperse across various software platforms, some of which are proprietary. In this work, we address this challenge by assembling a scalable CSM reconstruction pipeline capable of integrating transcriptomics data in CB models. We combined omics preprocessing methods inspired by previous efforts with in-house implementations of existing CSM algorithms and new model refinement and validation routines, all implemented in the Troppo Python-based open-source framework. The pipeline was validated with multi-omics datasets from the Cancer Cell Line Encyclopedia (CCLE), also including reference fluxomics measurements for the MCF7 cell line. We reconstructed over 6000 models based on the Human-GEM template model for 733 cell lines featured in the CCLE, using MCF7 models as reference to find the best parameter combinations. These reference models outperform earlier studies using the same template by comparing gene essentiality and fluxomics experiments. We also analysed the heterogeneity of breast cancer cell lines, identifying key changes in metabolism related to cancer aggressiveness. Despite the many challenges in CB modelling, we demonstrate using our pipeline that combining transcriptomics data in metabolic models can be used to investigate key metabolic shifts. Significant limitations were found on these models ability for reliable quantitative flux prediction, thus motivating further work in genome-wide phenotype prediction.

Genome-scale models of human metabolism are promising tools capable of contextualising large omics datasets within a framework that enables analysis and manipulation of metabolic phenotypes. Despite various successes in applying these methods to provide mechanistic hypotheses for deregulated metabolism in disease, there is no standardized workflow to extract these models using existing methods and the tools required to do so are mostly implemented using proprietary software. We have assembled a generic pipeline to extract and validate context-specific metabolic models using multi-omics datasets and implemented it using the troppo framework. We first validate our pipeline using MCF7 cell line models and assess their ability to predict lethal gene knockouts as well as flux activity using multi-omics data. We also demonstrate how this approach can be generalized for large-scale transcriptomics datasets and used to generate insights on the metabolic heterogeneity of cancer and relevant features for other data mining approaches. The pipeline is available as part of an open-source framework that is generic for a variety of applications.

Dissecting carbon metabolism of Yarrowia lipolytica type strain W29 using genome-scale metabolic modelling

Yarrowia lipolytica is a widely-used chassis cell in biotechnological applications. It has recently gained extensive research interest owing to its extraordinary ability of producing industrially valuable biochemicals from a variety of carbon sources. Genome-scale metabolic models (GSMMs) enable analyses of cellular metabolism for engineering various industrial hosts. In the present study, we developed a high-quality GSMM iYli21 for Y. lipolytica type strain W29 by extensive manual curation with Biolog experimental data. The model showed a high accuracy of 85.7% in predicting nutrient utilization. Transcriptomics data were integrated to delineate cellular metabolism of utilizing six individual metabolites as sole carbon sources. Comparisons showed that 302 reactions were commonly used, including those from TCA cycle, oxidative phosphorylation, and purine metabolism for energy and material supply. Whereas glycolytic reactions were employed only when glucose and glycerol used as sole carbon sources, gluconeogenesis and fatty acid oxidation reactions were specifically employed when fatty acid, alkane and glycerolipid were the sole carbon sources. Further test of 46 substrates for generating 5 products showed that hexanoate outcompeted other compounds in terms of maximum theoretical yield owing to the lowest carbon loss for energy supply. This newly generated model iYli21 will be a valuable tool in dissecting metabolic mechanism and guiding metabolic engineering of this important industrial cell factory.

Genome-scale modeling of yeast metabolism: retrospectives and perspectives

Yeasts have been widely used for production of bread, beer and wine, as well as for production of bioethanol, but they have also been designed as cell factories to produce various chemicals, advanced biofuels and recombinant proteins. To systematically understand and rationally engineer yeast metabolism, genome-scale metabolic models (GEMs) have been reconstructed for the model yeast Saccharomyces cerevisiae and nonconventional yeasts. Here, we review the historical development of yeast GEMs together with their recent applications, including metabolic flux prediction, cell factory design, culture condition optimization and multi-yeast comparative analysis. Furthermore, we present an emerging effort, namely the integration of proteome constraints into yeast GEMs, resulting in models with improved performance. At last, we discuss challenges and perspectives on the development of yeast GEMs and the integration of proteome constraints.

The yeastGemMap: A process diagram to assist yeast systems-metabolic studies

Visualization is a key aspect of the analysis of omics data. Although many tools can generate pathway maps for visualization of yeast metabolism, they fail in reconstructing genome-scale metabolic diagrams of compartmentalized metabolism. Here we report on the yeastGemMap, a process diagram of whole yeast metabolism created to assist data analysis in systems-metabolic studies. The map was manually reconstructed with reactions from a compartmentalized genome-scale metabolic model, based on biochemical process diagrams typically found in educational and specialized literature. The yeastGemMap consists of 3815 reactions representing 1150 genes, 2742 metabolites, and 14 compartments. Computational functions for adapting the graphical representation of the map are also reported. These functions modify the visual representation of the map to assist in three systems-metabolic tasks: illustrating reaction networks, interpreting metabolic flux data, and visualizing omics data. The versatility of the yeastGemMap and algorithms to assist visualization of systems-metabolic data was demonstrated in various tasks, including for single lethal reaction evaluation, flux balance analysis, and transcriptomic data analysis. For instance, visual interpretation of metabolic transcriptomes of thermally evolved and parental yeast strains allowed to demonstrate that evolved strains activate a preadaptation response at 30°C, which enabled thermotolerance. A quick interpretation of systems-metabolic data is promoted with yeastGemMap visualizations.

Path to improving the life cycle and quality of genome-scale models of metabolism

Genome-scale models of metabolism (GEMs) are key computational tools for the systems-level study of metabolic networks. Here, we describe the "GEM life cycle," which we subdivide into four stages: inception, maturation, specialization, and amalgamation. We show how different types of GEM reconstruction workflows fit in each stage and proceed to highlight two fundamental bottlenecks for GEM quality improvement: GEM maturation and content removal. We identify common characteristics contributing to increasing quality of maturing GEMs drawing from past independent GEM maturation efforts. We then shed some much-needed light on the latent and unrecognized but pervasive issue of content removal, demonstrating the substantial effects of model pruning on its solution space. Finally, we propose a novel framework for content removal and associated confidence-level assignment which will help guide future GEM development efforts, reduce duplication of effort across groups, potentially aid automated reconstruction platforms, and boost the reproducibility of model development.

Model-guided development of an evolutionarily stable yeast chassis

First-principle metabolic modelling holds potential for designing microbial chassis that are resilient against phenotype reversal due to adaptive mutations. Yet, the theory of model-based chassis design has rarely been put to rigorous experimental test. Here, we report the development of Saccharomyces cerevisiae chassis strains for dicarboxylic acid production using genome-scale metabolic modelling. The chassis strains, albeit geared for higher flux towards succinate, fumarate and malate, do not appreciably secrete these metabolites. As predicted by the model, introducing product-specific TCA cycle disruptions resulted in the secretion of the corresponding acid. Adaptive laboratory evolution further improved production of succinate and fumarate, demonstrating the evolutionary robustness of the engineered cells. In the case of malate, multi-omics analysis revealed a flux bypass at peroxisomal malate dehydrogenase that was missing in the yeast metabolic model. In all three cases, flux balance analysis integrating transcriptomics, proteomics and metabolomics data confirmed the flux re-routing predicted by the model. Taken together, our modelling and experimental results have implications for the computer-aided design of microbial cell factories.

Tolerance to NADH/NAD+ imbalance anticipates aging and anti-aging interventions

Redox couples coordinate cellular function, but the consequences of their imbalances are unclear. This is somewhat associated with the limitations of their experimental quantification. Here we circumvent these difficulties by presenting an approach that characterizes fitness-based tolerance profiles to redox couple imbalances using an in silico representation of metabolism. Focusing on the NADH/NAD⁺ redox couple in yeast, we demonstrate that reductive disequilibria generate metabolic syndromes comparable to those observed in cancer cells. The tolerance of yeast mutants to redox disequilibrium can also explain 30% of the variability in their experimentally measured chronological lifespan. Moreover, by predicting the significance of some metabolites to help stand imbalances, we correctly identify nutrients underlying mechanisms of pathology, lifespan-protecting molecules, or caloric restriction mimetics. Tolerance to redox imbalances becomes, in this way, a sound framework to recognize properties of the aging phenotype while providing a consistent biological rationale to assess anti-aging interventions.

•

We simulate how imbalances in NADH/NAD+ ratio modify cellular metabolic behavior

•

This reveals a mechanism to understand metabolic alterations at low growth rates

•

Tolerance to imbalance explains experimentally measured lifespan in yeast

•

We predict lifespan-protecting metabolites in yeast, animal, and human models

Machine learning applied for metabolic flux-based control of micro-aerated fermentations in bioreactors

Various bio-based processes depend on controlled micro-aerobic conditions to achieve a satisfactory product yield. However, the limiting oxygen concentration varies according to the micro-organism employed, while for industrial applications, there is no cost-effective way of measuring it at low levels. This study proposes a machine learning procedure within a metabolic flux-based control strategy (SUPERSYS_MCU) to address this issue. The control strategy used simulations of a genome-scale metabolic model to generate a surrogate model in the form of an artificial neural network, to be used in a micro-aerobic fermentation strategy (MF-ANN). The meta-model provided setpoints to the controller, allowing adjustment of the inlet air flow to control the oxygen uptake rate. The strategy was evaluated in micro-aerobic batch cultures employing industrial Saccharomyces cerevisiae yeast, with defined medium and glucose as the carbon source, as a case study. The performance of the proposed control scheme was compared with a conventional fermentation and with three previously reported micro-aeration strategies, including respiratory quotient-based control and constant air flow rate. Due to maintenance of the oxidative balance at the anaerobiosis threshold, the MF-ANN provided volumetric ethanol productivity of 4.16 g·L^-1 ·h^-1 and a yield of 0.48 g_ethanol .g_substrate ^-1 , which were higher than the values achieved for the other conditions studied (maximum of 3.4 g·L^-1 ·h^-1 and 0.35-0.40 g_ethanol ·g_substrate ^-1 , respectively). Due to its modular character, the MF-ANN strategy could be adapted to other micro-aerated bioprocesses.

Evaluating accessibility, usability and interoperability of genome-scale metabolic models for diverse yeasts species

Metabolic network reconstructions have become an important tool for probing cellular metabolism in the field of systems biology. They are used as tools for quantitative prediction but also as scaffolds for further knowledge contextualization. The yeast Saccharomyces cerevisiae was one of the first organisms for which a genome-scale metabolic model (GEM) was reconstructed, in 2003, and since then 45 metabolic models have been developed for a wide variety of relevant yeasts species. A systematic evaluation of these models revealed that-despite this long modeling history-the sequential process of tracing model files, setting them up for basic simulation purposes and comparing them across species and even different versions, is still not a generalizable task. These findings call the yeast modeling community to comply to standard practices on model development and sharing in order to make GEMs accessible and useful for a wider public.

Addressing uncertainty in genome-scale metabolic model reconstruction and analysis

The reconstruction and analysis of genome-scale metabolic models constitutes a powerful systems biology approach, with applications ranging from basic understanding of genotype-phenotype mapping to solving biomedical and environmental problems. However, the biological insight obtained from these models is limited by multiple heterogeneous sources of uncertainty, which are often difficult to quantify. Here we review the major sources of uncertainty and survey existing approaches developed for representing and addressing them. A unified formal characterization of these uncertainties through probabilistic approaches and ensemble modeling will facilitate convergence towards consistent reconstruction pipelines, improved data integration algorithms, and more accurate assessment of predictive capacity.

The Metano Modeling Toolbox MMTB: An Intuitive, Web-Based Toolbox Introduced by Two Use Cases

Genome-scale metabolic models are of high interest in a number of different research fields. Flux balance analysis (FBA) and other mathematical methods allow the prediction of the steady-state behavior of metabolic networks under different environmental conditions. However, many existing applications for flux optimizations do not provide a metabolite-centric view on fluxes. Metano is a standalone, open-source toolbox for the analysis and refinement of metabolic models. While flux distributions in metabolic networks are predominantly analyzed from a reaction-centric point of view, the Metano methods of split-ratio analysis and metabolite flux minimization also allow a metabolite-centric view on flux distributions. In addition, we present MMTB (Metano Modeling Toolbox), a web-based toolbox for metabolic modeling including a user-friendly interface to Metano methods. MMTB assists during bottom-up construction of metabolic models by integrating reaction and enzymatic annotation data from different databases. Furthermore, MMTB is especially designed for non-experienced users by providing an intuitive interface to the most commonly used modeling methods and offering novel visualizations. Additionally, MMTB allows users to upload their models, which can in turn be explored and analyzed by the community. We introduce MMTB by two use cases, involving a published model of Corynebacterium glutamicum and a newly created model of Phaeobacter inhibens.

Population-Based Parameter Identification for Dynamical Models of Biological Networks with an Application to Saccharomyces cerevisiae

One of the central elements in systems biology is the interaction between mathematical modeling and measured quantities. Typically, biological phenomena are represented as dynamical systems, and they are further analyzed and comprehended by identifying model parameters using experimental data. However, all model parameters cannot be found by gradient-based optimization methods by fitting the model to the experimental data due to the non-differentiable character of the problem. Here, we present POPI4SB, a Python-based framework for population-based parameter identification of dynamic models in systems biology. The code is built on top of PySCeS that provides an engine to run dynamic simulations. The idea behind the methodology is to provide a set of derivative-free optimization methods that utilize a population of candidate solutions to find a better solution iteratively. Additionally, we propose two surrogate-assisted population-based methods, namely, a combination of a k-nearest-neighbor regressor with the Reversible Differential Evolution and the Evolution of Distribution Algorithm, that speeds up convergence. We present the optimization framework on the example of the well-studied glycolytic pathway in Saccharomyces cerevisiae.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

Metagenomics for taxonomy profiling: tools and approaches

The study of metagenomics is an emerging field that identifies the total genetic materials in an organism along with the set of all genetic materials like deoxyribonucleic acid and ribose nucleic acid, which play a key role with the maintenance of cellular functions. The best part of this technology is that it gives more flexibility to environmental microbiologists to instantly pioneer the immense genetic variability of microbial communities. However, it is intensively complex to identify the suitable sequencing measures of any specific gene that can exclusively indicate the involvement of microbial metagenomes and be able to advance valuable results about these communities. This review provides an overview of the metagenomic advancement that has been advantageous for aggregation of more knowledge about specific genes, microbial communities and its metabolic pathways. More specific drawbacks of metagenomes technology mainly depend on sequence-based analysis. Therefore, this 'targeted based metagenomics' approach will give comprehensive knowledge about the ecological, evolutionary and functional sequence of significantly important genes that naturally exist in living beings either human, animal and microorganisms from distinctive ecosystems.

Fluxer: a web application to compute, analyze and visualize genome-scale metabolic flux networks

Next-generation sequencing has paved the way for the reconstruction of genome-scale metabolic networks as a powerful tool for understanding metabolic circuits in any organism. However, the visualization and extraction of knowledge from these large networks comprising thousands of reactions and metabolites is a current challenge in need of user-friendly tools. Here we present Fluxer (https://fluxer.umbc.edu), a free and open-access novel web application for the computation and visualization of genome-scale metabolic flux networks. Any genome-scale model based on the Systems Biology Markup Language can be uploaded to the tool, which automatically performs Flux Balance Analysis and computes different flux graphs for visualization and analysis. The major metabolic pathways for biomass growth or for biosynthesis of any metabolite can be interactively knocked-out, analyzed and visualized as a spanning tree, dendrogram or complete graph using different layouts. In addition, Fluxer can compute and visualize the k-shortest metabolic paths between any two metabolites or reactions to identify the main metabolic routes between two compounds of interest. The web application includes >80 whole-genome metabolic reconstructions of diverse organisms from bacteria to human, readily available for exploration. Fluxer enables the efficient analysis and visualization of genome-scale metabolic models toward the discovery of key metabolic pathways.

Systematically gap-filling the genome-scale metabolic model of CHO cells

OBJECTIVE

Chinese hamster ovary (CHO) cells are the leading cell factories for producing recombinant proteins in the biopharmaceutical industry. In this regard, constraint-based metabolic models are useful platforms to perform computational analysis of cell metabolism. These models need to be regularly updated in order to include the latest biochemical data of the cells, and to increase their predictive power. Here, we provide an update to iCHO1766, the metabolic model of CHO cells.

RESULTS

We expanded the existing model of Chinese hamster metabolism with the help of four gap-filling approaches, leading to the addition of 773 new reactions and 335 new genes. We incorporated these into an updated genome-scale metabolic network model of CHO cells, named iCHO2101. In this updated model, the number of reactions and pathways capable of carrying flux is substantially increased.

CONCLUSIONS

The present CHO model is an important step towards more complete metabolic models of CHO cells.

Gapsplit: efficient random sampling for non-convex constraint-based models

SUMMARY

Gapsplit generates random samples from convex and non-convex constraint-based models by targeting under-sampled regions of the solution space. Gapsplit provides uniform coverage of linear, mixed-integer and general non-linear models.

AVAILABILITY AND IMPLEMENTATION

Python and Matlab source code are freely available at http://jensenlab.net/tools.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-scale metabolic reconstruction of the non-model yeast Issatchenkia orientalis SD108 and its application to organic acids production

Many platform chemicals can be produced from renewable biomass by microorganisms, with organic acids making up a large fraction. Intolerance to the resulting low pH growth conditions, however, remains a challenge for the industrial production of organic acids by microorganisms. Issatchenkia orientalis SD108 is a promising host for industrial production because it is tolerant to acidic conditions as low as pH 2.0. With the goal to systematically assess the metabolic capabilities of this non-model yeast, we developed a genome-scale metabolic model for I. orientalis SD108 spanning 850 genes, 1826 reactions, and 1702 metabolites. In order to improve the model’s quantitative predictions, organism-specific macromolecular composition and ATP maintenance requirements were determined experimentally and implemented. We examined its network topology, including essential genes and flux coupling analysis and drew comparisons with the Yeast 8.3 model for Saccharomyces cerevisiae. We explored the carbon substrate utilization and examined the organism’s production potential for the industrially-relevant succinic acid, making use of the OptKnock framework to identify gene knockouts which couple production of the targeted chemical to biomass production. The genome-scale metabolic model iIsor850 is a data-supported curated model which can inform genetic interventions for overproduction.

•

Genome-scale metabolic model iIsor850 describes metabolism of I. orientalis SD108.

•

Customized biomass reaction highlights differences with S. cerevisiae.

•

Chemostat data elucidate growth-associated ATP maintenance.

•

Substrate utilization and CRISPR/Cas9 gene knockout phenotypes validate model.

•

Model pinpoints candidate gene deletions coupling succinic acid production to growth.

Reconstruction and analysis of genome-scale metabolic model of weak Crabtree positive yeast Lachancea kluyveri

Lachancea kluyveri, a weak Crabtree positive yeast, has been extensively studied for its unique URC pyrimidine catabolism pathway. It produces more biomass than Saccharomyces cerevisiae due to the underlying weak Crabtree effect and resorts to fermentation only in oxygen limiting conditions that renders it as a suitable industrial host. The yeast also produces ethyl acetate as a major overflow metabolite in aerobic conditions. Here, we report the first genome-scale metabolic model, iPN730, of L. kluyveri comprising of 1235 reactions, 1179 metabolites, and 730 genes distributed in 8 compartments. The in silico viability in different media conditions and the growth characteristics in various carbon sources show good agreement with experimental data. Dynamic flux balance analysis describes the growth dynamics, substrate utilization and product formation kinetics in various oxygen-limited conditions. We have also demonstrated the effect of switching carbon sources on the production of ethyl acetate under varying oxygen uptake rates. A phenotypic phase plane analysis described the energetic cost penalty of ethyl acetate and ethanol production on the specific growth rate of L. kluyveri. We generated the context specific models of L. kluyveri growing on uracil or ammonium salts as the sole nitrogen source. Differential flux calculated using flux variability analysis helped us in highlighting pathways like purine, histidine, riboflavin and pyrimidine metabolism associated with uracil degradation. The genome-scale metabolic construction of L. kluyveri will provide a better understanding of metabolism behind ethyl acetate production as well as uracil catabolism (pyrimidine degradation) pathway. iPN730 is an addition to genome-scale metabolic models of non-conventional yeasts that will facilitate system-wide omics analysis to understand fungal metabolic diversity.

Curation and Analysis of a Saccharomyces cerevisiae Genome-Scale Metabolic Model for Predicting Production of Sensory Impact Molecules under Enological Conditions

One approach for elucidating strain-to-strain metabolic differences is the use of genome-scale metabolic models (GSMMs). To date GSMMs have not focused on the industrially important area of flavor production and, as such; do not cover all the pathways relevant to flavor formation in yeast. Moreover, current models for Saccharomyces cerevisiae generally focus on carbon-limited and/or aerobic systems, which is not pertinent to enological conditions. Here, we curate a GSMM (iWS902) to expand on the existing Ehrlich pathway and ester formation pathways central to aroma formation in industrial winemaking, in addition to the existing sulfur metabolism and medium-chain fatty acid (MCFA) pathways that also contribute to production of sensory impact molecules. After validating the model using experimental data, we predict key differences in metabolism for a strain (EC 1118) in two distinct growth conditions, including differences for aroma impact molecules such as acetic acid, tryptophol, and hydrogen sulfide. Additionally, we propose novel targets for metabolic engineering for aroma profile modifications employing flux variability analysis with the expanded GSMM. The model provides mechanistic insights into the key metabolic pathways underlying aroma formation during alcoholic fermentation and provides a potential framework to contribute to new strategies to optimize the aroma of wines.

Genetic buffering and potentiation in metabolism

Cells adjust their metabolism in response to mutations, but how this reprogramming depends on the genetic context is not well known. Specifically, the absence of individual enzymes can affect reprogramming, and thus the impact of mutations in cell growth. Here, we examine this issue with an in silico model of Saccharomyces cerevisiae's metabolism. By quantifying the variability in the growth rate of 10000 different mutant metabolisms that accumulated changes in their reaction fluxes, in the presence, or absence, of a specific enzyme, we distinguish a subset of modifier genes serving as buffers or potentiators of variability. We notice that the most potent modifiers refer to the glycolysis pathway and that, more broadly, they show strong pleiotropy and epistasis. Moreover, the evidence that this subset depends on the specific growing condition strengthens its systemic underpinning, a feature only observed before in a toy model of a gene-regulatory network. Some of these enzymes also modulate the effect that biochemical noise and environmental fluctuations produce in growth. Thus, the reorganization of metabolism induced by mutations has not only direct physiological implications but also transforms the influence that other mutations have on growth. This is a general result with implications in the development of cancer therapies based on metabolic inhibitors.

Uptake of exogenous serine is important to maintain sphingolipid homeostasis in Saccharomyces cerevisiae

Sphingolipids are abundant and essential molecules in eukaryotes that have crucial functions as signaling molecules and as membrane components. Sphingolipid biosynthesis starts in the endoplasmic reticulum with the condensation of serine and palmitoyl-CoA. Sphingolipid biosynthesis is highly regulated to maintain sphingolipid homeostasis. Even though, serine is an essential component of the sphingolipid biosynthesis pathway, its role in maintaining sphingolipid homeostasis has not been precisely studied. Here we show that serine uptake is an important factor for the regulation of sphingolipid biosynthesis in Saccharomyces cerevisiae. Using genetic experiments, we find the broad-specificity amino acid permease Gnp1 to be important for serine uptake. We confirm these results with serine uptake assays in gnp1Δ cells. We further show that uptake of exogenous serine by Gnp1 is important to maintain cellular serine levels and observe a specific connection between serine uptake and the first step of sphingolipid biosynthesis. Using mass spectrometry-based flux analysis, we further observed imported serine as the main source for de novo sphingolipid biosynthesis. Our results demonstrate that yeast cells preferentially use the uptake of exogenous serine to regulate sphingolipid biosynthesis. Our study can also be a starting point to analyze the role of serine uptake in mammalian sphingolipid metabolism.

Sphingolipids (SPs) are membrane lipids globally required for eukaryotic life. In contrast to other lipid classes, SPs cannot be stored in the cell and therefore their levels have to be tightly regulated. Failure to maintain sphingolipid homeostasis can result in pathologies including neurodegeneration, childhood asthma and cancer. However, we are only starting to understand how SP biosynthesis is adjusted according to need. In this study, we use genetic and biochemical methods to show that the uptake of exogenous serine is necessary to maintain SP homeostasis in Saccharomyces cerevisiae. Serine is one of the precursors of long chain bases in cells, the first intermediate of SP metabolism. Our results suggest that the uptake of serine is directly coupled to SP biosynthesis at ER-plasma membrane contact sites. Overall, our study identifies serine uptake as a novel regulatory factor of SP homeostasis. While we use yeast as a discovery tool, these results also provide valuable insights into mammalian SP biology especially under pathological conditions.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

Quantitative Connection between Cell Size and Growth Rate by Phospholipid Metabolism

The processes involved in cell growth are extremely complicated even for a single cell organism such as Escherichia coli, while the relationship between growth rate and cell size is simple. We aimed to reveal the systematic link between them from the aspect of the genome-scale metabolic network. Since the growth rate reflects metabolic rates of bacteria and the cell size relates to phospholipid synthesis, a part of bacterial metabolic networks, we calculated the cell length from the cardiolipin synthesis rate, where the cardiolipin synthesis reaction is able to represent the phospholipid metabolism of Escherichia coli in the exponential growth phase. Combined with the flux balance analysis, it enables us to predict cell length and to examine the quantitative relationship between cell length and growth rate. By simulating bacteria growing in various nutrient media with the flux balance analysis and calculating the corresponding cell length, we found that the increase of the synthesis rate of phospholipid, the cell width, and the protein fraction in membranes caused the increase of cell length with growth rate. Different tendencies of phospholipid synthesis rate changing with growth rate result in different relationships between cell length and growth rate. The effects of gene deletions on cell size and growth rate are also examined. Knocking out the genes, such as Δ tktA, Δ tktB, Δ yqaB, Δ pgm, and Δ cysQ, affects growth rate largely while affecting cell length slightly. Results of this method are in good agreement with experiments.

CAMND: Comparative analysis of metabolic network decomposition based on previous and two new criteria, a web based application

Metabolic networks can model the behavior of metabolism in the cell. Since analyzing the whole metabolic networks is not easy, network modulation is an important issue to be investigated. Decomposing metabolic networks is a strategy to obtain better insight into metabolic functions. Additionally, decomposing these networks facilitates using computational methods, which are very slow when applied to the original genome-scale network. Several methods have been proposed for decomposing of the metabolic network. Therefore, it is necessary to evaluate these methods by suitable criteria. In this study, we introduce a web server package for decomposing of metabolic networks with 10 different methods, 9 datasets and the ability of computing 12 criteria, to evaluate and compare the results of different methods using ten previously defined and two new criteria which are based on Chebi ontology and Co-expression_of_Enzymes information. This package visualizes the obtained modules via "gephi" software. The ability of this package is that the user can examine whether two metabolites or reactions are in the same module or not. The functionality of the package can be easily extended to include new datasets and criteria. It also has the ability to compare the results of novel methods with the results of previously developed methods. The package is implemented in python and is available at http://eslahchilab.ir/softwares/dmn.

Evolution of C4 photosynthesis predicted by constraint-based modelling

Constraint-based modelling (CBM) is a powerful tool for the analysis of evolutionary trajectories. Evolution, especially evolution in the distant past, is not easily accessible to laboratory experimentation. Modelling can provide a window into evolutionary processes by allowing the examination of selective pressures which lead to particular optimal solutions in the model. To study the evolution of C4 photosynthesis from a ground state of C3 photosynthesis, we initially construct a C3 model. After duplication into two cells to reflect typical C4 leaf architecture, we allow the model to predict the optimal metabolic solution under various conditions. The model thus identifies resource limitation in conjunction with high photorespiratory flux as a selective pressure relevant to the evolution of C4. It also predicts that light availability and distribution play a role in guiding the evolutionary choice of possible decarboxylation enzymes. The data shows evolutionary CBM in eukaryotes predicts molecular evolution with precision.