BackgroundProstate cancer (PCa) is one of the most prevalent cancers among males, and its mortality rate is increasing due to biochemical recurrence (BCR). Glycolysis has been proven to play an important regulatory role in tumorigenesis. Although several key regulators or predictors involved in PCa progression have been found, the relationship between glycolysis and PCa is unclear; we aimed to develop a novel glycolysis-associated multifactor prediction model for better predicting the prognosis of PCa.MethodsDifferential mRNA expression profiles derived from the Cancer Genome Atlas (TCGA) PCa cohort were generated through the “edgeR” package. Glycolysis-related genes were obtained from the GSEA database. Univariate Cox and LASSO regression analyses were used to identify genes significantly associated with disease-free survival. ROC curves were applied to evaluate the predictive value of the model. An external dataset derived from Gene Expression Omnibus (GEO) was used to verify the predictive ability. Glucose consumption and lactic production assays were used to assess changes in metabolic capacity, and Transwell assays were used to assess the invasion and migration of PC3 cells.ResultsFive glycolysis-related genes were applied to construct a risk score prediction model. Patients with PCa derived from TCGA and GEO (GSE70770) were divided into high-risk and low-risk groups according to the median. In the TCGA cohort, the high-risk group had a poorer prognosis than the low-risk group, and the results were further verified in the GSE70770 cohort. In vitro experiments demonstrated that knocking down HMMR, KIF20A, PGM2L1, and ANKZF1 separately led to less glucose consumption, less lactic production, and inhibition of cell migration and invasion, and the results were the opposite with GPR87 knockdown.ConclusionThe risk score based on five glycolysis-related genes may serve as an accurate prognostic marker for PCa patients with BCR.