Mutations in the pncA and rpsA genes among 77 Mycobacterium tuberculosis isolates in Kazakhstan

Genetic diversities and drug resistance in Mycobacterium bovis isolates from zoonotic tuberculosis using whole genome sequencing

BACKGROUND

Zoonotic human tuberculosis (TB) caused by Mycobacterium bovis (M. bovis) is as vital as Mycobacterium tuberculosis, however with scarce available information. We aimed to use whole-genome sequencing (WGS) technology to take a deep insight into the circulating genotypes of human M. bovis and the genomic characteristics underlying virulence and drug resistance.

METHODS

The study included smear positive Ziehl-Neelsen samples from patients with suspected tuberculosis. Samples were cultured on Lowenstein-Jensen media and suspected colonies of M. bovis were selected to undergo DNA extraction and WGS. Data was analysed using the Bacterial and Viral Bioinformatics Resource Center (BV-BRC), and online bioinformatics tools. A phylogenetic tree was constructed for our sequenced strains, in addition to a set of 59 previously sequenced M. bovis genomes from different hosts and countries.

RESULTS

Out of total 112 mycobacterial positive cultures, five M. bovis were isolated and underwent WGS. All sequenced strains belonged to Mycobacterium tuberculosis var bovis, spoligotype BOV_1; BOV_11. Resistance gene mutations were determined in 100% of strains to pyrazinamide (pncA and rpsA), isoniazid (KatG and ahpC), ethambutol (embB, embC, embR and ubiA), streptomycin (rpsl) and fluoroquinolones (gyrA and gyrB). Rifampin (rpoB and rpoC) and delamanid (fbiC) resistance genes were found in 80% of strains. The major represented virulence classes were the secretion system, cell surface components and regulation system. The phylogenetic analysis revealed close genetic relatedness of three sequenced M. bovis strains to previous reported cow strains from Egypt and human strains from France, as well as relatedness of one M. bovis strain to four human Algerian strains. One sequenced strain was related to one cow strain from Egypt and a human strain from South Africa.

CONCLUSIONS

All sequenced M. bovis isolates showed the same spoligotype, but diverse phylogeny. Resistance gene mutations were detected for anti-TB drugs including pyrazinamide, isoniazid, streptomycin, ethambutol, fluoroquinolones, cycloserine, rifampin and delamanid. The virulence profile comprised genes assigned mainly to secretion system, cell surface components and regulation system. Phylogenetic analysis revealed genetic relatedness between our isolates and previously sequenced bovine strains from Egypt as well as human strains from other nearby countries in the region.

A robust computational quest: Discovering potential hits to improve the treatment of pyrazinamide-resistant Mycobacterium tuberculosis

The rise of pyrazinamide (PZA)-resistant strains of Mycobacterium tuberculosis (MTB) poses a major challenge to conventional tuberculosis (TB) treatments. PZA, a cornerstone of TB therapy, must be activated by the mycobacterial enzyme pyrazinamidase (PZase) to convert its active form, pyrazinoic acid, which targets the ribosomal protein S1. Resistance, often associated with mutations in the RpsA protein, complicates treatment and highlights a critical gap in the understanding of structural dynamics and mechanisms of resistance, particularly in the context of the G97D mutation. This study utilizes a novel integration of computational techniques, including multiscale biomolecular and molecular dynamics simulations, physicochemical and medicinal chemistry predictions, quantum computations and virtual screening from the ZINC and Chembridge databases, to elucidate the resistance mechanism and identify lead compounds that have the potential to improve treatment outcomes for PZA-resistant MTB, namely ZINC15913786, ZINC20735155, Chem10269711, Chem10279789 and Chem10295790. These computational methods offer a cost-effective, rapid alternative to traditional drug trials by bypassing the need for organic subjects while providing highly accurate insight into the binding sites and efficacy of new drug candidates. The need for rapid and appropriate drug development emphasizes the need for robust computational analysis to justify further validation through in vitro and in vivo experiments.

Prevalence of Beijing Central Asian/Russian Cluster 94-32 among Multidrug-Resistant M. tuberculosis in Kazakhstan

The Beijing genotype is the most distributed M. tuberculosis family in Kazakhstan. In this study, we identified dominant Beijing clusters in Kazakhstan and assessed their drug susceptibility profiles and association with the most widely spread mutation Ser531Leu of the rpoB gene and the mutation Ser315Thr of the katG gene associated with resistance to rifampicin and isoniazid, respectively. M. tuberculosis isolates (n = 540) from new TB cases were included in the study. MIRU-VNTR genotyping was performed for 540 clinical isolates to determine M. tuberculosis families using 24 loci. RD analysis was additionally performed for the Beijing isolates. The identification of mutations in the drug-resistance genes of M. tuberculosis was performed with allele-specific real-time PCR and Sanger sequencing. The Beijing genotype was identified in 60% (324/540) of the clinical isolates. Central Asian/Russian cluster 94-32 was the most distributed cluster among the Beijing isolates (50.3%; 163/324). Three other dominant Beijing clusters were identified as 94-33 (3.4%; 11/324), 100-32 (3.1%; 10/324) and 99-32 (3.1%; 10/324). The Beijing genotype was associated with drug-resistant TB (p < 0.0001), including multidrug-resistant TB (p < 0.0001), in our study. An association of the mutation Ser531Leu of the rpoB gene with the Beijing genotype was found (p < 0.0001; OR = 16.0000; 95%CI: 4.9161-52.0740). Among the Beijing isolates, cluster 94-32 showed an association with MDR-TB (p = 0.021). This is why the evaluation of the Beijing genotype and its clusters is needed to control MDR-TB in Kazakhstan.

Association between fatty acid metabolism gene mutations and Mycobacterium tuberculosis transmission revealed by whole genome sequencing

BACKGROUND

Fatty acid metabolism greatly promotes the virulence and pathogenicity of Mycobacterium tuberculosis (M.tb). However, the regulatory mechanism of fatty acid metabolism in M.tb remains to be elucidated, and limited evidence about the effects of gene mutations in fatty acid metabolism on the transmission of M.tb was reported.

RESULTS

Overall, a total of 3193 M.tb isolates were included in the study, of which 1596 (50%) were genomic clustered isolates. Most of the tuberculosis isolates belonged to lineage2(n = 2744,85.93%), followed by lineage4(n = 439,13.75%) and lineage3(n = 10,0.31%).Regression results showed that the mutations of gca (136,605, 317G > C, Arg106Pro; OR, 22.144; 95% CI, 2.591-189.272), ogt(1,477,346, 286G > C ,Gly96Arg; OR, 3.893; 95%CI, 1.432-10.583), and rpsA (1,834,776, 1235 C > T, Ala412Val; OR, 3.674; 95% CI, 1.217-11.091) were significantly associated with clustering; mutations in gca and rpsA were also significantly associated with clustering of lineage2. Mutation in arsA(3,001,498, 885 C > G, Thr295Thr; OR, 6.278; 95% CI, 2.508-15.711) was significantly associated with cross-regional clusters. We also found that 20 mutation sites were positively correlated with cluster size, while 11 fatty acid mutation sites were negatively correlated with cluster size.

CONCLUSION

Our research results suggested that mutations in genes related to fatty acid metabolism were related to the transmission of M.tb. This research could help in the future control of the transmission of M.tb.

Free energy landscape and thermodynamics properties of novel mutations in PncA of pyrazinamide resistance isolates of Mycobacterium tuberculosis

Pyrazinamide (PZA) is one of the first-line antituberculosis therapy, active against non-replicating Mycobacterium tuberculosis (Mtb). The conversion of PZA into pyrazinoic acid (POA), the active form, required the activity of pncA gene product pyrazinamidase (PZase) activity. Mutations occurred in pncA are the primary cause behind the PZA resistance. However, the resistance mechanism is important to explore using high throughput computational approaches. Here we aimed to explore the mechanism of PZA resistance behind novel P62T, L120R, and V130M mutations in PZase using 200 ns molecular dynamics (MD) simulations. MD simulations were performed to observe the structural changes for these three mutants (MTs) compared to the wild types (WT). Root means square fluctuation, the radius of gyration, free energy landscape, root means square deviation, dynamic cross-correlation motion, and pocket volume were found in variation between WT and MTs, revealing the effects of P62T, L120R, and V130M. The free energy conformational landscape of MTs differs significantly from the WT system, lowering the binding of PZA. The geometric shape complementarity of the drug (PZA) and target protein (PZase) further confirmed that P62T, L120R, and V130M affect the protein structure. These effects on PZase may cause vulnerability to convert PZA into POA.Communicated by Ramaswamy H. Sarma.

Whole-Genome Sequence-Based Characterization of Pre-XDR M. tuberculosis Clinical Isolates Collected in Kazakhstan

BACKGROUND

Kazakhstan has a high burden of multidrug-resistant tuberculosis in the Central Asian region. This study aimed to perform genomic characterization of Mycobacterium tuberculosis strains obtained from Kazakhstani patients with pre-extensively drug-resistant tuberculosis diagnosed in Kazakhstan.

METHODS

Whole-genome sequencing was performed on 10 pre-extensively drug-resistant M. tuberculosis strains from different regions of Kazakhstan. All strains had high-confidence resistance mutations according to the resistance grading system previously established by the World Health Organization. The genome analysis was performed using TB-Profiler, Mykrobe, CASTB, and ResFinder.

RESULTS

Valuable information for understanding the genetic diversity of tuberculosis in Kazakhstan can also be obtained from whole-genome sequencing. The results from the Phenotypic Drug Susceptibility Testing (DST) of bacterial strains were found to be consistent with the drug resistance information obtained from genomic data that characterized all isolates as pre-XDR. This information can help in developing targeted prevention and control strategies based on the local epidemiology of tuberculosis. Furthermore, the data obtained from whole-genome sequencing can help in tracing the transmission pathways of tuberculosis and facilitating early detection of outbreaks.

CONCLUSIONS

The results from whole-genome sequencing of tuberculosis clinical samples in Kazakhstan provide important insights into the drug resistance patterns and genetic diversity of tuberculosis in the country. These results can contribute to the improvement of tuberculosis control and management programs in Kazakhstan.

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.

Genomic Analysis of Multidrug-Resistant Mycobacterium tuberculosis Strains From Patients in Kazakhstan

Tuberculosis (TB) is an infectious disease that remains an essential public health problem in many countries. Despite decreasing numbers of new cases worldwide, the incidence of antibiotic-resistant forms (multidrug resistant and extensively drug-resistant) of TB is increasing. Next-generation sequencing technologies provide a high-throughput approach to identify known and novel potential genetic variants that are associated with drug resistance in Mycobacterium tuberculosis (Mtb). There are limited reports and data related to whole-genome characteristics of drug-resistant Mtb strains circulating in Kazakhstan. Here, we report whole-genome sequencing and analysis results of eight multidrug-resistant strains collected from TB patients in Kazakhstan. Genotyping and validation of all strains by MIRU-VNTR and spoligotyping methodologies revealed that these strains belong to the Beijing family. The spectrum of specific and potentially novel genomic variants (single-nucleotide polymorphisms, insertions, and deletions) related to drug resistance was identified and annotated. ResFinder, CARD, and CASTB antibiotic resistance databases were used for the characterization of genetic variants in genes associated with drug resistance. Our results provide reference data and genomic profiles of multidrug-resistant isolates for further comparative studies and investigations of genetic patterns in drug-resistant Mtb strains.

Comparative Performance of Genomic Methods for the Detection of Pyrazinamide Resistance and Heteroresistance in Mycobacterium tuberculosis

Pyrazinamide is an important component of both drug-susceptible and drug-resistant tuberculosis treatment regimens. Although approximately 50% of rifampin-resistant isolates are also resistant to pyrazinamide, pyrazinamide susceptibility testing is not routinely performed due to the challenging nature of the assay. We investigated the diagnostic accuracy of genotypic and phenotypic methods and explored the occurrence of pyrazinamide heteroresistance. We assessed pyrazinamide susceptibility among 358 individuals enrolled in the South African EXIT-RIF cohort using Sanger and targeted deep sequencing (TDS) of the pncA gene, whole-genome sequencing (WGS), and phenotypic drug susceptibility testing. We calculated the diagnostic accuracy of the different methods and investigated the prevalence and clinical impact of pncA heteroresistance. True pyrazinamide susceptibility status was assigned to each isolate using the Köser classification and expert rules. We observed 100% agreement across genotypic methods for detection of pncA fixed mutations; only TDS confidently identified three isolates (0.8%) with minor variants. For the 355 (99.2%) isolates that could be assigned true pyrazinamide status with confidence, phenotypic DST had a sensitivity of 96.5% (95% confidence interval [CI], 93.8 to 99.3%) and specificity of 100% (95% CI, 100 to 100%), both Sanger sequencing and WGS had a sensitivity of 97.1% (95% CI, 94.6 to 99.6%) and specificity of 97.8% (95% CI, 95.7 to 99.9%), and TDS had sensitivity of 98.8% (95% CI, 97.2 to 100%) and specificity of 97.8% (95% CI, 95.7 to 99.9%). We demonstrate high sensitivity and specificity for pyrazinamide susceptibility testing among all assessed genotypic methods. The prevalence of pyrazinamide heteroresistance in Mycobacterium tuberculosis isolates was lower than that identified for other first-line drugs.

Diagnostic Accuracy of Pyrazinamide Susceptibility Testing in Mycobacterium tuberculosis: A Systematic Review with Meta-Analysis

Introduction: Pyrazinamide (PZA) susceptibility testing plays a critical role in determining the appropriate treatment regimens for multidrug-resistant tuberculosis. We conducted a systematic review and meta-analysis to evaluate the diagnostic accuracy of sequencing PZA susceptibility tests against culture-based susceptibility testing methods as the reference standard. Methods: We searched the MEDLINE/PubMed, Embase, and Web of Science databases for the relevant records. The QUADAS-2 tool was used to assess the quality of the studies. Diagnostic accuracy measures (i.e., sensitivity and specificity) were pooled with a random-effects model. All statistical analyses were performed with Meta-DiSc (version 1.4, Cochrane Colloquium, Barcelona, Spain), STATA (version 14, Stata Corporation, College Station, TX), and RevMan (version 5.3, The Nordic Cochrane Centre, the Cochrane Collaboration, Copenhagen, Denmark) software. Results: A total of 72 articles, published between 2000 and 2019, comprising data for 8,701 isolates of Mycobacterium tuberculosis were included in the final analysis. The pooled sensitivity and specificity of the PZA sequencing test against all reference tests (the combination of BACTEC mycobacteria growth indicator tube 960 (MGIT 960), BACTEC 460, and proportion method) were 87% (95% CI: 85-88) and 94.7% (95% CI: 94-95). The positive likelihood ratio, negative likelihood ratio, diagnostic odds ratio, and the area under the curve estimates were found to be 12.0 (95% CI: 9.0-16.0), 0.17 (95% CI: 0.13-0.21), 106 (95% CI: 71-158), and 96%, respectively. Deek's test result indicated a low likelihood for publication bias (p = 0.01). Conclusions: Our analysis indicated that PZA sequencing may be used in combination with conventional tests due to the advantage of the time to result and in scenarios where culture tests are not feasible. Further work to improve molecular tests would benefit from the availability of standardized reference standards and improvements to the methodology.

Evaluating the Frequency of Resistance to Pyrazinamide Among Drug-resistant Strains of Mycobacterium tuberculosis in Isfahan, Iran

Background: Pyrazinamide is one of the most important first-line medications for the treatment of tuberculosis and an alternative intake for MDR-TB and XDR-TB patients. Objectives: The purpose of this study was to evaluate resistance to pyrazinamide in the isolates resistant to the Mycobacterium tuberculosis drug in patients in the city of Isfahan. Methods: In this study, the drug susceptibility test was performed with pyrazinamide using the proportion method and PZA assay on 47 isolates resistant to Mycobacterium tuberculosis. Then, the mutations of the pncA and rpsA genes of the isolates resistant to pyrazinamide were evaluated by the sequencing method. Results: According to the proportion method, 19 cases were resistant to pyrazinamide, 16 of which had mutations in their pncA and rpsA genes. Besides, five new mutations were recorded, and three isolates lacked mutations in the mentioned genes. Conclusions: Pyrazinamide resistance is high in MDR-TB and INH mono-resistant isolates. Therefore, evaluating the susceptibility to pyrazinamide in patients with MDR-TB before the initiation of treatment with pyrazinamide is considered essential.

Insights Into Mutations Induced Conformational Changes and Rearrangement of Fe2+ Ion in pncA Gene of Mycobacterium tuberculosis to Decipher the Mechanism of Resistance to Pyrazinamide

Pyrazinamide (PZA) is the first-line drug commonly used in treating Mycobacterium tuberculosis (Mtb) infections and reduces treatment time by 33%. This prodrug is activated and converted to an active form, Pyrazinoic acid (POA), by Pyrazinamidase (PZase) enzyme. Mtb resistance to PZA is the outcome of mutations frequently reported in pncA, rpsA, and panD genes. Among the mentioned genes, pncA mutations contribute to 72-99% of the total resistance to PZA. Thus, considering the vital importance of this gene in PZA resistance, its frequent mutations (D49N, Y64S, W68G, and F94A) were investigated through in-depth computational techniques to put conclusions that might be useful for new scaffolds design or structure optimization to improve the efficacy of the available drugs. Mutants and wild type PZase were used in extensive and long-run molecular dynamics simulations in triplicate to disclose the resistance mechanism induced by the above-mentioned point mutations. Our analysis suggests that these mutations alter the internal dynamics of PZase and hinder the correct orientation of PZA to the enzyme. Consequently, the PZA has a low binding energy score with the mutants compared with the wild type PZase. These mutations were also reported to affect the binding of Fe2+ ion and its coordinated residues. Conformational dynamics also revealed that β-strand two is flipped, which is significant in Fe2+ binding. MM-GBSA analysis confirmed that these mutations significantly decreased the binding of PZA. In conclusion, these mutations cause conformation alterations and deformities that lead to PZA resistance.

Drug resistance gene mutations and treatment outcomes in MDR-TB: A prospective study in Eastern China

Background

Multidrug-resistant tuberculosis (MDR-TB) poses a serious challenge to TB control. It is of great value to search for drug resistance mutation sites and explore the roles that they play in the diagnosis and prognosis of MDR-TB.

Methods

We consecutively enrolled MDR-TB patients from five cities in Jiangsu Province, China, between January 2013 and December 2014. Drug susceptibility tests of rifampin, isoniazid, ofloxacin, and kanamycin were routinely performed by proportion methods on Lowenstein–Jensen (LJ) medium. Drug resistance-related genes were sequenced, and the consistency of genetic mutations and phenotypic resistance was compared. The association between mutations and treatment outcomes was expressed as odds ratios (ORs) and 95% confidence intervals (CIs).

Results

Among 87 MDR-TB patients, 71 with treatment outcomes were involved in the analysis. The proportion of successful treatment was 50.7% (36/71). The rpoB gene exhibited the highest mutation rate (93.0%) followed by katG (70.4%), pncA (33.8%), gyrA (29.6%), eis (15.5%), rrs (12.7%), gyrB (9.9%) and rpsA (4.2%). Multivariable analysis demonstrated that patients with pncA gene mutations (adjusted OR: 19.69; 95% CI: 2.43–159.33), advanced age (adjusted OR: 13.53; 95% CI: 1.46–124.95), and nonstandard treatment (adjusted OR: 7.72; 95% CI: 1.35–44.35) had a significantly higher risk of poor treatment outcomes.

Conclusions

These results suggest that Mycobacterium tuberculosis gene mutations may be related to phenotypic drug susceptibility. The pncA gene mutation along with treatment regimen and age are associated with the treatment outcomes of MDR-TB.

Multidrug-resistant tuberculosis (MDR-TB) exacerbates the already serious tuberculosis epidemic, poses a notable threat to global tuberculosis control, and places a considerable burden on developing countries, as treatments for MDR-TB tend to be expensive, of limited efficacy, and toxic. Genotypic determinants of resistance to specific drugs or drug classes offer a rapid and highly specific alternative to phenotypic drug susceptibility testing. Although the relationship between gene mutations and drug resistance has been described previously, the strength of the association of mutations with the treatment outcomes of MDR tuberculosis have not been fully elucidated. The results of our study, which was conducted in a Chinese population, suggest that gene mutations in Mycobacterium tuberculosis may be related to phenotypic drug susceptibility. Mutation of the pncA gene contributes to a poor prognosis and can be applied to predict the treatment outcomes of MDR-TB.

Characterization of pncA Mutations and Prediction of PZA Resistance in Mycobacterium tuberculosis Clinical Isolates From Chongqing, China

Pyrazinamide (PZA) is widely used to treat drug-sensitive or multidrug resistance tuberculosis. However, conventional PZA susceptibility tests of clinical isolates are rather difficult because of the requirement of acid pH. Since resistance to pyrazinamide is primary mediated by mutation of pncA, an alternative way of PZA susceptibility test is to analyze the pyrazinamidase activities of Mycobacterium tuberculosis clinical isolates. Therefore, a database containing the full spectrum of pncA mutations along with pyrazinamidase activities will be beneficial. To characterize mutations of pncA in M. tuberculosis from Chongqing, China, the pncA gene was sequenced and analyzed in 465 clinical isolates. A total of 124 types of mutations were identified in 424 drug-resistant isolates, while no mutation was identified in the 31 pan-susceptible isolates. Ninety-four of the 124 mutations had previously been reported, and 30 new mutations were identified. Based on reported literatures, 294 isolates could be predicted resistant to pyrazinamide. Furthermore, pyrazinamidase activities of the 30 new mutations were tested using the Escherichia coli pncA gene knockout strain. The results showed that 24 of these new mutations (28 isolates) led to loss of pyrazinamidase activity and six (8 isolates) of them did not. Taken together, 322 isolates with pncA mutations could be predicted to be PZA resistant among the 424 drug-resistant isolates tested. Analysis of pncA mutations and their effects on pyrazinamidase activity will not only enrich our knowledge of comprehensive pncA mutations related with PZA resistance but also facilitate rapid molecular diagnosis of pyrazinamide resistance in M. tuberculosis.

A computational perspective on the dynamic behaviour of recurrent drug resistance mutations in the pncA gene from Mycobacterium tuberculosis

Tuberculosis is still one of the top 10 causes of death worldwide, particularly with the emergence of multidrug-resistant tuberculosis. As the most effective first-line anti-tuberculosis drug, pyrazinamide also develops resistance due to the mutation in the pncA gene. Among these mutations, seven mutations at positions F94L, F94S, K96N, K96R, G97C, G97D, and G97S are classified as high-level resistance mutations. However, the resistance mechanism of Mtb to PZA (pyrazinamide) caused by these mutations is still unclear. In this work, we combined molecular dynamics simulation, molecular mechanics/generalized-Born surface area calculation, principal component analysis, and free energy landscape analysis to explore the resistance mechanism of Mtb to PZA due to F94L, F94S, K96N, K96R, G97C, G97D, and G97S mutations, as well as compare interaction changes in wild-type and mutant PZA-bound complexes. The results of molecular mechanics/generalized-Born surface area calculations indicated that the binding free energy of PZA with seven mutants decreased. In mutant systems, the most significant interactions with His137 and Cys138 were lost. Besides, PCA and FEL confirmed significant variations in the protein dynamics during the simulation specifically by altering the Fe2+ binding and its destabilization. Furthermore, mutants also flipped the β-sheet 2, which also affects the binding of Fe2+ and PZA. In the G97D mutant, reported as most lethal, mutation causes the binding pocket to change considerably, so that the position of PZA has a large movement in the binding pocket. In this study, the resistance mechanism of PZA at the atomic level is proposed. The proposed drug-resistance mechanism will provide valuable guidance for the design of anti-tuberculosis drugs.

Whole genome sequence data of Mycobacterium tuberculosis XDR strain, isolated from patient in Kazakhstan

Drug-resistant tuberculosis (TB) is a major public health problem. Clinical Mycobacterium tuberculosis (MTB) isolate with Extensively drug-resistant tuberculosis (MTB-XDR) profile was subjected to whole-genome sequencing using a next-generation sequencing platform (NGS) Roche 454 GS FLX+ followed by bioinformatics sequence analysis. Quality of read was checked by FastQC, paired-end reads were trimmed using Trimmomatic. De novo genome assembly was conducted using Velvet v.1.2.10. The assembled genome of XDR-TB-1599 strain was functionally annotated using the PATRIC platform. Analysis of de novo assembled genome was performed using ResFinder, CARD, CASTB and TB-Profiler tools. MIRU_VNTR genotyping on 12 loci and spoligotyping have been performed for XDR-TB-1599 isolate. M. tuberculosis XDR-TB-1599 strain yielded an average read depth of 21-fold with overall 4 199 325 bp. The assembled genome contains 5528 protein-coding genes, including key drug resistance and virulence-associated genes and GC content of 65.4%. We identified that all proteins encoded by this strain contain conserved domains associated with the first-line anti-tuberculosis drugs such as rifampicin, isoniazid, streptomycin and ethionamide. TB-Profiler had higher average concordance results with phenotypic DST (drug susceptibility testing) in comparison with ResFinder, CARD, CASTB profiling to first-line (75% vs 50%) and second-line (25% vs 0%) of anti-TB drugs, correspondingly. To our knowledge, this is the first report of a highly annotated and characterized whole-genome sequence and de novo assembled XDR-TB M.tuberculosis strain isolated from a sputum of new TB case-patient from Kazakhstan performed on Roche 454 GS FLX+ platform. This report highlights an important role of whole-genome sequencing technology and analysis as an advanced approach for drug-resistance investigations of circulated TB isolates.

Mechanistic analysis of A46V, H57Y, and D129N in pyrazinamidase associated with pyrazinamide resistance

Pyrazinamide (PZA) is a component of first-line drugs, active against latent Mycobacterium tuberculosis (MTB) isolates. The prodrug is activated into the active form, pyrazinoic acid (POA) via pncA gene-encoded pyrazinamidase (PZase). Mutations in pncA have been reported, most commonly responsible for PZA-resistance in more than 70% of the resistant cases. In our previous study, we detected many mutations in PZase among PZA-resistance MTB isolates including A46V, H71Y, and D129N. The current study was aimed to investigate the molecular mechanism of PZA-resistance behind mutants (MTs) A46V, H71Y, and D129N in comparison with the wild type (WT) through molecular dynamic (MD) simulation. MTB positive samples were subjected to PZA drug susceptibility testing (DST) against critical concentration (100ug/ml). The resistant samples were subjected to pncA sequencing. Thirty-six various mutations have been observed in the coding region of pncA of PZA-resistant isolates (GenBank accession No. MH461111) including A46V, H71Y, and D129N. The post-simulation analysis revealed a significant variation in MTs structural dynamics as compared to the WT. Root means square deviations (RMSD) and Root means square fluctuation (RMSF) has been found in variation between WT and MTs. Folding effect and pocket volume were altered in MTs when compared with WT. Geometric matching supports the effect of mutation A46V, H71Y, and D129N on PZase structure that may have an insight effect on PZase dynamics, making them vulnerable to convert pro-PZA into active form, POA. In conclusion, the current analyses will provide useful information behind PZA-resistance for better management of drug-resistant TB.

Pyrazinamide resistance of novel mutations in pncA and their dynamic behavior

Pyrazinamide (PZA) is one of the essential anti-mycobacterium drugs, active against non-replicating Mycobacterium tuberculosis (MTB) isolates. PZA is converted into its active state, called pyrazinoic acid (POA), by action of pncA encoding pyrazinamidase (PZase). In the majority of PZA-resistance isolates, pncA harbored mutations in the coding region. In our recent report, we detected a number of novel variants in PZA-resistance (PZAR) MTB isolates, whose resistance mechanisms were yet to be determined. Here we performed several analyses to unveil the PZAR mechanism of R123P, T76P, G150A, and H71R mutants (MTs) through molecular dynamics (MD) simulations. In brief, culture positive MTB isolates were subjected to PZA susceptibility tests using the WHO recommended concentration of PZA (100 μg ml-1). The PZAR samples were screened for mutations in pncA along sensitive isolates through polymerase chain reactions and sequencing. A large number of variants (GeneBank accession no. MH461111), including R123P, T76P, G150A, and H71R, have been spotted in more than 70% of isolates. However, the mechanism of PZAR for mutants (MTs) R123P, T76P, G150A, and H71R was unknown. For the MTs and native PZase structures (WT), thermodynamic properties were compared using molecular dynamics simulations for 100 ns. The MTs structural activity was compared to the WT. Folding effect and pocket volume variations have been detected when comparing between WT and MTs. Geometric matching further confirmed the effect of R123P, T76P, G150A, and H71R mutations on PZase dynamics, making them vulnerable for activating the pro-drug into POA. This study offers a better understanding for management of PZAR TB. The results may be used as alternative diagnostic tools to infer PZA resistance at a structural dynamics level.

Mycobacterium tuberculosis ribosomal protein S1 (RpsA) and variants with truncated C-terminal end show absence of interaction with pyrazinoic acid

Pyrazinamide (PZA) is an antibiotic used in first- and second-line tuberculosis treatment regimens. Approximately 50% of multidrug-resistant tuberculosis and over 90% of extensively drug-resistant tuberculosis strains are also PZA resistant. Despite the key role played by PZA, its mechanisms of action are not yet fully understood. It has been postulated that pyrazinoic acid (POA), the hydrolyzed product of PZA, could inhibit trans-translation by binding to Ribosomal protein S1 (RpsA) and competing with tmRNA, the natural cofactor of RpsA. Subsequent data, however, indicate that these early findings resulted from experimental artifact. Hence, in this study we assess the capacity of POA to compete with tmRNA for RpsA. We evaluated RpsA wild type (WT), RpsA ∆A438, and RpsA ∆A438 variants with truncations towards the carboxy terminal end. Interactions were measured using Nuclear Magnetic Resonance spectroscopy (NMR), Isothermal Titration Calorimetry (ITC), Microscale Thermophoresis (MST), and Electrophoretic Mobility Shift Assay (EMSA). We found no measurable binding between POA and RpsA (WT or variants). This suggests that RpsA may not be involved in the mechanism of action of PZA in Mycobacterium tuberculosis, as previously thought. Interactions observed between tmRNA and RpsA WT, RpsA ∆A438, and each of the truncated variants of RpsA ∆A438, are reported.

Structural insights of catalytic mechanism in mutant pyrazinamidase of Mycobacterium tuberculosis

Pyrazinamidase (PZase) is a member of Fe-dependent amidohydrolases that activates pyrazinamide (PZA) into active pyrazinoic acid (POA). PZA, a nicotinamide analogue, is an essential first-line drug used in Mycobacterium tuberculosis (Mtb) treatment. The active form of PZA, POA, is toxic and potently inhibits the growth of latent Mtb, which makes it possible to shorten the conventional 9-month tuberculosis treatment to 6 months. In this study, an extensive molecular dynamics simulation was carried out to the study the resistance mechanism offered by the three mutations Q10P and D12A and G97D. Our results showed that two regions Gln10-His43, Phe50-Gly75 are profoundly affected by these mutations. Among the three mutations, Q10P and D12A mutations strongly disturb the communication among the catalytic triad (Asp8, Lys98 and Cys138). The oxyanion hole is formed between the backbone nitrogen atoms of A134 and C138 which stabilizes the hydroxyl anion of nicotinamide. The D12A mutation greatly disturbs the oxyanion hole formation followed by the Q10P and G97D. Our results also showed that these mutations destabilize the interaction between Fe²⁺ ion and Asp49, His51, His57 and His71. The binding pocket analysis showed that these mutations increase the cavity volume, which results in loose binding of PZA. MMGBSA analyzes have shown that these mutations reduce the binding affinity to the PZA drug. Our results may provide useful information for the design of new and effective PZase inhibitors based on structural information of WT and mutant PZases.Communicated by Ramaswamy H. Sarma.

Draft Genome Sequence of an Extensively Drug-Resistant Mycobacterium tuberculosis Clinical Isolate, 3485_MTB, from Nur-Sultan, Kazakhstan

Here, we report the draft genome sequence of an extensively drug-resistant Mycobacterium tuberculosis clinical isolate, 3485_MTB, from Nur-Sultan, Kazakhstan. The genome sequence is composed of 4,836,003 bp. The genome will provide more data on the genetic variations occurring in local drug-resistant isolates.

Here, we report the draft genome sequence of an extensively drug-resistant Mycobacterium tuberculosis clinical isolate, 3485_MTB, from Nur-Sultan, Kazakhstan. The genome sequence is composed of 4,836,003 bp. The genome will provide more data on the genetic variations occurring in local drug-resistant isolates.

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.

Dynamics Insights Into the Gain of Flexibility by Helix-12 in ESR1 as a Mechanism of Resistance to Drugs in Breast Cancer Cell Lines

Incidents of breast cancer (BC) are on the rise on a daily basis and have proven to be the most prevelant cause of death for women in both developed and developing countries. Among total BC cases diagnosed after menopause, 70% of cases are Estrogen Receptor (ER) positive (ER-positive or ER+). Mutations in the LBD (ligand-binding domain) of the ER have recently been reported to be the major cause of resistance to potent antagonists. In this study, the experimentally reported mutations K303R, E380Q, V392I, S463P, V524E, P535H, P536H, Y537C, Y537N, Y537S, and D538G were analyzed, and the most significant mutations were shortlisted based on multiple analyses. Initial analyses, such as mCSM stability, occluded depth analysis, mCSM-binding affinity, and FoldX energy changes shortlisted only six mutations as being highly resistant. Finally, simulations of force field-based molecular dynamics (MD on wild type (WT) ERα) on six mERα variants (E380Q, S463P, Y537S, Y537C, Y537N, and D538G) were carried out to justify mechanism of the resistance. It was observed that these mutations increased the flexibility of the H12. A bonding analysis suggested that previously reported important residue His524 lost bonding upon mutation. Other parameters, such as PCA (principal component analysis), DCCM (dynamics cross-correlation), and FEL (free energy landscape), verified that the shortlisted mutations affect the H12 helix, which opens up the co-activator binding conformation. These results provide deep insight into the mechanism of relative resistance posed to fulvestrant due to mutations in breast cancer. This study will facilitate further understanding of the important aspects of designing specific and more effective drugs.

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.

Value of pyrazinamide for composition of new treatment regimens for multidrug-resistant Mycobacterium tuberculosis in China

Background

Pyrazinamide still may be a useful drug for treatment of rifampin-resistant (RR-TB) or multidrug-resistant tuberculosis (MDR-TB) in China while awaiting scale up of new drugs and regimens including bedaquiline and linezolid. The level of pyrazinamide resistance among MDR-TB patients in China is not well established. Therefore, we assessed pyrazinamide resistance in a representative sample and explored determinants and patterns of pncA mutations.

Methods

MDR-TB isolates from the 2007 national drug resistance survey of China were sub-cultured and examined for pyrazinamide susceptibility by BACTEC MGIT 960 method. pncA mutations were identified by sequencing. Characteristics associated with pyrazinamide resistance were analyzed using univariable and multivariable log-binominal regression.

Results

Of 401 MDR-TB isolates, 324 were successfully sub-cultured and underwent drug susceptibility testing. Pyrazinamide resistance was prevalent in 40.7% of samples, similarly among new and previously treated MDR-TB patients. Pyrazinamide resistance in MDR-TB patients was associated with lower age (adjusted OR 0.54; 95% CI, 0.34–0.87 for those aged ≧60 years compared to < 40 years). Pyrazinamide resistance was not associated with gender, residential area, previous treatment history and Beijing genotype. Of 132 patients with pyrazinamide resistant MDR-TB, 97 (73.5%) had a mutation in the pncA gene; with 61 different point mutations causing amino acid change, and 11 frameshifts in the pncA gene. The mutations were scattered throughout the whole pncA gene and no hot spot region was identified.

Conclusions

Pyrazinamide resistance among MDR-TB patients in China is common, although less so in elderly patients. Therefore, pyrazinamide should only be used for treatment of RR/MDR-TB in China if susceptibility is confirmed. Molecular testing for detection of pyrazinamide resistance only based on pncA mutations has certain value for the rapid detection of pyrazinamide resistance in MDR-TB strains but other gene mutations conferring to pyrazinamide resistance still need to be explored to increase its predictive ability .

Structural Dynamics Behind Clinical Mutants of PncA-Asp12Ala, Pro54Leu, and His57Pro of Mycobacterium tuberculosis Associated With Pyrazinamide Resistance

Pyrazinamide (PZA) is one of the main FDA approved drugs to be used as the first line of defense against Mycobacterium Tuberculosis (MTB). It is activated into pyrazinoic acid (POA) via MTB's pncA gene-encoded pyrazinamidase (PZase). Mutations are most commonly responsible for PZA-resistance in nearly 70% of the resistant samples. In the present work, MTB positive samples were chosen for PZA drug susceptibility testing (DST) against critical concentration (100 ug/ml) of PZA. The resistant samples were subjected to pncA sequencing. As a result, 36 various mutations have been observed in the PZA resistant samples, uploaded to the NCBI (GeneBank accession no. MH461111). Here we report the mechanism of PZA resistance behind the three mutants (MTs), Asp12Ala, Pro54Leu, and His57Pro in comparison with the wild type (WT) through molecular dynamics simulation to unveil how these mutations affect the overall conformational stability. The post-simulation analyses revealed notable deviations as compared to the WT structure. Molecular docking studies of PZA with MTs and WT, pocket volume inspection and overall shape complementarity analysis confirmed the deleterious nature of these mutations and gave an insight into the mechanism behind PZA-resistance. These analyses provide vital information regarding MTB drug resistance and could be extremely useful in therapy management and overcoming its global burden.

Resistance To First-Line Antituberculosis Drugs And Prevalence Of pncA Mutations In Clinical Isolates Of Mycobacterium tuberculosis From Zunyi, Guizhou Province Of China

Background

China is one of the high-burden countries for multidrug-resistant tuberculosis (MDR-TB), and pyrazinamide is one of the anti-TB drugs used for the shorter MDR-TB treatment regimen. The aim of this study was to determine the correlation between pncA gene mutations and resistance to four first-line anti-TB drugs as well as treatment history in clinical isolates of Mycobacterium tuberculosis.

Patients and methods

M. tuberculosis clinical isolates were collected from 318 in-patients with smear-positive TB between October 2008 and September 2016 at a major hospital in Zunyi, Guizhou Province of China, and used for drug susceptibility testing against four first-line anti-TB drugs. Genomic DNA extracted from clinical isolates was used for PCR amplification and DNA sequencing of the pncA gene.

Results

Among 318 clinical isolates, 129 (40.6%), 170 (53.5%), 66 (20.8%) and 109 (34.3%) were resistant to rifampicin, isoniazid, ethambutol and streptomycin respectively. In addition, 124 clinical isolates were MDR-TB and 71.8% of them were previously treated cases. Sequencing results showed that 46.8% of MDR-TB and 2.2% of drug susceptible isolates harbored a pncA mutation, and 52 types of pncA mutations were detected from 64 isolates. The prevalence of pncA mutations in isolates resistant to first-line anti-TB drugs and previously treated TB cases was significantly higher than that in drug-susceptible isolates and new cases of TB.

Conclusion

High prevalence of pncA mutations in clinical isolates of M. tuberculosis from Zunyi, Guizhou Province of China, is correlated with resistance to four first-line anti-TB drugs, MDR-TB and previously treated TB cases.

Finding New Ways to Combat Multidrug-Resistant Tuberculosis

Tuberculosis (TB) is a major issue in global health and affects millions of people each year. Multidrug-resistant tuberculosis (MDR-TB) annually causes many deaths worldwide. Development of a way to diagnose and treat patients with MDR-TB can potentially reduce the incidence of the disease. The current study reviews the risk factors, pattern of progression, mechanism of resistance, and interaction between bacteria and the host immune system, which disrupts the immune response. It also targets the components of Mycobacterium tuberculosis (Mtb) and diagnosis and treatment options that could be available for clinical use in the near future. Mutations play an important role in development of MDR-TB and the selection of appropriate mutations can help to understand the type of resistance in patients to anti-TB drugs. In this way, they can be initially treated with proper and effective therapeutic choices, which can accelerate the course of treatment and improve patient health. Targeting the components and enzymes of Mtb is necessary for understanding bacterial survival and finding a way to destroy the pathogen and allow patients to recover faster and prevent the spread of disease, especially resistant strains.

Structural dynamics behind variants in pyrazinamidase and pyrazinamide resistance

Pyrazinamide (PZA) is an important component of first-line anti-tuberculosis (anti-TB) drugs. The anti-TB agent is activated into an active form, pyrazinoic acid (POA), by Mycobacterium tuberculosis (MTB) pncA gene encoding pyrazinamidase (PZase). The major cause of PZA-resistance has been associated with mutations in the pncA gene. We have detected several novel mutations including V131F, Q141P, R154T, A170P, and V180F (GeneBank Accession No. MH461111) in the pncA gene of PZA-resistant isolates during PZA drug susceptibility testing followed by pncA gene sequencing. Here, we investigated molecular mechanism of PZA-resistance by comparing the results of experimental and molecular dynamics. The mutants (MTs) and wild type (WT) PZase structures in apo and complex with PZA were subjected to molecular dynamic simulations (MD) at the 40 ns. Multiple factors, including root mean square deviations (RMSD), binding pocket, total energy, dynamic cross correlation, and root mean square fluctuations (RMSF) of MTs and WT were compared. The MTs attained a high deviation and fluctuation compared to WT. Binding pocket volumes of the MTs, were found, lower than the WT, and the docking scores were high than WT while shape complementarity scores were lower than that of the WT. Residual motion in MTs are seemed to be dominant in anti-correlated motion. Mutations at locations, V131F, Q141P, R154T, A170P, and V180F, might be involved in the structural changes, possibly affecting the catalytic property of PZase to convert PZA into POA. Our study provides useful information that will enhance the understanding for better management of TB. AbbreviationsDSTdrug susceptibility testingΔelecelectrostatic energyLJLowenstein-Jensen mediumMGITmycobacterium growth indicator tubesMTsmutantsMDmolecular dynamic simulationsMTBMycobacterium tuberculosisNALC-NaOHN-acetyl-l-cysteine-sodium hydroxideNIHNational Institutes of HealthNPTamount of substance (N), pressure (P) temperature (T)NVTmoles (N), volume (V) temperature (T)PZasepyrazinamidaseΔpspolar solvation energyPTRLProvincial Tuberculosis Reference LaboratoryRMSDroot mean square deviationsRMSFroot mean square fluctuationsΔSASAsolvent accessible surface area energyTBtuberculosisGTotaltotal binding free energyΔvdWVan der Waals energyWTwild typeCommunicated by Ramaswamy H. Sarma.

Structural and free energy landscape of novel mutations in ribosomal protein S1 (rpsA) associated with pyrazinamide resistance

Resistance to key first-line drugs is a major hurdle to achieve the global end tuberculosis (TB) targets. A prodrug, pyrazinamide (PZA) is the only drug, effective in latent TB, recommended in drug resistance and susceptible Mycobacterium tuberculosis (MTB) isolates. The prodrug conversion into active form, pyrazinoic acid (POA), required the activity of pncA gene encoded pyrazinamidase (PZase). Although pncA mutations have been commonly associated with PZA resistance but a small number of resistance cases have been associated with mutationss in RpsA protein. Here in this study a total of 69 PZA resistance isolates have been sequenced for pncA mutations. However, samples that were found PZA resistant but pncA wild type (pncAWT), have been sequenced for rpsA and panD genes mutation. We repeated a drug susceptibility testing according to the WHO guidelines on 18 pncAWT MTB isolates. The rpsA and panD genes were sequenced. Out of total 69 PZA resistant isolates, 51 harbored 36 mutations in pncA gene (GeneBank Accession No. MH46111) while, fifteen different mutations including seven novel, were detected in the fourth S1 domain of RpsA known as C-terminal (MtRpsACTD) end. We did not detect any mutations in panD gene. Among the rpsA mutations, we investigated the molecular mechanism of resistance behind mutations, D342N, D343N, A344P, and I351F, present in the MtRpsACTD through molecular dynamic simulations (MD). WT showed a good drug binding affinity as compared to mutants (MTs), D342N, D343N, A344P, and I351F. Binding pocket volume, stability, and fluctuations have been altered whereas the total energy, protein folding, and geometric shape analysis further explored a significant variation between WT and MTs. In conclusion, mutations in MtRpsACTD might be involved to alter the RpsA activity, resulting in drug resistance. Such molecular mechanism behind resistance may provide a better insight into the resistance mechanism to achieve the global TB control targets.

Pyrazinamide resistance and mutations L19R, R140H, and E144K in Pyrazinamidase of Mycobacterium tuberculosis

Pyrazinamide (PZA) is an important component of first-line antituberculosis drugs activated by Mycobacterium tuberculosis pyrazinamidase (PZase) into its active form pyrazinoic acid. Mutations in the pncA gene have been recognized as the major cause of PZA resistance. We detected some novel mutations, Leucine19Arginine (L19R), Arginine140Histidine (R140H), and Glutamic acid144 Lysine (E144K), in the pncA gene of PZA-resistant isolates in our wet lab PZA drug susceptibility testing and sequencing. As the molecular mechanism of resistance of these variants has not been reported earlier, we have performed multiple analyses to unveil different mechanisms of resistance because of PZase mutations L19R, R140H, and E144K. The mutants and native PZase structures were subjected to comprehensive computational molecular dynamics (MD) simulations at 100 nanoseconds in apo and drug-bound form. Mutants and native PZase binding pocket were compared to observe the consequence of mutations on the binding pocket size. Hydrogen bonding, Gibbs free energy, and natural ligand Fe ⁺² effect were also analyzed between native and mutants. A significant variation between native and mutant PZase structure activity was observed. The native PZase protein docking score was found to be the maximum, showing strong binding affinity in comparison with mutants. MD simulations explored the effect of the variants on the biological function of PZase. Hydrogen bonding, metal ion Fe ⁺² deviation, and fluctuation also seemed to be affected because of the mutations L19R, R140H, and E144K. The variants L19R, R140H, and E144K play a significant role in PZA resistance, altering the overall activity of native PZase, including metal ion Fe ⁺² displacement and free energy. This study offers valuable evidence for better management of drug-resistant tuberculosis.

Insights into the Mechanisms of the Pyrazinamide Resistance of Three Pyrazinamidase Mutants N11K, P69T, and D126N

In an effort to discover the mechanism of resistance offered by Mycobacterium tuberculosis (Mtb) toward the pyrazinamide (PZA) drug, an extensive molecular dynamics strategy was employed. PZA is a first-line prodrug that effectively cuts therapy time by 33% (from 9 to 6 months). Pyrazinamidase enzyme (PZase), encoded by the pncA gene, is responsible for the activation of prodrug PZA into pyrazinoic acid (POA). POA is toxic and potently inhibits the growth of latent Mtb even at low pH values. PZA resistance is caused by three genes pncA, rpsA, and panD. Among them, the pncA gene contributes 72-99% to the resistance. Hence, the present study focused on the novel mutations N11K, P69T, and D126N in the pncA gene. In the present study, the possible mechanism of these three mutations was studied through molecular dynamics simulation and docking techniques. Our in-depth analysis and results are in strong agreement with our experimental observation. The binding pocket analysis showed that mutations decrease the volume of the active site and hinder the correct orientation of PZA drug in the active site. Moreover, the Patchdock score was found to be low as compared to WT showing the disturbance of shape complementarity between PZase and PZA drug. These mutations were found to disturb the position of the Fe²⁺ ion. Among the mutations, D126N allosterically disturbed the position of the Fe²⁺ ion. MMGBSA analyses showed that these mutations decrease the binding affinity toward the PZA drug. In conclusion, mutations N11K, P69T, and D126N result in weak binding affinity with PZA and also cause significant structural deformations that lead to PZA resistance. This study provides useful information that mutations in other than active parts may also cause protein folding and ligand displacement effects, altering the biological functions.

Insight into novel clinical mutants of RpsA-S324F, E325K, and G341R of Mycobacterium tuberculosis associated with pyrazinamide resistance

Pyrazinamide (PZA) is an important component of first-line anti-tuberculosis drugs which is converted into active form, pyrazinoic acid (POA), by Mycobacterium tuberculosis (MTB) pncA gene encoded, pyrazinamidase (PZase). Mutations in pncA are detected in >70% of PZA resistant isolates but, noticeably, not in all. In this study, we selected 18 PZA-resistant but wild type pncA (pncAWT) MTB isolates. Drug susceptibility testing (DST) of all the isolates were repeated at the critical concentration of PZA drug. All these PZA-resistance but pncAWT isolates were subjected to RpsA sequencing. Fifteen different mutations were identified in eleven isolates, where seven were present in a conserved region including, Ser324Phe, Glu325Lys, Gly341Arg. As the molecular mechanism of resistance behind these variants has not been reported earlier, we have performed multiple analysis to unveil the mechanisms of resistance behind mutations S324F, E325K, and G341R. The mutant and wild type RpsA structures were subjected to comprehensive computational molecular dynamic simulations at 50 ns. Root mean square deviation (RMSD), Root mean square fluctuation (RMSF), and Gibbs free energy of mutants were analyzed in comparison with wild type. Docking score of wild type-RpsA has been found to be maximum, showing a strong binding affinity in comparison with mutants. Pocket volume, RMSD and RMSF have also been found to be altered, whereas total energy, folding effect (radius of gyration) and shape complimentarily analysis showed that variants S324F, E325K, and G341R have been playing a significant role behind PZA-resistance. The study offers valuable information for better management of drug resistance tuberculosis.

pncA Gene Mutations Associated with Pyrazinamide Resistance in Drug-Resistant Tuberculosis, South Africa and Georgia

Although pyrazinamide is commonly used for tuberculosis treatment, drug-susceptibility testing is not routinely available. We found polymorphisms in the pncA gene for 70% of multidrug-resistant and 96% of extensively drug-resistant Mycobacterium tuberculosis isolates from South Africa and Georgia. Assessment of pyrazinamide susceptibility may be prudent before using it in regimens for drug-resistant tuberculosis.

Pyrazinoic Acid Inhibits Mycobacterial Coenzyme A Biosynthesis by Binding to Aspartate Decarboxylase PanD

Previously, we showed that a major in vitro and in vivo mechanism of resistance to pyrazinoic acid (POA), the bioactive component of the critical tuberculosis (TB) prodrug pyrazinamide (PZA), involves missense mutations in the aspartate decarboxylase PanD, an enzyme required for coenzyme A biosynthesis. What is the mechanism of action of POA? Upon demonstrating that treatment of M. bovis BCG with POA resulted in a depletion of intracellular coenzyme A and confirming that this POA-mediated depletion is prevented by either missense mutations in PanD or exogenous supplementation of pantothenate, we hypothesized that POA binds to PanD and that this binding blocks the biosynthetic pathway. Here, we confirm both hypotheses. First, metabolomic analyses showed that POA treatment resulted in a reduction of the concentrations of all coenzyme A precursors downstream of the PanD-mediated catalytic step. Second, using isothermal titration calorimetry, we established that POA, but not its prodrug PZA, binds to PanD. Binding was abolished for mutant PanD proteins. Taken together, these findings support a mechanism of action of POA in which the bioactive component of PZA inhibits coenzyme A biosynthesis via binding to aspartate decarboxylase PanD. Together with previous works, these results establish PanD as a genetically, metabolically, and biophysically validated target of PZA.

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

Pyrazinamide (PZA) is a critical component of first- and second-line treatments of tuberculosis (TB), yet its mechanism of action largely remains an enigma. We carried out a genetic screen to isolate Mycobacterium bovis BCG mutants resistant to pyrazinoic acid (POA), the bioactive derivative of PZA, followed by whole genome sequencing of 26 POA resistant strains. Rather than finding mutations in the proposed candidate targets fatty acid synthase I and ribosomal protein S1, we found resistance conferring mutations in two pathways: missense mutations in aspartate decarboxylase panD, involved in the synthesis of the essential acyl carrier coenzyme A (CoA), and frameshift mutations in the vitro nonessential polyketide synthase genes mas and ppsA-E, involved in the synthesis of the virulence factor phthiocerol dimycocerosate (PDIM). Probing for cross resistance to two structural analogs of POA, nicotinic acid and benzoic acid, showed that the analogs share the PDIM- but not the CoA-related mechanism of action with POA. We demonstrated that POA depletes CoA in wild-type bacteria, which is prevented by mutations in panD. Sequencing 10 POA-resistant Mycobacterium tuberculosis H37Rv isolates confirmed the presence of at least 2 distinct mechanisms of resistance to the drug. The emergence of resistance through the loss of a virulence factor in vitro may explain the lack of clear molecular patterns in PZA-resistant clinical isolates, other than mutations in the prodrug-converting enzyme. The apparent interference of POA with virulence pathways may contribute to the drug's excellent in vivo efficacy compared to its modest in vitro potency.

Pyrazinamide resistance in Mycobacterium tuberculosis: Review and update

The global control and management of tuberculosis (TB) is faced with the formidable challenge of worsening scenarios of drug-resistant disease. Pyrazinamide (PZA) is an indispensable first-line drug used for the treatment of TB. It plays a key role in reducing TB relapse rates, shortening the course of the disease treatment from 9-12 months to 6 months, and the treatment of patients infected with bacillary strains that are resistant to at least isoniazid and rifampicin. Additionally, it is the only first-line anti-TB drug most likely to be maintained in all new regimens, which are aimed at reducing the treatment period of susceptible, multi-drug resistant and extensively drug-resistant TB. It has a preferential sterilizing activity against non-replicating persister bacilli with low metabolism at acid pH in vitro or in vivo during active inflammation where other drugs may not act so well. PZA seem to have a non-specific cellular target and instead, exerts its anti-mycobacterial effect by disrupting the membrane energetics, the trans-translation process, acidification of the cytoplasm and perhaps coenzyme A synthesis, which is required for survival of Mycobacterium tuberculosis (MTB) persisters. Indeed, the emergence of MTB strains resistant to PZA represents an important clinical and public health problem. The essential role of PZA in TB treatment underlines the need for accurate and rapid detection of its resistance. This article presents an updated review of the molecular mechanisms of drug action and resistance in MTB against PZA, commenting on the several research gaps and proposed drug targets for PZA.

Fighting an old disease with modern tools: characteristics and molecular detection methods of drug-resistant Mycobacterium tuberculosis

Tuberculosis (TB) is an ancient disease, but not a disease of the past. The increasing prevalence of drug-resistant strains of Mycobacterium tuberculosis, the causative agent of TB, demands new measures to combat the situation. Rapid and accurate detection of the pathogen, and its drug susceptibility pattern, is essential for timely initiation of treatment, and ultimately, control of the disease. Molecular-based methods offer a great chance to improve detection of drug-resistant TB; however, their development and usage should be accompanied with a profound understanding of drug resistance mechanisms and circulating M. tuberculosis strains in specific settings, as otherwise, the usefulness of such tests may be limited. This review gives an overview of the history of TB treatment and drug resistance, drug resistance mechanisms for the most commonly used drugs and molecular methods designed to detect drug-resistant strains.

Molecular snapshot of Mycobacterium tuberculosis population in Kazakhstan: a country-wide study

Republic of Kazakhstan is among the 27 high multidrug-resistant tuberculosis (MDR-TB) burden countries in the world. Here, we analyzed the population structure and phylogeography of Mycobacterium tuberculosis in Kazakhstan and impact of the identified genotypes on spread of drug resistant strains. A total of 159 M. tuberculosis isolates from different regions of Kazakhstan were typed using 24-MIRU-VNTR and spoligotyping, and the profiles were compared to the MIRU-VNTRplus and SITVIT_WEB databases. Eight isolates with double VNTR alleles were excluded from further analysis that was performed on 151 isolates. They were assigned to 10 families, Beijing (n = 109) being the largest and dominated by a single clonal cluster 94-32 and derived profiles (n = 101). The other families were represented mainly by LAM (n = 17), Ural (n = 8), NEW-1 (n = 3) and a new cluster named KAZ-1 (n = 8). Beijing, LAM and Ural isolates were detected in all parts of the country while Iran-specific family NEW-1 was found only in southern Kazakhstan (P = 0.001). A reduced scheme of 10 most polymorphic VNTR loci provided a discrimination similar to that achieved by 15-MIRU scheme and may be recommended for rapid preliminary screening of the clinical isolates in Kazakhstan. Multi-drug resistance was significantly more prevalent among Beijing (64/109) and LAM (7/17) strains compared to strains of other families (1/25; P = 0.0006 and 0.01, respectively). High prevalence of the genetically closely related MDR strains of the Beijing genotype found in different regions of Kazakhstan highlights their crucial impact on the current TB epidemic in this country.