Understanding the structural features of JAK2 inhibitors: a combined 3D-QSAR, DFT and molecular dynamics study

JAK2 plays a critical role in JAK/STAT signaling pathway and in patho-mechanism of myeloproliferative disorders and autoimmune diseases. Thus, effective JAK2 inhibitors provide a promising opportunity for the pharmaceutical intervention of many diseases. In this work, 3D-QSAR study was performed on a series of 1-amino-5H-pyrido-indole-4-carboxamide derivatives as JAK2 inhibitors to obtain reliable comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) models with three different alignment methods. Among the different alignment methods, ligand-based (CoMFA: q² = 0.676, r² = 0.979; CoMSIA: q² = 0.700, r² = 0.953) and pharmacophore-based alignment (CoMFA: q² = 0.710, r² = 0.982; CoMSIA: q² = 0.686, r² = 0.960) has produced better statistical results when compared to receptor-based alignment (CoMFA: q² = 0.507, r² = 0.979; CoMSIA: q² = 0.544, r² = 0.917). Statistical parameters indicated that data are well fitted and have high predictive ability. The presence of electrostatic and hydrophobic field is highly desirable for potent inhibitory activity, and the steric field plays a minor role in modulating the activity. The contour analysis indicates ARG980, ASN981, ASP939 and LEU937 have more possibility of interacting with bulky, hydrophobic groups in pyrido and positive and negative groups in pyrazole ring. Based on our findings, we have designed sixteen molecules and predicted its activity and drug-like properties. Subsequently, molecular docking, molecular dynamics and DFT calculations were performed to evaluate its potency.