Breast cancer (BC) ranks as the most prevalent malignancy affecting women globally, with apoptosis playing a pivotal role in its pathological progression. Despite the crucial role of apoptosis in BC development, there is limited research exploring the relationship between BC prognosis and apoptosis-related genes (ARGs). Therefore, this study aimed to establish a BC-specific risk model centered on apoptosis-related factors, presenting a novel approach for predicting prognosis and immune responses in BC patients. Utilizing data from The Cancer Gene Atlas (TCGA), Cox regression analysis was employed to identify differentially prognostic ARGs and construct prognostic models. The accuracy and clinical relevance of the model, along with its efficacy in predicting immunotherapy outcomes, were evaluated using independent datasets, Receiver Operator Characteristic (ROC) curves, and nomogram. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses were used to predict potential mechanical pathways. The CellMiner database is used to assess drug sensitivity of model genes. A survival risk model comprising eight prognostically relevant apoptotic genes (PMAIP1, TP53AIP1, TUBA3D, TUBA1C, BCL2A1, EMP1, GSN, F2) was established based on BC patient samples from TCGA. Calibration curves validated the ROC curve and nomogram, demonstrating excellent accuracy and clinical utility. In samples from the Gene Expression Omnibus (GEO) datasets and immunotherapy groups, the low-risk group (LRG) demonstrated enhanced immune cell infiltration and improved immunotherapy responses. Model genes also displayed positive associations with sensitivity to multiple drugs, including vemurafenib, dabrafenib, PD-98059, and palbociclib. This study successfully developed and validated a prognostic model based on ARGs, offering new insights into prognosis and immune response prediction in BC patients. These findings hold promise as valuable references for future research endeavors in this field.