51
|
Multiscale models quantifying yeast physiology: towards a whole-cell model. Trends Biotechnol 2021; 40:291-305. [PMID: 34303549 DOI: 10.1016/j.tibtech.2021.06.010] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2021] [Revised: 06/26/2021] [Accepted: 06/28/2021] [Indexed: 12/21/2022]
Abstract
The yeast Saccharomyces cerevisiae is widely used as a cell factory and as an important eukaryal model organism for studying cellular physiology related to human health and disease. Yeast was also the first eukaryal organism for which a genome-scale metabolic model (GEM) was developed. In recent years there has been interest in expanding the modeling framework for yeast by incorporating enzymatic parameters and other heterogeneous cellular networks to obtain a more comprehensive description of cellular physiology. We review the latest developments in multiscale models of yeast, and illustrate how a new generation of multiscale models could significantly enhance the predictive performance and expand the applications of classical GEMs in cell factory design and basic studies of yeast physiology.
Collapse
|
52
|
Kumar P, Adamczyk PA, Zhang X, Andrade RB, Romero PA, Ramanathan P, Reed JL. Active and machine learning-based approaches to rapidly enhance microbial chemical production. Metab Eng 2021; 67:216-226. [PMID: 34229079 DOI: 10.1016/j.ymben.2021.06.009] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2021] [Revised: 05/12/2021] [Accepted: 06/29/2021] [Indexed: 10/20/2022]
Abstract
In order to make renewable fuels and chemicals from microbes, new methods are required to engineer microbes more intelligently. Computational approaches, to engineer strains for enhanced chemical production typically rely on detailed mechanistic models (e.g., kinetic/stoichiometric models of metabolism)-requiring many experimental datasets for their parameterization-while experimental methods may require screening large mutant libraries to explore the design space for the few mutants with desired behaviors. To address these limitations, we developed an active and machine learning approach (ActiveOpt) to intelligently guide experiments to arrive at an optimal phenotype with minimal measured datasets. ActiveOpt was applied to two separate case studies to evaluate its potential to increase valine yields and neurosporene productivity in Escherichia coli. In both the cases, ActiveOpt identified the best performing strain in fewer experiments than the case studies used. This work demonstrates that machine and active learning approaches have the potential to greatly facilitate metabolic engineering efforts to rapidly achieve its objectives.
Collapse
Affiliation(s)
- Prashant Kumar
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, 1415 Engineering Dr., Madison, WI, 53706, USA; ZS Associates, 1560 Sherman Ave, Evanston, IL, 60201, USA
| | - Paul A Adamczyk
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, 1415 Engineering Dr., Madison, WI, 53706, USA
| | - Xiaolin Zhang
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, 1415 Engineering Dr., Madison, WI, 53706, USA
| | - Ramon Bonela Andrade
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, 1415 Engineering Dr., Madison, WI, 53706, USA
| | - Philip A Romero
- Department of Biochemistry, University of Wisconsin-Madison, 440 Henry Mall, Madison, WI, 53706, USA
| | - Parameswaran Ramanathan
- Department of Electrical and Computer Engineering, University of Wisconsin-Madison, 1415 Engineering Dr., Madison, WI, 53706, USA.
| | - Jennifer L Reed
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, 1415 Engineering Dr., Madison, WI, 53706, USA
| |
Collapse
|
53
|
Covert MW, Gillies TE, Kudo T, Agmon E. A forecast for large-scale, predictive biology: Lessons from meteorology. Cell Syst 2021; 12:488-496. [PMID: 34139161 PMCID: PMC8217727 DOI: 10.1016/j.cels.2021.05.014] [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/16/2020] [Revised: 04/01/2021] [Accepted: 05/18/2021] [Indexed: 11/19/2022]
Abstract
Quantitative systems biology, in which predictive mathematical models are constructed to guide the design of experiments and predict experimental outcomes, is at an exciting transition point, where the foundational scientific principles are becoming established, but the impact is not yet global. The next steps necessary for mathematical modeling to transform biological research and applications, in the same way it has already transformed other fields, is not completely clear. The purpose of this perspective is to forecast possible answers to this question-what needs to happen next-by drawing on the experience gained in another field, specifically meteorology. We review here a number of lessons learned in weather prediction that are directly relevant to biological systems modeling, and that we believe can enable the same kinds of global impact in our field as atmospheric modeling makes today.
Collapse
Affiliation(s)
- Markus W Covert
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA.
| | - Taryn E Gillies
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
| | - Takamasa Kudo
- Department of Chemical and Systems Biology, Stanford University, Stanford, CA 94305, USA
| | - Eran Agmon
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
54
|
|
55
|
Hameri T, Fengos G, Hatzimanikatis V. The effects of model complexity and size on metabolic flux distribution and control: case study in Escherichia coli. BMC Bioinformatics 2021; 22:134. [PMID: 33743594 PMCID: PMC7981984 DOI: 10.1186/s12859-021-04066-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2020] [Accepted: 03/08/2021] [Indexed: 12/31/2022] Open
Abstract
Background Significant efforts have been made in building large-scale kinetic models of cellular metabolism in the past two decades. However, most kinetic models published to date, remain focused around central carbon pathways or are built around ad hoc reduced models without clear justification on their derivation and usage. Systematic algorithms exist for reducing genome-scale metabolic reconstructions to build thermodynamically feasible and consistently reduced stoichiometric models. However, it is important to study how network complexity affects conclusions derived from large-scale kinetic models built around consistently reduced models before we can apply them to study biological systems. Results We reduced the iJO1366 Escherichia Coli genome-scale metabolic reconstruction systematically to build three stoichiometric models of different size. Since the reduced models are expansions around the core subsystems for which the reduction was performed, the models are nested. We present a method for scaling up the flux profile and the concentration vector reference steady-states from the smallest model to the larger ones, whilst preserving maximum equivalency. Populations of kinetic models, preserving similarity in kinetic parameters, were built around the reference steady-states and their metabolic sensitivity coefficients (MSCs) were computed. The MSCs were sensitive to the model complexity. We proposed a metric for measuring the sensitivity of MSCs to these structural changes. Conclusions We proposed for the first time a workflow for scaling up the size of kinetic models while preserving equivalency between the kinetic models. Using this workflow, we demonstrate that model complexity in terms of networks size has significant impact on sensitivity characteristics of kinetic models. Therefore, it is essential to account for the effects of network complexity when constructing kinetic models. The presented metric for measuring MSC sensitivity to structural changes can guide modelers and experimentalists in improving model quality and guide synthetic biology and metabolic engineering. Our proposed workflow enables the testing of the suitability of a kinetic model for answering certain study-specific questions. We argue that the model-based metabolic design targets that are common across models of different size are of higher confidence, while those that are different could be the objective of investigations for model improvement. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-021-04066-y.
Collapse
Affiliation(s)
- Tuure Hameri
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), 1015, Lausanne, Switzerland
| | - Georgios Fengos
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), 1015, Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), 1015, Lausanne, Switzerland.
| |
Collapse
|
56
|
Porokhin V, Amin SA, Nicks TB, Gopinarayanan VE, Nair NU, Hassoun S. Analysis of metabolic network disruption in engineered microbial hosts due to enzyme promiscuity. Metab Eng Commun 2021; 12:e00170. [PMID: 33850714 PMCID: PMC8039717 DOI: 10.1016/j.mec.2021.e00170] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2020] [Revised: 01/22/2021] [Accepted: 03/01/2021] [Indexed: 11/30/2022] Open
Abstract
Increasing understanding of metabolic and regulatory networks underlying microbial physiology has enabled creation of progressively more complex synthetic biological systems for biochemical, biomedical, agricultural, and environmental applications. However, despite best efforts, confounding phenotypes still emerge from unforeseen interplay between biological parts, and the design of robust and modular biological systems remains elusive. Such interactions are difficult to predict when designing synthetic systems and may manifest during experimental testing as inefficiencies that need to be overcome. Transforming organisms such as Escherichia coli into microbial factories is achieved via several engineering strategies, used individually or in combination, with the goal of maximizing the production of chosen target compounds. One technique relies on suppressing or overexpressing selected genes; another involves introducing heterologous enzymes into a microbial host. These modifications steer mass flux towards the set of desired metabolites but may create unexpected interactions. In this work, we develop a computational method, termed Metabolic Disruption Workflow (MDFlow), for discovering interactions and network disruptions arising from enzyme promiscuity – the ability of enzymes to act on a wide range of molecules that are structurally similar to their native substrates. We apply MDFlow to two experimentally verified cases where strains with essential genes knocked out are rescued by interactions resulting from overexpression of one or more other genes. We demonstrate how enzyme promiscuity may aid cells in adapting to disruptions of essential metabolic functions. We then apply MDFlow to predict and evaluate a number of putative promiscuous reactions that can interfere with two heterologous pathways designed for 3-hydroxypropionic acid (3-HP) production. Using MDFlow, we can identify putative enzyme promiscuity and the subsequent formation of unintended and undesirable byproducts that are not only disruptive to the host metabolism but also to the intended end-objective of high biosynthetic productivity and yield. As we demonstrate, MDFlow provides an innovative workflow to systematically identify incompatibilities between the native metabolism of the host and its engineered modifications due to enzyme promiscuity. Engineering modifications to cellular hosts result in undesirable byproducts. Metabolic Disruption: changes in engineered host due to enzyme promiscuity. Metabolic Disruption Workflow (MDFlow) uncovers metabolic disruption. MDFlow corroborates previously experimentally verified promiscuous interactions. MDFlow compares disruption due to heterologous pathways targeting 3-HP production.
Collapse
Affiliation(s)
| | - Sara A Amin
- Department of Computer Science, Tufts University, Medford, MA, USA
| | - Trevor B Nicks
- Department of Chemical and Biological Engineering, Tufts University, Medford, MA, USA
| | | | - Nikhil U Nair
- Department of Chemical and Biological Engineering, Tufts University, Medford, MA, USA
| | - Soha Hassoun
- Department of Computer Science, Tufts University, Medford, MA, USA.,Department of Chemical and Biological Engineering, Tufts University, Medford, MA, USA
| |
Collapse
|
57
|
Modelling Cell Metabolism: A Review on Constraint-Based Steady-State and Kinetic Approaches. Processes (Basel) 2021. [DOI: 10.3390/pr9020322] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.
Collapse
|
58
|
Seaver SMD, Liu F, Zhang Q, Jeffryes J, Faria JP, Edirisinghe JN, Mundy M, Chia N, Noor E, Beber M, Best AA, DeJongh M, Kimbrel JA, D’haeseleer P, McCorkle SR, Bolton JR, Pearson E, Canon S, Wood-Charlson EM, Cottingham RW, Arkin AP, Henry CS. The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes. Nucleic Acids Res 2021; 49:D575-D588. [PMID: 32986834 PMCID: PMC7778927 DOI: 10.1093/nar/gkaa746] [Citation(s) in RCA: 88] [Impact Index Per Article: 29.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2020] [Revised: 08/25/2020] [Accepted: 09/24/2020] [Indexed: 12/31/2022] Open
Abstract
For over 10 years, ModelSEED has been a primary resource for the construction of draft genome-scale metabolic models based on annotated microbial or plant genomes. Now being released, the biochemistry database serves as the foundation of biochemical data underlying ModelSEED and KBase. The biochemistry database embodies several properties that, taken together, distinguish it from other published biochemistry resources by: (i) including compartmentalization, transport reactions, charged molecules and proton balancing on reactions; (ii) being extensible by the user community, with all data stored in GitHub; and (iii) design as a biochemical 'Rosetta Stone' to facilitate comparison and integration of annotations from many different tools and databases. The database was constructed by combining chemical data from many resources, applying standard transformations, identifying redundancies and computing thermodynamic properties. The ModelSEED biochemistry is continually tested using flux balance analysis to ensure the biochemical network is modeling-ready and capable of simulating diverse phenotypes. Ontologies can be designed to aid in comparing and reconciling metabolic reconstructions that differ in how they represent various metabolic pathways. ModelSEED now includes 33,978 compounds and 36,645 reactions, available as a set of extensible files on GitHub, and available to search at https://modelseed.org/biochem and KBase.
Collapse
Affiliation(s)
- Samuel M D Seaver
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Filipe Liu
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Qizhi Zhang
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - James Jeffryes
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - José P Faria
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Janaka N Edirisinghe
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Michael Mundy
- Center for Individualized Medicine, Mayo Clinic, Rochester, MN 55905, USA
| | - Nicholas Chia
- Center for Individualized Medicine, Mayo Clinic, Rochester, MN 55905, USA
| | - Elad Noor
- Department of Biology, Institute of Molecular Systems Biology, Eidgenössische Technische Hochschule Zürich, CH-8093 Zürich, Switzerland
| | - Moritz E Beber
- Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Kongens Lyngby, 2800, Denmark
| | - Aaron A Best
- Department of Biology, Hope College, Holland, MI 49423, USA
| | - Matthew DeJongh
- Department of Computer Science, Hope College, Holland, MI 49423, USA
| | - Jeffrey A Kimbrel
- Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA
| | - Patrik D’haeseleer
- Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA
| | - Sean R McCorkle
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - Jay R Bolton
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Erik Pearson
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Shane Canon
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Elisha M Wood-Charlson
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Robert W Cottingham
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| | - Adam P Arkin
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Christopher S Henry
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| |
Collapse
|
59
|
Seaver SMD, Liu F, Zhang Q, Jeffryes J, Faria JP, Edirisinghe JN, Mundy M, Chia N, Noor E, Beber ME, Best AA, DeJongh M, Kimbrel JA, D'haeseleer P, McCorkle SR, Bolton JR, Pearson E, Canon S, Wood-Charlson EM, Cottingham RW, Arkin AP, Henry CS. The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes. Nucleic Acids Res 2021; 49:D1555. [PMID: 33179751 DOI: 10.1101/2020.03.31.018663] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023] Open
Abstract
ABSTRACTFor over ten years, ModelSEED has been a primary resource for the construction of draft genome-scale metabolic models based on annotated microbial or plant genomes. Now being released, the biochemistry database serves as the foundation of biochemical data underlying ModelSEED and KBase. The biochemistry database embodies several properties that, taken together, distinguish it from other published biochemistry resources by: (i) including compartmentalization, transport reactions, charged molecules and proton balancing on reactions;; (ii) being extensible by the user community, with all data stored in GitHub; and (iii) design as a biochemical “Rosetta Stone” to facilitate comparison and integration of annotations from many different tools and databases. The database was constructed by combining chemical data from many resources, applying standard transformations, identifying redundancies, and computing thermodynamic properties. The ModelSEED biochemistry is continually tested using flux balance analysis to ensure the biochemical network is modeling-ready and capable of simulating diverse phenotypes. Ontologies can be designed to aid in comparing and reconciling metabolic reconstructions that differ in how they represent various metabolic pathways. ModelSEED now includes 33,978 compounds and 36,645 reactions, available as a set of extensible files on GitHub, and available to search at https://modelseed.org and KBase.
Collapse
Affiliation(s)
- Samuel M D Seaver
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Filipe Liu
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Qizhi Zhang
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - James Jeffryes
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - José P Faria
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Janaka N Edirisinghe
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Michael Mundy
- Center for Individualized Medicine, Mayo Clinic, Rochester, MN 55905, USA
| | - Nicholas Chia
- Center for Individualized Medicine, Mayo Clinic, Rochester, MN 55905, USA
| | - Elad Noor
- Department of Biology, Institute of Molecular Systems Biology, Eidgenössische Technische Hochschule Zürich, CH-8093 Zürich, Switzerland
| | - Moritz E Beber
- Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Kongens Lyngby, 2800, Denmark
| | - Aaron A Best
- Department of Biology, Hope College, Holland, MI 49423, USA
| | - Matthew DeJongh
- Department of Computer Science, Hope College, Holland, MI 49423, USA
| | - Jeffrey A Kimbrel
- Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA
| | - Patrik D'haeseleer
- Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA
| | - Sean R McCorkle
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - Jay R Bolton
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Erik Pearson
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Shane Canon
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Elisha M Wood-Charlson
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Robert W Cottingham
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| | - Adam P Arkin
- Environmental Genomics and Systems Biology Division, E.O. Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Christopher S Henry
- Computing, Environment, and Life Sciences Division, Argonne National Laboratory, Lemont, IL 60439, USA
| |
Collapse
|
60
|
Alsiyabi A, Solis AG, Cahoon EB, Saha R. Dissecting the regulatory roles of ORM proteins in the sphingolipid pathway of plants. PLoS Comput Biol 2021; 17:e1008284. [PMID: 33507896 PMCID: PMC7872301 DOI: 10.1371/journal.pcbi.1008284] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Revised: 02/09/2021] [Accepted: 12/14/2020] [Indexed: 01/05/2023] Open
Abstract
Sphingolipids are a vital component of plant cellular endomembranes and carry out multiple functional and regulatory roles. Different sphingolipid species confer rigidity to the membrane structure, facilitate trafficking of secretory proteins, and initiate programmed cell death. Although the regulation of the sphingolipid pathway is yet to be uncovered, increasing evidence has pointed to orosomucoid proteins (ORMs) playing a major regulatory role and potentially interacting with a number of components in the pathway, including both enzymes and sphingolipids. However, experimental exploration of new regulatory interactions is time consuming and often infeasible. In this work, a computational approach was taken to address this challenge. A metabolic network of the sphingolipid pathway in plants was reconstructed. The steady-state rates of reactions in the network were then determined through measurements of growth and cellular composition of the different sphingolipids in Arabidopsis seedlings. The Ensemble modeling framework was modified to accurately account for activation mechanisms and subsequently used to generate sets of kinetic parameters that converge to the measured steady-state fluxes in a thermodynamically consistent manner. In addition, the framework was appended with an additional module to automate screening the parameters and to output models consistent with previously reported network responses to different perturbations. By analyzing the network's response in the presence of different combinations of regulatory mechanisms, the model captured the experimentally observed repressive effect of ORMs on serine palmitoyltransferase (SPT). Furthermore, predictions point to a second regulatory role of ORM proteins, namely as an activator of class II (or LOH1 and LOH3) ceramide synthases. This activating role was found to be modulated by the concentration of free ceramides, where an accumulation of these sphingolipid species dampened the activating effect of ORMs on ceramide synthase. The predictions pave the way for future guided experiments and have implications in engineering crops with higher biotic stress tolerance.
Collapse
Affiliation(s)
- Adil Alsiyabi
- Department of Chemical and Biomolecular Engineering, University of Nebraska-Lincoln, Lincoln, Nebraska, United States of America
| | - Ariadna Gonzalez Solis
- Center for Plant Science Innovation & Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, Nebraska, United States of America
| | - Edgar B. Cahoon
- Center for Plant Science Innovation & Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, Nebraska, United States of America
| | - Rajib Saha
- Department of Chemical and Biomolecular Engineering, University of Nebraska-Lincoln, Lincoln, Nebraska, United States of America
- Center for Plant Science Innovation, University of Nebraska-Lincoln, Lincoln, Nebraska, United States of America
| |
Collapse
|
61
|
Foster CJ, Wang L, Dinh HV, Suthers PF, Maranas CD. Building kinetic models for metabolic engineering. Curr Opin Biotechnol 2020; 67:35-41. [PMID: 33360621 DOI: 10.1016/j.copbio.2020.11.010] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2020] [Revised: 11/09/2020] [Accepted: 11/22/2020] [Indexed: 12/16/2022]
Abstract
Kinetic formalisms of metabolism link metabolic fluxes to enzyme levels, metabolite concentrations and their allosteric regulatory interactions. Though they require the identification of physiologically relevant values for numerous parameters, kinetic formalisms uniquely establish a mechanistic link across heterogeneous omics datasets and provide an overarching vantage point to effectively inform metabolic engineering strategies. Advances in computational power, gene annotation coverage, and formalism standardization have led to significant progress over the past few years. However, careful interpretation of model predictions, limited metabolic flux datasets, and assessment of parameter sensitivity remain as challenges. In this review we highlight fundamental considerations which influence model quality and prediction, advances in methodologies, and success stories of deploying kinetic models to guide metabolic engineering.
Collapse
Affiliation(s)
- Charles J Foster
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA
| | - Lin Wang
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA
| | - Hoang V Dinh
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA; DOE Center for Advanced Bioenergy and Bioproducts Innovation, The Pennsylvania State University, Univesity Park, PA, USA
| | - Patrick F Suthers
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA; DOE Center for Advanced Bioenergy and Bioproducts Innovation, The Pennsylvania State University, Univesity Park, PA, USA
| | - Costas D Maranas
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA.
| |
Collapse
|
62
|
Otero-Muras I, Carbonell P. Automated engineering of synthetic metabolic pathways for efficient biomanufacturing. Metab Eng 2020; 63:61-80. [PMID: 33316374 DOI: 10.1016/j.ymben.2020.11.012] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2020] [Revised: 11/15/2020] [Accepted: 11/20/2020] [Indexed: 12/19/2022]
Abstract
Metabolic engineering involves the engineering and optimization of processes from single-cell to fermentation in order to increase production of valuable chemicals for health, food, energy, materials and others. A systems approach to metabolic engineering has gained traction in recent years thanks to advances in strain engineering, leading to an accelerated scaling from rapid prototyping to industrial production. Metabolic engineering is nowadays on track towards a truly manufacturing technology, with reduced times from conception to production enabled by automated protocols for DNA assembly of metabolic pathways in engineered producer strains. In this review, we discuss how the success of the metabolic engineering pipeline often relies on retrobiosynthetic protocols able to identify promising production routes and dynamic regulation strategies through automated biodesign algorithms, which are subsequently assembled as embedded integrated genetic circuits in the host strain. Those approaches are orchestrated by an experimental design strategy that provides optimal scheduling planning of the DNA assembly, rapid prototyping and, ultimately, brings forward an accelerated Design-Build-Test-Learn cycle and the overall optimization of the biomanufacturing process. Achieving such a vision will address the increasingly compelling demand in our society for delivering valuable biomolecules in an affordable, inclusive and sustainable bioeconomy.
Collapse
Affiliation(s)
- Irene Otero-Muras
- BioProcess Engineering Group, IIM-CSIC, Spanish National Research Council, Vigo, 36208, Spain.
| | - Pablo Carbonell
- Institute of Industrial Control Systems and Computing (ai2), Universitat Politècnica de València, 46022, Spain.
| |
Collapse
|
63
|
Suthers PF, Foster CJ, Sarkar D, Wang L, Maranas CD. Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms. Metab Eng 2020; 63:13-33. [PMID: 33310118 DOI: 10.1016/j.ymben.2020.11.013] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2020] [Revised: 11/13/2020] [Accepted: 11/27/2020] [Indexed: 12/16/2022]
Abstract
Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.
Collapse
Affiliation(s)
- Patrick F Suthers
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA; DOE Center for Advanced Bioenergy and Bioproducts Innovation, The Pennsylvania State University, University Park, PA, USA
| | - Charles J Foster
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA
| | - Debolina Sarkar
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA
| | - Lin Wang
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA
| | - Costas D Maranas
- Department of Chemical Engineering, The Pennsylvania State University, University Park, PA, USA; DOE Center for Advanced Bioenergy and Bioproducts Innovation, The Pennsylvania State University, University Park, PA, USA.
| |
Collapse
|
64
|
Novel allosteric inhibition of phosphoribulokinase identified by ensemble kinetic modeling of Synechocystis sp. PCC 6803 metabolism. Metab Eng Commun 2020; 11:e00153. [PMID: 33312875 PMCID: PMC7721636 DOI: 10.1016/j.mec.2020.e00153] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2020] [Revised: 10/17/2020] [Accepted: 11/17/2020] [Indexed: 12/11/2022] Open
Abstract
The present study attempted a computer simulation of the metabolism of a model cyanobacteria, Synechocystis sp. PCC 6803 (PCC 6803) to predict allosteric inhibitions that are likely to occur in photoautotrophic and mixotrophic conditions as well as in a metabolically engineered strain. PCC 6803 is a promising host for direct biochemical production from CO2; however, further investigation of allosteric regulation is required for rational metabolic engineering to produce target compounds. Herein, ensemble modeling of microbial metabolism was applied to build accurate predictive models by synthesizing the results of multiple models with different parameter sets into a single score to identify plausible allosteric inhibitions. The data driven-computer simulation using metabolic flux, enzyme abundance, and metabolite concentration data successfully identified candidates for allosteric inhibition. The enzyme assay experiment using the recombinant protein confirmed isocitrate was a non-competitive inhibitor of phosphoribulokinase as a novel allosteric regulation of cyanobacteria metabolism.
Collapse
|
65
|
Antoniewicz MR. A guide to metabolic flux analysis in metabolic engineering: Methods, tools and applications. Metab Eng 2020; 63:2-12. [PMID: 33157225 DOI: 10.1016/j.ymben.2020.11.002] [Citation(s) in RCA: 60] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2020] [Revised: 10/28/2020] [Accepted: 11/01/2020] [Indexed: 12/22/2022]
Abstract
The field of metabolic engineering is primarily concerned with improving the biological production of value-added chemicals, fuels and pharmaceuticals through the design, construction and optimization of metabolic pathways, redirection of intracellular fluxes, and refinement of cellular properties relevant for industrial bioprocess implementation. Metabolic network models and metabolic fluxes are central concepts in metabolic engineering, as was emphasized in the first paper published in this journal, "Metabolic fluxes and metabolic engineering" (Metabolic Engineering, 1: 1-11, 1999). In the past two decades, a wide range of computational, analytical and experimental approaches have been developed to interrogate the capabilities of biological systems through analysis of metabolic network models using techniques such as flux balance analysis (FBA), and quantify metabolic fluxes using constrained-based modeling approaches such as metabolic flux analysis (MFA) and more advanced experimental techniques based on the use of stable-isotope tracers, i.e. 13C-metabolic flux analysis (13C-MFA). In this review, we describe the basic principles of metabolic flux analysis, discuss current best practices in flux quantification, highlight potential pitfalls and alternative approaches in the application of these tools, and give a broad overview of pragmatic applications of flux analysis in metabolic engineering practice.
Collapse
Affiliation(s)
- Maciek R Antoniewicz
- Department of Chemical Engineering, Metabolic Engineering and Systems Biology Laboratory, University of Michigan, Ann Arbor, MI, 48109, USA.
| |
Collapse
|
66
|
Liu Y, Su A, Li J, Ledesma-Amaro R, Xu P, Du G, Liu L. Towards next-generation model microorganism chassis for biomanufacturing. Appl Microbiol Biotechnol 2020; 104:9095-9108. [DOI: 10.1007/s00253-020-10902-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2020] [Revised: 09/03/2020] [Accepted: 09/10/2020] [Indexed: 11/29/2022]
|
67
|
Küken A, Gennermann K, Nikoloski Z. Characterization of maximal enzyme catalytic rates in central metabolism of Arabidopsis thaliana. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2020; 103:2168-2177. [PMID: 32656814 DOI: 10.1111/tpj.14890] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Revised: 05/06/2020] [Accepted: 06/09/2020] [Indexed: 06/11/2023]
Abstract
Availability of plant-specific enzyme kinetic data is scarce, limiting the predictive power of metabolic models and precluding identification of genetic factors of enzyme properties. Enzyme kinetic data are measured in vitro, often under non-physiological conditions, and conclusions elicited from modeling warrant caution. Here we estimate maximal in vivo catalytic rates for 168 plant enzymes, including photosystems I and II, cytochrome-b6f complex, ATP-citrate synthase, sucrose-phosphate synthase as well as enzymes from amino acid synthesis with previously undocumented enzyme kinetic data in BRENDA. The estimations are obtained by integrating condition-specific quantitative proteomics data, maximal rates of selected enzymes, growth measurements from Arabidopsis thaliana rosette with and fluxes through canonical pathways in a constraint-based model of leaf metabolism. In comparison to findings in Escherichia coli, we demonstrate weaker concordance between the plant-specific in vitro and in vivo enzyme catalytic rates due to a low degree of enzyme saturation. This is supported by the finding that concentrations of nicotinamide adenine dinucleotide (phosphate), adenosine triphosphate and uridine triphosphate, calculated based on our maximal in vivo catalytic rates, and available quantitative metabolomics data are below reported KM values and, therefore, indicate undersaturation of respective enzymes. Our findings show that genome-wide profiling of enzyme kinetic properties is feasible in plants, paving the way for understanding resource allocation.
Collapse
Affiliation(s)
- Anika Küken
- System Biology and Mathematical Modeling Group, Max Planck Institute of Molecular Plant Physiology, Am Mühlenberg 1, Potsdam-Golm, Germany
- Bioinformatics Group, Institute of Biochemistry and Biology, University of Potsdam, Karl-Liebknecht-Str. 24-25, Potsdam-Golm, Germany
| | - Kristin Gennermann
- Bioinformatics Group, Institute of Biochemistry and Biology, University of Potsdam, Karl-Liebknecht-Str. 24-25, Potsdam-Golm, Germany
| | - Zoran Nikoloski
- System Biology and Mathematical Modeling Group, Max Planck Institute of Molecular Plant Physiology, Am Mühlenberg 1, Potsdam-Golm, Germany
- Bioinformatics Group, Institute of Biochemistry and Biology, University of Potsdam, Karl-Liebknecht-Str. 24-25, Potsdam-Golm, Germany
| |
Collapse
|
68
|
Kinetic profiling of metabolic specialists demonstrates stability and consistency of in vivo enzyme turnover numbers. Proc Natl Acad Sci U S A 2020; 117:23182-23190. [PMID: 32873645 PMCID: PMC7502767 DOI: 10.1073/pnas.2001562117] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Enzyme turnover numbers (k cats) are essential for a quantitative understanding of cells. Because k cats are traditionally measured in low-throughput assays, they can be inconsistent, labor-intensive to obtain, and can miss in vivo effects. We use a data-driven approach to estimate in vivo k cats using metabolic specialist Escherichia coli strains that resulted from gene knockouts in central metabolism followed by metabolic optimization via laboratory evolution. By combining absolute proteomics with fluxomics data, we find that in vivo k cats are robust against genetic perturbations, suggesting that metabolic adaptation to gene loss is mostly achieved through other mechanisms, like gene-regulatory changes. Combining machine learning and genome-scale metabolic models, we show that the obtained in vivo k cats predict unseen proteomics data with much higher precision than in vitro k cats. The results demonstrate that in vivo k cats can solve the problem of inconsistent and low-coverage parameterizations of genome-scale cellular models.
Collapse
|
69
|
Gopalakrishnan S, Dash S, Maranas C. K-FIT: An accelerated kinetic parameterization algorithm using steady-state fluxomic data. Metab Eng 2020; 61:197-205. [DOI: 10.1016/j.ymben.2020.03.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Revised: 01/31/2020] [Accepted: 03/02/2020] [Indexed: 12/12/2022]
|
70
|
Lo-Thong O, Charton P, Cadet XF, Grondin-Perez B, Saavedra E, Damour C, Cadet F. Identification of flux checkpoints in a metabolic pathway through white-box, grey-box and black-box modeling approaches. Sci Rep 2020; 10:13446. [PMID: 32778715 PMCID: PMC7417601 DOI: 10.1038/s41598-020-70295-5] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2019] [Accepted: 07/27/2020] [Indexed: 11/29/2022] Open
Abstract
Metabolic pathway modeling plays an increasing role in drug design by allowing better understanding of the underlying regulation and controlling networks in the metabolism of living organisms. However, despite rapid progress in this area, pathway modeling can become a real nightmare for researchers, notably when few experimental data are available or when the pathway is highly complex. Here, three different approaches were developed to model the second part of glycolysis of E. histolytica as an application example, and have succeeded in predicting the final pathway flux: one including detailed kinetic information (white-box), another with an added adjustment term (grey-box) and the last one using an artificial neural network method (black-box). Afterwards, each model was used for metabolic control analysis and flux control coefficient determination. The first two enzymes of this pathway are identified as the key enzymes playing a role in flux control. This study revealed the significance of the three methods for building suitable models adjusted to the available data in the field of metabolic pathway modeling, and could be useful to biologists and modelers.
Collapse
Affiliation(s)
- Ophélie Lo-Thong
- University of Paris, UMR_S1134, BIGR, Inserm, 75015, Paris, France.,DSIMB, UMR_S1134, BIGR, Inserm, Laboratory of Excellence GR-Ex, Faculty of Sciences and Technology, University of La Reunion, 97715, Saint-Denis, France
| | - Philippe Charton
- University of Paris, UMR_S1134, BIGR, Inserm, 75015, Paris, France.,DSIMB, UMR_S1134, BIGR, Inserm, Laboratory of Excellence GR-Ex, Faculty of Sciences and Technology, University of La Reunion, 97715, Saint-Denis, France
| | - Xavier F Cadet
- PEACCEL, Artificial Intelligence Department, 6 square Albin Cachot, box 42, 75013, Paris, France
| | - Brigitte Grondin-Perez
- LE2P, Laboratory of Energy, Electronics and Processes EA 4079, Faculty of Sciences and Technology, University of La Reunion, 97444, St Denis cedex, France
| | - Emma Saavedra
- Departamento de Bioquímica, Instituto Nacional de Cardiología Ignacio Chávez, 14080, Mexico City, Mexico
| | - Cédric Damour
- LE2P, Laboratory of Energy, Electronics and Processes EA 4079, Faculty of Sciences and Technology, University of La Reunion, 97444, St Denis cedex, France
| | - Frédéric Cadet
- University of Paris, UMR_S1134, BIGR, Inserm, 75015, Paris, France. .,DSIMB, UMR_S1134, BIGR, Inserm, Laboratory of Excellence GR-Ex, Faculty of Sciences and Technology, University of La Reunion, 97715, Saint-Denis, France.
| |
Collapse
|
71
|
Volkova S, Matos MRA, Mattanovich M, Marín de Mas I. Metabolic Modelling as a Framework for Metabolomics Data Integration and Analysis. Metabolites 2020; 10:E303. [PMID: 32722118 PMCID: PMC7465778 DOI: 10.3390/metabo10080303] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2020] [Revised: 07/08/2020] [Accepted: 07/22/2020] [Indexed: 01/05/2023] Open
Abstract
Metabolic networks are regulated to ensure the dynamic adaptation of biochemical reaction fluxes to maintain cell homeostasis and optimal metabolic fitness in response to endogenous and exogenous perturbations. To this end, metabolism is tightly controlled by dynamic and intricate regulatory mechanisms involving allostery, enzyme abundance and post-translational modifications. The study of the molecular entities involved in these complex mechanisms has been boosted by the advent of high-throughput technologies. The so-called omics enable the quantification of the different molecular entities at different system layers, connecting the genotype with the phenotype. Therefore, the study of the overall behavior of a metabolic network and the omics data integration and analysis must be approached from a holistic perspective. Due to the close relationship between metabolism and cellular phenotype, metabolic modelling has emerged as a valuable tool to decipher the underlying mechanisms governing cell phenotype. Constraint-based modelling and kinetic modelling are among the most widely used methods to study cell metabolism at different scales, ranging from cells to tissues and organisms. These approaches enable integrating metabolomic data, among others, to enhance model predictive capabilities. In this review, we describe the current state of the art in metabolic modelling and discuss future perspectives and current challenges in the field.
Collapse
Affiliation(s)
| | | | | | - Igor Marín de Mas
- The Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark; (S.V.); (M.R.A.M.); (M.M.)
| |
Collapse
|
72
|
Macklin DN, Ahn-Horst TA, Choi H, Ruggero NA, Carrera J, Mason JC, Sun G, Agmon E, DeFelice MM, Maayan I, Lane K, Spangler RK, Gillies TE, Paull ML, Akhter S, Bray SR, Weaver DS, Keseler IM, Karp PD, Morrison JH, Covert MW. Simultaneous cross-evaluation of heterogeneous E. coli datasets via mechanistic simulation. Science 2020; 369:eaav3751. [PMID: 32703847 PMCID: PMC7990026 DOI: 10.1126/science.aav3751] [Citation(s) in RCA: 76] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2018] [Revised: 10/28/2019] [Accepted: 05/26/2020] [Indexed: 12/24/2022]
Abstract
The extensive heterogeneity of biological data poses challenges to analysis and interpretation. Construction of a large-scale mechanistic model of Escherichia coli enabled us to integrate and cross-evaluate a massive, heterogeneous dataset based on measurements reported by various groups over decades. We identified inconsistencies with functional consequences across the data, including that the total output of the ribosomes and RNA polymerases described by data are not sufficient for a cell to reproduce measured doubling times, that measured metabolic parameters are neither fully compatible with each other nor with overall growth, and that essential proteins are absent during the cell cycle-and the cell is robust to this absence. Finally, considering these data as a whole leads to successful predictions of new experimental outcomes, in this case protein half-lives.
Collapse
Affiliation(s)
- Derek N Macklin
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Travis A Ahn-Horst
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Heejo Choi
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Nicholas A Ruggero
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
- Department of Chemical Engineering, Stanford University, Stanford, CA 94305, USA
| | - Javier Carrera
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - John C Mason
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Gwanggyu Sun
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Eran Agmon
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Mialy M DeFelice
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Inbal Maayan
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Keara Lane
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Ryan K Spangler
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Taryn E Gillies
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Morgan L Paull
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
| | - Sajia Akhter
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
| | - Samuel R Bray
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
| | | | | | | | - Jerry H Morrison
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Markus W Covert
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA.
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
73
|
Ye C, Luo Q, Guo L, Gao C, Xu N, Zhang L, Liu L, Chen X. Improving lysine production through construction of an Escherichia coli enzyme-constrained model. Biotechnol Bioeng 2020; 117:3533-3544. [PMID: 32648933 DOI: 10.1002/bit.27485] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Revised: 05/28/2020] [Accepted: 07/09/2020] [Indexed: 12/28/2022]
Abstract
Microbial cell factories are widely used for the production of high-value chemicals. However, maximizing production titers is made difficult by the complicated regulatory mechanisms of these cell platforms. Here, kcat values were incorporated to construct an Escherichia coli enzyme-constrained model. The resulting ec_iML1515 model showed that the protein demand and protein synthesis rate were the key factors affecting lysine production. By optimizing the expression of the 20 top-demanded proteins, lysine titers reached 95.7 ± 0.7 g/L, with a 0.45 g/g glucose yield. Moreover, adjusting NH4 + and dissolved oxygen levels to regulate the synthesis rate of energy metabolism-related proteins caused lysine titers and glucose yields to increase to 193.6 ± 1.8 g/L and 0.74 g/g, respectively. The ec_iML1515 model provides insight into how enzymes required for the biosynthesis of certain products are distributed between and within metabolic pathways. This information can be used to accurately predict and rationally design lysine production.
Collapse
Affiliation(s)
- Chao Ye
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, China.,School of Food Science and Pharmaceutical Engineering, Nanjing Normal University, Nanjing, China.,Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi, China
| | - Qiuling Luo
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, China.,Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi, China
| | - Liang Guo
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, China.,Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi, China
| | - Cong Gao
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, China.,Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi, China
| | - Nan Xu
- College of Bioscience and Biotechnology, Yangzhou University, Yangzhou, China
| | - Li Zhang
- School of Marine and Bioengineering, Yancheng Institute of Technology, Yancheng, China
| | - Liming Liu
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, China.,Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi, China
| | - Xiulai Chen
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, China.,Key Laboratory of Industrial Biotechnology, Ministry of Education, Jiangnan University, Wuxi, China
| |
Collapse
|
74
|
Garcia S, Trinh CT. Harnessing Natural Modularity of Metabolism with Goal Attainment Optimization to Design a Modular Chassis Cell for Production of Diverse Chemicals. ACS Synth Biol 2020; 9:1665-1681. [PMID: 32470305 DOI: 10.1021/acssynbio.9b00518] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Modular design is key to achieve efficient and robust systems across engineering disciplines. Modular design potentially offers advantages to engineer microbial systems for biocatalysis, bioremediation, and biosensing, overcoming the slow and costly design-build-test-learn cycles in the conventional cell engineering approach. These systems consist of a modular (chassis) cell compatible with exchangeable modules that enable programmed functions such as overproduction of a desirable chemical. We previously proposed a multiobjective optimization framework coupled with metabolic flux models to design modular cells and solved it using multiobjective evolutionary algorithms. However, such approach might not achieve solution optimality and hence limits design options for experimental implementation. In this study, we developed the goal attainment formulation compatible with optimization algorithms that guarantee solution optimality. We applied goal attainment to design an Escherichia coli modular cell capable of synthesizing all molecules in a biochemically diverse library at high yields and rates with only a few genetic manipulations. To elucidate modular organization of the designed cells, we developed a flux variance clustering (FVC) method by identifying reactions with high flux variance and clustering them to identify metabolic modules. Using FVC, we identified reaction usage patterns for different modules in the modular cell, revealing that its broad pathway compatibility is enabled by the natural modularity and flexible flux capacity of endogenous core metabolism. Overall, this study not only sheds light on modularity in metabolic networks from their topology and metabolic functions but also presents a useful synthetic biology toolbox to design modular cells with broad applications.
Collapse
Affiliation(s)
- Sergio Garcia
- Department of Chemical and Biomolecular Engineering, The University of Tennessee, Knoxville, Tennessee 37996, United States
- Center for Bioenergy Innovation, Oak Ridge National Laboratory Oak Ridge, Tennessee 37830, United States
| | - Cong T. Trinh
- Department of Chemical and Biomolecular Engineering, The University of Tennessee, Knoxville, Tennessee 37996, United States
- Center for Bioenergy Innovation, Oak Ridge National Laboratory Oak Ridge, Tennessee 37830, United States
| |
Collapse
|
75
|
Volk MJ, Lourentzou I, Mishra S, Vo LT, Zhai C, Zhao H. Biosystems Design by Machine Learning. ACS Synth Biol 2020; 9:1514-1533. [PMID: 32485108 DOI: 10.1021/acssynbio.0c00129] [Citation(s) in RCA: 52] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Biosystems such as enzymes, pathways, and whole cells have been increasingly explored for biotechnological applications. However, the intricate connectivity and resulting complexity of biosystems poses a major hurdle in designing biosystems with desirable features. As -omics and other high throughput technologies have been rapidly developed, the promise of applying machine learning (ML) techniques in biosystems design has started to become a reality. ML models enable the identification of patterns within complicated biological data across multiple scales of analysis and can augment biosystems design applications by predicting new candidates for optimized performance. ML is being used at every stage of biosystems design to help find nonobvious engineering solutions with fewer design iterations. In this review, we first describe commonly used models and modeling paradigms within ML. We then discuss some applications of these models that have already shown success in biotechnological applications. Moreover, we discuss successful applications at all scales of biosystems design, including nucleic acids, genetic circuits, proteins, pathways, genomes, and bioprocesses. Finally, we discuss some limitations of these methods and potential solutions as well as prospects of the combination of ML and biosystems design.
Collapse
|
76
|
Post MJ, Levenberg S, Kaplan DL, Genovese N, Fu J, Bryant CJ, Negowetti N, Verzijden K, Moutsatsou P. Scientific, sustainability and regulatory challenges of cultured meat. ACTA ACUST UNITED AC 2020. [DOI: 10.1038/s43016-020-0112-z] [Citation(s) in RCA: 140] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
77
|
Bromig L, Kremling A, Marin-Sanguino A. Understanding biochemical design principles with ensembles of canonical non-linear models. PLoS One 2020; 15:e0230599. [PMID: 32353072 PMCID: PMC7192416 DOI: 10.1371/journal.pone.0230599] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2020] [Accepted: 03/03/2020] [Indexed: 12/22/2022] Open
Abstract
Systems biology applies concepts from engineering in order to understand biological networks. If such an understanding was complete, biologists would be able to design ad hoc biochemical components tailored for different purposes, which is the goal of synthetic biology. Needless to say that we are far away from creating biological subsystems as intricate and precise as those found in nature, but mathematical models and high throughput techniques have brought us a long way in this direction. One of the difficulties that still needs to be overcome is finding the right values for model parameters and dealing with uncertainty, which is proving to be an extremely difficult task. In this work, we take advantage of ensemble modeling techniques, where a large number of models with different parameter values are formulated and then tested according to some performance criteria. By finding features shared by successful models, the role of different components and the synergies between them can be better understood. We will address some of the difficulties often faced by ensemble modeling approaches, such as the need to sample a space whose size grows exponentially with the number of parameters, and establishing useful selection criteria. Some methods will be shown to reduce the predictions from many models into a set of understandable “design principles” that can guide us to improve or manufacture a biochemical network. Our proposed framework formulates models within standard formalisms in order to integrate information from different sources and minimize the dimension of the parameter space. Additionally, the mathematical properties of the formalism enable a partition of the parameter space into independent subspaces. Each of these subspaces can be paired with a set of criteria that depend exclusively on it, thus allowing a separate sampling/screening in spaces of lower dimension. By applying tests in a strict order where computationally cheaper tests are applied first to each subspace and applying computationally expensive tests to the remaining subset thereafter, the use of resources is optimized and a larger number of models can be examined. This can be compared to a complex database query where the order of the requests can make a huge difference in the processing time. The method will be illustrated by analyzing a classical model of a metabolic pathway with end-product inhibition. Even for such a simple model, the method provides novel insight.
Collapse
Affiliation(s)
- Lukas Bromig
- Specialty Division for Systems Biotechnology, Technische Universität München, Garching, Germany
| | - Andreas Kremling
- Specialty Division for Systems Biotechnology, Technische Universität München, Garching, Germany
| | - Alberto Marin-Sanguino
- Specialty Division for Systems Biotechnology, Technische Universität München, Garching, Germany
- * E-mail:
| |
Collapse
|
78
|
Clark TJ, Guo L, Morgan J, Schwender J. Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models. ANNUAL REVIEW OF PLANT BIOLOGY 2020; 71:303-326. [PMID: 32017600 DOI: 10.1146/annurev-arplant-050718-100221] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For 13C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.
Collapse
Affiliation(s)
- Teresa J Clark
- Biology Department, Brookhaven National Laboratory, Upton, New York 11973, USA; ,
| | - Longyun Guo
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907, USA; ,
| | - John Morgan
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907, USA; ,
| | - Jorg Schwender
- Biology Department, Brookhaven National Laboratory, Upton, New York 11973, USA; ,
| |
Collapse
|
79
|
Yeo HC, Hong J, Lakshmanan M, Lee DY. Enzyme capacity-based genome scale modelling of CHO cells. Metab Eng 2020; 60:138-147. [PMID: 32330653 DOI: 10.1016/j.ymben.2020.04.005] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2019] [Revised: 03/21/2020] [Accepted: 04/14/2020] [Indexed: 10/24/2022]
Abstract
Chinese hamster ovary (CHO) cells are most prevalently used for producing recombinant therapeutics in biomanufacturing. Recently, more rational and systems approaches have been increasingly exploited to identify key metabolic bottlenecks and engineering targets for cell line engineering and process development based on the CHO genome-scale metabolic model which mechanistically characterizes cell culture behaviours. However, it is still challenging to quantify plausible intracellular fluxes and discern metabolic pathway usages considering various clonal traits and bioprocessing conditions. Thus, we newly incorporated enzyme kinetic information into the updated CHO genome-scale model (iCHO2291) and added enzyme capacity constraints within the flux balance analysis framework (ecFBA) to significantly reduce the flux variability in biologically meaningful manner, as such improving the accuracy of intracellular flux prediction. Interestingly, ecFBA could capture the overflow metabolism under the glucose excess condition where the usage of oxidative phosphorylation is limited by the enzyme capacity. In addition, its applicability was successfully demonstrated via a case study where the clone- and media-specific lactate metabolism was deciphered, suggesting that the lactate-pyruvate cycling could be beneficial for CHO cells to efficiently utilize the mitochondrial redox capacity. In summary, iCHO2296 with ecFBA can be used to confidently elucidate cell cultures and effectively identify key engineering targets, thus guiding bioprocess optimization and cell engineering efforts as a part of digital twin model for advanced biomanufacturing in future.
Collapse
Affiliation(s)
- Hock Chuan Yeo
- Bioprocessing Technology Institute, Agency for Science, Technology and Research (A*STAR), 20 Biopolis Way, #06-01, 138668, Singapore; Department of Chemical and Biomolecular Engineering, National University of Singapore, 4 Engineering Drive 4, 117585, Singapore
| | - Jongkwang Hong
- Bioprocessing Technology Institute, Agency for Science, Technology and Research (A*STAR), 20 Biopolis Way, #06-01, 138668, Singapore
| | - Meiyappan Lakshmanan
- Bioprocessing Technology Institute, Agency for Science, Technology and Research (A*STAR), 20 Biopolis Way, #06-01, 138668, Singapore.
| | - Dong-Yup Lee
- Bioprocessing Technology Institute, Agency for Science, Technology and Research (A*STAR), 20 Biopolis Way, #06-01, 138668, Singapore; School of Chemical Engineering, Sungkyunkwan University, 2066 Seobu-ro, Jangan-gu, Suwon, Gyeonggi-do, 16419, Republic of Korea.
| |
Collapse
|
80
|
Li G, Cao H, Xu Y. Structural and functional analyses of microbial metabolic networks reveal novel insights into genome-scale metabolic fluxes. Brief Bioinform 2020; 20:1590-1603. [PMID: 29596572 DOI: 10.1093/bib/bby022] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2018] [Revised: 03/01/2018] [Indexed: 11/13/2022] Open
Abstract
We present here an integrated analysis of structures and functions of genome-scale metabolic networks of 17 microorganisms. Our structural analyses of these networks revealed that the node degree of each network, represented as a (simplified) reaction network, follows a power-law distribution, and the clustering coefficient of each network has a positive correlation with the corresponding node degree. Together, these properties imply that each network has exactly one large and densely connected subnetwork or core. Further analyses revealed that each network consists of three functionally distinct subnetworks: (i) a core, consisting of a large number of directed reaction cycles of enzymes for interconversions among intermediate metabolites; (ii) a catabolic module, with a largely layered structure consisting of mostly catabolic enzymes; (iii) an anabolic module with a similar structure consisting of virtually all anabolic genes; and (iv) the three subnetworks cover on average ∼56, ∼31 and ∼13% of a network's nodes across the 17 networks, respectively. Functional analyses suggest: (1) cellular metabolic fluxes generally go from the catabolic module to the core for substantial interconversions, then the flux directions to anabolic module appear to be determined by input nutrient levels as well as a set of precursors needed for macromolecule syntheses; and (2) enzymes in each subnetwork have characteristic ranges of kinetic parameters, suggesting optimized metabolic and regulatory relationships among the three subnetworks.
Collapse
Affiliation(s)
- Gaoyang Li
- College of Computer Science and Technology, Jilin University, Changchun, Jilin, China.,Computational Systems Biology Lab, Department of Biochemistry and Molecular Biology and Institute of Bioinformatics, University of Georgia, Athens, GA, USA
| | - Huansheng Cao
- Computational Systems Biology Lab, Department of Biochemistry and Molecular Biology and Institute of Bioinformatics, University of Georgia, Athens, GA, USA.,The BESC BioEnergy Research Center, Oak Ridge National Lab, Oak Ridge, TN, USA
| | - Ying Xu
- College of Computer Science and Technology, Jilin University, Changchun, Jilin, China.,Computational Systems Biology Lab, Department of Biochemistry and Molecular Biology and Institute of Bioinformatics, University of Georgia, Athens, GA, USA.,The BESC BioEnergy Research Center, Oak Ridge National Lab, Oak Ridge, TN, USA.,School of Public Health, Jilin University, Changchun, Jilin, China
| |
Collapse
|
81
|
Abstract
The biological fitness of microbes is largely determined by the rate with which they replicate their biomass composition. Mathematical models that maximize this balanced growth rate while accounting for mass conservation, reaction kinetics, and limits on dry mass per volume are inevitably non-linear. Here, we develop a general theory for such models, termed Growth Balance Analysis (GBA), which provides explicit expressions for protein concentrations, fluxes, and growth rates. These variables are functions of the concentrations of cellular components, for which we calculate marginal fitness costs and benefits that are related to metabolic control coefficients. At maximal growth rate, the net benefits of all concentrations are equal. Based solely on physicochemical constraints, GBA unveils fundamental quantitative principles of cellular resource allocation and growth; it accurately predicts the relationship between growth rates and ribosome concentrations in E. coli and yeast and between growth rate and dry mass density in E. coli.
Collapse
Affiliation(s)
- Hugo Dourado
- Institute for Computer Science & Department of Biology, Heinrich Heine University, 40221, Düsseldorf, Germany
| | - Martin J Lercher
- Institute for Computer Science & Department of Biology, Heinrich Heine University, 40221, Düsseldorf, Germany.
| |
Collapse
|
82
|
Dromms RA, Lee JY, Styczynski MP. LK-DFBA: a linear programming-based modeling strategy for capturing dynamics and metabolite-dependent regulation in metabolism. BMC Bioinformatics 2020; 21:93. [PMID: 32122331 PMCID: PMC7053146 DOI: 10.1186/s12859-020-3422-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2019] [Accepted: 02/17/2020] [Indexed: 01/11/2023] Open
Abstract
BACKGROUND The systems-scale analysis of cellular metabolites, "metabolomics," provides data ideal for applications in metabolic engineering. However, many of the computational tools for strain design are built around Flux Balance Analysis (FBA), which makes assumptions that preclude direct integration of metabolomics data into the underlying models. Finding a way to retain the advantages of FBA's linear structure while relaxing some of its assumptions could allow us to account for metabolite levels and metabolite-dependent regulation in strain design tools built from FBA, improving the accuracy of predictions made by these tools. We designed, implemented, and characterized a modeling strategy based on Dynamic FBA (DFBA), called Linear Kinetics-Dynamic Flux Balance Analysis (LK-DFBA), to satisfy these specifications. Our strategy adds constraints describing the dynamics and regulation of metabolism that are strictly linear. We evaluated LK-DFBA against alternative modeling frameworks using simulated noisy data from a small in silico model and a larger model of central carbon metabolism in E. coli, and compared each framework's ability to recapitulate the original system. RESULTS In the smaller model, we found that we could use regression from a dynamic flux estimation (DFE) with an optional non-linear parameter optimization to reproduce metabolite concentration dynamic trends more effectively than an ordinary differential equation model with generalized mass action rate laws when tested under realistic data sampling frequency and noise levels. We observed detrimental effects across all tested modeling approaches when metabolite time course data were missing, but found these effects to be smaller for LK-DFBA in most cases. With the E. coli model, we produced qualitatively reasonable results with similar properties to the smaller model and explored two different parameterization structures that yield trade-offs in computation time and accuracy. CONCLUSIONS LK-DFBA allows for calculation of metabolite concentrations and considers metabolite-dependent regulation while still retaining many computational advantages of FBA. This provides the proof-of-principle for a new metabolic modeling framework with the potential to create genome-scale dynamic models and the potential to be applied in strain engineering tools that currently use FBA.
Collapse
Affiliation(s)
- Robert A Dromms
- School of Chemical & Biomolecular Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Justin Y Lee
- School of Chemical & Biomolecular Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Mark P Styczynski
- School of Chemical & Biomolecular Engineering, Georgia Institute of Technology, Atlanta, GA, USA.
| |
Collapse
|
83
|
Ye C, Xu N, Gao C, Liu G, Xu J, Zhang W, Chen X, Nielsen J, Liu L. Comprehensive understanding of Saccharomyces cerevisiae phenotypes with whole-cell model WM_S288C. Biotechnol Bioeng 2020; 117:1562-1574. [PMID: 32022245 DOI: 10.1002/bit.27298] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Revised: 01/28/2020] [Accepted: 02/03/2020] [Indexed: 02/01/2023]
Abstract
Biological network construction for Saccharomyces cerevisiae is a widely used approach for simulating phenotypes and designing cell factories. However, due to a complicated regulatory mechanism governing the translation of genotype to phenotype, precise prediction of phenotypes remains challenging. Here, we present WM_S288C, a computational whole-cell model that includes 15 cellular states and 26 cellular processes and which enables integrated analyses of physiological functions of Saccharomyces cerevisiae. Using WM_S288C to predict phenotypes of S. cerevisiae, the functions of 1140 essential genes were characterized and linked to phenotypes at five levels. During the cell cycle, the dynamic allocation of intracellular molecules could be tracked in real-time to simulate cell activities. Additionally, one-third of non-essential genes were identified to affect cell growth via regulating nucleotide concentrations. These results demonstrated the value of WM_S288C as a tool for understanding and investigating the phenotypes of S. cerevisiae.
Collapse
Affiliation(s)
- Chao Ye
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, Jiangsu, China.,Key Laboratory of Industrial Biotechnology, Jiangnan University, Ministry of Education, Wuxi, Jiangsu, China.,National Engineering Laboratory for Cereal Fermentation Technology, Jiangnan University, Wuxi, Jiangsu, China
| | - Nan Xu
- College of Bioscience and Biotechnology, Yangzhou University, Yangzhou, Jiangsu, China
| | - Cong Gao
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, Jiangsu, China.,Key Laboratory of Industrial Biotechnology, Jiangnan University, Ministry of Education, Wuxi, Jiangsu, China.,National Engineering Laboratory for Cereal Fermentation Technology, Jiangnan University, Wuxi, Jiangsu, China
| | - Gaoqiang Liu
- Hunan Provincial Key Laboratory for Forestry Biotechnology, Central South University of Forestry and Technology, Changsha, Hunan, China
| | - Jianzhong Xu
- Key Laboratory of Industrial Biotechnology, Jiangnan University, Ministry of Education, Wuxi, Jiangsu, China
| | - Weiguo Zhang
- Key Laboratory of Industrial Biotechnology, Jiangnan University, Ministry of Education, Wuxi, Jiangsu, China
| | - Xiulai Chen
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, Jiangsu, China.,Key Laboratory of Industrial Biotechnology, Jiangnan University, Ministry of Education, Wuxi, Jiangsu, China.,National Engineering Laboratory for Cereal Fermentation Technology, Jiangnan University, Wuxi, Jiangsu, China
| | - Jens Nielsen
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
| | - Liming Liu
- State Key Laboratory of Food Science and Technology, Jiangnan University, Wuxi, Jiangsu, China.,Key Laboratory of Industrial Biotechnology, Jiangnan University, Ministry of Education, Wuxi, Jiangsu, China.,National Engineering Laboratory for Cereal Fermentation Technology, Jiangnan University, Wuxi, Jiangsu, China
| |
Collapse
|
84
|
Pecht T, Aschenbrenner AC, Ulas T, Succurro A. Modeling population heterogeneity from microbial communities to immune response in cells. Cell Mol Life Sci 2020; 77:415-432. [PMID: 31768606 PMCID: PMC7010691 DOI: 10.1007/s00018-019-03378-w] [Citation(s) in RCA: 2] [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: 11/05/2019] [Revised: 11/05/2019] [Accepted: 11/12/2019] [Indexed: 12/14/2022]
Abstract
Heterogeneity is universally observed in all natural systems and across multiple scales. Understanding population heterogeneity is an intriguing and attractive topic of research in different disciplines, including microbiology and immunology. Microbes and mammalian immune cells present obviously rather different system-specific biological features. Nevertheless, as typically occurs in science, similar methods can be used to study both types of cells. This is particularly true for mathematical modeling, in which key features of a system are translated into algorithms to challenge our mechanistic understanding of the underlying biology. In this review, we first present a broad overview of the experimental developments that allowed observing heterogeneity at the single cell level. We then highlight how this "data revolution" requires the parallel advancement of algorithms and computing infrastructure for data processing and analysis, and finally present representative examples of computational models of population heterogeneity, from microbial communities to immune response in cells.
Collapse
Affiliation(s)
- Tal Pecht
- Genomics and Immunoregulation, Life and Medical Sciences (LIMES) Institute, University of Bonn, Bonn, Germany
| | - Anna C Aschenbrenner
- Genomics and Immunoregulation, Life and Medical Sciences (LIMES) Institute, University of Bonn, Bonn, Germany
- Department of Internal Medicine and Radboud Center for Infectious Diseases (RCI), Radboud University Medical Center, 6525, Nijmegen, The Netherlands
| | - Thomas Ulas
- Genomics and Immunoregulation, Life and Medical Sciences (LIMES) Institute, University of Bonn, Bonn, Germany
| | - Antonella Succurro
- Genomics and Immunoregulation, Life and Medical Sciences (LIMES) Institute, University of Bonn, Bonn, Germany.
- West German Genome Center (WGGC), University of Bonn, Bonn, Germany.
| |
Collapse
|
85
|
Chen Y, Banerjee D, Mukhopadhyay A, Petzold CJ. Systems and synthetic biology tools for advanced bioproduction hosts. Curr Opin Biotechnol 2020; 64:101-109. [PMID: 31927061 DOI: 10.1016/j.copbio.2019.12.007] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Revised: 11/27/2019] [Accepted: 12/08/2019] [Indexed: 02/07/2023]
Abstract
The genomic revolution ushered in an era of discovery and characterization of enzymes from novel organisms that fueled engineering of microbes to produce commodity and high-value compounds. Over the past decade advances in synthetic biology tools in recent years contributed to significant progress in metabolic engineering efforts to produce both biofuels and bioproducts resulting in several such related items being brought to market. These successes represent a burgeoning bio-economy; however, significant resources and time are still necessary to progress a system from proof-of-concept to market. In order to fully realize this potential, methods that examine biological systems in a comprehensive, systematic and high-throughput manner are essential. Recent success in synthetic biology has coincided with the development of systems biology and analytical approaches that kept pace and scaled with technology development. Here, we review a selection of systems biology methods and their use in synthetic biology approaches for microbial biotechnology platforms.
Collapse
Affiliation(s)
- Yan Chen
- Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, CA, USA; Agile BioFoundry, Lawrence Berkeley National Laboratory, Emeryville, CA, USA; Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Deepanwita Banerjee
- Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, CA, USA; Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Aindrila Mukhopadhyay
- Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, CA, USA; Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA; Environmental Genomics and Systems Biology Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Christopher J Petzold
- Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, CA, USA; Agile BioFoundry, Lawrence Berkeley National Laboratory, Emeryville, CA, USA; Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA.
| |
Collapse
|
86
|
de Groot DH, Hulshof J, Teusink B, Bruggeman FJ, Planqué R. Elementary Growth Modes provide a molecular description of cellular self-fabrication. PLoS Comput Biol 2020; 16:e1007559. [PMID: 31986156 PMCID: PMC7004393 DOI: 10.1371/journal.pcbi.1007559] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2019] [Revised: 02/06/2020] [Accepted: 11/22/2019] [Indexed: 02/04/2023] Open
Abstract
In this paper we try to describe all possible molecular states (phenotypes) for a cell that fabricates itself at a constant rate, given its enzyme kinetics and the stoichiometry of all reactions. For this, we must understand the process of cellular growth: steady-state self-fabrication requires a cell to synthesize all of its components, including metabolites, enzymes and ribosomes, in proportions that match its own composition. Simultaneously, the concentrations of these components affect the rates of metabolism and biosynthesis, and hence the growth rate. We here derive a theory that describes all phenotypes that solve this circular problem. All phenotypes can be described as a combination of minimal building blocks, which we call Elementary Growth Modes (EGMs). EGMs can be used as the theoretical basis for all models that explicitly model self-fabrication, such as the currently popular Metabolism and Expression models. We then use our theory to make concrete biological predictions. We find that natural selection for maximal growth rate drives microorganisms to states of minimal phenotypic complexity: only one EGM will be active when growth rate is maximised. The phenotype of a cell is only extended with one more EGM whenever growth becomes limited by an additional biophysical constraint, such as a limited solvent capacity of a cellular compartment. The theory presented here extends recent results on Elementary Flux Modes: the minimal building blocks of cellular growth models that lack the self-fabrication aspect. Our theory starts from basic biochemical and evolutionary considerations, and describes unicellular life, both in growth-promoting and in stress-inducing environments, in terms of EGMs.
Collapse
Affiliation(s)
- Daan H. de Groot
- Systems Bioinformatics, Amsterdam Institute for Molecules, Medicines & Systems, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Josephus Hulshof
- Department of Mathematics, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Bas Teusink
- Systems Bioinformatics, Amsterdam Institute for Molecules, Medicines & Systems, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Frank J. Bruggeman
- Systems Bioinformatics, Amsterdam Institute for Molecules, Medicines & Systems, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Robert Planqué
- Systems Bioinformatics, Amsterdam Institute for Molecules, Medicines & Systems, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
- Department of Mathematics, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| |
Collapse
|
87
|
Wang G, Haringa C, Tang W, Noorman H, Chu J, Zhuang Y, Zhang S. Coupled metabolic-hydrodynamic modeling enabling rational scale-up of industrial bioprocesses. Biotechnol Bioeng 2019; 117:844-867. [PMID: 31814101 DOI: 10.1002/bit.27243] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Revised: 11/28/2019] [Accepted: 11/30/2019] [Indexed: 12/13/2022]
Abstract
Metabolomics aims to address what and how regulatory mechanisms are coordinated to achieve flux optimality, different metabolic objectives as well as appropriate adaptations to dynamic nutrient availability. Recent decades have witnessed that the integration of metabolomics and fluxomics within the goal of synthetic biology has arrived at generating the desired bioproducts with improved bioconversion efficiency. Absolute metabolite quantification by isotope dilution mass spectrometry represents a functional readout of cellular biochemistry and contributes to the establishment of metabolic (structured) models required in systems metabolic engineering. In industrial practices, population heterogeneity arising from fluctuating nutrient availability frequently leads to performance losses, that is reduced commercial metrics (titer, rate, and yield). Hence, the development of more stable producers and more predictable bioprocesses can benefit from a quantitative understanding of spatial and temporal cell-to-cell heterogeneity within industrial bioprocesses. Quantitative metabolomics analysis and metabolic modeling applied in computational fluid dynamics (CFD)-assisted scale-down simulators that mimic industrial heterogeneity such as fluctuations in nutrients, dissolved gases, and other stresses can procure informative clues for coping with issues during bioprocessing scale-up. In previous studies, only limited insights into the hydrodynamic conditions inside the industrial-scale bioreactor have been obtained, which makes case-by-case scale-up far from straightforward. Tracking the flow paths of cells circulating in large-scale bioreactors is a highly valuable tool for evaluating cellular performance in production tanks. The "lifelines" or "trajectories" of cells in industrial-scale bioreactors can be captured using Euler-Lagrange CFD simulation. This novel methodology can be further coupled with metabolic (structured) models to provide not only a statistical analysis of cell lifelines triggered by the environmental fluctuations but also a global assessment of the metabolic response to heterogeneity inside an industrial bioreactor. For the future, the industrial design should be dependent on the computational framework, and this integration work will allow bioprocess scale-up to the industrial scale with an end in mind.
Collapse
Affiliation(s)
- Guan Wang
- State Key Laboratory of Bioreactor Engineering, East China University of Science and Technology, Shanghai, People's Republic of China
| | - Cees Haringa
- Transport Phenomena, Chemical Engineering Department, Delft University of Technology, Delft, The Netherlands.,DSM Biotechnology Center, Delft, The Netherlands
| | - Wenjun Tang
- DSM Biotechnology Center, Delft, The Netherlands
| | - Henk Noorman
- DSM Biotechnology Center, Delft, The Netherlands.,Bioprocess Engineering, Department of Biotechnology, Delft University of Technology, Delft, The Netherlands
| | - Ju Chu
- State Key Laboratory of Bioreactor Engineering, East China University of Science and Technology, Shanghai, People's Republic of China
| | - Yingping Zhuang
- State Key Laboratory of Bioreactor Engineering, East China University of Science and Technology, Shanghai, People's Republic of China
| | - Siliang Zhang
- State Key Laboratory of Bioreactor Engineering, East China University of Science and Technology, Shanghai, People's Republic of China
| |
Collapse
|
88
|
Horvath N, Vilkhovoy M, Wayman JA, Calhoun K, Swartz J, Varner JD. Toward a genome scale sequence specific dynamic model of cell-free protein synthesis in Escherichia coli. Metab Eng Commun 2019; 10:e00113. [PMID: 32280586 PMCID: PMC7136494 DOI: 10.1016/j.mec.2019.e00113] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2018] [Revised: 10/15/2019] [Accepted: 11/19/2019] [Indexed: 11/09/2022] Open
Abstract
In this study, we developed a dynamic mathematical model of E. coli cell-free protein synthesis (CFPS). Model parameters were estimated from a dataset consisting of glucose, organic acids, energy species, amino acids, and protein product, chloramphenicol acetyltransferase (CAT) measurements. The model was successfully trained to simulate these measurements, especially those of the central carbon metabolism. We then used the trained model to evaluate the performance, e.g., the yield and rates of protein production. CAT was produced with an energy efficiency of 12%, suggesting that the process could be further optimized. Reaction group knockouts showed that protein productivity was most sensitive to the oxidative phosphorylation and glycolysis/gluconeogenesis pathways. Amino acid biosynthesis was also important for productivity, while overflow metabolism and TCA cycle affected the overall system state. In addition, translation was more important to productivity than transcription. Finally, CAT production was robust to allosteric control, as were most of the predicted metabolite concentrations; the exceptions to this were the concentrations of succinate and malate, and to a lesser extent pyruvate and acetate, which varied from the measured values when allosteric control was removed. This study is the first to use kinetic modeling to predict dynamic protein production in a cell-free E. coli system, and could provide a foundation for genome scale, dynamic modeling of cell-free E. coli protein synthesis. Protein production is biphasic, powered initially by glucose and later by pyruvate. Protein is produced with an energy efficiency of only 12%. Protein productivity is most sensitive to oxidative phosphorylation and glycolysis. Protein production is robust to allosteric control.
Collapse
Affiliation(s)
- Nicholas Horvath
- Robert Frederick Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, NY, 14853, USA
| | - Michael Vilkhovoy
- Robert Frederick Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, NY, 14853, USA
| | - Joseph A Wayman
- School of Applied and Engineering Physics, Cornell University, Ithaca, NY, 14853, USA
| | - Kara Calhoun
- School of Chemical Engineering, Stanford University, Stanford, CA, 94395, USA
| | - James Swartz
- School of Chemical Engineering, Stanford University, Stanford, CA, 94395, USA
| | - Jeffrey D Varner
- Robert Frederick Smith School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, NY, 14853, USA
| |
Collapse
|
89
|
Hameri T, Boldi MO, Hatzimanikatis V. Statistical inference in ensemble modeling of cellular metabolism. PLoS Comput Biol 2019; 15:e1007536. [PMID: 31815929 PMCID: PMC6922442 DOI: 10.1371/journal.pcbi.1007536] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2019] [Revised: 12/19/2019] [Accepted: 11/11/2019] [Indexed: 11/18/2022] Open
Abstract
Kinetic models of metabolism can be constructed to predict cellular regulation and devise metabolic engineering strategies, and various promising computational workflows have been developed in recent years for this. Due to the uncertainty in the kinetic parameter values required to build kinetic models, these workflows rely on ensemble modeling (EM) principles for sampling and building populations of models describing observed physiologies. Sensitivity coefficients from metabolic control analysis (MCA) of kinetic models can provide important insight about cellular control around a given physiological steady state. However, despite considering populations of kinetic models and their model outputs, current approaches do not provide adequate tools for statistical inference. To derive conclusions from model outputs, such as MCA sensitivity coefficients, it is necessary to rank/compare populations of variables with each other. Currently existing workflows consider confidence intervals (CIs) that are derived independently for each comparable variable. Hence, it is important to derive simultaneous CIs for the variables that we wish to rank/compare. Herein, we used an existing large-scale kinetic model of Escherichia Coli metabolism to present how univariate CIs can lead to incorrect conclusions, and we present a new workflow that applies three different multivariate statistical approaches. We use the Bonferroni and the exact normal methods to build symmetric CIs using the normality assumptions. We then suggest how bootstrapping can compute asymmetric CIs whilst relaxing this normality assumption. We conclude that the Bonferroni and the exact normal methods can provide simple and efficient ways for constructing reliable CIs, with the exact normal method favored over the Bonferroni when the compared variables present dependencies. Bootstrapping, despite its significantly higher computational cost, is recommended when comparing non-normal distributions of variables. Additionally, we show how the Bonferroni method can readily be used to estimate required sample numbers to attain a certain CI size.
Collapse
Affiliation(s)
- Tuure Hameri
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), Lausanne, Switzerland
| | - Marc-Olivier Boldi
- Department of Operations, Faculty of Business and Economics, Anthropole, University of Lausanne, Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), Lausanne, Switzerland
| |
Collapse
|
90
|
Dinh HV, Suthers PF, Chan SHJ, Shen Y, Xiao T, Deewan A, Jagtap SS, Zhao H, Rao CV, Rabinowitz JD, Maranas CD. A comprehensive genome-scale model for Rhodosporidium toruloides IFO0880 accounting for functional genomics and phenotypic data. Metab Eng Commun 2019; 9:e00101. [PMID: 31720216 PMCID: PMC6838544 DOI: 10.1016/j.mec.2019.e00101] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2019] [Revised: 08/19/2019] [Accepted: 08/20/2019] [Indexed: 12/21/2022] Open
Abstract
Rhodosporidium toruloides is a red, basidiomycetes yeast that can accumulate a large amount of lipids and produce carotenoids. To better assess this non-model yeast's metabolic capabilities, we reconstructed a genome-scale model of R. toruloides IFO0880's metabolic network (iRhto1108) accounting for 2204 reactions, 1985 metabolites and 1108 genes. In this work, we integrated and supplemented the current knowledge with in-house generated biomass composition and experimental measurements pertaining to the organism's metabolic capabilities. Predictions of genotype-phenotype relations were improved through manual curation of gene-protein-reaction rules for 543 reactions leading to correct recapitulations of 84.5% of gene essentiality data (sensitivity of 94.3% and specificity of 53.8%). Organism-specific macromolecular composition and ATP maintenance requirements were experimentally measured for two separate growth conditions: (i) carbon and (ii) nitrogen limitations. Overall, iRhto1108 reproduced R. toruloides's utilization capabilities for 18 alternate substrates, matched measured wild-type growth yield, and recapitulated the viability of 772 out of 819 deletion mutants. As a demonstration to the model's fidelity in guiding engineering interventions, the OptForce procedure was applied on iRhto1108 for triacylglycerol overproduction. Suggested interventions recapitulated many of the previous successful implementations of genetic modifications and put forth a few new ones.
Collapse
Affiliation(s)
- Hoang V. Dinh
- Department of Chemical Engineering, The Pennsylvania State University, University Park, 306 Chemical and Biomedical Engineering Building, PA, 16802-4400, USA
| | - Patrick F. Suthers
- Department of Chemical Engineering, The Pennsylvania State University, University Park, 306 Chemical and Biomedical Engineering Building, PA, 16802-4400, USA
| | - Siu Hung Joshua Chan
- Department of Chemical Engineering, The Pennsylvania State University, University Park, 306 Chemical and Biomedical Engineering Building, PA, 16802-4400, USA
| | - Yihui Shen
- Department of Chemistry, Princeton University, 285 Frick Laboratory, Princeton, NJ, 08544, USA
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, 08540, USA
| | - Tianxia Xiao
- Department of Chemistry, Princeton University, 285 Frick Laboratory, Princeton, NJ, 08544, USA
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, 08540, USA
| | - Anshu Deewan
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champagne, 114 Roger Adams Laboratory MC 712, Urbana, IL, 61801, USA
| | - Sujit S. Jagtap
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champagne, 114 Roger Adams Laboratory MC 712, Urbana, IL, 61801, USA
| | - Huimin Zhao
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champagne, 114 Roger Adams Laboratory MC 712, Urbana, IL, 61801, USA
- Carl R. Woese Institute for Genomic Biology, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA
| | - Christopher V. Rao
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champagne, 114 Roger Adams Laboratory MC 712, Urbana, IL, 61801, USA
- Carl R. Woese Institute for Genomic Biology, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA
| | - Joshua D. Rabinowitz
- Department of Chemistry, Princeton University, 285 Frick Laboratory, Princeton, NJ, 08544, USA
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, 08540, USA
| | - Costas D. Maranas
- Department of Chemical Engineering, The Pennsylvania State University, University Park, 306 Chemical and Biomedical Engineering Building, PA, 16802-4400, USA
| |
Collapse
|
91
|
Abstract
Abstract
Living organisms in analogy with chemical factories use simple molecules such as sugars to produce a variety of compounds which are necessary for sustaining life and some of which are also commercially valuable. The metabolisms of simple (such as bacteria) and higher organisms (such as plants) alike can be exploited to convert low value inputs into high value outputs. Unlike conventional chemical factories, microbial production chassis are not necessarily tuned for a single product overproduction. Despite the same end goal, metabolic and industrial engineers rely on different techniques for achieving productivity goals. Metabolic engineers cannot affect reaction rates by manipulating pressure and temperature, instead they have at their disposal a range of enzymes and transcriptional and translational processes to optimize accordingly. In this review, we first highlight how various analytical approaches used in metabolic engineering and synthetic biology are related to concepts developed in systems and control engineering. Specifically, how algorithmic concepts derived in operations research can help explain the structure and organization of metabolic networks. Finally, we consider the future directions and challenges faced by the field of metabolic network modeling and the possible contributions of concepts drawn from the classical fields of chemical and control engineering. The aim of the review is to offer a current perspective of metabolic engineering and all that it entails without requiring specialized knowledge of bioinformatics or systems biology.
Collapse
|
92
|
Foster CJ, Gopalakrishnan S, Antoniewicz MR, Maranas CD. From Escherichia coli mutant 13C labeling data to a core kinetic model: A kinetic model parameterization pipeline. PLoS Comput Biol 2019; 15:e1007319. [PMID: 31504032 PMCID: PMC6759195 DOI: 10.1371/journal.pcbi.1007319] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2018] [Revised: 09/24/2019] [Accepted: 08/02/2019] [Indexed: 12/02/2022] Open
Abstract
Kinetic models of metabolic networks offer the promise of quantitative phenotype prediction. The mechanistic characterization of enzyme catalyzed reactions allows for tracing the effect of perturbations in metabolite concentrations and reaction fluxes in response to genetic and environmental perturbation that are beyond the scope of stoichiometric models. In this study, we develop a two-step computational pipeline for the rapid parameterization of kinetic models of metabolic networks using a curated metabolic model and available 13C-labeling distributions under multiple genetic and environmental perturbations. The first step involves the elucidation of all intracellular fluxes in a core model of E. coli containing 74 reactions and 61 metabolites using 13C-Metabolic Flux Analysis (13C-MFA). Here, fluxes corresponding to the mid-exponential growth phase are elucidated for seven single gene deletion mutants from upper glycolysis, pentose phosphate pathway and the Entner-Doudoroff pathway. The computed flux ranges are then used to parameterize the same (i.e., k-ecoli74) core kinetic model for E. coli with 55 substrate-level regulations using the newly developed K-FIT parameterization algorithm. The K-FIT algorithm employs a combination of equation decomposition and iterative solution techniques to evaluate steady-state fluxes in response to genetic perturbations. k-ecoli74 predicted 86% of flux values for strains used during fitting within a single standard deviation of 13C-MFA estimated values. By performing both tasks using the same network, errors associated with lack of congruity between the two networks are avoided, allowing for seamless integration of data with model building. Product yield predictions and comparison with previously developed kinetic models indicate shifts in flux ranges and the presence or absence of mutant strains delivering flux towards pathways of interest from training data significantly impact predictive capabilities. Using this workflow, the impact of completeness of fluxomic datasets and the importance of specific genetic perturbations on uncertainties in kinetic parameter estimation are evaluated.
Collapse
Affiliation(s)
- Charles J. Foster
- Department of Chemical Engineering, The Pennsylvania State University, University Park, Pennsylvania, United States of America
| | - Saratram Gopalakrishnan
- Department of Chemical Engineering, The Pennsylvania State University, University Park, Pennsylvania, United States of America
| | - Maciek R. Antoniewicz
- Department of Chemical and Biomolecular Engineering, University of Delaware. Newark, Delaware, United States of America
| | - Costas D. Maranas
- Department of Chemical Engineering, The Pennsylvania State University, University Park, Pennsylvania, United States of America
| |
Collapse
|
93
|
Long CP, Antoniewicz MR. Metabolic flux responses to deletion of 20 core enzymes reveal flexibility and limits of E. coli metabolism. Metab Eng 2019; 55:249-257. [DOI: 10.1016/j.ymben.2019.08.003] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2019] [Revised: 08/03/2019] [Accepted: 08/03/2019] [Indexed: 02/08/2023]
|
94
|
Miskovic L, Béal J, Moret M, Hatzimanikatis V. Uncertainty reduction in biochemical kinetic models: Enforcing desired model properties. PLoS Comput Biol 2019; 15:e1007242. [PMID: 31430276 PMCID: PMC6716680 DOI: 10.1371/journal.pcbi.1007242] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2019] [Revised: 08/30/2019] [Accepted: 07/03/2019] [Indexed: 11/18/2022] Open
Abstract
A persistent obstacle for constructing kinetic models of metabolism is uncertainty in the kinetic properties of enzymes. Currently, available methods for building kinetic models can cope indirectly with uncertainties by integrating data from different biological levels and origins into models. In this study, we use the recently proposed computational approach iSCHRUNK (in Silico Approach to Characterization and Reduction of Uncertainty in the Kinetic Models), which combines Monte Carlo parameter sampling methods and machine learning techniques, in the context of Bayesian inference. Monte Carlo parameter sampling methods allow us to exploit synergies between different data sources and generate a population of kinetic models that are consistent with the available data and physicochemical laws. The machine learning allows us to data-mine the a priori generated kinetic parameters together with the integrated datasets and derive posterior distributions of kinetic parameters consistent with the observed physiology. In this work, we used iSCHRUNK to address a design question: can we identify which are the kinetic parameters and what are their values that give rise to a desired metabolic behavior? Such information is important for a wide variety of studies ranging from biotechnology to medicine. To illustrate the proposed methodology, we performed Metabolic Control Analysis, computed the flux control coefficients of the xylose uptake (XTR), and identified parameters that ensure a rate improvement of XTR in a glucose-xylose co-utilizing S. cerevisiae strain. Our results indicate that only three kinetic parameters need to be accurately characterized to describe the studied physiology, and ultimately to design and control the desired responses of the metabolism. This framework paves the way for a new generation of methods that will systematically integrate the wealth of available omics data and efficiently extract the information necessary for metabolic engineering and synthetic biology decisions. Kinetic models are the most promising tool for understanding the complex dynamic behavior of living cells. The primary goal of kinetic models is to capture the properties of the metabolic networks as a whole, and thus we need large-scale models for dependable in silico analyses of metabolism. However, uncertainty in kinetic parameters impedes the development of kinetic models, and uncertainty levels increase with the model size. Tools that will address the issues with parameter uncertainty and that will be able to reduce the uncertainty propagation through the system are therefore needed. In this work, we applied a method called iSCHRUNK that combines parameter sampling and machine learning techniques to characterize the uncertainties and uncover intricate relationships between the parameters of kinetic models and the responses of the metabolic network. The proposed method allowed us to identify a small number of parameters that determine the responses in the network regardless of the values of other parameters. As a consequence, in future studies of metabolism, it will be sufficient to explore a reduced kinetic space, and more comprehensive analyses of large-scale and genome-scale metabolic networks will be computationally tractable.
Collapse
Affiliation(s)
- Ljubisa Miskovic
- Laboratory of Computational Systems Biology (LCSB), EPFL, CH, Lausanne, Switzerland
| | - Jonas Béal
- Master's Program in Life Sciences and Technology, EPFL, CH, Lausanne, Switzerland
| | - Michael Moret
- Master's Program in Life Sciences and Technology, EPFL, CH, Lausanne, Switzerland
| | | |
Collapse
|
95
|
Antoniewicz MR. Synthetic methylotrophy: Strategies to assimilate methanol for growth and chemicals production. Curr Opin Biotechnol 2019; 59:165-174. [PMID: 31437746 DOI: 10.1016/j.copbio.2019.07.001] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2019] [Revised: 07/12/2019] [Accepted: 07/17/2019] [Indexed: 11/29/2022]
Abstract
Methanol is an attractive and broadly available substrate for large-scale bioproduction of fuels and chemicals. It contains more energy and electrons per carbon than carbohydrates and can be cheaply produced from natural gas. Synthetic methylotrophy refers to the development of non-native methylotrophs such as Escherichia coli and Corynebacterium glutamicum to utilize methanol as a carbon source. Here, we discuss recent advances in engineering these industrial hosts to assimilate methanol for growth and chemicals production through the introduction of the ribulose monophosphate (RuMP) cycle. In addition, we present novel strategies based on flux coupling and adaptive laboratory evolution to engineer new strains that can grow exclusively on methanol.
Collapse
Affiliation(s)
- Maciek R Antoniewicz
- Department of Chemical and Biomolecular Engineering, Metabolic Engineering and Systems Biology Laboratory, University of Delaware, Newark DE 19716, USA.
| |
Collapse
|
96
|
Greene J, Daniell J, Köpke M, Broadbelt L, Tyo KE. Kinetic ensemble model of gas fermenting Clostridium autoethanogenum for improved ethanol production. Biochem Eng J 2019. [DOI: 10.1016/j.bej.2019.04.021] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
97
|
Amann T, Schmieder V, Faustrup Kildegaard H, Borth N, Andersen MR. Genetic engineering approaches to improve posttranslational modification of biopharmaceuticals in different production platforms. Biotechnol Bioeng 2019; 116:2778-2796. [PMID: 31237682 DOI: 10.1002/bit.27101] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Revised: 03/27/2019] [Accepted: 06/18/2019] [Indexed: 12/18/2022]
Abstract
The number of approved biopharmaceuticals, where product quality attributes remain of major importance, is increasing steadily. Within the available variety of expression hosts, the production of biopharmaceuticals faces diverse limitations with respect to posttranslational modifications (PTM), while different biopharmaceuticals demand different forms and specifications of PTMs for proper functionality. With the growing toolbox of genetic engineering technologies, it is now possible to address general as well as host- or biopharmaceutical-specific product quality obstacles. In this review, we present diverse expression systems derived from mammalians, bacteria, yeast, plants, and insects as well as available genetic engineering tools. We focus on genes for knockout/knockdown and overexpression for meaningful approaches to improve biopharmaceutical PTMs and discuss their applicability as well as future trends in the field.
Collapse
Affiliation(s)
- Thomas Amann
- The Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Valerie Schmieder
- acib GmbH-Austrian Centre of Industrial Biotechnology, Graz, Austria.,Department of Biotechnology, BOKU University of Natural Resources and Life Sciences, Vienna, Austria
| | - Helene Faustrup Kildegaard
- The Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Nicole Borth
- Department of Biotechnology, BOKU University of Natural Resources and Life Sciences, Vienna, Austria
| | - Mikael Rørdam Andersen
- Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| |
Collapse
|
98
|
Human Systems Biology and Metabolic Modelling: A Review-From Disease Metabolism to Precision Medicine. BIOMED RESEARCH INTERNATIONAL 2019; 2019:8304260. [PMID: 31281846 PMCID: PMC6590590 DOI: 10.1155/2019/8304260] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/01/2018] [Revised: 02/07/2019] [Accepted: 05/20/2019] [Indexed: 01/06/2023]
Abstract
In cell and molecular biology, metabolism is the only system that can be fully simulated at genome scale. Metabolic systems biology offers powerful abstraction tools to simulate all known metabolic reactions in a cell, therefore providing a snapshot that is close to its observable phenotype. In this review, we cover the 15 years of human metabolic modelling. We show that, although the past five years have not experienced large improvements in the size of the gene and metabolite sets in human metabolic models, their accuracy is rapidly increasing. We also describe how condition-, tissue-, and patient-specific metabolic models shed light on cell-specific changes occurring in the metabolic network, therefore predicting biomarkers of disease metabolism. We finally discuss current challenges and future promising directions for this research field, including machine/deep learning and precision medicine. In the omics era, profiling patients and biological processes from a multiomic point of view is becoming more common and less expensive. Starting from multiomic data collected from patients and N-of-1 trials where individual patients constitute different case studies, methods for model-building and data integration are being used to generate patient-specific models. Coupled with state-of-the-art machine learning methods, this will allow characterizing each patient's disease phenotype and delivering precision medicine solutions, therefore leading to preventative medicine, reduced treatment, and in silico clinical trials.
Collapse
|
99
|
Garcia S, Trinh CT. Modular design: Implementing proven engineering principles in biotechnology. Biotechnol Adv 2019; 37:107403. [PMID: 31181317 DOI: 10.1016/j.biotechadv.2019.06.002] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2019] [Revised: 04/23/2019] [Accepted: 06/04/2019] [Indexed: 12/27/2022]
Abstract
Modular design is at the foundation of contemporary engineering, enabling rapid, efficient, and reproducible construction and maintenance of complex systems across applications. Remarkably, modularity has recently been discovered as a governing principle in natural biological systems from genes to proteins to complex networks within a cell and organism communities. The convergent knowledge of natural and engineered modular systems provides a key to drive modern biotechnology to address emergent challenges associated with health, food, energy, and the environment. Here, we first present the theory and application of modular design in traditional engineering fields. We then discuss the significance and impact of modular architectures on systems biology and biotechnology. Next, we focus on the very recent theoretical and experimental advances in modular cell engineering that seeks to enable rapid and systematic development of microbial catalysts capable of efficiently synthesizing a large space of useful chemicals. We conclude with an outlook towards theoretical and practical opportunities for a more systematic and effective application of modular cell engineering in biotechnology.
Collapse
Affiliation(s)
- Sergio Garcia
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, TN, United States of America; Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States of America
| | - Cong T Trinh
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, TN, United States of America; Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States of America.
| |
Collapse
|
100
|
|