1
|
Qian J, Ye C. Development and applications of genome-scale metabolic network models. ADVANCES IN APPLIED MICROBIOLOGY 2024; 126:1-26. [PMID: 38637105 DOI: 10.1016/bs.aambs.2024.02.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/20/2024]
Abstract
The genome-scale metabolic network model is an effective tool for characterizing the gene-protein-response relationship in the entire metabolic pathway of an organism. By combining various algorithms, the genome-scale metabolic network model can effectively simulate the influence of a specific environment on the physiological state of cells, optimize the culture conditions of strains, and predict the targets of genetic modification to achieve targeted modification of strains. In this review, we summarize the whole process of model building, sort out the various tools that may be involved in the model building process, and explain the role of various algorithms in model analysis. In addition, we also summarized the application of GSMM in network characteristics, cell phenotypes, metabolic engineering, etc. Finally, we discuss the current challenges facing GSMM.
Collapse
Affiliation(s)
- Jinyi Qian
- Ministry of Education Key Laboratory of NSLSCS, Nanjing Normal University, Nanjing, PR China
| | - Chao Ye
- Ministry of Education Key Laboratory of NSLSCS, Nanjing Normal University, Nanjing, PR China; School of Food Science and Pharmaceutical Engineering, Nanjing Normal University, Nanjing, PR China.
| |
Collapse
|
2
|
Kochen MA, Hellerstein JL, Sauro HM. First-order ultrasensitivity in phosphorylation cycles. Interface Focus 2024; 14:20230045. [PMID: 38344405 PMCID: PMC10853695 DOI: 10.1098/rsfs.2023.0045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Accepted: 01/24/2024] [Indexed: 05/09/2024] Open
Abstract
Cellular signal transduction takes place through a network of phosphorylation cycles. These pathways take the form of a multi-layered cascade of cycles. This work focuses on the sensitivity of single, double and n length cycles. Cycles that operate in the zero-order regime can become sensitive to changes in signal, resulting in zero-order ultrasensitivity (ZOU). Using frequency analysis, we confirm previous efforts that cascades can act as noise filters by computing the bandwidth. We show that n length cycles display what we term first-order ultrasensitivity which occurs even when the cycles are not operating in the zero-order regime. The magnitude of the sensitivity, however, has an upper bound equal to the number of cycles. It is known that ZOU can be significantly reduced in the presence of retroactivity. We show that the first-order ultrasensitivity is immune to retroactivity and that the ZOU and first-order ultrasensitivity can be blended to create systems with constant sensitivity over a wider range of signal. We show that the ZOU in a double cycle is only modestly higher compared with a single cycle. We therefore speculate that the double cycle has evolved to enable amplification even in the face of retroactivity.
Collapse
Affiliation(s)
- Michael A. Kochen
- Department of Bioengineering, University of Washington, Seattle, WA 98105, USA
| | | | - Herbert M. Sauro
- Department of Bioengineering, University of Washington, Seattle, WA 98105, USA
| |
Collapse
|
3
|
Kulyashov MA, Kolmykov SK, Khlebodarova TM, Akberdin IR. State-of the-Art Constraint-Based Modeling of Microbial Metabolism: From Basics to Context-Specific Models with a Focus on Methanotrophs. Microorganisms 2023; 11:2987. [PMID: 38138131 PMCID: PMC10745598 DOI: 10.3390/microorganisms11122987] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 12/09/2023] [Accepted: 12/13/2023] [Indexed: 12/24/2023] Open
Abstract
Methanotrophy is the ability of an organism to capture and utilize the greenhouse gas, methane, as a source of energy-rich carbon. Over the years, significant progress has been made in understanding of mechanisms for methane utilization, mostly in bacterial systems, including the key metabolic pathways, regulation and the impact of various factors (iron, copper, calcium, lanthanum, and tungsten) on cell growth and methane bioconversion. The implementation of -omics approaches provided vast amount of heterogeneous data that require the adaptation or development of computational tools for a system-wide interrogative analysis of methanotrophy. The genome-scale mathematical modeling of its metabolism has been envisioned as one of the most productive strategies for the integration of muti-scale data to better understand methane metabolism and enable its biotechnological implementation. Herein, we provide an overview of various computational strategies implemented for methanotrophic systems. We highlight functional capabilities as well as limitations of the most popular web resources for the reconstruction, modification and optimization of the genome-scale metabolic models for methane-utilizing bacteria.
Collapse
Affiliation(s)
- Mikhail A. Kulyashov
- Department of Computational Biology, Scientific Center for Information Technologies and Artificial Intelligence, Sirius University of Science and Technology, 354340 Sochi, Russia; (M.A.K.); (S.K.K.); (T.M.K.)
- Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
| | - Semyon K. Kolmykov
- Department of Computational Biology, Scientific Center for Information Technologies and Artificial Intelligence, Sirius University of Science and Technology, 354340 Sochi, Russia; (M.A.K.); (S.K.K.); (T.M.K.)
| | - Tamara M. Khlebodarova
- Department of Computational Biology, Scientific Center for Information Technologies and Artificial Intelligence, Sirius University of Science and Technology, 354340 Sochi, Russia; (M.A.K.); (S.K.K.); (T.M.K.)
- Department of Systems Biology, Institute of Cytology and Genetics SB RAS, 630090 Novosibirsk, Russia
- Kurchatov Genomics Center, Institute of Cytology and Genetics SB RAS, 630090 Novosibirsk, Russia
| | - Ilya R. Akberdin
- Department of Computational Biology, Scientific Center for Information Technologies and Artificial Intelligence, Sirius University of Science and Technology, 354340 Sochi, Russia; (M.A.K.); (S.K.K.); (T.M.K.)
- Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
| |
Collapse
|
4
|
Schneider P, Bekiaris PS, von Kamp A, Klamt S. StrainDesign: a comprehensive Python package for computational design of metabolic networks. Bioinformatics 2022; 38:4981-4983. [PMID: 36111857 PMCID: PMC9620819 DOI: 10.1093/bioinformatics/btac632] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Revised: 08/17/2022] [Accepted: 09/15/2022] [Indexed: 10/05/2023] Open
Abstract
SUMMARY Various constraint-based optimization approaches have been developed for the computational analysis and design of metabolic networks. Herein, we present StrainDesign, a comprehensive Python package that builds upon the COBRApy toolbox and integrates the most popular metabolic design algorithms, including nested strain optimization methods such as OptKnock, RobustKnock and OptCouple as well as the more general minimal cut sets approach. The optimization approaches are embedded in individual modules, which can also be combined for setting up more elaborate strain design problems. Advanced features, such as the efficient integration of GPR rules and the possibility to consider gene and reaction additions or regulatory interventions, have been generalized and are available for all modules. The package uses state-of-the-art preprocessing methods, supports multiple solvers and provides a number of enhanced tools for analyzing computed intervention strategies including 2D and 3D plots of user-selected metabolic fluxes or yields. Furthermore, a user-friendly interface for the StrainDesign package has been implemented in the GUI-based metabolic modeling software CNApy. StrainDesign provides thus a unique and rich framework for computational strain design in Python, uniting many algorithmic developments in the field and allowing modular extension in the future. AVAILABILITY AND IMPLEMENTATION The StrainDesign package can be retrieved from PyPi, Anaconda and GitHub (https://github.com/klamt-lab/straindesign) and is also part of the latest CNApy package.
Collapse
Affiliation(s)
- Philipp Schneider
- Analysis and Redesign of Biological Networks, Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg 39106, Germany
| | - Pavlos Stephanos Bekiaris
- Analysis and Redesign of Biological Networks, Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg 39106, Germany
| | - Axel von Kamp
- Analysis and Redesign of Biological Networks, Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg 39106, Germany
| | - Steffen Klamt
- Analysis and Redesign of Biological Networks, Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg 39106, Germany
| |
Collapse
|
5
|
Beura S, Kundu P, Das AK, Ghosh A. Metagenome-scale community metabolic modelling for understanding the role of gut microbiota in human health. Comput Biol Med 2022; 149:105997. [DOI: 10.1016/j.compbiomed.2022.105997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Revised: 07/03/2022] [Accepted: 08/14/2022] [Indexed: 11/03/2022]
|
6
|
Stoakes E, Savva GM, Coates R, Tejera N, Poolman MG, Grant AJ, Wain J, Singh D. Substrate Utilisation and Energy Metabolism in Non-Growing Campylobacter jejuni M1cam. Microorganisms 2022; 10:microorganisms10071355. [PMID: 35889074 PMCID: PMC9318392 DOI: 10.3390/microorganisms10071355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Revised: 06/23/2022] [Accepted: 06/30/2022] [Indexed: 02/04/2023] Open
Abstract
Campylobacter jejuni, the major cause of bacterial foodborne illness, is also a fastidious organism that requires strict growth requirements in the laboratory. Our aim was to study substrate utilisation and energy metabolism in non-growing C. jejuni to investigate the ability of these bacteria to survive so effectively in the food chain. We integrated phenotypic microarrays and genome-scale metabolic modelling (GSM) to investigate the survival of C. jejuni on 95 substrates. We further investigated the underlying metabolic re-adjustment associated with varying energy demands on each substrate. We identified amino acids, organic acids and H2, as single substrates supporting survival without growth. We identified several different mechanisms, which were used alone or in combination, for ATP production: substrate-level phosphorylation via acetate kinase, the TCA cycle, and oxidative phosphorylation via the electron transport chain that utilised alternative electron donors and acceptors. The benefit of ATP production through each of these mechanisms was associated with the cost of enzyme investment, nutrient availability and/or O2 utilisation. C. jejuni can utilise a wide range of substrates as energy sources, including organic acids commonly used for marination or preservation of ingredients, which might contribute to the success of their survival in changing environments.
Collapse
Affiliation(s)
- Emily Stoakes
- Department of Veterinary Medicine, University of Cambridge, Cambridge CB3 0ES, UK; (E.S.); (R.C.); (A.J.G.)
| | - George M. Savva
- Quadram Institute Bioscience, Norwich Research Park, Norwich NR4 7UQ, UK; (G.M.S.); (N.T.)
| | - Ruby Coates
- Department of Veterinary Medicine, University of Cambridge, Cambridge CB3 0ES, UK; (E.S.); (R.C.); (A.J.G.)
| | - Noemi Tejera
- Quadram Institute Bioscience, Norwich Research Park, Norwich NR4 7UQ, UK; (G.M.S.); (N.T.)
| | - Mark G. Poolman
- Cell System Modelling Group, Oxford Brookes University, Oxford OX3 0BP, UK;
| | - Andrew J. Grant
- Department of Veterinary Medicine, University of Cambridge, Cambridge CB3 0ES, UK; (E.S.); (R.C.); (A.J.G.)
| | - John Wain
- Quadram Institute Bioscience, Norwich Research Park, Norwich NR4 7UQ, UK; (G.M.S.); (N.T.)
- Correspondence: (J.W.); (D.S.)
| | - Dipali Singh
- Quadram Institute Bioscience, Norwich Research Park, Norwich NR4 7UQ, UK; (G.M.S.); (N.T.)
- Correspondence: (J.W.); (D.S.)
| |
Collapse
|
7
|
A genome-scale metabolic model of Cupriavidus necator H16 integrated with TraDIS and transcriptomic data reveals metabolic insights for biotechnological applications. PLoS Comput Biol 2022; 18:e1010106. [PMID: 35604933 PMCID: PMC9166356 DOI: 10.1371/journal.pcbi.1010106] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2021] [Revised: 06/03/2022] [Accepted: 04/14/2022] [Indexed: 11/29/2022] Open
Abstract
Exploiting biological processes to recycle renewable carbon into high value platform chemicals provides a sustainable and greener alternative to current reliance on petrochemicals. In this regard Cupriavidus necator H16 represents a particularly promising microbial chassis due to its ability to grow on a wide range of low-cost feedstocks, including the waste gas carbon dioxide, whilst also naturally producing large quantities of polyhydroxybutyrate (PHB) during nutrient-limited conditions. Understanding the complex metabolic behaviour of this bacterium is a prerequisite for the design of successful engineering strategies for optimising product yields. We present a genome-scale metabolic model (GSM) of C. necator H16 (denoted iCN1361), which is directly constructed from the BioCyc database to improve the readability and reusability of the model. After the initial automated construction, we have performed extensive curation and both theoretical and experimental validation. By carrying out a genome-wide essentiality screening using a Transposon-directed Insertion site Sequencing (TraDIS) approach, we showed that the model could predict gene knockout phenotypes with a high level of accuracy. Importantly, we indicate how experimental and computational predictions can be used to improve model structure and, thus, model accuracy as well as to evaluate potential false positives identified in the experiments. Finally, by integrating transcriptomics data with iCN1361 we create a condition-specific model, which, importantly, better reflects PHB production in C. necator H16. Observed changes in the omics data and in-silico-estimated alterations in fluxes were then used to predict the regulatory control of key cellular processes. The results presented demonstrate that iCN1361 is a valuable tool for unravelling the system-level metabolic behaviour of C. necator H16 and can provide useful insights for designing metabolic engineering strategies. Genome-scale metabolic models (GSMs) provide a tool for unravelling the complex metabolic behaviour of bacteria and how they adapt to changing environments and genetic perturbations, and thus offer invaluable insights for biotechnology applications. For a GSM to be used efficiently for strain development purposes, however, the model must be easily readable and reusable by other researchers, whilst being able to predict metabolic behaviour with a high level of accuracy. In this work, we developed a GSM for Cupriavidus necator H16 that is linked to the BioCyc database, which provides an efficient way of application, model update, integration of experimental data and network visualisation for other researchers. Using our model, we demonstrate how integrating experimental observations, including Transposon-directed Insertion site Sequencing (TraDIS) and omics data, can be used to compensate for the lack of regulatory, kinetic and thermodynamic information in GSMs, and thus improve model accuracy. Importantly, we found that TraDIS in vivo screening and GSM analysis are complementary approaches, which can be used in combination to provide reliable gene essentiality predictions. Overall, our results offer an informed strategy for the deliberate manipulation of C. necator H16 metabolic capabilities, towards its industrial application to convert greenhouse gases into biochemicals and biofuels.
Collapse
|
8
|
Díaz Calvo T, Tejera N, McNamara I, Langridge GC, Wain J, Poolman M, Singh D. Genome-Scale Metabolic Modelling Approach to Understand the Metabolism of the Opportunistic Human Pathogen Staphylococcus epidermidis RP62A. Metabolites 2022; 12:metabo12020136. [PMID: 35208211 PMCID: PMC8874387 DOI: 10.3390/metabo12020136] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Revised: 01/18/2022] [Accepted: 01/29/2022] [Indexed: 02/01/2023] Open
Abstract
Staphylococcus epidermidis is a common commensal of collagen-rich regions of the body, such as the skin, but also represents a threat to patients with medical implants (joints and heart), and to preterm babies. Far less studied than Staphylococcus aureus, the mechanisms behind this increasingly recognised pathogenicity are yet to be fully understood. Improving our knowledge of the metabolic processes that allow S. epidermidis to colonise different body sites is key to defining its pathogenic potential. Thus, we have constructed a fully curated, genome-scale metabolic model for S. epidermidis RP62A, and investigated its metabolic properties with a focus on substrate auxotrophies and its utilisation for energy and biomass production. Our results show that, although glucose is available in the medium, only a small portion of it enters the glycolytic pathways, whils most is utilised for the production of biofilm, storage and the structural components of biomass. Amino acids, proline, valine, alanine, glutamate and arginine, are preferred sources of energy and biomass production. In contrast to previous studies, we have shown that this strain has no real substrate auxotrophies, although removal of proline from the media has the highest impact on the model and the experimental growth characteristics. Further study is needed to determine the significance of proline, an abundant amino acid in collagen, in S. epidermidis colonisation.
Collapse
Affiliation(s)
- Teresa Díaz Calvo
- Quadram Institute Bioscience, Norwich Research Park, Norwich NR4 7UQ, UK;
| | - Noemi Tejera
- Microbes in the Food Chain, Quadram Institute Bioscience, Norwich Research Park, Norwich NR4 7UQ, UK; (N.T.); (G.C.L.); (J.W.)
| | - Iain McNamara
- Norwich Medical School, University of East Anglia, Norwich NR4 7UQ, UK;
- Department of Orthopaedics and Trauma, Norfolk and Norwich University Hospital NHS Foundation Trust, Norwich NR4 7UY, UK
| | - Gemma C. Langridge
- Microbes in the Food Chain, Quadram Institute Bioscience, Norwich Research Park, Norwich NR4 7UQ, UK; (N.T.); (G.C.L.); (J.W.)
| | - John Wain
- Microbes in the Food Chain, Quadram Institute Bioscience, Norwich Research Park, Norwich NR4 7UQ, UK; (N.T.); (G.C.L.); (J.W.)
- Norwich Medical School, University of East Anglia, Norwich NR4 7UQ, UK;
| | - Mark Poolman
- Cell System Modelling Group, Oxford Brookes University, Oxford OX3 OBP, UK;
| | - Dipali Singh
- Microbes in the Food Chain, Quadram Institute Bioscience, Norwich Research Park, Norwich NR4 7UQ, UK; (N.T.); (G.C.L.); (J.W.)
- Correspondence:
| |
Collapse
|
9
|
Lubbock AL, Lopez CF. Programmatic modeling for biological systems. CURRENT OPINION IN SYSTEMS BIOLOGY 2021; 27:100343. [PMID: 34485764 PMCID: PMC8411905 DOI: 10.1016/j.coisb.2021.05.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/30/2023]
Abstract
Computational modeling has become an established technique to encode mathematical representations of cellular processes and gain mechanistic insights that drive testable predictions. These models are often constructed using graphical user interfaces or domain-specific languages, with community standards used for interchange. Models undergo steady state or dynamic analysis, which can include simulation and calibration within a single application, or transfer across various tools. Here, we describe a novel programmatic modeling paradigm, whereby modeling is augmented with software engineering best practices. We focus on Python - a popular programming language with a large scientific package ecosystem. Models can be encoded as programs, adding benefits such as modularity, testing, and automated documentation generators, while still being extensible and exportable to standardized formats for use with external tools if desired. Programmatic modeling is a key technology to enable collaborative model development and enhance dissemination, transparency, and reproducibility.
Collapse
Affiliation(s)
- Alexander L.R. Lubbock
- Department of Biochemistry, Vanderbilt University, Nashville, Tennessee 37212, United States of America
- Vanderbilt-Ingram Cancer Center, Vanderbilt University, Nashville Tennessee 37212, United States of America
| | - Carlos F. Lopez
- Department of Biochemistry, Vanderbilt University, Nashville, Tennessee 37212, United States of America
- Vanderbilt-Ingram Cancer Center, Vanderbilt University, Nashville Tennessee 37212, United States of America
- Department of Biomedical Informatics, Vanderbilt University, Nashville, Tennessee 37212, United States of America
| |
Collapse
|
10
|
Fell DA. Metabolic Control Analysis. Metab Eng 2021. [DOI: 10.1002/9783527823468.ch6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
11
|
van Aalst M, Ebenhöh O, Matuszyńska A. Constructing and analysing dynamic models with modelbase v1.2.3: a software update. BMC Bioinformatics 2021; 22:203. [PMID: 33879053 PMCID: PMC8056244 DOI: 10.1186/s12859-021-04122-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Accepted: 04/07/2021] [Indexed: 11/12/2022] Open
Abstract
BACKGROUND Computational mathematical models of biological and biomedical systems have been successfully applied to advance our understanding of various regulatory processes, metabolic fluxes, effects of drug therapies, and disease evolution and transmission. Unfortunately, despite community efforts leading to the development of SBML and the BioModels database, many published models have not been fully exploited, largely due to a lack of proper documentation or the dependence on proprietary software. To facilitate the reuse and further development of systems biology and systems medicine models, an open-source toolbox that makes the overall process of model construction more consistent, understandable, transparent, and reproducible is desired. RESULTS AND DISCUSSION We provide an update on the development of modelbase, a free, expandable Python package for constructing and analysing ordinary differential equation-based mathematical models of dynamic systems. It provides intuitive and unified methods to construct and solve these systems. Significantly expanded visualisation methods allow for convenient analysis of the structural and dynamic properties of models. After specifying reaction stoichiometries and rate equations modelbase can automatically assemble the associated system of differential equations. A newly provided library of common kinetic rate laws reduces the repetitiveness of the computer programming code. modelbase is also fully compatible with SBML. Previous versions provided functions for the automatic construction of networks for isotope labelling studies. Now, using user-provided label maps, modelbase v1.2.3 streamlines the expansion of classic models to their isotope-specific versions. Finally, the library of previously published models implemented in modelbase is growing continuously. Ranging from photosynthesis to tumour cell growth to viral infection evolution, all these models are now available in a transparent, reusable and unified format through modelbase. CONCLUSION With this new Python software package, which is written in currently one of the most popular programming languages, the user can develop new models and actively profit from the work of others. modelbase enables reproducing and replicating models in a consistent, tractable and expandable manner. Moreover, the expansion of models to their isotopic label-specific versions enables simulating label propagation, thus providing quantitative information regarding network topology and metabolic fluxes.
Collapse
Affiliation(s)
- Marvin van Aalst
- Institute of Quantitative and Theoretical Biology, Heinrich Heine University, Universitätsstr. 1, 40225 Düsseldorf, Germany
| | - Oliver Ebenhöh
- Institute of Quantitative and Theoretical Biology, Heinrich Heine University, Universitätsstr. 1, 40225 Düsseldorf, Germany
- CEPLAS - Cluster of Excellence on Plant Sciences, Universitätsstr. 1, 40225 Düsseldorf, Germany
| | - Anna Matuszyńska
- Institute of Quantitative and Theoretical Biology, Heinrich Heine University, Universitätsstr. 1, 40225 Düsseldorf, Germany
- CEPLAS - Cluster of Excellence on Plant Sciences, Universitätsstr. 1, 40225 Düsseldorf, Germany
| |
Collapse
|
12
|
Villanova V, Singh D, Pagliardini J, Fell D, Le Monnier A, Finazzi G, Poolman M. Boosting Biomass Quantity and Quality by Improved Mixotrophic Culture of the Diatom Phaeodactylum tricornutum. FRONTIERS IN PLANT SCIENCE 2021; 12:642199. [PMID: 33897733 PMCID: PMC8063856 DOI: 10.3389/fpls.2021.642199] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Accepted: 02/22/2021] [Indexed: 06/12/2023]
Abstract
Diatoms are photoautotrophic unicellular algae and are among the most abundant, adaptable, and diverse marine phytoplankton. They are extremely interesting not only for their ecological role but also as potential feedstocks for sustainable biofuels and high-value commodities such as omega fatty acids, because of their capacity to accumulate lipids. However, the cultivation of microalgae on an industrial scale requires higher cell densities and lipid accumulation than those found in nature to make the process economically viable. One of the known ways to induce lipid accumulation in Phaeodactylum tricornutum is nitrogen deprivation, which comes at the expense of growth inhibition and lower cell density. Thus, alternative ways need to be explored to enhance the lipid production as well as biomass density to make them sustainable at industrial scale. In this study, we have used experimental and metabolic modeling approaches to optimize the media composition, in terms of elemental composition, organic and inorganic carbon sources, and light intensity, that boost both biomass quality and quantity of P. tricornutum. Eventually, the optimized conditions were scaled-up to 2 L photobioreactors, where a better system control (temperature, pH, light, aeration/mixing) allowed a further improvement of the biomass capacity of P. tricornutum to 12 g/L.
Collapse
Affiliation(s)
- Valeria Villanova
- Department of Biological and Environmental Sciences, University of Gothenburg, Gothenburg, Sweden
- Laboratoire de Physiologie Cellulaire et Végétale, Université Grenoble Alpes (UGA), Centre National de la Recherche Scientifique (CNRS), Commissariat á l'Énergie Atomique et aux Énergies Alternatives (CEA), Institut National de Recherche pour l'Agriculture, l'Alimentation et l'Environnement, Interdisciplinary Research Institute of Grenoble, CEA Grenoble, Grenoble, France
- Fermentalg SA, Libourne, France
| | - Dipali Singh
- Microbes in the Food Chain, Quadram Institute Biosciences, Norwich Research Park, Norwich, United Kingdom
- Cell System Modelling Group, Oxford Brookes University, Oxford, United Kingdom
| | | | - David Fell
- Cell System Modelling Group, Oxford Brookes University, Oxford, United Kingdom
| | | | - Giovanni Finazzi
- Laboratoire de Physiologie Cellulaire et Végétale, Université Grenoble Alpes (UGA), Centre National de la Recherche Scientifique (CNRS), Commissariat á l'Énergie Atomique et aux Énergies Alternatives (CEA), Institut National de Recherche pour l'Agriculture, l'Alimentation et l'Environnement, Interdisciplinary Research Institute of Grenoble, CEA Grenoble, Grenoble, France
| | - Mark Poolman
- Cell System Modelling Group, Oxford Brookes University, Oxford, United Kingdom
| |
Collapse
|
13
|
Ahmad A, Tiwari A, Srivastava S. A Genome-Scale Metabolic Model of Thalassiosira pseudonana CCMP 1335 for a Systems-Level Understanding of Its Metabolism and Biotechnological Potential. Microorganisms 2020; 8:microorganisms8091396. [PMID: 32932853 PMCID: PMC7563145 DOI: 10.3390/microorganisms8091396] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Revised: 07/31/2020] [Accepted: 08/07/2020] [Indexed: 01/27/2023] Open
Abstract
Thalassiosira pseudonana is a transformable and biotechnologically promising model diatom with an ability to synthesise nutraceuticals such as fucoxanthin and store a significant amount of polyglucans and lipids including omega-3 fatty acids. While it was the first diatom to be sequenced, a systems-level analysis of its metabolism has not been done yet. This work presents first comprehensive, compartmentalized, and functional genome-scale metabolic model of the marine diatom Thalassiosira pseudonana CCMP 1335, which we have termed iThaps987. The model includes 987 genes, 2477 reactions, and 2456 metabolites. Comparison with the model of another diatom Phaeodactylum tricornutum revealed presence of 183 unique enzymes (belonging primarily to amino acid, carbohydrate, and lipid metabolism) in iThaps987. Model simulations showed a typical C3-type photosynthetic carbon fixation and suggested a preference of violaxanthin-diadinoxanthin pathway over violaxanthin-neoxanthin pathway for the production of fucoxanthin. Linear electron flow was found be active and cyclic electron flow was inactive under normal phototrophic conditions (unlike green algae and plants), validating the model predictions with previous reports. Investigation of the model for the potential of Thalassiosira pseudonana CCMP 1335 to produce other industrially useful compounds suggest iso-butanol as a foreign compound that can be synthesized by a single-gene addition. This work provides novel insights about the metabolism and potential of the organism and will be helpful to further investigate its metabolism and devise metabolic engineering strategies for the production of various compounds.
Collapse
Affiliation(s)
- Ahmad Ahmad
- Systems Biology for Biofuel Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi 110067, India;
- Department of Biotechnology, Noida International University (NIU), Noida 203201, India
| | - Archana Tiwari
- Department of Biotechnology, Noida International University (NIU), Noida 203201, India
- Correspondence: (A.T.); (S.S.); Tel.: +91-958-264-9114 (A.T.); +91-11-2674-1361 (S.S.)
| | - Shireesh Srivastava
- Systems Biology for Biofuel Group, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi 110067, India;
- Correspondence: (A.T.); (S.S.); Tel.: +91-958-264-9114 (A.T.); +91-11-2674-1361 (S.S.)
| |
Collapse
|
14
|
Bommareddy RR, Wang Y, Pearcy N, Hayes M, Lester E, Minton NP, Conradie AV. A Sustainable Chemicals Manufacturing Paradigm Using CO 2 and Renewable H 2. iScience 2020; 23:101218. [PMID: 32559729 PMCID: PMC7303982 DOI: 10.1016/j.isci.2020.101218] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2020] [Revised: 05/07/2020] [Accepted: 05/28/2020] [Indexed: 12/01/2022] Open
Abstract
The chemical industry must decarbonize to align with UN Sustainable Development Goals. A shift toward circular economies makes CO2 an attractive feedstock for producing chemicals, provided renewable H2 is available through technologies such as supercritical water (scH2O) gasification. Furthermore, high carbon and energy efficiency is paramount to favorable techno-economics, which poses a challenge to chemo-catalysis. This study demonstrates continuous gas fermentation of CO2 and H2 by the cell factory, Cupriavidus necator, to (R,R)-2,3-butanediol and isopropanol as case studies. Although a high carbon efficiency of 0.75 [(C-mol product)/(C-mol CO2)] is exemplified, the poor energy efficiency of biological CO2 fixation requires ∼8 [(mol H2)/(mol CO2)], which is techno-economically infeasible for producing commodity chemicals. Heat integration between exothermic gas fermentation and endothermic scH2O gasification overcomes this energy inefficiency. This study unlocks the promise of sustainable manufacturing using renewable feedstocks by combining the carbon efficiency of bio-catalysis with energy efficiency enforced through process engineering.
Collapse
Affiliation(s)
- Rajesh Reddy Bommareddy
- BBSRC/EPSRC Synthetic Biology Research Centre, Biodiscovery Institute (BDI), School of Life Sciences, University of Nottingham, Nottingham NG7 2RD, UK.
| | - Yanming Wang
- BBSRC/EPSRC Synthetic Biology Research Centre, Biodiscovery Institute (BDI), School of Life Sciences, University of Nottingham, Nottingham NG7 2RD, UK
| | - Nicole Pearcy
- BBSRC/EPSRC Synthetic Biology Research Centre, Biodiscovery Institute (BDI), School of Life Sciences, University of Nottingham, Nottingham NG7 2RD, UK
| | - Martin Hayes
- Johnson Matthey Technology Centre, 28 Cambridge Science Park, Milton Road, Cambridge CB4 0 FP, UK
| | - Edward Lester
- Department of Chemical & Environmental Engineering, University of Nottingham, Nottingham NG7 2RD, UK
| | - Nigel P Minton
- BBSRC/EPSRC Synthetic Biology Research Centre, Biodiscovery Institute (BDI), School of Life Sciences, University of Nottingham, Nottingham NG7 2RD, UK
| | - Alex V Conradie
- Department of Chemical & Environmental Engineering, University of Nottingham, Nottingham NG7 2RD, UK.
| |
Collapse
|
15
|
Tejera N, Crossman L, Pearson B, Stoakes E, Nasher F, Djeghout B, Poolman M, Wain J, Singh D. Genome-Scale Metabolic Model Driven Design of a Defined Medium for Campylobacter jejuni M1cam. Front Microbiol 2020; 11:1072. [PMID: 32636809 PMCID: PMC7318876 DOI: 10.3389/fmicb.2020.01072] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2020] [Accepted: 04/29/2020] [Indexed: 12/17/2022] Open
Abstract
Campylobacter jejuni, the most frequent cause of food-borne bacterial gastroenteritis, is a fastidious organism when grown in the laboratory. Oxygen is required for growth, despite the presence of the metabolic mechanism for anaerobic respiration. Amino acid auxotrophies are variably reported and energy metabolism can occur through several electron donor/acceptor combinations. Overall, the picture is one of a flexible, but vulnerable metabolism. To understand Campylobacter metabolism, we have constructed a fully curated, metabolic model for the reference organism M1 (our variant is M1cam) and validated it through laboratory experiments. Our results show that M1cam is auxotrophic for methionine, niacinamide, and pantothenate. There are complete biosynthesis pathways for all amino acids except methionine and it can produce energy, but not biomass, in the absence of oxygen. M1cam will grow in DMEM/F-12 defined media but not in the previously published Campylobacter specific defined media tested. Using the model, we identified potential auxotrophies and substrates that may improve growth. With this information, we designed simple defined media containing inorganic salts, the auxotrophic substrates, L-methionine, niacinamide, and pantothenate, pyruvate and additional amino acids L-cysteine, L-serine, and L-glutamine for growth enhancement. Our defined media supports a 1.75-fold higher growth rate than Brucella broth after 48 h at 37°C and sustains the growth of other Campylobacter jejuni strains. This media can be used to design reproducible assays that can help in better understanding the adaptation, stress resistance, and the virulence mechanisms of this pathogen. We have shown that with a well-curated metabolic model it is possible to design a media to grow this fastidious organism. This has implications for the investigation of new Campylobacter species defined through metagenomics, such as C. infans.
Collapse
Affiliation(s)
- Noemi Tejera
- Microbes in Food Chain, Quadram Institute Biosciences, Norwich Research Park, Norwich, United Kingdom
| | - Lisa Crossman
- Microbes in Food Chain, Quadram Institute Biosciences, Norwich Research Park, Norwich, United Kingdom.,SequenceAnalysis.co.uk, NRP Innovation Centre, Norwich, United Kingdom.,University of East Anglia, Norwich, United Kingdom
| | - Bruce Pearson
- Microbes in Food Chain, Quadram Institute Biosciences, Norwich Research Park, Norwich, United Kingdom
| | - Emily Stoakes
- Department of Veterinary Medicine, University of Cambridge, Cambridge, United Kingdom
| | - Fauzy Nasher
- London School of Hygiene and Tropical Medicine, University of London, London, United Kingdom
| | - Bilal Djeghout
- Microbes in Food Chain, Quadram Institute Biosciences, Norwich Research Park, Norwich, United Kingdom
| | - Mark Poolman
- Cell Systems Modelling Group, Oxford Brookes University, Oxford, United Kingdom
| | - John Wain
- Microbes in Food Chain, Quadram Institute Biosciences, Norwich Research Park, Norwich, United Kingdom
| | - Dipali Singh
- Microbes in Food Chain, Quadram Institute Biosciences, Norwich Research Park, Norwich, United Kingdom
| |
Collapse
|
16
|
Ahmad A, Pathania R, Srivastava S. Biochemical Characteristics and a Genome-Scale Metabolic Model of an Indian Euryhaline Cyanobacterium with High Polyglucan Content. Metabolites 2020; 10:metabo10050177. [PMID: 32365713 PMCID: PMC7281201 DOI: 10.3390/metabo10050177] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2020] [Revised: 01/28/2020] [Accepted: 02/05/2020] [Indexed: 12/16/2022] Open
Abstract
Marine cyanobacteria are promising microbes to capture and convert atmospheric CO2 and light into biomass and valuable industrial bio-products. Yet, reports on metabolic characteristics of non-model cyanobacteria are scarce. In this report, we show that an Indian euryhaline Synechococcus sp. BDU 130192 has biomass accumulation comparable to a model marine cyanobacterium and contains approximately double the amount of total carbohydrates, but significantly lower protein levels compared to Synechococcus sp. PCC 7002 cells. Based on its annotated chromosomal genome sequence, we present a genome scale metabolic model (GSMM) of this cyanobacterium, which we have named as iSyn706. The model includes 706 genes, 908 reactions, and 900 metabolites. The difference in the flux balance analysis (FBA) predicted flux distributions between Synechococcus sp. PCC 7002 and Synechococcus sp. BDU130192 strains mimicked the differences in their biomass compositions. Model-predicted oxygen evolution rate for Synechococcus sp. BDU130192 was found to be close to the experimentally-measured value. The model was analyzed to determine the potential of the strain for the production of various industrially-useful products without affecting growth significantly. This model will be helpful to researchers interested in understanding the metabolism as well as to design metabolic engineering strategies for the production of industrially-relevant compounds.
Collapse
Affiliation(s)
- Ahmad Ahmad
- DBT-ICGEB Center for Advanced Bioenergy Research, International Centre for Genetic Engineering and Biotechnology, New Delhi 110067, India;
- Department of Biotechnology, Noida International University, Noida, U.P. 203201, India
| | - Ruchi Pathania
- Systems Biology for Biofuels Group, International Centre for Genetic Engineering and Biotechnology, New Delhi 110067, India;
| | - Shireesh Srivastava
- DBT-ICGEB Center for Advanced Bioenergy Research, International Centre for Genetic Engineering and Biotechnology, New Delhi 110067, India;
- Systems Biology for Biofuels Group, International Centre for Genetic Engineering and Biotechnology, New Delhi 110067, India;
- Correspondence: ; Tel.: +91-11-26741361 (ext. 450)
| |
Collapse
|
17
|
Gilbert J, Pearcy N, Norman R, Millat T, Winzer K, King J, Hodgman C, Minton N, Twycross J. Gsmodutils: a python based framework for test-driven genome scale metabolic model development. Bioinformatics 2019; 35:3397-3403. [PMID: 30759197 PMCID: PMC6748746 DOI: 10.1093/bioinformatics/btz088] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2018] [Revised: 01/29/2019] [Accepted: 02/12/2019] [Indexed: 12/13/2022] Open
Abstract
MOTIVATION Genome scale metabolic models (GSMMs) are increasingly important for systems biology and metabolic engineering research as they are capable of simulating complex steady-state behaviour. Constraints based models of this form can include thousands of reactions and metabolites, with many crucial pathways that only become activated in specific simulation settings. However, despite their widespread use, power and the availability of tools to aid with the construction and analysis of large scale models, little methodology is suggested for their continued management. For example, when genome annotations are updated or new understanding regarding behaviour is discovered, models often need to be altered to reflect this. This is quickly becoming an issue for industrial systems and synthetic biotechnology applications, which require good quality reusable models integral to the design, build, test and learn cycle. RESULTS As part of an ongoing effort to improve genome scale metabolic analysis, we have developed a test-driven development methodology for the continuous integration of validation data from different sources. Contributing to the open source technology based around COBRApy, we have developed the gsmodutils modelling framework placing an emphasis on test-driven design of models through defined test cases. Crucially, different conditions are configurable allowing users to examine how different designs or curation impact a wide range of system behaviours, minimizing error between model versions. AVAILABILITY AND IMPLEMENTATION The software framework described within this paper is open source and freely available from http://github.com/SBRCNottingham/gsmodutils. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- James Gilbert
- Synthetic Biology Research Centre, University of Nottingham, Nottingham, UK
| | - Nicole Pearcy
- Synthetic Biology Research Centre, University of Nottingham, Nottingham, UK
| | - Rupert Norman
- Synthetic Biology Research Centre, University of Nottingham, Nottingham, UK
- School of Biosciences, University of Nottingham, Sutton Bonington, Loughborough, UK
| | - Thomas Millat
- Synthetic Biology Research Centre, University of Nottingham, Nottingham, UK
| | - Klaus Winzer
- Synthetic Biology Research Centre, University of Nottingham, Nottingham, UK
| | - John King
- School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | - Charlie Hodgman
- Synthetic Biology Research Centre, University of Nottingham, Nottingham, UK
- School of Biosciences, University of Nottingham, Sutton Bonington, Loughborough, UK
| | - Nigel Minton
- Synthetic Biology Research Centre, University of Nottingham, Nottingham, UK
| | - Jamie Twycross
- School of Computer Science, University of Nottingham, Nottingham, UK
| |
Collapse
|
18
|
Workflow for Data Analysis in Experimental and Computational Systems Biology: Using Python as ‘Glue’. Processes (Basel) 2019. [DOI: 10.3390/pr7070460] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Bottom-up systems biology entails the construction of kinetic models of cellular pathways by collecting kinetic information on the pathway components (e.g., enzymes) and collating this into a kinetic model, based for example on ordinary differential equations. This requires integration and data transfer between a variety of tools, ranging from data acquisition in kinetics experiments, to fitting and parameter estimation, to model construction, evaluation and validation. Here, we present a workflow that uses the Python programming language, specifically the modules from the SciPy stack, to facilitate this task. Starting from raw kinetics data, acquired either from spectrophotometric assays with microtitre plates or from Nuclear Magnetic Resonance (NMR) spectroscopy time-courses, we demonstrate the fitting and construction of a kinetic model using scientific Python tools. The analysis takes place in a Jupyter notebook, which keeps all information related to a particular experiment together in one place and thus serves as an e-labbook, enhancing reproducibility and traceability. The Python programming language serves as an ideal foundation for this framework because it is powerful yet relatively easy to learn for the non-programmer, has a large library of scientific routines and active user community, is open-source and extensible, and many computational systems biology software tools are written in Python or have a Python Application Programming Interface (API). Our workflow thus enables investigators to focus on the scientific problem at hand rather than worrying about data integration between disparate platforms.
Collapse
|
19
|
Norman RO, Millat T, Schatschneider S, Henstra AM, Breitkopf R, Pander B, Annan FJ, Piatek P, Hartman HB, Poolman MG, Fell DA, Winzer K, Minton NP, Hodgman C. Genome‐scale model of
C. autoethanogenum
reveals optimal bioprocess conditions for high‐value chemical production from carbon monoxide. ENGINEERING BIOLOGY 2019. [DOI: 10.1049/enb.2018.5003] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Rupert O.J. Norman
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
- School of BiosciencesUniversity of NottinghamSutton Bonington Campus, Sutton BoningtonLeicestershireLE12 5RDUK
| | - Thomas Millat
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Sarah Schatschneider
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
- Evonik Nutrition and Care GmbHKantstr. 233798Halle‐KinsbeckGermany
| | - Anne M. Henstra
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Ronja Breitkopf
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Bart Pander
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Florence J. Annan
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Pawel Piatek
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Hassan B. Hartman
- Department of Biology and Medical SciencesOxford Brookes UniversityOxfordOX3 0BPUK
- Public Health England61 Colindale AvenueLondonNW9 5EQUK
| | - Mark G. Poolman
- Department of Biology and Medical SciencesOxford Brookes UniversityOxfordOX3 0BPUK
| | - David A. Fell
- Department of Biology and Medical SciencesOxford Brookes UniversityOxfordOX3 0BPUK
| | - Klaus Winzer
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Nigel P. Minton
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Charlie Hodgman
- Synthetic Biology Research CentreUniversity of Nottingham, University ParkNottinghamNG7 2RDUK
- School of BiosciencesUniversity of NottinghamSutton Bonington Campus, Sutton BoningtonLeicestershireLE12 5RDUK
| |
Collapse
|
20
|
Heirendt L, Arreckx S, Pfau T, Mendoza SN, Richelle A, Heinken A, Haraldsdóttir HS, Wachowiak J, Keating SM, Vlasov V, Magnusdóttir S, Ng CY, Preciat G, Žagare A, Chan SHJ, Aurich MK, Clancy CM, Modamio J, Sauls JT, Noronha A, Bordbar A, Cousins B, El Assal DC, Valcarcel LV, Apaolaza I, Ghaderi S, Ahookhosh M, Ben Guebila M, Kostromins A, Sompairac N, Le HM, Ma D, Sun Y, Wang L, Yurkovich JT, Oliveira MAP, Vuong PT, El Assal LP, Kuperstein I, Zinovyev A, Hinton HS, Bryant WA, Aragón Artacho FJ, Planes FJ, Stalidzans E, Maass A, Vempala S, Hucka M, Saunders MA, Maranas CD, Lewis NE, Sauter T, Palsson BØ, Thiele I, Fleming RMT. Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0. Nat Protoc 2019; 14:639-702. [PMID: 30787451 PMCID: PMC6635304 DOI: 10.1038/s41596-018-0098-2] [Citation(s) in RCA: 577] [Impact Index Per Article: 115.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Constraint-based reconstruction and analysis (COBRA) provides a molecular mechanistic framework for integrative analysis of experimental molecular systems biology data and quantitative prediction of physicochemically and biochemically feasible phenotypic states. The COBRA Toolbox is a comprehensive desktop software suite of interoperable COBRA methods. It has found widespread application in biology, biomedicine, and biotechnology because its functions can be flexibly combined to implement tailored COBRA protocols for any biochemical network. This protocol is an update to the COBRA Toolbox v.1.0 and v.2.0. Version 3.0 includes new methods for quality-controlled reconstruction, modeling, topological analysis, strain and experimental design, and network visualization, as well as network integration of chemoinformatic, metabolomic, transcriptomic, proteomic, and thermochemical data. New multi-lingual code integration also enables an expansion in COBRA application scope via high-precision, high-performance, and nonlinear numerical optimization solvers for multi-scale, multi-cellular, and reaction kinetic modeling, respectively. This protocol provides an overview of all these new features and can be adapted to generate and analyze constraint-based models in a wide variety of scenarios. The COBRA Toolbox v.3.0 provides an unparalleled depth of COBRA methods.
Collapse
Affiliation(s)
- Laurent Heirendt
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Sylvain Arreckx
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Thomas Pfau
- Life Sciences Research Unit, University of Luxembourg, Belvaux, Luxembourg
| | - Sebastián N Mendoza
- Center for Genome Regulation (Fondap 15090007), University of Chile, Santiago, Chile
- Mathomics, Center for Mathematical Modeling, University of Chile, Santiago, Chile
| | - Anne Richelle
- Department of Pediatrics, University of California, San Diego, School of Medicine, La Jolla, CA, USA
| | - Almut Heinken
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Hulda S Haraldsdóttir
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Jacek Wachowiak
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Sarah M Keating
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Hinxton, Cambridge, UK
| | - Vanja Vlasov
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Stefania Magnusdóttir
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Chiam Yu Ng
- Department of Chemical Engineering, The Pennsylvania State University, State College, PA, USA
| | - German Preciat
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Alise Žagare
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Siu H J Chan
- Department of Chemical Engineering, The Pennsylvania State University, State College, PA, USA
| | - Maike K Aurich
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Catherine M Clancy
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Jennifer Modamio
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - John T Sauls
- Department of Physics, and Bioinformatics and Systems Biology Program, University of California, San Diego, La Jolla, CA, USA
| | - Alberto Noronha
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | | | - Benjamin Cousins
- Algorithms and Randomness Center, School of Computer Science, Georgia Institute of Technology, Atlanta, GA, USA
| | - Diana C El Assal
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Luis V Valcarcel
- Biomedical Engineering and Sciences Department, TECNUN, University of Navarra, San Sebastián, Spain
| | - Iñigo Apaolaza
- Biomedical Engineering and Sciences Department, TECNUN, University of Navarra, San Sebastián, Spain
| | - Susan Ghaderi
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Masoud Ahookhosh
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Marouen Ben Guebila
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Andrejs Kostromins
- Institute of Microbiology and Biotechnology, University of Latvia, Riga, Latvia
| | - Nicolas Sompairac
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, Paris, France
| | - Hoai M Le
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Ding Ma
- Department of Management Science and Engineering, Stanford University, Stanford, CA, USA
| | - Yuekai Sun
- Department of Statistics, University of Michigan, Ann Arbor, MI, USA
| | - Lin Wang
- Department of Chemical Engineering, The Pennsylvania State University, State College, PA, USA
| | - James T Yurkovich
- Department of Bioengineering, University of California, San Diego, La Jolla, CA, USA
| | - Miguel A P Oliveira
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Phan T Vuong
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Lemmer P El Assal
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Inna Kuperstein
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, Paris, France
| | - Andrei Zinovyev
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, Paris, France
| | - H Scott Hinton
- Utah State University Research Foundation, North Logan, UT, USA
| | - William A Bryant
- Centre for Integrative Systems Biology and Bioinformatics, Department of Life Sciences, Imperial College London, London, UK
| | | | - Francisco J Planes
- Biomedical Engineering and Sciences Department, TECNUN, University of Navarra, San Sebastián, Spain
| | - Egils Stalidzans
- Institute of Microbiology and Biotechnology, University of Latvia, Riga, Latvia
| | - Alejandro Maass
- Center for Genome Regulation (Fondap 15090007), University of Chile, Santiago, Chile
- Mathomics, Center for Mathematical Modeling, University of Chile, Santiago, Chile
| | - Santosh Vempala
- Algorithms and Randomness Center, School of Computer Science, Georgia Institute of Technology, Atlanta, GA, USA
| | - Michael Hucka
- Department of Computing and Mathematical Sciences, California Institute of Technology, Pasadena, CA, USA
| | - Michael A Saunders
- Department of Management Science and Engineering, Stanford University, Stanford, CA, USA
| | - Costas D Maranas
- Department of Chemical Engineering, The Pennsylvania State University, State College, PA, USA
| | - Nathan E Lewis
- Department of Pediatrics, University of California, San Diego, School of Medicine, La Jolla, CA, USA
- Novo Nordisk Foundation Center for Biosustainability, University of California, San Diego, La Jolla, CA, USA
| | - Thomas Sauter
- Life Sciences Research Unit, University of Luxembourg, Belvaux, Luxembourg
| | - Bernhard Ø Palsson
- Department of Bioengineering, University of California, San Diego, La Jolla, CA, USA
- Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Kemitorvet, Lyngby, Denmark
| | - Ines Thiele
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Ronan M T Fleming
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg.
- Division of Systems Biomedicine and Pharmacology, Leiden Academic Centre for Drug Research, Faculty of Science, Leiden University, Leiden, The Netherlands.
| |
Collapse
|
21
|
Huma B, Kundu S, Poolman MG, Kruger NJ, Fell DA. Stoichiometric analysis of the energetics and metabolic impact of photorespiration in C3 plants. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2018; 96:1228-1241. [PMID: 30257035 DOI: 10.1111/tpj.14105] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/10/2018] [Revised: 09/10/2018] [Accepted: 09/17/2018] [Indexed: 06/08/2023]
Abstract
Analysis of the impact of photorespiration on plant metabolism is usually based on manual inspection of small network diagrams. Here we create a structural metabolic model that contains the reactions that participate in photorespiration in the plastid, peroxisome, mitochondrion and cytosol, and the metabolite exchanges between them. This model was subjected to elementary flux modes analysis, a technique that enumerates all the component, minimal pathways of a network. Any feasible photorespiratory metabolism in the plant will be some combination of the elementary flux modes (EFMs) that contain the Rubisco oxygenase reaction. Amongst the EFMs we obtained was the classic photorespiratory cycle, but there were also modes that involve photorespiration coupled with mitochondrial metabolism and ATP production, the glutathione-ascorbate cycle and nitrate reduction to ammonia. The modes analysis demonstrated the underlying basis of the metabolic linkages with photorespiration that have been inferred experimentally. The set of reactions common to all the elementary modes showed good agreement with the gene products of mutants that have been reported to have a defective phenotype in photorespiratory conditions. Finally, the set of modes provided a formal demonstration that photorespiration itself does not impact on the CO2 :O2 ratio (assimilation quotient), except in those modes associated with concomitant nitrate reduction.
Collapse
Affiliation(s)
- Benazir Huma
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92 APC Road, Kolkata, 700 009, West Bengal, India
| | - Sudip Kundu
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92 APC Road, Kolkata, 700 009, West Bengal, India
| | - Mark G Poolman
- Department of Biological and Medical Sciences, Oxford Brookes University, Gipsy Lane, Headington, Oxford, OX3 OBP, UK
| | - Nicholas J Kruger
- Department of Plant Sciences, University of Oxford, South Parks Road, Oxford, OX1 3RB, UK
| | - David A Fell
- Department of Biological and Medical Sciences, Oxford Brookes University, Gipsy Lane, Headington, Oxford, OX3 OBP, UK
| |
Collapse
|
22
|
Methods for automated genome-scale metabolic model reconstruction. Biochem Soc Trans 2018; 46:931-936. [PMID: 30065105 DOI: 10.1042/bst20170246] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Revised: 06/04/2018] [Accepted: 06/06/2018] [Indexed: 11/17/2022]
Abstract
In the era of next-generation sequencing and ubiquitous assembly and binning of metagenomes, new putative genome sequences are being produced from isolate and microbiome samples at ever-increasing rates. Genome-scale metabolic models have enormous utility for supporting the analysis and predictive characterization of these genomes based on sequence data. As a result, tools for rapid automated reconstruction of metabolic models are becoming critically important for supporting the analysis of new genome sequences. Many tools and algorithms have now emerged to support rapid model reconstruction and analysis. Here, we are comparing and contrasting the capabilities and output of a variety of these tools, including ModelSEED, Raven Toolbox, PathwayTools, SuBliMinal Toolbox and merlin.
Collapse
|
23
|
Chatterjee A, Huma B, Shaw R, Kundu S. Reconstruction of Oryza sativa indica Genome Scale Metabolic Model and Its Responses to Varying RuBisCO Activity, Light Intensity, and Enzymatic Cost Conditions. FRONTIERS IN PLANT SCIENCE 2017; 8:2060. [PMID: 29250098 PMCID: PMC5715477 DOI: 10.3389/fpls.2017.02060] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/23/2017] [Accepted: 11/17/2017] [Indexed: 05/12/2023]
Abstract
To combat decrease in rice productivity under different stresses, an understanding of rice metabolism is needed. Though there are different genome scale metabolic models (GSMs) of Oryza sativa japonica, no GSM with gene-protein-reaction association exist for Oryza sativa indica. Here, we report a GSM, OSI1136 of O.s. indica, which includes 3602 genes and 1136 metabolic reactions and transporters distributed across the cytosol, mitochondrion, peroxisome, and chloroplast compartments. Flux balance analysis of the model showed that for varying RuBisCO activity (Vc/Vo) (i) the activity of the chloroplastic malate valve increases to transport reducing equivalents out of the chloroplast under increased photorespiratory conditions and (ii) glyceraldehyde-3-phosphate dehydrogenase and phosphoglycerate kinase can act as source of cytosolic ATP under decreased photorespiration. Under increasing light conditions we observed metabolic flexibility, involving photorespiration, chloroplastic triose phosphate and the dicarboxylate transporters of the chloroplast and mitochondrion for redox and ATP exchanges across the intracellular compartments. Simulations under different enzymatic cost conditions revealed (i) participation of peroxisomal glutathione-ascorbate cycle in photorespiratory H2O2 metabolism (ii) different modes of the chloroplastic triose phosphate transporters and malate valve, and (iii) two possible modes of chloroplastic Glu-Gln transporter which were related with the activity of chloroplastic and cytosolic isoforms of glutamine synthetase. Altogether, our results provide new insights into plant metabolism.
Collapse
Affiliation(s)
| | | | | | - Sudip Kundu
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, Kolkata, India
| |
Collapse
|
24
|
Ahmad A, Hartman HB, Krishnakumar S, Fell DA, Poolman MG, Srivastava S. A Genome Scale Model of Geobacillus thermoglucosidasius (C56-YS93) reveals its biotechnological potential on rice straw hydrolysate. J Biotechnol 2017; 251:30-37. [DOI: 10.1016/j.jbiotec.2017.03.031] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2016] [Revised: 03/27/2017] [Accepted: 03/27/2017] [Indexed: 01/29/2023]
|
25
|
Pentjuss A, Stalidzans E, Liepins J, Kokina A, Martynova J, Zikmanis P, Mozga I, Scherbaka R, Hartman H, Poolman MG, Fell DA, Vigants A. Model-based biotechnological potential analysis of Kluyveromyces marxianus central metabolism. J Ind Microbiol Biotechnol 2017; 44:1177-1190. [PMID: 28444480 DOI: 10.1007/s10295-017-1946-8] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2016] [Accepted: 04/16/2017] [Indexed: 12/11/2022]
Abstract
The non-conventional yeast Kluyveromyces marxianus is an emerging industrial producer for many biotechnological processes. Here, we show the application of a biomass-linked stoichiometric model of central metabolism that is experimentally validated, and mass and charge balanced for assessing the carbon conversion efficiency of wild type and modified K. marxianus. Pairs of substrates (lactose, glucose, inulin, xylose) and products (ethanol, acetate, lactate, glycerol, ethyl acetate, succinate, glutamate, phenylethanol and phenylalanine) are examined by various modelling and optimisation methods. Our model reveals the organism's potential for industrial application and metabolic engineering. Modelling results imply that the aeration regime can be used as a tool to optimise product yield and flux distribution in K. marxianus. Also rebalancing NADH and NADPH utilisation can be used to improve the efficiency of substrate conversion. Xylose is identified as a biotechnologically promising substrate for K. marxianus.
Collapse
Affiliation(s)
- A Pentjuss
- Institute of Microbiology and Biotechnology, University of Latvia, Jelgavas str. 1, Riga, 1004, Latvia
| | - E Stalidzans
- Institute of Microbiology and Biotechnology, University of Latvia, Jelgavas str. 1, Riga, 1004, Latvia.
| | - J Liepins
- Institute of Microbiology and Biotechnology, University of Latvia, Jelgavas str. 1, Riga, 1004, Latvia
| | - A Kokina
- Institute of Microbiology and Biotechnology, University of Latvia, Jelgavas str. 1, Riga, 1004, Latvia
| | - J Martynova
- Institute of Microbiology and Biotechnology, University of Latvia, Jelgavas str. 1, Riga, 1004, Latvia
| | - P Zikmanis
- Institute of Microbiology and Biotechnology, University of Latvia, Jelgavas str. 1, Riga, 1004, Latvia
| | - I Mozga
- Institute of Microbiology and Biotechnology, University of Latvia, Jelgavas str. 1, Riga, 1004, Latvia
| | - R Scherbaka
- Institute of Microbiology and Biotechnology, University of Latvia, Jelgavas str. 1, Riga, 1004, Latvia
| | - H Hartman
- Department of Biological and Medical Sciences, Oxford Brookes University, Headington, OX, OX3 0BP, UK
| | - M G Poolman
- Department of Biological and Medical Sciences, Oxford Brookes University, Headington, OX, OX3 0BP, UK
| | - D A Fell
- Department of Biological and Medical Sciences, Oxford Brookes University, Headington, OX, OX3 0BP, UK
| | - A Vigants
- Institute of Microbiology and Biotechnology, University of Latvia, Jelgavas str. 1, Riga, 1004, Latvia
| |
Collapse
|
26
|
Pfau T, Pacheco MP, Sauter T. Towards improved genome-scale metabolic network reconstructions: unification, transcript specificity and beyond. Brief Bioinform 2016; 17:1060-1069. [PMID: 26615025 PMCID: PMC5142010 DOI: 10.1093/bib/bbv100] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2015] [Revised: 10/20/2015] [Indexed: 12/24/2022] Open
Abstract
Genome-scale metabolic network reconstructions provide a basis for the investigation of the metabolic properties of an organism. There are reconstructions available for multiple organisms, from prokaryotes to higher organisms and methods for the analysis of a reconstruction. One example is the use of flux balance analysis to improve the yields of a target chemical, which has been applied successfully. However, comparison of results between existing reconstructions and models presents a challenge because of the heterogeneity of the available reconstructions, for example, of standards for presenting gene-protein-reaction associations, nomenclature of metabolites and reactions or selection of protonation states. The lack of comparability for gene identifiers or model-specific reactions without annotated evidence often leads to the creation of a new model from scratch, as data cannot be properly matched otherwise. In this contribution, we propose to improve the predictive power of metabolic models by switching from gene-protein-reaction associations to transcript-isoform-reaction associations, thus taking advantage of the improvement of precision in gene expression measurements. To achieve this precision, we discuss available databases that can be used to retrieve this type of information and point at issues that can arise from their neglect. Further, we stress issues that arise from non-standardized building pipelines, like inconsistencies in protonation states. In addition, problems arising from the use of non-specific cofactors, e.g. artificial futile cycles, are discussed, and finally efforts of the metabolic modelling community to unify model reconstructions are highlighted.
Collapse
|
27
|
Flux balance analysis of genome-scale metabolic model of rice (Oryza sativa): aiming to increase biomass. J Biosci 2016; 40:819-28. [PMID: 26564982 DOI: 10.1007/s12038-015-9563-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
Due to socio-economic reasons, it is essential to design efficient stress-tolerant, more nutritious, high yielding rice varieties. A systematic understanding of the rice cellular metabolism is essential for this purpose. Here, we analyse a genome-scale metabolic model of rice leaf using Flux Balance Analysis to investigate whether it has potential metabolic flexibility to increase the biosynthesis of any of the biomass components. We initially simulate the metabolic responses under an objective to maximize the biomass components. Using the estimated maximum value of biomass synthesis as a constraint, we further simulate the metabolic responses optimizing the cellular economy. Depending on the physiological conditions of a cell, the transport capacities of intracellular transporters (ICTs) can vary. To mimic this physiological state, we randomly vary the ICTs' transport capacities and investigate their effects. The results show that the rice leaf has the potential to increase glycine and starch in a wide range depending on the ICTs' transport capacities. The predicted biosynthesis pathways vary slightly at the two different optimization conditions. With the constraint of biomass composition, the cell also has the metabolic plasticity to fix a wide range of carbon-nitrogen ratio.
Collapse
|
28
|
Yuan H, Cheung CYM, Poolman MG, Hilbers PAJ, van Riel NAW. A genome-scale metabolic network reconstruction of tomato (Solanum lycopersicum L.) and its application to photorespiratory metabolism. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2016; 85:289-304. [PMID: 26576489 DOI: 10.1111/tpj.13075] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2015] [Revised: 11/01/2015] [Accepted: 11/03/2015] [Indexed: 05/09/2023]
Abstract
Tomato (Solanum lycopersicum L.) has been studied extensively due to its high economic value in the market, and high content in health-promoting antioxidant compounds. Tomato is also considered as an excellent model organism for studying the development and metabolism of fleshy fruits. However, the growth, yield and fruit quality of tomatoes can be affected by drought stress, a common abiotic stress for tomato. To investigate the potential metabolic response of tomato plants to drought, we reconstructed iHY3410, a genome-scale metabolic model of tomato leaf, and used this metabolic network to simulate tomato leaf metabolism. The resulting model includes 3410 genes and 2143 biochemical and transport reactions distributed across five intracellular organelles including cytosol, plastid, mitochondrion, peroxisome and vacuole. The model successfully described the known metabolic behaviour of tomato leaf under heterotrophic and phototrophic conditions. The in silico investigation of the metabolic characteristics for photorespiration and other relevant metabolic processes under drought stress suggested that: (i) the flux distributions through the mevalonate (MVA) pathway under drought were distinct from that under normal conditions; and (ii) the changes in fluxes through core metabolic pathways with varying flux ratio of RubisCO carboxylase to oxygenase may contribute to the adaptive stress response of plants. In addition, we improved on previous studies of reaction essentiality analysis for leaf metabolism by including potential alternative routes for compensating reaction knockouts. Altogether, the genome-scale model provides a sound framework for investigating tomato metabolism and gives valuable insights into the functional consequences of abiotic stresses.
Collapse
Affiliation(s)
- Huili Yuan
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | | | - Mark G Poolman
- Cell Systems Modelling Group, Department of Biomedical and Medical Science, Oxford Brookes University, Oxford, UK
| | - Peter A J Hilbers
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
- Institute for Complex Molecular Systems, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Natal A W van Riel
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
- Institute for Complex Molecular Systems, Eindhoven University of Technology, Eindhoven, The Netherlands
| |
Collapse
|
29
|
Chatterjee A, Kundu S. Revisiting the chlorophyll biosynthesis pathway using genome scale metabolic model of Oryza sativa japonica. Sci Rep 2015; 5:14975. [PMID: 26443104 PMCID: PMC4595741 DOI: 10.1038/srep14975] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2015] [Accepted: 08/27/2015] [Indexed: 12/30/2022] Open
Abstract
Chlorophyll is one of the most important pigments present in green plants and rice is one of the major food crops consumed worldwide. We curated the existing genome scale metabolic model (GSM) of rice leaf by incorporating new compartment, reactions and transporters. We used this modified GSM to elucidate how the chlorophyll is synthesized in a leaf through a series of bio-chemical reactions spanned over different organelles using inorganic macronutrients and light energy. We predicted the essential reactions and the associated genes of chlorophyll synthesis and validated against the existing experimental evidences. Further, ammonia is known to be the preferred source of nitrogen in rice paddy fields. The ammonia entering into the plant is assimilated in the root and leaf. The focus of the present work is centered on rice leaf metabolism. We studied the relative importance of ammonia transporters through the chloroplast and the cytosol and their interlink with other intracellular transporters. Ammonia assimilation in the leaves takes place by the enzyme glutamine synthetase (GS) which is present in the cytosol (GS1) and chloroplast (GS2). Our results provided possible explanation why GS2 mutants show normal growth under minimum photorespiration and appear chlorotic when exposed to air.
Collapse
Affiliation(s)
- Ankita Chatterjee
- 1Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta India
| | - Sudip Kundu
- 1Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta India.,Center of Excellence in Systems Biology and Biomedical Engineering, TEQIP Phase-II, University of Calcutta India
| |
Collapse
|
30
|
Adlakha N, Pfau T, Ebenhöh O, Yazdani SS. Insight into metabolic pathways of the potential biofuel producer, Paenibacillus polymyxa ICGEB2008. BIOTECHNOLOGY FOR BIOFUELS 2015; 8:159. [PMID: 26413158 PMCID: PMC4583153 DOI: 10.1186/s13068-015-0338-4] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/23/2015] [Accepted: 09/09/2015] [Indexed: 06/05/2023]
Abstract
BACKGROUND Paenibacillus polymyxa is a facultative anaerobe known for production of hydrolytic enzymes and various important biofuel molecules. Despite its wide industrial use and the availability of its genome sequence, very little is known about metabolic pathways operative in the Paenibacillus system. Here, we report metabolic insights of an insect gut symbiont, Paenibacillus polymyxa ICGEB2008, and reveal pathways playing an important role in the production of 2,3-butanediol and ethanol. RESULT We developed a metabolic network model of P. polymyxa ICGEB2008 with 133 metabolites and 158 reactions. Flux balance analysis was employed to investigate the importance of redox balance in ICGEB2008. This led to the detection of the Bifid shunt, a pathway previously not described in Paenibacillus, which can uncouple the production of ATP from the generation of reducing equivalents. Using a combined experimental and modeling approach, we further studied pathways involved in 2,3-butanediol and ethanol production and also demonstrated the production of hydrogen by the organism. We could further show that the nitrogen source is critical for metabolite production by Paenibacillus, and correctly quantify the influence on the by-product metabolite profile of ICGEB2008. Both simulations and experiments showed that metabolic flux is diverted from ethanol to acetate production when an oxidized nitrogen source is utilized. CONCLUSION We have created a predictive model of the central carbon metabolism of P. polymyxa ICGEB2008 and could show the presence of the Bifid shunt and explain its role in ICGEB2008. An in-depth study has been performed to understand the metabolic pathways involved in ethanol, 2,3-butanediol and hydrogen production, which can be utilized as a basis for further metabolic engineering efforts to improve the efficiency of biofuel production by this P. polymyxa strain.
Collapse
Affiliation(s)
- Nidhi Adlakha
- />Synthetic Biology and Biofuels Group, International Centre for Genetic Engineering and Biotechnology, Aruna Asaf Ali Marg, New Delhi, India
| | - Thomas Pfau
- />Institute of Complex Systems and Mathematical Biology, University of Aberdeen, Aberdeen, UK
- />Life Sciences Research Unit, University of Luxembourg, Luxembourg City, Luxembourg
| | - Oliver Ebenhöh
- />Institute of Complex Systems and Mathematical Biology, University of Aberdeen, Aberdeen, UK
- />Cluster of Excellence on Plant Sciences (CEPLAS), Heinrich-Heine-University, Universitätsstraße 1, 40225 Düsseldorf, Germany
| | - Syed Shams Yazdani
- />Synthetic Biology and Biofuels Group, International Centre for Genetic Engineering and Biotechnology, Aruna Asaf Ali Marg, New Delhi, India
| |
Collapse
|
31
|
Shaw R, Kundu S. Metabolic Plasticity and Inter-Compartmental Interactions in Rice Metabolism: An Analysis from Reaction Deletion Study. PLoS One 2015. [PMID: 26222686 PMCID: PMC4519304 DOI: 10.1371/journal.pone.0133899] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
More than 20% of the total caloric intake of human population comes from rice. The expression of rice genes and hence, the concentration of enzymatic proteins might vary due to several biotic and abiotic stresses. It in turn, can influence the overall metabolism and survivability of rice plant. Thus, understanding the rice cellular metabolism, its plasticity and potential readjustments under different perturbations can help rice biotechnologists to design efficient rice cultivars. Here, using the flux balance analysis (FBA) method, with the help of in-silico reaction deletion strategy, we study the metabolic plasticity of genome-scale metabolic model of rice leaf. A set of 131 reactions, essential for the production of primary biomass precursors is identified; deletion of any of them can inhibit the overall biomass production. Usability Index (IU) for the rest of the reactions are estimated and based on this parameter, they are classified into three categories—maximally-favourable, quasi-favourable and unfavourable for the primary biomass production. The lower value of 1 − IU of a reaction suggests that the cell cannot easily bypass it for biomass production. While some of the alternative paths are energetically equally efficient, others demand for higher photon. The variations in (i) ATP/NADPH ratio, (ii) exchange of metabolites through chloroplastic transporters and (iii) total biomass production are also presented here. Mutual metabolic dependencies of different cellular compartments are also demonstrated.
Collapse
Affiliation(s)
- Rahul Shaw
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92 APC Road, Kolkata 700009, West Bengal, India
| | - Sudip Kundu
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92 APC Road, Kolkata 700009, West Bengal, India
- Centre of Excellence in Systems Biology & Biomedical Engineering, (TEQIP, Phase-II), University of Calcutta, 92 APC Road, Kolkata 700009, West Bengal, India
- * E-mail:
| |
Collapse
|
32
|
Stable isotope-labeling studies in metabolomics: new insights into structure and dynamics of metabolic networks. Bioanalysis 2014; 6:511-24. [PMID: 24568354 DOI: 10.4155/bio.13.348] [Citation(s) in RCA: 135] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
The rapid emergence of metabolomics has enabled system-wide measurements of metabolites in various organisms. However, advances in the mechanistic understanding of metabolic networks remain limited, as most metabolomics studies cannot routinely provide accurate metabolite identification, absolute quantification and flux measurement. Stable isotope labeling offers opportunities to overcome these limitations. Here we describe some current approaches to stable isotope-labeled metabolomics and provide examples of the significant impact that these studies have had on our understanding of cellular metabolism. Furthermore, we discuss recently developed software solutions for the analysis of stable isotope-labeled metabolomics data and propose the bioinformatics solutions that will pave the way for the broader application and optimal interpretation of system-scale labeling studies in metabolomics.
Collapse
|
33
|
Jelsbak L, Hartman H, Schroll C, Rosenkrantz JT, Lemire S, Wallrodt I, Thomsen LE, Poolman M, Kilstrup M, Jensen PR, Olsen JE. Identification of metabolic pathways essential for fitness of Salmonella Typhimurium in vivo. PLoS One 2014; 9:e101869. [PMID: 24992475 PMCID: PMC4081726 DOI: 10.1371/journal.pone.0101869] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2014] [Accepted: 06/12/2014] [Indexed: 01/27/2023] Open
Abstract
Bacterial infections remain a threat to human and animal health worldwide, and there is an urgent need to find novel targets for intervention. In the current study we used a computer model of the metabolic network of Salmonella enterica serovar Typhimurium and identified pairs of reactions (cut sets) predicted to be required for growth in vivo. We termed such cut sets synthetic auxotrophic pairs. We tested whether these would reveal possible combined targets for new antibiotics by analyzing the performance of selected single and double mutants in systemic mouse infections. One hundred and two cut sets were identified. Sixty-three of these included only pathways encoded by fully annotated genes, and from this sub-set we selected five cut sets involved in amino acid or polyamine biosynthesis. One cut set (asnA/asnB) demonstrated redundancy in vitro and in vivo and showed that asparagine is essential for S. Typhimurium during infection. trpB/trpA as well as single mutants were attenuated for growth in vitro, while only the double mutant was a cut set in vivo, underlining previous observations that tryptophan is essential for successful outcome of infection. speB/speF,speC was not affected in vitro but was attenuated during infection showing that polyamines are essential for virulence apparently in a growth independent manner. The serA/glyA cut-set was found to be growth attenuated as predicted by the model. However, not only the double mutant, but also the glyA mutant, were found to be attenuated for virulence. This adds glycine production or conversion of glycine to THF to the list of essential reactions during infection. One pair (thrC/kbl) showed true redundancy in vitro but not in vivo demonstrating that threonine is available to the bacterium during infection. These data add to the existing knowledge of available nutrients in the intra-host environment, and have identified possible new targets for antibiotics.
Collapse
Affiliation(s)
- Lotte Jelsbak
- Department of Veterinary Disease Biology, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg, Denmark
| | - Hassan Hartman
- Department of Medical and Biological Sciences, Faculty of Health and Life Science, Oxford Brookes University, Oxford, United Kingdom
| | - Casper Schroll
- Department of Veterinary Disease Biology, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg, Denmark
| | - Jesper T. Rosenkrantz
- Department of Veterinary Disease Biology, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg, Denmark
| | - Sebastien Lemire
- Center for Systems Microbiology, Department of Systems Biology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Inke Wallrodt
- Department of Veterinary Disease Biology, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg, Denmark
| | - Line E. Thomsen
- Department of Veterinary Disease Biology, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg, Denmark
| | - Mark Poolman
- Department of Medical and Biological Sciences, Faculty of Health and Life Science, Oxford Brookes University, Oxford, United Kingdom
| | - Mogens Kilstrup
- Center for Systems Microbiology, Department of Systems Biology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Peter R. Jensen
- Center for Systems Microbiology, Department of Systems Biology, Technical University of Denmark, Kongens Lyngby, Denmark
| | - John E. Olsen
- Department of Veterinary Disease Biology, Faculty of Health and Medical Sciences, University of Copenhagen, Frederiksberg, Denmark
- * E-mail:
| |
Collapse
|
34
|
Hartman HB, Fell DA, Rossell S, Jensen PR, Woodward MJ, Thorndahl L, Jelsbak L, Olsen JE, Raghunathan A, Daefler S, Poolman MG. Identification of potential drug targets in Salmonella enterica sv. Typhimurium using metabolic modelling and experimental validation. MICROBIOLOGY-SGM 2014; 160:1252-1266. [PMID: 24777662 DOI: 10.1099/mic.0.076091-0] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Salmonella enterica sv. Typhimurium is an established model organism for Gram-negative, intracellular pathogens. Owing to the rapid spread of resistance to antibiotics among this group of pathogens, new approaches to identify suitable target proteins are required. Based on the genome sequence of S. Typhimurium and associated databases, a genome-scale metabolic model was constructed. Output was based on an experimental determination of the biomass of Salmonella when growing in glucose minimal medium. Linear programming was used to simulate variations in the energy demand while growing in glucose minimal medium. By grouping reactions with similar flux responses, a subnetwork of 34 reactions responding to this variation was identified (the catabolic core). This network was used to identify sets of one and two reactions that when removed from the genome-scale model interfered with energy and biomass generation. Eleven such sets were found to be essential for the production of biomass precursors. Experimental investigation of seven of these showed that knockouts of the associated genes resulted in attenuated growth for four pairs of reactions, whilst three single reactions were shown to be essential for growth.
Collapse
Affiliation(s)
- Hassan B Hartman
- Department of Medical and Biological Sciences, Oxford Brookes University, Gipsy Lane, Headington, Oxford OX3 OBP, UK
| | - David A Fell
- Department of Medical and Biological Sciences, Oxford Brookes University, Gipsy Lane, Headington, Oxford OX3 OBP, UK
| | - Sergio Rossell
- Department of Systems Biology, Technical University of Denmark, Lyngby, Denmark
| | - Peter Ruhdal Jensen
- Department of Systems Biology, Technical University of Denmark, Lyngby, Denmark
| | - Martin J Woodward
- Department of Food and Nutritional Sciences, University of Reading, Reading, UK
| | - Lotte Thorndahl
- Department of Veterinary Disease Biology, University of Copenhagen, Copenhagen, Denmark
| | - Lotte Jelsbak
- Department of Veterinary Disease Biology, University of Copenhagen, Copenhagen, Denmark
| | - John Elmerdahl Olsen
- Department of Veterinary Disease Biology, University of Copenhagen, Copenhagen, Denmark
| | - Anu Raghunathan
- Department of Infectious Diseases, Mount Sinai School of Medicine, New York, NY, USA
| | - Simon Daefler
- Department of Infectious Diseases, Mount Sinai School of Medicine, New York, NY, USA
| | - Mark G Poolman
- Department of Medical and Biological Sciences, Oxford Brookes University, Gipsy Lane, Headington, Oxford OX3 OBP, UK
| |
Collapse
|
35
|
Micheloni A, Orsi G, De Maria C, Vozzi G. ADMET: ADipocyte METabolism mathematical model. Comput Methods Biomech Biomed Engin 2014; 18:1386-91. [DOI: 10.1080/10255842.2014.908855] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
|
36
|
Kalnenieks U, Pentjuss A, Rutkis R, Stalidzans E, Fell DA. Modeling of Zymomonas mobilis central metabolism for novel metabolic engineering strategies. Front Microbiol 2014; 5:42. [PMID: 24550906 PMCID: PMC3914154 DOI: 10.3389/fmicb.2014.00042] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2013] [Accepted: 01/21/2014] [Indexed: 12/21/2022] Open
Abstract
Mathematical modeling of metabolism is essential for rational metabolic engineering. The present work focuses on several types of modeling approach to quantitative understanding of central metabolic network and energetics in the bioethanol-producing bacterium Zymomonas mobilis. Combined use of Flux Balance, Elementary Flux Mode, and thermodynamic analysis of its central metabolism, together with dynamic modeling of the core catabolic pathways, can help to design novel substrate and product pathways by systematically analyzing the solution space for metabolic engineering, and yields insights into the function of metabolic network, hardly achievable without applying modeling tools.
Collapse
Affiliation(s)
- Uldis Kalnenieks
- Institute of Microbiology and Biotechnology, University of LatviaRiga, Latvia
| | - Agris Pentjuss
- Department of Computer Systems, Latvia University of AgricultureJelgava, Latvia
| | - Reinis Rutkis
- Institute of Microbiology and Biotechnology, University of LatviaRiga, Latvia
| | - Egils Stalidzans
- Institute of Microbiology and Biotechnology, University of LatviaRiga, Latvia
- Department of Computer Systems, Latvia University of AgricultureJelgava, Latvia
- SIA TIBITJelgava, Latvia
| | - David A. Fell
- Department of Biological and Medical Sciences, Oxford Brookes UniversityOxford, UK
| |
Collapse
|
37
|
Investigating host-pathogen behavior and their interaction using genome-scale metabolic network models. Methods Mol Biol 2014; 1184:523-62. [PMID: 25048144 DOI: 10.1007/978-1-4939-1115-8_29] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Genome Scale Metabolic Modeling methods represent one way to compute whole cell function starting from the genome sequence of an organism and contribute towards understanding and predicting the genotype-phenotype relationship. About 80 models spanning all the kingdoms of life from archaea to eukaryotes have been built till date and used to interrogate cell phenotype under varying conditions. These models have been used to not only understand the flux distribution in evolutionary conserved pathways like glycolysis and the Krebs cycle but also in applications ranging from value added product formation in Escherichia coli to predicting inborn errors of Homo sapiens metabolism. This chapter describes a protocol that delineates the process of genome scale metabolic modeling for analysing host-pathogen behavior and interaction using flux balance analysis (FBA). The steps discussed in the process include (1) reconstruction of a metabolic network from the genome sequence, (2) its representation in a precise mathematical framework, (3) its translation to a model, and (4) the analysis using linear algebra and optimization. The methods for biological interpretations of computed cell phenotypes in the context of individual host and pathogen models and their integration are also discussed.
Collapse
|
38
|
Poolman MG, Kundu S, Shaw R, Fell DA. Metabolic trade-offs between biomass synthesis and photosynthate export at different light intensities in a genome-scale metabolic model of rice. FRONTIERS IN PLANT SCIENCE 2014; 5:656. [PMID: 25506349 PMCID: PMC4246663 DOI: 10.3389/fpls.2014.00656] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/27/2014] [Accepted: 11/04/2014] [Indexed: 05/08/2023]
Abstract
Previously we have used a genome scale model of rice metabolism to describe how metabolism reconfigures at different light intensities in an expanding leaf of rice. Although this established that the metabolism of the leaf was adequately represented, in the model, the scenario was not that of the typical function of the leaf-to provide material for the rest of the plant. Here we extend our analysis to explore the transition to a source leaf as export of photosynthate increases at the expense of making leaf biomass precursors, again as a function of light intensity. In particular we investigate whether, when the leaf is making a smaller range of compounds for export to the phloem, the same changes occur in the interactions between mitochondrial and chloroplast metabolism as seen in biomass synthesis for growth when light intensity increases. Our results show that the same changes occur qualitatively, though there are slight quantitative differences reflecting differences in the energy and redox requirements for the different metabolic outputs.
Collapse
Affiliation(s)
- Mark G. Poolman
- Cell Systems Modelling Group, Department of Biological and Medical Science, Oxford Brookes UniversityOxford, UK
| | - Sudip Kundu
- Department of Biophysics, Molecular Biology, and Bioinformatics, Calcutta UniversityKolkata, India
| | - Rahul Shaw
- Department of Biophysics, Molecular Biology, and Bioinformatics, Calcutta UniversityKolkata, India
| | - David A. Fell
- Cell Systems Modelling Group, Department of Biological and Medical Science, Oxford Brookes UniversityOxford, UK
- *Correspondence: David A. Fell, Department of Biological and Medical Science, Oxford Brookes University, Oxford OX3 0BP, UK e-mail:
| |
Collapse
|
39
|
Cheung CYM, Williams TCR, Poolman MG, Fell DA, Ratcliffe RG, Sweetlove LJ. A method for accounting for maintenance costs in flux balance analysis improves the prediction of plant cell metabolic phenotypes under stress conditions. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2013; 75:1050-61. [PMID: 23738527 DOI: 10.1111/tpj.12252] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/03/2013] [Revised: 05/23/2013] [Accepted: 05/30/2013] [Indexed: 05/24/2023]
Abstract
Flux balance models of metabolism generally utilize synthesis of biomass as the main determinant of intracellular fluxes. However, the biomass constraint alone is not sufficient to predict realistic fluxes in central heterotrophic metabolism of plant cells because of the major demand on the energy budget due to transport costs and cell maintenance. This major limitation can be addressed by incorporating transport steps into the metabolic model and by implementing a procedure that uses Pareto optimality analysis to explore the trade-off between ATP and NADPH production for maintenance. This leads to a method for predicting cell maintenance costs on the basis of the measured flux ratio between the oxidative steps of the oxidative pentose phosphate pathway and glycolysis. We show that accounting for transport and maintenance costs substantially improves the accuracy of fluxes predicted from a flux balance model of heterotrophic Arabidopsis cells in culture, irrespective of the objective function used in the analysis. Moreover, when the new method was applied to cells under control, elevated temperature and hyper-osmotic conditions, only elevated temperature led to a substantial increase in cell maintenance costs. It is concluded that the hyper-osmotic conditions tested did not impose a metabolic stress, in as much as the metabolic network is not forced to devote more resources to cell maintenance.
Collapse
Affiliation(s)
- C Y Maurice Cheung
- Department of Plant Sciences, University of Oxford, South Parks Road, Oxford, OX1 3RB, UK
| | | | | | | | | | | |
Collapse
|
40
|
Poolman MG, Kundu S, Shaw R, Fell DA. Responses to light intensity in a genome-scale model of rice metabolism. PLANT PHYSIOLOGY 2013; 162:1060-72. [PMID: 23640755 PMCID: PMC3668040 DOI: 10.1104/pp.113.216762] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/22/2013] [Accepted: 04/30/2013] [Indexed: 05/08/2023]
Abstract
We describe the construction and analysis of a genome-scale metabolic model representing a developing leaf cell of rice (Oryza sativa) primarily derived from the annotations in the RiceCyc database. We used flux balance analysis to determine that the model represents a network capable of producing biomass precursors (amino acids, nucleotides, lipid, starch, cellulose, and lignin) in experimentally reported proportions, using carbon dioxide as the sole carbon source. We then repeated the analysis over a range of photon flux values to examine responses in the solutions. The resulting flux distributions show that (1) redox shuttles between the chloroplast, cytosol, and mitochondrion may play a significant role at low light levels, (2) photorespiration can act to dissipate excess energy at high light levels, and (3) the role of mitochondrial metabolism is likely to vary considerably according to the balance between energy demand and availability. It is notable that these organelle interactions, consistent with many experimental observations, arise solely as a result of the need for mass and energy balancing without any explicit assumptions concerning kinetic or other regulatory mechanisms.
Collapse
Affiliation(s)
- Mark G Poolman
- Department of Biology and Medical Science, Oxford Brookes University, Headington, Oxford OX3 OBP, United Kingdom.
| | | | | | | |
Collapse
|
41
|
Collakova E, Yen JY, Senger RS. Are we ready for genome-scale modeling in plants? PLANT SCIENCE : AN INTERNATIONAL JOURNAL OF EXPERIMENTAL PLANT BIOLOGY 2012; 191-192:53-70. [PMID: 22682565 DOI: 10.1016/j.plantsci.2012.04.010] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2012] [Revised: 04/17/2012] [Accepted: 04/18/2012] [Indexed: 05/02/2023]
Abstract
As it is becoming easier and faster to generate various types of high-throughput data, one would expect that by now we should have a comprehensive systems-level understanding of biology, biochemistry, and physiology at least in major prokaryotic and eukaryotic model systems. Despite the wealth of available data, we only get a glimpse of what is going on at the molecular level from the global perspective. The major reason is the high level of cellular complexity and our limited ability to identify all (or at least important) components and their interactions in virtually infinite number of internal and external conditions. Metabolism can be modeled mathematically by the use of genome-scale models (GEMs). GEMs are in silico metabolic flux models derived from available genome annotation. These models predict the combination of flux values of a defined metabolic network given the influence of internal and external signals. GEMs have been successfully implemented to model bacterial metabolism for over a decade. However, it was not until 2009 when the first GEM for Arabidopsis thaliana cell-suspension cultures was generated. Genome-scale modeling ("GEMing") in plants brings new challenges primarily due to the missing components and complexity of plant cells represented by the existence of: (i) photosynthesis; (ii) compartmentation; (iii) variety of cell and tissue types; and (iv) diverse metabolic responses to environmental and developmental cues as well as pathogens, insects, and competing weeds. This review presents a critical discussion of the advantages of existing plant GEMs, while identifies key targets for future improvements. Plant GEMs tend to be accurate in predicting qualitative changes in selected aspects of central carbon metabolism, while secondary metabolism is largely neglected mainly due to the missing (unknown) genes and metabolites. As such, these models are suitable for exploring metabolism in plants grown in favorable conditions, but not in field-grown plants that have to cope with environmental changes in complex ecosystems. AraGEM is the first GEM describing a photosynthetic and photorespiring plant cell (Arabidopsis thaliana). We demonstrate the use of AraGEM given the current (limited) knowledge of plant metabolism and reveal the unexpected robustness of AraGEM by a series of in silico simulations. The major focus of these simulations is on the assessment of the: (i) network connectivity; (ii) influence of CO₂ and photon uptake rates on cellular growth rates and production of individual biomass components; and (iii) stability of plant central carbon metabolism with internal pH changes.
Collapse
Affiliation(s)
- Eva Collakova
- Department of Plant Pathology, Physiology, and Weed Science, 308 Latham Hall, Virginia Tech, Blacksburg, VA, USA.
| | | | | |
Collapse
|
42
|
Quantification of metabolism in Saccharomyces cerevisiae under hyperosmotic conditions using elementary mode analysis. J Ind Microbiol Biotechnol 2012; 39:927-41. [PMID: 22354733 DOI: 10.1007/s10295-012-1090-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2010] [Accepted: 01/14/2012] [Indexed: 10/28/2022]
Abstract
Yeast metabolism under hyperosmotic stress conditions was quantified using elementary mode analysis to obtain insights into the metabolic status of the cell. The fluxes of elementary modes were determined as solutions to a linear program that used the stoichiometry of the elementary modes as constraints. The analysis demonstrated that distinctly different sets of elementary modes operate under normal and hyperosmotic conditions. During the adaptation phase, elementary modes that only produce glycerol are active, while elementary modes that yield biomass, ethanol, and glycerol become active after the adaptive phase. The flux distribution in the metabolic network, calculated using the fluxes in the elementary modes, was employed to obtain the flux ratio at key nodes. At the glucose 6-phosphate (G6P) node, 25% of the carbon influx was diverted towards the pentose phosphate pathway under normal growth conditions, while only 0.3% of the carbon flux was diverted towards the pentose phosphate pathway during growth at 1 M NaCl, indicating that cell growth is arrested under hyperosmotic conditions. Further, objective functions were used in the linear program to obtain optimal solution spaces corresponding to the different accumulation rates. The analysis demonstrated that while biomass formation was optimal under normal growth conditions, glycerol synthesis was closer to optimal during adaptation to osmotic shock.
Collapse
|
43
|
Trinh CT, Thompson RA. Elementary mode analysis: a useful metabolic pathway analysis tool for reprograming microbial metabolic pathways. Subcell Biochem 2012; 64:21-42. [PMID: 23080244 DOI: 10.1007/978-94-007-5055-5_2] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Elementary mode analysis is a useful metabolic pathway analysis tool to characterize cellular metabolism. It can identify all feasible metabolic pathways known as elementary modes that are inherent to a metabolic network. Each elementary mode contains a minimal and unique set of enzymatic reactions that can support cellular functions at steady state. Knowledge of all these pathway options enables systematic characterization of cellular phenotypes, analysis of metabolic network properties (e.g. structure, regulation, robustness, and fragility), phenotypic behavior discovery, and rational strain design for metabolic engineering application. This chapter focuses on the application of elementary mode analysis to reprogram microbial metabolic pathways for rational strain design and the metabolic pathway evolution of designed strains.
Collapse
Affiliation(s)
- Cong T Trinh
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, TN, USA,
| | | |
Collapse
|
44
|
Abstract
The detection and analysis of structural invariants in cellular reaction networks is of central importance to achieve a more comprehensive understanding of metabolism. In this work, we review different kinds of structural invariants in reaction networks and their Petri net-based representation. In particular, we discuss invariants that can be obtained from the left and right null spaces of the stoichiometric matrix which correspond to conserved moieties (P-invariants) and elementary flux modes (EFMs, minimal T-invariants). While conserved moieties can be used to detect stoichiometric inconsistencies in reaction networks, EFMs correspond to a mathematically rigorous definition of the concept of a biochemical pathway. As outlined here, EFMs allow to devise strategies for strain improvement, to assess the robustness of metabolic networks subject to perturbations, and to analyze the information flow in regulatory and signaling networks. Another important aspect addressed by this review is the limitation of metabolic pathway analysis using EFMs to small or medium-scale reaction networks. We discuss two recently introduced approaches to circumvent these limitations. The first is an algorithm to enumerate a subset of EFMs in genome-scale metabolic networks starting from the EFM with the least number of reactions. The second approach, elementary flux pattern analysis, allows to analyze pathways through specific subsystems of genome-scale metabolic networks. In contrast to EFMs, elementary flux patterns much more accurately reflect the metabolic capabilities of a subsystem of metabolism as well as its integration into the entire system.
Collapse
|
45
|
Xavier D, Vázquez S, Higuera C, Morán F, Montero F. Tools-4-Metatool (T4M): online suite of web-tools to process stoichiometric network analysis data from Metatool. Biosystems 2011; 105:169-72. [PMID: 21554926 DOI: 10.1016/j.biosystems.2011.04.004] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2010] [Revised: 03/04/2011] [Accepted: 04/06/2011] [Indexed: 11/15/2022]
Abstract
Tools-4-Metatool (T4M) is a suite of web-tools, implemented in PERL, which analyses, parses, and manipulates files related to Metatool. Its main goal is to assist the work with Metatool. T4M has two major sets of tools: Analysis and Compare. Analysis visualizes the results of Metatool (convex basis, elementary flux modes, and enzyme subsets) and facilitates the study of metabolic networks. It is composed of five tools: MDigraph, MetaMatrix, CBGraph, EMGraph, and SortEM. Compare was developed to compare different Metatool results from different networks. This set consists of: Compara and ComparaSub which compare network subsets providing outputs in different formats and ComparaEM that seeks for identical elementary modes in two metabolic networks. The suite T4M also includes one script that generates Metatool input: CBasis2Metatool, based on a Metatool output file that is filtered by a list of convex basis' metabolites. Finally, the utility CheckMIn checks the consistency of the Metatool input file. T4M is available at http://solea.quim.ucm.es/t4m.
Collapse
Affiliation(s)
- Daniela Xavier
- Departamento de Bioquímica y Biología Molecular I, Universidad Complutense Madrid, Avd. Complutense s/n, 28040 Madrid, Spain
| | | | | | | | | |
Collapse
|
46
|
Klamt S, von Kamp A. An application programming interface for CellNetAnalyzer. Biosystems 2011; 105:162-8. [PMID: 21315797 DOI: 10.1016/j.biosystems.2011.02.002] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2010] [Revised: 01/24/2011] [Accepted: 02/02/2011] [Indexed: 11/16/2022]
Abstract
CellNetAnalyzer (CNA) is a MATLAB toolbox providing computational methods for studying structure and function of metabolic and cellular signaling networks. In order to allow non-experts to use these methods easily, CNA provides GUI-based interactive network maps as a means of parameter input and result visualization. However, with the availability of high-throughput data, there is a need to make CNA's functionality also accessible in batch mode for automatic data processing. Furthermore, as some algorithms of CNA are of general relevance for network analysis it would be desirable if they could be called as sub-routines by other applications. For this purpose, we developed an API (application programming interface) for CNA allowing users (i) to access the content of network models in CNA, (ii) to use CNA's network analysis capabilities independent of the GUI, and (iii) to interact with the GUI to facilitate the development of graphical plugins. Here we describe the organization of network projects in CNA and the application of the new API functions to these projects. This includes the creation of network projects from scratch, loading and saving of projects and scenarios, and the application of the actual analysis methods. Furthermore, API functions for the import/export of metabolic models in SBML format and for accessing the GUI are described. Lastly, two example applications demonstrate the use and versatile applicability of CNA's API. CNA is freely available for academic use and can be downloaded from http://www.mpi-magdeburg.mpg.de/projects/cna/cna.html.
Collapse
Affiliation(s)
- Steffen Klamt
- Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse, Magdeburg, Germany.
| | | |
Collapse
|
47
|
Orsi G, De Maria C, Guzzardi M, Vozzi F, Vozzi G. HEMETβ: improvement of hepatocyte metabolism mathematical model. Comput Methods Biomech Biomed Engin 2011; 14:837-51. [DOI: 10.1080/10255842.2010.497145] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
48
|
Rajvanshi M, Venkatesh KV. Phenotypic characterization of Corynebacterium glutamicum under osmotic stress conditions using elementary mode analysis. J Ind Microbiol Biotechnol 2010; 38:1345-57. [PMID: 21132515 DOI: 10.1007/s10295-010-0918-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2010] [Accepted: 11/18/2010] [Indexed: 11/25/2022]
Abstract
Corynebacterium glutamicum, a soil bacterium, is used to produce amino acids such as lysine and glutamate. C. glutamicum is often exposed to osmolality changes in its medium, and the bacterium has therefore evolved several adaptive response mechanisms to overcome them. In this study we quantify the metabolic response of C. glutamicum under osmotic stress using elementary mode analysis (EMA). Further, we obtain the optimal phenotypic space for the synthesis of lysine and formation of biomass. The analysis demonstrated that with increasing osmotic stress, the flux towards trehalose formation and energy-generating pathways increased, while the flux of anabolic reactions diminished. Nodal analysis indicated that glucose-6-phosphate, phosphoenol pyruvate, and pyruvate nodes were capable of adapting to osmotic stress, whereas the oxaloacetic acid node was relatively unresponsive. Fewer elementary modes were active under stress indicating the rigid behavior of the metabolism in response to high osmolality. Optimal phenotypic space analysis revealed that under normal conditions the organism optimized growth during the initial log phase and lysine and trehalose formation during the stationary phase. However, under osmotic stress, the analysis demonstrated that the organism operates under suboptimal conditions for growth, and lysine and trehalose formation.
Collapse
Affiliation(s)
- Meghna Rajvanshi
- Department of Biosciences and Bioengineering, Indian Institute of Technology, Bombay, Powai, Mumbai, 400076, India
| | | |
Collapse
|
49
|
Radhakrishnan D, Rajvanshi M, Venkatesh KV. Phenotypic characterization of Corynebacterium glutamicum using elementary modes towards synthesis of amino acids. SYSTEMS AND SYNTHETIC BIOLOGY 2010; 4:281-91. [PMID: 22132055 PMCID: PMC3065593 DOI: 10.1007/s11693-011-9073-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2010] [Revised: 11/10/2010] [Accepted: 02/03/2011] [Indexed: 10/18/2022]
Abstract
UNLABELLED Elementary flux mode (EFM) analysis is a powerful tool to represent the metabolic network structure and can be further utilized for flux analysis. The method enables characterization and quantification of feasible phenotypes in microbes. EFM analysis was employed to characterize the phenotype of Corynebacterium glutamicum to yield various amino acids. The metabolic network of C. glutamicum yielded 62 elementary modes by incorporating the accumulation of amino acids namely, lysine, alanine, valine, glutamine and glutamate. The analysis also allowed us to compute the maximum theoretical yield for the synthesis of various amino acids. These 62 elementary modes were further used to obtain optimal phenotypic space towards accumulation of biomass and lysine. The study indicated that the optimal solution space from 62 elementary modes forms a super space which incorporates various mutants including lysine producing strain of C. glutamicum. The analysis was also extended to obtain sensitivity of the network to variation in the stoichiometry of NADP in the definition of biomass. ELECTRONIC SUPPLEMENTARY MATERIAL The online version of this article (doi:10.1007/s11693-011-9073-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Devesh Radhakrishnan
- Department of Chemical Engineering, Indian Institute of Technology, Bombay, Powai, Mumbai, 400076 India
| | - Meghna Rajvanshi
- Department of Biosciences & Bioengineering, Indian Institute of Technology, Bombay, Powai, Mumbai, 400076 India
| | - K. V. Venkatesh
- Department of Chemical Engineering, Indian Institute of Technology, Bombay, Powai, Mumbai, 400076 India
- Department of Biosciences & Bioengineering, Indian Institute of Technology, Bombay, Powai, Mumbai, 400076 India
| |
Collapse
|
50
|
Schallau K, Junker BH. Simulating plant metabolic pathways with enzyme-kinetic models. PLANT PHYSIOLOGY 2010; 152:1763-71. [PMID: 20118273 PMCID: PMC2850014 DOI: 10.1104/pp.109.149237] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/11/2009] [Accepted: 01/27/2010] [Indexed: 05/17/2023]
Affiliation(s)
| | - Björn H. Junker
- Department of Physiology and Cell Biology, Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), 06466 Gatersleben, Germany
| |
Collapse
|