This study focuses on understanding the structural features of 3,6-dihydroimidazo(4,5-d)pyrrolo(2,3-b)pyridin-2(1H)-one (dpp) derivatives to computationally identify new JAK inhibiting compounds. For the purpose, a novel virtual screening strategy, with 2D and 3D-QSAR (CoMFA and CoMSIA), data mining, pharmacophore modeling, ADMET prediction, multi-targeted protein-based docking and inverse QSAR, was employed. The 2D-QSAR equations developed for the JAK3, JAK2 and JAK1 involved five physicochemical descriptors. These descriptors correlate with the anti-RA activity with R2 values for JAK3, JAK2 and JAK1 are 0.9811, 0.8620 and 0.9740, respectively. The 3D-QSAR studies such as CoMFA and CoMSIA carried out through PLS analysis of the training set of JAK3, JAK2 and JAK1, gave Q2 values as 0.369, 0.476 and 0.490; [Formula: see text] values as 0.863, 0.684 and 0.724 and, F values as 23.098, 28.139 and 31.438, respectively. The contour maps produced by the CoMFA and CoMSIA models were used to understand the importance of hydrogen bond donor, acceptor, hydrophobic, steric and electrostatic interactions. The molecular docking studies of these selected compounds with various JAK proteins were carried out and the protein-ligand interactions were also studied. The study concluded that dpp15(s) is a highly potent JAK inhibitor with a very good predicted IC50 value.
Keywords: 2D and 3D QSAR; Janus kinase; Molecular docking; Rheumatoid arthritis; Weka.