Abstract High heterogeneity of lung adenocarcinoma (LUAD) is a major clinical challenge. This study aims to characterize the molecular features of LUAD through classification based on metabolism-related genes. A total of 500 LUAD samples from The Cancer Genome Atlas (TCGA) and 612 from Gene Expression Omnibus (GEO) were integrated with 2,753 metabolism-related genes to determine the molecular classification. Systematic bioinformatics analysis was used to conduct correlation analysis between metabolism-related classification and molecular characteristics of LUAD. LUAD patients were divided into three molecular clusters (C1-C3). Survival analysis revealed that C1 and C2 showed good and poor prognoses, respectively. Associational analysis of classification and molecular characteristics revealed that C1 was associated with low pathological stage, metabolic pathways, high metabolic process, active immune process and checkpoint, sensitive drug response, as well as a low genetic mutation. Nevertheless, C2 was associated with high pathological stage, carcinogenic pathways, low metabolic process, inactive immune signatures, resistant drug response, and frequent genetic mutation. Eventually, a classifier with 60 metabolic genes was constructed, confirming the robustness of molecular classification on LUAD. Our findings promote the understanding of LUAD molecular characteristics, and the research data may be used for providing information be helpful for clinical diagnosis and treatment.