Hanna Isaksson

Osteoblast-derived WNT16 represses osteoclastogenesis and prevents cortical bone fragility fractures.

The WNT16 locus is a major determinant of cortical bone thickness and nonvertebral fracture risk in humans. The disability, mortality and costs caused by osteoporosis-induced nonvertebral fractures are enormous. We demonstrate here that Wnt16-deficient mice develop spontaneous fractures as a result of low cortical thickness and high cortical porosity. In contrast, trabecular bone volume is not altered in these mice. Mechanistic studies revealed that WNT16 is osteoblast derived and inhibits human and mouse osteoclastogenesis both directly by acting on osteoclast progenitors and indirectly by increasing expression of osteoprotegerin (Opg) in osteoblasts. The signaling pathway activated by WNT16 in osteoclast progenitors is noncanonical, whereas the pathway activated in osteoblasts is both canonical and noncanonical. Conditional Wnt16 inactivation revealed that osteoblast-lineage cells are the principal source of WNT16, and its targeted deletion in osteoblasts increases fracture susceptibility. Thus, osteoblast-derived WNT16 is a previously unreported key regulator of osteoclastogenesis and fracture susceptibility. These findings open new avenues for the specific prevention or treatment of nonvertebral fractures, a substantial unmet medical need.

Sclerostin antibody treatment enhances metaphyseal bone healing in rats

Sclerostin is the product of the SOST gene. Loss‐of‐function mutations in the SOST gene result in a high‐bone‐mass phenotype, demonstrating that sclerostin is a negative regulator of bone mass. Primarily expressed by osteocytes in bone, sclerostin is reported to bind the LRP5/6 receptor, thereby antagonizing canonical Wnt signaling and negatively regulating bone formation. We therefore investigated whether systemic administration of a sclerostin‐neutralizing antibody would increase the regeneration of traumatized metaphyseal bone in rats. Young male rats had a screw inserted in the proximal tibia and were divided into six groups given 25 mg/kg of sclerostin antibody or control twice a week subcutaneously for 2 or 4 weeks. In four groups, the screws were tested for pull‐out strength. At the time of euthanasia, a similar screw also was inserted in the contralateral tibia and pull‐out tested immediately. Sclerostin antibody significantly increased the pull‐out force by almost 50% compared with controls after 2 and 4 weeks. Also, the screws inserted at the time of euthanasia showed increased pull‐out force. Micro–computed tomography (µCT) of the remaining two groups showed that the antibody led to a 30% increase in bone volume fraction in a region surrounding the screw. There also was a general increase in trabecular thickness in cancellous bone. Thus, as measured by the amount of bone and its mechanical resistance, the sclerostin antibody increased bone formation during metaphyseal repair but also in untraumatized bone.

A mechano-regulatory bone-healing model incorporating cell-phenotype specific activity

Phenomenological computational models of tissue regeneration and bone healing have been only partially successful in predicting experimental observations. This may be a result of simplistic modeling of cellular activity. Furthermore, phenomenological models are limited when considering the effects of combined physical and biological interventions. In this study, a new model of cell and tissue differentiation, using a more mechanistic approach, is presented and applied to fracture repair. The model directly couples cellular mechanisms to mechanical stimulation during bone healing and is based on the belief that the cells act as transducers during tissue regeneration. In the model, the cells within the matrix proliferate, differentiate, migrate, and produce extracellular matrix, all at cell-phenotype specific rates, based on the mechanical stimulation they experience. The model is assembled from coupled partial differentiation equations, which are solved using a newly developed finite element formulation. The evolution of four cell types, i.e. mesenchymal stem cells, fibroblasts, chondrocytes and osteoblasts, and the production of extracellular matrices of fibrous tissue, cartilage and bone are calculated. The material properties of the tissues are iteratively updated based on actual amounts of extracellular matrix in material elements at progressive time points. A two-dimensional finite element model of a long bone osteotomy was used to evaluate the models potential. The additional value of the presented model and the importance of including cell-phenotype specific activities when modeling tissue differentiation and bone healing, were demonstrated by comparing the predictions with phenomenological models. The models capacity was established by showing that it can correctly predict several aspects of bone healing, including cell and tissue distributions during normal fracture healing. Furthermore, it was able to predict experimentally established alterations due to excessive mechanical stimulation, periosteal stripping and impaired effects of cartilage remodeling.

Precision of nanoindentation protocols for measurement of viscoelasticity in cortical and trabecular bone

Nanoindentation has recently gained attention as a characterization technique for mechanical properties of biological tissues, such as bone, on the sub-micron level. However, optimal methods to characterize viscoelastic properties of bones are yet to be established. This study aimed to compare the time-dependent viscoelastic properties of bone tissue obtained with different nanoindentation methods. Bovine cortical and trabecular bone samples (n=8) from the distal femur and proximal tibia were dehydrated, embedded and polished. The material properties determined using nanoindentation were hardness and reduced modulus, as well as time-dependent parameters based on creep, loading-rate, dissipated energy and semi-dynamic testing under load control. Each loading protocol was repeated 160 times and the reproducibility was assessed based on the coefficient of variation (CV). Additionally, three well-characterized polymers were tested and CV values were calculated for reference. The employed methods were able to characterize time-dependent viscoelastic properties of bone. However, their reproducibility varied highly (CV 9-40%). The creep constant increased with increasing dwell time. The reproducibility was best with a 30s creep period (CV 18%). The dissipated energy was stable after three repeated load cycles, and the reproducibility improved with each cycle (CV 23%). The viscoelastic properties determined with semi-dynamic test increased with increase in frequency. These measurements were most reproducible at high frequencies (CV 9-10%). Our results indicate that several methods are feasible for the determination of viscoelastic properties of bone material. The high frequency semi-dynamic test showed the highest precision within the tested nanoindentation protocols.

Infrared spectroscopy indicates altered bone turnover and remodeling activity in renal osteodystrophy

Renal osteodystrophy alters metabolic activity and remodeling rate of bone and also may lead to different bone composition. The objective of this study was to characterize the composition of bone in high‐turnover renal osteodystrophy patients by means of Fourier transform infrared spectroscopic imaging (FTIRI). Iliac crest biopsies from healthy bone (n = 11) and patients with renal osteodystrophy (ROD, n = 11) were used in this study. The ROD samples were from patients with hyperparathyroid disease. By using FTIRI, phosphate‐to‐amide I ratio (mineral‐to‐matrix ratio), carbonate‐to‐phosphate ratio, and carbonate‐to‐amide I ratio (turnover rate/remodeling activity), as well as the collagen cross‐link ratio (collagen maturity), were quantified. Histomorphometric analyses were conducted for comparison. The ROD samples showed significantly lower carbonate‐to‐phosphate (p < .01) and carbonate‐to‐amide I (p < .001) ratios. The spatial variation across the trabeculae highlighted a significantly lower degree of mineralization (p < .05) at the edges of the trabeculae in the ROD samples than in normal bone. Statistically significant linear correlations were found between histomorphometric parameters related to bone‐remodeling activity and number of bone cells and FTIRI‐calculated parameters based on carbonate‐to‐phosphate and carbonate‐to‐amide I ratios. Hence the results suggested that FTIRI parameters related to carbonate may be indicative of turnover and remodeling rate of bone.

Sensitivity of tissue differentiation and bone healing predictions to tissue properties

Computational models are employed as tools to investigate possible mechano-regulation pathways for tissue differentiation and bone healing. However, current models do not account for the uncertainty in input parameters, and often include assumptions about parameter values that are not yet established. The aim was to clarify the importance of the assumed tissue material properties in a computational model of tissue differentiation during bone healing. An established mechano-biological model was employed together with a statistical approach. The model included an adaptive 2D finite element model of a fractured long bone. Four outcome criteria were quantified: (1) ability to predict sequential healing events, (2) amount of bone formation at specific time points, (3) total time until healing, and (4) mechanical stability at specific time points. Statistical analysis based on fractional factorial designs first involved a screening experiment to identify the most significant tissue material properties. These seven properties were studied further with response surface methodology in a three-level Box-Behnken design. Generally, the sequential events were not significantly influenced by any properties, whereas rate-dependent outcome criteria and mechanical stability were significantly influenced by Youngs modulus and permeability. Poissons ratio and porosity had minor effects. The amount of bone formation at early, mid and late phases of healing, the time until complete healing and the mechanical stability were all mostly dependent on three material properties; permeability of granulation tissue, Youngs modulus of cartilage and permeability of immature bone. The consistency between effects of the most influential parameters was high. To increase accuracy and predictive capacity of computational models of bone healing, the most influential tissue mechanical properties should be accurately quantified.

The effects of Dickkopf-1 antibody on metaphyseal bone and implant fixation under different loading conditions

The secreted protein Dickkopf-1 (Dkk1) is an antagonist of canonical Wnt signaling, expressed during fracture healing. It is unclear how it is involved in the mechanical control of bone maintenance. We investigated the response to administration of a Dkk1 neutralizing antibody (Dkk1-ab) in metaphyseal bone under different loading conditions, with or without trauma. In this three part experiment, 120 rats had a screw or bone chamber inserted either unilaterally or bilaterally in the proximal tibia. Mechanical (pull-out) testing, μCT and histology were used for evaluation. The animals were injected with either 10mg/kg Dkk1-ab or saline every 14days for 14, 28, or 42days. Antibody treatment increased bone formation around the screws and improved their fixation. After 28days, the pull-out force was increased by over 100%. In cancellous bone, the bone volume fraction was increased by 50%. In some animals, one hind limb was paralyzed with Botulinum toxin A (Botox) to create a mechanically unloaded environment. This did not increase the response to antibody treatment with regard to screw fixation, but in cancellous bone, the bone volume fraction increased by 233%. Thus, the response in unloaded, untraumatized bone was proportionally larger, suggesting that Dkk1 may be up-regulated in unloaded bone. There was also an increase in thickness of the metaphyseal cortex. In bone chambers, the antibody treatment increased the bone volume fraction. The results suggest that antibodies blocking Dkk1 might be used to stimulate bone formation especially during implant fixation, fracture repair, or bone disuse. It also seems that Dkk1 is up-regulated both after metaphyseal trauma and after unloading, and that Dkk1 is involved in mechano-transduction.

Longitudinal elastic properties and porosity of cortical bone tissue vary with age in human proximal femur

Tissue level structural and mechanical properties are important determinants of bone strength. As an individual ages, microstructural changes occur in bone, e.g., trabeculae and cortex become thinner and porosity increases. However, it is not known how the elastic properties of bone change during aging. Bone tissue may lose its elasticity and become more brittle and prone to fractures as it ages. In the present study the age-dependent variation in the spatial distributions of microstructural and microelastic properties of the human femoral neck and shaft were evaluated by using acoustic microscopy. Although these properties may not be directly measured in vivo, there is a major interest to investigate their relationships with the linear elastic measurements obtained by diagnostic ultrasound at the most severe fracture sites, e.g., the femoral neck. However, before the validity of novel in vivo techniques can be established, it is essential to understand the age-dependent variation in tissue elastic properties and porosity at different skeletal sites. A total of 42 transverse cross-sectional bone samples were obtained from the femoral neck (Fn) and proximal femoral shaft (Ps) of 21 men (mean±SD age 47.1±17.8, range 17-82years). Samples were quantitatively imaged using a scanning acoustic microscope (SAM) equipped with a 50MHz ultrasound transducer. Distributions of the elastic coefficient (c33) of cortical (Ct) and trabecular (Tr) tissues and microstructure of cortex (cortical thickness Ct.Th and porosity Ct.Po) were determined. Variations in c33 were observed with respect to tissue type (c33Trc33(Ct.Fn)=35.3GPa>c33(Tr.Ps)=33.8GPa>c33(Tr.Fn)=31.9GPa), and cadaver age (R(2)=0.28-0.46, p<0.05). Regional variations in porosity were found in the neck (superior 13.1%; inferior 6.1%; anterior 10.1%; posterior 8.6%) and in the shaft (medial 9.5%; lateral 7.7%; anterior 8.6%; posterior 12.0%). In conclusion, significant variations in elastic coefficients were detected between femoral neck and shaft as well as between the quadrants of the cross-sections of neck and shaft. Moreover, an age-related increase in cortical porosity and a stiffening of the bone tissue were observed. These findings may explain in part the increase in susceptibility to suffer low energy fractures during aging and highlight the potential of ultrasound in clinical osteoporosis diagnostics.

Remodeling of fracture callus in mice is consistent with mechanical loading and bone remodeling theory

During the remodeling phase of fracture healing in mice, the callus gradually transforms into a double cortex, which thereafter merges into one cortex. In large animals, a double cortex normally does not form. We investigated whether these patterns of remodeling of the fracture callus in mice can be explained by mechanical loading. Morphologies of fractures after 21, 28, and 42 days of healing were determined from an in vivo mid‐diaphyseal femoral osteotomy healing experiment in mice. Bone density distributions from microCT at 21 days were converted into adaptive finite element models. To assess the effect of loading mode on bone remodeling, a well‐established remodeling algorithm was used to examine the effect of axial force or bending moment on bone structure. All simulations predicted that under axial loading, the callus remodeled to form a single cortex. When a bending moment was applied, dual concentric cortices developed in all simulations, corresponding well to the progression of remodeling observed experimentally and resulting in quantitatively comparable callus areas of woven and lamellar bone. Effects of biological differences between species or other reasons cannot be excluded, but this study demonstrates how a difference in loading mode could explain the differences between the remodeling phase in small rodents and larger mammals.

Determining the most important cellular characteristics for fracture healing, using design of experiments methods

Computational models are employed as tools to investigate possible mechanoregulation pathways for tissue differentiation and bone healing. However, current models do not account for the uncertainty in input parameters, and often include assumptions about parameter values that are not yet established. The objective of this study was to determine the most important cellular characteristics of a mechanoregulatory model describing both cell phenotype-specific and mechanobiological processes that are active during bone healing using a statistical approach. The computational model included an adaptive two-dimensional finite element model of a fractured long bone. Three different outcome criteria were quantified: (1) ability to predict sequential healing events, (2) amount of bone formation at early, mid and late stages of healing and (3) the total time until complete healing. For the statistical analysis, first a resolution IV fractional factorial design (L(64)) was used to identify the most significant factors. Thereafter, a three-level Taguchi orthogonal array (L(27)) was employed to study the curvature (non-linearity) of the 10 identified most important parameters. The results show that the ability of the model to predict the sequences of normal fracture healing was predominantly influenced by the rate of matrix production of bone, followed by cartilage degradation (replacement). The amount of bone formation at early stages was solely dependent on matrix production of bone and the proliferation rate of osteoblasts. However, the amount of bone formation at mid and late phases had the rate of matrix production of cartilage as the most influential parameter. The time to complete healing was primarily dependent on the rate of cartilage degradation during endochondral ossification, followed by the rate of cartilage formation. The analyses of the curvature revealed a linear response for parameters related to bone, where higher rates of formation were more beneficial to healing. In contrast, parameters related to fibrous tissue and cartilage showed optimum levels. Some fibrous connective tissue- and cartilage formation was beneficial to bone healing, but too much of either tissue delayed bone formation. The identified significant parameters and processes are further confirmed by in vivo animal experiments in the literature. This study illustrates the potential of design of experiments methods for evaluating computational mechanobiological model parameters and suggests that further experiments should preferably focus at establishing values of parameters related to cartilage formation and degradation.