The purpose of this study is to couple a cellular bone homeostasis model with the pharmacokinetics (PK) and mechanism of action of denosumab, an inhibitor of receptor activator of nuclear factor-κB ligand, to characterize the time course of serum N-telopeptide (NTX), a bone resorption biomarker, following single escalating doses in multiple myeloma (MM) patients. Mean PK and median serum NTX temporal profiles were extracted from a previously conducted randomized, double-blind, double-dummy, active-controlled, multicenter study including 25 MM patients receiving escalating denosumab doses.
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.
Long-term use of denosumab and its association with skeletal-related events and osteonecrosis of the jaw
Denosumab, an inhibitor of receptor activator of nuclear factor kappa-B ligand, reduces skeletal-related events (SREs) and is approved for solid tumors with bone metastases. We studied long-term denosumab efficacy and safety because real-world data is scarce. This single-arm, single-center retrospective study included denosumab-treated breast cancer patients with bone metastases. Kaplan–Meier survival curves assessed exposure, SREs, osteonecrosis of the jaw (ONJ), and death. 132 patients were enrolled. The median denosumab exposure was 28.3 months (range 1.0–84.9). In the first year, 11.1% experienced SREs. This increased to 18.6% in the second, 21% in the third, and 35.1% in the fourth year and beyond. The median time to first on-study SRE has not been reached. 10 denosumab users (7.6%) developed ONJ. ONJ incidence was 0.9% in the first year, 6.2% in the second, 13.6% in the third, and 16.2% in subsequent years. The median time to first on-study ONJ has not been reached yet. Seven patients resumed denosumab after careful management of ONJ. Our data suggest that long-term treatment with denosumab may further prevent or postpone SREs at the cost of an increased risk of ONJ. The majority of patients who resumed denosumab did not experience a recurrence of ONJ.
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
Ten‐Year Simulation of the Effects of Denosumab on Bone Remodeling in Human Biopsies
ABSTRACT Postmenopausal osteoporosis is a disease manifesting in degradation of bone mass and microarchitecture, leading to weakening and increased risk of fracture. Clinical trials are an essential tool for evaluating new treatments and may provide further mechanistic understanding of their effects in vivo. However, the histomorphometry from clinical trials is limited to 2D images and reflects single time points. Biochemical markers of bone turnover give global insight into a drug's action, but not the local dynamics of the bone remodeling process and the cells involved. Additionally, comparative trials necessitate separate treatment groups, meaning only aggregated measures can be compared. In this study, in silico modeling based on histomorphometry and pharmacokinetic data was used to assess the effects of treatment versus control on μCT scans of the same biopsy samples over time, matching the changes in bone volume fraction observed in biopsies from denosumab and placebo groups through year 10 of the FREEDOM Extension trial. In the simulation, treatment decreased osteoclast number, which led to a modest increase in trabecular thickness and osteocyte stress shielding. Long‐term bone turnover suppression led to increased RANKL production, followed by a small increase in osteoclast number at the end of the 6‐month–dosing interval, especially at the end of the Extension study. Lack of treatment led to a significant loss of bone mass and structure. The study's results show how in silico models can generate predictions of denosumab cellular action over a 10‐year period, matching static and dynamic morphometric measures assessed in clinical biopsies. The use of in silico models with clinical trial data can be a method to gain further insight into fundamental bone biology and how treatments can perturb this. With rigorous validation, such models could be used for informing the design of clinical trials, such that the number of participants could be reduced to a minimum to show efficacy. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.
Randomized, Double-Blind Study of Denosumab Versus Zoledronic Acid in the Treatment of Bone Metastases in Patients With Advanced Cancer (Excluding Breast and Prostate Cancer) or Multiple Myeloma
Purpose This study compared denosumab, a fully human monoclonal anti-receptor activator of nuclear factor kappa-B ligand antibody, with zoledronic acid (ZA) for delaying or preventing skeletal-related events (SRE) in patients with advanced cancer and bone metastases (excluding breast and prostate) or myeloma. Patients and Methods Eligible patients were randomly assigned in a double-blind, double-dummy design to receive monthly subcutaneous denosumab 120 mg (n = 886) or intravenous ZA 4 mg (dose adjusted for renal impairment; n = 890). Daily supplemental calcium and vitamin D were strongly recommended. The primary end point was time to first on-study SRE (pathologic fracture, radiation or surgery to bone, or spinal cord compression). Results Denosumab was noninferior to ZA in delaying time to first on-study SRE (hazard ratio, 0.84; 95% CI, 0.71 to 0.98; P = .0007). Although directionally favorable, denosumab was not statistically superior to ZA in delaying time to first on-study SRE (P = .03 unadjusted; P = .06 adjusted for multiplicity) or time to first-and-subsequent (multiple) SRE (rate ratio, 0.90; 95% CI, 0.77 to 1.04; P = .14). Overall survival and disease progression were similar between groups. Hypocalcemia occurred more frequently with denosumab. Osteonecrosis of the jaw occurred at similarly low rates in both groups. Acute-phase reactions after the first dose occurred more frequently with ZA, as did renal adverse events and elevations in serum creatinine based on National Cancer Institute Common Toxicity Criteria for Adverse Events grading. Conclusion Denosumab was noninferior (trending to superiority) to ZA in preventing or delaying first on-study SRE in patients with advanced cancer metastatic to bone or myeloma. Denosumab represents a potential novel treatment option with the convenience of subcutaneous administration and no requirement for renal monitoring or dose adjustment.
Denosumab versus zoledronic acid in bone disease treatment of newly diagnosed multiple myeloma: an international, double-blind, double-dummy, randomised, controlled, phase 3 study
In patients with newly diagnosed multiple myeloma, denosumab was non-inferior to zoledronic acid for time to skeletal-related events. The results from this study suggest denosumab could be an additional option for the standard of care for patients with multiple myeloma with bone disease.
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.
Systems Modeling of Bortezomib and Dexamethasone Combinatorial Effects on Bone Homeostasis in Multiple Myeloma Patients
Osteolytic bone disease is one of the most debilitating manifestations of multiple myeloma (MM). Bortezomib is a proteasome inhibitor that shows both anticancer and bone anabolic properties and is being evaluated for its positive effects in MM patients with skeletal complications. Dexamethasone is a potent corticosteroid that is often given in combination with bortezomib for its antineoplastic effects; however, bone loss and osteoporosis are major adverse effects of long-term steroid-based therapies.
Physiological and pathophysiological bone turnover — role of the immune system
Receptor activator of NF-κB ligand (RANKL) and osteoprotegerin (OPG) are key downstream effectors of bone resorption; tumour necrosis factor (TNF) might synergize with RANKL to superinduce osteoclastic bone resorptionB cells, regulated by T cells, are a key source of basal OPG whereas activated T cells and B cells are key sources of TNF and RANKL in inflammatory conditionsShort-term antiresorptive therapies might safely prevent bone loss associated with combination antiretroviral therapyAnergic and parathyroid-hormone-treated T cells secrete the protein Wnt-10b, which promotes bone formationNovel therapeutic strategies targeting the immune system might promote bone formation and decrease bone resorption to manage osteoporotic bone loss and prevent fracture
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.
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.
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.
Osteoimmunology: The Conceptual Framework Unifying the Immune and Skeletal Systems
The immune and skeletal systems share a variety of molecules, including cytokines, chemokines, hormones, receptors, and transcription factors. Bone cells interact with immune cells under physiological and pathological conditions. Osteoimmunology was created as a new interdisciplinary field in large part to highlight the shared molecules and reciprocal interactions between the two systems in both heath and disease. Receptor activator of NF-κB ligand (RANKL) plays an essential role not only in the development of immune organs and bones, but also in autoimmune diseases affecting bone, thus effectively comprising the molecule that links the two systems. Here we review the function, gene regulation, and signal transduction of osteoimmune molecules, including RANKL, in the context of osteoclastogenesis as well as multiple other regulatory functions. Osteoimmunology has become indispensable for understanding the pathogenesis of a number of diseases such as rheumatoid arthritis (RA). We review the various osteoimmune pathologies, including the bone destruction in RA, in which pathogenic helper T cell subsets [such as IL-17-expressing helper T (Th17) cells] induce bone erosion through aberrant RANKL expression. We also focus on cellular interactions and the identification of the communication factors in the bone marrow, discussing the contribution of bone cells to the maintenance and regulation of hematopoietic stem and progenitors cells. Thus the time has come for a basic reappraisal of the framework for understanding both the immune and bone systems. The concept of a unified osteoimmune system will be absolutely indispensable for basic and translational approaches to diseases related to bone and/or the immune system.
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.
Osteoimmunology: shared mechanisms and crosstalk between the immune and bone systems
Osteoimmunology is an interdisciplinary field, that covers the shared mechanisms and interactions between bone cells and immune cellsReceptor activator of nuclear factor-κB (NF-κB) ligand (RANKL) is an osteoclast-differentiation factor that links the activated immune system and bone loss. In addition, abnormal bone homeostasis has been observed in various mice deficient in immunomodulatory molecules.Osteoclast differentiation is dependent on the transcription factor nuclear factor of activated T cells, cytoplasmic 1 (NFATc1), which is induced and activated by RANKL and its co-stimulatory (immunoglobulin-like) receptors.Interleukin-17 (IL-17)-producing T helper cells (TH17 cells) are the key T-cell subset that links T-cell activation and bone destruction in autoimmune arthritis.Bone cells are involved in the maintenance and mobilization of haematopoietic stem cells.Osteoimmunology is becoming increasingly important for understanding the pathogenesis of, and developing new therapeutic strategies for, diseases that affect both systems.
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.
