3D-QSAR and Relative Binding Affinity Estimation of Focal Adhesion Kinase Inhibitors

Precise binding affinity predictions are essential for structure-based drug discovery (SBDD). Focal adhesion kinase (FAK) is a member of the tyrosine kinase protein family and is overexpressed in a variety of human malignancies. Inhibition of FAK using small molecules is a promising therapeutic option for several types of cancer. Here, we conducted computational modeling of FAK targeting inhibitors using 3-dimensional structure-activity relationship (3D-QSAR), molecular dynamics (MD), and hybrid topology-based free energy perturbation (FEP) methods. The structure-activity relationship (SAR) studies between the physicochemical descriptors and inhibitory activities of the chemical compounds were performed with reasonable statistical accuracy using CoMFA and CoMSIA. These are two well-known 3D-QSAR methods based on the principle of supervised machine learning (ML). Essential information regarding residue-specific binding interactions was determined using the MD and MM-PB/GBSA methods. Finally, physics-based relative binding free energy (〖∆∆G〗_RBFE^(A→B)) values of analogous ligands were estimated using the alchemical FEP simulation. An acceptable agreement was observed between the experimental and computed relative binding free energies. The overall results using ML and physics-based hybrid approaches could be useful for the rational optimization of accessible lead compounds with similar scaffolds targeting the FAK receptor.