Multiple myeloma (MM) is a malignant plasma cell tumor with high heterogeneity, characterized by anemia, hypercalcemia, renal failure, and lytic bone lesions. Although various powerful prognostic factors and models have been exploited, the development of more accurate prognosis and treatment for MM patients is still facing many challenges. Given the essential roles of super-enhancer (SE) associated genes in the tumorigenesis of MM, we tried to initially screen and identify the significant prognostic factors from SE associated genes in MM by the least absolute shrinkage and selection operator (Lasso) penalized Cox regression, univariate and multivariate Cox regression analysis using GSE24080 and GSE9782 datasets. Risk score model of five genes including CSGALNACT1, FAM53B, TAPBPL, REPIN1, and DDX11, was further constructed and the Kaplan-Meier (K-M) curves showed that the low-risk group seems to have better clinical outcome of survival compared to the high-risk group. Time-dependent receiver operating characteristic (ROC) curves presented the favorable performance of the model. An interactive nomogram consisting of the five-gene risk group and eleven clinical traits was established and identified by calibration curves. Therefore, the risk score model of SE associated five genes developed here could be used to predict the prognosis of MM patients, which may assist the clinical treatment of MM patients in the future.