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

Fangfang Wang,Wei Yang,Ran Li,Zhihai Sui,Guijuan Cheng,Bo Zhou
DOI: https://doi.org/10.1016/j.arabjc.2021.103144
IF: 6.212
2021-06-01
Arabian Journal of Chemistry
Abstract: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 (R<sup>2</sup><sub>cv</sub>=0.528, R<sup>2</sup><sub>pred</sub>=0.7557) and CoMSIA model (R<sup>2</sup><sub>cv</sub>=0.757, R<sup>2</sup><sub>pred</sub>=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.
chemistry, multidisciplinary
What problem does this paper attempt to address?
This paper aims to solve the problem of how to design more effective focal adhesion kinase (FAK) inhibitors through molecular modeling methods. Specifically, the researchers used methods such as three - dimensional quantitative structure - activity relationship analysis (3D - QSAR), molecular docking, molecular dynamics simulations, and binding free - energy analysis to better understand the structure - activity relationship of FAK inhibitors and their mechanisms of action. ### Main problems addressed in the paper 1. **Understanding the structure - activity relationship**: Through 3D - QSAR models (including CoMFA and CoMSIA models), the researchers hope to reveal the influence of different structural features on the activity of FAK inhibitors. 2. **Optimizing inhibitor design**: Through molecular docking and molecular dynamics simulations, the researchers hope to find key binding sites and interaction patterns, thereby guiding the design of new FAK inhibitors. 3. **Verifying the predictive ability of the model**: Through the verification of an external test set, evaluate the predictive ability and reliability of the established 3D - QSAR model. ### Main methods - **3D - QSAR analysis**: Generated CoMFA and CoMSIA models, which are respectively used to predict the steric field and electrostatic field effects of inhibitors. - **Molecular docking**: Determine the binding mode of inhibitors to FAK. - **Molecular dynamics simulations**: Verify the reliability of molecular docking results and provide a higher - precision structure of the binding complex. - **Binding free - energy analysis**: Evaluate the thermodynamic stability of the binding of inhibitors to FAK. ### Main findings - **3D - QSAR models**: Two satisfactory 3D - QSAR models were generated. Among them, the cross - validation correlation coefficient \( R^2_{cv} \) of the CoMFA model is 0.528, and the prediction correlation coefficient \( R^2_{\text{pred}} \) is 0.7557; the \( R^2_{cv} \) of the CoMSIA model is 0.757, and the \( R^2_{\text{pred}} \) is 0.8362. - **Contour map analysis**: Through CoMFA and CoMSIA contour maps, key structural features that affect the activity of inhibitors were revealed, such as the introduction of small groups, large groups, and electronegative groups at specific positions. - **Molecular dynamics simulations**: Verified the reliability of molecular docking results and found the important roles of key residues (such as Arg14, Glu88, Cys90, Arg138, Asn139, Leu141, Leu155) in ligand - receptor binding. These results will be helpful for optimizing the design of this series of FAK inhibitors and improving their inhibitory activity.