Predictive computational modeling to define effective treatment strategies for bone metastatic prostate cancer

The ability to rapidly assess the efficacy of therapeutic strategies for incurable bone metastatic prostate cancer is an urgent need. Pre-clinical in vivo models are limited in their ability to define the temporal effects of therapies on simultaneous multicellular interactions in the cancer-bone microenvironment. Integrating biological and computational modeling approaches can overcome this limitation. Here, we generated a biologically driven discrete hybrid cellular automaton (HCA) model of bone metastatic prostate cancer to identify the optimal therapeutic window for putative targeted therapies. As proof of principle, we focused on TGFβ because of its known pleiotropic cellular effects. HCA simulations predict an optimal effect for TGFβ inhibition in a pre-metastatic setting with quantitative outputs indicating a significant impact on prostate cancer cell viability, osteoclast formation and osteoblast differentiation. In silico predictions were validated in vivo with models of bone metastatic prostate cancer (PAIII and C4-2B). Analysis of human bone metastatic prostate cancer specimens reveals heterogeneous cancer cell use of TGFβ. Patient specific information was seeded into the HCA model to predict the effect of TGFβ inhibitor treatment on disease evolution. Collectively, we demonstrate how an integrated computational/biological approach can rapidly optimize the efficacy of potential targeted therapies on bone metastatic prostate cancer.