Identification of EGFR inhibitors as potential agents for cancer therapy: pharmacophore-based modeling, molecular docking, and molecular dynamics investigations

Context As a member of a large family of proteins that together regulate various aspects of cell growth and development, the epidermal growth factor receptor (EGFR) is a validated target for the development of new drugs. Herein, we compiled a library of 62 compounds from the PubChem database with similar pharmacophores as osimertinib, which to our knowledge represents the only drug capable of overcoming EGFR-T790M-mutated NSCLC until date. Subsequently, we launched a docking-based virtual screening campaign against the EGFR kinase with the compiled chemical entities. The virtual screen identified 3 hit candidates (CID_126667097, CID_137660592, and CID_137659061) with lower binding energy/higher affinity (− 8.7 kcal/mol, − 8.6 kcal/mol, and − 8.5 kcal/mol, respectively) than the standard osimertinib (− 8.4 kcal/mol). Molecular dynamics metrics such as RMSD, RMSF, ROG, and intermolecular H-bond were used to substantiate the stability of the promising drug candidates at the binding pocket of EGFR after 100,000 ps production run. Overall, our molecular modeling study portrayed CID_126667097, CID_137660592, and CID_137659061 as lead-like drug candidates that may be further developed for the treatment of EGFR-associated cancer disease. Methods Molecular docking was conducted with Autodock Vina. A total of 62 compounds were compiled for the docking screen, which were then downloaded in SMILE format and converted to Protein Data Bank (PDB) format using the Openbabel online server. Finally, Gromacs 2022.3 was used to perform MD simulation to substantiate the stability of the hit candidates.