An integrative model of prostate cancer interaction with the bone microenvironment

Tumor growth for remodeling process: A 2D approach

This paper aims to present a comprehensive framework for coupling tumor-bone remodeling processes in a 2-dimensional geometry. This is achieved by introducing a bio-inspired damage that represents the growing tumor, which subsequently affects the main populations involved in the remodeling process, namely, osteoclasts, osteoblasts, and bone tissue. The model is constructed using a set of differential equations based on the Komarova's and Ayati's models, modified to incorporate the bio-inspired damage that may result in tumor mass formation. Three distinct models were developed. The first two models are based on the Komarova's governing equations, with one demonstrating an osteolytic behavior and the second one an osteoblastic model. The third model is a variation of Ayati's model, where the bio-inspired damage is induced through the paracrine and autocrine parameters, exhibiting an osteolytic behavior. The obtained results are consistent with existing literature, leading us to believe that our in-silico experiments will serve as a cornerstone for paving the way towards targeted interventions and personalized treatment strategies, ultimately improving the quality of life for those affected by these conditions.

Mathematical modeling of the effects of Wnt-10b on bone metabolism

Bone health is determined by factors including bone metabolism or remodeling. Wnt-10b alters osteoblastogenesis through pre-osteoblast proliferation and differentiation and osteoblast apoptosis rate, which collectively lead to the increase of bone density. To model this, we adapted a previously published model of bone remodeling. The resulting model for the bone compartment includes differential equations for active osteoclasts, pre-osteoblasts, osteoblasts, osteocytes, and the amount of bone present at the remodeling site. Our alterations to the original model consist of extending it past a single remodeling cycle and implementing a direct relationship to Wnt-10b. Four new parameters were estimated and validated using normalized data from mice. The model connects Wnt-10b to bone metabolism and predicts the change in trabecular bone volume caused by a change in Wnt-10b input. We find that this model predicts the expected increase in pre-osteoblasts and osteoblasts while also pointing to a decrease in osteoclasts when Wnt-10b is increased.

Experimental-based mechanobiological modeling of the anabolic and catabolic effects of breast cancer on bone remodeling

Bone is a biological tissue characterized by its hierarchical organization. This material has the ability to be continually renewed, which makes it highly adaptative to external loadings. Bone renewing is managed by a dynamic biological process called bone remodeling (BR), where continuous resorption of old bone and formation of new bone permits to change the bone composition and microstructure. Unfortunately, because of several factors, such as age, hormonal imbalance, and a variety of pathologies including cancer metastases, this process can be disturbed leading to various bone diseases. In this study, we have investigated the effect of breast cancer (BC) metastases causing osteolytic bone loss. BC has the ability to affect bone quantity in different ways in each of its primary and secondary stages. Based on a BR mathematical model, we modeled the BC cells' interaction with bone cells to assess their effect on bone volume fraction (BV/TV) evolution during the remodeling process. Some of the parameters used in our model have been determined experimentally using the enzyme-linked immune-sorbent assay (ELISA) and the MTT assay. Our numerical simulations show that primary BC plays a significant role in enhancing bone-forming cells' activity leading to a 6.22% increase in BV/TV over 1 year. On the other hand, secondary BC causes a noticeable decrease in BV/TV reaching 15.74% over 2 years.

Antigen receptor therapy in bone metastasis via optimal control for different human life stages

In this work we propose a bone metastasis model using power law growth functions in order to describe the biochemical interactions between bone cells and cancer cells. Experimental studies indicate that bone remodeling cycles are different for human life stages: childhood, young adulthood, and adulthood. In order to include such differences in our study, we estimate the model parameter values for each human life stage via bifurcation analysis. Results reveal an intrinsic relationship between the active period of remodeling cycles and the proliferation of cancer cells. Subsequently, using optimal control theory we analyze a possible antigen receptor therapy as a new treatment for bone metastasis. Theoretical results such as existence of optimal solutions are proved. Numerical simulations for late stages of bone metastasis are presented and a discussion of our results is carried out.

Developing and studying the dynamical behavior of a nonlinear mathematical model for cancers with tumor by considering immune system role

We are constrained by widespread cancerous diseases to improve treatment methods which save patients and provide better living conditions during and after the treatment period. Because of the complexity of the treatment process, mathematical models need to be used in order to have a better understanding of the process. However, deriving an adequate complex model that can capture the disease pattern which could be confirmed by simulations and experiments has its own barriers. In this paper, a new mathematical model is developed concerning immune system effect on cancer. The model is introduced using nonlinear ordinary differential equations. Also, the qualitative behavior of the proposed system is studied in order to examine the extent of the model with respect to the nature of tumor evolution. Thus, number and status of equilibria points in line with the existence of limit cycles are obtained for sub-systems and the whole system. Meanwhile, possible bifurcations are mentioned, and the consequent evolutions are described. It is shown that the model conforms well to natural possibilities, cancer growth or remission. Thus, the model would be fit for further studies for prediction and contemplating treatment method, especially for immune stimulating drugs and immunotherapy.

Toward a Mathematical Modeling of Diseases' Impact on Bone Remodeling: Technical Review

A wide variety of bone diseases have hitherto been discovered, such as osteoporosis, Paget's disease, osteopetrosis, and metastatic bone disease, which are not well defined in terms of changes in biochemical and mechanobiological regulatory factors. Some of these diseases are secondary to other pathologies, including cancer, or to some clinical treatments. To better understand bone behavior and prevent its deterioration, bone biomechanics have been the subject of mathematical modeling that exponentially increased over the last years. These models are becoming increasingly complex. The current paper provides a timely and critical analysis of previously developed bone remodeling mathematical models, particularly those addressing bone diseases. Besides, mechanistic pharmacokinetic/pharmacodynamic (PK/PD) models, which englobe bone disease and its treatment's effect on bone health. Therefore, the review starts by presenting bone remodeling cycle and mathematical models describing this process, followed by introducing some bone diseases and discussing models of pathological mechanisms affecting bone, and concludes with exhibiting the available bone treatment procedures considered in the PK/PD models.

Review: Mathematical Modeling of Prostate Cancer and Clinical Application

We review and synthesize key findings and limitations of mathematical models for prostate cancer, both from theoretical work and data-validated approaches, especially concerning clinical applications. Our focus is on models of prostate cancer dynamics under treatment, particularly with a view toward optimizing hormone-based treatment schedules and estimating the onset of treatment resistance under various assumptions. Population models suggest that intermittent or adaptive therapy is more beneficial to delay cancer relapse as compared to the standard continuous therapy if treatment resistance comes at a competitive cost for cancer cells. Another consensus among existing work is that the standard biomarker for cancer growth, prostate-specific antigen, may not always correlate well with cancer progression. Instead, its doubling rate appears to be a better indicator of tumor growth. Much of the existing work utilizes simple ordinary differential equations due to difficulty in collecting spatial data and due to the early success of using prostate-specific antigen in mathematical modeling. However, a shift toward more complex and realistic models is taking place, which leaves many of the theoretical and mathematical questions unexplored. Furthermore, as adaptive therapy displays better potential than existing treatment protocols, an increasing number of studies incorporate this treatment into modeling efforts. Although existing modeling work has explored and yielded useful insights on the treatment of prostate cancer, the road to clinical application is still elusive. Among the pertinent issues needed to be addressed to bridge the gap from modeling work to clinical application are (1) real-time data validation and model identification, (2) sensitivity analysis and uncertainty quantification for model prediction, and (3) optimal treatment/schedule while considering drug properties, interactions, and toxicity. To address these issues, we suggest in-depth studies on various aspects of the parameters in dynamical models such as the evolution of parameters over time. We hope this review will assist future attempts at studying prostate cancer.

Mathematical modelling of the role of Endo180 network in the development of metastatic bone disease in prostate cancer

Metastatic bone disease (MBD) is a common complication of advanced cancer and recent research suggests that Endo180 expression is dysregulated through the TGFβ-TGFβR-SMAD2/3 signalling pathway during the invasion of tumour cells in the development of MBD. We here provide a model for the dysregulation of the Endo180 network to demonstrate its vital contribution to bone destruction as well as tumour cell growth. The model consisted of a set of ordinary differential equations and reconstructed variations in the bone cells, resultant bone volume, and biochemical factors involved in the TGFβ-TGFβR-SMAD2/3 signalling pathway over time. The model also investigated the underlying mechanism in which the change of TGFβ affects the TGFβ-TGFβR-SMAD2/3 signalling pathway and the resultant Endo180 expression in osteoblastic and tumour cells. The model links the appearance of tumour cells with the inhibition of TGFβ binding to its receptors on osteoblastic cells, to affect TGFβ-TGFβR-SMAD2/3 signalling and Endo180 expression. Temporal variation in bone cells, bone volume, and the biochemical factors involved in the TGFβ-TGFβR-SMAD2/3 pathway as demonstrated in the model simulations agree with published experimental data. The model can be refined based on further discoveries but allows the influence of Endo180 network dysregulation on bone remodelling in MBD to be established. This model could aid in the development of Endo180 targeted therapies for MBD in the future.

Mechanistic modelling of prostate-specific antigen dynamics shows potential for personalized prediction of radiation therapy outcome

External beam radiation therapy is a widespread treatment for prostate cancer. The ensuing patient follow-up is based on the evolution of the prostate-specific antigen (PSA). Serum levels of PSA decay due to the radiation-induced death of tumour cells and cancer recurrence usually manifest as a rising PSA. The current definition of biochemical relapse requires that PSA reaches nadir and starts increasing, which delays the use of further treatments. Also, these methods do not account for the post-radiation tumour dynamics that may contain early information on cancer recurrence. Here, we develop three mechanistic models of post-radiation PSA evolution. Our models render superior fits of PSA data in a patient cohort and provide a biological justification for the most common empirical formulation of PSA dynamics. We also found three model-based prognostic variables: the proliferation rate of the survival fraction, the ratio of radiation-induced cell death rate to the survival proliferation rate, and the time to PSA nadir since treatment termination. We argue that these markers may enable the early identification of biochemical relapse, which would permit physicians to subsequently adapt patient monitoring to optimize the detection and treatment of cancer recurrence.

Pathological Bases and Clinical Impact of Intratumor Heterogeneity in Clear Cell Renal Cell Carcinoma

PURPOSE OF REVIEW

Intratumor heterogeneity is an inherent event in tumor development that is receiving much attention in the last years since it is responsible for most failures of current targeted therapies. The purpose of this review is to offer clinicians an updated insight of the multiple manifestations of a complex event that impacts significantly patient's life.

RECENT FINDINGS

Clear cell renal cell carcinoma is the most common renal tumor and a paradigmatic example of a heterogeneous neoplasm. Next-generation sequencing has demonstrated that intratumor heterogeneity encompasses genetic, epigenetic, and microenvironmental variability. Currently accepted protocols of tumor sampling seem insufficient in unveiling intratumor heterogeneity with reliability and need to be updated. This variability challenges the precise morphological diagnosis, its molecular characterization, and the selection of optimal personalized therapies in clear cell renal cell carcinoma, a neoplasm traditionally considered chemo- and radio-resistant. We review the state of the art of the different approaches to intratumor heterogeneity in clear cell renal cell carcinomas, from the simple morphology to the most sophisticated massive sequencing tools.