How circulating tumor cluster biology contributes to the metastatic cascade: from invasion to dissemination and dormancy.

Circulating tumor cells (CTCs) are known to be prognostic for metastatic relapse and are detected in patients as solitary cells or cell clusters. Circulating tumor cell clusters (CTC clusters) have been observed clinically for decades and are of significantly higher metastatic potential compared to solitary CTCs. Recent studies suggest distinct differences in CTC cluster biology regarding invasion and survival in circulation. However, differences regarding dissemination, dormancy, and reawakening require more investigations compared to solitary CTCs. Here, we review the current state of CTC cluster research and consider their clinical significance. In addition, we discuss the concept of collective invasion by CTC clusters and molecular evidence as to how cluster survival in circulation compares to that of solitary CTCs. Molecular differences between solitary and clustered CTCs during dormancy and reawakening programs will also be discussed. We also highlight future directions to advance our current understanding of CTC cluster biology.

Integrated computational and in vivo models reveal Key Insights into macrophage behavior during bone healing.

Myeloid-derived monocyte and macrophages are key cells in the bone that contribute to remodeling and injury repair. However, their temporal polarization status and control of bone-resorbing osteoclasts and bone-forming osteoblasts responses is largely unknown. In this study, we focused on two aspects of monocyte/macrophage dynamics and polarization states over time: 1) the injury-triggered pro- and anti-inflammatory monocytes/macrophages temporal profiles, 2) the contributions of pro- versus anti-inflammatory monocytes/macrophages in coordinating healing response. Bone healing is a complex multicellular dynamic process. While traditional in vitro and in vivo experimentation may capture the behavior of select populations with high resolution, they cannot simultaneously track the behavior of multiple populations. To address this, we have used an integrated coupled ordinary differential equations (ODEs)-based framework describing multiple cellular species to in vivo bone injury data in order to identify and test various hypotheses regarding bone cell populations dynamics. Our approach allowed us to infer several biological insights including, but not limited to,: 1) anti-inflammatory macrophages are key for early osteoclast inhibition and pro-inflammatory macrophage suppression, 2) pro-inflammatory macrophages are involved in osteoclast bone resorptive activity, whereas osteoblasts promote osteoclast differentiation, 3) Pro-inflammatory monocytes/macrophages rise during two expansion waves, which can be explained by the anti-inflammatory macrophages-mediated inhibition phase between the two waves. In addition, we further tested the robustness of the mathematical model by comparing simulation results to an independent experimental dataset. Taken together, this novel comprehensive mathematical framework allowed us to identify biological mechanisms that best recapitulate bone injury data and that explain the coupled cellular population dynamics involved in the process. Furthermore, our hypothesis testing methodology could be used in other contexts to decipher mechanisms in complex multicellular processes.

Histone deacetylase inhibition prevents the growth of primary and metastatic osteosarcoma.

Overall survival rates for patients with advanced osteosarcoma have remained static for over three decades. An in vitro analysis of osteosarcoma cell lines for sensitivity to an array of approved cancer therapies revealed that panobinostat, a broad spectrum histone deacetalyase (HDAC) inhibitor, is highly effective at triggering osteosarcoma cell death. Using in vivo models of orthotopic and metastatic osteosarcoma, here we report that panobinostat impairs the growth of primary osteosarcoma in bone and spontaneous metastasis to the lung, the most common site of metastasis for this disease. Further, pretreatment of mice with panobinostat prior to tail vein inoculation of osteosarcoma prevents the seeding and growth of lung metastases. Additionally, panobinostat impaired the growth of established lung metastases and improved overall survival, and these effects were also manifest in the lung metastatic SAOS2-LM7 model. Mechanistically, the efficacy of panobinostat was linked to high expression of HDAC1 and HDAC2 in osteosarcoma, and silencing of HDAC1 and 2 greatly reduced osteosarcoma growth in vitro. In accordance with these findings, treatment with the HDAC1/2 selective inhibitor romidepsin compromised the growth of osteosarcoma in vitro and in vivo. Analysis of patient-derived xenograft osteosarcoma cell lines further demonstrated the sensitivity of the disease to panobinostat or romidepsin. Collectively, these studies provide rationale for clinical trials in osteosarcoma patients using the approved therapies panobinostat or romidepsin.

Proteolytic Regulation of Parathyroid Hormone-Related Protein: Functional Implications for Skeletal Malignancy.

Parathyroid hormone-related protein (PTHrP), with isoforms ranging from 139 to 173 amino acids, has long been implicated in the development and regulation of multiple tissues, including that of the skeleton, via paracrine and autocrine signaling. PTHrP is also known as a potent mediator of cancer-induced bone disease, contributing to a vicious cycle between tumor cells and the bone microenvironment that drives the formation and progression of metastatic lesions. The abundance of roles ascribed to PTHrP have largely been attributed to the N-terminal 1-36 amino acid region, however, activities for mid-region and C-terminal products as well as additional shorter N-terminal species have also been described. Studies of the protein sequence have indicated that PTHrP is susceptible to post-translational proteolytic cleavage by multiple classes of proteases with emerging evidence pointing to novel functional roles for these PTHrP products in regulating cell behavior in homeostatic and pathological contexts. As a consequence, PTHrP products are also being explored as potential biomarkers of disease. Taken together, our enhanced understanding of the post-translational regulation of PTHrP bioactivity could assist in developing new therapeutic approaches that can effectively treat skeletal malignancies.

Size Matters: Metastatic Cluster Size and Stromal Recruitment in the Establishment of Successful Prostate Cancer to Bone Metastases.

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.

New approaches to selectively target cancer-associated matrix metalloproteinase activity.

Heightened matrix metalloproteinase (MMP) activity has been noted in the context of the tumor microenvironment for many years, and causal roles for MMPs have been defined across the spectrum of cancer progression. This is primarily due to the ability of the MMPs to process extracellular matrix (ECM) components and to regulate the bioavailability/activity of a large repertoire of cytokines and growth factors. These characteristics made MMPs an attractive target for therapeutic intervention but notably clinical trials performed in the 1990s did not fulfill the promise of preclinical studies. The reason for the failure of early MMP inhibitor (MMPI) clinical trials that are multifold but arguably principal among them was the inability of early MMP-based inhibitors to selectively target individual MMPs and to distinguish between MMPs and other members of the metzincin family. In the decades that have followed the MMP inhibitor trials, innovations in chemical design, antibody-based strategies, and nanotechnologies have greatly enhanced our ability to specifically target and measure the activity of MMPs. These advances provide us with the opportunity to generate new lines of highly selective MMPIs that will not only extend the overall survival of cancer patients, but will also afford us the ability to utilize heightened MMP activity in the tumor microenvironment as a means by which to deliver MMPIs or MMP activatable prodrugs.

Integrating new discoveries into the "vicious cycle" paradigm of prostate to bone metastases.

In prostate to bone metastases, the "vicious cycle" paradigm has been traditionally used to illustrate how metastases manipulate the bone forming osteoblasts and resorbing osteoclasts in order to yield factors that facilitate growth and establishment. However, recent advances have illustrated that the cycle is far more complex than this simple interpretation. In this review, we will discuss the role of exosomes and hematopoietic/mesenchymal stem/stromal cells (MSC) that facilitate the establishment and activation of prostate metastases and how cells including myeloid-derived suppressor cells, macrophages, T cells, and nerve cells contribute to the momentum of the vicious cycle. The increased complexity of the tumor-bone microenvironment requires a system level approach. The evolution of computational models to interrogate the tumor-bone microenvironment is also discussed, and the application of this integrated approach should allow for the development of effective therapies to treat and cure prostate to bone metastases.

Current and emerging therapies for bone metastatic castration-resistant prostate cancer.

BACKGROUND: A paucity of therapeutic options is available to treat men with metastatic castration-resistant prostate cancer (mCRPC). However, recent developments in our understanding of the disease have resulted in several new therapies that show promise in improving overall survival rates in this patient population. METHODS: Agents approved for use in the United States and those undergoing clinical trials for the treatment of mCRPC are reviewed. Recent contributions to the understanding of prostate biology and bone metastasis are discussed as well as how the underlying mechanisms may represent opportunities for therapeutic intervention. New challenges to delivering effective mCRPC treatment will also be examined. RESULTS: New and emerging treatments that target androgen synthesis and utilization or the microenvironment may improve overall survival rates for men diagnosed with mCRPC. Determining how factors derived from the primary tumor can promote the development of premetastatic niches and how prostate cancer cells parasitize niches in the bone microenvironment, thus remaining dormant and protected from systemic therapy, could yield new therapies to treat mCRPC. Challenges such as intratumoral heterogeneity and patient selection can potentially be circumvented via computational biology approaches. CONCLUSIONS: The emergence of novel treatments for mCRPC, combined with improved patient stratification and optimized therapy sequencing, suggests that significant gains may be made in terms of overall survival rates for men diagnosed with this form of cancer.

Rapid ex vivo imaging of PAIII prostate to bone tumor with SWIFT-MRI.

PURPOSE: The limiting factor for MRI of skeletal/mineralized tissue is fast transverse relaxation. A recent advancement in MRI technology, SWIFT (Sweep Imaging with Fourier Transform), is emerging as a new approach to overcome this difficulty. Among other techniques like UTE, ZTE, and WASPI, the application of SWIFT technology has the strong potential to impact preclinical and clinical imaging, particularly in the context of primary or metastatic bone cancers because it has the added advantage of imaging water in mineralized tissues of bone allowing MRI images to be obtained of tissues previously visible only with modalities such as computed tomography (CT). The goal of the current study is to examine the feasibility of SWIFT for the assessment of the prostate cancer induced changes in bone formation (osteogenesis) and destruction (osteolysis) in ex vivo specimens. METHODS: A luciferase expressing prostate cancer cell line (PAIII) or saline control was inoculated directly into the tibia of 6-week-old immunocompromised male mice. Tumor growth was assessed weekly for 3 weeks before euthanasia and dissection of the tumor bearing and sham tibias. The ex vivo mouse tibia specimens were imaged with a 9.4 Tesla (T) and 7T MRI systems. SWIFT images are compared with traditional gradient-echo and spin-echo MRI images as well as CT and histological sections. RESULTS: SWIFT images with nominal resolution of 78 µm are obtained with the tumor and different bone structures identified. Prostate cancer induced changes in the bone microstructure are visible in SWIFT images, which is supported by spin-echo, high resolution CT and histological analysis. CONCLUSION: SWIFT MRI is capable of high-quality high-resolution ex vivo imaging of bone tumor and surrounding bone and soft tissues. Furthermore, SWIFT MRI shows promise for in vivo bone tumor imaging, with the added benefits of nonexposure to ionizing radiation, quietness, and speed.

The bone ecosystem facilitates multiple myeloma relapse and the evolution of heterogeneous drug resistant disease.

Multiple myeloma (MM) is an osteolytic malignancy that is incurable due to the emergence of treatment resistant disease. Defining how, when and where myeloma cell intrinsic and extrinsic bone microenvironmental mechanisms cause relapse is challenging with current biological approaches. Here, we report a biology-driven spatiotemporal hybrid agent-based model of the MM-bone microenvironment. Results indicate MM intrinsic mechanisms drive the evolution of treatment resistant disease but that the protective effects of bone microenvironment mediated drug resistance (EMDR) significantly enhances the probability and heterogeneity of resistant clones arising under treatment. Further, the model predicts that targeting of EMDR deepens therapy response by eliminating sensitive clones proximal to stroma and bone, a finding supported by in vivo studies. Altogether, our model allows for the study of MM clonal evolution over time in the bone microenvironment and will be beneficial for optimizing treatment efficacy so as to significantly delay disease relapse.

MMP-7 promotes prostate cancer-induced osteolysis via the solubilization of RANKL.

We developed a rodent model that mimics the osteoblastic and osteolytic changes associated with human metastatic prostate cancer. Microarray analysis identified MMP-7, cathepsin-K, and apolipoprotein D as being upregulated at the tumor-bone interface. MMP-7, which was produced by osteoclasts at the tumor-bone interface, was capable of processing RANKL to a soluble form that promoted osteoclast activation. MMP-7-deficient mice demonstrated reduced prostate tumor-induced osteolysis and RANKL processing. This study suggests that inhibition of MMP-7 will have therapeutic benefit in the treatment of prostate cancer-induced osteolysis.

γδ-Enriched CAR-T cell therapy for bone metastatic castrate-resistant prostate cancer.

Immune checkpoint blockade has been largely unsuccessful for the treatment of bone metastatic castrate-resistant prostate cancer (mCRPC). Here, we report a combinatorial strategy to treat mCRPC using γδ-enriched chimeric antigen receptor (CAR) T cells and zoledronate (ZOL). In a preclinical murine model of bone mCRPC, γδ CAR-T cells targeting prostate stem cell antigen (PSCA) induced a rapid and significant regression of established tumors, combined with increased survival and reduced cancer-associated bone disease. Pretreatment with ZOL, a U.S. Food and Drug Administration-approved bisphosphonate prescribed to mitigate pathological fracture in mCRPC patients, resulted in CAR-independent activation of γδ CAR-T cells, increased cytokine secretion, and enhanced antitumor efficacy. These data show that the activity of the endogenous Vγ9Vδ2 T cell receptor is preserved in CAR-T cells, allowing for dual-receptor recognition of tumor cells. Collectively, our findings support the use of γδ CAR-T cell therapy for mCRPC treatment.

MMP-2 contributes to the development of the mouse ventral prostate by impacting epithelial growth and morphogenesis.

Epithelial growth, branching, and canalization are important morphogenetic events of the rodent ventral prostate (VP) that take place during the first postnatal week. In this study, we evaluated the effect of knocking out MMP-2 (MMP-2(-/-)), by examining developmental and structural aspects of the VP in MMP-2(-/-) mice. Neonate (day 6) MMP-2(-/-) mice showed fewer epithelial tips, a lower epithelial cell proliferation rate, and also reticulin fiber accumulation. The VP of adult MMP-2(-/-) mice showed lower relative weight, smaller epithelial and smooth-muscle cell volume, and a larger amount of thicker reticulin fibers. No differences in cell proliferation or apoptotic index were noted between adult MMP-2(-/-) and wild-type mice. MMP-9 was found in the adult MMP-2(-/-), but not in the wild-type. In conclusion, MMP-2 function is essential for the epithelial morphogenesis of the mouse VP, and expression of MMP-9 is not sufficient for acquisition of the normal adult histology.

Computational modeling reveals a key role for polarized myeloid cells in controlling osteoclast activity during bone injury repair.

Bone-forming osteoblasts and -resorbing osteoclasts control bone injury repair, and myeloid-derived cells such as monocytes and macrophages are known to influence their behavior. However, precisely how these multiple cell types coordinate and regulate each other over time within the bone marrow to restore bone is difficult to dissect using biological approaches. Conversely, mathematical modeling lends itself well to this challenge. Therefore, we generated an ordinary differential equation (ODE) model powered by experimental data (osteoblast, osteoclast, bone volume, pro- and anti-inflammatory myeloid cells) obtained from intra-tibially injured mice. Initial ODE results using only osteoblast/osteoclast populations demonstrated that bone homeostasis could not be recovered after injury, but this issue was resolved upon integration of pro- and anti-inflammatory myeloid population dynamics. Surprisingly, the ODE revealed temporal disconnects between the peak of total bone mineralization/resorption, and osteoblast/osteoclast numbers. Specifically, the model indicated that osteoclast activity must vary greatly (> 17-fold) to return the bone volume to baseline after injury and suggest that osteoblast/osteoclast number alone is insufficient to predict bone the trajectory of bone repair. Importantly, the values of osteoclast activity fall within those published previously. These data underscore the value of mathematical modeling approaches to understand and reveal new insights into complex biological processes.

EphA2 Is a Clinically Relevant Target for Breast Cancer Bone Metastatic Disease.

EphA2 receptor tyrosine kinase (RTK) is highly expressed in breast tumor cells across multiple molecular subtypes and correlates with poor patient prognosis. In this study, the potential role of EphA2 in this clinically relevant phenomenon is investigated as metastasis of breast cancer to bone is a major cause of morbidity and mortality in patients. It was found that the EphA2 function in breast cancer cells promotes osteoclast activation and the development of osteolytic bone disease. Blocking EphA2 function molecularly and pharmacologically in breast tumors reduced the number and size of bone lesions and the degree of osteolytic disease in intratibial and intracardiac mouse models, which correlated with a significant decrease in the number of osteoclasts at the tumor-bone interface. EphA2 loss of function in tumor cells impaired osteoclast progenitor differentiation in coculture, which is mediated, at least in part, by reduced expression of IL-6. EPHA2 transcript levels are enriched in human breast cancer bone metastatic lesions relative to visceral metastatic sites; EphA2 protein expression was detected in breast tumor cells in bone metastases in patient samples, supporting the clinical relevance of the study's findings. These data provide a strong rationale for the development and application of molecularly targeted therapies against EphA2 for the treatment of breast cancer bone metastatic disease. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC. on behalf of American Society for Bone and Mineral Research.

Host-Derived Matrix Metalloproteinase-13 Activity Promotes Multiple Myeloma-Induced Osteolysis and Reduces Overall Survival.

Multiple myeloma promotes systemic skeletal bone disease that greatly contributes to patient morbidity. Resorption of type I collagen-rich bone matrix by activated osteoclasts results in the release of sequestered growth factors that can drive progression of the disease. Matrix metalloproteinase-13 (MMP13) is a collagenase expressed predominantly in the skeleton by mesenchymal stromal cells (MSC) and MSC-derived osteoblasts. Histochemical analysis of human multiple myeloma specimens also demonstrated that MMP13 largely localizes to the stromal compartment compared with CD138+ myeloma cells. In this study, we further identified that multiple myeloma induces MMP13 expression in bone stromal cells. Because of its ability to degrade type I collagen, we examined whether bone stromal-derived MMP13 contributed to myeloma progression. Multiple myeloma cells were inoculated into wild-type or MMP13-null mice. In independent in vivo studies, MMP13-null mice demonstrated significantly higher overall survival rates and lower levels of bone destruction compared with wild-type controls. Unexpectedly, no differences in type I collagen processing between the groups were observed. Ex vivo stromal coculture assays showed reduced formation and activity in MMP13-null osteoclasts. Analysis of soluble factors from wild-type and MMP13-null MSCs revealed decreased bioavailability of various osteoclastogenic factors including CXCL7. CXCL7 was identified as a novel MMP13 substrate and regulator of osteoclastogenesis. Underscoring the importance of host MMP13 catalytic activity in multiple myeloma progression, we demonstrate the in vivo efficacy of a novel and highly selective MMP13 inhibitor that provides a translational opportunity for the treatment of this incurable disease. SIGNIFICANCE: Genetic and pharmacologic approaches show that bone stromal-derived MMP13 catalytic activity is critical for osteoclastogenesis, bone destruction, and disease progression. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/81/9/2415/F1.large.jpg.

Mesenchymal stem cell-derived interleukin-28 drives the selection of apoptosis resistant bone metastatic prostate cancer.

Bone metastatic prostate cancer (PCa) promotes mesenchymal stem cell (MSC) recruitment and their differentiation into osteoblasts. However, the effects of bone-marrow derived MSCs on PCa cells are less explored. Here, we report MSC-derived interleukin-28 (IL-28) triggers prostate cancer cell apoptosis via IL-28 receptor alpha (IL-28Rα)-STAT1 signaling. However, chronic exposure to MSCs drives the selection of prostate cancer cells that are resistant to IL-28-induced apoptosis and therapeutics such as docetaxel. Further, MSC-selected/IL-28-resistant prostate cancer cells grow at accelerated rates in bone. Acquired resistance to apoptosis is PCa cell intrinsic, and is associated with a shift in IL-28Rα signaling via STAT1 to STAT3. Notably, STAT3 ablation or inhibition impairs MSC-selected prostate cancer cell growth and survival. Thus, bone marrow MSCs drive the emergence of therapy-resistant bone metastatic prostate cancer yet this can be disabled by targeting STAT3.

Quantification and Optimization of Standard-of-Care Therapy to Delay the Emergence of Resistant Bone Metastatic Prostate Cancer.

BACKGROUND: Bone metastatic prostate cancer (BMPCa), despite the initial responsiveness to androgen deprivation therapy (ADT), inevitably becomes resistant. Recent clinical trials with upfront treatment of ADT combined with chemotherapy or novel hormonal therapies (NHTs) have extended overall patient survival. These results indicate that there is significant potential for the optimization of standard-of-care therapies to delay the emergence of progressive metastatic disease. METHODS: Here, we used data extracted from human bone metastatic biopsies pre- and post-abiraterone acetate/prednisone to generate a mathematical model of bone metastatic prostate cancer that can unravel the treatment impact on disease progression. Intra-tumor heterogeneity in regard to ADT and chemotherapy resistance was derived from biopsy data at a cellular level, permitting the model to track the dynamics of resistant phenotypes in response to treatment from biological first-principles without relying on data fitting. These cellular data were mathematically correlated with a clinical proxy for tumor burden, utilizing prostate-specific antigen (PSA) production as an example. RESULTS: Using this correlation, our model recapitulated the individual patient response to applied treatments in a separate and independent cohort of patients (n = 24), and was able to estimate the initial resistance to the ADT of each patient. Combined with an intervention-decision algorithm informed by patient-specific prediction of initial resistance, we propose to optimize the sequence of treatments for each patient with the goal of delaying the evolution of resistant disease and limit cancer cell growth, offering evidence for an improvement against retrospective data. CONCLUSIONS: Our results show how minimal but widely available patient information can be used to model and track the progression of BMPCa in real time, offering a clinically relevant insight into the patient-specific evolutionary dynamics of the disease and suggesting new therapeutic options for intervention. TRIAL REGISTRATION: NCT # 01953640. FUNDING: Funded by an NCI U01 (NCI) U01CA202958-01 and a Moffitt Team Science Award. CCL and DB were partly funded by an NCI PSON U01 (U01CA244101). AA was partly funded by a Department of Defense Prostate Cancer Research Program (W81XWH-15-1-0184) fellowship. LC was partly funded by a postdoctoral fellowship (PF-13-175-01-CSM) from the American Cancer Society.

Prostate cancer-derived MMP-3 controls intrinsic cell growth and extrinsic angiogenesis.

Bone metastatic prostate cancer significantly impacts patient quality of life and overall survival, and despite available therapies, it is presently incurable with an unmet need for improved treatment options. As mediators of tumor progression, matrix metalloproteinases (MMPs) can degrade extracellular matrix components and regulate growth factor and cytokine bioactivity. Depending on tissue context, MMPs can either promote or inhibit tumorigenesis. Therefore, it is essential to study individual MMPs in specific cancer contexts and microenvironments to support the design and application of selective MMP inhibitors. Here we report that tumor-derived MMP-3 contributes to bone metastatic prostate cancer progression via intrinsic and extrinsic routes. MMP-3 ablation in prostate cancer cell lines significantly reduced in vitro growth combined with lowered AKT and ERK phosphorylation and total VEGFR1 and FGFR3 protein levels. In vivo, MMP-3 ablated tumors grew at a slower rate and were significantly less vascularized. Quantitative PCR analyses of wild type and MMP-3 silenced prostate cancer cells also demonstrate downregulation of a wide array of angiogenic factors. The extrinsic role for MMP-3 in angiogenesis was supported by in vitro endothelial tube formation assays where the lack of MMP-3 in prostate cancer conditioned media resulted in slower rates of tube formation. Taken together, our results suggest that tumor-derived MMP-3 contributes to prostate cancer growth in bone. These data indicate that selective inhibition of MMP-3 and/or targeting MMP generated products could be efficacious for the treatment of prostate to bone metastases.

Modeling Osteocyte Network Formation: Healthy and Cancerous Environments.

Advanced cancers, such as prostate and breast cancers, commonly metastasize to bone. In the bone matrix, dendritic osteocytes form a spatial network allowing communication between osteocytes and the osteoblasts located on the bone surface. This communication network facilitates coordinated bone remodeling. In the presence of a cancerous microenvironment, the topology of this network changes. In those situations, osteocytes often appear to be either overdifferentiated (i.e., there are more dendrites than healthy bone) or underdeveloped (i.e., dendrites do not fully form). In addition to structural changes, histological sections from metastatic breast cancer xenografted mice show that number of osteocytes per unit area is different between healthy bone and cancerous bone. We present a stochastic agent-based model for bone formation incorporating osteoblasts and osteocytes that allows us to probe both network structure and density of osteocytes in bone. Our model both allows for the simulation of our spatial network model and analysis of mean-field equations in the form of integro-partial differential equations. We considered variations of our model to study specific physiological hypotheses related to osteoblast differentiation; for example predicting how changing biological parameters, such as rates of bone secretion, rates of cancer formation, and rates of osteoblast differentiation can allow for qualitatively different network topologies. We then used our model to explore how commonly applied therapies such as bisphosphonates (e.g., zoledronic acid) impact osteocyte network formation.