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.
Synopsis
Sex‐specific, whole‐body human metabolic models were developed and constrained with physiological, dietary, and metabolomic data. They recapitulate known whole‐body metabolic functions and enable mechanistic exploration of host‐microbiome co‐metabolism.
Sex‐specific whole‐body metabolic reconstructions represent the integrated function of 26 organs and six blood cell types.
Stoichiometric reconstructions of metabolism can be constrained with whole‐body physiological and metabolomic data to generate personalized models.
Whole‐body metabolic models recapitulate known inter‐organ metabolic cycles and energy use, successfully predict known biomarkers of inherited metabolic diseases, and explore potential host‐microbiome co‐metabolism.
Sex‐specific, whole‐body human metabolic models were developed and constrained with physiological, dietary, and metabolomic data. They recapitulate known whole‐body metabolic functions and enable mechanistic exploration of host‐microbiome co‐metabolism.
The human gut microbiome performs important functions in human health and disease. A classic example for host-gut microbial co-metabolism is host biosynthesis of primary bile acids and their ...subsequent deconjugation and transformation by the gut microbiome. To understand these system-level host-microbe interactions, a mechanistic, multi-scale computational systems biology approach that integrates the different types of omic data is needed. Here, we use a systematic workflow to computationally model bile acid metabolism in gut microbes and microbial communities.
Therefore, we first performed a comparative genomic analysis of bile acid deconjugation and biotransformation pathways in 693 human gut microbial genomes and expanded 232 curated genome-scale microbial metabolic reconstructions with the corresponding reactions (available at https://vmh.life ). We then predicted the bile acid biotransformation potential of each microbe and in combination with other microbes. We found that each microbe could produce maximally six of the 13 secondary bile acids in silico, while microbial pairs could produce up to 12 bile acids, suggesting bile acid biotransformation being a microbial community task. To investigate the metabolic potential of a given microbiome, publicly available metagenomics data from healthy Western individuals, as well as inflammatory bowel disease patients and healthy controls, were mapped onto the genomes of the reconstructed strains. We constructed for each individual a large-scale personalized microbial community model that takes into account strain-level abundances. Using flux balance analysis, we found considerable variation in the potential to deconjugate and transform primary bile acids between the gut microbiomes of healthy individuals. Moreover, the microbiomes of pediatric inflammatory bowel disease patients were significantly depleted in their bile acid production potential compared with that of controls. The contributions of each strain to overall bile acid production potential across individuals were found to be distinct between inflammatory bowel disease patients and controls. Finally, bottlenecks limiting secondary bile acid production potential were identified in each microbiome model.
This large-scale modeling approach provides a novel way of analyzing metagenomics data to accelerate our understanding of the metabolic interactions between the host and gut microbiomes in health and diseases states. Our models and tools are freely available to the scientific community.
Omics experiments are ubiquitous in biological studies, leading to a deluge of data. However, it is still challenging to connect changes in these data to changes in cell functions because of complex ...interdependencies between genes, proteins, and metabolites. Here, we present a framework allowing researchers to infer how metabolic functions change on the basis of omics data. To enable this, we curated and standardized lists of metabolic tasks that mammalian cells can accomplish. Genome-scale metabolic networks were used to define gene sets associated with each metabolic task. We further developed a framework to overlay omics data on these sets and predict pathway usage for each metabolic task. We demonstrated how this approach can be used to quantify metabolic functions of diverse biological samples from the single cell to whole tissues and organs by using multiple transcriptomic datasets. To facilitate its adoption, we integrated the approach into GenePattern (www.genepattern.org—CellFie).
Display omitted
•An alternative approach for the interpretation of omics data is proposed•A collection of tasks covering major metabolic activities in mammalian cells is presented•The framework predicts how metabolic functions change on the basis of omics data•Metabolic functions from single cells to tissues and organs can be quantified
The existence of complex interdependencies between genes, proteins, and metabolites challenge the interpretation of omics experiments. Data-driven approaches have been particularly useful for identifying gene sets of interest. However, it remains difficult to gain a mechanistic understanding of and to quantify a cell's functions from enriched ontology terms. Genome-scale systems biology models can be used to analyze these datasets, but they require specialized training and can take extensive effort to deploy. Here, we developed a framework to directly predict how changes in omics experiments correspond to cell or tissue functions, thereby facilitating phenotype-relevant interpretation of these complex datum types.
Richelle et al. present a framework to comprehensively quantify the propensity of a mammalian cell to be responsible for a metabolic function. This approach can be used to facilitate phenotype-relevant interpretation of transcriptomic datasets from the single-cell level of these biological entities to their organization in tissues and organs.
Investigation of molecular mechanisms of human disorders, especially rare diseases, require exploration of various knowledge repositories for building precise hypotheses and complex data ...interpretation. Recently, increasingly more resources offer diagrammatic representation of such mechanisms, including disease-dedicated schematics in pathway databases and disease maps. However, collection of knowledge across them is challenging, especially for research projects with limited manpower.
In this article we present an automated workflow for construction of maps of molecular mechanisms for rare diseases. The workflow requires a standardized definition of a disease using Orphanet or HPO identifiers to collect relevant genes and variants, and to assemble a functional, visual repository of related mechanisms, including data overlays. The diagrams composing the final map are unified to a common systems biology format from CellDesigner SBML, GPML and SBML+layout+render. The constructed resource contains disease-relevant genes and variants as data overlays for immediate visual exploration, including embedded genetic variant browser and protein structure viewer.
We demonstrate the functionality of our workflow on two examples of rare diseases: Kawasaki disease and retinitis pigmentosa. Two maps are constructed based on their corresponding identifiers. Moreover, for the retinitis pigmentosa use-case, we include a list of differentially expressed genes to demonstrate how to tailor the workflow using omics datasets.
In summary, our work allows for an ad-hoc construction of molecular diagrams combined from different sources, preserving their layout and graphical style, but integrating them into a single resource. This allows to reduce time consuming tasks of prototyping of a molecular disease map, enabling visual exploration, hypothesis building, data visualization and further refinement. The code of the workflow is open and accessible at https://gitlab.lcsb.uni.lu/minerva/automap/.
The application of constraint-based modeling to functionally analyze metagenomic data has been limited so far, partially due to the absence of suitable toolboxes.
To address this gap, we created a ...comprehensive toolbox to model (i) microbe-microbe and host-microbe metabolic interactions, and (ii) microbial communities using microbial genome-scale metabolic reconstructions and metagenomic data. The Microbiome Modeling Toolbox extends the functionality of the constraint-based reconstruction and analysis toolbox.
The Microbiome Modeling Toolbox and the tutorials at https://git.io/microbiomeModelingToolbox.
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.
Abstract
A multitude of factors contribute to complex diseases and can be measured with ‘omics’ methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the ...Virtual Metabolic Human (VMH, www.vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources ‘Human metabolism’, ‘Gut microbiome’, ‘Disease’, ‘Nutrition’, and ‘ReconMaps’. The VMH captures 5180 unique metabolites, 17 730 unique reactions, 3695 human genes, 255 Mendelian diseases, 818 microbes, 632 685 microbial genes and 8790 food items. The VMH’s unique features are (i) the hosting of the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; (ii) seven human metabolic maps for data visualization; (iii) a nutrition designer; (iv) a user-friendly webpage and application-programming interface to access its content; (v) user feedback option for community engagement and (vi) the connection of its entities to 57 other web resources. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community.
Flux balance analysis, and its variants, are widely used methods for predicting steady-state reaction rates in biochemical reaction networks. The exploration of high dimensional networks with such ...methods is currently hampered by software performance limitations.
DistributedFBA.jl is a high-level, high-performance, open-source implementation of flux balance analysis in Julia. It is tailored to solve multiple flux balance analyses on a subset or all the reactions of large and huge-scale networks, on any number of threads or nodes.
The code is freely available on github.com/opencobra/COBRA.jl. The documentation can be found at opencobra.github.io/COBRA.jl.
ronan.mt.fleming@gmail.com.
Flux balance analysis and its variants are widely used methods for predicting steady-state reaction rates in biochemical reaction networks. The exploration of high dimensional networks with such ...methods is currently hampered by software performance limitations.
DistributedFBA.jl is a high-level, high-performance, open-source implementation of flux balance analysis in Julia. It is tailored to solve multiple flux balance analyses on a subset or all the reactions of large and huge-scale networks, on any number of threads or nodes.
The code is freely available on github.com/opencobra/COBRA.jl. The documentation can be found at opencobra.github.io/COBRA.jl.
ronan.mt.fleming@gmail.com.
Supplementary data are available at Bioinformatics online.
We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS‐CoV‐2 infection, improving data interpretation and predicting key ...targets of intervention. Here, we describe a large‐scale community effort to build an open access, interoperable and computable repository of COVID‐19 molecular mechanisms. The COVID‐19 Disease Map (C19DMap) is a graphical, interactive representation of disease‐relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph‐based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS‐CoV‐2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID‐19 or similar pandemics in the long‐term perspective.
SYNOPSIS
COVID‐19 Disease Map is a large‐scale collection of curated computational models and diagrams of molecular mechanisms involved in SARS‐CoV‐2 infection. The map supports the computational exploration of pathways affected by the virus.
COVID‐19 Disease Map was built by over 20 independent biocuration teams and harmonised using systems biology standards.
Biocuration efforts were assisted by the systematic use of text‐ and AI‐assisted mining of relevant bioinformatic databases and platforms.
Case studies illustrate the applications of the map for visual exploration and computational analysis of SARS‐CoV‐2 pathways in combination with omic data.
The map is an open‐access effort, with all content and code shared in public repositories.
COVID‐19 Disease Map is a large‐scale collection of curated computational models and diagrams of molecular mechanisms involved in SARS‐CoV‐2 infection. The map supports the computational exploration of pathways affected by the virus.