The treatment of patients with nonsmall cell lung cancer (NSCLC) using epidermal growth factor receptor (EGFR) inhibitors is complicated by drug-sensitive activating L858R/T790M and L858R/T790M/C797S mutations. To overcome drug resistance, a series of furopyridine (PD) compounds were virtually screened to identify potent EGFR inhibitors using molecular docking and molecular dynamics (MD) simulations based on the solvated interaction energy (SIE) method. Several PD compounds identified from virtual screening demonstrated the potential to suppress both wild-type and mutant forms of EGFR, with IC50 values in the nanomolar range. Among these, PD18 and PD56 exhibited highly potent inhibitory activity against both wild-type and mutant forms of EGFR, surpassing the efficacy of known drugs. Additionally, both PD compounds were cytotoxic to NSCLC cell lines (A549 and H1975) while being nontoxic to normal cell lines (Vero). The interaction mechanisms of both PD compounds complexed with wild-type and mutant forms of EGFR were elucidated through 500 ns molecular dynamics simulations. The predicted binding affinity from molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) correlated well with the experimental binding affinity derived from IC50 values. Furthermore, it was observed that van der Waals interactions, rather than electrostatic interactions, played a significant role in interacting with EGFR's active site. The strong inhibitory activity against EGFR was attributed to two key residues, M793 and S797, via hydrogen bonding, corresponding with lower solvent accessibility and a higher number of atomic contacts. Therefore, these potent compounds could be developed as promising drugs targeting both wild-type and mutant EGFR for the treatment of NSCLC.
Read full abstract