A multiscale model to predict current absolute risk of femoral fracture in a postmenopausal population

Osteoporotic hip fractures are a major healthcare problem. Fall severity and bone strength are important risk factors of hip fracture. This study aims to obtain a mechanistic explanation for fracture risk in dependence of these risk factors. A novel modelling approach is developed that combines models at different scales to overcome the challenge of a large space–time domain of interest and considers the variability of impact forces between potential falls in a subject. The multiscale model and its component models are verified with respect to numerical approximations made therein, the propagation of measurement uncertainties of model inputs is quantified, and model predictions are validated against experimental and clinical data. The main results are model predicted absolute risk of current fracture (ARF0) that ranged from 1.93 to 81.6% (median 36.1%) for subjects in a retrospective cohort of 98 postmenopausal British women (49 fracture cases and 49 controls); ARF0 was computed up to a precision of 1.92 percentage points (pp) due to numerical approximations made in the model; ARF0 possessed an uncertainty of 4.00 pp due to uncertainties in measuring model inputs; ARF0 classified observed fracture status in the above cohort with AUC = 0.852 (95% CI 0.753–0.918), 77.6% specificity (95% CI 63.4–86.5%) and 81.6% sensitivity (95% CI 68.3–91.1%). These results demonstrate that ARF0 can be computed using the model with sufficient precision to distinguish between subjects and that the novel mechanism of fracture risk determination based on fall dynamics, hip impact and bone strength can be considered validated

Effect of trabecular bone loss on cortical strain rate during impact in an in vitro model of avian femur

BACKGROUND: Osteoporotic hip fractures occur due to loss of cortical and trabecular bone mass and consequent degradation in whole bone strength. The direct cause of most fractures is a fall, and hence, characterizing the mechanical behavior of a whole osteopenic bone under impact is important. However, very little is known about the mechanical interactions between cortical and trabecular bone during impact, and it is specifically unclear to what extent epiphyseal trabecular bone contributes to impact resistance of whole bones. We hypothesized that trabecular bone serves as a structural support to the cortex during impact, and hence, loss of a critical mass of trabecular bone reduces internal constraining of the cortex, and, thereby, decreases the impact tolerance of the whole bone. METHODS: To test this hypothesis, we conducted cortical strain rate measurements in adult chicken's proximal femora subjected to a Charpy impact test, after removing different trabecular bone core masses to simulate different osteopenic severities. RESULTS: We found that removal of core trabecular bone decreased by ~10-fold the cortical strain rate at the side opposite to impact (p < 0.01), i.e. from 359,815 ± 1799 μm/m per second (mean ± standard error) for an intact (control) specimen down to 35,997 ± 180 μm/m per second where 67% of the total trabecular bone mass (~0.7 grams in adult chicken) were removed. After normalizing the strain rate by the initial weight of bone specimens, a sigmoid relation emerged between normalized strain rate and removed mass of trabecular bone, showing very little effect on the cortex strain rate if below 10% of the trabecular mass is removed, but most of the effect was already apparent for less than 30% trabecular bone loss. An analytical model of the experiments supported this behavior. CONCLUSION: We conclude that in our in vitro avian model, loss of over 10% of core trabecular bone substantially altered the deformation response of whole bone to impact, which supports the above hypothesis and indicates that integrity of trabecular bone is critical for resisting impact loads

The application of digital volume correlation (DVC) to evaluate strain predictions generated by finite element models of the osteoarthritic humeral head

Continuum-level finite element models (FEMs) of the humerus offer the ability to evaluate joint replacement designs preclinically; however, experimental validation of these models is critical to ensure accuracy. The objective of the current study was to quantify experimental full-field strain magnitudes within osteoarthritic (OA) humeral heads by combining mechanical loading with volumetric microCT imaging and digital volume correlation (DVC). The experimental data was used to evaluate the accuracy of corresponding FEMs. Six OA humeral head osteotomies were harvested from patients being treated with total shoulder arthroplasty and mechanical testing was performed within a microCT scanner. MicroCT images (33.5 µm isotropic voxels) were obtained in a pre- and post-loaded state and BoneDVC was used to quantify full-field experimental strains (≈ 1 mm nodal spacing, accuracy = 351 µstrain, precision = 518 µstrain). Continuum-level FEMs with two types of boundary conditions (BCs) were simulated: DVC-driven and force-driven. Accuracy of the FEMs was found to be sensitive to the BC simulated with better agreement found with the use of DVC-driven BCs (slope = 0.83, r2 = 0.80) compared to force-driven BCs (slope = 0.22, r2 = 0.12). This study quantified mechanical strain distributions within OA trabecular bone and demonstrated the importance of BCs to ensure the accuracy of predictions generated by corresponding FEMs

Are CT-Based Finite Element Model Predictions of Femoral Bone Strengthening Clinically Useful?

Purpose of Review: This study reviews the available literature to compare the accuracy of areal bone mineral density derived from dual X-ray absorptiometry (DXA-aBMD) and of subject-specific finite element models derived from quantitative computed tomography (QCT-SSFE) in predicting bone strength measured experimentally on cadaver bones, as well as their clinical accuracy both in terms of discrimination and prediction. Based on this information, some basic cost-effectiveness calculations are performed to explore the use of QCT-SSFE instead of DXA-aBMD in (a) clinical studies with femoral strength as endpoint, (b) predictor of the risk of hip fracture in low bone mass patients. Recent Findings: Recent improvements involving the use of smooth-boundary meshes, better anatomical referencing for proximal-only scans, multiple side-fall directions, and refined boundary conditions increase the predictive accuracy of QCT-SSFE. Summary: If these improvements are adopted, QCT-SSFE is always preferable over DXA-aBMD in clinical studies with femoral strength as the endpoint, while it is not yet cost-effective as a hip fracture risk predictor, although pathways that combine both QCT-SSFE and DXA-aBMD are promising

Patient-specific finite element estimated femur strength as a predictor of the risk of hip fracture: the effect of methodological determinants

Summary: A finite element modelling pipeline was adopted to predict femur strength in a retrospective cohort of 100 women. The effects of the imaging protocol and the meshing technique on the ability of the femur strength to classify the fracture and the control groups were analysed. Introduction: The clinical standard to estimate the risk of osteoporotic hip fracture is based on the areal bone mineral density (aBMD). A few retrospective studies have concluded that finite element (FE)-based femoral strength is a better classifier of fracture and control groups than the aBMD, while others could not find significant differences. We investigated the effect of the imaging protocol and of the FE modelling techniques on the discriminatory power of femoral strength. Methods: A retrospective cohort of 100 post-menopausal women (50 with hip fracture, 50 controls) was examined. Each subject received a dual-energy absorptiometry (DXA) exam and a computed tomography (CT) scan of the proximal femur region. Each case was modelled a number of times, using different modelling pipelines, and the results were compared in terms of accuracy in discriminating the fracture and the control cases. The baseline pipeline involved local anatomical orientation and mesh morphing. Revised pipelines involved global anatomical orientation using a full-femur atlas registration and an optimised meshing algorithm. Minimum physiological (MPhyS) and pathological (MPatS) strengths were estimated for each subject. Area under the receiver operating characteristic (ROC) curve (AUC) was calculated to compare the ability of MPhyS, MPatS and aBMD to classify the control and the cases. Results: Differences in the modelling protocol were found to considerably affect the accuracy of the FE predictors. For the most optimised protocol, logistic regression showed aBMD Neck , MPhyS and MPatS to be significantly associated with the facture status, with AUC of 0.75, 0.75 and 0.79, respectively. Conclusion: The study emphasized the necessity of modelling the whole femur anatomy to develop a robust FE-based tool for hip fracture risk assessment. FE-strength performed only slightly better than the aBMD in discriminating the fracture and control cases. Differences between the published studies can be explained in terms of differences in the modelling protocol and cohort design

Analysis of miRNA and mRNA Expression Profiles Highlights Alterations in Ionizing Radiation Response of Human Lymphocytes under Modeled Microgravity

BACKGROUND: Ionizing radiation (IR) can be extremely harmful for human cells since an improper DNA-damage response (DDR) to IR can contribute to carcinogenesis initiation. Perturbations in DDR pathway can originate from alteration in the functionality of the microRNA-mediated gene regulation, being microRNAs (miRNAs) small noncoding RNA that act as post-transcriptional regulators of gene expression. In this study we gained insight into the role of miRNAs in the regulation of DDR to IR under microgravity, a condition of weightlessness experienced by astronauts during space missions, which could have a synergistic action on cells, increasing the risk of radiation exposure. METHODOLOGY/PRINCIPAL FINDINGS: We analyzed miRNA expression profile of human peripheral blood lymphocytes (PBL) incubated for 4 and 24 h in normal gravity (1 g) and in modeled microgravity (MMG) during the repair time after irradiation with 0.2 and 2Gy of \u3b3-rays. Our results show that MMG alters miRNA expression signature of irradiated PBL by decreasing the number of radio-responsive miRNAs. Moreover, let-7i*, miR-7, miR-7-1*, miR-27a, miR-144, miR-200a, miR-598, miR-650 are deregulated by the combined action of radiation and MMG. Integrated analyses of miRNA and mRNA expression profiles, carried out on PBL of the same donors, identified significant miRNA-mRNA anti-correlations of DDR pathway. Gene Ontology analysis reports that the biological category of "Response to DNA damage" is enriched when PBL are incubated in 1 g but not in MMG. Moreover, some anti-correlated genes of p53-pathway show a different expression level between 1 g and MMG. Functional validation assays using luciferase reporter constructs confirmed miRNA-mRNA interactions derived from target prediction analyses. CONCLUSIONS/SIGNIFICANCE: On the whole, by integrating the transcriptome and microRNome, we provide evidence that modeled microgravity can affects the DNA-damage response to IR in human PBL