1
|
Fleming RMT, Haraldsdottir HS, Minh LH, Vuong PT, Hankemeier T, Thiele I. Cardinality optimization in constraint-based modelling: application to human metabolism. Bioinformatics 2023; 39:btad450. [PMID: 37697651 PMCID: PMC10495685 DOI: 10.1093/bioinformatics/btad450] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Revised: 05/12/2023] [Indexed: 09/13/2023] Open
Abstract
MOTIVATION Several applications in constraint-based modelling can be mathematically formulated as cardinality optimization problems involving the minimization or maximization of the number of nonzeros in a vector. These problems include testing for stoichiometric consistency, testing for flux consistency, testing for thermodynamic flux consistency, computing sparse solutions to flux balance analysis problems and computing the minimum number of constraints to relax to render an infeasible flux balance analysis problem feasible. Such cardinality optimization problems are computationally complex, with no known polynomial time algorithms capable of returning an exact and globally optimal solution. RESULTS By approximating the zero-norm with nonconvex continuous functions, we reformulate a set of cardinality optimization problems in constraint-based modelling into a difference of convex functions. We implemented and numerically tested novel algorithms that approximately solve the reformulated problems using a sequence of convex programs. We applied these algorithms to various biochemical networks and demonstrate that our algorithms match or outperform existing related approaches. In particular, we illustrate the efficiency and practical utility of our algorithms for cardinality optimization problems that arise when extracting a model ready for thermodynamic flux balance analysis given a human metabolic reconstruction. AVAILABILITY AND IMPLEMENTATION Open source scripts to reproduce the results are here https://github.com/opencobra/COBRA.papers/2023_cardOpt with general purpose functions integrated within the COnstraint-Based Reconstruction and Analysis toolbox: https://github.com/opencobra/cobratoolbox.
Collapse
Affiliation(s)
- Ronan M T Fleming
- Metabolomics and Analytics Center, Leiden Academic Centre for Drug Research, Leiden University, Wassenaarseweg 76, Leiden 2333 CC, The Netherlands
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6 avenue du Swing, Belvaux L-4362, Luxembourg
- School of Medicine, National University of Ireland, University Rd, Galway H91 TK33, Ireland
| | - Hulda S Haraldsdottir
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6 avenue du Swing, Belvaux L-4362, Luxembourg
| | - Le Hoai Minh
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6 avenue du Swing, Belvaux L-4362, Luxembourg
| | - Phan Tu Vuong
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6 avenue du Swing, Belvaux L-4362, Luxembourg
- Mathematical Sciences School, University of Southampton, University Road, Southampton SO17 1BJ, United Kingdom
| | - Thomas Hankemeier
- Metabolomics and Analytics Center, Leiden Academic Centre for Drug Research, Leiden University, Wassenaarseweg 76, Leiden 2333 CC, The Netherlands
| | - Ines Thiele
- School of Medicine, National University of Ireland, University Rd, Galway H91 TK33, Ireland
| |
Collapse
|
2
|
Heinken A, Hertel J, Acharya G, Ravcheev DA, Nyga M, Okpala OE, Hogan M, Magnúsdóttir S, Martinelli F, Nap B, Preciat G, Edirisinghe JN, Henry CS, Fleming RMT, Thiele I. Genome-scale metabolic reconstruction of 7,302 human microorganisms for personalized medicine. Nat Biotechnol 2023; 41:1320-1331. [PMID: 36658342 PMCID: PMC10497413 DOI: 10.1038/s41587-022-01628-0] [Citation(s) in RCA: 42] [Impact Index Per Article: 42.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Accepted: 11/30/2022] [Indexed: 01/21/2023]
Abstract
The human microbiome influences the efficacy and safety of a wide variety of commonly prescribed drugs. Designing precision medicine approaches that incorporate microbial metabolism would require strain- and molecule-resolved, scalable computational modeling. Here, we extend our previous resource of genome-scale metabolic reconstructions of human gut microorganisms with a greatly expanded version. AGORA2 (assembly of gut organisms through reconstruction and analysis, version 2) accounts for 7,302 strains, includes strain-resolved drug degradation and biotransformation capabilities for 98 drugs, and was extensively curated based on comparative genomics and literature searches. The microbial reconstructions performed very well against three independently assembled experimental datasets with an accuracy of 0.72 to 0.84, surpassing other reconstruction resources and predicted known microbial drug transformations with an accuracy of 0.81. We demonstrate that AGORA2 enables personalized, strain-resolved modeling by predicting the drug conversion potential of the gut microbiomes from 616 patients with colorectal cancer and controls, which greatly varied between individuals and correlated with age, sex, body mass index and disease stages. AGORA2 serves as a knowledge base for the human microbiome and paves the way to personalized, predictive analysis of host-microbiome metabolic interactions.
Collapse
Affiliation(s)
- Almut Heinken
- School of Medicine, University of Galway, Galway, Ireland
- Ryan Institute, University of Galway, Galway, Ireland
- INSERM UMRS 1256, Nutrition, Genetics, and Environmental Risk Exposure (NGERE), University of Lorraine, Nancy, France
| | - Johannes Hertel
- School of Medicine, University of Galway, Galway, Ireland
- Department of Psychiatry and Psychotherapy, University Medicine Greifswald, Greifswald, Germany
| | - Geeta Acharya
- Integrated BioBank of Luxembourg, Dudelange, Luxembourg
| | - Dmitry A Ravcheev
- School of Medicine, University of Galway, Galway, Ireland
- Ryan Institute, University of Galway, Galway, Ireland
| | | | | | - Marcus Hogan
- School of Medicine, University of Galway, Galway, Ireland
- Ryan Institute, University of Galway, Galway, Ireland
| | - Stefanía Magnúsdóttir
- Center for Molecular Medicine, University Medical Center Utrecht, Utrecht, the Netherlands
| | - Filippo Martinelli
- School of Medicine, University of Galway, Galway, Ireland
- Ryan Institute, University of Galway, Galway, Ireland
| | - Bram Nap
- School of Medicine, University of Galway, Galway, Ireland
- Ryan Institute, University of Galway, Galway, Ireland
| | - German Preciat
- Leiden Academic Centre for Drug Research, Leiden University, Leiden, the Netherlands
| | - Janaka N Edirisinghe
- Computation Institute, University of Chicago, Chicago, IL, USA
- Mathematics and Computer Science Division, Argonne National Laboratory, Argonne, IL, USA
| | - Christopher S Henry
- Mathematics and Computer Science Division, Argonne National Laboratory, Argonne, IL, USA
| | - Ronan M T Fleming
- School of Medicine, University of Galway, Galway, Ireland
- Leiden Academic Centre for Drug Research, Leiden University, Leiden, the Netherlands
| | - Ines Thiele
- School of Medicine, University of Galway, Galway, Ireland.
- Ryan Institute, University of Galway, Galway, Ireland.
- Division of Microbiology, University of Galway, Galway, Ireland.
- APC Microbiome Ireland, Cork, Ireland.
| |
Collapse
|
3
|
Eicher T, Kinnebrew G, Patt A, Spencer K, Ying K, Ma Q, Machiraju R, Mathé EA. Metabolomics and Multi-Omics Integration: A Survey of Computational Methods and Resources. Metabolites 2020; 10:E202. [PMID: 32429287 PMCID: PMC7281435 DOI: 10.3390/metabo10050202] [Citation(s) in RCA: 61] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2020] [Revised: 05/07/2020] [Accepted: 05/13/2020] [Indexed: 02/06/2023] Open
Abstract
As researchers are increasingly able to collect data on a large scale from multiple clinical and omics modalities, multi-omics integration is becoming a critical component of metabolomics research. This introduces a need for increased understanding by the metabolomics researcher of computational and statistical analysis methods relevant to multi-omics studies. In this review, we discuss common types of analyses performed in multi-omics studies and the computational and statistical methods that can be used for each type of analysis. We pinpoint the caveats and considerations for analysis methods, including required parameters, sample size and data distribution requirements, sources of a priori knowledge, and techniques for the evaluation of model accuracy. Finally, for the types of analyses discussed, we provide examples of the applications of corresponding methods to clinical and basic research. We intend that our review may be used as a guide for metabolomics researchers to choose effective techniques for multi-omics analyses relevant to their field of study.
Collapse
Affiliation(s)
- Tara Eicher
- Biomedical Informatics Department, The Ohio State University College of Medicine, Columbus, OH 43210, USA; (T.E.); (G.K.); (K.S.); (Q.M.); (R.M.)
- Computer Science and Engineering Department, The Ohio State University College of Engineering, Columbus, OH 43210, USA
| | - Garrett Kinnebrew
- Biomedical Informatics Department, The Ohio State University College of Medicine, Columbus, OH 43210, USA; (T.E.); (G.K.); (K.S.); (Q.M.); (R.M.)
- Comprehensive Cancer Center, The Ohio State University and James Cancer Hospital, Columbus, OH 43210, USA;
- Bioinformatics Shared Resource Group, The Ohio State University, Columbus, OH 43210, USA
| | - Andrew Patt
- Division of Preclinical Innovation, National Center for Advancing Translational Sciences, NIH, 9800 Medical Center Dr., Rockville, MD, 20892, USA;
- Biomedical Sciences Graduate Program, The Ohio State University, Columbus, OH 43210, USA
| | - Kyle Spencer
- Biomedical Informatics Department, The Ohio State University College of Medicine, Columbus, OH 43210, USA; (T.E.); (G.K.); (K.S.); (Q.M.); (R.M.)
- Biomedical Sciences Graduate Program, The Ohio State University, Columbus, OH 43210, USA
- Nationwide Children’s Research Hospital, Columbus, OH 43210, USA
| | - Kevin Ying
- Comprehensive Cancer Center, The Ohio State University and James Cancer Hospital, Columbus, OH 43210, USA;
- Molecular, Cellular and Developmental Biology Program, The Ohio State University, Columbus, OH 43210, USA
| | - Qin Ma
- Biomedical Informatics Department, The Ohio State University College of Medicine, Columbus, OH 43210, USA; (T.E.); (G.K.); (K.S.); (Q.M.); (R.M.)
| | - Raghu Machiraju
- Biomedical Informatics Department, The Ohio State University College of Medicine, Columbus, OH 43210, USA; (T.E.); (G.K.); (K.S.); (Q.M.); (R.M.)
- Computer Science and Engineering Department, The Ohio State University College of Engineering, Columbus, OH 43210, USA
- Department of Pathology, Wexner Medical Center, The Ohio State University, Columbus, OH 43210, USA
- Translational Data Analytics Institute, The Ohio State University, Columbus, OH 43210, USA
| | - Ewy A. Mathé
- Biomedical Informatics Department, The Ohio State University College of Medicine, Columbus, OH 43210, USA; (T.E.); (G.K.); (K.S.); (Q.M.); (R.M.)
- Division of Preclinical Innovation, National Center for Advancing Translational Sciences, NIH, 9800 Medical Center Dr., Rockville, MD, 20892, USA;
| |
Collapse
|
4
|
Thiele I, Sahoo S, Heinken A, Hertel J, Heirendt L, Aurich MK, Fleming RMT. Personalized whole-body models integrate metabolism, physiology, and the gut microbiome. Mol Syst Biol 2020; 16:e8982. [PMID: 32463598 PMCID: PMC7285886 DOI: 10.15252/msb.20198982] [Citation(s) in RCA: 103] [Impact Index Per Article: 25.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2019] [Revised: 12/18/2019] [Accepted: 12/23/2019] [Indexed: 12/12/2022] Open
Abstract
Comprehensive molecular-level models of human metabolism have been generated on a cellular level. However, models of whole-body metabolism have not been established as they require new methodological approaches to integrate molecular and physiological data. We developed a new metabolic network reconstruction approach that used organ-specific information from literature and omics data to generate two sex-specific whole-body metabolic (WBM) reconstructions. These reconstructions capture the metabolism of 26 organs and six blood cell types. Each WBM reconstruction represents whole-body organ-resolved metabolism with over 80,000 biochemical reactions in an anatomically and physiologically consistent manner. We parameterized the WBM reconstructions with physiological, dietary, and metabolomic data. The resulting WBM models could recapitulate known inter-organ metabolic cycles and energy use. We also illustrate that the WBM models can predict known biomarkers of inherited metabolic diseases in different biofluids. Predictions of basal metabolic rates, by WBM models personalized with physiological data, outperformed current phenomenological models. Finally, integrating microbiome data allowed the exploration of host-microbiome co-metabolism. Overall, the WBM reconstructions, and their derived computational models, represent an important step toward virtual physiological humans.
Collapse
Affiliation(s)
- Ines Thiele
- School of MedicineNational University of IrelandGalwayIreland
- Discipline of MicrobiologySchool of Natural SciencesNational University of IrelandGalwayIreland
- APC MicrobiomeCorkIreland
- Luxembourg Centre for Systems BiomedicineUniversity of LuxembourgEsch‐sur‐AlzetteLuxembourg
| | - Swagatika Sahoo
- Luxembourg Centre for Systems BiomedicineUniversity of LuxembourgEsch‐sur‐AlzetteLuxembourg
- Present address:
Department of Chemical Engineering, and Initiative for Biological Systems Engineering (IBSE)Indian Institute of TechnologyChennaiIndia
| | - Almut Heinken
- School of MedicineNational University of IrelandGalwayIreland
| | - Johannes Hertel
- School of MedicineNational University of IrelandGalwayIreland
- Department of Psychiatry and PsychotherapyUniversity Medicine GreifswaldGreifswaldGermany
| | - Laurent Heirendt
- Luxembourg Centre for Systems BiomedicineUniversity of LuxembourgEsch‐sur‐AlzetteLuxembourg
| | - Maike K Aurich
- Luxembourg Centre for Systems BiomedicineUniversity of LuxembourgEsch‐sur‐AlzetteLuxembourg
| | - Ronan MT Fleming
- School of MedicineNational University of IrelandGalwayIreland
- Division of Analytical BiosciencesLeiden Academic Centre for Drug ResearchFaculty of ScienceUniversity of LeidenLeidenThe Netherlands
| |
Collapse
|
5
|
Structural conserved moiety splitting of a stoichiometric matrix. J Theor Biol 2020; 499:110276. [PMID: 32333975 DOI: 10.1016/j.jtbi.2020.110276] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2019] [Revised: 04/03/2020] [Accepted: 04/06/2020] [Indexed: 11/20/2022]
Abstract
Characterising biochemical reaction network structure in mathematical terms enables the inference of functional biochemical consequences from network structure with existing mathematical techniques and spurs the development of new mathematics that exploits the peculiarities of biochemical network structure. The structure of a biochemical network may be specified by reaction stoichiometry, that is, the relative quantities of each molecule produced and consumed in each reaction of the network. A biochemical network may also be specified at a higher level of resolution in terms of the internal structure of each molecule and how molecular structures are transformed by each reaction in a network. The stoichiometry for a set of reactions can be compiled into a stoichiometric matrix N∈Zm×n, where each row corresponds to a molecule and each column corresponds to a reaction. We demonstrate that a stoichiometric matrix may be split into the sum of m-rank(N) moiety transition matrices, each of which corresponds to a subnetwork accessible to a structurally identifiable conserved moiety. The existence of this moiety matrix splitting is a property that distinguishes a stoichiometric matrix from an arbitrary rectangular matrix.
Collapse
|
6
|
Magnúsdóttir S, Heinken A, Fleming RMT, Thiele I. Reply to "Challenges in modeling the human gut microbiome". Nat Biotechnol 2019; 36:686-691. [PMID: 30080835 DOI: 10.1038/nbt.4212] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Stefanía Magnúsdóttir
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Esch-sur-Alzette, Luxembourg
| | - Almut Heinken
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Esch-sur-Alzette, Luxembourg
| | - Ronan M T Fleming
- Division of Analytical Biosciences, Leiden Academic Centre for Drug Research, Faculty of Science, University of Leiden, Leiden, The Netherlands
| | - Ines Thiele
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Esch-sur-Alzette, Luxembourg
| |
Collapse
|
7
|
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: 609] [Impact Index Per Article: 121.8] [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
|
8
|
Magnúsdóttir S, Thiele I. Modeling metabolism of the human gut microbiome. Curr Opin Biotechnol 2017; 51:90-96. [PMID: 29258014 DOI: 10.1016/j.copbio.2017.12.005] [Citation(s) in RCA: 77] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2017] [Revised: 12/05/2017] [Accepted: 12/06/2017] [Indexed: 12/15/2022]
Abstract
The human gut microbiome plays an important part in human health. The complexity of the microbiome makes it difficult to determine the detailed metabolic functions and cross-talk occurs between the individual species. In silico systems biology studies of the microbiome can help to identify metabolite exchanges among gut microbes. Constraint-based reconstruction and analysis methods use biochemically accurate genome-scale metabolic networks of microorganisms to simulate metabolism between species in a given microbiome and help generate novel hypotheses on microbial interactions. Here, we review metabolic modeling studies that have investigated metabolic functions of the gut microbiome.
Collapse
Affiliation(s)
- Stefanía Magnúsdóttir
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Esch-sur-Alzette, Luxembourg
| | - Ines Thiele
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Esch-sur-Alzette, Luxembourg.
| |
Collapse
|
9
|
Cannon WR, Baker SE. Non-steady state mass action dynamics without rate constants: dynamics of coupled reactions using chemical potentials. Phys Biol 2017; 14:055003. [PMID: 28675379 DOI: 10.1088/1478-3975/aa7d80] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Comprehensive and predictive simulation of coupled reaction networks has long been a goal of biology and other fields. Currently, metabolic network models that utilize enzyme mass action kinetics have predictive power but are limited in scope and application by the fact that the determination of enzyme rate constants is laborious and low throughput. We present a statistical thermodynamic formulation of the law of mass action for coupled reactions at both steady states and non-stationary states. The formulation uses chemical potentials instead of rate constants. When used to model deterministic systems, the method corresponds to a rescaling of the time dependent reactions in such a way that steady states can be reached on the same time scale but with significantly fewer computational steps. The relationships between reaction affinities, free energy changes and generalized detailed balance are central to the discussion. The significance for applications in systems biology are discussed as is the concept and assumption of maximum entropy production rate as a biological principle that links thermodynamics to natural selection.
Collapse
Affiliation(s)
- William R Cannon
- Biological Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Blvd, Richland, WA 99352, United States of America. Author to whom any correspondence should be addressed
| | | |
Collapse
|
10
|
Ma D, Yang L, Fleming RMT, Thiele I, Palsson BO, Saunders MA. Reliable and efficient solution of genome-scale models of Metabolism and macromolecular Expression. Sci Rep 2017; 7:40863. [PMID: 28098205 PMCID: PMC5241643 DOI: 10.1038/srep40863] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2016] [Accepted: 12/13/2016] [Indexed: 12/25/2022] Open
Abstract
Constraint-Based Reconstruction and Analysis (COBRA) is currently the only methodology that permits integrated modeling of Metabolism and macromolecular Expression (ME) at genome-scale. Linear optimization computes steady-state flux solutions to ME models, but flux values are spread over many orders of magnitude. Data values also have greatly varying magnitudes. Standard double-precision solvers may return inaccurate solutions or report that no solution exists. Exact simplex solvers based on rational arithmetic require a near-optimal warm start to be practical on large problems (current ME models have 70,000 constraints and variables and will grow larger). We have developed a quadruple-precision version of our linear and nonlinear optimizer MINOS, and a solution procedure (DQQ) involving Double and Quad MINOS that achieves reliability and efficiency for ME models and other challenging problems tested here. DQQ will enable extensive use of large linear and nonlinear models in systems biology and other applications involving multiscale data.
Collapse
Affiliation(s)
- Ding Ma
- Stanford University, Dept of Management Science and Engineering, Stanford, CA 94305, USA
| | - Laurence Yang
- University of California, San Diego, Dept of Bioengineering, La Jolla, CA 92093, USA
| | - Ronan M T Fleming
- University of Luxembourg, Luxembourg Centre for Systems Biomedicine, L-4365 Esch-sur-Alzette, Luxembourg
| | - Ines Thiele
- University of Luxembourg, Luxembourg Centre for Systems Biomedicine, L-4365 Esch-sur-Alzette, Luxembourg
| | - Bernhard O Palsson
- University of California, San Diego, Dept of Bioengineering, La Jolla, CA 92093, USA.,Technical University of Denmark, Novo Nordisk Foundation Center for Biosustainability, 2970 Hørsholm, Denmark
| | - Michael A Saunders
- Stanford University, Dept of Management Science and Engineering, Stanford, CA 94305, USA
| |
Collapse
|
11
|
Magnúsdóttir S, Heinken A, Kutt L, Ravcheev DA, Bauer E, Noronha A, Greenhalgh K, Jäger C, Baginska J, Wilmes P, Fleming RMT, Thiele I. Generation of genome-scale metabolic reconstructions for 773 members of the human gut microbiota. Nat Biotechnol 2016; 35:81-89. [DOI: 10.1038/nbt.3703] [Citation(s) in RCA: 434] [Impact Index Per Article: 54.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2016] [Accepted: 09/20/2016] [Indexed: 02/06/2023]
|