Metabolic modeling and experimental studies on carbon and energy metabolism of Mycobacterium tuberculosis

Tuberculosis is an ancient human disease affecting millions of people every year. The common causative agent, Mycobacterium tuberculosis, is a peculiar pathogen that can persist in the host tissues for decades. Many of the key attributes of M. tuberculosis can be traced to its metabolism. Therefore, novel insights into the metabolism of this bacterium could lead to the identification of new drug targets as well as better methods to understand the mechanism of action of anti-tuberculosis drugs. Analysis of the metabolic network of this bacterium using a genome-scale model offers an unconventional approach to understanding utilization of metabolic pathways of this pathogen during growth and non-replicating persistence. Metabolic adaptation of a pathogen is central to its life style in vivo. To understand the capacity and limitations of M. tuberculosis, we developed an updated and thermodynamically curated genome scale model – iOS783. This model has been verified with respect to the observed phenotypes for targeted gene deletion studies reported previously. Evaluation of iOS783 with thermodynamics based flux balance analysis (TFBA) improved overall reliability of the model results and led to identification of possible bottleneck reactions. In addition, several new hypotheses were generated with iOS783, involving the role of isocitrate lyase enzymes, propionate toxicity, and energy metabolism of this pathogen. Analysis of iOS783 for conditional gene essentiality emphasizes the importance of the carbon source available to the pathogen. Several genes show a conditional impact on biomass production rate of the model, depending on the availability of carbohydrates or fatty acids as the carbon substrate. Conditional essentiality is even more pronounced when ATP replenishment capacity is used as a measure of a pathogen’s ability to survive under stress in a non-growing state. As expected, ATP synthase plays a major role in estimations of ATP replenishment capacity based on iOS783. Curiously, computational analysis of the metabolic network with ATP synthase activity blocked indicates that the damaged network can use substrate-level phosphorylation reactions to sustain ATP production to fulfill non-growth-associated ATP requirements when glucose is available. In order to establish the dependence of M. tuberculosis on ATP synthase, we attempted to construct a deletion of the chromosomal atp operon. Our results show that the native copy of these genes can be deleted only when a second set of genes are introduced elsewhere in the genome, indicating essentiality of the ATP synthase for growth. Although our results indicate that ATP synthase is essential for growth, we observed carbon source specific efficacy of an ATP synthase inhibitor, bedaquiline, on M. tuberculosis. These experimental results are consistent with ATP turnover estimations based on modeling results. Furthermore, analysis of iOS783 suggests a novel role in energy metabolism for recently described glycine dehydrogenase (GDH) on energy metabolism. Deletion of ald gene responsible for GDH activity reduces the ATP replenishment capacity of iOS783, especially during metabolism of fatty acid derivative carbon sources. Consistent with these modeling results, wetlab experiments designed to test this hypothesis show that bedaquiline has a stronger bactericidal effect on an M. tuberculosis ∆ald mutant, compared to the wild-type parental strain, when cells are cultured in medium containing acetate as the carbon source. Metabolic modelling is already well established as an invaluable tool for engineering of strains with desired characteristics. For pathogenic species, however, the application of this approach has so far been quite limited. The work presented in this thesis exemplifies an interdisciplinary approach combining computational modeling and experimental studies of M. tuberculosis metabolism. First, iOS783 was utilized to interpret existing experimental results and to generate novel hypotheses, thereby providing the impetus for new wetlab experiments. Conversely, results from iOS783-inspired wetlab experiments were used to further refine and expand the iOS783 model. This cyclic exchange between computational modeling and wetlab experimentation illustrates how a close interdisciplinary partnership can create a positive feedback loop that enriches both disciplines.

Quantitative Mass Spectrometry Reveals Plasticity of Metabolic Networks in Mycobacterium smegmatis

Florence Armand, Diego Chiappe, Tarun Chopra, Romain Hamelin, John McKinney, Marc Moniatte

Mycobacterium tuberculosis has a remarkable ability to persist within the human host as a clinically inapparent or chronically active infection. Fatty acids are thought to be an important carbon source used by the bacteria during long term infection. Catabolism of fatty acids requires reprogramming of metabolic networks, and enzymes central to this reprogramming have been targeted for drug discovery. Mycobacterium smegmatis, a nonpathogenic relative of M. tuberculosis, is often used as a model system because of the similarity of basic cellular processes in these two species. Here, we take a quantitative proteomics-based approach to achieve a global view of how the M. smegmatis metabolic network adjusts to utilization of fatty acids as a carbon source. Two-dimensional liquid chromatography and mass spectrometry of isotopically labeled proteins identified a total of 3,067 proteins with high confidence. This number corresponds to 44% of the predicted M. smegmatis proteome and includes most of the predicted metabolic enzymes. Compared with glucose-grown cells, 162 proteins showed differential abundance in acetate- or propionate-grown cells. Among these, acetate-grown cells showed a higher abundance of proteins that could constitute a functional glycerate pathway. Gene inactivation experiments confirmed that both the glyoxylate shunt and the glycerate pathway are operational in M. smegmatis. In addition to proteins with annotated functions, we demonstrate carbon source-dependent differential abundance of proteins that have not been functionally characterized. These proteins might play as-yet-unidentified roles in mycobacterial carbon metabolism. This study reveals several novel features of carbon assimilation in M. smegmatis, which suggests significant functional plasticity of metabolic networks in this organism.

Amer Soc Biochemistry Molecular Biology Inc2014

Differential expression of iron-, carbon-, and oxygen-responsive mycobacterial genes in the lungs of chronically infected mice and tuberculosis patients

John McKinney

Pathogenetic processes that facilitate the entry, replication, and persistence of Mycobacterium tuberculosis (MTB) in the mammalian host likely include the regulated expression of specific sets of genes at different stages of infection. Identification of genes that are differentially expressed in vivo would provide insights into host-pathogen interactions in tuberculosis (TB); this approach might be particularly valuable for the study of human TB, where experimental opportunities are limited. In this study, the levels of selected MTB mRNAs were quantified in vitro in axenic culture, in vivo in the lungs of mice, and in lung specimens obtained from TB patients with active disease. We report the differential expression of MTB mRNAs associated with iron limitation, alternative carbon metabolism, and cellular hypoxia, conditions that are thought to exist within the granulomatous lesions of TB, in the lungs of wild-type C57BL/6 mice as compared with bacteria grown in vitro. Analysis of the same set of mRNAs in lung specimens obtained from TB patients revealed differences in MTB gene expression in humans as compared with mice.

2003

Single-cell studies of Mycobacterium smegmatis cell cycle using time-lapse fluorescence microscopy

Joëlle Xiao Yuan Ven

Cell division is one of the most primordial cellular processes and is essential for life propagation. It is composed of crucial events, namely genome duplication, chromosomes segregation and membrane separation. Failure of these processes will have dramatic impacts on the development of any organisms, including cell death, polyploidy, or cell cycle arrest. In the case of bacterial pathogens, blocking cell division efficiently using targeted drugs is of interest to cure bacterial infections. However, despite its ancestral origin, the mechanisms of cell division have diverged amongst organisms. Given its crucial role, as well as the potential therapeutic applications, cell division is an intensive field of research and many aspects are yet to be fully understood. Amongst bacteria, the phyla Proteobacteria and Firmicutes have been intensively studied using their respective model organisms Escherichia coli and Bacillus subtilis. However, there is a paucity of research on cell division in Actinobacteria, the phylum that contains the worldâs leading pathogen, Mycobacterium tuberculosis. In particular, the molecular players involved in the process of chromosome segregation and the timing of cell division have not been fully investigated. In the well-studied bacteria E. coli and B. subtilis the Min system determines the division site localization, while the nucleoid occlusion system (Noc or SlmA, respectively) determines the timing of division. In mycobacteria, however, no homologs for the Min or Noc systems have been identified thus far. However, the partitioning (par) genes are currently thought to be responsible for symmetric chromosome segregation. In addition, FtsK has been identified as a homolog of the B. subtilis protein SpoIIIE, which is responsible for rescuing DNA in the forming septum and pulling it in the forespore during starvation-induced sporulation. Overall, the role and dynamics of DNA, where FtsK might play a role, during mycobacterial cell division, still remain to be uncovered. Current methods to study DNA localization using snapshot microscopy typically utilise DNA dyes, such as DAPI, Hoechst, or SYTO. However, these compounds impair cell growth, preventing the dynamic study of mycobacterial division. To circumvent this shortcoming, I constructed a Mycobacterium smegmatis reporter strain expressing a Dendra2 fusion protein in-frame to the histone-like protein Hlp, and compared it to conventional DNA dyes. Additionally, I characterized FtsK and, for the first time, demonstrated its essentiality in M. smegmatis. As FtsK is essential I adapted a dual repression system to knock down FtsK, which proved to be a synthetic rescue with the par genes deletion. Using a custom-made microfluidic device and long-term time-lapse microscopy, the distribution and dynamics of both Hlp and FtsK fluorescent reporter proteins were followed at the single-cell level. This approach allowed the spatial and temporal characterization of these proteins with respect to other septal genes in wild-type cells as well as in Îpar mutants. Finally, I propose a new model where the par genes are important for chromosome orientation and FtsK is important for chromosome segregation. Altogether, this thesis represents a step forward towards a better understanding of the mechanisms underlying cell division in mycobacteria and in the long-term, a potential source of novel treatments to control the division of M. tuberculosis, for instance by targeting FtsK.