An integrated computational investigation to unveil the structural impacts of mutation on the InhA structural gene of Mycobacterium tuberculosis

Abstract

Growing concern about the difficulty in diagnosis and treatments of drug-resistant tuberculosis falls under the major global health issues. There is an urgent need for finding novel strategies to develop drugs or bioactive molecules against the global threat of Mycobacterium tuberculosis (MTB). Isoniazid (INH) is a front line drug against tuberculosis; it primarily targets the enoyl-acyl carrier protein reductase (InhA), a potent drug target in the mycolic acid pathway of MTB. To gain deeper insight into the impact of INH resistant mutation and its influence on the structural dynamics of InhA, combined conformational dynamics and residue interaction network (RIN) studies were performed. The molecular dynamics investigation provided a hint about the structural changes altering protein activity. The principal component analysis (PCA) based free energy landscape plot highlighted the highest stable part of wild-type (WT) and mutant structures. Intriguingly, the mutation at the 78th position of InhA from its native residue valine to alanine increases the structural stability with higher NADH binding affinity. The MM-PBSA based binding energy calculations confirm that electrostatic interactions played a critical role in the binding of NADH at the binding site of InhA. The calculated binding energy score, as well as potential hydrogen bonds and salt bridge networks, proved the strong binding of mutant InhA as compared to WT. Further, the mutation potentially altered the protein network topology, thereby subsequently affected the landscape of NADH binding. The present study is an attempt to understand the structural and functional impact associated with a drug-resistant mutation (V78A) thus it will be helpful in designing potent inhibitors against drug-resistant tuberculosis.