A coupled mechano-biochemical model for bone adaptation

A review of the current status and future prospects of the bone remodeling process: Biological and mathematical perspectives

This review dives into the complex dynamics of bone remodeling, combining biological insights with mathematical perspectives to better understand this fundamental aspect of skeletal health. Bone, being a crucial part of our body, constantly renews itself, and with the growing number of individuals facing bone-related issues, research in this field is vital. In this review, we categorized and classified most common mathematical models used to simulate the mechanical behavior of bone under different loading and health conditions, shedding light on the evolving landscape of bone biology. While current models have effectively captured the essence of healthy bone remodeling, the ever-expanding knowledge in bone biology suggests an update in mathematical methods. Knowing the role of the skeleton in whole-body physiology, and looking at the recent discoveries about activities of bone cells emphasize the urgency of refining our mathematical descriptions of the bone remodeling process. The underexplored impact of bone diseases like osteoporosis, Paget's disease, or breast cancer on bone remodeling also points to the need for intensified research into diverse disease types and their unique effects on bone health. By reviewing a range of bone remodeling models, we show the necessity for tailor-made mathematical models to decipher their roots and enhance patient treatment strategies. Collaboration among scientists from various domains is pivotal to surmount these challenges, ensuring improved accuracy and applicability of mathematical models. Ultimately, this effort aims to deepen our understanding of bone remodeling processes and their broader implications for diverse health conditions.

Tumor growth for remodeling process: A 2D approach

This paper aims to present a comprehensive framework for coupling tumor-bone remodeling processes in a 2-dimensional geometry. This is achieved by introducing a bio-inspired damage that represents the growing tumor, which subsequently affects the main populations involved in the remodeling process, namely, osteoclasts, osteoblasts, and bone tissue. The model is constructed using a set of differential equations based on the Komarova's and Ayati's models, modified to incorporate the bio-inspired damage that may result in tumor mass formation. Three distinct models were developed. The first two models are based on the Komarova's governing equations, with one demonstrating an osteolytic behavior and the second one an osteoblastic model. The third model is a variation of Ayati's model, where the bio-inspired damage is induced through the paracrine and autocrine parameters, exhibiting an osteolytic behavior. The obtained results are consistent with existing literature, leading us to believe that our in-silico experiments will serve as a cornerstone for paving the way towards targeted interventions and personalized treatment strategies, ultimately improving the quality of life for those affected by these conditions.

Computational model of articular cartilage regeneration induced by scaffold implantation in vivo

Computational models allow to explain phenomena that cannot be observed through an animal model, such as the strain and stress states which can highly influence regeneration of the tissue. For this purpose, we have developed a simulation tool to determine the mechanical conditions provided by the polymeric scaffold. The computational model considered the articular cartilage, the subchondral bone, and the scaffold. All materials were modeled as poroelastic, and the cartilage had linear-elastic oriented collagen fibers. This model was able to explain the remodeling process that subchondral bone goes through, and how the scaffold allowed the conditions for cartilage regeneration. These results suggest that the use of scaffolds might lead the cartilaginous tissue growth in vivo by providing a better mechanical environment. Moreover, the developed computational model demonstrated to be useful as a tool prior experimental in vivo studies, by predicting the possible outcome of newly proposed treatments allowing to discard approaches that might not bring good results.

Computational prediction of the long-term behavior of the femoral density after THR using the Silent Hip stem

Aseptic loosening due to the progressive periprosthetic bone resorption following total hip replacement is a crucial concern, that causes complications and failure of the arthroplasty surgery. The mismatch in stiffness between the hip implant and the surrounding femoral bone is one of the key factors leading to bone density resorption. This paper aimed to investigate the long-term response of the femoral bone after THR using the Silent Hip stem. For this purpose, a validated thermodynamic-based computational model was used to compute the change in bone density before and after THR. This model incorporated essential factors involved in bone remodeling process, such as mechanical loading, and biochemical affinities. The results of the numerical simulations using 3D finite element analysis were analyzed in five zones of interest qualitatively and quantitatively. Bone density predictions showed notable bone resorption in cervical areas, specifically in zone 1 and zone 5 of -18.7% and -14%, respectively. Conversely, bone formation was observed in the greater trochanter area (zone 2) of +25%. Stress shielding seemed to occur at cervical area due to the reduction in the mechanical loading in this region. Based on the quantitative analysis of the bone density distribution throughout the femoral bone, it appears that the Silent Hip stem achieved less bone resorption compared to conventional hip stem designs reported in the literature, which could be used for active patients.

Multiscale modeling of bone tissue mechanobiology

Mechanical environment has a crucial role in our organism at the different levels, ranging from cells to tissues and our own organs. This regulatory role is especially relevant for bones, given their importance as load-transmitting elements that allow the movement of our body as well as the protection of vital organs from load impacts. Therefore bone, as living tissue, is continuously adapting its properties, shape and repairing itself, being the mechanical loads one of the main regulatory stimuli that modulate this adaptive behavior. Here we review some key results of bone mechanobiology from computational models, describing the effect that changes associated to the mechanical environment induce in bone response, implant design and scaffold-driven bone regeneration.

A General Mechano-Pharmaco-Biological Model for Bone Remodeling Including Cortisol Variation

The process of bone remodeling requires a strict coordination of bone resorption and formation in time and space in order to maintain consistent bone quality and quantity. Bone-resorbing osteoclasts and bone-forming osteoblasts are the two major players in the remodeling process. Their coordination is achieved by generating the appropriate number of osteoblasts since osteoblastic-lineage cells govern the bone mass variation and regulate a corresponding number of osteoclasts. Furthermore, diverse hormones, cytokines and growth factors that strongly link osteoblasts to osteoclasts coordinated these two cell populations. The understanding of this complex remodeling process and predicting its evolution is crucial to manage bone strength under physiologic and pathologic conditions. Several mathematical models have been suggested to clarify this remodeling process, from the earliest purely phenomenological to the latest biomechanical and mechanobiological models. In this current article, a general mathematical model is proposed to fill the gaps identified in former bone remodeling models. The proposed model is the result of combining existing bone remodeling models to present an updated model, which also incorporates several important parameters affecting bone remodeling under various physiologic and pathologic conditions. Furthermore, the proposed model can be extended to include additional parameters in the future. These parameters are divided into four groups according to their origin, whether endogenous or exogenous, and the cell population they affect, whether osteoclasts or osteoblasts. The model also enables easy coupling of biological models to pharmacological and/or mechanical models in the future.

On the effect of antiresorptive drugs on the bone remodeling of the mandible after dental implantation: a mathematical model

Bone remodeling identifies the process of permanent bone change with new bone formation and old bone resorption. Understanding this process is essential in many applications, such as optimizing the treatment of diseases like osteoporosis, maintaining bone density in long-term periods of disuse, or assessing the long-term evolution of the bone surrounding prostheses after implantation. A particular case of study is the bone remodeling process after dental implantation. Despite the overall success of this type of implants, the increasing life expectancy in developed countries has boosted the demand for dental implants in patients with osteoporosis. Although several studies demonstrate a high success rate of dental implants in osteoporotic patients, it is also known that the healing time and the failure rate increase, necessitating the adoption of pharmacological measures to improve bone quality in those patients. However, the general efficacy of these antiresorptive drugs for osteoporotic patients is still controversial, requiring more experimental and clinical studies. In this work, we investigate the effect of different doses of several drugs, used nowadays in osteoporotic patients, on the evolution of bone density after dental implantation. With this aim, we use a pharmacokinetic-pharmacodynamic (PK/PD) mathematical model that includes the effect of antiresorptive drugs on the RANK/RANK-L/OPG pathway, as well as the mechano-chemical coupling with external mechanical loads. This mechano-PK/PD model is then used to analyze the evolution of bone in normal and osteoporotic mandibles after dental implantation with different drug dosages. We show that using antiresorptive agents such as bisphosphonates or denosumab increases bone density and the associated mechanical properties, but at the same time, it also increases bone brittleness. We conclude that, despite the many limitations of these very complex models, the one presented here is capable of predicting qualitatively the evolution of some of the main biological and chemical variables associated with the process of bone remodeling in patients receiving drugs for osteoporosis, so it could be used to optimize dental implant design and coating for osteoporotic patients, as well as the drug dosage protocol for patient-specific treatments.

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, Paget's disease, osteopetrosis, and metastatic bone disease, which are not well defined in terms of changes in biochemical and mechanobiological regulatory factors. 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 becoming increasingly complex. The current paper provides a timely and critical analysis of previously developed bone remodeling mathematical models, particularly those addressing bone diseases. Besides, mechanistic pharmacokinetic/pharmacodynamic (PK/PD) models, which englobe bone disease and its treatment's 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 concludes with exhibiting the available bone treatment procedures considered in the PK/PD models.

Simulation of bone remodeling around a femoral prosthesis using a model that accounts for biological and mechanical interactions

The present study focuses on a model for three-dimensional bone remodeling of the human femur that considers cellular dynamics to determine the volume fraction of new bone. The model considers the interaction among responsive osteoblasts, active osteoblasts, and osteoclasts, as well as signaling molecules and parathyroid hormone (PTH). The stimulus of the model has a systemic origin due to the PTH effect, and a local origin due to the action of cytokines, growth factors, and mechanical stimuli near the site of the bone cells. The present work considers that the mechanical stimulus that activates cellular activity is obtained from stresses acting on the bone tissue and the number of daily loading cycles. In addition to simulating the bone modeling process in an intact femur, the numerical model is used to simulate bone adaptation in relation to the stress shielding phenomenon after the implantation of a femoral prosthesis. The results showed that the simulations provide a distribution of bone density that is similar to a radiograph and, in addition, allows the visualization of osteoblast and osteoclast dynamics in bone adaptation response after prosthesis implantation.

A mechano-chemo-biological model for bone remodeling with a new mechano-chemo-transduction approach

Bone remodeling is a fundamental biological process that develops in bone tissue along its whole lifetime. It refers to a continuous bone transformation with new bone formation and old bone resorption that changes the internal microstructure and composition of the tissue. The main objectives of bone remodeling are: repair of the internal microcracks; adaptation of the macroscopic stiffness and strength to the actual changing mechanical demands; and control of the calcium homeostasis. Understanding this process and predicting its evolution is critical to reduce the effects of long-term disuse as happens during periods of reduced mobility. It is also important in the design of bone implants to avoid long-term stress shielding. Many mathematical models have been proposed from the earliest purely phenomenological to the latest that include biological knowledge. However, there still exists a lack of connection between the mechanical driving force and the biochemical and cell processes it triggers. Here, and following previous works that model independently the mechanobiological and biochemical processes in bone remodeling, we present a more complete model, useful for both cortical and trabecular bone, that uses a new mechanotransduction approach based on the effect of strains onto the bonding-unbonding rate of RANK/RANKL/OPG receptor-ligand reactions. We compare the results of this model with previous ones, showing a good agreement in similar conditions. We also apply it to realistic situations such as a femoral bone after implantation of a hip prosthesis, getting similar results to the clinical ones in the final bone density distribution. Finally, we extend this approach to the anisotropic case, getting not only the mean density, but also the directional homogenization of the microstructure. This biochemical approach permits, not only to predict the bone evolution under changes in the mechanical loads, but also, to consider anabolic and catabolic drugs to control bone density, such as those used in osteoporosis.

BONE ORTHOTROPIC REMODELING AS A THERMODYNAMICALLY-DRIVEN EVOLUTION

In this paper, we present and discuss a model of bone remodeling set up in the framework of the theory of generalized continuum mechanics which was first introduced by DiCarlo et al. [Sur le remodelage des tissus osseux anisotropes, Comptes Rendus Mécanique 334(11):651–661, 2006]. Bone is described as an orthotropic body experiencing remodeling as a rotation of its microstructure. Thus, the complete kinematic description of a material point is provided by its position in space and a rotation tensor describing the orientation of its microstructure. Material motion is driven by energetic considerations, namely by the application of the Clausius–Duhem inequality to the microstructured material. Within this framework of orthotropic remodeling, some key features of the remodeling equilibrium configurations are deduced in the case of homogeneous strain or stress loading conditions. First, it is shown that remodeling equilibrium configurations correspond to energy extrema. Second, stability of the remodeling equilibrium configurations is assessed in terms of the local convexity of the strain and complementary energy functionals hence recovering some classical energy theorems. Eventually, it is shown that the remodeling equilibrium configurations are not only highly dependent on the loading conditions, but also on the material properties.

Salvianolic acid B and danshensu induce osteogenic differentiation of rat bone marrow stromal stem cells by upregulating the nitric oxide pathway

The aim of the present study was to investigate the effect of salvianolic acid B (Sal B) and danshensu (DSU) on the osteogenic differentiation of bone marrow mesenchymal stem cells (MSCs) and the mechanisms of the effects. The osteogenic differentiation of MSCs in culture was assessed by measuring alkaline phosphatase (ALP) activity, osteocalcin (OCN) production, nitric oxide (NO) production and the mRNA expression levels of osteoprotegerin (OPG) and its ligand by MSCs. MSCs were successfully induced to differentiate into osteoblasts and adipocytes. Sal B and DSU increased the ALP activity and the production of OCN in the absence of an ossification inducer. The increase in ALP activity was more pronounced when induction was combined with the osteogenic inducer, Sal B, which enhanced the expression of OPG; however, Sal B reduced the expression of receptor activator of nuclear factor-κB ligand (RANKL) by MSCs. Sal B reversed the inhibitory effect of N-nitro L-arginine methylester on the MSCs and increased ALP activity, OCN content and the OPG/RANKL ratio. Based on these results, it was concluded that Sal B increases the osteogenic differentiation of MSCs, most likely by regulating the nitric oxide pathway.

Subject-specific musculoskeletal loading of the tibia: Computational load estimation

The systematic development of subject-specific computer models for the analysis of personalized treatments is currently a reality. In fact, many advances have recently been developed for creating virtual finite element-based models. These models accurately recreate subject-specific geometries and material properties from recent techniques based on quantitative image analysis. However, to determine the subject-specific forces, we need a full gait analysis, typically in combination with an inverse dynamics simulation study. In this work, we aim to determine the subject-specific forces from the computer tomography images used to evaluate bone density. In fact, we propose a methodology that combines these images with bone remodelling simulations and artificial neural networks. To test the capability of this novel technique, we quantify the personalized forces for five subject-specific tibias using our technique and a gait analysis. We compare both results, finding that similar vertical loads are estimated by both methods and that the dominant part of the load can be reliably computed. Therefore, we can conclude that the numerical-based technique proposed in this work has great potential for estimating the main forces that define the mechanical behaviour of subject-specific bone.

On the Use of Bone Remodelling Models to Estimate the Density Distribution of Bones. Uniqueness of the Solution

Bone remodelling models are widely used in a phenomenological manner to estimate numerically the distribution of apparent density in bones from the loads they are daily subjected to. These simulations start from an arbitrary initial distribution, usually homogeneous, and the density changes locally until a bone remodelling equilibrium is achieved. The bone response to mechanical stimulus is traditionally formulated with a mathematical relation that considers the existence of a range of stimulus, called dead or lazy zone, for which no net bone mass change occurs. Implementing a relation like that leads to different solutions depending on the starting density. The non-uniqueness of the solution has been shown in this paper using two different bone remodelling models: one isotropic and another anisotropic. It has also been shown that the problem of non-uniqueness is only mitigated by removing the dead zone, but it is not completely solved unless the bone formation and bone resorption rates are limited to certain maximum values.

Mechano-chemical coupling in Belousov-Zhabotinskii reactions

Mechano-chemical coupling has been recently recognised as an important effect in various systems as chemical reactivity can be controlled through an applied mechanical loading. Namely, Belousov-Zhabotinskii reactions in polymer gels exhibit self-sustained oscillations and have been identified to be reasonably controllable and definable to the extent that they can be harnessed to perform mechanical work at specific locations. In this paper, we use our theoretical work of nonlinear mechano-chemical coupling and investigate the possibility of providing an explanation of phenomena found in experimental research by means of this theory. We show that mechanotransduction occurs as a response to both static and dynamic mechanical stimulation, e.g., volume change and its rate, as observed experimentally and discuss the difference of their effects on oscillations. Plausible values of the quasi-stoichiometric parameter f of Oregonator model are estimated together with its dependence on mechanical stimulation. An increase in static loading, e.g., pressure, is predicted to have stimulatory effect whereas dynamic loading, e.g., rate of volume change, is predicted to be stimulatory only up to a certain threshold. Further, we offer a physically consistent explanation of the observed phenomena why some Belousov-Zhabotinskii gels require an additional mechanical stimulation to show emergence of oscillation or why "revival" of oscillations in Belousov-Zhabotinskii reactions is possible together with indications for further experimental setups.