Ligand-Based Pharmacophore Modeling, Molecular Docking, and Molecular Dynamic Studies of Dual Tyrosine Kinase Inhibitor of EGFR and VEGFR2

Epidermal growth factor receptor (EGFR) and vascular endothelial growth factor receptor 2 (VEGFR2) play an important role in cancer growth. Both of them have close relationships. Expression of EGFR will induce an angiogenic factor (VEGF) release for binding with VEGFR2. However, the existence of VEGF up-regulation independent of EGFR leads to cancer cell resistance to anti-EGFR. Therefore, a therapeutic approach targeting EGFR and VEGFR2 simultaneously may improve the outcome of cancer treatment. The present study was designed to identify potential compounds as a dual inhibitor of EGFR and VEGFR2 by the computational method. Firstly, the ligand-based pharmacophore model for each target was setup to screen of ZINC database of purchasable compounds. The hit compounds obtained by pharmacophore screening were then further screened by molecular docking studies. Taking erlotinib (EGFR inhibitor) and axitinib (VEGFR2 inhibitor) as reference drugs, six potential compounds (ZINC08398597, ZINC12047553, ZINC16525481, ZINC17418102, ZINC21942954, and ZINC38484632) were selected based on their docking scores and binding interaction. However, molecular dynamics simulations demonstrated that only ZINC16525481 and ZINC38484632 which have good binding free energy and stable hydrogen bonding interactions with EGFR and VEGFR2. The result represents a promising starting point for developing potent dual tyrosine kinases inhibitor of EGFR and VEGFR2.