Reconstruction of a catalogue of genome-scale metabolic models with enzymatic constraints using GECKO 2.0

New insights in metabolism modelling to decipher plant-microbe interactions

Plant disease outbreaks, exacerbated by climate change, threaten food security and environmental sustainability world-wide. Plants interact with a wide range of microorganisms. The quest for resilient agriculture requires a deep insight into the molecular and ecological interplays between plants and their associated microbial communities. Omics methods, by profiling entire molecular sets, have shed light on these complex interactions. Nonetheless, deciphering the relationships among thousands of molecular components remains a formidable challenge, and studies that integrate these components into cohesive biological networks involving plants and associated microbes are still limited. Systems biology has the potential to predict the effects of biotic and abiotic perturbations on these networks. It is therefore a promising framework for addressing the full complexity of plant-microbiome interactions.

Precise metabolic modeling in post-omics era: accomplishments and perspectives

Microbes have been extensively utilized for their sustainable and scalable properties in synthesizing desired bio-products. However, insufficient knowledge about intracellular metabolism has impeded further microbial applications. The genome-scale metabolic models (GEMs) play a pivotal role in facilitating a global understanding of cellular metabolic mechanisms. These models enable rational modification by exploring metabolic pathways and predicting potential targets in microorganisms, enabling precise cell regulation without experimental costs. Nonetheless, simplified GEM only considers genome information and network stoichiometry while neglecting other important bio-information, such as enzyme functions, thermodynamic properties, and kinetic parameters. Consequently, uncertainties persist particularly when predicting microbial behaviors in complex and fluctuant systems. The advent of the omics era with its massive quantification of genes, proteins, and metabolites under various conditions has led to the flourishing of multi-constrained models and updated algorithms with improved predicting power and broadened dimension. Meanwhile, machine learning (ML) has demonstrated exceptional analytical and predictive capacities when applied to training sets of biological big data. Incorporating the discriminant strength of ML with GEM facilitates mechanistic modeling efficiency and improves predictive accuracy. This paper provides an overview of research innovations in the GEM, including multi-constrained modeling, analytical approaches, and the latest applications of ML, which may contribute comprehensive knowledge toward genetic refinement, strain development, and yield enhancement for a broad range of biomolecules.

Coarse-Grained Modeling Elucidates Differential Metabolism of Saccharomyces cerevisiae under Varied Nutrient Limitations

Microorganisms such as Saccharomyces cerevisiae have a native ability to adapt their metabolism to varying nutrient conditions. Understanding their responses to nutrient limitations is critical for decoding cellular physiology and designing strategies for metabolic engineering. While the influence of carbon availability on yeast metabolism has been extensively studied, the role of nitrogen availability remains relatively underexplored. In this study, we utilized a coarse-grained kinetic model to systematically analyze and compare the effects of carbon and nitrogen limitations on yeast metabolism. Our model successfully revealed the differential metabolic characteristics of S. cerevisiae under carbon- and nitrogen-limited chemostat conditions. It also highlighted the significance of protein activity regulation at varying carbon-to-nitrogen ratios, and elucidated distinct strategies employed to maintain ATP homeostasis. This study provides a computational tool for investigating yeast physiology under nutrient limitations and offers quantitative and mechanistic insights into yeast metabolism.

The Dawn of High-Throughput and Genome-Scale Kinetic Modeling: Recent Advances and Future Directions

Researchers have invested much effort into developing kinetic models due to their ability to capture dynamic behaviors, transient states, and regulatory mechanisms of metabolism, providing a detailed and realistic representation of cellular processes. Historically, the requirements for detailed parametrization and significant computational resources created barriers to their development and adoption for high-throughput studies. However, recent advancements, including the integration of machine learning with mechanistic metabolic models, the development of novel kinetic parameter databases, and the use of tailor-made parametrization strategies, are reshaping the field of kinetic modeling. In this Review, we discuss these developments and offer future directions, highlighting the potential of these advances to drive progress in systems and synthetic biology, metabolic engineering, and medical research at an unprecedented scale and pace.

Influence of Amino Acids on Quorum Sensing-Related Pathways in Pseudomonas aeruginosa PAO1: Insights from the GEM iJD1249

BACKGROUND/OBJECTIVES

Amino acids (AAs) play a critical role in diseases such as cystic fibrosis where Pseudomonas aeruginosa PAO1 adapts its metabolism in response to host-derived nutrients. The adaptation influences virulence and complicates antibiotic treatment mainly for the antimicrobial resistance context. D- and L-AAs have been analyzed for their impact on quorum sensing (QS), a mechanism that regulates virulence factors. This research aimed to reconstruct the genome-scale metabolic model (GEM) of P. aeruginosa PAO1 to investigate the metabolic roles of D- and L-AAs in QS-related pathways.

METHODS

The updated GEM, iJD1249, was reconstructed by using protocols to integrate data from previous models and refined with well-standardized in silico media (LB, M9, and SCFM) to improve flux balance analysis accuracy. The model was used to explore the metabolic impact of D-Met, D-Ala, D-Glu, D-Ser, L-His, L-Glu, L-Arg, and L-Ornithine (L-Orn) at 5 and 50 mM in QS-related pathways, focusing on the effects on bacterial growth and carbon flux distributions.

RESULTS

Among the tested AAs, D-Met was the only one that did not enhance the growth rate of P. aeruginosa PAO1, while L-Arg and L-Orn increased fluxes in the L-methionine biosynthesis pathway, influencing the metH gene. These findings suggest a differential metabolic role for D-and L-AAs in QS-related pathways.

CONCLUSIONS

Our results shed some light on the metabolic impact of AAs on QS-related pathways and their potential role in P. aeruginosa virulence. Future studies should assess D-Met as a potential adjuvant in antimicrobial strategies, optimizing the concentration in combination with antibiotics to maximize its therapeutic effectiveness.

Gene network centrality analysis identifies key regulators coordinating day-night metabolic transitions in Synechococcus elongatus PCC 7942 despite limited accuracy in predicting direct regulator-gene interactions

Synechococcus elongatus PCC 7942 is a model organism for studying circadian regulation and bioproduction, where precise temporal control of metabolism significantly impacts photosynthetic efficiency and CO2-to-bioproduct conversion. Despite extensive research on core clock components, our understanding of the broader regulatory network orchestrating genome-wide metabolic transitions remains incomplete. We address this gap by applying machine learning tools and network analysis to investigate the transcriptional architecture governing circadian-controlled gene expression. While our approach showed moderate accuracy in predicting individual transcription factor-gene interactions - a common challenge with real expression data - network-level topological analysis successfully revealed the organizational principles of circadian regulation. Our analysis identified distinct regulatory modules coordinating day-night metabolic transitions, with photosynthesis and carbon/nitrogen metabolism controlled by day-phase regulators, while nighttime modules orchestrate glycogen mobilization and redox metabolism. Through network centrality analysis, we identified potentially significant but previously understudied transcriptional regulators: HimA as a putative DNA architecture regulator, and TetR and SrrB as potential coordinators of nighttime metabolism, working alongside established global regulators RpaA and RpaB. This work demonstrates how network-level analysis can extract biologically meaningful insights despite limitations in predicting direct regulatory interactions. The regulatory principles uncovered here advance our understanding of how cyanobacteria coordinate complex metabolic transitions and may inform metabolic engineering strategies for enhanced photosynthetic bioproduction from CO2.

Comprehensive evaluation of the capacities of microbial cell factories

Systems metabolic engineering is facilitating the development of high-performing microbial cell factories for producing chemicals and materials. However, constructing an efficient microbial cell factory still requires exploring and selecting various host strains, as well as identifying the best-suited metabolic engineering strategies, which demand significant time, effort, and costs. Here, we comprehensively evaluate the capacities of various microbial cell factories and propose strategies for systems metabolic engineering steps, including host strain selection, metabolic pathway reconstruction, and metabolic flux optimization. We analyze the metabolic capacities of five representative industrial microorganisms as cell factories for the production of 235 different bio-based chemicals and suggest the most suitable host strain for the corresponding chemical production. To improve the innate metabolic capacity by constructing more efficient metabolic pathways, heterologous metabolic reactions, and cofactor exchanges are systematically analyzed. Additionally, we present metabolic engineering strategies, which include up- and down-regulation target reactions, for the improved production of chemicals. Altogether, this study will serve as a comprehensive resource for the systems metabolic engineering of microorganisms in the bio-based production of chemicals.

Protein-constrained models pinpoints the role of underground metabolism in robustness of metabolic phenotypes

Integrating enzyme parameters into constraint-based models have significantly improved the prediction of physiological and molecular traits. To further improve these models, we integrated promiscuous enzyme activities that jointly comprise the so-called underground metabolism by developing the CORAL toolbox, which increases the resolution of modeled enzyme resource allocation. Applying CORAL to a protein-constrained model of Escherichia coli revealed that underground metabolism resulted in larger flexibility of metabolic fluxes and enzyme usage. Simulating metabolic defects where the main activity of a promiscuous enzyme was blocked but promiscuous activities remained functional showed a small enzyme redistribution to the side activities. Further, blocking pairs of main activities showed that non-promiscuous enzymes exhibited larger impact on growth than promiscuous enzymes. These simulations showed that promiscuous enzymes can compensate for these defects, in line with experimental evidence. Together, our results indicated that promiscuous enzyme activities are vital to maintain robust metabolic function and growth.

NEXT-FBA: A hybrid stoichiometric/data-driven approach to improve intracellular flux predictions

Genome-scale metabolic models (GEMs) have been widely utilized to understand cellular metabolism. The application of GEMs has been advanced by computational methods that enable the prediction and analysis of intracellular metabolic states. However, the accuracy and biological relevance of these predictions often suffer from the many degrees of freedom and scarcity of available data to constrain the models adequately. Here, we introduce Neural-net EXtracellular Trained Flux Balance Analysis, (NEXT-FBA), a novel computational methodology that addresses these limitations by utilizing exometabolomic data to derive biologically relevant constraints for intracellular fluxes in GEMs. We achieve this by training artificial neural networks (ANNs) with exometabolomic data from Chinese hamster ovary (CHO) cells and correlating it with 13C-labeled intracellular fluxomic data. By capturing the underlying relationships between exometabolomics and cell metabolism, NEXT-FBA predicts upper and lower bounds for intracellular reaction fluxes to constrain GEMs. We demonstrate the efficacy of NEXT-FBA across several validation experiments, where it outperforms existing methods in predicting intracellular flux distributions that align closely with experimental observations. Furthermore, a case study demonstrates how NEXT-FBA can guide bioprocess optimization by identifying key metabolic shifts and refining flux predictions to yield actionable process and metabolic engineering targets. Overall, NEXT-FBA aims to improve the accuracy and biological relevance of intracellular flux predictions in metabolic modelling, with minimal input data requirements for pre-trained models.

Computational biology predicts metabolic engineering targets for increased production of 103 valuable chemicals in yeast

Development of efficient cell factories that can compete with traditional chemical production processes is complex and generally driven by case-specific strategies, based on the product and microbial host of interest. Despite major advancements in the field of metabolic modeling in recent years, prediction of genetic modifications for increased production remains challenging. Here, we present a computational pipeline that leverages the concept of protein limitations in metabolism for prediction of optimal combinations of gene engineering targets for enhanced chemical bioproduction. We used our pipeline for prediction of engineering targets for 103 different chemicals using Saccharomyces cerevisiae as a host. Furthermore, we identified sets of gene targets predicted for groups of multiple chemicals, suggesting the possibility of rational model-driven design of platform strains for diversified chemical production.

CatPred: a comprehensive framework for deep learning in vitro enzyme kinetic parameters

Estimation of enzymatic activities still heavily relies on experimental assays, which can be cost and time-intensive. We present CatPred, a deep learning framework for predicting in vitro enzyme kinetic parameters, including turnover numbers (kcat), Michaelis constants (Km), and inhibition constants (Ki). CatPred addresses key challenges such as the lack of standardized datasets, performance evaluation on enzyme sequences that are dissimilar to those used during training, and model uncertainty quantification. We explore diverse learning architectures and feature representations, including pretrained protein language models and three-dimensional structural features, to enable robust predictions. CatPred provides accurate predictions with query-specific uncertainty estimates, with lower predicted variances correlating with higher accuracy. Pretrained protein language model features particularly enhance performance on out-of-distribution samples. CatPred also introduces benchmark datasets with extensive coverage (~23 k, 41 k, and 12 k data points for kcat, Km, and Ki respectively). Our framework performs competitively with existing methods while offering reliable uncertainty quantification.

Modeling for understanding and engineering metabolism

Metabolism is at the core of all functions of living cells as it provides Gibbs free energy and building blocks for synthesis of macromolecules, which are necessary for structures, growth, and proliferation. Metabolism is a complex network composed of thousands of reactions catalyzed by enzymes involving many co-factors and metabolites. Traditionally it has been difficult to study metabolism as a whole network and most traditional efforts were therefore focused on specific metabolic pathways, enzymes, and metabolites. By using engineering principles of mathematical modeling to analyze and study metabolism, as well as engineer it, that is, design and build, new metabolic features, it is possible to gain many new fundamental insights as well as applications in biotechnology. Here, we present the history and basic principles of engineering metabolism, as well as the newest developments in the field. We are using examples of applications in: (1) production of protein pharmaceuticals and chemicals; (2) basic studies of metabolism; and (3) impacting health care. We will end by discussing how engineering metabolism can benefit from advances in artificial intelligence (AI)-based models.

COBREXA 2: tidy and scalable construction of complex metabolic models

SUMMARY

Constraint-based metabolic models offer a scalable framework to investigate biological systems using optimality principles. Construction and simulation of detailed models that utilize multiple kinds of constraint systems pose a significant coding overhead, complicating implementation of new types of analyses. We present an improved version of the constraint-based metabolic modeling package COBREXA, which utilizes a hierarchical model construction framework that decouples the implemented analysis algorithms into independent, yet re-combinable, building blocks. By removing the need to re-implement modeling components, assembly of complex metabolic models is simplified, which we demonstrate on use-cases of resource-balanced models, and enzyme-constrained flux balance models of interacting bacterial communities. Notably, these models show improved predictive capabilities in both monoculture and community settings. In perspective, the re-usable model-building components in COBREXA 2 provide a sustainable way to handle increasingly complex models in constraint-based modeling.

AVAILABILITY AND IMPLEMENTATION

COBREXA 2 is available from https://github.com/COBREXA/COBREXA.jl, and from Julia package repositories. COBREXA 2 works on all major operating systems and computer architectures. Documentation is available at https://cobrexa.github.io/COBREXA.jl/.

Artificial design of the genome: from sequences to the 3D structure of chromosomes

Genome design is the foundation of genome synthesis, which provides a new platform for deepening our understanding of biological systems by exploring the fundamental components and structure of the genome. Artificial genome designs can endow unnatural genomes with desired functions. We provide a comprehensive overview of genome design principles ranging from DNA sequences to the 3D structure of chromosomes. Furthermore, we highlight applications of genome design in gene expression, genome structure, genome function, and biocontainment, and discuss the potential of artificial intelligence (AI) in genome design.

Genome-scale models in human metabologenomics

Metabologenomics integrates metabolomics with other omics data types to comprehensively study the genetic and environmental factors that influence metabolism. These multi-omics data can be incorporated into genome-scale metabolic models (GEMs), which are highly curated knowledge bases that explicitly account for genes, transcripts, proteins and metabolites. By including all known biochemical reactions catalysed by enzymes and transporters encoded in the human genome, GEMs analyse and predict the behaviour of complex metabolic networks. Continued advancements to the scale and scope of GEMs - from cells and tissues to microbiomes and the whole body - have helped to design effective treatments and develop better diagnostic tools for metabolic diseases. Furthermore, increasing amounts of multi-omics data are incorporated into GEMs to better identify the underlying mechanisms, biomarkers and potential drug targets of metabolic diseases.

Metabolic modeling identifies determinants of thermal growth responses in Arabidopsis thaliana

Temperature is a critical environmental factor affecting nearly all plant processes, including growth, development, and yield. Yet, despite decades of research, we lack the ability to predict plant performance at different temperatures, limiting the development of climate-resilient crops. Further, there is a pressing need to bridge the gap between the prediction of physiological and molecular traits to improve our understanding and manipulation of plant temperature responses. Here, we developed the first enzyme-constrained model of Arabidopsis thaliana's metabolism, facilitating predictions of growth-related phenotypes at different temperatures. We showed that the model can be employed for in silico identification of genes that affect plant growth at suboptimal growth temperature. Using mutant lines, we validated the genes predicted to affect plant growth, demonstrating the potential of metabolic modeling in accurately predicting plant thermal responses. The temperature-dependent enzyme-constrained metabolic model provides a template that can be used for developing sophisticated strategies to engineer climate-resilient crops.

Thermodynamic tools for more efficient biotechnological processes: an example in poly-(3-hydroxybutyrate) production from carbon monoxide

Modern biotechnology requires the integration of several disciplines, with thermodynamics being a crucial one. Experimental approaches frequently used in biotechnology, such as rewiring of metabolic networks or culturing of micro-organisms in engineered environments, can benefit from the application of thermodynamic tools. In this paper, we provide an overview of several thermodynamic tools that are useful for the design and optimization of biotechnological processes, and we demonstrate their potential application in the production of poly-(3-hydroxybutyrate) (PHB) from carbon monoxide (CO). We discuss how these tools can aid in the design of metabolic engineering strategies, the calculation of expected yields, the assessment of the thermodynamic feasibility of the targeted conversions, the identification of potential thermodynamic bottlenecks, and the selection of genetic engineering targets. Although we illustrate these tools using the specific example of PHB production from CO, they can be applied to other substrates and products.

Harnessing the optimization of enzyme catalytic rates in engineering of metabolic phenotypes

The increasing availability of enzyme turnover number measurements from experiments and of turnover number predictions from deep learning models prompts the use of these enzyme parameters in precise metabolic engineering. Yet, there is no computational approach that allows the prediction of metabolic engineering strategies that rely on the modification of turnover numbers. It is also unclear if modifications of turnover numbers without alterations in the host's transcriptional regulatory machinery suffice to increase the production of chemicals of interest. Here, we present a constraint-based modeling approach, termed Overcoming Kinetic rate Obstacles (OKO), that uses enzyme-constrained metabolic models to predict in silico strategies to increase the production of a given chemical, while ensuring specified cell growth. We demonstrate that the application of OKO to enzyme-constrained metabolic models of Escherichia coli and Saccharomyces cerevisiae results in strategies that can at least double the production of over 40 compounds with little penalty to growth. Interestingly, we show that the overproduction of compounds of interest does not entail only an increase in the values of turnover numbers. Lastly, we demonstrate that a refinement of OKO, allowing also for manipulation of enzyme abundance, facilitates the usage of the available compendia and deep learning models of turnover numbers in the design of precise metabolic engineering strategies. Our results expand the usage of genome-scale metabolic models toward the identification of targets for protein engineering, allowing their direct usage in the generation of innovative metabolic engineering designs for various biotechnological applications.

Integrating Molecular Perspectives: Strategies for Comprehensive Multi-Omics Integrative Data Analysis and Machine Learning Applications in Transcriptomics, Proteomics, and Metabolomics

With the advent of high-throughput technologies, the field of omics has made significant strides in characterizing biological systems at various levels of complexity. Transcriptomics, proteomics, and metabolomics are the three most widely used omics technologies, each providing unique insights into different layers of a biological system. However, analyzing each omics data set separately may not provide a comprehensive understanding of the subject under study. Therefore, integrating multi-omics data has become increasingly important in bioinformatics research. In this article, we review strategies for integrating transcriptomics, proteomics, and metabolomics data, including co-expression analysis, metabolite-gene networks, constraint-based models, pathway enrichment analysis, and interactome analysis. We discuss combined omics integration approaches, correlation-based strategies, and machine learning techniques that utilize one or more types of omics data. By presenting these methods, we aim to provide researchers with a better understanding of how to integrate omics data to gain a more comprehensive view of a biological system, facilitating the identification of complex patterns and interactions that might be missed by single-omics analyses.

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.

Integration of proteomic data with genome-scale metabolic models: A methodological overview

The integration of proteomics data with constraint-based reconstruction and analysis (COBRA) models plays a pivotal role in understanding the relationship between genotype and phenotype and bridges the gap between genome-level phenomena and functional adaptations. Integrating a generic genome-scale model with information on proteins enables generation of a context-specific metabolic model which improves the accuracy of model prediction. This review explores methodologies for incorporating proteomics data into genome-scale models. Available methods are grouped into four distinct categories based on their approach to integrate proteomics data and their depth of modeling. Within each category section various methods are introduced in chronological order of publication demonstrating the progress of this field. Furthermore, challenges and potential solutions to further progress are outlined, including the limited availability of appropriate in vitro data, experimental enzyme turnover rates, and the trade-off between model accuracy, computational tractability, and data scarcity. In conclusion, methods employing simpler approaches demand fewer kinetic and omics data, consequently leading to a less complex mathematical problem and reduced computational expenses. On the other hand, approaches that delve deeper into cellular mechanisms and aim to create detailed mathematical models necessitate more extensive kinetic and omics data, resulting in a more complex and computationally demanding problem. However, in some cases, this increased cost can be justified by the potential for more precise predictions.

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.

Design and construction of artificial metabolic pathways for the bioproduction of useful compounds

To efficiently produce useful compounds using biological cells, it is essential to optimally design all metabolic reactions and pathways, including not only the flow of carbon within the cell but also the production and consumption of energy and the balance of oxidation-reduction. Computational scientific methods are effective for the rational design of metabolic pathways and the optimization of metabolic fluxes. Based on this blueprint, it is crucial to accurately construct the cell, test and analyze whether it conforms to the design, and learn from the results to redesign the system in an effective cycle. This review introduces essential metabolic design techniques in synthetic biology and discusses the potential of using plant cells or plant genes effectively in synthetic biology for the production of useful compounds.

Relieving metabolic burden to improve robustness and bioproduction by industrial microorganisms

Metabolic burden is defined by the influence of genetic manipulation and environmental perturbations on the distribution of cellular resources. The rewiring of microbial metabolism for bio-based chemical production often leads to a metabolic burden, followed by adverse physiological effects, such as impaired cell growth and low product yields. Alleviating the burden imposed by undesirable metabolic changes has become an increasingly attractive approach for constructing robust microbial cell factories. In this review, we provide a brief overview of metabolic burden engineering, focusing specifically on recent developments and strategies for diminishing the burden while improving robustness and yield. A variety of examples are presented to showcase the promise of metabolic burden engineering in facilitating the design and construction of robust microbial cell factories. Finally, challenges and limitations encountered in metabolic burden engineering are discussed.

ECMpy 2.0: A Python package for automated construction and analysis of enzyme-constrained models

Genome-scale metabolic models (GEMs) have been widely employed to predict microorganism behaviors. However, GEMs only consider stoichiometric constraints, leading to a linear increase in simulated growth and product yields as substrate uptake rates rise. This divergence from experimental measurements prompted the creation of enzyme-constrained models (ecModels) for various species, successfully enhancing chemical production. Building upon studies that allocate macromolecule resources, we developed a Python-based workflow (ECMpy) that constructs an enzyme-constrained model. This involves directly imposing an enzyme amount constraint in GEM and accounting for protein subunit composition in reactions. However, this procedure demands manual collection of enzyme kinetic parameter information and subunit composition details, making it rather user-unfriendly. In this work, we've enhanced the ECMpy toolbox to version 2.0, broadening its scope to automatically generate ecGEMs for a wider array of organisms. ECMpy 2.0 automates the retrieval of enzyme kinetic parameters and employs machine learning for predicting these parameters, which significantly enhances parameter coverage. Additionally, ECMpy 2.0 introduces common analytical and visualization features for ecModels, rendering computational results more user accessible. Furthermore, ECMpy 2.0 seamlessly integrates three published algorithms that exploit ecModels to uncover potential targets for metabolic engineering. ECMpy 2.0 is available at https://github.com/tibbdc/ECMpy or as a pip package (https://pypi.org/project/ECMpy/).

Cell factory design with advanced metabolic modelling empowered by artificial intelligence

Advances in synthetic biology and artificial intelligence (AI) have provided new opportunities for modern biotechnology. High-performance cell factories, the backbone of industrial biotechnology, are ultimately responsible for determining whether a bio-based product succeeds or fails in the fierce competition with petroleum-based products. To date, one of the greatest challenges in synthetic biology is the creation of high-performance cell factories in a consistent and efficient manner. As so-called white-box models, numerous metabolic network models have been developed and used in computational strain design. Moreover, great progress has been made in AI-powered strain engineering in recent years. Both approaches have advantages and disadvantages. Therefore, the deep integration of AI with metabolic models is crucial for the construction of superior cell factories with higher titres, yields and production rates. The detailed applications of the latest advanced metabolic models and AI in computational strain design are summarized in this review. Additionally, approaches for the deep integration of AI and metabolic models are discussed. It is anticipated that advanced mechanistic metabolic models powered by AI will pave the way for the efficient construction of powerful industrial chassis strains in the coming years.

Molecular dynamics of DNA repair and carcinogen interaction: Implications for cancer initiation, progression, and therapeutic strategies

The integrity of the genetic material in human cells is continuously challenged by environmental agents and endogenous stresses. Among these, environmental carcinogens are pivotal in initiating complex DNA lesions that can lead to malignant transformations if not properly repaired. This review synthesizes current knowledge on the molecular dynamics of DNA repair mechanisms and their interplay with various environmental carcinogens, providing a comprehensive overview of how these interactions contribute to cancer initiation and progression. We examine key DNA repair pathways including base excision repair, nucleotide excision repair, and double-strand break repair and their regulatory networks, highlighting how defects in these pathways can exacerbate carcinogen-induced damage. Further, we discuss how understanding these molecular interactions offers novel insights into potential therapeutic strategies. This includes leveraging synthetic lethality concepts and designing targeted therapies that exploit specific DNA repair vulnerabilities in cancer cells. By integrating recent advances in molecular biology, genetics, and oncology, this review aims to illuminate the complex landscape of DNA repair and carcinogen-induced carcinogenesis, setting the stage for future research and therapeutic innovations.

Current State, Challenges, and Opportunities in Genome-Scale Resource Allocation Models: A Mathematical Perspective

Stoichiometric genome-scale metabolic models (generally abbreviated GSM, GSMM, or GEM) have had many applications in exploring phenotypes and guiding metabolic engineering interventions. Nevertheless, these models and predictions thereof can become limited as they do not directly account for protein cost, enzyme kinetics, and cell surface or volume proteome limitations. Lack of such mechanistic detail could lead to overly optimistic predictions and engineered strains. Initial efforts to correct these deficiencies were by the application of precursor tools for GSMs, such as flux balance analysis with molecular crowding. In the past decade, several frameworks have been introduced to incorporate proteome-related limitations using a genome-scale stoichiometric model as the reconstruction basis, which herein are called resource allocation models (RAMs). This review provides a broad overview of representative or commonly used existing RAM frameworks. This review discusses increasingly complex models, beginning with stoichiometric models to precursor to RAM frameworks to existing RAM frameworks. RAM frameworks are broadly divided into two categories: coarse-grained and fine-grained, with different strengths and challenges. Discussion includes pinpointing their utility, data needs, highlighting framework strengths and limitations, and appropriateness to various research endeavors, largely through contrasting their mathematical frameworks. Finally, promising future applications of RAMs are discussed.

Pan-Cancer, Genome-Scale Metabolic Network Analysis of over 10,000 Patients Elucidates Relationship between Metabolism and Survival

Despite the high variability in cancer biology, cancers nevertheless exhibit cohesive hallmarks across multiple cancer types, notably dysregulated metabolism. Metabolism plays a central role in cancer biology, and shifts in metabolic pathways have been linked to tumor aggressiveness and likelihood of response to therapy. We therefore sought to interrogate metabolism across cancer types and understand how intrinsic modes of metabolism vary within and across indications and how they relate to patient prognosis. We used context specific genome-scale metabolic modeling to simulate metabolism across 10,915 patients from 34 cancer types from The Cancer Genome Atlas and the MMRF-COMMPASS study. We found that cancer metabolism clustered into modes characterized by differential glycolysis, oxidative phosphorylation, and growth rate. We also found that the simulated activities of metabolic pathways are intrinsically prognostic across cancer types, especially tumor growth rate, fatty acid biosynthesis, folate metabolism, oxidative phosphorylation, steroid metabolism, and glutathione metabolism. This work shows the prognostic power of individual patient metabolic modeling across multiple cancer types. Additionally, it shows that analyzing large-scale models of cancer metabolism with survival information provides unique insights into underlying relationships across cancer types and suggests how therapies designed for one cancer type may be repurposed for use in others.

Microbiome modeling: a beginner's guide

Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with each other, their environment, and hosts via ecological interactions. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding between microbiologists and modelers/bioinformaticians, stemming from a lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored for microbiologists, researchers new to microbiome modeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metaomics methods, common modeling formalisms, and how models facilitate microbiome control. It concludes with guidelines and repositories for modeling. Each section provides entry-level information, example applications, and important references, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling.

Development of a genome-scale metabolic model for the lager hybrid yeast S. pastorianus to understand the evolution of metabolic pathways in industrial settings

In silico tools such as genome-scale metabolic models have shown to be powerful for metabolic engineering of microorganisms. Saccharomyces pastorianus is a complex aneuploid hybrid between the mesophilic Saccharomyces cerevisiae and the cold-tolerant Saccharomyces eubayanus. This species is of biotechnological importance because it is the primary yeast used in lager beer fermentation and is also a key model for studying the evolution of hybrid genomes, including expression pattern of ortholog genes, composition of protein complexes, and phenotypic plasticity. Here, we created the iSP_1513 GSMM for S. pastorianus CBS1513 to allow top-down computational approaches to predict the evolution of metabolic pathways and to aid strain optimization in production processes. The iSP_1513 comprises 4,062 reactions, 1,808 alleles, and 2,747 metabolites, and takes into account the functional redundancy in the gene-protein-reaction rule caused by the presence of orthologous genes. Moreover, a universal algorithm to constrain GSMM reactions using transcriptome data was developed as a python library and enabled the integration of temperature as parameter. Essentiality data sets, growth data on various carbohydrates and volatile metabolites secretion were used to validate the model and showed the potential of media engineering to improve specific flavor compounds. The iSP_1513 also highlighted the different contributions of the parental sub-genomes to the oxidative and non-oxidative parts of the pentose phosphate pathway. Overall, the iSP_1513 GSMM represent an important step toward understanding the metabolic capabilities, evolutionary trajectories, and adaptation potential of S. pastorianus in different industrial settings.

IMPORTANCE

Genome-scale metabolic models (GSMM) have been successfully applied to predict cellular behavior and design cell factories in several model organisms, but no models to date are currently available for hybrid species due to their more complex genetics and general lack of molecular data. In this study, we generated a bespoke GSMM, iSP_1513, for this industrial aneuploid hybrid Saccharomyces pastorianus, which takes into account the aneuploidy and functional redundancy from orthologous parental alleles. This model will (i) help understand the metabolic capabilities and adaptive potential of S. pastorianus (domestication processes), (ii) aid top-down predictions for strain development (industrial biotechnology), and (iii) allow predictions of evolutionary trajectories of metabolic pathways in aneuploid hybrids (evolutionary genetics).

Rewiring Saccharomyces cerevisiae metabolism for optimised Taxol® precursors production

Saccharomyces cerevisiae has been conveniently used to produce Taxol® anticancer drug early precursors. However, the harmful impact of oxidative stress by the first cytochrome P450-reductase enzymes (CYP725A4-POR) of Taxol® pathway has hampered sufficient progress in yeast. Here, we evolved an oxidative stress-resistant yeast strain with three-fold higher titre of their substrate, taxadiene. The performance of the evolved and parent strains were then evaluated in galactose-limited chemostats before and under the oxidative stress by an oxidising agent. The interaction of evolution and oxidative stress was comprehensively evaluated through transcriptomics and metabolite profiles integration in yeast enzyme-constrained genome scale model. Overall, the evolved strain showed improved respiration, reduced overflow metabolites production and oxidative stress re-induction tolerance. The cross-protection mechanism also potentially contributed to better heme, flavin and NADPH availability, essential for CYP725A4 and POR optimal activity in yeast. The results imply that the evolved strain is a robust cell factory for future efforts towards Taxol© production.

Emergence of Orchestrated and Dynamic Metabolism of Saccharomyces cerevisiae

Microbial metabolism is a fundamental cellular process that involves many biochemical events and is distinguished by its emergent properties. While the molecular details of individual reactions have been increasingly elucidated, it is not well understood how these reactions are quantitatively orchestrated to produce collective cellular behaviors. Here we developed a coarse-grained, systems, and dynamic mathematical framework, which integrates metabolic reactions with signal transduction and gene regulation to dissect the emergent metabolic traits of Saccharomyces cerevisiae. Our framework mechanistically captures a set of characteristic cellular behaviors, including the Crabtree effect, diauxic shift, diauxic lag time, and differential growth under nutrient-altered environments. It also allows modular expansion for zooming in on specific pathways for detailed metabolic profiles. This study provides a systems mathematical framework for yeast metabolic behaviors, providing insights into yeast physiology and metabolic engineering.

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.

Protein constraints in genome-scale metabolic models: Data integration, parameter estimation, and prediction of metabolic phenotypes

Genome-scale metabolic models provide a valuable resource to study metabolism and cell physiology. These models are employed with approaches from the constraint-based modeling framework to predict metabolic and physiological phenotypes. The prediction performance of genome-scale metabolic models can be improved by including protein constraints. The resulting protein-constrained models consider data on turnover numbers (kcat ) and facilitate the integration of protein abundances. In this systematic review, we present and discuss the current state-of-the-art regarding the estimation of kinetic parameters used in protein-constrained models. We also highlight how data-driven and constraint-based approaches can aid the estimation of turnover numbers and their usage in improving predictions of cellular phenotypes. Finally, we identify standing challenges in protein-constrained metabolic models and provide a perspective regarding future approaches to improve the predictive performance.

Resource allocation in mammalian systems

Cells execute biological functions to support phenotypes such as growth, migration, and secretion. Complementarily, each function of a cell has resource costs that constrain phenotype. Resource allocation by a cell allows it to manage these costs and optimize their phenotypes. In fact, the management of resource constraints (e.g., nutrient availability, bioenergetic capacity, and macromolecular machinery production) shape activity and ultimately impact phenotype. In mammalian systems, quantification of resource allocation provides important insights into higher-order multicellular functions; it shapes intercellular interactions and relays environmental cues for tissues to coordinate individual cells to overcome resource constraints and achieve population-level behavior. Furthermore, these constraints, objectives, and phenotypes are context-dependent, with cells adapting their behavior according to their microenvironment, resulting in distinct steady-states. This review will highlight the biological insights gained from probing resource allocation in mammalian cells and tissues.

Reconstruction, simulation and analysis of enzyme-constrained metabolic models using GECKO Toolbox 3.0

Genome-scale metabolic models (GEMs) are computational representations that enable mathematical exploration of metabolic behaviors within cellular and environmental constraints. Despite their wide usage in biotechnology, biomedicine and fundamental studies, there are many phenotypes that GEMs are unable to correctly predict. GECKO is a method to improve the predictive power of a GEM by incorporating enzymatic constraints using kinetic and omics data. GECKO has enabled reconstruction of enzyme-constrained metabolic models (ecModels) for diverse organisms, which show better predictive performance than conventional GEMs. In this protocol, we describe how to use the latest version GECKO 3.0; the procedure has five stages: (1) expansion from a starting metabolic model to an ecModel structure, (2) integration of enzyme turnover numbers into the ecModel structure, (3) model tuning, (4) integration of proteomics data into the ecModel and (5) simulation and analysis of ecModels. GECKO 3.0 incorporates deep learning-predicted enzyme kinetics, paving the way for improved metabolic models for virtually any organism and cell line in the absence of experimental data. The time of running the whole protocol is organism dependent, e.g., ~5 h for yeast.

Brain energy metabolism is optimized to minimize the cost of enzyme synthesis and transport

The energy metabolism of the brain is poorly understood partly due to the complex morphology of neurons and fluctuations in ATP demand over time. To investigate this, we used metabolic models that estimate enzyme usage per pathway, enzyme utilization over time, and enzyme transportation to evaluate how these parameters and processes affect ATP costs for enzyme synthesis and transportation. Our models show that the total enzyme maintenance energy expenditure of the human body depends on how glycolysis and mitochondrial respiration are distributed both across and within cell types in the brain. We suggest that brain metabolism is optimized to minimize the ATP maintenance cost by distributing the different ATP generation pathways in an advantageous way across cell types and potentially also across synapses within the same cell. Our models support this hypothesis by predicting export of lactate from both neurons and astrocytes during peak ATP demand, reproducing results from experimental measurements reported in the literature. Furthermore, our models provide potential explanation for parts of the astrocyte-neuron lactate shuttle theory, which is recapitulated under some conditions in the brain, while contradicting other aspects of the theory. We conclude that enzyme usage per pathway, enzyme utilization over time, and enzyme transportation are important factors for defining the optimal distribution of ATP production pathways, opening a broad avenue to explore in brain metabolism.

Simultaneous application of enzyme and thermodynamic constraints to metabolic models using an updated Python implementation of GECKO

IMPORTANCE

The metabolism of biological cells is an intricate network of reactions that interconvert chemical compounds, gathering energy, and using that energy to grow. The static analysis of these metabolic networks can be turned into a computational model that can efficiently output the distribution of fluxes in the network. With the inclusion of enzymes in the network, we can also interpret the role and concentrations of the metabolic proteins. However, the models and the experimental data often clash, resulting in a network that cannot grow. Here, we tackle this situation with a suite of relaxation algorithms in a package called geckopy. Geckopy also integrates with other software to allow for adding thermodynamic and metabolomic constraints. In addition, to ensure that enzyme-constrained models follow the community standards, a format for the proteins is postulated. We hope that the package and algorithms presented here will be useful for the constraint-based modeling community.

Advances in genome-scale metabolic models of industrially important fungi

Many fungal species have been used industrially for production of biofuels and bioproducts. Developing strains with better performance in biomanufacturing contexts requires a systematic understanding of cellular metabolism. Genome-scale metabolic models (GEMs) offer a comprehensive view of interconnected pathways and a mathematical framework for downstream analysis. Recently, GEMs have been developed or updated for several industrially important fungi. Some of them incorporate enzyme constraints, enabling improved predictions of cell states and proteome allocation. Here, we provide an overview of these newly developed GEMs and computational methods that facilitate construction of enzyme-constrained GEMs and utilize flux predictions from GEMs. Furthermore, we highlight the pivotal roles of these GEMs in iterative design-build-test-learn cycles, ultimately advancing the field of fungal biomanufacturing.

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.

Accurate prediction of in vivo protein abundances by coupling constraint-based modelling and machine learning

Quantification of how different environmental cues affect protein allocation can provide important insights for understanding cell physiology. While absolute quantification of proteins can be obtained by resource-intensive mass-spectrometry-based technologies, prediction of protein abundances offers another way to obtain insights into protein allocation. Here we present CAMEL, a framework that couples constraint-based modelling with machine learning to predict protein abundance for any environmental condition. This is achieved by building machine learning models that leverage static features, derived from protein sequences, and condition-dependent features predicted from protein-constrained metabolic models. Our findings demonstrate that CAMEL results in excellent prediction of protein allocation in E. coli (average Pearson correlation of at least 0.9), and moderate performance in S. cerevisiae (average Pearson correlation of at least 0.5). Therefore, CAMEL outperformed contending approaches without using molecular read-outs from unseen conditions and provides a valuable tool for using protein allocation in biotechnological applications.

PARROT: Prediction of enzyme abundances using protein-constrained metabolic models

Protein allocation determines the activity of cellular pathways and affects growth across all organisms. Therefore, different experimental and machine learning approaches have been developed to quantify and predict protein abundance and how they are allocated to different cellular functions, respectively. Yet, despite advances in protein quantification, it remains challenging to predict condition-specific allocation of enzymes in metabolic networks. Here, using protein-constrained metabolic models, we propose a family of constrained-based approaches, termed PARROT, to predict how much of each enzyme is used based on the principle of minimizing the difference between a reference and an alternative growth condition. To this end, PARROT variants model the minimization of enzyme reallocation using four different (combinations of) distance functions. We demonstrate that the PARROT variant that minimizes the Manhattan distance between the enzyme allocation of a reference and an alternative condition outperforms existing approaches based on the parsimonious distribution of fluxes or enzymes for both Escherichia coli and Saccharomyces cerevisiae. Further, we show that the combined minimization of flux and enzyme allocation adjustment leads to inconsistent predictions. Together, our findings indicate that minimization of protein allocation rather than flux redistribution is a governing principle determining steady-state pathway activity for microorganism grown in alternative growth conditions.

Model-driven insights into the effects of temperature on metabolism

Temperature affects cellular processes at different spatiotemporal scales, and identifying the genetic and molecular mechanisms underlying temperature responses paves the way to develop approaches for mitigating the effects of future climate scenarios. A systems view of the effects of temperature on cellular physiology can be obtained by focusing on metabolism since: (i) its functions depend on transcription and translation and (ii) its outcomes support organisms' development, growth, and reproduction. Here we provide a systematic review of modelling efforts directed at investigating temperature effects on properties of single biochemical reactions, system-level traits, metabolic subsystems, and whole-cell metabolism across different prokaryotes and eukaryotes. We compare and contrast computational approaches and theories that facilitate modelling of temperature effects on key properties of enzymes and their consideration in constraint-based as well as kinetic models of metabolism. In addition, we provide a summary of insights from computational approaches, facilitating integration of omics data from temperature-modulated experiments with models of metabolic networks, and review the resulting biotechnological applications. Lastly, we provide a perspective on how different types of metabolic modelling can profit from developments in machine learning and models of different cellular layers to improve model-driven insights into the effects of temperature relevant for biotechnological applications.

Enzyme-constrained metabolic model and in silico metabolic engineering of Clostridium ljungdahlii for the development of sustainable production processes

Constraint-based genome-scale models (GEMs) of microorganisms provide a powerful tool for predicting and analyzing microbial phenotypes as well as for understanding how these are affected by genetic and environmental perturbations. Recently, MATLAB and Python-based tools have been developed to incorporate enzymatic constraints into GEMs. These constraints enhance phenotype predictions by accounting for the enzyme cost of catalyzed model´s reactions, thereby reducing the space of possible metabolic flux distributions. In this study, enzymatic constraints were added to an existing GEM of Clostridium ljungdahlii, a model acetogenic bacterium, by including its enzyme turnover numbers (kcats) and molecular masses, using the Python-based AutoPACMEN approach. When compared to the metabolic model iHN637, the enzyme cost-constrained model (ec_iHN637) obtained in our study showed an improved predictive ability of growth rate and product profile. The model ec_iHN637 was then employed to perform in silico metabolic engineering of C. ljungdahlii, by using the OptKnock computational framework to identify knockouts to enhance the production of desired fermentation products. The in silico metabolic engineering was geared towards increasing the production of fermentation products by C. ljungdahlii, with a focus on the utilization of synthesis gas and CO2. This resulted in different engineering strategies for overproduction of valuable metabolites under different feeding conditions, without redundant knockouts for different products. Importantly, the results of the in silico engineering results indicated that the mixotrophic growth of C. ljungdahlii is a promising approach to coupling improved cell growth and acetate and ethanol productivity with net CO2 fixation.

Personalised modelling of clinical heterogeneity between medium-chain acyl-CoA dehydrogenase patients

BACKGROUND

Monogenetic inborn errors of metabolism cause a wide phenotypic heterogeneity that may even differ between family members carrying the same genetic variant. Computational modelling of metabolic networks may identify putative sources of this inter-patient heterogeneity. Here, we mainly focus on medium-chain acyl-CoA dehydrogenase deficiency (MCADD), the most common inborn error of the mitochondrial fatty acid oxidation (mFAO). It is an enigma why some MCADD patients-if untreated-are at risk to develop severe metabolic decompensations, whereas others remain asymptomatic throughout life. We hypothesised that an ability to maintain an increased free mitochondrial CoA (CoASH) and pathway flux might distinguish asymptomatic from symptomatic patients.

RESULTS

We built and experimentally validated, for the first time, a kinetic model of the human liver mFAO. Metabolites were partitioned according to their water solubility between the bulk aqueous matrix and the inner membrane. Enzymes are also either membrane-bound or in the matrix. This metabolite partitioning is a novel model attribute and improved predictions. MCADD substantially reduced pathway flux and CoASH, the latter due to the sequestration of CoA as medium-chain acyl-CoA esters. Analysis of urine from MCADD patients obtained during a metabolic decompensation showed an accumulation of medium- and short-chain acylcarnitines, just like the acyl-CoA pool in the MCADD model. The model suggested some rescues that increased flux and CoASH, notably increasing short-chain acyl-CoA dehydrogenase (SCAD) levels. Proteome analysis of MCADD patient-derived fibroblasts indeed revealed elevated levels of SCAD in a patient with a clinically asymptomatic state. This is a rescue for MCADD that has not been explored before. Personalised models based on these proteomics data confirmed an increased pathway flux and CoASH in the model of an asymptomatic patient compared to those of symptomatic MCADD patients.

CONCLUSIONS

We present a detailed, validated kinetic model of mFAO in human liver, with solubility-dependent metabolite partitioning. Personalised modelling of individual patients provides a novel explanation for phenotypic heterogeneity among MCADD patients. Further development of personalised metabolic models is a promising direction to improve individualised risk assessment, management and monitoring for inborn errors of metabolism.

Control analysis in the identification of key enzymes driving metabolic adaptations: Towards drug target discovery

Metabolic Control Analysis (MCA) marked a turning point in understanding the design principles of metabolic network control by establishing control coefficients as a means to quantify the degree of control that an enzyme exerts on flux or metabolite concentrations. MCA has demonstrated that control of metabolic pathways is distributed among many enzymes rather than depending on a single rate-limiting step. MCA also proved that this distribution depends not only on the stoichiometric structure of the network but also on other kinetic determinants, such as the degree of saturation of the enzyme active site, the distance to thermodynamic equilibrium, and metabolite feedback regulatory loops. Consequently, predicting the alterations that occur during metabolic adaptation in response to strong changes involving a redistribution in such control distribution can be challenging. Here, using the framework provided by MCA, we illustrate how control distribution in a metabolic pathway/network depends on enzyme kinetic determinants and to what extent the redistribution of control affects our predictions on candidate enzymes suitable as targets for small molecule inhibition in the drug discovery process. Our results uncover that kinetic determinants can lead to unexpected control distribution and outcomes that cannot be predicted solely from stoichiometric determinants. We also unveil that the inference of key enzyme-drivers of an observed metabolic adaptation can be dramatically improved using mean control coefficients and ruling out those enzyme activities that are associated with low control coefficients. As the use of constraint-based stoichiometric genome-scale metabolic models (GSMMs) becomes increasingly prevalent for identifying genes/enzymes that could be potential drug targets, we anticipate that incorporating kinetic determinants and ruling out enzymes with low control coefficients into GSMM workflows will facilitate more accurate predictions and reveal novel therapeutic targets.

Proteomics and constraint-based modelling reveal enzyme kinetic properties of Chlamydomonas reinhardtii on a genome scale

Metabolic engineering of microalgae offers a promising solution for sustainable biofuel production, and rational design of engineering strategies can be improved by employing metabolic models that integrate enzyme turnover numbers. However, the coverage of turnover numbers for Chlamydomonas reinhardtii, a model eukaryotic microalga accessible to metabolic engineering, is 17-fold smaller compared to the heterotrophic cell factory Saccharomyces cerevisiae. Here we generate quantitative protein abundance data of Chlamydomonas covering 2337 to 3708 proteins in various growth conditions to estimate in vivo maximum apparent turnover numbers. Using constrained-based modeling we provide proxies for in vivo turnover numbers of 568 reactions, representing a 10-fold increase over the in vitro data for Chlamydomonas. Integration of the in vivo estimates instead of in vitro values in a metabolic model of Chlamydomonas improved the accuracy of enzyme usage predictions. Our results help in extending the knowledge on uncharacterized enzymes and improve biotechnological applications of Chlamydomonas.

Practical sampling of constraint-based models: Optimized thinning boosts CHRR performance

Thinning is a sub-sampling technique to reduce the memory footprint of Markov chain Monte Carlo. Despite being commonly used, thinning is rarely considered efficient. For sampling constraint-based models, a highly relevant use-case in systems biology, we here demonstrate that thinning boosts computational and, thereby, sampling efficiencies of the widely used Coordinate Hit-and-Run with Rounding (CHRR) algorithm. By benchmarking CHRR with thinning with simplices and genome-scale metabolic networks of up to thousands of dimensions, we find a substantial increase in computational efficiency compared to unthinned CHRR, in our examples by orders of magnitude, as measured by the effective sample size per time (ESS/t), with performance gains growing with polytope (effective network) dimension. Using a set of benchmark models we derive a ready-to-apply guideline for tuning thinning to efficient and effective use of compute resources without requiring additional coding effort. Our guideline is validated using three (out-of-sample) large-scale networks and we show that it allows sampling convex polytopes uniformly to convergence in a fraction of time, thereby unlocking the rigorous investigation of hitherto intractable models. The derivation of our guideline is explained in detail, allowing future researchers to update it as needed as new model classes and more training data becomes available. CHRR with deliberate utilization of thinning thereby paves the way to keep pace with progressing model sizes derived with the constraint-based reconstruction and analysis (COBRA) tool set. Sampling and evaluation pipelines are available at https://jugit.fz-juelich.de/IBG-1/ModSim/fluxomics/chrrt.

Functional decomposition of metabolism allows a system-level quantification of fluxes and protein allocation towards specific metabolic functions

Quantifying the contribution of individual molecular components to complex cellular processes is a grand challenge in systems biology. Here we establish a general theoretical framework (Functional Decomposition of Metabolism, FDM) to quantify the contribution of every metabolic reaction to metabolic functions, e.g. the synthesis of biomass building blocks. FDM allowed for a detailed quantification of the energy and biosynthesis budget for growing Escherichia coli cells. Surprisingly, the ATP generated during the biosynthesis of building blocks from glucose almost balances the demand from protein synthesis, the largest energy expenditure known for growing cells. This leaves the bulk of the energy generated by fermentation and respiration unaccounted for, thus challenging the common notion that energy is a key growth-limiting resource. Moreover, FDM together with proteomics enables the quantification of enzymes contributing towards each metabolic function, allowing for a first-principle formulation of a coarse-grained model of global protein allocation based on the structure of the metabolic network.