Metabolic modeling predicts metabolite changes in Mycobacterium tuberculosis

BackgroundMycobacterium tuberculosis (MTB) is the causal agent of the disease tuberculosis (TB). Metabolic adaptations are thought to be critical to the survival of MTB during pathogenesis. Computational tools that can be used to study MTB metabolism in silico and prioritize resource-intensive experimental work could significantly accelerate research.ResultsWe have developed E-Flux-MFC, an enhancement of our original E-Flux method that enables the prediction of changes in the production of external and internal metabolites corresponding to gene expression measurements. We have used this method to simulate the changes in the metabolic state of Mycobacterium tuberculosis (MTB). We have validated the accuracy of E-Flux-MFC for predicting changes in lipids and metabolites during a hypoxia time course using previously published metabolomics and transcriptomics data. We have further validated the accuracy of the method for predicting changes in MTB lipids following the deletion and induction of two well-studied transcription factors (TFs). We have applied the method to predict the metabolic impact of the induction of each of the approximately 180 MTB TFs using a previously generated and publically available expression data set.ConclusionsE-flux-MFC can be used to study global changes in MTB metabolites from gene expression data associated with environmental and genetic perturbations. The application of this method to a data set of MTB TF perturbations provides a resource for studying the large number of TFs whose functions remain unknown. Most TFs impact metabolites indirectly through the propagation of gene expression changes through the regulatory network rather than through their direct regulons. E-Flux-MFC is also applicable to any organism for which accurate metabolic models are available.

[1]  B. Gicquel,et al.  An essential role for phoP in Mycobacterium tuberculosis virulence , 2001, Molecular microbiology.

[2]  T. Myers,et al.  The Transcriptional Responses of Mycobacterium tuberculosis to Inhibitors of Metabolism , 2004, Journal of Biological Chemistry.

[3]  B. Abomoelak,et al.  A Novel Lipase Belonging to the Hormone-sensitive Lipase Family Induced under Starvation to Utilize Stored Triacylglycerol in Mycobacterium tuberculosis* , 2006, Journal of Biological Chemistry.

[4]  Gordon K Smyth,et al.  Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments , 2004, Statistical applications in genetics and molecular biology.

[5]  H. Maamar,et al.  Mycobacterium tuberculosis Uses Host Triacylglycerol to Accumulate Lipid Droplets and Acquires a Dormancy-Like Phenotype in Lipid-Loaded Macrophages , 2011, PLoS pathogens.

[6]  W. Bishai,et al.  Mechanisms of latency in Mycobacterium tuberculosis. , 1998, Trends in microbiology.

[7]  Peter M. Small,et al.  The W-Beijing Lineage of Mycobacterium tuberculosis Overproduces Triglycerides and Has the DosR Dormancy Regulon Constitutively Upregulated , 2007, Journal of bacteriology.

[8]  Marisa Klopper,et al.  Emergence and Spread of Extensively and Totally Drug-Resistant Tuberculosis, South Africa , 2013, Emerging Infectious Diseases.

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

[10]  David G. Russell,et al.  Tuberculosis: What We Don’t Know Can, and Does, Hurt Us , 2010, Science.

[11]  V. Kapur,et al.  Characterization of a two-component system, devR-devS, of Mycobacterium tuberculosis. , 2000, Tubercle and lung disease : the official journal of the International Union against Tuberculosis and Lung Disease.

[12]  Carolyn R Bertozzi,et al.  Cholesterol catabolism by Mycobacterium tuberculosis requires transcriptional and metabolic adaptations. , 2012, Chemistry & biology.

[13]  D. Sherman,et al.  Role of cholesterol in Mycobacterium tuberculosis infection. , 2009, Indian journal of experimental biology.

[14]  Brigitte Gicquel,et al.  The Virulence-associated Two-component PhoP-PhoR System Controls the Biosynthesis of Polyketide-derived Lipids in Mycobacterium tuberculosis* , 2006, Journal of Biological Chemistry.

[15]  William R. Jacobs,et al.  Complex lipid determines tissue-specific replication of Mycobacterium tuberculosis in mice , 1999, Nature.

[16]  Stefan Niemann,et al.  Mycobacterium tuberculosis wears what it eats. , 2010, Cell host & microbe.

[17]  David G. Russell,et al.  Who puts the tubercle in tuberculosis? , 2007, Nature Reviews Microbiology.

[18]  E. Chan,et al.  Multidrug-resistant and extensively drug-resistant tuberculosis: a review , 2008, Current opinion in infectious diseases.

[19]  C. Sohaskey,et al.  Nonreplicating persistence of mycobacterium tuberculosis. , 2001, Annual review of microbiology.

[20]  Markus J. Herrgård,et al.  Network-based prediction of human tissue-specific metabolism , 2008, Nature Biotechnology.

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

[22]  G. Migliori,et al.  First tuberculosis cases in Italy resistant to all tested drugs. , 2007, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[23]  Yves Van de Peer,et al.  The Mycobacterium tuberculosis regulatory network and hypoxia , 2013, Nature.

[24]  Hao Wu,et al.  MAANOVA: A Software Package for the Analysis of Spotted cDNA Microarray Experiments , 2003 .

[25]  E. Ruppin,et al.  Computational reconstruction of tissue-specific metabolic models: application to human liver metabolism , 2010, Molecular systems biology.

[26]  Jens Nielsen,et al.  Architecture of transcriptional regulatory circuits is knitted over the topology of bio-molecular interaction networks , 2008, BMC Systems Biology.

[27]  Amit Singh,et al.  Mycobacterium tuberculosis WhiB3 responds to O2 and nitric oxide via its [4Fe-4S] cluster and is essential for nutrient starvation survival , 2007, Proceedings of the National Academy of Sciences.

[28]  E. Ruppin,et al.  Predicting metabolic biomarkers of human inborn errors of metabolism , 2009, Molecular systems biology.

[29]  Melanie I. Stefan,et al.  Molecules for memory: modelling CaMKII , 2007, BMC Systems Biology.

[30]  Dirk Schnappinger,et al.  Regulation of the Mycobacterium tuberculosis hypoxic response gene encoding α-crystallin , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[31]  B. Palsson,et al.  Constraining the metabolic genotype–phenotype relationship using a phylogeny of in silico methods , 2012, Nature Reviews Microbiology.

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

[33]  Thomas Dick,et al.  Mycobacterium bovis BCG Response Regulator Essential for Hypoxic Dormancy , 2002, Journal of bacteriology.

[34]  B. Palsson,et al.  Regulation of gene expression in flux balance models of metabolism. , 2001, Journal of theoretical biology.

[35]  B. Palsson,et al.  Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110 , 1994, Applied and environmental microbiology.

[36]  Kyle J. Minch,et al.  Mapping and manipulating the Mycobacterium tuberculosis transcriptome using a transcription factor overexpression-derived regulatory network , 2014, Genome Biology.

[37]  Kellen L. Olszewski,et al.  Reconstruction and flux-balance analysis of the Plasmodium falciparum metabolic network , 2010, Molecular systems biology.

[38]  R. Heinrich,et al.  Metabolic Pathway Analysis: Basic Concepts and Scientific Applications in the Post‐genomic Era , 1999, Biotechnology progress.

[39]  Jason A. Papin,et al.  Metabolic pathways in the post-genome era. , 2003, Trends in biochemical sciences.

[40]  T. Hughes,et al.  Mapping pathways and phenotypes by systematic gene overexpression. , 2006, Molecular cell.

[41]  Jonathan M. Dreyfuss,et al.  Reconstruction and Validation of a Genome-Scale Metabolic Model for the Filamentous Fungus Neurospora crassa Using FARM , 2013, PLoS Comput. Biol..

[42]  Yang Liu,et al.  Transcriptional Adaptation of Mycobacterium tuberculosis within Macrophages , 2003, The Journal of experimental medicine.

[43]  R. Brosch,et al.  Evolutionary history of tuberculosis shaped by conserved mutations in the PhoPR virulence regulator , 2014, Proceedings of the National Academy of Sciences.

[44]  Gary K. Schoolnik,et al.  The Response of Mycobacterium Tuberculosis to Reactive Oxygen and Nitrogen Species , 2011, Front. Microbio..

[45]  I. Smith,et al.  Cholesterol metabolism increases the metabolic pool of propionate in Mycobacterium tuberculosis. , 2009, Biochemistry.

[46]  L. Gabbasova,et al.  Global tuberculosis report (2014) , 2014 .

[47]  E. Rubin,et al.  Genes required for mycobacterial growth defined by high density mutagenesis , 2003, Molecular microbiology.

[48]  Yukari C. Manabe,et al.  Latent Mycobacterium tuberculosis–persistence, patience, and winning by waiting , 2000, Nature Medicine.

[49]  O. Demin,et al.  The Edinburgh human metabolic network reconstruction and its functional analysis , 2007, Molecular systems biology.

[50]  F. Doyle,et al.  Dynamic flux balance analysis of diauxic growth in Escherichia coli. , 2002, Biophysical journal.

[51]  Andreas Hoppe,et al.  Antimalarial drug targets in Plasmodium falciparum predicted by stage-specific metabolic network analysis , 2010, BMC Systems Biology.

[52]  Bernhard O. Palsson,et al.  Investigating the metabolic capabilities of Mycobacterium tuberculosis H37Rv using the in silico strain iNJ661 and proposing alternative drug targets , 2007 .

[53]  Z. Udwadia,et al.  Totally drug-resistant tuberculosis (TDR-TB) in India: every dark cloud has a silver lining , 2012, Journal of Epidemiology & Community Health.

[54]  Irina Kolesnikova,et al.  A Thiolase of Mycobacterium tuberculosis Is Required for Virulence and Production of Androstenedione and Androstadienedione from Cholesterol , 2009, Infection and Immunity.

[55]  Carl Nathan,et al.  Taming tuberculosis: a challenge for science and society. , 2009, Cell host & microbe.

[56]  John Chan,et al.  Tuberculosis: Latency and Reactivation , 2001, Infection and Immunity.

[57]  A M Dannenberg,et al.  Immunopathogenesis of pulmonary tuberculosis. , 1993, Hospital practice.

[58]  Jaques Reifman,et al.  Development and analysis of an in vivo-compatible metabolic network of Mycobacterium tuberculosis , 2010, BMC Systems Biology.

[59]  Amit Singh,et al.  Mycobacterium tuberculosis WhiB3 Maintains Redox Homeostasis by Regulating Virulence Lipid Anabolism to Modulate Macrophage Response , 2009, PLoS pathogens.

[60]  S. Kaufmann,et al.  How can immunology contribute to the control of tuberculosis? , 2001, Nature Reviews Immunology.

[61]  David G. Russell,et al.  Mycobacterium and the coat of many lipids , 2002, The Journal of cell biology.

[62]  Martin Tompa,et al.  Rv3133c/dosR is a transcription factor that mediates the hypoxic response of Mycobacterium tuberculosis , 2003, Molecular microbiology.

[63]  Johnjoe McFadden,et al.  Carbon flux rerouting during Mycobacterium tuberculosis growth arrest , 2010, Molecular microbiology.

[64]  Marisa Klopper,et al.  Emergence and Spread of Extensively and Totally Drug-Resistant Tuberculosis, South Africa , 2013, Emerging infectious diseases.

[65]  C. Gille,et al.  HepatoNet1: a comprehensive metabolic reconstruction of the human hepatocyte for the analysis of liver physiology , 2010, Molecular systems biology.

[66]  B. Palsson,et al.  Toward Metabolic Phenomics: Analysis of Genomic Data Using Flux Balances , 1999, Biotechnology progress.

[67]  Sharon L. Kendall,et al.  Cholesterol utilization in mycobacteria is controlled by two TetR-type transcriptional regulators: kstR and kstR2 , 2010, Microbiology.

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

[69]  James C. Sacchettini,et al.  Persistence of Mycobacterium tuberculosis in macrophages and mice requires the glyoxylate shunt enzyme isocitrate lyase , 2000, Nature.

[70]  Parissa Farnia,et al.  Emergence of new forms of totally drug-resistant tuberculosis bacilli: super extensively drug-resistant tuberculosis or totally drug-resistant strains in iran. , 2009, Chest.

[71]  David G. Russell,et al.  Intracellular Mycobacterium tuberculosis Exploits Host-derived Fatty Acids to Limit Metabolic Stress* , 2013, The Journal of Biological Chemistry.

[72]  Fong-Fu Hsu,et al.  aprABC: a Mycobacterium tuberculosis complex‐specific locus that modulates pH‐driven adaptation to the macrophage phagosome , 2011, Molecular microbiology.

[73]  Christopher J Petzold,et al.  Lipidomics reveals control of Mycobacterium tuberculosis virulence lipids via metabolic coupling , 2007, Proceedings of the National Academy of Sciences.

[74]  B. Abomoelak,et al.  Induction of a Novel Class of Diacylglycerol Acyltransferases and Triacylglycerol Accumulation in Mycobacterium tuberculosis as It Goes into a Dormancy-Like State in Culture , 2004, Journal of bacteriology.

[75]  Scott L. Zeger,et al.  The Analysis of Gene Expression Data: Methods and Software , 2013 .

[76]  Tige R. Rustad,et al.  The Enduring Hypoxic Response of Mycobacterium tuberculosis , 2008, PloS one.

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

[78]  D. Sherman,et al.  The secret lives of the pathogenic mycobacteria. , 2003, Annual review of microbiology.

[79]  Ernesto S. Nakayasu,et al.  Model-driven multi-omic data analysis elucidates metabolic immunomodulators of macrophage activation , 2012, Molecular systems biology.

[80]  Beatriz Galán,et al.  Characterization of the KstR-dependent promoter of the gene for the first step of the cholesterol degradative pathway in Mycobacterium smegmatis. , 2011, Microbiology.

[81]  Robert Urbanczik,et al.  The geometry of the flux cone of a metabolic network. , 2005, Biophysical journal.

[82]  S. Klamt,et al.  GSMN-TB: a web-based genome-scale network model of Mycobacterium tuberculosis metabolism , 2007, Genome Biology.

[83]  D. Schnappinger,et al.  Biosynthesis and Recycling of Nicotinamide Cofactors in Mycobacterium tuberculosis , 2008, Journal of Biological Chemistry.

[84]  B. Palsson,et al.  Saccharomyces cerevisiae phenotypes can be predicted by using constraint-based analysis of a genome-scale reconstructed metabolic network , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[85]  Bernhard O. Palsson,et al.  Constraint-based analysis of metabolic capacity of Salmonella typhimurium during host-pathogen interaction , 2009, BMC Systems Biology.

[86]  David G. Russell,et al.  Mycobacterium tuberculosis: here today, and here tomorrow , 2001, Nature Reviews Molecular Cell Biology.

[87]  Christopher M. Sassetti,et al.  igr Genes and Mycobacterium tuberculosis Cholesterol Metabolism , 2009, Journal of bacteriology.

[88]  Bonnie Berger,et al.  MetaMerge: scaling up genome-scale metabolic reconstructions with application to Mycobacterium tuberculosis , 2012, Genome Biology.

[89]  Stephen S Fong,et al.  Metabolic gene–deletion strains of Escherichia coli evolve to computationally predicted growth phenotypes , 2004, Nature Genetics.

[90]  Tom M. Conrad,et al.  Omic data from evolved E. coli are consistent with computed optimal growth from genome-scale models , 2010, Molecular systems biology.

[91]  Monica L. Mo,et al.  Global reconstruction of the human metabolic network based on genomic and bibliomic data , 2007, Proceedings of the National Academy of Sciences.

[92]  C. Sassetti,et al.  Metabolic Regulation of Mycobacterial Growth and Antibiotic Sensitivity , 2011, PLoS biology.

[93]  Neema Jamshidi,et al.  Updating and curating metabolic pathways of TB. , 2013, Tuberculosis.

[94]  C. Lima,et al.  Crystal Structure of PapA5, a Phthiocerol Dimycocerosyl Transferase from Mycobacterium tuberculosis* , 2004, Journal of Biological Chemistry.

[95]  B. Palsson,et al.  In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data , 2001, Nature Biotechnology.

[96]  B. Palsson,et al.  Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions , 2010, Molecular systems biology.

[97]  Sara B. Collins,et al.  Temporal Expression-based Analysis of Metabolism , 2012, PLoS Comput. Biol..

[98]  V. Dubey,et al.  Disruption of msl3 abolishes the synthesis of mycolipanoic and mycolipenic acids required for polyacyltrehalose synthesis in Mycobacterium tuberculosis H37Rv and causes cell aggregation , 2002, Molecular microbiology.

[99]  N. Price,et al.  Probabilistic integrative modeling of genome-scale metabolic and regulatory networks in Escherichia coli and Mycobacterium tuberculosis , 2010, Proceedings of the National Academy of Sciences.

[100]  Irina Kolesnikova,et al.  The Mycobacterium tuberculosis PhoPR two‐component system regulates genes essential for virulence and complex lipid biosynthesis , 2006, Molecular microbiology.

[101]  A. Steyn,et al.  Mycobacterium tuberculosis DosS is a redox sensor and DosT is a hypoxia sensor , 2007, Proceedings of the National Academy of Sciences.

[102]  J. Nielsen,et al.  Uncovering transcriptional regulation of metabolism by using metabolic network topology. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[103]  James E Gomez,et al.  M. tuberculosis persistence, latency, and drug tolerance. , 2004, Tuberculosis.

[104]  Andrzej M. Kierzek,et al.  Systems-Based Approaches to Probing Metabolic Variation within the Mycobacterium tuberculosis Complex , 2013, PloS one.

[105]  D. Sherman,et al.  Regulation of the Mycobacterium tuberculosis hypoxic response gene encoding alpha -crystallin. , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[107]  Markus R. Wenk,et al.  Triacylglycerol Utilization Is Required for Regrowth of In Vitro Hypoxic Nonreplicating Mycobacterium bovis Bacillus Calmette-Guerin , 2009, Journal of bacteriology.

[108]  H. Sprecher,et al.  The Mycobacterium tuberculosis pks2 Gene Encodes the Synthase for the Hepta- and Octamethyl-branched Fatty Acids Required for Sulfolipid Synthesis* , 2001, The Journal of Biological Chemistry.

[109]  Ed Reznik,et al.  Flux Imbalance Analysis and the Sensitivity of Cellular Growth to Changes in Metabolite Pools , 2013, PLoS Comput. Biol..

[110]  Andrzej M. Kierzek,et al.  Differential Producibility Analysis (DPA) of Transcriptomic Data with Metabolic Networks: Deconstructing the Metabolic Response of M. tuberculosis , 2011, PLoS Comput. Biol..

[111]  Bernhard O. Palsson,et al.  Metabolic flux balance analysis and the in silico analysis of Escherichia coli K-12 gene deletions , 2000, BMC Bioinformatics.

[112]  Ben Sidders,et al.  A highly conserved transcriptional repressor controls a large regulon involved in lipid degradation in Mycobacterium smegmatis and Mycobacterium tuberculosis , 2007, Molecular microbiology.