Overexpression of polo-like kinase 1 (PLK1) has been found in many different types of cancers. With its essential role in cell proliferation, PLK1 has been determined to be a broad-spectrum anti-cancer target. In this study, 3D-QSAR, molecular docking, and molecular dynamics (MD) simulations were applied on a series of novel pteridinone derivatives as PLK1 inhibitors to discover anti-cancer drug candidates. In this work, three models—CoMFA (Q² = 0.67, R² = 0.992), CoMSIA/SHE (Q² = 0.69, R² = 0.974), and CoMSIA/SEAH (Q² = 0.66, R² = 0.975)—of pteridinone derivatives were established. The three models that were established gave Rpred2 = 0.683, Rpred 2= 0.758, and Rpred 2= 0.767, respectively. Thus, the predictive abilities of the three proposed models were successfully evaluated. The relations between the different champs and activities were well-demonstrated by the contour chart of the CoMFA and CoMSIA/SEAH models. The results of molecular docking indicated that residues R136, R57, Y133, L69, L82, and Y139 were the active sites of the PLK1 protein (PDB code: 2RKU), in which the more active ligands can inhibit the enzyme of PLK1. The results of the molecular dynamic MD simulation diagram were obtained to reinforce the previous molecular docking results, which showed that both inhibitors remained stable in the active sites of the PLK1 protein (PDB code: 2RKU) for 50 ns. Finally, a check of the ADME-Tox properties of the two most active molecules showed that molecular N° 28 could represent a good drug candidate for the therapy of prostate cancer diseases.
Keywords: 3D-QSAR; PLK1 inhibitors; anti-cancer; dynamic simulation; molecular docking.