Molecular description of pyrimidine-based inhibitors with activity against FAK combining 3D-QSAR analysis, molecular docking and molecular dynamics

Focal adhesion kinase (FAK) is a promising target for developing more effective anticancer drugs. To better understand the structure-activity relationships and mechanism of actions of FAK inhibitors, a molecular modeling study using 3D-QSAR, molecular docking, molecular dynamics simulations, and binding free energy analysis were conducted. Two types of satisfactory 3D-QSAR models were generated, comprising the CoMFA model (R2cv = 0.528, R2pred = 0.7557) and CoMSIA model (R2cv = 0.757, R2pred = 0.8362), for predicting the inhibitory activities of novel inhibitors. The derived contour maps indicate structural characteristics for substituents on the template. Molecular docking, molecular dynamic simulations and binding free energy calculations further reveal that the binding of inhibitors to FAK is mainly contributed from hydrophobic, electrostatic and hydrogen bonding interactions. In addition, some key residues (Arg14, Glu88, Cys90, Arg138, Asn139, Leu141, and Leu155) responsible for ligand-receptor binding are highlighted. All structural information obtained from 3D-QSAR models and molecular dynamics is consist with the available experimental activities. All the results will facilitate the optimization of this series of FAK inhibitors with higher inhibitory activities.