Quantitative analysis of imprint shape and its relation to mechanical properties measured by microindentation in bone

Microindentation in bone is a micromechanical testing technique routinely used to extract material properties related to bone quality. As the analysis of microindentation data is based on assumptions about the contact between sample and surface, the aim of this study was to quantify the topological variability of indentations in bone and examine its relationship with mechanical properties. Indentations were performed in dry human and ovine bone in axial and transverse directions and their topology was measured by atomic force microscopy. Statistical shape modeling of the residual imprint allowed to define a mean shape and to describe the variability in terms of 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was found to be highly consistent and free of any pile up while differing mostly by depth between species and direction. A few of the topological parameters, in particular depth, showed significant but rather weak and inconsistent correlations to variations in mechanical properties. The mechanical response of bone as well as the residual imprint shape was highly consistent within each category. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small deviations from an ideally flat surface.

The Initial Slope of the Variogram, Foundation of the Trabecular Bone Score, Is Not or Poorly Associated With Vertebral Strength

Trabecular bone score (TBS) rests on the textural analysis of DXA to reflect the decay in trabecular structure characterising osteoporosis. Yet, its discriminative power in fracture studies remains incomprehensible as prior biomechanical tests found no correlation with vertebral strength. To verify this result possibly due to an unrealistic set-up and to cover a wide range of loading scenarios, the data from three previous biomechanical studies using different experimental settings was used. They involved the compressive failure of 62 human lumbar vertebrae loaded 1) via intervertebral discs to mimic the in vivo situation (“full vertebra”), 2) via the classical endplate embedding (“vertebral body”) or 3) via a ball joint to induce anterior wedge failure (“vertebral section”). HR-pQCT scans acquired prior testing were used to simulate anterior-posterior DXA from which areal bone mineral density (aBMD) and the initial slope of the variogram (ISV), the early definition of TBS, were evaluated. Finally, the relation of aBMD and ISV with failure load (Fexp) and apparent failure stress (σexp) was assessed and their relative contribution to a multi-linear model was quantified via ANOVA. We found that, unlike aBMD, ISV did not significantly correlate with Fexp and σexp, except for the “vertebral body” case (r2 = 0.396, p = 0.028). Aside from the “vertebra section” set-up where it explained only 6.4% of σexp (p = 0.037), it brought no significant improvement to aBMD. These results indicate that ISV, a replica of TBS, is a poor surrogate for vertebral strength no matter the testing set-up, which supports the prior observations and raises a fortiori the question of the deterministic factors underlying the statistical relationship between TBS and vertebral fracture risk.

Temperature Prediction Model for Bone Drilling Based on Density Distribution and In Vivo Experiments for Minimally Invasive Robotic Cochlear Implantation

Surgical robots have been proposed ex vivo to drill precise holes in the temporal bone for minimally invasive cochlear implantation. The main risk of the procedure is damage of the facial nerve due to mechanical interaction or due to temperature elevation during the drilling process. To evaluate the thermal risk of the drilling process, a simplified model is proposed which aims to enable an assessment of risk posed to the facial nerve for a given set of constant process parameters for different mastoid bone densities. The model uses the bone density distribution along the drilling trajectory in the mastoid bone to calculate a time dependent heat production function at the tip of the drill bit. Using a time dependent moving point source Green's function, the heat equation can be solved at a certain point in space so that the resulting temperatures can be calculated over time. The model was calibrated and initially verified with in vivo temperature data. The data was collected in minimally invasive robotic drilling of 12 holes in four different sheep. The sheep were anesthetized and the temperature elevations were measured with a thermocouple which was inserted in a previously drilled hole next to the planned drilling trajectory. Bone density distributions were extracted from pre-operative CT data by averaging Hounsfield values over the drill bit diameter. Post-operative [Formula: see text]CT data was used to verify the drilling accuracy of the trajectories. The comparison of measured and calculated temperatures shows a very good match for both heating and cooling phases. The average prediction error of the maximum temperature was less than 0.7 °C and the average root mean square error was approximately 0.5 °C. To analyze potential thermal damage, the model was used to calculate temperature profiles and cumulative equivalent minutes at 43 °C at a minimal distance to the facial nerve. For the selected drilling parameters, temperature elevation profiles and cumulative equivalent minutes suggest that thermal elevation of this minimally invasive cochlear implantation surgery may pose a risk to the facial nerve, especially in sclerotic or high density mastoid bones. Optimized drilling parameters need to be evaluated and the model could be used for future risk evaluation.

Molecular and cellular basis of bone resorption.

Osteoclast research has an exciting history and a challenging future. More than 3 decades ago, it became evident that bone-resorbing osteoclasts are of hematopoietic origin and are ultimately linked to the "basic multicellular unit," where they team up with the other cell types, including bone-forming osteoblasts. Since 2 decades, we have learned about the signaling pathways controlling genes relevant for osteoclastogenesis and bone resorption. It took another decade until the hypothesized "osteoclast differentiation" factor was discovered and was translated into an approved pharmacologic strategy. Here, the focus is on another molecular target, cathepsin K, a cysteine protease being released by the osteoclast into the resorption compartment. Genetic deletion and pharmacological blocking of cathepsin K reduces bone resorption but with ongoing bone formation. This observation not only holds great promise to become a new pharmacologic strategy, but it also provides new insights into the coordinated work of cells in the "basic multicellular unit" and thus, bridges the history and future of osteoclast research. This article is a short primer on osteoclast biology for readers of the special issue on odanacatib, a cathepsin K inhibitor.

Loss of bone strength in HLA-B27 transgenic rats is characterized by a high bone turnover and is mainly osteoclast-driven

Objective Although osteopenia is frequent in spondyloarthritis (SpA), the underlying cellular mechanisms and association with other symptoms are poorly understood. This study aimed to characterize bone loss during disease progression, determine cellular alterations, and assess the contribution of inflammatory bowel disease (IBD) to bone loss in HLA-B27 transgenic rats. Methods Bones of 2-, 6-, and 12-month-old non-transgenic, disease-free HLA-B7 and disease-associated HLA-B27 transgenic rats were examined using peripheral quantitative computed tomography, μCT, and nanoindentation. Cellular characteristics were determined by histomorphometry and ex vivo cultures. The impact of IBD was determined using [21-3 x 283-2]F1 rats, which develop arthritis and spondylitis, but not IBD. Results HLA-B27 transgenic rats continuously lost bone mass with increasing age and had impaired bone material properties, leading to a 3-fold decrease in bone strength at 12 months of age. Bone turnover was increased in HLA-B27 transgenic rats, as evidenced by a 3-fold increase in bone formation and a 6-fold increase in bone resorption parameters. Enhanced osteoclastic markers were associated with a larger number of precursors in the bone marrow and a stronger osteoclastogenic response to RANKL or TNFα. Further, IBD-free [21-3 x 283-2]F1 rats also displayed decreased total and trabecular bone density. Conclusions HLA-B27 transgenic rats lose an increasing amount of bone density and strength with progressing age, which is primarily mediated via increased bone remodeling in favor of bone resorption. Moreover, IBD and bone loss seem to be independent features of SpA in HLA-B27 transgenic rats.

Continuum damage interactions between tension and compression in osteonal bone

Skeletal diseases such as osteoporosis impose a severe socio-economic burden to ageing societies. Decreasing mechanical competence causes a rise in bone fracture incidence and mortality especially after the age of 65 y. The mechanisms of how bone damage is accumulated under different loading modes and its impact on bone strength are unclear. We hypothesise that damage accumulated in one loading mode increases the fracture risk in another. This study aimed at identifying continuum damage interactions between tensile and compressive loading modes. We propose and identify the material constants of a novel piecewise 1D constitutive model capable of describing the mechanical response of bone in combined tensile and compressive loading histories. We performed several sets of loading–reloading experiments to compute stiffness, plastic strains, and stress-strain curves. For tensile overloading, a stiffness reduction (damage) of 60% at 0.65% accumulated plastic strain was detectable as stiffness reduction of 20% under compression. For compressive overloading, 60% damage at 0.75% plastic strain was detectable as a stiffness reduction of 50% in tension. Plastic strain at ultimate stress was the same in tension and compression. Compression showed softening and tension exponential hardening in the post-yield regime. The hardening behaviour in compression is unaffected by a previous overload in tension but the hardening behaviour in tension is affected by a previous overload in compression as tensile reloading strength is significantly reduced. This paper demonstrates how damage accumulated under one loading mode affects the mechanical behaviour in another loading mode. To explain this and to illustrate a possible implementation we proposed a theoretical model. Including such loading mode dependent damage and plasticity behaviour in finite element models will help to improve fracture risk analysis of whole bones and bone implant structures.

Experimental validation of a nonlinear μ FE model based on cohesive-frictional plasticity for trabecular bone

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of μFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

Comparison of Mixed and Kinematic Uniform Boundary Conditions in Homogenized Elasticity of Femoral Trabecular Bone Using Microfinite Element Analyses

Mechanical properties of human trabecular bone play an important role in age-related bone fragility and implant stability. Micro-finite element (microFE) analysis allows computing the apparent elastic properties of trabecular bone biopsies, but the results depend on the type of applied boundary conditions (BCs). In this study, 167 femoral trabecular cubic biopsies with a side length of 5.3 mm were analyzed using microFE analysis to compare their stiffness systematically with kinematic uniform boundary conditions (KUBCs) and periodicity-compatible mixed uniform boundary conditions (PMUBCs). The obtained elastic constants were then used in the volume fraction and fabric-based orthotropic Zysset-Curnier model to identify their respective model parameters. As expected, PMUBCs lead to more compliant apparent elastic properties than KUBCs, especially in shear. The differences in stiffness decreased with bone volume fraction and mean intercept length. Unlike KUBCs, PMUBCs were sensitive to heterogeneity of the biopsies. The Zysset-Curnier model predicted apparent elastic constants successfully in both cases with adjusted coefficients of determination of 0.986 for KUBCs and 0.975 for PMUBCs. The role of these boundary conditions in finite element analyses of whole bones and bone-implant systems will need to be investigated in future work.

Effect of boundary conditions on yield properties of human femoral trabecular bone

Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7–72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient R2≥0.93. As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.

Investigation of bone inhibition in non-union spinal fusion caused by intervertebral disc cells in vitro

Introduction: Discectomy and spinal fusion is the gold standard for spinal surgery to relieve pain. However, fusion can be hindered for yet unknown reasons that lead to non-fusions with pseudo-arthrose. It is hence appealing to develop biomaterials that can enhance bone formation. Clinical observations indicate that presence of residual intervertebral disc (IVD) tissue might hinder the ossification. We hypothesize that BMP-antagonists are constantly secreted by IVD cells and potentially prevent the ossification process. Furthermore, L51P, the engineered BMP2 variant, stimulates osteoinduction of bone marrow-derived mesenchymal stem cells (MSC) by antagonizing BMP-inhibitors. Methods: Human MSCs, primary nucleus pulposus (NPC) and annulus pulposus cells (AFC) were isolated and expanded in monolayer cultures up to passage 3. IVD cells were seeded in 1.2% alginate beads (4Mio/mL) and separated by culture inserts from MSCs in a co-culture set-up. MSCs were kept in 1:control medium, 2:osteogenic medium+alginate control, 3:osteogenic medium+NPC (±L51P) and 4:osteogenic medium+AFC (±L51P) for 21 days. Relative gene expression of bone-related genes, Alkaline Phosphatase (ALP) assay and histological staining were performed. Results: Osteogenesis of MSCs was hindered as shown by reduced alizarin red staining in the presence of NPC. No such inhibition was observed if co-cultured with alginate only or in the presence of AFC. The results were confirmed on the RNA and protein level. Addition of L51P to the co-cultures induced mineralization of MSCs, however a reduced ALP was observed. Conclusion: We demonstrated that NPC secrete BMP-antagonists that prevent osteogenesis of MSCs and L51P can antagonize BMP-antagonists and induce bone formation.

Prediction of Trabecular Bone Anisotropy from Quantitative Computed Tomography using Supervised Learning and a Novel Morphometric Feature Descriptor

Patient-specific biomechanical models including local bone mineral density and anisotropy have gained importance for assessing musculoskeletal disorders. However the trabecular bone anisotropy captured by high-resolution imaging is only available at the peripheral skeleton in clinical practice. In this work, we propose a supervised learning approach to predict trabecular bone anisotropy that builds on a novel set of pose invariant feature descriptors. The statistical relationship between trabecular bone anisotropy and feature descriptors were learned from a database of pairs of high resolution QCT and clinical QCT reconstructions. On a set of leave-one-out experiments, we compared the accuracy of the proposed approach to previous ones, and report a mean prediction error of 6% for the tensor norm, 6% for the degree of anisotropy and 19◦ for the principal tensor direction. These findings show the potential of the proposed approach to predict trabecular bone anisotropy from clinically available QCT images.

Biomechanical Role of Bone Anisotropy Estimated on Clinical CT Scans by Image Registration

Image-based modeling is a popular approach to perform patient-specific biomechanical simulations. Accurate modeling is critical for orthopedic application to evaluate implant design and surgical planning. It has been shown that bone strength can be estimated from the bone mineral density (BMD) and trabecular bone architecture. However, these findings cannot be directly and fully transferred to patient-specific modeling since only BMD can be derived from clinical CT. Therefore, the objective of this study was to propose a method to predict the trabecular bone structure using a µCT atlas and an image registration technique. The approach has been evaluated on femurs and patellae under physiological loading. The displacement and ultimate force for femurs loaded in stance position were predicted with an error of 2.5% and 3.7%, respectively, while predictions obtained with an isotropic material resulted in errors of 7.3% and 6.9%. Similar results were obtained for the patella, where the strain predicted using the registration approach resulted in an improved mean squared error compared to the isotropic model. We conclude that the registration of anisotropic information from of a single template bone enables more accurate patient-specific simulations from clinical image datasets than isotropic model.

Osteogenic differentiation of bone marrow stromal cells is hindered by the presence of intervertebral disc cells.

BACKGROUND Clinical observations indicate that the presence of nucleus pulposus (NP) tissue during spinal fusion hinders the rate of disc ossification. While the underlying mechanism remains unknown, this observation could be due to incomplete removal of NP cells (NPCs) that secrete factors preventing disc calcification, such as bone morphogenetic protein (BMP) antagonists including noggin and members of the DAN (differential screening selected gene aberrative in neuroblastoma) family. METHODS Monolayer human bone marrow-derived mesenchymal stem cells (MSCs) were cocultured withNPCs and annulus fibrosus cells (AFCs) embedded in alginate for 21 days. At the end of coculture, MSCs were stained for mineral deposition by alizarin red, and relative expression of bone-related genes [Runt-related transcription factor 2, (RUNX2), Osteopontin (OPN), and Alkaline phosphatase (ALP)] and ALP activity were analyzed. Relative expression of three BMP antagonists, chordin (CHRD), gremlin (GREM1), and noggin (NOG), was determined in primary human NPCs and AFCs. These cells were also stained for Gremlin and Noggin by immunocytochemistry. RESULTS Alizarin red staining showed that MSC osteogenesis in monolayer cultures was inhibited by coculture with NPCs or AFCs. ALP activity and RT-PCR analyses confirmed these results and demonstrated inhibition of osteogenesis of MSC in the presence of disc cells. NOG was significantly up-regulated in MSCs after coculture. Relative gene expression of intervertebral disc (IVD) cells showed higher expression of GREM1 in NPCs than in AFCs. CONCLUSIONS We show that primary IVD cells inhibit osteogenesis of MSCs. BMP inhibitors NOG, GREM1 and CHRD were expressed in IVD cells. GREM1 appears to be differentially expressed in NPCs and AFCs. Our results have implications for the design and development of treatments for non-union in spinal fusion.