Oxygen Availability and Metabolic Dynamics During Mycobacterium tuberculosis Latency

Objective: In vitro models of Mycobacterium tuberculosis (Mtb) nonreplicating persistence (NRP) suggest the rate of oxygen (O2) depletion is a significant determinant in persistence. However, few studies have characterized the metabolic association between slow and rapid O2 depletion rates and successful versus failed persistence. Methods: We developed a theoretical model of Mtb metabolic adaptation that includes O2 driven genetic modulation of enzymes in the tricarboxylic acid cycle, energy, and redox recycling pathways. We conducted an in silico study of Mtb adaptation, and investigated the metabolic dynamics that enable persistence during slow versus rapid anaerobiosis. Results: Consistent with in vitro studies, during rapid anaerobiosis the most significant enzymatic changes occurred during the active growth period while the majority of the slow anaerobic system's adaptive response occurred during the first phase of aerobic shiftdown (NRP1). The characteristic response of the two conditions differed, with the slow anaerobic system exhibiting positive adaptation response, upregulating six enzymes during NRP1, while the rapid system exhibited a negative adaptation response, decreasing the levels of seven enzymes during active growth. Conclusion: Results of our study illustrate the intricate metabolic balance Mtb achieves during adaptation to anaerobic conditions, and how failure in redox recycling correlates to failure in persistence. Significance: We have provided the first theoretical description of the genetic and metabolic death profile for Mtb during anaerobic growth. Insight into Mtbs dynamic adaptation response provides a unique systems framework for the development of targeted therapies to reduce the likelihood of mycobacterial persistence and the related incidence of latent tuberculosis infection.

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

[2]  Sabine Ehrt,et al.  Expression profiling of host pathogen interactions: how Mycobacterium tuberculosis and the macrophage adapt to one another. , 2006, Microbes and infection.

[3]  G. Schoolnik,et al.  Mycobacterium tuberculosis gene expression during adaptation to stationary phase and low-oxygen dormancy. , 2004, Tuberculosis.

[4]  E E May,et al.  BioXyce: an engineering platform for the study of cellular systems. , 2009, IET systems biology.

[5]  Steven J. Plimpton,et al.  Oxygen Modulates the Effectiveness of Granuloma Mediated Host Response to Mycobacterium tuberculosis: A Multiscale Computational Biology Approach , 2016, Front. Cell. Infect. Microbiol..

[6]  Nicola Mulder,et al.  A mathematical representation of the development of Mycobacterium tuberculosis active, latent and dormant stages. , 2012, Journal of theoretical biology.

[7]  Tudor I. Oprea,et al.  A systems chemical biology study of malate synthase and isocitrate lyase inhibition in Mycobacterium tuberculosis during active and NRP growth , 2013, Comput. Biol. Chem..

[8]  Eric R. Keiter,et al.  Xyce Parallel Electronic Simulator : reference guide, version 2.0. , 2004 .

[9]  L. Vogler.  Lehninger Principles Of Biochemistry 5th Edition , 2016 .

[10]  Elebeoba May,et al.  Circuit-Based Models of Biomolecular System Dynamics , 2011 .

[11]  Michael S. Eldred,et al.  DAKOTA : a multilevel parallel object-oriented framework for design optimization, parameter estimation, uncertainty quantification, and sensitivity analysis. Version 5.0, user's reference manual. , 2010 .

[12]  H. Sauro.  Enzyme Kinetics for Systems Biology , 2012 .

[13]  Nicholas Fisher,et al.  Chapter 17 Type II NADH: quinone oxidoreductases of Plasmodium falciparum and Mycobacterium tuberculosis kinetic and high-throughput assays. , 2009, Methods in enzymology.

[14]  E. A. Dawes,et al.  Biochemical calculations, 2nd Edition By Irwin H. Segel. John Wiley & Sons, London, New York, Sydney, Toronto. xvii + 441 pp. £5.35 , 1976 .

[15]  Harvey Rubin,et al.  The role of RelMtb-mediated adaptation to stationary phase in long-term persistence of Mycobacterium tuberculosis in mice , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Shin-Il Kim,et al.  Mycobacterial granulomas: keys to a long-lasting host-pathogen relationship. , 2004, Clinical immunology.

[17]  T. Marshall,et al.  Common angiotensin receptor blockers may directly modulate the immune system via VDR, PPAR and CCR2b , 2006, Theoretical Biology and Medical Modelling.

[18]  E. May,et al.  Addressing biological circuit simulation accuracy: Reachability for parameter identification and initial conditions , 2007, 2007 IEEE/NIH Life Science Systems and Applications Workshop.

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

[20]  Rick Lyons,et al.  The temporal expression profile of Mycobacterium tuberculosis infection in mice. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

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

[22]  Sophia Lefantzi,et al.  DAKOTA : a multilevel parallel object-oriented framework for design optimization, parameter estimation, uncertainty quantification, and sensitivity analysis. , 2011 .

[23]  Thomas Dick,et al.  The protonmotive force is required for maintaining ATP homeostasis and viability of hypoxic, nonreplicating Mycobacterium tuberculosis , 2008, Proceedings of the National Academy of Sciences.

[24]  Antje Chang,et al.  BRENDA, the enzyme information system in 2011 , 2010, Nucleic Acids Res..

[25]  James C Sacchettini,et al.  Biochemical and Structural Studies of Malate Synthase fromMycobacterium tuberculosis * , 2002, The Journal of Biological Chemistry.

[26]  V. Singh,et al.  Kinetic modeling of tricarboxylic acid cycle and glyoxylate bypass in Mycobacterium tuberculosis, and its application to assessment of drug targets , 2006, Theoretical Biology and Medical Modelling.

[27]  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.

[28]  L G Wayne,et al.  Glyoxylate metabolism and adaptation of Mycobacterium tuberculosis to survival under anaerobic conditions , 1982, Infection and immunity.

[29]  D. Kirschner,et al.  A methodology for performing global uncertainty and sensitivity analysis in systems biology. , 2008, Journal of theoretical biology.

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

[31]  B. Abomoelak,et al.  A Novel In Vitro Multiple-Stress Dormancy Model for Mycobacterium tuberculosis Generates a Lipid-Loaded, Drug-Tolerant, Dormant Pathogen , 2009, PloS one.

[32]  M. Shiloh,et al.  Mycobacterium tuberculosis senses host-derived carbon monoxide during macrophage infection. , 2008, Cell host & microbe.

[33]  Eric R. Keiter,et al.  Xyce parallel electronic simulator. , 2010 .

[34]  Nicola Mulder,et al.  Understanding TB latency using computational and dynamic modelling procedures. , 2013, Infection, genetics and evolution : journal of molecular epidemiology and evolutionary genetics in infectious diseases.

[35]  L. Wayne,et al.  An in vitro model for sequential study of shiftdown of Mycobacterium tuberculosis through two stages of nonreplicating persistence , 1996, Infection and immunity.