Abstract

The enoyl-ACP reductase enzyme (InhA) from M. tuberculosis is recognized as the primary target of isoniazid (INH), a first-line antibiotic for tuberculosis treatment. To identify the specific interactions of INH-NAD adduct and its derivative adducts in InhA binding pocket, molecular docking calculations and quantum chemical calculations were performed on a set of INH derivative adducts. Reliable binding modes of INH derivative adducts in the InhA pocket were established using the Autodock 3.05 program, which shows a good ability to reproduce the X-ray bound conformation with rmsd of less than 1.0 Å. The interaction energies of the INH-NAD adduct and its derivative adducts with individual amino acids in the InhA binding pocket were computed based on quantum chemical calculations at the MP2/6-31G (d) level. The molecular docking and quantum chemical calculation results reveal that hydrogen bond interactions are the main interactions for adduct binding. To clearly delineate the linear relationship between structure and activity of these adducts, CoMFA and CoMSIA models were set up based on molecular docking alignment. The resulting CoMFA and CoMSIA models are in conformity with the best statistical qualities, in which r2cv is 0.67 and 0.74, respectively. Structural requirements of isoniazid derivatives that can be incorporated into the isoniazid framework to improve the activity have been identified through CoMFA and CoMSIA steric and electrostatic contour maps. The integrated results from structure-based, ligand-based design approaches and quantum chemical calculations provide useful structural information facilitating the design of new and more potentially effective antitubercular agents as follow: the R substituents of isoniazid derivatives should contain a large plane and both sides of the plane should contain an electropositive group. Moreover, the steric and electrostatic fields of the 4-pyridyl ring are optimal for greater potency.

Highlights

  • Tuberculosis (TB) caused by Mycobacterium tuberculosis (M. tuberculosis) still remains a major cause of illness and death worldwide, especially in Asia and Africa

  • The docked conformation of INH-NAD adduct is close to the binding mode found in the X-ray crystal structure with rmsd of 0.44 Å, indicating that the docking parameters are reasonable to generate the binding mode of INH-NAD adduct in the inhibits a 2-trans-enoyl-acyl carrier protein reductase (InhA) binding pocket

  • Based on the interaction energy calculation results, all high attraction energies correlate with the hydrogen bond interactions of the INH-NAD adduct in the InhA binding pocket observed in the X-ray crystal structure

Read more

Summary

Introduction

Tuberculosis (TB) caused by Mycobacterium tuberculosis (M. tuberculosis) still remains a major cause of illness and death worldwide, especially in Asia and Africa. First-line drugs are mainly bactericidal and combine a high degree of efficacy with a relatively low toxicity to the patient during treatment. These drugs include isoniazid, rifampicin, streptomycin, ethambutol, pyrazinamide, and fluoroquinolones. Isoniazid (INH) has the greatest bactericidal activity and is used almost from the outset of tuberculosis chemotherapy [11,12,13] This antibiotic inhibits a 2-trans-enoyl-acyl carrier protein reductase (InhA) displaying a long-chain fatty acid elongation activity. The reactive species generated from the activation process forms a covalent adduct with NAD+ that is a potent inhibitor of InhA [20,23,24].

Validation of the molecular docking calculations
Molecular docking analysis of the highly active compounds
Docking analysis of the moderately active compounds
Docking analysis of the weakly active compounds
The favorable interactions for binding of INH derivative adducts
Interaction energy of the INH-NAD adduct in the InhA binding pocket
Interaction energy of the highly active compounds in the InhA binding pocket
Interaction energy of the weakly active compounds in the InhA binding pocket
CoMFA and CoMSIA models
2.10. CoMFA and CoMSIA contour analysis of INH derivative adducts
2.11. The steric contour analysis of INH derivative adducts
2.12. The electrostatic contour analysis of INH derivative adducts
2.13. The structural requirement of the R substituent of INH derivative adducts
Biological activity data
Geometry optimization
Molecular docking calculations
Interaction energy calculations
Training and test sets
Molecular alignment rules for CoMFA and CoMSIA modeling
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call