1
|
Wu S, Qu Z, Chen D, Wu H, Caiyin Q, Qiao J. Deciphering and designing microbial communities by genome-scale metabolic modelling. Comput Struct Biotechnol J 2024; 23:1990-2000. [PMID: 38765607 PMCID: PMC11098673 DOI: 10.1016/j.csbj.2024.04.055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2024] [Revised: 04/21/2024] [Accepted: 04/21/2024] [Indexed: 05/22/2024] Open
Abstract
Microbial communities are shaped by the complex interactions among organisms and the environment. Genome-scale metabolic models (GEMs) can provide deeper insights into the complexity and ecological properties of various microbial communities, revealing their intricate interactions. Many researchers have modified GEMs for the microbial communities based on specific needs. Thus, GEMs need to be comprehensively summarized to better understand the trends in their development. In this review, we summarized the key developments in deciphering and designing microbial communities using different GEMs. A timeline of selected highlights in GEMs indicated that this area is evolving from the single-strain level to the microbial community level. Then, we outlined a framework for constructing GEMs of microbial communities. We also summarized the models and resources of static and dynamic community-level GEMs. We focused on the role of external environmental and intracellular resources in shaping the assembly of microbial communities. Finally, we discussed the key challenges and future directions of GEMs, focusing on the integration of GEMs with quorum sensing mechanisms, microbial ecology interactions, machine learning algorithms, and automatic modeling, all of which contribute to consortia-based applications in different fields.
Collapse
Affiliation(s)
- Shengbo Wu
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
- Zhejiang Shaoxing Research Institute of Tianjin University, Shaoxing 312300, China
| | - Zheping Qu
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
| | - Danlei Chen
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
- Zhejiang Shaoxing Research Institute of Tianjin University, Shaoxing 312300, China
| | - Hao Wu
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
- Zhejiang Shaoxing Research Institute of Tianjin University, Shaoxing 312300, China
| | - Qinggele Caiyin
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
- Zhejiang Shaoxing Research Institute of Tianjin University, Shaoxing 312300, China
- Key Laboratory of Systems Bioengineering, Ministry of Education (Tianjin University), Tianjin 300072, China
| | - Jianjun Qiao
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China
- Zhejiang Shaoxing Research Institute of Tianjin University, Shaoxing 312300, China
- Key Laboratory of Systems Bioengineering, Ministry of Education (Tianjin University), Tianjin 300072, China
- Frontiers Science Center for Synthetic Biology (Ministry of Education), Tianjin University, Tianjin 300072, China
| |
Collapse
|
2
|
Zare F, Fleming RMT. Integration of proteomic data with genome-scale metabolic models: A methodological overview. Protein Sci 2024; 33:e5150. [PMID: 39275997 PMCID: PMC11400636 DOI: 10.1002/pro.5150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Revised: 06/29/2024] [Accepted: 08/06/2024] [Indexed: 09/16/2024]
Abstract
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.
Collapse
Affiliation(s)
- Farid Zare
- School of Medicine, University of Galway, Galway, Ireland
| | | |
Collapse
|
3
|
Kong Y, Chen H, Huang X, Chang L, Yang B, Chen W. Precise metabolic modeling in post-omics era: accomplishments and perspectives. Crit Rev Biotechnol 2024:1-19. [PMID: 39198033 DOI: 10.1080/07388551.2024.2390089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Revised: 07/18/2024] [Accepted: 07/23/2024] [Indexed: 09/01/2024]
Abstract
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.
Collapse
Affiliation(s)
- Yawen Kong
- State Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi, P. R. China
- School of Food Science and Technology, Jiangnan University, Wuxi, P. R. China
| | - Haiqin Chen
- State Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi, P. R. China
- School of Food Science and Technology, Jiangnan University, Wuxi, P. R. China
| | - Xinlei Huang
- The Key Laboratory of Industrial Biotechnology, School of Biotechnology, Jiangnan University, Wuxi, P. R. China
| | - Lulu Chang
- State Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi, P. R. China
- School of Food Science and Technology, Jiangnan University, Wuxi, P. R. China
| | - Bo Yang
- State Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi, P. R. China
- School of Food Science and Technology, Jiangnan University, Wuxi, P. R. China
| | - Wei Chen
- State Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi, P. R. China
- School of Food Science and Technology, Jiangnan University, Wuxi, P. R. China
- National Engineering Research Center for Functional Food, Jiangnan University, Wuxi, P. R. China
| |
Collapse
|
4
|
Oftadeh O, Hatzimanikatis V. Genome-scale models of metabolism and expression predict the metabolic burden of recombinant protein expression. Metab Eng 2024; 84:109-116. [PMID: 38880390 DOI: 10.1016/j.ymben.2024.06.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2023] [Revised: 05/10/2024] [Accepted: 06/10/2024] [Indexed: 06/18/2024]
Abstract
The production of recombinant proteins in a host using synthetic constructs such as plasmids comes at the cost of detrimental effects such as reduced growth, energetic inefficiencies, and other stress responses, collectively known as metabolic burden. Increasing the number of copies of the foreign gene increases the metabolic load but increases the expression of the foreign protein. Thus, there is a trade-off between biomass and product yield in response to changes in heterologous gene copy number. This work proposes a computational method, rETFL (recombinant Expression and Thermodynamic Flux), for analyzing and predicting the responses of recombinant organisms to the introduction of synthetic constructs. rETFL is an extension to the ETFL formulations designed to reconstruct models of metabolism and expression (ME-models). We have illustrated the capabilities of the method in four studies to (i) capture the growth reduction in plasmid-containing E. coli and recombinant protein production; (ii) explore the trade-off between biomass and product yield as plasmid copy number is varied; (iii) predict the emergence of overflow metabolism in recombinant E. coli in agreement with experimental data; and (iv) investigate the individual pathways and enzymes affected by the presence of the plasmid. We anticipate that rETFL will serve as a comprehensive platform for integrating available omics data for recombinant organisms and making context-specific predictions that can help optimize recombinant expression systems for biopharmaceutical production and gene therapy.
Collapse
Affiliation(s)
- Omid Oftadeh
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), CH, 1015, Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), CH, 1015, Lausanne, Switzerland.
| |
Collapse
|
5
|
Schroeder WL, Suthers PF, Willis TC, Mooney EJ, Maranas CD. Current State, Challenges, and Opportunities in Genome-Scale Resource Allocation Models: A Mathematical Perspective. Metabolites 2024; 14:365. [PMID: 39057688 PMCID: PMC11278519 DOI: 10.3390/metabo14070365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2024] [Revised: 06/24/2024] [Accepted: 06/26/2024] [Indexed: 07/28/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Wheaton L. Schroeder
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA 16802, USA
- The Center for Bioenergy Innovation, Oak Ridge, TN 37830, USA
| | - Patrick F. Suthers
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA 16802, USA
- The Center for Bioenergy Innovation, Oak Ridge, TN 37830, USA
- DOE Center for Advanced Bioenergy and Bioproducts Innovation, The Pennsylvania State University, University Park, PA 16802, USA
| | - Thomas C. Willis
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA 16802, USA
- The Center for Bioenergy Innovation, Oak Ridge, TN 37830, USA
| | - Eric J. Mooney
- DOE Center for Advanced Bioenergy and Bioproducts Innovation, The Pennsylvania State University, University Park, PA 16802, USA
- Department of Biochemistry, Microbiology and Molecular Biology, The Pennsylvania State University, University Park, PA 16802, USA
| | - Costas D. Maranas
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA 16802, USA
- The Center for Bioenergy Innovation, Oak Ridge, TN 37830, USA
- DOE Center for Advanced Bioenergy and Bioproducts Innovation, The Pennsylvania State University, University Park, PA 16802, USA
| |
Collapse
|
6
|
Munro V, Kelly V, Messner CB, Kustatscher G. Cellular control of protein levels: A systems biology perspective. Proteomics 2024; 24:e2200220. [PMID: 38012370 DOI: 10.1002/pmic.202200220] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2023] [Revised: 11/13/2023] [Accepted: 11/15/2023] [Indexed: 11/29/2023]
Abstract
How cells regulate protein levels is a central question of biology. Over the past decades, molecular biology research has provided profound insights into the mechanisms and the molecular machinery governing each step of the gene expression process, from transcription to protein degradation. Recent advances in transcriptomics and proteomics have complemented our understanding of these fundamental cellular processes with a quantitative, systems-level perspective. Multi-omic studies revealed significant quantitative, kinetic and functional differences between the genome, transcriptome and proteome. While protein levels often correlate with mRNA levels, quantitative investigations have demonstrated a substantial impact of translation and protein degradation on protein expression control. In addition, protein-level regulation appears to play a crucial role in buffering protein abundances against undesirable mRNA expression variation. These findings have practical implications for many fields, including gene function prediction and precision medicine.
Collapse
Affiliation(s)
- Victoria Munro
- Wellcome Centre for Cell Biology, University of Edinburgh, Edinburgh, UK
| | - Van Kelly
- Wellcome Centre for Cell Biology, University of Edinburgh, Edinburgh, UK
| | - Christoph B Messner
- Precision Proteomics Center, Swiss Institute of Allergy and Asthma Research (SIAF), University of Zurich, Davos, Switzerland
| | - Georg Kustatscher
- Wellcome Centre for Cell Biology, University of Edinburgh, Edinburgh, UK
| |
Collapse
|
7
|
Goelzer A, Rajjou L, Chardon F, Loudet O, Fromion V. Resource allocation modeling for autonomous prediction of plant cell phenotypes. Metab Eng 2024; 83:86-101. [PMID: 38561149 DOI: 10.1016/j.ymben.2024.03.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2023] [Revised: 02/19/2024] [Accepted: 03/29/2024] [Indexed: 04/04/2024]
Abstract
Predicting the plant cell response in complex environmental conditions is a challenge in plant biology. Here we developed a resource allocation model of cellular and molecular scale for the leaf photosynthetic cell of Arabidopsis thaliana, based on the Resource Balance Analysis (RBA) constraint-based modeling framework. The RBA model contains the metabolic network and the major macromolecular processes involved in the plant cell growth and survival and localized in cellular compartments. We simulated the model for varying environmental conditions of temperature, irradiance, partial pressure of CO2 and O2, and compared RBA predictions to known resource distributions and quantitative phenotypic traits such as the relative growth rate, the C:N ratio, and finally to the empirical characteristics of CO2 fixation given by the well-established Farquhar model. In comparison to other standard constraint-based modeling methods like Flux Balance Analysis, the RBA model makes accurate quantitative predictions without the need for empirical constraints. Altogether, we show that RBA significantly improves the autonomous prediction of plant cell phenotypes in complex environmental conditions, and provides mechanistic links between the genotype and the phenotype of the plant cell.
Collapse
Affiliation(s)
- Anne Goelzer
- Université Paris-Saclay, INRAE, MaIAGE, 78350, Jouy-en-Josas, France.
| | - Loïc Rajjou
- Université Paris-Saclay, INRAE, AgroParisTech, Institut Jean-Pierre Bourgin (IJPB), 78000, Versailles, France
| | - Fabien Chardon
- Université Paris-Saclay, INRAE, AgroParisTech, Institut Jean-Pierre Bourgin (IJPB), 78000, Versailles, France
| | - Olivier Loudet
- Université Paris-Saclay, INRAE, AgroParisTech, Institut Jean-Pierre Bourgin (IJPB), 78000, Versailles, France
| | - Vincent Fromion
- Université Paris-Saclay, INRAE, MaIAGE, 78350, Jouy-en-Josas, France.
| |
Collapse
|
8
|
Ferreira MADM, Silveira WBD, Nikoloski Z. Protein constraints in genome-scale metabolic models: Data integration, parameter estimation, and prediction of metabolic phenotypes. Biotechnol Bioeng 2024; 121:915-930. [PMID: 38178617 DOI: 10.1002/bit.28650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 10/24/2023] [Accepted: 12/18/2023] [Indexed: 01/06/2024]
Abstract
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.
Collapse
Affiliation(s)
| | | | - Zoran Nikoloski
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, Germany
- Systems Biology and Mathematical Modeling, Max Planck Institute of Molecular Plant Physiology, Potsdam, Germany
| |
Collapse
|
9
|
Chen Y, Gustafsson J, Tafur Rangel A, Anton M, Domenzain I, Kittikunapong C, Li F, Yuan L, Nielsen J, Kerkhoven EJ. Reconstruction, simulation and analysis of enzyme-constrained metabolic models using GECKO Toolbox 3.0. Nat Protoc 2024; 19:629-667. [PMID: 38238583 DOI: 10.1038/s41596-023-00931-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2023] [Accepted: 10/13/2023] [Indexed: 03/10/2024]
Abstract
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.
Collapse
Affiliation(s)
- Yu Chen
- Department of Life Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Key Laboratory of Quantitative Synthetic Biology, Shenzhen Institute of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
| | - Johan Gustafsson
- Department of Life Sciences, Chalmers University of Technology, Gothenburg, Sweden
| | - Albert Tafur Rangel
- Department of Life Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Novo Nordisk Foundation Center for Biosustainability, Technology University of Denmark, Lyngby, Denmark
| | - Mihail Anton
- Department of Life Sciences, National Bioinformatics Infrastructure Sweden, Science for Life Laboratory, Chalmers University of Technology, Gothenburg, Sweden
| | - Iván Domenzain
- Department of Life Sciences, Chalmers University of Technology, Gothenburg, Sweden
| | | | - Feiran Li
- Department of Life Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Institute of Biopharmaceutical and Health Engineering, Tsinghua Shenzhen International Graduate School, Tsinghua University, Shenzhen, China
| | - Le Yuan
- Department of Life Sciences, Chalmers University of Technology, Gothenburg, Sweden
| | - Jens Nielsen
- Department of Life Sciences, Chalmers University of Technology, Gothenburg, Sweden
- BioInnovation Institute, Copenhagen, Denmark
| | - Eduard J Kerkhoven
- Department of Life Sciences, Chalmers University of Technology, Gothenburg, Sweden.
- Novo Nordisk Foundation Center for Biosustainability, Technology University of Denmark, Lyngby, Denmark.
- SciLifeLab, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden.
| |
Collapse
|
10
|
Hu M, Suthers PF, Maranas CD. KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states. Metab Eng 2024; 82:123-133. [PMID: 38336004 DOI: 10.1016/j.ymben.2024.02.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2023] [Revised: 01/24/2024] [Accepted: 02/06/2024] [Indexed: 02/12/2024]
Abstract
Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. This versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.
Collapse
Affiliation(s)
- Mengqi Hu
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, 16802, USA; DOE Center for Advanced Bioenergy and Bioproducts Innovation, USA
| | - Patrick F Suthers
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, 16802, USA; DOE Center for Advanced Bioenergy and Bioproducts Innovation, USA
| | - Costas D Maranas
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, 16802, USA; DOE Center for Advanced Bioenergy and Bioproducts Innovation, USA.
| |
Collapse
|
11
|
Széliová D, Müller S, Zanghellini J. Costs of ribosomal RNA stabilization affect ribosome composition at maximum growth rate. Commun Biol 2024; 7:196. [PMID: 38368456 PMCID: PMC10874399 DOI: 10.1038/s42003-024-05815-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2023] [Accepted: 01/12/2024] [Indexed: 02/19/2024] Open
Abstract
Ribosomes are key to cellular self-fabrication and limit growth rate. While most enzymes are proteins, ribosomes consist of 1/3 protein and 2/3 ribonucleic acid (RNA) (in E. coli).Here, we develop a mechanistic model of a self-fabricating cell, validated across diverse growth conditions. Through resource balance analysis (RBA), we explore the variation in maximum growth rate with ribosome composition, assuming constant kinetic parameters.Our model highlights the importance of RNA instability. If we neglect it, RNA synthesis is always cheaper than protein synthesis, leading to an RNA-only ribosome at maximum growth rate. Upon accounting for RNA turnover, we find that a mixed ribosome composed of RNA and proteins maximizes growth rate. To account for RNA turnover, we explore two scenarios regarding the activity of RNases. In (a) degradation is proportional to RNA content. In (b) ribosomal proteins cooperatively mitigate RNA instability by protecting it from misfolding and subsequent degradation. In both cases, higher protein content elevates protein synthesis costs and simultaneously lowers RNA turnover expenses, resulting in mixed RNA-protein ribosomes. Only scenario (b) aligns qualitatively with experimental data across varied growth conditions.Our research provides fresh insights into ribosome biogenesis and evolution, paving the way for understanding protein-rich ribosomes in archaea and mitochondria.
Collapse
Affiliation(s)
- Diana Széliová
- Department of Analytical Chemistry, University of Vienna, Vienna, 1090, Austria
| | - Stefan Müller
- Faculty of Mathematics, University of Vienna, Vienna, 1090, Austria
| | - Jürgen Zanghellini
- Department of Analytical Chemistry, University of Vienna, Vienna, 1090, Austria.
| |
Collapse
|
12
|
Flamholz AI, Goyal A, Fischer WW, Newman DK, Phillips R. The proteome is a terminal electron acceptor. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.31.578293. [PMID: 38352589 PMCID: PMC10862836 DOI: 10.1101/2024.01.31.578293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/25/2024]
Abstract
Microbial metabolism is impressively flexible, enabling growth even when available nutrients differ greatly from biomass in redox state. E. coli, for example, rearranges its physiology to grow on reduced and oxidized carbon sources through several forms of fermentation and respiration. To understand the limits on and evolutionary consequences of metabolic flexibility, we developed a mathematical model coupling redox chemistry with principles of cellular resource allocation. Our integrated model clarifies key phenomena, including demonstrating that autotrophs grow slower than heterotrophs because of constraints imposed by intracellular production of reduced carbon. Our model further indicates that growth is improved by adapting the redox state of biomass to nutrients, revealing an unexpected mode of evolution where proteins accumulate mutations benefiting organismal redox balance.
Collapse
Affiliation(s)
- Avi I. Flamholz
- Division of Biology and Biological Engineering, California Institute of Technology; Pasadena, CA 91125
| | - Akshit Goyal
- Physics of Living Systems, Department of Physics, Massachusetts Institute of Technology; Cambridge, MA 02139
- International Centre for Theoretical Sciences, Tata Institute of Fundamental Research; Bengaluru 560089
| | - Woodward W. Fischer
- Division of Geological & Planetary Sciences, California Institute of Technology; Pasadena, CA 91125
| | - Dianne K. Newman
- Division of Biology and Biological Engineering, California Institute of Technology; Pasadena, CA 91125
- Division of Geological & Planetary Sciences, California Institute of Technology; Pasadena, CA 91125
| | - Rob Phillips
- Division of Biology and Biological Engineering, California Institute of Technology; Pasadena, CA 91125
- Department of Physics, California Institute of Technology; Pasadena, CA 91125, USA
| |
Collapse
|
13
|
Joseph C, Zafeiropoulos H, Bernaerts K, Faust K. Predicting microbial interactions with approaches based on flux balance analysis: an evaluation. BMC Bioinformatics 2024; 25:36. [PMID: 38262921 PMCID: PMC10804772 DOI: 10.1186/s12859-024-05651-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2023] [Accepted: 01/11/2024] [Indexed: 01/25/2024] Open
Abstract
BACKGROUND Given a genome-scale metabolic model (GEM) of a microorganism and criteria for optimization, flux balance analysis (FBA) predicts the optimal growth rate and its corresponding flux distribution for a specific medium. FBA has been extended to microbial consortia and thus can be used to predict interactions by comparing in-silico growth rates for co- and monocultures. Although FBA-based methods for microbial interaction prediction are becoming popular, a systematic evaluation of their accuracy has not yet been performed. RESULTS Here, we evaluate the accuracy of FBA-based predictions of human and mouse gut bacterial interactions using growth data from the literature. For this, we collected 26 GEMs from the semi-curated AGORA database as well as four previously published curated GEMs. We tested the accuracy of three tools (COMETS, Microbiome Modeling Toolbox and MICOM) by comparing growth rates predicted in mono- and co-culture to growth rates extracted from the literature and also investigated the impact of different tool settings and media. We found that except for curated GEMs, predicted growth rates and their ratios (i.e. interaction strengths) do not correlate with growth rates and interaction strengths obtained from in vitro data. CONCLUSIONS Prediction of growth rates with FBA using semi-curated GEMs is currently not sufficiently accurate to predict interaction strengths reliably.
Collapse
Affiliation(s)
- Clémence Joseph
- Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Laboratory of Molecular Bacteriology, KU Leuven, 3000, Leuven, Belgium
| | - Haris Zafeiropoulos
- Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Laboratory of Molecular Bacteriology, KU Leuven, 3000, Leuven, Belgium
| | - Kristel Bernaerts
- Department of Chemical Engineering, Chemical and Biochemical Reactor Engineering and Safety (CREaS), KU Leuven, 3001, Leuven, Belgium
| | - Karoline Faust
- Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Laboratory of Molecular Bacteriology, KU Leuven, 3000, Leuven, Belgium.
| |
Collapse
|
14
|
Tummler K, Klipp E. Data integration strategies for whole-cell modeling. FEMS Yeast Res 2024; 24:foae011. [PMID: 38544322 PMCID: PMC11042497 DOI: 10.1093/femsyr/foae011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2023] [Revised: 03/15/2024] [Accepted: 03/26/2024] [Indexed: 04/25/2024] Open
Abstract
Data makes the world go round-and high quality data is a prerequisite for precise models, especially for whole-cell models (WCM). Data for WCM must be reusable, contain information about the exact experimental background, and should-in its entirety-cover all relevant processes in the cell. Here, we review basic requirements to data for WCM and strategies how to combine them. As a species-specific resource, we introduce the Yeast Cell Model Data Base (YCMDB) to illustrate requirements and solutions. We discuss recent standards for data as well as for computational models including the modeling process as data to be reported. We outline strategies for constructions of WCM despite their inherent complexity.
Collapse
Affiliation(s)
- Katja Tummler
- Humboldt-Universität zu Berlin, Faculty of Life Sciences, Institute of Biology, Theoretical Biophysics,, Invalidenstr. 42, 10115 Berlin, Germany
| | - Edda Klipp
- Humboldt-Universität zu Berlin, Faculty of Life Sciences, Institute of Biology, Theoretical Biophysics,, Invalidenstr. 42, 10115 Berlin, Germany
| |
Collapse
|
15
|
Cerk K, Ugalde‐Salas P, Nedjad CG, Lecomte M, Muller C, Sherman DJ, Hildebrand F, Labarthe S, Frioux C. Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing. Microb Biotechnol 2024; 17:e14396. [PMID: 38243750 PMCID: PMC10832553 DOI: 10.1111/1751-7915.14396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2023] [Revised: 11/27/2023] [Accepted: 12/20/2023] [Indexed: 01/21/2024] Open
Abstract
Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.
Collapse
Affiliation(s)
- Klara Cerk
- Quadram Institute BioscienceNorwichUK
- Earlham InstituteNorwichUK
| | | | - Chabname Ghassemi Nedjad
- Inria, University of Bordeaux, INRAETalenceFrance
- University of Bordeaux, CNRS, Bordeaux INP, LaBRI, UMR 5800TalenceFrance
| | - Maxime Lecomte
- Inria, University of Bordeaux, INRAETalenceFrance
- INRAE STLO¸University of RennesRennesFrance
| | | | | | - Falk Hildebrand
- Quadram Institute BioscienceNorwichUK
- Earlham InstituteNorwichUK
| | - Simon Labarthe
- Inria, University of Bordeaux, INRAETalenceFrance
- INRAE, University of Bordeaux, BIOGECO, UMR 1202CestasFrance
| | | |
Collapse
|
16
|
Qiu S, Zhao S, Yang A. DLTKcat: deep learning-based prediction of temperature-dependent enzyme turnover rates. Brief Bioinform 2023; 25:bbad506. [PMID: 38189538 PMCID: PMC10772988 DOI: 10.1093/bib/bbad506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2023] [Revised: 11/29/2023] [Accepted: 12/08/2023] [Indexed: 01/09/2024] Open
Abstract
The enzyme turnover rate, ${k}_{cat}$, quantifies enzyme kinetics by indicating the maximum efficiency of enzyme catalysis. Despite its importance, ${k}_{cat}$ values remain scarce in databases for most organisms, primarily because of the cost of experimental measurements. To predict ${k}_{cat}$ and account for its strong temperature dependence, DLTKcat was developed in this study and demonstrated superior performance (log10-scale root mean squared error = 0.88, R-squared = 0.66) than previously published models. Through two case studies, DLTKcat showed its ability to predict the effects of protein sequence mutations and temperature changes on ${k}_{cat}$ values. Although its quantitative accuracy is not high enough yet to model the responses of cellular metabolism to temperature changes, DLTKcat has the potential to eventually become a computational tool to describe the temperature dependence of biological systems.
Collapse
Affiliation(s)
- Sizhe Qiu
- Department of Engineering Science, University of Oxford, OX1 3PJ, United Kingdom
| | - Simiao Zhao
- Radcliffe Department of Medicine, University of Oxford, OX3 9DU, United Kingdom
| | - Aidong Yang
- Department of Engineering Science, University of Oxford, OX1 3PJ, United Kingdom
| |
Collapse
|
17
|
Moura Ferreira MAD, Wendering P, Arend M, Batista da Silveira W, Nikoloski Z. Accurate prediction of in vivo protein abundances by coupling constraint-based modelling and machine learning. Metab Eng 2023; 80:184-192. [PMID: 37802292 DOI: 10.1016/j.ymben.2023.09.014] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Revised: 09/10/2023] [Accepted: 09/25/2023] [Indexed: 10/08/2023]
Abstract
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.
Collapse
Affiliation(s)
| | - Philipp Wendering
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, 14476, Germany; Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, Potsdam, 14476, Germany
| | - Marius Arend
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, 14476, Germany; Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, Potsdam, 14476, Germany
| | | | - Zoran Nikoloski
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, 14476, Germany; Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, Potsdam, 14476, Germany.
| |
Collapse
|
18
|
Hu XP, Schroeder S, Lercher MJ. Proteome efficiency of metabolic pathways in Escherichia coli increases along the nutrient flow. mSystems 2023; 8:e0076023. [PMID: 37795991 PMCID: PMC10654084 DOI: 10.1128/msystems.00760-23] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Accepted: 08/24/2023] [Indexed: 10/06/2023] Open
Abstract
IMPORTANCE Protein translation is the most expensive cellular process in fast-growing bacteria, and efficient proteome usage should thus be under strong natural selection. However, recent studies show that a considerable part of the proteome is unneeded for instantaneous cell growth in Escherichia coli. We still lack a systematic understanding of how this excess proteome is distributed across different pathways as a function of the growth conditions. We estimated the minimal required proteome across growth conditions in E. coli and compared the predictions with experimental data. We found that the proteome allocated to the most expensive internal pathways, including translation and the synthesis of amino acids and cofactors, is near the minimally required levels. In contrast, transporters and central carbon metabolism show much higher proteome levels than the predicted minimal abundance. Our analyses show that the proteome fraction unneeded for instantaneous cell growth decreases along the nutrient flow in E. coli.
Collapse
Affiliation(s)
- Xiao-Pan Hu
- Institute for Computer Science, Heinrich Heine University, Düsseldorf, Germany
- Department of Biology, Heinrich Heine University, Düsseldorf, Germany
| | - Stefan Schroeder
- Institute for Computer Science, Heinrich Heine University, Düsseldorf, Germany
- Department of Biology, Heinrich Heine University, Düsseldorf, Germany
| | - Martin J. Lercher
- Institute for Computer Science, Heinrich Heine University, Düsseldorf, Germany
- Department of Biology, Heinrich Heine University, Düsseldorf, Germany
| |
Collapse
|
19
|
Ferreira MADM, da Silveira WB, Nikoloski Z. PARROT: Prediction of enzyme abundances using protein-constrained metabolic models. PLoS Comput Biol 2023; 19:e1011549. [PMID: 37856550 PMCID: PMC10617714 DOI: 10.1371/journal.pcbi.1011549] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Revised: 10/31/2023] [Accepted: 09/29/2023] [Indexed: 10/21/2023] Open
Abstract
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.
Collapse
Affiliation(s)
| | | | - Zoran Nikoloski
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, Germany
- Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, Potsdam, Germany
| |
Collapse
|
20
|
Wendering P, Nikoloski Z. Model-driven insights into the effects of temperature on metabolism. Biotechnol Adv 2023; 67:108203. [PMID: 37348662 DOI: 10.1016/j.biotechadv.2023.108203] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2023] [Revised: 05/22/2023] [Accepted: 06/18/2023] [Indexed: 06/24/2023]
Abstract
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.
Collapse
Affiliation(s)
- Philipp Wendering
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, 14476 Potsdam, Germany; Systems Biology and Mathematical Modeling, Max Planck Institute of Molecular Plant Physiology, 14476 Potsdam, Germany
| | - Zoran Nikoloski
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, 14476 Potsdam, Germany; Systems Biology and Mathematical Modeling, Max Planck Institute of Molecular Plant Physiology, 14476 Potsdam, Germany.
| |
Collapse
|
21
|
Bruggeman FJ, Teusink B, Steuer R. Trade-offs between the instantaneous growth rate and long-term fitness: Consequences for microbial physiology and predictive computational models. Bioessays 2023; 45:e2300015. [PMID: 37559168 DOI: 10.1002/bies.202300015] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Revised: 07/14/2023] [Accepted: 07/19/2023] [Indexed: 08/11/2023]
Abstract
Microbial systems biology has made enormous advances in relating microbial physiology to the underlying biochemistry and molecular biology. By meticulously studying model microorganisms, in particular Escherichia coli and Saccharomyces cerevisiae, increasingly comprehensive computational models predict metabolic fluxes, protein expression, and growth. The modeling rationale is that cells are constrained by a limited pool of resources that they allocate optimally to maximize fitness. As a consequence, the expression of particular proteins is at the expense of others, causing trade-offs between cellular objectives such as instantaneous growth, stress tolerance, and capacity to adapt to new environments. While current computational models are remarkably predictive for E. coli and S. cerevisiae when grown in laboratory environments, this may not hold for other growth conditions and other microorganisms. In this contribution, we therefore discuss the relationship between the instantaneous growth rate, limited resources, and long-term fitness. We discuss uses and limitations of current computational models, in particular for rapidly changing and adverse environments, and propose to classify microbial growth strategies based on Grimes's CSR framework.
Collapse
Affiliation(s)
- Frank J Bruggeman
- Systems Biology Lab/AIMMS, VU University, Amsterdam, The Netherlands
| | - Bas Teusink
- Systems Biology Lab/AIMMS, VU University, Amsterdam, The Netherlands
| | - Ralf Steuer
- Institute for Theoretical Biology (ITB), Institute for Biology, Humboldt-University of Berlin, Berlin, Germany
| |
Collapse
|
22
|
Arend M, Zimmer D, Xu R, Sommer F, Mühlhaus T, Nikoloski Z. Proteomics and constraint-based modelling reveal enzyme kinetic properties of Chlamydomonas reinhardtii on a genome scale. Nat Commun 2023; 14:4781. [PMID: 37553325 PMCID: PMC10409818 DOI: 10.1038/s41467-023-40498-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Accepted: 08/01/2023] [Indexed: 08/10/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Marius Arend
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, 14476, Potsdam, Germany
- Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, 14476, Potsdam, Germany
- Bioinformatics and Mathematical Modeling Department, Center of Plant Systems Biology and Biotechnology, 4000, Plovdiv, Bulgaria
| | - David Zimmer
- Computational Systems Biology, TU Kaiserslautern, 67663, Kaiserslautern, Germany
| | - Rudan Xu
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, 14476, Potsdam, Germany
- Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, 14476, Potsdam, Germany
| | - Frederik Sommer
- Molecular Biotechnology & Systems Biology, TU Kaiserslautern, 67663, Kaiserslautern, Germany
| | - Timo Mühlhaus
- Computational Systems Biology, TU Kaiserslautern, 67663, Kaiserslautern, Germany
| | - Zoran Nikoloski
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, 14476, Potsdam, Germany.
- Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, 14476, Potsdam, Germany.
- Bioinformatics and Mathematical Modeling Department, Center of Plant Systems Biology and Biotechnology, 4000, Plovdiv, Bulgaria.
| |
Collapse
|
23
|
Faure L, Mollet B, Liebermeister W, Faulon JL. A neural-mechanistic hybrid approach improving the predictive power of genome-scale metabolic models. Nat Commun 2023; 14:4669. [PMID: 37537192 PMCID: PMC10400647 DOI: 10.1038/s41467-023-40380-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Accepted: 07/19/2023] [Indexed: 08/05/2023] Open
Abstract
Constraint-based metabolic models have been used for decades to predict the phenotype of microorganisms in different environments. However, quantitative predictions are limited unless labor-intensive measurements of media uptake fluxes are performed. We show how hybrid neural-mechanistic models can serve as an architecture for machine learning providing a way to improve phenotype predictions. We illustrate our hybrid models with growth rate predictions of Escherichia coli and Pseudomonas putida grown in different media and with phenotype predictions of gene knocked-out Escherichia coli mutants. Our neural-mechanistic models systematically outperform constraint-based models and require training set sizes orders of magnitude smaller than classical machine learning methods. Our hybrid approach opens a doorway to enhancing constraint-based modeling: instead of constraining mechanistic models with additional experimental measurements, our hybrid models grasp the power of machine learning while fulfilling mechanistic constrains, thus saving time and resources in typical systems biology or biological engineering projects.
Collapse
Affiliation(s)
- Léon Faure
- MICALIS Institute, INRAE, AgroParisTech, University of Paris-Saclay, 78350, Jouy-en-Josas, France
| | - Bastien Mollet
- Ecole Normale Supérieure of Lyon, 69342, Lyon, France
- UMR MIA, INRAE, AgroParisTech, University of Paris-Saclay, 91120, Palaiseau, France
| | | | - Jean-Loup Faulon
- MICALIS Institute, INRAE, AgroParisTech, University of Paris-Saclay, 78350, Jouy-en-Josas, France.
- Manchester Institute of Biotechnology, University of Manchester, Manchester, M1 7DN, UK.
| |
Collapse
|
24
|
Dourado H, Liebermeister W, Ebenhöh O, Lercher MJ. Mathematical properties of optimal fluxes in cellular reaction networks at balanced growth. PLoS Comput Biol 2023; 19:e1011156. [PMID: 37279246 DOI: 10.1371/journal.pcbi.1011156] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2022] [Accepted: 05/04/2023] [Indexed: 06/08/2023] Open
Abstract
The physiology of biological cells evolved under physical and chemical constraints, such as mass conservation across the network of biochemical reactions, nonlinear reaction kinetics, and limits on cell density. For unicellular organisms, the fitness that governs this evolution is mainly determined by the balanced cellular growth rate. We previously introduced growth balance analysis (GBA) as a general framework to model and analyze such nonlinear systems, revealing important analytical properties of optimal balanced growth states. It has been shown that at optimality, only a minimal subset of reactions can have nonzero flux. However, no general principles have been established to determine if a specific reaction is active at optimality. Here, we extend the GBA framework to study the optimality of each biochemical reaction, and we identify the mathematical conditions determining whether a reaction is active or not at optimal growth in a given environment. We reformulate the mathematical problem in terms of a minimal number of dimensionless variables and use the Karush-Kuhn-Tucker (KKT) conditions to identify fundamental principles of optimal resource allocation in GBA models of any size and complexity. Our approach helps to identify from first principles the economic values of biochemical reactions, expressed as marginal changes in cellular growth rate; these economic values can be related to the costs and benefits of proteome allocation into the reactions' catalysts. Our formulation also generalizes the concepts of Metabolic Control Analysis to models of growing cells. We show how the extended GBA framework unifies and extends previous approaches of cellular modeling and analysis, putting forward a program to analyze cellular growth through the stationarity conditions of a Lagrangian function. GBA thereby provides a general theoretical toolbox for the study of fundamental mathematical properties of balanced cellular growth.
Collapse
Affiliation(s)
- Hugo Dourado
- Institute for Computer Science and Department of Biology, Heinrich-Heine Universität, Düsseldorf, Germany
| | | | - Oliver Ebenhöh
- Quantitative and Theoretical Biology, Heinrich-Heine Universität, Düsseldorf, Germany
| | - Martin J Lercher
- Institute for Computer Science and Department of Biology, Heinrich-Heine Universität, Düsseldorf, Germany
| |
Collapse
|
25
|
Wagner L, Jules M, Borkowski O. What remains from living cells in bacterial lysate-based cell-free systems. Comput Struct Biotechnol J 2023; 21:3173-3182. [PMID: 37333859 PMCID: PMC10275740 DOI: 10.1016/j.csbj.2023.05.025] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 05/23/2023] [Accepted: 05/23/2023] [Indexed: 06/20/2023] Open
Abstract
Because they mimic cells while offering an accessible and controllable environment, lysate-based cell-free systems (CFS) have emerged as valuable biotechnology tools for synthetic biology. Historically used to uncover fundamental mechanisms of life, CFS are nowadays used for a multitude of purposes, including protein production and prototyping of synthetic circuits. Despite the conservation of fundamental functions in CFS like transcription and translation, RNAs and certain membrane-embedded or membrane-bound proteins of the host cell are lost when preparing the lysate. As a result, CFS largely lack some essential properties of living cells, such as the ability to adapt to changing conditions, to maintain homeostasis and spatial organization. Regardless of the application, shedding light on the black-box of the bacterial lysate is necessary to fully exploit the potential of CFS. Most measurements of the activity of synthetic circuits in CFS and in vivo show significant correlations because these only require processes that are preserved in CFS, like transcription and translation. However, prototyping circuits of higher complexity that require functions that are lost in CFS (cell adaptation, homeostasis, spatial organization) will not show such a good correlation with in vivo conditions. Both for prototyping circuits of higher complexity and for building artificial cells, the cell-free community has developed devices to reconstruct cellular functions. This mini-review compares bacterial CFS to living cells, focusing on functional and cellular process differences and the latest developments in restoring lost functions through complementation of the lysate or device engineering.
Collapse
|
26
|
Bodeit O, Ben Samir I, Karr JR, Goelzer A, Liebermeister W. RBAtools: a programming interface for Resource Balance Analysis models. BIOINFORMATICS ADVANCES 2023; 3:vbad056. [PMID: 37179703 PMCID: PMC10168585 DOI: 10.1093/bioadv/vbad056] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Revised: 03/21/2023] [Accepted: 04/21/2023] [Indexed: 05/15/2023]
Abstract
Motivation Efficient resource allocation can contribute to an organism's fitness and can improve evolutionary success. Resource Balance Analysis (RBA) is a computational framework that models an organism's growth-optimal proteome configurations in various environments. RBA software enables the construction of RBA models on genome scale and the calculation of medium-specific, growth-optimal cell states including metabolic fluxes and the abundance of macromolecular machines. However, existing software lacks a simple programming interface for non-expert users, easy to use and interoperable with other software. Results The python package RBAtools provides convenient access to RBA models. As a flexible programming interface, it enables the implementation of custom workflows and the modification of existing genome-scale RBA models. Its high-level functions comprise simulation, model fitting, parameter screens, sensitivity analysis, variability analysis and the construction of Pareto fronts. Models and data are represented as structured tables and can be exported to common data formats for fluxomics and proteomics visualization. Availability and implementation RBAtools documentation, installation instructions and tutorials are available at https://sysbioinra.github.io/rbatools/. General information about RBA and related software can be found at rba.inrae.fr.
Collapse
Affiliation(s)
- Oliver Bodeit
- MaIAGE, Université Paris-Saclay, INRAE, 78350 Jouy-en-Josas, France
- Institute of Quantitative and Theoretical Biology, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
- Institute of Biology, Theoretical Biophysics, Humboldt-Universität zu Berlin, 10115 Berlin, Germany
- Institute of Biochemistry, Charité—Universitätsmedizin Berlin, 10117 Berlin, Germany
| | - Inès Ben Samir
- MaIAGE, Université Paris-Saclay, INRAE, 78350 Jouy-en-Josas, France
| | - Jonathan R Karr
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
| | - Anne Goelzer
- MaIAGE, Université Paris-Saclay, INRAE, 78350 Jouy-en-Josas, France
| | | |
Collapse
|
27
|
Wendering P, Arend M, Razaghi-Moghadam Z, Nikoloski Z. Data integration across conditions improves turnover number estimates and metabolic predictions. Nat Commun 2023; 14:1485. [PMID: 36932067 PMCID: PMC10023748 DOI: 10.1038/s41467-023-37151-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Accepted: 03/03/2023] [Indexed: 03/19/2023] Open
Abstract
Turnover numbers characterize a key property of enzymes, and their usage in constraint-based metabolic modeling is expected to increase the prediction accuracy of diverse cellular phenotypes. In vivo turnover numbers can be obtained by integrating reaction rate and enzyme abundance measurements from individual experiments. Yet, their contribution to improving predictions of condition-specific cellular phenotypes remains elusive. Here, we show that available in vitro and in vivo turnover numbers lead to poor prediction of condition-specific growth rates with protein-constrained models of Escherichia coli and Saccharomyces cerevisiae, particularly when protein abundances are considered. We demonstrate that correction of turnover numbers by simultaneous consideration of proteomics and physiological data leads to improved predictions of condition-specific growth rates. Moreover, the obtained estimates are more precise than corresponding in vitro turnover numbers. Therefore, our approach provides the means to correct turnover numbers and paves the way towards cataloguing kcatomes of other organisms.
Collapse
Affiliation(s)
- Philipp Wendering
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, Germany
- Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, Potsdam, Germany
| | - Marius Arend
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, Germany
- Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, Potsdam, Germany
| | - Zahra Razaghi-Moghadam
- Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, Potsdam, Germany
| | - Zoran Nikoloski
- Bioinformatics, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, Germany.
- Systems Biology and Mathematical Modelling, Max Planck Institute of Molecular Plant Physiology, Potsdam, Germany.
| |
Collapse
|
28
|
The choice of the objective function in flux balance analysis is crucial for predicting replicative lifespans in yeast. PLoS One 2022; 17:e0276112. [PMID: 36227951 PMCID: PMC9560524 DOI: 10.1371/journal.pone.0276112] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 09/28/2022] [Indexed: 11/10/2022] Open
Abstract
Flux balance analysis (FBA) is a powerful tool to study genome-scale models of the cellular metabolism, based on finding the optimal flux distributions over the network. While the objective function is crucial for the outcome, its choice, even though motivated by evolutionary arguments, has not been directly connected to related measures. Here, we used an available multi-scale mathematical model of yeast replicative ageing, integrating cellular metabolism, nutrient sensing and damage accumulation, to systematically test the effect of commonly used objective functions on features of replicative ageing in budding yeast, such as the number of cell divisions and the corresponding time between divisions. The simulations confirmed that assuming maximal growth is essential for reaching realistic lifespans. The usage of the parsimonious solution or the additional maximisation of a growth-independent energy cost can improve lifespan predictions, explained by either increased respiratory activity using resources otherwise allocated to cellular growth or by enhancing antioxidative activity, specifically in early life. Our work provides a new perspective on choosing the objective function in FBA by connecting it to replicative ageing.
Collapse
|
29
|
Anderson ME, Smith JL, Grossman AD. Multiple mechanisms for overcoming lethal over-initiation of DNA replication. Mol Microbiol 2022; 118:426-442. [PMID: 36053906 PMCID: PMC9825946 DOI: 10.1111/mmi.14976] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2022] [Revised: 08/14/2022] [Accepted: 08/25/2022] [Indexed: 01/12/2023]
Abstract
DNA replication is highly regulated and primarily controlled at the step of initiation. In bacteria, the replication initiator DnaA and the origin of replication oriC are the primary targets of regulation. Perturbations that increase or decrease replication initiation can cause a decrease in cell fitness. We found that multiple mechanisms, including an increase in replication elongation and a decrease in replication initiation, can compensate for lethal over-initiation. We found that in Bacillus subtilis, under conditions of rapid growth, loss of yabA, a negative regulator of replication initiation, caused a synthetic lethal phenotype when combined with the dnaA1 mutation that also causes replication over-initiation. We isolated several classes of suppressors that restored viability to dnaA1 ∆yabA double mutants. Some suppressors (relA, nrdR) stimulated replication elongation. Others (dnaC, cshA) caused a decrease in replication initiation. One class of suppressors decreased replication initiation in the dnaA1 ∆yabA mutant by causing a decrease in the amount of the replicative helicase, DnaC. We found that decreased levels of helicase in otherwise wild-type cells were sufficient to decrease replication initiation during rapid growth, indicating that the replicative helicase is limiting for replication initiation. Our results highlight the multiple mechanisms cells use to regulate DNA replication.
Collapse
Affiliation(s)
- Mary E. Anderson
- Department of BiologyMassachusetts Institute of TechnologyCambridgeMassachusettsUSA
| | - Janet L. Smith
- Department of BiologyMassachusetts Institute of TechnologyCambridgeMassachusettsUSA
| | - Alan D. Grossman
- Department of BiologyMassachusetts Institute of TechnologyCambridgeMassachusettsUSA
| |
Collapse
|
30
|
Wilken SE, Besançon M, Kratochvíl M, Foko Kuate CA, Trefois C, Gu W, Ebenhöh O. Interrogating the effect of enzyme kinetics on metabolism using differentiable constraint-based models. Metab Eng 2022; 74:72-82. [PMID: 36152931 DOI: 10.1016/j.ymben.2022.09.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Revised: 09/08/2022] [Accepted: 09/10/2022] [Indexed: 10/31/2022]
Abstract
Metabolic models are typically characterized by a large number of parameters. Traditionally, metabolic control analysis is applied to differential equation-based models to investigate the sensitivity of predictions to parameters. A corresponding theory for constraint-based models is lacking, due to their formulation as optimization problems. Here, we show that optimal solutions of optimization problems can be efficiently differentiated using constrained optimization duality and implicit differentiation. We use this to calculate the sensitivities of predicted reaction fluxes and enzyme concentrations to turnover numbers in an enzyme-constrained metabolic model of Escherichia coli. The sensitivities quantitatively identify rate limiting enzymes and are mathematically precise, unlike current finite difference based approaches used for sensitivity analysis. Further, efficient differentiation of constraint-based models unlocks the ability to use gradient information for parameter estimation. We demonstrate this by improving, genome-wide, the state-of-the-art turnover number estimates for E. coli. Finally, we show that this technique can be generalized to arbitrarily complex models. By differentiating the optimal solution of a model incorporating both thermodynamic and kinetic rate equations, the effect of metabolite concentrations on biomass growth can be elucidated. We benchmark these metabolite sensitivities against a large experimental gene knockdown study, and find good alignment between the predicted sensitivities and in vivo metabolome changes. In sum, we demonstrate several applications of differentiating optimal solutions of constraint-based metabolic models, and show how it connects to classic metabolic control analysis.
Collapse
Affiliation(s)
- St Elmo Wilken
- Institute of Quantitative and Theoretical Biology, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225, Düsseldorf, Germany; Cluster of Excellence on Plant Sciences, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225, Düsseldorf, Germany.
| | - Mathieu Besançon
- Department for AI in Society, Science, and Technology, Zuse Institute Berlin, Takustraße 7, 14195, Berlin, Germany
| | - Miroslav Kratochvíl
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Campus Belval, L-4367, Belvaux, Luxembourg
| | - Chilperic Armel Foko Kuate
- Institute of Quantitative and Theoretical Biology, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225, Düsseldorf, Germany
| | - Christophe Trefois
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Campus Belval, L-4367, Belvaux, Luxembourg
| | - Wei Gu
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Campus Belval, L-4367, Belvaux, Luxembourg
| | - Oliver Ebenhöh
- Institute of Quantitative and Theoretical Biology, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225, Düsseldorf, Germany; Cluster of Excellence on Plant Sciences, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225, Düsseldorf, Germany
| |
Collapse
|
31
|
Thuillier K, Baroukh C, Bockmayr A, Cottret L, Paulevé L, Siegel A. MERRIN: MEtabolic regulation rule INference from time series data. Bioinformatics 2022; 38:ii127-ii133. [PMID: 36124795 DOI: 10.1093/bioinformatics/btac479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
MOTIVATION Many techniques have been developed to infer Boolean regulations from a prior knowledge network (PKN) and experimental data. Existing methods are able to reverse-engineer Boolean regulations for transcriptional and signaling networks, but they fail to infer regulations that control metabolic networks. RESULTS We present a novel approach to infer Boolean rules for metabolic regulation from time-series data and a PKN. Our method is based on a combination of answer set programming and linear programming. By solving both combinatorial and linear arithmetic constraints, we generate candidate Boolean regulations that can reproduce the given data when coupled to the metabolic network. We evaluate our approach on a core regulated metabolic network and show how the quality of the predictions depends on the available kinetic, fluxomics or transcriptomics time-series data. AVAILABILITY AND IMPLEMENTATION Software available at https://github.com/bioasp/merrin. SUPPLEMENTARY INFORMATION Supplementary data are available at https://doi.org/10.5281/zenodo.6670164.
Collapse
Affiliation(s)
- Kerian Thuillier
- INRIA, CNRS, IRISA, University of Rennes, Rennes F-35000, France
| | - Caroline Baroukh
- LIPME, INRAE, CNRS, Université de Toulouse, Castanet-Tolosan F-31326, France
| | - Alexander Bockmayr
- Institute of Mathematics, Freie Universität Berlin, Berlin D-14195, Germany
| | - Ludovic Cottret
- LIPME, INRAE, CNRS, Université de Toulouse, Castanet-Tolosan F-31326, France
| | - Loïc Paulevé
- Univ. Bordeaux, Bordeaux INP, CNRS, LaBRI, UMR5800, F-33400 Talence, France
| | - Anne Siegel
- INRIA, CNRS, IRISA, University of Rennes, Rennes F-35000, France
| |
Collapse
|
32
|
Domenzain I, Sánchez B, Anton M, Kerkhoven EJ, Millán-Oropeza A, Henry C, Siewers V, Morrissey JP, Sonnenschein N, Nielsen J. Reconstruction of a catalogue of genome-scale metabolic models with enzymatic constraints using GECKO 2.0. Nat Commun 2022; 13:3766. [PMID: 35773252 PMCID: PMC9246944 DOI: 10.1038/s41467-022-31421-1] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2021] [Accepted: 06/16/2022] [Indexed: 01/08/2023] Open
Abstract
Genome-scale metabolic models (GEMs) have been widely used for quantitative exploration of the relation between genotype and phenotype. Streamlined integration of enzyme constraints and proteomics data into such models was first enabled by the GECKO toolbox, allowing the study of phenotypes constrained by protein limitations. Here, we upgrade the toolbox in order to enhance models with enzyme and proteomics constraints for any organism with a compatible GEM reconstruction. With this, enzyme-constrained models for the budding yeasts Saccharomyces cerevisiae, Yarrowia lipolytica and Kluyveromyces marxianus are generated to study their long-term adaptation to several stress factors by incorporation of proteomics data. Predictions reveal that upregulation and high saturation of enzymes in amino acid metabolism are common across organisms and conditions, suggesting the relevance of metabolic robustness in contrast to optimal protein utilization as a cellular objective for microbial growth under stress and nutrient-limited conditions. The functionality of GECKO is expanded with an automated framework for continuous and version-controlled update of enzyme-constrained GEMs, also producing such models for Escherichia coli and Homo sapiens. In this work, we facilitate the utilization of enzyme-constrained GEMs in basic science, metabolic engineering and synthetic biology purposes.
Collapse
Affiliation(s)
- Iván Domenzain
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden
- Novo Nordisk Foundation Center for Biosustainability, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden
| | - Benjamín Sánchez
- Department of Biotechnology and Biomedicine, Technical University of Denmark, 2800, Kongens Lyngby, Denmark
- Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, 2800, Kongens Lyngby, Denmark
| | - Mihail Anton
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden
- Department of Biology and Biological Engineering, National Bioinformatics Infrastructure Sweden, Science for Life Laboratory, Chalmers University of Technology, Kemivägen 10, SE-412 58, Gothenburg, Sweden
| | - Eduard J Kerkhoven
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden
- Novo Nordisk Foundation Center for Biosustainability, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden
| | - Aarón Millán-Oropeza
- Plateforme d'analyse protéomique Paris Sud-Ouest (PAPPSO), INRAE, MICALIS Institute, Université Paris-Saclay, 78350, Jouy-en-Josas, France
| | - Céline Henry
- Plateforme d'analyse protéomique Paris Sud-Ouest (PAPPSO), INRAE, MICALIS Institute, Université Paris-Saclay, 78350, Jouy-en-Josas, France
| | - Verena Siewers
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden
- Novo Nordisk Foundation Center for Biosustainability, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden
| | - John P Morrissey
- School of Microbiology, Environmental Research Institute and APC Microbiome Ireland, University College Cork, T12 K8AF, Cork, Ireland
| | - Nikolaus Sonnenschein
- Department of Biotechnology and Biomedicine, Technical University of Denmark, 2800, Kongens Lyngby, Denmark
| | - Jens Nielsen
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden.
- Novo Nordisk Foundation Center for Biosustainability, Chalmers University of Technology, SE-412 96, Gothenburg, Sweden.
- BioInnovation Institute, Ole Maaløes Vej 3, 2200, Copenhagen, Denmark.
| |
Collapse
|
33
|
Kerkhoven EJ. Advances in constraint-based models: methods for improved predictive power based on resource allocation constraints. Curr Opin Microbiol 2022; 68:102168. [PMID: 35691074 DOI: 10.1016/j.mib.2022.102168] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Revised: 05/03/2022] [Accepted: 05/05/2022] [Indexed: 11/18/2022]
Abstract
The concept of metabolic models with resource allocation constraints has been around for over a decade and has clear advantages even when implementation is relatively rudimentary. Nonetheless, the number of organisms for which such a model is reconstructed is low. Various approaches exist, from coarse-grained consideration of enzyme usage to fine-grained description of protein translation. These approaches are reviewed here, with a particular focus on user-friendly solutions that can introduce resource allocation constraints to metabolic models of any organism. The availability of kcat data is a major hurdle, where recent advances might help to fill in the numerous gaps that exist for this data, especially for nonmodel organisms.
Collapse
Affiliation(s)
- Eduard J Kerkhoven
- Department of Biology and Biological Engineering, Chalmers University of Technology, Kemivägen 10, SE412 96 Gothenburg, Sweden; Novo Nordisk Foundation Center for Biosustainability, Chalmers University of Technology, Kemivägen 10, SE-412 96 Gothenburg, Sweden.
| |
Collapse
|
34
|
Yunus IS, Lee TS. Applications of targeted proteomics in metabolic engineering: advances and opportunities. Curr Opin Biotechnol 2022; 75:102709. [DOI: 10.1016/j.copbio.2022.102709] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2021] [Revised: 02/15/2022] [Accepted: 02/23/2022] [Indexed: 12/22/2022]
|
35
|
Verhagen KJA, Eerden SA, Sikkema BJ, Wahl SA. Predicting Metabolic Adaptation Under Dynamic Substrate Conditions Using a Resource-Dependent Kinetic Model: A Case Study Using Saccharomyces cerevisiae. Front Mol Biosci 2022; 9:863470. [PMID: 35651815 PMCID: PMC9149170 DOI: 10.3389/fmolb.2022.863470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Accepted: 03/29/2022] [Indexed: 11/26/2022] Open
Abstract
Exposed to changes in their environment, microorganisms will adapt their phenotype, including metabolism, to ensure survival. To understand the adaptation principles, resource allocation-based approaches were successfully applied to predict an optimal proteome allocation under (quasi) steady-state conditions. Nevertheless, for a general, dynamic environment, enzyme kinetics will have to be taken into account which was not included in the linear resource allocation models. To this end, a resource-dependent kinetic model was developed and applied to the model organism Saccharomyces cerevisiae by combining published kinetic models and calibrating the model parameters to published proteomics and fluxomics datasets. Using this approach, we were able to predict specific proteomes at different dilution rates under chemostat conditions. Interestingly, the approach suggests that the occurrence of aerobic fermentation (Crabtree effect) in S. cerevisiae is not caused by space limitation in the total proteome but rather an effect of constraints on the mitochondria. When exposing the approach to repetitive, dynamic substrate conditions, the proteome space was allocated differently. Less space was predicted to be available for non-essential enzymes (reserve space). This could indicate that the perceived “overcapacity” present in experimentally measured proteomes may very likely serve a purpose in increasing the robustness of a cell to dynamic conditions, especially an increase of proteome space for the growth reaction as well as of the trehalose cycle that was shown to be essential in providing robustness upon stronger substrate perturbations. The model predictions of proteome adaptation to dynamic conditions were additionally evaluated against respective experimentally measured proteomes, which highlighted the model’s ability to accurately predict major proteome adaptation trends. This proof of principle for the approach can be extended to production organisms and applied for both understanding metabolic adaptation and improving industrial process design.
Collapse
Affiliation(s)
- K. J. A. Verhagen
- Department of Biotechnology, Delft University of Technology, Delft, Netherlands
| | - S. A. Eerden
- Department of Biotechnology, Delft University of Technology, Delft, Netherlands
| | - B. J. Sikkema
- Department of Biotechnology, Delft University of Technology, Delft, Netherlands
| | - S. A. Wahl
- Department of Biotechnology, Delft University of Technology, Delft, Netherlands
- Lehrstuhl für Bioverfahrenstechnik, FAU Erlangen-Nürnberg, Erlangen, Germany
- *Correspondence: S. A. Wahl,
| |
Collapse
|
36
|
Kleijn IT, Martínez-Segura A, Bertaux F, Saint M, Kramer H, Shahrezaei V, Marguerat S. Growth-rate-dependent and nutrient-specific gene expression resource allocation in fission yeast. Life Sci Alliance 2022; 5:e202101223. [PMID: 35228260 PMCID: PMC8886410 DOI: 10.26508/lsa.202101223] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2021] [Revised: 01/25/2022] [Accepted: 01/25/2022] [Indexed: 12/20/2022] Open
Abstract
Cellular resources are limited and their relative allocation to gene expression programmes determines physiological states and global properties such as the growth rate. Here, we determined the importance of the growth rate in explaining relative changes in protein and mRNA levels in the simple eukaryote Schizosaccharomyces pombe grown on non-limiting nitrogen sources. Although expression of half of fission yeast genes was significantly correlated with the growth rate, this came alongside wide-spread nutrient-specific regulation. Proteome and transcriptome often showed coordinated regulation but with notable exceptions, such as metabolic enzymes. Genes positively correlated with growth rate participated in every level of protein production apart from RNA polymerase II-dependent transcription. Negatively correlated genes belonged mainly to the environmental stress response programme. Critically, metabolic enzymes, which represent ∼55-70% of the proteome by mass, showed mostly condition-specific regulation. In summary, we provide a rich account of resource allocation to gene expression in a simple eukaryote, advancing our basic understanding of the interplay between growth-rate-dependent and nutrient-specific gene expression.
Collapse
Affiliation(s)
- Istvan T Kleijn
- Medical Research Council London Institute of Medical Sciences (MRC LMS), London, UK
- Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London, UK
| | - Amalia Martínez-Segura
- Medical Research Council London Institute of Medical Sciences (MRC LMS), London, UK
- Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK
| | - François Bertaux
- Medical Research Council London Institute of Medical Sciences (MRC LMS), London, UK
- Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London, UK
| | - Malika Saint
- Medical Research Council London Institute of Medical Sciences (MRC LMS), London, UK
- Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK
| | - Holger Kramer
- Medical Research Council London Institute of Medical Sciences (MRC LMS), London, UK
- Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK
| | - Vahid Shahrezaei
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London, UK
| | - Samuel Marguerat
- Medical Research Council London Institute of Medical Sciences (MRC LMS), London, UK
- Institute of Clinical Sciences (ICS), Faculty of Medicine, Imperial College London, London, UK
| |
Collapse
|
37
|
Muntoni AP, Braunstein A, Pagnani A, De Martino D, De Martino A. Relationship between fitness and heterogeneity in exponentially growing microbial populations. Biophys J 2022; 121:1919-1930. [DOI: 10.1016/j.bpj.2022.04.012] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2021] [Revised: 09/13/2021] [Accepted: 04/08/2022] [Indexed: 11/02/2022] Open
|
38
|
Quantitative metabolic fluxes regulated by trans-omic networks. Biochem J 2022; 479:787-804. [PMID: 35356967 PMCID: PMC9022981 DOI: 10.1042/bcj20210596] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2022] [Revised: 03/11/2022] [Accepted: 03/15/2022] [Indexed: 12/21/2022]
Abstract
Cells change their metabolism in response to internal and external conditions by regulating the trans-omic network, which is a global biochemical network with multiple omic layers. Metabolic flux is a direct measure of the activity of a metabolic reaction that provides valuable information for understanding complex trans-omic networks. Over the past decades, techniques to determine metabolic fluxes, including 13C-metabolic flux analysis (13C-MFA), flux balance analysis (FBA), and kinetic modeling, have been developed. Recent studies that acquire quantitative metabolic flux and multi-omic data have greatly advanced the quantitative understanding and prediction of metabolism-centric trans-omic networks. In this review, we present an overview of 13C-MFA, FBA, and kinetic modeling as the main techniques to determine quantitative metabolic fluxes, and discuss their advantages and disadvantages. We also introduce case studies with the aim of understanding complex metabolism-centric trans-omic networks based on the determination of metabolic fluxes.
Collapse
|
39
|
Baroukh C, Zemouri M, Genin S. Trophic preferences of the pathogen Ralstonia solanacearum and consequences on its growth in xylem sap. Microbiologyopen 2022; 11:e1240. [PMID: 35212480 PMCID: PMC8770891 DOI: 10.1002/mbo3.1240] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2021] [Accepted: 09/16/2021] [Indexed: 11/16/2022] Open
Abstract
Ralstonia solanacearum is one of the most destructive pathogens worldwide. In the last 30 years, the molecular mechanisms at the origin of R. solanacearum pathogenicity have been studied in depth. However, the nutrition status of the pathogen once inside the plant has been poorly investigated. Yet, the pathogen needs substrates to sustain a fast-enough growth, maintain its virulence and subvert the host immunity. This study aimed to explore in-depth the xylem environment where the pathogen is abundant, and its trophic preferences. First, we determined the composition of tomato xylem sap, where fast multiplication of the pathogen occurs. Then, kinetic growth on single and mixtures of carbon sources in relation to this environment was performed to fully quantify growth. Finally, we calculated the concentration of available metabolites in the xylem sap flux to assess how much it can support bacterial growth in planta. Overall, the study underlines the adaptation of R. solanacearum to the xylem environment and the fact that the pathogen assimilates several substrates at the same time in media composed of several carbon sources. It also provides metrics on key physiological parameters governing the growth of this major pathogen, which will be instrumental in the future to better understand its metabolic behavior during infection.
Collapse
Affiliation(s)
| | - Meriem Zemouri
- LIPMEINRACNRSUniversité de ToulouseCastanet‐TolosanFrance
| | - Stéphane Genin
- LIPMEINRACNRSUniversité de ToulouseCastanet‐TolosanFrance
| |
Collapse
|
40
|
Kim M, Sung J, Chia N. Resource-allocation constraint governs structure and function of microbial communities in metabolic modeling. Metab Eng 2022; 70:12-22. [PMID: 34990848 DOI: 10.1016/j.ymben.2021.12.011] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Revised: 12/01/2021] [Accepted: 12/29/2021] [Indexed: 10/19/2022]
Abstract
Predictive modeling tools for assessing microbial communities are important for realizing transformative capabilities of microbiomes in agriculture, ecology, and medicine. Constraint-based community-scale metabolic modeling is unique in its potential for making mechanistic predictions regarding both the structure and function of microbial communities. However, accessing this potential requires an understanding of key physicochemical constraints, which are typically considered on a per-species basis. What is needed is a means of incorporating global constraints relevant to microbial ecology into community models. Resource-allocation constraint, which describes how limited resources should be distributed to different cellular processes, sets limits on the efficiency of metabolic and ecological processes. In this study, we investigate the implications of resource-allocation constraints in community-scale metabolic modeling through a simple mechanism-agnostic implementation of resource-allocation constraints directly at the flux level. By systematically performing single-, two-, and multi-species growth simulations, we show that resource-allocation constraints are indispensable for predicting the structure and function of microbial communities. Our findings call for a scalable workflow for implementing a mechanistic version of resource-allocation constraints to ultimately harness the full potential of community-scale metabolic modeling tools.
Collapse
Affiliation(s)
- Minsuk Kim
- Microbiome Program, Center for Individualized Medicine, Mayo Clinic, Rochester, MN, 55905, USA; Division of Surgical Research, Department of Surgery, Mayo Clinic, Rochester, MN, 55905, USA
| | - Jaeyun Sung
- Microbiome Program, Center for Individualized Medicine, Mayo Clinic, Rochester, MN, 55905, USA; Division of Surgical Research, Department of Surgery, Mayo Clinic, Rochester, MN, 55905, USA; Division of Rheumatology, Department of Medicine, Mayo Clinic, Rochester, MN, 55905, USA
| | - Nicholas Chia
- Microbiome Program, Center for Individualized Medicine, Mayo Clinic, Rochester, MN, 55905, USA; Division of Surgical Research, Department of Surgery, Mayo Clinic, Rochester, MN, 55905, USA.
| |
Collapse
|
41
|
Jahn M, Crang N, Janasch M, Hober A, Forsström B, Kimler K, Mattausch A, Chen Q, Asplund-Samuelsson J, Hudson EP. Protein allocation and utilization in the versatile chemolithoautotroph Cupriavidus necator. eLife 2021; 10:69019. [PMID: 34723797 PMCID: PMC8591527 DOI: 10.7554/elife.69019] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2021] [Accepted: 10/30/2021] [Indexed: 12/12/2022] Open
Abstract
Bacteria must balance the different needs for substrate assimilation, growth
functions, and resilience in order to thrive in their environment. Of all
cellular macromolecules, the bacterial proteome is by far the most important
resource and its size is limited. Here, we investigated how the highly versatile
'knallgas' bacterium Cupriavidus necator reallocates protein
resources when grown on different limiting substrates and with different growth
rates. We determined protein quantity by mass spectrometry and estimated enzyme
utilization by resource balance analysis modeling. We found that C.
necator invests a large fraction of its proteome in functions that
are hardly utilized. Of the enzymes that are utilized, many are present in
excess abundance. One prominent example is the strong expression of CBB cycle
genes such as Rubisco during growth on fructose. Modeling and mutant competition
experiments suggest that CO2-reassimilation through Rubisco does not
provide a fitness benefit for heterotrophic growth, but is rather an investment
in readiness for autotrophy.
Collapse
Affiliation(s)
- Michael Jahn
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| | - Nick Crang
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| | - Markus Janasch
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| | - Andreas Hober
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| | - Björn Forsström
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| | - Kyle Kimler
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| | - Alexander Mattausch
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| | - Qi Chen
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| | - Johannes Asplund-Samuelsson
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| | - Elton Paul Hudson
- School of Engineering Sciences in Chemistry, Biotechnology and Health, Science for Life Laboratory, KTH - Royal Institute of Technology, Stockholm, Sweden
| |
Collapse
|
42
|
Zeng H, Rohani R, Huang WE, Yang A. Understanding and mathematical modelling of cellular resource allocation in microorganisms: a comparative synthesis. BMC Bioinformatics 2021; 22:467. [PMID: 34583645 PMCID: PMC8479906 DOI: 10.1186/s12859-021-04382-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Accepted: 09/20/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND The rising consensus that the cell can dynamically allocate its resources provides an interesting angle for discovering the governing principles of cell growth and metabolism. Extensive efforts have been made in the past decade to elucidate the relationship between resource allocation and phenotypic patterns of microorganisms. Despite these exciting developments, there is still a lack of explicit comparison between potentially competing propositions and a lack of synthesis of inter-related proposals and findings. RESULTS In this work, we have reviewed resource allocation-derived principles, hypotheses and mathematical models to recapitulate important achievements in this area. In particular, the emergence of resource allocation phenomena is deciphered by the putative tug of war between the cellular objectives, demands and the supply capability. Competing hypotheses for explaining the most-studied phenomenon arising from resource allocation, i.e. the overflow metabolism, have been re-examined towards uncovering the potential physiological root cause. The possible link between proteome fractions and the partition of the ribosomal machinery has been analysed through mathematical derivations. Finally, open questions are highlighted and an outlook on the practical applications is provided. It is the authors' intention that this review contributes to a clearer understanding of the role of resource allocation in resolving bacterial growth strategies, one of the central questions in microbiology. CONCLUSIONS We have shown the importance of resource allocation in understanding various aspects of cellular systems. Several important questions such as the physiological root cause of overflow metabolism and the correct interpretation of 'protein costs' are shown to remain open. As the understanding of the mechanisms and utility of resource application in cellular systems further develops, we anticipate that mathematical modelling tools incorporating resource allocation will facilitate the circuit-host design in synthetic biology.
Collapse
Affiliation(s)
- Hong Zeng
- Beijing Advanced Innovation Center for Food Nutrition and Human Health, Beijing Technology and Business University, Beijing, 100048, China
| | - Reza Rohani
- Department of Engineering Science, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK
| | - Wei E Huang
- Department of Engineering Science, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK
| | - Aidong Yang
- Department of Engineering Science, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK.
| |
Collapse
|
43
|
Combining Kinetic and Constraint-Based Modelling to Better Understand Metabolism Dynamics. Processes (Basel) 2021. [DOI: 10.3390/pr9101701] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
To understand the phenotypic capabilities of organisms, it is useful to characterise cellular metabolism through the analysis of its pathways. Dynamic mathematical modelling of metabolic networks is of high interest as it provides the time evolution of the metabolic components. However, it also has limitations, such as the necessary mechanistic details and kinetic parameters are not always available. On the other hand, large metabolic networks exhibit a complex topological structure which can be studied rather efficiently in their stationary regime by constraint-based methods. These methods produce useful predictions on pathway operations. In this review, we present both modelling techniques and we show how they bring complementary views of metabolism. In particular, we show on a simple example how both approaches can be used in conjunction to shed some light on the dynamics of metabolic networks.
Collapse
|
44
|
Alsiyabi A, Chowdhury NB, Long D, Saha R. Enhancing in silico strain design predictions through next generation metabolic modeling approaches. Biotechnol Adv 2021; 54:107806. [PMID: 34298108 DOI: 10.1016/j.biotechadv.2021.107806] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2021] [Revised: 06/22/2021] [Accepted: 07/15/2021] [Indexed: 02/06/2023]
Abstract
The reconstruction and analysis of metabolic models has garnered increasing attention due to the multitude of applications in which these have proven to be practical. The growing number of generated metabolic models has been accompanied by an exponentially expanding arsenal of tools used to analyze them. In this work, we discussed the biological relevance of a number of promising modeling frameworks, focusing on the questions and hypotheses each method is equipped to address. To this end, we critically analyzed the steady-state modeling approaches focusing on resource allocation and incorporation of thermodynamic considerations which produce promising results and aid in the generation and experimental validation of numerous predictions. For smaller networks involving more complex regulation, we addressed kinetic modeling techniques which show encouraging results in addressing questions outside the scope of steady-state modeling. Finally, we discussed the potential application of the discussed frameworks within the field of strain design. Adoption of such methodologies is believed to significantly enhance the accuracy of in silico predictions and hence decrease the number of design-build-test cycles required.
Collapse
Affiliation(s)
- Adil Alsiyabi
- Chemical and Biomolecular Engineering, University of Nebraska-Lincoln, United States of America
| | - Niaz Bahar Chowdhury
- Chemical and Biomolecular Engineering, University of Nebraska-Lincoln, United States of America
| | - Dianna Long
- Complex Biosystems, University of Nebraska-Lincoln, United States of America
| | - Rajib Saha
- Chemical and Biomolecular Engineering, University of Nebraska-Lincoln, United States of America; Complex Biosystems, University of Nebraska-Lincoln, United States of America.
| |
Collapse
|
45
|
Petrizzelli MS, de Vienne D, Nidelet T, Noûs C, Dillmann C. Data integration uncovers the metabolic bases of phenotypic variation in yeast. PLoS Comput Biol 2021; 17:e1009157. [PMID: 34264947 PMCID: PMC8315545 DOI: 10.1371/journal.pcbi.1009157] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2021] [Revised: 07/27/2021] [Accepted: 06/07/2021] [Indexed: 12/13/2022] Open
Abstract
The relationship between different levels of integration is a key feature for understanding the genotype-phenotype map. Here, we describe a novel method of integrated data analysis that incorporates protein abundance data into constraint-based modeling to elucidate the biological mechanisms underlying phenotypic variation. Specifically, we studied yeast genetic diversity at three levels of phenotypic complexity in a population of yeast obtained by pairwise crosses of eleven strains belonging to two species, Saccharomyces cerevisiae and S. uvarum. The data included protein abundances, integrated traits (life-history/fermentation) and computational estimates of metabolic fluxes. Results highlighted that the negative correlation between production traits such as population carrying capacity (K) and traits associated with growth and fermentation rates (Jmax) is explained by a differential usage of energy production pathways: a high K was associated with high TCA fluxes, while a high Jmax was associated with high glycolytic fluxes. Enrichment analysis of protein sets confirmed our results. This powerful approach allowed us to identify the molecular and metabolic bases of integrated trait variation, and therefore has a broad applicability domain. The integration of data at different levels of cellular organization is an important goal in computational biology for understanding the way the genotypic variation translates into phenotypic variation. Novel profiling technologies and accurate high-throughput phenotyping now allows genomic, transcriptomic, metabolic and proteomic characterization of a large number of individuals under various environmental conditions. However, the metabolic fluxes remain difficult to measure. In this work, we take advantage of recent advances in genome-scale functional annotation and constraint-based metabolic modeling to provide a mathematical framework that allows to estimate internal cellular fluxes from protein abundances and elucidate the biological mechanisms underlying phenotypic variation. Applied to yeast as a model system, this approach highlights that the negative correlation between production traits such as maximum population size and growth and fermentation traits is explained by a differential usage of energy production pathways. The ability to identify molecular and metabolic bases of the variation of integrated traits through population studies has a broad applicability domain.
Collapse
Affiliation(s)
- Marianyela Sabina Petrizzelli
- Université Paris-Saclay, INRAE, CNRS, AgroParisTech, GQE–Le Moulon, Gif-sur-Yvette, France
- Institut Curie, PSL Research University, Paris, France
- INSERM, U900, Paris, France
- CBIO-Centre for Computational Biology, MINES ParisTech, PSL Research University, Paris, France
- * E-mail:
| | - Dominique de Vienne
- Université Paris-Saclay, INRAE, CNRS, AgroParisTech, GQE–Le Moulon, Gif-sur-Yvette, France
| | - Thibault Nidelet
- SPO, INRAE, Montpellier SupAgro, Université de Montpellier, Montpellier, France
| | | | - Christine Dillmann
- Université Paris-Saclay, INRAE, CNRS, AgroParisTech, GQE–Le Moulon, Gif-sur-Yvette, France
| |
Collapse
|
46
|
Zhou J, Zhuang Y, Xia J. Integration of enzyme constraints in a genome-scale metabolic model of Aspergillus niger improves phenotype predictions. Microb Cell Fact 2021; 20:125. [PMID: 34193117 PMCID: PMC8247156 DOI: 10.1186/s12934-021-01614-2] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2021] [Accepted: 06/16/2021] [Indexed: 11/26/2022] Open
Abstract
Background Genome-scale metabolic model (GSMM) is a powerful tool for the study of cellular metabolic characteristics. With the development of multi-omics measurement techniques in recent years, new methods that integrating multi-omics data into the GSMM show promising effects on the predicted results. It does not only improve the accuracy of phenotype prediction but also enhances the reliability of the model for simulating complex biochemical phenomena, which can promote theoretical breakthroughs for specific gene target identification or better understanding the cell metabolism on the system level. Results Based on the basic GSMM model iHL1210 of Aspergillus niger, we integrated large-scale enzyme kinetics and proteomics data to establish a GSMM based on enzyme constraints, termed a GEM with Enzymatic Constraints using Kinetic and Omics data (GECKO). The results show that enzyme constraints effectively improve the model’s phenotype prediction ability, and extended the model’s potential to guide target gene identification through predicting metabolic phenotype changes of A. niger by simulating gene knockout. In addition, enzyme constraints significantly reduced the solution space of the model, i.e., flux variability over 40.10% metabolic reactions were significantly reduced. The new model showed also versatility in other aspects, like estimating large-scale \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$k_{{cat}}$$\end{document}kcat values, predicting the differential expression of enzymes under different growth conditions. Conclusions This study shows that incorporating enzymes’ abundance information into GSMM is very effective for improving model performance with A. niger. Enzyme-constrained model can be used as a powerful tool for predicting the metabolic phenotype of A. niger by incorporating proteome data. In the foreseeable future, with the fast development of measurement techniques, and more precise and rich proteomics quantitative data being obtained for A. niger, the enzyme-constrained GSMM model will show greater application space on the system level. Supplementary Information The online version contains supplementary material available at 10.1186/s12934-021-01614-2.
Collapse
Affiliation(s)
- Jingru Zhou
- State Key Laboratory of Bioreactor Engineering, East China University of Science and Technology, Shanghai, 200237, China
| | - Yingping Zhuang
- State Key Laboratory of Bioreactor Engineering, East China University of Science and Technology, Shanghai, 200237, China
| | - Jianye Xia
- State Key Laboratory of Bioreactor Engineering, East China University of Science and Technology, Shanghai, 200237, China. .,Tianjin Institute of Industrial Biotechnology, Chinese Academy of Science, Tianjin, 300308, China.
| |
Collapse
|
47
|
Lloyd CJ, Monk J, Yang L, Ebrahim A, Palsson BO. Computation of condition-dependent proteome allocation reveals variability in the macro and micro nutrient requirements for growth. PLoS Comput Biol 2021; 17:e1007817. [PMID: 34161321 PMCID: PMC8259983 DOI: 10.1371/journal.pcbi.1007817] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2020] [Revised: 07/06/2021] [Accepted: 05/31/2021] [Indexed: 11/21/2022] Open
Abstract
Sustaining a robust metabolic network requires a balanced and fully functioning proteome. In addition to amino acids, many enzymes require cofactors (coenzymes and engrafted prosthetic groups) to function properly. Extensively validated resource allocation models, such as genome-scale models of metabolism and gene expression (ME-models), have the ability to compute an optimal proteome composition underlying a metabolic phenotype, including the provision of all required cofactors. Here we apply the ME-model for Escherichia coli K-12 MG1655 to computationally examine how environmental conditions change the proteome and its accompanying cofactor usage. We found that: (1) The cofactor requirements computed by the ME-model mostly agree with the standard biomass objective function used in models of metabolism alone (M-models); (2) ME-model computations reveal non-intuitive variability in cofactor use under different growth conditions; (3) An analysis of ME-model predicted protein use in aerobic and anaerobic conditions suggests an enrichment in the use of peroxyl scavenging acids in the proteins used to sustain aerobic growth; (4) The ME-model could describe how limitation in key protein components affect the metabolic state of E. coli. Genome-scale models have thus reached a level of sophistication where they reveal intricate properties of functional proteomes and how they support different E. coli lifestyles. Escherichia coli is capable of growing in many environments, each of which requires a different collection of enzymes to metabolize the nutrients within that environment. Each individual enzyme requires its own set of amino acids and oftentimes cofactors, which are accessory molecules essential for the enzyme to function. Thus, the composition of the micronutrients (amino acids, cofactors, etc.) within a cell will differ depending on its metabolic needs. The presented work is the first effort to employ metabolic models to probe the connection between E. coli’s diverse growth environments and its biomass composition. We first show how differences in model-predicted enzyme use for aerobic or anaerobic growth results in distinct amino acid and cofactor usage. Alternatively, we show that the metabolic models can predict how modifying the cell’s biomass composition will affect growth. For example, by modeling the exposure of E. coli to trimethoprim or sulfamethoxazole—two antibiotics that target folate (vitamin B9) synthesis—we predicted how E. coli could adapt to grow under folate-limited conditions. This work demonstrates how models can be used to study antibiotic resistance of drugs that target amino acid or cofactor synthesis.
Collapse
Affiliation(s)
- Colton J. Lloyd
- Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America
| | - Jonathan Monk
- Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America
| | - Laurence Yang
- Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America
| | - Ali Ebrahim
- Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America
| | - Bernhard O. Palsson
- Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America
- Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Lyngby, Denmark
- Department of Pediatrics, University of California, San Diego, La Jolla, California, United States of America
- * E-mail:
| |
Collapse
|
48
|
Bruggeman FJ, Planqué R, Molenaar D, Teusink B. Searching for principles of microbial physiology. FEMS Microbiol Rev 2021; 44:821-844. [PMID: 33099619 PMCID: PMC7685786 DOI: 10.1093/femsre/fuaa034] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Accepted: 08/02/2020] [Indexed: 12/27/2022] Open
Abstract
Why do evolutionarily distinct microorganisms display similar physiological behaviours? Why are transitions from high-ATP yield to low(er)-ATP yield metabolisms so widespread across species? Why is fast growth generally accompanied with low stress tolerance? Do these regularities occur because most microbial species are subject to the same selective pressures and physicochemical constraints? If so, a broadly-applicable theory might be developed that predicts common microbiological behaviours. Microbial systems biologists have been working out the contours of this theory for the last two decades, guided by experimental data. At its foundations lie basic principles from evolutionary biology, enzyme biochemistry, metabolism, cell composition and steady-state growth. The theory makes predictions about fitness costs and benefits of protein expression, physicochemical constraints on cell growth and characteristics of optimal metabolisms that maximise growth rate. Comparisons of the theory with experimental data indicates that microorganisms often aim for maximisation of growth rate, also in the presence of stresses; they often express optimal metabolisms and metabolic proteins at optimal concentrations. This review explains the current status of the theory for microbiologists; its roots, predictions, experimental evidence and future directions.
Collapse
Affiliation(s)
- Frank J Bruggeman
- Systems Biology Lab, AIMMS, De Boelelaan 1108, 1081 HZ, VU University, Amsterdam, The Netherlands
| | - Robert Planqué
- Department of Mathematics, De Boelelaan 1111, 1081 HV, VU University, Amsterdam, The Netherlands
| | - Douwe Molenaar
- Systems Biology Lab, AIMMS, De Boelelaan 1108, 1081 HZ, VU University, Amsterdam, The Netherlands
| | - Bas Teusink
- Systems Biology Lab, AIMMS, De Boelelaan 1108, 1081 HZ, VU University, Amsterdam, The Netherlands
| |
Collapse
|
49
|
Proteomic Adaptation of Clostridioides difficile to Treatment with the Antimicrobial Peptide Nisin. Cells 2021; 10:cells10020372. [PMID: 33670309 PMCID: PMC7918085 DOI: 10.3390/cells10020372] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2020] [Revised: 02/04/2021] [Accepted: 02/08/2021] [Indexed: 01/05/2023] Open
Abstract
Clostridioides difficile is the leading cause of antibiotic-associated diarrhea but can also result in more serious, life-threatening conditions. The incidence of C. difficile infections in hospitals is increasing, both in frequency and severity, and antibiotic-resistant C. difficile strains are advancing. Against this background antimicrobial peptides (AMPs) are an interesting alternative to classic antibiotics. Information on the effects of AMPs on C. difficile will not only enhance the knowledge for possible biomedical application but may also provide insights into mechanisms of C. difficile to adapt or counteract AMPs. This study applies state-of-the-art mass spectrometry methods to quantitatively investigate the proteomic response of C. difficile 630∆erm to sublethal concentrations of the AMP nisin allowing to follow the cellular stress adaptation in a time-resolved manner. The results do not only point at a heavy reorganization of the cellular envelope but also resulted in pronounced changes in central cellular processes such as carbohydrate metabolism. Further, the number of flagella per cell was increased during the adaptation process. The potential involvement of flagella in nisin adaptation was supported by a more resistant phenotype exhibited by a non-motile but hyper-flagellated mutant.
Collapse
|
50
|
Regueira A, Lema JM, Mauricio-Iglesias M. Microbial inefficient substrate use through the perspective of resource allocation models. Curr Opin Biotechnol 2021; 67:130-140. [PMID: 33540363 DOI: 10.1016/j.copbio.2021.01.015] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Revised: 01/19/2021] [Accepted: 01/20/2021] [Indexed: 01/15/2023]
Abstract
Microorganisms extract energy from substrates following strategies that may seem suboptimal at first glance. Beyond the so-called yield-rate trade-off, resource allocation models, which focus on assigning different functional roles to the limited number of enzymes that a cell can support, offer a framework to interpret the inefficient substrate use by microorganisms. We review here relevant examples of substrate conversions where a significant part of the available energy is not utilised and how resource allocation models offer a mechanistic interpretation thereof, notably for open mixed cultures. Future developments are identified, in particular, the challenge of considering metabolic flexibility towards uncertain environmental changes instead of strict fixed optimality objectives, with the final goal of increasing the prediction capabilities of resource allocation models. Finally, we highlight the relevance of resource allocation to understand and enable a promising biorefinery platform revolving around lactate, which would increase the flexibility of waste-to-chemical biorefinery schemes.
Collapse
Affiliation(s)
- Alberte Regueira
- CRETUS Institute, Department of Chemical Engineering, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain.
| | - Juan M Lema
- CRETUS Institute, Department of Chemical Engineering, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | - Miguel Mauricio-Iglesias
- CRETUS Institute, Department of Chemical Engineering, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| |
Collapse
|