Application of a mechanism-based disease systems model for osteoporosis to clinical data

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.

Integrated model for denosumab and ibandronate pharmacodynamics in postmenopausal women

ABSTRACT This study aims to characterize the pharmacodynamic properties of denosumab, a RANK ligand inhibitor, and ibandronate, a bisphosphonate, using an integrated bone homeostasis model in postmenopausal women. Mean temporal profiles of denosumab, serum and urine N‐telopeptide (sNTX, uNTX), lumbar spine bone mineral density (BMD) following denosumab administration, and urine C‐telopeptide (uCTX) and lumbar spine BMD upon ibandronate administration were extracted from the literature. A mechanistic model was developed that integrates denosumab pharmacokinetics with binding to RANK ligand and ibandronate inhibition of osteoclast precursor differentiation to active osteoclasts (AOC). Biomarker concentrations were linked to the AOC pool. The BMD was characterized by a turnover model with stimulation of bone formation and degradation by AOB (active osteoblasts) and AOC pools. The estimated basal sNTX, uNTX and uCTX concentrations were 7.24 n m , 14.4 nmol/mmolCr and 31 µg/mmolCr. The BMD degradation rate was 0.00161 day ‐1 with stimulation constants associated with AOB and AOC of 1214 and 790 p m ‐1 . The plasma ibandronate concentration producing 50% of maximum inhibition of osteoclast differentiation was 522 ng/l. The integrated model, which incorporates multiple pathways of therapeutic intervention, quantitatively describes changes in clinical biomarkers of bone turnover and BMD after denosumab and ibandronate exposures in postmenopausal women. Copyright © 2011 John Wiley & Sons, Ltd.

Bone Remodeling in Post-menopausal Osteoporosis

Bone mass in the skeleton is dependent on the coordinated activities of bone-forming osteoblasts and bone-resorbing osteoclasts in discrete bone multi-cellular units. Remodeling of bone in these units is important not only for maintaining bone mass, but also to repair microdamage, to prevent accumulation of too much old bone, and for mineral homeostasis. The activities of osteoblasts and osteoclasts are controlled by a variety of hormones and cytokines, as well as by mechanical loading. Most importantly, sex hormones are very crucial for keeping bone mass in balance, and the lack of either estrogen or testosterone leads to decreased bone mass and increased risk for osteoporosis. The prevalence of osteoporotic fractures is increasing dramatically in the Western part of the world and is a major health problem in many countries. In the present review, the cellular and molecular mechanisms controlling bone remodeling and the influence of sex hormones on these processes are summarized. In a separate paper in this issue, the pathogenesis of post-menopausal osteoporosis will be compared with that of inflammation-induced bone remodeling, including the evidence for and against the hypothesis that concomitant post-menopausal osteoporotic disease influences the progression of periodontal disease.

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.

Predicting Nonlinear Changes in Bone Mineral Density Over Time Using a Multiscale Systems Pharmacology Model

A mathematical model component that extends an existing physiologically based multiscale systems pharmacology model (MSPM) of calcium and bone homeostasis was developed, enabling prediction of nonlinear changes in lumbar spine bone mineral density (LSBMD). Data for denosumab, a monoclonal antibody osteoporosis treatment, dosed at several levels and regimens, was used for fitting the BMD component. Bone marker and LSBMD data extracted from the literature described on/off‐treatment effects of denosumab over 48 months [Miller, P.D. et al . Effect of denosumab on bone density and turnover in postmenopausal women with low bone mass after long‐term continued, discontinued, and restarting of therapy: a randomized blinded phase 2 clinical trial. Bone 43, 222–229 (2008)]. An indirect model linking bone markers to LSBMD was embedded in the existing MSPM, reasonably predicting nonlinear increases in LSBMD during treatment (24 months); LSBMD declines following discontinuation and increases upon treatment reinstitution. This study demonstrates the utility of MSPM extension to describe a phenomena of interest not originally in a model, and the ability of this updated MSPM to predict nonlinear longitudinal changes in the clinically relevant endpoint, LSBMD, with denosumab treatment. CPT: Pharmacometrics & Systems Pharmacology (2012) 1, e14; doi: 10.1038/psp.2012.15 ; advance online publication 14 November 2012

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.

Critical Role of Estrogens on Bone Homeostasis in Both Male and Female: From Physiology to Medical Implications

Bone is a multi-skilled tissue, protecting major organs, regulating calcium phosphate balance and producing hormones. Its development during childhood determines height and stature as well as resistance against fracture in advanced age. Estrogens are key regulators of bone turnover in both females and males. These hormones play a major role in longitudinal and width growth throughout puberty as well as in the regulation of bone turnover. In women, estrogen deficiency is one of the major causes of postmenopausal osteoporosis. In this review, we will summarize the main clinical and experimental studies reporting the effects of estrogens not only in females but also in males, during different life stages. Effects of estrogens on bone involve either Estrogen Receptor (ER)α or ERβ depending on the type of bone (femur, vertebrae, tibia, mandible), the compartment (trabecular or cortical), cell types involved (osteoclasts, osteoblasts and osteocytes) and sex. Finally, we will discuss new ongoing strategies to increase the benefit/risk ratio of the hormonal treatment of menopause.

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.

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.

The Development of Molecular Biology of Osteoporosis

Osteoporosis is one of the major bone disorders that affects both women and men, and causes bone deterioration and bone strength. Bone remodeling maintains bone mass and mineral homeostasis through the balanced action of osteoblasts and osteoclasts, which are responsible for bone formation and bone resorption, respectively. The imbalance in bone remodeling is known to be the main cause of osteoporosis. The imbalance can be the result of the action of various molecules produced by one bone cell that acts on other bone cells and influence cell activity. The understanding of the effect of these molecules on bone can help identify new targets and therapeutics to prevent and treat bone disorders. In this article, we have focused on molecules that are produced by osteoblasts, osteocytes, and osteoclasts and their mechanism of action on these cells. We have also summarized the different pharmacological osteoporosis treatments that target different molecular aspects of these bone cells to minimize osteoporosis.

Coping with time scales in disease systems analysis: application to bone remodeling

In this study we demonstrate the added value of mathematical model reduction for characterizing complex dynamic systems using bone remodeling as an example. We show that for the given parameter values, the mechanistic RANK-RANKL-OPG pathway model proposed by Lemaire et al. (J Theor Biol 229:293–309, 2004) can be reduced to a simpler model, which can describe the dynamics of the full Lemaire model to very good approximation. The response of both models to changes in the underlying physiology and therapeutic interventions was evaluated in four physiologically meaningful scenarios: (i) estrogen deficiency/estrogen replacement therapy, (ii) Vitamin D deficiency, (iii) ageing, and (iv) chronic glucocorticoid treatment and its cessation. It was found that on the time scale of disease progression and therapeutic intervention, the models showed negligible differences in their dynamic properties and were both suitable for characterizing the impact of estrogen deficiency and estrogen replacement therapy, Vitamin D deficiency, ageing, and chronic glucocorticoid treatment and its cessation on bone forming (osteoblasts) and bone resorbing (osteoclasts) cells. It was also demonstrated how the simpler model could help in elucidating qualitative properties of the observed dynamics, such as the absence of overshoot and rebound, and the different dynamics of onset and washout.

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.

Application of a Systems Pharmacology‐Based Placebo Population Model to Analyze Long‐Term Data of Postmenopausal Osteoporosis

Osteoporosis is a progressive bone disease characterized by decreased bone mass resulting in increased fracture risk. The objective of this investigation was to test whether a recently developed dise...

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.

Simulated Interventions to Ameliorate Age-Related Bone Loss Indicate the Importance of Timing

Bone remodeling is the continuous process of bone resorption by osteoclasts and bone formation by osteoblasts, in order to maintain homeostasis. The activity of osteoclasts and osteoblasts is regulated by a network of signaling pathways, including Wnt, parathyroid hormone (PTH), RANKL/OPG and TGF-β, in response to stimuli such as mechanical loading. During aging there is a gradual loss of bone mass due to dysregulation of signaling pathways. This may be due to a decline in physical activity with age and/or changes in hormones and other signaling molecules. In particular, hormones such as PTH have a circadian rhythm which may be disrupted in aging. Due to the complexity of the molecular and cellular networks involved in bone remodeling, several mathematical models have been proposed to aid understanding of the processes involved. However, to date there are no models which explicitly consider the effects of mechanical loading, the circadian rhythm of PTH and the dynamics of signaling molecules on bone remodeling. Therefore, we have constructed a network model of the system using a modular approach which will allow further modifications as required in future research. The model was used to simulate the effects of mechanical loading and also the effects of different interventions such as continuous or intermittent administration of PTH. Our model predicts that the absence of regular mechanical loading and/or an impaired PTH circadian rhythm leads to a gradual decrease in bone mass over time which can be restored by simulated interventions and that the effectiveness of some interventions may depend on their timing.

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.