Comparison of microCT and an inverse finite element approach for biomechanical analysis: results in a mesenchymal stem cell therapeutic system for fracture healing

Novel approaches to correlate computerized tomography imaging of bone fracture callus to callus structural mechanics

Estimating the mechanical properties of bone in vivo without destructive testing would be useful for research and clinical orthopedic applications. Micro-computerized tomography (μCT) imaging can provide quantitative, high-resolution 3D representations of bone morphology and is generally the basis from which bone mechanical properties are non-destructively estimated. The goal of this study was to develop metrics using qualitative and quantitative aspects of bone microarchitecture derived from μCT imaging to estimate the mechanical integrity of bone fracture calluses. Mechanical testing data (peak torque) and μCT image data from 12 rat femur fractures were collected at 4 weeks after fracture. MATLAB was used to analyze the callus μCT imaging data which were then correlated to the empirically determined peak torque of the callus. One metric correlated Z-rays, linear contiguities of voxels running parallel to the neutral axis of the femur and through the fracture callus, to peak torque. Other metrics were based on voxel linkage values (LVs), which is a novel measurement defined by the number of voxels surrounding a given voxel (ranging from 1 to 27) that are all above a specified threshold. Linkage values were utilized to segment the callus and compute healing scores (termed eRUST) based on the modified Radiographic Union Score for Tibial fractures (mRUST). Linkage values were also used to calculate linked bone areas (LBAs). All metrics positively correlated with peak torque, yielding correlations of determination (R2) of 0.863 for eRUST, 0.792 for Z-ray scoring, and 0.764 for a normalized Linked Bone Area metric. These novel metrics appear to be promising approaches for extrapolating fracture callus structural properties from bone microarchitecture using objective analytical methods and without resorting to computationally complex finite element analyses.

The association between mineralised tissue formation and the mechanical local in vivo environment: Time-lapsed quantification of a mouse defect healing model

An improved understanding of how local mechanical stimuli guide the fracture healing process has the potential to enhance clinical treatment of bone injury. Recent preclinical studies of bone defect in animal models have used cross-sectional data to examine this phenomenon indirectly. In this study, a direct time-lapsed imaging approach was used to investigate the local mechanical strains that precede the formation of mineralised tissue at the tissue scale. The goal was to test two hypotheses: 1) the local mechanical signal that precedes the onset of tissue mineralisation is higher in areas which mineralise, and 2) this local mechanical signal is independent of the magnitude of global mechanical loading of the tissue in the defect. Two groups of mice with femoral defects of length 0.85 mm (n = 10) and 1.45 mm (n = 9) were studied, allowing for distinct distributions of tissue scale strains in the defects. The regeneration and (re)modelling of mineralised tissue was observed weekly using in vivo micro-computed tomography (micro-CT), which served as a ground truth for resolving areas of mineralised tissue formation. The mechanical environment was determined using micro-finite element analysis (micro-FE) on baseline images. The formation of mineralised tissue showed strong association with areas of higher mechanical strain (area-under-the-curve: 0.91 ± 0.04, true positive rate: 0.85 ± 0.05) while surface based strains could correctly classify 43% of remodelling events. These findings support our hypotheses by showing a direct association between the local mechanical strains and the formation of mineralised tissue.

Fracture healing in mice lacking Pten in osteoblasts: a micro-computed tomography image-based analysis of the mechanical properties of the femur

In the United States, approximately eight million osseous fractures are reported annually, of which 5-10% fail to create a bony union. Osteoblast-specific deletion of the gene Pten in mice has been found to stimulate bone growth and accelerate fracture healing. Healing rates at four weeks increased in femurs from Pten osteoblast conditional knock-out mice (Pten-CKO) compared to wild-type mice (WT) of the same genetic strain as measured by an increase in mechanical stiffness and failure load in four-point bending tests. Preceding mechanical testing, each femur was imaged using a Skyscan 1172 micro-computed tomography (μCT) scanner (Skyscan, Kontich, Belgium). The present study used µCT image-based analysis to test the hypothesis that the increased femoral fracture force and stiffness in Pten-CKO were due to greater section properties with the same effective material properties as that of the WT. The second moment of area and section modulus were computed in ImageJ 1.46 (National Institutes of Health) and used to predict the effective flexural modulus and the stress at failure for fourteen pairs of intact and callus WT and twelve pairs of intact and callus Pten-CKO femurs. For callus and intact femurs, the failure stress and tissue mineral density of the Pten-CKO and WT were not different; however, the section properties of the Pten-CKO were more than twice as large 28 days post-fracture. It was therefore concluded, when the gene Pten was conditionally knocked-out in osteoblasts, the resulting increased bending stiffness and force to fracture were due to increased section properties.

Finite element analysis of a bone healing model: 1-year follow-up after internal fixation surgery for femoral fracture

OBJECTIVE

Finite element analysis was used to compare preoperative and postoperative stress distribution of a bone healing model of femur fracture, to identify whether broken ends of fractured bone would break or not after fixation dislodgement one year after intramedullary nailing. Method s: Using fast, personalized imaging, bone healing models of femur fracture were constructed based on data from multi-slice spiral computed tomography using Mimics, Geomagic Studio, and Abaqus software packages. The intramedullary pin was removed by Boolean operations before fixation was dislodged. Loads were applied on each model to simulate a person standing on one leg. The von Mises stress distribution, maximum stress, and its location was observed. Results : According to 10 kinds of display groups based on material assignment, the nodes of maximum and minimum von Mises stress were the same before and after dislodgement, and all nodes of maximum von Mises stress were outside the fracture line. The maximum von Mises stress node was situated at the bottom quarter of the femur. The von Mises stress distribution was identical before and after surgery. Conclusion : Fast, personalized model establishment can simulate fixation dislodgement before operation, and personalized finite element analysis was performed to successfully predict whether nail dislodgement would disrupt femur fracture or not.