Biomechanical assessment of vertebrae with lytic metastases with subject-specific finite element models

Simulation of cancer progression in bone by a virtual thermal flux with a case study on lumbar vertebrae with multiple myeloma

BACKGROUND

Two main problems examining the mechanism of cancer progression in the tissues using the computational models are lack of enough knowledge on the effective factors for such events in vivo environments and lack of specific parameters in the available computational models to simulate such complicated reactions.

METHODS

In this study, it was tried to simulate the progression of cancerous lesions in the bone tissues by an independent parameter from the anatomical and physiological characteristics of the tissues, so to degrade the orthotropic mechanical properties of the bone tissues, a virtual temperature was determined to be used by a well-known framework for simulation of damages in the composite materials. First, the reliability of the FE model to simulate hyperelastic response in the intervertebral discs (IVDs) and progressive failure in the bony components were verified by simulation of some In-Vitro tests, available in the literature. Then, the progression of the osteolytic damage was simulated in a clinical case with multiple myeloma in the lumbar vertebrae.

RESULTS

The FE model could simulate stress-shielding and diffusion of lesion in the posterior elements of the damaged vertebra which led to spinal stenosis. The load carrying shares associated with the anterior half and the posterior half of the examined vertebral body and the posterior elements were estimated equal to 41 %, 47 % and 12 %, respectively for the intact condition, that changed to 14 %, 16 % and 70 %, when lesion occupied one third of the vertebral body.

CONCLUSION

Correlation of the FE results with the deformation shapes, observed in the MRIs for the clinical case study, indicated appropriateness of the procedure, proposed for simulation of the progressive osteolytic damage in the vertebral segments. The future studies may follow simulation of tumor growth for various metastatic tissues using the method, established here.

The effect of deep learning-based lesion segmentation on failure load calculations of metastatic femurs using finite element analysis

Bone ranks as the third most frequent tissue affected by cancer metastases, following the lung and liver. Bone metastases are often painful and may result in pathological fracture, which is a major cause of morbidity and mortality in cancer patients. To quantify fracture risk, finite element (FE) analysis has shown to be a promising tool, but metastatic lesions are typically not specifically segmented and therefore their mechanical properties may not be represented adequately. Deep learning methods potentially provide the opportunity to automatically segment these lesions and change the mechanical properties more adequately. In this study, our primary focus was to gain insight into the performance of an automatic segmentation algorithm for femoral metastatic lesions using deep learning methods and the subsequent effects on FE outcomes. The aims were to determine the similarity between manual segmentation and automatic segmentation; the differences in predicted failure load between FE models with automatically segmented osteolytic and mixed lesions and the models with CT-based lesion values (the gold standard); and the effect on the BOne Strength (BOS) score (failure load adjusted for body weight) and subsequent fracture risk assessments. From two patient cohorts, a total number of 50 femurs with osteolytic and mixed metastatic lesions were included in this study. The femurs were segmented from CT images and transferred into FE meshes. The material behavior was implemented as non-linear isotropic. These FE models were considered as gold standard (Finite Element no Segmented Lesion: FE-no-SL), whereby the local calcium equivalent density of both femur and metastatic lesion was extracted from CT-values. Lesions in the femur were manually segmented by two biomechanical experts after which final lesion segmentation for each femur was obtained based on consensus of opinions between two observers. Subsequently, a self-configuring variant of the popular deep learning model U-Net known as nnU-Net was used to automatically segment metastatic lesions within the femur. For these models with segmented lesions (Finite Element with Segmented Lesion: FE-with-SL), the calcium equivalent density within the metastatic lesions was set to zero after being segmented by the neural network, simulating absence of load-bearing capacity of these lesions. The models (either with or without automatically segmented lesions) were loaded incrementally in axial direction until failure was simulated. Dice coefficient was used to evaluate the similarity of the manual and automatic segmentation. Mean calcium equivalent density values within the automatically segmented lesions were calculated. Failure loads and patterns were determined. Furthermore, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were calculated for both groups by comparing the predictions to the occurrence or absence of actual fracture within the patient cohorts. The automatic segmentation algorithm performed in a none-robust manner. Dice coefficients describing the similarity between consented manual and automatic segmentations were relatively low (mean 0.45 ± standard deviation 0.33, median 0.54). Failure load difference between the FE-no-SL and FE-with-SL groups varied from 0 % to 48 % (mean 6.6 %). Correlation analysis of failure loads between the two groups showed a strong relationship (R2 > 0.9). From the 50 cases, four cases showed clear deviations for which models with automatic lesion segmentation (FE-with-SL) showed considerably lower failure loads. In the whole database including osteolytic and mixed lesions, sensitivity and NPV remained the same, but specificity and PPV decreased from 94 % to 83 %, and from 78 % to 54 % respectively from FE-no-SL to FE-with-SL. This study indicates that the nnU-Net yielded none-robust outcomes in femoral lesion segmentation and that other segmentation algorithms should be considered. However, the difference in failure pattern and failure load between FE models with automatically segmented osteolytic and mixed lesions were relatively small in most cases with a few exceptions. On the other hand, the accuracy of fracture risk assessment using the BOS score was lower compared to the FE-no-SL. In conclusion, this study showed that automatic lesion segmentation is a none-solved issue and therefore, quantifying lesion characteristics and the subsequent effect on the fracture risk using deep learning will remain challenging.

The effect of periprosthetic bone loss on the failure risk of tibial total knee arthroplasty

The effect of long-term periprosthetic bone loss on the process of aseptic loosening of tibial total knee arthroplasty (TKA) is subject to debate. Contradicting studies can be found in literature, reporting either bone resorption or bone formation before failure of the tibial tray. The aim of the current study was to investigate the effects of bone resorption on failure of tibial TKA, by simulating clinical postoperative bone density changes in finite element analysis (FEA) models and FEA models were created of two tibiae representing cases with good and poor initial bone quality which were subjected to a walking configuration and subsequently to a traumatic stumbling load. Bone failure was simulated using a crushable foam model incorporating progressive yielding. Repetitive loading under a level walking load did not result in failure of the periprosthetic bone in neither the good nor poor bone quality tibia at the baseline bone densities. When applying a stumble load, a collapse of the tibial reconstruction was noticed in the poor bone quality model. Incorporating postoperative bone loss led to a significant increase of the failure risk, particularly for the poor bone quality model in which subsidence of the tibial component was substantial. Our results suggest bone loss can lead to an increased risk of a collapse of the tibial component, particularly in case of poor bone quality at the time of surgery. The study also examined the probability of medial or lateral subsidence of the implant and aimed to improve clinical implications. The FEA model simulated plastic deformation of the bone and implant subsidence, with further validation required via mechanical experiments.

Patient-specific bone material modelling can improve the predicted biomechanical outcomes of sacral fracture fixation techniques: A comparative finite element study

OBJECTIVE

To evaluate and compare the biomechanical efficacy of six iliosacral screw fixation techniques for treating unilateral AO Type B2 (Denis Type II) sacral fractures using literature-based and QCT-based bone material properties in finite element (FE) models.

METHODS

Two FE models of the intact pelvis were constructed: the literature-based model (LBM) with bone material properties taken from the literature, and the patient-specific model (PSM) with QCT-derived bone material properties. Unilateral transforaminal sacral fracture was modelled to assess different fixation techniques: iliosacral screw (ISS) at the first sacral vertebra (S1) (ISS1), ISS at the second sacral vertebra (S2) (ISS2), ISS at S1 and S2 (ISS12), transverse iliosacral screws (TISS) at S1 (TISS1), TISS at S2 (TISS2), and TISS at S1 and S2 (TISS12). A 600 N vertical load with both acetabula fixed was applied. Vertical stiffness (VS), relative interfragmentary displacement (RID), and the von Mises stress values in the screws and fracture interface were analysed.

RESULTS

The lowest and highest normalised VS was given by ISS1 and TISS12 techniques for LBM and PSM, with 137 % and 149 %, and 375 % and 472 %, respectively. In comparison with the LBM, the patient-specific bone modelling increased the maximum screw stress values by 19.3, 16.3, 27.8, 2.3, 24.4 and 7.8 % for ISS1, ISS2, ISS12, TISS1, TISS2 and TISS12, respectively. The maximum RID values were between 0.10 mm and 0.47 mm for all fixation techniques in both models. The maximum von Mises stress results on the fracture interface show a substantial difference between the two models, as PSM (mean ± SD of 15.76 ± 8.26 MPa) gave lower stress values for all fixation techniques than LBM (mean ± SD of 28.95 ± 6.91 MPa).

CONCLUSION

The differences in stress distribution underline the importance of considering locally defined bone material properties when investigating internal mechanical parameters. Based on the results, all techniques demonstrated clinically sufficient stability, with TISS12 being superior from a biomechanical standpoint. Both LBM and PSM models indicated a consistent trend in ranking the fixation techniques based on stability. However, long-term clinical trials are recommended to confirm the findings of the study.

The role of bone metastases on the mechanical competence of human vertebrae

Spine is the most common site for bone metastases. The evaluation of the mechanical competence and failure location in metastatic vertebrae is a biomechanical and clinical challenge. Little is known about the failure behaviour of vertebrae with metastatic lesions. The aim of this study was to use combined micro-Computed Tomography (microCT) and time-lapsed mechanical testing to reveal the failure location in metastatic vertebrae. Fifteen spine segments, each including a metastatic and a radiologically healthy vertebra, were tested in compression up to failure within a microCT. Volumetric strains were measured using Digital Volume Correlation. The images of undeformed and deformed specimens were overlapped to identify the failure location. Vertebrae with lytic metastases experienced the largest average compressive strains (median ± standard deviation: -8506 ± 4748microstrain), followed by the vertebrae with mixed metastases (-7035 ± 15605microstrain), the radiologically healthy vertebrae (-5743 ± 5697microstrain), and the vertebrae with blastic metastases (-3150 ± 4641microstrain). Strain peaks were localised within and nearby the lytic lesions or around the blastic tissue. Failure between the endplate and the metastasis was identified in vertebrae with lytic metastases, whereas failure localised around the metastasis in vertebrae with blastic lesions. This study showed for the first time the role of metastases on the vertebral internal deformations. While lytic lesions lead to failure of the metastatic vertebra, vertebrae with blastic metastases are more likely to induce failure in the adjacent vertebrae. Nevertheless, every metastatic lesion affects the vertebral deformation differently, making it essential to assess how the lesion affects the bone microstructure. These results suggest that the properties of the lesion (type, size, location within the vertebral body) should be considered when developing clinical tools to predict the risk of fracture in patients with metastatic lesions.

Mechanoregulation may drive osteolysis during bone metastasis: A finite element analysis of the mechanical environment within bone tissue during bone metastasis and osteolytic resorption

Metastatic bone disease occurs in 70-80% of advanced breast cancer patients and bone tissue is accepted to have attractive physical properties that facilitate cancer cell attraction, adhesion, and invasion. Bone cells also facilitate tumour invasion by biochemical signalling and through resorption of the bone matrix (osteolysis), which releases factors that further stimulate tumour cell activity. The evolving mechanical environment during tumour invasion might play an important role in these processes, as the activity of both bone and cancer cells is regulated by mechanical cues. In particular bone loss and altered mineralisation have been reported, yet how these alter the mechanical environment local to bone and tumour cells is unknown. The objective of this study is to quantify changes in the mechanical environment within bone tissue, during bone metastasis and osteolytic resorption, using finite element analysis (FEA) models reconstructed from high-resolution μCT images of metastatic mouse bone. In particular, we quantify time-dependent changes in mechanical stimuli, local to and distant from an invading tumour mass, to investigate putative mechanobiological cues for osteolysis during bone metastasis. We report here that in early metastasis (3 weeks after tumour inoculation), there was a decrease in strain distribution within the proximal femur trabecular and distal cortical bone tissue. These changes in the mechanical environment preceded extensive osteolytic destruction, but coincided with the onset of early osteolysis, cortical thickening and mineralisation of proximal and distal femur bone. We propose that early changes in the mechanical environment within bone tissue may activate resorption by osteoclast cells and thereby contribute to the extensive osteolytic bone loss at later stage (6 weeks) bone metastasis.

Evaluation of Load-To-Strength Ratios in Metastatic Vertebrae and Comparison With Age- and Sex-Matched Healthy Individuals

Vertebrae containing osteolytic and osteosclerotic bone metastases undergo pathologic vertebral fracture (PVF) when the lesioned vertebrae fail to carry daily loads. We hypothesize that task-specific spinal loading patterns amplify the risk of PVF, with a higher degree of risk in osteolytic than in osteosclerotic vertebrae. To test this hypothesis, we obtained clinical CT images of 11 cadaveric spines with bone metastases, estimated the individual vertebral strength from the CT data, and created spine-specific musculoskeletal models from the CT data. We established a musculoskeletal model for each spine to compute vertebral loading for natural standing, natural standing + weights, forward flexion + weights, and lateral bending + weights and derived the individual vertebral load-to-strength ratio (LSR). For each activity, we compared the metastatic spines' predicted LSRs with the normative LSRs generated from a population-based sample of 250 men and women of comparable ages. Bone metastases classification significantly affected the CT-estimated vertebral strength (Kruskal-Wallis, p < 0.0001). Post-test analysis showed that the estimated vertebral strength of osteosclerotic and mixed metastases vertebrae was significantly higher than that of osteolytic vertebrae (p = 0.0016 and p = 0.0003) or vertebrae without radiographic evidence of bone metastasis (p = 0.0010 and p = 0.0003). Compared with the median (50%) LSRs of the normative dataset, osteolytic vertebrae had higher median (50%) LSRs under natural standing (p = 0.0375), natural standing + weights (p = 0.0118), and lateral bending + weights (p = 0.0111). Surprisingly, vertebrae showing minimal radiographic evidence of bone metastasis presented significantly higher median (50%) LSRs under natural standing (p < 0.0001) and lateral bending + weights (p = 0.0009) than the normative dataset. Osteosclerotic vertebrae had lower median (50%) LSRs under natural standing (p < 0.0001), natural standing + weights (p = 0.0005), forward flexion + weights (p < 0.0001), and lateral bending + weights (p = 0.0002), a trend shared by vertebrae with mixed lesions. This study is the first to apply musculoskeletal modeling to estimate individual vertebral loading in pathologic spines and highlights the role of task-specific loading in augmenting PVF risk associated with specific bone metastatic types. Our finding of high LSRs in vertebrae without radiologically observed bone metastasis highlights that patients with metastatic spine disease could be at an increased risk of vertebral fractures even at levels where lesions have not been identified radiologically.

Biomechanical analysis of vertebral wedge deformity in elderly women with quantitative CT-based finite element analysis

BACKGROUND

To explore the vertebral deformity angle (VD angle) of 1^st lumbar vertebral body (L1) in elderly women, investigate the influence of VD on vertebral stiffness (VS) by biomechanical analysis using quantitative computed tomography-based finite element analysis (QCT-FEA).

METHODS

Two hundred seventy eight participants were recruited, and underwent QCT scan. Measured VD angles of L1, and constructed QCT-FEA models of L1 with the minimum (0.59°), median (5.79°) and maximum (11.15°) VD angles, respectively. Loads in two directions were applied on the upper edge of L1 with a force of 700 N, and vertebral stiffness (VS) was defined as the ratio of 700 N and displacement at the superior reference point: (1) perpendicular to the upper edge of L1 (defined as VS-U); (2) perpendicular to the lower edge of L1(defined as VS-L).

RESULTS

Age was very weak positively correlated with VD angle, moderate negatively correlated with vBMD, and moderate negatively correlated with VS (P < 0.05). VS-U was significantly different among three VD angles, so was VS-L (P < 0.001). VS-U was higher than VS-L in 5.79° and 11.15° VD angles (P < 0.05), however no difference in 0.59° VD angles (P > 0.10).

CONCLUSIONS

VD angle of L1 was slightly increased with age and not correlated with vBMD, and VS was moderate negatively correlated with age, showing that the vertebral body was more likely to fracture with aging. VS-U and VS-L were gradually decreased with the increase of VD angle, and VS-L was lower than VS-U with the increase of VD angle, which showed that vertebral body was more prone to fracture when the load was perpendicular to the lower edge of the vertebral body as the VD angle increasing.

Osteolytic vs. Osteoblastic Metastatic Lesion: Computational Modeling of the Mechanical Behavior in the Human Vertebra after Screws Fixation Procedure

Metastatic lesions compromise the mechanical integrity of vertebrae, increasing the fracture risk. Screw fixation is usually performed to guarantee spinal stability and prevent dramatic fracture events. Accordingly, predicting the overall mechanical response in such conditions is critical to planning and optimizing surgical treatment. This work proposes an image-based finite element computational approach describing the mechanical behavior of a patient-specific instrumented metastatic vertebra by assessing the effect of lesion size, location, type, and shape on the fracture load and fracture patterns under physiological loading conditions. A specific constitutive model for metastasis is integrated to account for the effect of the diseased tissue on the bone material properties. Computational results demonstrate that size, location, and type of metastasis significantly affect the overall vertebral mechanical response and suggest a better way to account for these parameters in estimating the fracture risk. Combining multiple osteolytic lesions to account for the irregular shape of the overall metastatic tissue does not significantly affect the vertebra fracture load. In addition, the combination of loading mode and metastasis type is shown for the first time as a critical modeling parameter in determining fracture risk. The proposed computational approach moves toward defining a clinically integrated tool to improve the management of metastatic vertebrae and quantitatively evaluate fracture risk.

A Three-Dimensional Cement Quantification Method for Decision Prediction of Vertebral Recompression after Vertebroplasty

Objective

Proposing parameters to quantify cement distribution and increasing accuracy for decision prediction of vertebroplasty postoperative complication.

Methods

Finite element analysis was used to biomechanically assess vertebral mechanics (n = 51) after percutaneous vertebroplasty (PVP) or kyphoplasty (PKP). The vertebral space was divided into 27 portions. The numbers of cement occupied portions and numbers of cement-endplate contact portions were defined as overall distribution number (oDN) and overall endplate contact number (oEP), respectively. And cement distribution was parametrized by oDN and oEP. The determination coefficients of vertebral mechanics and parameters (R 2) can validate the correlation of proposed parameters with vertebral mechanics.

Results

oDN and oEP were mainly correlated with failure load (R 2 = 0.729) and stiffness (R 2 = 0.684), respectively. oDN, oEP, failure load, and stiffness had obvious difference between the PVP group and the PKP group (P < 0.05). The regional endplate contact number in the front column is most correlated with vertebral stiffness (R 2 = 0.59) among all regional parameters. Cement volume and volume fraction are not dominant factors of vertebral augmentation, and they are not suitable for postoperative fracture risk prediction.

Conclusions

Proposed parameters with high correlation on vertebral mechanics are promising for clinical utility. The oDN and oEP can strongly affect augmented vertebral mechanics thus is suitable for postoperative fracture risk prediction. The parameters are beneficial for decision-making process of revision surgery necessity. Parametrized methods are also favorable for surgeon's preoperative planning. The methods can be inspirational for clinical image recognition development and auxiliary diagnosis.

Type, size, and position of metastatic lesions explain the deformation of the vertebrae under complex loading conditions

BACKGROUND

Bone metastases may lead to spine instability and increase the risk of fracture. Scoring systems are available to assess critical metastases, but they lack specificity, and provide uncertain indications over a wide range, where most cases fall. The aim of this work was to use a novel biomechanical approach to evaluate the effect of lesion type, size, and location on the deformation of the metastatic vertebra.

METHOD

Vertebrae with metastases were identified from 16 human spines from a donation programme. The size and position of the metastases, and the Spine Instability Neoplastic Score (SINS) were evaluated from clinical Quantitative Computed Tomography images. Thirty-five spine segments consisting of metastatic vertebrae and adjacent healthy controls were biomechanically tested in four different loading conditions. The strain distribution over the entire vertebral bodies was measured with Digital Image Correlation. Correlations between the features of the metastasis (type, size, position and SINS) and the deformation of the metastatic vertebrae were statistically explored.

RESULTS

The metastatic type (lytic, blastic, mixed) characterizes the vertebral behaviour (Kruskal-Wallis, p = 0.04). In fact, the lytic metastases showed more critical deformation compared to the control vertebrae (average: 2-fold increase, with peaks of 14-fold increase). By contrast, the vertebrae with mixed or blastic metastases did not show a clear trend, with deformations similar or lower than the controls. Once the position of the lytic lesion with respect to the loading direction was taken into account, the size of the lesion was significantly correlated with the perturbation to the strain distribution (r² = 0.72, p < 0.001). Conversely, the SINS poorly correlated with the mechanical evidence, and only in case of lytic lesions (r² = 0.25, p < 0.0001).

CONCLUSION

These results highlight the relevance of the size and location of the lytic lesion, which are marginally considered in the current clinical scoring systems, in driving the spinal biomechanical instability. The strong correlation with the biomechanical evidence indicates that these parameters are representative of the mechanical competence of the vertebra. The improved explanatory power compared to the SINS suggests including them in future guidelines for the clinical practice.

Development of a crushable foam model for human trabecular bone

Finite element (FE) simulations can be used to evaluate the mechanical behavior of human bone and allow for quantitative prediction of press-fit implant fixation. An adequate material model that captures post-yield behavior is essential for a realistic simulation. The crushable foam (CF) model is a constitutive model that has recently been proposed in this regard. Compression tests under uniaxial and confined loading conditions were performed on 59 human trabecular bone specimens. Three essential material parameters were obtained as a function of bone mineral density (BMD) to develop the isotropic CF model. The related constitutive rule was implemented in FE models and the results were compared to the experimental data. The CF model provided an accurate simulation of uniaxial compression tests and the post-yield behavior of the stress-strain was well-matched with the experimental results. The model was able to reproduce the confined response of the bone up to 15% of strain. This model allows for simulation of the mechanical behavior of the cellular structure of human bone and adequately predicts the post-yield response of trabecular bone, particularly under uniaxial loading conditions. The model can be further improved to simulate bone collapse due to local overload around orthopaedic implants.

A novel approach to evaluate the effects of artificial bone focal lesion on the three-dimensional strain distributions within the vertebral body

The spine is the first site for incidence of bone metastasis. Thus, the vertebrae have a high potential risk of being weakened by metastatic tissues. The evaluation of strength of the bone affected by the presence of metastases is fundamental to assess the fracture risk. This work proposes a robust method to evaluate the variations of strain distributions due to artificial lesions within the vertebral body, based on in situ mechanical testing and digital volume correlation. Five porcine vertebrae were tested in compression up to 6500N inside a micro computed tomography scanner. For each specimen, images were acquired before and after the application of the load, before and after the introduction of the artificial lesions. Principal strains were computed within the bone by means of digital volume correlation (DVC). All intact specimens showed a consistent strain distribution, with peak minimum principal strain in the range -1.8% to -0.7% in the middle of the vertebra, demonstrating the robustness of the method. Similar distributions of strains were found for the intact vertebrae in the different regions. The artificial lesion generally doubled the strain in the middle portion of the specimen, probably due to stress concentrations close to the defect. In conclusion, a robust method to evaluate the redistribution of the strain due to artificial lesions within the vertebral body was developed and will be used in the future to improve current clinical assessment of fracture risk in metastatic spines.

Patient-Specific Finite Element Models of Posterior Pedicle Screw Fixation: Effect of Screw's Size and Geometry

Pedicle screw fixation is extensively performed to treat spine injuries or diseases and it is common for thoracolumbar fractures. Post-operative complications may arise from this surgery leading to back pain or revisions. Finite element (FE) models could be used to predict the outcomes of surgeries but should be verified when both simplified and realistic designs of screws are used. The aim of this study was to generate patient-specific Computed Tomography (CT)-based FE models of human vertebrae with two pedicle screws, verify the models, and use them to evaluate the effect of the screws' size and geometry on the mechanical properties of the screws-vertebra structure. FE models of the lumbar vertebra implanted with two pedicle screws were created from anonymized CT-scans of three patients. Compressive loads were applied to the head of the screws. The mesh size was optimized for realistic and simplified geometry of the screws with a mesh refinement study. Finally, the optimal mesh size was used to evaluate the sensitivity of the model to changes in screw's size (diameter and length) and geometry (realistic or simplified). For both simplified and realistic models, element sizes of 0.6 mm in the screw and 1.0 mm in the bone allowed to obtain relative differences of approximately 5% or lower. Changes in screw's length resulted in 4-10% differences in maximum deflection, 1-6% differences in peak stress in the screws, 10-22% differences in mean strain in the bone around the screw; changes in screw's diameter resulted in 28-36% differences in maximum deflection, 6-27% differences in peak stress in the screws, and 30-47% differences in mean strain in the bone around the screw. The maximum deflection predicted with realistic or simplified screws correlated very well (R ² = 0.99). The peak stress in screws with realistic or simplified design correlated well (R ² = 0.82) but simplified models underestimated the peak stress. In conclusion, the results showed that the diameter of the screw has a major role on the mechanics of the screw-vertebral structure for each patient. Simplified screws can be used to estimate the mechanical properties of the implanted vertebrae, but the systematic underestimation of the peak stress should be considered when interpreting the results from the FE analyses.

Finite element analysis of wedge and biconcave deformity in four different height restoration after augmentation of osteoporotic vertebral compression fractures

Purpose

Biomechanical comparison of wedge and biconcave deformity of different height restoration after augmentation of osteoporotic vertebral compression fractures was analyzed by three-dimensional finite element analysis (FEA).

Methods

Three-dimensional finite element model (FEM) of T11-L2 segment was constructed from CT scan of elderly osteoporosis patient. The von Mises stresses of vertebrae, intervertebral disc, facet joints, displacement, and range of motion (ROM) of wedge and biconcave deformity were compared at four different heights (Genant 0–3 grade) after T12 vertebral augmentation.

Results

In wedge deformity, the stress of T12 decreased as the vertebral height in neutral position, flexion, extension, and left axial rotation, whereas increased sharply in bending at Genant 0; L1 and L2 decreased in all positions excluding flexion of L2, and T11 increased in neutral position, flexion, extension, and right axial rotation at Genant 0. No significant changes in biconcave deformity. The stress of T11-T12, T12-L1, and L1-L2 intervertebral disc gradually increased or decreased under other positions in wedge fracture, whereas L1-L2 no significant change in biconcave fracture. The utmost overall facet joint stress is at Genant 3, whereas there is no significant change under the same position in biconcave fracture. The displacement and ROM of the wedge fracture had ups and downs, while a decline in all positions excluding extension in biconcave fracture.

Conclusions

The vertebral restoration height after augmentation to Genant 0 affects the von Mises stress, displacement, and ROM in wedge deformity, which may increase the risk of fracture, whereas restored or not in biconcave deformity.

Computational Challenges in Tissue Engineering for the Spine

This paper deals with a brief review of the recent developments in computational modelling applied to innovative treatments of spine diseases. Additionally, it provides a perspective on the research directions expected for the forthcoming years. The spine is composed of distinct and complex tissues that require specific modelling approaches. With the advent of additive manufacturing and increasing computational power, patient-specific treatments have moved from being a research trend to a reality in clinical practice, but there are many issues to be addressed before such approaches become universal. Here, it is identified that the major setback resides in validation of these computational techniques prior to approval by regulatory agencies. Nevertheless, there are very promising indicators in terms of optimised scaffold modelling for both disc arthroplasty and vertebroplasty, powered by a decisive contribution from imaging methods.

Effect of pedicle screw angles on the fracture risk of the human vertebra: A patient-specific computational model

The assessment of a human vertebra's stability after a screws fixation procedure and its fracture risk is still an open clinical problem. The accurate evaluation of fracture risk requires that all fracture mechanical determinants such as geometry, constitutive behavior, loading modes, and screws angulation are accounted for, which requires biomechanics-based analyses. As such, in the present work we investigate the effect of pedicle screws angulation on a patient-specific model of non osteoporotic lumbar vertebra, derived from clinical CT images. We propose a novel computational approach of fracture analysis and compare the effects of fixation stability in the lumbar spine. We considered a CT-based three-dimensional FE model of bilaterally instrumented L4 vertebra virtually implanting pedicle screws according to clinical guidelines. Nine screws trajectories were selected combining three craniocaudal and mediolateral angles, thus investigated through a parametric computational analysis. Bone was modeled as an elastic material with element-wise inhomogeneous properties fine-tuned on CT data. We implemented a custom algorithm to identify the thin cortical layer correctly from CT images ensuring reliable material properties in the computational model. Physiological motion (i.e., flexion, extension, axial rotation, lateral bending) was further accomplished by simultaneously loading the vertebra and the implant. We simulated local progressive damage of the bone by using a quasi-static force-driven incremental approach and considering a stress-based fracture criterion. Ductile-like and brittle-like fractures were found. Statistical analyses show significant differences comparing screws trajectories and averaging the results among six loading modes. In particular, we identified the caudomedial trajectory as the least critical case, thus safer from a clinical perspective. Instead, medial and craniolaterally oriented screws entailed higher peak and average stresses, though no statistical evidence classified such loads as the most critical scenarios. A quantitative validation procedure will be required in the future to translate our findings into clinical practice. Besides, to apply the results to the target osteoporotic population, new studies will be needed, including a specimen from an osteoporotic patient and the effect of osteoporosis on the constitutive model of bone.

Conventional finite element models estimate the strength of metastatic human vertebrae despite alterations of the bone's tissue and structure

INTRODUCTION

Pathologic vertebral fractures are a major clinical concern in the management of cancer patients with metastatic spine disease. These fractures are a direct consequence of the effect of bone metastases on the anatomy and structure of the vertebral bone. The goals of this study were twofold. First, we evaluated the effect of lytic, blastic and mixed (both lytic and blastic) metastases on the bone structure, on its material properties, and on the overall vertebral strength. Second, we tested the ability of bone mineral content (BMC) measurements and standard FE methodologies to predict the strength of real metastatic vertebral bodies.

METHODS

Fifty-seven vertebral bodies from eleven cadaver spines containing lytic, blastic, and mixed metastatic lesions from donors with breast, esophageal, kidney, lung, or prostate cancer were scanned using micro-computed tomography (μCT). Based on radiographic review, twelve vertebrae were selected for nanoindentation testing, while the remaining forty-five vertebrae were used for assessing their compressive strength. The μCT reconstruction was exploited to measure the vertebral BMC and to establish two finite element models. 1) a micro finite element (μFE) model derived at an image resolution of 24.5 μm and 2) homogenized FE (hFE) model derived at a resolution of 0.98 mm. Statistical analyses were conducted to measure the effect of the bone metastases on BV/TV, indentation modulus (E_it), ratio of plastic/total work (W_Pl/W_tot), and in vitro vertebral strength (F_exp). The predictive value of BMC, μFE stiffness, and hFE strength were evaluated against the in vitro measurements.

RESULTS

Blastic vertebral bodies exhibit significantly higher BV/TV compared to the mixed (p = 0.0205) and lytic (p = 0.0216) vertebral bodies. No significant differences were found between lytic and mixed vertebrae (p = 0.7584). Blastic bone tissue exhibited a 5.8% lower median E_it (p< 0.001) and a 3.3% lower median W_pl/W_tot (p<0.001) compared to non-involved bone tissue. No significant differences were measured between lytic and non-involved bone tissues. F_exp ranged from 1.9 to 13.8 kN, was strongly associated with hFE strength (R²=0.78, p< 0.001) and moderately associated with BMC (R²=0.66, p< 0.001) and μFE stiffness (R²=0.66, p< 0.001), independently of the lesion type.

DISCUSSION

Our findings show that tumour-induced osteoblastic metastases lead to slightly, but significantly lower bone tissue properties compared to controls, while osteolytic lesions appear to have a negligible impact. These effects may be attributed to the lower mineralization and woven nature of bone forming in blastic lesions whilst the material properties of bone in osteolytic vertebrae appeared little changed. The moderate association between BMC- and FE-based predictions to fracture strength suggest that vertebral strength is affected by the changes of bone mass induced by the metastatic lesions, rather than altered tissue properties. In a broader context, standard hFE approaches generated from CTs at clinical resolution are robust to the lesion type when predicting vertebral strength. These findings open the door for the development of FE-based prediction tools that overcomes the limitations of BMC in accounting for shape and size of the metastatic lesions. Such tools may help clinicians to decide whether a patient needs the prophylactic fixation of an impending fracture.

Numerical implementation of an osmo-poro-visco-hyperelastic finite element solver: application to the intervertebral disc

This work deals with the finite element (FE) implementation of a biphasic poroelastic formulation specifically developed to address the intricate behaviour of the Intervertebral Disc (IVD) and other highly hydrated soft tissues. This formulation is implemented in custom FE solver V-Biomech, being the validation performed with a lumbar IVD model, which was compared against the analogous FE model of Williams et al. and the experiments of Tyrrell et al. Good agreement with these benchmarks was achieved, meaning that V-Biomech and its novel poroelastic formulation are a viable alternative for simulation of biphasic soft tissues.

Assessing the Mechanical Weakness of Vertebrae Affected by Primary Tumors: A Feasibility Study

Patients spend months between the primary spinal tumor diagnosis and the surgical treatment, due to the need for performing chemotherapy and/or radiotherapy. During this period, they are exposed to an unknown risk of fracture. The aim of this study was to assess if it is possible to measure the mechanical strain in vertebrae affected by primary tumors, so as to open the way to an evidence-based scoring or prediction tool. We performed biomechanical tests on three vertebrae with bone tumor removed from patients. The tests were designed so as not to compromise the standard surgical and diagnostic procedures. Non-destructive mechanical tests in combination with state-of-the-art digital image correlation allowed to measure the distribution of strain on the surface of the vertebra. Our study has shown that the strains in the tumor region is circa 3 times higher than in the healthy bones, with principal strain peaks of 40,000/-20,000 microstrain, indicating a stress concentration potentially triggering vertebral fracture. This study has proven it is possible to analyze the mechanical behavior of primary tumor vertebrae as part of the clinical treatment protocol. This will allow building a tool for quantifying the risk of fracture and improving decision making in spine tumors.

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.