Advances in Continuous and Discrete Models - In this paper, we propose an impulsive mathematical model of bone formation and resorption accounting for the number of active osteoclastic cells, bone...
Parathyroid hormone temporal effects on bone formation and resorption
Parathyroid hormone (PTH) paradoxically causes net bone loss (resorption) when administered in a continuous fashion, and net bone formation (deposition) when administered intermittently. Currently no pharmacological formulations are available to promote bone formation, as needed for the treatment of osteoporosis. The paradoxical behavior of PTH confuses endocrinologists, thus, a model bone resorption or deposition dependent on the timing of PTH administration would de-mystify this behavior and provide the basis for logical drug formulation. We developed a mathematical model that accounts for net bone loss with continuous PTH administration and net bone formation with intermittent PTH administration, based on the differential effects of PTH on the osteoblastic and osteoclastic populations of cells.
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.
Mathematical Model of Bone Remodeling Captures the Antiresorptive and Anabolic Actions of Various Therapies
A better understanding of the molecular pathways regulating the bone remodeling process should help in the development of new antiresorptive regulators and anabolic regulators, that is, regulators of bone resorption and of bone formation. Understanding the mechanisms by which parathyroid hormone (PTH) influences bone formation and how it switches from anabolic to catabolic action is important for treating osteoporosis (Poole and Reeve in Curr Opin Pharmacol 5:612–617, 2005). In this paper we describe a mathematical model of bone remodeling that incorporates, extends, and integrates several models of particular aspects of this biochemical system (Cabal et al. in J Bone Miner Res 28(8):1830–1836, 2013; Lemaire et al. in J Theor Biol 229:293–309, 2004; Peterson and Riggs in Bone 46:49–63, 2010; Raposo et al. in J Clin Endocrinol Metab 87(9):4330–4340, 2002; Ross et al. in J Disc Cont Dyn Sys Series B 17(6):2185–2200, 2012). We plan to use this model as a bone homeostasis platform to develop anabolic and antiresorptive compounds. The model will allow us to test hypotheses about the dynamics of compounds and to test the potential benefits of combination therapies. At the core of the model is the idealized account of osteoclast and osteoblast signaling given by Lemaire et al. (J Theor Biol 229:293–309, 2004). We have relaxed some of their assumptions about the roles of osteoprotegerin, transforming growth factor $$\upbeta $$, and receptor activator of nuclear factor $$\upkappa $$B ligand; we have devised more detailed models of the interactions of these species. We have incorporated a model of the effect of calcium sensing receptor antagonists on remodeling (Cabal et al. in J Bone Miner Res 28(8):1830–1836, 2013). We have also incorporated a basic model of the effects of vitamin D on calcium homeostasis. We have included a simple model of the mechanism proposed by Bellido et al. (2003), Ross et al. (J Disc Cont Dyn Sys Series B 17(6):2185–2200, 2012), of the influence of PTH on osteoblast apoptosis, a mechanism that accounts for the anabolic response to pulsatile PTH administration. Finally, we have devised a simple model of the administration and effects of bisphosphonates. The biomarkers in the model are procollagen type 1 amino-terminal propeptide and C-terminal telopeptide. Bone mineral density is the model’s principal endpoint.
Analytic approach to explore dynamical osteoporotic bone turnover
The dynamics of the osteoporotic bone turnover is studied in this paper with the aid of stability analysis of the associated mathematical model. Osteoporosis, which is a common bone disorder, is studied in this papper in detail with an emphasis on the relative threshold values. We examine the expository signaling among the bone cells named osteoclast and osteoblast. Main functioning of osteoblasts is bone formation, whereas osteoclasts are bone removal cells.
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 Model of Paracrine Interactions between Osteoclasts and Osteoblasts Predicts Anabolic Action of Parathyroid Hormone on Bone
Abstract. To restore falling plasma calcium levels, PTH promotes calcium liberation from bone. PTH targets bone-forming cells, osteoblasts, to increase exp
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 signaling and PTH treatments of osteoporosis
In this paper we analyze and generalize the dynamical system introduced by Lemaire and co-workers [14] as a model of cell signaling in bone remodeling. We show that for large classes of parameter values, including the physically-realistic baseline values, the system has a unique physically relevant equilibrium which is a global attractor. We generalize that model, minimally, to incorporate a mechanism by which parathyroid hormone (PTH) retards osteoblast apoptosis. We show that with this mechanism, which is a simplified version of that proposed by Bellido and co-workers [4], the model exhibits a well-known phenomenon that has puzzled researchers: the system responds catabolically to the continuous administration of PTH andally to appropriately pulsed administration of PTH.
Application of a mechanism-based disease systems model for osteoporosis to clinical data
A recently proposed mechanism-based disease systems model for osteoporosis (Schmidt et al., J Pharmacokinet Pharmacodyn 38:873–900, 2011) was applied to clinical data from post-menopausal women (n = 767) receiving various doses of the selective estrogen receptor modulator tibolone. Plasma bone-specific alkaline phosphatase activity and urinary N-telopeptide were used as biomarkers reflecting the activity of osteoblasts (bone forming cells) and osteoclasts (bone removing cells), respectively. These data were analyzed in conjunction with data on osteocalcin and on bone mineral density (BMD) (both lumbar spine and total hip), which reflect the activity of both cell types. While the dynamics of bone turnover markers changes rapidly, closely following changes in the activity of bone cells, changes in BMD are slower and have their own dynamics. Application of the mechanism-based disease systems model to the clinical data allowed for an adequate description of the data and yielded parameter estimates that are consistent with physiological values reported in the literature (Lemaire et al., J Theor Biol 229:293–309, 2004). The fitted model enabled characterization of (i) the critical time scales involved in disease progression, (ii) the dynamics of the system during onset and offset of the therapeutic intervention, and (iii) the distinction between responders and low-responders to tibolone treatment.
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.
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.
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.
Studying Bone Remodelling and Tumour Growth for Therapy Predictive Control
Bone remodelling consists of cycles of bone resorption and formation executed mainly by osteoclasts and osteoblasts. Healthy bone remodelling is disrupted by diseases such as Multiple Myeloma and bone metastatic diseases. In this paper, a simple mathematical model with differential equations, which takes into account the evolution of osteoclasts, osteoblasts, bone mass and bone metastasis growth, is improved with a pharmacokinetic and pharmacodynamic (PK/PD) scheme of the drugs denosumab, bisphosphonates, proteasome inhibitors and paclitaxel. The major novelty is the inclusion of drug resistance phenomena, which resulted in two variations of the model, corresponding to different paradigms of the origin and development of the tumourous cell resistance condition. These models are then used as basis for an optimization of the drug dose applied, paving the way for personalized medicine. A Nonlinear Model Predictive Control scheme is used, which takes advantage of the convenient properties of a suggested adaptive and democratic variant of Particle Swarm Optimization. Drug prescriptions obtained in this way provide useful insights into dose administration strategies. They also show how results may change depending on which of the two very different paradigms of drug resistance is used to model the behaviour of the tumour.
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.
Modeling the Macrophage-Mediated Inflammation Involved in the Bone Fracture Healing Process
A new mathematical model is presented to study the effects of macrophages on the bone fracture healing process. The model consists of a system of nonlinear ordinary differential equations that represents the interactions among classically and alternatively activated macrophages, mesenchymal stem cells, osteoblasts, and pro- and anti-inflammatory cytokines. A qualitative analysis of the model is performed to determine the equilibria and their corresponding stability properties. Numerical simulations are also presented to support the theoretical results, and to monitor the evolution of a broken bone for different types of fractures under various medical interventions. The model can be used to guide clinical experiments and to explore possible medical treatments that accelerate the bone fracture healing process, either by surgical interventions or drug administrations.
Integrative modeling of drug therapy and the bone turnover
Denosumab has been successfully used for the postmenopausal osteoporosis treatment. This research is focused on the computational analysis of the effect of denosumab on bone remodeling.
