This research explored the immune characteristics of natural killer (NK) cells in lung adenocarcinoma (LUAD) and their predictive role on patient survival and immunotherapy response. Molecular subtyping of LUAD samples was performed by evaluating NK cell-associated pathways and genes in The Cancer Genome Atlas (TCGA) dataset using consistent clustering. 12 programmed cell death (PCD) patterns were acquired from previous study. Riskscore prognostic models were constructed using Least absolute shrinkage and selection operator (Lasso) and Cox regression. The model stability was validated in Gene Expression Omnibus database (GEO). We classified LUAD into three different molecular subgroups based on NK cell-related genes, with the worst prognosis in C1 patients and the optimal in C3. Homologous Recombination Defects, purity and ploidy, TMB, LOH, Aneuploidy Score, were the most high-expressed in C1 and the least expressed in C3. ImmuneScore was the highest in C3 type, suggesting greater immune infiltration in C3 subtype. C1 subtypes had higher TIDE scores, indicating that C1 subtypes may benefit less from immunotherapy. Generally, C3 subtype presented highest PCD patterns scores. With four genes, ANLN, FAM83A, RHOV and PARP15, we constructed a LUAD risk prediction model with significant differences in immune cell composition, cell cycle related pathways between the two risk groups. Samples in C1 and high group were more sensitive to chemotherapy drug. The score of PCD were differences in high- and low-groups. Finally, we combined Riskscore and clinical features to improve the performance of the prediction model, and the calibration curve and decision curve verified that the great robustness of the model. We identified three stable molecular subtypes of LUAD and constructed a prognostic model based on NK cell-related genes, maybe have a greater potential for application in predicting immunotherapy response and patient prognosis.