Mathematical modelling of the pathogenesis of multiple myeloma-induced bone disease

Osteoporosis induced by cellular senescence: A mathematical model

Osteoporosis is a disease characterized by loss of bone mass, where bones become fragile and more likely to fracture. Bone density begins to decrease at age 50, and a state of osteoporosis is defined by loss of more than 25%. Cellular senescence is a permanent arrest of normal cell cycle, while maintaining cell viability. The number of senescent cells increase with age. Since osteoporosis is an aging disease, it is natural to consider the question to what extend senescent cells induce bone density loss and osteoporosis. In this paper we use a mathematical model to address this question. We determine the percent of bone loss for men and women during age 50 to 100 years, and the results depend on the rate η of net formation of senescent cell, with η = 1 being the average rate. In the case η = 1, the model simulations are in agreement with empirical data. We also consider senolytic drugs, like fisetin and quercetin, that selectively eliminate senescent cells, and assess their efficacy in terms of reducing bone loss. For example, at η = 1, with estrogen hormonal therapy and early treatment with fisetin, bone density loss for women by age 75 is 23.4% (below osteoporosis), while with no treatment with fisetin it is 25.8% (osteoporosis); without even a treatment with estrogen hormonal therapy, bone loss of 25.3% occurs already at age 65.

A review of mathematical modeling of bone remodeling from a systems biology perspective

Bone remodeling is an essential, delicately balanced physiological process of coordinated activity of bone cells that remove and deposit new bone tissue in the adult skeleton. Due to the complex nature of this process, many mathematical models of bone remodeling have been developed. Each of these models has unique features, but they have underlying patterns. In this review, the authors highlight the important aspects frequently found in mathematical models for bone remodeling and discuss how and why these aspects are included when considering the physiology of the bone basic multicellular unit, which is the term used for the collection of cells responsible for bone remodeling. The review also emphasizes the view of bone remodeling from a systems biology perspective. Understanding the systemic mechanisms involved in remodeling will help provide information on bone pathology associated with aging, endocrine disorders, cancers, and inflammatory conditions and enhance systems pharmacology. Furthermore, some features of the bone remodeling cycle and interactions with other organ systems that have not yet been modeled mathematically are discussed as promising future directions in the field.

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.

Breast Cancer Exosomal microRNAs Facilitate Pre-Metastatic Niche Formation in the Bone: A Mathematical Model

Pre-metastatic niche is a location where cancer cells, separating from a primary tumor, find "fertile soil" for growth and proliferation, ensuring successful metastasis. Exosomal miRNAs of breast cancer are known to enter the bone and degrade it, which facilitates cancer cells invasion into the bone interior and ensures its successful colonization. In this paper, we use a mathematical model to first describe, in health, the continuous remodeling of the bone by bone-forming osteoblasts, bone-resorbing osteoclasts and the RANKL-OPG-RANK signaling system, which keeps the balance between bone formation and bone resorption. We next demonstrate how breast cancer exosomal miRNAs disrupt this balance, either by increasing or by decreasing the ratio of osteoclasts/osteoblasts, which results in abnormal high bone resorption or abnormal high bone forming, respectively, and in bone weakening in both cases. Finally we consider the case of abnormally high resorption and evaluate the effect of drugs, which may increase bone density to normal level, thus protecting the bone from invasion by cancer cells.

Spatial extension of a bone remodeling dynamics model and its finite element analysis

There are many works dealing with the dynamics of bone remodeling, proposing increasingly complex and complete models. In the recent years, the efforts started to focus on developing models that not only reproduce the temporal evolution, but also include the spatial aspects of this phenomenon. In this work, we propose the spatial extension of an existing model that includes the dynamics of osteocytes. The spatial dependence is modeled in terms of a linear diffusion, as proposed in previous works dealing with related problems. The resulting model is then written in its variational form, and fully discretized using the well-known finite element method and a combination of the implicit and explicit Euler schemes. The numerical algorithm is then analyzed, proving some a priori error estimates and its linear convergence. Finally, we extend the examples already published for the temporal model to one and two dimensions, showing the dynamics of the solution in the spatial domain.

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.

Mathematical modelling of the role of GADD45β in the pathogenesis of multiple myeloma

Multiple myeloma (MM) is an incurable disease with relatively high morbidity and mortality rates. Great efforts were made to develop nuclear factor-kappa B (NF-κB)-targeted therapies against MM disease. However, these treatments influence MM cells as well as normal cells, inevitably causing serious side effects. Further research showed that NF-κB signalling promotes the survival of MM cells by interacting with JNK signalling through growth arrest and DNA damage-inducible beta (GADD45β), the downstream module of NF-κB signalling. The GADD45β-targeted intervention was suggested to be an effective and MM cell-specific treatment. However, the underlying mechanism through which GADD45β promotes the survival of MM cells is usually ignored in the previous models. A mathematical model of MM is built in this paper to investigate how NF-κB signalling acts along with JNK signalling through GADD45β and MKK7 to promote the survival of MM cells. The model cannot only mimic the variations in bone cells, the bone volume and MM cells with time, but it can also examine how the NF-κB pathway acts with the JNK pathway to promote the development of MM cells. In addition, the model also investigates the efficacies of GADD45β- and NF-κB-targeted treatments, suggesting that GADD45β-targeted therapy is more effective but has no apparent side effects. The simulation results match the experimental observations. It is anticipated that this model could be employed as a useful tool to initially investigate and even explore potential therapies involving the NF-κB and JNK pathways in the future.

Modeling the Human Bone-Tumor Niche: Reducing and Replacing the Need for Animal Data

Bone is the most common site for cancer metastasis. Understanding the interactions within the complex, heterogeneous bone-tumor microenvironment is essential for the development of new therapeutics. Various animal models of tumor-induced bone disease are routinely used to provide valuable information on the relationship between cancer cells and the skeleton. However, new model systems exist that offer an alternative approach to the use of animals and might more accurately reveal the cellular interactions occurring within the human bone-tumor niche. This review highlights replacement models that mimic the bone microenvironment and where cancer metastases and tumor growth might be assessed alongside bone turnover. Such culture models include the use of calcified regions of animal tissue and scaffolds made from bone mineral hydroxyapatite, synthetic polymers that can be manipulated during manufacture to create structures resembling trabecular bone surfaces, gel composites that can be modified for stiffness and porosity to resemble conditions in the tumor-bone microenvironment. Possibly the most accurate model system involves the use of fresh human bone samples, which can be cultured ex vivo in the presence of human tumor cells and demonstrate similar cancer cell-bone cell interactions as described in vivo. In addition, the use of mathematical modeling and computational biology approaches provide an alternative to preliminary animal testing. The use of such models offers the capacity to mimic significant elements of the human bone-tumor environment, and complement, refine, or replace the use of preclinical models. © 2020 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

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.

A Retrospective Analysis: A Novel Index Predicts Survival and Risk-Stratification for Bone Destruction in 419 Newly Diagnosed Multiple Myelomas

Objective

Multiple myeloma (MM) patients with bone destruction are difficult to restore, so it is of great clinical significance to further explore the factors affecting MM bone destruction.

Methods and results

This study retrospectively analyzed 419 cases with MM. Multiple linear regression analysis showed that those MM patients with a higher concentration of Ca²⁺ in serum, higher positive rate of CD138 immuno-phenotype and advanced in stage with 13q34 deletion in cytogenetics would be more prone to bone destruction, while total bile acid (TBA) and kappa chain isotope negatively correlated with bone destruction in MM patients. The Kaplan-Meier analysis indicated that Ca²⁺, serum β2-microglobulin (β2-MG), hemoglobin (HGB), creatinine (CREA), uric acid (UA) and age correlated with the survival of bone destruction in MM patients. Cox regression analysis further showed that the independent prognostic factors of β2-MG and CREA had a higher risk for early mortality in bone destruction patients. Moreover, an index was calculated based on β2-MG and globulin (GLB) to white blood cell (WBC) ratio to predict the poor survival of bone destruction patients.

Conclusion

We provide a novel marker to predict the prognosis of myeloma patients using routine examination method instead of bone marrow aspiration, and provide a reference for clinical evaluation.

Modeling the Effects of Multiple Myeloma on Kidney Function

Multiple myeloma (MM), a plasma cell cancer, is associated with many health challenges, including damage to the kidney by tubulointerstitial fibrosis. We develop a mathematical model which captures the qualitative behavior of the cell and protein populations involved. Specifically, we model the interaction between cells in the proximal tubule of the kidney, free light chains, renal fibroblasts, and myeloma cells. We analyze the model for steady-state solutions to find a mathematically and biologically relevant stable steady-state solution. This foundational model provides a representation of dynamics between key populations in tubulointerstitial fibrosis that demonstrates how these populations interact to affect patient prognosis in patients with MM and renal impairment.

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.

A virtual approach to evaluate therapies for management of multiple myeloma induced bone disease

Multiple myeloma bone disease is devastating for patients and a major cause of morbidity. The disease leads to bone destruction by inhibiting osteoblast activity while stimulating osteoclast activity. Recent advances in multiple myeloma research have improved our understanding of the pathogenesis of multiple myeloma-induced bone disease and suggest several potential therapeutic strategies. However, the effectiveness of some potential therapeutic strategies still requires further investigation and optimization. In this paper, a recently developed mathematical model is extended to mimic and then evaluate three therapies of the disease, namely: bisphosphonates, bortezomib and TGF-β inhibition. The model suggests that bisphosphonates and bortezomib treatments not only inhibit bone destruction, but also reduce the viability of myeloma cells. This contributes to the current debate as to whether bisphosphonate therapy has an anti-tumour effect. On the other hand, the analyses indicate that treatments designed to inhibit TGF-β do not reduce bone destruction, although it appears that they might reduce the viability of myeloma cells, which again contributes to the current controversy regarding the efficacy of TGF-β inhibition in multiple myeloma-induced bone disease.

A predator-prey based mathematical model of the bone remodelling cycle: exploring the relationship between the model parameters and biochemical factors

Bone remodelling is a vital process which enables bone to repair, renew and optimize itself. Disorders in the bone remodelling process are inevitably manifested in bone-related diseases, such as hypothyroidism, primary hyperparathyroidism and osteoporosis. In our previous work, a predator-prey based mathematical model was developed to simulate bone remodelling cycles under normal and two pathological conditions, hypothyroidism and primary hyperparathyroidism, for trabecular bone at a fixed point. However, the biochemical meanings of the model parameters were not fully explored. This article first extends the previous work by proposing relationships between the model parameters and biochemical factors involved in the bone remodelling process and by examining whether those relationships do predict the behaviours observed in vivo. The model is then applied to the simulation and investigation of bone remodelling of postmenopausal osteoporosis. The proposed connections are supported by good agreement between the model simulations and published experimental observations for the normal condition and all three pathological variations in bone remodelling.