Pharmacophore modeling of JAK1: A target infested with activity-cliffs

Janus kinase 1 (JAK1) is protein kinase involved in autoimmune diseases (AIDs). JAK1 inhibitors have shown promising results in treating AIDs. JAK1 inhibitors are known to exhibit regions of SAR discontinuity or activity cliffs (ACs). ACs represent fundamental challenge to successful QSAR/pharmacophore modeling because QSAR modeling rely on the basic premise that activity is a smooth continuous function of structure. We propose that ACs exist because active ACs members exhibit subtle, albeit critical, enthalpic features absent from their inactive twins. In this context we compared the performances of two computational modeling workflows in extracting valid pharmacophores from 151 diverse JAK1 inhibitors that include ACs: QSAR-guided pharmacophore selection versus docking-based comparative intermolecular contacts analysis (db-CICA). The two methods were judged based on the receiver operating characteristic (ROC) curves of their corresponding pharmacophore models and their abilities to distinguish active members among established JAK1 ACs. db-CICA modeling significantly outperformed ligand-based pharmacophore modeling. The resulting optimal db-CICA pharmacophore was used as virtual search query to scan the National Cancer Institute (NCI) database for novel JAK1 inhibitory leads. The most active hit showed IC₅₀ of 1.04 μM. This study proposes the use of db-CICA modeling as means to extract valid pharmacophores from SAR data infested with ACs.