Single-Cell Untargeted Lipidomics Using Liquid Chromatography and Data-Dependent Acquisition after Live Cell Selection.

We report the development and validation of an untargeted single-cell lipidomics method based on microflow chromatography coupled to a data-dependent mass spectrometry method for fragmentation-based identification of lipids. Given the absence of single-cell lipid standards, we show how the methodology should be optimized and validated using a dilute cell extract. The methodology is applied to dilute pancreatic cancer and macrophage cell extracts and standards to demonstrate the sensitivity requirements for confident assignment of lipids and classification of the cell type at the single-cell level. The method is then coupled to a system that can provide automated sampling of live, single cells into capillaries under microscope observation. This workflow retains the spatial information and morphology of cells during sampling and highlights the heterogeneity in lipid profiles observed at the single-cell level. The workflow is applied to show changes in single-cell lipid profiles as a response to oxidative stress, coinciding with expanded lipid droplets. This demonstrates that the workflow is sufficiently sensitive to observing changes in lipid profiles in response to a biological stimulus. Understanding how lipids vary in single cells will inform future research into a multitude of biological processes as lipids play important roles in structural, biophysical, energy storage, and signaling functions.

Nanocapillary sampling coupled to liquid chromatography mass spectrometry delivers single cell drug measurement and lipid fingerprints.

This work describes the development of a new approach to measure drug levels and lipid fingerprints in single living mammalian cells. Nanocapillary sampling is an approach that enables the selection and isolation of single living cells under microscope observation. Here, live single cell nanocapillary sampling is coupled to liquid chromatography for the first time. This allows molecular species to be separated prior to ionisation and improves measurement precision of drug analytes. The efficiency of transferring analytes from the sampling capillary into a vial was optimised in this work. The analysis was carried out using standard flow liquid chromatography coupled to widely available mass spectrometry instrumentation, highlighting opportunities for widespread adoption. The method was applied to 30 living cells, revealing cell-to-cell heterogeneity in the uptake of different drug molecules. Using this system, we detected 14-158 lipid features per single cell, revealing the association between bedaquiline uptake and lipid fingerprints.

One-shot 13 C15 N-metabolic flux analysis for simultaneous quantification of carbon and nitrogen flux.

Metabolic flux is the final output of cellular regulation and has been extensively studied for carbon but much less is known about nitrogen, which is another important building block for living organisms. For the tuberculosis pathogen, this is particularly important in informing the development of effective drugs targeting the pathogen's metabolism. Here we performed 13 C15 N dual isotopic labeling of Mycobacterium bovis BCG steady state cultures, quantified intracellular carbon and nitrogen fluxes and inferred reaction bidirectionalities. This was achieved by model scope extension and refinement, implemented in a multi-atom transition model, within the statistical framework of Bayesian model averaging (BMA). Using BMA-based 13 C15 N-metabolic flux analysis, we jointly resolve carbon and nitrogen fluxes quantitatively. We provide the first nitrogen flux distributions for amino acid and nucleotide biosynthesis in mycobacteria and establish glutamate as the central node for nitrogen metabolism. We improved resolution of the notoriously elusive anaplerotic node in central carbon metabolism and revealed possible operation modes. Our study provides a powerful and statistically rigorous platform to simultaneously infer carbon and nitrogen metabolism in any biological system.

Dual RNA Sequencing of Mycobacterium tuberculosis-Infected Human Splenic Macrophages Reveals a Strain-Dependent Host-Pathogen Response to Infection.

Tuberculosis (TB) is caused by Mycobacterium tuberculosis (Mtb), leading to pulmonary and extrapulmonary TB, whereby Mtb is disseminated to many other organs and tissues. Dissemination occurs early during the disease, and bacteria can be found first in the lymph nodes adjacent to the lungs and then later in the extrapulmonary organs, including the spleen. The early global gene expression response of human tissue macrophages and intracellular clinical isolates of Mtb has been poorly studied. Using dual RNA-seq, we have explored the mRNA profiles of two closely related clinical strains of the Latin American and Mediterranean (LAM) family of Mtb in infected human splenic macrophages (hSMs). This work shows that these pathogens mediate a distinct host response despite their genetic similarity. Using a genome-scale host-pathogen metabolic reconstruction to analyze the data further, we highlight that the infecting Mtb strain also determines the metabolic response of both the host and pathogen. Thus, macrophage ontogeny and the genetic-derived program of Mtb direct the host-pathogen interaction.

Murine Norovirus Infection Results in Anti-inflammatory Response Downstream of Amino Acid Depletion in Macrophages.

Murine norovirus (MNV) infection results in a late translation shutoff that is proposed to contribute to the attenuated and delayed innate immune response observed both in vitro and in vivo. Recently, we further demonstrated the activation of the α subunit of eukaryotic initiation factor 2 (eIF2α) kinase GCN2 during MNV infection, which has been previously linked to immunomodulation and resistance to inflammatory signaling during metabolic stress. While viral infection is usually associated with activation of double-stranded RNA (dsRNA) binding pattern recognition receptor PKR, we hypothesized that the establishment of a metabolic stress in infected cells is a proviral event, exploited by MNV to promote replication through weakening the activation of the innate immune response. In this study, we used multi-omics approaches to characterize cellular responses during MNV replication. We demonstrate the activation of pathways related to the integrated stress response, a known driver of anti-inflammatory phenotypes in macrophages. In particular, MNV infection causes an amino acid imbalance that is associated with GCN2 and ATF2 signaling. Importantly, this reprogramming lacks the features of a typical innate immune response, with the ATF/CHOP target GDF15 contributing to the lack of antiviral responses. We propose that MNV-induced metabolic stress supports the establishment of host tolerance to viral replication and propagation. IMPORTANCE During viral infection, host defenses are typically characterized by the secretion of proinflammatory autocrine and paracrine cytokines, potentiation of the interferon (IFN) response, and induction of the antiviral response via activation of JAK and Stat signaling. To avoid these and propagate, viruses have evolved strategies to evade or counteract host sensing. In this study, we demonstrate that murine norovirus controls the antiviral response by activating a metabolic stress response that activates the amino acid response and impairs inflammatory signaling. This highlights novel tools in the viral countermeasures arsenal and demonstrates the importance of the currently poorly understood metabolic reprogramming occurring during viral infections.

Metabolic fluxes for nutritional flexibility of Mycobacterium tuberculosis.

The co-catabolism of multiple host-derived carbon substrates is required by Mycobacterium tuberculosis (Mtb) to successfully sustain a tuberculosis infection. However, the metabolic plasticity of this pathogen and the complexity of the metabolic networks present a major obstacle in identifying those nodes most amenable to therapeutic interventions. It is therefore critical that we define the metabolic phenotypes of Mtb in different conditions. We applied metabolic flux analysis using stable isotopes and lipid fingerprinting to investigate the metabolic network of Mtb growing slowly in our steady-state chemostat system. We demonstrate that Mtb efficiently co-metabolises either cholesterol or glycerol, in combination with two-carbon generating substrates without any compartmentalisation of metabolism. We discovered that partitioning of flux between the TCA cycle and the glyoxylate shunt combined with a reversible methyl citrate cycle is the critical metabolic nodes which underlie the nutritional flexibility of Mtb. These findings provide novel insights into the metabolic architecture that affords adaptability of bacteria to divergent carbon substrates and expand our fundamental knowledge about the methyl citrate cycle and the glyoxylate shunt.

Structural and Molecular Dynamics of Mycobacterium tuberculosis Malic Enzyme, a Potential Anti-TB Drug Target.

Tuberculosis (TB) is the most lethal bacterial infectious disease worldwide. It is notoriously difficult to treat, requiring a cocktail of antibiotics administered over many months. The dense, waxy outer membrane of the TB-causing agent, Mycobacterium tuberculosis (Mtb), acts as a formidable barrier against uptake of antibiotics. Subsequently, enzymes involved in maintaining the integrity of the Mtb cell wall are promising drug targets. Recently, we demonstrated that Mtb lacking malic enzyme (MEZ) has altered cell wall lipid composition and attenuated uptake by macrophages. These results suggest that MEZ contributes to lipid biosynthesis by providing reductants in the form of NAD(P)H. Here, we present the X-ray crystal structure of MEZ to 3.6 Å. We use biochemical assays to demonstrate MEZ is dimeric in solution and to evaluate the effects of pH and allosteric regulators on its kinetics and thermal stability. To assess the interactions between MEZ and its substrate malate and cofactors, Mn2+ and NAD(P)+, we ran a series of molecular dynamics (MD) simulations. First, the MD analysis corroborates our empirical observations that MEZ is unusually flexible, which persists even with the addition of substrate and cofactors. Second, the MD simulations reveal that dimeric MEZ subunits alternate between open and closed states, and that MEZ can stably bind its NAD(P)+ cofactor in multiple conformations, including an inactive, compact NAD+ form. Together the structure of MEZ and insights from its dynamics can be harnessed to inform the design of MEZ inhibitors that target Mtb and not human malic enzyme homologues.

Bedaquiline reprograms central metabolism to reveal glycolytic vulnerability in Mycobacterium tuberculosis.

The approval of bedaquiline (BDQ) for the treatment of tuberculosis has generated substantial interest in inhibiting energy metabolism as a therapeutic paradigm. However, it is not known precisely how BDQ triggers cell death in Mycobacterium tuberculosis (Mtb). Using 13C isotopomer analysis, we show that BDQ-treated Mtb redirects central carbon metabolism to induce a metabolically vulnerable state susceptible to genetic disruption of glycolysis and gluconeogenesis. Metabolic flux profiles indicate that BDQ-treated Mtb is dependent on glycolysis for ATP production, operates a bifurcated TCA cycle by increasing flux through the glyoxylate shunt, and requires enzymes of the anaplerotic node and methylcitrate cycle. Targeting oxidative phosphorylation (OXPHOS) with BDQ and simultaneously inhibiting substrate level phosphorylation via genetic disruption of glycolysis leads to rapid sterilization. Our findings provide insight into the metabolic mechanism of BDQ-induced cell death and establish a paradigm for the development of combination therapies that target OXPHOS and glycolysis.

A systematic evaluation of Mycobacterium tuberculosis Genome-Scale Metabolic Networks.

Metabolism underpins the pathogenic strategy of the causative agent of TB, Mycobacterium tuberculosis (Mtb), and therefore metabolic pathways have recently re-emerged as attractive drug targets. A powerful approach to study Mtb metabolism as a whole, rather than just individual enzymatic components, is to use a systems biology framework, such as a Genome-Scale Metabolic Network (GSMN) that allows the dynamic interactions of all the components of metabolism to be interrogated together. Several GSMNs networks have been constructed for Mtb and used to study the complex relationship between the Mtb genotype and its phenotype. However, the utility of this approach is hampered by the existence of multiple models, each with varying properties and performances. Here we systematically evaluate eight recently published metabolic models of Mtb-H37Rv to facilitate model choice. The best performing models, sMtb2018 and iEK1011, were refined and improved for use in future studies by the TB research community.

Intracellular Mycobacterium leprae Utilizes Host Glucose as a Carbon Source in Schwann Cells.

New approaches are needed to control leprosy, but understanding of the biology of the causative agent Mycobacterium leprae remains rudimentary, principally because the pathogen cannot be grown in axenic culture. Here, we applied 13C isotopomer analysis to measure carbon metabolism of M. leprae in its primary host cell, the Schwann cell. We compared the results of this analysis with those of a related pathogen, Mycobacterium tuberculosis, growing in its primary host cell, the macrophage. Using 13C isotopomer analysis with glucose as the tracer, we show that whereas M. tuberculosis imports most of its amino acids directly from the host macrophage, M. leprae utilizes host glucose pools as the carbon source to biosynthesize the majority of its amino acids. Our analysis highlights the anaplerotic enzyme phosphoenolpyruvate carboxylase required for this intracellular diet of M. leprae, identifying this enzyme as a potential antileprosy drug target.IMPORTANCE Leprosy remains a major problem in the world today, particularly affecting the poorest and most disadvantaged sections of society in the least developed countries of the world. The long-term aim of research is to develop new treatments and vaccines, and these aims are currently hampered by our inability to grow the pathogen in axenic culture. In this study, we probed the metabolism of M. leprae while it is surviving and replicating inside its primary host cell, the Schwann cell, and compared it to a related pathogen, M. tuberculosis, replicating in macrophages. Our analysis revealed that unlike M. tuberculosis, M. leprae utilized host glucose as a carbon source and that it biosynthesized its own amino acids, rather than importing them from its host cell. We demonstrated that the enzyme phosphoenolpyruvate carboxylase plays a crucial role in glucose catabolism in M. leprae Our findings provide the first metabolic signature of M. leprae in the host Schwann cell and identify novel avenues for the development of antileprosy drugs.

Intracellular Mycobacterium tuberculosis Exploits Multiple Host Nitrogen Sources during Growth in Human Macrophages.

Nitrogen metabolism of Mycobacterium tuberculosis (Mtb) is crucial for the survival of this important pathogen in its primary human host cell, the macrophage, but little is known about the source(s) and their assimilation within this intracellular niche. Here, we have developed 15N-flux spectral ratio analysis (15N-FSRA) to explore Mtb's nitrogen metabolism; we demonstrate that intracellular Mtb has access to multiple amino acids in the macrophage, including glutamate, glutamine, aspartate, alanine, glycine, and valine; and we identify glutamine as the predominant nitrogen donor. Each nitrogen source is uniquely assimilated into specific amino acid pools, indicating compartmentalized metabolism during intracellular growth. We have discovered that serine is not available to intracellular Mtb, and we show that a serine auxotroph is attenuated in macrophages. This work provides a systems-based tool for exploring the nitrogen metabolism of intracellular pathogens and highlights the enzyme phosphoserine transaminase as an attractive target for the development of novel anti-tuberculosis therapies.

The anaplerotic node is essential for the intracellular survival of Mycobacterium tuberculosis.

Enzymes at the phosphoenolpyruvate (PEP)-pyruvate-oxaloacetate or anaplerotic (ANA) node control the metabolic flux to glycolysis, gluconeogenesis, and anaplerosis. Here we used genetic, biochemical, and 13C isotopomer analysis to characterize the role of the enzymes at the ANA node in intracellular survival of the world's most successful bacterial pathogen, Mycobacterium tuberculosis (Mtb). We show that each of the four ANA enzymes, pyruvate carboxylase (PCA), PEP carboxykinase (PCK), malic enzyme (MEZ), and pyruvate phosphate dikinase (PPDK), performs a unique and essential metabolic function during the intracellular survival of Mtb. We show that in addition to PCK, intracellular Mtb requires PPDK as an alternative gateway into gluconeogenesis. Propionate and cholesterol detoxification was also identified as an essential function of PPDK revealing an unexpected role for the ANA node in the metabolism of these physiologically important intracellular substrates and highlighting this enzyme as a tuberculosis (TB)-specific drug target. We show that anaplerotic fixation of CO2 through the ANA node is essential for intracellular survival of Mtb and that Mtb possesses three enzymes (PCA, PCK, and MEZ) capable of fulfilling this function. In addition to providing a back-up role in anaplerosis we show that MEZ also has a role in lipid biosynthesis. MEZ knockout strains have an altered cell wall and were deficient in the initial entry into macrophages. This work reveals that the ANA node is a focal point for controlling the intracellular replication of Mtb, which goes beyond canonical gluconeogenesis and represents a promising target for designing novel anti-TB drugs.

Systems-based approaches to probing metabolic variation within the Mycobacterium tuberculosis complex.

The Mycobacterium tuberculosis complex includes bovine and human strains of the tuberculosis bacillus, including Mycobacterium tuberculosis, Mycobacterium bovis and the Mycobacterium bovis BCG vaccine strain. M. bovis has evolved from a M. tuberculosis-like ancestor and is the ancestor of the BCG vaccine. The pathogens demonstrate distinct differences in virulence, host range and metabolism, but the role of metabolic differences in pathogenicity is poorly understood. Systems biology approaches have been used to investigate the metabolism of M. tuberculosis, but not to probe differences between tuberculosis strains. In this study genome scale metabolic networks of M. bovis and M. bovis BCG were constructed and interrogated, along with a M. tuberculosis network, to predict substrate utilisation, gene essentiality and growth rates. The models correctly predicted 87-88% of high-throughput phenotype data, 75-76% of gene essentiality data and in silico-predicted growth rates matched measured rates. However, analysis of the metabolic networks identified discrepancies between in silico predictions and in vitro data, highlighting areas of incomplete metabolic knowledge. Additional experimental studies carried out to probe these inconsistencies revealed novel insights into the metabolism of these strains. For instance, that the reduction in metabolic capability observed in bovine tuberculosis strains, as compared to M. tuberculosis, is not reflected by current genetic or enzymatic knowledge. Hence, the in silico networks not only successfully simulate many aspects of the growth and physiology of these mycobacteria, but also provide an invaluable tool for future metabolic studies.

13C-flux spectral analysis of host-pathogen metabolism reveals a mixed diet for intracellular Mycobacterium tuberculosis.

Whereas intracellular carbon metabolism has emerged as an attractive drug target, the carbon sources of intracellularly replicating pathogens, such as the tuberculosis bacillus Mycobacterium tuberculosis, which causes long-term infections in one-third of the world's population, remain mostly unknown. We used a systems-based approach--(13)C-flux spectral analysis (FSA) complemented with manual analysis-to measure the metabolic interaction between M. tuberculosis and its macrophage host cell. (13)C-FSA analysis of experimental data showed that M. tuberculosis obtains a mixture of amino acids, C1 and C2 substrates from its host cell. We experimentally confirmed that the C1 substrate was derived from CO2. (13)C labeling experiments performed on a phosphoenolpyruvate carboxykinase mutant revealed that intracellular M. tuberculosis has access to glycolytic C3 substrates. These findings provide constraints for developing novel chemotherapeutics.

¹³C metabolic flux analysis identifies an unusual route for pyruvate dissimilation in mycobacteria which requires isocitrate lyase and carbon dioxide fixation.

Mycobacterium tuberculosis requires the enzyme isocitrate lyase (ICL) for growth and virulence in vivo. The demonstration that M. tuberculosis also requires ICL for survival during nutrient starvation and has a role during steady state growth in a glycerol limited chemostat indicates a function for this enzyme which extends beyond fat metabolism. As isocitrate lyase is a potential drug target elucidating the role of this enzyme is of importance; however, the role of isocitrate lyase has never been investigated at the level of in vivo fluxes. Here we show that deletion of one of the two icl genes impairs the replication of Mycobacterium bovis BCG at slow growth rate in a carbon limited chemostat. In order to further understand the role of isocitrate lyase in the central metabolism of mycobacteria the effect of growth rate on the in vivo fluxes was studied for the first time using ¹³C-metabolic flux analysis (MFA). Tracer experiments were performed with steady state chemostat cultures of BCG or M. tuberculosis supplied with ¹³C labeled glycerol or sodium bicarbonate. Through measurements of the ¹³C isotopomer labeling patterns in protein-derived amino acids and enzymatic activity assays we have identified the activity of a novel pathway for pyruvate dissimilation. We named this the GAS pathway because it utilizes the Glyoxylate shunt and Anapleurotic reactions for oxidation of pyruvate, and Succinyl CoA synthetase for the generation of succinyl CoA combined with a very low flux through the succinate--oxaloacetate segment of the tricarboxylic acid cycle. We confirm that M. tuberculosis can fix carbon from CO2 into biomass. As the human host is abundant in CO2 this finding requires further investigation in vivo as CO2 fixation may provide a point of vulnerability that could be targeted with novel drugs. This study also provides a platform for further studies into the metabolism of M. tuberculosis using ¹³C-MFA.

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

A general paucity of knowledge about the metabolic state of Mycobacterium tuberculosis within the host environment is a major factor impeding development of novel drugs against tuberculosis. Current experimental methods do not allow direct determination of the global metabolic state of a bacterial pathogen in vivo, but the transcriptional activity of all encoded genes has been investigated in numerous microarray studies. We describe a novel algorithm, Differential Producibility Analysis (DPA) that uses a metabolic network to extract metabolic signals from transcriptome data. The method utilizes Flux Balance Analysis (FBA) to identify the set of genes that affect the ability to produce each metabolite in the network. Subsequently, Rank Product Analysis is used to identify those metabolites predicted to be most affected by a transcriptional signal. We first apply DPA to investigate the metabolic response of E. coli to both anaerobic growth and inactivation of the FNR global regulator. DPA successfully extracts metabolic signals that correspond to experimental data and provides novel metabolic insights. We next apply DPA to investigate the metabolic response of M. tuberculosis to the macrophage environment, human sputum and a range of in vitro environmental perturbations. The analysis revealed a previously unrecognized feature of the response of M. tuberculosis to the macrophage environment: a down-regulation of genes influencing metabolites in central metabolism and concomitant up-regulation of genes that influence synthesis of cell wall components and virulence factors. DPA suggests that a significant feature of the response of the tubercle bacillus to the intracellular environment is a channeling of resources towards remodeling of its cell envelope, possibly in preparation for attack by host defenses. DPA may be used to unravel the mechanisms of virulence and persistence of M. tuberculosis and other pathogens and may have general application for extracting metabolic signals from other "-omics" data.

System-level strategies for studying the metabolism of Mycobacterium tuberculosis.

Despite decades of research many aspects of the biology of Mycobacterium tuberculosis remain unclear and this is reflected in the antiquated tools available to treat and prevent tuberculosis and consequently this disease remains a serious public health problem. Important discoveries linking M. tuberculosis's metabolism and pathogenesis have renewed interest in this area of research. Previous experimental studies were limited to the analysis of individual genes or enzymes whereas recent advances in computational systems biology and high throughput experimental technologies now allow metabolism to be studied on a genome scale. Here we discuss the progress being made in applying system level approaches to studying the metabolism of this important pathogen. The information from these studies will fundamentally change our approach to tuberculosis research and lead to new targets for therapeutic drugs and vaccines.

Systems biology of the metabolism of Mycobacterium tuberculosis.

Despite decades of research, many aspects of the biology of Mycobacterium tuberculosis remain unclear, and this is reflected in the antiquated tools available to treat and prevent tuberculosis and consequently this disease remains a serious public health problem. Important discoveries linking the metabolism of M. tuberculosis and pathogenesis has renewed interest in this area of research. Previous experimental studies were limited to the analysis of individual genes or enzymes, whereas recent advances in computational systems biology and high-throughput experimental technologies now allows metabolism to be studied on a genome scale. In the present article, we discuss the progress being made in applying system-level approaches to study the metabolism of this important pathogen.

The genetic requirements for fast and slow growth in mycobacteria.

Mycobacterium tuberculosis infects a third of the world's population. Primary tuberculosis involving active fast bacterial replication is often followed by asymptomatic latent tuberculosis, which is characterised by slow or non-replicating bacteria. Reactivation of the latent infection involving a switch back to active bacterial replication can lead to post-primary transmissible tuberculosis. Mycobacterial mechanisms involved in slow growth or switching growth rate provide rational targets for the development of new drugs against persistent mycobacterial infection. Using chemostat culture to control growth rate, we screened a transposon mutant library by Transposon site hybridization (TraSH) selection to define the genetic requirements for slow and fast growth of Mycobacterium bovis (BCG) and for the requirements of switching growth rate. We identified 84 genes that are exclusively required for slow growth (69 hours doubling time) and 256 genes required for switching from slow to fast growth. To validate these findings we performed experiments using individual M. tuberculosis and M. bovis BCG knock out mutants. We have demonstrated that growth rate control is a carefully orchestrated process which requires a distinct set of genes encoding several virulence determinants, gene regulators, and metabolic enzymes. The mce1 locus appears to be a component of the switch to slow growth rate, which is consistent with the proposed role in virulence of M. tuberculosis. These results suggest novel perspectives for unravelling the mechanisms involved in the switch between acute and persistent TB infections and provide a means to study aspects of this important phenomenon in vitro.