Modeling the Macrophage-Mediated Inflammation Involved in the Bone Fracture Healing Process

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.

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, Patget’s disease, osteopetrosis, and metastatic bone disease, which is not well defined. 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 more and more involving the impact of bone diseases on its quality. The current paper provides a meaningful and critical analysis of previously developed bone remodeling mathematical models, particularly those addressing bone diseases. Besides, mechanistic pharmacokinetic/pharmacodynamic (PK/PD) models, which englobe either bone disease and its treatments’ 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 ends by exhibiting the available bone treatment procedures and the PK/PD models considering these lasts.

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

SUMMARY Multiple myeloma (MM) is the second most common haematological malignancy and results in destructive bone lesions. The interaction between MM cells and the bone microenvironment plays an important role in the development of the tumour cells and MM‐induced bone disease and forms a ‘vicious cycle’ of tumour development and bone destruction, intensified by suppression of osteoblast activity and promotion of osteoclast activity. In this paper, a mathematical model is proposed to simulate how the interaction between MM cells and the bone microenvironment facilitates the development of the tumour cells and the resultant bone destruction. It includes both the roles of inhibited osteoblast activity and stimulated osteoclast activity. The model is able to mimic the temporal variation of bone cell concentrations and resultant bone volume after the invasion and then removal of the tumour cells and explains why MM‐induced bone lesions rarely heal even after the complete removal of MM cells. The behaviour of the model compares well with published experimental data. The model serves as a first step to understand the development of MM‐induced bone disease and could be applied further to evaluate the current therapies against MM‐induced bone disease and even suggests new potential therapeutic targets. © 2014 The Authors. International Journal for Numerical Methods in Biomedical Engineering published by John Wiley & Sons Ltd

Mathematical modeling of postmenopausal osteoporosis and its treatment by the anti‐catabolic drug denosumab

SUMMARY Denosumab, a fully human monoclonal antibody, has been approved for the treatment of postmenopausal osteoporosis. The therapeutic effect of denosumab rests on its ability to inhibit osteoclast differentiation. Here, we present a computational approach on the basis of coupling a pharmacokinetics model of denosumab with a pharmacodynamics model for quantifying the effect of denosumab on bone remodeling. The pharmacodynamics model comprises an integrated systems biology‐continuum micromechanics approach, including a bone cell population model, considering the governing biochemical factors of bone remodeling (including the action of denosumab), and a multiscale micromechanics‐based bone mechanics model, for implementing the mechanobiology of bone remodeling in our model. Numerical studies of postmenopausal osteoporosis show that denosumab suppresses osteoclast differentiation, thus strongly curtailing bone resorption. Simulation results also suggest that denosumab may trigger a short‐term bone volume gain, which is, however, followed by constant or decreasing bone volume. This evolution is accompanied by a dramatic decrease of the bone turnover rate by more than one order of magnitude. The latter proposes dominant occurrence of secondary mineralization (which is not anymore impeded through cellular activity), leading to higher mineral concentration per bone volume. This explains the overall higher bone mineral density observed in denosumab‐related clinical studies. Copyright © 2013 John Wiley & Sons, Ltd.

Mathematical Modeling of Spatio-Temporal Dynamics of a Single Bone Multicellular Unit

Abstract During bone remodeling, bone-resorbing osteoclasts and bone-forming osteoblasts are organized in bone multicellular units (BMUs), which travel at a rate of 20–40 μm/d for 6–12 mo, maintaining a cylindrical structure. However, the interplay of local BMU geometry with biochemical regulation is poorly understood. We developed a mathematical model of BMU describing changes in time and space of the concentrations of proresorptive cytokine RANKL and its inhibitor osteoprotegerin (OPG), in osteoclast and osteoblast numbers, and in bone mass. We assumed that osteocytes surrounding a microfracture produce RANKL, which attracted osteoclasts. OPG and RANKL were produced by osteoblasts and diffused through bone, RANKL was eliminated by binding to OPG and RANK. Osteoblasts were coupled to osteoclasts through paracrine factors. The evolution of the BMU arising from this model was studied using numerical simulations. Our model recapitulated the spatio-temporal dynamics observed in vivo in a cross-section of bone. In response to a RANKL field, osteoclasts moved as a well-confined cutting cone. The coupling of osteoclasts to osteoblasts allowed for sufficient recruitment of osteoblasts to the resorbed surfaces. The RANKL field was the highest at the microfracture in front of the BMU, whereas the OPG field peaked at the back of the BMU, resulting in the formation of a RANKL/OPG gradient, which strongly affected the rate of BMU progression and its size. Thus, the spatial organization of a BMU provides important constraints on the roles of RANKL and OPG as well as possibly other regulators in determining the outcome of remodeling in the BMU.

A review of recent developments in mathematical modeling of bone remodeling

In this article, we summarize the developments in the mathematical modeling of the mechanics of bone and related biological phenomena. We will devote special attention to the results of the last 10–15 years, although we will cover some relevant classical work to better frame the more recent researches. We will propose a division of the literature based on the main aim of the model (mechanical/biomathematical) and the type of biological phenomena considered (stimulus, growth, cell population dynamics). Finally, we will suggest some possible directions for future investigations.

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

Bone remodeling is an essential physiological process 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.

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.

Dynamic Biochemical and Cellular Models of Bone Physiology: Integrating Remodeling Processes, Tumor Growth, and Therapy

Bone is an activeCoelho, Rui M.  Neto, Joana P. connective tissue composed of different types of cells. The dynamic behavior of bone remodeling processes is typically represented throughValério, Duarte differential equations, which represent the physiological phenomena occurring in this organ. These models takeVinga, Susana into account the tight biochemical regulation between osteoclasts and osteoblasts and have also been enriched with variables and parameters related to bone pathologies and treatment. This chapter reviews some of the more recent models describing bone physiology, focusing on those that include the main cellular processes, along the biochemical control, and also the pharmacokinetic/pharmacodynamic (PK/PD) of the most common treatments for diseases such as cancer. These models are then compared in terms of the simulations obtained and, finally, some highlights on integrating them with the biomechanical component of the system which will be given. These models are expected to provide a valuable insight into this complex system and to support the development of clinical decision systems for bone pathologies.

Dynamics of Bone Cell Interactions and Differential Responses to PTH and Antibody-Based Therapies

We propose a mathematical model describing the dynamics of osteoblasts and osteoclasts in bone remodeling. The goal of this work is to develop an integrated modeling framework for bone remodeling and bone cell signaling dynamics that could be used to explore qualitatively combination treatments for osteoporosis in humans. The model has been calibrated using 57 checks from the literature. Specific global optimization methods based on qualitative objectives have been developed to perform the model calibration. We also added pharmacokinetics representations of three drugs to the model, which are teriparatide (PTH(1–34)), denosumab (a RANKL antibody) and romosozumab (a sclerostin antibody), achieving excellent goodness-of-fit of human clinical data. The model reproduces the paradoxical effects of PTH on the bone mass, where continuous administration of PTH results in bone loss but intermittent administration of PTH leads to bone gain, thus proposing an explanation of this phenomenon. We used the model to simulate different categories of osteoporosis. The main attributes of each disease are qualitatively well captured by the model, for example changes in bone turnover in the disease states. We explored dosing regimens for each disease based on the combination of denosumab and romosozumab, identifying adequate ratios and doses of both drugs for subpopulations of patients in function of categories of osteoporosis and the degree of severity of the disease.

Computational Modeling of Interactions between Multiple Myeloma and the Bone Microenvironment

Multiple Myeloma (MM) is a B-cell malignancy that is characterized by osteolytic bone lesions. It has been postulated that positive feedback loops in the interactions between MM cells and the bone microenvironment form reinforcing ‘vicious cycles’, resulting in more bone resorption and MM cell population growth in the bone microenvironment. Despite many identified MM-bone interactions, the combined effect of these interactions and their relative importance are unknown. In this paper, we develop a computational model of MM-bone interactions and clarify whether the intercellular signaling mechanisms implemented in this model appropriately drive MM disease progression. This new computational model is based on the previous bone remodeling model of Pivonka et al. [1], and explicitly considers IL-6 and MM-BMSC (bone marrow stromal cell) adhesion related pathways, leading to formation of two positive feedback cycles in this model. The progression of MM disease is simulated numerically, from normal bone physiology to a well established MM disease state. Our simulations are consistent with known behaviors and data reported for both normal bone physiology and for MM disease. The model results suggest that the two positive feedback cycles identified for this model are sufficient to jointly drive the MM disease progression. Furthermore, quantitative analysis performed on the two positive feedback cycles clarifies the relative importance of the two positive feedback cycles, and identifies the dominant processes that govern the behavior of the two positive feedback cycles. Using our proposed quantitative criteria, we identify which of the positive feedback cycles in this model may be considered to be ‘vicious cycles’. Finally, key points at which to block the positive feedback cycles in MM-bone interactions are identified, suggesting potential drug targets.

Complex Dynamics of Osteoclast Formation and Death in Long-Term Cultures

BackgroundOsteoclasts, cells responsible for bone resorption, contribute to the development of degenerative, metabolic and neoplastic bone diseases, which are often characterized by persistent changes in bone microenvironment. We aimed to investigate the dynamics of osteoclast formation and death in cultures that considerably exceeded the length of standard protocol and to design a mathematical model describing osteoclastogenesis.Methodology/Principal FindingsRAW 264.7 monocytic cells fuse to form multinucleated osteoclasts upon treatment with pro-resorptive cytokine RANKL. We have found that in long-term experiments (15–26 days), the dynamics of changes in osteoclast numbers was remarkably complex and qualitatively variable in different experiments. Whereas 19 of 46 experiments exhibited single peak of osteoclast formation, in 27 experiments we observed development of successive waves of osteoclast formation and death. Periodic changes in osteoclast numbers were confirmed in long-term cultures of mouse bone marrow cells treated with M-CSF and RANKL. Because the dynamics of changes in osteoclast numbers was found to be largely independent of monocytes, a two-species model of ordinary differential equations describing the changes in osteoclasts and monocytes was ineffective in recapitulating the oscillations in osteoclast numbers. Following experimental observation that medium collected from mature osteoclasts inhibited osteoclastogenesis in fresh cultures, we introduced a third variable, factor f, to describe osteoclast-derived inhibitor. This model allowed us to simulate the oscillatory changes in osteoclasts, which were coupled to oscillatory changes in the factor f, whereas monocytes changed exponentially. Importantly, to achieve the experimentally observed oscillations with increasing amplitude, we also had to assume that osteoclast presence stimulates osteoclast formation.Conclusions/SignificanceThis study identifies the critical role for osteoclast autocrine regulation in controlling long-term dynamic of osteoclast formation and death and describes the complementary roles for negative and positive feedback mediators in determining the sharp dynamics of activation and inactivation of osteoclasts.

<i>In silico</i> biology of bone modelling and remodelling: regeneration

Bone regeneration is the process whereby bone is able to (scarlessly) repair itself from trauma, such as fractures or implant placement. Despite extensive experimental research, many of the mechanisms involved still remain to be elucidated. Over the last decade, many mathematical models have been established to investigate the regeneration process in silico . The first models considered only the influence of the mechanical environment as a regulator of the healing process. These models were followed by the development of bioregulatory models where mechanics was neglected and regeneration was regulated only by biological stimuli such as growth factors. The most recent mathematical models couple the influences of both biological and mechanical stimuli. Examples are given to illustrate the added value of mathematical regeneration research, specifically in the in silico design of treatment strategies for non-unions. Drawbacks of the current continuum-type models, together with possible solutions in extending the models towards other time and length scales are discussed. Finally, the demands for dedicated and more quantitative experimental research are presented.

Optimizing the dynamics of bone turnover with genetic algorithm

Control systems and the modeling strategies are not only limited to engineering problems. These approaches can be used in the field of bio-mathematics as well and modern studies have promoted this approach to a great extent. The computational modeling and simulation of bone metastasis is painful yet critical after cancer invades the body. This vicious cycle is complex, and several research centers worldwide are devoted to understanding the dynamics and setting up a treatment strategy for this life-threatening behavior of cancer. Cancerous cells activation and the corresponding process of metastasis is reported to boost during the periodic waves of COVID-19, due to the inflammatory nature of the infection associated with SARS-2 and its variants. The bone cells are comprised of two types of cells responsible for bone formation and resorption. The computational framework of such cells, in spatial form, can help the researchers forecast the bone dynamics in a robust manner where the impact of cancer is incorporated into the computational model as a source of perturbation. A series of computational models are presented to explore the complex behavior of bone metastasis with COVID-19 induced infection. The finite difference algorithm is used to simulate the nonlinear computational model. The results obtained are in close agreement with the experimental findings. The computational results can help explore the vicious cycle’s fate and help set up control strategies through drug therapies.

Osteoimmunology: A Current Update of the Interplay Between Bone and the Immune System

Immunology, already a discipline in its own right, has become a major part of many different medical fields. However, its relationship to orthopedics and trauma surgery has unfortunately, and perhaps unjustly, been developing rather slowly. Discover-ies in recent years have emphasized the immense breadth of communication and connection between both systems and, importantly, the highly promising therapeu-tic opportunities.Recent discoveries of factors originally assigned to the immune system have now also been shown to have a significant impact on bone health and disease, which has greatly changed how we approach treatment of bone pathologies. In case of bone fracture, immune cells, especially macrophages, are present through-out the whole healing process, assure defense against pathogens and discharge a complex variety of effectors to regulate bone modelling. In rheumatoid arthritis and osteoporosis, the immune system contributes to the formation of the pathological and chronic conditions. Fascinatingly, prosthesis failure is not at all solely a me-chanical problem of improper strain but works in conjunction with an active contri-bution of the immune system as a reaction to irritant debris from material wear.Unravelling conjoined mechanisms of the immune and osseous systems heralds therapeutic possibilities for ailments of both. Contemplation of the bone as merely an unchanging support pillar is outdated and obsolete. Instead it is manda-tory that this highly diverse network be incorporated in our understanding of the immune system and hematopoiesis.

A mechano-chemo-biological model for bone remodeling with a new mechano-chemo-transduction approach

Bone remodeling is a fundamental biological process that develops in bone tissue along its whole lifetime. It refers to a continuous bone transformation with new bone formation and old bone resorption that changes the internal microstructure and composition of the tissue. The main objectives of bone remodeling are: repair of the internal microcracks; adaptation of the macroscopic stiffness and strength to the actual changing mechanical demands; and control of the calcium homeostasis. Understanding this process and predicting its evolution is critical to reduce the effects of long-term disuse as happens during periods of reduced mobility. It is also important in the design of bone implants to avoid long-term stress shielding. Many mathematical models have been proposed from the earliest purely phenomenological to the latest that include biological knowledge. However, there still exists a lack of connection between the mechanical driving force and the biochemical and cell processes it triggers. Here, and following previous works that model independently the mechanobiological and biochemical processes in bone remodeling, we present a more complete model, useful for both cortical and trabecular bone, that uses a new mechanotransduction approach based on the effect of strains onto the bonding–unbonding rate of RANK/RANKL/OPG receptor–ligand reactions. We compare the results of this model with previous ones, showing a good agreement in similar conditions. We also apply it to realistic situations such as a femoral bone after implantation of a hip prosthesis, getting similar results to the clinical ones in the final bone density distribution. Finally, we extend this approach to the anisotropic case, getting not only the mean density, but also the directional homogenization of the microstructure. This biochemical approach permits, not only to predict the bone evolution under changes in the mechanical loads, but also, to consider anabolic and catabolic drugs to control bone density, such as those used in osteoporosis.