Sulfite reduction in mycobacteria

Divergent downstream biosynthetic pathways are supported by <sc>L</sc>-cysteine synthases of Mycobacterium tuberculosis

Mycobacterium tuberculosis's (Mtb) autarkic lifestyle within the host involves rewiring its transcriptional networks to combat host-induced stresses. With the help of RNA sequencing performed under various stress conditions, we identified that genes belonging to Mtb sulfur metabolism pathways are significantly upregulated during oxidative stress. Using an integrated approach of microbial genetics, transcriptomics, metabolomics, animal experiments, chemical inhibition, and rescue studies, we investigated the biological role of non-canonical L-cysteine synthases, CysM and CysK2. While transcriptome signatures of RvΔcysM and RvΔcysK2 appear similar under regular growth conditions, we observed unique transcriptional signatures when subjected to oxidative stress. We followed pool size and labelling (34S) of key downstream metabolites, viz. mycothiol and ergothioneine, to monitor L-cysteine biosynthesis and utilization. This revealed the significant role of distinct L-cysteine biosynthetic routes on redox stress and homeostasis. CysM and CysK2 independently facilitate Mtb survival by alleviating host-induced redox stress, suggesting they are not fully redundant during infection. With the help of genetic mutants and chemical inhibitors, we show that CysM and CysK2 serve as unique, attractive targets for adjunct therapy to combat mycobacterial infection.

Dimerization of assimilatory NADPH-dependent sulfite reductase reveals elements for diflavin reductase binding at a minimal interface

Escherichia coli NADPH-dependent assimilatory sulfite reductase is responsible for fixing sulfur for incorporation into sulfur-containing biomolecules. The oxidoreductase is composed of two subunits, an NADPH, FMN, and FAD-binding diflavin reductase and an iron siroheme and Fe4S4-containing oxidase. How they interact has been an unknown for over 50 years because the complex is highly flexible, thus has been intransigent for traditional X-ray or cryo-EM structural analysis. Using a combination of the chameleon plunging system with a fluorinated lipid we overcame the challenge of preserving the minimal dimer between the subunits for high-resolution cryo-EM analysis. Here, we report the first structure of the complex between the reductase and oxidase, revealing how they interact in a minimal interface. Further, we determined the structural elements that discriminate between the pairing of a siroheme-containing oxidase with a diflavin reductase or a ferredoxin partner to channel the six electrons that reduce sulfite to sulfide.

Genome Sequence of Eubacterium limosum B2 and Evolution for Growth on a Mineral Medium with Methanol and CO2 as Sole Carbon Sources

Eubacterium limosum is an acetogen that can produce butyrate along with acetate as the main fermentation end-product from methanol, a promising C1 feedstock. Although physiological characterization of E. limosum B2 during methylotrophy was previously performed, the strain was cultured in a semi-defined medium, limiting the scope for further metabolic insights. Here, we sequenced the complete genome of the native strain and performed adaptive laboratory evolution to sustain growth on methanol mineral medium. The evolved population significantly improved its maximal growth rate by 3.45-fold. Furthermore, three clones from the evolved population were isolated on methanol mineral medium without cysteine by the addition of sodium thiosulfate. To identify mutations related to growth improvement, the whole genomes of wild-type E. limosum B2, the 10th, 25th, 50th, and 75th generations, and the three clones were sequenced. We explored the total proteomes of the native and the best evolved clone (n°2) and noticed significant differences in proteins involved in gluconeogenesis, anaplerotic reactions, and sulphate metabolism. Furthermore, a homologous recombination was found in subunit S of the type I restriction-modification system between both strains, changing the structure of the subunit, its sequence recognition and the methylome of the evolved clone. Taken together, the genomic, proteomic and methylomic data suggest a possible epigenetic mechanism of metabolic regulation.

Biodesulfurization Induces Reprogramming of Sulfur Metabolism in Rhodococcus qingshengii IGTS8: Proteomics and Untargeted Metabolomics

Sulfur metabolism in fuel-biodesulfurizing bacteria and the underlying physiological adaptations are not understood, which has impeded the development of a commercially viable bioprocess for fuel desulfurization. To fill these knowledge gaps, we performed comparative proteomics and untargeted metabolomics in cultures of the biodesulfurization reference strain Rhodococcus qingshengii IGTS8 grown on either inorganic sulfate or the diesel-borne organosulfur compound dibenzothiophene as a sole sulfur source. Dibenzothiophene significantly altered the biosynthesis of many sulfur metabolism proteins and metabolites in a growth phase-dependent manner, which enabled us to reconstruct the first experimental model for sulfur metabolism in a fuel-biodesulfurizing bacterium. All key pathways related to assimilatory sulfur metabolism were represented in the sulfur proteome, including uptake of the sulfur sources, sulfur acquisition, and assimilatory sulfate reduction, in addition to biosynthesis of key sulfur-containing metabolites such as S-adenosylmethionine, coenzyme A, biotin, thiamin, molybdenum cofactor, mycothiol, and ergothioneine (low-molecular weight thiols). Fifty-two proteins exhibited significantly different abundance during at least one growth phase. Sixteen proteins were uniquely detected and 47 proteins were significantly more abundant in the dibenzothiophene culture during at least one growth phase. The sulfate-free dibenzothiophene-containing culture reacted to sulfate starvation by restricting sulfur assimilation, enforcing sulfur-sparing, and maintaining redox homeostasis. Biodesulfurization triggered alternative pathways for sulfur assimilation different from those operating in the inorganic sulfate culture. Sulfur metabolism reprogramming and metabolic switches in the dibenzothiophene culture were manifested in limiting sulfite reduction and biosynthesis of cysteine, while boosting the production of methionine via the cobalamin-independent pathway, as well as the biosynthesis of the redox buffers mycothiol and ergothioneine. The omics data underscore the key role of sulfur metabolism in shaping the biodesulfurization phenotype and highlight potential targets for improving the biodesulfurization catalytic activity via metabolic engineering.

IMPORTANCE For many decades, research on biodesulfurization of fossil fuels was conducted amid a large gap in knowledge of sulfur metabolism and its regulation in fuel-biodesulfurizing bacteria, which has impeded the development of a commercially viable bioprocess. In addition, lack of understanding of biodesulfurization-associated metabolic and physiological adaptations prohibited the development of efficient biodesulfurizers. Our integrated omics-based findings reveal the assimilatory sulfur metabolism in the biodesulfurization reference strain Rhodococcus qingshengii IGTS8 and show how sulfur metabolism and oxidative stress response were remodeled and orchestrated to shape the biodesulfurization phenotype. Our findings not only explain the frequently encountered low catalytic activity of native fuel-biodesulfurizing bacteria but also uncover unprecedented potential targets in sulfur metabolism that could be exploited via metabolic engineering to boost the biodesulfurization catalytic activity, a prerequisite for commercial application.

Mycobacterium tuberculosis H2S Functions as a Sink to Modulate Central Metabolism, Bioenergetics, and Drug Susceptibility

H2S is a potent gasotransmitter in eukaryotes and bacteria. Host-derived H2S has been shown to profoundly alter M. tuberculosis (Mtb) energy metabolism and growth. However, compelling evidence for endogenous production of H2S and its role in Mtb physiology is lacking. We show that multidrug-resistant and drug-susceptible clinical Mtb strains produce H2S, whereas H2S production in non-pathogenic M. smegmatis is barely detectable. We identified Rv3684 (Cds1) as an H2S-producing enzyme in Mtb and show that cds1 disruption reduces, but does not eliminate, H2S production, suggesting the involvement of multiple genes in H2S production. We identified endogenous H2S to be an effector molecule that maintains bioenergetic homeostasis by stimulating respiration primarily via cytochrome bd. Importantly, H2S plays a key role in central metabolism by modulating the balance between oxidative phosphorylation and glycolysis, and it functions as a sink to recycle sulfur atoms back to cysteine to maintain sulfur homeostasis. Lastly, Mtb-generated H2S regulates redox homeostasis and susceptibility to anti-TB drugs clofazimine and rifampicin. These findings reveal previously unknown facets of Mtb physiology and have implications for routine laboratory culturing, understanding drug susceptibility, and improved diagnostics.

C-terminal Redox Domain of Arabidopsis APR1 is a Non-Canonical Thioredoxin Domain with Glutaredoxin Function

Sulfur is an essential nutrient that can be converted into utilizable metabolic forms to produce sulfur-containing metabolites in plant. Adenosine 5'-phosphosulfate (APS) reductase (APR) plays a vital role in catalyzing the reduction of activated sulfate to sulfite, which requires glutathione. Previous studies have shown that the C-terminal domain of APR acts as a glutathione-dependent reductase. The crystal structure of the C-terminal redox domain of Arabidopsis APR1 (AtAPR1) shows a conserved α/β thioredoxin fold, but not a glutaredoxin fold. Further biochemical studies of the redox domain from AtAPR1 provided evidence to support the structural observation. Collectively, our results provide structural and biochemical information to explain how the thioredoxin fold exerts the glutaredoxin function in APR.

Transcriptional Approach for Decoding the Mechanism of rpoC Compensatory Mutations for the Fitness Cost in Rifampicin-Resistant Mycobacterium tuberculosis

Multidrug-resistant tuberculosis (TB), defined as TB resistant to the two first-line drugs, isoniazid and rifampin, is a serious challenge to global TB eradication efforts. Although mutations in rpoA or rpoC have been proposed to compensate for this fitness cost due to rpoB mutation in rifampicin-resistant Mycobacterium tuberculosis mutants, whether the compensatory effect exists and the underlying mechanisms of compensation remain unclear. Here, we used RNA sequencing to investigate the global transcriptional profiles of 6 rifampin-resistant clinical isolates with either single mutation in rpoB or dual mutations in rpoB/rpoC, as well as 3 rifampin-susceptible clinical isolates, trying to prove the potential compensatory effect of rpoC by transcriptomic alteration. In rifampin-free conditions, rpoC mutation was associated with M. tuberculosis upregulation of ribosomal protein-coding genes, dysregulation of growth-related essential genes and balancing the expression of arginine and glutamate synthesis-associated genes. Upon rifampin exposure of M. tuberculosis isolates, rpoC mutations were associated with the upregulation of the oxidative phosphorylation machinery, which was inhibited in the rpoB single mutants, as well as stabilization of the expression of rifampin-regulated essential genes and balancing the expression of genes involved in metabolism of sulfur-containing amino acids. Taken together, our data suggest that rpoC mutation may compensate for the fitness defect of rifampicin-resistant M. tuberculosis by altering gene expression in response to rifampin exposure.

Biosynthesis of Orthogonal Molecules Using Ferredoxin and Ferredoxin-NADP+ Reductase Systems Enables Genetically Encoded PhyB Optogenetics

Transplanting metabolic reactions from one species into another has many uses as a research tool with applications ranging from optogenetics to crop production. Ferredoxin (Fd), the enzyme that most often supplies electrons to these reactions, is often overlooked when transplanting enzymes from one species to another because most cells already contain endogenous Fd. However, we have shown that the production of chromophores used in Phytochrome B (PhyB) optogenetics is greatly enhanced in mammalian cells by expressing bacterial and plant Fds with ferredoxin-NADP+ reductases (FNR). We delineated the rate limiting factors and found that the main metabolic precursor, heme, was not the primary limiting factor for producing either the cyanobacterial or plant chromophores, phycocyanobilin or phytochromobilin, respectively. In fact, Fd is limiting, followed by Fd+FNR and finally heme. Using these findings, we optimized the PCB production system and combined it with a tissue penetrating red/far-red sensing PhyB optogenetic gene switch in animal cells. We further characterized this system in several mammalian cell lines using red and far-red light. Importantly, we found that the light-switchable gene system remains active for several hours upon illumination, even with a short light pulse, and requires very small amounts of light for maximal activation. Boosting chromophore production by matching metabolic pathways with specific ferredoxin systems will enable the unparalleled use of the many PhyB optogenetic tools and has broader implications for optimizing synthetic metabolic pathways.

The transcriptome of Mycobacterium tuberculosis in a lipid-rich dormancy model through RNAseq analysis

Tuberculosis (TB) is currently the number one killer among infectious diseases worldwide. Lipids are abundant molecules during the infectious cycle of Mycobacterium tuberculosis (Mtb) and studies better mimicking its actual metabolic state during pathogenesis are needed. Though most studies have focused on the mycobacterial lipid metabolism under standard culture conditions, little is known about the transcriptome of Mtb in a lipid environment. Here we determined the transcriptome of Mtb H37Rv in a lipid-rich environment (cholesterol and fatty acid) under aerobic and hypoxic conditions, using RNAseq. Lipids significantly induced the expression of 368 genes. A main core lipid response was observed involving efflux systems, iron caption and sulfur reduction. In co-expression with ncRNAs and other genes discussed below, may act coordinately to prepare the machinery conferring drug tolerance and increasing a persistent population. Our findings could be useful to tag relevant pathways for the development of new drugs, vaccines and new strategies to control TB.

Energetics of Respiration and Oxidative Phosphorylation in Mycobacteria

Mycobacteria inhabit a wide range of intracellular and extracellular environments. Many of these environments are highly dynamic and therefore mycobacteria are faced with the constant challenge of redirecting their metabolic activity to be commensurate with either replicative growth or a non-replicative quiescence. A fundamental feature in this adaptation is the ability of mycobacteria to respire, regenerate reducing equivalents and generate ATP via oxidative phosphorylation. Mycobacteria harbor multiple primary dehydrogenases to fuel the electron transport chain and two terminal respiratory oxidases, an aa3 -type cytochrome c oxidase and cytochrome bd-type menaquinol oxidase, are present for dioxygen reduction coupled to the generation of a protonmotive force. Hypoxia leads to the downregulation of key respiratory complexes, but the molecular mechanisms regulating this expression are unknown. Despite being obligate aerobes, mycobacteria have the ability to metabolize in the absence of oxygen and a number of reductases are present to facilitate the turnover of reducing equivalents under these conditions (e.g. nitrate reductase, succinate dehydrogenase/fumarate reductase). Hydrogenases and ferredoxins are also present in the genomes of mycobacteria suggesting the ability of these bacteria to adapt to an anaerobic-type of metabolism in the absence of oxygen. ATP synthesis by the membrane-bound F1FO-ATP synthase is essential for growing and non-growing mycobacteria and the enzyme is able to function over a wide range of protonmotive force values (aerobic to hypoxic). The discovery of lead compounds that target respiration and oxidative phosphorylation in Mycobacterium tuberculosis highlights the importance of this area for the generation of new front line drugs to combat tuberculosis.

Pyridoxal-phosphate dependent mycobacterial cysteine synthases: Structure, mechanism and potential as drug targets

The alarming increase of drug resistance in Mycobacterium tuberculosis strains poses a severe threat to human health. Chemotherapy is particularly challenging because M. tuberculosis can persist in the lungs of infected individuals; estimates of the WHO indicate that about 1/3 of the world population is infected with latent tuberculosis providing a large reservoir for relapse and subsequent spread of the disease. Persistent M. tuberculosis shows considerable tolerance towards conventional antibiotics making treatment particularly difficult. In this phase the bacilli are exposed to oxygen and nitrogen radicals generated as part of the host response and redox-defense mechanisms are thus vital for the survival of the pathogen. Sulfur metabolism and de novo cysteine biosynthesis have been shown to be important for the redox homeostasis in persistent M. tuberculosis and these pathways could provide promising targets for novel antibiotics for the treatment of the latent form of the disease. Recent research has provided evidence for three de novo metabolic routes of cysteine biosynthesis in M. tuberculosis, each with a specific PLP dependent cysteine synthase with distinct substrate specificities. In this review we summarize our present understanding of these pathways, with a focus on the advances on functional and mechanistic characterization of mycobacterial PLP dependent cysteine synthases, their role in the various pathways to cysteine, and first attempts to develop specific inhibitors of mycobacterial cysteine biosynthesis. This article is part of a Special Issue entitled: Cofactor-dependent proteins: evolution, chemical diversity and bio-applications.

Energetics of pathogenic bacteria and opportunities for drug development

The emergence and spread of drug-resistant pathogens and our inability to develop new antimicrobials to overcome resistance has inspired scientists to consider new targets for drug development. Cellular bioenergetics is an area showing promise for the development of new antimicrobials, particularly in the discovery of new anti-tuberculosis drugs where several new compounds have entered clinical trials. In this review, we have examined the bioenergetics of various bacterial pathogens, highlighting the versatility of electron donor and acceptor utilisation and the modularity of electron transport chain components in bacteria. In addition to re-examining classical concepts, we explore new literature that reveals the intricacies of pathogen energetics, for example, how Salmonella enterica and Campylobacter jejuni exploit host and microbiota to derive powerful electron donors and sinks; the strategies Mycobacterium tuberculosis and Pseudomonas aeruginosa use to persist in lung tissues; and the importance of sodium energetics and electron bifurcation in the chemiosmotic anaerobe Fusobacterium nucleatum. A combination of physiological, biochemical, and pharmacological data suggests that, in addition to the clinically-approved target F1Fo-ATP synthase, NADH dehydrogenase type II, succinate dehydrogenase, hydrogenase, cytochrome bd oxidase, and menaquinone biosynthesis pathways are particularly promising next-generation drug targets. The realisation of cellular energetics as a rich target space for the development of new antimicrobials will be dependent upon gaining increased understanding of the energetic processes utilised by pathogens in host environments and the ability to design bacterial-specific inhibitors of these processes.

Crystallization of the C-terminal redox domain of the sulfur-assimilatory enzyme APR1 from Arabidopsis thaliana

Plant-type APS reductase (APR), which catalyzes the reduction of activated sulfate to sulfite in plants, consists of a reductase domain and a C-terminal redox domain showing sequence homology to thioredoxin but possessing the activity of glutaredoxin. In order to understand the structural and biochemical properties of the redox domain of plant-type APS reductase, the C-terminal domain of APR1 (APR1C) from Arabidopsis thaliana was crystallized using the sitting-drop vapour-diffusion method. X-ray diffraction data were collected to a resolution of 2.70 Å on the SPXF beamline BL13B1 at the NSRRC, Taiwan. The crystals belonged to space group P43212 or P41212, with unit-cell parameters a = b = 58.2, c = 86.7 Å. With one molecule per asymmetric unit, the crystal volume per unit protein weight (VM) is 2.64 Å(3) Da(-1), which corresponds to a solvent content of approximately 53.49%. Further structure-based functional studies of APR1C would extend knowledge of the molecular mechanism and regulation of APR.

An obligately aerobic soil bacterium activates fermentative hydrogen production to survive reductive stress during hypoxia

Oxygen availability is a major factor and evolutionary force determining the metabolic strategy of bacteria colonizing an environmental niche. In the soil, conditions can switch rapidly between oxia and anoxia, forcing soil bacteria to remodel their energy metabolism accordingly. Mycobacterium is a dominant genus in the soil, and all its species are obligate aerobes. Here we show that an obligate aerobe, the soil actinomycete Mycobacterium smegmatis, adopts an anaerobe-type strategy by activating fermentative hydrogen production to adapt to hypoxia. This process is controlled by the two-component system DosR-DosS/DosT, an oxygen and redox sensor that is well conserved in mycobacteria. We show that DosR tightly regulates the two [NiFe]-hydrogenases: Hyd3 (MSMEG_3931-3928) and Hyd2 (MSMEG_2719-2718). Using genetic manipulation and high-sensitivity GC, we demonstrate that Hyd3 facilitates the evolution of H2 when oxygen is depleted. Combined activity of Hyd2 and Hyd3 was necessary to maintain an optimal NAD(+)/NADH ratio and enhanced adaptation to and survival of hypoxia. We demonstrate that fermentatively-produced hydrogen can be recycled when fumarate or oxygen become available, suggesting Mycobacterium smegmatis can switch between fermentation, anaerobic respiration, and aerobic respiration. Hydrogen metabolism enables this obligate aerobe to rapidly meet its energetic needs when switching between microoxic and anoxic conditions and provides a competitive advantage in low oxygen environments.

Vitamin B(12) metabolism in Mycobacterium tuberculosis

Mycobacterium tuberculosis is included among a select group of bacteria possessing the capacity for de novo biosynthesis of vitamin B12, the largest and most complex natural organometallic cofactor. The bacillus is also able to scavenge B12 and related corrinoids utilizing an ATP-binding cassette-type protein that is distinct from the only known bacterial B12-specific transporter, BtuFCD. Consistent with the inferred requirement for vitamin B12 for metabolic function, the M. tuberculosis genome encodes two B12 riboswitches and three B12-dependent enzymes. Two of these enzymes have been shown to operate in methionine biosynthesis (MetH) and propionate utilization (MutAB), while the function of the putative nrdZ-encoded ribonucleotide reductase remains unknown. Taken together, these observations suggest that M. tuberculosis has the capacity to regulate core metabolic functions according to B12 availability - whether acquired via endogenous synthesis or through uptake from the host environment - and, therefore, imply that there is a role for vitamin B12 in pathogenesis, which remains poorly understood.

Perturbation of cytochrome c maturation reveals adaptability of the respiratory chain in Mycobacterium tuberculosis

Mycobacterium tuberculosis depends on aerobic respiration for growth and utilizes an aa3-type cytochrome c oxidase for terminal electron transfer. Cytochrome c maturation in bacteria requires covalent attachment of heme to apocytochrome c, which occurs outside the cytoplasmic membrane. We demonstrate that in M. tuberculosis the thioredoxin-like protein Rv3673c, which we named CcsX, is required for heme insertion in cytochrome c. Inactivation of CcsX resulted in loss of c-type heme absorbance, impaired growth and virulence of M. tuberculosis, and induced cytochrome bd oxidase. This suggests that the bioenergetically less efficient bd oxidase can compensate for deficient cytochrome c oxidase activity, highlighting the flexibility of the M. tuberculosis respiratory chain. A spontaneous mutation in the active site of vitamin K epoxide reductase (VKOR) suppressed phenotypes of the CcsX mutant and abrogated the activity of the disulfide bond-dependent alkaline phosphatase, which shows that VKOR is the major disulfide bond catalyzing protein in the periplasm of M. tuberculosis.

IMPORTANCE

Mycobacterium tuberculosis requires oxygen for growth; however, the biogenesis of respiratory chain components in mycobacteria has not been explored. Here, we identified a periplasmic thioredoxin, CcsX, necessary for heme insertion into cytochrome c. We investigated the consequences of disrupting cytochrome c maturation (CCM) for growth and survival of M. tuberculosis in vitro and for its pathogenesis. Appearance of a second-site suppressor mutation in the periplasmic disulfide bond catalyzing protein VKOR indicates the strong selective pressure for a functional cytochrome c oxidase. The observation that M. tuberculosis is able to partially compensate for defective CCM by upregulation of the cytochrome bd oxidase exposes a functional role of this alternative terminal oxidase under normal aerobic conditions and during pathogenesis. This suggests that targeting both oxidases simultaneously might be required to effectively disrupt respiration in M. tuberculosis.

The regulation of sulfur metabolism in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) has evolved into a highly successful human pathogen. It deftly subverts the bactericidal mechanisms of alveolar macrophages, ultimately inducing granuloma formation and establishing long-term residence in the host. These hallmarks of Mtb infection are facilitated by the metabolic adaptation of the pathogen to its surrounding environment and the biosynthesis of molecules that mediate its interactions with host immune cells. The sulfate assimilation pathway of Mtb produces a number of sulfur-containing metabolites with important contributions to pathogenesis and survival. This pathway is regulated by diverse environmental cues and regulatory proteins that mediate sulfur transactions in the cell. Here, we discuss the transcriptional and biochemical mechanisms of sulfur metabolism regulation in Mtb and potential small molecule regulators of the sulfate assimilation pathway that are collectively poised to aid this intracellular pathogen in its expert manipulation of the host. From this global analysis, we have identified a subset of sulfur-metabolizing enzymes that are sensitive to multiple regulatory cues and may be strong candidates for therapeutic intervention.

aprABC: a Mycobacterium tuberculosis complex-specific locus that modulates pH-driven adaptation to the macrophage phagosome

Following phagocytosis by macrophages, Mycobacterium tuberculosis (Mtb) senses the intracellular environment and remodels its gene expression for growth in the phagosome. We have identified an acid and phagosome regulated (aprABC) locus that is unique to the Mtb complex and whose gene expression is induced during growth in acidic environments in vitro and in macrophages. Using the aprA promoter, we generated a strain that exhibits high levels of inducible fluorescence in response to growth in acidic medium in vitro and in macrophages. aprABC expression is dependent on the two-component regulator phoPR, linking phoPR signalling to pH sensing. Deletion of the aprABC locus causes defects in gene expression that impact aggregation, intracellular growth, and the relative levels of storage and cell wall lipids. We propose a model where phoPR senses the acidic pH of the phagosome and induces aprABC expression to fine-tune processes unique for intracellular adaptation of Mtb complex bacteria.

Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions

A human alveolar macrophage genome-scale metabolic reconstruction was reconstructed from tailoring a global human metabolic network, Recon 1, by using computational algorithms and manual curation.

A genome-scale host–pathogen network of the human alveolar macrophage and Mycobacterium tuberculosis is presented. This involved integrating two genome-scale network reconstructions.

The reaction activity and gene essentiality predictions of the host–pathogen model represent a more accurate depiction of infection.

Integration of high-throughput data into a host-pathogen model followed by systems analysis was performed in order to elucidate major metabolic differences under different types of M. tuberculosis infection.

Mycobacterium tuberculosis (M. tb) is an insidious and highly persistent pathogen that affects one-third of the world's population (WHO, 2009). Metabolism is foundational to M. tb's infection ability and the ensuing host–pathogen interactions. In addition, M. tb has a heterogeneous clinical presentation and can infect virtually every tissue. Depending on the location of the infection, different metabolic pathways are active and inactive in both the host and pathogen cells. In this study, we sought to model the host–pathogen interactions of the human alveolar macrophage and M. tb as well as detail the metabolic differences in specific infection types using genome-scale metabolic reconstructions (Figure 4A).

Genome-scale metabolic reconstructions are knowledge bases of all known metabolic reactions of a given organism. Reconstructions have been shown to elucidate the mechanistic genotype-to-phenotype relationship through the integration of high-throughput and physiological data (Oberhardt et al, 2009). Genome-scale reconstructions are converted into mathematical models under the constraints-based reconstruction and analysis (COBRA) platform (Becker et al, 2007). COBRA models use network stoichiometry and steady-state mass balances to define a solution space of potential flux states that a network can take. Thus, the COBRA approach does not require kinetic parameters.

Recently, the global human metabolic network, Recon 1, has been reconstructed (Duarte et al, 2007). To understand the metabolic host–pathogen integrations of M. tb with its human host, we first tailored the global human metabolic network into a cell-specific metabolic reconstruction of the human alveolar macrophage. This was carried out using established computational algorithms (Becker and Palsson, 2008; Shlomi et al, 2008) and manual curation to confirm the included and excluded reactions. The human alveolar macrophage reconstruction, iAB-AMØ-1410, accounts for 1410 genes, 3012 intracellular reactions, and 2572 metabolites (Figure 4C). iAB-AMØ-1410 was able to accurately predict maximum ATP and NO production rates obtained from experimental data (Griscavage et al, 1993; Newsholme et al, 1999).

The second step to studying host–pathogen interactions was integration of the human alveolar macrophage reconstruction with an existing genome-scale metabolic model of M. tb, iNJ661 (Jamshidi and Palsson, 2007). Interfacial constraints were set to create a phagosomal environment that was hypoxic, nitrosative, rich in fatty acids, and poor in carbohydrates. From the onset, it was apparent that some oxygen (<15% of in vitro uptake) was required for proper simulations. In addition, algorithmic tailoring of the M. tb biomass objective function was performed to better represent an infectious state. The integrated host–pathogen metabolic reconstruction was dubbed iAB-AMØ-1410-Mt-661.

Analysis of the integrated host–pathogen metabolic reconstruction resulted in three main findings. First, by setting interfacial constraints and tailoring the biomass objective function, the solution space better represents an infectious state. Without adding artificial constraints to the host portion of the integrated model, the iAB-AMØ-1410 solution space is greatly reduced (Figure 4B). Macrophage glycolysis and nitric oxide production are up-regulated and macrophage ATP production, nucleotide synthesis, and amino-acid metabolism are suppressed. In addition, M. tb glycolysis is suppressed and isocitrate lyase is up-regulated for generation of acetyl-CoA. Fatty acid oxidation pathways and production of mycolic acids are increased, while production of nucleotides, peptidoglycans, and phenolic glycolipids are reduced. The modified solution space of the alveolar macrophage and M. tb better represents the infectious state.

Second, the host-pathogen model more accurately predicts M. tb gene deletion tests than the current in vitro model, iNJ661. The host-pathogen model predicted 11 essential genes and 37 unessential genes differently than iNJ661. A total of 22 of the differentially predicted genes have been experimentally characterized (Sassetti and Rubin, 2003; Sohaskey, 2008). The host-pathogen model correctly predicted 18 of the 22 genes. Thus, iAB-AMØ-1410-Mt-661 is a more accurate platform for studying infectious states of M. tb.

Finally, we sought to determine metabolic differences in both the macrophage and M. tb between three different types of infection: latent, pulmonary, and meningeal. Transcription profiling data of the macrophage for the three infections (Thuong et al, 2008) were integrated in the context of the host–pathogen network to elucidate the reaction activity of the three infections. There was wide heterogeneity in the three infection states; some of these differences are highlighted. Macrophage hyaluronan synthase and export were only active in the pulmonary infection. This is potentially interesting from a pharmaceutical viewpoint as hyaluronan has been implicated as a potential carbon source for extracellular M. tb (Hirayama et al, 2009). In addition, we detected metabolic activity differences in M. tb pathways that have been previously discussed as potential drug targets (Eoh et al, 2007; Boshoff et al, 2008). Polyprenyl metabolic reactions were only active in the latent state infection, while de novo synthesis of nicotinamide cofactors was only active in latent and meningeal M. tb infections.

Host-pathogen modeling represents a novel approach for studying metabolic interactions during infection. iAB-AMØ-1410-Mt-661 is a more accurate platform for understanding the biology and pathophysiology of M. tb infection. Most importantly, genome-scale metabolic reconstructions can act as scaffolds for integrating high-throughput data. Particularly, in this study we were able to discern reaction activity differences between different infection types.

Metabolic coupling of Mycobacterium tuberculosis to its host is foundational to its pathogenesis. Computational genome-scale metabolic models have shown utility in integrating -omic as well as physiologic data for systemic, mechanistic analysis of metabolism. To date, integrative analysis of host–pathogen interactions using in silico mass-balanced, genome-scale models has not been performed. We, therefore, constructed a cell-specific alveolar macrophage model, iAB-AMØ-1410, from the global human metabolic reconstruction, Recon 1. The model successfully predicted experimentally verified ATP and nitric oxide production rates in macrophages. This model was then integrated with an M. tuberculosis H37Rv model, iNJ661, to build an integrated host–pathogen genome-scale reconstruction, iAB-AMØ-1410-Mt-661. The integrated host–pathogen network enables simulation of the metabolic changes during infection. The resulting reaction activity and gene essentiality targets of the integrated model represent an altered infectious state. High-throughput data from infected macrophages were mapped onto the host–pathogen network and were able to describe three distinct pathological states. Integrated host–pathogen reconstructions thus form a foundation upon which understanding the biology and pathophysiology of infections can be developed.

Structural enzymology of sulphur metabolism in Mycobacterium tuberculosis

The emergence of multidrug-resistant strains of Mycobacterium tuberculosis poses a serious threat to human health and has led to world-wide efforts focusing on the development of novel vaccines and antibiotics against this pathogen. Sulphur metabolism in this organism has been linked to essential processes such as virulence and redox defence. The cysteine biosynthetic pathway is up-regulated in models of persistent M. tuberculosis infections and provides potential targets for novel anti-mycobacterial agents, directed specifically toward the pathogen in its persistent phase. Functional and structural characterization of enzymes from sulfur metabolism establishes a necessary framework for the design of strong binding inhibitors that might be developed into new drugs. This review summarizes recent progress in the elucidation of the structural enzymology of the sulphate reduction and cysteine biosynthesis pathways.

3'-Phosphoadenosine-5'-phosphate phosphatase activity is required for superoxide stress tolerance in Streptococcus mutans

Aerobic microorganisms have evolved different strategies to withstand environmental oxidative stresses generated by various reactive oxygen species (ROS). For the facultative anaerobic human oral pathogen Streptococcus mutans, the mechanisms used to protect against ROS are not fully understood, since it does not possess catalase, an enzyme that degrades hydrogen peroxide. In order to elucidate the genes that are essential for superoxide stress response, methyl viologen (MV)-sensitive mutants of S. mutans were generated via ISS1 mutagenesis. Screening of approximately 2,500 mutants revealed six MV-sensitive mutants, each containing an insertion in one of five genes, including a highly conserved hypothetical gene, SMU.1297. Sequence analysis suggests that SMU.1297 encodes a hypothetical protein with a high degree of homology to the Bacillus subtilis YtqI protein, which possesses an oligoribonuclease activity that cleaves nano-RNAs and a phosphatase activity that degrades 3'-phosphoadenosine-5'-phosphate (pAp) and 3'-phosphoadenosine-5'-phosphosulfate (pApS) to produce AMP; the latter activity is similar to the activity of the Escherichia coli CysQ protein, which is required for sulfur assimilation. SMU.1297 was deleted using a markerless Cre-loxP-based strategy; the SMU.1297 deletion mutant was just as sensitive to MV as the ISS1 insertion mutant. Complementation of the deletion mutant with wild-type SMU.1297, in trans, restored the parental phenotype. Biochemical analyses with purified SMU.1297 protein demonstrated that it has pAp phosphatase activity similar to that of YtqI but apparently lacks an oligoribonuclease activity. The ability of SMU.1297 to dephosphorylate pApS in vivo was confirmed by complementation of an E. coli cysQ mutant with SMU.1297 in trans. Thus, our results suggest that SMU.1297 is involved in superoxide stress tolerance in S. mutans. Furthermore, the distribution of homologs of SMU.1297 in streptococci indicates that this protein is essential for superoxide stress tolerance in these organisms.

Rv2131c from Mycobacterium tuberculosis is a CysQ 3'-phosphoadenosine-5'-phosphatase

Mycobacterium tuberculosis (Mtb) produces a number of sulfur-containing metabolites that contribute to its pathogenesis and ability to survive in the host. These metabolites are products of the sulfate assimilation pathway. CysQ, a 3′-phosphoadenosine-5′-phosphatase, is considered an important regulator of this pathway in plants, yeast, and other bacteria. By controlling the pools of 3′-phosphoadenosine 5′-phosphate (PAP) and 3′-phosphoadenosine 5′-phosphosulfate (PAPS), CysQ has the potential to modulate flux in the biosynthesis of essential sulfur-containing metabolites. Bioinformatic analysis of the Mtb genome suggests the presence of a CysQ homologue encoded by the gene Rv2131c. However, a recent biochemical study assigned the protein’s function as a class IV fructose-1,6-bisphosphatase. In the present study, we expressed Rv2131c heterologously and found that the protein dephosphorylates PAP in a magnesium-dependent manner, with optimal activity observed between pH 8.5 and pH 9.5 using 0.5 mM MgCl₂. A sensitive electrospray ionization mass spectrometry-based assay was used to extract the kinetic parameters for PAP, revealing a K_m (8.1 ± 3.1 μM) and k_cat (5.4 ± 1.1 s⁻¹) comparable to those reported for other CysQ enzymes. The second-order rate constant for PAP was determined to be over 3 orders of magnitude greater than those determined for myo-inositol 1-phosphate (IMP) and fructose 1,6-bisphosphate (FBP), previously considered to be the primary substrates of this enzyme. Moreover, the ability of the Rv2131c-encoded enzyme to dephosphorylate PAP and PAPS in vivo was confirmed by functional complementation of an Escherichia coli ΔcysQ mutant. Taken together, these studies indicate that Rv2131c encodes a CysQ enzyme that may play a role in mycobacterial sulfur metabolism.

Significance analysis of microarray for relative quantitation of LC/MS data in proteomics

Background

Although fold change is a commonly used criterion in quantitative proteomics for differentiating regulated proteins, it does not provide an estimation of false positive and false negative rates that is often desirable in a large-scale quantitative proteomic analysis. We explore the possibility of applying the Significance Analysis of Microarray (SAM) method (PNAS 98:5116-5121) to a differential proteomics problem of two samples with replicates. The quantitative proteomic analysis was carried out with nanoliquid chromatography/linear iron trap-Fourier transform mass spectrometry. The biological sample model included two Mycobacterium smegmatis unlabeled cell cultures grown at pH 5 and pH 7. The objective was to compare the protein relative abundance between the two unlabeled cell cultures, with an emphasis on significance analysis of protein differential expression using the SAM method. Results using the SAM method are compared with those obtained by fold change and the conventional t-test.

Results

We have applied the SAM method to solve the two-sample significance analysis problem in liquid chromatography/mass spectrometry (LC/MS) based quantitative proteomics. We grew the pH5 and pH7 unlabelled cell cultures in triplicate resulting in 6 biological replicates. Each biological replicate was mixed with a common ¹⁵N-labeled reference culture cells for normalization prior to SDS/PAGE fractionation and LC/MS analysis. For each biological replicate, one center SDS/PAGE gel fraction was selected for triplicate LC/MS analysis. There were 121 proteins quantified in at least 5 of the 6 biological replicates. Of these 121 proteins, 106 were significant in differential expression by the t-test (p < 0.05) based on peptide-level replicates, 54 were significant in differential expression by SAM with Δ = 0.68 cutoff and false positive rate at 5%, and 29 were significant in differential expression by the t-test (p < 0.05) based on protein-level replicates. The results indicate that SAM appears to overcome the false positives one encounters using the peptide-based t-test while allowing for identification of a greater number of differentially expressed proteins than the protein-based t-test.

Conclusion

We demonstrate that the SAM method can be adapted for effective significance analysis of proteomic data. It provides much richer information about the protein differential expression profiles and is particularly useful in the estimation of false discovery rates and miss rates.