Background/Objectives: Despite current treatments extending the lifespan of Glioblastoma (GBM) patients, the average survival time is around 15–18 months, underscoring the fatality of GBM. This study aims to investigate the impact of sample heterogeneity on gene expression in GBM, identify key metabolic pathways and gene modules, and explore potential therapeutic targets. Methods: In this study, we analysed GBM transcriptome data derived from The Cancer Genome Atlas (TCGA) using genome-scale metabolic models (GEMs) and co-expression networks. We examine transcriptome data incorporating tumour purity scores (TPSs), allowing us to assess the impact of sample heterogeneity on gene expression profiles. We analysed the metabolic profile of GBM by generating condition-specific GEMs based on the TPS group. Results: Our findings revealed that over 90% of genes showing brain and glioma specificity in RNA expression demonstrate a high positive correlation, underscoring their expression is dominated by glioma cells. Conversely, negatively correlated genes are strongly associated with immune responses, indicating a complex interaction between glioma and immune pathways and non-tumorigenic cell dominance on gene expression. TPS-based metabolic profile analysis was supported by reporter metabolite analysis, highlighting several metabolic pathways, including arachidonic acid, kynurenine and NAD pathway. Through co-expression network analysis, we identified modules that significantly overlap with TPS-correlated genes. Notably, SOX11 and GSX1 are upregulated in High TPS, show a high correlation with TPS, and emerged as promising therapeutic targets. Additionally, NCAM1 exhibits a high centrality score within the co-expression module, which shows a positive correlation with TPS. Moreover, LILRB4, an immune-related gene expressed in the brain, showed a negative correlation and upregulated in Low TPS, highlighting the importance of modulating immune responses in the GBM mechanism. Conclusions: Our study uncovers sample heterogeneity’s impact on gene expression and the molecular mechanisms driving GBM, and it identifies potential therapeutic targets for developing effective treatments for GBM patients.