Abstract Background Neuroblastoma (NB), a prevalent pediatric solid tumor, presents formidable challenges due to its high malignancy and intricate pathogenesis. The role of disulfidptosis, a novel form of programmed cell death, remains poorly understood in the context of NB. Methods Gaussian mixture model (GMM)-identified disulfidptosis-related molecular subtypes in NB, differential gene analysis, survival analysis, and gene set variation analysis were conducted subsequently. Weighted gene co-expression network analysis (WGCNA) selected modular genes most relevant to the disulfidptosis core pathways. Integration of machine learning approaches revealed the combination of the Least absolute shrinkage and selection operator (LASSO) and Random Survival Forest (RSF) provided optimal dimensionality reduction of the modular genes. The resulting model was validated, and a nomogram assessed disulfidptosis characteristics in NB. Core genes were filtered and subjected to tumor phenotype and disulfidptosis-related experiments. Results GMM clustering revealed three distinct subtypes with diverse prognoses, showing significant variations in glucose metabolism, cytoskeletal structure, and tumor-related pathways. WGCNA highlighted the red module of genes highly correlated with disulfide isomerase activity, cytoskeleton formation, and glucose metabolism. The LASSO and RSF combination yielded the most accurate and stable prognostic model, with a significantly worse prognosis for high-scoring patients. Cytological experiments targeting core genes (CYFIP1, EMILIN1) revealed decreased cell proliferation, migration, invasion abilities, and evident cytoskeletal deformation upon core gene knockdown. Conclusions This study showcases the utility of disulfidptosis-related gene scores for predicting prognosis and molecular subtypes of NB. The identified core genes, CYFIP1 and EMILIN1, hold promise as potential therapeutic targets and diagnostic markers for NB.