Integrated computational approaches for designing potent pyrimidine-based CDK9 inhibitors: 3D-QSAR, docking, and molecular dynamics simulations

Abstract

CDK9 is an emerging target for the development of anticancer drugs. The development of CDK9 inhibitors with significant potency had consistently posed a formidable challenge. In the current research, a number of computational methodologies, such as, 3D-QSAR, molecular docking, fingerprint analysis, molecular dynamic (MD) simulations followed by MMGB/PBSA and ADMET studies were used systemically to uncover the binding mechanism of pyrimidine derivatives against CDK9. The CoMFA and CoMSIA models having high q² (0.53, 0.54) and r² values (0.96, 0.93) respectively indicating that model could accurately predict the bioactivities of CDK9 inhibitors. Using the R-group exploration technique implemented by the Spark™ by Cresset group, the structural requirements revealed by the contour maps of model were utilized strategically to create an in-house library of 100 new CDK9 inhibitors. Additionally, the compounds from the in-house library were mapped into 3D-QSAR model which predicted pIC₅₀ values comparable to the experimental values. A comparison between 3D-QSAR generated contours and molecular docking conformation of ligands was performed to elucidate the essentials of CDK9 inhibitor design. MD simulations (100 ns) were performed on the selected docked complexes A21, A14 and D98 which contributed in validating the binding interactions. According to the findings of binding free energy analysis (MMGB/PBSA), It was observed that residues CYS106 and GLU107 had a considerable tendency to facilitate ligand-protein interactions via H-bond interactions. The aforementioned findings have the potential to enhance researchers comprehension of the mechanism underlying CDK9 inhibition and may be utilized in the development of innovative and efficacious CDK9 inhibitors.

Keywords: Binding free energy; CDK9; CoMFA/CoMSIA modeling; Distance analysis; Residual decomposition.