Dysregulated epidermal growth factor receptor (EGFR) expression is frequently observed in non-small cell lung cancer (NSCLC) growth and metastasis. Despite recent successes in the development of tyrosine kinase inhibitors (TKIs), inevitable resistance to TKIs has led to urgent calls for novel EGFR inhibitors. Herein, we report a rational workflow used to identify novel EGFR-TKIs by combining hybrid ligand- and structure-based pharmacophore models. Three types of models were developed in this workflow, including 3D QSAR-, common feature-, and structure-based EGFR-TK domain-containing pharmacophores. A National Cancer Institute (NCI) compound dataset was adopted for multiple-stage pharmacophore-based virtual screening (PBVS) of various pharmacophore models. The six top-scoring compounds were identified through the PBVS pipeline coupled with molecular docking. Among these compounds, NSC609077 exerted a significant inhibitory effect on EGFR activity in gefitinib-resistant H1975 cells, as determined by an enzyme-linked immunosorbent assay (ELISA). Further investigations showed that NSC609077 inhibited the anchorage-dependent growth and migration of lung cancer cells. Furthermore, NSC609077 exerted a suppressive effect on the EGFR/PI3K/AKT pathway in H1975 cells. In conclusion, these findings suggest that hybrid virtual screening may accelerate the development of targeted drugs for lung cancer treatment.