Ten-Year Simulation of the Effects of Denosumab on Bone Remodeling in Human Biopsies

Advances in mechanobiological pharmacokinetic-pharmacodynamic models of osteoporosis treatment - Pathways to optimise and exploit existing therapies

Osteoporosis (OP) is a chronic progressive bone disease which is characterised by reduction of bone matrix volume and changes in the bone matrix properties which can ultimately lead to bone fracture. The two major forms of OP are related to aging and/or menopause. With the worldwide increase of the elderly population, particularly age-related OP poses a serious health issue which puts large pressure on health care systems. A major challenge for development of new drug treatments for OP and comparison of drug efficacy with existing treatments is due to current regulatory requirements which demand testing of drugs based on bone mineral density (BMD) in phase 2 trials and fracture risk in phase 3 trials. This requires large clinical trials to be conducted and to be run for long time periods, which is very costly. This, together with the fact that there are already many drugs available for treatment of OP, makes the development of new drugs inhibitive. Furthermore, an increased trend of the use of different sequential drug therapies has been observed in OP management, such as sequential anabolic-anticatabolic drug treatment or switching from one anticatabolic drug to another. Running clinical trials for concurrent and sequential therapies is neither feasible nor practical due to large number of combinatorial possibilities. In silico mechanobiological pharmacokinetic-pharmacodynamic (PK-PD) models of OP treatments allow predictions beyond BMD, i.e. bone microdamage and degree of mineralisation can also be monitored. This will help to inform clinical drug usage and development by identifying the most promising scenarios to be tested clinically (confirmatory trials rather than exploratory only trials), optimise trial design and identify subgroups of the population that show benefit-risk profiles (both good and bad) that are different from the average patient. In this review, we provide examples of the predictive capabilities of mechanobiological PK-PD models. These include simulation results of PMO treatment with denosumab, implications of denosumab drug holidays and coupling of bone remodelling models with calcium and phosphate systems models that allows to investigate the effects of co-morbidities such as hyperparathyroidism and chronic kidney disease together with calcium and vitamin D status on drug efficacy.

An in silico micro-multiphysics agent-based approach for simulating bone regeneration in a mouse femur defect model

Bone defects represent a challenging clinical problem as they can lead to non-union. In silico models are well suited to study bone regeneration under varying conditions by linking both cellular and systems scales. This paper presents an in silico micro-multiphysics agent-based (micro-MPA) model for bone regeneration following an osteotomy. The model includes vasculature, bone, and immune cells, as well as their interaction with the local environment. The model was calibrated by time-lapsed micro-computed tomography data of femoral osteotomies in C57Bl/6J mice, and the differences between predicted bone volume fractions and the longitudinal in vivo measurements were quantitatively evaluated using root mean square error (RMSE). The model performed well in simulating bone regeneration across the osteotomy gap, with no difference (5.5% RMSE, p = 0.68) between the in silico and in vivo groups for the 5-week healing period - from the inflammatory phase to the remodelling phase - in the volume spanning the osteotomy gap. Overall, the proposed micro-MPA model was able to simulate the influence of the local mechanical environment on bone regeneration, and both this environment and cytokine concentrations were found to be key factors in promoting bone regeneration. Further, the validated model matched clinical observations that larger gap sizes correlate with worse healing outcomes and ultimately simulated non-union. This model could help design and guide future experimental studies in bone repair, by identifying which are the most critical in vivo experiments to perform.

Megakaryocyte Secreted Factors Regulate Bone Marrow Niche Cells During Skeletal Homeostasis, Aging, and Disease

The bone marrow microenvironment contains a diverse array of cell types under extensive regulatory control and provides for a novel and complex mechanism for bone regulation. Megakaryocytes (MKs) are one such cell type that potentially acts as a master regulator of the bone marrow microenvironment due to its effects on hematopoiesis, osteoblastogenesis, and osteoclastogenesis. While several of these processes are induced/inhibited through MK secreted factors, others are primarily regulated by direct cell-cell contact. Notably, the regulatory effects that MKs exert on these different cell populations has been found to change with aging and disease states. Overall, MKs are a critical component of the bone marrow that should be considered when examining regulation of the skeletal microenvironment. An increased understanding of the role of MKs in these physiological processes may provide insight into novel therapies that can be used to target specific pathways important in hematopoietic and skeletal disorders.

Changes in RANKL and TRAcP 5b after discontinuation of denosumab suggest RANKL mediated formation of osteoclasts results in the increased bone resorption

In patients discontinuing long-term denosumab, RANKL levels are high 6 months after the last denosumab injection. Nine and 12 months after the last denosumab injection RANKL levels are lower, but TRAcP 5b levels are higher, suggesting that accumulated RANKL increases the number of active osteoclasts.

PURPOSE

The rapid increase in bone turnover occurring when discontinuing long-term treatment with denosumab (DMAB), an antibody that neutralizes receptor activator of nuclear factor kappa-B ligand (RANKL) is not fully understood. We aimed to investigate the mechanisms underlying the rebound activation of bone resorption by measuring tartrate-resistant acid phosphatase 5b (TRAcP 5b), RANKL, osteoprotegerin (OPG), C-terminal collagen crosslinks (CTX), and procollagen type I N-propeptide (P1NP) in patients discontinuing long-term DMAB.

METHODS

Sixty-one patients with BMD T-score > - 2.5 at the spine and hip discontinuing long-term DMAB were randomized to treatment with zoledronate (ZOL) 6 months (6 M group, n = 20), 9 months (9 M group, n = 20) or 12 months after the last DMAB injection or when bone turnover was high (12 M group, n = 21). Bone turnover markers were measured immediately before initiation of ZOL treatment.

RESULTS

We found higher CTX and PINP in the 9 M and 12 M groups compared to the 6 M group (p < 0.001). In the 6 M group, TRAcP 5b was lower and RANKL higher than in the other two groups (p < 0.001). TRAcP 5b correlated negatively with RANKL (R = - 0.54), and time since the last DMAB injection correlated positively with CTX (R = 0.56), PINP (R = 0.72), TRAcP 5b (R = 0.51) and negatively with RANKL (R = - 0.70) (p < 0.001 for all). We found no difference in OPG between groups.

CONCLUSION

Following discontinuation of long-term DMAB, we find high levels of RANKL, which most likely result in an increase in the number of active osteoclasts (illustrated by TRAcP5b) causing an increased bone turnover.

Spatio-temporal simulations of bone remodelling using a bone cell population model based on cell availability

Here we developed a spatio-temporal bone remodeling model to simulate the action of Basic Multicelluar Units (BMUs). This model is based on two major extensions of a temporal-only bone cell population model (BCPM). First, the differentiation into mature resorbing osteoclasts and mature forming osteoblasts from their respective precursor cells was modelled as an intermittent process based on precursor cells availability. Second, the interaction between neighbouring BMUs was considered based on a "metabolic cost" argument which warrants that no new BMU will be activated in the neighbourhood of an existing BMU. With the proposed model we have simulated the phases of the remodelling process obtaining average periods similar to those found in the literature: resorption ( ∼ 22 days)-reversal (∼8 days)-formation (∼65 days)-quiescence (560-600 days) and an average BMU activation frequency of ∼1.6 BMUs/year/mm³. We further show here that the resorption and formation phases of the BMU become coordinated only by the presence of TGF-β (transforming growth factor β), i.e., a major coupling factor stored in the bone matrix. TGF-β is released through resorption so upregulating osteoclast apoptosis and accumulation of osteoblast precursors, i.e., facilitating the transition from the resorption to the formation phase at a given remodelling site. Finally, we demonstrate that this model can explain targeted bone remodelling as the BMUs are steered towards damaged bone areas in order to commence bone matrix repair.

Trabecular bone remodeling in the aging mouse: A micro-multiphysics agent-based in silico model using single-cell mechanomics

Bone remodeling is regulated by the interaction between different cells and tissues across many spatial and temporal scales. Notably, in silico models are regarded as powerful tools to further understand the signaling pathways that regulate this intricate spatial cellular interplay. To this end, we have established a 3D multiscale micro-multiphysics agent-based (micro-MPA) in silico model of trabecular bone remodeling using longitudinal in vivo data from the sixth caudal vertebra (CV6) of PolgA^{(D257A/D257A)} mice, a mouse model of premature aging. Our in silico model includes a variety of cells as single agents and receptor-ligand kinetics, mechanomics, diffusion and decay of cytokines which regulate the cells' behavior. We highlighted its capabilities by simulating trabecular bone remodeling in the CV6 of five mice over 4 weeks and we evaluated the static and dynamic morphometry of the trabecular bone microarchitecture. Based on the progression of the average trabecular bone volume fraction (BV/TV), we identified a configuration of the model parameters to simulate homeostatic trabecular bone remodeling, here named basal. Crucially, we also produced anabolic, anti-anabolic, catabolic and anti-catabolic responses with an increase or decrease by one standard deviation in the levels of osteoprotegerin (OPG), receptor activator of nuclear factor kB ligand (RANKL), and sclerostin (Scl) produced by the osteocytes. Our results showed that changes in the levels of OPG and RANKL were positively and negatively correlated with the BV/TV values after 4 weeks in comparison to basal levels, respectively. Conversely, changes in Scl levels produced small fluctuations in BV/TV in comparison to the basal state. From these results, Scl was deemed to be the main driver of equilibrium while RANKL and OPG were shown to be involved in changes in bone volume fraction with potential relevance for age-related bone features. Ultimately, this micro-MPA model provides valuable insights into how cells respond to their local mechanical environment and can help to identify critical pathways affected by degenerative conditions and ageing.

Denosumab Can Prevent Collapse in Patients with Early-Stage Steroid-Induced Osteonecrosis of the Femoral Head by Inhibiting Osteoclasts and Autophagy

OBJECTIVE

The osteoclastic bone resorption inhibitors might have positive effect in preventing femoral head collapse in patients with osteonecrosis of the femoral head (ONFH). However, as a novel osteoclastic inhibitor, whether denosumab can prevent collapse in steroid-induced ONFH remains unknown. This study aims to evaluate the treatment effect of denosumab and the potential protective mechanism.

METHODS

This was a retrospective study. A total of 161 patients with steroid-induced ONFH who underwent denosumab treatment were reviewed, and 209 untreated patients were selected as controls. Their clinical characteristics and radiological exam results were obtained. Patients were treated with 60 mg denosumab every 6 months for 2 years. The primary outcome was the incidence of femoral head collapse at 2 years after the initial diagnosis of ONFH. Secondary outcomes included the Harris hip score, progression of osteosclerosis, increase in necrotic area, bone marrow oedema relief, and bone mineral density increase in the femoral head. The Mann-Whitney U test and chi-square tests were performed to identify the differences between the continuous and categorical variables, respectively. A multivariate logistic regression model was built to identify the factors associated with the treatment effect of denosumab.

RESULTS

The incidence of femoral head collapse was 42.24% (68/161) in the denosumab group and 54.07% (113/209) in the control group (χ²  = 5.094, p = 0.024; relative risk = 0.787, 95% CI = 0.627-0.973). The excellent-good rates of the Harris hip score were 63.98% (103/161) in the denosumab group and 44.98% (94/209) in the control group (χ²  = 13.186, p < 0.001). The incidence of osteosclerosis progression in the denosumab group was 55.28% (89/161), which was significantly higher than that in the control group (43.54%, 91/209, χ²  = 5.016, p = 0.025). Meanwhile, a significant increase in bone mineral density was identified in 29.19% (47/161) and 7.18% (15/209) of patients in the denosumab and control groups, respectively (χ²  = 31.600, p < 0.001). The osteoclastic cytoplasm expression of LC3-II was more positive in the control group than in the denosumab group (immunohistochemistry scoring: 3.58 ± 2.27 vs 6.33 ± 2.64, Z = -2.684, p = 0.007). A total of three independent factors were considered to be associated with the positive treatment effect of denosumab, the time of first denosumab administration (OR = 2.010, 95% CI = 1.272-3.177), osteosclerosis (OR = 1.583, 95% CI = 1.024-2.445), and the necrotic area before denosumab administration (medium necrotic area: OR = 2.084, 95% CI = 1.245-3.487; large necrotic area: OR = 2.211, 95% CI = 1.255-3.893).

CONCLUSIONS

The current study demonstrated that denosumab had a positive effect on preventing femoral head collapse in patients with steroid ONFH. This effect might be closely associated with the inhibition of osteoclasts and their autophagy.

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.

The effect of osteoporosis treatment on bone mass

Over the last two decades there have been significant developments in the pharmacotherapy of osteoporosis. The therapeutic arsenal has expanded with monoclonal antibodies which have been developed based on discoveries of the molecular mechanisms underlying bone resorption and bone formation. Denosumab, the antibody binding RANKL, inhibits bone resorption, and romosozumab, the antibody binding sclerostin, inhibits bone resorption and stimulates bone formation as well. Both antibodies have shown potent anti-fracture efficacy in randomized clinical trials and this review will discuss the preclinical and clinical studies focusing on the effects on bone mass. After discontinuation of these antibodies, bone mineral density quickly returns to baseline and in the case of denosumab, discontinuation can not only induce rebound bone loss, but also the occurrence of vertebral fractures. Therefore, sequential antiresorptive therapy to maintain bone mass gains and anti-fracture efficacy is of utmost importance and will also be discussed in this review.