Chemical-genetic interaction mapping links carbon metabolism and cell wall structure to tuberculosis drug efficacy.

Current chemotherapy against Mycobacterium tuberculosis (Mtb), an important human pathogen, requires a multidrug regimen lasting several months. While efforts have been made to optimize therapy by exploiting drug­drug synergies, testing new drug combinations in relevant host environments remains arduous. In particular, host environments profoundly affect the bacterial metabolic state and drug efficacy, limiting the accuracy of predictions based on in vitro assays alone. In this study, we utilized conditional Mtb knockdown mutants of essential genes as an experimentally tractable surrogate for drug treatment and probe the relationship between Mtb carbon metabolism and chemical­genetic interactions (CGIs). We examined the antitubercular drugs isoniazid, rifampicin, and moxifloxacin and found that CGIs are differentially responsive to the metabolic state, defining both environment-independent and -dependent interactions. Specifically, growth on the in vivo­relevant carbon source, cholesterol, reduced rifampicin efficacy by altering mycobacterial cell surface lipid composition. We report that a variety of perturbations in cell wall synthesis pathways restore rifampicin efficacy during growth on cholesterol, and that both environment-independent and cholesterol-dependent in vitro CGIs could be leveraged to enhance bacterial clearance in the mouse infection model. Our findings present an atlas of chemical­genetic­environmental interactions that can be used to optimize drug­drug interactions, as well as provide a framework for understanding in vitro correlates of in vivo efficacy.

An improved statistical method to identify chemical-genetic interactions by exploiting concentration-dependence.

Chemical-genetics (C-G) experiments can be used to identify interactions between inhibitory compounds and bacterial genes, potentially revealing the targets of drugs, or other functionally interacting genes and pathways. C-G experiments involve constructing a library of hypomorphic strains with essential genes that can be knocked-down, treating it with an inhibitory compound, and using high-throughput sequencing to quantify changes in relative abundance of individual mutants. The hypothesis is that, if the target of a drug or other genes in the same pathway are present in the library, such genes will display an excessive fitness defect due to the synergy between the dual stresses of protein depletion and antibiotic exposure. While assays at a single drug concentration are susceptible to noise and can yield false-positive interactions, improved detection can be achieved by requiring that the synergy between gene and drug be concentration-dependent. We present a novel statistical method based on Linear Mixed Models, called CGA-LMM, for analyzing C-G data. The approach is designed to capture the dependence of the abundance of each gene in the hypomorph library on increasing concentrations of drug through slope coefficients. To determine which genes represent candidate interactions, CGA-LMM uses a conservative population-based approach in which genes with negative slopes are considered significant only if they are outliers with respect to the rest of the population (assuming that most genes in the library do not interact with a given inhibitor). We applied the method to analyze 3 independent hypomorph libraries of M. tuberculosis for interactions with antibiotics with anti-tubercular activity, and we identify known target genes or expected interactions for 7 out of 9 drugs where relevant interacting genes are known.

Genome-wide gene expression tuning reveals diverse vulnerabilities of M. tuberculosis.

Antibacterial agents target the products of essential genes but rarely achieve complete target inhibition. Thus, the all-or-none definition of essentiality afforded by traditional genetic approaches fails to discern the most attractive bacterial targets: those whose incomplete inhibition results in major fitness costs. In contrast, gene "vulnerability" is a continuous, quantifiable trait that relates the magnitude of gene inhibition to the effect on bacterial fitness. We developed a CRISPR interference-based functional genomics method to systematically titrate gene expression in Mycobacterium tuberculosis (Mtb) and monitor fitness outcomes. We identified highly vulnerable genes in various processes, including novel targets unexplored for drug discovery. Equally important, we identified invulnerable essential genes, potentially explaining failed drug discovery efforts. Comparison of vulnerability between the reference and a hypervirulent Mtb isolate revealed incomplete conservation of vulnerability and that differential vulnerability can predict differential antibacterial susceptibility. Our results quantitatively redefine essential bacterial processes and identify high-value targets for drug development.

Mycobacterial STAND adenylyl cyclases: The HTH domain binds DNA to form biocrystallized nucleoids.

Mycobacteria harbor a unique class of adenylyl cyclases with a complex domain organization consisting of an N-terminal putative adenylyl cyclase domain fused to a nucleotide-binding adaptor shared by apoptotic protease-activating factor-1, plant resistance proteins, and CED-4 (NB-ARC) domain, a tetratricopeptide repeat (TPR) domain, and a C-terminal helix-turn-helix (HTH) domain. The products of the rv0891c-rv0890c genes represent a split gene pair, where Rv0891c has sequence similarity to adenylyl cyclases, and Rv0890c harbors the NB-ARC-TPR-HTH domains. Rv0891c had very low adenylyl cyclase activity so it could represent a pseudoenzyme. By analyzing the genomic locus, we could express and purify Rv0890c and find that the NB-ARC domain binds ATP and ADP, but does not hydrolyze these nucleotides. Using systematic evolution of ligands by exponential enrichment (SELEX), we identified DNA sequences that bound to the HTH domain of Rv0890c. Uniquely, the HTH domain could also bind RNA. Atomic force microscopy revealed that binding of Rv0890c to DNA was sequence independent, and binding of adenine nucleotides to the protein induced the formation of higher order structures that may represent biocrystalline nucleoids. This represents the first characterization of this group of proteins and their unusual biochemical properties warrant further studies into their physiological roles in future.

Depletion of the DarG antitoxin in Mycobacterium tuberculosis triggers the DNA-damage response and leads to cell death.

Of the ~80 putative toxin-antitoxin (TA) modules encoded by the bacterial pathogen Mycobacterium tuberculosis (Mtb), three contain antitoxins essential for bacterial viability. One of these, Rv0060 (DNA ADP-ribosyl glycohydrolase, DarGMtb ), functions along with its cognate toxin Rv0059 (DNA ADP-ribosyl transferase, DarTMtb ), to mediate reversible DNA ADP-ribosylation (Jankevicius et al., 2016). We demonstrate that DarTMtb -DarGMtb form a functional TA pair and essentiality of darGMtb is dependent on the presence of darTMtb , but simultaneous deletion of both darTMtb -darGMtb does not alter viability of Mtb in vitro or in mice. The antitoxin, DarGMtb , forms a cytosolic complex with DNA-repair proteins that assembles independently of either DarTMtb or interaction with DNA. Depletion of DarGMtb alone is bactericidal, a phenotype that is rescued by expression of an orthologous antitoxin, DarGTaq , from Thermus aquaticus. Partial depletion of DarGMtb triggers a DNA-damage response and sensitizes Mtb to drugs targeting DNA metabolism and respiration. Induction of the DNA-damage response is essential for Mtb to survive partial DarGMtb -depletion and leads to a hypermutable phenotype.

Statistical analysis of variability in TnSeq data across conditions using zero-inflated negative binomial regression.

BACKGROUND: Deep sequencing of transposon mutant libraries (or TnSeq) is a powerful method for probing essentiality of genomic loci under different environmental conditions. Various analytical methods have been described for identifying conditionally essential genes whose tolerance for insertions varies between two conditions. However, for large-scale experiments involving many conditions, a method is needed for identifying genes that exhibit significant variability in insertions across multiple conditions. RESULTS: In this paper, we introduce a novel statistical method for identifying genes with significant variability of insertion counts across multiple conditions based on Zero-Inflated Negative Binomial (ZINB) regression. Using likelihood ratio tests, we show that the ZINB distribution fits TnSeq data better than either ANOVA or a Negative Binomial (in a generalized linear model). We use ZINB regression to identify genes required for infection of M. tuberculosis H37Rv in C57BL/6 mice. We also use ZINB to perform a analysis of genes conditionally essential in H37Rv cultures exposed to multiple antibiotics. CONCLUSIONS: Our results show that, not only does ZINB generally identify most of the genes found by pairwise resampling (and vastly out-performs ANOVA), but it also identifies additional genes where variability is detectable only when the magnitudes of insertion counts are treated separately from local differences in saturation, as in the ZINB model.

Autoinhibitory mechanism and activity-related structural changes in a mycobacterial adenylyl cyclase.

An adenylyl cyclase from Mycobacterium avium, Ma1120, is a functional orthologue of a pseudogene Rv1120c from Mycobacterium tuberculosis. We report the crystal structure of Ma1120 in a monomeric form and its truncated construct as a dimer. Ma1120 exists as a monomer in solution and crystallized as a monomer in the absence of substrate or inhibitor. An additional α-helix present at the N-terminus of the monomeric structure blocks the active site by interacting with the substrate binding residues and occupying the dimer interface region. However, the enzyme has been found to be active in solution, indicating the movement of the helix away from the interface to facilitate the formation of active dimers in conditions favourable for catalysis. Thus, the N-terminal helix of Ma1120 keeps the enzyme in an autoinhibited state when it is not active. Deletion of this helix enabled us to crystallize the molecule as an active homodimer in the presence of a P-site inhibitor 2',5'-dideoxy-3'-ATP, or pyrophosphate along with metal ions. The substrate specifying lysine residue plays a dual role of interacting with the substrate and stabilizing the dimer. The dimerization loop region harbouring the second substrate specifying residue, an aspartate, shows significant differences in conformation and position between the monomeric and dimeric structures. Thus, this study has not only revealed that significant structural transitions are required for the interconversion of the inactive and the active forms of the enzyme, but also provided precise nature of these transitions.

The adenylyl cyclase Rv2212 modifies the proteome and infectivity of Mycobacterium bovis BCG.

All organisms have the capacity to sense and respond to environmental changes. These signals often involve the use of second messengers such as cyclic adenosine monophosphate (cAMP). This second messenger is widely distributed among organisms and coordinates gene expression related with pathogenesis, virulence, and environmental adaptation. Genomic analysis in Mycobacterium tuberculosis has identified 16 adenylyl cyclases (AC) and one phosphodiesterase, which produce and degrade cAMP, respectively. To date, ten AC have been biochemically characterized and only one (Rv0386) has been found to be important during murine infection with M. tuberculosis. Here, we investigated the impact of hsp60-driven Rv2212 gene expression in Mycobacterium bovis Bacillus Calmette-Guerin (BCG) during growth in vitro, and during macrophage and mice infection. We found that hsp60-driven expression of Rv2212 resulted in an increased capacity of replication in murine macrophages but an attenuated phenotype in lungs and spleen when administered intravenously in mice. Furthermore, this strain displayed an altered proteome mainly affecting proteins associated with stress conditions (bfrB, groEL-2, DnaK) that could contribute to the attenuated phenotype observed in mice.

Paralogous cAMP receptor proteins in Mycobacterium smegmatis show biochemical and functional divergence.

The cyclic AMP receptor protein (CRP) family of transcription factors consists of global regulators of bacterial gene expression. Here, we identify two paralogous CRPs in the genome of Mycobacterium smegmatis that have 78% identical sequences and characterize them biochemically and functionally. The two proteins (MSMEG_0539 and MSMEG_6189) show differences in cAMP binding affinity, trypsin sensitivity, and binding to a CRP site that we have identified upstream of the msmeg_3781 gene. MSMEG_6189 binds to the CRP site readily in the absence of cAMP, while MSMEG_0539 binds in the presence of cAMP, albeit weakly. msmeg_6189 appears to be an essential gene, while the Δmsmeg_0539 strain was readily obtained. Using promoter-reporter constructs, we show that msmeg_3781 is regulated by CRP binding, and its transcription is repressed by MSMEG_6189. Our results are the first to characterize two paralogous and functional CRPs in a single bacterial genome. This gene duplication event has subsequently led to the evolution of two proteins whose biochemical differences translate to differential gene regulation, thus catering to the specific needs of the organism.

A fluorescent nucleic acid nanodevice quantitatively images elevated cyclic adenosine monophosphate in membrane-bound compartments.

cAMPhor: In the presence of cAMP, cAMPhor folds into a structure that binds DFHBI (green), increasing its fluorescence, while Alexa 647 (red) functions as a normalizing dye. It can thus be used to spatially image cAMP quantitatively in membrane-bound compartments.