Compartmentalization of the Edinburgh Human Metabolic Network

Unraveling Connective Tissue Growth Factor as a Therapeutic Target and Assessing Kahweol as a Potential Drug Candidate in Triple-Negative Breast Cancer Treatment

Triple-negative breast cancer (TNBC) is characterized by aggressive behavior and limited treatment options, necessitating the identification of novel therapeutic targets. In this study, we investigated the clinical significance of connective tissue growth factor (CTGF) as a prognostic marker and explored the potential therapeutic effects of kahweol, a coffee diterpene molecule, in TNBC treatment. Initially, through a survival analysis on breast cancer patients from The Cancer Genome Atlas (TCGA) database, we found that CTGF exhibited significant prognostic effects exclusively in TNBC patients. To gain mechanistic insights, we performed the functional annotation and gene set enrichment analyses, revealing the involvement of CTGF in migratory pathways relevant to TNBC treatment. Subsequently, in vitro experiments using MDA-MB 231 cells, a representative TNBC cell line, demonstrated that recombinant CTGF (rCTGF) administration enhanced cell motility, whereas CTGF knockdown using CTGF siRNA resulted in reduced motility. Notably, rCTGF restored kahweol-reduced cell motility, providing compelling evidence for the role of CTGF in mediating kahweol's effects. At the molecular level, kahweol downregulated the protein expression of CTGF as well as critical signaling molecules, such as p-ERK, p-P38, p-PI3K/AKT, and p-FAK, associated with cell motility. In summary, our findings propose CTGF as a potential prognostic marker for guiding TNBC treatment and suggest kahweol as a promising antitumor compound capable of regulating CTGF expression to suppress cell motility in TNBC. These insights hold promise for the development of targeted therapies and improved clinical outcomes for TNBC patients.

A new metabolic model of Drosophila melanogaster and the integrative analysis of Parkinson's disease

High conservation of the disease-associated genes between flies and humans facilitates the common use of Drosophila melanogaster to study metabolic disorders under controlled laboratory conditions. However, metabolic modeling studies are highly limited for this organism. We here report a comprehensively curated genome-scale metabolic network model of Drosophila using an orthology-based approach. The gene coverage and metabolic information of the draft model derived from a reference human model were expanded via Drosophila-specific KEGG and MetaCyc databases, with several curation steps to avoid metabolic redundancy and stoichiometric inconsistency. Furthermore, we performed literature-based curations to improve gene-reaction associations, subcellular metabolite locations, and various metabolic pathways. The performance of the resulting Drosophila model (8,230 reactions, 6,990 metabolites, and 2,388 genes), iDrosophila1 (https://github.com/SysBioGTU/iDrosophila), was assessed using flux balance analysis in comparison with the other currently available fly models leading to superior or comparable results. We also evaluated the transcriptome-based prediction capacity of iDrosophila1, where differential metabolic pathways during Parkinson's disease could be successfully elucidated. Overall, iDrosophila1 is promising to investigate system-level metabolic alterations in response to genetic and environmental perturbations.

Principles and functions of metabolic compartmentalization

Metabolism has historically been studied at the levels of whole cells, whole tissues and whole organisms. As a result, our understanding of how compartmentalization-the spatial and temporal separation of pathways and components-shapes organismal metabolism remains limited. At its essence, metabolic compartmentalization fulfils three important functions or 'pillars': establishing unique chemical environments, providing protection from reactive metabolites and enabling the regulation of metabolic pathways. However, how these pillars are established, regulated and maintained at both the cellular and systemic levels remains unclear. Here we discuss how the three pillars are established, maintained and regulated within the cell and discuss the consequences of dysregulation of metabolic compartmentalization in human disease. Organelles are increasingly emerging as 'command-and-control centres' and the increased understanding of metabolic compartmentalization is revealing new aspects of metabolic homeostasis, with this knowledge being translated into therapies for the treatment of cancer and certain neurodegenerative diseases.

Reconstruction of Litopenaeus vannamei Genome-Scale Metabolic Network Model and Nutritional Requirements Analysis of Different Shrimp Commercial Varieties

As an important tool for systematic analysis, genome-scale metabolic network (GSMN) model has been widely used in various organisms. However, there are few reports on the GSMNs of aquatic crustaceans. Litopenaeus vannamei is the largest and most productive shrimp species. Feed improvement is one of the important methods to improve the yield of L. vannamei and control water pollution caused by the inadequate absorption of feed. In this work, the first L. vannamei GSMN named iGH3005 was reconstructed and applied to the optimization of feed. iGH3005 was reconstructed based on the genomic data. The model includes 2,292 reactions and 3,005 genes. iGH3005 was used to analyze the nutritional requirements of five different L. vannamei commercial varieties and the genes influencing the metabolism of the nutrients. Based on the simulation, we found that tyrosine-protein kinase src64b like may catalyze different reactions in different commercial varieties. The preference of carbohydrate utilization is different in various commercial varieties, which may due to the different expressions of some genes. In addition, this investigation suggests that a rational and targeted modification in the macronutrient content of shrimp feed would lead to an increase in growth and feed conversion rate. The feed for different commercial varieties should be adjusted accordingly, and possible adjustment schemes were provided. The results of this work provided important information for physiological research and optimization of the components in feed of L. vannamei.

Reconstruction of Eriocheir sinensis Y-organ Genome-Scale Metabolic Network and Differential Analysis After Eyestalk Ablation

Genome-scale metabolic network (GSMN) has been proven to be a useful tool for the system analysis of organism metabolism and applied to deeply explore the metabolic functions or mechanisms in many organisms, including model or non-model organisms. However, the systematic studies on the metabolisms of aquatic animals are seldom reported, especially the aquatic crustaceans. In this work, we reconstructed an Eriocheir sinensis Y-organ GSMN based on the transcriptome sequencing of Y-organ, which includes 1,645 reactions, 1,885 unigenes, and 1,524 metabolites distributed in 100 pathways and 11 subsystems. Functional module and centrality analysis of the GSMN show the main metabolic functions of Y-organ. Further analysis of the differentially expressed unigenes in Y-organ after eyestalk ablation reveals that 191 genes in the network were up-regulated and 283 were down-regulated. The unigenes associated with the ecdysone synthetic pathway were all up-regulated, which is consistent with the report on the increasing secretion of ecdysone after eyestalk ablation. Besides, we compared the Y-organ GSMN with that of E. sinensis eyestalk and hepatopancreas, and we analyzed the specific metabolisms in each organ. The specific metabolisms and pathways of these three networks are closely related to their corresponding metabolic functions. The GSMN reconstructed in this work provides a new method and many novel clues for further understanding the physiological function of Y-organ. It also supplies a new platform for the investigation of the interactions among different organs in the growth process of E. sinensis.

Pathway Analysis for Targeted and Untargeted Metabolomics

Recent advances in analytical techniques, particularly LC-MS, generate increasingly large and complex metabolomics datasets. Pathway analysis tools help place the experimental observations into relevant biological or disease context. This chapter provides an overview of the general concepts and common tools for pathway analysis, including Mummichog for untargeted metabolomics. Examples of pathway mapping, MetScape, and Mummichog are explained. This serves as both a practical tutorial and a timely survey of pathway analysis for label-free metabolomics data.

Compartment and hub definitions tune metabolic networks for metabolomic interpretations

BACKGROUND

Metabolic networks represent all chemical reactions that occur between molecular metabolites in an organism's cells. They offer biological context in which to integrate, analyze, and interpret omic measurements, but their large scale and extensive connectivity present unique challenges. While it is practical to simplify these networks by placing constraints on compartments and hubs, it is unclear how these simplifications alter the structure of metabolic networks and the interpretation of metabolomic experiments.

RESULTS

We curated and adapted the latest systemic model of human metabolism and developed customizable tools to define metabolic networks with and without compartmentalization in subcellular organelles and with or without inclusion of prolific metabolite hubs. Compartmentalization made networks larger, less dense, and more modular, whereas hubs made networks larger, more dense, and less modular. When present, these hubs also dominated shortest paths in the network, yet their exclusion exposed the subtler prominence of other metabolites that are typically more relevant to metabolomic experiments. We applied the non-compartmental network without metabolite hubs in a retrospective, exploratory analysis of metabolomic measurements from 5 studies on human tissues. Network clusters identified individual reactions that might experience differential regulation between experimental conditions, several of which were not apparent in the original publications.

CONCLUSIONS

Exclusion of specific metabolite hubs exposes modularity in both compartmental and non-compartmental metabolic networks, improving detection of relevant clusters in omic measurements. Better computational detection of metabolic network clusters in large data sets has potential to identify differential regulation of individual genes, transcripts, and proteins.

Integration of Metabolomic and Other Omics Data in Population-Based Study Designs: An Epidemiological Perspective

It is not controversial that study design considerations and challenges must be addressed when investigating the linkage between single omic measurements and human phenotypes. It follows that such considerations are just as critical, if not more so, in the context of multi-omic studies. In this review, we discuss (1) epidemiologic principles of study design, including selection of biospecimen source(s) and the implications of the timing of sample collection, in the context of a multi-omic investigation, and (2) the strengths and limitations of various techniques of data integration across multi-omic data types that may arise in population-based studies utilizing metabolomic data.

Human Systems Biology and Metabolic Modelling: A Review-From Disease Metabolism to Precision Medicine

In cell and molecular biology, metabolism is the only system that can be fully simulated at genome scale. Metabolic systems biology offers powerful abstraction tools to simulate all known metabolic reactions in a cell, therefore providing a snapshot that is close to its observable phenotype. In this review, we cover the 15 years of human metabolic modelling. We show that, although the past five years have not experienced large improvements in the size of the gene and metabolite sets in human metabolic models, their accuracy is rapidly increasing. We also describe how condition-, tissue-, and patient-specific metabolic models shed light on cell-specific changes occurring in the metabolic network, therefore predicting biomarkers of disease metabolism. We finally discuss current challenges and future promising directions for this research field, including machine/deep learning and precision medicine. In the omics era, profiling patients and biological processes from a multiomic point of view is becoming more common and less expensive. Starting from multiomic data collected from patients and N-of-1 trials where individual patients constitute different case studies, methods for model-building and data integration are being used to generate patient-specific models. Coupled with state-of-the-art machine learning methods, this will allow characterizing each patient's disease phenotype and delivering precision medicine solutions, therefore leading to preventative medicine, reduced treatment, and in silico clinical trials.

Elucidation of Novel Therapeutic Targets for Acute Myeloid Leukemias with RUNX1-RUNX1T1 Fusion

The RUNX1-RUNX1T1 fusion is a frequent chromosomal alteration in acute myeloid leukemias (AMLs). Although RUNX1-RUNX1T1 fusion protein has pivotal roles in the development of AMLs with the fusion, RUNX1-RUNX1T1, fusion protein is difficult to target, as it lacks kinase activities. Here, we used bioinformatic tools to elucidate targetable signaling pathways in AMLs with RUNX1-RUNX1T1 fusion. After analysis of 93 AML cases from The Cancer Genome Atlas (TCGA) database, we found expression of 293 genes that correlated to the expression of the RUNX1-RUNX1T1 fusion gene. Based on these 293 genes, the cyclooxygenase (COX), vascular endothelial growth factor receptor (VEGFR), platelet-derived growth factor receptor (PDGFR), and fibroblast growth factor receptor (FGFR) pathways were predicted to be specifically activated in AMLs with RUNX1-RUNX1T1 fusion. Moreover, the in vitro proliferation of AML cells with RUNX1-RUNX1T1 fusion decreased significantly more than that of AML cells without the fusion, when the pathways were inhibited pharmacologically. The results indicate that novel targetable signaling pathways could be identified by the analysis of the gene expression features of AMLs with non-targetable genetic alterations. The elucidation of specific molecular targets for AMLs that have a specific genetic alteration would promote personalized treatment of AMLs and improve clinical outcomes.

Treatment Avenues in Myalgic Encephalomyelitis/Chronic Fatigue Syndrome: A Split-gender Pharmacogenomic Study of Gene-expression Modules

PURPOSE

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a debilitating multisymptom illness impacting up to 1 million people in the United States. As the pathogenesis and etiology of this complex condition are unclear, prospective treatments are limited. Identifying US Food and Drug Administration-approved drugs that may be repositioned as treatments for ME/CFS may offer a rapid and cost-effective solution.

METHODS

Here we used gene-expression data from 33 patients with Fukuda-defined ME/CFS (23 females, 10 males) and 21 healthy demographically comparable controls (15 females, 6 males) to identify differential expression of predefined gene-module sets based on nonparametric statistics. Differentially expressed gene modules were then annotated via over-representation analysis using the Consensus Pathway database. Differentially expressed modules were then regressed onto measures of fatigue and cross-referenced with drug atlas and pharmacogenomics databases to identify putative treatment agents.

FINDINGS

The top 1% of modules identified in males indicated small effect sizes in modules associated with immune regulation and mitochondrial dysfunction. In females, modules identified included those related to immune factors and cardiac/blood factors, returning effect sizes ranging from very small to intermediate (0.147 < Cohen δ < 0.532). Regression analysis indicated that B-cell receptors, T-cell receptors, tumor necrosis factor α, transforming growth factor β, and metabolic and cardiac modules were strongly correlated with multiple composite measures of fatigue. Cross-referencing identified genes with pharmacogenomics data indicated immunosuppressants as potential treatments of ME/CFS symptoms.

IMPLICATIONS

The findings from our analysis suggest that ME/CFS symptoms are perpetuated by immune dysregulation that may be approached via immune modulation-based treatment strategies.

Genetic Predisposition for Immune System, Hormone, and Metabolic Dysfunction in Myalgic Encephalomyelitis/Chronic Fatigue Syndrome: A Pilot Study

Introduction: Myalgic Encephalomyelitis/ Chronic Fatigue Syndrome (ME/CFS) is a multifactorial illness of unknown etiology with considerable social and economic impact. To investigate a putative genetic predisposition to ME/CFS we conducted genome-wide single-nucleotide polymorphism (SNP) analysis to identify possible variants. Methods: 383 ME/CFS participants underwent DNA testing using the commercial company 23andMe. The deidentified genetic data was then filtered to include only non-synonymous and nonsense SNPs from exons and microRNAs, and SNPs close to splice sites. The frequencies of each SNP were calculated within our cohort and compared to frequencies from the Kaviar reference database. Functional annotation of pathway sets containing SNP genes with high frequency in ME/CFS was performed using over-representation analysis via ConsensusPathDB. Furthermore, these SNPs were also scored using the Combined Annotation Dependent Depletion (CADD) algorithm to gauge their deleteriousness. Results: 5693 SNPs were found to have at least 10% frequency in at least one cohort (ME/CFS or reference) and at least two-fold absolute difference for ME/CFS. Functional analysis identified the majority of SNPs as related to immune system, hormone, metabolic, and extracellular matrix organization. CADD scoring identified 517 SNPs in these pathways that are among the 10% most deleteriousness substitutions to the human genome.

Sparse network modeling and metscape-based visualization methods for the analysis of large-scale metabolomics data

Motivation

Recent technological advances in mass spectrometry, development of richer mass spectral libraries and data processing tools have enabled large scale metabolic profiling. Biological interpretation of metabolomics studies heavily relies on knowledge-based tools that contain information about metabolic pathways. Incomplete coverage of different areas of metabolism and lack of information about non-canonical connections between metabolites limits the scope of applications of such tools. Furthermore, the presence of a large number of unknown features, which cannot be readily identified, but nonetheless can represent bona fide compounds, also considerably complicates biological interpretation of the data.

Results

Leveraging recent developments in the statistical analysis of high-dimensional data, we developed a new Debiased Sparse Partial Correlation algorithm (DSPC) for estimating partial correlation networks and implemented it as a Java-based CorrelationCalculator program. We also introduce a new version of our previously developed tool Metscape that enables building and visualization of correlation networks. We demonstrate the utility of these tools by constructing biologically relevant networks and in aiding identification of unknown compounds.

Availability and Implementation

http://metscape.med.umich.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Modelling pyruvate dehydrogenase under hypoxia and its role in cancer metabolism

Metabolism is the only biological system that can be fully modelled at genome scale. As a result, metabolic models have been increasingly used to study the molecular mechanisms of various diseases. Hypoxia, a low-oxygen tension, is a well-known characteristic of many cancer cells. Pyruvate dehydrogenase (PDH) controls the flux of metabolites between glycolysis and the tricarboxylic acid cycle and is a key enzyme in metabolic reprogramming in cancer metabolism. Here, we develop and manually curate a constraint-based metabolic model to investigate the mechanism of pyruvate dehydrogenase under hypoxia. Our results characterize the activity of pyruvate dehydrogenase and its decline during hypoxia. This results in lactate accumulation, consistent with recent hypoxia studies and a well-known feature in cancer metabolism. We apply machine-learning techniques on the flux datasets to identify reactions that drive these variations. We also identify distinct features on the structure of the variables and individual metabolic components in the switch from normoxia to hypoxia. Our results provide a framework for future studies by integrating multi-omics data to predict condition-specific metabolic phenotypes under hypoxia.

Compartmentalized metabolic network reconstruction of microbial communities to determine the effect of agricultural intervention on soils

Soil microbial communities are responsible for a wide range of ecological processes and have an important economic impact in agriculture. Determining the metabolic processes performed by microbial communities is crucial for understanding and managing ecosystem properties. Metagenomic approaches allow the elucidation of the main metabolic processes that determine the performance of microbial communities under different environmental conditions and perturbations. Here we present the first compartmentalized metabolic reconstruction at a metagenomics scale of a microbial ecosystem. This systematic approach conceives a meta-organism without boundaries between individual organisms and allows the in silico evaluation of the effect of agricultural intervention on soils at a metagenomics level. To characterize the microbial ecosystems, topological properties, taxonomic and metabolic profiles, as well as a Flux Balance Analysis (FBA) were considered. Furthermore, topological and optimization algorithms were implemented to carry out the curation of the models, to ensure the continuity of the fluxes between the metabolic pathways, and to confirm the metabolite exchange between subcellular compartments. The proposed models provide specific information about ecosystems that are generally overlooked in non-compartmentalized or non-curated networks, like the influence of transport reactions in the metabolic processes, especially the important effect on mitochondrial processes, as well as provide more accurate results of the fluxes used to optimize the metabolic processes within the microbial community.

Analysing Algorithms and Data Sources for the Tissue-Specific Reconstruction of Liver Healthy and Cancer Cells

Genome-Scale Metabolic Models (GSMMs), mathematical representations of the cell metabolism in different organisms including humans, are resourceful tools to simulate metabolic phenotypes and understand associated diseases, such as obesity, diabetes and cancer. In the last years, different algorithms have been developed to generate tissue-specific metabolic models that simulate different phenotypes for distinct cell types. Hepatocyte cells are one of the main sites of metabolic conversions, mainly due to their diverse physiological functions. Most of the liver's tissue is formed by hepatocytes, being one of the largest and most important organs regarding its biological functions. Hepatocellular carcinoma is, also, one of the most important human cancers with high mortality rates. In this study, we will analyze four different algorithms (MBA, mCADRE, tINIT and FASTCORE) for tissue-specific model reconstruction, based on a template model and two types of data sources: transcriptomics and proteomics. These methods will be applied to the reconstruction of metabolic models for hepatocyte cells and HepG2 cancer cell line. The models will be analyzed and compared under different perspectives, emphasizing their functional analysis considering a set of metabolic liver tasks. The results show that there is no "ideal" algorithm. However, with the current analysis, we were able to retrieve knowledge about the metabolism of the liver.

Reconstruction of genome-scale human metabolic models using omics data

The impact of genome-scale human metabolic models on human systems biology and medical sciences is becoming greater, thanks to increasing volumes of model building platforms and publicly available omics data. The genome-scale human metabolic models started with Recon 1 in 2007, and have since been used to describe metabolic phenotypes of healthy and diseased human tissues and cells, and to predict therapeutic targets. Here we review recent trends in genome-scale human metabolic modeling, including various generic and tissue/cell type-specific human metabolic models developed to date, and methods, databases and platforms used to construct them. For generic human metabolic models, we pay attention to Recon 2 and HMR 2.0 with emphasis on data sources used to construct them. Draft and high-quality tissue/cell type-specific human metabolic models have been generated using these generic human metabolic models. Integration of tissue/cell type-specific omics data with the generic human metabolic models is the key step, and we discuss omics data and their integration methods to achieve this task. The initial version of the tissue/cell type-specific human metabolic models can further be computationally refined through gap filling, reaction directionality assignment and the subcellular localization of metabolic reactions. We review relevant tools for this model refinement procedure as well. Finally, we suggest the direction of further studies on reconstructing an improved human metabolic model.

Zooming-in on cancer metabolic rewiring with tissue specific constraint-based models

The metabolic rearrangements occurring in cancer cells can be effectively investigated with a Systems Biology approach supported by metabolic network modeling. We here present tissue-specific constraint-based core models for three different types of tumors (liver, breast and lung) that serve this purpose. The core models were extracted and manually curated from the corresponding genome-scale metabolic models in the Human Metabolic Atlas database with a focus on the pathways that are known to play a key role in cancer growth and proliferation. Along similar lines, we also reconstructed a core model from the original general human metabolic network to be used as a reference model. A comparative Flux Balance Analysis between the reference and the cancer models highlighted both a clear distinction between the two conditions and a heterogeneity within the three different cancer types in terms of metabolic flux distribution. These results emphasize the need for modeling approaches able to keep up with this tumoral heterogeneity in order to identify more suitable drug targets and develop effective treatments. According to this perspective, we identified key points able to reverse the tumoral phenotype toward the reference one or vice-versa.

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.

Genome-scale metabolic modelling of hepatocytes reveals serine deficiency in patients with non-alcoholic fatty liver disease

Several liver disorders result from perturbations in the metabolism of hepatocytes, and their underlying mechanisms can be outlined through the use of genome-scale metabolic models (GEMs). Here we reconstruct a consensus GEM for hepatocytes, which we call iHepatocytes2322, that extends previous models by including an extensive description of lipid metabolism. We build iHepatocytes2322 using Human Metabolic Reaction 2.0 database and proteomics data in Human Protein Atlas, which experimentally validates the incorporated reactions. The reconstruction process enables improved annotation of the proteomics data using the network centric view of iHepatocytes2322. We then use iHepatocytes2322 to analyse transcriptomics data obtained from patients with non-alcoholic fatty liver disease. We show that blood concentrations of chondroitin and heparan sulphates are suitable for diagnosing non-alcoholic steatohepatitis and for the staging of non-alcoholic fatty liver disease. Furthermore, we observe serine deficiency in patients with NASH and identify PSPH, SHMT1 and BCAT1 as potential therapeutic targets for the treatment of non-alcoholic steatohepatitis.

Reconstruction and analysis of a genome-scale metabolic model for Eriocheir sinensis eyestalks

The eyestalk of Eriocheir sinensis has significant biological functions with many nerve peptide hormones expressed in the X-organ which exists in the eyestalk. A metabolic network model is an effective tool for the systematic study of E. sinensis eyestalks. In this work, we reconstructed a metabolic network model for E. sinensis eyestalks based on transcriptome sequencing. The model contains 1304 reactions, 1381 unigenes and 1243 metabolites distributing in 98 pathways. The reconstructed metabolic network model was used for the functional module and block metabolite analysis of eyestalks, which reveals that the function of the eyestalk network agrees with its function as the centre of the E. sinensis endocrine system. The difference expression analysis of reactions in the model indicates that the eyestalk mainly functions in the regulation of amino acids, carbohydrate and nucleotide metabolism.

Dissecting Germ Cell Metabolism through Network Modeling

Metabolic pathways are increasingly postulated to be vital in programming cell fate, including stemness, differentiation, proliferation, and apoptosis. The commitment to meiosis is a critical fate decision for mammalian germ cells, and requires a metabolic derivative of vitamin A, retinoic acid (RA). Recent evidence showed that a pulse of RA is generated in the testis of male mice thereby triggering meiotic commitment. However, enzymes and reactions that regulate this RA pulse have yet to be identified. We developed a mouse germ cell-specific metabolic network with a curated vitamin A pathway. Using this network, we implemented flux balance analysis throughout the initial wave of spermatogenesis to elucidate important reactions and enzymes for the generation and degradation of RA. Our results indicate that primary RA sources in the germ cell include RA import from the extracellular region, release of RA from binding proteins, and metabolism of retinal to RA. Further, in silico knockouts of genes and reactions in the vitamin A pathway predict that deletion of Lipe, hormone-sensitive lipase, disrupts the RA pulse thereby causing spermatogenic defects. Examination of other metabolic pathways reveals that the citric acid cycle is the most active pathway. In addition, we discover that fatty acid synthesis/oxidation are the primary energy sources in the germ cell. In summary, this study predicts enzymes, reactions, and pathways important for germ cell commitment to meiosis. These findings enhance our understanding of the metabolic control of germ cell differentiation and will help guide future experiments to improve reproductive health.

Using gene expression signatures to identify novel treatment strategies in gulf war illness

Background

Gulf War Illness (GWI) is a complex multi-symptom disorder that affects up to one in three veterans of this 1991 conflict and for which no effective treatment has been found. Discovering novel treatment strategies for such a complex chronic illness is extremely expensive, carries a high probability of failure and a lengthy cycle time. Repurposing Food and Drug Administration approved drugs offers a cost-effective solution with a significantly abbreviated timeline.

Methods

Here, we explore drug re-purposing opportunities in GWI by combining systems biology and bioinformatics techniques with pharmacogenomic information to find overlapping elements in gene expression linking GWI to successfully treated diseases. Gene modules were defined based on cellular function and their activation estimated from the differential expression of each module’s constituent genes. These gene modules were then cross-referenced with drug atlas and pharmacogenomic databases to identify agents currently used successfully for treatment in other diseases. To explore the clinical use of these drugs in illnesses similar to GWI we compared gene expression patterns in modules that were significantly expressed in GWI with expression patterns in those same modules in other illnesses.

Results

We found 19 functional modules with significantly altered gene expression patterns in GWI. Within these modules, 45 genes were documented drug targets. Illnesses with highly correlated gene expression patterns overlapping considerably with GWI were found in 18 of the disease conditions studied. Brain, muscular and autoimmune disorders composed the bulk of these.

Conclusion

Of the associated drugs, immunosuppressants currently used in treating rheumatoid arthritis, and hormone based therapies were identified as the best available candidates for treating GWI symptoms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-015-0111-3) contains supplementary material, which is available to authorized users.

Alterations in cancer cell metabolism: the Warburg effect and metabolic adaptation

The Warburg effect means higher glucose uptake of cancer cells compared to normal tissues, whereas a smaller fraction of this glucose is employed for oxidative phosphorylation. With the advent of high throughput technologies and computational systems biology, cancer cell metabolism has been reinvestigated over the last decades toward identifying various events underlying "how" and "why" a cancer cell employs aerobic glycolysis. Significant progress has been shaped to revise the Warburg effect. In this study, we have integrated the gene expression of 13 different cancer cells with the genome-scale metabolic network of human (Recon1) based on the E-Flux method, and analyzed them based on constraint-based modeling. Results show that regardless of significant up- and down-regulated metabolic genes, the distribution of metabolic changes is similar in different cancer types. These findings support the theory that the Warburg effect is a consequence of metabolic adaptation in cancer cells.

Metabolomics and diabetes: analytical and computational approaches

Diabetes is characterized by altered metabolism of key molecules and regulatory pathways. The phenotypic expression of diabetes and associated complications encompasses complex interactions between genetic, environmental, and tissue-specific factors that require an integrated understanding of perturbations in the network of genes, proteins, and metabolites. Metabolomics attempts to systematically identify and quantitate small molecule metabolites from biological systems. The recent rapid development of a variety of analytical platforms based on mass spectrometry and nuclear magnetic resonance have enabled identification of complex metabolic phenotypes. Continued development of bioinformatics and analytical strategies has facilitated the discovery of causal links in understanding the pathophysiology of diabetes and its complications. Here, we summarize the metabolomics workflow, including analytical, statistical, and computational tools, highlight recent applications of metabolomics in diabetes research, and discuss the challenges in the field.

Comparison of human cell signaling pathway databases--evolution, drawbacks and challenges

Elucidating the complexities of cell signaling pathways is of immense importance to gain understanding about various biological phenomenon, such as dynamics of gene/protein expression regulation, cell fate determination, embryogenesis and disease progression. The successful completion of human genome project has also helped experimental and theoretical biologists to analyze various important pathways. To advance this study, during the past two decades, systematic collections of pathway data from experimental studies have been compiled and distributed freely by several databases, which also integrate various computational tools for further analysis. Despite significant advancements, there exist several drawbacks and challenges, such as pathway data heterogeneity, annotation, regular update and automated image reconstructions, which motivated us to perform a thorough review on popular and actively functioning 24 cell signaling databases. Based on two major characteristics, pathway information and technical details, freely accessible data from commercial and academic databases are examined to understand their evolution and enrichment. This review not only helps to identify some novel and useful features, which are not yet included in any of the databases but also highlights their current limitations and subsequently propose the reasonable solutions for future database development, which could be useful to the whole scientific community.

Computational strategies for a system-level understanding of metabolism

Cell metabolism is the biochemical machinery that provides energy and building blocks to sustain life. Understanding its fine regulation is of pivotal relevance in several fields, from metabolic engineering applications to the treatment of metabolic disorders and cancer. Sophisticated computational approaches are needed to unravel the complexity of metabolism. To this aim, a plethora of methods have been developed, yet it is generally hard to identify which computational strategy is most suited for the investigation of a specific aspect of metabolism. This review provides an up-to-date description of the computational methods available for the analysis of metabolic pathways, discussing their main advantages and drawbacks. In particular, attention is devoted to the identification of the appropriate scale and level of accuracy in the reconstruction of metabolic networks, and to the inference of model structure and parameters, especially when dealing with a shortage of experimental measurements. The choice of the proper computational methods to derive in silico data is then addressed, including topological analyses, constraint-based modeling and simulation of the system dynamics. A description of some computational approaches to gain new biological knowledge or to formulate hypotheses is finally provided.

Systems biology approaches for studying the pathogenesis of non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) is a progressive disease of increasing public health concern. In western populations the disease has an estimated prevalence of 20%-40%, rising to 70%-90% in obese and type II diabetic individuals. Simplistically, NAFLD is the macroscopic accumulation of lipid in the liver, and is viewed as the hepatic manifestation of the metabolic syndrome. However, the molecular mechanisms mediating both the initial development of steatosis and its progression through non-alcoholic steatohepatitis to debilitating and potentially fatal fibrosis and cirrhosis are only partially understood. Despite increased research in this field, the development of non-invasive clinical diagnostic tools and the discovery of novel therapeutic targets has been frustratingly slow. We note that, to date, NAFLD research has been dominated by in vivo experiments in animal models and human clinical studies. Systems biology tools and novel computational simulation techniques allow the study of large-scale metabolic networks and the impact of their dysregulation on health. Here we review current systems biology tools and discuss the benefits to their application to the study of NAFLD. We propose that a systems approach utilising novel in silico modelling and simulation techniques is key to a more comprehensive, better targeted NAFLD research strategy. Such an approach will accelerate the progress of research and vital translation into clinic.

Anomalies in network bridges involved in bile Acid metabolism predict outcomes of colorectal cancer patients

Biomarkers prognostic for colorectal cancer (CRC) would be highly desirable in clinical practice. Proteins that regulate bile acid (BA) homeostasis, by linking metabolic sensors and metabolic enzymes, also called bridge proteins, may be reliable prognostic biomarkers for CRC. Based on a devised metric, "bridgeness," we identified bridge proteins involved in the regulation of BA homeostasis and identified their prognostic potentials. The expression patterns of these bridge proteins could distinguish between normal and diseased tissues, suggesting that these proteins are associated with CRC pathogenesis. Using a supervised classification system, we found that these bridge proteins were reproducibly prognostic, with high prognostic ability compared to other known markers.

Prioritizing candidate disease metabolites based on global functional relationships between metabolites in the context of metabolic pathways

Identification of key metabolites for complex diseases is a challenging task in today's medicine and biology. A special disease is usually caused by the alteration of a series of functional related metabolites having a global influence on the metabolic network. Moreover, the metabolites in the same metabolic pathway are often associated with the same or similar disease. Based on these functional relationships between metabolites in the context of metabolic pathways, we here presented a pathway-based random walk method called PROFANCY for prioritization of candidate disease metabolites. Our strategy not only takes advantage of the global functional relationships between metabolites but also sufficiently exploits the functionally modular nature of metabolic networks. Our approach proved successful in prioritizing known metabolites for 71 diseases with an AUC value of 0.895. We also assessed the performance of PROFANCY on 16 disease classes and found that 4 classes achieved an AUC value over 0.95. To investigate the robustness of the PROFANCY, we repeated all the analyses in two metabolic networks and obtained similar results. Then we applied our approach to Alzheimer's disease (AD) and found that a top ranked candidate was potentially related to AD but had not been reported previously. Furthermore, our method was applicable to prioritize the metabolites from metabolomic profiles of prostate cancer. The PROFANCY could identify prostate cancer related-metabolites that are supported by literatures but not considered to be significantly differential by traditional differential analysis. We also developed a freely accessible web-based and R-based tool at http://bioinfo.hrbmu.edu.cn/PROFANCY.

An atlas of genetic influences on human blood metabolites

Genome-wide association scans with high-throughput metabolic profiling provide unprecedented insights into how genetic variation influences metabolism and complex disease. Here we report the most comprehensive exploration of genetic loci influencing human metabolism thus far, comprising 7,824 adult individuals from 2 European population studies. We report genome-wide significant associations at 145 metabolic loci and their biochemical connectivity with more than 400 metabolites in human blood. We extensively characterize the resulting in vivo blueprint of metabolism in human blood by integrating it with information on gene expression, heritability and overlap with known loci for complex disorders, inborn errors of metabolism and pharmacological targets. We further developed a database and web-based resources for data mining and results visualization. Our findings provide new insights into the role of inherited variation in blood metabolic diversity and identify potential new opportunities for drug development and for understanding disease.

MetDisease--connecting metabolites to diseases via literature

MOTIVATION

In recent years, metabolomics has emerged as an approach to perform large-scale characterization of small molecules in biological systems. Metabolomics posed a number of bioinformatics challenges associated in data analysis and interpretation. Genome-based metabolic reconstructions have established a powerful framework for connecting metabolites to genes through metabolic reactions and enzymes that catalyze them. Pathway databases and bioinformatics tools that use this framework have proven to be useful for annotating experimental metabolomics data. This framework can be used to infer connections between metabolites and diseases through annotated disease genes. However, only about half of experimentally detected metabolites can be mapped to canonical metabolic pathways. We present a new Cytoscape 3 plug-in, MetDisease, which uses an alternative approach to link metabolites to disease information. MetDisease uses Medical Subject Headings (MeSH) disease terms mapped to PubChem compounds through literature to annotate compound networks.

AVAILABILITY AND IMPLEMENTATION

MetDisease can be downloaded from http://apps.cytoscape.org/apps/metdisease or installed via the Cytoscape app manager. Further information about MetDisease can be found at http://metdisease.ncibi.org

CONTACT

akarnovs@med.umich.edu

SUPPLEMENTARY INFORMATION

Supplementary Data are available at Bioinformatics online.

Metabolic cancer biology: structural-based analysis of cancer as a metabolic disease, new sights and opportunities for disease treatment

The cancer cell metabolism or the Warburg effect discovery goes back to 1924 when, for the first time Otto Warburg observed, in contrast to the normal cells, cancer cells have different metabolism. With the initiation of high throughput technologies and computational systems biology, cancer cell metabolism renaissances and many attempts were performed to revise the Warburg effect. The development of experimental and analytical tools which generate high-throughput biological data including lots of information could lead to application of computational models in biological discovery and clinical medicine especially for cancer. Due to the recent availability of tissue-specific reconstructed models, new opportunities in studying metabolic alteration in various kinds of cancers open up. Structural approaches at genome-scale levels seem to be suitable for developing diagnostic and prognostic molecular signatures, as well as in identifying new drug targets. In this review, we have considered these recent advances in structural-based analysis of cancer as a metabolic disease view. Two different structural approaches have been described here: topological and constraint-based methods. The ultimate goal of this type of systems analysis is not only the discovery of novel drug targets but also the development of new systems-based therapy strategies.

Metabolomics and systems pharmacology: why and how to model the human metabolic network for drug discovery

Metabolism represents the 'sharp end' of systems biology, because changes in metabolite concentrations are necessarily amplified relative to changes in the transcriptome, proteome and enzyme activities, which can be modulated by drugs. To understand such behaviour, we therefore need (and increasingly have) reliable consensus (community) models of the human metabolic network that include the important transporters. Small molecule 'drug' transporters are in fact metabolite transporters, because drugs bear structural similarities to metabolites known from the network reconstructions and from measurements of the metabolome. Recon2 represents the present state-of-the-art human metabolic network reconstruction; it can predict inter alia: (i) the effects of inborn errors of metabolism; (ii) which metabolites are exometabolites, and (iii) how metabolism varies between tissues and cellular compartments. However, even these qualitative network models are not yet complete. As our understanding improves so do we recognise more clearly the need for a systems (poly)pharmacology.

Gene duplication and phenotypic changes in the evolution of mammalian metabolic networks

Metabolic networks attempt to describe the complete suite of biochemical reactions available to an organism. One notable feature of these networks in mammals is the large number of distinct proteins that catalyze the same reaction. While the existence of these isoenzymes has long been known, their evolutionary significance is still unclear. Using a phylogenetically-aware comparative genomics approach, we infer enzyme orthology networks for sixteen mammals as well as for their common ancestors. We find that the pattern of isoenzymes copy-number alterations (CNAs) in these networks is suggestive of natural selection acting on the retention of certain gene duplications. When further analyzing these data with a machine-learning approach, we found that that the pattern of CNAs is also predictive of several important phenotypic traits, including milk composition and geographic range. Integrating tools from network analyses, phylogenetics and comparative genomics both allows the prediction of phenotypes from genetic data and represents a means of unifying distinct biological disciplines.

Systems medicine: a new approach to clinical practice

Most respiratory diseases are considered complex diseases as their susceptibility and outcomes are determined by the interaction between host-dependent factors (genetic factors, comorbidities, etc.) and environmental factors (exposure to microorganisms or allergens, treatments received, etc.) The reductionist approach in the study of diseases has been of fundamental importance for the understanding of the different components of a system. Systems biology or systems medicine is a complementary approach aimed at analyzing the interactions between the different components within one organizational level (genome, transcriptome, proteome), and then between the different levels. Systems medicine is currently used for the interpretation and understanding of the pathogenesis and pathophysiology of different diseases, biomarker discovery, design of innovative therapeutic targets, and the drawing up of computational models for different biological processes. In this review we discuss the most relevant concepts of the theory underlying systems medicine, as well as its applications in the various biological processes in humans.

Fast reconstruction of compact context-specific metabolic network models

Systemic approaches to the study of a biological cell or tissue rely increasingly on the use of context-specific metabolic network models. The reconstruction of such a model from high-throughput data can routinely involve large numbers of tests under different conditions and extensive parameter tuning, which calls for fast algorithms. We present fastcore, a generic algorithm for reconstructing context-specific metabolic network models from global genome-wide metabolic network models such as Recon X. fastcore takes as input a core set of reactions that are known to be active in the context of interest (e.g., cell or tissue), and it searches for a flux consistent subnetwork of the global network that contains all reactions from the core set and a minimal set of additional reactions. Our key observation is that a minimal consistent reconstruction can be defined via a set of sparse modes of the global network, and fastcore iteratively computes such a set via a series of linear programs. Experiments on liver data demonstrate speedups of several orders of magnitude, and significantly more compact reconstructions, over a rival method. Given its simplicity and its excellent performance, fastcore can form the backbone of many future metabolic network reconstruction algorithms.

Network pharmacology strategies toward multi-target anticancer therapies: from computational models to experimental design principles

Polypharmacology has emerged as novel means in drug discovery for improving treatment response in clinical use. However, to really capitalize on the polypharmacological effects of drugs, there is a critical need to better model and understand how the complex interactions between drugs and their cellular targets contribute to drug efficacy and possible side effects. Network graphs provide a convenient modeling framework for dealing with the fact that most drugs act on cellular systems through targeting multiple proteins both through on-target and off-target binding. Network pharmacology models aim at addressing questions such as how and where in the disease network should one target to inhibit disease phenotypes, such as cancer growth, ideally leading to therapies that are less vulnerable to drug resistance and side effects by means of attacking the disease network at the systems level through synergistic and synthetic lethal interactions. Since the exponentially increasing number of potential drug target combinations makes pure experimental approach quickly unfeasible, this review depicts a number of computational models and algorithms that can effectively reduce the search space for determining the most promising combinations for experimental evaluation. Such computational-experimental strategies are geared toward realizing the full potential of multi-target treatments in different disease phenotypes. Our specific focus is on system-level network approaches to polypharmacology designs in anticancer drug discovery, where we give representative examples of how network-centric modeling may offer systematic strategies toward better understanding and even predicting the phenotypic responses to multi-target therapies.

Investigating host-pathogen behavior and their interaction using genome-scale metabolic network models

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.

Controllability in cancer metabolic networks according to drug targets as driver nodes

Networks are employed to represent many nonlinear complex systems in the real world. The topological aspects and relationships between the structure and function of biological networks have been widely studied in the past few decades. However dynamic and control features of complex networks have not been widely researched, in comparison to topological network features. In this study, we explore the relationship between network controllability, topological parameters, and network medicine (metabolic drug targets). Considering the assumption that targets of approved anticancer metabolic drugs are driver nodes (which control cancer metabolic networks), we have applied topological analysis to genome-scale metabolic models of 15 normal and corresponding cancer cell types. The results show that besides primary network parameters, more complex network metrics such as motifs and clusters may also be appropriate for controlling the systems providing the controllability relationship between topological parameters and drug targets. Consequently, this study reveals the possibilities of following a set of driver nodes in network clusters instead of considering them individually according to their centralities. This outcome suggests considering distributed control systems instead of nodal control for cancer metabolic networks, leading to a new strategy in the field of network medicine.

Software platforms to facilitate reconstructing genome-scale metabolic networks

System-level analyses of microbial metabolism are facilitated by genome-scale reconstructions of microbial biochemical networks. A reconstruction provides a structured representation of the biochemical transformations occurring within an organism, as well as the genes necessary to carry out these transformations, as determined by the annotated genome sequence and experimental data. Network reconstructions also serve as platforms for constraint-based computational techniques, which facilitate biological studies in a variety of applications, including evaluation of network properties, metabolic engineering and drug discovery. Bottom-up metabolic network reconstructions have been developed for dozens of organisms, but until recently, the pace of reconstruction has failed to keep up with advances in genome sequencing. To address this problem, a number of software platforms have been developed to automate parts of the reconstruction process, thereby alleviating much of the manual effort previously required. Here, we review four such platforms in the context of established guidelines for network reconstruction. While many steps of the reconstruction process have been successfully automated, some manual evaluation of the results is still required to ensure a high-quality reconstruction. Widespread adoption of these platforms by the scientific community is underway and will be further enabled by exchangeable formats across platforms.

A novel untargeted metabolomics correlation-based network analysis incorporating human metabolic reconstructions

BACKGROUND

Metabolomics has become increasingly popular in the study of disease phenotypes and molecular pathophysiology. One branch of metabolomics that encompasses the high-throughput screening of cellular metabolism is metabolic profiling. In the present study, the metabolic profiles of different tumour cells from colorectal carcinoma and breast adenocarcinoma were exposed to hypoxic and normoxic conditions and these have been compared to reveal the potential metabolic effects of hypoxia on the biochemistry of the tumour cells; this may contribute to their survival in oxygen compromised environments. In an attempt to analyse the complex interactions between metabolites beyond routine univariate and multivariate data analysis methods, correlation analysis has been integrated with a human metabolic reconstruction to reveal connections between pathways that are associated with normoxic or hypoxic oxygen environments.

RESULTS

Correlation analysis has revealed statistically significant connections between metabolites, where differences in correlations between cells exposed to different oxygen levels have been highlighted as markers of hypoxic metabolism in cancer. Network mapping onto reconstructed human metabolic models is a novel addition to correlation analysis. Correlated metabolites have been mapped onto the Edinburgh human metabolic network (EHMN) with the aim of interlinking metabolites found to be regulated in a similar fashion in response to oxygen. This revealed novel pathways within the metabolic network that may be key to tumour cell survival at low oxygen. Results show that the metabolic responses to lowering oxygen availability can be conserved or specific to a particular cell line. Network-based correlation analysis identified conserved metabolites including malate, pyruvate, 2-oxoglutarate, glutamate and fructose-6-phosphate. In this way, this method has revealed metabolites not previously linked, or less well recognised, with respect to hypoxia before. Lactate fermentation is one of the key themes discussed in the field of hypoxia; however, malate, pyruvate, 2-oxoglutarate, glutamate and fructose-6-phosphate, which are connected by a single pathway, may provide a more significant marker of hypoxia in cancer.

CONCLUSIONS

Metabolic networks generated for each cell line were compared to identify conserved metabolite pathway responses to low oxygen environments. Furthermore, we believe this methodology will have general application within metabolomics.

Consensus and conflict cards for metabolic pathway databases

Background

The metabolic network of H. sapiens and many other organisms is described in multiple pathway databases. The level of agreement between these descriptions, however, has proven to be low. We can use these different descriptions to our advantage by identifying conflicting information and combining their knowledge into a single, more accurate, and more complete description. This task is, however, far from trivial.

Results

We introduce the concept of Consensus and Conflict Cards (C₂Cards) to provide concise overviews of what the databases do or do not agree on. Each card is centered at a single gene, EC number or reaction. These three complementary perspectives make it possible to distinguish disagreements on the underlying biology of a metabolic process from differences that can be explained by different decisions on how and in what detail to represent knowledge. As a proof-of-concept, we implemented C₂Cards^Human, as a web application http://www.molgenis.org/c2cards, covering five human pathway databases.

Conclusions

C₂Cards can contribute to ongoing reconciliation efforts by simplifying the identification of consensus and conflicts between pathway databases and lowering the threshold for experts to contribute. Several case studies illustrate the potential of the C₂Cards in identifying disagreements on the underlying biology of a metabolic process. The overviews may also point out controversial biological knowledge that should be subject of further research. Finally, the examples provided emphasize the importance of manual curation and the need for a broad community involvement.

Novel insights into obesity and diabetes through genome-scale metabolic modeling

The growing prevalence of metabolic diseases, such as obesity and diabetes, are putting a high strain on global healthcare systems as well as increasing the demand for efficient treatment strategies. More than 360 million people worldwide are suffering from type 2 diabetes (T2D) and, with the current trends, the projection is that 10% of the global adult population will be affected by 2030. In light of the systemic properties of metabolic diseases as well as the interconnected nature of metabolism, it is necessary to begin taking a holistic approach to study these diseases. Human genome-scale metabolic models (GEMs) are topological and mathematical representations of cell metabolism and have proven to be valuable tools in the area of systems biology. Successful applications of GEMs include the process of gaining further biological and mechanistic understanding of diseases, finding potential biomarkers, and identifying new drug targets. This review will focus on the modeling of human metabolism in the field of obesity and diabetes, showing its vast range of applications of clinical importance as well as point out future challenges.