Prostate cancer (PCa) impacts over 180,000 men every year in the USA alone, with 26,000 patients expected to succumb to the disease ( cancer.gov ). The primary cause of death is metastasis, with secondary lesions most commonly occurring in the skeleton. Prostate cancer to bone metastasis is an important, yet poorly understood, process that is difficult to explore with experimental techniques alone. To this end we have utilized a hybrid (discrete-continuum) cellular automaton model of normal bone matrix homeostasis that allowed us to investigate how metastatic PCa can disrupt the bone microenvironment. Our previously published results showed that PCa cells can recruit mesenchymal stem cells (MSCs) that give rise to bone-building osteoblasts. MSCs are also thought to be complicit in the establishment of successful bone metastases (Lu, in Mol Cancer Res 4(4):221-233, 2006). Here we have explored the aspects of early metastatic colonization and shown that the size of PCa clusters needs to be within a specific range to become successfully established: sufficiently large to maximize success, but not too large to risk failure through competition among cancer and stromal cells for scarce resources. Furthermore, we show that MSC recruitment can promote the establishment of a metastasis and compensate for relatively low numbers of PCa cells seeding the bone microenvironment. Combined, our results highlight the utility of biologically driven computational models that capture the complex and dynamic dialogue between cells during the initiation of active metastases.