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

Abstract

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, costeffective, 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 is associated with the Sphere Exclusion method (PLS-SE), having r² and q² 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 q²= 0.8509 and a non-crossvalidated correlation coefficient pred_r²= 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.

Keywords: 2D; 3D QSAR; ADME-T.; Molecular Docking; Pharmacophore Modelling; Reverse transcriptase inhibitors; anti-HIV.