# Structure-based identification of novel FAK1 inhibitors using pharmacophore modeling, molecular dynamics, and MM/PBSA calculations

**Authors:** Amirhossein Hajipasha, Nilofar Ghaffari Cherati, Mohammad Darzi, Seyedeh Sana Nateghi, Seyed Arshia Sadat Mohsenian, Mahla Noorzaei, Fatemeh Kashiri Dinaki, Mohammad Halimi

PMC · DOI: 10.1038/s41598-025-23203-8 · Scientific Reports · 2025-11-11

## TL;DR

This paper uses computational methods to identify new inhibitors for FAK1, a protein linked to cancer progression, and highlights a promising candidate for further study.

## Contribution

The study introduces a novel computational workflow combining pharmacophore modeling, molecular dynamics, and MM/PBSA to discover FAK1 inhibitors.

## Key findings

- Pharmacophore modeling identified potential FAK1 inhibitors from the ZINC database.
- Molecular dynamics simulations confirmed the stability of four compounds, with ZINC23845603 showing the strongest binding.
- ZINC23845603 exhibited favorable pharmacokinetic properties and binding energy, making it a top candidate for experimental validation.

## Abstract

Focal adhesion kinase 1 (FAK1) is a non-receptor tyrosine kinase involved in cancer metastasis and tumor progression. Due to its role in regulating cell migration and survival, FAK1 is considered a promising target for cancer therapy. While large chemical databases such as ZINC offer a wide variety of molecules, many of them have not yet been explored as potential FAK1 inhibitors.In this study, we applied several computational methods to identify new inhibitors of the FAK1 kinase domain. Pharmacophore models were built based on the FAK1–P4N complex (PDB ID: 6YOJ), and the most statistically reliable model was used to screen compounds from the ZINC database. Hits were first docked using AutoDock Vina in PyRx, and seventeen compounds with acceptable pharmacokinetic properties and low predicted toxicity were selected for more precise docking via SwissDock. Four promising candidates—ZINC23845603, ZINC44851809, ZINC266691666, and ZINC20267780—were chosen for molecular dynamics (MD) simulations using GROMACS. The stability and behavior of each protein–ligand complex were examined, and binding free energies were calculated using the MM/PBSA method. Among them, ZINC23845603 showed strong binding and interaction features similar to the known ligand P4N. Given its favorable binding energy and pharmacokinetic profile, ZINC23845603 may be a good candidate for further experimental studies targeting FAK1.

The online version contains supplementary material available at 10.1038/s41598-025-23203-8.

## Linked entities

- **Proteins:** PTK2 (protein tyrosine kinase 2)
- **Chemicals:** P4N (PubChem CID 11612883)
- **Diseases:** cancer (MONDO:0004992)

## Full-text entities

- **Genes:** PTK2 (protein tyrosine kinase 2) [NCBI Gene 5747] {aka FADK, FADK 1, FAK, FAK1, FRNK, PPP1R71}, TXK (TXK tyrosine kinase) [NCBI Gene 7294] {aka BTKL, PSCTK5, PTK4, RLK, TKL}
- **Diseases:** cancer (MESH:D009369), toxicity (MESH:D064420)
- **Chemicals:** ZINC (MESH:D015032), P4N (-)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12606203/full.md

## Figures

16 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12606203/full.md

## References

7 references — full list in the complete paper: https://tomesphere.com/paper/PMC12606203/full.md

---
Source: https://tomesphere.com/paper/PMC12606203