Differential producibility analysis (DPA) of transcriptomic data with metabolic networks: deconstructing the metabolic response of M. tuberculosis

Differential producibility analysis reveals drug-associated carbon and nitrogen metabolite expressions in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) is one of the world's successful pathogens that flexibly adapts its metabolic nature during infection of the host, and in response to drugs. Here we used genome scale metabolic modelling coupled with differential producibility analysis (DPA) to translate RNA-seq datasets into metabolite signals and identified drug-associated metabolic response profiles. We tested four tuberculosis (TB) drugs bedaquiline (BDQ), isoniazid (INH), rifampicin (RIF), and clarithromycin (CLA); conducted RNA-seq experiments of Mtb exposed to the individual drugs at subinhibitory concentrations, followed by DPA of gene expression data to map up and downregulated metabolites. Here we highlight those metabolic pathways that were flexibly used by Mtb to tolerate stress generated upon drug exposure. BDQ and INH upregulated maximum number of central carbon metabolites in glycolysis, pentose phosphate pathway and tri-carboxylic acid cycle with concomitant downregulation of lipid and amino acid metabolite classes. Oxaloacetate was significantly upregulated in all four drug-treated Mtb cells highlighting it as an important metabolite in Mtb's metabolism. Amino acid metabolism was selectively induced by different drugs. We have enhanced our knowledge on Mtb's carbon and nitrogen metabolic adaptations in the presence of drugs and identify metabolic nodes for therapeutic development against TB. Our work also provides DPA omics platform to interrogate RNA-seq datasets of any organism that can be reconstructed as a genome scale metabolic network.

Dissecting Host-Pathogen Interactions in TB Using Systems-Based Omic Approaches

Tuberculosis (TB) is a devastating infectious disease that kills over a million people every year. There is an increasing burden of multi drug resistance (MDR) and extensively drug resistance (XDR) TB. New and improved therapies are urgently needed to overcome the limitations of current treatment. The causative agent, Mycobacterium tuberculosis (Mtb) is one of the most successful pathogens that can manipulate host cell environment for adaptation, evading immune defences, virulence, and pathogenesis of TB infection. Host-pathogen interaction is important to establish infection and it involves a complex set of processes. Metabolic cross talk between the host and pathogen is a facet of TB infection and has been an important topic of research where there is growing interest in developing therapies and drugs that target these interactions and metabolism of the pathogen in the host. Mtb scavenges multiple nutrient sources from the host and has adapted its metabolism to survive in the intracellular niche. Advancements in systems-based omic technologies have been successful to unravel host-pathogen interactions in TB. In this review we discuss the application and usefulness of omics in TB research that provides promising interventions for developing anti-TB therapies.

GSMN-ML- a genome scale metabolic network reconstruction of the obligate human pathogen Mycobacterium leprae

Leprosy, caused by Mycobacterium leprae, has plagued humanity for thousands of years and continues to cause morbidity, disability and stigmatization in two to three million people today. Although effective treatment is available, the disease incidence has remained approximately constant for decades so new approaches, such as vaccine or new drugs, are urgently needed for control. Research is however hampered by the pathogen's obligate intracellular lifestyle and the fact that it has never been grown in vitro. Consequently, despite the availability of its complete genome sequence, fundamental questions regarding the biology of the pathogen, such as its metabolism, remain largely unexplored. In order to explore the metabolism of the leprosy bacillus with a long-term aim of developing a medium to grow the pathogen in vitro, we reconstructed an in silico genome scale metabolic model of the bacillus, GSMN-ML. The model was used to explore the growth and biomass production capabilities of the pathogen with a range of nutrient sources, such as amino acids, glucose, glycerol and metabolic intermediates. We also used the model to analyze RNA-seq data from M. leprae grown in mouse foot pads, and performed Differential Producibility Analysis to identify metabolic pathways that appear to be active during intracellular growth of the pathogen, which included pathways for central carbon metabolism, co-factor, lipids, amino acids, nucleotides and cell wall synthesis. The GSMN-ML model is thereby a useful in silico tool that can be used to explore the metabolism of the leprosy bacillus, analyze functional genomic experimental data, generate predictions of nutrients required for growth of the bacillus in vitro and identify novel drug targets.

Network Analyses in Plant Pathogens

Even in the age of big data in Biology, studying the connections between the biological processes and the molecular mechanisms behind them is a challenging task. Systems biology arose as a transversal discipline between biology, chemistry, computer science, mathematics, and physics to facilitate the elucidation of such connections. A scenario, where the application of systems biology constitutes a very powerful tool, is the study of interactions between hosts and pathogens using network approaches. Interactions between pathogenic bacteria and their hosts, both in agricultural and human health contexts are of great interest to researchers worldwide. Large amounts of data have been generated in the last few years within this area of research. However, studies have been relatively limited to simple interactions. This has left great amounts of data that remain to be utilized. Here, we review the main techniques in network analysis and their complementary experimental assays used to investigate bacterial-plant interactions. Other host-pathogen interactions are presented in those cases where few or no examples of plant pathogens exist. Furthermore, we present key results that have been obtained with these techniques and how these can help in the design of new strategies to control bacterial pathogens. The review comprises metabolic simulation, protein-protein interactions, regulatory control of gene expression, host-pathogen modeling, and genome evolution in bacteria. The aim of this review is to offer scientists working on plant-pathogen interactions basic concepts around network biology, as well as an array of techniques that will be useful for a better and more complete interpretation of their data.

Metabolic modeling predicts metabolite changes in Mycobacterium tuberculosis

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Mycobacterium tuberculosis metabolism

Metabolism underpins the physiology and pathogenesis of Mycobacterium tuberculosis. However, although experimental mycobacteriology has provided key insights into the metabolic pathways that are essential for survival and pathogenesis, determining the metabolic status of bacilli during different stages of infection and in different cellular compartments remains challenging. Recent advances-in particular, the development of systems biology tools such as metabolomics-have enabled key insights into the biochemical state of M. tuberculosis in experimental models of infection. In addition, their use to elucidate mechanisms of action of new and existing antituberculosis drugs is critical for the development of improved interventions to counter tuberculosis. This review provides a broad summary of mycobacterial metabolism, highlighting the adaptation of M. tuberculosis as specialist human pathogen, and discusses recent insights into the strategies used by the host and infecting bacillus to influence the outcomes of the host-pathogen interaction through modulation of metabolic functions.

Systems-level modeling of mycobacterial metabolism for the identification of new (multi-)drug targets

Systems-level metabolic network reconstructions and the derived constraint-based (CB) mathematical models are efficient tools to explore bacterial metabolism. Approximately one-fourth of the Mycobacterium tuberculosis (Mtb) genome contains genes that encode proteins directly involved in its metabolism. These represent potential drug targets that can be systematically probed with CB models through the prediction of genes essential (or the combination thereof) for the pathogen to grow. However, gene essentiality depends on the growth conditions and, so far, no in vitro model precisely mimics the host at the different stages of mycobacterial infection, limiting model predictions. These limitations can be circumvented by combining expression data from in vivo samples with a validated CB model, creating an accurate description of pathogen metabolism in the host. To this end, we present here a thoroughly curated and extended genome-scale CB metabolic model of Mtb quantitatively validated using 13C measurements. We describe some of the efforts made in integrating CB models and high-throughput data to generate condition specific models, and we will discuss challenges ahead. This knowledge and the framework herein presented will enable to identify potential new drug targets, and will foster the development of optimal therapeutic strategies.

Promise and reality in the expanding field of network interaction analysis: metabolic networks

In the last few decades, metabolic networks revealed their capabilities as powerful tools to analyze the cellular metabolism. Many research fields (eg, metabolic engineering, diagnostic medicine, pharmacology, biochemistry, biology and physiology) improved the understanding of the cell combining experimental assays and metabolic network-based computations. This process led to the rise of the “systems biology” approach, where the theory meets experiments and where two complementary perspectives cooperate in the study of biological phenomena. Here, the reconstruction of metabolic networks is presented, along with established and new algorithms to improve the description of cellular metabolism. Then, advantages and limitations of modeling algorithms and network reconstruction are discussed.

Mycobacterial Lipidomics

Lipidomics is a distinct subspecialty of metabolomics concerned with hydrophobic molecules that organize into membranes. Most of the lipid classes present in Mycobacterium tuberculosis are found only in Actinobacteria and show extreme structural diversity. This article highlights the conceptual basis and the practical challenges associated with the mass spectrometry–based lipidomic study of M. tuberculosis to solve basic questions about the virulence of this lipid-laden organism.

Deciphering the response of Mycobacterium smegmatis to nitrogen stress using bipartite active modules

BACKGROUND

The ability to adapt to environments with fluctuating nutrient availability is vital for bacterial survival. Although essential for growth, few nitrogen metabolism genes have been identified or fully characterised in mycobacteria and nitrogen stress survival mechanisms are unknown.

RESULTS

A global transcriptional analysis of the mycobacterial response to nitrogen stress, showed a significant change in the differential expression of 16% of the Mycobacterium smegmatis genome. Gene expression changes were mapped onto the metabolic network using Active Modules for Bipartite Networks (AMBIENT) to identify metabolic pathways showing coordinated transcriptional responses to the stress. AMBIENT revealed several key features of the metabolic response not identified by KEGG enrichment alone. Down regulated reactions were associated with the general reduction in cellular metabolism as a consequence of reduced growth rate. Up-regulated modules highlighted metabolic changes in nitrogen assimilation and scavenging, as well as reactions involved in hydrogen peroxide metabolism, carbon scavenging and energy generation.

CONCLUSIONS

Application of an Active Modules algorithm to transcriptomic data identified key metabolic reactions and pathways altered in response to nitrogen stress, which are central to survival under nitrogen limiting environments.

Parallel exploitation of diverse host nutrients enhances Salmonella virulence

Pathogen access to host nutrients in infected tissues is fundamental for pathogen growth and virulence, disease progression, and infection control. However, our understanding of this crucial process is still rather limited because of experimental and conceptual challenges. Here, we used proteomics, microbial genetics, competitive infections, and computational approaches to obtain a comprehensive overview of Salmonella nutrition and growth in a mouse typhoid fever model. The data revealed that Salmonella accessed an unexpectedly diverse set of at least 31 different host nutrients in infected tissues but the individual nutrients were available in only scarce amounts. Salmonella adapted to this situation by expressing versatile catabolic pathways to simultaneously exploit multiple host nutrients. A genome-scale computational model of Salmonella in vivo metabolism based on these data was fully consistent with independent large-scale experimental data on Salmonella enzyme quantities, and correctly predicted 92% of 738 reported experimental mutant virulence phenotypes, suggesting that our analysis provided a comprehensive overview of host nutrient supply, Salmonella metabolism, and Salmonella growth during infection. Comparison of metabolic networks of other pathogens suggested that complex host/pathogen nutritional interfaces are a common feature underlying many infectious diseases.

Modeling phenotypic metabolic adaptations of Mycobacterium tuberculosis H37Rv under hypoxia

The ability to adapt to different conditions is key for Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), to successfully infect human hosts. Adaptations allow the organism to evade the host immune responses during acute infections and persist for an extended period of time during the latent infectious stage. In latently infected individuals, estimated to include one-third of the human population, the organism exists in a variety of metabolic states, which impedes the development of a simple strategy for controlling or eradicating this disease. Direct knowledge of the metabolic states of M. tuberculosis in patients would aid in the management of the disease as well as in forming the basis for developing new drugs and designing more efficacious drug cocktails. Here, we propose an in silico approach to create state-specific models based on readily available gene expression data. The coupling of differential gene expression data with a metabolic network model allowed us to characterize the metabolic adaptations of M. tuberculosis H37Rv to hypoxia. Given the microarray data for the alterations in gene expression, our model predicted reduced oxygen uptake, ATP production changes, and a global change from an oxidative to a reductive tricarboxylic acid (TCA) program. Alterations in the biomass composition indicated an increase in the cell wall metabolites required for cell-wall growth, as well as heightened accumulation of triacylglycerol in preparation for a low-nutrient, low metabolic activity life style. In contrast, the gene expression program in the deletion mutant of dosR, which encodes the immediate hypoxic response regulator, failed to adapt to low-oxygen stress. Our predictions were compatible with recent experimental observations of M. tuberculosis activity under hypoxic and anaerobic conditions. Importantly, alterations in the flow and accumulation of a particular metabolite were not necessarily directly linked to differential gene expression of the enzymes catalyzing the related metabolic reactions.

Mycobacterium tuberculosis latently infects one-third of the human population and is responsible for millions of deaths worldwide every year. The ability of the pathogen to persist in the human population stems from its capacity to adapt to host-induced stresses and adjust its metabolism to different host environments. We have developed a novel model to interpret M. tuberculosis H37Rv metabolic adjustment by combining gene transcription data with a genome-scale metabolic network model. Using our model, we were able to identify the changes in the metabolic program associated with hypoxia, predict phenotypic change, and determine the critical metabolic enzymes and pathways that are required for pathogen survival. In particular, we predicted the switch in the tricarboxylic acid cycle from an oxidative to a reductive path. The altered importance of different metabolites and pathways under hypoxic conditions may provide guidance for designing novel, adjuvant drug therapies for clearing persistent and latent infections.

Network analysis has diverse roles in drug discovery

Computational biologists use network analysis to uncover relationships between various data types of interest for drug discovery. For example, signalling and metabolic pathways are commonly used to understand disease states and drug mechanisms. However, several other flavours of network analysis techniques are also applicable in a drug discovery context. Recent advances include networks that encompass relationships between diseases, molecular mechanisms and gene targets. Even social networks that mirror interactions within the scientific community are helping to foster collaborations and novel research. We review how these different types of network analysis approaches facilitate drug discovery and their associated challenges.

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.