Combined Effects of Exercise and Denosumab Treatment on Local Failure in Post-menopausal Osteoporosis-Insights from Bone Remodelling Simulations Accounting for Mineralisation and Damage

Wnt/β-catenin signaling pathway: proteins' roles in osteoporosis and cancer diseases and the regulatory effects of natural compounds on osteoporosis

Osteoblasts are mainly derived from mesenchymal stem cells in the bone marrow. These stem cells can differentiate into osteoblasts, which have the functions of secreting bone matrix, promoting bone formation, and participating in bone remodeling. The abnormality of osteoblasts can cause a variety of bone-related diseases, including osteoporosis, delayed fracture healing, and skeletal deformities. In recent years, with the side effects caused by the application of PTH drugs, biphosphonate drugs, and calmodulin drugs, people have carried out more in-depth research on the mechanism of osteoblast differentiation, and are actively looking for natural compounds for the treatment of osteoporosis. The Wnt/β-catenin signaling pathway is considered to be one of the important pathways of osteoblast differentiation, and has become an important target for the treatment of osteoporosis. The Wnt/β-catenin signaling pathway, whether its activation is enhanced or its expression is weakened, will cause a variety of diseases including tumors. This review will summarize the effect of Wnt/β-catenin signaling pathway on osteoblast differentiation and the correlation between the related proteins in the pathway and human diseases. At the same time, the latest research progress of natural compounds targeting Wnt/β-catenin signaling pathway against osteoporosis is summarized.

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 approach to elucidate the pathways leading to primary osteoporosis: age-related vs. postmenopausal

Numerical models of bone remodelling have traditionally been used to perform in silico tests of bone loss in postmenopausal women and also to simulate the response to different drug treatments. These models simulate the menopausal oestrogen decline by altering certain signalling pathways. However, they do not consider the simultaneous effect that ageing can have on cell function and bone remodelling, and thus on bone loss. Considering ageing and oestrogen decline together is important for designing osteoporosis treatments that can selectively counteract one or the other disease mechanism. A previously developed bone cell population model was adapted to consider the effect of ageing through: (1) the decrease of TGF- β contained in the bone matrix and (2) an increased production of sclerostin by non-skeletal cells. Oestrogen deficiency is simulated in three different ways: (a) an increase in RANKL expression, (b) a decrease in OPG production, and (c) an increase in the responsiveness of osteoclasts to RANKL. The effect of ageing was validated using the cross-sectional study of (Riggs et al. in J Bone Miner Res 19: 1945-1954, 2004) on BMD of trabecular bone of the vertebral body of men. The joint effect of ageing and oestrogen deficiency was validated using these same clinical results but in women. In ageing, the effect of the increasing production of sclerostin is more important than the decrease of TGF- β , while the three mechanisms used to simulate the effect of oestrogen deficiency produce almost identical responses. The results show that an early menopause leads to a lower average density in the fifth decade, but after the sixth decade the average density is independent of the age at menopause. Treatment of osteoporosis with denosumab was also simulated to conclude that the drug is not very effective if started before 10 years after menopause or before age 60.

Bone turnover and mineralisation kinetics control trabecular BMDD and apparent bone density: insights from a discrete statistical bone remodelling model

The mechanical quality of trabecular bone is influenced by its mineral content and spatial distribution, which is controlled by bone remodelling and mineralisation. Mineralisation kinetics occur in two phases: a fast primary mineralisation and a secondary mineralisation that can last from several months to years. Variations in bone turnover and mineralisation kinetics can be observed in the bone mineral density distribution (BMDD). Here, we propose a statistical spatio-temporal bone remodelling model to study the effects of bone turnover (associated with the activation frequency Ac . f ) and mineralisation kinetics (associated with secondary mineralisation T sec ) on BMDD. In this model, individual basic multicellular units (BMUs) are activated discretely on trabecular surfaces that undergo typical bone remodelling periods. Our results highlight that trabecular BMDD is strongly regulated by Ac . f and T sec in a coupled way. Ca wt% increases with lower Ac . f and short T sec . For example, aAc . f = 4 BMU/year/mm3 and T sec = 8 years result in a mean Ca wt% of 25, which is in accordance with Ca wt% values reported in quantitative backscattered electron imaging (qBEI) experiments. However, for lower Ac . f and shorter T sec (from 0.5 to 4 years) one obtains a high Ca wt% and a very narrow skew BMDD to the right. This close link between Ac . f and T sec highlights the importance of considering both characteristics to draw meaningful conclusion about bone quality. Overall, this model represents a new approach to modelling healthy and diseased bone and can aid in developing deeper insights into disease states like osteoporosis.

Computational Analysis of the Influence of Menopause and Ageing on Bone Mineral Density, Exploring the Impact of Bone Turnover and Focal Bone Balance-A Study on Overload and Underload Scenarios

This study aims to investigate the impact of hormonal imbalances during menopause, compounded by the natural ageing process, on bone health. Specifically, it examines the effects of increased bone turnover and focal bone balance on bone mass. A three-dimensional computational bone remodeling model was employed to simulate the response of the femur to habitual loads over a 19-year period, spanning premenopause, menopause, and postmenopause. The model was calibrated using experimental bone mineral density data from the literature to ensure accurate simulations. The study reveals that individual alterations in bone turnover or focal bone balance do not fully account for the observed experimental outcomes. Instead, simultaneous changes in both factors provide a more comprehensive explanation, leading to increased porosity while maintaining the material-to-apparent density ratio. Additionally, different load scenarios were tested, demonstrating that reaching the clinical osteoporosis threshold is independent of the timing of load changes. However, underload scenarios resulted in the threshold being reached approximately 6 years earlier than overload scenarios. These findings hold significant implications for strategies aimed at delaying the onset of osteoporosis and minimizing fracture risks through targeted mechanical stimulation during the early stages of menopause.

Effect of Aerobic/Strength Training on RANKL Gene DNA Methylation Levels

BACKGROUND

The osteoclastogenesis RANKL gene plays a key role in bone remodeling. The hypomethylation of its promoter region may cause osteoporosis. The present study aimed to elucidate the influence of physical activity on DNA methylation changes of RANKL promoter cytosine-phosphate-guanine (CpG)-rich region in active and sedentary adults and to assess the effect of aerobic and strength training on RANKL DNA methylation changes among Tunisian-North African adults.

METHODS

A total of 104 participants including 52 adults (58% males and 42% females) and 52 adults (31% males and 69% females) were recruited for the observational and interventional part of the study, respectively. The intervention consisted of 12 weeks of aerobic training (30 min/session) followed by 10 minutes of strengthening exercises. All participants completed the International Physical Activity Questionnaire and provided blood samples for quantitative methylation-specific polymerase chain reaction (PCR) analysis.

RESULTS

The study revealed a significant difference (P = 6 × 10-10) in the methylation level of the RANKL promoter region between active and sedentary adults, with a 6.68-fold increase observed in the active group. After the intervention, both the trained (P = 41 × 10-5) and untrained (P = .002) groups displayed high methylation levels in the RANKL promoter region. In addition, the trained group exhibited significant improvements in heart rate (P = 2.2 × 10-16), blood pressure (P = 39 × 10-3), maximal oxygen uptake (P = 1.5 × 10-7), and fat mass (P = 7 × 10-4).

CONCLUSION

Exploring epigenetic modifications in the RANKL promoter region may contribute to a more comprehensive understanding of the complexity of osteoporosis. This suggests that aerobic/strength training could potentially improve the bone system, reducing its vulnerability to osteoporosis by increasing RANKL DNA methylation.

Mechanisms underlying the long-term and withdrawal effects of denosumab therapy on bone

Denosumab, a human monoclonal antibody against receptor activator of nuclear factor-κB ligand (RANKL), is a potent inhibitor of osteoclast differentiation and activity. As the first biologic drug used to treat osteoporosis, denosumab has shown potent anti-resorptive properties and anti-fracture efficacy. The effects of this drug are also unique compared with the effects of bisphosphonates: namely, long-term treatment with this drug results in a continuous gain of bone mineral density, whereas withdrawal of the drug results in a transient overshoot in bone turnover and rapid bone loss. Although the mechanisms for these specific effects remain incompletely understood, emerging experimental and clinical data have started to highlight potential biological and pharmacological mechanisms by which denosumab might affect osteoclasts, as well as osteoblasts, and cause both sustained bone gain and bone loss upon treatment cessation. This Perspective discusses those potential mechanisms and the future studies and clinical implications that might ensue from these findings.

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.

Mechanistic PK-PD model of alendronate treatment of postmenopausal osteoporosis predicts bone site-specific response

Alendronate is the most widely used drug for postmenopausal osteoporosis (PMO). It inhibits bone resorption, affecting osteoclasts. Pharmacokinetics (PK) and pharmacodynamics (PD) of alendronate have been widely studied, but few mathematical models exist to simulate its effect. In this work, we have developed a PK model for alendronate, valid for short- and long-term treatments, and a mechanistic PK-PD model for the treatment of PMO to predict bone density gain (BDG) at the hip and lumbar spine. According to our results, at least three compartments are required in the PK model to predict the effect of alendronate in both the short and long terms. Clinical data of a 2-year treatment of alendronate, reproduced by our PK-PD model, demonstrate that bone response is site specific (hip: 7% BDG, lumbar spine: 4% BDG). We identified that this BDG is mainly due to an increase in tissue mineralization and a decrease in porosity. The difference in BDG between sites is linked to the different loading and dependence of the released alendronate on the bone-specific surface and porosity. Osteoclast population diminishes quickly within the first month of alendronate treatment. Osteoblast population lags behind but also falls due to coupling of resorption and formation. Two dosing regimens were studied (70 mg weekly and 10 mg daily), and both showed very similar BDG evolution, indicating that alendronate accumulates quickly in bone and saturates. The proposed PK-PD model could provide a valuable tool to analyze the effect of alendronate and to design patient-specific treatments, including drug combinations.

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.

Assessment of Strategies for Safe Drug Discontinuation and Transition of Denosumab Treatment in PMO—Insights From a Mechanistic PK/PD Model of Bone Turnover

Denosumab (Dmab) treatment against postmenopausal osteoporosis (PMO) has proven very efficient in increasing bone mineral density (BMD) and reducing the risk of bone fractures. However, concerns have been recently raised regarding safety when drug treatment is discontinued. Mechanistic pharmacokinetic-pharmacodynamic (PK-PD) models are the most sophisticated tools to develop patient specific drug treatments of PMO to restore bone mass. However, only a few PK-PD models have addressed the effect of Dmab drug holidays on changes in BMD. We showed that using a standard bone cell population model (BCPM) of bone remodelling it is not possible to account for the spike in osteoclast numbers observed after Dmab discontinuation. We show that inclusion of a variable osteoclast precursor pool in BCPMs is essential to predict the experimentally observed rapid rise in osteoclast numbers and the associated increases in bone resorption. This new model also showed that Dmab withdrawal leads to a rapid increase of damage in the bone matrix, which in turn decreases the local safety factor for fatigue failure. Our simulation results show that changes in BMD strongly depend on Dmab concentration in the central compartment. Consequently, bone weight (BW) might play an important factor in calculating effective Dmab doses. The currently clinically prescribed constant Dmab dose of 60 mg injected every 6 months is less effective in increasing BMD for patients with high BW (2.5% for 80 kg in contrast to 8% for 60 kg after 6 years of treatment). However, bone loss observed 24 months after Dmab withdrawal is less pronounced in patients with high BW (3.5% for 80kg and 8.5% for 60 kg). Finally, we studied how to safely discontinue Dmab treatment by exploring several transitional and combined drug treatment strategies. Our simulation results indicate that using transitional reduced Dmab doses are not effective in reducing rapid bone loss. However, we identify that use of a bisphosphonate (BP) is highly effective in avoiding rapid bone loss and increase in bone tissue damage compared to abrupt withdrawal of Dmab. Furthermore, the final values of BMD and damage were not sensitive to the time of administration of the BP.

Advances in pathogenesis and therapeutic strategies for osteoporosis

Osteoporosis, is the most common bone disorder worldwide characterized by low bone mineral density, leaving affected bones vulnerable to fracture. Bone homeostasis depends on the precise balance between bone resorption by osteoclasts and bone matrix formation by mesenchymal lineage osteoblasts, and involves a series of complex and highly regulated steps. Bone homeostasis will be disrupted when the speed of bone resorption is faster than bone formation. Based on various regulatory mechanisms of bone homeostasis, a series of drugs targeting osteoporosis have emerged in clinical practice, including bisphosphonates, selective estrogen receptor modulators, calcitonin, molecular-targeted drugs and so on. However, many drugs have major adverse effects or are unsuitable for long-term use. Therefore, it is very urgent to find more effective therapeutic drugs based on the new pathogenesis of osteoporosis. In this review, we summarize novel mechanisms involved in the pathological process of osteoporosis, including the roles of gut microbiome, autophagy, iron balance and cellular senescence. Based on the above pathological mechanism, we found promising drugs for osteoporosis treatment, such as: probiotics, alpha-ketoglutarate, senolytics and hydrogen sulfide. This new finding may provide an important basis for elucidating the complex pathological mechanisms of osteoporosis and provide promising drugs for clinical osteoporosis treatment.