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.