3D-QSAR, docking, molecular dynamics simulation and free energy calculation studies of some pyrimidine derivatives as novel JAK3 inhibitors

Janus kinase 3 (JAK3) is a promising drug target for the treatment of inflammatory diseases, autoimmune disorders, organ transplant rejection and various cancers. In the present study, 3D-QSAR, docking, MD simulation and MM/PBSA studies were performed on a series of pyrimidine-based JAK3 inhibitors. A reliable COMSIA (q2 = 0.717 and r2 = 0.986) model was developed and validated using external validation test set, bootstrapping, progressive scrambling and rm2 metrics analyses. Structural requirements identified through contour maps of the model were strategically utilized to computationally design 170 novel JAK3 inhibitors with improved potency. Docking studies were performed on the selected data set and newly designed compounds to show their binding mode and to identify important interacting residues inside the active site of JAK3. In addition, docking results of the selected designed compounds inside the active sites of JAK1, JAK2 and TYK2 indicated their JAK3 selectivity. MD simulation (100 ns) on the docked complex of compound 28 (one of highly active compounds of the data set) assisted in the further exploration of the binding interactions. Some crucial residues like Lys830 (glycine-rich loop), Val836, Ala853, Leu905 (hinge region), Cys909, Asn954, Leu956 and Ala966 were identified. Hydrogen bond interactions with hinge residue Leu905 were critical for the binding of JAK3 inhibitors. Additionally, MM/PBSA calculation provided the binding free energy of the compound 28. Newly designed molecules showed promising results in the preliminary in silico ADMET evaluations. Outcomes of the study can further be exploited to develop potent JAK3 inhibitors.