Influence of different mechanical stimuli in a multi-scale mechanobiological isotropic model for bone remodelling

A 3D in Silico Multi-Tissue Evolution Model Highlights the Relevance of Local Strain Accumulation in Bone Fracture Remodeling

Since 5–10% of all bone fractures result in non-healing situations, a thorough understanding of the various bone fracture healing phases is necessary to propose adequate therapeutic strategies. In silico models have greatly contributed to the understanding of the influence of mechanics on tissue formation and resorption during the soft and hard callus phases. However, the late-stage remodeling phase has not been investigated from a mechanobiological viewpoint so far. Here, we propose an in silico multi-tissue evolution model based on mechanical strain accumulation to investigate the mechanobiological regulation of bone remodeling during the late phase of healing. Computer model predictions are compared to histological data of two different pre-clinical studies of bone healing. The model predicted the bone marrow cavity re-opening and the resorption of the external callus. Our results suggest that the local strain accumulation can explain the fracture remodeling process and that this mechanobiological response is conserved among different mammal species. Our study paves the way for further understanding of non-healing situations that could help adapting therapeutic strategies to foster bone healing.

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.

NUMERICAL SIMULATION OF APPARENT DENSITY EVOLUTION OF TRABECULAR BONE UNDER FATIGUE LOADING: EFFECT OF BONE INITIAL PROPERTIES

Bone remodeling is a physiological phenomenon coupling resorption and formation processes that are mainly mediated by osteoclasts and osteoblasts, in response to mechanical stimuli transduced by osteocytes to biochemical signals activating the bone multicellular unit. Under normal loading conditions, bone resorption and formation are balanced by a homeostasis process. When bone is subjected to overstress, microdamaging occurs, which induces a modification of the structural integrity and microarchitecture. This has drawn significant attention to the mechanical properties of bone. In this context, the current study has been carried out with the aim of numerically investigating the impact of the mechanical properties on the remodeling process of the trabecular bone under cyclic loading, highlighting the effects of different values of the mineral density and the Young’s modulus. This was performed using a mechanobiological model, coupling mechanical and biological approaches, allowing to numerically simulate the effect of the selected parameters for a 20-year-period of cyclic loading for 2D and 3D models of a human femur head. The current work is an explorative numerical study, and the obtained results revealed the changes in the overall stiffness of the bone according to the mechanical properties.

The Influence of Geometry of Implants for Direct Skeletal Attachment of Limb Prosthesis on Rehabilitation Program and Stress-Shielding Intensity

The purpose of the research was to evaluate the influence of selected parameters of the implants for bone anchored prostheses on possibility of conducting static load bearing exercises and stress-shielding intensity. A press-fit implant, a threaded implant, and the proposed design were compared using the finite element method. For the analyses two features were examined: diameter (19.0 – 21.0 mm) and length (75.0 – 130.0 mm). To define the possibility of conducting rehabilitation exercises the micromotion of implants while axial loading with a force up to 1000 N was examined to evaluate the changes at implant-bone interface. The stress-shielding intensity was estimated by bone mass loss over 60 months. The results suggest that, in terms of micromotion generated during rehabilitation exercises, the threaded (max. micromotion of 16.00 μm) and the proposed (max. micromotion of 45.43 μm) implants ensure low and appropriate micromotion. In the case of the press-fit solution the load values should be selected with care, as there is a risk of losing primary stabilisation. The allowed forces (that do not stimulate the organism to generate fibrous tissue) were approx. 140 N in the case of the length of 75 mm, increasing up to 560 N, while using the length of 130 mm. Moreover, obtained stress-shielding intensities suggest that the proposed implant should provide appropriate secondary stability, similar to the threaded solution, due to the low bone mass loss during long-term use (improving at the same time more bone remodelling in distal Gruen zones, by providing lower bone mass loss by approx. 13% to 20% in dependency of the length and diameter used). On this basis it can be concluded that the proposed design can be an appropriate alternative to commercially used implants.

A novel mechanobiological model can predict how physiologically relevant dynamic loading causes proteoglycan loss in mechanically injured articular cartilage

Cartilage provides low-friction properties and plays an essential role in diarthrodial joints. A hydrated ground substance composed mainly of proteoglycans (PGs) and a fibrillar collagen network are the main constituents of cartilage. Unfortunately, traumatic joint loading can destroy this complex structure and produce lesions in tissue, leading later to changes in tissue composition and, ultimately, to post-traumatic osteoarthritis (PTOA). Consequently, the fixed charge density (FCD) of PGs may decrease near the lesion. However, the underlying mechanisms leading to these tissue changes are unknown. Here, knee cartilage disks from bovine calves were injuriously compressed, followed by a physiologically relevant dynamic compression for twelve days. FCD content at different follow-up time points was assessed using digital densitometry. A novel cartilage degeneration model was developed by implementing deviatoric and maximum shear strain, as well as fluid velocity controlled algorithms to simulate the FCD loss as a function of time. Predicted loss of FCD was quite uniform around the cartilage lesions when the degeneration algorithm was driven by the fluid velocity, while the deviatoric and shear strain driven mechanisms exhibited slightly discontinuous FCD loss around cracks. Our degeneration algorithm predictions fitted well with the FCD content measured from the experiments. The developed model could subsequently be applied for prediction of FCD depletion around different cartilage lesions and for suggesting optimal rehabilitation protocols.

A comparative analysis of internal bone remodelling concepts in a novel implant for direct skeletal attachment of limb prosthesis evaluation: A finite element analysis

Nowadays, numerous internal bone remodelling concepts are under development, in order to estimate long-term functionality of implants by evaluating the intensity of stress-shielding effect. This effect is also analysed for the implants for direct skeletal attachment, considered as a better exoprosthesis fixation method than prosthetic sockets. Most of bone remodelling approaches are based on basic concepts, differing with certain assumptions, which may affect the accuracy of the results. This article compares commonly used internal bone remodelling concepts and evaluates the functionality of the proposed Limb Prosthesis Osseointegrated Fixation System for direct skeletal attachment of limb prosthesis in comparison with two currently available implants: the Intraosseous Transcutaneous Amputation Prosthesis and the Osseointegrated Prostheses for the Rehabilitation of Amputees. Three concepts were chosen: without and with lazy zone and with the use of quadratic formula which considers bone overloading. Therefore, three finite element models were created with identical femur implanted with each of analysed implants. The implants were loaded with loads that refer to two stages of gait cycle (heel strike and toe-off). The analysed concepts have given similar results, allowing to assume that each of them can be successfully used to estimate internal bone remodelling around analysed implants for direct skeletal attachment of limb prosthesis. The results also present higher functionality of the proposed implant for direct skeletal attachment of limb prosthesis due to a significant reduction in stress-shielding in the analysed areas around implant in comparison with the Intraosseous Transcutaneous Amputation Prosthesis and the Osseointegrated Prostheses for the Rehabilitation of Amputees. It suggests that the proposed design is a better alternative to the currently used solutions.