







It is concluded that intermittent administration of calcitonin is an effective therapy for the treatment of osteoporosis.
Forecasting the critical role of intermittent therapies for the control of bone resorption
Osteoporosis is a chronic metabolic disease characterized by an imbalance of bone resorption and formation, leading to bone fragility and increased susceptibility to fracture. Parathyroid hormone is approved therapy for the treatment of osteoporosis.

Modeling osteoporosis to design and optimize pharmacological therapies comprising multiple drug types
For the treatment of postmenopausal osteoporosis, several drug classes with different mechanisms of action are available. Since only a limited set of dosing regimens and drug combinations can be tested in clinical trials, it is currently unclear whether common medication strategies achieve optimal bone mineral density gains or are outperformed by alternative dosing schemes and combination therapies that have not been explored so far. Here, we develop a mathematical framework of drug interventions for postmenopausal osteoporosis that unifies fundamental mechanisms of bone remodeling and the mechanisms of action of four drug classes: bisphosphonates, parathyroid hormone analogs, sclerostin inhibitors, and receptor activator of NF-κB ligand inhibitors. Using data from several clinical trials, we calibrate and validate the model, demonstrating its predictive capacity for complex medication scenarios, including sequential and parallel drug combinations. Via simulations, we reveal that there is a large potential to improve gains in bone mineral density by exploiting synergistic interactions between different drug classes, without increasing the total amount of drug administered.
Biology of Bone Tissue: Structure, Function, and Factors That Influence Bone Cells
Bone tissue is continuously remodeled through the concerted actions of bone cells, which include bone resorption by osteoclasts and bone formation by osteoblasts, whereas osteocytes act as mechanosensors and orchestrators of the bone remodeling process. This process is under the control of local (e.g., growth factors and cytokines) and systemic (e.g., calcitonin and estrogens) factors that all together contribute for bone homeostasis. An imbalance between bone resorption and formation can result in bone diseases including osteoporosis. Recently, it has been recognized that, during bone remodeling, there are an intricate communication among bone cells. For instance, the coupling from bone resorption to bone formation is achieved by interaction between osteoclasts and osteoblasts. Moreover, osteocytes produce factors that influence osteoblast and osteoclast activities, whereas osteocyte apoptosis is followed by osteoclastic bone resorption. The increasing knowledge about the structure and functions of bone cells contributed to a better understanding of bone biology. It has been suggested that there is a complex communication between bone cells and other organs, indicating the dynamic nature of bone tissue. In this review, we discuss the current data about the structure and functions of bone cells and the factors that influence bone remodeling.

Estrogen therapy for osteoporosis in the modern era
Menopause predisposes women to osteoporosis due to declining estrogen levels. This results in a decrease in bone mineral density (BMD) and an increase in fractures. Osteoporotic fractures lead to substantial morbidity and mortality, and are considered one of the largest public health priorities by the World Health Organization (WHO). It is therefore essential for menopausal women to receive appropriate guidance for the prevention and management of osteoporosis. The Women’s Health Initiative (WHI) randomized controlled trial first proved hormonal therapy (HT) reduces the incidence of all osteoporosis-related fractures in postmenopausal women. However, the study concluded that the adverse effects outweighed the potential benefits on bone, leading to a significant decrease in HT use for menopausal symptoms. Additionally, HT was not used as first-line therapy for osteoporosis and fractures. Subsequent studies have challenged these initial conclusions and have shown significant efficacy of HT in various doses, durations, regimens, and routes of administration. These studies support that HT improves BMD and reduces fracture risk in women with and without osteoporosis. Furthermore, the studies suggest that low-dose and transdermal HT are less likely associated with the adverse effects of breast cancer, endometrial hyperplasia, coronary artery disease (CAD), and venous thromboembolism (VTE) previously observed in standard-dose oral HT regimens. Given the need for estrogen in menopausal women and evidence supporting the cost effectiveness, safety, and efficacy of HT, we propose that HT should be considered for the primary prevention and treatment of osteoporosis in appropriate candidates. HT should be individualized and the once “lowest dose for shortest period of time” concept should no longer be used. This review will focus on the prior and current studies for various HT formulations used for the prevention and treatment of osteoporosis, exploring the safety profile of low-dose and transdermal HT that have been shown to be safer than oral standard-dose HT.
Parathyroid Hormone and Teriparatide for the Treatment of Osteoporosis: A Review of the Evidence and Suggested Guidelines for Its Use
All therapies currently recommended for the management of osteoporosis act mainly to inhibit bone resorption and reduce bone remodeling. PTH and its analog

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.
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.

Current Status of the Diagnosis and Management of Osteoporosis
Osteoporosis has been defined as the silent disease of the 21st century, becoming a public health risk due to its severity, chronicity and progression and affecting mainly postmenopausal women and older adults. Osteoporosis is characterized by an imbalance between bone resorption and bone production. It is diagnosed through different methods such as bone densitometry and dual X-rays. The treatment of this pathology focuses on different aspects. On the one hand, pharmacological treatments are characterized by the use of anti-resorptive drugs, as well as emerging regenerative medicine treatments such as cell therapies and the use of bioactive hydrogels. On the other hand, non-pharmacological treatments are associated with lifestyle habits that should be incorporated, such as physical activity, diet and the cessation of harmful habits such as a high consumption of alcohol or smoking. This review seeks to provide an overview of the theoretical basis in relation to bone biology, the existing methods for diagnosis and the treatments of osteoporosis, including the development of new strategies.

Bisphosphonates: Mechanism of Action and Role in Clinical Practice
Bisphosphonates are primary agents in the current pharmacological arsenal against osteoclast-mediated bone loss due to osteoporosis, Paget disease of bone, malignancies metastatic to bone, multiple myeloma, and hypercalcemia of malignancy. In addition to currently approved uses, bisphosphonates are commonly prescribed for prevention and treatment of a variety of other skeletal conditions, such as low bone density and osteogenesis imperfecta. However, the recent recognition that bisphosphonate use is associated with pathologic conditions including osteonecrosis of the jaw has sharpened the level of scrutiny of the current widespread use of bisphosphonate therapy.

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.

Clinical Data for Parametrization of In Silico Bone Models Incorporating Cell-Cytokine Dynamics: A Systematic Review of Literature
In silico simulations aim to provide fast, inexpensive, and ethical alternatives to years of costly experimentation on animals and humans for studying bone remodeling, its deregulation during osteoporosis and the effect of therapeutics. Within the varied spectrum of in silico modeling techniques, bone cell population dynamics and agent-based multiphysics simulations have recently emerged as useful tools to simulate the effect of specific signaling pathways. In these models, parameters for cell and cytokine behavior are set based on experimental values found in literature; however, their use is currently limited by the lack of clinical in vivo data on cell numbers and their behavior as well as cytokine concentrations, diffusion, decay and reaction rates. Further, the settings used for these parameters vary across research groups, prohibiting effective cross-comparisons. This review summarizes and evaluates the clinical trial literature that can serve as input or validation for in silico models of bone remodeling incorporating cells and cytokine dynamics in post-menopausal women in treatment, and control scenarios. The GRADE system was used to determine the level of confidence in the reported data, and areas lacking in reported measures such as binding site occupancy, reaction rates and cell proliferation, differentiation and apoptosis rates were highlighted as targets for further research. We propose a consensus for the range of values that can be used for the cell and cytokine settings related to the RANKL-RANK-OPG, TGF-β and sclerostin pathways and a Levels of Evidence-based method to estimate parameters missing from clinical trial literature.

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.

Disease Systems Analysis of Bone Mineral Density and Bone Turnover Markers in Response to Alendronate, Placebo, and Washout in Postmenopausal Women
A previously established mechanism‐based disease systems model for osteoporosis that is based on a mathematically reduced version of a model describing the interactions between osteoclast (bone removing) and osteoblast (bone forming) cells in bone remodeling has been applied to clinical data from women ( n = 1,379) receiving different doses and treatment regimens of alendronate, placebo, and washout. The changes in the biomarkers, plasma bone‐specific alkaline phosphatase activity (BSAP), urinary N‐telopeptide (NTX), lumbar spine bone mineral density (BMD), and total hip BMD, were linked to the underlying mechanistic core of the model. The final model gave an accurate description of all four biomarkers for the different treatments. Simulations were used to visualize the dynamics of the underlying network and the natural disease progression upon alendronate treatment and discontinuation. These results complement the previous applications of this mechanism‐based disease systems model to data from various treatments for osteoporosis.

An impulsive mathematical model of bone formation and resorption: effects of parathyroid hormone, calcitonin and impulsive estrogen supplement
In this paper, we propose an impulsive mathematical model of bone formation and resorption accounting for the number of active osteoclastic cells, bone resorbing cells, and the number of active osteoblastic cells, bone forming cells, based on the effects of parathyroid hormone and calcitonin with impulsive estrogen supplement. The model is then analyzed theoretically in terms of its stability, permanence and oscillatory behavior. The conditions on the model parameters, for which the desirable behaviors of the solution of the system can be expected, are derived. Numerical simulations are also carried out in order to support our theoretical predictions. The results indicate that the frequency and dosage of the estrogen supplements are important since the behavior of the solution of the system depends on the frequency and dosage of the estrogen supplements.

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.

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.
