Biomechanical study of the spinal cord in thoracic ossification of the posterior longitudinal ligament

Effect of Posterior Decompression on the Spinal Cord of Thoracic Ossification of the Posterior Longitudinal Ligament: A Finite Element Analysis

OBJECTIVE

Thoracic ossification of the posterior longitudinal ligament (T-OPLL) causes myelopathy. Although posterior decompression for T-OPLL has shown positive results, patients with kyphotic curvatures often endure poor outcomes. Posterior decompression with fusion (PDF) has demonstrated better results compared to posterior decompression alone. This study aims to evaluate the effects of the posterior procedures for T-OPLL.

METHODS

A 3-dimensional finite element model of the C2-T12 spine, created from medical images, was used to develop the following T3-T4 OPLL compression models: an intact model (no surgery), 25% canal occupancy ratio (COR) OPLL, a discontinuous 25% COR OPLL, a continuous 50% COR OPLL, and a discontinuous 50% COR OPLL. These models were analyzed to evaluate the effects of posterior decompression (laminectomy [LN]) with varied fixation lengths (LN T3-T4, PDF T3-T4, LN T2-T5, and PDF T2-T5) in neutral, flexion, and extension positions.

RESULTS

Increased discontinuity in OPLL led to increased stress on the spinal cord. Posterior decompression reduced spinal cord stress in the neutral posture. However, in flexion and extension, spinal cord stress increased for LN T3-T4, LN T2-T5, and PDF T3-T4 compared to the neutral posture. Notably, PDF T2-T5 prevented an increase in spinal cord stress during these motions.

CONCLUSIONS

Effective management of intervertebral mobility and the appropriate length of decompression are crucial for addressing the thickness and mobility of T-OPLL.

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.

Effects of a contusion load on spinal cord with different curvatures

The cervical spine injury is a complicated procedure in the combination of different injury loads and postures. The aim of this study is to investigate the injury mechanism considering different types of cervical curvatures subjected to contusion loads. A finite element model of a cervical spinal cord was constructed. Gray matter, white matter and pia matter were modeled and hyperelastic material properties were assigned. Convergence analysis and validation analysis were carried out. The model was simulated in 3 different spinal curvatures and loaded by 2 directions with 4 compression degrees. The maximum von Mises stress in the whole model was concentrated in the pia matter in all loading cases. When investigating spinal cord injury, the pia matter must be considered. For all three curvatures, the stress in the gray matter and white matter was higher in front-to-back loading condition than that in back-to-front loading condition. The front-to-back impact may cause a larger damage. A back-to-front load damaged the structure around the central canal and a front-to-back contusion load damaged the anterior horn of the spinal cord at most time. From the view of the maximum stress, the lordotic curvature did not show significant buffering effect. However, the pathological curvature had large areas affected and the lordotic curvature showed some benefits to some degree from the view of stress distribution.

Compression analysis of the gray and white matter of the spinal cord

The spinal cord is composed of gray matter and white matter. It is well known that the properties of these two tissues differ considerably. Spinal diseases often present with symptoms that are caused by spinal cord compression. Understanding the mechanical properties of gray and white matter would allow us to gain a deep understanding of the injuries caused to the spinal cord and provide information on the pathological changes to these distinct tissues in several disorders. Previous studies have reported on the physical properties of gray and white matter, however, these were focused on longitudinal tension tests. Little is known about the differences between gray and white matter in terms of their response to compression. We therefore performed mechanical compression test of the gray and white matter of spinal cords harvested from cows and analyzed the differences between them in response to compression. We conducted compression testing of gray matter and white matter to detect possible differences in the collapse rate. We found that increased compression (especially more than 50% compression) resulted in more severe injuries to both the gray and white matter. The present results on the mechanical differences between gray and white matter in response to compression will be useful when interpreting findings from medical imaging in patients with spinal conditions.

Analysis of stress application at the thoracolumbar junction and influence of vertebral body collapse on the spinal cord and cauda equina

The thoracolumbar junction comprises the spinal cord, nerve roots and the cauda equina, exhibiting unique anatomical features that may give rise to a diverse array of symptoms under conditions of injury, thus complicating the diagnosis of compressive disorders. The present study aimed to examine varying degrees and forms of compression at this level of the spinal cord using a two-dimensional model to calculate the relationship of these variables to injury. The degree of compression was expressed as a percentage of the spinal canal that was occupied. Results were compared with findings from clinical observations to assess the validity of the model. Analysis revealed that higher levels of compression/spinal canal occupation are associated with the presence of neurological symptoms. This finding was consistent with clinical data. Results of the present analysis warrant further research involving evaluation of compression with respect to other parameters, such as blood flow, as well as more anatomically accurate three-dimensional analysis.

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.

Risk Factors for Ineffectiveness of Posterior Decompression and Dekyphotic Corrective Fusion with Instrumentation for Beak-Type Thoracic Ossification of the Posterior Longitudinal Ligament: A Single Institute Study

BACKGROUND: Thoracic ossification of the posterior longitudinal ligament (T-OPLL) is treated surgically with instrumented posterior decompression and fusion. However, the factors determining the outcome of this approach and the efficacy of additional resection of T-OPLL are unknown.

OBJECTIVE: To identify these factors in a prospective study at a single institution.

METHODS: The subjects were 70 consecutive patients with beak-type T-OPLL who underwent posterior decompression and dekyphotic fusion and had an average of 4.8 years of follow-up (minimum of 2 years). Of these patients, 4 (6%; group R) had no improvement or aggravation, were not ambulatory for 3 weeks postoperatively, and required additional T-OPLL resection; while 66 (group N) required no further T-OPLL resection. Clinical records, gait status, intraoperative ultrasonography, intraoperative neurophysiological monitoring (IONM), plain radiography, computed tomography and magnetic resonance imaging findings, and Japanese Orthopaedic Association (JOA) score were compared between the groups.

RESULTS: Preoperatively, patients in group R had significantly higher rates of severe motor paralysis, nonambulatory status, positive prone and supine position test, no spinal cord floating in intraoperative ultrasonography, and deterioration of IONM at the end of surgery (P < .05). In preoperative radiography, the OPLL spinal cord kyphotic angle difference in fused area, OPLL length, and OPLL canal stenosis were significantly higher in group R (P < .05). At final follow-up, JOA scores improved similarly in both groups.

CONCLUSION: Preoperative severe motor paralysis, nonambulatory status, positive prone and supine position test, radiographic spinal cord compression due to beak-type T-OPLL, and intraoperative residual spinal cord compression and deterioration of IONM were associated with ineffectiveness of posterior decompression and fusion with instrumentation. Our 2-stage strategy may be appropriate for beak-type T-OPLL surgery.

Cervical ossification of the posterior longitudinal ligament: factors affecting the effect of posterior decompression

OBJECTIVE

Decompression procedures for cervical myelopathy of ossification of the posterior longitudinal ligament (OPLL) are anterior decompression with fusion, laminoplasty, and posterior decompression with fusion. Preoperative and postoperative stress analyses were performed for compression from hill-shaped cervical OPLL using 3-dimensional finite element method (FEM) spinal cord models.

METHODS

Three FEM models of vertebral arch, OPLL, and spinal cord were used to develop preoperative compression models of the spinal cord to which 10%, 20%, and 30% compression was applied; a posterior compression with fusion model of the posteriorly shifted vertebral arch; an advanced kyphosis model following posterior decompression with the spinal cord stretched in the kyphotic direction; and a combined model of advanced kyphosis following posterior decompression and intervertebral mobility. The combined model had discontinuity in the middle of OPLL, assuming the presence of residual intervertebral mobility at the level of maximum cord compression, and the spinal cord was mobile according to flexion of vertebral bodies by 5°, 10°, and 15°.

RESULTS

In the preoperative compression model, intraspinal stress increased as compression increased. In the posterior decompression with fusion model, intraspinal stress decreased, but partially persisted under 30% compression. In the advanced kyphosis model, intraspinal stress increased again. As anterior compression was higher, the stress increased more. In the advanced kyphosis  + intervertebral mobility model, intraspinal stress increased more than in the only advanced kyphosis model following decompression. Intraspinal stress increased more as intervertebral mobility increased.

CONCLUSION

In high residual compression or instability after posterior decompression, anterior decompression with fusion or posterior decompression with instrumented fusion should be considered.

Biomechanical Behaviors in Three Types of Spinal Cord Injury Mechanisms

Clinically, spinal cord injuries (SCIs) are radiographically evaluated and diagnosed from plain radiographs, computed tomography (CT), and magnetic resonance imaging. However, it is difficult to conclude that radiographic evaluation of SCI can directly explain the fundamental mechanism of spinal cord damage. The von-Mises stress and maximum principal strain are directly associated with neurological damage in the spinal cord from a biomechanical viewpoint. In this study, the von-Mises stress and maximum principal strain in the spinal cord as well as the cord cross-sectional area (CSA) were analyzed under various magnitudes for contusion, dislocation, and distraction SCI mechanisms, using a finite-element (FE) model of the cervical spine with spinal cord including white matter, gray matter, dura mater with nerve roots, and cerebrospinal fluid (CSF). A regression analysis was performed to find correlation between peak von-Mises stress/peak maximum principal strain at the cross section of the highest reduction in CSA and corresponding reduction in CSA of the cord. Dislocation and contusion showed greater peak stress and strain values in the cord than distraction. The substantial increases in von-Mises stress as well as CSA reduction similar to or more than 30% were produced at a 60% contusion and a 60% dislocation, while the maximum principal strain was gradually increased as injury severity elevated. In addition, the CSA reduction had a strong correlation with peak von-Mises stress/peak maximum principal strain for the three injury mechanisms, which might be fundamental information in elucidating the relationship between radiographic and mechanical parameters related to SCI.

Stress analysis of the cervical spinal cord: Impact of the morphology of spinal cord segments on stress

OBJECTIVE

Although there are several classifications for cervical myelopathy, these do not take differences between spinal cord segments into account. Moreover, there has been no report of stress analyses for individual segments to date.

METHODS

By using the finite element method, we constructed 3-dimensional spinal cord models comprised of gray matter, white matter, and pia mater of the second to eighth cervical vertebrae (C2-C8). We placed compression components (disc and yellow ligament) at the front and back of these models, and applied compression to the posterior section covering 10%, 20%, 30%, or 40% of the anteroposterior diameter of each cervical spinal cord segment.

RESULTS

Our results revealed that, under compression applied to an area covering 10%, 20%, or 30% of the anteroposterior diameter of the cervical spinal cord segment, sites of increased stress varied depending on the morphology of each cervical spinal cord segment. Under 40% compression, stress was increased in the gray matter, lateral funiculus, and posterior funiculus of all spinal cord segments, and stress differences between the segments were smaller.

CONCLUSION

These results indicate that, under moderate compression, sites of increased stress vary depending on the morphology of each spinal cord segment or the shape of compression components, and also that the variability of symptoms may depend on the direction of compression. However, under severe compression, the differences among the cervical spinal segments are smaller, which may facilitate diagnosis.

Effect of posterior decompression extent on biomechanical parameters of the spinal cord in cervical ossification of the posterior longitudinal ligament

Ossification of the posterior longitudinal ligament is a common cause of the cervical myelopathy due to compression of the spinal cord. Patients with ossification of the posterior longitudinal ligament usually require the decompression surgery, and there is a need to better understand the optimal surgical extent with which sufficient decompression without excessive posterior shifting can be achieved. However, few quantitative studies have clarified this optimal extent for decompression of cervical ossification of the posterior longitudinal ligament. We used finite element modeling of the cervical spine and spinal cord to investigate the effect of posterior decompression extent for continuous-type cervical ossification of the posterior longitudinal ligament on changes in stress, strain, and posterior shifting that occur with three different surgical methods (laminectomy, laminoplasty, and hemilaminectomy). As posterior decompression extended, stress and strain in the spinal cord decreased and posterior shifting of the cord increased. The location of the decompression extent also influenced shifting. Laminectomy and laminoplasty were very similar in terms of decompression results, and both were superior to hemilaminectomy in all parameters tested. Decompression to the extents of C3-C6 and C3-C7 of laminectomy and laminoplasty could be considered sufficient with respect to decompression itself. Our findings provide fundamental information regarding the treatment of cervical ossification of the posterior longitudinal ligament and can be applied to patient-specific surgical planning.

Mechanical properties of nerve roots and rami radiculares isolated from fresh pig spinal cords

No reports have described experiments designed to determine the strength characteristics of spinal nerve roots and rami radiculares for the purpose of explaining the complexity of symptoms of medullary cone lesions and cauda equina syndrome. In this study, to explain the pathogenesis of cauda equina syndrome, monoaxial tensile tests were performed to determine the strength characteristics of spinal nerve roots and rami radiculares, and analysis was conducted to evaluate the stress-strain relationship and strength characteristics. Using the same tensile test device, the nerve root and ramus radiculares isolated from the spinal cords of pigs were subjected to the tensile test and stress relaxation test at load strain rates of 0.1, 1, 10, and 100 s^-1 under identical settings. The tensile strength of the nerve root was not rate dependent, while the ramus radiculares tensile strength tended to decrease as the strain rate increased. These findings provide important insights into cauda equina symptoms, radiculopathy, and clinical symptoms of the medullary cone.

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

BACKGROUND

There are few studies focusing on the prediction of stress distribution according to the types of ossification of the posterior longitudinal ligament, which can be fundamental information associated with clinical aspects such as the relationship between stress level and neurological symptom severity. In this study, the influence of sagittal and axial types of ossification of the posterior longitudinal ligament on mechanical stress in the cervical spinal cord was investigated.

METHODS

A three-dimensional finite element model of the cervical spine with spinal cord was developed and validated. The von Mises stresses in the cord and the reduction in cross-sectional areas and volume of the cord were investigated for various axial and sagittal types according to the occupying ratio of ossification of the posterior longitudinal ligament in the spinal canal.

FINDINGS

The influence of axial type was less than that of the sagittal type, even though the central type showed higher maximum stresses in the cord, especially for the continuous type. With a 60% occupying ratio of ossification of the posterior longitudinal ligament, the maximum stress was significantly high and the cross-sectional area of the spinal cord was reduced by more than 30% of the intact area regardless of sagittal or axial types. Finally, a higher level of sagittal extension would increase the peak cord tissue stress, which would be related to the neurological dysfunction and tissue damage.

INTERPRETATION

Quantitative investigation of biomechanical characteristics such as mechanical stress may provide fundamental information for pre-operative planning of treatment for ossification of the posterior longitudinal ligament.

Cervical ossification of the posterior longitudinal ligament: Biomechanical analysis of the influence of static and dynamic factors

OBJECTIVE

Cervical myelopathy due to ossification of the posterior longitudinal ligament (OPLL) is induced by static factors, dynamic factors, or a combination of both. We used a three-dimensional finite element method (3D-FEM) to analyze the stress distributions in the cervical spinal cord under static compression, dynamic compression, or a combination of both in the context of OPLL.

METHODS

Experimental conditions were established for the 3D-FEM spinal cord, lamina, and hill-shaped OPLL. To simulate static compression of the spinal cord, anterior compression at 10, 20, and 30% of the anterior-posterior diameter of the spinal cord was applied by the OPLL. To simulate dynamic compression, the OPLL was rotated 5°, 10°, and 15° in the flexion direction. To simulate combined static and dynamic compression under 10 and 20% anterior static compression, the OPLL was rotated 5°, 10°, and 15° in the flexion direction.

RESULTS

The stress distribution in the spinal cord increased following static and dynamic compression by cervical OPLL. However, the stress distribution did not increase throughout the entire spinal cord. For combined static and dynamic compression, the stress distribution increased as the static compression increased, even for a mild range of motion (ROM).

CONCLUSION

Symptoms may appear under static or dynamic compression only. However, under static compression, the stress distribution increases with the ROM of the responsible level and this makes it very likely that symptoms will worsen. We conclude that cervical OPLL myelopathy is induced by static factors, dynamic factors, and a combination of both.

Biomechanical analysis of cervical myelopathy due to ossification of the posterior longitudinal ligament: Effects of posterior decompression and kyphosis following decompression

Cervical ossification of the posterior longitudinal ligament (OPLL) results in myelopathy. Conservative treatment is usually ineffective, thus, surgical treatment is required. One of the reasons for the poor surgical outcome following laminoplasty for cervical OPLL is kyphosis. In the present study, a 3-dimensional finite element method (3D-FEM) was used to analyze the stress distribution in preoperative, posterior decompression and kyphosis models of OPLL. The 3D-FEM spinal cord model established in this study consisted of gray and white matter, as well as pia mater. For the preoperative model, 30% anterior static compression was applied to OPLL. For the posterior decompression model, the lamina was shifted backwards and for the kyphosis model, the spinal cord was studied at 10, 20, 30, 40 and 50° kyphosis. In the preoperative model, high stress distributions were observed in the spinal cord. In the posterior decompression model, stresses were lower than those observed in the preoperative model. In the kyphosis model, an increase in the angle of kyphosis resulted in augmented stress on the spinal cord. Therefore, the results of the present study indicated that posterior decompression was effective, but stress distribution increased with the progression of kyphosis. In cases where kyphosis progresses following surgery, detailed follow-ups are required in case the symptoms worsen.

Biomechanical analysis of cervical spondylotic myelopathy: the influence of dynamic factors and morphometry of the spinal cord

OBJECTIVE

Patients with cervical spondylotic myelopathy (CSM) have the same clinical symptoms that vary according to the degree of spinal cord compression and the cross-sectional cord shape. We used a three-dimensional finite element method (3D-FEM) to analyze the stress distributions of the spinal cord with neck extension under three cross-sectional cord shapes.

METHODS

Experimental condition for the 3D-FEM spinal cord, ligamentum flavum, and anterior compression shape (central, lateral, and diffuse types) was established. To simulate neck extension, the spinal cord was extended by 20° and the ligamentum flavum was shifted distally according to movement of the cephalad lamina.

RESULTS

The stress distribution in the spinal cord increased due to invagination of the ligamentum flavum into the neck extension. The range of stress distribution observed for the diffuse type was wider than for the central and lateral types. In addition, the stress distribution in the spinal cord was increased by the pincer movement of the ligamentum flavum and by the anterior compression of the spinal cord. The range of stress distribution observed for the diffuse type under antero-posterior compression was also wider than for the central and lateral types.

CONCLUSION

This simulation model showed that the clinical symptoms of CSM due to compression of the diffuse type may be stronger than for the central and lateral types. Therefore, careful follow-up is recommended for anterior compression of the spinal cord of diffuse type.