A patient-specific finite element methodology to predict damage accumulation in vertebral bodies under axial compression, sagittal flexion and combined loads

Due to the inherent limitations of DXA, assessment of the biomechanical properties of vertebral bodies relies increasingly on CT-based finite element (FE) models, but these often use simplistic material behaviour and/or single loading cases. In this study, we applied a novel constitutive law for bone elasticity, plasticity and damage to FE models created from coarsened pQCT images of human vertebrae, and compared vertebral stiffness, strength and damage accumulation for axial compression, anterior flexion and a combination of these two cases. FE axial stiffness and strength correlated with experiments and were linearly related to flexion properties. In all loading modes, damage localised preferentially in the trabecular compartment. Damage for the combined loading was higher than cumulated damage produced by individual compression and flexion. In conclusion, this FE method predicts stiffness and strength of vertebral bodies from CT images with clinical resolution and provides insight into damage accumulation in various loading modes.

Embedding of human vertebral bodies leads to higher ultimate load and altered damage localisation under axial compression

Computer tomography (CT)-based finite element (FE) models of vertebral bodies assess fracture load in vitro better than dual energy X-ray absorptiometry, but boundary conditions affect stress distribution under the endplates that may influence ultimate load and damage localisation under post-yield strains. Therefore, HRpQCT-based homogenised FE models of 12 vertebral bodies were subjected to axial compression with two distinct boundary conditions: embedding in polymethylmethalcrylate (PMMA) and bonding to a healthy intervertebral disc (IVD) with distinct hyperelastic properties for nucleus and annulus. Bone volume fraction and fabric assessed from HRpQCT data were used to determine the elastic, plastic and damage behaviour of bone. Ultimate forces obtained with PMMA were 22% higher than with IVD but correlated highly (R2 = 0.99). At ultimate force, distinct fractions of damage were computed in the endplates (PMMA: 6%, IVD: 70%), cortex and trabecular sub-regions, which confirms previous observations that in contrast to PMMA embedding, failure initiated underneath the nuclei in healthy IVDs. In conclusion, axial loading of vertebral bodies via PMMA embedding versus healthy IVD overestimates ultimate load and leads to distinct damage localisation and failure pattern.

Region Specific Response of Intervertebral Disc Cells to Complex Dynamic Loading: An Organ Culture Study Using a Dynamic Torsion-Compression Bioreactor

The spine is routinely subjected to repetitive complex loading consisting of axial compression, torsion, flexion and extension. Mechanical loading is one of the important causes of spinal diseases, including disc herniation and disc degeneration. It is known that static and dynamic compression can lead to progressive disc degeneration, but little is known about the mechanobiology of the disc subjected to combined dynamic compression and torsion. Therefore, the purpose of this study was to compare the mechanobiology of the intervertebral disc when subjected to combined dynamic compression and axial torsion or pure dynamic compression or axial torsion using organ culture. We applied four different loading modalities 1. control: no loading (NL), 2. cyclic compression (CC), 3. cyclic torsion (CT), and 4. combined cyclic compression and torsion (CCT) on bovine caudal disc explants using our custom made dynamic loading bioreactor for disc organ culture. Loads were applied for 8 h/day and continued for 14 days, all at a physiological magnitude and frequency. Our results provided strong evidence that complex loading induced a stronger degree of disc degeneration compared to one degree of freedom loading. In the CCT group, less than 10\% nucleus pulposus (NP) cells survived the 14 days of loading, while cell viabilities were maintained above 70\% in the NP of all the other three groups and in the annulus fibrosus (AF) of all the groups. Gene expression analysis revealed a strong up-regulation in matrix genes and matrix remodeling genes in the AF of the CCT group. Cell apoptotic activity and glycosaminoglycan content were also quantified but there were no statistically significant differences found. Cell morphology in the NP of the CCT was changed, as shown by histological evaluation. Our results stress the importance of complex loading on the initiation and progression of disc degeneration.

Classification of thoracolumbar fractures and dislocations

A classification of injuries is necessary in order to develop a common language for treatment indications and outcomes. Several classification systems have been proposed, the most frequently used is the Denis classification. The problem of this classification system is that it is based on an assumption, which is anatomically unidentifiable: the so-called middle column. For this reason, few years ago, a group of spine surgeons has developed a new classification system, which is based on the severity of the injury. The severity is defined by the pathomorphological findings, the prognosis in terms of healing and potential of neurological damage. This classification is based on three major groups: A = isolated anterior column injuries by axial compression, B = disruption of the posterior ligament complex by distraction posteriorly, and group C = corresponding to group B but with rotation. There is an increasing severity from A to C, and within each group, the severity usually increases within the subgroups from .1, .2, .3. All these pathomorphologies are supported by a mechanism of injury, which is responsible for the extent of the injury. The type of injury with its groups and subgroups is able to suggest the treatment modality.

The effects of dynamic loading on the intervertebral disc

Loading is important to maintain the balance of matrix turnover in the intervertebral disc (IVD). Daily cyclic diurnal assists in the transport of large soluble factors across the IVD and its surrounding circulation and applies direct and indirect stimulus to disc cells. Acute mechanical injury and accumulated overloading, however, could induce disc degeneration. Recently, there is more information available on how cyclic loading, especially axial compression and hydrostatic pressure, affects IVD cell biology. This review summarises recent studies on the response of the IVD and stem cells to applied cyclic compression and hydrostatic pressure. These studies investigate the possible role of loading in the initiation and progression of disc degeneration as well as quantifying a physiological loading condition for the study of disc degeneration biological therapy. Subsequently, a possible physiological/beneficial loading range is proposed. This physiological/beneficial loading could provide insight into how to design loading regimes in specific system for the testing of various biological therapies such as cell therapy, chemical therapy or tissue engineering constructs to achieve a better final outcome. In addition, the parameter space of 'physiological' loading may also be an important factor for the differentiation of stem cells towards most ideally 'discogenic' cells for tissue engineering purpose.

Effect of reduced oxygen tension and long-term mechanical stimulation on chondrocyte-polymer constructs

We have investigated the influence of long-term confined dynamic compression and surface motion under low oxygen tension on tissue-engineered cell-scaffold constructs. Porous polyurethane scaffolds (8 mm x 4 mm) were seeded with bovine articular chondrocytes and cultured under normoxic (21% O(2)) or hypoxic (5% O(2)) conditions for up to 4 weeks. By means of our joint-simulating bioreactor, cyclic axial compression (10-20%; 0.5 Hz) was applied for 1 h daily with a ceramic ball, which simultaneously oscillated over the construct surface (+/-25 degrees; 0.5 Hz). Culture under reduced oxygen tension resulted in an increase in mRNA levels of type II collagen and aggrecan, whereas the expression of type I collagen was down-regulated at early time points. A higher glycosaminoglycan content was found in hypoxic than in normoxic constructs. Immunohistochemical analysis showed more intense type II and weaker type I collagen staining in hypoxic than in normoxic cultures. Type II collagen gene expression was slightly elevated after short-term loading, whereas aggrecan mRNA levels were not influenced by the applied mechanical stimuli. Of importance, the combination of loading and low oxygen tension resulted in a further down-regulation of collagen type I mRNA expression, contributing to the stabilization of the chondrocytic phenotype. Histological results confirmed the beneficial effect of mechanical loading on chondrocyte matrix synthesis. Thus, mechanical stimulation combined with low oxygen tension is an effective tool for modulating the chondrocytic phenotype and should be considered when chondrocytes or mesenchymal stem cells are cultured and differentiated with the aim of generating cartilage-like tissue in vitro.

Experimental folate and vitamin B12 deficiency does not alter bone quality in rats

Hyperhomocysteinemia (HHCY) has been linked to fragility fractures and osteoporosis. Folate and vitamin B(12) deficiencies are among the main causes of HHCY. However, the impact of these vitamins on bone health has been poorly studied. This study analyzed the effect of folate and vitamin B(12) deficiency on bone in rats. We used two groups of rats: a control group (Co, n = 10) and a vitamin-deficient group (VitDef, n = 10). VitDef animals were fed for 12 wk with a folate- and vitamin B(12)-free diet. Co animals received an equicaloric control diet. Tissue and plasma concentrations of homocysteine (HCY), S-adenosyl-homocysteine (SAH), and S-adenosyl-methionine (SAM) were measured. Bone quality was assessed by biomechanical testing (maximum force of an axial compression test; F(max)), histomorphometry (bone area/total area; B.Ar./T.Ar.], and the measurement of biochemical bone turnover markers (osteocalcin, collagen I C-terminal cross-laps [CTX]). VitDef animals developed significant HHCY (Co versus VitDef: 6.8 +/- 2.7 versus 61.1 +/- 12.8 microM, p < 0.001) that was accompanied by a high plasma concentration of SAH (Co versus VitDef: 24.1 +/- 5.9 versus 86.4 +/- 44.3 nM, p < 0.001). However, bone tissue concentrations of HCY, SAH, and SAM were similar in the two groups. Fmax, B.Ar./T.Ar., OC, and CTX did not differ between VitDef and Co animals, indicating that bone quality was not affected. Folate and vitamin B(12) deficiency induces distinct HHCY but has no effect on bone health in otherwise healthy adult rats. The unchanged HCY metabolism in bone is the most probable explanation for the missing effect of the vitamin-free diet on bone.

Augmentation of mechanical properties in osteoporotic vertebral bones--a biomechanical investigation of vertebroplasty efficacy with different bone cements

Recent clinical trials have reported favorable early results for transpedicular vertebral cement reinforcement of osteoporotic vertebral insufficiencies. There is, however, a lack of basic data on the application, safety and biomechanical efficacy of materials such as polymethyl-methacrylate (PMMA) and calciumphospate (CaP) cements. The present study analyzed 33 vertebral pairs from five human cadaver spines. Thirty-nine vertebrae were osteoporotic (bone mineral density < 0.75 g/cm2), 27 showed nearly normal values. The cranial vertebra of each pair was augmented with either PMMA (Palacos E-Flow) or experimental brushite cement (EBC), with the caudal vertebra as a control. PMMA and EBC were easy to inject, and vertebral fillings of 20-50% were achieved. The maximal possible filling was inversely correlated to the bone mineral density (BMD) values. Cement extrusion into the spinal canal was observed in 12% of cases. All specimens were subjected to axial compression tests in a displacement-controlled mode. From load-displacement curves, the stiffness, S, and the maximal force before failure, Fmax, were determined. Compared with the native control vertebrae, a statistically significant increase in vertebral stiffness and Fmax was observed by the augmentation. With PMMA the stiffness increased by 174% (P = 0.018) and Fmax by 195% (P = 0.001); the corresponding augmentation with EBC was 120% (P = 0.03) and 113% (P = 0.002). The lower the initial BMD, the more pronounced was the augmentation effect. Both PMMA and EBC augmentation reliably and significantly raised the stiffness and maximal tolerable force until failure in osteoporotic vertebral bodies. In non-porotic specimens, no significant increase was achieved.

Compressive strength of interbody cages in the lumbar spine: the effect of cage shape, posterior instrumentation and bone density

One goal of interbody fusion is to increase the height of the degenerated disc space. Interbody cages in particular have been promoted with the claim that they can maintain the disc space better than other methods. There are many factors that can affect the disc height maintenance, including graft or cage design, the quality of the surrounding bone and the presence of supplementary posterior fixation. The present study is an in vitro biomechanical investigation of the compressive behaviour of three different interbody cage designs in a human cadaveric model. The effect of bone density and posterior instrumentation were assessed. Thirty-six lumbar functional spinal units were instrumented with one of three interbody cages: (1) a porous titanium implant with endplate fit (Stratec), (2) a porous, rectangular carbon-fibre implant (Brantigan) and (3) a porous, cylindrical threaded implant (Ray). Posterior instrumentation (USS) was applied to half of the specimens. All specimens were subjected to axial compression displacement until failure. Correlations between both the failure load and the load at 3 mm displacement with the bone density measurements were observed. Neither the cage design nor the presence of posterior instrumentation had a significant effect on the failure load. The loads at 3 mm were slightly less for the Stratec cage, implying lower axial stiffness, but were not different with posterior instrumentation. The large range of observed failure loads overlaps the potential in vivo compressive loads, implying that failure of the bone-implant interface may occur clinically. Preoperative measurements of bone density may be an effective tool to predict settling around interbody cages.

In vitro biomechanical testing of a micro external skeletal fixator

OBJECTIVE To biomechanically test the properties of three different Universal Micro External Fixator (UMEX™) configurations with regard to their use in very small animals (<5kg) and compare the UMEX system to the widely used IMEX External Skeletal Fixation (SK™) system in terms of stiffness, space needed for pin placement and weight. METHODS Three different UMEX configurations (type Ia, type Ib, and type II modified) and one SK configuration type Ia were used to stabilize Delrin plastic rods in a 1 cm fracture gap model. These constructs were tested in axial compression, craniocaudal bending, mediolateral bending, and torsion. Testing was conducted within the elastic range and mean stiffness in each mode was determined from the slope of the linear portion of the load-deformation curve. A Kruskal Wallis one-way analysis of variance on ranks test was utilized to assess differences between constructs (p <0.05). RESULTS The UMEX type II modified configuration was significantly stiffer than the other UMEX configurations and the SK type Ia, except in craniocaudal bending, where the SK type Ia configuration was stiffer than all UMEX constructs. The UMEX type Ia configuration was significantly the weakest of those frames. The UMEX constructs were lighter and smaller than the SK, thus facilitating closer pin placement. CONCLUSIONS Results supported previous reports concerning the superiority of more complex constructs regarding stiffness. The UMEX system appears to be a valid alternative for the treatment of long-bone fractures in very small animals.

Finite element analysis for prediction of bone strength

Article preview View full access options BoneKEy Reports | Review Print Email Share/bookmark Finite element analysis for prediction of bone strength Philippe K Zysset, Enrico Dall'Ara, Peter Varga & Dieter H Pahr Affiliations Corresponding author BoneKEy Reports (2013) 2, Article number: 386 (2013) doi:10.1038/bonekey.2013.120 Received 03 January 2013 Accepted 25 June 2013 Published online 07 August 2013 Article tools Citation Reprints Rights & permissions Abstract Abstract• References• Author information Finite element (FE) analysis has been applied for the past 40 years to simulate the mechanical behavior of bone. Although several validation studies have been performed on specific anatomical sites and load cases, this study aims to review the predictability of human bone strength at the three major osteoporotic fracture sites quantified in recently completed in vitro studies at our former institute. Specifically, the performance of FE analysis based on clinical computer tomography (QCT) is compared with the ones of the current densitometric standards, bone mineral content, bone mineral density (BMD) and areal BMD (aBMD). Clinical fractures were produced in monotonic axial compression of the distal radii, vertebral sections and in side loading of the proximal femora. QCT-based FE models of the three bones were developed to simulate as closely as possible the boundary conditions of each experiment. For all sites, the FE methodology exhibited the lowest errors and the highest correlations in predicting the experimental bone strength. Likely due to the improved CT image resolution, the quality of the FE prediction in the peripheral skeleton using high-resolution peripheral CT was superior to that in the axial skeleton with whole-body QCT. Because of its projective and scalar nature, the performance of aBMD in predicting bone strength depended on loading mode and was significantly inferior to FE in axial compression of radial or vertebral sections but not significantly inferior to FE in side loading of the femur. Considering the cumulated evidence from the published validation studies, it is concluded that FE models provide the most reliable surrogates of bone strength at any of the three fracture sites.

Early changes in biochemical markers of bone formation during teriparatide therapy correlate with improvements in vertebral strength in men with glucocorticoid-induced osteoporosis

Summary Changes of the bone formation marker PINP correlated positively with improvements in vertebral strength in men with glucocorticoid-induced osteoporosis (GIO) who received 18-month treatment with teriparatide, but not with risedronate. These results support the use of PINP as a surrogate marker of bone strength in GIO patients treated with teriparatide. Introduction To investigate the correlations between biochemical markers of bone turnover and vertebral strength estimated by finite element analysis (FEA) in men with GIO. Methods A total of 92 men with GIO were included in an 18-month, randomized, open-label trial of teriparatide (20 μg/day, n = 45) and risedronate (35 mg/week, n = 47). High-resolution quantitative computed tomography images of the 12th thoracic vertebra obtained at baseline, 6 and 18 months were converted into digital nonlinear FE models and subjected to anterior bending, axial compression and torsion. Stiffness and strength were computed for each model and loading mode. Serum biochemical markers of bone formation (amino-terminal-propeptide of type I collagen [PINP]) and bone resorption (type I collagen cross-linked C-telopeptide degradation fragments [CTx]) were measured at baseline, 3 months, 6 months and 18 months. A mixed-model of repeated measures analysed changes from baseline and between-group differences. Spearman correlations assessed the relationship between changes from baseline of bone markers with FEA variables. Results PINP and CTx levels increased in the teriparatide group and decreased in the risedronate group. FEA-derived parameters increased in both groups, but were significantly higher at 18 months in the teriparatide group. Significant positive correlations were found between changes from baseline of PINP at 3, 6 and 18 months with changes in FE strength in the teriparatide-treated group, but not in the risedronate group. Conclusions Positive correlations between changes in a biochemical marker of bone formation and improvement of biomechanical properties support the use of PINP as a surrogate marker of bone strength in teriparatide-treated GIO patients.

High resolution quantitative computed tomography-based assessment of trabecular microstructure and strength estimates by finite-element analysis of the spine, but not DXA, reflects vertebral fracture status in men with glucocorticoid-induced osteoporosis

High-resolution quantitative computed tomography (HRQCT)-based analysis of spinal bone density and microstructure, finite element analysis (FEA), and DXA were used to investigate the vertebral bone status of men with glucocorticoid-induced osteoporosis (GIO). DXA of L1–L3 and total hip, QCT of L1–L3, and HRQCT of T12 were available for 73 men (54.6±14.0years) with GIO. Prevalent vertebral fracture status was evaluated on radiographs using a semi-quantitative (SQ) score (normal=0 to severe fracture=3), and the spinal deformity index (SDI) score (sum of SQ scores of T4 to L4 vertebrae). Thirty-one (42.4%) subjects had prevalent vertebral fractures. Cortical BMD (Ct.BMD) and thickness (Ct.Th), trabecular BMD (Tb.BMD), apparent trabecular bone volume fraction (app.BV/TV), and apparent trabecular separation (app.Tb.Sp) were analyzed by HRQCT. Stiffness and strength of T12 were computed by HRQCT-based nonlinear FEA for axial compression, anterior bending and axial torsion. In logistic regressions adjusted for age, glucocorticoid dose and osteoporosis treatment, Tb.BMD was most closely associated with vertebral fracture status (standardized odds ratio [sOR]: Tb.BMD T12: 4.05 [95% CI: 1.8–9.0], Tb.BMD L1–L3: 3.95 [1.8–8.9]). Strength divided by cross-sectional area for axial compression showed the most significant association with spine fracture status among FEA variables (2.56 [1.29–5.07]). SDI was best predicted by a microstructural model using Ct.Th and app.Tb.Sp (r2=0.57, p<0.001). Spinal or hip DXA measurements did not show significant associations with fracture status or severity. In this cross-sectional study of males with GIO, QCT, HRQCT-based measurements and FEA variables were superior to DXA in discriminating between patients of differing prevalent vertebral fracture status. A microstructural model combining aspects of cortical and trabecular bone reflected fracture severity most accurately.

Finite element analyses of human vertebral bodies embedded in polymethylmethalcrylate or loaded via the hyperelastic intervertebral disc models provide equivalent predictions of experimental strength

Quantitative computer tomography (QCT)-based finite element (FE) models of vertebral body provide better prediction of vertebral strength than dual energy X-ray absorptiometry. However, most models were validated against compression of vertebral bodies with endplates embedded in polymethylmethalcrylate (PMMA). Yet, loading being as important as bone density, the absence of intervertebral disc (IVD) affects the strength. Accordingly, the aim was to assess the strength predictions of the classic FE models (vertebral body embedded) against the in vitro and in silico strengths of vertebral bodies loaded via IVDs. High resolution peripheral QCT (HR-pQCT) were performed on 13 segments (T11/T12/L1). T11 and L1 were augmented with PMMA and the samples were tested under a 4° wedge compression until failure of T12. Specimen-specific model was generated for each T12 from the HR-pQCT data. Two FE sets were created: FE-PMMA refers to the classical vertebral body embedded model under axial compression; FE-IVD to their loading via hyperelastic IVD model under the wedge compression as conducted experimentally. Results showed that FE-PMMA models overestimated the experimental strength and their strength prediction was satisfactory considering the different experimental set-up. On the other hand, the FE-IVD models did not prove significantly better (Exp/FE-PMMA: R²=0.68; Exp/FE-IVD: R²=0.71, p=0.84). In conclusion, FE-PMMA correlates well with in vitro strength of human vertebral bodies loaded via real IVDs and FE-IVD with hyperelastic IVDs do not significantly improve this correlation. Therefore, it seems not worth adding the IVDs to vertebral body models until fully validated patient-specific IVD models become available.