2D-QSAR and docking study of a series of coumarin derivatives as inhibitors of CDK (anticancer activity) with an application of the molecular docking method

Quantitative Structure Activity Relationship (QSAR) analysis techniques are tools largely utilized in many research fields, including drug discovery processes. In this work electronic descriptors are calculated with the Gaussian 03W software using the DFT method with the BecKe 3-parameters exchange functional and Lee-Yang-Parr correlation functional, with Kohn and Sham orbitals (KS) developed on a Gaussian Basis of type 6-31G (d), in combination with five Lipinski parameters that have been calculated with ChemOffice software, in order to develop a statistically verified 2D-QSAR model able to predict the biological activity of new molecules belonging to the same range of coumarins rather than chemical synthesis and biological evaluations that require more time and resources. Two QSAR models against both MCF-7 and HepG-2 cell lines are obtained using the multiple linear regression method. The predictive power of these models has been confirmed by internal and external validation. The Leverage method was used to determine the domain of applicability of the 2D-QSAR models developed. The results indicate that the best QSAR model is the one that links the 2D descriptors with the CDK inhibitory activity of the cell line (HepG-2) R(2) = 0.748, R(2)cv = 0.618, MSE = 0.03 for the learning series and R(2) = 0.73, MSE = 0.18 for the test series. This model implies that coumarin inhibitory activity is strongly related to dipole moment and the number of hydrogen bond donors. The results obtained suggest the importance of studying structure-activity relationships as a principal axis in drug design. The docking procedure using AutoDOCK Tools was also used to understand the mechanisms of molecular interactions and consequently, to develop new inhibitors.