The effect of non-fusion dynamic stabilization on biomechanical responses of the implanted lumbar spine during whole-body vibration

Stepwise reduction of bony density in patients induces a higher risk of annular tears by deteriorating the local biomechanical environment

BACKGROUND CONTEXT

The relationship between osteoporosis and intervertebral disc degeneration (IDD) remains unclear. Considering that annular tear is the primary phenotype of IDD in the lumbar spine, the deteriorating local biomechanical environment may be the main trigger for annular tears.

PURPOSE

To investigate whether poor bone mineral density (BMD) in the vertebral bodies may increase the risk of annular tears via the degradation of the local biomechanical environment.

STUDY DESIGN

This study was a retrospective investigation with relevant numerical mechanical simulations.

PATIENT SAMPLE

A total of 64 patients with low back pain (LBP) and the most severe IDD in the L4-L5 motion segment were enrolled.

OUTCOME MEASURES

Annulus integration status was assessed using diffusion tensor fibre tractography (DTT). Hounsfield unit (HU) values of adjacent vertebral bodies were employed to determine BMD. Numerical simulations were conducted to compute stress values in the annulus of models with different BMDs and body positions.

METHODS

The clinical data of the 64 patients with low back pain were collected retrospectively. The BMD of the vertebral bodies was measured using the HU values, and the annulus integration status was determined according to DTT. The data of the patients with and without annular tears were compared, and regression analysis was used to identify the independent risk factors for annular tears. Furthermore, finite element models of the L4-L5 motion segment were constructed and validated, followed by estimating the maximum stress on the post and postlateral interfaces between the superior and inferior bony endplates (BEPs) and the annulus.

RESULTS

Patients with lower HU values in their vertebral bodies had significantly higher incidence rates of annular tears, with decreased HU values being an independent risk factor for annular tears. Moreover, increased stress on the BEP-annulus interfaces was associated with a stepwise reduction of bony density (ie, elastic modulus) in the numerical models.

CONCLUSIONS

The stepwise reduction of bony density in patients results in a higher risk of annular tears by deteriorating the local biomechanical environment. Thus, osteoporosis should be considered to be a potential risk factor for IDD biomechanically.

Customized design and biomechanical property analysis of 3D-printed tantalum intervertebral cages

BACKGROUND

Intervertebral cages used in clinical applications were often general products with standard specifications, which were challenging to match with the cervical vertebra and prone to cause stress shielding and subsidence.

OBJECTIVE

To design and fabricate customized tantalum (Ta) intervertebral fusion cages that meets the biomechanical requirements of the cervical segment.

METHODS

The lattice intervertebral cages were customized designed and fabricated by the selective laser melting. The joint and muscle forces of the cervical segment under different movements were analyzed using reverse dynamics method. The stress characteristics of cage, plate, screws and vertebral endplate were analyzed by finite element analysis. The fluid flow behaviors and permeability of three lattice structures were simulated by computational fluid dynamics. Compression tests were executed to investigate the biomechanical properties of the cages.

RESULTS

Compared with the solid cages, the lattice-filled structures significantly reduced the stress of cages and anterior fixation system. In comparison to the octahedroid and quaddiametral lattice-filled cages, the bitriangle lattice-filled cage had a lower stress shielding rate, higher permeability, and superior subsidence resistance ability.

CONCLUSION

The inverse dynamics simulation combined with finite element analysis is an effective method to investigate the biomechanical properties of the cervical vertebra during movements.

Treatment of Lumbar Degenerative Disease with a Novel Interlaminar Screw Elastic Spacer Technique: A Finite Element Analysis

A novel interlaminar elastic screw spacer technique was designed to maintain lumbar mobility in treating lumbar degenerative diseases. A validated finite element model of L4/5 was used to establish an ISES-1/2 model and an ISES-1/3 model based on different insertion points, a unilateral fixation model and a bilateral fixed model based on different fixation methods, and a Coflex-F model based on different implants. The elastic rods were used to fix screws. Under the same mechanical conditions, we compared the biomechanical characteristics to investigate the optimal entry point for ISES technology, demonstrate the effectiveness of unilateral fixation, and validate the feasibility of the ISES technique. Compared to ISES-1/3, the ISES-1/2 model had lower intradiscal pressure, facet cartilage stress, and posterior structural stress. Compared to the ISES-BF model, the ISES-UF model had lower intervertebral pressure, larger mobility, and smaller stress on the posterior structures. The ISES model had a similar intervertebral pressure and limitation of extension as the Coflex-F model. The ISES model retained greater mobility and reduced the stress on the facet cartilage and posterior structure compared with the Coflex-F model. Our study suggests that the ISES technique is a promising treatment of lumbar degenerative diseases, especially those with osteoporosis.

Biomechanical Evaluation of Rigid Interspinous Process Fixation Combined With Lumbar Interbody Fusion Using Hybrid Testing Protocol

Rigid interspinous process fixation (RIPF) has been recently discussed as an alternative to pedicle screw fixation (PSF) for reducing trauma in lumbar interbody fusion (LIF) surgery. This study aimed to investigate biomechanics of the lumbar spine with RIPF, and also to compare biomechanical differences between two postoperative stages (before and after bony fusion). Based on an intact finite-element model of lumbosacral spine, the models of single-level LIF with RIPF or conventional PSF were developed and were computed for biomechanical responses to the moments of four physiological motions using hybrid testing protocol. It was found that compared with PSF, range of motion (ROM), intradiscal pressure (IDP), and facet joint forces (FJF) at adjacent segments of the surgical level for RIPF were decreased by up to 8.4%, 2.3%, and 16.8%, respectively, but ROM and endplate stress at the surgical segment were increased by up to 285.3% and 174.3%, respectively. The results of comparison between lumbar spine with RIPF before and after bony fusion showed that ROM and endplate stress at the surgical segment were decreased by up to 62.6% and 40.4%, respectively, when achieved to bony fusion. These findings suggest that lumbar spine with RIPF as compared to PSF has potential to decrease the risk of adjacent segment degeneration but might have lower stability of surgical segment and an increased risk of cage subsidence; When achieved bony fusion, it might be helpful for the lumbar spine with RIPF in increasing stability of surgical segment and reducing failure of bone contact with cage.

A swelling-based biphasic analysis on the quasi-static biomechanical behaviors of healthy and degenerative intervertebral discs

BACKGROUND AND OBJECTIVE

The degeneration of intervertebral discs is significantly dependent of the changes in tissue composition ratio and tissue structure. Up to the present, the effects of degeneration on the quasi-static biomechanical responses of discs have not been well understood. The goal of this study is to quantitatively analyze the quasi-static responses of healthy and degenerative discs.

METHODS

Four biphasic swelling-based finite element models are developed and quantitatively validated. Four quasi-static test protocols, including the free-swelling, slow-ramp, creep and stress-relaxation, are implemented. The double Voigt and double Maxwell models are further used to extract the immediate (or residual), short-term and long-term responses of these tests.

RESULTS

Simulation results show that both the swelling-induced pressure in the nucleus pulposus and the initial modulus decrease with degeneration. In the free-swelling test of discs possessing healthy cartilage endplates, simulation results show that over 80% of the total strain is contributed by the short-term response. The long-term response is dominant for discs with degenerated permeability in cartilage endplates. For the creep test, over 50% of the deformation is contributed by the long-term response. In the stress-relaxation test, the long-term stress contribution occupies approximately 31% of total response and is independent of degeneration. Both the residual and short-term responses vary monotonically with degeneration. In addition, both the glycosaminoglycan content and permeability affect the engineering equilibrium time constants of the rheologic models, in which the determining factor is the permeability.

CONCLUSIONS

The content of glycosaminoglycan in intervertebral soft tissues and the permeability of cartilage endplates are two critical factors that affect the fluid-dependent viscoelastic responses of intervertebral discs. The component proportions of the fluid-dependent viscoelastic responses depend also strongly on test protocols. In the slow-ramp test, the glycosaminoglycan content is responsible for the changes of the initial modulus. Since existing computational models simulate disc degenerations only by altering disc height, boundary conditions and material stiffness, the current work highlights the significance of biochemical composition and cartilage endplates permeability in the biomechanical behaviors of degenerated discs.

Patterns of Vertebral Bone Marrow Edema in the Normal Healing Process of Lumbar Interbody Fusion: Baseline Data for Diagnosis of Pathological Events

STUDY DESIGN

Retrospective investigation using a prospectively collected database.

OBJECTIVE

To examine the appearance and characteristics of vertebral bone marrow edema (BME) in the normal healing of lumbar interbody fusion.

SUMMARY OF BACKGROUND DATA

Although BME in pathological spinal conditions has been well-documented, the patterns and characteristics of BME in the normal healing process of spinal fusion remains unexplored.

MATERIALS AND METHODS

We reviewed imaging from 225 patients with normal healing following posterior lumbar interbody fusion or transforaminal lumbar interbody fusion. BME was identified on magnetic resonance imaging at the third postoperative week and categorized with respect to its appearance, including assessment of area and extension within the relevant vertebrae.

RESULTS

Three hundred eighty-nine of the 450 instrumented vertebrae (86.4%) displayed evidence BME. All instances of BME were associated with the area of contact with the endplate. The average extent of BME was 32.7±1.0%. BME within normal healing following interbody fusion could be categorized into four types: no edema (13.6%), anterior corner (36.6%), around-the-cage focal (48.0%), and diffuse (1.8%). Anterior corner BME was significantly associated with instances of single cage placement than in dual cages (42.6% vs. 24.7%, P =0.0002). Single cages had a significantly higher rate of BME than dual cages (92.0% vs. 75.3%, P <0.0001). The extent of BME was significantly greater in the single cage cohort (36.9% vs. 24.2% in dual cages, P <0.0001).

CONCLUSIONS

This serves as the first study demonstrating the patterns of BME associated with normal healing following lumbar interbody fusion procedures. Anterior corner BME and around-the-cage focal BME were the most common patterns encountered, with diffuse BME a relatively rare pattern. These findings might contribute to the better differentiation of postoperative pathological events from normal healing following lumbar interbody fusion.

LEVEL OF EVIDENCE

4.

Biomechanical analysis of lumbar nonfusion dynamic stabilization using a pedicle screw-based dynamic stabilizer or an interspinous process spacer

This study aimed to investigate and compare the effects of two widely used nonfusion posterior dynamic stabilization (NPDS) devices, pedicle screw-based dynamic stabilizer (PSDS) and interspinous process spacer (IPS), on biomechanics of the implanted lumbar spine under static and vibration loadings. The finite element model of healthy human lumbosacral segment was modified to incorporate NPDS device insertion at L4-L5 segment. Bioflex and DIAM were used as PSDS-based and IPS-based NPDS devices, respectively. As a comparison, lumbar interbody fusion with rigid stabilization was also simulated at L4-L5. For static loading, segmental range of motion (ROM) of the models under moments of four physiological motions was computed using hybrid testing protocol. For vibration loading, resonant modes and dynamic stress of the models under vertical excitation were extracted through random response analysis. The results showed that compared with the rigid fusion model, ROM of the nonfusion models was higher at L4-L5 level but lower at adjacent levels (L1- L2, L2-L3, L3-L4, L5-S1). Compared with the Bioflex model, the DIAM model produced higher ROM at L4-L5 level but lower ROM at adjacent levels, especially under lateral bending and axial rotation; resonant frequency of the DIAM model was slightly lower; dynamic response of nucleus stress at L4-L5 level was slightly higher for the DIAM model, and the dynamic stress at adjacent levels was no obvious difference between the nonfusion models. This study reveals biomechanical differences between the Bioflex and DIAM systems, which may provide references for selecting surgical approaches in clinical practice.

Single-level cervical disc arthroplasty in the spine with reversible kyphosis: A finite element study

Background

Our previous studies found the single‐level cervical disc arthroplasty (CDA) might be a feasible treatment for the patients with reversible kyphosis (RK). Theoretically, the change of cervical alignment from lordosis to RK comes with the biomechanical alteration of prostheses and cervical spine. However, the biomechanical data of CDA in the spine with RK have not been reported. This study aimed at establishing finite element (FE) models to (1) explore the effects of RK on the biomechanics of artificial cervical disc; (2) investigate the biomechanical differences of single‐level anterior cervical discectomy and fusion (ACDF) and CDA in the cervical spine with RK.

Methods

The FE models of the cervical spine with lordosis and RK were constructed, then three single‐level surgical models were developed: (1) RK + ACDF; (2) RK + CDA; (3) lordosis + CDA. A 73.6‐N follower load combined with 1 N·m moment was applied at the C2 vertebra to produce cervical motion.

Results

At the surgical level, “lordosis + CDA” had the greatest ROM (except for flexion) while “RK + ACDF” had the minimum ROM. However, at adjacent levels, the ROM of “RK + ACDF” increased by 4.05% to 38.04% in comparison to “RK + CDA.” “RK + ACDF” had the greatest prosthesis interface stress, while the maximum prosthesis interface stress of “RK + CDA” was at least 2.15 times higher than “lordosis + CDA.” Similarly, “RK + ACDF” had the greatest intradiscal pressure (IDP) at adjacent levels, while the IDP of “RK + CDA” was 1.6 to 6.7 times higher than “lordosis + CDA.” At the surgical level, “RK + CDA” had the greatest facet joint stress (except for extension), which was 1.9 to 11.2 times higher than “lordosis + CDA.” At the adjacent levels, “RK + CDA” had the greatest facet joint stress (except for extension), followed by “RK + ACDF” and “lordosis + CDA” in descending order.

Conclusions

RK significantly changed the biomechanics of CDA, which is demonstrated by the decreased ROM and the significantly increased prosthesis interface stress, IDP, and facet joint stress in the “RK + CDA” model. Compared with ACDF, CDA overall exhibited a better biomechanical performance in the cervical spine with RK, with the increased ROM of surgical level and facet joint stress and the decreased ROM of adjacent levels, prosthesis interface stress, and IDP.

Biomaterials for Interbody Fusion in Bone Tissue Engineering

In recent years, interbody fusion cages have played an important role in interbody fusion surgery for treating diseases like disc protrusion and spondylolisthesis. However, traditional cages cannot achieve satisfactory results due to their unreasonable design, poor material biocompatibility, and induced osteogenesis ability, limiting their application. There are currently 3 ways to improve the fusion effect, as follows. First, the interbody fusion cage is designed to facilitate bone ingrowth through the preliminary design. Second, choose interbody fusion cages made of different materials to meet the variable needs of interbody fusion. Finally, complete post-processing steps, such as coating the designed cage, to achieve a suitable osseointegration microstructure, and add other bioactive materials to achieve the most suitable biological microenvironment of bone tissue and improve the fusion effect. The focus of this review is on the design methods of interbody fusion cages, a comparison of the advantages and disadvantages of various materials, the influence of post-processing techniques and additional materials on interbody fusion, and the prospects for the future development of interbody fusion cages.

Biomechanical Investigation of Lumbar Interbody Fusion Supplemented with Topping-off Instrumentation Using Different Dynamic Stabilization Devices

STUDY DESIGN

A biomechanical comparison study using finite element method.

OBJECTIVE

The aim of this study was to investigate effects of different dynamic stabilization devices, including pedicle-based dynamic stabilization system (PBDSS) and interspinous process spacer (ISP), used for topping-off implants on biomechanical responses of human spine after lumbar interbody fusion.

SUMMARY OF BACKGROUND DATA

Topping-off stabilization technique has been proposed to prevent adjacent segment degeneration following lumbar spine fusion. PBDSS and ISP are the most used dynamic stabilizers for topping-off instrumentation. However, biomechanical differences between them still remain unclear.

METHODS

A validated, normal FE model of human lumbosacral spine was employed. Based on this model, rigid fusion at L4-L5 and moderately disc degeneration at L3-L4 were simulated and used as a comparison baseline. Subsequently, Bioflex and DIAM systems were instrumented at L3-L4 segment to construct PBDSS-based and ISP-based topping-off models. Biomechanical responses of the models to bending moments and vertical vibrational excitation were computed using FE static and random response analyses, respectively.

RESULTS

Results from static analysis showed that at L3-L4, the response parameters including annulus stress and range of motion were decreased by 41.6% to 85.2% for PBDSS-based model and by 6.3% to 67% for ISP-based model compared with rigid fusion model. At L2-L3, these parameters were lower in ISP-based model than in PBDSS-based model. Results from random response analysis showed that topping-off instrumentation increased resonant frequency of spine system but decreased dynamic response of annulus stress at L3-L4. PBDSS-based model generated lower dynamic stress than ISP-based model at L3-L4, but the dynamic stress was higher at L2-L3 for PBDSSbased model.

CONCLUSION

Under static and vibration loadings, the PBDSSbased topping-off device (Bioflex) provided a better protection for transition segment, and likelihood of degeneration of supraadjacent segment might be relatively lower when using the ISPbased topping-off device (DIAM).Level of Evidence: 5.

Prediction of the natural frequencies of different degrees of degenerated human lumbar segments L2-L3 using dynamic finite element analysis

BACKGROUND AND OBJECTIVE

Chronic exposure to resonant environment may cause more serious injuries to human lumbar spine than other vibrations. On the condition that the natural frequency of human lumbar spine is known, excitation frequency from an external vibration source can be optimized to keep away from the natural frequency and thus avoid lumbar resonance. Therefore, this study aimed to present an approach to predict the natural frequency of the human lumbar spine.

METHODS

Four poroelastic finite element models of human L2-L3 spinal motion segments with different degrees of degeneration were established. Dynamic finite element analyses of these models during 1 h of vibration were then conducted. The mechanical parameters of these models under vibrations at different excitation frequencies were predicted. The excitation frequencies that resulted in the greatest changes in the lumbar mechanical parameters were identified as the natural frequencies of the established L2-L3 spinal motion segments.

RESULTS

Simulation results showed that the natural frequencies of the healthy and mildly degenerated L2-L3 spinal motion segments, moderately degenerated L2-L3 spinal motion segments, and seriously degenerated L2-L3 spinal motion segments were in the range of 5-7, 3-5, and 1-3 Hz, respectively.

CONCLUSIONS

The predicted results indicated that the natural frequencies of the human L2-L3 spinal motion segments gradually decreased with the severity of degeneration. These phenomena may be related to changes in the lumbar structures and materials because of degeneration. This study provided a feasible method to predict the lumbar natural frequencies for different populations, which may be helpful in optimizing external vibration sources to avoid lumbar resonance.

Stress analysis of intervertebral disc during occupational activities

BACKGROUND AND OBJECTIVE

Manual material handling activities cause large compression of the intervertebral disc of the lumbar spine. Intradiscal pressure (IDP) has generally been employed to predict the risk of low back injury. As an alternative to in vivo measurements, either motion analysis or finite element (FE) analysis has been used to estimate IDP. The purpose of this study is to propose a new biomechanical method that integrates FE analysis with motion analysis, in order to estimate the stresses and deformations of the intervertebral disc of the lumbar spine during occupational activities.

METHODS

In the proposed method, motion analysis is performed first by using motion capture data, and the results are employed as input data to FE analysis at specific times of interest during motion. In this method, an in-house interface program is used to scale an initial reference FE model to the subject of experiment, and transformed to the corresponding posture at a specific time during motion. The muscle forces and GRF obtained from motion analysis are applied to FE analysis as boundary and loading conditions. For a total of eighteen occupational activities, the IDP, shear stress, and strain of the L4-L5 segment are estimated.

RESULTS

Under each in vivo activity, the predicted IDP was in overall agreement with the available in vivo data. For lifting activities according to lift origin position, the maximum IDP occurred in the far-knee position immediately after lifting. As the lift origin position moved away from the spine, the stresses and strains in the disc increased.

CONCLUSIONS

This new proposed method is expected to allow the estimation of the stresses and deformations in the intervertebral disc during various occupational activities.

Prediction of the influence of vertical whole-body vibration on biomechanics of spinal segments after lumbar interbody fusion surgery

BACKGROUND

Previous studies have shown that for healthy spine, cyclic loading encountered due to whole-body vibration exposure generated higher responses in spinal tissues than static loading. However, how whole-body vibration affects spine biomechanics after interbody fusion surgery is poorly understood. This study aimed at comparing the effects of vibration loading on spinal segments between postsurgical and healthy lumbar spines.

METHODS

A validated finite element model of healthy human lumbosacral spine was modified to simulate interbody fusion at L4-L5 level considering the statuses immediately after surgery (before bony fusion) and after bony fusion. Biomechanical responses at its adjacent levels for the healthy and fusion models to a sinusoidal axial vibration load of ±40 N and the corresponding static axal loads (-40 N and 40 N) were computed using transient dynamic and static analyses, respectively.

FINDINGS

For both healthy and fusion models, vibration amplitudes of the predicted responses were significantly higher than the corresponding changing amplitudes under static loads. Specifically, the increasing effect of vibration load in disc bulge, disc stress and intradiscal pressure at L3-L4 level reached 255.9%, 215.0% and 224.4% for the healthy model, 157.4%, 177.8% and 171.8% for the fusion model (before bony fusion), 141.9%, 152.6% and 160.1% for the fusion model (after bony fusion).

INTERPRETATION

Although whole-body vibration is still more dangerous for the lumbar spine after interbody fusion surgery than static loading, the sensitivity of adjacent segment in postsurgical spine to vibration loading is decreased compared with healthy spine, especially when reaching to bony fusion.

Biomechanical analysis of lumbar interbody fusion supplemented with various posterior stabilization systems

PURPOSE

Biomechanical comparison between rigid and non-rigid posterior stabilization systems following lumbar interbody fusion has been conducted in several studies. However, most of these previous studies mainly focused on investigating biomechanics of adjacent spinal segments or spine stability. The objective of the present study was to compare biomechanical responses of the fusion devices when using different posterior instrumentations.

METHODS

Finite-element model of the intact human lumbar spine (L1-sacrum) was modified to simulate implantation of the fusion cage at L4-L5 level supplemented with different posterior stabilization systems including (i) pedicle screw-based fixation using rigid connecting rods (titanium rods), (ii) pedicle screw-based fixation using flexible connecting rods (PEEK rods) and (iii) dynamic interspinous spacer (DIAM). Stress responses were compared among these various models under bending moments.

RESULTS

The highest and lowest stresses in endplate, fusion cage and bone graft were found at the fused L4-L5 level with DIAM and titanium rod stabilization systems, respectively. When using PEEK rod for the pedicle screw fixation, peak stress in the pedicle screw was lower but the ratio of peak stress in the rods to yield stress of the rod material was higher than using titanium rod.

CONCLUSIONS

Compared with conventional rigid posterior stabilization system, the use of non-rigid stabilization system (i.e., the PEEK rod system and DIAM system) following lumbar interbody fusion might increase the risks of cage subsidence and cage damage, but promote bony fusion due to higher stress in the bone graft. For the pedicle screw-based rod stabilization system, using PEEK rod might reduce the risk of screw breakage but increased breakage risk of the rod itself.

Biomechanical analysis of lumbar spine with interbody fusion surgery and U-shaped lumbar interspinous spacers

Previous research indicates whole-body vibration may lead to low back pain. The aim of this study is assessing the dynamic characteristics of a lumbar spine with Coflex and Coflex-F (commercial implants used as lumbar interspinous spacers) and effect of lumbar interbody fusion surgery. A transient dynamic analysis is performed on three numerical lumbar spine models under the loading condition of a vertical sinusoidal force of ±40 N with a compressive follower preload of 400 N. Also, Coflex-F model with and without interbody fusion surgery is analyzed under the same loading condition. The results show that the maximum value and vibration amplitude of von Mises stress in annulus ground substance (AGS) and intradiscal pressure (IDP) at implanted segment decrease from healthy model to Coflex model, and Coflex-F model. By contrast, for adjacent segments the maximum value of implanted models are larger than that of healthy model. The maximum value of endplates with and without cage are 2.44 MPa and 1.73 MPa (L4 inferior endplate), 1.94 MPa and 1.42 MPa (L5 superior endplate), respectively. The vibration amplitude of Coflex-F model with fusion surgery is smaller than that without fusion surgery. Coflex and Coflex-F not only protect implanted segment but also have a negative effect on adjacent segments. Inserting cage for Coflex-F model can absorb vibration energy at adjacent segments.