Personalized whole-body models integrate metabolism, physiology, and the gut microbiome.

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.

Arterio-venous metabolomics exploration reveals major changes across liver and intestine in the obese Yucatan minipig.

Blood circulation mainly aims at distributing the nutrients required for tissue metabolism and collecting safely the by-products of all tissues to be further metabolized or eliminated. The simultaneous study of arterial (A) and venous (V) specific metabolites therefore has appeared to be a more relevant approach to understand and study the metabolism of a given organ. We propose to implement this approach by applying a metabolomics (NMR) strategy on paired AV blood across the intestine and liver on high fat/high sugar (HFHS)-fed minipigs. Our objective was to unravel kinetically and sequentially the metabolic adaptations to early obesity/insulin resistance onset specifically on these two tissues. After two months of HFHS feeding our study of AV ratios of the metabolome highlighted three major features. First, the hepatic metabolism switched from carbohydrate to lipid utilization. Second, the energy demand of the intestine increased, resulting in an enhanced uptake of glutamine, glutamate, and the recruitment of novel energy substrates (choline and creatine). Third, the uptake of methionine and threonine was considered to be driven by an increased intestine turnover to cope with the new high-density diet. Finally, the unique combination of experimental data and modelling predictions suggested that HFHS feeding was associated with changes in tryptophan metabolism and fatty acid ß-oxidation, which may play an important role in lipid hepatic accumulation and insulin sensitivity.

Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0.

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.

A systems approach reveals distinct metabolic strategies among the NCI-60 cancer cell lines.

The metabolic phenotype of cancer cells is reflected by the metabolites they consume and by the byproducts they release. Here, we use quantitative, extracellular metabolomic data of the NCI-60 panel and a novel computational method to generate 120 condition-specific cancer cell line metabolic models. These condition-specific cancer models used distinct metabolic strategies to generate energy and cofactors. The analysis of the models' capability to deal with environmental perturbations revealed three oxotypes, differing in the range of allowable oxygen uptake rates. Interestingly, models based on metabolomic profiles of melanoma cells were distinguished from other models through their low oxygen uptake rates, which were associated with a glycolytic phenotype. A subset of the melanoma cell models required reductive carboxylation. The analysis of protein and RNA expression levels from the Human Protein Atlas showed that IDH2, which was an essential gene in the melanoma models, but not IDH1 protein, was detected in normal skin cell types and melanoma. Moreover, the von Hippel-Lindau tumor suppressor (VHL) protein, whose loss is associated with non-hypoxic HIF-stabilization, reductive carboxylation, and promotion of glycolysis, was uniformly absent in melanoma. Thus, the experimental data supported the predicted role of IDH2 and the absence of VHL protein supported the glycolytic and low oxygen phenotype predicted for melanoma. Taken together, our approach of integrating extracellular metabolomic data with metabolic modeling and the combination of different network interrogation methods allowed insights into the metabolism of cells.

MetaboTools: A Comprehensive Toolbox for Analysis of Genome-Scale Metabolic Models.

Metabolomic data sets provide a direct read-out of cellular phenotypes and are increasingly generated to study biological questions. Previous work, by us and others, revealed the potential of analyzing extracellular metabolomic data in the context of the metabolic model using constraint-based modeling. With the MetaboTools, we make our methods available to the broader scientific community. The MetaboTools consist of a protocol, a toolbox, and tutorials of two use cases. The protocol describes, in a step-wise manner, the workflow of data integration, and computational analysis. The MetaboTools comprise the Matlab code required to complete the workflow described in the protocol. Tutorials explain the computational steps for integration of two different data sets and demonstrate a comprehensive set of methods for the computational analysis of metabolic models and stratification thereof into different phenotypes. The presented workflow supports integrative analysis of multiple omics data sets. Importantly, all analysis tools can be applied to metabolic models without performing the entire workflow. Taken together, the MetaboTools constitute a comprehensive guide to the intra-model analysis of extracellular metabolomic data from microbial, plant, or human cells. This computational modeling resource offers a broad set of computational analysis tools for a wide biomedical and non-biomedical research community.

Computational Modeling of Human Metabolism and Its Application to Systems Biomedicine.

Modern high-throughput techniques offer immense opportunities to investigate whole-systems behavior, such as those underlying human diseases. However, the complexity of the data presents challenges in interpretation, and new avenues are needed to address the complexity of both diseases and data. Constraint-based modeling is one formalism applied in systems biology. It relies on a genome-scale reconstruction that captures extensive biochemical knowledge regarding an organism. The human genome-scale metabolic reconstruction is increasingly used to understand normal cellular and disease states because metabolism is an important factor in many human diseases. The application of human genome-scale reconstruction ranges from mere querying of the model as a knowledge base to studies that take advantage of the model's topology and, most notably, to functional predictions based on cell- and condition-specific metabolic models built based on omics data.An increasing number and diversity of biomedical questions are being addressed using constraint-based modeling and metabolic models. One of the most successful biomedical applications to date is cancer metabolism, but constraint-based modeling also holds great potential for inborn errors of metabolism or obesity. In addition, it offers great prospects for individualized approaches to diagnostics and the design of disease prevention and intervention strategies. Metabolic models support this endeavor by providing easy access to complex high-throughput datasets. Personalized metabolic models have been introduced. Finally, constraint-based modeling can be used to model whole-body metabolism, which will enable the elucidation of metabolic interactions between organs and disturbances of these interactions as either causes or consequence of metabolic diseases. This chapter introduces constraint-based modeling and describes some of its contributions to systems biomedicine.

Prediction of intracellular metabolic states from extracellular metabolomic data.

Metabolic models can provide a mechanistic framework to analyze information-rich omics data sets, and are increasingly being used to investigate metabolic alternations in human diseases. An expression of the altered metabolic pathway utilization is the selection of metabolites consumed and released by cells. However, methods for the inference of intracellular metabolic states from extracellular measurements in the context of metabolic models remain underdeveloped compared to methods for other omics data. Herein, we describe a workflow for such an integrative analysis emphasizing on extracellular metabolomics data. We demonstrate, using the lymphoblastic leukemia cell lines Molt-4 and CCRF-CEM, how our methods can reveal differences in cell metabolism. Our models explain metabolite uptake and secretion by predicting a more glycolytic phenotype for the CCRF-CEM model and a more oxidative phenotype for the Molt-4 model, which was supported by our experimental data. Gene expression analysis revealed altered expression of gene products at key regulatory steps in those central metabolic pathways, and literature query emphasized the role of these genes in cancer metabolism. Moreover, in silico gene knock-outs identified unique control points for each cell line model, e.g., phosphoglycerate dehydrogenase for the Molt-4 model. Thus, our workflow is well-suited to the characterization of cellular metabolic traits based on extracellular metabolomic data, and it allows the integration of multiple omics data sets into a cohesive picture based on a defined model context.

Membrane transporters in a human genome-scale metabolic knowledgebase and their implications for disease.

Membrane transporters enable efficient cellular metabolism, aid in nutrient sensing, and have been associated with various diseases, such as obesity and cancer. Genome-scale metabolic network reconstructions capture genomic, physiological, and biochemical knowledge of a target organism, along with a detailed representation of the cellular metabolite transport mechanisms. Since the first reconstruction of human metabolism, Recon 1, published in 2007, progress has been made in the field of metabolite transport. Recently, we published an updated reconstruction, Recon 2, which significantly improved the metabolic coverage and functionality. Human metabolic reconstructions have been used to investigate the role of metabolism in disease and to predict biomarkers and drug targets. Given the importance of cellular transport systems in understanding human metabolism in health and disease, we analyzed the coverage of transport systems for various metabolite classes in Recon 2. We will review the current knowledge on transporters (i.e., their preferred substrates, transport mechanisms, metabolic relevance, and disease association for each metabolite class). We will assess missing coverage and propose modifications and additions through a transport module that is functional when combined with Recon 2. This information will be valuable for further refinements. These data will also provide starting points for further experiments by highlighting areas of incomplete knowledge. This review represents the first comprehensive overview of the transporters involved in central metabolism and their transport mechanisms, thus serving as a compendium of metabolite transporters specific for human metabolic reconstructions.

A community-driven global reconstruction of human metabolism.

Multiple models of human metabolism have been reconstructed, but each represents only a subset of our knowledge. Here we describe Recon 2, a community-driven, consensus 'metabolic reconstruction', which is the most comprehensive representation of human metabolism that is applicable to computational modeling. Compared with its predecessors, the reconstruction has improved topological and functional features, including ∼2× more reactions and ∼1.7× more unique metabolites. Using Recon 2 we predicted changes in metabolite biomarkers for 49 inborn errors of metabolism with 77% accuracy when compared to experimental data. Mapping metabolomic data and drug information onto Recon 2 demonstrates its potential for integrating and analyzing diverse data types. Using protein expression data, we automatically generated a compendium of 65 cell type-specific models, providing a basis for manual curation or investigation of cell-specific metabolic properties. Recon 2 will facilitate many future biomedical studies and is freely available at http://humanmetabolism.org/.

Contextualization procedure and modeling of monocyte specific TLR signaling.

Innate immunity is the first line of defense against invasion of pathogens. Toll-like receptor (TLR) signaling is involved in a variety of human diseases extending far beyond immune system-related diseases, affecting a number of different tissues and cell-types. Computational models often do not account for cell-type specific differences in signaling networks. Investigation of these differences and its phenotypic implications could increase understanding of cell signaling and processes such as inflammation. The wealth of knowledge for TLR signaling has been recently summarized in a stoichiometric signaling network applicable for constraint-based modeling and analysis (COBRA). COBRA methods have been applied to investigate tissue-specific metabolism using omics data integration. Comparable approaches have not been conducted using signaling networks. In this study, we present ihsTLRv2, an updated TLR signaling network accounting for the association of 314 genes with 558 network reactions. We present a mapping procedure for transcriptomic data onto signaling networks and demonstrate the generation of a monocyte-specific TLR network. The generated monocyte network is characterized through expression of a specific set of isozymes rather than reduction of pathway contents. While further tailoring the network to a specific stimulation condition, we observed that the quantitative changes in gene expression due to LPS stimulation affected the tightly connected set of genes. Differential expression influenced about one third of the entire TLR signaling network, in particular, NF-κB activation. Thus, a cell-type and condition-specific signaling network can provide functional insight into signaling cascades. Furthermore, we demonstrate the energy dependence of TLR signaling pathways in monocytes.

Feature-based attention and the suppression of non-relevant object features.

Feature-directed attention has been recently studied in various psychophysical, electrophysiological, and imaging studies. Convincing evidence has been obtained for its global effectiveness, but there is a debate about the processing fate of non-attended target features. A number of studies demonstrated feature-directed attention being associated with co-selection of non-relevant object features, thus resulting in selection of the entire object, whereas most other studies did not examine the extent to which processing of non-attended features was affected. Here, we present the results of two psychophysical experiments consisting of a Posner-like paradigm in which subjects were cued either to an individual feature or the entire object. We measured reaction times to changes in speed or colour of one of two simultaneously presented gratings. Our results strongly support the view that feature-based selection is a unique selection process different from object-based selection in that it can be associated with active suppression of non-relevant features.