Influence of sagittal and axial types of ossification of posterior longitudinal ligament on mechanical stress in cervical spinal cord: A finite element analysis

Finite Element Analysis for Degenerative Cervical Myelopathy: Scoping Review of the Current Findings and Design Approaches, Including Recommendations on the Choice of Material Properties

BACKGROUND

Degenerative cervical myelopathy (DCM) is a slow-motion spinal cord injury caused via chronic mechanical loading by spinal degenerative changes. A range of different degenerative changes can occur. Finite element analysis (FEA) can predict the distribution of mechanical stress and strain on the spinal cord to help understand the implications of any mechanical loading. One of the critical assumptions for FEA is the behavior of each anatomical element under loading (ie, its material properties).

OBJECTIVE

This scoping review aims to undertake a structured process to select the most appropriate material properties for use in DCM FEA. In doing so, it also provides an overview of existing modeling approaches in spinal cord disease and clinical insights into DCM.

METHODS

We conducted a scoping review using qualitative synthesis. Observational studies that discussed the use of FEA models involving the spinal cord in either health or disease (including DCM) were eligible for inclusion in the review. We followed the PRISMA-ScR (Preferred Reporting Items for Systematic Reviews and Meta-Analyses extension for Scoping Reviews) guidelines. The MEDLINE and Embase databases were searched to September 1, 2021. This was supplemented with citation searching to retrieve the literature used to define material properties. Duplicate title and abstract screening and data extraction were performed. The quality of evidence was appraised using the quality assessment tool we developed, adapted from the Newcastle-Ottawa Scale, and shortlisted with respect to DCM material properties, with a final recommendation provided. A qualitative synthesis of the literature is presented according to the Synthesis Without Meta-Analysis reporting guidelines.

RESULTS

A total of 60 papers were included: 41 (68%) "FEA articles" and 19 (32%) "source articles." Most FEA articles (33/41, 80%) modeled the gray matter and white matter separately, with models typically based on tabulated data or, less frequently, a hyperelastic Ogden variant or linear elastic function. Of the 19 source articles, 14 (74%) were identified as describing the material properties of the spinal cord, of which 3 (21%) were considered most relevant to DCM. Of the 41 FEA articles, 15 (37%) focused on DCM, of which 9 (60%) focused on ossification of the posterior longitudinal ligament. Our aggregated results of DCM FEA indicate that spinal cord loading is influenced by the pattern of degenerative changes, with decompression alone (eg, laminectomy) sufficient to address this as opposed to decompression combined with other procedures (eg, laminectomy and fusion).

CONCLUSIONS

FEA is a promising technique for exploring the pathobiology of DCM and informing clinical care. This review describes a structured approach to help future investigators deploy FEA for DCM. However, there are limitations to these recommendations and wider uncertainties. It is likely that these will need to be overcome to support the clinical translation of FEA to DCM.

Effect of Different Types of Ossification of the Posterior Longitudinal Ligament on the Dynamic Biomechanical Response of the Spinal Cord: A Finite Element Analysis

Ossification of the posterior longitudinal ligament (OPLL) has been identified as an important cause of cervical myelopathy. However, the biomechanical mechanism between the OPLL type and the clinical characteristics of myelopathy remains unclear. The aim of this study was to evaluate the effect of different types of OPLL on the dynamic biomechanical response of the spinal cord. A three-dimensional finite element model of the fluid-structure interaction of the cervical spine with spinal cord was established and validated. The spinal cord stress and strain, cervical range of motion (ROM) in different types of OPLL models were predicted during dynamic flexion and extension activity. Different types of OPLL models showed varying degrees of increase in stress and strain under the process of flexion and extension, and there was a surge toward the end of extension. Larger spinal cord stress was observed in segmental OPLL. For continuous and mixed types of OPLL, the adjacent segments of OPLL showed a dramatic increase in ROM, while the ROM of affected segments was limited. As a dynamic factor, flexion and extension of the cervical spine play an amplifying role in OPLL-related myelopathy, while appropriate spine motion is safe and permitted. Segmental OPLL patients are more concerned about the spinal cord injury induced by large stress, and patients with continuous OPLL should be noted to progressive injuries of adjacent level.

Finite element modeling of the human cervical spinal cord and its applications: A systematic review

Background Context

Finite element modeling (FEM) is an established tool to analyze the biomechanics of complex systems. Advances in computational techniques have led to the increasing use of spinal cord FEMs to study cervical spinal cord pathology. There is considerable variability in the creation of cervical spinal cord FEMs and to date there has been no systematic review of the technique. The aim of this study was to review the uses, techniques, limitations, and applications of FEMs of the human cervical spinal cord.

Methods

A literature search was performed through PubMed and Scopus using the words finite element analysis, spinal cord, and biomechanics. Studies were selected based on the following inclusion criteria: (1) use of human spinal cord modeling at the cervical level; (2) model the cervical spinal cord with or without the osteoligamentous spine; and (3) the study should describe an application of the spinal cord FEM.

Results

Our search resulted in 369 total publications, 49 underwent reviews of the abstract and full text, and 23 were included in the study. Spinal cord FEMs are used to study spinal cord injury and trauma, pathologic processes, and spine surgery. Considerable variation exists in the derivation of spinal cord geometries, mathematical models, and material properties. Less than 50% of the FEMs incorporate the dura mater, cerebrospinal fluid, nerve roots, and denticulate ligaments. Von Mises stress, and strain of the spinal cord are the most common outputs studied. FEM offers the opportunity for dynamic simulation, but this has been used in only four studies.

Conclusions

Spinal cord FEM provides unique insight into the stress and strain of the cervical spinal cord in various pathological conditions and allows for the simulation of surgical procedures. Standardization of modeling parameters, anatomical structures and inclusion of patient-specific data are necessary to improve the clinical translation.

Effect of Cervical Stenosis and Rate of Impact on Risk of Spinal Cord Injury During Whiplash Injury

STUDY DESIGN

Finite Element Study.

OBJECTIVE

To determine the risk of spinal cord injury with pre-existing cervical stenosis during a whiplash injury.

SUMMARY OF BACKGROUND DATA

Patients with cervical spinal stenosis are often cautioned on the potential increased risk of spinal cord injury (SCI) from minor trauma such as rear impact whiplash injuries. However, there is no consensus on the degree of canal stenosis or the rate of impact that predisposes cervical SCI from minor trauma.

METHODS

A previously validated three-dimensional finite element model of the human head-neck complex with the spinal cord and activated cervical musculature was used. Rear impact acceleration was applied at 1.8 m/s and 2.6 m/s. Progressive spinal stenosis was simulated at the C5 to C6 segment, from 14 mm to 6 mm, at 2 mm intervals of ventral disk protrusion. Spinal cord von Mises stress and maximum principal strain were extracted and normalized with respect to the 14 mm spine at each cervical spine level from C2 to C7.

RESULTS

The mean segmental range of motion was 7.3 degrees at 1.8 m/s and 9.3 degrees at 2.6 m/s. Spinal cord stress above the threshold for SCI was noted at C5 to C6 for 6 mm stenosis at 1.8 m/s and 2.6 m/s. The segment (C6-C7) inferior to the level of maximum stenosis also showed increasing stress and strain with a higher rate of impact. For 8 mm stenosis, spinal cord stress exceeded SCI thresholds only at 2.6 m/s. Spinal cord strain above SCI thresholds were only noted in the 6 mm stenosis model at 2.6 m/s.

CONCLUSION

Increased spinal stenosis and rate of impact are associated with greater magnitude and spatial distribution of spinal cord stress and strain during a whiplash injury. Spinal canal stenosis of 6 mm was associated with consistent elevation of spinal cord stress and strain above SCI thresholds at 2.6 m/s.

DLX5 regulates the osteogenic differentiation of spinal ligaments cells derived from ossification of the posterior longitudinal ligament patients via NOTCH signaling

Background

Ossification of the posterior longitudinal ligaments (OPLL) is common disorder characterized by heterotopic ossification of the spinal ligaments. Mechanical stimulation (MS) plays an important role in OPLL. DLX5 is an essential transcription factor required for osteoblast differentiation. However, the role of DLX5 during in OPLL is unclear. This study aims to investigate whether DLX5 is associated with OPLL progression under MS.

Methods

Stretch stimulation was applied to spinal ligaments cells derived from OPLL (OPLL cells) and non-OPLL (non-OPLL cells) patients. Expression of DLX5 and osteogenesis-related genes were determined by quantitative real-time polymerase chain reaction and Western blot. The osteogenic differentiation ability of the cells was measured using alkaline phosphatase (ALP) staining and alizarin red staining. The protein expression of DLX5 in the tissues and the nuclear translocation of NOTCH intracellular domain (NICD) was examined by immunofluorescence.

Results

Compared with non-OPLL cells, OPLL cells expressed higher levels of DLX5 in vitro and vivo (p < 0.01). Upregulated expression of DLX5 and osteogenesis-related genes (OSX, RUNX2, and OCN) were observed in OPLL cells induced with stretch stimulation and osteogenic medium, whereas there was no change in the non-OPLL cells (p < 0.01). Cytoplasmic NICD protein translocated from the cytoplasm to the nucleus inducing DLX5 under stretch stimulation, which was reduced by the NOTCH signaling inhibitors (DAPT) (p < 0.01).

Conclusions

These data suggest that DLX5 play a critical role in MS-induced progression of OPLL through NOTCH signaling, which provides a new insight into the pathogenesis of OPLL.

Effects of cervical rotatory manipulation on the cervical spinal cord complex with ossification of the posterior longitudinal ligament in the vertebral canal: A finite element study

Background: There are few studies focusing on biomechanism of spinal cord injury according to the ossification of the posterior longitudinal ligament (OPLL) during cervical rotatory manipulation (CRM). This study aimed to explore the biomechanical effects of CRM on the spinal cord, dura matter and nerve roots with OPLL in the cervical vertebral canal. Methods: Three validated FE models of the craniocervical spine and spinal cord complex were constructed by adding mild, moderate, and severe OPLL to the healthy FE model, respectively. We simulated the static compression of the spinal cord by OPLL and the dynamic compression during CRM in the flexion position. The stress distribution of the spinal cord complex was investigated. Results: The cervical spinal cord experienced higher von Mises stress under static compression by the severe OPLL. A higher von Mises stress was observed on the spinal cord in the moderate and severe OPLL models during CRM. The dura matter and nerve roots had a higher von Mises stress in all three models during CRM. Conclusion: The results show a high risk in performing CRM in the flexion position on patients with OPLL, in that different occupying ratios in the vertebral canal due to OPLL could significantly increase the stress on the spinal cord complex.

Stress Distribution on Spinal Cord According to Type of Laminectomy for Large Focal Cervical Ossification of Posterior Longitudinal Ligament Based on Finite Element Method

Most studies on the ossification of the posterior longitudinal ligament (OPLL) using the finite element method were conducted in the neutral state, and the resulting decompression was judged to be good. As these studies do not reflect the actual behavior of the cervical spine, this study conducted an analysis in the neutral state and a biomechanical analysis during flexion and extension behaviors. After validation via the construction of an intact cervical spine model, the focal OPLL model was inserted into the C4-C5 segment and a simulation was performed. The neutral state was shown by applying a fixed condition to the lower part of the T1 and Y-axis fixed condition of the spinal cord and simulating spinal cord compression with OPLL. For flexion and extension simulation, a ±30-degree displacement was additionally applied to the top of the C2 dens. Accordingly, it was confirmed that spinal cord decompression did not work well during the flexion and extension behaviors, but rather increased. Thus, if patients with focal OPLL inevitably need to undergo posterior decompression, additional surgery using an anterior approach should be considered.

Biomechanical comparison of spinal cord compression types occurring in Degenerative Cervical Myelopathy

BACKGROUND

Degenerative Cervical Myelopathy results from spine degenerations narrowing the spinal canal and inducing cord compressions. Prognosis is challenging. This study aimed at simulating typical spinal cord compressions observed in patients with a realistic model to better understand pathogenesis for later prediction of patients' evolution.

METHODS

A 30% reduction in cord cross-sectional area at C5-C6 was defined as myelopathy threshold based on Degenerative Cervical Myelopathy features from literature and MRI measurements in 20 patients. Four main compression types were extracted from MRIs and simulated with a comprehensive three-dimensional finite element spine model. Median diffuse, median focal and lateral types were modelled as disk herniation while circumferential type additionally involved ligamentum flavum hypertrophy. All stresses were quantified along inferior-superior axis, compression development and across atlas-defined spinal cord regions.

FINDINGS

Anterior gray and white matter globally received the highest stress while lateral pathways were the least affected. Median diffuse compression induced the highest stresses. Circumferential type focused stresses in posterior gray matter. Along inferior-superior axis, those two types showed a peak of constraints at compression site while median focal and lateral types showed lower values but extending further.

INTERPRETATION

Median diffuse type would be the most detrimental based on stress amplitude. Anterior regions would be the most at risk, except for circumferential type where posterior regions would be equally affected. In addition to applying constraints, ischemia could be a significant component explaining the early demyelination reported in lateral pathways. Moving towards patient-specific simulations, biomechanical models could become strong predictors for degenerative changes.

Biomechanical Effects on Cervical Spinal Cord and Nerve Root Following Laminoplasty for Ossification of the Posterior Longitudinal Ligament in the Cervical Spine: A Comparison Between Open-Door and Double-Door Laminoplasty Using Finite Element Analysis

Many clinical case series have reported the predisposing factors for C5 palsy and have presented comparisons of the two types of laminoplasty. However, there have been no biomechanical studies focusing on cervical spinal cord and nerve root following laminoplasty. The purpose of this study is to investigate biomechanical changes in the spinal cord and nerve roots following the two most common types of laminoplasty, open-door and double-door laminoplasty, for cervical ossification of the posterior longitudinal ligament (OPLL). A finite element (FE) model of the cervical spine and spinal cord with nerve root complex structures was developed. Stress changes in the spinal cord and nerve roots, posterior shift of the spinal cord, and displacement of the cervical nerve roots were analyzed with two types of cervical laminoplasty models for variations in the degree of canal occupying ratio and shape of the OPLL. The shape and degree of spinal cord compression caused by the OPLL had more influence on the changes in stress, posterior shift of the spinal cord, and displacement of the nerve root than the type of laminoplasty. The lateral-type OPLL resulted in imbalanced stress on the nerve roots and the highest nerve root displacement. Type of laminoplasty and shape and degree of spinal cord compression caused by OPLL were found to influence the changes in stress and posterior displacement of the cervical spinal cord and nerve roots. Lateral-type OPLL might contribute to the development of C5 palsy due to the imbalanced stress and tension on the nerve roots after laminoplasty.

Increased stress and strain on the spinal cord due to ossification of the posterior longitudinal ligament in the cervical spine under flexion after laminectomy

Myelopathy in the cervical spine due to cervical ossification of the posterior longitudinal ligament could be induced by static compression and/or dynamic factors. It has been suggested that dynamic factors need to be considered when planning and performing the decompression surgery on patients with the ossification of the posterior longitudinal ligament. A finite element model of the C2-C7 cervical spine in the neutral position was developed and used to generate flexion and extension of the cervical spine. The segmental ossification of the posterior longitudinal ligament on the C5 was assumed, and laminectomy was performed on C4-C6 according to a conventional surgical technique. For various occupying ratios of the ossified ligament between 20% and 60%, von-Mises stresses, maximum principal strains in the spinal cord, and cross-sectional area of the cord were investigated in the pre-operative and laminectomy models under flexion, neutral position, and extension. The results were consistent with previous experimental and computational studies in terms of stress, strain, and cross-sectional area. Flexion leads to higher stresses and strains in the cord than the neutral position and extension, even after decompression surgery. These higher stresses and strains might be generated by residual compression occurring at the segment with the ossification of the posterior longitudinal ligament. This study provides fundamental information under different neck positions regarding biomechanical characteristics of the spinal cord in cervical ossification of the posterior longitudinal ligament.