Zwahlen, Alexander

A bone loading estimation algorithm was previously developed that provides in vivo loading conditions required for in vivo bone remodelling simulations. The algorithm derives a bone's loading history from its microstructure as assessed by high-resolution (HR) computed tomography (CT). This reverse engineering approach showed accurate and realistic results based on micro-CT and HR-peripheral quantitative CT images. However, its voxel size dependency, reproducibility and sensitivity still need to be investigated, which is the purpose of this study. Voxel size dependency was tested on cadaveric distal radii with micro-CT images scanned at 25 µm and downscaled to 50, 61, 75, 82, 100, 125 and 150 µm. Reproducibility was calculated with repeated in vitro as well as in vivo HR-pQCT measurements at 82 µm. Sensitivity was defined using HR-pQCT images from women with fracture versus non-fracture, and low versus high bone volume fraction, expecting similar and different loading histories, respectively. Our results indicate that the algorithm is voxel size independent within an average (maximum) error of 8.2% (32.9%) at 61 µm, but that the dependency increases considerably at voxel sizes bigger than 82 µm. In vitro and in vivo reproducibility are up to 4.5% and 10.2%, respectively, which is comparable to other in vitro studies and slightly higher than in other in vivo studies. Subjects with different bone volume fraction were clearly distinguished but not subjects with and without fracture. This is in agreement with bone adapting to customary loading but not to fall loads. We conclude that the in vivo bone loading estimation algorithm provides reproducible, sensitive and fairly voxel size independent results at up to 82 µm, but that smaller voxel sizes would be advantageous.

ABSTRACT More accurate techniques to estimate fracture risk could help reduce the burden of fractures in postmenopausal women. Although micro-finite element (µFE) simulations allow a direct assessment of bone mechanical performance, in this first clinical study we investigated whether the additional information obtained using geometrically and materially nonlinear µFE simulations allows a better discrimination between fracture cases and controls. We used patient data and high-resolution peripheral quantitative computed tomography (HRpQCT) measurements from our previous clinical study on fracture risk, which compared 100 postmenopausal women with a distal forearm fracture to 105 controls. Analyzing these data with the nonlinear µFE simulations, the odds ratio (OR) for the factor-of-risk (yield load divided by the expected fall load) was marginally higher (1.99; 95% confidence interval [CI], 1.41–2.77) than for the factor-of-risk computed from linear µFE (1.89; 95% CI, 1.37–2.69). The yield load and the energy absorbed up to the yield point as computed from nonlinear µFE were highly correlated with the initial stiffness (R2 = 0.97 and 0.94, respectively) and could therefore be derived from linear simulations with little loss in precision. However, yield deformation was not related to any other measurement performed and was itself a good predictor of fracture risk (OR, 1.89; 95% CI, 1.39–2.63). Moreover, a combined risk score integrating information on relative bone strength (yield load-based factor-of-risk), bone ductility (yield deformation), and the structural integrity of the bone under critical loads (cortical plastic volume) improved the separation of cases and controls by one-third (OR, 2.66; 95% CI, 1.84–4.02). We therefore conclude that nonlinear µFE simulations provide important additional information on the risk of distal forearm fractures not accessible from linear µFE nor from other techniques assessing bone microstructure, density, or mass. © 2013 American Society for Bone and Mineral Research.