Computational reconstruction of tissue-specific metabolic models: application to human liver metabolism

Metabolic network modeling approaches for investigating the "hungry cancer"

Metabolism is the functional phenotype of a cell, at a given condition, resulting from an intricate interplay of various regulatory processes. The study of these dynamic metabolic processes and their capabilities help to identify the fundamental properties of living systems. Metabolic deregulation is an emerging hallmark of cancer cells. This deregulation results in rewiring of the metabolic circuitry conferring an exploitative metabolic advantage for the tumor cells which leads to a distinct benefit in survival and lays the basis for unbound progression. Metabolism can be considered as a thermodynamic open-system in which source substrates of high value are being processed through a well established interconnected biochemical conversion system, strictly obeying physiochemical principles, generating useful intermediates and finally resulting in the release of byproducts. Based on this basic principle of an input-output balance, various models have been developed to interrogate metabolism elucidating its underlying functional properties. However, only a few modeling approaches have proved computationally feasible in elucidating the metabolic nature of cancer at a systems level. Besides this, statistical approaches have been set up to identify biochemical pathways being more relevant for specific types of tumor cells. In this review, we are briefly introducing the basic statistical approaches followed by the major modeling concepts. We have put an emphasis on the methods and their applications that have been used to a greater extent in understanding the metabolic remodeling of cancer.

Metabolic states following accumulation of intracellular aggregates: implications for neurodegenerative diseases

The formation of intracellular aggregates is a common etiology of several neurodegenerative diseases. Mitochondrial defects and oxidative stress has been pointed as the major mechanistic links between the accumulation of intracellular aggregates and cell death. In this work we propose a “metabolic cell death by overcrowding” as an alternative hypothesis. Using a model of neuron metabolism, we predict that as the concentration of protein aggregates increases the neurons transit through three different metabolic phases. The first phase (0–6 mM) corresponds with the normal neuron state, where the neuronal activity is sustained by the oxidative phosphorylation of lactate. The second phase (6–8.6 mM) is characterized by a mixed utilization of lactate and glucose as energy substrates and a switch from ammonia uptake to ammonia release by neurons. In the third phase (8.6–9.3 mM) neurons are predicted to support their energy demands from glycolysis and an alternative pathway for energy generation, involving reactions from serine synthesis, one carbon metabolism and the glycine cleavage system. The model also predicts a decrease in the maximum neuronal capacity for energy generation with increasing the concentration of protein aggregates. Ultimately this maximum capacity becomes zero when the protein aggregates reach a concentration of about 9.3 mM, predicting the cessation of neuronal activity.

Genome-scale modeling of human metabolism - a systems biology approach

Altered metabolism is linked to the appearance of various human diseases and a better understanding of disease-associated metabolic changes may lead to the identification of novel prognostic biomarkers and the development of new therapies. Genome-scale metabolic models (GEMs) have been employed for studying human metabolism in a systematic manner, as well as for understanding complex human diseases. In the past decade, such metabolic models - one of the fundamental aspects of systems biology - have started contributing to the understanding of the mechanistic relationship between genotype and phenotype. In this review, we focus on the construction of the Human Metabolic Reaction database, the generation of healthy cell type- and cancer-specific GEMs using different procedures, and the potential applications of these developments in the study of human metabolism and in the identification of metabolic changes associated with various disorders. We further examine how in silico genome-scale reconstructions can be employed to simulate metabolic flux distributions and how high-throughput omics data can be analyzed in a context-dependent fashion. Insights yielded from this mechanistic modeling approach can be used for identifying new therapeutic agents and drug targets as well as for the discovery of novel biomarkers. Finally, recent advancements in genome-scale modeling and the future challenge of developing a model of whole-body metabolism are presented. The emergent contribution of GEMs to personalized and translational medicine is also discussed.

Flux-coupled genes and their use in metabolic flux analysis

As large volumes of omics data have become available, systems biology is playing increasingly important roles in elucidating new biological phenomena, especially through genome-scale metabolic network modeling and simulation. Much effort has been exerted on integrating omics data with metabolic flux simulation, but further development is necessary for more accurate flux estimation. To move one step forward, we adopted the concept of flux-coupled genes (FCGs), which show that their expression transition patterns upon perturbations are correlated with their corresponding flux values, as additional constraints in metabolic flux analysis. It was found that gnd, pfkB, rpe, sdhB, sdhD, sucA, and zwf genes, mostly associated with pentose phosphate pathway and TCA cycle, were the most consistent FCGs in Escherichia coli based on its transcriptome and (13) C-flux data obtained from the chemostat cultivation at five different dilution rates. Consequently, constraints-based flux analyses with FCGs as additional constraints were conducted for the seven single-gene knockout mutants, compared with those obtained without using FCGs. This strategy of constraining the metabolic flux analysis with FCGs is expected to be useful due to the relative ease in obtaining transcriptional information in the functional genomics era.

Inferring metabolic states in uncharacterized environments using gene-expression measurements

The large size of metabolic networks entails an overwhelming multiplicity in the possible steady-state flux distributions that are compatible with stoichiometric constraints. This space of possibilities is largest in the frequent situation where the nutrients available to the cells are unknown. These two factors: network size and lack of knowledge of nutrient availability, challenge the identification of the actual metabolic state of living cells among the myriad possibilities. Here we address this challenge by developing a method that integrates gene-expression measurements with genome-scale models of metabolism as a means of inferring metabolic states. Our method explores the space of alternative flux distributions that maximize the agreement between gene expression and metabolic fluxes, and thereby identifies reactions that are likely to be active in the culture from which the gene-expression measurements were taken. These active reactions are used to build environment-specific metabolic models and to predict actual metabolic states. We applied our method to model the metabolic states of Saccharomyces cerevisiae growing in rich media supplemented with either glucose or ethanol as the main energy source. The resulting models comprise about 50% of the reactions in the original model, and predict environment-specific essential genes with high sensitivity. By minimizing the sum of fluxes while forcing our predicted active reactions to carry flux, we predicted the metabolic states of these yeast cultures that are in large agreement with what is known about yeast physiology. Most notably, our method predicts the Crabtree effect in yeast cells growing in excess glucose, a long-known phenomenon that could not have been predicted by traditional constraint-based modeling approaches. Our method is of immediate practical relevance for medical and industrial applications, such as the identification of novel drug targets, and the development of biotechnological processes that use complex, largely uncharacterized media, such as biofuel production.

Metabolic fluxes are steady-state rates of metabolite interconversion within living cells. They determine the rates of growth and product formation, and are of biotechnological and medical importance. An important and pressing question is how to identify the actual distribution of fluxes in living cells among the manifold possibilities that complex metabolic networks allow. One way to address this question is to constrain the space of possibilities using gene-expression measurements. Here we present a method that uses gene-expression measurements to infer the metabolic state of cells growing in uncharacterized environments. We applied this method to model the metabolism of Saccharomyces cerevisiae grown with glucose or ethanol as main energy source. Our modeling enables the prediction of genes that are essential for growth in either environment. We also show that our method predicts aspects of the energy metabolism of these cultures that are in large agreement with what is known about yeast physiology. Our method is of direct practical importance in the fields of biotechnology and medicine, such as in vivo drug target identification, where nutrient conditions are largely unknown.

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/.

Hepatocellular carcinoma: a systems biology perspective

Hepatocellular carcinomas (HCCs) have different etiology and heterogenic genomic alterations lead to high complexity. The molecular features of HCC have largely been studied by gene expression and proteome profiling focusing on the correlations between the expression of specific markers and clinical data. Integration of the increasing amounts of data in databases has facilitated the link of genomic and proteomic profiles of HCC to disease state and clinical outcome. Despite the current knowledge, specific molecular markers remain to be identified and new strategies are required to establish novel-targeted therapies. In the last years, mathematical models reconstructing gene and protein networks based on experimental data of HCC have been developed providing powerful tools to predict candidate interactions and potential targets for therapy. Furthermore, the combination of dynamic and logical mathematical models with quantitative data allows detailed mechanistic insights into system properties. To address effects at the organ level, mathematical models reconstructing the three-dimensional organization of liver lobules were developed. In the future, integration of different modeling approaches capturing the effects at the cellular up to the organ level is required to address the complex properties of HCC and to enable the discovery of new targets for HCC prevention or treatment.

p53 promotes the expression of gluconeogenesis-related genes and enhances hepatic glucose production

Background

The p53 tumor suppressor protein is a transcription factor that initiates transcriptional programs aimed at inhibiting carcinogenesis. p53 represses metabolic pathways that support tumor development (such as glycolysis and the pentose phosphate pathway (PPP)) and enhances metabolic pathways that are considered counter-tumorigenic such as fatty acid oxidation.

Findings

In an attempt to comprehensively define metabolic pathways regulated by p53, we performed two consecutive high-throughput analyses in human liver-derived cells with varying p53 statuses. A gene expression microarray screen followed by constraint-based modeling (CBM) predicting metabolic changes imposed by the transcriptomic changes suggested a role for p53 in enhancing gluconeogenesis (de novo synthesis of glucose). Examining glucogenic gene expression revealed a p53-dependent induction of genes involved in both gluconeogenesis (G6PC, PCK2) and in supplying glucogenic precursors (glycerol kinase (GK), aquaporin 3 (AQP3), aquaporin 9 (AQP9) and glutamic-oxaloacetic transaminase 1 (GOT1)). Accordingly, p53 augmented hepatic glucose production (HGP) in both human liver cells and primary mouse hepatocytes.

Conclusions

These findings portray p53 as a novel regulator of glucose production. By facilitating glucose export, p53 may prevent it from being shunted to pro-cancerous pathways such as glycolysis and the PPP. Thus, our findings suggest a metabolic pathway through which p53 may inhibit tumorigenesis.

Application of Top-Down and Bottom-up Systems Approaches in Ruminant Physiology and Metabolism

Systems biology is a computational field that has been used for several years across different scientific areas of biological research to uncover the complex interactions occurring in living organisms. Applications of systems concepts at the mammalian genome level are quite challenging, and new complimentary computational/experimental techniques are being introduced. Most recent work applying modern systems biology techniques has been conducted on bacteria, yeast, mouse, and human genomes. However, these concepts and tools are equally applicable to other species including ruminants (e.g., livestock). In systems biology, both bottom-up and top-down approaches are central to assemble information from all levels of biological pathways that must coordinate physiological processes. A bottom-up approach encompasses draft reconstruction, manual curation, network reconstruction through mathematical methods, and validation of these models through literature analysis (i.e., bibliomics). Whereas top-down approach encompasses metabolic network reconstructions using ‘omics’ data (e.g., transcriptomics, proteomics) generated through DNA microarrays, RNA-Seq or other modern high-throughput genomic techniques using appropriate statistical and bioinformatics methodologies. In this review we focus on top-down approach as a means to improve our knowledge of underlying metabolic processes in ruminants in the context of nutrition. We also explore the usefulness of tissue specific reconstructions (e.g., liver and adipose tissue) in cattle as a means to enhance productive efficiency.

Analysis of omics data with genome-scale models of metabolism

Over the past decade a massive amount of research has been dedicated to generating omics data to gain insight into a variety of biological phenomena, including cancer, obesity, biofuel production, and infection. Although most of these omics data are available publicly, there is a growing concern that much of these data sit in databases without being used or fully analyzed. Statistical inference methods have been widely applied to gain insight into which genes may influence the activities of others in a given omics data set, however, they do not provide information on the underlying mechanisms or whether the interactions are direct or distal. Biochemically, genetically, and genomically consistent knowledge bases are increasingly being used to extract deeper biological knowledge and understanding from these data sets than possible by inferential methods. This improvement is largely due to knowledge bases providing a validated biological context for interpreting the data.

Future of medicine: models in predictive diagnostics and personalized medicine

Molecular medicine is undergoing fundamental changes driving the whole area towards a revolution in modern medicine. The breakthrough was generated the fast-developing technologies in molecular biology since the first draft sequence of the human genome was published. The technological advances enabled the analysis of biological samples from cells and organs to whole organisms in a depth that was not possible before. These technologies are increasingly implemented in the medical and health care system to study diseases and refine diagnostics. As a consequence, the understanding of diseases and the health status of an individual patient is now based on an enormous amount of data that can only be interpreted in the context of the body as a whole. Systems biology as a new field in the life sciences develops new approaches for data integration and interpretation. Systems medicine as a specialized aspect of systems biology combines in an interdisciplinary approach all expertise necessary to decipher the human body in all its complexity. This created new challenges in the area of information and communication technologies to provide the infrastructure and technology needed to cope with the data flood that will accompany the next generation of medicine. The new initiative 'IT Future of Medicine' aims at driving this development even further and integrates not only molecular data (especially genomic information), but also anatomical, physiological, environmental, and lifestyle data in a predictive model approach-the 'virtual patient'-that will allow the clinician or the general practitioner to predict and anticipate the optimal treatment for the individual patient. The application of the virtual patient model will allow truly personalized medicine.

Integration of clinical data with a genome-scale metabolic model of the human adipocyte

We evaluated the presence/absence of proteins encoded by 14 077 genes in adipocytes obtained from different tissue samples using immunohistochemistry. By combining this with previously published adipocyte-specific proteome data, we identified proteins associated with 7340 genes in human adipocytes. This information was used to reconstruct a comprehensive and functional genome-scale metabolic model of adipocyte metabolism. The resulting metabolic model, iAdipocytes1809, enables mechanistic insights into adipocyte metabolism on a genome-wide level, and can serve as a scaffold for integration of omics data to understand the genotype-phenotype relationship in obese subjects. By integrating human transcriptome and fluxome data, we found an increase in the metabolic activity around androsterone, ganglioside GM2 and degradation products of heparan sulfate and keratan sulfate, and a decrease in mitochondrial metabolic activities in obese subjects compared with lean subjects. Our study hereby shows a path to identify new therapeutic targets for treating obesity through combination of high throughput patient data and metabolic modeling.

Novel methodologies in metabolic profiling with a focus on molecular diagnostic applications

The metabolome contains all the biological end points of genomic, transcriptomic and proteomic perturbations, also including the influence of gut microbiota and the environment, giving a direct picture of an organism's ongoing metabolic state. Metabolomics thus has the potential to be an effective tool for early diagnosis of disease, and also to be a predictor of treatment response and survival. In recent years, the development of instrumental systems has enabled more comprehensive coverage of the metabolome. Advances in mass spectrometry and chromatography have particularly improved both the efficiency of nontargeted metabolic profiling as well as the sensitivity and reliability of targeted analyses. Mass spectrometric techniques are also increasingly becoming accepted as a routine diagnostic tool in clinical laboratories. This review summarizes the most recent advances and current challenges in metabolomics, with a focus on mass spectrometric methods utilized in biomarker research, highlighted with selected examples.

Reconstruction of genome-scale metabolic models for 126 human tissues using mCADRE

Background

Human tissues perform diverse metabolic functions. Mapping out these tissue-specific functions in genome-scale models will advance our understanding of the metabolic basis of various physiological and pathological processes. The global knowledgebase of metabolic functions categorized for the human genome (Human Recon 1) coupled with abundant high-throughput data now makes possible the reconstruction of tissue-specific metabolic models. However, the number of available tissue-specific models remains incomplete compared with the large diversity of human tissues.

Results

We developed a method called metabolic Context-specificity Assessed by Deterministic Reaction Evaluation (mCADRE). mCADRE is able to infer a tissue-specific network based on gene expression data and metabolic network topology, along with evaluation of functional capabilities during model building. mCADRE produces models with similar or better functionality and achieves dramatic computational speed up over existing methods. Using our method, we reconstructed draft genome-scale metabolic models for 126 human tissue and cell types. Among these, there are models for 26 tumor tissues along with their normal counterparts, and 30 different brain tissues. We performed pathway-level analyses of this large collection of tissue-specific models and identified the eicosanoid metabolic pathway, especially reactions catalyzing the production of leukotrienes from arachidnoic acid, as potential drug targets that selectively affect tumor tissues.

Conclusions

This large collection of 126 genome-scale draft metabolic models provides a useful resource for studying the metabolic basis for a variety of human diseases across many tissues. The functionality of the resulting models and the fast computational speed of the mCADRE algorithm make it a useful tool to build and update tissue-specific metabolic models.

Genome-level transcription data of Yersinia pestis analyzed with a new metabolic constraint-based approach

Background

Constraint-based computational approaches, such as flux balance analysis (FBA), have proven successful in modeling genome-level metabolic behavior for conditions where a set of simple cellular objectives can be clearly articulated. Recently, the necessity to expand the current range of constraint-based methods to incorporate high-throughput experimental data has been acknowledged by the proposal of several methods. However, these methods have rarely been used to address cellular metabolic responses to some relevant perturbations such as antimicrobial or temperature-induced stress. Here, we present a new method for combining gene-expression data with FBA (GX-FBA) that allows modeling of genome-level metabolic response to a broad range of environmental perturbations within a constraint-based framework. The method uses mRNA expression data to guide hierarchical regulation of cellular metabolism subject to the interconnectivity of the metabolic network.

Results

We applied GX-FBA to a genome-scale model of metabolism in the gram negative bacterium Yersinia pestis and analyzed its metabolic response to (i) variations in temperature known to induce virulence, and (ii) antibiotic stress. Without imposition of any a priori behavioral constraints, our results show strong agreement with reported phenotypes. Our analyses also lead to novel insights into how Y. pestis uses metabolic adjustments to counter different forms of stress.

Conclusions

Comparisons of GX-FBA predicted metabolic states with fluxomic measurements and different reported post-stress phenotypes suggest that mass conservation constraints and network connectivity can be an effective representative of metabolic flux regulation in constraint-based models. We believe that our approach will be of aid in the in silico evaluation of cellular goals under different conditions and can be used for a variety of analyses such as identification of potential drug targets and their action.

Towards multiscale plant models: integrating cellular networks

One of the ambitions of 'crop systems biology' is to combine information from molecular biology with a broader view of plant development and growth. In the context of modeling, this calls for a multiscale perspective that focuses on the interplay between cellular and macroscopic studies. With this in mind, in this review we aim to draw attention to a panel of approaches that were developed in the context of systems biology and are used for analyzing and describing the behavior of cellular networks. Ultimately, insights obtained from these methods can be exploited to refine the description of plant processes, leading to integrated plant-cellular models.

Towards multiscale plant models: integrating cellular networks

One of the ambitions of 'crop systems biology' is to combine information from molecular biology with a broader view of plant development and growth. In the context of modeling, this calls for a multiscale perspective that focuses on the interplay between cellular and macroscopic studies. With this in mind, in this review we aim to draw attention to a panel of approaches that were developed in the context of systems biology and are used for analyzing and describing the behavior of cellular networks. Ultimately, insights obtained from these methods can be exploited to refine the description of plant processes, leading to integrated plant-cellular models.

MIRAGE: a functional genomics-based approach for metabolic network model reconstruction and its application to cyanobacteria networks

Genome-scale metabolic network reconstructions are considered a key step in quantifying the genotype-phenotype relationship. We present a novel gap-filling approach, MetabolIc Reconstruction via functionAl GEnomics (MIRAGE), which identifies missing network reactions by integrating metabolic flux analysis and functional genomics data. MIRAGE's performance is demonstrated on the reconstruction of metabolic network models of E. coli and Synechocystis sp. and validated via existing networks for these species. Then, it is applied to reconstruct genome-scale metabolic network models for 36 sequenced cyanobacteria amenable for constraint-based modeling analysis and specifically for metabolic engineering. The reconstructed network models are supplied via standard SBML files.

Metabolic Consequences of TGFb Stimulation in CulturedPrimary Mouse Hepatocytes Screened from Transcript Data with ModeScore

TGFβ signaling plays a major role in the reorganization of liver tissue upon injury and is an important driver of chronic liver disease. This is achieved by a deep impact on a cohort of cellular functions. To comprehensively assess the full range of affected metabolic functions, transcript changes of cultured mouse hepatocytes were analyzed with a novel method (ModeScore), which predicts the activity of metabolic functions by scoring transcript expression changes with 987 reference flux distributions, which yielded the following hypotheses. TGFβ multiplies down-regulation of most metabolic functions occurring in culture stressed controls. This is especially pronounced for tyrosine degradation, urea synthesis, glucuronization capacity, and cholesterol synthesis. Ethanol degradation and creatine synthesis are down-regulated only in TGFβ treated hepatocytes, but not in the control. Among the few TGFβ dependently up-regulated functions, synthesis of various collagens is most pronounced. Further interesting findings include: down-regulation of glucose export is postponed by TGFβ, TGFβ up-regulates the synthesis capacity of ketone bodies only as an early response, TGFβ suppresses the strong up-regulation of Vanin, and TGFβ induces re-formation of ceramides and sphingomyelin.

Reconstruction of Danio rerio metabolic model accounting for subcellular compartmentalisation

Plant and microbial metabolic engineering is commonly used in the production of functional foods and quality trait improvement. Computational model-based approaches have been used in this important endeavour. However, to date, fish metabolic models have only been scarcely and partially developed, in marked contrast to their prominent success in metabolic engineering. In this study we present the reconstruction of fully compartmentalised models of the Danio rerio (zebrafish) on a global scale. This reconstruction involves extraction of known biochemical reactions in D. rerio for both primary and secondary metabolism and the implementation of methods for determining subcellular localisation and assignment of enzymes. The reconstructed model (ZebraGEM) is amenable for constraint-based modelling analysis, and accounts for 4,988 genes coding for 2,406 gene-associated reactions and only 418 non-gene-associated reactions. A set of computational validations (i.e., simulations of known metabolic functionalities and experimental data) strongly testifies to the predictive ability of the model. Overall, the reconstructed model is expected to lay down the foundations for computational-based rational design of fish metabolic engineering in aquaculture.

Impact of stoichiometry representation on simulation of genotype-phenotype relationships in metabolic networks

Genome-scale metabolic networks provide a comprehensive structural framework for modeling genotype-phenotype relationships through flux simulations. The solution space for the metabolic flux state of the cell is typically very large and optimization-based approaches are often necessary for predicting the active metabolic state under specific environmental conditions. The objective function to be used in such optimization algorithms is directly linked with the biological hypothesis underlying the model and therefore it is one of the most relevant parameters for successful modeling. Although linear combination of selected fluxes is widely used for formulating metabolic objective functions, we show that the resulting optimization problem is sensitive towards stoichiometry representation of the metabolic network. This undesirable sensitivity leads to different simulation results when using numerically different but biochemically equivalent stoichiometry representations and thereby makes biological interpretation intrinsically subjective and ambiguous. We hereby propose a new method, Minimization of Metabolites Balance (MiMBl), which decouples the artifacts of stoichiometry representation from the formulation of the desired objective functions, by casting objective functions using metabolite turnovers rather than fluxes. By simulating perturbed metabolic networks, we demonstrate that the use of stoichiometry representation independent algorithms is fundamental for unambiguously linking modeling results with biological interpretation. For example, MiMBl allowed us to expand the scope of metabolic modeling in elucidating the mechanistic basis of several genetic interactions in Saccharomyces cerevisiae.

One of the challenging tasks in systems biology is to quantitatively predict the metabolic behavior of the cell under given genetic and environmental constraints. To this end, genome-scale metabolic reconstructions and simulation tools are indispensable. The choice of the objective function to be used for simulating genome-scale metabolic models is dependent on the biological context and one of the most relevant parameters for successful modeling. Formulation of the intended objective function often requires the use of multiple fluxes, e.g. the sum of fluxes through ATP-producing reactions. We demonstrate that the existing tools confound biological interpretation of the simulations due to undesired dependence on the representation of stoichiometry and propose a new tool – Minimization of Metabolites Balance (MiMBl). MiMBl allows casting of the desired biological objective functions into linear optimization models and gives consistent simulation results when using numerically different but biochemically equivalent stoichiometry representations. We demonstrate relevance of MiMBl for addressing biological questions through improved predictions of genetic interactions within the yeast metabolic network. Genetic interactions imply functional relationship between the genes and therefore allow assessing different hypotheses for the underlying biological principles. MiMBl explains several of the genetic interactions as outcome of flux re-routing for minimal metabolite turnover adjustments.

Mechanistic modeling of aberrant energy metabolism in human disease

Dysfunction in energy metabolism—including in pathways localized to the mitochondria—has been implicated in the pathogenesis of a wide array of disorders, ranging from cancer to neurodegenerative diseases to type II diabetes. The inherent complexities of energy and mitochondrial metabolism present a significant obstacle in the effort to understand the role that these molecular processes play in the development of disease. To help unravel these complexities, systems biology methods have been applied to develop an array of computational metabolic models, ranging from mitochondria-specific processes to genome-scale cellular networks. These constraint-based (CB) models can efficiently simulate aspects of normal and aberrant metabolism in various genetic and environmental conditions. Development of these models leverages—and also provides a powerful means to integrate and interpret—information from a wide range of sources including genomics, proteomics, metabolomics, and enzyme kinetics. Here, we review a variety of mechanistic modeling studies that explore metabolic functions, deficiency disorders, and aberrant biochemical pathways in mitochondria and related regions in the cell.

Integrating cellular metabolism into a multiscale whole-body model

Cellular metabolism continuously processes an enormous range of external compounds into endogenous metabolites and is as such a key element in human physiology. The multifaceted physiological role of the metabolic network fulfilling the catalytic conversions can only be fully understood from a whole-body perspective where the causal interplay of the metabolic states of individual cells, the surrounding tissue and the whole organism are simultaneously considered. We here present an approach relying on dynamic flux balance analysis that allows the integration of metabolic networks at the cellular scale into standardized physiologically-based pharmacokinetic models at the whole-body level. To evaluate our approach we integrated a genome-scale network reconstruction of a human hepatocyte into the liver tissue of a physiologically-based pharmacokinetic model of a human adult. The resulting multiscale model was used to investigate hyperuricemia therapy, ammonia detoxification and paracetamol-induced toxication at a systems level. The specific models simultaneously integrate multiple layers of biological organization and offer mechanistic insights into pathology and medication. The approach presented may in future support a mechanistic understanding in diagnostics and drug development.

Cellular metabolism is a key element in human physiology. Ideally the metabolic network needs to be considered within the context of the surrounding tissue and organism since the various levels of biological organization are mutually influencing each other. To mechanistically describe the interplay between intracellular space and extracellular environment, we here integrate the genome-scale metabolic network model HepatoNet1 at the cellular scale into physiologically-based pharmacokinetic models at the whole-body level. The resulting multiscale model allows the quantitative description of metabolic behavior in the context of time-resolved metabolite concentration profiles in the body and the surrounding liver tissue. The model has been applied to three case studies covering fundamental aspects of medicine and pharmacology: drug administration, biomarker identification and drug-induced toxication. Most notably, our multiscale approach fosters an improved quantitative understanding of drug action and the impact of metabolic disorders at an organism level, based on a genome-scale representation of cellular metabolism. Computational models such as the one presented include various aspects of human physiology and may therefore significantly support rational approaches in medical diagnostics and pharmaceutical drug development in the future.

Enhancing the prioritization of disease-causing genes through tissue specific protein interaction networks

The prioritization of candidate disease-causing genes is a fundamental challenge in the post-genomic era. Current state of the art methods exploit a protein-protein interaction (PPI) network for this task. They are based on the observation that genes causing phenotypically-similar diseases tend to lie close to one another in a PPI network. However, to date, these methods have used a static picture of human PPIs, while diseases impact specific tissues in which the PPI networks may be dramatically different. Here, for the first time, we perform a large-scale assessment of the contribution of tissue-specific information to gene prioritization. By integrating tissue-specific gene expression data with PPI information, we construct tissue-specific PPI networks for 60 tissues and investigate their prioritization power. We find that tissue-specific PPI networks considerably improve the prioritization results compared to those obtained using a generic PPI network. Furthermore, they allow predicting novel disease-tissue associations, pointing to sub-clinical tissue effects that may escape early detection.

Mechanistic systems modeling to guide drug discovery and development

A crucial question that must be addressed in the drug development process is whether the proposed therapeutic target will yield the desired effect in the clinical population. Pharmaceutical and biotechnology companies place a large investment on research and development, long before confirmatory data are available from human trials. Basic science has greatly expanded the computable knowledge of disease processes, both through the generation of large omics data sets and a compendium of studies assessing cellular and systemic responses to physiologic and pathophysiologic stimuli. Given inherent uncertainties in drug development, mechanistic systems models can better inform target selection and the decision process for advancing compounds through preclinical and clinical research.

What mRNA Abundances Can Tell us about Metabolism

Inferring decreased or increased metabolic functions from transcript profiles is at first sight a bold and speculative attempt because of the functional layers in between: proteins, enzymatic activities, and reaction fluxes. However, the growing interest in this field can easily be explained by two facts: the high quality of genome-scale metabolic network reconstructions and the highly developed technology to obtain genome-covering RNA profiles. Here, an overview of important algorithmic approaches is given by means of criteria by which published procedures can be classified. The frontiers of the methods are sketched and critical voices are being heard. Finally, an outlook for the prospects of the field is given.

Metabolomics view on gut microbiome modulation by polyphenol-rich foods

Health is influenced by genetic, lifestyle, and diet determinants; therefore, nutrition plays an essential role in health management. Still, the substantiation of nutritional health benefits is challenged by the intrinsic macro- and micronutrient complexity of foods and individual responses. Evidence of healthy effects of food requires new strategies not only to stratify populations according to their metabolic requirements but also to predict and measure individual responses to dietary intakes. The influence of the gut microbiome and its interaction with the host is pivotal to understand nutrition and metabolism. Thus, the modulation of the gut microbiome composition by alteration of food habits has potentialities in health improvement or even disease prevention. Dietary polyphenols are naturally occurring constituents in vegetables and fruits, including coffee and cocoa. They are commonly associated to health benefits, although mechanistic evidence in vivo is not yet fully understood. Polyphenols are extensively metabolized by gut bacteria into a complex series of end-products that support a significant effect on the functional ecology of symbiotic partners that can affect the host physiology. This review reports recent nutritional metabolomics inspections of gut microbiota-host metabolic interactions with a particular focus on the cometabolism of cocoa and coffee polyphenols.

CardioNet: a human metabolic network suited for the study of cardiomyocyte metabolism

Background

Availability of oxygen and nutrients in the coronary circulation is a crucial determinant of cardiac performance. Nutrient composition of coronary blood may significantly vary in specific physiological and pathological conditions, for example, administration of special diets, long-term starvation, physical exercise or diabetes. Quantitative analysis of cardiac metabolism from a systems biology perspective may help to a better understanding of the relationship between nutrient supply and efficiency of metabolic processes required for an adequate cardiac output.

Results

Here we present CardioNet, the first large-scale reconstruction of the metabolic network of the human cardiomyocyte comprising 1793 metabolic reactions, including 560 transport processes in six compartments. We use flux-balance analysis to demonstrate the capability of the network to accomplish a set of 368 metabolic functions required for maintaining the structural and functional integrity of the cell. Taking the maintenance of ATP, biosynthesis of ceramide, cardiolipin and further important phospholipids as examples, we analyse how a changed supply of glucose, lactate, fatty acids and ketone bodies may influence the efficiency of these essential processes.

Conclusions

CardioNet is a functionally validated metabolic network of the human cardiomyocyte that enables theorectical studies of cellular metabolic processes crucial for the accomplishment of an adequate cardiac output.

Systems biology of stored blood cells: can it help to extend the expiration date?

With increasingly stringent regulations regarding deferral and elimination of blood donors it will become increasingly important to extend the expiration date of blood components beyond the current allowed storage periods. One reason for the storage time limit for blood components is that platelets and red blood cells develop a condition called storage lesions during their storage in plastic blood containers. Systems biology provides comprehensive bio-chemical descriptions of organisms through quantitative measurements and data integration in mathematical models. The biological knowledge for a target organism can be translated in a mathematical format and used to compute physiological properties. The use of systems biology represents a concrete solution in the study of blood cell storage lesions, and it may open up new avenues towards developing better storage methods and better storage media, thereby extending the storage period of blood components. This article is part of a Special Issue entitled: Integrated omics.

Detecting and investigating substrate cycles in a genome-scale human metabolic network

Substrate cycles, also known as futile cycles, are cyclic metabolic routes that dissipate energy by hydrolysing cofactors such as ATP. They were first described to occur in the muscles of bumblebees and brown adipose tissue in the 1970s. A popular example is the conversion of fructose 6-phosphate to fructose 1,6-bisphosphate and back. In the present study, we analyze a large number of substrate cycles in human metabolism that consume ATP and discuss their statistics. For this purpose, we use two recently published methods (i.e. EFMEvolver and the K-shortest EFM method) to calculate samples of 100,000 and 15,000 substrate cycles, respectively. We find an unexpectedly high number of substrate cycles in human metabolism, with up to 100 reactions per cycle, utilizing reactions from up to six different compartments. An analysis of tissue-specific models of liver and brain metabolism shows that there is selective pressure that acts against the uncontrolled dissipation of energy by avoiding the coexpression of enzymes belonging to the same substrate cycle. This selective force is particularly strong against futile cycles that have a high flux as a result of thermodynamic principles.

Integration of expression data in genome-scale metabolic network reconstructions

With the advent of high-throughput technologies, the field of systems biology has amassed an abundance of "omics" data, quantifying thousands of cellular components across a variety of scales, ranging from mRNA transcript levels to metabolite quantities. Methods are needed to not only integrate this omics data but to also use this data to heighten the predictive capabilities of computational models. Several recent studies have successfully demonstrated how flux balance analysis (FBA), a constraint-based modeling approach, can be used to integrate transcriptomic data into genome-scale metabolic network reconstructions to generate predictive computational models. In this review, we summarize such FBA-based methods for integrating expression data into genome-scale metabolic network reconstructions, highlighting their advantages as well as their limitations.

Identification of aberrant pathways and network activities from high-throughput data

Many complex diseases such as cancer are associated with changes in biological pathways and molecular networks rather than being caused by single gene alterations. A major challenge in the diagnosis and treatment of such diseases is to identify characteristic aberrancies in the biological pathways and molecular network activities and elucidate their relationship to the disease. This review presents recent progress in using high-throughput biological assays to decipher aberrant pathways and network activities. In particular, this review provides specific examples in which high-throughput data have been applied to identify relationships between diseases and aberrant pathways and network activities. The achievements in this field have been remarkable, but many challenges have yet to be addressed.

Network biology methods integrating biological data for translational science

The explosion of biomedical data, both on the genomic and proteomic side as well as clinical data, will require complex integration and analysis to provide new molecular variables to better understand the molecular basis of phenotype. Currently, much data exist in silos and is not analyzed in frameworks where all data are brought to bear in the development of biomarkers and novel functional targets. This is beginning to change. Network biology approaches, which emphasize the interactions between genes, proteins and metabolites provide a framework for data integration such that genome, proteome, metabolome and other -omics data can be jointly analyzed to understand and predict disease phenotypes. In this review, recent advances in network biology approaches and results are identified. A common theme is the potential for network analysis to provide multiplexed and functionally connected biomarkers for analyzing the molecular basis of disease, thus changing our approaches to analyzing and modeling genome- and proteome-wide data.

Mathematical and statistical modeling in cancer systems biology

Cancer is a major health problem with high mortality rates. In the post-genome era, investigators have access to massive amounts of rapidly accumulating high-throughput data in publicly available databases, some of which are exclusively devoted to housing Cancer data. However, data interpretation efforts have not kept pace with data collection, and gained knowledge is not necessarily translating into better diagnoses and treatments. A fundamental problem is to integrate and interpret data to further our understanding in Cancer Systems Biology. Viewing cancer as a network provides insights into the complex mechanisms underlying the disease. Mathematical and statistical models provide an avenue for cancer network modeling. In this article, we review two widely used modeling paradigms: deterministic metabolic models and statistical graphical models. The strength of these approaches lies in their flexibility and predictive power. Once a model has been validated, it can be used to make predictions and generate hypotheses. We describe a number of diverse applications to Cancer Biology, including, the system-wide effects of drug-treatments, disease prognosis, tumor classification, forecasting treatment outcomes, and survival predictions.

Model-driven multi-omic data analysis elucidates metabolic immunomodulators of macrophage activation

Macrophages are central players in immune response, manifesting divergent phenotypes to control inflammation and innate immunity through release of cytokines and other signaling factors. Recently, the focus on metabolism has been reemphasized as critical signaling and regulatory pathways of human pathophysiology, ranging from cancer to aging, often converge on metabolic responses. Here, we used genome-scale modeling and multi-omics (transcriptomics, proteomics, and metabolomics) analysis to assess metabolic features that are critical for macrophage activation. We constructed a genome-scale metabolic network for the RAW 264.7 cell line to determine metabolic modulators of activation. Metabolites well-known to be associated with immunoactivation (glucose and arginine) and immunosuppression (tryptophan and vitamin D3) were among the most critical effectors. Intracellular metabolic mechanisms were assessed, identifying a suppressive role for de-novo nucleotide synthesis. Finally, underlying metabolic mechanisms of macrophage activation are identified by analyzing multi-omic data obtained from LPS-stimulated RAW cells in the context of our flux-based predictions. Our study demonstrates metabolism's role in regulating activation may be greater than previously anticipated and elucidates underlying connections between activation and metabolic effectors.

A linearized constraint-based approach for modeling signaling networks

With the unparalleled increase in the availability of biological data over the last couple of decades, accurate and computable models are becoming increasingly important for unraveling complex biological phenomena. Past efforts to model signaling networks have utilized various computational methods, including Boolean and constraint-based modeling (CBM) approaches. These approaches are based on solving mixed integer linear programs; hence, they may not scale up for the analysis of large networks and are not amenable for applications based on sampling the full spectrum of the solution space. Here we propose a new CBM approach that is fully linear and does not involve integer variables, thereby overcoming the aforementioned limitations. We describe a novel optimization procedure for model construction and demonstrate the utility of our approach on a reconstructed model of the human epidermal growth factor receptor (EGFR) pathway, spanning 322 species and 211 connections. We compare our model's predictions to experimental phosphorylation data and to the predictions inferred via an additional Boolean-based EGFR signaling model. Our results show high prediction accuracy (75%) and high similarity to the Boolean model. Considering the marked computational advantages in terms of scalability and sampling utilization obtained by having a linear mode, these results demonstrate the potential promise of this framework for the study of cellular signaling.

Reconstruction of genome-scale active metabolic networks for 69 human cell types and 16 cancer types using INIT

Development of high throughput analytical methods has given physicians the potential access to extensive and patient-specific data sets, such as gene sequences, gene expression profiles or metabolite footprints. This opens for a new approach in health care, which is both personalized and based on system-level analysis. Genome-scale metabolic networks provide a mechanistic description of the relationships between different genes, which is valuable for the analysis and interpretation of large experimental data-sets. Here we describe the generation of genome-scale active metabolic networks for 69 different cell types and 16 cancer types using the INIT (Integrative Network Inference for Tissues) algorithm. The INIT algorithm uses cell type specific information about protein abundances contained in the Human Proteome Atlas as the main source of evidence. The generated models constitute the first step towards establishing a Human Metabolic Atlas, which will be a comprehensive description (accessible online) of the metabolism of different human cell types, and will allow for tissue-level and organism-level simulations in order to achieve a better understanding of complex diseases. A comparative analysis between the active metabolic networks of cancer types and healthy cell types allowed for identification of cancer-specific metabolic features that constitute generic potential drug targets for cancer treatment.

A network perspective on metabolic inconsistency

BACKGROUND

Integrating gene expression profiles and metabolic pathways under different experimental conditions is essential for understanding the coherence of these two layers of cellular organization. The network character of metabolic systems can be instrumental in developing concepts of agreement between expression data and pathways. A network-driven interpretation of gene expression data has the potential of suggesting novel classifiers for pathological cellular states and of contributing to a general theoretical understanding of gene regulation.

RESULTS

Here, we analyze the coherence of gene expression patterns and a reconstruction of human metabolism, using consistency scores obtained from network and constraint-based analysis methods. We find a surprisingly strong correlation between the two measures, demonstrating that a substantial part of inconsistencies between metabolic processes and gene expression can be understood from a network perspective alone. Prompted by this finding, we investigate the topological context of the individual biochemical reactions responsible for the observed inconsistencies. On this basis, we are able to separate the differential contributions that bear physiological information about the system, from the unspecific contributions that unravel gaps in the metabolic reconstruction. We demonstrate the biological potential of our network-driven approach by analyzing transcriptome profiles of aldosterone producing adenomas that have been obtained from a cohort of Primary Aldosteronism patients. We unravel systematics in the data that could not have been resolved by conventional microarray data analysis. In particular, we discover two distinct metabolic states in the adenoma expression patterns.

CONCLUSIONS

The methodology presented here can help understand metabolic inconsistencies from a network perspective. It thus serves as a mediator between the topology of metabolic systems and their dynamical function. Finally, we demonstrate how physiologically relevant insights into the structure and dynamics of metabolic networks can be obtained using this novel approach.

Applications of system-level models of metabolism for analysis of bacterial physiology and identification of new drug targets

For nearly all of the 20th century, biologists gained considerable insights into the fundamental principles of cellular dynamics by examining select modules of biochemical processes. This form of analysis provides detailed information about the workings of the examined pathways. However, any attempt to alter the normal function of bacteria (perhaps for industrial or medicinal goals) requires a detailed global understanding of cellular mechanisms. The reductionist mode of analysis cannot provide the required information for developing the needed perspective on the complex interactions of biochemical pathways. Thankfully, the increasing availability of microbial genomic, transcriptomic, proteomic and other high-throughput data permits system-level analyses of microbiology. During the past two decades, systems biologists have developed constraint-based genome-scale models (GSM) of metabolism for a variety of pathogens. These models are important tools for assessing the metabolic capabilities of various genotypes. Simultaneously, new computational methods have been developed that use these network reconstructions to answer an array of important immunological questions. The objective of this article is to briefly review some of the uses of GSMs for studying bacterial metabolism under different conditions and to discuss how the calculated solutions can be used for rational design of drugs.

Advanced stoichiometric analysis of metabolic networks of mammalian systems

Metabolic engineering tools have been widely applied to living organisms to gain a comprehensive understanding about cellular networks and to improve cellular properties. Metabolic flux analysis (MFA), flux balance analysis (FBA), and metabolic pathway analysis (MPA) are among the most popular tools in stoichiometric network analysis. Although application of these tools into well-known microbial systems is extensive in the literature, various barriers prevent them from being utilized in mammalian cells. Limited experimental data, complex regulatory mechanisms, and the requirement of more complex nutrient media are some major obstacles in mammalian cell systems. However, mammalian cells have been used to produce therapeutic proteins, to characterize disease states or related abnormal metabolic conditions, and to analyze the toxicological effects of some medicinally important drugs. Therefore, there is a growing need for extending metabolic engineering principles to mammalian cells in order to understand their underlying metabolic functions. In this review article, advanced metabolic engineering tools developed for stoichiometric analysis including MFA, FBA, and MPA are described. Applications of these tools in mammalian cells are discussed in detail, and the challenges and opportunities are highlighted.

Integrative genomics strategies to elucidate the complexity of drug response

Pharmacogenomic investigation from both genome-wide association studies and experiments focused on candidate loci involved in drug mechanism and metabolism has yielded a substantial and increasing list of robust genetic effects on drug therapy in humans. At the same time, reasonably comprehensive molecular data such as gene expression, proteomic and metabolomic data are now available for collections of hundreds to thousands of individuals. If these data are structured in a statistically robust and computationally tractable way, such as a network model, they can aid in the analysis of new pharmacogenomics studies by suggesting novel hypotheses for the regulation of genes involved in drug metabolism and response. Similarly, hypotheses taken from these same models can direct genome-wide association studies by focusing the genome-wide association studies analysis on a number of specific hypotheses informed by the relationships customarily seen between a gene's expression or protein activity and genetic variation at a particular locus. Network models based on other sorts of systematic biological data such as cell-based surveys of drug effect on gene expression and mining of literature and electronic medical records for associations between clinical and molecular phenotypes also promise similar utility. Although surely primitive in comparison with what will be developed, these model-based approaches to leveraging the increasing volume of data generated in the course of patient care and medical research nevertheless suggest a huge opportunity to improve our understanding of biological systems involved in pharmacogenomics and apply them to questions of medical relevance.

Stoichiometry based steady-state hepatic flux analysis: computational and experimental aspects

: The liver has many complex physiological functions, including lipid, protein and carbohydrate metabolism, as well as bile and urea production. It detoxifies toxic substances and medicinal products. It also plays a key role in the onset and maintenance of abnormal metabolic patterns associated with various disease states, such as burns, infections and major traumas. Liver cells have been commonly used in in vitro experiments to elucidate the toxic effects of drugs and metabolic changes caused by aberrant metabolic conditions, and to improve the functions of existing systems, such as bioartificial liver. More recently, isolated liver perfusion systems have been increasingly used to characterize intrinsic metabolic changes in the liver caused by various perturbations, including systemic injury, hepatotoxin exposure and warm ischemia. Metabolic engineering tools have been widely applied to these systems to identify metabolic flux distributions using metabolic flux analysis or flux balance analysis and to characterize the topology of the networks using metabolic pathway analysis. In this context, hepatic metabolic models, together with experimental methodologies where hepatocytes or perfused livers are mainly investigated, are described in detail in this review. The challenges and opportunities are also discussed extensively.

Human metabolic network: reconstruction, simulation, and applications in systems biology

Metabolism is crucial to cell growth and proliferation. Deficiency or alterations in metabolic functions are known to be involved in many human diseases. Therefore, understanding the human metabolic system is important for the study and treatment of complex diseases. Current reconstructions of the global human metabolic network provide a computational platform to integrate genome-scale information on metabolism. The platform enables a systematic study of the regulation and is applicable to a wide variety of cases, wherein one could rely on in silico perturbations to predict novel targets, interpret systemic effects, and identify alterations in the metabolic states to better understand the genotype-phenotype relationships. In this review, we describe the reconstruction of the human metabolic network, introduce the constraint based modeling approach to analyze metabolic networks, and discuss systems biology applications to study human physiology and pathology. We highlight the challenges and opportunities in network reconstruction and systems modeling of the human metabolic system.

Frontiers in metabolic reconstruction and modeling of plant genomes

A major goal of post-genomic biology is to reconstruct and model in silico the metabolic networks of entire organisms. Work on bacteria is well advanced, and is now under way for plants and other eukaryotes. Genome-scale modelling in plants is much more challenging than in bacteria. The challenges come from features characteristic of higher organisms (subcellular compartmentation, tissue differentiation) and also from the particular severity in plants of a general problem: genome content whose functions remain undiscovered. This problem results in thousands of genes for which no function is known ('undiscovered genome content') and hundreds of enzymatic and transport functions for which no gene is yet identified. The severity of the undiscovered genome content problem in plants reflects their genome size and complexity. To bring the challenges of plant genome-scale modelling into focus, we first summarize the current status of plant genome-scale models. We then highlight the challenges - and ways to address them - in three areas: identifying genes for missing processes, modelling tissues as opposed to single cells, and finding metabolic functions encoded by undiscovered genome content. We also discuss the emerging view that a significant fraction of undiscovered genome content encodes functions that counter damage to metabolites inflicted by spontaneous chemical reactions or enzymatic mistakes.

Constraining the metabolic genotype-phenotype relationship using a phylogeny of in silico methods

Reconstructed microbial metabolic networks facilitate a mechanistic description of the genotype-phenotype relationship through the deployment of constraint-based reconstruction and analysis (COBRA) methods. As reconstructed networks leverage genomic data for insight and phenotype prediction, the development of COBRA methods has accelerated following the advent of whole-genome sequencing. Here, we describe a phylogeny of COBRA methods that has rapidly evolved from the few early methods, such as flux balance analysis and elementary flux mode analysis, into a repertoire of more than 100 methods. These methods have enabled genome-scale analysis of microbial metabolism for numerous basic and applied uses, including antibiotic discovery, metabolic engineering and modelling of microbial community behaviour.

Using the reconstructed genome-scale human metabolic network to study physiology and pathology

Metabolism plays a key role in many major human diseases. Generation of high-throughput omics data has ushered in a new era of systems biology. Genome-scale metabolic network reconstructions provide a platform to interpret omics data in a biochemically meaningful manner. The release of the global human metabolic network, Recon 1, in 2007 has enabled new systems biology approaches to study human physiology, pathology and pharmacology. There are currently more than 20 publications that utilize Recon 1, including studies of cancer, diabetes, host-pathogen interactions, heritable metabolic disorders and off-target drug binding effects. In this mini-review, we focus on the reconstruction of the global human metabolic network and four classes of its application. We show that computational simulations for numerous pathologies have yielded clinically relevant results, many corroborated by existing or newly generated experimental data.

Systems medicine and metabolic modelling

Several complex diseases are caused by the malfunction of human metabolism, and deciphering the underlying molecular mechanisms can elucidate their aetiology. Systems biology is an integrative approach combining experimental and computational biology to identify and describe the molecular mechanisms of complex biological systems. Systems medicine has the potential to elucidate the onset and progression of complex metabolic diseases through the use of computational approaches. Advances in biotechnology have resulted in the provision of high-throughput data, which provide information about different metabolic processes. The systems medicine approach can utilize such data to reconstruct genome-scale metabolic models that can be used to study the function of specific enzymes and pathways in the context of the complete metabolic network. In this review, we outline the importance of genome-scale models in systems medicine and discuss how they may contribute towards the development of personalized medicine.

Individualized therapy of HHT driven by network analysis of metabolomic profiles

Background

Hereditary Hemorrhagic Telangiectasia (HHT) is an autosomal dominant disease with a varying range of phenotypes involving abnormal vasculature primarily manifested as arteriovenous malformations in various organs, including the nose, brain, liver, and lungs. The varied presentation and involvement of different organ systems makes the choice of potential treatment medications difficult.

Results

A patient with a mixed-clinical presentation and presumed diagnosis of HHT, severe exertional dyspnea, and diffuse pulmonary shunting at the microscopic level presented for treatment. We sought to analyze her metabolomic plasma profile to assist with pharmacologic treatment selection. Fasting serum samples from 5 individuals (4 healthy and 1 with HHT) were metabolomically profiled.

A global metabolic network reconstruction, Recon 1, was used to help guide the choice of medication via analysis of the differential metabolism between the patient and healthy controls using metabolomic data. Flux Balance Analysis highlighted changes in metabolic pathway activity, notably in nitric oxide synthase (NOS), which suggested a potential link between changes in vascular endothelial function and metabolism. This finding supported the use of an already approved medication, bevacizumab (Avastin). Following 2 months of treatment, the patient's metabolic profile shifted, becoming more similar to the control subject profiles, suggesting that the treatment was addressing at least part of the pathophysiological state.

Conclusions

In this 'individualized case study' of personalized medicine, we carry out untargeted metabolomic profiling of a patient and healthy controls. Rather than filtering the data down to a single value, these data are analyzed in the context of a network model of metabolism, in order to simulate the biochemical phenotypic differences between healthy and disease states; the results then guide the therapy. This presents one approach to achieving the goals of individualized medicine through Systems Biology and causal models analysis.