Using Bioconductor Package BiGGR for Metabolic Flux Estimation Based on Gene Expression Changes in Brain

Predicting the distribution of metabolic fluxes in biochemical networks is of major interest in systems biology. Several databases provide metabolic reconstructions for different organisms. Software to analyze flux distributions exists, among others for the proprietary MATLAB environment. Given the large user community for the R computing environment, a simple implementation of flux analysis in R appears desirable and will facilitate easy interaction with computational tools to handle gene expression data. We extended the R software package BiGGR, an implementation of metabolic flux analysis in R. BiGGR makes use of public metabolic reconstruction databases, and contains the BiGG database and the reconstruction of human metabolism Recon2 as Systems Biology Markup Language (SBML) objects. Models can be assembled by querying the databases for pathways, genes or reactions of interest. Fluxes can then be estimated by maximization or minimization of an objective function using linear inverse modeling algorithms. Furthermore, BiGGR provides functionality to quantify the uncertainty in flux estimates by sampling the constrained multidimensional flux space. As a result, ensembles of possible flux configurations are constructed that agree with measured data within precision limits. BiGGR also features automatic visualization of selected parts of metabolic networks using hypergraphs, with hyperedge widths proportional to estimated flux values. BiGGR supports import and export of models encoded in SBML and is therefore interoperable with different modeling and analysis tools. As an application example, we calculated the flux distribution in healthy human brain using a model of central carbon metabolism. We introduce a new algorithm termed Least-squares with equalities and inequalities Flux Balance Analysis (Lsei-FBA) to predict flux changes from gene expression changes, for instance during disease. Our estimates of brain metabolic flux pattern with Lsei-FBA for Alzheimer’s disease agree with independent measurements of cerebral metabolism in patients. This second version of BiGGR is available from Bioconductor.

[1]  M A Holliday,et al.  Metabolic rate and organ size during growth from infancy to maturity and during late gastation and early infancy. , 1971, Pediatrics.

[2]  Holliday Ma,et al.  Metabolic rate and organ size during growth from infancy to maturity and during late gastation and early infancy. , 1971, Pediatrics.

[3]  B. Persson,et al.  Cerebral blood flow and metabolic rate of oxygen, glucose, lactate, pyruvate, ketone bodies and amino acids , 1980, Acta neurologica Scandinavica.

[4]  K. Jellinger,et al.  Buspirone and dystonia , 1991, Neurology.

[5]  S. Hoyer.  Oxidative energy metabolism in Alzheimer brain. Studies in early-onset and late-onset cases. , 1992, Molecular and chemical neuropathology.

[6]  K. Ishii,et al.  Decreased medial temporal oxygen metabolism in Alzheimer's disease shown by PET. , 1996, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[7]  Peter Green,et al.  Markov chain Monte Carlo in Practice , 1996 .

[8]  K. Ishii,et al.  Reduction of cerebellar glucose metabolism in advanced Alzheimer's disease. , 1997, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[9]  Hubert Blain,et al.  Abnormal Substrate Levels That Depend upon Mitochondrial Function in Cerebrospinal Fluid from Alzheimer Patients , 1998, Gerontology.

[10]  R. Mahadevan,et al.  The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. , 2003, Metabolic engineering.

[11]  Hiroaki Kitano,et al.  The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models , 2003, Bioinform..

[12]  Jean YH Yang,et al.  Bioconductor: open software development for computational biology and bioinformatics , 2004, Genome Biology.

[13]  I. Hozo,et al.  Estimating the mean and variance from the median, range, and the size of a sample , 2005, BMC medical research methodology.

[14]  Gordon K. Smyth,et al.  limma: Linear Models for Microarray Data , 2005 .

[15]  R. Shulman,et al.  The contribution of GABA to glutamate/glutamine cycling and energy metabolism in the rat cortex in vivo. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Steffen Klamt,et al.  Structural and functional analysis of cellular networks with CellNetAnalyzer , 2007, BMC Systems Biology.

[17]  D. Hovda,et al.  Increased Pentose Phosphate Pathway Flux after Clinical Traumatic Brain Injury: A [1,2-13C2]glucose Labeling Study in Humans , 2007, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[18]  Daniela Calvetti,et al.  Bayesian flux balance analysis applied to a skeletal muscle metabolic model. , 2007, Journal of theoretical biology.

[19]  Christopher R. Myers,et al.  Universally Sloppy Parameter Sensitivities in Systems Biology Models , 2007, PLoS Comput. Biol..

[20]  Eric M Reiman,et al.  Gene expression profiles in anatomically and functionally distinct regions of the normal aged human brain. , 2007, Physiological genomics.

[21]  K. Ülgen,et al.  Reconstruction and flux analysis of coupling between metabolic pathways of astrocytes and neurons: application to cerebral hypoxia , 2007, Theoretical Biology and Medical Modelling.

[22]  Bernhard O. Palsson,et al.  Context-Specific Metabolic Networks Are Consistent with Experiments , 2008, PLoS Comput. Biol..

[23]  Winnie S. Liang,et al.  Alzheimer's disease is associated with reduced expression of energy metabolism genes in posterior cingulate neurons , 2008, Proceedings of the National Academy of Sciences.

[24]  Karline Soetaert,et al.  Quantifying Food Web Flows Using Linear Inverse Models , 2009, Ecosystems.

[25]  John P. Lowry,et al.  An integrative dynamic model of brain energy metabolism using in vivo neurochemical measurements , 2009, Journal of Computational Neuroscience.

[26]  Bernhard O. Palsson,et al.  BiGG: a Biochemical Genetic and Genomic knowledgebase of large scale metabolic reconstructions , 2010, BMC Bioinformatics.

[27]  Karline Soetaert,et al.  xsample(): An R Function for Sampling Linear Inverse Problems , 2009 .

[28]  Miguel Rocha,et al.  OptFlux: an open-source software platform for in silico metabolic engineering , 2010, BMC Systems Biology.

[29]  Daniela Calvetti,et al.  Astrocytes as the glucose shunt for glutamatergic neurons at high activity: an in silico study. , 2009, Journal of neurophysiology.

[30]  Steffen Klamt,et al.  Hypergraphs and Cellular Networks , 2009, PLoS Comput. Biol..

[31]  Falk Schreiber,et al.  FBA-SimVis: interactive visualization of constraint-based metabolic models , 2009, Bioinform..

[32]  Desmond S. Lun,et al.  Interpreting Expression Data with Metabolic Flux Models: Predicting Mycobacterium tuberculosis Mycolic Acid Production , 2009, PLoS Comput. Biol..

[33]  Eytan Ruppin,et al.  iMAT: an integrative metabolic analysis tool , 2010, Bioinform..

[34]  Ines Thiele,et al.  Computationally efficient flux variability analysis , 2010, BMC Bioinformatics.

[35]  Daniela Calvetti,et al.  Metabolica: A statistical research tool for analyzing metabolic networks , 2010, Comput. Methods Programs Biomed..

[36]  Neil Swainston,et al.  Integration of metabolic databases for the reconstruction of genome-scale metabolic networks , 2010, BMC Systems Biology.

[37]  Jeffrey D Orth,et al.  What is flux balance analysis? , 2010, Nature Biotechnology.

[38]  B. Palsson,et al.  Large-scale in silico modeling of metabolic interactions between cell types in the human brain , 2010, Nature Biotechnology.

[39]  Jason A. Papin,et al.  Functional integration of a metabolic network model and expression data without arbitrary thresholding , 2011, Bioinform..

[40]  Michèle Allard,et al.  Brain fuel metabolism, aging, and Alzheimer's disease. , 2011, Nutrition (Burbank, Los Angeles County, Calif.).

[41]  Andrzej M. Kierzek,et al.  SurreyFBA: a command line tool and graphics user interface for constraint-based modeling of genome-scale metabolic reaction networks , 2011, Bioinform..

[42]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.

[43]  Joost Boele,et al.  FAME, the Flux Analysis and Modeling Environment , 2012, BMC Systems Biology.

[44]  Hannes Hettling,et al.  Simulating the physiology of athletes during endurance sports events: modelling human energy conversion and metabolism , 2011, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[45]  Neil Swainston,et al.  Improving metabolic flux predictions using absolute gene expression data , 2012, BMC Systems Biology.

[46]  Joshua A. Lerman,et al.  COBRApy: COnstraints-Based Reconstruction and Analysis for Python , 2013, BMC Systems Biology.

[47]  Jaques Reifman,et al.  Modeling Phenotypic Metabolic Adaptations of Mycobacterium tuberculosis H37Rv under Hypoxia , 2012, PLoS Comput. Biol..

[48]  Ali Navid,et al.  Genome-level transcription data of Yersinia pestis analyzed with a New metabolic constraint-based approach , 2012, BMC Systems Biology.

[49]  Jaap Heringa,et al.  Computational estimation of tricarboxylic acid cycle fluxes using noisy NMR data from cardiac biopsies , 2013, BMC Systems Biology.

[50]  Martin J. Lercher,et al.  sybil – Efficient constraint-based modelling in R , 2013, BMC Systems Biology.

[51]  Intawat Nookaew,et al.  The RAVEN Toolbox and Its Use for Generating a Genome-scale Metabolic Model for Penicillium chrysogenum , 2013, PLoS Comput. Biol..

[52]  Nicolas Le Novère,et al.  BioModels Database: a repository of mathematical models of biological processes. , 2013, Methods in molecular biology.

[53]  Ronan M. T. Fleming,et al.  A community-driven global reconstruction of human metabolism , 2013, Nature Biotechnology.

[54]  Daniel Machado,et al.  Systematic Evaluation of Methods for Integration of Transcriptomic Data into Constraint-Based Models of Metabolism , 2014, PLoS Comput. Biol..

[55]  Dong-Yup Lee,et al.  Software applications for flux balance analysis , 2014, Briefings Bioinform..

[56]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[57]  A. Stefani,et al.  CSF lactate levels, τ proteins, cognitive decline: a dynamic relationship in Alzheimer's disease , 2014, Journal of Neurology, Neurosurgery & Psychiatry.