Computational modeling of long-term effects of prophylactic vertebroplasty on bone adaptation

A Density-Dependent Target Stimulus for Inverse Bone (Re)modeling with Homogenized Finite Element Models

Inverse bone (re)modeling (IBR) can infer physiological loading conditions from the bone microstructure. IBR scales unit loads, imposed on finite element (FE) models of a bone, such that the trabecular microstructure is homogeneously loaded and the difference to a target stimulus is minimized. Micro-FE (µFE) analyses are typically used to model the microstructure, but computationally more efficient, homogenized FE (hFE) models, where the microstructure is replaced by an equivalent continuum, could be used instead. However, also the target stimulus has to be translated from the tissue to the continuum level. In this study, a new continuum-level target stimulus relating relative bone density and strain energy density is proposed. It was applied using different types of hFE models to predict the physiological loading of 21 distal radii sections, which was subsequently compared to µFE-based IBR. The hFE models were able to correctly identify the dominant load direction and showed a high correlation of the predicted forces, but mean magnitude errors ranged from − 14.7 to 26.6% even for the best models. While µFE-based IBR can still be regarded as a gold standard, hFE-based IBR enables faster predictions, the usage of more sophisticated boundary conditions, and the usage of clinical images.

A 20-Year Review of Biomechanical Experimental Studies on Spine Implants Used for Percutaneous Surgical Repair of Vertebral Compression Fractures

A vertebral compression fracture (VCF) is an injury to a vertebra of the spine affecting the cortical walls and/or middle cancellous section. The most common risk factor for a VCF is osteoporosis, thus predisposing the elderly and postmenopausal women to this injury. Clinical consequences include loss of vertebral height, kyphotic deformity, altered stance, back pain, reduced mobility, reduced abdominal space, and reduced thoracic space, as well as early mortality. To restore vertebral mechanical stability, overall spine function, and patient quality of life, the original percutaneous surgical intervention has been vertebroplasty, whereby bone cement is injected into the affected vertebra. Because vertebroplasty cannot fully restore vertebral height, newer surgical techniques have been developed, such as kyphoplasty, stents, jacks, coils, and cubes. But, relatively few studies have experimentally assessed the biomechanical performance of these newer procedures. This article reviews over 20 years of scientific literature that has experimentally evaluated the biomechanics of percutaneous VCF repair methods. Specifically, this article describes the basic operating principles of the repair methods, the study protocols used to experimentally assess their biomechanical performance, and the actual biomechanical data measured, as well as giving a number of recommendations for future research directions.

Bone Mechanoregulation Allows Subject-Specific Load Estimation Based on Time-Lapsed Micro-CT and HR-pQCT in Vivo

Patients at high risk of fracture due to metabolic diseases frequently undergo long-term antiresorptive therapy. However, in some patients, treatment is unsuccessful in preventing fractures or causes severe adverse health outcomes. Understanding load-driven bone remodelling, i.e., mechanoregulation, is critical to understand which patients are at risk for progressive bone degeneration and may enable better patient selection or adaptive therapeutic intervention strategies. Bone microarchitecture assessment using high-resolution peripheral quantitative computed tomography (HR-pQCT) combined with computed mechanical loads has successfully been used to investigate bone mechanoregulation at the trabecular level. To obtain the required mechanical loads that induce local variances in mechanical strain and cause bone remodelling, estimation of physiological loading is essential. Current models homogenise strain patterns throughout the bone to estimate load distribution in vivo, assuming that the bone structure is in biomechanical homoeostasis. Yet, this assumption may be flawed for investigating alterations in bone mechanoregulation. By further utilising available spatiotemporal information of time-lapsed bone imaging studies, we developed a mechanoregulation-based load estimation (MR) algorithm. MR calculates organ-scale loads by scaling and superimposing a set of predefined independent unit loads to optimise measured bone formation in high-, quiescence in medium-, and resorption in low-strain regions. We benchmarked our algorithm against a previously published load history (LH) algorithm using synthetic data, micro-CT images of murine vertebrae under defined experimental in vivo loadings, and HR-pQCT images from seven patients. Our algorithm consistently outperformed LH in all three datasets. In silico-generated time evolutions of distal radius geometries (n = 5) indicated significantly higher sensitivity, specificity, and accuracy for MR than LH (p < 0.01). This increased performance led to substantially better discrimination between physiological and extra-physiological loading in mice (n = 8). Moreover, a significantly (p < 0.01) higher association between remodelling events and computed local mechanical signals was found using MR [correct classification rate (CCR) = 0.42] than LH (CCR = 0.38) to estimate human distal radius loading. Future applications of MR may enable clinicians to link subtle changes in bone strength to changes in day-to-day loading, identifying weak spots in the bone microstructure for local intervention and personalised treatment approaches.

How the clinical dosage of bone cement biomechanically affects adjacent vertebrae

Objective

This study evaluated the biomechanical changes in the adjacent vertebrae under a physiological load (500 N) when the clinically relevant amount of bone cement was injected into fractured cadaver vertebral bodies.

Methods

The embalmed cadaver thoracolumbar specimens in which each vertebral body (T12–L2) had a BMD of < 0.75 g/cm² were used for the experiment. For establishing a fracture model, the upper one third of the L1 vertebra was performed wedge osteotomy and the superior endplate was kept complete. Stiffness of specimens was measured in different states. Strain of the adjacent vertebral body and intervertebral disc were measured in pre-fracture, post-fracture, and after augmentation by non-contact optical strain measurement system.

Results

The average amount of bone cement was 4.4 ml (3.8–5.0 ml). The stiffness of after augmentation was significantly higher than the stiffness of post-fracture (p < 0.05), but still lower than pre-fracture stiffness (p < 0.05). After augmentation, the adjacent upper vertebral strain showed no significant difference (p > 0.05) with pre-fracture, while the strain of adjacent lower vertebral body was significantly higher than that before fracture (p < 0.05). In flexion, T12/L1 intervertebral disc strain was significantly greater after augmentation than after the fracture (p < 0.05), but there was no significant difference from that before the fracture (p > 0.05); L1/2 vertebral strain after augmentation was significantly less than that after the fracture (p < 0.05), but there was no significant difference from that before the fracture (p > 0.05).

Conclusions

PVP may therefore have partially reversed the abnormal strain state of adjacent vertebral bodies which was caused by fracture.

Effect of size and location of simulated lytic lesions on the structural properties of human vertebral bodies, a micro-finite element study

Currently, the Spinal Instability Neoplastic Score system is used in clinics to evaluate the risk of fracture in patients with spinal metastases. This method, however, does not always provide a clear guideline due to the complexity in accounting for the effect of metastatic lesions on vertebral stability. The aim of this study was to use a validated micro Finite Element (microFE) modelling approach to analyse the effect of the size and location of lytic metastases on the mechanical properties of human vertebral bodies. Micro Computed Tomography based microFE models were generated with and without lytic lesions simulated as holes within a human vertebral body. Single and multiple lytic lesions were simulated with four different sizes and in five different locations. Bone was assumed homogenous, isotropic and linear elastic, and each vertebra was loaded in axial compression. It was observed that the size of lytic lesions was linearly related with the reduction in structural properties of the vertebral body (reduction of stiffness between 3% and 30% for lesion volume between 4% and 35%). The location of lytic lesions did not show a clear effect on predicted structural properties. Single or multiple lesions with the same volume provided similar results. Locally, there was a homogeneous distribution of axial principal strains among the models with and without lytic lesions. This study highlights the potential of microFE models to study the effect of lesions on the mechanical properties of the human vertebral body.

•

MicroFE models can show the effect of lytic lesions on vertebral properties.

•

The size of the lesions was more critical than the location of the lesions.

•

Lesions affecting the cortical shell had a larger effect on the local strains.

•

Multiple lesions showed a similar effect to single lesions.

Numerical investigation of the effect of bone cement porosity on osteoporotic femoral augmentation

Femoroplasty is the injection of bone cement into the proximal femur, enhances the bone load capacity, and is typically applied to osteoporotic femora. To minimize the required injected volume of bone cement and maximize the load capacity enhancement, an optimization problem must be solved, where the modulus of elasticity of the augmented bone is a key element. This paper, through the numerical investigation of a fall on the greater trochanter of an osteoporotic femur, compares different ways to calculate this modulus and introduces an approach, based on the concept of bone cement porosity, which provides results statistically similar to those obtained with other considerations. Based on this approach, the present paper quantifies the correlation between degree of osteoporosis and optimum volume of bone cement. It concludes with an exhaustive search that reveals the effect of the bone cement porosity on the optimum volume of PMMA, for various combinations of the frontal and transverse angles of the fall on the greater trochanter.

Biomechanical effects of different vertebral heights after augmentation of osteoporotic vertebral compression fracture: a three-dimensional finite element analysis

Background

Clinical results have shown that different vertebral heights have been restored post-augmentation of osteoporotic vertebral compression fractures (OVCFs) and the treatment results are consistent. However, no significant results regarding biomechanical effects post-augmentation have been found with different types of vertebral deformity or vertebral heights by biomechanical analysis. Therefore, the present study aimed to investigate the biomechanical effects between different vertebral heights of OVCFs before and after augmentation using three-dimensional finite element analysis.

Methods

Four patients with OVCFs of T12 underwent computed tomography (CT) of the T11-L1 levels. The CT images were reconstructed as simulated three-dimensional finite-element models of the T11-L1 levels (before and after the T12 vertebra was augmented with cement). Four different kinds of vertebral height models included Genant semi-quantitative grades 0, 1, 2, and 3, which simulated unilateral augmentation. These models were assumed to represent vertical compression and flexion, left flexion, and right flexion loads, and the von Mises stresses of the T12 vertebral body were assessed under different vertebral heights before and after bone cement augmentation.

Results

Data showed that the von Mises stresses significantly increased under four loads of OVCFs of the T12 vertebral body before the operation from grade 0 to grade 3 vertebral heights. The maximum stress of grade 3 vertebral height pre-augmentation was produced at approximately 200%, and at more than 200% for grade 0. The von Mises stresses were significantly different between different vertebral heights preoperatively. The von Mises stresses of the T12 vertebral body significantly decreased in four different loads and at different vertebral body heights (grades 0–3) after augmentation. There was no significant difference between the von Mises stresses of grade 0, 1, and 3 vertebral heights postoperatively. The von Mises stress significantly decreased between pre-augmentation and post-augmentation in T12 OVCF models of grade 0–3 vertebral heights.

Conclusion

Vertebral augmentation can sufficiently reduce von Mises stresses at different heights of OVCFs of the vertebral body, although this technique does not completely restore vertebral height to the anatomical criteria.