Ligand- and structure-based identification of novel CDK9 inhibitors for the potential treatment of leukemia

Cyclin-dependent kinase 9 (CDK9) plays a vital role in controlling cell transcription and has been an attractive target for cancer treatment. Herein, ten predictive models derived from 1330 unique molecules against CDK9 were constructed based on molecular fingerprints and graphs using two conventional machine learning and four deep learning methods. The evaluation results showed that FP-GNN deep learning architecture performed best for CDK9 inhibitors prediction with the highest BA and F1 values of 0.681 and 0.912 for testing set. We then performed virtual screening to identify new CDK9 inhibitors by incorporating the optimal established predictive model and molecular docking. Five compounds were identified to show broad anticancer activity against various cancer cell lines through bioassays. For example, C9 exhibited antiproliferative activities against HeLa, MOLM-13 and MDA-MB-231 with IC50 values of 2.53, 3.92 and 11.65 μM. Kinase inhibition assay results demonstrated that these compounds displayed submicromolar (214 ∼ 504 nM) inhibitory activities against CDK9. Further cellular mechanism evaluation revealed that C9 suppressed the activity of CDK9 and interfered with the expression of Mcl-1 and cleaved PARP in MOLM-13 cells, resulting in the induction of cellular apoptosis. In addition, C9 displayed a good stability in rat liver microsomes, artificial gastrointestinal fluid and plasm. An online platform (called DEEPCDK9Pred) was developed based on the FP-GNN models to predict or design new CDK9 inhibitors. Collectively, our findings demonstrated that FP-GNN algorithm can achieve accurate prediction of CDK9 inhibitors and the subsequent discovery of C9 as a new potential CDK9 inhibitor deserves further structural modification for the treatment of leukemia.