Mechano-Regulation of Trabecular Bone Adaptation Is Controlled by the Local in vivo Environment and Logarithmically Dependent on Loading Frequency

Effect of different running protocols on bone morphology and microarchitecture of the forelimbs in a male Wistar rat model

BACKGROUND

It is accepted that the metabolic response of bone tissue depends on the intensity of the mechanical loads, but also on the type and frequency of stress applied to it. Physical exercise such as running involves stresses which, under certain conditions, have been shown to have the best osteogenic effects. However, at high intensity, it can be deleterious for bone tissue. Consequently, there is no clear consensus as to which running modality would have the best osteogenic effects.

AIM

Our objective was to compare the effects of three running modalities on morphological and micro-architectural parameters on forelimb bones.

METHODS

Forty male Wistar rats were randomly divided into four groups: high intensity interval training (HIIT), continuous running, combined running ((alternating HIIT and continuous modalities) and sedentary (control). The morphometry, trabecular microarchitecture and cortical porosity of the ulna, radius and humerus were analyzed using micro-tomography.

RESULTS

All three running modalities resulted in bone adaptation, with an increase in the diaphyseal diameter of all three bones. The combined running protocol had positive effects on the trabecular thickness in the distal ulna. The HIIT protocol resulted in an increase in both medio-lateral diameter and cortical bone area over total area (Ct.Ar/Tt.Ar) at the ulnar shaft compared with sedentary condition. Moreover, the HIIT protocol decreased the mean surface area of the medulla (Ma.Ar) according to sedentary condition at the ulnar shaft.

CONCLUSION

This study has shown that HIIT resulted in a decrease in trabecular bone fraction in favor of cortical bone area at the ulna.

Elucidating the mechano-molecular dynamics of TRAP activity using CRISPR/Cas9 mediated fluorescent reporter mice

Osteoclasts are essential for bone remodeling by adapting their resorptive activity in response to their mechanical in vivo environment. However, the molecular mechanisms underlying this process remain unclear. Here, we demonstrated the role of tartrate-resistant acid phosphatase (TRAP, Acp5), a key enzyme secreted by osteoclasts, in bone remodeling and mechanosensitivity. Using CRISPR/Cas9 reporter mice, we demonstrated bone cell reporter (BCRIbsp/Acp5) mice feature fluorescent TRAP-deficient osteoclasts and examined their activity during mechanically driven trabecular bone remodeling. Although BCRIbsp/Acp5 mice exhibited trabecular bone impairments and reduced resorption capacity in vitro, RNA sequencing revealed unchanged levels of key osteoclast-associated genes such as Ctsk, Mmp9, and Calcr. These findings, in conjunction with serum carboxy-terminal collagen crosslinks (CTX) and in vivo mechanical loading outcomes collectively indicated an unaltered bone resorption capacity of osteoclasts in vivo. Furthermore, we demonstrated similar mechanoregulation during trabecular bone remodeling in BCRIbsp/Acp5 and wild-type (WT) mice. Hence, this study provides valuable insights into the dynamics of TRAP activity in the context of bone remodeling and mechanosensation.

Mechanoregulation analysis of bone formation in tissue engineered constructs requires a volumetric method using time-lapsed micro-computed tomography

Bone can adapt its microstructure to mechanical loads through mechanoregulation of the (re)modeling process. This process has been investigated in vivo using time-lapsed micro-computed tomography (micro-CT) and micro-finite element (FE) analysis using surface-based methods, which are highly influenced by surface curvature. Consequently, when trying to investigate mechanoregulation in tissue engineered bone constructs, their concave surfaces make the detection of mechanoregulation impossible when using surface-based methods. In this study, we aimed at developing and applying a volumetric method to non-invasively quantify mechanoregulation of bone formation in tissue engineered bone constructs using micro-CT images and FE analysis. We first investigated hydroxyapatite scaffolds seeded with human mesenchymal stem cells that were incubated over 8 weeks with one mechanically loaded and one control group. Higher mechanoregulation of bone formation was measured in loaded samples with an area under the curve for the receiver operating curve (AUCformation) of 0.633-0.637 compared to non-loaded controls (AUCformation: 0.592-0.604) during culture in osteogenic medium (p < 0.05). Furthermore, we applied the method to an in vivo mouse study investigating the effect of loading frequencies on bone adaptation. The volumetric method detected differences in mechanoregulation of bone formation between loading conditions (p < 0.05). Mechanoregulation in bone formation was more pronounced (AUCformation: 0.609-0.642) compared to the surface-based method (AUCformation: 0.565-0.569, p < 0.05). Our results show that mechanoregulation of formation in bone tissue engineered constructs takes place and its extent can be quantified with a volumetric mechanoregulation method using time-lapsed micro-CT and FE analysis. STATEMENT OF SIGNIFICANCE: Many efforts have been directed towards optimizing bone scaffolds for tissue growth. However, the impact of the scaffolds mechanical environment on bone growth is still poorly understood, requiring accurate assessment of its mechanoregulation. Existing surface-based methods were unable to detect mechanoregulation in tissue engineered constructs, due to predominantly concave surfaces in scaffolds. We present a volumetric approach to enable the precise and non-invasive quantification and analysis of mechanoregulation in bone tissue engineered constructs by leveraging time-lapsed micro-CT imaging, image registration, and finite element analysis. The implications of this research extend to diverse experimental setups, encompassing culture conditions, and material optimization, and investigations into bone diseases, enabling a significant stride towards comprehensive advancements in bone tissue engineering and regenerative medicine.

The loading direction dramatically affects the mechanical properties of the mouse tibia

Introduction: The in vivo tibial loading mouse model has been extensively used to evaluate bone adaptation in the tibia after mechanical loading treatment. However, there is a prevailing assumption that the load is applied axially to the tibia. The aim of this in silico study was to evaluate how much the apparent mechanical properties of the mouse tibia are affected by the loading direction, by using a validated micro-finite element (micro-FE) model of mice which have been ovariectomized and exposed to external mechanical loading over a two-week period. Methods: Longitudinal micro-computed tomography (micro-CT) images were taken of the tibiae of eleven ovariectomized mice at ages 18 and 20 weeks. Six of the mice underwent a mechanical loading treatment at age 19 weeks. Micro-FE models were generated, based on the segmented micro-CT images. Three models using unitary loads were linearly combined to simulate a range of loading directions, generated as a function of the angle from the inferior-superior axis (θ, 0°-30° range, 5° steps) and the angle from the anterior-posterior axis (ϕ, 0°: anterior axis, positive anticlockwise, 0°-355° range, 5° steps). The minimum principal strain was calculated and used to estimate the failure load, by linearly scaling the strain until 10% of the nodes reached the critical strain level of -14,420 με. The apparent bone stiffness was calculated as the ratio between the axial applied force and the average displacement along the longitudinal direction, for the loaded nodes. Results: The results demonstrated a high sensitivity of the mouse tibia to the loading direction across all groups and time points. Higher failure loads were found for several loading directions (θ = 10°, ϕ 205°-210°) than for the nominal axial case (θ = 0°, ϕ = 0°), highlighting adaptation of the bone for loading directions far from the nominal axial one. Conclusion: These results suggest that in studies which use mouse tibia, the loading direction can significantly impact the failure load. Thus, the magnitude and direction of the applied load should be well controlled during the experiments.

Pathophysiological mechanism of acute bone loss after fracture

BACKGROUND

Acute bone loss after fracture is associated with various effects on the complete recovery process and a risk of secondary fractures among patients. Studies have reported similarities in pathophysiological mechanisms involved in acute bone loss after fractures and osteoporosis. However, given the silence nature of bone loss and bone metabolism complexities, the actual underlying pathophysiological mechanisms have yet to be fully elucidated.

AIM OF REVIEW

To elaborate the latest findings in basic research with a focus on acute bone loss after fracture. To briefly highlight potential therapeutic targets and current representative drugs. To arouse researchers' attention and discussion on acute bone loss after fracture.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Bone loss after fracture is associated with immobilization, mechanical unloading, blood supply damage, sympathetic nerve regulation, and crosstalk between musculoskeletals among other factors. Current treatment strategies rely on regulation of osteoblasts and osteoclasts, therefore, there is a need to elucidate on the underlying mechanisms of acute bone loss after fractures to inform the development of efficacious and safe drugs. In addition, attention should be paid towards ensuring long-term skeletal health.

Mechanostat parameters estimated from time-lapsed in vivo micro-computed tomography data of mechanically driven bone adaptation are logarithmically dependent on loading frequency

Mechanical loading is a key factor governing bone adaptation. Both preclinical and clinical studies have demonstrated its effects on bone tissue, which were also notably predicted in the mechanostat theory. Indeed, existing methods to quantify bone mechanoregulation have successfully associated the frequency of (re)modeling events with local mechanical signals, combining time-lapsed in vivo micro-computed tomography (micro-CT) imaging and micro-finite element (micro-FE) analysis. However, a correlation between the local surface velocity of (re)modeling events and mechanical signals has not been shown. As many degenerative bone diseases have also been linked to impaired bone (re)modeling, this relationship could provide an advantage in detecting the effects of such conditions and advance our understanding of the underlying mechanisms. Therefore, in this study, we introduce a novel method to estimate (re)modeling velocity curves from time-lapsed in vivo mouse caudal vertebrae data under static and cyclic mechanical loading. These curves can be fitted with piecewise linear functions as proposed in the mechanostat theory. Accordingly, new (re)modeling parameters can be derived from such data, including formation saturation levels, resorption velocity moduli, and (re)modeling thresholds. Our results revealed that the norm of the gradient of strain energy density yielded the highest accuracy in quantifying mechanoregulation data using micro-finite element analysis with homogeneous material properties, while effective strain was the best predictor for micro-finite element analysis with heterogeneous material properties. Furthermore, (re)modeling velocity curves could be accurately described with piecewise linear and hyperbola functions (root mean square error below 0.2 µm/day for weekly analysis), and several (re)modeling parameters determined from these curves followed a logarithmic relationship with loading frequency. Crucially, (re)modeling velocity curves and derived parameters could detect differences in mechanically driven bone adaptation, which complemented previous results showing a logarithmic relationship between loading frequency and net change in bone volume fraction over 4 weeks. Together, we expect this data to support the calibration of in silico models of bone adaptation and the characterization of the effects of mechanical loading and pharmaceutical treatment interventions in vivo.

Using Finite Element Modeling in Bone Mechanoadaptation

PURPOSE OF THE REVIEW

Bone adapts structure and material properties in response to its mechanical environment, a process called mechanoadpatation. For the past 50 years, finite element modeling has been used to investigate the relationships between bone geometry, material properties, and mechanical loading conditions. This review examines how we use finite element modeling in the context of bone mechanoadpatation.

RECENT FINDINGS

Finite element models estimate complex mechanical stimuli at the tissue and cellular levels, help explain experimental results, and inform the design of loading protocols and prosthetics. FE modeling is a powerful tool to study bone adaptation as it complements experimental approaches. Before using FE models, researchers should determine whether simulation results will provide complementary information to experimental or clinical observations and should establish the level of complexity required. As imaging technics and computational capacity continue increasing, we expect FE models to help in designing treatments of bone pathologies that take advantage of mechanoadaptation of bone.

Effect of stretch frequency on osteogenesis of periodontium during periodontal ligament distraction

OBJECTIVES

Periodontal ligament distraction (PDLD) can accelerate orthodontic tooth movement (OTM). However, the effect of stretch frequency on osseous formation during PDLD remains unclear. Here, we sought to identify the effect of PDLD frequency on the osteogenic remodelling of the periodontium.

MATERIALS AND METHODS

(i) In vitro, five human periodontal ligament stem cell (PDLSC) cultures were randomized to either static conditions or exposure to a cyclic stretch force involving 12% deformation at frequencies of 0.3, 0.5, 0.7 or 1.0 Hz for 12 h, and the osteogenic differentiation of PDLSCs was assessed using Western blotting. (ii) In vivo, 18 beagle dogs underwent orthodontic distalization of bilateral maxillary first premolars. In the test groups, PDLD was performed at a frequency of two or six times/day, while Ni-Ti coil springs were applied to mimic traditional OTM in the control group. The amount of OTM and histological staining was estimated after force loading for 5, 10 and 15 days.

RESULTS

(i) In vitro, the expression of osteogenic-specific markers (runt-related transcription factor 2 [Runx2], type I collagen [COL-I] and osteocalcin [OCN]) increased with the frequency of tensile force, to a peak at 0.7 Hz. (ii) In vivo, both PDLD groups displayed a greater rate of OTM and a higher bone metabolism than the control group. The expression of COL-I and OCN was significantly reinforced in the six times/day-PDLD group in comparison to the two times/day-PDLD group.

CONCLUSIONS

The cyclic stretch force enhances osteogenesis of the periodontium in a frequency-dependent manner.

Increased osteoclastogenesis contributes to bone loss in the Costello syndrome Hras G12V mouse model

RAS GTPases are ubiquitous GDP/GTP-binding proteins that function as molecular switches in cellular signalling and control numerous signalling pathways and biological processes. Pathogenic mutations in RAS genes severely affect cellular homeostasis, leading to cancer when occurring in somatic cells and developmental disorders when the germline is affected. These disorders are generally termed as RASopathies and among them Costello syndrome (CS) is a distinctive entity that is caused by specific HRAS germline mutations. The majority of these mutations affect residues 12 and 13, the same sites as somatic oncogenic HRAS mutations. The hallmarks of the disease include congenital cardiac anomalies, impaired thriving and growth, neurocognitive impairments, distinctive craniofacial anomalies, and susceptibility to cancer. Adult patients often present signs of premature aging including reduced bone mineral density and osteoporosis. Using a CS mouse model harbouring a Hras G12V germline mutation, we aimed at determining whether this model recapitulates the patients’ bone phenotype and which bone cells are driving the phenotype when mutated. Our data revealed that Hras G12V mutation induces bone loss in mice at certain ages. In addition, we identified that bone loss correlated with an increased number of osteoclasts in vivo and Hras G12V mutations increased osteoclastogenesis in vitro. Last, but not least, mutant osteoclast differentiation was reduced by treatment in vitro with MEK and PI3K inhibitors, respectively. These results indicate that Hras is a novel regulator of bone homeostasis and an increased osteoclastogenesis due to Hras G12V mutation contributes to bone loss in the Costello syndrome.

Controlled mechanical early loads improve bone quality and quantity around implants: An in vivo experimental study

OBJECTIVES

This study aimed to investigate the effects of early loads on bone quality and quantity around implants and to compare the effects of early loads on bone quality and quantity with the effects of conventional loads.

MATERIALS AND METHODS

Grade IV-titanium implants with buttress threads were placed in rat maxillary bone 4 weeks after extraction of first molars. A controlled mechanical load (10 N, 3 Hz, 1,800 cycles, 2 days/week) was started via the implants 1 and 3 weeks after implant placement for 2 weeks (early and conventional loads, respectively). Bone quality, defined as distribution of bone cells, types and orientation of collagen fibers, and production of semaphorin3A, its receptor neuropilin-1, and sclerostin, were quantitatively evaluated.

RESULTS

Early loads substantially and positively affected bone quality by changing the preferential alignment of collagen fibers with increased production of type I and III collagens, semaphorin3A, and neuropilin-1, increased osteoblast numbers, decreased production of sclerostin, and decreased osteoclast numbers both inside and outside the implant threads, when compared with non-loaded conditions. Conventional loads changed bone quality around implants slightly. Interestingly, early loads had significantly stronger effects on bone quality and quantity based on the evaluation parameters than conventional loads.

CONCLUSIONS

This is the first report to provide scientific evidence for load initiation time based on both bone quality and quantity around implants. These new findings show that implants with buttress threads transmitted early loads optimally to bone tissue by improving bone quality and quantity inside and outside the implant threads.

Mechanoregulated trabecular bone adaptation: Progress report on in silico approaches

Adaptation is the process by which bone responds to changes in loading environment and modulates its properties and spatial organization to meet the mechanical demands. Adaptation in trabecular bone is achieved through increase in bone mass and alignment of trabecular-bone morphology along the loading direction. This transformation of internal microstructure is governed by mechanical stimuli sensed by mechanosensory cells in the bone matrix. Realisation of adaptation in the form of local bone-resorption and -formation activities as a function of mechanical stimuli is still debated. In silico modelling is a useful tool for simulation of various scenarios that cannot be investigated in vivo and particularly well suited for prediction of trabecular bone adaptation. This progress report presents the recent advances in in silico modelling of mechanoregulated adaptation at the scale of trabecular bone tissue. Four well-established bone-adaptation models are reviewed in terms of their recent improvements and validation. They consider various mechanical factors: (i) strain energy density, (ii) strain and damage, (iii) stress nonuniformity and (iv) daily stress. Contradictions of these models are discussed and their ability to describe adequately a real-life mechanoregulation process in bone is compared.

Effect of Tensile Frequency on the Osteogenic Differentiation of Periodontal Ligament Stem Cells

Purpose

The role of periodontal ligament stem cells (PDLSCs) in mediating osteogenesis involved in orthodontic tooth movement (OTM) is well established. However, various relevant in vitro studies vary in the frequency of tension. The effect of tensile frequency on the mechanotransduction of PDLSCs is not clear. The current study aimed to determine the effect of different tensile frequencies on the osteogenic differentiation of PDLSCs and to identify important mechano-sensitivity genes.

Methods

Human PDLSCs were isolated, identified, and subjected to cyclic equibiaxial tensile strain of 12% at different frequencies of 0.1 Hz, 0.5 Hz, 0.7 Hz, or static cultures. Osteogenic differentiation of PDLSCs was assessed by using Western blotting. High-throughput sequencing was used to identify differential mRNA expression. Short time-series expression miner (STEM) was utilized to describe the frequency patterns of the mRNAs. The functions and enriched pathways were identified, and the hub genes were identified and validated.

Results

We found that the osteoblastic differentiation capacity of PDLSCs increased with tensile frequency in the range of 0.1–0.7 Hz. Eight frequency-tendency gene expression profiles were identified to be statistically significant. Tensile frequency-specific expressed genes, such as SALL1 and EYA1, which decreased with the increase in tensile frequency, were found.

Conclusion

The osteoblastic differentiation of PDLSCs under mechanical tensile force is frequency dependent. EYA1 and SALL1 were identified as potential important tensile frequency-sensitive genes, which may contribute to the cyclic tension-induced osteogenic differentiation of PDLSCs in a frequency-dependent manner.

Image-based radiodensity profilometry measures early remodeling at the bone-callus interface in sheep

Bone healing has been traditionally described as a four-phase process: inflammatory response, soft callus formation, hard callus development, and remodeling. The remodeling phase has been largely neglected in most numerical mechanoregulation models of fracture repair in favor of capturing early healing using a pre-defined callus domain. However, in vivo evidence suggests that remodeling occurs concurrently with repair and causes changes in cortical bone adjacent to callus that are typically neglected in numerical models of bone healing. The objective of this study was to use image processing techniques to quantify this early-stage remodeling in ovine osteotomies. To accomplish this, we developed a numerical method for radiodensity profilometry with optimization-based curve fitting to mathematically model the bone density gradients in the radial direction across the cortical wall and callus. After assessing data from 26 sheep, we defined a dimensionless density fitting function that revealed significant remodeling occurring in the cortical wall adjacent to callus during early healing, a 23% average reduction in density compared to intact. This fitting function is robust for modeling radial density gradients in both intact bone and fracture repair scenarios and can capture a wide variety of the healing responses. The fitting function can also be scaled easily for comparison to numerical model predictions and may be useful for validating future mechanoregulatory models of coupled fracture repair and remodeling.

Effect of External Mechanical Stimuli on Human Bone: a narrative review

Bone is a living composite material that has the capacity to adapt and respond to both internal and external stimuli. This capacity allows bone to adapt its structure to habitual loads and repair microdamage. Although human bone evolved to adapt to normal physiologic loading (for example from gravitational and muscle forces), these same biological pathways can potentially be activated through other types of external stimuli such as pulsed electromagnetic fields, mechanical vibration, and others. This review summarizes what is currently known about how human bone adapts to various types of external stimuli. We highlight how studies on sports-specific athletes and other exercise interventions have clarified the role of mechanical loading on bone structure. We also discuss clinical scenarios, such as spinal cord injury, where mechanical loading is drastically reduced, leading to rapid bone loss and permanent alterations to bone structure. Finally, we highlight areas of emerging research and unmet clinical need.

Positive interactions of mechanical loading and PTH treatments on spatio-temporal bone remodelling

Osteoporosis is one of the most common skeletal diseases, but current therapies are limited to generalized antiresorptive or anabolic interventions, which do not target regions that would benefit from improvements to skeletal health. To improve the evaluation of treatment plans, we used a spatio-temporal multiscale approach that combines longitudinal in vivo micro-computed tomography (micro-CT) and in silico subject-specific finite element modeling to quantitatively map bone adaptation changes due to disease and treatment at high resolution. Our findings show time and region-dependent modifications in bone remodelling following one and two sets of mechanical loading and/or pharmacological interventions. The multiscale results highlighted that the distal section was unaffected by mechanical loading alone but the proximal tibia had the greatest gain from positive interactions of combined therapies. Mechanical loading abated the catabolic effect of PTH, but the main benefit of combined treatments occurred from the additive interactions of the two therapies in periosteal apposition. These results provide detailed insight into the efficacy of combined treatments, facilitating the optimisation of dosage and treatment duration in preclinical mouse studies, and the development of novel interventions for skeletal diseases. STATEMENT OF SIGNIFICANCE: Combined mechanical loading and pharmacotherapy have the potential to slow osteoporosis-induced bone loss but current therapies do not target the regions in need of strengthening. We show for the first time spatial region-dependant interactions between PTH and mechanical loading treatment in OVX mouse tibiae, highlighting local regions in the tibia that benefitted from separate and combined treatments. Combined experimental-computational analysis also detailed the lasting period of each treatment per location in the tibia, the extent of positive (or negative) interactions of the combined therapies, and the impact of each treatment on the regulation of bone adaptation spatio-temporally. This approach can be used to create hypothesis about the interactions of different treatments to optimise the design of biomaterials and medical interventions.

The micro and ultrastructural anatomy of bone spicules found in the osteochondral junction of bovine patellae with early joint degeneration

The structural changes in the tissues of the osteochondral junction are a topic of interest, especially considering how bone changes are involved in the initiation and progression of osteoarthritis (OA). Our research group has previously demonstrated that at the cement line boundary between the zone of calcified cartilage (ZCC) and the subchondral bone, in mature bovine patellae with early OA, there are numerous bone spicules that have emerged from the underlying bone. These spicules contain a central vascular canal and a bone cuff. In this study, we use high-resolution differential interference contrast optical microscopy and scanning electron microscopy to compare the cartilage-bone junction of three groups of mature bovine patellae showing healthy to mild to moderately degenerate cartilage. The ZCC and bone junction was carefully examined to estimate the frequency of marrow spaces, bone spicules and fully formed bone bulges. The results reveal that bone spicules are associated with all grades of cartilage tissue studied, with the most occurring in the intermediate stages of tissue health. The micro and ultrastructure of the bone spicule are consistent with that of an osteon, especially those found in compression zones in long bones. Also considering the coexistence of marrow spaces and fully formed bone, this study suggests that these bone spicules arise similar to the formation of osteons in the bone remodelling process. The significance of this conclusion is in the way researchers approach the bone formation issue in the early degenerative joint. Instead of endochondral ossification, we propose that bone formation in OA is more akin to a combination of primary bone remodelling and de novo bone formation.

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.

Bone remodelling patterns around orthodontic mini-implants migrating in bone: an experimental study in rat vertebrae

Background

Orthodontic implant migration has been clinically observed in presence of continuous loading forces. Recent studies indicate that osteocytes play a crucial role in this phenomenon.

Objectives

Aim of this study was to investigate local osteocytic gene expression, protein expression, and bone micro-structure in peri-implant regions of pressure and tension.

Material and methods

The present work reports a complementary analysis to a previous micro-computed tomography study. Two customized mini-implants were placed in one caudal rat vertebra and connected by a nickel–titanium contraction spring generating different forces (i.e. 0, 0.5, 1.0, and 1.5 N). Either at 2 or 8 weeks, the vertebrae were harvested and utilized for 1. osteocytic gene expression using laser capture micro-dissection on frozen sections coupled with qPCR, 2. haematoxylin–eosin staining for qualitative and quantitative analyses, 3. immunofluorescence staining and analysis, and 4. bone-to-implant contact on undecalcified samples.

Results

At the two time points for all the performed analyses no significant differences were observed with respect to the applied force magnitudes and cell harvesting localization. However, descriptive histological analysis revealed remarkable bone remodelling at 2 weeks of loading. At 8 weeks the implants were osseointegrated and, especially in 1.0 and 1.5 N groups, newly formed bone presented a characteristic load bearing architecture with trabecula oriented in the direction of the loading.

Conclusions

The present study confirmed that stress-induced bone remodelling is the biological mechanism of orthodontic implant migration. Bone apposition was found at ‘tension’ and ‘pressure’ sites thus limiting implant migration over time.

The mechanical influence of bone spicules in the osteochondral junction: A finite element modelling study

While much has been done to study how cartilage responds to mechanical loading, as well as modelling such responses, arguably less has been accomplished around the mechanics of the cartilage-bone junction. Previously, it has been reported that the presence of bony spicules invading the zone of calcified cartilage, preceded the formation of new subchondral bone and the advancing of the cement line (Thambyah and Broom in Osteoarthr Cartil 17:456-463, 2009). In this study, the morphology and frequency of bone spicules in the cartilage-bone interface of osteochondral beams subjected to three-point bending were modelled, and the results are discussed within the context of biomechanical theories on bone formation. It was found that the stress and strain magnitudes, and their distribution were sensitive to the presence and number of spicules. Spicule numbers and shape were shown to affect the strain energy density (SED) distribution in the areas of the cement line adjacent to spicules. Stresses, strains and SED analyses thus provided evidence that the mechanical environment with the addition of spicules promotes bone formation in the cartilage-bone junction.

TRPV4 mediates cell damage induced by hyperphysiological compression and regulates COX2/PGE2 in intervertebral discs

BACKGROUND

Aberrant mechanical loading of the spine causes intervertebral disc (IVD) degeneration and low back pain. Current therapies do not target the mediators of the underlying mechanosensing and mechanotransduction pathways, as these are poorly understood. This study investigated the role of the mechanosensitive transient receptor potential vanilloid 4 (TRPV4) ion channel in dynamic compression of bovine nucleus pulposus (NP) cells in vitro and mouse IVDs in vivo.

METHODS

Degenerative changes and the expression of the inflammatory mediator cyclooxygenase 2 (COX2) were examined histologically in the IVDs of mouse tails that were dynamically compressed at a short repetitive hyperphysiological regime (vs sham). Bovine NP cells embedded in an agarose-collagen hydrogel were dynamically compressed at a hyperphysiological regime in the presence or absence of the selective TRPV4 antagonist GSK2193874. Lactate dehydrogenase (LDH) and prostaglandin E2 (PGE2) release, as well as phosphorylation of mitogen-activated protein kinases (MAPKs), were analyzed. Degenerative changes and COX2 expression were further evaluated in the IVDs of trpv4-deficient mice (vs wild-type; WT).

RESULTS

Dynamic compression caused IVD degeneration in vivo as previously shown but did not affect COX2 expression. Dynamic compression significantly augmented LDH and PGE2 releases in vitro, which were significantly reduced by TRPV4 inhibition. Moreover, TRPV4 inhibition during dynamic compression increased the activation of the extracellular signal-regulated kinases 1/2 (ERK) MAPK pathway by 3.13-fold compared to non-compressed samples. Trpv4-deficient mice displayed mild IVD degeneration and decreased COX2 expression compared to WT mice.

CONCLUSIONS

TRPV4 therefore regulates COX2/PGE2 and mediates cell damage induced by hyperphysiological dynamic compression, possibly via ERK. Targeted TRPV4 inhibition or knockdown might thus constitute promising therapeutic approaches to treat patients suffering from IVD pathologies caused by aberrant mechanical stress.

Bone Mechanoregulation Allows Subject-Specific Load Estimation Based on Time-Lapsed Micro-CT and HR-pQCT in Vivo

Patients at high risk of fracture due to metabolic diseases frequently undergo long-term antiresorptive therapy. However, in some patients, treatment is unsuccessful in preventing fractures or causes severe adverse health outcomes. Understanding load-driven bone remodelling, i.e., mechanoregulation, is critical to understand which patients are at risk for progressive bone degeneration and may enable better patient selection or adaptive therapeutic intervention strategies. Bone microarchitecture assessment using high-resolution peripheral quantitative computed tomography (HR-pQCT) combined with computed mechanical loads has successfully been used to investigate bone mechanoregulation at the trabecular level. To obtain the required mechanical loads that induce local variances in mechanical strain and cause bone remodelling, estimation of physiological loading is essential. Current models homogenise strain patterns throughout the bone to estimate load distribution in vivo, assuming that the bone structure is in biomechanical homoeostasis. Yet, this assumption may be flawed for investigating alterations in bone mechanoregulation. By further utilising available spatiotemporal information of time-lapsed bone imaging studies, we developed a mechanoregulation-based load estimation (MR) algorithm. MR calculates organ-scale loads by scaling and superimposing a set of predefined independent unit loads to optimise measured bone formation in high-, quiescence in medium-, and resorption in low-strain regions. We benchmarked our algorithm against a previously published load history (LH) algorithm using synthetic data, micro-CT images of murine vertebrae under defined experimental in vivo loadings, and HR-pQCT images from seven patients. Our algorithm consistently outperformed LH in all three datasets. In silico-generated time evolutions of distal radius geometries (n = 5) indicated significantly higher sensitivity, specificity, and accuracy for MR than LH (p < 0.01). This increased performance led to substantially better discrimination between physiological and extra-physiological loading in mice (n = 8). Moreover, a significantly (p < 0.01) higher association between remodelling events and computed local mechanical signals was found using MR [correct classification rate (CCR) = 0.42] than LH (CCR = 0.38) to estimate human distal radius loading. Future applications of MR may enable clinicians to link subtle changes in bone strength to changes in day-to-day loading, identifying weak spots in the bone microstructure for local intervention and personalised treatment approaches.

Cortical Thickness Adaptive Response to Mechanical Loading Depends on Periosteal Position and Varies Linearly With Loading Magnitude

The aim of the current study was to quantify the local effect of mechanical loading on cortical bone formation response at the periosteal surface using previously obtained μCT data from a mouse tibia mechanical loading study. A novel image analysis algorithm was developed to quantify local cortical thickness changes (ΔCt.Th) along the periosteal surface due to different peak loads (0N ≤ F ≤ 12N) applied to right-neurectomised mature female C57BL/6 mice. Furthermore, beam analysis was performed to analyse the local strain distribution including regions of tensile, compressive, and low strain magnitudes. Student's paired t-test showed that ΔCt.Th in the proximal (25%), proximal/middle (37%), and middle (50%) cross-sections (along the z-axis of tibia) is strongly associated with the peak applied loads. These changes are significant in a majority of periosteal positions, in particular those experiencing high compressive or tensile strains. No association between F and ΔCt.Th was found in regions around the neutral axis. For the most distal cross-section (75%), the association of loading magnitude and ΔCt.Th was not as pronounced as the more proximal cross-sections. Also, bone formation responses along the periosteum did not occur in regions of highest compressive and tensile strains predicted by beam theory. This could be due to complex experimental loading conditions which were not explicitly accounted for in the mechanical analysis. Our results show that the bone formation response depends on the load magnitude and the periosteal position. Bone resorption due to the neurectomy of the loaded tibia occurs throughout the entire cross-sectional region for all investigated cortical sections 25, 37, 50, and 75%. For peak applied loads higher than 4 N, compressive and tensile regions show bone formation; however, regions around the neutral axis show constant resorption. The 50% cross-section showed the most regular ΔCt.Th response with increased loading when compared to 25 and 37% cross-sections. Relative thickness gains of approximately 70, 60, and 55% were observed for F = 12 N in the 25, 37, and 50% cross-sections. ΔCt.Th at selected points of the periosteum follow a linear response with increased peak load; no lazy zone was observed at these positions.

Application of subject-specific adaptive mechanical loading for bone healing in a mouse tail vertebral defect

Methods to repair bone defects arising from trauma, resection, or disease, continue to be sought after. Cyclic mechanical loading is well established to influence bone (re)modelling activity, in which bone formation and resorption are correlated to micro-scale strain. Based on this, the application of mechanical stimulation across a bone defect could improve healing. However, if ignoring the mechanical integrity of defected bone, loading regimes have a high potential to either cause damage or be ineffective. This study explores real-time finite element (rtFE) methods that use three-dimensional structural analyses from micro-computed tomography images to estimate effective peak cyclic loads in a subject-specific and time-dependent manner. It demonstrates the concept in a cyclically loaded mouse caudal vertebral bone defect model. Using rtFE analysis combined with adaptive mechanical loading, mouse bone healing was significantly improved over non-loaded controls, with no incidence of vertebral fractures. Such rtFE-driven adaptive loading regimes demonstrated here could be relevant to clinical bone defect healing scenarios, where mechanical loading can become patient-specific and more efficacious. This is achieved by accounting for initial bone defect conditions and spatio-temporal healing, both being factors that are always unique to the patient.

Finite element analysis of human skull bone adaptation to mechanical loading

Bones self-optimize their mechanical behavior in response to mechanical stimulus. The objective of this research was to develop an integrated bone remodeling and stress binning algorithms into a finite element environment to elucidate the evolution of the bone properties as a function of loading. The bone remodeling algorithm was used to calculate the change in the density and elastic modulus based on the strain energy stimulus. The stress-binning procedure seeks to assign the properties to each element based on the levels of stress from the previous cycle, eliminating pseudo-lazy-zoning and stress dilation effects. The developed algorithms were used to analyze the response skull to loading associated with orthodontic devices. Specifically, a load was applied between the roots of the canine teeth and the first premolars while constraining the foramen magnum. Full-field contours of the displacement, strain, and strain energy were extracted after each remodeling cycle at nine commonly cephalometric landmarks. The results indicate that the overall mechanical response and the associated properties reached a steady-state behavior after nearly 50 cycles of applying the algorithm, where different zones within the skull exhibited unique evolution based on the locations from the loading and boundary sites. When approaching this steady-state condition, it was found that the upper incisor displacement is reduced by 72%, and the density is reduced by almost 7.5%. The finite element approach can be used in defining the treatment process by dynamically changing the loads. Future research will focus on integrating the time-dependent behavior of the bone.