Limitations of current in vitro test protocols for investigation of instrumented adjacent segment biomechanics: critical analysis of the literature

Determining a relative total lumbar range of motion to alleviate adjacent segment degeneration after transforaminal lumbar interbody fusion: a finite element analysis

BACKGROUND

A reduction in total lumbar range of motion (ROM) after lumbar fusion may offset the increase in intradiscal pressure (IDP) and facet joint force (FJF) caused by the abnormally increased ROM at adjacent segments. This study aimed to determine a relative total lumbar ROM rather than an ideal adjacent segment ROM to guide postoperative waist activities and further delay adjacent segment degeneration (ASD).

METHODS

An intact L1-S1 finite element model was constructed and validated. Based on this, a surgical model was created to allow the simulation of L4/5 transforaminal lumbar interbody fusion (TLIF). Under the maximum total L1-S1 ROM, the ROM, IDP, and FJF of each adjacent segment between the intact and TLIF models were compared to explore the biomechanical influence of lumbar fusion on adjacent segments. Subsequently, the functional relationship between total L1-S1 ROM and IDP or total L1-S1 ROM and FJF was fitted in the TLIF model to calculate the relative total L1-S1 ROMs without an increase in IDP and FJF.

RESULTS

Compared with those of the intact model, the ROM, IDP, and FJF of the adjacent segments in the TLIF model increased by 12.6-28.9%, 0.1-6.8%, and 0-134.2%, respectively. As the total L1-S1 ROM increased, the IDP and FJF of each adjacent segment increased by varying degrees. The relative total L1-S1 ROMs in the TLIF model were 11.03°, 12.50°, 12.14°, and 9.82° in flexion, extension, lateral bending, and axial rotation, respectively.

CONCLUSIONS

The relative total L1-S1 ROMs after TLIF were determined, which decreased by 19.6-29.3% compared to the preoperative ones. Guiding the patients to perform postoperative waist activities within these specific ROMs, an increase in the IDP and FJF of adjacent segments may be effectively offset, thereby alleviating ASD.

Development of a novel FE model for investigation of interactions of multi-motion segments of the lumbar spine

The spinal anatomy is composed of a series of motion segments (MSs). Although finite element (FE) analysis has been extensively used to investigate the spinal biomechanics with various simplifications of the spinal structures, it is still a challenge to investigate the interactions of different MSs. Anatomical studies have shown that there are major spine ligaments connecting not only single-MS (i.e., two consecutive vertebrae) but also spanning multi-vertebral bones or multi-MSs. However, the effects of the multi-MS spanning ligaments on the spine biomechanics have not been investigated previously. This study developed an FE model of the lumbar spine by simulating the anterior longitudinal ligaments (ALLs) in two portions, one connecting a single-MS and the other spanning two MSs, with varying physiological cross-sectional area (PCSA) ratios of the two portions. The spine biomechanics during extension motion were investigated. The results showed that on average, the constraining forces by the two-MS spanning elements were ∼18% of those of the single-MS ALL elements when the PCSA ratio was 50%, but the two-MS ALL elements also applied compressive forces on the anterior surfaces of the vertebrae. Decreases in intradiscal pressure were also calculated when the two-MS spanning ALL elements were included in the spine model. The multi-MS spanning ligaments were shown to synergistically function with the single-MS elements in spine biomechanics, especially in the interactions of different MSs. The novel lumbar FE model could therefore provide a useful analysis tool for investigation of physiological functions of the spine.

Off-Axis Loading Fixture for Spine Biomechanics: Combined Compression and Bending

The spine is a multi-tissue musculoskeletal system that supports large multi-axial loads and motions during physiological activities. The healthy and pathological biomechanical function of the spine and its subtissues are generally studied using cadaveric specimens that often require multi-axis biomechanical test systems to mimic the complex loading environment of the spine. Unfortunately, an off-the-shelf device can easily exceed 200,000 USD, while a custom device requires extensive time and experience in mechatronics. Our goal was to develop a cost-appropriate compression and bending (flexion-extension and lateral bending) spine testing system that requires little time and minimal technical knowledge. Our solution was an off-axis loading fixture (OLaF) that mounts to an existing uni-axial test frame and requires no additional actuators. OLaF requires little machining, with most components purchased off-the-shelf, and costs less than 10,000 USD. The only external transducer required is a six-axis load cell. Furthermore, OLaF is controlled using the existing uni-axial test frame's software, while the load data is collected using the software included with the six-axis load cell. Here we provide the design rationale for how OLaF develops primary motions and loads and minimizes off-axis secondary constraints, verify the primary kinematics using motion capture, and demonstrate that the system is capable of applying physiologically relevant, noninjurious, axial compression and bending. While OLaF is limited to compression and bending studies it produces repeatable physiologically relevant biomechanics, with high quality data, and minimal startup costs.

Development and Validation of a Novel In Vitro Joint Testing System for Reproduction of In Vivo Dynamic Muscle Force

In vitro biomechanical experiments utilizing cadaveric specimens are one of the most effective methods for rehearsing surgical procedures, testing implants, and guiding postoperative rehabilitation. Applying dynamic physiological muscle force to the specimens is a challenge to reconstructing the environment of bionic mechanics in vivo, which is often ignored in the in vitro experiment. The current work aims to establish a hardware platform and numerical computation methods to reproduce dynamic muscle forces that can be applied to mechanical testing on in vitro specimens. Dynamic muscle loading is simulated through numerical computation, and the inputs of the platform will be derived. Then, the accuracy and robustness of the platform will be evaluated through actual muscle loading tests in vitro. The tests were run on three muscles (gastrocnemius lateralis, the rectus femoris, and the semitendinosus) around the knee joint and the results showed that the platform can accurately reproduce the magnitude of muscle strength (errors range from -6.2% to 1.81%) and changing pattern (goodness-of-fit range coefficient ranges from 0.00 to 0.06) of target muscle forces. The robustness of the platform is mainly manifested in that the platform can still accurately reproduce muscle force after changing the hardware combination. Additionally, the standard deviation of repeated test results is very small (standard ranges of hardware combination 1: 0.34 N~2.79 N vs. hardware combination 2: 0.68 N~2.93 N). Thus, the platform can stably and accurately reproduce muscle forces in vitro, and it has great potential to be applied in the future musculoskeletal loading system.

Biomechanics after spinal decompression and posterior instrumentation

PURPOSE

The aim of this study was to elucidate segmental range of motion (ROM) before and after common decompression and fusion procedures on the lumbar spine.

METHODS

ROM of fourteen fresh-frozen human cadaver lumbar segments (L1/2: 4, L3/4: 5, L5/S1: 5) was evaluated in six loading directions: flexion/extension (FE), lateral bending (LB), lateral shear (LS), anterior shear (AS), axial rotation (AR), and axial compression/distraction (AC). ROM was tested with and without posterior instrumentation under the following conditions: 1) native 2) after unilateral laminotomy, 3) after midline decompression, and 4) after nucleotomy.

RESULTS

Median native ROM was FE 6.8°, LB 5.6°, and AR 1.7°, AS 1.8 mm, LS 1.4 mm, AC 0.3 mm. Unilateral laminotomy significantly increased ROM by 6% (FE), 3% (LB), 12% (AR), 11% (AS), and 8% (LS). Midline decompression significantly increased these numbers to 15%, 5%, 21%, 20%, and 19%, respectively. Nucleotomy further increased ROM in all directions, most substantially in AC of 153%. Pedicle screw fixation led to ROM decreases of 82% in FE, 72% in LB, 42% in AR, 31% in AS, and 17% in LS. In instrumented segments, decompression only irrelevantly affected ROM.

CONCLUSIONS

The amount of posterior decompression significantly impacts ROM of the lumbar spine. The here performed biomechanical study allows creation of a simplified rule of thumb: Increases in segmental ROM of approximately 10%, 20%, and 50% can be expected after unilateral laminotomy, midline decompression, and nucleotomy, respectively. Instrumentation decreases ROM by approximately 80% in bending moments and accompanied decompression procedures only minorly destabilize the instrumentation construct.

Biomechanical analysis of single- and double-level cervical disc arthroplasty using finite element modeling

Recently, many different types of artificial discs have been introduced to persevere the biomechanical behavior of the cervical spine. This study compares the biomechanical behavior of single- and double-level cervical disc arthroplasty, that is "Prestige LP and Mobi-C" on the index and adjacent segment. A three-dimension finite element model of C2-C7 was developed and validated. In single-level prostheses, the Prestige LP or Mobi-C was implanted in the segment C5-C6, while the double-level arthroplasty was integrated at both segments C4-C5 and C5-C6 in the FE model. The intact FE and prosthesis-modified models were constrained from the inferior endplate of the vertebra C7 and applied a compressive load of 73.6 N with a moment load of 1 Nm on the odontoid process of the vertebra C2 to produce flexion/extension, lateral bending, and axial rotation. The prosthesis-modified model's range of motion and intradiscal pressure were determined and compared to the intact model. Also examined the impact of the prostheses on the stress at the bone-implant interface. The range of motion of the implanted segments in both single- and double-levels arthroplasty was increased while that of the adjacent segment of implanted segments decreased. The intradiscal pressure in both levels of arthroplasty was greater than in the intact model. In conclusion, Mobi-C's cervical prostheses could better preserve the normal range of motion and maintain intradiscal pressure than the Prestige LP.

Lumbosacral transitional vertebrae alter the distribution of lumbar mobility–Preliminary results of a radiographic evaluation

Background

Lumbo-sacral transitional vertebrae (LSTV) are one of the most common congenital variances of the spine. They are associated with an increased frequency of degeneration in the cranial adjacent segment. Hypermobility and concomitant increased loads are discussed as a possible reason for segmental degeneration. We therefore examined the lumbar and segmental motion distribution in patients with LSTV with flexion-extension radiographs.

Methods

A retrospective study of 51 patients with osteochondrosis L5/S1 with flexion and extension radiographs was performed. Of these, 17 patients had LSTV and were matched 1:1 for age and sex with patients without LSTV out of the collective of the remaining 34 patients. The lumbar and segmental range of motion (RoM) by segmental lordosis angle and the segmental wedge angle were determined. Normal distribution of parameters was observed by Kolmogorov-Smirnov-test. Parametric data were compared by paired T-test. Non-parametric data were compared by Wilcoxon-rank-sum-test. Correlations were observed using Spearman’s Rank correlation coefficient. A p-value <0.05 was stated as statistically significant.

Results

Patients with LSTV had mean age of 52.2±10.9, control group of 48.9±10.3. Both groups included 7 females and 10 males. Patients with LSTV presented with reduced RoM of the lumbar spine (LSTV 37.3°±19.2°, control 52.1°±20.5°, p = 0.065), however effects were statistically insignificant. LSTV significantly decreased segmental RoM in the transitional segment (LSTV 1.8°±2.7°, control 6.7°±6.0°, p = 0.003). Lumbar motion distribution differed significantly; while RoM was decreased in the transitional segment, (LSTV 5.7%, control 16.2%, p = 0.002), the distribution of lumbar motion to the cranial adjacent segment was increased (LSTV 30.7%, control 21.6%, p = 0.007).

Conclusion

Patients with LSTV show a reduced RoM in the transitional segment and a significantly increased motion distribution to the cranial adjacent segment in flexion-extension radiographs. The increased proportion of mobility in the cranial adjacent segment possibly explain the higher rates of degeneration within the segment.

In Silico Meta-Analysis of Boundary Conditions for Experimental Tests on the Lumbar Spine

The study of the spine range of motion under given external load has been the object of many studies in literature, finalised to a better understanding of the spine biomechanics, its physiology, eventual pathologic conditions and possible rehabilitation strategies. However, the huge amount of experimental work performed so far cannot be straightforwardly analysed due to significant differences among loading set-ups. This work performs a meta-analysis of various boundary conditions in literature, focusing on the flexion/extension behaviour of the lumbar spine. The comparison among range of motions is performed virtually through a validated multibody model. Results clearly illustrated the effect of various boundary conditions which can be met in literature, so justifying differences of biomechanical behaviours reported by authors implementing different set-up: for example, a higher value of the follower load can indeed result in a stiffer behaviour; the application of force producing spurious moments results in an apparently more deformable behaviour, however the respective effects change at various segments along the spine due to its natural curvature. These outcomes are reported not only in qualitative, but also in quantitative terms. The numerical approach here followed to perform the meta-analysis is original and it proved to be effective thanks to the bypass of the natural variability among specimens which might completely or partially hinder the effect of some boundary conditions. In addition, it can provide very complete information since the behaviour of each functional spinal unit can be recorded. On the whole, the work provided an extensive review of lumbar spine loading in flexion/extension.

Changes of adjacent segment biomechanics after anterior cervical interbody fusion with different profile design plate: single- versus double-level

Low-profile angle-stable spacer Zero-P is claimed to reduce the morbidity associated with traditional plate and cage construct (PCC). Both Zero-P and PCC could achieve comparable mid- and long-term clinical and radiological outcomes in anterior cervical discectomy and fusion (ACDF). It is not clear whether Zero-P can reduce the incidence of adjacent segment degeneration (ASD), especially in multi-segmental fusion. This study aimed to test the effect of fusion level with Zero-P versus with PCC on adjacent-segment biomechanics in ACDF. A three-dimensional finite element (FE) model of an intact C2-T1 segment was built and validated. Six single- or double-level instrumented conditions were modeled from this intact FE model using Zero-P or the standard PCC. The biomechanical responses of adjacent segments at the cephalad and caudal levels of the operation level were assessed in terms of range of motion (ROM), stresses in the endplate and disc, loads in the facets. When comparing the increase of adjacent-segment motion in single-level PCC fusion versus Zero-P fusion, a significantly larger increase was found in double-level fusion condition. The fold changes of PCC versus Zero-P of intradiscal and endplate stress, and facet load at adjacent levels in the double-level fusion spine were significantly larger than that in the single-level fusion spine during the sagittal, the transverse, and the frontal plane motion. The increased value of biomechanical features was greater at above segment than that at below. The fold changes of PCC versus Zero-P at adjacent segment were most notable in flexion and extension movement. Low-profile device could decrease adjacent segment biomechanical burden compared to traditional PCC in ACDF, especially in double-level surgery. Zero-P could be a good alternative for traditional PCC in ACDF. Further clinical/in vivo studies will be necessary to explore the approaches selected for this study is warranted.

Novel force-displacement control passive finite element models of the spine to simulate intact and pathological conditions; comparisons with traditional passive and detailed musculoskeletal models

Passive finite element (FE) models of the spine are commonly used to simulate intact and various pre- and postoperative pathological conditions. Being devoid of muscles, these traditional models are driven by simplistic loading scenarios, e.g., a constant moment and compressive follower load (FL) that do not properly mimic the complex in vivo loading condition under muscle exertions. We aim to develop novel passive FE models that are driven by more realistic yet simple loading scenarios, i.e., in vivo vertebral rotations and pathological-condition dependent FLs (estimated based on detailed musculoskeletal finite element (MS-FE) models). In these novel force-displacement control FE models, unlike the traditional passive FE models, FLs vary not only at different spine segments (T12-S1) but between intact, pre- and postoperative conditions. Intact, preoperative degenerated, and postoperative fused conditions at the L4-L5 segment for five static in vivo activities in upright and flexed postures were simulated by the traditional passive FE, novel force-displacement control FE, and gold-standard detailed MS-FE spine models. Our findings indicate that, when compared to the MS-FE models, the force-displacement control passive FE models could accurately predict the magnitude of disc compression force, intradiscal pressure, annulus maximal von Mises stress, and vector sum of all ligament forces at adjacent segments (L3-L4 and L5-S1) but failed to predict disc shear and facet joint forces. In this regard, the force-displacement control passive FE models were much more accurate than the traditional passive FE models. Clinical recommendations made based on traditional passive FE models should, therefore, be interpreted with caution.

Adjacent segments biomechanics following lumbar fusion surgery: a musculoskeletal finite element model study

PURPOSE

This study exploits a novel musculoskeletal finite element (MS-FE) spine model to evaluate the post-fusion (L4-L5) alterations in adjacent segment kinetics.

METHODS

Unlike the existing MS models with idealized representation of spinal joints, this model predicts stress/strain distributions in all passive tissues while organically coupled to a MS model. This generic (in terms of musculature and material properties) model uses population-based in vivo vertebral sagittal rotations, gravity loads, and an optimization algorithm to calculate muscle forces. Simulations represent individuals with an intact L4-L5, a preoperative severely degenerated L4-L5 (by reducing the disc height by ~ 60% and removing the nucleus incompressibility), and a postoperative fused L4-L5 segment with either a fixed or an altered lumbopelvic rhythm with respect to the intact condition (based on clinical observations). Changes in spine kinematics and back muscle cross-sectional areas (due to intraoperative injuries) are considered based on in vivo data while simulating three activities in upright/flexed postures.

RESULTS

Postoperative changes in some adjacent segment kinetics were found considerable (i.e., larger than 25%) that depended on the postoperative lumbopelvic kinematics and preoperative L4-L5 disc condition. Postoperative alterations in adjacent disc shear, facet/ligament forces, and annulus stresses/strains were greater (> 25%) than those found in intradiscal pressure and compression (< 25%). Kinetics of the lower (L5-S1) and upper (L3-L4) adjacent segments were altered to different degrees.

CONCLUSION

Alterations in segmental rotations mainly affected adjacent disc shear forces, facet/ligament forces, and annulus/collagen fibers stresses/strains. An altered lumbopelvic rhythm (increased pelvis rotation) tends to mitigate some of these surgically induced changes.

Biomechanical changes at the adjacent segments induced by a lordotic porous interbody fusion cage

Biomechanical changes at the adjacent segments after interbody fusion are common instigators of adjacent segment degeneration (ASD). This study aims to investigate how the presence of a lordotic porous cage affects the biomechanical performance of the adjacent segments. A finite element model (FEM) of a lumbar spine implanted with a lordotic cage at L3-L4 was validated by in-vitro testing. The stress distribution on the cage and range of motion (ROM) of L3-L4 were used to assess the stability of the implant. Three angles of cage (0° = non-restoration, 7° = normal restoration and 11° = over-restoration) were modelled with different porosities (0%, 30% and 60%) and evaluated in the motions of flexion, extension, lateral bending and rotation. The ROM, intervertebral disc pressure (IDP) and facet joint force (FJF) were used to evaluate biomechanical changes at the adjacent segments in each model. The results indicated that porous cages produced more uniform stress distribution, but cage porosity did not influence the ROM, IDP and FJF at L2-L3 and L4-L5. Increasing the cage lordotic angle acted to decrease the ROM and IDP, and increase the FJF of L4-L5, but did not alter the ROM of L2-L3. In conclusion, changes in ROM, IDP and FJF at the adjacent segments were mainly influenced by the lordotic angle of the cage and not by the porosity. A larger angle of lordotic cage was shown to reduce the ROM and IDP, and increase the FJF of the lower segment (L4-L5), but had little effect on the ROM of the upper segment (L2-L3).

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.

Lumbar Stabilization with DSS-HPS® System: Radiological Outcomes and Correlation with Adjacent Segment Degeneration

Arthrodesis has always been considered the main treatment of degenerative lumbar disease. Adjacent segment degeneration is one of the major topics related to fusion surgery. Non-fusion surgery may prevent this because of the protective effect of persisting segmental motion. The aims of the study were (1) to describe the radiological outcomes in the adjacent vertebral segment after lumbar stabilization with DSS-HPS^® system and (2) to verify the hypothesis that this system prevents the degeneration of the adjacent segment. This is a retrospective monocentric analysis of twenty-seven patients affected by degenerative lumbar disease underwent spinal hybrid stabilization with the DSS-HPS^® system between January 2016 and January 2019. All patients completed 1-year radiological follow-up. Preoperative X-rays and magnetic resonance images, as well as postoperative radiographs at 1, 6 and 12 months, were evaluated by one single observer. Pre- and post-operative anterior and posterior disc height at the dynamic (DL) and adjacent level (AL) were measured; segmental angle (SA) of the dynamized level were measured. There was a statistically significant decrease of both anterior (p = 0.0003 for the DL, p = 0.036 for the AL) and posterior disc height (p = 0.00000 for the DL, p = 0.00032 for the AL); there were a statistically significant variations of the segmental angle (p = 0.00000). Eleven cases (40.7%) of radiological progression of disc degeneration were found. The DSS-HPS^® system does not seem to reduce progression of lumbar disc degeneration in a radiologic evaluation, both in the dynamized and adjacent level.

Biomechanical effects of uncinate process excision in cervical disc arthroplasty

BACKGROUND

Studies on the role of uncinate process have been limited to responses of the intact spine and patient's outcomes, and procedures to perform the excision. The aim of this study was to determine the role of uncinate process on the biomechanical response at the index and adjacent levels in three artificial discs used in cervical disc arthroplasty.

METHODS

A validated finite element model of cervical spine was used. Flexion, extension, and lateral moments and follower load were applied to Bryan, Mobi-C, and Prestige LP artificial discs at C5-C6 level with and without uncinate process. Ranges of motion at index level and adjacent caudal and cranial segments, intradiscal pressures at adjacent segments, and facet loads at index level and adjacent segments were obtained. Data were normalized with respect to the preservation of uncinate process.

FINDINGS

Uncinate process removal increased motions up to 27% at index and decreased up to 10% at adjacent levels, decreased disc pressures up to 14% at adjacent segments, decreased facet loads at adjacent segments up to 14%, while at index level, change in loads depended on mode and arthroplasty, with Mobi-C responding with up to 51% increase and Bryan disc up to 11% decrease, while Prestige LP increased loads by 17% in extension and decreased by 9%% in lateral bending.

INTERPRETATION

As surgical selection is based on morphology and surgeon's experience, the present computational findings provide quantitative information for an optimal choice of the device and procedure, while further studies (in vitro/clinical) would be required.

Biomechanical effects of lumbar fusion surgery on adjacent segments using musculoskeletal models of the intact, degenerated and fused spine

Adjacent segment disorders are prevalent in patients following a spinal fusion surgery. Postoperative alterations in the adjacent segment biomechanics play a role in the etiology of these conditions. While experimental approaches fail to directly quantify spinal loads, previous modeling studies have numerous shortcomings when simulating the complex structures of the spine and the pre/postoperative mechanobiology of the patient. The biomechanical effects of the L4–L5 fusion surgery on muscle forces and adjacent segment kinetics (compression, shear, and moment) were investigated using a validated musculoskeletal model. The model was driven by in vivo kinematics for both preoperative (intact or severely degenerated L4–L5) and postoperative conditions while accounting for muscle atrophies. Results indicated marked changes in the kinetics of adjacent L3–L4 and L5–S1 segments (e.g., by up to 115% and 73% in shear loads and passive moments, respectively) that depended on the preoperative L4–L5 disc condition, postoperative lumbopelvic kinematics and, to a lesser extent, postoperative changes in the L4–L5 segmental lordosis and muscle injuries. Upper adjacent segment was more affected post-fusion than the lower one. While these findings identify risk factors for adjacent segment disorders, they indicate that surgical and postoperative rehabilitation interventions should focus on the preservation/restoration of patient’s normal segmental kinematics.

Biomechanical analysis of lumbar fusion with proximal interspinous process device implantation

Lumbar spinal fusion may cause adjacent segment degeneration (ASD) in the long term. Recently, inserting an interspinous process device (IPD) proximal to the fusion has been proposed to prevent ASD. The aim of this study was to investigate the biomechanics of lumbar fusion with proximal IPD implantation (LFPI) under both static loads and whole body vibration (WBV). A previously validated finite element (FE) model of the L1-5 lumbar spine was modified to simulate L4-5 fusion. Three different IPDs (Coflex-F, Wallis and DIAM) were inserted at the L3-4 segment of the fusion model to construct the LFPI models. The intact and surgical FE models were analyzed under static loads and WBV, respectively. Under static loading conditions, LFPI decreased range of motion (ROM) and intradiscal pressure (IDP) at the transition segment L3-4 compared with the fusion case. At the segment (L2-3) adjacent to the transition level, LFPI induced higher motion and IDP than rigid fusion. Under WBV, vibration amplitudes of the L3-4 IDP and L4-5 facet joint force (FJF) decreased by more than 54.3% after surgery. The LFPI model with the DIAM system offered the most comparable biomechanics to the intact model under static loads, and decreased the dynamic responses of the L4-5 FJF under WBV. The LFPI model with the Wallis and Coflex-F systems could stabilize the transition segment, and decrease dynamic responses of the L3-4 IDP. The DIAM system may be more suitable in LFPI.

Influence of structural and material property uncertainties on biomechanics of intervertebral discs - Implications for disc tissue engineering

This study investigated how variations of structural and material properties of human intervertebral discs (IVDs) affect the biomechanical responses of the IVDs under simulated physiological loading conditions using a stochastic finite element (SFE) model. An SFE method, which combined an anatomic FE model of human lumbar L3-4 segment and probabilistic analysis of its structural and material properties, was used to generate a dataset of 500 random disc samples with varying structural and material properties. The sensitivity of the biomechanical responses, including intervertebral displacements/rotations, intradiscal pressures (IDP), fiber stresses and matrix strains of annulus fibrosus (AF), were systematically quantified under various physiological loading conditions, including a 500N compression and 7.5Nm moments in the 3 primary rotations. Significant variations of the IDPs, IVD displacements/rotations, and stress/strain distributions were found using the dataset of 500 ramdom disc samples. Under all the loading conditions, the IDPs were positively correlated with the Poisson's ratio of the NP (r = 0.46 to 0.75, p = 0.004-0.001) and negatively with the Young's modulus of the annulus matrix (r = -0.48 to -0.65, p = 0.003-0.001). The primary intervertebral rotations were significantly affected by the Young's modulus of the annulus matrix (r = -0.44 to -0.71, p = 0.001-0.032) and the orientations of the annular fibers (r = -0.45 to -0.69, p = 0.001-0.029). The heterogeneity of structures and material properties of the IVD had distinct effects on the biomechanical performances of the IVD. These data could help improve our understanding of the intrinsic biomechanics of the IVD and provide references for optimal design of tissue engineered discs by controlling structural and material properties of the disc components.

Prediction of biomechanical responses of human lumbar discs - a stochastic finite element model analysis

BACKGROUND

Accurate biomechanical investigation of human intervertebral discs (IVDs) is difficult because of their complicated structural and material features.

AIM

To investigate probabilistic distributions of the biomechanical responses of the IVD by considering varying nonlinear structural and material properties using a stochastic finite element (FE) model.

METHODS

A FE model of a L3-4 disc was reconstructed, including the nucleus pulposus (NP), annular matrix and fibers. A Monte Carlo method was used to randomly generate 500 sets of the nonlinear material properties and fiber orientations of the disc that were implemented into the FE model. The FE model was analyzed under seven loading conditions: a 500 N compressive force, a 7.5Nm moment simulating flexion, extension, left-right lateral bending, and left-right axial rotation, respectively. The distributions of the ranges of motion (ROMs), intradiscal pressures (IDP), fiber stresses and matrix strains of the disc were analyzed.

RESULTS

Under the compressive load, the displacement varied between 0.29 mm and 0.76 mm. Under the 7.5Nm moment, the ROMs varied between 3.0° and 6.0° in primary rotations. The IDPs varied within 0.3 MPa under all the loading conditions. The maximal fiber stress (3.22 ± 0.64 MPa) and matrix strain (0.27 ± 0.12%) were observed under the flexion and extension moments, respectively.

CONCLUSION

The IVD biomechanics could be dramatically affected by the structural and material parameters used to construct the FE model. The stochastic FE model that includes the probabilistic distributions of the structural and material parameters provides a useful approach to analyze the statistical ranges of the biomechanical responses of the IVDs.

Assessing the biofidelity of in vitro biomechanical testing of the human cervical spine

In vitro biomechanical studies of the osteoligamentous spine are widely used to characterize normal biomechanics, identify injury mechanisms, and assess the effects of degeneration and surgical instrumentation on spine mechanics. The objective of this study was to determine how well four standards in vitro loading paradigms replicate in vivo kinematics with regards to the instantaneous center of rotation and arthrokinematics in relation to disc deformation. In vivo data were previously collected from 20 asymptomatic participants (45.5 ± 5.8 years) who performed full range of motion neck flexion-extension (FE) within a biplane x-ray system. Intervertebral kinematics were determined with sub-millimeter precision using a validated model-based tracking process. Ten cadaveric spines (51.8 ± 7.3 years) were tested in FE within a robotic testing system. Each specimen was tested under four loading conditions: pure moment, axial loading, follower loading, and combined loading. The in vivo and in vitro bone motion data were directly compared. The average in vitro instant center of rotation was significantly more anterior in all four loading paradigms for all levels. In general, the anterior and posterior disc heights were larger in the in vitro models than in vivo. However, after adjusting for gender, the observed differences in disc height were not statistically significant. This data suggests that in vitro biomechanical testing alone may fail to replicate in vivo conditions, with significant implications for novel motion preservation devices such as cervical disc arthroplasty implants.

A comprehensive approach for the validation of lumbar spine finite element models investigating post-fusion adjacent segment effects

Spinal fusion surgery is usually followed by accelerated degenerative changes in the unfused segments above and below the treated segment(s), i.e., adjacent segment disease (ASD). While a number of risk factors for ASD have been suggested, its exact pathogenesis remains to be identified. Finite element (FE) models are indispensable tools to investigate mechanical effects of fusion surgeries on post-fusion changes in the adjacent segment kinematics and kinetics. Existing modeling studies validate only their intact FE model against in vitro data and subsequently simulate post-fusion in vivo conditions. The present study provides a novel approach for the comprehensive validation of a lumbar (T12-S1) FE model in post-fusion conditions. Sixteen simulated fusion surgeries, performed on cadaveric specimens using various testing and loading conditions, were modeled by this FE model. Predictions for adjacent segment range of motion (RoM) and intradiscal pressure (IDP) were compared with those obtained from the corresponding in vitro tests. Overall, 70% of the predicted adjacent segment RoMs were within the range of in vitro data for both intact and post-fusion conditions. Correlation (r) values between model and in vitro findings for the adjacent segment RoMs were positive and greater than 0.84. Most of the predicted IDPs were, however, out of the narrow range of in vitro IDPs at the adjacent segments but with great positive correlations (r ≥ 0.89). FE modeling studies investigating the effect of fusion surgery on in vivo adjacent segment biomechanics are encouraged to use post-surgery in vitro data to validate their FE model.

A minimum 8-year follow-up comparative study of decompression and coflex stabilization with decompression and fusion

The current study aimed to compare the outcomes of decompression and interlaminar stabilisation with those of decompression and fusion for the treatment of lumbar degenerative disease (LDD) at a minimum 8-year follow-up. The current study also aimed to analyse the risk factors of radiographic adjacent segment degeneration (ASD). A total of 82 consecutive patients with LDD who underwent surgery between June 2007 and February 2011 were retrospectively reviewed. Of these patients, 39 underwent decompression and Coflex interspinous stabilisation (Coflex group) and 43 underwent decompression and posterior lumbar interbody fusion (PLIF) (PLIF group). All patients had a minimum of 8-years of follow-up data. Radiographic and clinical outcomes were compared between the groups, and the risk factors of developing radiographic ASD were also evaluated. The Oswestry disability index and visual analogue scale leg and back pain scores of both groups significantly improved compared with the baseline (all P<0.05), and no difference were indicated between the two groups at each follow-up time point (P>0.05). The Coflex group exhibited preserved mobility (P<0.001), which was associated with a decreased amount of blood loss (P<0.001), shorter duration of surgery (P=0.001), shorter duration of hospital stay and a lower incidence of ASD (12.8 vs. 32.56%; P=0.040) compared with the fusion group. The current study indicated that coflex and fusion technologies are safe and effective for the treatment of LDD, based on long-term follow-up data. However, Coflex interspinous stabilisation was revealed to reduce ASD incidence. Under strict indications, Coflex interspinous stabilisation is an effective and safe treatment method.

Loading of the lumbar spine during transition from standing to sitting: effect of fusion versus motion preservation at L4-L5 and L5-S1

BACKGROUND CONTEXT

Transition from standing to sitting significantly decreases lumbar lordosis with the greatest lordosis-loss occurring at L4-S1. Fusing L4-S1 eliminates motion and thus the proximal mobile segments maybe recruited during transition from standing to sitting to compensate for the loss of L4-S1 mobility. This may subject proximal segments to supra-physiologic flexion loading.

PURPOSE

Assess effects of instrumented fusion versus motion preservation at L4-L5 and L5-S1 on lumbar spine loads and proximal segment motions during transition from standing to sitting.

STUDY DESIGN

Biomechanical study using human thoracolumbar spine specimens.

METHODS

A novel laboratory model was used to simulate lumbosacral alignment changes caused by a person's transition from standing to sitting in eight T10-sacrum spine specimens. The sacrum was tilted in the sagittal plane while constraining anterior-posterior translation of T10. Continuous loading-data and segmental motion-data were collected over a range of sacral slope values, which represented transition from standing to different sitting postures. We compared different constructs involving fusions and motion preserving prostheses across L4-S1.

RESULTS

After L4-S1 fusion, the sacrum could not be tilted as far posteriorly compared to the intact spine for the same applied moment (p<.001). For the same reduction in sacral slope, L4-S1 fusion induced 2.9 times the flexion moment in the lumbar spine and required 2.4 times the flexion motion of the proximal segments as the intact condition (p<.001). Conversely, motion preservation at L4-S1 restored lumbar spine loads and proximal segment motions to intact specimen levels during transition from standing to sitting.

CONCLUSIONS

In general, sitting requires lower lumbar segments to undergo flexion, thereby increasing load on the lumbar disks. L4-S1 fusion induced greater moments and increased flexion of proximal segments to attain a comparable seated posture. Motion preservation using a total joint replacement prosthesis at L4-S1 restored the lumbar spine loads and proximal segment motion to intact specimen levels during transition from standing to sitting.

CLINICAL SIGNIFICANCE

After L4-S1 fusion, increased proximal segment loading during sitting may cause discomfort in some patients and may lead to junctional breakdown over time. Preserving motion at L4-S1 may improve patient comfort and function during activities of daily living, and potentially decrease the need for adjacent level surgery.

Biomechanical models: key considerations in study design

This manuscript summarizes presentations of a symposium on key considerations in design of biomechanical models at the 2019 Basic Science Focus Forum of the Orthopaedic Trauma Association. The first section outlines the most important characteristics of a high-quality biomechanical study. The second section considers choices associated with designing experiments using finite element modeling versus synthetic bones versus human specimens. The third section discusses appropriate selection of experimental protocols and finite element analyses. The fourth section considers the pros and cons of use of biomechanical research for implant design. Finally, the fifth section examines how results from biomechanical studies can be used when clinical evidence is lacking or contradictory. When taken together, these presentations emphasize the critical importance of biomechanical research and the need to carefully consider and optimize models when designing a biomechanical study.

The impact of different artificial disc heights during total cervical disc replacement: an in vitro biomechanical study

BACKGROUND

The principles of choosing an appropriate implant height remain controversial in total cervical disc replacement (TDR). By performing an in vitro biomechanical study and exploring the biomechanical impact of implant height on facet joint and motion function, the study aimed to offer valid proposals regarding implant height selection during TDR.

METHODS

A total of 6 fresh-frozen male cadaveric cervical spines (C2-C7) with 5 mm intervertebral disc height at C5/6 level were enrolled in the study. Specimens with the intact condition and with different height artificial discs were tested. Facet joint pressures and range of motion under each condition were recorded using a specialized machine.

RESULTS

The artificial disc heights that were involved in this study were 5 mm, 6 mm, and 7 mm. The range of motion decreased along with the increment of implant height, while facet joint pressure showed an opposite trend. Specimens with a 5 mm implant height could provide a similar range of motion (11.8° vs. 12.2° in flexion-extension, 8.7° vs. 9.0° in rotation, 7.9° vs. 8.2° in lateral bending) and facet joint pressure (27.8 psi vs. 25.2 psi in flexion, 59.7 psi vs. 58.9 psi in extension, 24.0 psi vs. 22.7 psi in rotation, 32.0 psi vs. 28.8 psi in lateral bending) compared with intact specimens. Facet joint pressure of specimens with 6 mm implant height (≥ 1 mm in height) increased during flexion at the C5-6 segment (30.4 psi vs. 25.2 psi, P = 0.076). However, specimens with 7 mm implant height (≥ 2 mm in height) showed a significant reduction in motion (9.5° vs. 12.2° in flexion-extension, P < 0.001) and increment of facet joint pressure at C5-6 segment (44.6 psi vs. 25.2 psi in flexion, 90.3 psi vs. 58.9 psi in extension, P < 0.0001) and adjacent segments.

CONCLUSIONS

This study suggested that an appropriate artificial disc height can achieve near-normal biomechanical properties and is recommended. We should be very cautious when using artificial discs ≥ 1 mm in height compared to normal. However, implants ≥ 2 mm in height compared to normal significantly increased the facet joint pressure and decreased the range of motion; therefore, it should not be used in clinical practice.

Effect of different designs of interspinous process devices on the instrumented and adjacent levels after double-level lumbar decompression surgery: A finite element analysis

Recently, various designs and material manufactured interspinous process devices (IPDs) are on the market in managing symptomatic lumbar spinal stenosis (LSS). However, atraumatic fracture of the intervening spinous process has been reported in patients, particularly, double or multiple level lumbar decompression surgery with IPDs. This study aimed to biomechanically investigate the effects of few commercial IPDs, namely DIAM^TM, Coflex^TM, and M-PEEK, which were implanted into the L2-3, L3-4 double-level lumbar spinal processes. A validated finite element model of musculoskeletal intact lumbar spinal column was modified to accommodate the numerical analysis of different implants. The range of motion (ROM) between each vertebra, stiffness of the implanted level, intra stress on the intervertebral discs and facet joints, and the contact forces on spinous processes were compared. Among the three implants, the Coflex system showed the largest ROM restriction in extension and caused the highest stress over the disc annulus at the adjacent levels, as well as the sandwich phenomenon on the spinous process at the instrumented levels. Further, the DIAM device provided a superior loading-sharing between the two bridge supports, and the M-PEEK system offered a superior load-sharing from the superior spinous process to the lower pedicle screw. The limited motion at the instrumented segments were compensated by the upper and lower adjacent functional units, however, this increasing ROM and stress would accelerate the degeneration of un-instrumented segments.

Meta-analyses comparing spine simulators with cadavers and finite element models by analysing range-of-motion data before and after lumbar total disc replacement

•

Data from two different in vitro test methods for the same evaluation subjects were compared. It was investigated whether spine simulators with real human cadavers (SSCs) and finite element models (FEMs) provide the same data exemplarily for range of motion (ROM) before and after insertion of motion-retaining devices.

•

Only fifty-nine percent of SSC meta-analyses show restored ROM after insertion of the device compared to the intact spinal segment. In FEM meta-analyses, ROM is restored in ninety percent.

•

Ten percent of ROM analyses show significantly different data between SSCs and FEMs.

•

With regard to the included studies, data generated by SSCs and FEMs cannot be used unrestricted as alternative and complementary data.

•

Our analysis provides a new approach to compare data from associated test methods.

Background

Range-of-motion (ROM) data generated by the in vitro test methods of spine simulators with cadavers (SSCs) and finite element models (FEMs) are used alternatively and complementarily for in vitro evaluations.

Aim of Review

Our purpose is to compare exemplary segmental ROM data from SSCs and FEMs before and after ball-and-socket total disc replacement (bsTDR) to determine whether the two test methods provide the same data for the same evaluation subjects.

Key Scientific Concepts of Review

We performed 70 meta-analyses (MAs) and 20 additional comparative analyses based on data from 21 SSC studies used for 39 MAs and 16 FEM studies used for 31 MAs. Only fifty-nine percent (n = 23/39) of SSC MAs show a restored ROM after bsTDR, whereas in FEM MAs, the ROM is restored in ninety percent (n = 28/31). Among the analyses comparing data from the same spinal segments, motion directions and bsTDR, SSC and FEM data are significantly different in ten percent (n = 2/20). According to our results, data generated by SSCs and FEMs cannot be used as alternative and complementary data without restriction. The quality of the evaluation methods itself as well as potential technical reasons for the discrepant results were not our evaluation target. Further SSC and FEM data should be compared using the same approach.

Risk of adjacent segment disease after 'topping-off' multi-level lumbar fusions with posterior dynamic stabilisers: an observational cohort study

PURPOSE

To determine whether 'topping-off' lumbar fusions, using posterior dynamic stabilising devices (PDSs) with specific biomechanical parameters, reduces the risk of adjacent segment disease (ASD).

METHODS

Survival analysis of two non-randomised cohorts, with or without 'topping-off' (T/O or NoT/O), compared the risk of further surgery for ASD following multi-level posterior lumbar interbody fusion (PLIF). The study sample comprised consecutive patients, aged 55 + years, with degenerative pathology at 2, 3 or 4 levels. The NoT/O cohort underwent surgery between August 1993 and September 2019 (n = 425) and the T/O cohort between September 2011 and September 2019 (n = 146). Comparison of ASD risk between cohorts used Cox proportional hazards (CPH) modelling and Kaplan-Meier survivorship analysis.

RESULTS

Analysis was completed on 571 operations across 507 patients. Median follow-up was 63 months (range 0.3-196) and 37 months (range 1.7-98) for the NoT/O and T/O cohorts, respectively. Of 423 patients, 125 (29.6%) patients in the NoT/O cohort underwent further surgery for ASD and 16/145 (11.03%) in the T/O cohort. The hazard ratio (T/O: NoT/O) from the CPH model was 0.42 (95% CL: 0.24-0.74, P = 0.003). Mean annual incidence across the first 5 years was 5.0% in the NoT/O cohort compared with 2.8% in the T/O cohort (P = 0.029). No patient required surgery or developed ASD at a 'topped-off' level. Two patients developed asymptomatic pedicle screw loosening at the level of the PDS device. PROMs were similar between cohorts.

CONCLUSION

This large, non-randomised, observational study found an approximately 60% reduction in further surgery for ASD with the use of the PDS to 'top-off' PLIF fusions. PDS device-related complications were very low.

Investigation of Alterations in the Lumbar Disc Biomechanics at the Adjacent Segments After Spinal Fusion Using a Combined In Vivo and In Silico Approach

The development of adjacent segment degeneration (ASD) is a major concern after lumbar spinal fusion surgery, but the causative mechanisms remain unclear. This study used a combined in vivo and in silico method to investigate the changes of anatomical dimensions and biomechanical responses of the adjacent segment (L3-4) after spinal fusion (L4-S1) in five patients under weight-bearing upright standing conditions. The in vivo adjacent disc height changes before and after fusion were measured using a dual fluoroscopic imaging system (DFIS), and the measured in vivo intervertebral positions and orientations were used as displacement boundary conditions of the patient-specific three-dimensional (3D) finite element (FE) disc models to simulate the biomechanical responses of adjacent discs to fusion of the diseased segments. Our data (represented by medians and 95% confidence intervals) showed that a significant decrease by - 0.8 (- 1.2, - 0.4) mm (p < 0.05) in the adjacent disc heights occurred at the posterior region after fusion. The significant increases in disc tissue strains and stresses, 0.32 (0.21, 0.43) mm/mm (p < 0.05) and 1.70 (1.07, 3.60) MPa (p < 0.05), respectively, after fusion were found in the posterolateral portions of the outermost annular lamella. The intradiscal pressure of the adjacent disc was significantly increased by 0.29 (0.13, 0.47) MPa after fusion (p < 0.05). This study demonstrated that fusion could cause alterations in adjacent disc biomechanics, and the combined in vivo and in silico method could be a valuable tool for the quantitative assessment of ASD after fusion.

Dynamic foraminal dimensions during neck motion 6.5 years after fusion and artificial disc replacement

Objective

To compare changes in foraminal motion at two time points post-surgery between artificial disc replacement (ADR) and anterior cervical discectomy and fusion (ACDF).

Methods

Eight ACDF and 6 ADR patients (all single-level C5-6) were tested at 2 years (T1) and 6.5 years (T2) post-surgery. The minimum foraminal height (FH.Min) and width (FW.Min) achieved during neck axial rotation and extension, and the range of these dimensions during motion (FH.Rn and FW.Rn, respectively) were measured using a biplane dynamic x-ray system, CT imaging and model-based tracking while patients performed neck axial rotation and extension tasks. Two-way mixed ANOVA was employed for analysis.

Results

In neck extension, significant interactions were found between year post-surgery and type of surgery for FW.Rn at C5-6 (p<0.006) and C6-7 (p<0.005), and for FH.Rn at C6-7 (p<0.01). Post-hoc analysis indicated decreases over time in FW.Rn for ACDF (p<0.01) and increases in FH.Rn for ADR (p<0.03) at the C6-7 adjacent level. At index level, FW.Rn was comparable between ACDF and ADR at T1, but was smaller for ACDF than for ADR at T2 (p<0.002). In axial rotation, differences were found between T1 and T2 but did not depend on type of surgery (p>0.7).

Conclusions

Changes were observed in the range of foraminal geometry at adjacent levels from 2 years to 6.5 years post-surgery that were different between ACDF and ADR for neck extension. These changes are contrary to the notion that motion at adjacent levels continue to increase following ACDF as compared to ADR over the long term.

ICR in human cadaveric specimens: An essential parameter to consider in a new lumbar disc prosthesis design

STUDY DESIGN

Biomechanical study in cadaveric specimens.

BACKGROUND

The commercially available lumbar disc prostheses do not reproduce the intact disc's Instantaneous centre of Rotation (ICR), thus inducing an overload on adjacent anatomical structures, promoting secondary degeneration.

AIM

To examine biomechanical testing of cadaveric lumbar spine specimens in order to evaluate and define the ICR of intact lumbar discs.

MATERIAL AND METHODS

Twelve cold preserved fresh human cadaveric lumbosacral spine specimens were subjected to computerized tomography (CT), magnetic resonance imaging (MRI) and biomechanical testing. Kinematic studies were performed to analyse range of movements in order to determine ICR.

RESULTS

Flexoextension and lateral bending tests showed a positive linear correlation between the angle rotated and the displacement of the ICR in different axes.

DISCUSSION

ICR has not been taken into account in any of the available literature regarding lumbar disc prosthesis. Considering our results, neither the actual ball-and-socket nor the withdrawn elastomeric nucleus models fit the biomechanics of the lumbar spine, which could at least in part explain the failure rates of the implants in terms of postoperative failed back syndrome (low back pain). It is reasonable to consider then that an implant should also adapt the equations of the movement of the intact ICR of the joint to the post-surgical ICR.

CONCLUSIONS

This is the first cadaveric study on the ICR of the human lumbar spine. We have shown that it is feasible to calculate and consider this parameter in order to design future prosthesis with improved clinical and biomechanical characteristics.

Multiobjective Design Optimization of a Biconcave Mobile-Bearing Lumbar Total Artificial Disk Considering Spinal Kinematics, Facet Joint Loading, and Metal-on-Polyethylene Contact Mechanics

Total disk arthroplasty (TDA) using an artificial disk (AD) is an attractive surgical technique for the treatment of spinal disorders, since it can maintain or restore spinal motion (unlike interbody fusion). However, adverse surgical outcomes of contemporary lumbar TDAs have been reported. We previously proposed a new mobile-bearing AD design concept featuring a biconcave ultrahigh-molecular-weight polyethylene (UHMWPE) mobile core. The objective of this study was to develop an artificial neural network (NN) based multiobjective optimization framework to refine the biconcave-core AD design considering multiple TDA performance metrics, simultaneously. We hypothesized that there is a tradeoff relationship between the performance metrics in terms of range of motion (ROM), facet joint force (FJF), and polyethylene contact pressure (PCP). By searching the resulting three-dimensional (3D) Pareto frontier after multiobjective optimization, it was found that there was a "best-tradeoff" AD design, which could balance all the three metrics, without excessively sacrificing each metric. However, for each single-objective optimum AD design, only one metric was optimal, and distinct sacrifices were observed in the other two metrics. For a commercially available biconvex-core AD design, the metrics were even worse than the poorest outcomes of the single-objective optimum AD designs. Therefore, multiobjective design optimization could be useful for achieving native lumbar segment biomechanics and minimal PCPs, as well as for improving the existing lumbar motion-preserving surgical treatments.

Total disc arthroplasties change the kinematics of functional spinal units during lateral bending

BACKGROUND

Information about kinematics in different functional spinal units before and after total disc arthroplasties is necessary to improve prostheses and determine indications. There is little information about the nonstationary instantaneous helical axis of rotation under lateral bending in the cervical spine before and after total disc arthroplasty.

METHODS

Kinematic analyses were performed with an established measuring apparatus on 8 human functional spinal units (C3/C4, C5/C6) under intact conditions and after total disc arthroplasty with two different types of prostheses: Bryan and Prestige. The instantaneous helical axis, migration, and stiffness of the segments were calculated.

FINDINGS

The instantaneous helical axis direction was always inclined ventrally. Ventral inclination was significantly higher in segment C3/C4 than in segment C5/C6 under all conditions (p < 0.001). Both types of arthroplasties significantly increased ventral inclination compared to intact conditions. In both segments, the path length of the instantaneous helical axis' migration was significantly longer after total disc arthroplasty with Bryan (p = 0.001) and shorter after Prestige (p < 0.001) prostheses than under intact conditions. After both types of arthroplasties, the migration path length was significantly longer and the stiffness was significantly lower in segment C3/C4 than in segment C5/C6.

INTERPRETATION

Both types of arthroplasties changed the kinematics of both segments during lateral bending. Altered instantaneous helical axis migration, greater ventral inclination and less stiffness after both arthroplasties indicate unphysiological motion. Both arthroplasties had greater impact on segment C3/C4 than on segment C5/C6 in terms of hypermobility. Increased translational motion after total disc arthroplasty with a Bryan prosthesis might be caused by the prosthetic design.

Biomechanical analysis of single-level interbody fusion with different internal fixation rod materials: a finite element analysis

Background

Lumbar spinal fusion with rigid spinal fixators as one of the high risk factors related to adjacent-segment failure. The purpose of this study is to investigate how the material properties of spinal fixation rods influence the biomechanical behavior at the instrumented and adjacent levels through the use of the finite element method.

Methods

Five finite element models were constructed in our study to simulate the human spine pre- and post-surgery. For the four post-surgical models, the spines were implanted with rods made of three different materials: (i) titanium rod, (ii) PEEK rod with interbody PEEK cage, (iii) Biodegradable rod with interbody PEEK cage, and (iv) PEEK cage without pedicle screw fixation (no rods).

Results

Fusion of the lumbar spine using PEEK or biodegradable rods allowed a similar ROM at both the fusion and adjacent levels under all conditions. The models with PEEK and biodegradable rods also showed a similar increase in contact forces at adjacent facet joints, but both were less than the model with a titanium rod.

Conclusions

Flexible rods or cages with non-instrumented fusion can mitigate the increased contact forces on adjacent facet joints typically found following spinal fixation, and could also reduce the level of stress shielding at the bone graft.

The strain distribution in the lumbar anterior longitudinal ligament is affected by the loading condition and bony features: An in vitro full-field analysis

The role of the ligaments is fundamental in determining the spine biomechanics in physiological and pathological conditions. The anterior longitudinal ligament (ALL) is fundamental in constraining motions especially in the sagittal plane. The ALL also confines the intervertebral discs, preventing herniation. The specific contribution of the ALL has indirectly been investigated in the past as a part of whole spine segments where the structural flexibility was measured. The mechanical properties of isolated ALL have been measured as well. The strain distribution in the ALL has never been measured under pseudo-physiological conditions, as part of multi-vertebra spine segments. This would help elucidate the biomechanical function of the ALL. The aim of this study was to investigate in depth the biomechanical function of the ALL in front of the lumbar vertebrae and of the intervertebral disc. Five lumbar cadaveric spine specimens were subjected to different loading scenarios (flexion-extension, lateral bending, axial torsion) using a state-of-the-art spine tester. The full-field strain distribution on the anterior surface was measured using digital image correlation (DIC) adapted and validated for application to spine segments. The measured strain maps were highly inhomogeneous: the ALL was generally more strained in front of the discs than in front of the vertebrae, with some locally higher strains both imputable to ligament fibers and related to local bony defects. The strain distributions were significantly different among the loading configurations, but also between opposite directions of loading (flexion vs. extension, right vs. left lateral bending, clockwise vs. counterclockwise torsion). This study allowed for the first time to assess the biomechanical behaviour of the anterior longitudinal ligament for the different loading of the spine. We were able to identify both the average trends, and the local effects related to osteophytes, a key feature indicative of spine degeneration.

Total disc arthroplasties alter the characteristics of the instantaneous helical axis of the cervical functional spinal units C3/C4 and C5/C6 during flexion and extension in in vitro conditions

Total disc arthroplasty (TDA) increases the risk of adjacent segment disease (ASD). Kinematic analyses are necessary to compare the intact condition (IC) with alterations after TDA to develop better prostheses. A well-established 6D measuring apparatus (resolution < 2.4 μm; 400 positions/cycle) was used. Kinematics of the flexion and extension of 8 human cervical spine segments (cFSU) C3/C4 and C5/C6 (67.9 ± 13.2 y) were analyzed in the IC and after TDA (Bryan® Cervical Disc [B-TDA], Prestige LP® Cervical Disc [P-TDA]). The migration of the instantaneous helical axis (IHA) and the stiffness of the segments were calculated. Analyses demonstrated a stretched U-curved IHA migration in the sagittal plane. The IHA positions were significantly more cranial in cFSU C5/C6 than in C3/C4 in IC and after either TDA (IC: p < 0.001; B-TDA: p = 0.001; P-TDA: p = 0.045). In cFSU C3/C4 IHA positions shifted anteriocranially after either TDA (p < 0.001). In cFSU C5/C6, the IHA positions were significantly more anterocranial after B-TDA than in IC and after P-TDA (anterior: p < 0.001; cranial: p = 0.005). After B-TDA, the IHA migration path length was significantly longer in cFSU C3/C4 than in C5/C6 (p = 0.007) and longer than in IC in both cFSU (C3/C4: p = 0.047; C5/C6: p < 0.001). Stiffness was increased after both TDA. Various kinematic alterations were observed after both TDA. Increased translation and IHA position shifting after both TDA might indicate abnormal strain and a derogated benefit of TDA. These results imply the most abnormal strain after B-TDA. The lower cFSU might be more susceptible to alterations after TDA than the upper cFSU.

Biomechanical In Vitro Test of a Novel Dynamic Spinal Stabilization System Incorporating Polycarbonate Urethane Material Under Physiological Conditions

Posterior dynamic stabilization systems (PDSS) were developed to provide stabilization to pathologic or hypermobile spinal segments while maintaining the healthy biomechanics of the spine. Numerous novel dynamic devices incorporate the temperature and moisture dependent material polycarbonate urethane (PCU) due to its mechanical properties and biocompatibility. In this study, standardized pure moment in vitro tests were carried out on human lumbar spines to evaluate the performance of a device containing PCU. An environmental chamber with controlled moisture and temperature was included in the setup to meet the requirements of testing under physiological conditions. Three test conditions were compared: (1) native spine, (2) dynamic instrumentation, and (3) dynamic instrumentation with decompression. The ranges of motion, centers of rotation, and relative pedicle screw motions were evaluated. The device displayed significant stiffening in flexion-extension, lateral bending, and axial rotation load directions. A reduction of the native range of motion diminished the stiffening effect along the spinal column and has the potential to reduce the risk of the onset of degeneration of an adjacent segment. In combination with decompression, the implant decreased the native range of motion for flexion-extension and skew bending, but not for lateral bending and axial rotation. Curve fittings using the sigmoid function were performed to parameterize all load-deflection curves in order to enhance accurate numerical model calibrations and comparisons. The device caused a shift of the center of rotation (COR) in the posterior and caudal direction during flexion-extension loading.

Effects on hip stress following sacroiliac joint fixation: A finite element study

For those patients who suffer from low back pain generated by the sacroiliac (SI) joint, fusion of the SI joint has proven to be an effective means of stabilizing it and reducing pain. Though it has shown promise, SI joint fusion raises clinical questions regarding its effect on neighboring joints such as the hip. As such, the purpose of this study was to determine the effects of SI joint fixation on the hip. A finite element spine-sacroiliac-hip (SSIH) model was developed and its functionality was verified against SI joint range of motion (ROM) and hip contact stress, respectively. The intact model was fixed in double leg stance at the distal femora, and loading was applied at the lumbar spine to simulate stance, flexion, extension, right and left lateral bending, and right and left axial rotation. Functionality was confirmed by measuring and comparing SI joint ROM and contact stress and area at the hip with data from the literature. Following verification of the intact SSIH model, both unilateral and bilateral SI joint fixation were modeled; hip contact stress and area were compared to the intact state. Average hip contact stress was ~2 MPa, with most motions resulting in changes less than 5% relative to intact; contact area changed less than 10% for any motion. Clinical significance: these results demonstrated that SI joint fixation with triangular titanium implants imparted little change in stress at the hip, which suggests that the risk of developing adjacent segment disease is likely low. Future clinical studies may be executed to confirm the results of this computational study.

In vivo changes in adjacent segment kinematics after lumbar decompression and fusion

The pathogenesis of lumbar adjacent segment disease is thought to be secondary to altered biomechanics resulting from fusion. Direct in vivo evidence for altered biomechanics following lumbar fusion is lacking. This study's aim was to describe in vivo kinematics of the superior adjacent segment relative to the fused segment before and after lumbar fusion. This study analyzed seven patients with symptomatic lumbar degenerative spondylolisthesis (5 M, 2F; age 65 ± 5.1 years) using a biplane radiographic imaging system. Each subject performed two to three trials of continuous flexion of their torso according to established protocols. Synchronized biplane radiographs were acquired at 20 images per second one month before and six months after single-level fusion at L4-L5 or L5-S1, or two-level fusion at L3-L5 or L4-S1. A previously validated volumetric model-based tracking process was used to track the position and orientation of vertebrae in the radiographic images. Intervertebral flexion/extension and AP translation (slip) at the superior adjacent segment were calculated over the entire dynamic flexion activity. Skin-mounted surface markers were tracked using conventional motion analysis and used to determine torso flexion. Change in adjacent segment kinematics after fusion was determined at corresponding angles of dynamic torso flexion. Changes in adjacent segment motion varied across patients, however, all patients maintained or increased the amount of adjacent segment slip or intervertebral flexion/extension. No patients demonstrated both decreased adjacent segment slip and decreased rotation. This study suggests that short-term changes in kinematics at the superior adjacent segment after lumbar fusion appear to be patient-specific.

Development of a Biconcave Mobile-Bearing Lumbar Total Disc Arthroplasty Concept Using Finite Element Analysis and Design Optimization

Total disc arthroplasty (TDA) is a motion-preserving surgical treatment for spinal disorders. However, adverse surgical outcomes, such as abnormal kinematics, facet joint (FJ) overloading, and polyethylene (PE) failures, have limited wide application of lumbar TDAs. The objectives of this computational study were to elucidate how implant design and FJ articulation both influence metal-on-polyethylene (MoP) motion and contact mechanics, as well as to propose and refine a new mobile-bearing TDA concept which enhanced postoperative performance. Simulation results show that abnormal motions (lift-off and/or unsymmetrical motion) are alleviated in fixed-/mobile-bearing TDA-treated segments, as the FJ gap increases. It clearly demonstrates that FJ articulation guides segmental motion and interferes with intended MoP articulation. For an existing biconvex mobile-bearing design, component impingement leads to a peak PE stress of 20.8 MPa (yield stress: 13 MPa), indicating a high risk of PE creep/fracture. Therefore, we proposed a new TDA concept featuring a biconcave PE core with a smooth shape, in order to strengthen the PE rim and mitigate edge-loading. Furthermore, the biconcave-core TDA was optimally designed to promote normal segmental range of motion (ROM), or to minimize polyethylene contact pressure (PCP). In extension (the severest loading scenario), the biconvex-core TDA design caused a ROM 3.6° (+88%) greater than the intact segment and a peak PCP of 116.5 MPa. In contrast, ROM-optimal or PCP-optimal biconcave-core TDA designs decreased the ROM difference to 0.0° or the peak PCP to 24.3 MPa. Therefore, this new TDA design can potentially reduce the incidence of hypermotion and PE damage. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:1805-1816, 2019.

An upper bound computational model for investigation of fusion effects on adjacent segment biomechanics of the lumbar spine

Prediction of the biomechanical effects of fusion surgery on adjacent segments is a challenge in computational biomechanics of the spine. In this study, a two-segment L3-L4-L5 computational model was developed to simulate the effects of spinal fusion on adjacent segment biomechanical responses under a follower load condition. The interaction between the degenerative segment (L4-5) and the adjacent segment (L3-4) was simulated using an equivalent follower spring. The spring stiffness was calibrated using a rigid fusion of a completely degenerated disc model at the L4-5 level, resulting in an upper bound response at the adjacent (L3-4) segment. The obtained upper bound equivalent follower spring was used to simulate the upper bound biomechanical responses of fusion of the disc with different degeneration grades. It was predicted that as the disc degeneration grade at the degenerative segment decreased, the effect on the adjacent segment responses decreased accordingly after fusion. The data indicated that the upper bound computational model can be a useful computational tool for evaluation of the interaction between segments and for investigation of the biomechanical mechanisms of adjacent segment degeneration after fusion.

Optimization of compressive loading parameters to mimic in vivo cervical spine kinematics in vitro

The human cervical spine supports substantial compressive load in vivo. However, the traditional in vitro testing methods rarely include compressive loads, especially in investigations of multi-segment cervical spine constructs. Previously, a systematic comparison was performed between the standard pure moment with no compressive loading and published compressive loading techniques (follower load - FL, axial load - AL, and combined load - CL). The systematic comparison was structured a priori using a statistical design of experiments and the desirability function approach, which was chosen based on the goal of determining the optimal compressive loading parameters necessary to mimic the segmental contribution patterns exhibited in vivo. The optimized set of compressive loading parameters resulted in in vitro segmental rotations that were within one standard deviation and 10% of average percent error of the in vivo mean throughout the entire motion path. As hypothesized, the values for the optimized independent variables of FL and AL varied dynamically throughout the motion path. FL was not necessary at the extremes of the flexion-extension (FE) motion path but peaked through the neutral position, whereas, a large negative value of AL was necessary in extension and increased linearly to a large positive value in flexion. Although further validation is required, the long-term goal is to develop a "physiologic" in vitro testing method, which will be valuable for evaluating adjacent segment effect following spinal fusion surgery, disc arthroplasty instrumentation testing and design, as well as mechanobiology experiments where correct kinematics and arthrokinematics are critical.

Material failure in dynamic spine implants: are the standardized implant tests before market launch sufficient?

PURPOSE

International Standards Organization (ISO) 12189 and American Society for Testing and Materials F2624 are two standard material specification and test methods for spinal implant devices. The aim of this study was to assess whether the existing and required tests before market launch are sufficient.

METHODS

In three prospective studies, patients were treated due to degenerative disease of the lumbar spine or spondylolisthesis with lumbar interbody fusion and dynamic stabilization of the cranial adjacent level. The CD HORIZON BalanC rod and S⁴ Dynamic rod were implanted in 45 and 11 patients, respectively.

RESULTS

A fatigue fracture of the material of the topping off system has been found in five cases (11%) for the group fitted with the CD HORIZON BalanC rod. In the group using the S⁴ Dynamic rod group, a material failure of the dynamic part was demonstrated in seven patients (64%). All three studies were interrupted due to these results, and a report to the Federal Institute for Drugs and Medical Devices was generated.

CONCLUSION

Spinal implants have to be checked by a notified body before market launch. The notified body verifies whether the implants fulfil the requirements of the current standards. These declared studies suggest that the current standards for the testing of load bearing capacity and stand ability of dynamic spine implants might be insufficient. Revised standards depicting sufficient deformation and load pattern have to be developed and counted as a requirement for the market launch of an implant. These slides can be retrieved under Electronic Supplementary Material.

Bone-Preserving Decompression Procedures Have a Minor Effect on the Flexibility of the Lumbar Spine

Objective

To mitigate the risk of iatrogenic instability, new posterior decompression techniques able to preserve musculoskeletal structures have been introduced but never extensively investigated from a biomechanical point of view. This study was aimed to investigate the impact on spinal flexibility caused by a unilateral laminotomy for bilateral decompression, in comparison to the intact condition and a laminectomy with preservation of a bony bridge at the vertebral arch. Secondary aims were to investigate the biomechanical effects of two-level decompression and the quantification of the restoration of stability after posterior fixation.

Methods

A universal spine tester was used to measure the flexibility of six L2–L5 human spine specimens in intact conditions and after decompression and fixation surgeries. An incremental damage protocol was applied : 1) unilateral laminotomy for bilateral decompression at L3–L4; 2) on three specimens, the unilateral laminotomy was extended to L4–L5; 3) laminectomy with preservation of a bony bridge at the vertebral arch (at L3–L4 in the first three specimens and at L4–L5 in the rest); and 4) pedicle screw fixation at the involved levels.

Results

Unilateral laminotomy for bilateral decompression had a minor influence on the lumbar flexibility. In flexion-extension, the median range of motion increased by 8%. The bone-preserving laminectomy did not cause major changes in spinal flexibility. Two-level decompression approximately induced a twofold destabilization compared to the single-level treatment, with greater effect on the lower level. Posterior fixation reduced the flexibility to values lower than in the intact conditions in all cases.

Conclusion

In vitro testing of human lumbar specimens revealed that unilateral laminotomy for bilateral decompression and bone-preserving laminectomy induced a minor destabilization at the operated level. In absence of other pathological factors (e.g., clinical instability, spondylolisthesis), both techniques appear to be safe from a biomechanical point of view.

The shape and mobility of the thoracic spine in asymptomatic adults - A systematic review of in vivo studies

A comprehensive knowledge of the thoracic shape and kinematics is essential for effective risk prevention, diagnose and proper management of thoracic disorders and assessment of treatment or rehabilitation strategies as well as for in silico and in vitro models for realistic applications of boundary conditions. After an extensive search of the existing literature, this study summarizes 45 studies on in vivo thoracic kyphosis and kinematics and creates a systematic and detailed database. The thoracic kyphosis over T1-12 determined using non-radiological devices (34°) was relatively less than measured using radiological devices (40°) during standing. The majority of kinematical measurements are based on non-radiological devices. The thoracic range of motion (RoM) was greatest during axial rotation (40°), followed by lateral bending (26°), and flexion (21°) when determined using non-radiological devices during standing. The smallest RoM was identified during extension (13°). The lower thoracic level (T8-12) contributed more to the RoM than the upper (T1-4) and middle (T4-8) levels during flexion and lateral bending. During axial rotation and extension, the middle level (T4-8) contributed the most. Coupled motion was evident, mostly during lateral bending and axial rotation. With aging, the thoracic kyphosis increased by about 3° per decade, whereas the RoM decreased by about 5° per decade for all load directions. These changes with aging mainly occurred in the lower region (T6-12). The influence of sex on thoracic kyphosis and the RoM has been described as partly contradictory. Obesity was found to decrease the thoracic RoM. Studies comparing standing, sitting and lying reported the effect of posture as significant.

A new in vitro spine test rig to track multiple vertebral motions under physiological conditions

In vitro pure moment spine tests are commonly used to analyse surgical implants in cadaveric models. Most of the tests are performed at room temperature. However, some new dynamic instrumentation devices and soft tissues show temperature-dependent material properties. Therefore, the aim of this study is to develop a new test rig, which allows applying pure moments on lumbar spine specimens in a vapour-filled chamber at body temperature. As no direct sight is given in the vapour-filled closed chamber, a magnetic tracking (MT) system with implantable receivers was used. Four human cadaveric lumbar spines (L2–L5) were tested in a vapour atmosphere at body temperature with a native and rigid instrumented group. In conclusion, the experimental set-up allows vertebral motion tracking of multiple functional spinal units (FSUs) in a moisture environment at body temperature.

The influence of spinal fusion length on proximal junction biomechanics: a parametric computational study

PURPOSE

Proximal junctional kyphosis and failure are frequent complications in adult spinal deformity surgery with long fusion constructs. The aim of this study was to assess the biomechanics of the proximal segment for fusions of various lengths.

METHODS

A previously established musculoskeletal model of thoracolumbar spine was used to simulate full-range flexion task for fusions (modeled by introduction of rigid constraints) with lower instrumented vertebrae at L5 or S1 and upper instrumented vertebrae (UIV) at any level above, up to T2. Inverse dynamics simulations with force-dependent kinematics were performed for gradually increasing spinal flexion in order to predict global and segmental range of flexion, maximum passive moment, segmental compression and shear forces, which were compared to the uninstrumented case.

RESULTS

For long fusions, with the UIV at T11 or higher, the model predicted an increase in segmental flexion (by 33-860%, or 1.6°-4.7°) and passive moment (by 39-1370%, or 13-31 Nm) at the proximal junction-generally increasing with fusion length. While the maximum shear force was 57-239% (135-283 N) higher for the proximal junction at the upper thorax (UIV at T6 or above), the compression forces were reduced by up to 44% (375 N).

CONCLUSIONS

The length of the instrumentation has an important effect on the proximal segment biomechanics. Despite the limitations of the current model, musculoskeletal modeling appears to be a promising and versatile method to support planning of spinal instrumentation surgeries in the future. These slides can be retrieved under Electronic Supplementary Material.