Elucidating the structural basis of diphenyl ether derivatives as highly potent enoyl-ACP reductase inhibitors through molecular dynamics simulations and 3D-QSAR study

Identifying the Structural Features of Diphenyl Ether Analogues for InhA Inhibition: A 2D and 3D QSAR Based Study

Background:

Enoyl acyl carrier protein reductase (InhA) is a validated target for Mycobacterium. It is an enzyme which is associated with the biosynthesis of mycolic acids in type II fatty acid synthase system. Mycobacterial cell wall majorly comprises mycolic acids, which are responsible for virulence of the microorganism. Several diphenyl ether derivatives have been known to be direct inhibitors of InhA.

Objective:

In the present work, a Quantitative Structure Activity Relationship (QSAR) study was performed to identify the structural features of diphenyl ether analogues which contribute to InhA inhibitory activity in a favourable way.

Method:

Both 2D and 3D QSAR models were built and compared. Several fingerprint based 2D QSAR models were generated and their relationship with the structural features was studied. Models which corroborated the inhibitory activity of the molecules with their structural features were selected and studied in detail.

Results:

A 2D-QSAR model, with dendritic fingerprints having regression coefficient, for test set molecules Q2 =0.8132 and for the training set molecules, R2 =0.9607 was obtained. Additionally, an atom-based 3D-QSAR model with Q2 =0.7697 and R2 =0.9159 was also constructed.

Conclusion:

The data reported by various models provides guidance for the designing of structurally new diphenyl ether inhibitors with potential activity against InhA of M. tuberculosis.

Aqueous Molecular Dynamics Simulations of the M. tuberculosis Enoyl-ACP Reductase-NADH System and Its Complex with a Substrate Mimic or Diphenyl Ethers Inhibitors

Molecular dynamics (MD) simulations of 12 aqueous systems of the NADH-dependent enoyl-ACP reductase from Mycobacterium tuberculosis (InhA) were carried out for up to 20–40 ns using the GROMACS 4.5 package. Simulations of the holoenzyme, holoenzyme-substrate, and 10 holoenzyme-inhibitor complexes were conducted in order to gain more insight about the secondary structure motifs of the InhA substrate-binding pocket. We monitored the lifetime of the main intermolecular interactions: hydrogen bonds and hydrophobic contacts. Our MD simulations demonstrate the importance of evaluating the conformational changes that occur close to the active site of the enzyme-cofactor complex before and after binding of the ligand and the influence of the water molecules. Moreover, the protein-inhibitor total steric (E_LJ) and electrostatic (E_C) interaction energies, related to Gly96 and Tyr158, are able to explain 80% of the biological response variance according to the best linear equation, pK_i = 7.772 − 0.1885 × Gly96 + 0.0517 × Tyr158 (R² = 0.80; n = 10), where interactions with Gly96, mainly electrostatic, increase the biological response, while those with Tyr158 decrease. These results will help to understand the structure-activity relationships and to design new and more potent anti-TB drugs.