Recent breakthrough therapies have improved survival rates in non-small cell lung cancer (NSCLC), but a paradigm for prospective confirmation is still lacking. Patientdatasets were mainly downloaded from TCGA, CPTAC and GEO. We conducted downstream analysis by collecting metagenes and generated 42-gene subtype classifiers to elucidate biological pathways. Subsequently, scRNA, eRNA, methylation, mutation, and copy number variation were depicted from a phenotype perspective. Enhancing the clinical translatability of molecular subtypes, preclinical models including CMAP, CCLE, and GDSC were utilized for drug repositioning. Importantly, we verified the presence of previously described three phenotypes including bronchioid, neuroendocrine, and squamoid. Poor prognosis was seen in squamoid and neuroendocrine clusters for treatment-naive and immunotherapy populations. The neuroendocrine cluster was dominated by STK11 mutations and 14q13.3 amplifications, whose related methylated loci are predictive of immunotherapy. And the greatest therapeutic potential lies in the bronchioid cluster. We further estimated the relative cell abundance of the tumor microenvironment (TME), specific cell types could be reflected among three clusters. Meanwhile, the higher portion of immune cell infiltration belonged to bronchioid and squamoid, not the neuroendocrine cluster. In drug repositioning, MEK inhibitors resisted bronchioid but were squamoid-sensitive. To conceptually validate compounds/targets, we employed RNA-seq and CCK-8/western blot assays. Our results indicated that dinaciclib and alvocidib exhibited similar activity and sensitivity in the neuroendocrine cluster. Also, a lineage factor named KLF5 recognized by inferred transcriptional factors activity could be suppressed by verteporfin.