Three-Dimensional Osseointegration Patterns of Cementless Femoral Stems: An ex Vivo Study with High-Resolution Imaging and Histological Evaluation

BACKGROUND

Osseointegration is essential for the long-term survival of cementless femoral stems and is dependent on periprosthetic bone quality and correct implantation technique. The aim of this study was to evaluate the 3-dimensional long-term fixation patterns of, and bone microarchitecture around, cementless hip stems.

METHODS

Four specimens with varying degrees of bone quality and fixation characteristics from body donors who had received Alloclassic Zweymüller hip stems during their lifetime (mean time in situ at the time of death: 12.73 years) were evaluated with use of radiographs, high-resolution computed tomography (CT) scans, and hard-tissue histology. The CT voxel size was 85 µm, and the following parameters were calculated: total bone volume, total bone volume fraction, trabecular bone volume, trabecular bone volume fraction, cortical bone volume, cortical bone volume fraction, and cortical thickness. Bone-implant contact and canal fill index values for each Gruen zone of the specimens were calculated with use of histological samples.

RESULTS

Femoral stems with apparently good cortical contact on clinical radiographs showed higher values for cortical bone volume, trabecular bone volume, and cortical thickness in the high-resolution CT analysis than femoral stems with apparently weak cortical contact on clinical radiographs. Based on the histological evaluation, the mean bone-implant contact ranged from 22.94% to 57.24% and the mean canal fill index ranged from 52.33% to 69.67% among the specimens.

CONCLUSIONS

This study demonstrated different osseointegration patterns of cementless femoral stems on the basis of radiographs, high-resolution CT scans, and histological evaluation. Femora with high cortical bone volume and cortical thickness were associated with higher canal fill indices, whereas femora with low cortical bone volume and cortical thickness had lower canal fill indices and showed a characteristic corner-anchorage pattern.

CLINICAL RELEVANCE

Osseointegration patterns and thus the long-term survival of cementless femoral stems are dependent on cortical bone volume and cortical thickness.

Development of a density-based topology optimization of homogenized lattice structures for individualized hip endoprostheses and validation using micro-FE

Prosthetic implants, particularly hip endoprostheses, often lead to stress shielding because of a mismatch in compliance between the bone and the implant material, adversely affecting the implant's longevity and effectiveness. Therefore, this work aimed to demonstrate a computationally efficient method for density-based topology optimization of homogenized lattice structures in a patient-specific hip endoprosthesis. Thus, the root mean square error (RMSE) of the stress deviations between the physiological femur model and the optimized total hip arthroplasty (THA) model compared to an unoptimized-THA model could be reduced by 81 % and 66 % in Gruen zone (GZ) 6 and 7. However, the method relies on homogenized finite element (FE) models that only use a simplified representation of the microstructural geometry of the bone and implant. The topology-optimized hip endoprosthesis with graded lattice structures was synthesized using algorithmic design and analyzed in a virtual implanted state using micro-finite element (micro-FE) analysis to validate the optimization method. Homogenized FE and micro-FE models were compared based on averaged von Mises stresses in multiple regions of interest. A strong correlation (CCC > 0.97) was observed, indicating that optimizing homogenized lattice structures yields reliable outcomes. The graded implant was additively manufactured to ensure the topology-optimized result's feasibility.

The deep trabecular structure of first metacarpals in extant hominids

OBJECTIVES

Recent studies have associated subarticular trabecular bone distribution in the extant hominid first metacarpal (Mc1) with observed thumb use, to infer fossil hominin thumb use. Here, we analyze the entire Mc1 to test for interspecific differences in: (1) the absolute volume of trabecular volume fraction, (2) the distribution of the deeper trabecular network, and (3) the distribution of trabeculae in the medullary cavity, especially beneath the Mc1 disto-radial flange.

MATERIALS AND METHODS

Trabecular bone was imaged using micro-computed tomography in a sample of Homo sapiens (n = 11), Pan paniscus (n = 10), Pan troglodytes (n = 11), Gorilla gorilla (n = 10) and Pongo sp., (n = 7). Using Canonical Holistic Morphometric Analysis (cHMA), we tested for interspecific differences in the trabecular bone volume fraction (BV/TV) and its relative distribution (rBV/TV) throughout the Mc1, including within the head, medullary cavity, and base.

RESULTS

P. paniscus had the highest, and H. sapiens the lowest, BV/TV relative to other species. rBV/TV distribution statistically distinguished the radial concentrations and lack of medullary trabecular bone in the H. sapiens Mc1 from all other hominids. H. sapiens and, to a lesser extent, G. gorilla also had a significantly higher trabecular volume beneath the disto-radial flange relative to other hominids.

DISCUSSION

These results are consistent with differences in observed thumb use in these species and may also reflect systemic differences in bone volume fraction. The trabecular bone extension into the medullary cavity and concentrations beneath the disto-radial flange may represent crucial biomechanical signals that will aid in the inference of fossil hominin thumb use.

Trabecular architecture of the distal femur in extant hominids

Extant great apes are characterized by a wide range of locomotor, postural and manipulative behaviours that each require the limbs to be used in different ways. In addition to external bone morphology, comparative investigation of trabecular bone, which (re-)models to reflect loads incurred during life, can provide novel insights into bone functional adaptation. Here, we use canonical holistic morphometric analysis (cHMA) to analyse the trabecular morphology in the distal femoral epiphysis of Homo sapiens (n = 26), Gorilla gorilla (n = 14), Pan troglodytes (n = 15) and Pongo sp. (n = 9). We test two predictions: (1) that differing locomotor behaviours will be reflected in differing trabecular architecture of the distal femur across Homo, Pan, Gorilla and Pongo; (2) that trabecular architecture will significantly differ between male and female Gorilla due to their different levels of arboreality but not between male and female Pan or Homo based on previous studies of locomotor behaviours. Results indicate that trabecular architecture differs among extant great apes based on their locomotor repertoires. The relative bone volume and degree of anisotropy patterns found reflect habitual use of extended knee postures during bipedalism in Homo, and habitual use of flexed knee posture during terrestrial and arboreal locomotion in Pan and Gorilla. Trabecular architecture in Pongo is consistent with a highly mobile knee joint that may vary in posture from extension to full flexion. Within Gorilla, trabecular architecture suggests a different loading of knee in extension/flexion between females and males, but no sex differences were found in Pan or Homo, supporting our predictions. Inter- and intra-specific variation in trabecular architecture of distal femur provides a comparative context to interpret knee postures and, in turn, locomotor behaviours in fossil hominins.

Fracture pattern analysis of fractures to the diaphysis of the fifth metatarsal

BACKGROUND

Fractures to the fifth's metatarsal (MT-V) diaphysis are common. These are inconsistently referred to as diaphyseal-, shaft-, or Dancer's fractures. A comprehensive analysis of the MT-V fracture morphology is missing. The aim was to qualitatively and quantitatively analyze fracture patterns of MT-V diaphyseal fractures.

HYPOTHESIS

Fractures to the shaft of the fifth metatarsal feature specific fracture morphologies.

MATERIALS AND METHODS

Retrospective, radiologic database study. Included were all acute, isolated MT-V shaft fractures (including the proximal [Lawrence and Botte (L&B) III] and distal meta-diaphysis). Demographics and fracture characteristics were assessed. Each proximal fracture line was drawn, scaled, and a qualitative and quantitative fracture line analysis was conducted. The quantitative fracture line analysis aimed at identifying dens clusters with arbitrary shape using the DBSCAN algorithm. Data are presented as mean±standard deviation.

RESULTS

Out of 704 eligible MT-V fractures, 156 met the inclusion criteria. Patient's mean age was 46±19 years and 94% suffered a low energy trauma. Qualitative and quantitative fracture line analysis revealed three distinct fracture patterns. The proximal (30%) and distal (5%) meta-diaphyseal clusters showed a predominant transverse fracture pattern. The vast majority of diaphyseal fractures (56%) were spiral/oblique fractures, progressing from the proximal lateral meta-diaphyseal region in an oblique course at 61±9° to the medial distal diaphyseal cortex. Seven percent of diaphyseal fractures showed a transverse fracture pattern.

DISCUSSION

Based on a qualitative and quantitative analysis of all MT-V shaft fractures, three distinct fracture clusters were identified with homogeneous fracture patterns. MT-V shaft fractures should therefore be classified as proximal meta-diaphyseal (L&B Type III), diaphyseal (oblique or transverse) and distal meta-diaphyseal.

LEVEL OF PROOF

IV; retrospective database study.

Feasibility of aluminum phantom radiography for osteoporosis detection in postmenopausal women with a fragility fracture of the distal radius compared to DXA and HR-pQCT

Recently, promising results have been reported for detection of osteoporosis with use of an aluminum phantom. Therefore, the aim of this study was to evaluate the feasibility of radiography-based bone mineral density (BMD) measurement using a graded aluminum phantom. This study included 27 postmenopausal women with a distal radius fracture. Aluminum phantom radiography of the healthy radius was conducted as well as high-resolution peripheral quantitative computed tomography (HR-pQCT) measurement of the ultradistal radius and dual energy X-ray absorptiometry (DXA) of the radius, spine, and hip. A strong correlation was observed between aluminum phantom radiography-based mean gray value (mGV) and DXA-derived BMD, especially for the ultradistal radius (ρ = 0.75; p < 0.001). A moderate correlation for the femoral neck (ρ = 0.61 and p < 0.001) between modalities was found. Radius mGV and HR-pQCT-derived BMD only showed a moderate correlation (ρ = 0.48; p < 0.09). Aluminum phantom radiography might serve as a cost efficient, highly available, low-radiation dose screening, and diagnostic method for osteoporosis additively to DXA measurements. Especially, an application in areas with constrained DXA availability and such as preoperative trauma settings would be beneficial. However, further investigation and assessment of specificity and sensitivity is needed.

Hip joint load prediction using inverse bone remodeling with homogenized FE models: Comparison to micro-FE and influence of material modeling strategy

BACKGROUND AND OBJECTIVE

Measuring physiological loading conditions in vivo can be challenging, as methods are invasive or pose a high modeling effort. However, the physiological loading of bones is also imprinted in the bone microstructure due to bone (re)modeling. This information can be retrieved by inverse bone remodeling (IBR). Recently, an IBR method based on micro-finite-element (µFE) modeling was translated to homogenized-FE (hFE) to decrease computational effort and tested on the distal radius. However, this bone has a relatively simple geometry and homogeneous microstructure. Therefore, the objective of this study was to assess the agreement of hFE-based IBR with µFE-based IBR to predict hip joint loading from the head of the femur; a bone with more complex loading as well as more heterogeneous microstructure.

METHODS

hFE-based IBR was applied to a set of 19 femoral heads using four different material mapping laws. One model with a single homogeneous material for both trabecular and cortical volume and three models with a separated cortex and either homogeneous, density-dependent inhomogeneous, or density and fabric-dependent orthotropic material. Three different evaluation regions (full bone, trabecular bone only, head region only) were defined, in which IBR was applied. µFE models were created for the same bones, and the agreement of the predicted hip joint loading history obtained from hFE and µFE models was evaluated. The loading history was discretized using four unit load cases.

RESULTS

The computational time for FE solving was decreased on average from 500 h to under 1 min (CPU time) when using hFE models instead of µFE models. Using more information in the material model in the hFE models led to a better prediction of hip joint loading history. Inhomogeneous and inhomogeneous orthotropic models gave the best agreement to µFE-based IBR (RMSE% <14%). The evaluation region only played a minor role.

CONCLUSIONS

hFE-based IBR was able to reconstruct the dominant joint loading of the femoral head in agreement with µFE-based IBR and required considerably lower computational effort. Results indicate that cortical and trabecular bone should be modeled separately and at least density-dependent inhomogeneous material properties should be used with hFE models of the femoral head to predict joint loading.

Accuracy of osseointegrated screw-bone construct stiffness and peri-implant loading predicted by homogenized FE models relative to micro-FE models

Computational predictions of stiffness and peri-implant loading of screw-bone constructs are highly relevant to investigate and improve bone fracture fixations. Homogenized finite element (hFE) models have been used for this purpose in the past, but their accuracy has been questioned given the numerous simplifications, such as neglecting screw threads and modelling the trabecular bone structure as a continuum. This study aimed to investigate the accuracy of hFE models of an osseointegrated screw-bone construct when compared to micro-FE models considering the simplified screw geometry and different trabecular bone material models. Micro-FE and hFE models were created from 15 cylindrical bone samples with a virtually inserted, osseointegrated screw (fully bonded interface). Micro-FE models were created including the screw with threads (=reference models) and without threads to quantify the error due to screw geometry simplification. In the hFE models, the screws were modelled without threads and four different trabecular bone material models were used, including orthotropic and isotropic material derived from homogenization with kinematic uniform boundary conditions (KUBC), as well as from periodicity-compatible mixed uniform boundary conditions (PMUBC). Three load cases were simulated (pullout, shear in two directions) and errors in the construct stiffness and the volume average strain energy density (SED) in the peri-implant region were evaluated relative to the micro-FE model with a threaded screw. The pooled error caused by only omitting screw threads was low (max: 8.0%) compared to the pooled error additionally including homogenized trabecular bone material (max: 92.2%). Stiffness was predicted most accurately using PMUBC-derived orthotropic material (error: -0.7 ± 8.0%) and least accurately using KUBC-derived isotropic material (error: +23.1 ± 24.4%). Peri-implant SED averages were generally well correlated (R² ≥ 0.76), but slightly over- or underestimated by the hFE models and SED distributions were qualitatively different between hFE and micro-FE models. This study suggests that osseointegrated screw-bone construct stiffness can be predicted accurately using hFE models when compared to micro-FE models and that volume average peri-implant SEDs are well correlated. However, the hFE models are highly sensitive to the choice of trabecular bone material properties. PMUBC-derived isotropic material properties represented the best trade-off between model accuracy and complexity in this study.

The predictive ability of a QCT-FE model of the proximal femoral stiffness under multiple load cases is strongly influenced by experimental uncertainties

Despite significant improvements in terms of the predictive ability of Quantitative Computed Tomography based Finite Element (QCT-FE) models in estimating femoral strength (fracture load and stiffness), no substantial clinical adoption of this method has taken place to date. Narrowing the wide variability of FE results by standardizing the methodology and validation protocols, as well as reducing the uncertainties in the FEA process have been proposed as routes towards improved reliability. The aim of this study was to: First, validate a QCT-FE model of proximal femoral stiffness in multiple stance load cases, and second, using a parametric approach, determine the influence of select experimental and modeling parameters on the predictive ability of our model. Ten fresh frozen human femoral samples were tested in neutral stance, 15° adducted and 15° abducted load cases. Voxel-based linear-elastic QCT-FE models of the samples were generated to predict the models' stiffness values in all load cases. The base FE models were validated against the experimental results using linear regression. Thirty six deviated models were created using the minimum and maximum values of experiment-based "plausible range" for 18 parameters in 4 categories of embedding, loading, material, and segmentation. The predictive ability of the models were compared in terms of the coefficient of determination (R²) of the linear regression between the measured and predicted stiffness values in all load cases. Our model was capable of capturing 90% of the variation in the experimental stiffness of the samples in neutral stance position (R² = 0.9, concordance correlation coefficient (CCC) = 0.93, percent root mean squared error (RMSE%) = 8.4%, slope and intercept not significantly different from unity and zero, respectively). Embedding and loading categories strongly affected the predictive ability of the models with an average percent difference in R² of 4.36% ± 2.77 and 2.96% ± 1.69 for the stance-neutral load case, respectively. The performance of the models were significantly different in adducted and abducted load cases with their R² dropping to 71% and 70%, respectively. Similarly, off-axes load cases were affected by the parameters differently compared to the neutral load case, with the loading parameter category imposing more than 10% difference on their R², larger than all other categories. We also showed that automatically selecting the best performing plausible value for each parameter and each sample would result in a perfectly linear correlation (R²> 0.99) between the "tuned" model's predicted stiffness and experimental results. Based on our results, high sensitivity of the model performance to experimental parameters requires extra diligence in modeling the embedding geometry and the loading angles since these sources of uncertainty could dwarf the effects of material modeling and image processing parameters. The results of this study could help in improving the robustness of the QCT-FE models of proximal femur by limiting the uncertainties in the experimental and modeling steps.

Comparison of linear and nonlinear stepwise μFE displacement predictions to digital volume correlation measurements of trabecular bone biopsies

Digital volume correlation (DVC) enables to evaluate the ability of μFE models in predicting experimental results on the mesoscale. In this study predicted displacement fields of three different linear and materially nonlinear μFE simulation methods were compared to DVC measured displacement fields at specific load steps in the elastic regime (Step_El) and after yield (Step_Ult). Five human trabecular bone biopsies from a previous study were compressed in several displacement steps until failure. At every compression step, μCT images (resolution: 36 μm) were recorded. A global DVC algorithm was applied to compute the displacement fields at all loading steps. The unloaded 3D images were then used to generate homogeneous, isotropic, linear and materially nonlinear μFE models. Three different μFE simulation methods were used: linear (L), nonlinear (NL), and nonlinear stepwise (NLS). Regarding L and NL, the boundary conditions were derived from the interpolated displacement fields at Step_El and Step_Ult, while for the NLS method nonlinear changes of the boundary conditions of the experiments were captured using the DVC displacement field of every available load step until Step_El and Step_Ult. The predicted displacement fields of all μFE simulation methods were in good agreement with the DVC measured displacement fields (individual specimens: R²>0.83 at Step_El and R²>0.59 at Step_Ult; pooled data: R²>0.97 at Step_El and R²>0.92 at Step_Ult). At Step_El, all three simulation methods showed similar intercepts, slopes, and coefficients of determination while the nonlinear μFE models improved the prediction of the displacement fields slightly in all Cartesian directions at Step_Ult (individual specimens: L: R²>0.59 and NL, NLS: R²>0.68; pooled data: L: R²>0.92 and NL, NLS: R²>0.94). Damaged/overstrained elements in L, NL, and NLS occurred at similar locations but the number of overstrained elements was overestimated when using the L simulation method. Considering the increased solving time of the nonlinear μFE models as well as the acceptable performance in displacement prediction of the linear μFE models, one can conclude that for similar use cases linear μFE models represent the best compromise between computational effort and accuracy of the displacement field predictions.

Comparison of different microCT-based morphology assessment tools using human trabecular bone

MicroCT-based morphological parameters are often used to quantify the structural properties of trabecular bone. Various software tools are available for calculating these parameters. Studies that examine the comparability of their results are rare. Four different software tools were used to analyse a set of 701 microCT images from human trabecular bone samples. Bone volume to total volume (BV/TV), bone surface (BS), trabecular thickness (Tb. Th.) and degree of anisotropy (DA) were evaluated. BV/TV shows very low difference (-0.18 ± 0.15%). The difference in BS could be reduced below 5% if artificial cut surfaces are not included. Tb. Th. and Tb. Sp. show differences of maximal -12% although the same theoretical background is used. DA is most critical with differences from 4.75 ± 3.70% (medtool vs. Scanco), over -38.61 ± 13.15% (BoneJ vs. Scanco), up to 80.52 ± 50.04% (medtool vs. BoneJ). Quantitative results should be considered with caution, especially when comparing different studies. Introducing standardization procedures and the disclosure of underlying algorithms and their respective implementations could improve this issue.

Metacarpophalangeal joint loads during bonobo locomotion: model predictions versus proxies

The analysis of internal trabecular and cortical bone has been an informative tool for drawing inferences about behaviour in extant and fossil primate taxa. Within the hand, metacarpal bone architecture has been shown to correlate well with primate locomotion; however, the extent of morphological differences across taxa is unexpectedly small given the variability in hand use. One explanation for this observation is that the activity-related differences in the joint loads acting on the bone are simply smaller than estimated based on commonly used proxies (i.e. external loading and joint posture), which neglect the influence of muscle forces. In this study, experimental data and a musculoskeletal finger model are used to test this hypothesis by comparing differences between climbing and knuckle-walking locomotion of captive bonobos (Pan paniscus) based on (i) joint load magnitude and direction predicted by the models and (ii) proxy estimations. The results showed that the activity-related differences in predicted joint loads are indeed much smaller than the proxies would suggest, with joint load magnitudes being almost identical between the two locomotor modes. Differences in joint load directions were smaller but still evident, indicating that joint load directions might be a more robust indicator of variation in hand use than joint load magnitudes. Overall, this study emphasizes the importance of including muscular forces in the interpretation of skeletal remains and promotes the use of musculoskeletal models for correct functional interpretations.

Musculoskeletal models of a human and bonobo finger: parameter identification and comparison to in vitro experiments

Introduction

Knowledge of internal finger loading during human and non-human primate activities such as tool use or knuckle-walking has become increasingly important to reconstruct the behaviour of fossil hominins based on bone morphology. Musculoskeletal models have proven useful for predicting these internal loads during human activities, but load predictions for non-human primate activities are missing due to a lack of suitable finger models. The main goal of this study was to implement both a human and a representative non-human primate finger model to facilitate comparative studies on metacarpal bone loading. To ensure that the model predictions are sufficiently accurate, the specific goals were: (1) to identify species-specific model parameters based on in vitro measured fingertip forces resulting from single tendon loading and (2) to evaluate the model accuracy of predicted fingertip forces and net metacarpal bone loading in a different loading scenario.

Materials & Methods

Three human and one bonobo (Pan paniscus) fingers were tested in vitro using a previously developed experimental setup. The cadaveric fingers were positioned in four static postures and load was applied by attaching weights to the tendons of the finger muscles. For parameter identification, fingertip forces were measured by loading each tendon individually in each posture. For the evaluation of model accuracy, the extrinsic flexor muscles were loaded simultaneously and both the fingertip force and net metacarpal bone force were measured. The finger models were implemented using custom Python scripts. Initial parameters were taken from literature for the human model and own dissection data for the bonobo model. Optimized model parameters were identified by minimizing the error between predicted and experimentally measured fingertip forces. Fingertip forces and net metacarpal bone loading in the combined loading scenario were predicted using the optimized models and the remaining error with respect to the experimental data was evaluated.

Results

The parameter identification procedure led to minor model adjustments but considerably reduced the error in the predicted fingertip forces (root mean square error reduced from 0.53/0.69 N to 0.11/0.20 N for the human/bonobo model). Both models remained physiologically plausible after the parameter identification. In the combined loading scenario, fingertip and net metacarpal forces were predicted with average directional errors below 6° and magnitude errors below 12%.

Conclusions

This study presents the first attempt to implement both a human and non-human primate finger model for comparative palaeoanthropological studies. The good agreement between predicted and experimental forces involving the action of extrinsic flexors—which are most relevant for forceful grasping—shows that the models are likely sufficiently accurate for comparisons of internal loads occurring during human and non-human primate manual activities.

QCT-based finite element prediction of pathologic fractures in proximal femora with metastatic lesions

Predicting pathologic fractures in femora with metastatic lesions remains a clinical challenge. Currently used guidelines are inaccurate, especially to predict non-impeding fractures. This study evaluated the ability of a nonlinear quantitative computed tomography (QCT)-based homogenized voxel finite element (hvFE) model to predict patient-specific pathologic fractures. The hvFE model was generated highly automated from QCT images of human femora. The femora were previously loaded in a one-legged stance setup in order to assess stiffness, failure load, and fracture location. One femur of each pair was tested in its intact state, while the contralateral femur included a simulated lesion on either the superolateral- or the inferomedial femoral neck. The hvFE model predictions of the stiffness (0.47 < R² < 0.94), failure load (0.77 < R² < 0.98), and exact fracture location (68%) were in good agreement with the experimental data. However, the model underestimated the failure load by a factor of two. The hvFE models predicted significant differences in stiffness and failure load for femora with superolateral- and inferomedial lesions. In contrast, standard clinical guidelines predicted identical fracture risk for both lesion sites. This study showed that the subject-specific QCT-based hvFE model could predict the effect of metastatic lesions on the biomechanical behaviour of the proximal femur with moderate computational time and high level of automation and could support treatment strategy in patients with metastatic bone disease.

Inverse remodelling algorithm identifies habitual manual activities of primates based on metacarpal bone architecture

Previously, a micro-finite element (micro-FE)-based inverse remodelling method was presented in the literature that reconstructs the loading history of a bone based on its architecture alone. Despite promising preliminary results, it remains unclear whether this method is sensitive enough to detect differences of bone loading related to pathologies or habitual activities. The goal of this study was to test the sensitivity of the inverse remodelling method by predicting joint loading histories of metacarpal bones of species with similar anatomy but clearly distinct habitual hand use. Three groups of habitual hand use were defined using the most representative primate species: manipulation (human), suspensory locomotion (orangutan), and knuckle-walking locomotion (bonobo, chimpanzee, gorilla). Nine to ten micro-computed tomography scans of each species ([Formula: see text] in total) were used to create micro-FE models of the metacarpal head region. The most probable joint loading history was predicted by optimally scaling six load cases representing joint postures ranging from [Formula: see text] (extension) to [Formula: see text] (flexion). Predicted mean joint load directions were significantly different between knuckle-walking and non-knuckle-walking groups ([Formula: see text]) and in line with expected primary hand postures. Mean joint load magnitudes tended to be larger in species using their hands for locomotion compared to species using them for manipulation. In conclusion, this study shows that the micro-FE-based inverse remodelling method is sensitive enough to detect differences of joint loading related to habitual manual activities of primates and might, therefore, be useful for palaeoanthropologists to reconstruct the behaviour of extinct species and for biomedical applications such as detecting pathological joint loading.

A novel experimental design for the measurement of metacarpal bone loading and deformation and fingertip force

Background

Musculoskeletal and finite element modelling are often used to predict joint loading and bone strength within the human hand, but there is a lack of in vitro evidence of the force and strain experienced by hand bones.

Methods

This study presents a novel experimental setup that allows the positioning of a cadaveric digit in a variety of postures with the measurement of force and strain experienced by the third metacarpal. The setup allows for the measurement of fingertip force as well. We tested this experimental setup using three cadaveric human third digits in which the flexor tendons were loaded in two tendon pathways: (1) parallel to the metacarpal bone shaft, with bowstringing; (2) a semi-physiological condition in which the tendons were positioned closer to the bone shaft.

Results

There is substantial variation in metacarpal net force, metacarpal strain and fingertip force between the two tendon pathways. The net force acting on the metacarpal bone is oriented palmarly in the parallel tendon condition, causing tension along the dorsum of the metacarpal shaft, while the force increases and is oriented dorsally in the semi-physiological condition, causing compression of the dorsal metacarpal shaft. Fingertip force is also greater in the semi-physiological condition, implying a more efficient grip function. Inter-individual variation is observed in the radioulnar orientation of the force experienced by the metacarpal bone, the fingertip force, and the strain patterns on the metacarpal shaft.

Conclusion

This study demonstrates a new method for measuring force and strain experienced by the metacarpal, and fingertip force in cadaveric digits that can, in turn, inform computation models. Inter-individual variation in loads experienced by the third digit suggest that there are differences in joint contact and/or internal bone structure across individuals that are important to consider in clinical and evolutionary contexts.

Plausibility and parameter sensitivity of micro-finite element-based joint load prediction at the proximal femur

A micro-finite element-based method to estimate the bone loading history based on bone architecture was recently presented in the literature. However, a thorough investigation of the parameter sensitivity and plausibility of this method to predict joint loads is still missing. The goals of this study were (1) to analyse the parameter sensitivity of the joint load predictions at one proximal femur and (2) to assess the plausibility of the results by comparing load predictions of ten proximal femora to in vivo hip joint forces measured with instrumented prostheses (available from www.orthoload.com ). Joint loads were predicted by optimally scaling the magnitude of four unit loads (inclined [Formula: see text] to [Formula: see text] with respect to the vertical axis) applied to micro-finite element models created from high-resolution computed tomography scans ([Formula: see text]m voxel size). Parameter sensitivity analysis was performed by varying a total of nine parameters and showed that predictions of the peak load directions (range 10[Formula: see text]-[Formula: see text]) are more robust than the predicted peak load magnitudes (range 2344.8-4689.5 N). Comparing the results of all ten femora with the in vivo loading data of ten subjects showed that peak loads are plausible both in terms of the load direction (in vivo: [Formula: see text], predicted: [Formula: see text]) and magnitude (in vivo: [Formula: see text], predicted: [Formula: see text]). Overall, this study suggests that micro-finite element-based joint load predictions are both plausible and robust in terms of the predicted peak load direction, but predicted load magnitudes should be interpreted with caution.

Effect of simulated metastatic lesions on the biomechanical behavior of the proximal femur

Pathologic fractures of femora in patients with metastatic cancer are associated with high morbidity and mortality. Prediction of impending fractures is based on unspecific clinical criteria or past clinician's experience, which leads to underestimation or overtreatment. The aim of this study was to investigate the effect of the site of metastatic lesions on biomechanical behavior of the proximal femur. Sixteen pairs of human femora were scanned with quantitative computed tomography (QCT) to asses bone mineral density. One femur of each pair remained intact while a defined lesion was reamed out in either the superolateral or inferomedial portion of the femoral neck of the contralateral femur. All femora were loaded in a mechanical test setup mimicking one-legged stance and stiffness, failure load, and fracture location were determined. In the biomechanical experiments the superolateral lesion and the inferomedial lesion caused a stiffness reduction of 19% and 66%, respectively. The average failure load was 40% and 75% lower for specimens with the superolateral (4.53 ± 1.56 kN) and inferomedial (1.89 ± 1.73 kN) lesions, respectively, compared to intact specimens (7.66 ± 3.34 kN). Lesions in the femoral neck led to reduction in both stiffness and failure load of the proximal femur. Furthermore, the site of the lesion had a large effect on the magnitude of the reduction in biomechanical properties. The presented data emphasize the importance of differentiating between locations of the lesion in pathologic fracture prediction of the metastatic femur and underline the insufficient accuracy of current predictive guidelines. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:2407-2414, 2017.

Analysis of the three-dimensional anatomical variance of the distal radius using 3D shape models

BACKGROUND

Various medical fields rely on detailed anatomical knowledge of the distal radius. Current studies are limited to two-dimensional analysis and biased by varying measurement locations. The aims were to 1) generate 3D shape models of the distal radius and investigate variations in the 3D shape, 2) generate and assess morphometrics in standardized cut planes, and 3) test the model's classification accuracy.

METHODS

The local radiographic database was screened for CT-scans of intact radii. 1) The data sets were segmented and 3D surface models generated. Statistical 3D shape models were computed (overall, gender and side separate) and the 3D shape variation assessed by evaluating the number of modes. 2) Anatomical landmarks were assigned and used to define three standardized cross-sectional cut planes perpendicular to the main axis. Cut planes were generated for the mean shape models and each individual radius. For each cut plane, the following morphometric parameters were calculated and compared: maximum width and depth, perimeter and area. 3) The overall shape model was utilized to evaluate the predictive value (leave one out cross validation) for gender and side identification within the study population.

RESULTS

Eighty-six radii (45 left, 44% female, 40 ± 18 years) were included. 1) Overall, side and gender specific statistical 3D models were successfully generated. The first mode explained 37% of the overall variance. Left radii had a higher shape variance (number of modes: 20 female / 23 male) compared to right radii (number of modes: 6 female / 6 male). 2) Standardized cut planes could be defined using anatomical landmarks. All morphometric parameters decreased from distal to proximal. Male radii were larger than female radii with no significant side difference. 3) The overall shape model had a combined median classification probability for side and gender of 80%.

CONCLUSIONS

Statistical 3D shape models of the distal radius can be generated using clinical CT-data sets. These models can be used to assess overall bone variance, define and analyze standardized cut-planes, and identify the gender of an unknown sample. These data highlight the potential of shape models to assess the 3D anatomy and anatomical variance of human bones.

The effect of the extensor mechanism on maximum isometric fingertip forces: A numerical study on the index finger

The extensor mechanism is a tendinous network connecting intrinsic and extrinsic muscles of the finger and its function has not yet been fully understood. The goal of this study was to assess the effect of the extensor mechanism on the maximum isometric fingertip forces - a parameter which is essential for grasping. For this purpose, maximum fingertip forces in all directions (i.e. feasible force sets) of two musculoskeletal models of the index finger were compared: the wEM model included a full representation of the extensor mechanism, whereas in the noEM model the extensor mechanism was replaced by a single extensor tendon without connectivity to intrinsic muscles. The feasible force sets were computed in the flexion-extension plane for nine postures. Forces in four predefined directions (palmar, proximal, dorsal, and distal), and the peak resultant forces were evaluated. Averaged forces in all four predefined directions were considerably larger in the wEM model (+187.6%). However, peak resultant forces were slightly lower in the wEM model (-4.3% on average). The general advantage of the wEM model could be explained by co-contraction of intrinsic and extrinsic extensor muscles which allowed reaching larger activation levels of the extrinsic flexors. Only within a narrow range of force directions the co-contraction of intrinsic muscles limited the fingertip forces and lead to lower peak resultant forces in the wEM model. Rather than maximizing peak resultant forces, it appears that the extensor mechanism is a sophisticated tool for increasing maximum fingertip forces over a broad range of postures and force directions - making the finger more versatile during grasping.

The influence of distal screw length on the primary stability of volar plate osteosynthesis--a biomechanical study

Background

Extensor tendon irritation is one of the most common complications following volar locking plate osteosynthesis (VLPO) for distal radius fractures. It is most likely caused by distal screws protruding the dorsal cortex. Shorter distal screws could avoid this, yet the influence of distal screw length on the primary stability in VLPO is unknown. The aim of this study was to compare 75 to 100 % distal screw lengths in VLPO.

Methods

A biomechanical study was conducted on 11 paired fresh-frozen radii. HRpQCT scans were performed to assess bone mineral density (BMD) and bone mineral content (BMC). The specimens were randomized pair-wise into two groups: 100 % (group A) and 75 % (group B) unicortical distal screw lengths. A validated fracture model for extra-articular distal radius fractures (AO-23 A3) was used. Polyaxial volar locking plates were mounted, and distal screws was inserted using a drill guide block. For group A, the distal screw tips were intended to be flush or just short of the dorsal cortex. In group B, a target screw length of 75 % was calculated. The specimens were tested to failure using a displacement-controlled axial compression test. Primary biomechanical stability was assessed by stiffness, elastic limit, and maximum force as well as with residual tilt, which quantified plastic deformation.

Results

Nine specimens were tested successfully. BMD and BMC did not differ between the two groups. The mean distal screw length of group A was 21.7 ± 2.6 mm (range: 16 to 26 mm), for group B 16.9 ± 1.9 mm (range: 12 to 20 mm). Distal screws in group B were on average 5.6 ± 0.9 mm (range: 3 to 7 mm) shorter than measured. No significant differences were found for stiffness (706 ± 103 N/mm vs. 660 ± 124 N/mm), elastic limit (177 ± 25 N vs. 167 ± 36 N), maximum force (493 ± 139 N vs. 471 ± 149 N), or residual tilt (7.3° ± 0.7° vs. 7.1° ± 1.3°).

Conclusion

The 75 % distal screw length in VLPO provides similar primary stability to 100 % unicortical screw length. This study, for the first time, provides the biomechanical basis to choose distal screws significantly shorter then measured.

Electronic supplementary material

The online version of this article (doi:10.1186/s13018-015-0283-8) contains supplementary material, which is available to authorized users.