Cancer stem cells (CSCs) possess capacity to both self-renew and generate all cells within a tumor, and are thought to drive tumor recurrence. Targeting the stem cell niche to eradicate CSCs represents an important area of therapeutic development. The complex nature of many interacting elements of the stem cell niche, including both intracellular signals and microenvironmental growth factors and cytokines, creates a challenge in choosing which elements to target, alone or in combination. Stochastic stimulation techniques allow for the careful study of complex systems in biology and medicine and are ideal for the investigation of strategies aimed at CSC eradication. We present a mathematical model of the breast cancer stem cell (BCSC) niche to predict population dynamics during carcinogenesis and in response to treatment. Using data from cell line and mouse xenograft experiments, we estimate rates of interconversion between mesenchymal and epithelial states in BCSCs and find that EMT/MET transitions occur frequently. We examine bulk tumor growth dynamics in response to alterations in the rate of symmetric self-renewal of BCSCs and find that small changes in BCSC behavior can give rise to the Gompertzian growth pattern observed in breast tumors. Finally, we examine stochastic reaction kinetic simulations in which elements of the breast cancer stem cell niche are inhibited individually and in combination. We find that slowing self-renewal and disrupting the positive feedback loop between IL-6, Stat3 activation, and NF-[kappa]B signaling by simultaneous inhibition of IL-6 and HER2 is the most effective combination to eliminate both mesenchymal and epithelial populations of BCSCs. Predictions from our model and simulations show excellent agreement with experimental data showing the efficacy of combined HER2 and Il-6 blockade in reducing BCSC populations. Our findings will be directly examined in a planned clinical trial of combined HER2 and IL-6 targeted therapy in HER2-positive breast cancer.