Mechanical Loading of the Femoral Neck in Human Locomotion

The effect of regular running on the bone tissue of middle-aged men and women

BACKGROUND

Many authors consider running to be a protective physical activity (PA) in bone health. However, many studies also show inconsistencies in their results. The objective of the study is to analyze the effect of cumulative loading rate (TCL) on the bone mass of middle-aged runners and non-runners is assessed.

METHODS

This cross-sectional study included 322 individuals. There were 212 runners (109 male, 103 female) and those were individuals who did >10 km of running per week. There were 110 non-runners (54 male, 56 female). This group included individuals who did not adhere to the WHO (2020) recommendations for PA. The average age in the individual groups ranged from 40.9±4.1 to 42.3±4.8 years. Bone parameters were measured on the lower extremities and vertebral spine using the DXA method (Hologic QDR Horizon A). Multi-regression dependencies analysis was used to assess the results.

RESULTS

The results of the multi-regression dependencies analysis showed that the bone mineral content (BMC) and bone mineral density (BMD) are significantly influenced by the TCL and gender.

CONCLUSIONS

Therefore, we can conclude that running could be a suitable PA for preventing the reduction of BMD in the middle-aged population, especially in the lower limbs.

Influence of a change in activity regime on femoral bone architecture and failure behaviour

The incidence and morbidity of femoral fractures increases drastically with age. Femoral architecture and associated fracture risk are strongly influenced by loading during physical activities and it has been shown that the rate of loss of bone mineral density is significantly lower for active individuals than inactive. The objective of this work is to evaluate the impact of a cessation of some physical activities on elderly femoral structure and fracture behaviour. The authors previously established a biofidelic finite element model of the femur considered as a structure optimised to loading associated with daily activities. The same structural optimisation algorithm was used here to quantify the changes in bone architecture following cessation of stair climbing and sit-to-stand. Side fall fracture simulations were run on the adapted bone structures using a damage elasticity formulation. Total cortical and trabecular bone volume and failure load reduced in all cases of activity cessation. Bone loss distribution was strongly heterogeneous, with some locations even showing increased bone volume. This work suggests that maintaining the physical activities involved in the daily routine of a young healthy adult would help reduce the risk of femoral fracture in the elderly population by preventing bone loss.

Macroscopic differences in adult human femora are linked to body mass index

Bone functional adaptation is routinely invoked to interpret skeletal morphology despite ongoing debate regarding the limits of the bone response to mechanical stimuli. The wide variation in human body mass presents an opportunity to explore the relationship between mechanical load and skeletal response in weight-bearing elements. Here, we examine variation in femoral macroscopic morphology as a function of body mass index (BMI), which is used as a metric of load history. A sample of 80 femora (40 female; 40 male) from recent modern humans was selected from the Texas State University Donated Skeletal Collection. Femora were imaged using x-ray computed tomography (voxel size ~0.5 mm), and segmented to produce surface models. Landmark-based geometric morphometric analyses based on the Coherent Point Drift algorithm were conducted to quantify shape. Principal components analyses were used to summarize shape variation, and component scores were regressed on BMI. Within the male sample, increased BMI was associated with a mediolaterally expanded femoral shaft, as well as increased neck-shaft angle and decreased femoral neck anteversion angle. No statistically significant relationships between shape and BMI were found in the female sample. While mechanical stimulus has traditionally been applied to changes in long bong diaphyseal shape it appears that bone functional adaptation may also result in fundamental changes in the shape of skeletal elements.

Effect of Muscle Forces on Femur During Level Walking Using a Virtual Population of Older Women

Aging is associated with a greater risk of muscle and bone disorders such as sarcopenia and osteoporosis. These conditions substantially affect one's mobility and quality of life. In the past, muscles and bones are often studied separately using generic or scaled information that are not personal-specific, nor are they representative of the large variations seen in the elderly population. Consequently, the mechanical interaction between the aged muscle and bone is not well understood, especially when carrying out daily activities. This study presents a coupling approach across the body and the organ level, using fully personal-specific musculoskeletal and finite element models in order to study femoral loading during level walking. Variations in lower limb muscle volume/force were examined using a virtual population. These muscle forces were then applied to the finite element model of the femur to study the variations in predicted strains. The study shows that effective coupling across two scales can be carried out to study the muscle-bone interaction in elderly women. The generation of a virtual population is a feasible approach to augment anatomical variations based on a small population that could mimic variations seen in a larger cohort. This is a valuable alternative to overcome the limitation or the need to collect dataset from a large population, which is both time and resource consuming.

Using tissue clearing and light sheet fluorescence microscopy for the three-dimensional analysis of sensory and sympathetic nerve endings that innervate bone and dental tissue of mice

Bone and dental tissues are richly innervated by sensory and sympathetic neurons. However, the characterization of the morphology, molecular phenotype, and distribution of nerves that innervate hard tissue has so far mostly been limited to thin histological sections. This approach does not adequately capture dispersed neuronal projections due to the loss of important structural information during three-dimensional (3D) reconstruction. In this study, we modified the immunolabeling-enabled imaging of solvent-cleared organs (iDISCO/iDISCO+) clearing protocol to image high-resolution neuronal structures in whole femurs and mandibles collected from perfused C57Bl/6 mice. Axons and their nerve terminal endings were immunolabeled with antibodies directed against protein gene product 9.5 (pan-neuronal marker), calcitonin gene-related peptide (peptidergic nociceptor marker), or tyrosine hydroxylase (sympathetic neuron marker). Volume imaging was performed using light sheet fluorescence microscopy. We report high-quality immunolabeling of the axons and nerve terminal endings for both sensory and sympathetic neurons that innervate the mouse femur and mandible. Importantly, we are able to follow their projections through full 3D volumes, highlight how extensive their distribution is, and show regional differences in innervation patterns for different parts of each bone (and surrounding tissues). Mapping the distribution of sensory and sympathetic axons, and their nerve terminal endings, in different bony compartments may be important in further elucidating their roles in health and disease.

Composition and micromechanical properties of the femoral neck compact bone in relation to patient age, sex and hip fracture occurrence

Current clinical methods of bone health assessment depend to a great extent on bone mineral density (BMD) measurements. However, these methods only act as a proxy for bone strength and are often only carried out after the fracture occurs. Besides BMD, composition and tissue-level mechanical properties are expected to affect the whole bone's strength and toughness. While the elastic properties of the bone extracellular matrix (ECM) have been extensively investigated over the past two decades, there is still limited knowledge of the yield properties and their relationship to composition and architecture. In the present study, morphological, compositional and micropillar compression bone data was collected from patients who underwent hip arthroplasty. Femoral neck samples from 42 patients were collected together with anonymous clinical information about age, sex and primary diagnosis (coxarthrosis or hip fracture). The femoral neck cortex from the inferomedial region was analyzed in a site-matched manner using a combination of micromechanical testing (nanoindentation, micropillar compression) together with micro-CT and quantitative polarized Raman spectroscopy for both morphological and compositional characterization. Mechanical properties, as well as the sample-level mineral density, were constant over age. Only compositional properties demonstrate weak dependence on patient age: decreasing mineral to matrix ratio (p = 0.02, R2 = 0.13, 2.6 % per decade) and increasing amide I sub-peak ratio I∼1660/I∼1683 (p = 0.04, R2 = 0.11, 1.5 % per decade). The patient's sex and diagnosis did not seem to influence investigated bone properties. A clear zonal dependence between interstitial and osteonal cortical zones was observed for compositional and elastic bone properties (p < 0.0001). Site-matched microscale analysis confirmed that all investigated mechanical properties except yield strain demonstrate a positive correlation with the mineral fraction of bone. The output database is the first to integrate the experimentally assessed microscale yield properties, local tissue composition and morphology with the available patient clinical information. The final dataset was used for bone fracture risk prediction in-silico through the principal component analysis and the Naïve Bayes classification algorithm. The analysis showed that the mineral to matrix ratio, indentation hardness and micropillar yield stress are the most relevant parameters for bone fracture risk prediction at 70 % model accuracy (0.71 AUC). Due to the low number of samples, further studies to build a universal fracture prediction algorithm are anticipated with the higher number of patients (N > 200). The proposed classification algorithm together with the output dataset of bone tissue properties can be used for the future comparison of existing methods to evaluate bone quality as well as to form a better understanding of the mechanisms through which bone tissue is affected by aging or disease.

Validity of evaluating spinal kinetics without participant-specific kinematics

Musculoskeletal models are commonly used to estimate in vivo spinal loads under various loading conditions. Typically, participant-specific measured kinematics (PSMK) are coupled with participant-specific models, but obtaining PSMK data can be costly and infeasible in large studies or clinical practice. Thus, we evaluated two alternative methods to estimate spinal loads without PSMK: 1) ensemble average kinematics (EAK) based on kinematics from all participants; and 2) using separately measured individual kinematics (SMIK) from multiple other participants as inputs, then averaging the resulting loads. This study compares the dynamic spine loading patterns and peak loads in older adults performing five lifting tasks using PSMK, EAK and SMIK. Median root mean square errors of EAK and SMIK methods versus PSMK ranged from 18 to 72% body weight for compressive loads and from 2 to 25% body weight for shear loads, with median cross-correlations ranging from 0.931 to 0.991. The root mean square errors and cross-correlations between repeated PSMK trials fell within similar ranges. Compressive peak loads evaluated by EAK and SMIK were not different than PSMK in 12 of 15 cases, while by comparison repeated PSMK trials were not different in 13 of 15 cases. Overall, the resulting spine loading magnitudes and profiles using EAK or SMIK were not notably different than using a PSMK approach, and differences were not greater than between two PSMK trials. Thus, these findings indicate that these approaches may be used to make reasonable estimates of dynamic spinal loading without direct measurement of participant kinematics.

Effect of gait types and external weight carrying strategies on the femoral neck strains during stair descent

Gait and weight carrying method may change the femoral neck load during stair descent. Applying specific gait and weight carrying methods may reduce the femoral neck load during stair descent, which may reduce hip pain, hip pain related falls and fall related fractures for the older population. The purpose of this study was to test the effect of different gait types (step-over-step v.s. step-by-step) and external weight carrying strategies (ipsilateral v.s. contralateral side) on the femoral neck load, discover which method could reduce the femoral neck load effectively. Seventeen healthy adults from 50 to 70 yrs old were recruited. The kinematic and kinetic analysis, musculoskeletal modelling method were used to estimate the joint and muscle loads for the lower extremities. Finite element analysis was used with the femur model to calculate the femoral neck strains during stair descent with different gait types and weight carrying strategies. The compressive strains were reduced for step-by-step gait method than step-over-step (p<0.015, 12.3-17.4% decrease of strains), the tensile strains were significantly increased for the trailing leg of step-by-step than the leading leg (p<0.001, 24.7% increase of strains). Contralateral weight carrying increased compressive and tensile strains than ipsilateral (p<0.001, 9.9-24.5% increase of strains) in most conditions. Applying step-by-step method and avoiding contralateral side weight carrying could be effective to reduce femoral neck strains. These outcomes could be helpful for the older population to reduce the risks of hip pain, femoral neck pain or pain related falls and fractures.

Elastic Modulus and Its Relation to Apparent Mineral Density in Juvenile Equine Bones of the Lower Limb

Density-modulus relationships are necessary to develop finite element models of bones that may be used to evaluate local tissue response to different physical activities. It is unknown if juvenile equine trabecular bone may be described by the same density-modulus as adult equine bone, and how the density-modulus relationship varies with anatomical location and loading direction. To answer these questions, trabecular bone cores from the third metacarpal (MC3) and proximal phalanx (P1) bones of juvenile horses (age <1 yr) were machined in the longitudinal (n = 134) and transverse (n = 90) directions and mechanically tested in compression. Elastic modulus was related to apparent computed tomography density of each sample using power law regressions. We found that density-modulus relationships for juvenile equine trabecular bone were significantly different for each anatomical location (MC3 versus P1) and orientation (longitudinal versus transverse). Use of the incorrect density-modulus relationship resulted in increased root mean squared percent error of the modulus prediction by 8-17%. When our juvenile density-modulus relationship was compared to one of an equivalent location in adult horses, the adult relationship resulted in an approximately 80% increase in error of the modulus prediction. Moving forward, more accurate models of young bone can be developed and used to evaluate potential exercise regimens designed to encourage bone adaptation.

Exploration of the synergistic role of cortical thickness asymmetry ("Trabecular Eccentricity" concept) in reducing fracture risk in the human femoral neck and a control bone (Artiodactyl Calcaneus)

The mechanobiology of the human femoral neck is a focus of research for many reasons including studies that aim to curb age-related bone loss that contributes to a near-exponential rate of hip fractures. Many believe that the femoral neck is often loaded in rather simple bending, which causes net tension stress in the upper (superior) femoral neck and net compression stress in its inferior aspect ("T/C paradigm"). This T/C loading regime lacks in vivo proof. The "C/C paradigm" is a plausible alternative simplified load history that is characterized by a gradient of net compression across the entire femoral neck; action of the gluteus medius and external rotators of the hip are important in this context. It is unclear which paradigm is at play in natural loading due to lack of in vivo bone strain data and deficiencies in understanding mechanisms and manifestations of bone adaptation in tension vs. compression. For these reasons, studies of the femoral neck would benefit from being compared to a 'control bone' that has been proven, by strain data, to be habitually loaded in bending. The artiodactyl (sheep and deer) calcaneus model has been shown to be a very suitable control in this context. However, the application of this control in understanding the load history of the femoral neck has only been attempted in two prior studies, which did not examine the interplay between cortical and trabecular bone, or potential load-sharing influences of tendons and ligaments. Our first goal is to compare fracture risk factors of the femoral neck in both paradigms. Our second goal is to compare and contrast the deer calcaneus to the human femoral neck in terms of fracture risk factors in the T/C paradigm (the C/C paradigm is not applicable in the artiodactyl calcaneus due to its highly constrained loading). Our third goal explores interplay between dorsal/compression and plantar/tension regions of the deer calcaneus and the load-sharing roles of a nearby ligament and tendon, with insights for translation to the femoral neck. These goals were achieved by employing the analytical model of Fox and Keaveny (J. Theoretical Biology 2001, 2003) that estimates fracture risk factors of the femoral neck. This model focuses on biomechanical advantages of the asymmetric distribution of cortical bone in the direction of habitual loading. The cortical thickness asymmetry of the femoral neck (thin superior cortex, thick inferior cortex) reflects the superior-inferior placement of trabecular bone (i.e., "trabecular eccentricity," TE). TE helps the femoral neck adapt to typical stresses and strains through load-sharing between superior and inferior cortices. Our goals were evaluated in the context of TE. Results showed the C/C paradigm has lower risk factors for the superior cortex and for the overall femoral neck, which is clinically relevant. TE analyses of the deer calcaneus revealed important synergism in load-sharing between the plantar/tension cortex and adjacent ligament/tendon, which challenges conventional understanding of how this control bone achieves functional adaptation. Comparisons with the control bone also exposed important deficiencies in current understanding of human femoral neck loading and its potential histocompositional adaptations.

Spatial Variation in Young Ovine Cortical Bone Properties

Significant effort continues to be made to understand whether differences exist in the structural, compositional, and mechanical properties of cortical bone subjected to different strain modes or magnitudes. We evaluated juvenile sheep femora (age = 4 months) from the anterior and posterior quadrants at three points along the diaphysis as a model system for variability in loading. Micro-CT scans (50 micron) were used to measure cortical thickness and mineral density. Three point bending tests were performed to measure the flexural modulus, strength, and post-yield displacement. There was no difference in cortical thickness or density between anterior or posterior quadrants; however, density was consistently higher in the middle diaphysis. Interestingly, bending modulus and strength were higher in anterior quadrants compared to posterior quadrants. Together, our results suggest that there is a differential spatial response of bone in terms of elastic bending modulus and mechanical strength. The origins of this difference may lie within the variation in ongoing mineralization, in combination with the collagen-rich plexiform structure, and whether this is related to strain mode remains to be explored. These data suggest that in young ovine cortical bone, modulation of strength occurs via potentially complex interactions of both mineral and collagen-components that may be different in regions of bone exposed to variable amounts of strain. Further work is needed to confirm the physiological load state of bone during growth to better elucidate the degree to which these variations are a function of the local mechanical environment.

Understanding compressive viscoelastic properties of additively manufactured PLA for bone-mimetic scaffold design

Bone tissue engineering has been recognized as a promising strategy to repair or replace damaged bone tissues. The mechanical properties of bone scaffolds play a critical role in successful bone regeneration, as it is essential to match the mechanical properties of the scaffold with the surrounding bone tissue. In this study, we investigated the effects of fused deposition modeling (FDM) process parameters, including printing speed, printing temperature, and layer thickness, on the compressive viscoelastic properties of polylactic acid (PLA) scaffolds. The compressive viscoelastic properties of bulk PLA specimens were characterized using a Zhu-Wang-Tang (ZWT) constitutive model under different compressive strain rates. A comprehensive statistical analysis comprising multivariate and univariate analysis of variance (MANOVA and ANOVA) and Tukey's post hoc analysis was utilized to quantify the effect of each FDM parameter on the viscoelastic mechanical properties of the PLA specimens. Subsequently, we fabricated modified face-centered cubic (MFCC) scaffolds using FDM and varied the FDM process parameters to achieve a compressive viscoelastic response that matched the natural trabecular bone tissue. The viscoelastic performance of the MFCC scaffolds was compared with traditional orthogonal cylindrical struts (OCS) scaffolds. Our methodology contributes to the design of bone-mimetic scaffolds with optimized mechanical properties by controlling FDM process parameters.

Above knee socket prosthesis use changes proximal femur morphology

Patients with transfemoral amputation (TFA) are up to six times more likely to develop hip osteoarthritis (OA) in either or both the intact and residual limb, which is primarily attributed to habitually altered joint loading due to compensatory movement patterns. However, joint loading patterns differ between limbs, which confounds the understanding of loading-induced OA etiology across limbs. It remains unknown if altered loading due to amputation results in bony shape changes at the hip, which is a known etiological factor in the development of hip OA. Retrospective computed tomography images were collected of the residual limb for 31 patients with unilateral TFA (13F/18M; age: 51.7 ± 9.9 y/o; time since amputation: 13.7 ± 12.4 years) and proximal femur for a control group of 29 patients (13F/16M; age: 42.0 ± 12.27 years) and used to create 3D geometries of the proximal femur. Femoral 3D geometric variation was quantified using statistical shape modeling (SSM), a computational tool which placed 2048 corresponding particles on each geometry. Independent modes of variation were created using principal component analysis. 2D radiographic measures of the proximal femur, including common measures such as α-angle, head neck offset, and neck shaft angle, were quantified on digitally reconstructed radiographs (DRRs). SSM results were then compared to 2D measures using Pearson correlation coefficients (r). Two-sample t-tests were used to determine if there were significant differences between the TFA and control group means of 2D radiographic measurements (p < 0.05). Patients with TFA had greater femoral head asphericity within the SSM, which was moderately correlated to head-neck offset (r = -0.54) and α-angle (r = 0.63), as well as greater trochanteric torsion, which was strongly correlated to the novel radiographic measure of trochanteric torsion (r = -0.78), compared to controls. For 2D measures, the neck-shaft angle was smaller in the TFA group compared to the control group (p = 0.01) while greater trochanter height was larger in the TFA group compared to the control group (p = 0.04). These results indicate altered loading from transfemoral prosthesis use changes proximal femur bony morphology, including femoral head asphericity and greater trochanter changes. Greater trochanter morphologic changes, though not a known factor to OA, affect moment arm and line of action of the primary hip abductors, the major muscles which contribute to joint loading and hip stability. Thus, chronic altered loading of the amputated limb hip, whether under- or overloading, results in bony changes to the proximal femur which may contribute to the etiological progression and development of OA.

Using Finite Element Modeling in Bone Mechanoadaptation

PURPOSE OF THE REVIEW

Bone adapts structure and material properties in response to its mechanical environment, a process called mechanoadpatation. For the past 50 years, finite element modeling has been used to investigate the relationships between bone geometry, material properties, and mechanical loading conditions. This review examines how we use finite element modeling in the context of bone mechanoadpatation.

RECENT FINDINGS

Finite element models estimate complex mechanical stimuli at the tissue and cellular levels, help explain experimental results, and inform the design of loading protocols and prosthetics. FE modeling is a powerful tool to study bone adaptation as it complements experimental approaches. Before using FE models, researchers should determine whether simulation results will provide complementary information to experimental or clinical observations and should establish the level of complexity required. As imaging technics and computational capacity continue increasing, we expect FE models to help in designing treatments of bone pathologies that take advantage of mechanoadaptation of bone.

Effect of Yi Jin Jing exercise plus Elastic Band Resistance exercise on overall bone mineral density in postmenopausal women

OBJECTIVES

This work aimed to study the effects of Yi Jin Jing plus Elastic Band Resistance exercise on bone mineral density at all parts of the body and bone metabolism index levels in postmenopausal women.

DESIGN

Randomized controlled trial.

METHODS

Forty postmenopausal women were randomly assigned equally to the exercise or to the control group. The control group maintained their lifestyle behaviors unaltered, whereas the exercise group received  Yi Jin Jing plus Elastic Band Resistance exercise. The primary outcome was overall bone mineral density at each part, and the secondary one was bone metabolism indicator levels and bone mineral density on both sides.

RESULTS

The results after six months showed increased bone mineral density at all parts of the body in the exercise group (spine, P = 0.002; thighs, lumbar, and whole body, P < 0.05) and decreased bone mineral density in the control group (trunk, pelvis, and spine, P < 0.01). In particular, the decrease and increase were greater on the non-preferred (left) side than on the right side. As for bone metabolism indexes, β-Crosslaps levels reduced (P = 0.016) and a significant increase in 1,25-(OH)2-D3 (P < 0.001) can be observed in the exercise group.

CONCLUSIONS

The results suggested that Yi Jin Jing plus Elastic Band Resistance exercise could delay the overall decrease of bone mineral density in postmenopausal women, especially on the non-preferred side. It also increased bone formation metabolite levels and inhibited bone resorption metabolite levels.

Modification to Mirels scoring system location component improves fracture prediction for metastatic disease of the proximal femur

BACKGROUND

Correctly identifying patients at risk of femoral fracture due to metastatic bone disease remains a clinical challenge. Mirels criteria remains the most widely referenced method with the advantage of being easily calculated but it suffers from poor specificity. The purpose of this study was to develop and evaluate a modified Mirels scoring system through scoring modification of the original Mirels location component within the proximal femur.

METHODS

Computational (finite element) experiments were performed to quantify strength reduction in the proximal femur caused by simulated lytic lesions at defined locations. Virtual spherical defects representing lytic lesions were placed at 32 defined locations based on axial (4 axial positions: neck, intertrochanteric, subtrochanteric or diaphyseal) and circumferential (8 circumferential: 45-degree intervals) positions. Finite element meshes were created, material property assignment was based on CT mineral density, and femoral head/greater trochanter loading consistent with stair ascent was applied. The strength of each femur with a simulated lesion divided by the strength of the intact femur was used to calculate the Location-Based Strength Fraction (LBSF). A modified Mirels location score was next defined for each of the 32 lesion locations with an assignment of 1 (LBSF > 75%), 2 (LBSF: 51-75%), and 3 (LBSF: 0-50%). To test the new scoring system, data from 48 patients with metastatic disease to the femur, previously enrolled in a Musculoskeletal Tumor Society (MSTS) cross-sectional study was used. The lesion location was identified for each case based on axial and circumferential location from the CT images and assigned an original (2 or 3) and modified (1,2, or 3) Mirels location score. The total score for each was then calculated. Eight patients had a fracture of the femur and 40 did not over a 4-month follow-up period. Logistic regression and decision curve analysis were used to explore relationships between clinical outcome (Fracture/No Fracture) and the two Mirels scoring methods.

RESULTS

The location-based strength fraction (LBSF) was lowest for lesions in the subtrochanteric and diaphyseal regions on the lateral side of the femur; lesions in these regions would be at greatest risk of fracture. Neck lesions located at the anterior and antero-medial positions were at the lowest risk of fracture. When grouped, neck lesions had the highest LBSF (83%), followed by intertrochanteric (72%), with subtrochanteric (50%) and diaphyseal lesions (49%) having the lowest LBSF. There was a significant difference (p < 0.0001) in LBSF between each axial location, except subtrochanteric and diaphyseal which were not different from each other (p = 0.96). The area under the receiver operator characteristic (ROC) curve using logistic regression was greatest for modified Mirels Score using site specific location of the lesion (Modified Mirels-ss, AUC = 0.950), followed by a modified Mirels Score using axial location of lesion (Modified Mirels-ax, AUC = 0.941). Both were an improvement over the original Mirels score (AUC = 0.853). Decision curve analysis was used to quantify the relative risks of identifying patients that would fracture (TP, true positives) and those erroneously predicted to fracture (FP, false positives) for the original and modified Mirels scoring systems. The net benefit of the scoring system weighed the benefits (TP) and harms (FP) on the same scale. At a threshold probability of fracture of 10%, use of the modified Mirels scoring reduced the number of false positives by 17-20% compared to Mirels scoring.

CONCLUSIONS

A modified Mirels scoring system, informed by detailed analysis of the influence of lesion location, improved the ability to predict impending pathological fractures of the proximal femur for patients with metastatic bone disease. Decision curve analysis is a useful tool to weigh costs and benefits concerning fracture risk and could be combined with other patient/clinical factors that contribute to clinical decision making.

Potential of Neuromuscular Electrical Stimulation as a Bone Loss Countermeasure in Microgravity

INTRODUCTION: For future long-duration spaceflight missions, additional methods of loading the skeleton may be required to supplement exercise to minimize bone loss. Neuromuscular electrical stimulation (NMES) can elicit muscular contractions that create strain on bone. However, the potential effectiveness of NMES on the proximal femur during disuse is not known.METHODS: We measured the maximum isometric force of NMES-induced contractions of the rectus femoris and the hamstrings of 10 subjects (5 male, 5 female), sitting with the hips and knees at 90 degrees of flexion. We employed 2-D biomechanical models of the knee and hip to estimate the hip joint reaction forces, applied these forces to a generic femur finite element analysis model, and qualitatively compared the peak principal strains of the proximal femoral neck to the peak strains modeled in previous studies for other forms of exercise.RESULTS: The average peak tensile/compressive strains were 1380 ± 719 µε/-2179 ± 1130 µε and 573 ± 345 µε/-900 ± 543 µε for the male and female subjects, respectively. While results varied between studies, the strains achieved during NMES generally were comparable to those achieved during walking or stairs, with some individuals matching higher intensity activities.DISCUSSION: This study demonstrated that isometric NMES contractions of the thigh muscles can create strain in the proximal femoral neck similar to that achieved during low impact activities. While NMES alone will unlikely create a sufficient daily strain stimulus to prevent bone loss, it will likely improve the current spaceflight countermeasures by adding more frequent loading throughout the day.Abitante TJ, Bouxsein ML, Duda KR, Newman DJ. Potential of neuromuscular electrical stimulation as a bone loss countermeasure in microgravity. Aerosp Med Hum Perform. 2022; 93(11):774–782.

Social stratification reflected in bone mineral density and stature: Spectral imaging and osteoarchaeological findings from medieval Norway

This study presents skeletal material from five medieval burial sites in Eastern Norway, confined to one royal burial church, one Dominican monastery, and three burial sites representing parish populations. We combine osteological analysis and Dual Energy X-Ray Absorptiometry, studying the remains of 227 individuals (102 females and 125 males) employing young, middle, and old adult age categories. The aim is to assess bone mineral density as a skeletal indicator of socioeconomic status including stature as a variable. We detected that socioeconomic status significantly affected bone mineral density and stature. Individuals of high status had higher bone mineral density (0.07 g/cm2, p = 0.003) and taller stature (1.85 cm, p = 0.017) than individuals from the parish population. We detected no significant relationship between young adult bone mineral density and socioeconomic status (p = 0.127 and 0.059 for females and males, respectively). For males, high young adult bone mineral density and stature varied concordantly in both status groups. In contrast, females of high status were significantly taller than females in the parish population (p = 0.011). Our findings indicate quite different conditions during growth and puberty for the two groups of females. The age-related pattern of bone variation also portrayed quite different trajectories for the two socioeconomic status groups of both sexes. We discuss sociocultural practices (living conditions during childhood and puberty, as well as nutritional and lifestyle factors in adult life), possibly explaining the differences in bone mineral density between the high-status and parish population groups. The observation of greater differences in bone mineral density and stature for females than males in the medieval society of Norway is also further discussed.

Multidirectional basketball activities load different regions of the tibia: A subject-specific muscle-driven finite element study

The tibia is a common site for bone stress injuries, which are believed to develop from microdamage accumulation to repetitive sub-yield strains. There is a need to understand how the tibia is loaded in vivo to understand how bone stress injuries develop and design exercises to build a more robust bone. Here, we use subject-specific, muscle-driven, finite element simulations of 11 basketball players to calculate strain and strain rate distributions at the midshaft and distal tibia during six activities: walking, sprinting, lateral cut, jumping after landing, changing direction from forward-to-backward sprinting, and changing direction while side shuffling. Maximum compressive strains were at least double maximum tensile strains during the stance phase of all activities. Sprinting and lateral cut had the highest compressive (-2,862 ± 662 με and -2,697 ± 495 με, respectively) and tensile (973 ± 208 με and 942 ± 223 με, respectively) strains. These activities also had the highest strains rates (peak compressive strain rate = 64,602 ± 19,068 με/s and 37,961 ± 14,210 με/s, respectively). Compressive strains principally occurred in the posterior tibia for all activities; however, tensile strain location varied. Activities involving a change in direction increased tensile loads in the anterior tibia. These observations may guide preventative and management strategies for tibial bone stress injuries. In terms of prevention, the strain distributions suggest individuals should perform activities involving changes in direction during growth to adapt different parts of the tibia and develop a more fatigue resistant bone. In terms of management, the greater strain and strain rates during sprinting than jumping suggests jumping activities may be commenced earlier than full pace running. The greater anterior tensile strains during changes in direction suggest introduction of these types of activities should be delayed during recovery from an anterior tibial bone stress injury, which have a high-risk of healing complications.

A Comparative Study on the Multiscale Mechanical Responses of Human Femoral Neck Between the Young and the Elderly Using Finite Element Method

Background: Femoral neck fracture (FNF) is the most serious bone disease in the elderly population. The multiscale mechanical response is a key to predicting the strength of the femoral neck, assessing the risk of FNF, and exploring the role of mechanosensation and mechanotransmission in bone remodeling, especially in the context of aging bone.

Methods: Multiscale finite element (FE) models of the proximal femur for both young and elderly people were developed. The models included organ scale (proximal femur), tissue scale (cortical bone), tissue element scale (osteon), and cell scale [osteocyte lacuna-canalicular network (LCN) and extracellular matrix (ECM), OLCEM]. The mechanical responses of cortical bone and osteocytes in the mid-femoral neck and the differences in mechanical responses between these two scales were investigated.

Results: The mechanical responses of cortical bone and osteocyte showed significant differences between the elderly and the young. The minimum principal strains and mean SEDs of cortical bone in the elderly were 2.067–4.708 times and 3.093–14.385 times of the values in the young, respectively; the minimum principal strains and mean SEDs of osteocyte in the elderly were 1.497–3.246 times and 3.044–12 times of the values in the young, respectively; the amplification factors of minimum principal strain in the inferior (Inf), anterior (Ant), and posterior (Post) quadrants in the young were 1.241–1.804 times of the values in the elderly, but the amplification factor of minimum principal strain in the superior (Sup) quadrant was 87.4% of the value in the elderly; the amplification factors of mean SED in the young were 1.124–9.637 times of the values in the elderly.

Conclusion: The mass and bone mineral density (BMD) of cortical bone in the femoral neck is closely related to the mechanical response of osteocytes, which provides a new idea for improving cortical bone quality. Perhaps cortical bone quality could be improved by stimulating osteocytes. Quadrantal differences of bone quality in the mid-femoral neck should be considered to improve fracture risk prediction in the future.

Effect of fall direction on the lower hip fracture risk in athletes with different loading histories: A finite element modeling study in multiple sideways fall configurations

Physical loading makes bones stronger through structural adaptation. Finding effective modes of exercise to improve proximal femur strength has the potential to decrease hip fracture risk. Previous proximal femur finite element (FE) modeling studies have indicated that the loading history comprising impact exercises is associated with substantially higher fracture load. However, those results were limited only to one specified fall direction. It remains thus unclear whether exercise-induced higher fracture load depends on the fall direction. To address this, using magnetic resonance images of proximal femora from 91 female athletes (mean age 24.7 years with >8 years competitive career) and their 20 non-athletic but physically active controls (mean age 23.7 years), proximal femur FE models were created in 12 different sideways fall configurations. The athletes were divided into five groups by typical loading patterns of their sports: high-impact (H-I: 9 triple- and 10 high-jumpers), odd-impact (O-I: 9 soccer and 10 squash players), high-magnitude (H-M: 17 powerlifters), repetitive-impact (R-I: 18 endurance runners), and repetitive non-impact (R-NI: 18 swimmers). Compared to the controls, the FE models showed that the H-I and R-I groups had significantly (p < 0.05) higher fracture loads, 11-17% and 22-28% respectively, in all fall directions while the O-I group had significantly 10-11% higher fracture loads in four fall directions. The H-M and R-NI groups did not show significant benefit in any direction. Also, the analyses of the minimum fall strength (MFS) among these multiple fall configurations confirmed significantly 15%, 11%, and 14% higher MFSs in these impact groups, respectively, compared to the controls. These results suggest that the lower hip fracture risk indicated by higher fracture loads in athletes engaged in high impact or repetitive impact sports is independent of fall direction whereas the lower fracture risk attributed to odd-impact exercise is more modest and specific to the fall direction. Moreover, in concordance with the literature, the present study also confirmed that the fracture risk increases if the impact is imposed on the more posterolateral aspect of the hip. The present results highlight the importance of engaging in the impact exercises to prevent hip fractures and call for retrospective studies to investigate whether specific impact exercise history in adolescence and young adulthood is also associated with lower incidence of hip fractures in later life.

Finite Element Assessment of Bone Fragility from Clinical Images

PURPOSE OF REVIEW

We re-evaluated clinical applications of image-to-FE models to understand if clinical advantages are already evident, which proposals are promising, and which questions are still open.

RECENT FINDINGS

CT-to-FE is useful in longitudinal treatment evaluation and groups discrimination. In metastatic lesions, CT-to-FE strength alone accurately predicts impending femoral fractures. In osteoporosis, strength from CT-to-FE or DXA-to-FE predicts incident fractures similarly to DXA-aBMD. Coupling loads and strength (possibly in dynamic models) may improve prediction. One promising MRI-to-FE workflow may now be tested on clinical data. Evidence of artificial intelligence usefulness is appearing. CT-to-FE is already clinical in opportunistic CT screening for osteoporosis, and risk of metastasis-related impending fractures. Short-term keys to improve image-to-FE in osteoporosis may be coupling FE with fall risk estimates, pool FE results with other parameters through robust artificial intelligence approaches, and increase reproducibility and cross-validation of models. Modeling bone modifications over time and bone fracture mechanics are still open issues.

Physical activity induced adaptation can increase proximal femur strength under loading from a fall onto the greater trochanter

Physical activity enhances proximal femur bone mass, but it remains unclear whether the benefits translate into an enhanced ability to resist fracture related loading. We recently used baseball pitchers as a within-subject controlled model to demonstrate physical activity induced proximal femur adaptation in regions associated with weight bearing and increased strength under single-leg stance loading. However, there was no measurable benefit to resisting common injurious loading (e.g. a fall onto the greater trochanter). A lack of power and a small physical activity effect size may have contributed to the latter null finding. Softball pitchers represent an alternative within-subject controlled model to explore adaptation of the proximal femur to physical activity, exhibiting greater dominant-to-nondominant (D-to-ND) leg differences than baseball pitchers. The current study used quantitative computed tomography, statistical parametric mapping, and subject-specific finite element (FE) modeling to explore adaptation of the proximal femur to physical activity in female softball pitchers (n = 25). Female cross-country runners (n = 15) were included as symmetrically loaded controls, showing very limited D-to-ND leg differences. Softball pitchers had D-to-ND leg differences in proximal femur, femoral neck, and trochanteric volumetric bone mineral density and content, and femoral neck volume. Voxel-based morphometry analyses and cortical bone mapping showed D-to-ND leg differences within a large region connecting the superior femoral head, inferior femoral neck and medial intertrochanteric region, and within the greater trochanter. FE modeling revealed pitchers had 19.4% (95%CI, 15.0 to 23.9%) and 4.9% (95%CI, 1.7 to 8.2%) D-to-ND leg differences in predicted ultimate strength under single-leg stance loading and a fall onto the greater trochanter, respectively. These data affirm the spatial and strength adaptation of the proximal femur to weight bearing directed loading and demonstrate that the changes can also have benefits, albeit smaller, on resisting loads associated with a sideways fall onto the greater trochanter.

Modeling the musculoskeletal loading in bone remodeling at the hip of a child

BACKGROUND AND OBJECTIVES

The mechanical load associated with physical activity affects the bone adaptation process. The bone adaptationeffect varies with age, being more effective during childhood and adolescence, particularly during pre-pubertal years. Bone-strengthening physical activity is recommended for children and adolescents. The number of time periods (bouts) per day of vigorous physical activity seems to be more important than the total cumulative time for optimal bone strength. So, the aim of this study was to evaluate the effects of weight-bearing physical activity on bone mineral density (BMD) of the proximal femur through computational simulation considering the intensity, exposure time (bouts) and regionalization of the results.

METHODS

For this purpose, a finite element model of a 7 year-old child femur was developed based on computed tomography images. Musculoskeletal loads were obtained from experimental kinematic data of weight-bearing physical activity performed by children of the same age (standing, walking, running, jumping). The effects of physical activity on BMD of several regions of interest of the femur were analyzed using a bone remodeling model. A daily accumulation of 400 min of physical activity (200 min walking and 200 min standing) was considered as reference, against with which the effects of additional 10 min loading bouts were compared: 10 min bouts of vigorous intensity physical activity vs. 10 min bouts of light to moderate intensity physical activity.

RESULTS

The simulations revealed greater increases in BMD associated with higher intensity and longer duration of physical activity. The largest BMD increases occurs during the first 10 min bout compared to longer durations and in less mineralized central regions compared to regions far from the neutral axis of the bone.

CONCLUSION

Weight bearing physical activity is more effective in bone remodeling when the musculoskeletal loading is more intense and of short duration and, under these conditions, less mineralized regions are more positively impacted.

Real-time replication of three-dimensional and time-varying physiological loading cycles for bone and implant testing: A novel protocol demonstrated for the proximal human femur while walking

In vitro real-time replication of three-dimensional, time-varying load profiles acting on human bones during physical activity can advance bone and implant testing protocols. This study aimed to develop a novel protocol for applying the three-dimensional, time-varying hip contact force while walking to a human femur specimen. The target force profile was obtained from the literature. A proximal femur from an elderly female donor was instrumented using ten rosette strain gages and tested using a custom-made hexapod robot. A load-control algorithm determined the robot position generating the target force at low frequency (0.0004 Hz). Five cycles of the robot position were played back at five intermediate frequencies up to real-time (0.04, 0.08, 0.16, 0.4, and 0.8 Hz). The hip reaction force, the length of the actuators (position), and cortical strains were compared. The error in the load-control force was 0.3 ± 4.2 N (mean ± SD). The last three force, position, and strain cycles varied by less than 1.1% for every frequency analyzed. Across frequencies, the force increased by 28% at 0.8 Hz as a logarithmic function of frequency (R² = 0.98). The position and strain error linearly increased with frequency up to 0.4 Hz. The median position error and the interquartile range of the strain error reached 15% and 13% at 0.8 Hz. Changes of force and cortical strain at increasing frequencies were linearly related (R² = 0.99). Therefore, the protocol developed can provide repeatable three-dimensional time-varying load profiles, although the comparison of the specimen deformation obtained across frequencies should be considered with care, particularly in the higher frequency range. This information supports the design of dynamic tests of bone and implants.

Finite element analysis of femoral neck strains during stair ascent and descent

For older population, a better understanding of the hip joint loading environment is needed for the prevention of hip pain, and the reduction of the stress fractures and fall risks. Using the motion analysis and inverse dynamics methods, combined with musculoskeletal modelling, static optimization, and finite element (FE) femur model, the difference of femoral neck strains between stair ascent vs. descent, young vs. older populations was compared. A two-way repeated-measures MANOVA was applied to test the effect of age and stair direction on the femoral neck strains. The strains at the femoral neck cross-section were greater for stair descent than ascent for both age groups (mostly P = 0.001 to 0006) but there was no difference between age groups. In this study, femoral neck strains represented bone loading environment in more direct ways than joint reaction forces/moments or joint contact forces, the risk of hip pain, falls and stress fractures might be greater during stair descent than ascent. Possible preventative methods to reduce these risks should be developed in the future studies.

Development of axial compression and combined axial compression and torque loading configurations to reproduce strain in the implanted femur during activities of daily living

Femoral strain is indicative of the potential for bone remodeling (strain energy density, SED) and periprosthetic femoral fracture (magnitude of principal strains) after total hip arthroplasty (THA). Previous modeling studies have evaluated femoral strains in THA-implanted femurs under gait loads including both physiological hip contact force and femoral muscle forces. However, experimental replication of the complex muscle forces during activities of daily living (ADLs) is difficult for in vitro assessment of femoral implant or fixation hardware. Alternatively, cadaveric tests using simplified loading configurations have been developed to assess post-THA bone mechanics, although no current studies have demonstrated simplified loading configurations used in mechanical tests may simulate the physiological femoral strains under ADL loads. Using an optimization approach integrated with finite element analysis, this study developed axial compression and combined axial compression and torque testing configurations for three common ADLs (gait, stair-descent and sit-to-stand) via matching the SED profile of the femur in THA-implanted models of three specimens. The optimized simplified-loading models showed good agreement in predicting bone remodeling stimuli (post-THA change in SED per unit mass) and fatigue regions as compared with the ADL-loading models, as well as other modeling and clinical studies. The optimized simplified test configurations can provide a physiological-loading based pre-clinical platform for the evaluation of implant/fixation devices of the femur.

Paggiosi M, McCloskey EV, Viceconti M, Mazzà C, Li X. Femoral neck strain prediction during level walking using a combined musculoskeletal and finite element model approach

Recently, coupled musculoskeletal-finite element modelling approaches have emerged as a way to investigate femoral neck loading during various daily activities. Combining personalised gait data with finite element models will not only allow us to study changes in motion/movement, but also their effects on critical internal structures, such as the femur. However, previous studies have been hampered by the small sample size and the lack of fully personalised data in order to construct the coupled model. Therefore, the aim of this study was to build a pipeline for a fully personalised multiscale (body-organ level) model to investigate the strain levels at the femoral neck during a normal gait cycle. Five postmenopausal women were included in this study. The CT and MRI scans of the lower limb, and gait data were collected for all participants. Muscle forces derived from the body level musculoskeletal models were used as boundary constraints on the finite element femur models. Principal strains were estimated at the femoral neck region during a full gait cycle. Considerable variation was found in the predicted peak strain among individuals with mean peak first principal strain of 0.24% ± 0.11% and mean third principal strain of -0.29% ± 0.24%. For four individuals, two overall peaks of the maximum strains were found to occur when both feet were in contact with the floor, while one individual had one peak at the toe-off phase. Both the joint contact forces and the muscular forces were found to substantially influence the loading at the femoral neck. A higher correlation was found between the predicted peak strains and the gluteus medius (R² ranged between 0.95 and 0.99) than the hip joint contact forces (R² ranged between 0.63 and 0.96). Therefore, the current findings suggest that personal variations are substantial, and hence it is important to consider multiple subjects before deriving general conclusions for a target population.

Postoperative lower limb joint kinematics following tibial plateau fracture: A 2-year longitudinal study

BACKGROUND

The goal of postoperative tibial plateau fracture (TPF) management is to ensure surgical fixation is maintained while returning patients to normal function as soon as possible, allowing patients to resume their normal activities of daily living. The aim of this study was to investigate longitudinal changes in lower limb joint kinematics following TPF and determine how these kinematics relate to self-reported function.

METHODS

Patients presenting with a TPF were recruited (n = 18) and undertook gait analysis at six postoperative time points (two weeks, six weeks, three months, six months, one and two years). Lower limb joint kinematics were assessed at each time point based on gait data. Statistical parametric mapping (SPM) was undertaken to investigate the change in joint kinematic traces with time. The Knee Injury and Osteoarthritis Outcome Score (KOOS) was assessed at each time point to obtain self-reported outcomes. A healthy reference was also analyzed and used for qualitative comparison of joint kinematics.

RESULTS AND SIGNIFICANCE

Knee kinematics showed improvements with time, however only minor changes were noted after six weeks at the hip, and six months at the knee and ankle relative to two weeks postoperative. SPM identified significant improvements with time in hip (p < 0.001) and knee (p = 0.003) flexion. No significant changes were observed with time at the ankle however, when compared to the healthy reference, participants showed reduced plantarflexion at two years. Lower limb joint ROM showed significant weak to moderate correlation with the ADL sub-scale of the KOOS (hip r = 0.442, knee r = 0.303, ankle r = 0.367). The lack of significant changes with time and overall reduced plantarflexion at the ankle potentially reduces propulsive capacity during gait up to two years postoperative. In this study, we see a deficiency in joint kinematics in TPF patients up to two years when compared to a healthy reference, especially at the ankle.

Real-Time Prediction of Joint Forces by Motion Capture and Machine Learning

Conventional biomechanical modelling approaches involve the solution of large systems of equations that encode the complex mathematical representation of human motion and skeletal structure. To improve stability and computational speed, being a common bottleneck in current approaches, we apply machine learning to train surrogate models and to predict in near real-time, previously calculated medial and lateral knee contact forces (KCFs) of 54 young and elderly participants during treadmill walking in a speed range of 3 to 7 km/h. Predictions are obtained by fusing optical motion capture and musculoskeletal modeling-derived kinematic and force variables, into regression models using artificial neural networks (ANNs) and support vector regression (SVR). Training schemes included either data from all subjects (LeaveTrialsOut) or only from a portion of them (LeaveSubjectsOut), in combination with inclusion of ground reaction forces (GRFs) in the dataset or not. Results identify ANNs as the best-performing predictor of KCFs, both in terms of Pearson R (0.89-0.98 for LeaveTrialsOut and 0.45-0.85 for LeaveSubjectsOut) and percentage normalized root mean square error (0.67-2.35 for LeaveTrialsOut and 1.6-5.39 for LeaveSubjectsOut). When GRFs were omitted from the dataset, no substantial decrease in prediction power of both models was observed. Our findings showcase the strength of ANNs to predict simultaneously multi-component KCF during walking at different speeds-even in the absence of GRFs-particularly applicable in real-time applications that make use of knee loading conditions to guide and treat patients.

Play During Growth: the Effect of Sports on Bone Adaptation

PURPOSE OF REVIEW

The development of exercise interventions for bone health requires an understanding of normative growth trends. Here, we summarize changes in bone during growth and the effect of participating in sports on structural and compositional measures in different bones in males and females.

RECENT FINDINGS

Growing females and males have similar normalized density and bone area fraction until age 16, after which males continue increasing at a faster rate than females. All metrics for both sexes tend to plateau or decline in the early 20s. Areal BMD measures indicate significant heterogeneity in adaptation to sport between regions of the body. High-resolution CT data indicate changes in structure are more readily apparent than changes in density. While adaptation to sport is spatially heterogeneous, participation in weight-bearing activities that involve dynamic muscle contractions tends to result in increased bone adaptation.

Cortical thinning and accumulation of large cortical pores in the tibia reflect local structural deterioration of the femoral neck

INTRODUCTION

Cortical bone thinning and a rarefaction of the trabecular architecture represent possible causes of increased femoral neck (FN) fracture risk. Due to X-ray exposure limits, the bone microstructure is rarely measurable in the FN of subjects but can be assessed at the tibia. Here, we studied whether changes of the tibial cortical microstructure, which were previously reported to be associated with femur strength, are also associated with structural deteriorations of the femoral neck.

METHODS

The cortical and trabecular architectures in the FN of 19 humans were analyzed ex vivo on 3D microcomputed tomography images with 30.3 μm voxel size. Cortical thickness (Ct.Th_tibia), porosity (Ct.Po_tibia) and pore size distribution in the tibiae of the same subjects were measured using scanning acoustic microscopy (12 μm pixel size). Femur strength during sideways falls was simulated with homogenized voxel finite element models.

RESULTS

Femur strength was associated with the total (vBMD_tot; R² = 0.23, p < 0.01) and trabecular (vBMD_trab; R² = 0.26, p < 0.01) volumetric bone mineral density (vBMD), with the cortical thickness (Ct.Th_FN; R² = 0.29, p < 0.001) and with the trabecular bone volume fraction (Tb.BV/TV_FN; R² = 0.34, p < 0.001), separation (Tb.Sp_FN; R² = 0.25, p < 0.01) and number (Tb.N_FN; R² = 0.32, p < 0.001) of the femoral neck. Moreover, smaller Ct.Th_tibia was associated with smaller Ct.Th_FN (R² = 0.31, p < 0.05), lower Tb.BV/TV_FN (R² = 0.29, p < 0.05), higher Tb.Sp_FN (R² = 0.33, p < 0.05) and lower Tb.N_FN (R² = 0.42, p < 0.01). A higher prevalence of pores with diameter > 100 μm in tibial cortical bone (relCt.Po_100μm-tibia) indicated higher Tb.Sp_FN (R² = 0.36, p < 0.01) and lower Tb.N_FN (R² = 0.45, p < 0.01).

CONCLUSION

Bone resorption and structural decline of the femoral neck may be identified in vivo by measuring cortical bone thickness and large pores in the tibia.

The relationship between tibiofemoral geometry and musculoskeletal function during normal activity

BACKGROUND

The effect of tibiofemoral geometry on musculoskeletal function is important to movement biomechanics.

RESEARCH QUESTION

We hypothesised that tibiofemoral geometry determines tibiofemoral motion and musculoskeletal function. We then aimed at 1) modelling tibiofemoral motion during normal activity as a function of tibiofemoral geometry in healthy adults; and 2) quantifying the effect of tibiofemoral geometry on musculoskeletal function.

METHODS

We used motion data for six activity types and CT images of the knee from 12 healthy adults. Geometrical variation of the tibia and femoral articular surfaces were measured in the CT images. The geometry-based tibiofemoral motion was calculated by fitting a parallel mechanism to geometrical variation in the cohort. Matched musculoskeletal models embedding the geometry-based tibiofemoral joint motion and a common generic tibiofemoral motion of reference were generated and used to calculate joint angles, net joint moments, muscle and joint forces for the six activities analysed. The tibiofemoral model was validated against bi-planar fluoroscopy measurements for walking for all the six planes of motion. The effect of tibiofemoral geometry on musculoskeletal function was the difference between the geometry-based model and the model of reference.

RESULTS

The geometry-based tibiofemoral motion described the pattern and the variation during walking for all six motion components, except the pattern of anterior tibial translation. Tibiofemoral geometry had moderate effect on cohort-averages of musculoskeletal function (R² = 0.60-1), although its effect was high in specific instances of the model, outputs and activities analysed, reaching 2.94 BW for the ankle reaction force during stair descent. In conclusion, tibiofemoral geometry is a major determinant of tibiofemoral motion during walking.

SIGNIFICANCE

Geometrical variations of the tibiofemoral joint are important for studying musculoskeletal function during normal activity in specific individuals but not for studying cohort averages of musculoskeletal function. This finding expands current knowledge of movement biomechanics.

Cortical bone mapping improves finite element strain prediction accuracy at the proximal femur

Despite evidence of the biomechanical role of cortical bone, current state of the art finite element models of the proximal femur built from clinical CT data lack a subject-specific representation of the bone cortex. Our main research hypothesis is that the subject-specific modelling of cortical bone layer from CT images, through a deconvolution procedure known as Cortical Bone Mapping (CBM, validated for cortical thickness and density estimates) can improve the accuracy of CT-based FE models of the proximal femur, currently limited by partial volume artefacts. Our secondary hypothesis is that a careful choice of cortical-specific density-elasticity relationship may improve model accuracy. We therefore: (i) implemented a procedure to include subject-specific CBM estimates of both cortical thickness and density in CT-based FE models. (ii) defined alternative models that included CBM estimates and featured a cortical-specific or an independently optimised density-elasticity relationship. (iii) tested our hypotheses in terms of elastic strain estimates and failure load and location prediction, by comparing with a published cohort of 14 femurs, where strain and strength in stance and fall loading configuration were experimentally measured, and estimated through reference FE models that did not explicitly model the cortical compartment. Our findings support the main hypothesis: an explicit modelling of the proximal femur cortical bone layer including CBM estimates of cortical bone thickness and density increased the FE strains prediction, mostly by reducing peak errors (average error reduced by 30%, maximum error and 95th percentile of error distribution halved) and especially when focusing on the femoral neck locations (all error metrics at least halved). We instead rejected the secondary hypothesis: changes in cortical density-elasticity relationship could not improve validation performances. From these improved baseline strain estimates, further work is needed to achieve accurate strength predictions, as models incorporating cortical thickness and density produced worse estimates of failure load and equivalent estimates of failure location when compared to reference models. In summary, we recommend including local estimates of cortical thickness and density in FE models to estimate bone strains in physiological conditions, and especially when designing exercise studies to promote bone strength.

Heterogeneous Spatial and Strength Adaptation of the Proximal Femur to Physical Activity: A Within-Subject Controlled Cross-Sectional Study

Physical activity (PA) enhances proximal femur bone mass, as assessed using projectional imaging techniques. However, these techniques average data over large volumes, obscuring spatially heterogeneous adaptations. The current study used quantitative computed tomography, statistical parameter mapping, and subject-specific finite element (FE) modeling to explore spatial adaptation of the proximal femur to PA. In particular, we were interested in adaptation occurring at the superior femoral neck and improving strength under loading from a fall onto the greater trochanter. High/long jump athletes (n = 16) and baseball pitchers (n = 16) were utilized as within-subject controlled models as they preferentially load their take-off leg and leg contralateral to their throwing arm, respectively. Controls (n = 15) were included but did not show any dominant-to-nondominant (D-to-ND) leg differences. Jumping athletes showed some D-to-ND leg differences but less than pitchers. Pitchers had 5.8% (95% confidence interval [CI] 3.9%-7.6%) D-to-ND leg differences in total hip volumetric bone mineral density (vBMD), with increased vBMD in the cortical compartment of the femoral neck and trochanteric cortical and trabecular compartments. Voxel-based morphometry analyses and cortical bone mapping showed pitchers had D-to-ND leg differences within the regions of the primary compressive trabeculae, inferior femoral neck, and greater trochanter but not the superior femoral neck. FE modeling revealed pitchers had 4.1% (95% CI 1.4%-6.7%) D-to-ND leg differences in ultimate strength under single-leg stance loading but no differences in ultimate strength to a fall onto the greater trochanter. These data indicate the asymmetrical loading associated with baseball pitching induces proximal femur adaptation in regions associated with weight bearing and muscle contractile forces and increases strength under single-leg stance loading. However, there were no benefits evident at the superior femoral neck and no measurable improvement in ultimate strength to common injurious loading during aging (ie, fall onto the greater trochanter), raising questions as to how to better target these variables with PA. © 2019 American Society for Bone and Mineral Research.

Population-Based Bone Strain During Physical Activity: A Novel Method Demonstrated for the Human Femur

Statistical methods are increasingly used in biomechanics for studying bone geometry, bone density distribution and function in the population. However, relating population-based bone variation to strain during activity is computationally challenging. Here, we describe a new method for calculating strain in a population, using the Superposition Principle Method Squared (SPM²), and we demonstrate the method for calculating strain in human femurs. Computed-tomography images and motion capture while walking in 21 healthy adult women were obtained earlier. Variation of femur geometry and bone distribution were modelled using active shape and appearance modelling (ASAM). Femoral strain was modelled as the weighted sum of strain generated by each force in the model plus a strain variation assumed a quadratic function of the ASAM scores. The quadratic coefficients were fitted to 35 instances drawn from the ASAM model by varying each eigenmode by ± 2 SD. The equivalent strain in matched finite-element and SPM² calculations was obtained for 40 frames of walking for three independent cases and 50 ASAM instances. Finite-element and SPM2 solutions for walking were obtained in 44 and 3 min respectively. The SPM2 model accurately predicted strain for the three independent instances (R-squared 0.83-0.94) and the 50 ASAM instances (R-squared 0.95-1.00). The method developed enables fast and accurate calculation of population-based femoral strain.

Sex-Specific Muscular Mediation of the Relationship Between Physical Activity and Cortical Bone in Young Adults

Muscle mass is a commonly cited mediator of the relationship between physical activity (PA) and bone, representing the mechanical forces generated during PA. However, neuromuscular properties (eg, peak force) also account for unique portions of variance in skeletal outcomes. We used serial multiple mediation to explore the intermediary role of muscle mass and force in the relationships between cortical bone and moderate-to-vigorous intensity PA (MVPA). In a cross-sectional sample of young adults (n = 147, 19.7 ± 0.7 years old, 52.4% female) cortical diaphyseal bone was assessed via peripheral quantitative computed tomography at the mid-tibia. Peak isokinetic torque in knee extension was assessed via Biodex dynamometer. Thigh fat-free soft tissue (FFST) mass, assessed via dual-energy X-ray absorptiometry, represented the muscular aspect of tibial mechanical forces. Habitual MVPA was assessed objectively over 7 days using Actigraph GT3X+ accelerometers. Participants exceeded MVPA guidelines (89.14 ± 27.29 min/day), with males performing 44.5% more vigorous-intensity activity relative to females (p < 0.05). Males had greater knee extension torque and thigh FFST mass compared to females (55.3%, and 34.2%, respectively, all p < 0.05). In combined-sex models, controlling for tibia length and age, MVPA was associated with strength strain index (pSSI) through two indirect pathways: (i) thigh FFST mass (b = 1.11 ± 0.37; 95% CI, 0.47 to 1.93), and (i) thigh FFST mass and knee extensor torque in sequence (b = 0.30 ± 0.16; 95% CI, 0.09 to 0.73). However, in sex-specific models MVPA was associated with pSSI indirectly through its relationship with knee extensor torque in males (b = 0.78 ± 0.48; 95% CI, 0.04 to 2.02) and thigh FFST mass in females (b = 1.12 ± 0.50; 95% CI, 0.37 to 2.46). Bootstrapped CIs confirmed these mediation pathways. The relationship between MVPA and cortical structure appears to be mediated by muscle in young adults, with potential sex-differences in the muscular pathway. If confirmed, these findings may highlight novel avenues for the promotion of bone strength in young adults. © 2019 American Society for Bone and Mineral Research.

Temporal and spatial changes in bone accrual, density, and strain energy density in growing foals

Bone adaptation is in part driven by mechanical loading, and exercise during youth has been shown to have life-long benefits for bone health. However, the development of early exercise-based interventions that reduce the incidence of fractures in racing horses is limited by the lack of characterization of normal development in growing bone. Previous efforts to quantify bone development in the horse have relied on repeated radiographs or peripheral quantitative computed tomography scans, which are limited in their assessment of the entire bone. In this study, we acquired computed tomography scans of three Standardbred trotting colts longitudinally between 2 and 12 months of age. Finite-element models were constructed of the left forelimb proximal phalanx and used to assess strain energy density during quiet standing. Growth related changes in mineral density and bone area fraction in the distal epiphysis, mid-diaphysis, and proximal epiphysis were evaluated. Mineral density and bone area fraction uniformly increased in the diaphysis and strain energy density was constant during growth, indicating adaptation to quiet standing. Bone mineral density and bone area fraction increased in the medial quadrant of the proximal epiphysis but not in the fracture-prone lateral quadrant. The data presented provides a benchmark of normal growth trajectories that can be used to evaluate the effect of training regimens during growth.

Rat bone properties and their relationship to gait during growth

Allometric relationships have been studied over different Orders of mammals to understand how bone accommodates the mechanical demands associated with increasing mass. However, less attention has been given to the scaling of bone within a single lifetime. We aimed to determine how bone morphology and tissue density are related to (1) bending and compressive strength, and (2) gait dynamics. Longitudinal in vivo computed tomography of the hindlimbs and gait data were collected from female rats (n=5, age 8-20 weeks). Cross-sectional properties and tissue density were measured at the diaphysis, distal and proximal regions of the tibia and scaling exponents were calculated. Finite element models of the tibia were used to simulate loading during walking using joint forces from inverse dynamics calculation to determine the strain energy density and longitudinal strain at the midshaft. Second moment of area at the diaphysis followed strain similarity-based allometry, while bone area trended towards positive allometry. Strain energy in the diaphysis under transverse loading was lower than axial loading throughout growth. While both axial and transverse loading resulted in bending, tensile strains were mitigated by a change in the neutral axis and resulted in overall lower longitudinal tensile strains. The tissue density and cross-sectional properties initially increased and converged by 11 weeks of age and were correlated with changes in ground reaction forces. The scaling analyses imply that rodent tibia is (re)modeled in order to sustain bending at the midshaft during growth. The finite element results and relatively constant density after 10 weeks of age indicate that structural parameters may be the primary determinant of bone strength in the growing rodent tibia. The correlations between bone properties and joint angles imply that the changes in posture may affect bone growth in specific regions.

Adaptation of the proximal humerus to physical activity: A within-subject controlled study in baseball players

The proximal humerus is a common, yet understudied site for osteoporotic fracture. The current study explored the impact of prolonged physical activity on proximal humerus bone health by comparing bone properties between the throwing and nonthrowing arms within professional baseball players. The proximal humerus in throwing arms had 28.1% (95% CI, 17.8 to 38.3%) greater bone mass compared to nonthrowing arms, as assessed using dual-energy x-ray absorptiometry. At the level of the surgical neck, computed tomography revealed 12.0% (95% CI, 8.2 to 15.8%) greater total cross-sectional area and 31.0% (95% CI, 17.8 to 44.2%) greater cortical thickness within throwing arms, which contributed to 56.8% (95% CI, 44.9 to 68.8%) greater polar moment of inertia (i.e., estimated ability to resist torsional forces) compared to nonthrowing arms. Within the humeral head and greater tubercle regions, throwing arms had 3.1% (95% CI, 1.1 to 5.1%) more trabecular bone, as assessed using high-resolution magnetic resonance imaging. Three-dimensional mapping of voxel- and vertex-wise differences between arms using statistical parametric mapping techniques revealed throwing arms had adaptation within much of the proximal diaphysis, especially the posterolateral cortex. The pattern of proximal diaphysis adaptation approximated the pattern of strain energy distribution within the proximal humerus during a fastball pitch derived from a musculoskeletal and finite element model in a representative player. These data demonstrate the adaptive ability of the proximal humerus to physical activity-related mechanical loads. It remains to be established how they translate to exercise prescription to improve bone health within the proximal humerus; however, they provide unique insight into the relationship between prolonged loading and skeletal adaptation at a clinically relevant osteoporotic site.

A novel training-free method for real-time prediction of femoral strain

Surrogate methods for rapid calculation of femoral strain are limited by the scope of the training data. We compared a newly developed training-free method based on the superposition principle (Superposition Principle Method, SPM) and popular surrogate methods for calculating femoral strain during activity. Finite-element calculations of femoral strain, muscle, and joint forces for five different activity types were obtained previously. Multi-linear regression, multivariate adaptive regression splines, and Gaussian process were trained for 50, 100, 200, and 300 random samples generated using Latin Hypercube (LH) and Design of Experiment (DOE) sampling. The SPM method used weighted linear combinations of 173 activity-independent finite-element analyses accounting for each muscle and hip contact force. Across the surrogate methods, we found that 200 DOE samples consistently provided low error (RMSE < 100 µε), with model construction time ranging from 3.8 to 63.3 h and prediction time ranging from 6 to 1236 s per activity. The SPM method provided the lowest error (RMSE = 40 µε), the fastest model construction time (3.2 h) and the second fastest prediction time per activity (36 s) after Multi-linear Regression (6 s). The SPM method will enable large numerical studies of femoral strain and will narrow the gap between bone strain prediction and real-time clinical applications.

Biomechanics of Femoral Neck Fractures and Implications for Fixation

Fractures of the femoral neck can occur in young healthy individuals due to high loads occurring during motor vehicle accidents, impacts, or falls. Failure forces are lower if impacts occur sideways onto the greater trochanter as compared with vertical loading of the hip. Bone density, bone geometry, and thickness of cortical bone at the femoral neck contribute to its mechanical strength. Femoral neck fractures in young adults require accurate reduction and stable internal fixation. The available techniques for fracture fixation at the femoral neck (cannulated screws, hip screw systems, proximal femur plates, and cephallomedullary nails) are reviewed with respect to their competence to provide biomechanical stability. Mechanically unstable fractures require a load-bearing implant, such as hip screws, with antirotational screws or intramedullary nails. Subcapital or transcervical fracture patterns and noncomminuted fractures enable load sharing and can be securely fixed with cannulated screws or solitary hip screw systems without compromising fixation stability.

Efficacy and efficiency of multivariate linear regression for rapid prediction of femoral strain fields during activity

Multivariate Linear Regression-based (MLR) surrogate models were explored to reduce the computational cost of predicting femoral strains during normal activity in comparison with finite element analysis. The musculoskeletal model of one individual, the finite-element model of the right femur, and experimental force and motion data for normal walking, fast walking, stair ascent, stair descent, and rising from a chair were obtained from a previous study. Equivalent Von Mises strain was calculated for 1000 frames uniformly distributed across activities. MLR surrogate models were generated using training sets of 50, 100, 200 and 300 samples. The finite-element and MLR analyses were compared using linear regression. The Root Mean Square Error (RMSE) and the 95th percentile of the strain error distribution were used as indicators of average and peak error. The MLR model trained using 200 samples (RMSE < 108 µε; peak error < 228 µε) was used as a reference. The finite-element method required 66 s per frame on a standard desktop computer. The MLR model required 0.1 s per frame plus 1848 s of training time. RMSE ranged from 1.2% to 1.3% while peak error ranged from 2.2% to 3.6% of the maximum micro-strain (5020 µε). Performance within an activity was lower during early and late stance, with RMSE of 4.1% and peak error of 8.6% of the maximum computed micro-strain. These results show that MLR surrogate models may be used to rapidly and accurately estimate strain fields in long bones during daily physical activity.

Exercise for the prevention of osteoporosis in postmenopausal women: an evidence-based guide to the optimal prescription

BACKGROUND

Osteoporosis and related fragility fractures are a global public health problem in which pharmaceutical agents targeting bone mineral density (BMD) are the first line of treatment. However, pharmaceuticals have no effect on improving other key fracture risk factors, including low muscle strength, power and functional capacity, all of which are associated with an increased risk for falls and fracture, independent of BMD. Targeted exercise training is the only strategy that can simultaneously improve multiple skeletal and fall-related risk factors, but it must be appropriately prescribed and tailored to the desired outcome(s) and the specified target group.

OBJECTIVES

In this review, we provide an overview of the general principles of training and specific loading characteristics underlying current exercise guidelines for the prevention of osteoporosis, and an update on the latest scientific evidence with regard to the type and dose of exercise shown to positively influence bone mass, structure and strength and reduce fracture risk in postmenopausal women.