Combined multi-pharmacophore, molecular docking and molecular dynamic study for discovery of promising MTH1 inhibitors

A Combined approach of Pharmacophore Modeling, QSAR Study, Molecular Docking and in silico ADME/Tox prediction of 4-Arylthio & 4- Aryloxy-3- Iodopyridine-2(1H)-one analogs to identify potential Reverse Transcriptase inhibitor: Anti-HIV agents

BACKGROUND

Reverse transcriptase is an important therapeutic target to treat AIDS caused by the Human Immunodeficiency Virus (HIV). Despite many effective anti-HIV drugs, reverse transcriptase (RT) inhibitors remain the cornerstone of the drug regimen to treat AIDS. In the present work, we have expedited the use of different computational modules and presented an easy, cost-effective and high throughput screening method to identify potential reverse transcriptase inhibitors.

METHODS

A congeneric series of 4-Arylthio & 4-Aryloxy-3- Iodopyridine-2(1H)-one analogs having anti-HIV activity were subjected to structure-based 2D, 3D QSAR, Pharmacophore Modeling, and Molecular Docking to elucidate the structural properties required for the design of potent HIV-RT inhibitors. Prediction of preliminary Pharmacokinetic and the Drug Likeliness profile was performed for these compounds by in silico ADME study.

RESULTS

The 2D and 3D- QSAR models were developed by correlating two and three-dimensional descriptors with activity (pIC50) by sphere exclusion method and k-nearest neighbor molecular field analysis approach, respectively. The significant 2D- QSAR model developed by Partial Least Square associated with the Sphere Exclusion method (PLS-SE) having r2 and q2 values 0.9509 and 0.8038 respectively. The 3D-QSAR model by Step Wise variable selection method (SW-kNN MFA) is more significant which has a cross-validated squared correlation coefficient q2= 0.8509 and a non-cross-validated correlation coefficient pred_r2= 0.8102. The pharmacophore hypothesis was developed which comprised 5 features includes 3 aliphatic regions (Ala), 1 H-bond donor (HDr) and 1 H-bond acceptor (HAc). Docking studies of the selected inhibitors with the active site of reverse transcriptase enzyme showed hydrogen bond and π - π interaction with LYS-101, LYS-103, TYR- 181, TYR-188 and TRP-229 residues present at the active site. All the candidates with good bioavailability and ADMET drug likeliness properties.

CONCLUSION

The results of the present work provide more useful information and important structural insights for the discovery, design of novel and potent reverse transcriptase inhibitors with high therapeutic windows in the future.