Fast Prediction of Femoral Biomechanics Using Supervised Machine Learning and Statistical Shape Modeling

Finite element (FE) analysis is an important computational tool in biomechanics. However, its adoption into clinical practice has been hampered by its computational complexity and required high technical competences for clinicians. In this paper we propose a supervised learning approach to predict the outcome of the FE analysis. We demonstrate our approach on clinical CT and X-ray femur images for FE predictions ( FEP), with features extracted, respectively, from a statistical shape model and from 2D-based morphometric and density information. Using leave-one-out experiments and sensitivity analysis, comprising a database of 89 clinical cases, our method is capable of predicting the distribution of stress values for a walking loading condition with an average correlation coefficient of 0.984 and 0.976, for CT and X-ray images, respectively. These findings suggest that supervised learning approaches have the potential to leverage the clinical integration of mechanical simulations for the treatment of musculoskeletal conditions.

Femoral Biomechanic and microtomography from male rats submitted to dietary restricition supplemented with sucrose

Food restriction reduces body weight and influence bone mass and also is correlated with bone mineral density (BMD). Mechanisms have been proposed for the loss of BMD after body weight reduction, including reduced energy intake. Growing 8 wk-old Wistar male rats were randomly divided into Control and Calorie restriction associated with sucrose 30% (CRS). These animals were subjected to intermittent food restriction during 8 weeks and had free access to tap water and sucrose30% in distilled water. The rats were euthanized at the end of week 8, blood collected from abdominal aorta artery, femurs cleaned of adherent soft tissues, scanned using dual energy X-ray absorptiometry, structural and material properties determined by three-point bending testing in the mid-diaphyseal region, bone surface tested in a microhardness tester and microstructure was assessed in a microcomputer tomography. In CRS animals body weight decreased significantly relative to the Control animals. There was a clear option for high-sucrose beverage in CRS animals. No difference was observed in biochemical, densitometric and biomechanical analyzes. Results from micro CT showed only significant difference in connectivity of trabecular bone. It has been suggested that rats submitted to food restriction consumed sugar not because of its inherent palatability, but in order to alter their macronutrient balance and animals need to meet energy demands in high-sucrose.

Sensitivity of proximal femoral stiffness and areal bone mineral density to changes in bone geometry and density

Areal bone mineral density (aBMD) is the most common surrogate measurement for assessing the bone strength of the proximal femur associated with osteoporosis. Additional factors, however, contribute to the overall strength of the proximal femur, primarily the anatomical geometry. Finite element analysis (FEA) is an effective and widely used computerbased simulation technique for modeling mechanical loading of various engineering structures, providing predictions of displacement and induced stress distribution due to the applied load. FEA is therefore inherently dependent upon both density and anatomical geometry. FEA may be performed on both three-dimensional and two-dimensional models of the proximal femur derived from radiographic images, from which the mechanical stiffness may be redicted. It is examined whether the outcome measures of two-dimensional FEA, two-dimensional, finite element analysis of X-ray images (FEXI), and three-dimensional FEA computed stiffness of the proximal femur were more sensitive than aBMD to changes in trabecular bone density and femur geometry. It is assumed that if an outcome measure follows known trends with changes in density and geometric parameters, then an increased sensitivity will be indicative of an improved prediction of bone strength. All three outcome measures increased non-linearly with trabecular bone density, increased linearly with cortical shell thickness and neck width, decreased linearly with neck length, and were relatively insensitive to neck-shaft angle. For femoral head radius, aBMD was relatively insensitive, with two-dimensional FEXI and threedimensional FEA demonstrating a non-linear increase and decrease in sensitivity, respectively. For neck anteversion, aBMD decreased non-linearly, whereas both two-dimensional FEXI and three dimensional FEA demonstrated a parabolic-type relationship, with maximum stiffness achieved at an angle of approximately 15o. Multi-parameter analysis showed that all three outcome measures demonstrated their highest sensitivity to a change in cortical thickness. When changes in all input parameters were considered simultaneously, three and twodimensional FEA had statistically equal sensitivities (0.41±0.20 and 0.42±0.16 respectively, p = ns) that were significantly higher than the sensitivity of aBMD (0.24±0.07; p = 0.014 and 0.002 for three-dimensional and two-dimensional FEA respectively). This simulation study suggests that since mechanical integrity and FEA are inherently dependent upon anatomical geometry, FEXI stiffness, being derived from conventional two-dimensional radiographic images, may provide an improvement in the prediction of bone strength of the proximal femur than currently provided by aBMD.

Shortening cemented femoral implants. An in vitro investigation to quantify exeter femoral implant rotational stability vs simulated implant length

The Exeter stems vary in length from 90 to 150 mm. The shorter stems generally have lower offsets. The purpose of this study was to determine if length of stem, with fixed offset, affected rotational stability. Mechanical testing was carried out on 10 implant-cement constructs with 2 loading profiles, rising from chair and stair climbing, at different simulated implant lengths using purpose-built apparatus. This paper presents a mechanism for clinically observed rotational stability and explains the mechanical characteristics required for rotational stability in Exeter femoral stems. © 2012.

The mechanical effect of the existing cement mantle on the in-cement femoral revision

Cement-in-cement revision hip arthroplasty is an increasingly popular technique to replace a loose femoral stem which retains much of the original cement mantle. However, some concern exists regarding the retention of the existing fatigued and aged cement in such cement-in-cement revisions. This study investigates whether leaving an existing fatigued and aged cement mantle degrades the mechanical performance of a cement-in-cement revision construct.

Three-dimensional primary stability of cementless femoral stems

OBJECTIVE: This study investigates by means of a new bone-prosthesis interface motion detector whether conceptual design differences of femoral stems are reflected in their primary stability pattern. DESIGN: An in vitro experiment using a biaxial materials testing machine in combination with three-dimensional motion measurement devices was performed. BACKGROUND: Primary stability of uncemented total hip replacements is considered to be a prerequisite for the quality of bony ongrowth to the femoral stem. Dynamic motion as a response to loading as well as total motion of the prosthesis have to be considered under quasi-physiological cyclic loading conditions. METHODS: Seven paired fresh cadaveric femora were used for the testing of two types of uncemented femoral stems with different anchoring concepts: CLS stem (Spotorno) and Cone Prosthesis (Wagner). Under sinusoidal cyclic loading mimicking in vivo hip joint forces a new measurement technique was applied allowing for the analysis of the three-dimensional interface motion. RESULTS: Considerable differences between the two prostheses could be detected both in their dynamic motion and total motion behaviour. Whereas the CLS stem, due to the wedge-shaped concept, provides smaller total motions, the longitudinal ribs of the Cone prostheses result in a substantially smaller dynamic motion. CONCLUSIONS: The measuring technique provided reliable and accurate data illustrating the three-dimensional interface motion of uncemented femoral stems.

Measurement of structural anisotropy in femoral trabecular bone using clinical-resolution CT images

Discrepancies in finite-element model predictions of bone strength may be attributed to the simplified modeling of bone as an isotropic structure due to the resolution limitations of clinical-level Computed Tomography (CT) data. The aim of this study is to calculate the preferential orientations of bone (the principal directions) and the extent to which bone is deposited more in one direction compared to another (degree of anisotropy). Using 100 femoral trabecular samples, the principal directions and degree of anisotropy were calculated with a Gradient Structure Tensor (GST) and a Sobel Structure Tensor (SST) using clinical-level CT. The results were compared against those calculated with the gold standard Mean-Intercept-Length (MIL) fabric tensor using micro-CT. There was no significant difference between the GST and SST in the calculation of the main principal direction (median error=28°), and the error was inversely correlated to the degree of transverse isotropy (r=−0.34, p<0.01). The degree of anisotropy measured using the structure tensors was weakly correlated with the MIL-based measurements (r=0.2, p<0.001). Combining the principal directions with the degree of anisotropy resulted in a significant increase in the correlation of the tensor distributions (r=0.79, p<0.001). Both structure tensors were robust against simulated noise, kernel sizes, and bone volume fraction. We recommend the use of the GST because of its computational efficiency and ease of implementation. This methodology has the promise to predict the structural anisotropy of bone in areas with a high degree of anisotropy, and may improve the in vivo characterization of bone.

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.

Rebound-associated vertebral fractures after discontinuation of denosumab-from clinic and biomechanics.

Rebound-associated vertebral fractures may follow treatment discontinuation of highly potent reversible bone antiresorptives, resulting from the synergy of rapid bone resorption and accelerated microdamage accumulation in trabecular bone. INTRODUCTION The purposes of this study are to characterize rebound-associated vertebral fractures following the discontinuation of a highly potent reversible antiresorptive therapy based on clinical observation and propose a pathophysiological rationale. METHODS This study is a case report of multiple vertebral fractures early after discontinuation of denosumab therapy in a patient with hormone receptor-positive non-metastatic breast cancer treated with an aromatase inhibitor. RESULTS Discontinuation of highly potent reversible bone antiresorptives such as denosumab may expose patients to an increased fracture risk due to the joined effects of absent microdamage repair during therapy followed by synchronous excess activation of multiple bone remodelling units at the time of loss-of-effect. We suggest the term rebound-associated vertebral fractures (RVF) for this phenomenon characterized by the presence of multiple new clinical vertebral fractures, associated with either no or low trauma, in a context consistent with the presence of high bone turnover and rapid loss of lumbar spine bone mineral density (BMD) occurring within 3 to 12 months after discontinuation (loss-of-effect) of a reversible antiresorptive therapy in the absence of secondary causes of bone loss or fractures. Unlike atypical femoral fractures that emerge from failure of microdamage repair in cortical bone with long-term antiresorptive treatment, RVF originate from the synergy of rapid bone resorption and accelerated microdamage accumulation in trabecular bone triggered by the discontinuation of highly potent reversible antiresorptives. CONCLUSIONS Studies are urgently needed to i) prove the underlying pathophysiological processes suggested above, ii) establish the predictive criteria exposing patients to an increased risk of RVF, and iii) determine appropriate treatment regimens to be applied in such patients.