Pyrazinamide Resistance Is Caused by Two Distinct Mechanisms: Prevention of Coenzyme A Depletion and Loss of Virulence Factor Synthesis

Propionate prevents loss of the PDIM virulence lipid in Mycobacterium tuberculosis

Phthiocerol dimycocerosate (PDIM) is an essential virulence lipid of Mycobacterium tuberculosis. In vitro culturing rapidly selects for spontaneous PDIM-negative mutants that have attenuated virulence and increased cell wall permeability, thus impacting the relevance of experimental findings. PDIM loss can also reduce the efficacy of the BCG Pasteur vaccine. Here we show that vancomycin susceptibility can rapidly screen for M. tuberculosis PDIM production. We find that metabolic deficiency of methylmalonyl-CoA impedes the growth of PDIM-producing bacilli, selecting for PDIM-negative variants. Supplementation with odd-chain fatty acids, cholesterol or vitamin B12 restores PDIM-positive bacterial growth. Specifically, we show that propionate supplementation enhances PDIM-producing bacterial growth and selects against PDIM-negative mutants, analogous to in vivo conditions. Our study provides a simple approach to screen for and maintain PDIM production, and reveals how discrepancies between the host and in vitro nutrient environments can attenuate bacterial pathogenicity.

Prediction of pyrazinamide resistance in Mycobacterium tuberculosis using structure-based machine-learning approaches

Background

Pyrazinamide is one of four first-line antibiotics used to treat tuberculosis; however, antibiotic susceptibility testing for pyrazinamide is challenging. Resistance to pyrazinamide is primarily driven by genetic variation in pncA, encoding an enzyme that converts pyrazinamide into its active form.

Methods

We curated a dataset of 664 non-redundant, missense amino acid mutations in PncA with associated high-confidence phenotypes from published studies and then trained three different machine-learning models to predict pyrazinamide resistance. All models had access to a range of protein structural-, chemical- and sequence-based features.

Results

The best model, a gradient-boosted decision tree, achieved a sensitivity of 80.2% and a specificity of 76.9% on the hold-out test dataset. The clinical performance of the models was then estimated by predicting the binary pyrazinamide resistance phenotype of 4027 samples harbouring 367 unique missense mutations in pncA derived from 24 231 clinical isolates.

Conclusions

This work demonstrates how machine learning can enhance the sensitivity/specificity of pyrazinamide resistance prediction in genetics-based clinical microbiology workflows, highlights novel mutations for future biochemical investigation, and is a proof of concept for using this approach in other drugs.

Recombinant expression and functional characterization of FadD2 protein in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) is a crucial and highly destructive intracellular pathogen responsible for causing tuberculosis (TB). The emergence and dissemination of multi-drug resistant Mtb has further aggravated the TB crisis, leading to high mortality. Mtb FadD2 is a fatty acyl-coenzyme A (CoA) synthetase that modifies the cell envelope and plays an important role in reducing Mtb susceptibility to pyrazinoic acid (POA). However, the functional mechanism of Mtb FadD2 remains poorly understood. Here, we successfully expressed, purified and obtained monomeric FadD2 by using buffer (500 mM NaCl, 20 mM Tris-HCl, pH7.4 and 5 % glycerol). Palmitate was found to be the optimal substrate for FadD2. Fatty acyl-CoA synthetase activity reached maximum at 450 μM palmitate, and the Km value was 318.2 μM for palmitate. The results of mutation experiments indicated the critical role of T370 and K551 in the enzymatic activity of FadD2. Our work provides a guideline and concept for the development of novel drugs against Mtb.

Cholesterol-Mediated Coenzyme A Depletion in Catabolic Mutants of Mycobacteria Leads to Toxicity

Cholesterol is a critical growth substrate for Mycobacterium tuberculosis (Mtb) during infection, and the cholesterol catabolic pathway has been targeted for the development of new antimycobacterial agents. A key metabolite in cholesterol catabolism is 3aα-H-4α(3'-propanoate)-7aβ-methylhexahydro-1,5-indanedione (HIP). Many of the HIP metabolites are acyl-coenzyme A (CoA) thioesters, whose accumulation in deletion mutants can cause cholesterol-mediated toxicity. We used LC-MS/MS analysis to demonstrate that deletion of genes involved in HIP catabolism leads to acyl-CoA accumulation with concomitant depletion of free CoASH, leading to dysregulation of central metabolic pathways. CoASH and acyl-CoAs inhibited PanK, the enzyme that catalyzes the first step in the transformation of pantothenate to CoASH. Inhibition was competitive with respect to ATP with Kic values ranging from 9 μM for CoASH to 57 μM for small acyl-CoAs and 180 ± 30 μM for cholesterol-derived acyl-CoA. These findings link two critical metabolic pathways and suggest that therapeutics targeting cholesterol catabolic enzymes could both prevent the utilization of an important growth substrate and simultaneously sequester CoA from essential cellular processes, leading to bacterial toxicity.

Inhibition of host PARP1 contributes to the anti-inflammatory and antitubercular activity of pyrazinamide

The antibiotic pyrazinamide (PZA) is a cornerstone of tuberculosis (TB) therapy that shortens treatment durations by several months despite being only weakly bactericidal. Intriguingly, PZA is also an anti-inflammatory molecule shown to specifically reduce inflammatory cytokine signaling and lesion activity in TB patients. However, the target and clinical importance of PZA's host-directed activity during TB therapy remain unclear. Here, we identify the host enzyme Poly(ADP-ribose) Polymerase 1 (PARP1), a pro-inflammatory master regulator strongly activated in TB, as a functionally relevant host target of PZA. We show that PZA inhibits PARP1 enzymatic activity in macrophages and in mice where it reverses TB-induced PARP1 activity in lungs to uninfected levels. Utilizing a PZA-resistant mutant, we demonstrate that PZA's immune-modulatory effects are PARP1-dependent but independent of its bactericidal activity. Importantly, PZA's bactericidal efficacy is impaired in PARP1-deficient mice, suggesting that immune modulation may be an integral component of PZA's antitubercular activity. In addition, adjunctive PARP1 inhibition dramatically reduces inflammation and lesion size in mice and may be a means to reduce lung damage and shorten TB treatment duration. Together, these findings provide insight into PZA's mechanism of action and the therapeutic potential of PARP1 inhibition in the treatment of TB.

The PDIM paradox of Mycobacterium tuberculosis: new solutions to a persistent problem

Phthiocerol dimycocerosate (PDIM) is an essential virulence lipid of Mycobacterium tuberculosis. In vitro culturing rapidly selects for spontaneous mutations that cause PDIM loss leading to virulence attenuation and increased cell wall permeability. We discovered that PDIM loss is due to a metabolic deficiency of methylmalonyl-CoA that impedes the growth of PDIM-producing bacilli. This can be remedied by supplementation with odd-chain fatty acids, cholesterol, or vitamin B12. We developed a much-needed facile and scalable routine assay for PDIM production and show that propionate supplementation enhances the growth of PDIM-producing bacilli and selects against PDIM-negative mutants, analogous to in vivo conditions. Our results solve a major issue in tuberculosis research and exemplify how discrepancies between the host and in vitro nutrient environments can attenuate bacterial pathogenicity.

Pharmacokinetic Considerations for Optimizing Inhaled Spray-Dried Pyrazinoic Acid Formulations

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), remains a leading cause of death with 1.6 million deaths worldwide reported in 2021. Oral pyrazinamide (PZA) is an integral part of anti-TB regimens, but its prolonged use has the potential to drive the development of PZA-resistant Mtb. PZA is converted to the active moiety pyrazinoic acid (POA) by the Mtb pyrazinamidase encoded by pncA, and mutations in pncA are associated with the majority of PZA resistance. Conventional oral and parenteral therapies may result in subtherapeutic exposure in the lung; hence, direct pulmonary administration of POA may provide an approach to rescue PZA efficacy for treating pncA-mutant PZA-resistant Mtb. The objectives of the current study were to (i) develop novel dry powder POA formulations, (ii) assess their feasibility for pulmonary delivery using physicochemical characterization, (iii) evaluate their pharmacokinetics (PK) in the guinea pig model, and (iv) develop a mechanism-based pharmacokinetic model (MBM) using in vivo PK data to select a formulation providing adequate exposure in epithelial lining fluid (ELF) and lung tissue. We developed three POA formulations for pulmonary delivery and characterized their PK in plasma, ELF, and lung tissue following passive inhalation in guinea pigs. Additionally, the PK of POA following oral, intravenous, and intratracheal administration was characterized in guinea pigs. The MBM was used to simultaneously model PK data following administration of POA and its formulations via the different routes. The MBM described POA PK well in plasma, ELF, and lung tissue. Physicochemical analyses and MBM predictions suggested that POA maltodextrin was the best among the three formulations and an excellent candidate for further development as it has: (i) the highest ELF-to-plasma exposure ratio (203) and lung tissue-to-plasma exposure ratio (30.4) compared with POA maltodextrin and leucine (75.7/16.2) and POA leucine salt (64.2/19.3) and (ii) the highest concentration in ELF (CmaxELF: 171 nM) within 15.5 min, correlating with a fast transfer into ELF after pulmonary administration (KPM: 22.6 1/h). The data from the guinea pig allowed scaling, using the MBM to a human dose of POA maltodextrin powder demonstrating the potential feasibility of an inhaled product.

PROTAC antibiotics: the time is now

INTRODUCTION

Novel antibiotics are needed to keep antibiotic resistance at bay and to improve treatment of the many drug-susceptible infections for which current therapies achieve poor cure rates. While revolutionizing human therapeutics, the concept of targeted protein degradation (TPD) by bifunctional proteolysis targeting chimeras (PROTACs) has not yet been applied to the discovery of antibiotics. A major obstacle precluding successful translation of this strategy to antibiotic development is that bacteria lack the E3 ligase-proteasome system exploited by human PROTACs to facilitate target degradation.

AREAS COVERED

The authors describe the serendipitous discovery of the first monofunctional target-degrading antibiotic pyrazinamide, supporting TPD as a viable and novel approach in antibiotic discovery. They then discuss the rational design, mechanism, and activity of the first bifunctional antibacterial target degrader BacPROTAC, enabling a generalizable approach to TPD in bacteria.

EXPERT OPINION

BacPROTACs demonstrate that linking a target directly to a bacterial protease complex can promote target degradation. BacPROTACs successfully bypass the 'middleman' E3 ligase, providing an entry strategy for the generation of antibacterial PROTACs. We speculate that antibacterial PROTACs will not only expand the target space but may also improve treatment by allowing dosage reduction, stronger bactericidal activity and activity against drug-tolerant 'persisters.'

Structure activity relationship of pyrazinoic acid analogs as potential antimycobacterial agents

Tuberculosis (TB) remains a leading cause of infectious disease-related mortality and morbidity. Pyrazinamide (PZA) is a critical component of the first-line TB treatment regimen because of its sterilizing activity against non-replicating Mycobacterium tuberculosis (Mtb), but its mechanism of action has remained enigmatic. PZA is a prodrug converted by pyrazinamidase encoded by pncA within Mtb to the active moiety, pyrazinoic acid (POA) and PZA resistance is caused by loss-of-function mutations to pyrazinamidase. We have recently shown that POA induces targeted protein degradation of the enzyme PanD, a crucial component of the coenzyme A biosynthetic pathway essential in Mtb. Based on the newly identified mechanism of action of POA, along with the crystal structure of PanD bound to POA, we designed several POA analogs using structure for interpretation to improve potency and overcome PZA resistance. We prepared and tested ring and carboxylic acid bioisosteres as well as 3, 5, 6 substitutions on the ring to study the structure activity relationships of the POA scaffold. All the analogs were evaluated for their whole cell antimycobacterial activity, and a few representative molecules were evaluated for their binding affinity, towards PanD, through isothermal titration calorimetry. We report that analogs with ring and carboxylic acid bioisosteres did not significantly enhance the antimicrobial activity, whereas the alkylamino-group substitutions at the 3 and 5 position of POA were found to be up to 5 to 10-fold more potent than POA. Further development and mechanistic analysis of these analogs may lead to a next generation POA analog for treating TB.

Mutations Associated with Pyrazinamide Resistance in Mycobacterium tuberculosis: A Review and Update

Pyrazinamide (PZA) has remained a keystone of tuberculosis (TB) therapy, and it possesses high imperative sterilizing action that can facilitate reduction in the present chemotherapy regimen. The combination of PZA works both with first- and second-line TB drugs, notably fluoroquinolones, clofazimine, bedaquiline, delamanid and pretomanid. Pyrazinamide inhibits various targets that are involved in different cellular processes like energy production (pncA), trans-translation (rpsA) and pantothenate/coenzyme A (panD) which are required for persistence of the pathogen. It is well known that pncA gene encoding pyrazinamidase is involved in the transition of PZA into the active form of pyrazinoic acid, which implies that mutation in the pncA gene can develop PZA resistance in Mycobacterium tuberculosis (M. tuberculosis) strain leading to a major clinical and public health concern. Therefore, it is very crucial to understand its resistance mechanism and to detect it precisely to help in the management of the disease. Scope of this review is to have a deep understanding of molecular mechanism of PZA resistance with its multiple targets which would help study the association of mutations and its resistance in M. tuberculosis. This will in turn help learn about the resistance of PZA and develop more accurate molecular diagnostic tool for drug-resistant TB in future TB therapy.

Structural and Mechanistic Insights into Mycobacterium abscessus Aspartate Decarboxylase PanD and a Pyrazinoic Acid-Derived Inhibitor

Mycobacterium tuberculosis (Mtb) aspartate decarboxylase PanD is required for biosynthesis of the essential cofactor coenzyme A and targeted by the first line drug pyrazinamide (PZA). PZA is a prodrug that is converted by a bacterial amidase into its bioactive form pyrazinoic acid (POA). Employing structure-function analyses we previously identified POA-based inhibitors of Mtb PanD showing much improved inhibitory activity against the enzyme. Here, we performed the first structure-function studies on PanD encoded by the nontuberculous mycobacterial lung pathogen Mycobacterium abscessus (Mab), shedding light on the differences and similarities of Mab and Mtb PanD. Solution X-ray scattering data provided the solution structure of the entire tetrameric Mab PanD, which in comparison to the structure of the derived C-terminal truncated Mab PanD1-114 mutant revealed the orientation of the four flexible C-termini relative to the catalytic core. Enzymatic studies of Mab PanD1-114 explored the essentiality of the C-terminus for catalysis. A library of recombinant Mab PanD mutants based on structural information and PZA/POA resistant PanD mutations in Mtb illuminated critical residues involved in the substrate tunnel and enzymatic activity. Using our library of POA analogues, we identified (3-(1-naphthamido)pyrazine-2-carboxylic acid) (analogue 2) as the first potent inhibitor of Mab PanD. The inhibitor shows mainly electrostatic- and hydrogen bonding interaction with the target enzyme as explored by isothermal titration calorimetry and confirmed by docking studies. The observed unfavorable entropy indicates that significant conformational changes are involved in the binding process of analogue 2 to Mab PanD. In contrast to PZA and POA, which are whole-cell inactive, analogue 2 exerts appreciable antibacterial activity against the three subspecies of Mab.

Visualizing Pyrazinamide Action by Live Single-Cell Imaging of Phagosome Acidification and Mycobacterium tuberculosis pH Homeostasis

Mycobacterium tuberculosis segregates within multiple subcellular niches with different biochemical and biophysical properties that, upon treatment, may impact antibiotic distribution, accumulation, and efficacy. However, it remains unclear whether fluctuating intracellular microenvironments alter mycobacterial homeostasis and contribute to antibiotic enrichment and efficacy. Here, we describe a live dual-imaging approach to monitor host subcellular acidification and M. tuberculosis intrabacterial pH. By combining this approach with pharmacological and genetic perturbations, we show that M. tuberculosis can maintain its intracellular pH independently of the surrounding pH in human macrophages. Importantly, unlike bedaquiline (BDQ), isoniazid (INH), or rifampicin (RIF), the drug pyrazinamide (PZA) displays antibacterial efficacy by disrupting M. tuberculosis intrabacterial pH homeostasis in cellulo. By using M. tuberculosis mutants, we confirmed that intracellular acidification is a prerequisite for PZA efficacy in cellulo. We anticipate this imaging approach will be useful to identify host cellular environments that affect antibiotic efficacy against intracellular pathogens. IMPORTANCE We still do not completely understand why tuberculosis (TB) treatment requires the combination of several antibiotics for up to 6 months. M. tuberculosis is a facultative intracellular pathogen, and it is still unknown whether heterogenous and dynamic intracellular populations of bacteria in different cellular environments affect antibiotic efficacy. By developing a dual live imaging approach to monitor mycobacterial pH homeostasis, host cell environment, and antibiotic action, we show here that intracellular localization of M. tuberculosis affects the efficacy of one first-line anti-TB drug. Our observations can be applicable to the treatment of other intracellular pathogens and help to inform the development of more effective combined therapies for tuberculosis that target heterogenous bacterial populations within the host.

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

Efforts to improve tuberculosis therapy include optimizing multidrug regimens to take advantage of drug–drug synergies. However, the complex host environment has a profound effect on bacterial metabolic state and drug activity, making predictions of optimal drug combinations difficult. In this study, we leverage a newly developed library of conditional knockdown Mycobacterium tuberculosis mutants in which genetic depletion of essential genes mimics the effect of drug therapy. This tractable system allowed us to assess the effect of growth condition on predicted drug–drug interactions. We found that these interactions can be differentially sensitive to the metabolic state, and select in vitro–defined interactions can be leveraged to accelerate bacterial killing during infection. These findings suggest strategies for optimizing tuberculosis therapy.

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.

Experimental confirmation that an uncommon rrs gene mutation (g878a) of Mycobacterium tuberculosis confers resistance to streptomycin

The effective treatment of patients diagnosed with drug resistant tuberculosis is highly dependent upon the ability to rapidly and accurately determine the antibiotic susceptibility profile of the Mycobacterium tuberculosis isolate(s) involved. Thus, as more clinical microbiology laboratories advance towards the use of DNA sequence-based diagnostics, it is imperative that their predictive functions extend beyond the well-known resistance mutations, in order to also encompass as many of the lower-frequency mutations as possible. However, in most cases, the fundamental experimental proof that links these uncommon mutations with phenotypic resistance is lacking. One such example is the g878a polymorphism within the rrs 16s rRNA gene. We, and others, have identified this mutation within a small number of drug-resistant isolates, although a consensus regarding exactly which aminoglycoside antibiotic(s) it confers resistance toward has not previously been reached. Here we have employed oligo-mediated recombineering to introduce the g878a polymorphism into the rrs gene of M. bovis BCG - a close relative of M. tuberculosis - and demonstrate that it confers low-level resistance to streptomycin alone. It does not confer cross-resistance towards amikacin, capreomycin, nor kanamycin. We also demonstrate that the rrs g878a mutation exerts a substantial fitness defect in vitro , that may at least in part explain why clinical isolates bearing this mutation appear to be quite rare. Overall, this study provides clarity to the phenotype attributable to the rrs g878a mutation and is relevant to the future implementation of genomics-based diagnostics, as well as the clinical management of patients where this particular polymorphism is encountered.

In Vitro Resistance against DNA Gyrase Inhibitor SPR719 in Mycobacterium avium and Mycobacterium abscessus

The aminobenzimidazole SPR719 targets DNA gyrase in Mycobacterium tuberculosis. The molecule acts as inhibitor of the enzyme’s ATPase located on the Gyrase B subunit of the tetrameric Gyrase A₂B₂ protein. SPR719 is also active against non-tuberculous mycobacteria (NTM) and recently entered clinical development for lung disease caused by these bacteria. Resistance against SPR719 in NTM has not been characterized. Here, we determined spontaneous in vitro resistance frequencies in single step resistance development studies, MICs of resistant strains, and resistance associated DNA sequence polymorphisms in two major NTM pathogens Mycobacterium avium and Mycobacterium abscessus. A low-frequency resistance (10^−8/CFU) was associated with missense mutations in the ATPase domain of the Gyrase B subunit in both bacteria, consistent with inhibition of DNA gyrase as the mechanism of action of SPR719 against NTM. For M. abscessus, but not for M. avium, a second, high-frequency (10^−6/CFU) resistance mechanism was observed. High-frequency SPR719 resistance was associated with frameshift mutations in the transcriptional repressor MAB_4384 previously shown to regulate expression of the drug efflux pump system MmpS5/MmpL5. Our results confirm DNA gyrase as target of SPR719 in NTM and reveal differential resistance development in the two NTM species, with M. abscessus displaying high-frequency indirect resistance possibly involving drug efflux.

IMPORTANCE Clinical emergence of resistance to new antibiotics affects their utility. Characterization of in vitro resistance is a first step in the profiling of resistance properties of novel drug candidates. Here, we characterized in vitro resistance against SPR719, a drug candidate for the treatment of lung disease caused by non-tuberculous mycobacteria (NTM). The identified resistance associated mutations and the observed differential resistance behavior of the two characterized NTM species provide a basis for follow-up studies of resistance in vivo to further inform clinical development of SPR719.

Pantothenate and CoA biosynthesis in Apicomplexa and their promise as antiparasitic drug targets

The Apicomplexa phylum comprises thousands of distinct intracellular parasite species, including coccidians, haemosporidians, piroplasms, and cryptosporidia. These parasites are characterized by complex and divergent life cycles occupying a variety of host niches. Consequently, they exhibit distinct adaptations to the differences in nutritional availabilities, either relying on biosynthetic pathways or by salvaging metabolites from their host. Pantothenate (Pan, vitamin B5) is the precursor for the synthesis of an essential cofactor, coenzyme A (CoA), but among the apicomplexans, only the coccidian subgroup has the ability to synthesize Pan. While the pathway to synthesize CoA from Pan is largely conserved across all branches of life, there are differences in the redundancy of enzymes and possible alternative pathways to generate CoA from Pan. Impeding the scavenge of Pan and synthesis of Pan and CoA have been long recognized as potential targets for antimicrobial drug development, but in order to fully exploit these critical pathways, it is important to understand such differences. Recently, a potent class of pantothenamides (PanAms), Pan analogs, which target CoA-utilizing enzymes, has entered antimalarial preclinical development. The potential of PanAms to target multiple downstream pathways make them a promising compound class as broad antiparasitic drugs against other apicomplexans. In this review, we summarize the recent advances in understanding the Pan and CoA biosynthesis pathways, and the suitability of these pathways as drug targets in Apicomplexa, with a particular focus on the cyst-forming coccidian, Toxoplasma gondii, and the haemosporidian, Plasmodium falciparum.

Structure-Aware Mycobacterium tuberculosis Functional Annotation Uncloaks Resistance, Metabolic, and Virulence Genes

Accurate and timely functional genome annotation is essential for translating basic pathogen research into clinically impactful advances. Here, through literature curation and structure-function inference, we systematically update the functional genome annotation of Mycobacterium tuberculosis virulent type strain H37Rv. First, we systematically curated annotations for 589 genes from 662 publications, including 282 gene products absent from leading databases. Second, we modeled 1,711 underannotated proteins and developed a semiautomated pipeline that captured shared function between 400 protein models and structural matches of known function on Protein Data Bank, including drug efflux proteins, metabolic enzymes, and virulence factors. In aggregate, these structure- and literature-derived annotations update 940/1,725 underannotated H37Rv genes and generate hundreds of functional hypotheses. Retrospectively applying the annotation to a recent whole-genome transposon mutant screen provided missing function for 48% (13/27) of underannotated genes altering antibiotic efficacy and 33% (23/69) required for persistence during mouse tuberculosis (TB) infection. Prospective application of the protein models enabled us to functionally interpret novel laboratory generated pyrazinamide (PZA)-resistant mutants of unknown function, which implicated the emerging coenzyme A depletion model of PZA action in the mutants’ PZA resistance. Our findings demonstrate the functional insight gained by integrating structural modeling and systematic literature curation, even for widely studied microorganisms. Functional annotations and protein structure models are available at https://tuberculosis.sdsu.edu/H37Rv in human- and machine-readable formats.

IMPORTANCEMycobacterium tuberculosis, the primary causative agent of tuberculosis, kills more humans than any other infectious bacterium. Yet 40% of its genome is functionally uncharacterized, leaving much about the genetic basis of its resistance to antibiotics, capacity to withstand host immunity, and basic metabolism yet undiscovered. Irregular literature curation for functional annotation contributes to this gap. We systematically curated functions from literature and structural similarity for over half of poorly characterized genes, expanding the functionally annotated Mycobacterium tuberculosis proteome. Applying this updated annotation to recent in vivo functional screens added functional information to dozens of clinically pertinent proteins described as having unknown function. Integrating the annotations with a prospective functional screen identified new mutants resistant to a first-line TB drug, supporting an emerging hypothesis for its mode of action. These improvements in functional interpretation of clinically informative studies underscore the translational value of this functional knowledge. Structure-derived annotations identify hundreds of high-confidence candidates for mechanisms of antibiotic resistance, virulence factors, and basic metabolism and other functions key in clinical and basic tuberculosis research. More broadly, they provide a systematic framework for improving prokaryotic reference annotations.

Whole-Genome SNP Analysis Identifies Putative Mycobacterium bovis Transmission Clusters in Livestock and Wildlife in Catalonia, Spain

The high-resolution WGS analyses of MTBC strains have provided useful insight for determining sources of infection for animal tuberculosis. In Spain, tuberculosis in livestock is caused by Mycobacterium bovis and Mycobacterium caprae, where wildlife reservoirs play an important role. We analyzed a set of 125 M. bovis isolates obtained from livestock and wildlife from Catalonia to investigate strain diversity and identify possible sources and/or causes of infection. Whole-genome SNP profiles were used for phylogenetic reconstruction and pairwise SNP distance analysis. Additionally, SNPs were investigated to identify virulence and antimicrobial resistance factors to investigate clade-specific associations. Putative transmission clusters (≤12 SNPs) were identified, and associated epidemiological metadata were used to determine possible explanatory factors for transmission. M. bovis distribution was heterogeneous, with 7 major clades and 21 putative transmission clusters. In order of importance, the explanatory factors associated were proximity and neighborhood, residual infection, livestock-wildlife interaction, shared pasture, and movement. Genes related to lipid transport and metabolism showed the highest number of SNPs. All isolates were pyrazinamide resistant, and five were additionally resistant to isoniazid, but no clade-specific associations could be determined. Our findings highlight the importance of high-resolution molecular surveillance to monitor bovine tuberculosis dynamics in a low-prevalence setting.

Activity of Pyrazinamide against Mycobacterium tuberculosis at Neutral pH in PZA-S1 Minimal Medium

Susceptibility testing of tuberculosis (TB) drugs on Mycobacterium tuberculosis is essential for the rapid detection of strains resistant to the drugs, providing the patient with effective treatment, and preventing the spread of drug-resistant TB strains. Pyrazinamide (PZA) is one of the first-line agents used for the treatment of TB. However, current phenotypic PZA susceptibility testing is unreliable due to its performance in acidic pH conditions. The aims of this study were to develop minimal media to determine the activity of PZA at a neutral pH at 37 °C to avoid problems caused by an acidic pH, which is currently used in PZA susceptibility tests, and to identify PZA-resistant M. tuberculosis in media with reproducibility and accuracy. Different minimal media were used to determine the activity of PZA using the broth microdilution method with M. tuberculosis H37Ra as the reference strain. The PZA-S1 minimal medium was proposed as the most suitable medium. PZA inhibited the growth of M. tuberculosis in PZA-S1 at a neutral pH of 6.8, which is the optimal pH for M. tuberculosis growth. Moreover, PZA showed activity at a neutral pH on a PZA-S1 agar plate when using the disk diffusion method. PZA-resistant M. tuberculosis could be identified at a neutral pH in PZA-S1 minimal medium. This study establishes valuable information regarding the testing of PZA’s susceptibility in relation to M. tuberculosis at a neutral pH of 6.8 with reliability and accuracy in clinical settings.

Structural and Genomic Insights Into Pyrazinamide Resistance in Mycobacterium tuberculosis Underlie Differences Between Ancient and Modern Lineages

Resistance to drugs used to treat tuberculosis disease (TB) continues to remain a public health burden, with missense point mutations in the underlying Mycobacterium tuberculosis bacteria described for nearly all anti-TB drugs. The post-genomics era along with advances in computational and structural biology provide opportunities to understand the interrelationships between the genetic basis and the structural consequences of M. tuberculosis mutations linked to drug resistance. Pyrazinamide (PZA) is a crucial first line antibiotic currently used in TB treatment regimens. The mutational promiscuity exhibited by the pncA gene (target for PZA) necessitates computational approaches to investigate the genetic and structural basis for PZA resistance development. We analysed 424 missense point mutations linked to PZA resistance derived from ∼35K M. tuberculosis clinical isolates sourced globally, which comprised the four main M. tuberculosis lineages (Lineage 1-4). Mutations were annotated to reflect their association with PZA resistance. Genomic measures (minor allele frequency and odds ratio), structural features (surface area, residue depth and hydrophobicity) and biophysical effects (change in stability and ligand affinity) of point mutations on pncA protein stability and ligand affinity were assessed. Missense point mutations within pncA were distributed throughout the gene, with the majority (>80%) of mutations with a destabilising effect on protomer stability and on ligand affinity. Active site residues involved in PZA binding were associated with multiple point mutations highlighting mutational diversity due to selection pressures at these functionally important sites. There were weak associations between genomic measures and biophysical effect of mutations. However, mutations associated with PZA resistance showed statistically significant differences between structural features (surface area and residue depth), but not hydrophobicity score for mutational sites. Most interestingly M. tuberculosis lineage 1 (ancient lineage) exhibited a distinct protein stability profile for mutations associated with PZA resistance, compared to modern lineages.

Intracellular localisation of Mycobacterium tuberculosis affects efficacy of the antibiotic pyrazinamide

To be effective, chemotherapy against tuberculosis (TB) must kill the intracellular population of the pathogen, Mycobacterium tuberculosis. However, how host cell microenvironments affect antibiotic accumulation and efficacy remains unclear. Here, we use correlative light, electron, and ion microscopy to investigate how various microenvironments within human macrophages affect the activity of pyrazinamide (PZA), a key antibiotic against TB. We show that PZA accumulates heterogeneously among individual bacteria in multiple host cell environments. Crucially, PZA accumulation and efficacy is maximal within acidified phagosomes. Bedaquiline, another antibiotic commonly used in combined TB therapy, enhances PZA accumulation via a host cell-mediated mechanism. Thus, intracellular localisation and specific microenvironments affect PZA accumulation and efficacy. Our results may explain the potent in vivo efficacy of PZA, compared to its modest in vitro activity, and its critical contribution to TB combination chemotherapy.

Mycobacterium tuberculosis PanD Structure-Function Analysis and Identification of a Potent Pyrazinoic Acid-Derived Enzyme Inhibitor

A common strategy employed in antibacterial drug discovery is the targeting of biosynthetic processes that are essential and specific for the pathogen. Specificity in particular avoids undesirable interactions with potential enzymatic counterparts in the human host, and it ensures on-target toxicity. Synthesis of pantothenate (Vitamine B5), which is a precursor of the acyl carrier coenzyme A, is an example of such a pathway. In Mycobacterium tuberculosis (Mtb), which is the causative agent of tuberculosis (TB), pantothenate is formed by pantothenate synthase, utilizing D-pantoate and β-Ala as substrates. β-Ala is mainly formed by the decarboxylation of l-aspartate, generated by the decarboxylase PanD, which is a homo-oliogomer in solution. Pyrazinoic acid (POA), which is the bioactive form of the TB prodrug pyrazinamide, binds and inhibits PanD activity weakly. Here, we generated a library of recombinant Mtb PanD mutants based on structural information and PZA/POA resistance mutants. Alterations in oligomer formation, enzyme activity, and/or POA binding were observed in respective mutants, providing insights into essential amino acids for Mtb PanD's proper structural assembly, decarboxylation activity and drug interaction. This information provided the platform for the design of novel POA analogues with modifications at position 3 of the pyrazine ring. Analogue 2, which incorporates a bulky naphthamido group at this position, displayed a 1000-fold increase in enzyme inhibition, compared to POA, along with moderately improved antimycobacterial activity. The data demonstrate that an improved understanding of mechanistic and enzymatic features of key metabolic enzymes can stimulate design of more-potent PanD inhibitors.

Potency boost of a Mycobacterium tuberculosis dihydrofolate reductase inhibitor by multienzyme F420H2-dependent reduction

Bacterial metabolism can cause intrinsic drug resistance but can also convert inactive parent drugs into bioactive derivatives, as is the case for several antimycobacterial prodrugs. Here, we show that the intrabacterial metabolism of a Mtb dihydrofolate reductase (DHFR) inhibitor with moderate affinity for its target boosts its on-target activity by two orders of magnitude. This is a “prodrug-like” antimycobacterial that possesses baseline activity in the absence of intracellular bioactivation. By elucidating the metabolic enhancement mechanism, we have provided the basis for the rational optimization of a class of DHFR inhibitors and uncovered an antibacterial drug discovery concept.

Triaza-coumarin (TA-C) is a Mycobacterium tuberculosis (Mtb) dihydrofolate reductase (DHFR) inhibitor with an IC₅₀ (half maximal inhibitory concentration) of ∼1 µM against the enzyme. Despite this moderate target inhibition, TA-C shows exquisite antimycobacterial activity (MIC₅₀, concentration inhibiting growth by 50% = 10 to 20 nM). Here, we investigated the mechanism underlying this potency disconnect. To confirm that TA-C targets DHFR and investigate its unusual potency pattern, we focused on resistance mechanisms. In Mtb, resistance to DHFR inhibitors is frequently associated with mutations in thymidylate synthase thyA, which sensitizes Mtb to DHFR inhibition, rather than in DHFR itself. We observed thyA mutations, consistent with TA-C interfering with the folate pathway. A second resistance mechanism involved biosynthesis of the redox coenzyme F₄₂₀. Thus, we hypothesized that TA-C may be metabolized by Mtb F₄₂₀–dependent oxidoreductases (FDORs). By chemically blocking the putative site of FDOR-mediated reduction in TA-C, we reproduced the F₄₂₀-dependent resistance phenotype, suggesting that F₄₂₀H₂-dependent reduction is required for TA-C to exert its potent antibacterial activity. Indeed, chemically synthesized TA-C-Acid, the putative product of TA-C reduction, displayed a 100-fold lower IC₅₀ against DHFR. Screening seven recombinant Mtb FDORs revealed that at least two of these enzymes reduce TA-C. This redundancy in activation explains why no mutations in the activating enzymes were identified in the resistance screen. Analysis of the reaction products confirmed that FDORs reduce TA-C at the predicted site, yielding TA-C-Acid. This work demonstrates that intrabacterial metabolism converts TA-C, a moderately active “prodrug,” into a 100-fold-more-potent DHFR inhibitor, thus explaining the disconnect between enzymatic and whole-cell activity.

Pyrazinamide Susceptibility Is Driven by Activation of the SigE-Dependent Cell Envelope Stress Response in Mycobacterium tuberculosis

Pyrazinamide (PZA) plays a crucial role in first-line tuberculosis drug therapy. Unlike other antimicrobial agents, PZA is active against Mycobacterium tuberculosis only at low pH. The basis for this conditional drug susceptibility remains undefined. In this study, we utilized a genome-wide approach to interrogate potentiation of PZA action. We found that mutations in numerous genes involved in central metabolism as well as cell envelope maintenance and stress response are associated with PZA resistance. Further, we demonstrate that constitutive activation of the cell envelope stress response can drive PZA susceptibility independent of environmental pH. Consequently, exposure to peptidoglycan synthesis inhibitors, such as beta-lactams and d-cycloserine, potentiate PZA action through triggering this response. These findings illuminate a regulatory mechanism for conditional PZA susceptibility and reveal new avenues for enhancing potency of this important drug through targeting activation of the cell envelope stress response. IMPORTANCE For decades, pyrazinamide has served as a cornerstone of tuberculosis therapy. Unlike any other antitubercular drug, pyrazinamide requires an acidic environment to exert its action. Despite its importance, the driver of this conditional susceptibility has remained unknown. In this study, a genome-wide approach revealed that pyrazinamide action is governed by the cell envelope stress response. This observation was validated by orthologous approaches that demonstrate that a central player of this response, SigE, is both necessary and sufficient for potentiation of pyrazinamide action. Moreover, constitutive activation of this response through deletion of the anti-sigma factor gene rseA or exposure of bacilli to drugs that target the cell wall was found to potently drive pyrazinamide susceptibility independent of environmental pH. These findings force a paradigm shift in our understanding of pyrazinamide action and open new avenues for improving diagnostic and therapeutic tools for tuberculosis.

Bioenergetic Inhibitors: Antibiotic Efficacy and Mechanisms of Action in Mycobacterium tuberculosis

Development of novel anti-tuberculosis combination regimens that increase efficacy and reduce treatment timelines will improve patient compliance, limit side-effects, reduce costs, and enhance cure rates. Such advancements would significantly improve the global TB burden and reduce drug resistance acquisition. Bioenergetics has received considerable attention in recent years as a fertile area for anti-tuberculosis drug discovery. Targeting the electron transport chain (ETC) and oxidative phosphorylation machinery promises not only to kill growing cells but also metabolically dormant bacilli that are inherently more drug tolerant. Over the last two decades, a broad array of drugs targeting various ETC components have been developed. Here, we provide a focused review of the current state of art of bioenergetic inhibitors of Mtb with an in-depth analysis of the metabolic and bioenergetic disruptions caused by specific target inhibition as well as their synergistic and antagonistic interactions with other drugs. This foundation is then used to explore the reigning theories on the mechanisms of antibiotic-induced cell death and we discuss how bioenergetic inhibitors in particular fail to be adequately described by these models. These discussions lead us to develop a clear roadmap for new lines of investigation to better understand the mechanisms of action of these drugs with complex mechanisms as well as how to leverage that knowledge for the development of novel, rationally-designed combination therapies to cure TB.

Vitamin in the Crosshairs: Targeting Pantothenate and Coenzyme A Biosynthesis for New Antituberculosis Agents

Despite decades of dedicated research, there remains a dire need for new drugs against tuberculosis (TB). Current therapies are generations old and problematic. Resistance to these existing therapies results in an ever-increasing burden of patients with disease that is difficult or impossible to treat. Novel chemical entities with new mechanisms of action are therefore earnestly required. The biosynthesis of coenzyme A (CoA) has long been known to be essential in Mycobacterium tuberculosis (Mtb), the causative agent of TB. The pathway has been genetically validated by seminal studies in vitro and in vivo. In Mtb, the CoA biosynthetic pathway is comprised of nine enzymes: four to synthesize pantothenate (Pan) from l-aspartate and α-ketoisovalerate; five to synthesize CoA from Pan and pantetheine (PantSH). This review gathers literature reports on the structure/mechanism, inhibitors, and vulnerability of each enzyme in the CoA pathway. In addition to traditional inhibition of a single enzyme, the CoA pathway offers an antimetabolite strategy as a promising alternative. In this review, we provide our assessment of what appear to be the best targets, and, thus, which CoA pathway enzymes present the best opportunities for antitubercular drug discovery moving forward.

Exploiting Homoplasy in Genome-Wide Association Studies to Enhance Identification of Antibiotic-Resistance Mutations in Bacterial Genomes

Many antibacterial drugs have multiple mechanisms of resistance, which are often represented simultaneously by a mixture of resistance mutations (some more frequent than others) in a clinical population. This presents a challenge for Genome-Wide Association Studies (GWAS) methods, making it difficult to detect less prevalent resistance mechanisms purely through (weak) statistical associations. Homoplasy, or the occurrence of multiple independent mutations at the same site, is often observed with drug resistance mutations and can be a strong indicator of positive selection. However, traditional GWAS methods, such as those based on allele counting or linear regression, are not designed to take homoplasy into account. In this article, we present a new method, called ECAT (for Evolutionary Cluster-based Association Test), that extends traditional regression-based GWAS methods with the ability to take advantage of homoplasy. This is achieved through a preprocessing step which identifies hypervariable regions in the genome exhibiting statistically significant clusters of distinct evolutionary changes, to which association testing by a linear mixed model (LMM) is applied using GEMMA (a well-established LMM-based GWAS tool). Thus, the approach can be viewed as extending GEMMA from the usual site- or gene-level analysis to focusing on clustered regions of mutations. This approach was evaluated on a large collection of more than 600 clinical isolates of multidrug-resistant (MDR) Mycobacterium tuberculosis from Lima, Peru. We show that ECAT does a better job of detecting known resistance mutations for several antitubercular drugs (including less prevalent mutations with weaker associations), compared with (site- or gene-based) GEMMA, as representative of existing GWAS methods. The power of the multiphase approach in ECAT comes from focusing association testing on the hypervariable regions of the genome, which reduces complexity in the model and increases statistical power.

A biochemically-interpretable machine learning classifier for microbial GWAS

Current machine learning classifiers have successfully been applied to whole-genome sequencing data to identify genetic determinants of antimicrobial resistance (AMR), but they lack causal interpretation. Here we present a metabolic model-based machine learning classifier, named Metabolic Allele Classifier (MAC), that uses flux balance analysis to estimate the biochemical effects of alleles. We apply the MAC to a dataset of 1595 drug-tested Mycobacterium tuberculosis strains and show that MACs predict AMR phenotypes with accuracy on par with mechanism-agnostic machine learning models (isoniazid AUC = 0.93) while enabling a biochemical interpretation of the genotype-phenotype map. Interpretation of MACs for three antibiotics (pyrazinamide, para-aminosalicylic acid, and isoniazid) recapitulates known AMR mechanisms and suggest a biochemical basis for how the identified alleles cause AMR. Extending flux balance analysis to identify accurate sequence classifiers thus contributes mechanistic insights to GWAS, a field thus far dominated by mechanism-agnostic results.

Current machine learning classifiers have been applied to whole-genome sequencing data to identify determinants of antimicrobial resistance, but they lack interpretability. Here the authors present a metabolic machine learning classifier that uses flux balance analysis to estimate the biochemical effects of alleles.

Mycobacterium smegmatis msmeg_3314 is involved in pyrazinamide and fluoroquinolones susceptibility via NAD+/NADH dysregulation

Aim: To identify and characterize new mycobacterium pyrazinamide (PZA) resistance genes in addition to pncA, rpsA and panD. Materials & methods: To screen a Tn7 M. smegmatis mc²155 transposon library using 50 μM PZA and a PZA hypersensitive mutant (M492) was obtained. MIC was further used to confirm the hypersensitivity of M492 mutant by culturing the mutant in Middlebrook 7H9 liquid medium at 37°C. Results: msmeg_3314 is the gene underlying the hypersensitive phenotype of mutant M492. The observed resistance to PZA and fluoroquinolones involved the alteration of Mycobacterium cell wall permeability and the dissipation of the proton motive force. NAD⁺/NADH dysregulation and attenuated glyoxylate shunt might underlie the declined scavenging capacity of reactive oxygen species in the msmeg_3314-deficient mutants. Conclusion: msmeg_ 3314 is a novel gene involved in pyrazinamide resistance and might be a new candidate for drugs target.

Pyrazinamide triggers degradation of its target aspartate decarboxylase

Pyrazinamide is a sterilizing first-line tuberculosis drug. Genetic, metabolomic and biophysical analyses previously demonstrated that pyrazinoic acid, the bioactive form of the prodrug pyrazinamide (PZA), interrupts biosynthesis of coenzyme A in Mycobacterium tuberculosis by binding to aspartate decarboxylase PanD. While most drugs act by inhibiting protein function upon target binding, we find here that pyrazinoic acid is only a weak enzyme inhibitor. We show that binding of pyrazinoic acid to PanD triggers degradation of the protein by the caseinolytic protease ClpC1-ClpP. Thus, the old tuberculosis drug pyrazinamide exerts antibacterial activity by acting as a target degrader, a mechanism of action that has recently emerged as a successful strategy in drug discovery across disease indications. Our findings provide the basis for the rational discovery of next generation PZA.

It has been shown that the bioactive component of pyrazinamide, pyrazinoic acid (POA), blocks coenzyme A biosynthesis in M. tuberculosis by binding to the aspartate decarboxylase PanD. Here the authors show that pyrazinamide triggers degradation of PanD by stimulating its degradation by the caseinolytic protease Clp.

The Bewildering Antitubercular Action of Pyrazinamide

Pyrazinamide (PZA) is a cornerstone antimicrobial drug used exclusively for the treatment of tuberculosis (TB). Due to its ability to shorten drug therapy by 3 months and reduce disease relapse rates, PZA is considered an irreplaceable component of standard first-line short-course therapy for drug-susceptible TB and second-line treatment regimens for multidrug-resistant TB. Despite over 60 years of research on PZA and its crucial role in current and future TB treatment regimens, the mode of action of this unique drug remains unclear. Defining the mode of action for PZA will open new avenues for rational design of novel therapeutic approaches for the treatment of TB. In this review, we discuss the four prevailing models for PZA action, recent developments in modulation of PZA susceptibility and resistance, and outlooks for future research and drug development.

Detection of Novel Gene Mutations Associated with Pyrazinamide Resistance in Multidrug-Resistant Mycobacterium tuberculosis Clinical Isolates in Southern China

Objective

Pyrazinamide (PZA) is a cornerstone of modern tuberculosis regimens. This study aimed to investigate the performance of genotypic testing of pncA ⁺ upstream region, rpsA, panD, Rv2783c, and clpC1 genes to add insights for more accurate molecular diagnosis of PZA-resistant (R) Mycobacterium tuberculosis.

Methods

Drug susceptibility testing, sequencing analysis of PZA-related genes including the entire operon of pncA (Rv2044c-pncA-Rv2042c) and PZase assay were performed for 448 M. tuberculosis clinical isolates.

Results

Our data showed that among 448 M. tuberculosis clinical isolates, 113 were MDR, 195 pre-XDR and 70 XDR TB, while the remaining 70 strains had other combinations of drug-resistance. A total of 60.04% (269/448) M. tuberculosis clinical isolates were resistant to PZA, of which 78/113 were MDR, 119/195 pre-XDR and 29/70 XDR TB strains. PZA^R isolates have predominance (83.3%) of Beijing genotype. Genotypic characterization of Rv2044c-pncA-Rv2042c revealed novel nonsynonymous mutations in Rv2044c with negative PZase activity which led to confer PZA^R. Compared with phenotypic data, 84.38% (227/269) PZA^R strains with mutations in pncA ⁺ upstream region exhibited 83.64% sensitivity but the combined evaluation of the mutations in rpsA 2.60% (7/269), panD 1.48% (4/269), Rv2783c 1.11% (3/269) and Rv2044c 0.74% (2/269) increased the sensitivity to 89.59%. Fifty-seven novel mutations were identified in this study. Interestingly, a frameshift deletion (C-114del) in upstream of pncA^wt nullified the effect of A-11G mutation and induced positive PZase activity, divergent from five PZase negative A-11G PZA^R mutants. Twenty-six PZA^R strains having wild-type-sequenced genes with positive or negative PZase suggest the existence of unknown resistance mechanisms.

Conclusion

Our study revealed that PZA^R rate in MDR and pre-XDR TB was markedly higher in southern China. The concomitant evaluation of pncA ⁺ UFR, rpsA, panD, Rv2783c, and Rv2044c provides more dependable genotypic results of PZA resistance. Fifty-seven novel mutations/indels in this study may play a vital role as diagnostic markers. The upstream region of pncA and PZase regulation are valuable to explore the unknown mechanism of PZA-resistance.

High-Throughput Screen for Cell Wall Synthesis Network Module in Mycobacterium tuberculosis Based on Integrated Bioinformatics Strategy

BACKGROUND

Mycobacterium tuberculosis is one of the deadliest pathogens in humans. Co-infection of M. tuberculosis with HIV and the emergence of multi-drug-resistant tuberculosis (TB) constitute a serious global threat. However, no effective anti-TB drugs are available, with the exception of first-line drugs such as isoniazid. The cell wall of M. tuberculosis, which is primarily responsible for the lack of effective anti-TB drugs and the escape of the bacteria from host immunity, is an important drug target. The core components of the cell wall of M. tuberculosis are peptidoglycan, arabinogalactan, and mycotic acid. However, the functional genome and metabolic regulation pathways for the M. tuberculosis cell wall are still unknown. In this study, we used the biclustering algorithm integrated into cMonkey, sequence alignment, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and other bioinformatics methods to scan the whole genome of M. tuberculosis as well as to identify and statistically analyze the genes related to the synthesis of the M. tuberculosis cell wall.

METHOD

We performed high-throughput genome-wide screening for M. tuberculosis using Biocarta, KEGG, National Cancer Institute Pathway Interaction Database (NCI-PID), HumanCyc, and Reactome. We then used the Database of Origin and Registration (DOOR) established in our laboratory to classify the collection of operons for M. tuberculosis cell wall synthetic genes. We used the cMonkey double clustering algorithm to perform clustering analysis on the gene expression profile of M. tuberculosis for cell wall synthesis. Finally, we visualized the results using Cytoscape.

RESULT AND CONCLUSION

Through bioinformatics and statistical analyses, we identified 893 M. tuberculosis H37Rv cell wall synthesis genes, distributed in 20 pathways, involved in 46 different functions related to cell wall synthesis, and clustered in 386 modules. We identified important pivotal genes and proteins in the cell wall synthesis pathway such as murA, a class of operons containing genes involved in cell wall synthesis such as ID6951, and a class of operons indispensable for the survival of the bacteria. In addition, we found 41 co-regulatory modules for cell wall synthesis and five co-expression networks of molecular complexes involved in peptidoglycan biosynthesis, membrane transporter synthesis, and other cell wall processes.

Pharmacological and Molecular Mechanisms Behind the Sterilizing Activity of Pyrazinamide

Inclusion of pyrazinamide (PZA) in the tuberculosis (TB) drug regimen during the 1970s enabled a reduction in treatment duration from 12 to 6 months. PZA has this remarkable effect in patients despite displaying poor potency against Mycobacterium tuberculosis (Mtb) in vitro. The pharmacological basis for the in vivo sterilizing activity of the drug has remained obscure and its bacterial target controversial. Recently it was shown that PZA penetrates necrotic caseous TB lung lesions and kills nongrowing, drug-tolerant bacilli. Furthermore, it was uncovered that PZA inhibits bacterial Coenzyme A biosynthesis. It may block this pathway by triggering degradation of its target, aspartate decarboxylase. The elucidation of the pharmacological and molecular mechanisms of PZA provides the basis for the rational discovery of the next-generation PZA with improved in vitro potency while maintaining attractive pharmacological properties.

Galleria mellonella: An Infection Model for Screening Compounds Against the Mycobacterium tuberculosis Complex

Drug screening models have a vital role in the development of novel antimycobacterial agents which are urgently needed to tackle drug-resistant tuberculosis (TB). We recently established the larvae of the insect Galleria mellonella (greater wax moth) as a novel infection model for the Mycobacterium tuberculosis complex. Here we demonstrate its use as a rapid and reproducible screen to evaluate antimycobacterial drug efficacy using larvae infected with bioluminescent Mycobacterium bovis BCG lux. Treatment improved larval survival outcome and, with the exception of pyrazinamide, was associated with a significant reduction in in vivo mycobacterial bioluminescence over a 96 h period compared to the untreated controls. Isoniazid and rifampicin displayed the greatest in vivo efficacy and survival outcome. Thus G. mellonella, infected with bioluminescent mycobacteria, can rapidly determine in vivo drug efficacy, and has the potential to significantly reduce and/or replace the number of animals used in TB research.

Targeted protein degradation in antibacterial drug discovery?

Drug induced degradation of a target protein is a novel concept in drug discovery. Traditionally drugs modulate activity, as opposed to abundance, of their targets. Degradation inducing ligands act catalytically. Thus, one advantage of target degradation over the classical on-target mechanism is that lower drug concentration may be sufficient to cause the desired cellular effects. The first promoters of target degradation were discovered unintentionally: it turned out that some drugs 'accidently' promote degradation of their target by the cellular proteolytic machinery. Elegant methods were developed to target specific proteins of interest for degradation, thus enabling the rational discovery of degradation inducers. The application of targeted degradation has so far been limited to human cells. Recently, we discovered that an antibacterial drug, the anti-tuberculosis antibiotic pyrazinamide, functions as a promotor of degradation of its bacterial target. Increasing antimicrobial resistance makes the discovery of novel antibiotics more urgent than ever. Can rational target degradation be applied for the discovery of anti-bacterials? Here, we first discuss briefly some historic examples and then recent approaches in rational target degradation for human diseases. Then, we describe how the first anti-bacterial target degradation promoter pyrazinamide triggers removal of its target. Efforts are under way to exploit this specific mechanistic knowledge for the discovery of next generation pyrazinamide. We end with the big - and open - question whether targeted protein degradation as an approach to anti-bacterial drug discovery can be generalized, similar to what has been achieved in the area of drug discovery for human diseases.

Impact of the host environment on the antitubercular action of pyrazinamide

Pyrazinamide remains the only drug in the tuberculosis pharmacopeia to drastically shorten first-line therapy from nine to six months. Due to its unparalleled ability to sterilize non-replicating bacilli and reduce relapse rates, PZA is expected to be irreplaceable in future therapies against tuberculosis. While the molecular target of PZA is unclear, recent pharmacokinetic studies using small animal models and patient samples have highlighted the importance of host metabolism and immune responses in PZA efficacy. Delineating which host factors are important for PZA action will be integral to the design of next-generation therapies to shorten current TB drug regimens as well as to overcome treatment limitations in some patients. In this review, we discuss evidence for influence of the host environment on PZA activity, targets for PZA mechanism of action, recent studies in PZA pharmacokinetics, PZA antagonism and synergy with other first-line anti-TB drugs, and implications for future research.

'Happy the man, who, studying nature's laws, Thro' known effects can trace the secret cause.' Do we have enough pieces to solve the pyrazinamide puzzle?

A low pH was assumed to be required for the activity of pyrazinoic acid (the active form of pyrazinamide) against Mycobacterium tuberculosis, but recently activity has been demonstrated at neutral pH. Renewed interest in pyrazinamide has led to an increasing number of potential targets and the suspicion that pyrazinamide is a 'dirty drug'. However, it is our opinion that the recent demonstration that pyrazinoic acid is active against PanD provides an alternative explanation for the secret of pyrazinamide's unusual activity. In this article we propose that PanD is the primary target of pyrazinoic acid but expression of pyrazinoic acid susceptibility requires an intact stress response. As the mycobacterial stress response requires the interaction of a number of genes, disruption of any could result in an inability to enter the susceptible phenotype. We believe this model can explain most of the recent observations of the seemingly diverse spectrum of activity of pyrazinamide.

Phase variation in Mycobacterium tuberculosis glpK produces transiently heritable drug tolerance

The length and complexity of tuberculosis (TB) therapy, as well as the propensity of Mycobacterium tuberculosis to develop drug resistance, are major barriers to global TB control efforts. M. tuberculosis is known to have the ability to enter into a drug-tolerant state, which may explain many of these impediments to TB treatment. We have identified a mechanism of genetically encoded but rapidly reversible drug tolerance in M. tuberculosis caused by transient frameshift mutations in a homopolymeric tract (HT) of 7 cytosines (7C) in the glpK gene. Inactivating frameshift mutations associated with the 7C HT in glpK produce small colonies that exhibit heritable multidrug increases in minimal inhibitory concentrations and decreases in drug-dependent killing; however, reversion back to a fully drug-susceptible large-colony phenotype occurs rapidly through the introduction of additional insertions or deletions in the same glpK HT region. These reversible frameshift mutations in the 7C HT of M. tuberculosis glpK occur in clinical isolates, accumulate in M. tuberculosis-infected mice with further accumulation during drug treatment, and exhibit a reversible transcriptional profile including induction of dosR and sigH and repression of kstR regulons, similar to that observed in other in vitro models of M. tuberculosis tolerance. These results suggest that GlpK phase variation may contribute to drug tolerance, treatment failure, and relapse in human TB. Drugs effective against phase-variant M. tuberculosis may hasten TB treatment and improve cure rates.

Comparative fitness analysis of D-cycloserine resistant mutants reveals both fitness-neutral and high-fitness cost genotypes

Drug resistant infections represent one of the most challenging medical problems of our time. D-cycloserine is an antibiotic used for six decades without significant appearance and dissemination of antibiotic resistant strains, making it an ideal model compound to understand what drives resistance evasion. We therefore investigated why Mycobacterium tuberculosis fails to become resistant to D-cycloserine. To address this question, we employed a combination of bacterial genetics, genomics, biochemistry and fitness analysis in vitro, in macrophages and in mice. Altogether, our results suggest that the ultra-low rate of emergence of D-cycloserine resistance mutations is the dominant biological factor delaying the appearance of clinical resistance to this antibiotic. Furthermore, we also identified potential compensatory mechanisms able to minimize the severe fitness costs of primary D-cycloserine resistance conferring mutations.

D-cycloserine (DCS) has been used for decades to treat Mycobacterium tuberculosis (Mtb) but resistance is rarely observed in clinical isolates. Here, the authors report ultra-low rate of emergence of resistance mutations as the underlying mechanism of DCS resistance evasion in Mtb.

Methyl arachidonyl fluorophosphonate inhibits Mycobacterium tuberculosis thioesterase TesA and globally affects vancomycin susceptibility

Phthiocerol dimycocerosates and phenolic glycolipids (PGL) are considered as major virulence elements of Mycobacterium tuberculosis, in particular because of their involvement in cell wall impermeability and drug resistance. The biosynthesis of these waxy lipids involves multiple enzymes, including thioesterase A (TesA). We observed that purified recombinant M. tuberculosis TesA is able to dimerize in the presence of palmitoyl-CoA and our 3D structure model of TesA with this acyl-CoA suggests hydrophobic interaction requirement for dimerization. Furthermore, we identified that methyl arachidonyl fluorophosphonate, which inhibits TesA by covalently modifying the catalytic serine, also displays a synergistic antimicrobial activity with vancomycin further warranting the development of TesA inhibitors as valuable antituberculous drug candidates.

Common Variants in the Glycerol Kinase Gene Reduce Tuberculosis Drug Efficacy

Despite the administration of multiple drugs that are highly effective in vitro, tuberculosis (TB) treatment requires prolonged drug administration and is confounded by the emergence of drug-resistant strains. To understand the mechanisms that limit antibiotic efficacy, we performed a comprehensive genetic study to identify Mycobacterium tuberculosis genes that alter the rate of bacterial clearance in drug-treated mice. Several functionally distinct bacterial genes were found to alter bacterial clearance, and prominent among these was the glpK gene that encodes the glycerol-3-kinase enzyme that is necessary for glycerol catabolism. Growth on glycerol generally increased the sensitivity of M. tuberculosis to antibiotics in vitro, and glpK-deficient bacteria persisted during antibiotic treatment in vivo, particularly during exposure to pyrazinamide-containing regimens. Frameshift mutations in a hypervariable homopolymeric region of the glpK gene were found to be a specific marker of multidrug resistance in clinical M. tuberculosis isolates, and these loss-of-function alleles were also enriched in extensively drug-resistant clones. These data indicate that frequently observed variation in the glpK coding sequence produces a drug-tolerant phenotype that can reduce antibiotic efficacy and may contribute to the evolution of resistance.IMPORTANCE TB control is limited in part by the length of antibiotic treatment needed to prevent recurrent disease. To probe mechanisms underlying survival under antibiotic pressure, we performed a genetic screen for M. tuberculosis mutants with altered susceptibility to treatment using the mouse model of TB. We identified multiple genes involved in a range of functions which alter sensitivity to antibiotics. In particular, we found glycerol catabolism mutants were less susceptible to treatment and that common variation in a homopolymeric region in the glpK gene was associated with drug resistance in clinical isolates. These studies indicate that reversible high-frequency variation in carbon metabolic pathways can produce phenotypically drug-tolerant clones and have a role in the development of resistance.

Acid Fasting: Modulation of Mycobacterium tuberculosis Metabolism at Acidic pH

Mycobacterium tuberculosis (Mtb) senses and adapts to acidic host environments during the course of pathogenesis. Mutants defective in acidic pH-dependent adaptations are often attenuated during macrophage or animal infections, supporting that these pathways are essential for pathogenesis and represent important new targets for drug discovery. This review examines a confluence of findings supporting that Mtb has restricted metabolism at acidic pH that results in the slowing of bacterial growth and changes in redox homeostasis. It is proposed that induction of the PhoPR regulon and anaplerotic metabolism, in concert with the restricted use of specific carbon sources, functions to counter reductive stress associated with acidic pH.

Gut Microbiota Metabolite Indole Propionic Acid Targets Tryptophan Biosynthesis in Mycobacterium tuberculosis

New drugs against tuberculosis are urgently needed. The tryptophan (Trp) analog indole propionic acid (IPA) is the first antitubercular metabolite produced by human gut bacteria. Here, we show that this antibiotic blocks Trp synthesis, an in vivo essential biosynthetic pathway in M. tuberculosis. Intriguingly, IPA acts by decoupling a bacterial feedback regulatory mechanism: it mimics Trp as allosteric inhibitor of anthranilate synthase, thereby switching off Trp synthesis regardless of intracellular Trp levels. The identification of IPA’s target paves the way for the discovery of more potent TrpE ligands employing rational, target-based lead optimization.

Indole propionic acid (IPA), produced by the gut microbiota, is active against Mycobacterium tuberculosisin vitro and in vivo. However, its mechanism of action is unknown. IPA is the deamination product of tryptophan (Trp) and thus a close structural analog of this essential aromatic amino acid. De novo Trp biosynthesis in M. tuberculosis is regulated through feedback inhibition: Trp acts as an allosteric inhibitor of anthranilate synthase TrpE, which catalyzes the first committed step in the Trp biosynthesis pathway. Hence, we hypothesized that IPA may mimic Trp as an allosteric inhibitor of TrpE and exert its antimicrobial effect by blocking synthesis of Trp at the TrpE catalytic step. To test our hypothesis, we carried out metabolic, chemical rescue, genetic, and biochemical analyses. Treatment of mycobacteria with IPA inhibited growth and reduced the intracellular level of Trp, an effect abrogated upon supplementation of Trp in the medium. Missense mutations at the allosteric Trp binding site of TrpE eliminated Trp inhibition and caused IPA resistance. In conclusion, we have shown that IPA blocks Trp biosynthesis in M. tuberculosis via inhibition of TrpE by mimicking the physiological allosteric inhibitor of this enzyme.

Machine learning and structural analysis of Mycobacterium tuberculosis pan-genome identifies genetic signatures of antibiotic resistance

Mycobacterium tuberculosis is a serious human pathogen threat exhibiting complex evolution of antimicrobial resistance (AMR). Accordingly, the many publicly available datasets describing its AMR characteristics demand disparate data-type analyses. Here, we develop a reference strain-agnostic computational platform that uses machine learning approaches, complemented by both genetic interaction analysis and 3D structural mutation-mapping, to identify signatures of AMR evolution to 13 antibiotics. This platform is applied to 1595 sequenced strains to yield four key results. First, a pan-genome analysis shows that M. tuberculosis is highly conserved with sequenced variation concentrated in PE/PPE/PGRS genes. Second, the platform corroborates 33 genes known to confer resistance and identifies 24 new genetic signatures of AMR. Third, 97 epistatic interactions across 10 resistance classes are revealed. Fourth, detailed structural analysis of these genes yields mechanistic bases for their selection. The platform can be used to study other human pathogens.

Comparing and Validating Machine Learning Models for Mycobacterium tuberculosis Drug Discovery

Tuberculosis is a global health dilemma. In 2016, the WHO reported 10.4 million incidences and 1.7 million deaths. The need to develop new treatments for those infected with Mycobacterium tuberculosis ( Mtb) has led to many large-scale phenotypic screens and many thousands of new active compounds identified in vitro. However, with limited funding, efforts to discover new active molecules against Mtb needs to be more efficient. Several computational machine learning approaches have been shown to have good enrichment and hit rates. We have curated small molecule Mtb data and developed new models with a total of 18,886 molecules with activity cutoffs of 10 μM, 1 μM, and 100 nM. These data sets were used to evaluate different machine learning methods (including deep learning) and metrics and to generate predictions for additional molecules published in 2017. One Mtb model, a combined in vitro and in vivo data Bayesian model at a 100 nM activity yielded the following metrics for 5-fold cross validation: accuracy = 0.88, precision = 0.22, recall = 0.91, specificity = 0.88, kappa = 0.31, and MCC = 0.41. We have also curated an evaluation set ( n = 153 compounds) published in 2017, and when used to test our model, it showed the comparable statistics (accuracy = 0.83, precision = 0.27, recall = 1.00, specificity = 0.81, kappa = 0.36, and MCC = 0.47). We have also compared these models with additional machine learning algorithms showing Bayesian machine learning models constructed with literature Mtb data generated by different laboratories generally were equivalent to or outperformed deep neural networks with external test sets. Finally, we have also compared our training and test sets to show they were suitably diverse and different in order to represent useful evaluation sets. Such Mtb machine learning models could help prioritize compounds for testing in vitro and in vivo.

Genetics and roadblocks of drug resistant tuberculosis

Considering the extensive evolutionary history of Mycobacterium tuberculosis, anti-Tuberculosis (TB) drug therapy exerts a recent selective pressure. However, in a microorganism devoid of horizontal gene transfer and with a strictly clonal populational structure such as M. tuberculosis the usual, but not sole, path to overcome drug susceptibility is through de novo mutations on a relatively strict set of genes. The possible allelic diversity that can be associated with drug resistance through several mechanisms such as target alteration or target overexpression, will dictate how these genes can become associated with drug resistance. The success demonstrated by this pathogenic microbe in this latter process and its ability to spread is currently one of the major obstacles to an effective TB elimination. This article reviews the action mechanism of the more important anti-TB drugs, including bedaquiline and delamanid, along with new findings on specific resistance mechanisms. With the development, validation and endorsement of new in vitro molecular tests for drug resistance, knowledge on these resistance mechanisms and microevolutionary dynamics leading to the emergence and fixation of drug resistance mutations within the host is highly important. Additionally, the fitness toll imposed by resistance development is also herein discussed together with known compensatory mechanisms. By elucidating the possible mechanisms that enable one strain to reacquire the original fitness levels, it will be theoretically possible to make more informed decisions and develop novel strategies that can force M. tuberculosis microevolutionary trajectory down through a path of decreasing fitness levels.

Impact of immunopathology on the antituberculous activity of pyrazinamide

In the 1970s, inclusion of pyrazinamide (PZA) in the drug regimen of tuberculosis (TB) patients for the first 2 mo achieved a drastic reduction of therapy duration. Until now, however, the mechanisms underlying PZA's unique contribution to efficacy have remained controversial, and animal efficacy data vary across species. To understand how PZA kills bacterial populations present in critical lung lesion compartments, we first characterized a rabbit model of active TB, showing striking similarities in lesion types and fates to nonhuman primate models deemed the most appropriate surrogates of human TB. We next employed this model with lesion-centric molecular and bacteriology readouts to demonstrate that PZA exhibits potent activity against Mycobacterium tuberculosis residing in difficult-to-sterilize necrotic lesions. Our data also indicate that PZA is slow acting, suggesting that PZA administration beyond the first 2 mo may accelerate the cure. In conclusion, we provide a pharmacodynamic explanation for PZA's treatment-shortening effect and deliver new tools to dissect the contribution of immune response versus drug at the lesion level.