The effect of bone fragment size and cerebrospinal fluid on spinal cord deformation during trauma: an ex vivo study

Finite element modeling and analysis of effect of preexisting cervical degenerative disease on the spinal cord during flexion and extension

Recent studies have emphasized the importance of dynamic activity in the development of myelopathy. However, current knowledge of how degenerative factors affect the spinal cord during motion is still limited. This study aimed to investigate the effect of various types of preexisting herniated cervical disc and the ligamentum flavum ossification on the spinal cord during cervical flexion and extension. A detailed dynamic fluid-structure interaction finite element model of the cervical spine with the spinal cord was developed and validated. The changes of von Mises stress and maximum principal strain within the spinal cord in the period of normal, hyperflexion, and hyperextension were investigated, considering various types and grades of disc herniation and ossification of the ligamentum flavum. The flexion and extension of the cervical spine with spinal canal encroachment induced high stress and strain inside the spinal cord, and this effect was also amplified by increased canal encroachments and cervical hypermobility. The spinal cord might evade lateral encroachment, leading to a reduction in the maximum stress and principal strain within the spinal cord in local-type herniation. Although the impact was limited in the case of diffuse type, the maximum stress tended to appear in the white matter near the encroachment site while compression from both ventral and dorsal was essential to make maximum stress appear in the grey matter. The existence of canal encroachment can reduce the safe range for spinal cord activities, and hypermobility activities may induce spinal cord injury. Besides, the ligamentum flavum plays an important role in the development of central canal syndrome.Significance. This model will enable researchers to have a better understanding of the influence of cervical degenerative diseases on the spinal cord during extension and flexion.

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.

The biomechanical effect of different types of ossification of the ligamentum flavum on the spinal cord during cervical dynamic activities

Ossification of the ligamentum flavum (OLF) is thought to be an influential etiology of myelopathy, as thickened ligamentum flavum causes the stenosis of the vertebral canal, which could subsequently compress the spinal cord. Unfortunately, there was little information available on the effects of cervical OLF on spinal cord compression, such as the relationship between the progression of cervical OLF and nervous system symptoms during dynamic cervical spine activities. In this research, a finite element model of C1-C7 including the spinal cord featured by dynamic fluid-structure interaction was reconstructed and utilized to analyze how different types of cervical OLF affect principal strain and stress distribution in spinal cord during spinal activities towards six directions. For patients with cervical OLF, cervical extension induces higher stress within the spinal cord among all directions. From the perspective of biomechanics, extension leads to stress concentration in the lateral corticospinal tracts or the posterior of gray matter. Low energy damage to the spinal cord would be caused by the high and fluctuating stresses during cervical movements to the affected side for patients with unilateral OLF at lower grades.

The biomechanical implications of neck position in cervical contusion animal models of SCI

Large animal contusion models of spinal cord injury are an essential precursor to developing and evaluating treatment options for human spinal cord injury. Reducing variability in these experiments has been a recent focus as it increases the sensitivity with which treatment effects can be detected while simultaneously decreasing the number of animals required in a study. Here, we conducted a detailed review to explore if head and neck positioning in a cervical contusion model of spinal cord injury could be a factor impacting the biomechanics of a spinal cord injury, and thus, the resulting outcomes. By reviewing existing literature, we found evidence that animal head/neck positioning affects the exposed level of the spinal cord, morphology of the spinal cord, tissue mechanics and as a result the biomechanics of a cervical spinal cord injury. We posited that neck position could be a hidden factor contributing to variability. Our results indicate that neck positioning is an important factor in studying biomechanics, and that reporting these values can improve inter-study consistency and comparability and that further work needs to be done to standardize positioning for cervical spinal cord contusion injury models.

Effect of degenerative factors on cervical spinal cord during flexion and extension: a dynamic finite element analysis

Spinal cord injury (SCI) is a global problem that brings a heavy burden to both patients and society. Recent investigations indicated degenerative disease is taking an increasing part in SCI with the growth of the aging population. However, little insight has been gained about the effect of cervical degenerative disease on the spinal cord during dynamic activities. In this work, a dynamic fluid-structure interaction model was developed and validated to investigate the effect of anterior and posterior encroachment caused by degenerative disease on the spinal cord during normal extension and flexion. Maximum von-Mises stress and maximum principal strain were observed at the end of extension and flexion. The abnormal stress distribution caused by degenerative factors was concentrated in the descending tracts of the spinal cord. Our finding indicates that the excessive motion of the cervical spine could potentially exacerbate spinal cord injury and enlarge injury areas. Stress and strain remained low compared to extension during moderate flexion. This suggests that patients with cervical degenerative disease should avoid frequent or excessive flexion and extension which could result in motor function impairment, whereas moderate flexion is safe. Besides, encroachment caused by degenerative factors that are not significant in static imaging could also cause cord compression during normal activities.

Comparison of numerical methods for cerebrospinal fluid representation and fluid-structure interaction during transverse impact of a finite element spinal cord model

Spinal cord impacts can have devastating consequences. Computational models can investigate such impacts but require biofidelic numerical representations of the neural tissues and fluid-structure interaction with cerebrospinal fluid. Achieving this biofidelity is challenging, particularly for efficient implementation of the cerebrospinal fluid in full computational human body models. The goal of this study was to assess the biofidelity and computational efficiency of fluid-structure interaction methods representing the cerebrospinal fluid interacting with the spinal cord, dura, and pia mater using experimental pellet impact test data from bovine spinal cords. Building on an existing finite element model of the spinal cord and pia mater, an orthotropic hyperelastic constitutive model was proposed for the dura mater and fit to literature data. The dura mater and cerebrospinal fluid were integrated with the existing finite element model to assess four fluid-structure interaction methods under transverse impact: Lagrange, pressurized volume, smoothed particle hydrodynamics, and arbitrary Lagrangian-Eulerian. The Lagrange method resulted in an overly stiff mechanical response, whereas the pressurized volume method over-predicted compression of the neural tissues. Both the smoothed particle hydrodynamics and arbitrary Lagrangian-Eulerian methods were able to effectively model the impact response of the pellet on the dura mater, outflow of the cerebrospinal fluid, and compression of the spinal cord; however, the arbitrary Lagrangian-Eulerian compute time was approximately five times higher than smoothed particle hydrodynamics. Crucial to implementation in human body models, the smoothed particle hydrodynamics method provided a computationally efficient and representative approach to model spinal cord fluid-structure interaction during transverse impact.

A Hyper-Viscoelastic Continuum-Level Finite Element Model of the Spinal Cord Assessed for Transverse Indentation and Impact Loading

Finite Element (FE) modelling of spinal cord response to impact can provide unique insights into the neural tissue response and injury risk potential. Yet, contemporary human body models (HBMs) used to examine injury risk and prevention across a wide range of impact scenarios often lack detailed integration of the spinal cord and surrounding tissues. The integration of a spinal cord in contemporary HBMs has been limited by the need for a continuum-level model owing to the relatively large element size required to be compatible with HBM, and the requirement for model development based on published material properties and validation using relevant non-linear material data. The goals of this study were to develop and assess non-linear material model parameters for the spinal cord parenchyma and pia mater, and incorporate these models into a continuum-level model of the spinal cord with a mesh size conducive to integration in HBM. First, hyper-viscoelastic material properties based on tissue-level mechanical test data for the spinal cord and hyperelastic material properties for the pia mater were determined. Secondly, the constitutive models were integrated in a spinal cord segment FE model validated against independent experimental data representing transverse compression of the spinal cord-pia mater complex (SCP) under quasi-static indentation and dynamic impact loading. The constitutive model parameters were fit to a quasi-linear viscoelastic model with an Ogden hyperelastic function, and then verified using single element test cases corresponding to the experimental strain rates for the spinal cord (0.32–77.22 s⁻¹) and pia mater (0.05 s⁻¹). Validation of the spinal cord model was then performed by re-creating, in an explicit FE code, two independent ex-vivo experimental setups: 1) transverse indentation of a porcine spinal cord-pia mater complex and 2) dynamic transverse impact of a bovine SCP. The indentation model accurately matched the experimental results up to 60% compression of the SCP, while the impact model predicted the loading phase and the maximum deformation (within 7%) of the SCP experimental data. This study quantified the important biomechanical contribution of the pia mater tissue during spinal cord deformation. The validated material models established in this study can be implemented in computational HBM.

A novel CX3CR1 inhibitor AZD8797 facilitates early recovery of rat acute spinal cord injury by inhibiting inflammation and apoptosis

The present study aimed to evaluate the effect of the CX3CR1 inhibitor AZD8797 in early recovery after acute SCI and elucidate its potential mechanism in blocking inflammation and apoptosis. Adult rats were sacrificed after 3, 7, 10, or 14 days of SCI. The injured spinal tissues were collected for assessing C-X3-C motif chemokine ligand 1(CX3CL1)/C-X3-C motif chemokine receptor 1 (CX3CR1) expression at each time point via western blotting (WB) and quantitative PCR. The cellular localization of the proteins was detected by immunofluorescence. Another batch of rats (subdivided into sham, injury model, AZD8797 and methylprednisolone groups) were used to evaluate locomotive recovery with a Basso Beattie Bresnahan score. Based on the expression level of CX3CR1, these rats were sacrificed at the most prominent stage of CX3CR1 expression (10 days after SCI), for assessing the serum levels of tumor necrosis factor-α/interleukin (IL)-6/IL-1β and the expression of CX3CL1/CX3CR1/caspase 3/Bcl-2/Bax in the spinal cord tissues through WB and ELISA. Additionally, apoptosis and necrosis in the injured spinal cord were evaluated by terminal deoxynucleotidyl transferase-mediated dUTP nick-end labeling staining/fluoro-jade B staining. Expression levels of both CX3CR1 and CX3CL1 reached their peak 10 days after the injury, followed by a dramatic downward trend at 14 days. The enhanced expression of CX3CR1 was detected in astrocytes and microglia of the injured spinal cord. AZD8797 improved locomotive recovery after 10 days of SCI and was as effective as methylprednisolone. The effect of AZD8797 was mediated by suppressing apoptosis, necrosis and inflammatory responses, as assessed by WB/ELISA and morphological examinations. The current study has demonstrated that AZD8797 can effectively block overwhelming inflammation, apoptosis and necrosis after SCI and facilitate early recovery of locomotive function.

Engineering approaches to understanding mechanisms of spinal column injury leading to spinal cord injury

BACKGROUND

The mechanical interactions occurring between the spinal column and spinal cord during an injury event are complex and variable, and likely have implications for the clinical presentation and prognosis of the individual.

METHODS

The engineering approaches that have been developed to better understand spinal column and cord interactions during an injury event are discussed. These include injury models utilising human and animal cadaveric specimens, in vivo anaesthetised animals, finite element models, inanimate physical systems and combinations thereof.

FINDINGS

The paper describes the development of these modelling approaches, discusses the advantages and disadvantages of the various models, and the major outcomes that have had implications for spinal cord injury research and clinical practice.

INTERPRETATION

The contribution of these four engineering approaches to understanding the interaction between the biomechanics and biology of spinal cord injury is substantial; they have improved our understanding of the factors contributing to the spinal column disruption, the degree of spinal cord deformation or motion, and the resultant neurological deficit and imaging features. Models of the injury event are challenging to produce, but technological advances are likely to improve these models and, consequently, our understanding of the mechanical context in which the biological injury occurs.

Strain rate dependent behavior of the porcine spinal cord under transverse dynamic compression

The accurate description of the mechanical properties of spinal cord tissue benefits to clinical evaluation of spinal cord injuries and is a required input for analysis tools such as finite element models. Unfortunately, available data in the literature generally relate mechanical properties of the spinal cord under quasi-static loading conditions, which is not adapted to the study of traumatic behavior, as neurological tissue adopts a viscoelastic behavior. Thus, the objective of this study is to describe mechanical properties of the spinal cord up to mechanical damage, under dynamic loading conditions. A total of 192 porcine cervical to lumbar spinal cord samples were compressed in a transverse direction. Loading conditions included ramp tests at 0.5, 5 or 50 s−1 and cyclic loading at 1, 10 or 20 Hz. Results showed that spinal cord behavior was significantly influenced by strain rate. Mechanical damage occurred at 0.64, 0.68 and 0.73 strains for 0.5, 5 or 50 s−1 loadings, respectively. Variations of behavior between the tested strain rates were explained by cyclic loading results, which revealed behavior more or less viscous depending on strain rate. Also, a parameter (stress multiplication factor) was introduced to allow transcription of a stress–strain behavior curve to different strain rates. This factor was described and was significantly different for cervical, thoracic and lumbar vertebral heights, and for the strain rates evaluated in this study.

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.

Analysis of the Role of CX3CL1 (Fractalkine) and Its Receptor CX3CR1 in Traumatic Brain and Spinal Cord Injury: Insight into Recent Advances in Actions of Neurochemokine Agents

CX3CL1 (fractalkine) is the only member of the CX3C (delta) subfamily of chemokines which is unique and combines the properties of both chemoattractant and adhesion molecules. The two-form ligand can exist either in a soluble form, like all other chemokines, and as a membrane-anchored molecule. CX3CL1 discloses its biological properties through interaction with one dedicated CX3CR1 receptor which belongs to a family of G protein-coupled receptors (GPCR). The CX3CL1/CX3CR1 axis acts in many physiological phenomena including those occurring in the central nervous system (CNS), by regulating the interactions between neurons, microglia, and immune cells. Apart from the role under physiological conditions, the CX3CL1/CX3CR1 axis was implied to have a role in different neuropathologies such as traumatic brain injury (TBI) and spinal cord injury (SCI). CNS injuries represent a serious public health problem, despite improvements in therapeutic management. To date, no effective treatment has been determined, so they constitute a leading cause of death and severe disability. The course of TBI and SCI has two consecutive poorly demarcated phases: the initial, primary injury and secondary injury. Recent evidence has implicated the role of the CX3CL1/CX3CR1 axis in neuroinflammatory processes occurring after CNS injuries. The importance of the CX3CL1/CX3CR1 axis in the pathophysiology of TBI and SCI in the context of systemic and direct local immune response is still under investigation. This paper, based on a review of the literature, updates and summarizes the current knowledge about CX3CL1/CX3CR1 axis involvement in TBI and SCI pathogenesis, indicating possible molecular and cellular mechanisms with a potential target for therapeutic intervention.

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.

EXPERIMENTAL METHODS FOR THE BIOMECHANICAL INVESTIGATION OF THE HUMAN SPINE: A REVIEW

In vitro mechanical testing of spinal specimens is extremely important to better understand the biomechanics of the healthy and diseased spine, fracture, and to test/optimize surgical treatment. While spinal testing has extensively been carried out in the past four decades, testing methods are quite diverse. This paper aims to provide a critical overview of the in vitro methods for mechanical testing the human spine at different scales. Specimens of different type are used, according to the aim of the study: spine segments (two or more adjacent vertebrae) are used both to investigate the spine kinematics, and the mechanical properties of the spine components (vertebrae, ligaments, discs); single vertebrae (whole vertebra, isolated vertebral body, or vertebral body without endplates) are used to investigate the structural properties of the vertebra itself; core specimens are extracted to test the mechanical properties of the trabecular bone at the tissue-level; mechanical properties of spine soft tissue (discs, ligaments, spinal cord) are measured on isolated elements, or on tissue specimens. Identification of consistent reference frames is still a debated issue. Testing conditions feature different pre-conditioning and loading rates, depending on the simulated action. Tissue specimen preservation is a very critical issue, affecting test results. Animal models are often used as a surrogate. However, because of different structure and anatomy, extreme caution is required when extrapolating to the human spine. In vitro loading conditions should be based on reliable in vivo data. Because of the high complexity of the spine, such information (either through instrumented implants or through numerical modeling) is currently unsatisfactory. Because of the increasing ability of computational models in predicting biomechanical properties of musculoskeletal structures, a synergy is possible (and desirable) between in vitro experiments and numerical modeling. Future perspectives in spine testing include integration of mechanical and structural properties at different dimensional scales (from the whole-body-level down to the tissue-level) so that organ-level models (which are used to predict the most relevant phenomena such as fracture) include information from all dimensional scales.

Cerebrospinal Fluid Pressures Resulting From Experimental Traumatic Spinal Cord Injuries in a Pig Model

Despite considerable effort over the last four decades, research has failed to translate into consistently effective treatment options for spinal cord injury (SCI). This is partly attributed to differences between the injury response of humans and rodent models. Some of this difference could be because the cerebrospinal fluid (CSF) layer of the human spine is relatively large, while that of the rodents is extremely thin. We sought to characterize the fluid impulse induced in the CSF by experimental SCIs of moderate and high human-like severity, and to compare this with previous studies in which fluid impulse has been associated with neural tissue injury. We used a new in vivo pig model (n = 6 per injury group, mean age 124.5 days, 20.9 kg) incorporating four miniature pressure transducers that were implanted in pairs in the subarachnoid space, cranial, and caudal to the injury at 30 mm and 100 mm. Tissue sparing was assessed with Eriochrome Cyanine and Neutral Red staining. The median peak pressures near the injury were 522.5 and 868.8 mmHg (range 96.7–1430.0) and far from the injury were 7.6 and 36.3 mmHg (range 3.8–83.7), for the moderate and high injury severities, respectively. Pressure impulse (mmHg.ms), apparent wave speed, and apparent attenuation factor were also evaluated. The data indicates that the fluid pressure wave may be sufficient to affect the severity and extent of primary tissue damage close to the injury site. However, the CSF pressure was close to normal physiologic values at 100 mm from the injury. The high injury severity animals had less tissue sparing than the moderate injury severity animals; this difference was statistically significant only within 1.6 mm of the epicenter. These results indicate that future research seeking to elucidate the mechanical origins of primary tissue damage in SCI should consider the effects of CSF. This pig model provides advantages for basic and preclinical SCI research due to its similarities to human scale, including the existence of a human-like CSF fluid layer.

Development of a large-animal model to measure dynamic cerebrospinal fluid pressure during spinal cord injury

Object

Spinal cord injury (SCI) often results in considerable permanent neurological impairment, and unfortunately, the successful translation of effective treatments from laboratory models to human patients is lacking. This may be partially attributed to differences in anatomy, physiology, and scale between humans and rodent models. One potentially important difference between the rodent and human spinal cord is the presence of a significant CSF volume within the intrathecal space around the human cord. While the CSF may “cushion” the spinal cord, pressure waves within the CSF at the time of injury may contribute to the extent and severity of the primary injury. The objective of this study was to develop a model of contusion SCI in a miniature pig and establish the feasibility of measuring spinal CSF pressure during injury.

Methods

A custom weight-drop device was used to apply thoracic contusion SCI to 17 Yucatan miniature pigs. Impact load and velocity were measured. Using fiber optic pressure transducers implanted in the thecal sac, CSF pressures resulting from 2 injury severities (caused by 50-g and 100-g weights released from a 50-cm height) were measured.

Results

The median peak impact loads were 54 N and 132 N for the 50-g and 100-g injuries, respectively. At a nominal 100 mm from the injury epicenter, the authors observed a small negative pressure peak (median −4.6 mm Hg [cranial] and −5.8 mm Hg [caudal] for 50 g; −27.6 mm Hg [cranial] and −27.2 mm Hg [caudal] for 100 g) followed by a larger positive pressure peak (median 110.5 mm Hg [cranial] and 77.1 mm Hg [caudal] for 50 g; 88.4 mm Hg [cranial] and 67.2 mm Hg [caudal] for 100 g) relative to the preinjury pressure. There were no significant differences in peak pressure between the 2 injury severities or the caudal and cranial transducer locations.

Conclusions

A new model of contusion SCI was developed to measure spinal CSF pressures during the SCI event. The results suggest that the Yucatan miniature pig is an appropriate model for studying CSF, spinal cord, and dura interactions during injury. With further development and characterization it may be an appropriate in vivo largeanimal model of SCI to answer questions regarding pathological changes, therapeutic safety, or treatment efficacy, particularly where humanlike dimensions and physiology are important.

Maximum principal strain correlates with spinal cord tissue damage in contusion and dislocation injuries in the rat cervical spine

The heterogeneity of the primary mechanical mechanism of spinal cord injury (SCI) is not currently used to tailor treatment strategies because the effects of these distinct patterns of acute mechanical damage on long-term neuropathology have not been fully investigated. A computational model of SCI enables the dynamic analysis of mechanical forces and deformations within the spinal cord tissue that would otherwise not be visible from histological tissue sections. We created a dynamic, three-dimensional finite element (FE) model of the rat cervical spine and simulated contusion and dislocation SCI mechanisms. We investigated the relationship between maximum principal strain and tissue damage, and compared primary injury patterns between mechanisms. The model incorporated the spinal cord white and gray matter, the dura mater, cerebrospinal fluid, spinal ligaments, intervertebral discs, a rigid indenter and vertebrae, and failure criteria for ligaments and vertebral endplates. High-speed (∼ 1 m/sec) contusion and dislocation injuries were simulated between vertebral levels C3 and C6 to match previous animal experiments, and average peak maximum principal strains were calculated for several regions at the injury epicenter and at 1-mm intervals from +5 mm rostral to -5 mm caudal to the lesion. Average peak principal strains were compared to tissue damage measured previously in the same regions via axonal permeability to 10-kD fluorescein-dextran. Linear regression of tissue damage against peak maximum principal strain for pooled data within all white matter regions yielded similar and significant (p<0.0001) correlations for both contusion (R(2)=0.86) and dislocation (R(2)=0.52). The model enhances our understanding of the differences in injury patterns between SCI mechanisms, and provides further evidence for the link between principal strain and tissue damage.

The importance of fluid-structure interaction in spinal trauma models

While recent studies have demonstrated the importance of the initial mechanical insult in the severity of spinal cord injury, there is a lack of information on the detailed cord-column interaction during such events. In vitro models have demonstrated the protective properties of the cerebrospinal fluid, but visualization of the impact is difficult. In this study a computational model was developed in order to clarify the role of the cerebrospinal fluid and provide a more detailed picture of the cord-column interaction. The study was validated against a parallel in vitro study on bovine tissue. Previous assumptions about complete subdural collapse before any cord deformation were found to be incorrect. Both the presence of the dura mater and the cerebrospinal fluid led to a reduction in the longitudinal strains within the cord. The division of the spinal cord into white and grey matter perturbed the bone fragment trajectory only marginally. In conclusion, the cerebrospinal fluid had a significant effect on the deformation pattern of the cord during impact and should be included in future models. The type of material models used for the spinal cord and the dura mater were found to be important to the stress and strain values within the components, but less important to the fragment trajectory.