Identification of novel CDK 9 inhibitors based on virtual screening, molecular dynamics simulation, and biological evaluation

Cyclin-dependent kinase 9 (CDK9) is a member of the CDK subfamily and plays a major role in the regulation of transcriptional elongation. It has attracted widespread attention as a therapeutic target for cancer. Here, we aimed to explore novel CDK 9 inhibitors by using a hybrid virtual screening strategy. A hybrid virtual screening strategy was constructed with computer-aided drug design (CADD). First, compounds were filtered in accordance with Lipinski's rule of five and adsorption, distribution, metabolism, excretion, and toxicity (ADMET) properties. Second, a 3D-QSAR pharmacophore model was built and used as a 3D query to screen the obtained hit compounds. Third, the hit compounds were subjected to molecular docking studies. Fourth, molecular dynamics (MD) simulations were performed on CDK9 in complex with the final hits to examine the structural stability. Finally, CDK9 kinase biochemical assay was performed to identify the biological activity of the hit compounds. Seven hit compounds were screened out. These hit compounds showed drug-like properties in accordance with Lipinski's rule of five and ADMET. Complexes involving the six hit compounds bound to CDK9 exhibited good structural stability in the MD simulation. Furthermore, these six hit compounds had strong inhibitory activity against CDK9 kinase. In particular, hit 3 showed the most promising activity with the percentage of 71%. The six hit compounds may be promising novel CDK9 inhibitors, and the hybrid virtual screening strategy designed in this study provides an important reference for the design and synthesis of novel CDK9 inhibitors.