The mechanical properties of rat spinal cord in vitro

Bioprinted Microchannel Scaffolds Modulate Neuronal Differentiation of Encapsulated Human Spinal Cord Progenitor Cells

A potential approach for treating spinal cord injuries is the implantation of human induced pluripotent stem cells-derived spinal cord progenitor cells (SCPCs) encapsulated in hydrogels. Digital light processing (DLP) enables the fabrication of scaffolds with high microchannel packing density, which are essential for neurofilament infiltration. In this study, SCPCs were encapsulated in gelatin methacrylate (GelMA)-based bioinks for single-layer printing via DLP bioprinting to incorporate human SCPCs within microchannel scaffolds at a reduced printing time. Mechanical properties were evaluated through degradation studies and compression testing, revealing that while the presence of poly(ethylene glycol) diacrylate improved printability and scaffold stability, it adversely affected cell survival. Scaffolds with higher GelMA concentration (10%) induced greater extent of motor neuronal differentiation as compared to those with 7.5% GelMA concentration (9.4 ± 5.1% vs 3.70 ± 2.6%, p < 0.001). In contrast, the scaffolds with lower GelMA concentration increased interneuron differentiation compared to those with higher GelMA concentration (7.3 ± 1.7% vs 1.6 ± 1.8%, p < 0.01), indicating that stiffness and GelMA content may modulate SCPC differentiation to specific neural subtypes. Overall, the encapsulation of SCPCs within the GelMA microchannel scaffold highlights the significance of material composition and stiffness in 3D printability and neuronal differentiation for spinal cord injury treatment.

Mechanical characteristics of spinal cord tissue by indentation

The mechanical properties of brain and spinal cord tissue have proven to be extremely complex and difficult to assess. Due to the heterogeneous and ultra-soft nature of the tissue, the available literature shows a large variance in mechanical parameters derived from experiments. In this study, we performed a series of indentation experiments to systematically investigate the mechanical properties of porcine spinal cord tissue in terms of their sensitivity to indentation tip diameter, loading rate, holding time, ambient temperature along with cyclic and oscillatory dynamic loading. Our results show that spinal cord white matter tissue is more compliant than grey matter tissue with apparent moduli of 128.7 and 403.8 Pa, respectively. They show similar viscoelastic behavior with stress relaxation time constants of τ1=1.38s and τ2=36.29s for grey matter and τ1=1.46s and τ2=46.10s for white matter, while the initial peak force decreased by 54 % for grey and 59 % for white matter tissue. An increase of the applied loading rate by two orders of magnitude led to an approximate doubling of the apparent modulus for both tissue types. Thermal variations showed a decrease in apparent modulus of up to 30 % after heating from 20 to 37.0 °C. Our dynamic tests revealed a significant influence of cyclic preload on the stiffness, with a drop of up to 20 % and a relative decrease of up to 60 % after the first cycle compared to the total modulus drop after five cycles for spinal cord grey matter tissue. Oscillatory indentation experiments identified similar loss moduli for spinal cord grey and white matter tissue and a higher storage modulus for white matter tissue. This work provides systematic insights into the mechanical properties of spinal cord tissue under different loading scenarios using nanoindentation.

Regional biomechanical characterization of the spinal cord tissue: dynamic mechanical response

Characterizing the dynamic mechanical properties of spinal cord tissue is deemed important for developing a comprehensive knowledge of the mechanisms underlying spinal cord injury. However, complex viscoelastic properties are vastly underexplored due to the spinal cord shows heterogeneous properties. To investigate regional differences in the biomechanical properties of spinal cord, we provide a mechanical characterization method (i.e., dynamic mechanical analysis) that facilitates robust measurement of spinal cord ex vivo, at small deformations, in the dynamic regimes. Load-unload cycles were applied to the tissue surface at sinusoidal frequencies of 0.05, 0.10, 0.50 and 1.00 Hz ex vivo within 2 h post mortem. We report the main response features (e.g., nonlinearities, rate dependencies, hysteresis and conditioning) of spinal cord tissue dependent on anatomical origin, and quantify the viscoelastic properties through the measurement of peak force, moduli, and hysteresis and energy loss. For all three anatomical areas (cervical, thoracic, and lumbar spinal cord tissues), the compound, storage, and loss moduli responded similarly to increasing strain rates. Notably, the complex modulus values of ex vivo spinal cord tissue rose nonlinearly with rising test frequency. Additionally, at every strain rate, it was shown that the tissue in the thoracic spinal cord was significantly more rigid than the tissue in the cervical or lumbar spinal cord, with compound modulus values roughly 1.5-times that of the lumbar region. At strain rates between 0.05 and 0.50 Hz, tan δ values for thoracic (that is, 0.26, 0.25, 0.06, respectively) and lumbar (that is, 0.27, 0.25, 0.07, respectively) spinal cord regions were similar, respectively, which were higher than cervical (that is, 0.21, 0.21, 0.04, respectively) region. The conditioning effects tend to be greater at relative higher deformation rates. Interestingly, no marked difference of conditioning ratios is observed among all three anatomical regions, regardless of loading rate. These findings lay a foundation for further comparison between healthy and diseased spinal cord to the future development of spinal cord scaffold and helps to advance our knowledge of neuroscience.

Biomimetic Dual-Network Collagen Fibers with Porous and Mechanical Cues Reconstruct Neural Stem Cell Niche via AKT/YAP Mechanotransduction after Spinal Cord Injury

Tissue engineering scaffolds can mediate the maneuverability of neural stem cell (NSC) niche to influence NSC behavior, such as cell self-renewal, proliferation, and differentiation direction, showing the promising application in spinal cord injury (SCI) repair. Here, dual-network porous collagen fibers (PCFS) are developed as neurogenesis scaffolds by employing biomimetic plasma ammonia oxidase catalysis and conventional amidation cross-linking. Following optimizing the mechanical parameters of PCFS, the well-matched Young's modulus and physiological dynamic adaptability of PCFS (4.0 wt%) have been identified as a neurogenetic exciter after SCI. Remarkably, porous topographies and curving wall-like protrusions are generated on the surface of PCFS by simple and non-toxic CO2 bubble-water replacement. As expected, PCFS with porous and matched mechanical properties can considerably activate the cadherin receptor of NSCs and induce a series of serine-threonine kinase/yes-associated protein mechanotransduction signal pathways, encouraging cellular orientation, neuron differentiation, and adhesion. In SCI rats, implanted PCFS with matched mechanical properties further integrated into the injured spinal cords, inhibited the inflammatory progression and decreased glial and fibrous scar formation. Wall-like protrusions of PCFS drive multiple neuron subtypes formation and even functional neural circuits, suggesting a viable therapeutic strategy for nerve regeneration and functional recovery after SCI.

The mechanical properties of the spinal cord: a systematic review

BACKGROUND CONTENT

Spinal cord compression is a source of pathology routinely seen in clinical practice. However, there remain unanswered questions surrounding both the understanding of pathogenesis and the best method of treatment. This arises from limited real-life testing of the mechanical properties of the spinal cord, either through cadaveric human specimens or animal testing, both of which suffer from methodological, as well as ethical, issues.

PURPOSE

To conduct a review of the literature on the mechanical properties of the spinal cord.

STUDY DESIGN/SETTING

A systematic review of the literature on the mechanical properties of the spinal cord is undertaken.

PATIENT SAMPLE

All literature reporting the testing of the mechanical properties of the spinal cord.

OUTCOME MEASURES

Reported physiological mechanical properties of the spinal cord.

METHODS

The methodological quality of the studies has been assessed within the ARRIVE guidelines using the CAMARADES framework and SYRCLE's risk of bias tool. This paper details the methodologies and results of the reported testing.

RESULTS

We show that (1) the research quality of previous work does not follow published guidelines on animal treatment or risk of bias, (2) no standard protocol has been employed for sample preparation or mechanical testing, (3) this leads to a wide distribution of results for the tested mechanical properties, not applicable to the living human or animal, and (4) animal testing is not a good proxy for human application.

CONCLUSIONS

The findings summarize the sum of current knowledge inherent to the mechanical properties of the spinal cord and may contribute to the development of a physical model which is applicable to the living human for analysis and testing in a controlled and repeatable fashion. Such a model would be the basis for further clinical research to improve outcomes from spinal cord compression.

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.

The mechanical properties of the spinal cord: a protocol for a systematic review of previous testing procedures and results

BACKGROUND

Spinal cord compression is a pathology seen in routine clinical practice. However, there remain a number of unanswered questions around both the understanding of the pathogenesis and the best method of treatment of the condition. This is partly due to the issues of the real-life testing of the physical properties of the spinal cord, either through the use of cadaveric human specimens or through animal testing, both of which have methodological, as well as ethical, issues.

DESIGN AND METHODS

This paper details a protocol for a systematic review of the literature on the mechanical properties of the spinal cord. We will conduct a literature search of a number of electronic databases, along with the grey literature, as a single-stage search. All literature will be screened for appropriate studies which will then be reviewed fully to extract relevant information on the methodology and mechanics of the reported testing along with the results. Two reviewers will separately screen and extract the data, with a comparison of results to ensure concordance. Conflicts will be resolved through discussion and independent arbitration as required. The methodological quality of the studies will be assessed within the ARRIVE guidelines using the CAMARADES framework and SYRCLE risk of bias tool. A narrative synthesis will be created with the appropriate tables to describe the demographics and findings of the included studies.

DISCUSSION

The systematic review described here will form the basis of an understanding of the current literature around the physical properties of the spinal cord. This will allow future work to develop a physical model of the spinal cord, which is translatable to patients for analysis and testing in a controlled and repeatable fashion. Such a model would be the basis for further clinical research to improve outcomes from this condition.

TRIAL REGISTRATION

Prospero registration number: CRD42022361933.

A finite element model of contusion spinal cord injury in rodents

Traumatic spinal cord injuries result from high impact forces acting on the spine and are proceeded by an extensive secondary inflammatory response resulting in motor, sensory, and autonomic dysfunction. Experimental in vivo traumatic spinal cord injuries in rodents using a contusion model have been extremely useful in elucidating the underlying pathophysiology of these injuries. However, the relationship between the pathophysiology and the biomechanical factors is still not well understood. Therefore, the aim of this research is to provide a comprehensive analysis of the biomechanics of traumatic spinal cord injury in a rat contusion model. This is achieved through the development and validation of a finite element model of the thoracic rat spinal cord and subsequently simulating controlled cortical impact-induced traumatic spinal cord injury. The effects of impactor velocity, depth, and geometry on the resulting stresses and strains within the spinal cord are investigated. Our results show that increasing impactor depth results in larger stresses and strains within the spinal cord tissue as expected. Further, for the first time ever our results show that impactor geometry (spherical versus cylindrical) plays an important role in the distribution and magnitude of stresses and strains within the cord. Therefore, finite element modelling can be a powerful tool used to predict stresses and strains that occur in spinal cord tissue during trauma.

The mechanical behavior of bovine spinal cord white matter under various strain rate conditions: tensile testing and visco-hyperelastic constitutive modeling

The mechanical behavior of the white matter is important for estimating the damage of the spinal cord during accidents. In this study, we conducted uniaxial tension testing in vitro of bovine spinal cord white matter under extremely high strain rate conditions (up to 100 s^-1). A visco-hyperelastic constitutive law for modeling the strain rate-dependent behavior of the bovine spinal cord white matter was developed. A set of material constants was obtained using a Levenberg-Marquardt fitting algorithm to match the uniaxial tension experimental data with various strain rates. Our experimental data confirmed that the modulus and tensile strength increased when the strain rate is higher. For the extremely high strain rate condition (100 s^-1), we found that both the modulus and failure stress significantly increased compared with the low strain rate case. These new data in terms of mechanical response at high strain rate provide insight into the spine injury mechanism caused by high-speed impact. Moreover, the developed constitutive model will allow researchers to perform more realistic finite element modeling and simulation of spinal cord injury damage under various complicated conditions.

Dynamic changes in mechanical properties of the adult rat spinal cord after injury

Spinal cord injury (SCI), a debilitating medical condition that can cause irreversible loss of neurons and permanent paralysis, currently has no cure. However, regenerative medicine may offer a promising treatment. Given that numerous regenerative strategies aim to deliver cells and materials in the form of tissue-engineered therapies, understanding and characterising the mechanical properties of the spinal cord tissue is very important. In this study, we have systematically characterised the spatiotemporal changes in elastic stiffness (elastic modulus, Pa) and viscosity (drop in peak force, %) of injured rat thoracic spinal cord tissues at distinct time points after crush injury using the indentation technique. Our results demonstrate that in comparison with uninjured spinal cord tissue, the injured tissues exhibited lower stiffness (median 3281 Pa versus 9632 Pa; P < 0.001) but demonstrated elevated viscosity (median 80% versus 57%; P < 0.001) at 3 days postinjury. Between 4 and 6 weeks after SCI, the overall viscoelastic properties of injured tissues returned to baseline values. At 12 weeks after SCI, in comparison with uninjured tissue, the injured spinal cord tissues displayed a significant increase in both elasticity (median 13698 Pa versus 9920 Pa; P < 0.001) and viscosity (median 64% versus 58%; P < 0.001). This work constitutes the first quantitative mapping of spatiotemporal changes in spinal cord tissue elasticity and viscosity in injured rats, providing a mechanical basis of the tissue for future studies on the development of biomaterials for SCI repair. STATEMENT OF SIGNIFICANCE: Spinal cord injury (SCI) is a devastating disease often leading to permanent paralysis. While enormous progress in understanding the molecular pathomechanisms of SCI has been made, the mechanical properties of injured spinal cord tissue have received considerably less attention. This study provides systematic characterization of the biomechanical evolution of rat spinal cord tissue after SCI using a microindentation test method. We find spinal cord tissue behaves significantly softer but more viscous immediately postinjury. As time passes, the lesion site gradually returns to baseline values and then displays pronounced increased viscoelastic properties. As host tissue mechanical properties are a crucial consideration for any biomaterial implanted into central nervous system, our results may have important implications for further studies of SCI repair.

Composition, host responses and clinical applications of bioadhesives

Bioadhesives are medical devices used to join or seal tissues that have been injured or incised. They have been classified into tissue adhesives, sealants, and hemostatic agents. Bioadhesives such as FloSeal®, CoSeal®, BioGlue®, Evicel®, Tisseel®, Progel™ PALS, and TissuGlu® have been commercialized and used in clinical setting. They can be formulated with natural or synthetic components or a combination of both including albumin, glutaraldehyde, chitosan, cyanoacrylate, fibrin and thrombin, gelatin, polyethylene glycol (PEG), along with urethanes. Each formulation has intrinsic properties and has been developed and validated for a specific application. This review article briefs the mechanisms by which bioadhesives forms adhesion to tissues and highlights the correlation between bioadhesives composition and their potential host responses. Furthermore, clinical applications of bioadhesives and their application-driven requirements are outlined.

Directional Submicrofiber Hydrogel Composite Scaffolds Supporting Neuron Differentiation and Enabling Neurite Alignment

Cell cultures aiming at tissue regeneration benefit from scaffolds with physiologically relevant elastic moduli to optimally trigger cell attachment, proliferation and promote differentiation, guidance and tissue maturation. Complex scaffolds designed with guiding cues can mimic the anisotropic nature of neural tissues, such as spinal cord or brain, and recall the ability of human neural progenitor cells to differentiate and align. This work introduces a cost-efficient gelatin-based submicron patterned hydrogel–fiber composite with tuned stiffness, able to support cell attachment, differentiation and alignment of neurons derived from human progenitor cells. The enzymatically crosslinked gelatin-based hydrogels were generated with stiffnesses from 8 to 80 kPa, onto which poly(ε-caprolactone) (PCL) alignment cues were electrospun such that the fibers had a preferential alignment. The fiber–hydrogel composites with a modulus of about 20 kPa showed the strongest cell attachment and highest cell proliferation, rendering them an ideal differentiation support. Differentiated neurons aligned and bundled their neurites along the aligned PCL filaments, which is unique to this cell type on a fiber–hydrogel composite. This novel scaffold relies on robust and inexpensive technology and is suitable for neural tissue engineering where directional neuron alignment is required, such as in the spinal cord.

Effect of Velocity and Contact Stress Area on the Dynamic Behavior of the Spinal Cord Under Different Testing Conditions

Knowledge of the dynamic behavior of the spinal cord under different testing conditions is critical for our understanding of biomechanical mechanisms of spinal cord injury. Although velocity and contact stress area are known to affect external mechanical stress or energy upon sudden traumatic injury, quantitative investigation of the two clinically relevant biomechanical variables is limited. Here, freshly excised rat spinal-cord–pia-arachnoid constructs were tested through indentation using indenters of different sizes (radii: 0.25, 0.50, and 1.00 mm) at various loading rates ranging from 0.04 to 0.20 mm/s. This analysis found that the ex vivo specimen displayed significant nonlinear viscoelasticity at <10% of specimen thickness depth magnitudes. At higher velocity and larger contact stress area, the cord withstood a higher peak load and exhibited more sensitive mechanical relaxation responses (i.e., increasing amplitude and speed of the drop in peak load). Additionally, the cord became stiffer (i.e., increasing elastic modulus) and softer (i.e., decreasing elastic modulus) at a higher velocity and larger contact stress area, respectively. These findings will improve our understanding of the real-time complex biomechanics involved in traumatic spinal cord injury.

Correlating Tissue Mechanics and Spinal Cord Injury: Patient-Specific Finite Element Models of Unilateral Cervical Contusion Spinal Cord Injury in Non-Human Primates

Non-human primate (NHP) models are the closest approximation of human spinal cord injury (SCI) available for pre-clinical trials. The NHP models, however, include broader morphological variability that can confound experimental outcomes. We developed subject-specific finite element (FE) models to quantify the relationship between impact mechanics and SCI, including the correlations between FE outcomes and tissue damage. Subject-specific models of cervical unilateral contusion SCI were generated from pre-injury MRIs of six NHPs. Stress and strain outcomes were compared with lesion histology using logit analysis. A parallel generic model was constructed to compare the outcomes of subject-specific and generic models. The FE outcomes were correlated more strongly with gray matter damage (0.29 < R2 < 0.76) than white matter (0.18 < R2 < 0.58). Maximum/minimum principal strain, Von-Mises and Tresca stresses showed the strongest correlations (0.31 < R2 < 0.76) with tissue damage in the gray matter while minimum principal strain, Von-Mises stress, and Tresca stress best predicted white matter damage (0.23 < R2 < 0.58). Tissue damage thresholds varied for each subject. The generic FE model captured the impact biomechanics in two of the four models; however, the correlations between FE outcomes and tissue damage were weaker than the subject-specific models (gray matter [0.25 < R2 < 0.69] and white matter [R2 < 0.06] except for one subject [0.26 < R2 < 0.48]). The FE mechanical outputs correlated with tissue damage in spinal cord white and gray matters, and the subject-specific models accurately mimicked the biomechanics of NHP cervical contusion impacts.

Mechanical properties of spinal cord grey matter and white matter in confined compression

To better understand the link between spinal cord impact and the resulting tissue damage, computational models are often used. These models typically simulate the spinal cord as a homogeneous and isotropic material. Recent research suggests that grey and white matter tissue differences and directional differences, i.e. anisotropy, are important to predict spinal cord damage. The objective of this research was to characterize the mechanical properties of spinal cord grey and white matter tissue in confined compression. Spinal cords (n = 12) were harvested immediately following euthanasia from Yorkshire and Yucatan pigs. The spinal cords were flash frozen (60 s at -80 °C) and prepared into four types of test samples: grey matter axial, grey matter transverse, white matter axial, and white matter transverse. Each sample type was thawed, and subsequently tested in confined compression within 6 h of euthanasia. Samples were compressed to 10% strain at a quasi-static strain rate (0.001/sec) and allowed to relax for 120 s. A quasi-linear viscoelastic model combining a first-order exponential with a 1-term Prony series characterized the loading and relaxation responses respectively. The effect of tissue type (grey matter vs. white matter), direction (axial vs. transverse), and their interaction were evaluated with a two-way ANOVA (p < 0.05) with peak stress, aggregate modulus, and relaxation time as dependent variables. This study found grey matter to be 1.6-2 times stiffer than white matter and both grey and white matter were isotropic in compression. These findings should be emphasized when studying SCI biomechanics using computational models.

Fabrication of a biomimetic spinal cord tissue construct with heterogenous mechanical properties using intrascaffold cell assembly

In tissue engineering studies, scaffolds play a very important role in offering both physical and chemical cues for cell growth and tissue regeneration. However, in some cases, tissue regeneration requires scaffolds with high mechanical properties (e.g., bone and cartilage), while cells need a soft mechanical microenvironment. In this study, to mimic the heterogenous mechanical properties of a spinal cord tissue, a biomimetic rat tissue construct is fabricated. A collagen-coated poly(lactic-co-glycolic acid) scaffold is manufactured using thermally induced phase separation casting. Primary rat neural cells (P01 Wistar rat cortex) with soft hydrogels are later printed within the scaffold using an image-guided intrascaffold cell assembly technique. The scaffolds have unidirectional microporous structure with parallel axial macrochannels (260 ± 4 µm in diameter). Scaffolds showed mechanical properties similar to rat spine (ultimate tensile strength: 0.085 MPa, Young's modulus [stretch]: 0.31 MPa). The bioink composed of gelatin/alginate/fibrinogen is precisely printed into the macrochannels and showed mechanical properties suitable for neural cells (Young's modulus [compressive]: 3.814 kPa). Scaffold interface, cell viability, and immunostaining analyses show uniform distribution of stable, healthy, and elongated neural cells and neurites over 14 culture days in vitro. The results demonstrated that this method can serve as a valuable tool to aid manufacturing of tissue constructs requiring heterogenous mechanical properties for complex cell and/or biomolecule assembly.

Basic biomechanics of spinal cord injury - How injuries happen in people and how animal models have informed our understanding

The wide variability, or heterogeneity, in human spinal cord injury is due partially to biomechanical factors. This review summarizes our current knowledge surrounding the patterns of human spinal column injury and the biomechanical factors affecting injury. The biomechanics of human spinal injury is studied most frequently with human cadaveric models and the features of the two most common injury patterns, burst fracture and fracture dislocation, are outlined. The biology of spinal cord injury is typically studied with animal models and the effects of the most relevant biomechanical factors - injury mechanism, injury velocity, and residual compression, are described. Tissue damage patterns and behavioural outcomes following dislocation or distraction injury mechanisms differ from the more commonly used contusion mechanism. The velocity of injury affects spinal cord damage, principally in the white matter. Ongoing, or residual compression after the initial impact does affect spinal cord damage, but few models exist that replicate the clinical scenario. Future research should focus on the effects of these biomechanical factors in different preclinical animal models as recent data suggests that treatment outcomes may vary between models.

The influence of microenvironment and extracellular matrix molecules in driving neural stem cell fate within biomaterials

Transplantation of stem cells is a promising potential therapy for central nervous system disease and injury. The capacity for self-renewal, proliferation of progenitor cells, and multi-lineage potential underscores the need for controlling stem cell fate. Furthermore, transplantation within a hostile environment can lead to significant cell death and limited therapeutic potential. Tissue-engineered materials have been developed to both regulate stem cell fate, increase transplanted cell viability, and improve therapeutic outcomes. Traditionally, regulation of stem cell differentiation has been driven through soluble signals, such as growth factors. While these signals are important, insoluble factors from the local microenvironment or extracellular matrix (ECM) molecules also contribute to stem cell activity and fate. Understanding the microenvironment factors that influence stem cell fate, such as mechanical properties, topography, and presentation of specific ECM ligands, is necessary for designing improved biomaterials. Here we review some of the microenvironment factors that regulate stem cell fate and how they can be incorporated into biomaterials as part of potential CNS therapies.

Comparing Predictive Accuracy and Computational Costs for Viscoelastic Modeling of Spinal Cord Tissues

The constitutive equation used to characterize and model spinal tissues can significantly influence the conclusions from experimental and computational studies. Therefore, researchers must make critical judgements regarding the balance of computational efficiency and predictive accuracy necessary for their purposes. The objective of this study is to quantitatively compare the fitting and prediction accuracy of linear viscoelastic (LV), quasi-linear viscoelastic (QLV), and (fully) non-linear viscoelastic (NLV) modeling of spinal-cord-pia-arachnoid-construct (SCPC), isolated cord parenchyma, and isolated pia-arachnoid-complex (PAC) mechanics in order to better inform these judgements. Experimental data collected during dynamic cyclic testing of each tissue condition were used to fit each viscoelastic formulation. These fitted models were then used to predict independent experimental data from stress-relaxation testing. Relative fitting accuracy was found not to directly reflect relative predictive accuracy, emphasizing the need for material model validation through predictions of independent data. For the SCPC and isolated cord, the NLV formulation best predicted the mechanical response to arbitrary loading conditions, but required significantly greater computational run time. The mechanical response of the PAC under arbitrary loading conditions was best predicted by the QLV formulation.

High-Speed Fluoroscopy to Measure Dynamic Spinal Cord Deformation in an In Vivo Rat Model

Although spinal cord deformation is thought to be a predictor of injury severity, few researchers have investigated dynamic cord deformation, in vivo, during impact. This is needed to establish correlations among impact parameters, internal cord deformation, and histological and functional outcomes. Relying on surface deformations alone may not sufficiently represent spinal cord deformation. The objective of this study was to develop a high-speed fluoroscopic method of tracking the surface and internal cord deformations of rat spinal cord during experimental cord injury. Two radio-opaque beads were injected into the cord at C5/6 in the dorsal and ventral white matter. Four additional beads were glued to the surface of the cord. Dynamic bead displacement was tracked during a dorsal impact (130 mm/sec, 1 mm depth) by high-speed radiographic imaging at 3000 FPS, laterally. The internal spinal cord beads displaced significantly more than the surface beads in the ventral direction (1.1-1.9 times) and more than most surface beads in the cranial direction (1.2-1.5 times). The dorsal beads (internal and surface) displaced more than the ventral beads during all impacts. The bead displacement pattern implies that the spinal cord undergoes complex internal and surface deformations during impact. Residual displacement of the internal beads was significantly greater than that of the surface beads in the cranial-caudal direction but not the dorsoventral direction. Finite element simulation confirmed that the additional bead mass likely had little effect on the internal cord deformations. These results support the merit of this technique for measuring in vivo spinal cord deformation.

Viscoelasticity of spinal cord and meningeal tissues

Compared to the outer dura mater, the mechanical behavior of spinal pia and arachnoid meningeal layers has received very little attention in the literature. This is despite experimental evidence of their importance with respect to the overall spinal cord stiffness and recovery following compression. Accordingly, inclusion of the mechanical contribution of the pia and arachnoid maters would improve the predictive accuracy of finite element models of the spine, especially in the distribution of stresses and strain through the cord's cross-section. However, to-date, only linearly elastic moduli for what has been previously identified as spinal pia mater is available in the literature. This study is the first to quantitatively compare the viscoelastic behavior of isolated spinal pia-arachnoid-complex, neural tissue of the spinal cord parenchyma, and intact construct of the two. The results show that while it only makes up 5.5% of the overall cross-sectional area, the thin membranes of the innermost meninges significantly affect both the elastic and viscous response of the intact construct. Without the contribution of the pia and arachnoid maters, the spinal cord has very little inherent stiffness and experiences significant relaxation when strained. The ability of the fitted non-linear viscoelastic material models of each condition to predict independent data within experimental variability supports their implementation into future finite element computational studies of the spine.

STATEMENT OF SIGNIFICANCE

The neural tissue of the spinal cord is surrounded by three fibrous layers called meninges which are important in the behavior of the overall spinal-cord-meningeal construct. While the mechanical properties of the outermost layer have been reported, the pia mater and arachnoid mater have received considerably less attention. This study is the first to directly compare the behavior of the isolated neural tissue of the cord, the isolated pia-arachnoid complex, and the construct of these individual components. The results show that, despite being very thin, the inner meninges significantly affect the elastic and time-dependent response of the spinal cord, which may have important implications for studies of spinal cord injury.

A ternary nanofibrous scaffold potential for central nerve system tissue engineering

In the present research, a ternary polycaprolactone (PCL)/gelatin/fibrinogen nanofibrous scaffold for tissue engineering application was developed. Through this combination, PCL improved the scaffold mechanical properties; meanwhile, gelatin and fibrinogen provided more hydrophilicity and cell proliferation. Three types of nanofibrous scaffolds containing different fibrinogen contents were prepared and characterized. Morphological study of the nanofibers showed that the prepared nanofibers were smooth, uniform without any formation of beads with a significant reduction in nanofiber diameter after incorporation of fibrinogen. The chemical characterization of the scaffolds confirmed that no chemical reaction occurred between the scaffold components. The tensile test results of the scaffolds showed that increasing in fibrinogen content led to a decrease in mechanical properties. Furthermore, adipose-derived stem cells were employed to evaluate cell-scaffold interaction. Cell culture results indicated that higher cell proliferation occurred for the higher amount of fibrinogen. Statistical analysis was also carried out to evaluate the significant difference for the obtained results of water droplet contact angle and cell culture. Therefore, the results confirmed that PCL/gel/fibrinogen scaffold has a good potential for tissue engineering applications including central nerve system tissue engineering. © 2018 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 106A:2394-2401, 2018.

Radiography used to measure internal spinal cord deformation in an in vivo rat model

Little is known about the internal mechanics of the in vivo spinal cord during injury. The objective of this study was to develop a method of tracking internal and surface deformation of in vivo rat spinal cord during compression using radiography. Since neural tissue is radio-translucent, radio-opaque markers were injected into the spinal cord. Two tantalum beads (260 µm) were injected into the cord (dorsal and ventral) at C5 of nine anesthetized rats. Four beads were glued to the lateral surface of the cord, caudal and cranial to the injection site. A compression plate was displaced 0.5 mm, 2 mm, and 3 mm into the spinal cord and lateral X-ray images were taken before, during, and after each compression for measuring bead displacements. Potential bead migration was monitored for by comparing displacements of the internal and glued surface beads. Dorsal beads moved significantly more than ventral beads with a range in averages of 0.57-0.71 mm and 0.31-0.35 mm respectively. Bead displacements during 0.5 mm compressions were significantly lower than 2 mm and 3 mm compressions. There was no statistically significant migration of the internal beads. The results indicate the merit of this technique for measuring in vivo spinal cord deformation. The pattern of bead displacements illustrates the complex internal and surface deformations of the spinal cord during transverse compression. This information is needed for validating physical and finite element spinal cord surrogates and to define relationships between loading parameters, internal cord deformation, and biological and functional outcomes.

Comparison of in vivo and ex vivo viscoelastic behavior of the spinal cord

Despite efforts to simulate the in vivo environment, post-mortem degradation and lack of blood perfusion complicate the use of ex vivo derived material models in computational studies of spinal cord injury. In order to quantify the mechanical changes that manifest ex vivo, the viscoelastic behavior of in vivo and ex vivo porcine spinal cord samples were compared. Stress-relaxation data from each condition were fit to a non-linear viscoelastic model using a novel characterization technique called the direct fit method. To validate the presented material models, the parameters obtained for each condition were used to predict the respective dynamic cyclic response. Both ex vivo and in vivo samples displayed non-linear viscoelastic behavior with a significant increase in relaxation with applied strain. However, at all three strain magnitudes compared, ex vivo samples experienced a higher stress and greater relaxation than in vivo samples. Significant differences between model parameters also showed distinct relaxation behaviors, especially in non-linear relaxation modulus components associated with the short-term response (0.1-1 s). The results of this study underscore the necessity of utilizing material models developed from in vivo experimental data for studies of spinal cord injury, where the time-dependent properties are critical. The ability of each material model to accurately predict the dynamic cyclic response validates the presented methodology and supports the use of the in vivo model in future high-resolution finite element modeling efforts.

STATEMENT OF SIGNIFICANCE

Neural tissues (such as the brain and spinal cord) display time-dependent, or viscoelastic, mechanical behavior making it difficult to model how they respond to various loading conditions, including injury. Methods that aim to characterize the behavior of the spinal cord almost exclusively use ex vivo cadaveric or animal samples, despite evidence that time after death affects the behavior compared to that in a living animal (in vivo response). Therefore, this study directly compared the mechanical response of ex vivo and in vivo samples to quantify these differences for the first time. This will allow researchers to draw more accurate conclusions about spinal cord injuries based on ex vivo data (which are easier to obtain) and emphasizes the importance of future in vivo experimental animal work.

The development and validation of a numerical integration method for non-linear viscoelastic modeling

Compelling evidence that many biological soft tissues display both strain- and time-dependent behavior has led to the development of fully non-linear viscoelastic modeling techniques to represent the tissue’s mechanical response under dynamic conditions. Since the current stress state of a viscoelastic material is dependent on all previous loading events, numerical analyses are complicated by the requirement of computing and storing the stress at each step throughout the load history. This requirement quickly becomes computationally expensive, and in some cases intractable, for finite element models. Therefore, we have developed a strain-dependent numerical integration approach for capturing non-linear viscoelasticity that enables calculation of the current stress from a strain-dependent history state variable stored from the preceding time step only, which improves both fitting efficiency and computational tractability. This methodology was validated based on its ability to recover non-linear viscoelastic coefficients from simulated stress-relaxation (six strain levels) and dynamic cyclic (three frequencies) experimental stress-strain data. The model successfully fit each data set with average errors in recovered coefficients of 0.3% for stress-relaxation fits and 0.1% for cyclic. The results support the use of the presented methodology to develop linear or non-linear viscoelastic models from stress-relaxation or cyclic experimental data of biological soft tissues.

An Elastomeric Polymer Matrix, PEUU-Tac, Delivers Bioactive Tacrolimus Transdurally to the CNS in Rat

Central nervous system (CNS) neurons fail to regrow injured axons, often resulting in permanently lost neurologic function. Tacrolimus is an FDA-approved immunosuppressive drug with known neuroprotective and neuroregenerative properties in the CNS. However, tacrolimus is typically administered systemically and blood levels required to effectively treat CNS injuries can lead to lethal, off-target organ toxicity. Thus, delivering tacrolimus locally to CNS tissues may provide therapeutic control over tacrolimus levels in CNS tissues while minimizing off-target toxicity. Herein we show an electrospun poly(ester urethane) urea and tacrolimus elastomeric matrix (PEUU-Tac) can deliver tacrolimus trans-durally to CNS tissues. In an acute CNS ischemia model in rat, the optic nerve (ON) was clamped for 10s and then PEUU-Tac was used as an ON wrap and sutured around the injury site. Tacrolimus was detected in PEUU-Tac wrapped ONs at 24 h and 14 days, without significant increases in tacrolimus blood levels. Similar to systemically administered tacrolimus, PEUU-Tac locally decreased glial fibrillary acidic protein (GFAP) at the injury site and increased growth associated protein-43 (GAP-43) expression in ischemic ONs from the globe to the chiasm, consistent with decreased astrogliosis and increased retinal ganglion cell (RGC) axon growth signaling pathways. These initial results suggest PEUU-Tac is a biocompatible elastic matrix that delivers bioactive tacrolimus trans-durally to CNS tissues without significantly increasing tacrolimus blood levels and off-target toxicity.

•

PEUU-Tac locally delivers tacrolimus to CNS tissues

•

PEUU-Tac positively modulates CNS tissue remodeling

•

PEUU-Tac minimizes off-target tacrolimus toxicity

Central nervous system (CNS) injury typically results in permanently lost neurological function. Tacrolimus is an FDA-approved drug used during organ transplantation that also has CNS neuroprotective and neuroregenerative properties. However, tacrolimus is typically delivered systemically in the blood and delivering effective concentrations to CNS tissues requires tacrolimus blood levels that can lead to adverse side effects in multiple organs. Herein we show that PEUU-Tac, a tacrolimus-eluting matrix, can locally deliver tacrolimus to injured CNS tissues without increasing blood levels, suggesting PEUU-Tac can be used to treat CNS injuries locally while minimizing adverse side effects.

Damage Accumulation Modeling and Rate Dependency of Spinal Dura Mater

As the strongest of the meningeal tissues, the spinal dura mater plays an important role in the overall behavior of the spinal cord-meningeal complex (SCM). It follows that the accumulation of damage affects the dura mater's ability to protect the cord from excessive mechanical loads. Unfortunately, current computational investigations of spinal cord injury (SCI) etiology typically do not include postyield behavior. Therefore, a more detailed description of the material behavior of the spinal dura mater, including characterization of damage accumulation, is required to comprehensively study SCI. Continuum mechanics-based viscoelastic damage theories have been previously applied to other biological tissues; however, the current work is the first to report damage accumulation modeling in a tissue of the SCM complex. Longitudinal (i.e., cranial-to-caudal long-axis) samples of ovine cervical dura mater were tensioned-to-failure at one of three strain rates (quasi-static, 0.05/s, and 0.3/s). The resulting stress–strain data were fit to a hyperelastic continuum damage model to characterize the strain-rate-dependent subfailure and failure behavior. The results show that the damage behavior of the fibrous and matrix components of the dura mater are strain-rate dependent, with distinct behaviors when exposed to strain rates above that experienced during normal voluntary neck motion suggesting the possible existence of a protective mechanism.

Protection of cortex by overlying meninges tissue during dynamic indentation of the adolescent brain

Traumatic brain injury (TBI) has become a recent focus of biomedical research with a growing international effort targeting material characterization of brain tissue and simulations of trauma using computer models of the head and brain to try to elucidate the mechanisms and pathogenesis of TBI. The meninges, a collagenous protective tri-layer, which encloses the entire brain and spinal cord has been largely overlooked in these material characterization studies. This has resulted in a lack of accurate constitutive data for the cranial meninges, particularly under dynamic conditions such as those experienced during head impacts. The work presented here addresses this lack of data by providing for the first time, in situ large deformation material properties of the porcine dura-arachnoid mater composite under dynamic indentation. It is demonstrated that this tissue is substantially stiffer (shear modulus, μ=19.10±8.55kPa) and relaxes at a slower rate (τ₁=0.034±0.008s, τ₂=0.336±0.077s) than the underlying brain tissue (μ=6.97±2.26kPa, τ₁=0.021±0.007s, τ₂=0.199±0.036s), reducing the magnitudes of stress by 250% and 65% for strains that arise during indentation-type deformations in adolescent brains.

STATEMENT OF SIGNIFICANCE

We present the first mechanical analysis of the protective capacity of the cranial meninges using in situ micro-indentation techniques. Force-relaxation tests are performed on in situ meninges and cortex tissue, under large strain dynamic micro-indentation. A quasi-linear viscoelastic model is used subsequently, providing time-dependent mechanical properties of these neural tissues under loading conditions comparable to what is experienced in TBI. The reported data highlights the large differences in mechanical properties between these two tissues. Finite element simulations of the indentation experiments are also performed to investigate the protective capacity of the meninges. These simulations show that the meninges protect the underlying brain tissue by reducing the overall magnitude of stress by 250% and up to 65% for strains.

Biomechanical properties of the spinal cord: implications for tissue engineering and clinical translation

Spinal cord injury is a severely debilitating condition which can leave individuals paralyzed and suffering from autonomic dysfunction. Regenerative medicine may offer a promising solution to this problem. Previous research has focused primarily on exploring the cellular and biological aspects of the spinal cord, yet relatively little remains known about the biomechanical properties of spinal cord tissue. Given that a number of regenerative strategies aim to deliver cells and materials in the form of tissue-engineered therapies, understanding the biomechanical properties of host spinal cord tissue is important. We review the relevant biomechanical properties of spinal cord tissue and provide the baseline knowledge required to apply these important physical concepts to spinal cord tissue engineering.

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.

The Transverse Isotropy of Spinal Cord White Matter Under Dynamic Load

The rostral-caudally aligned fiber-reinforced structure of spinal cord white matter (WM) gives rise to transverse isotropy in the material. Stress and strain patterns generated in the spinal cord parenchyma following spinal cord injury (SCI) are multidirectional and dependent on the mechanism of the injury. Our objective was to develop a WM constitutive model that captures the material transverse isotropy under dynamic loading. The WM mechanical behavior was extracted from the published tensile and compressive experiments. Combinations of isotropic and fiber-reinforcing models were examined in a conditional quasi-linear viscoelastic (QLV) formulation to capture the WM mechanical behavior. The effect of WM transverse isotropy on SCI model outcomes was evaluated by simulating a nonhuman primate (NHP) contusion injury experiment. A second-order reduced polynomial hyperelastic energy potential conditionally combined with a quadratic reinforcing function in a four-term Prony series QLV model best captured the WM mechanical behavior (0.89 < R2 < 0.99). WM isotropic and transversely isotropic material models combined with discrete modeling of the pia mater resulted in peak impact forces that matched the experimental outcomes. The transversely isotropic WM with discrete pia mater resulted in maximum principal strain (MPS) distributions which effectively captured the combination of ipsilateral peripheral WM sparing, ipsilateral injury and contralateral sparing, and the rostral/caudal spread of damage observed in in vivo injuries. The results suggest that the WM transverse isotropy could have an important role in correlating tissue damage with mechanical measures and explaining the directional sensitivity of the spinal cord to injury.

CNS cell distribution and axon orientation determine local spinal cord mechanical properties

Mechanical signaling plays an important role in cell physiology and pathology. Many cell types, including neurons and glial cells, respond to the mechanical properties of their environment. Yet, for spinal cord tissue, data on tissue stiffness are sparse. To investigate the regional and direction-dependent mechanical properties of spinal cord tissue at a spatial resolution relevant to individual cells, we conducted atomic force microscopy (AFM) indentation and tensile measurements on acutely isolated mouse spinal cord tissue sectioned along the three major anatomical planes, and correlated local mechanical properties with the underlying cellular structures. Stiffness maps revealed that gray matter is significantly stiffer than white matter irrespective of directionality (transverse, coronal, and sagittal planes) and force direction (compression or tension) (K_g= ∼130 Pa vs. K_w= ∼70 Pa); both matters stiffened with increasing strain. When all data were pooled for each plane, gray matter behaved like an isotropic material under compression; however, subregions of the gray matter were rather heterogeneous and anisotropic. For example, in sagittal sections the dorsal horn was significantly stiffer than the ventral horn. In contrast, white matter behaved transversely isotropic, with the elastic stiffness along the craniocaudal (i.e., longitudinal) axis being lower than perpendicular to it. The stiffness distributions we found under compression strongly correlated with the orientation of axons, the areas of cell nuclei, and cellular in plane proximity. Based on these morphological parameters, we developed a phenomenological model to estimate local mechanical properties of central nervous system (CNS) tissue. Our study may thus ultimately help predicting local tissue stiffness, and hence cell behavior in response to mechanical signaling under physiological and pathological conditions, purely based on histological data.

Development of an integrated optical coherence tomography-gas nozzle system for surgical laser ablation applications: preliminary findings of in situ spinal cord deformation due to gas flow effects

Gas assisted laser machining of materials is a common practice in the manufacturing industry. Advantages in using gas assistance include reducing the likelihood of flare-ups in flammable materials and clearing away ablated material in the cutting path. Current surgical procedures and research do not take advantage of this and in the case for resecting osseous tissue, gas assisted ablation can help minimize charring and clear away debris from the surgical site. In the context of neurosurgery, the objective is to cut through osseous tissue without damaging the underlying neural structures. Different inert gas flow rates used in laser machining could cause deformations in compliant materials. Complications may arise during surgical procedures if the dura and spinal cord are damaged by these deformations. We present preliminary spinal deformation findings for various gas flow rates by using optical coherence tomography to measure the depression depth at the site of gas delivery.

Impact depth and the interaction with impact speed affect the severity of contusion spinal cord injury in rats

Spinal cord injury (SCI) biomechanics suggest that the mechanical factors of impact depth and speed affect the severity of contusion injury, but their interaction is not well understood. The primary aim of this work was to examine both the individual and combined effects of impact depth and speed in contusion SCI on the cervical spinal cord. Spinal cord contusions between C5 and C6 were produced in anesthetized rats at impact speeds of 8, 80, or 800 mm/s with displacements of 0.9 or 1.5 mm (n=8/group). After 7 days postinjury, rats were assessed for open-field behavior, euthanized, and spinal cords were harvested. Spinal cord tissue sections were stained for demyelination (myelin-based protein) and tissue sparing (Luxol fast blue). In parallel, a finite element model of rat spinal cord was used to examine the resulting maximum principal strain in the spinal cord during impact. Increasing impact depth from 0.9 to 1.5 mm reduced open-field scores (p<0.01) above 80 mm/s, reduced gray (GM) and white matter (WM) sparing (p<0.01), and increased the amount of demyelination (p<0.01). Increasing impact speed showed similar results at the 1.5-mm impact depth, but not the 0.9-mm impact depth. Linear correlation analysis with finite element analysis strain showed correlations (p<0.001) with nerve fiber damage in the ventral (R(2)=0.86) and lateral (R(2)=0.74) regions of the spinal cord and with WM (R(2)=0.90) and GM (R(2)=0.76) sparing. The results demonstrate that impact depth is more important in determining the severity of SCI and that threshold interactions exist between impact depth and speed.

Characterization of a Novel, Magnetic Resonance Imaging-Compatible Rodent Model Spinal Cord Injury Device

Rodent models of acute spinal cord injury (SCI) are often used to investigate the effects of injury mechanism, injury speed, and cord displacement magnitude, on the ensuing cascade of biological damage in the cord. However, due to its small size, experimental observations have largely been limited to the gross response of the cord. To properly understand the relationship between mechanical stimulus and biological damage, more information is needed about how the constituent tissues of the cord (i.e., gray and white matter) respond to injurious stimuli. To address this limitation, we developed a novel magnetic resonance imaging (MRI)-compatible test apparatus that can impose either a contusion-type or dislocation-type acute cervical SCI in a rodent model and facilitate MR-imaging of the cervical spinal cord in a 7 T MR scanner. In this study, we present the experimental performance parameters of the MR rig. Utilizing cadaveric specimens and static radiographs, we report contusion magnitude accuracy that for a desired 1.8 mm injury, a nominal 1.78 mm injury (SD = 0.12 mm) was achieved. High-speed video analysis was employed to determine the injury speeds for both mechanisms and were found to be 1147 mm/s (SD = 240 mm/s) and 184 mm/s (SD = 101 mm/s) for contusion and dislocation injuries, respectively. Furthermore, we present qualitative pilot data from a cadaveric trial, employing the MR rig, to show the expected results from future studies.

Nonlinear viscoelastic characterization of the porcine spinal cord

Although quasi-static and quasi-linear viscoelastic properties of the spinal cord have been reported previously, there are no published studies that have investigated the fully (strain-dependent) nonlinear viscoelastic properties of the spinal cord. In this study, stress relaxation experiments and dynamic cycling were performed on six fresh porcine lumbar cord specimens to examine their viscoelastic mechanical properties. The stress relaxation data were fitted to a modified superposition formulation and a novel finite ramp time correction technique was applied. The parameters obtained from this fitting methodology were used to predict the average dynamic cyclic viscoelastic behavior of the porcine cord. The data indicate that the porcine spinal cord exhibited fully nonlinear viscoelastic behavior. The average weighted root mean squared error for a Heaviside ramp fit was 2.8 kPa, which was significantly greater (p<0.001) than that of the nonlinear (comprehensive viscoelastic characterization method) fit (0.365 kPa). Further, the nonlinear mechanical parameters obtained were able to accurately predict the dynamic behavior, thus exemplifying the reliability of the obtained nonlinear parameters. These parameters will be important for future studies investigating various damage mechanisms of the spinal cord and studies developing high-resolution finite elements models of the spine.

Finite Element Modeling of CNS White Matter Kinematics: Use of a 3D RVE to Determine Material Properties

Axonal injury represents a critical target area for the prevention and treatment of traumatic brain and spinal cord injuries. Finite element (FE) models of the head and/or brain are often used to predict brain injury caused by external mechanical loadings, such as explosive waves and direct impact. The accuracy of these numerical models depends on correctly determining the material properties and on the precise depiction of the tissues' microstructure (microscopic level). Moreover, since the axonal microstructure for specific regions of the brain white matter is locally oriented, the stress, and strain fields are highly anisotropic and axon orientation dependent. Additionally, mechanical strain has been identified as the proximal cause of axonal injury, which further demonstrates the importance of this multi-scale relationship. In this study, our previously developed FE and kinematic axonal models are coupled and applied to a pseudo 3-dimensional representative volume element of central nervous system white matter to investigate the multi-scale mechanical behavior. An inverse FE procedure was developed to identify material parameters of spinal cord white matter by combining the results of uniaxial testing with FE modeling. A satisfactory balance between simulation and experiment was achieved via optimization by minimizing the squared error between the simulated and experimental force-stretch curve. The combination of experimental testing and FE analysis provides a useful analysis tool for soft biological tissues in general, and specifically enables evaluations of the axonal response to tissue-level loading and subsequent predictions of axonal damage.

The effects of substrate elastic modulus on neural precursor cell behavior

The spinal cord has a limited capacity to self-repair. After injury, endogenous stem cells are activated and migrate, proliferate, and differentiate into glial cells. The absence of neuronal differentiation has been partly attributed to the interaction between the injured microenvironment and neural stem cells. In order to improve post-injury neuronal differentiation and/or maturation potential, cell-cell and cell-biochemical interactions have been investigated. However, little is known about the role of stem cell-matrix interactions on stem cell-mediated repair. Here, we specifically examined the effects of matrix elasticity on stem cell-mediated repair in the spinal cord, since spinal cord injury results in drastic changes in parenchyma elasticity and viscosity. Spinal cord-derived neural precursor cells (NPCs) were grown on bis-acrylamide substrates with various rigidities. NPC growth, proliferation, and differentiation were examined and optimal in the range of normal spinal cord elasticity. In conclusion, limitations in NPC growth, proliferation, and neuronal differentiation were encountered when substrate elasticity was not within normal spinal cord tissue elasticity ranges. These studies elucidate the effect injury mediated mechanical changes may have on tissue repair by stem cells. Furthermore, this information can be applied to the development of future neuroregenerative biomaterials for spinal cord repair.

Gross morphological changes of the spinal cord immediately after surgical decompression in a large animal model of traumatic spinal cord injury

STUDY DESIGN

Quantitative in vivo ultrasound imaging study of spinal cord and dura morphology after acute experimental spinal cord injury (SCI) and decompression in a pig model.

OBJECTIVE

To study the morphological changes of the spinal cord and dura immediately after surgical decompression for acute SCI.

SUMMARY OF BACKGROUND DATA

Surgical decompression for traumatic SCI is currently a topic of debate. After decompression, relief of bony impingement on the thecal sac and spinal cord can be confirmed intraoperatively. However, postoperative imaging often reveals that the cord has swollen to fill the subarachnoid space. Little is known about the extent and timing of this morphological response.

METHODS

Yucatan miniature pigs received sham surgery (N = 1) or a moderate (N = 6, 20 g, 2.3 m/s) or high (N = 6, 20 g, 4.7 m/s) severity weight-drop SCI followed by 8 hours of sustained compression (100 g) and 6 hours of postdecompression monitoring. Sagittal-plane ultrasound images were used to quantify spinal cord, dura, and subarachnoid space dimensions preinjury and once per hour after decompression.

RESULTS

Animals with a moderate SCI exhibited a residual cord deformation of up to 0.64 mm within 10 minutes of decompression, which tended to resolve during 6 hours because of tissue relaxation and swelling. For animals with high-severity SCIs, cord swelling was immediate and resulted in occlusion of the subarachnoid space within 10 minutes to 5 hours, whereas this occurred for only half of the moderate injury group.

CONCLUSION

Decompression of an acute SCI may result in residual cord deformation followed by gradual swelling or immediate swelling leading to subarachnoid occlusion. The response is dependent on initial injury severity. These observations may partly explain the lack of benefit of decompression in some patients and suggest a need to reduce cord swelling to optimize the clinical outcome after acute SCI.

Mechanical characterization of the injured spinal cord after lateral spinal hemisection injury in the rat

The glial scar formed at the site of traumatic spinal cord injury (SCI) has been classically hypothesized to be a potent physical and biochemical barrier to nerve regeneration. One longstanding hypothesis is that the scar acts as a physical barrier due to its increased stiffness in comparison to uninjured spinal cord tissue. However, the information regarding the mechanical properties of the glial scar in the current literature is mostly anecdotal and not well quantified. We monitored the mechanical relaxation behavior of injured rat spinal cord tissue at the site of mid-thoracic spinal hemisection 2 weeks and 8 weeks post-injury using a microindentation test method. Elastic moduli were calculated and a modified standard linear model (mSLM) was fit to the data to estimate the relaxation time constant and viscosity. The SLM was modified to account for a spectrum of relaxation times, a phenomenon common to biological tissues, by incorporating a stretched exponential term. Injured tissue exhibited significantly lower stiffness and elastic modulus in comparison to uninjured control tissue, and the results from the model parameters indicated that the relaxation time constant and viscosity of injured tissue were significantly higher than controls. This study presents direct micromechanical measurements of injured spinal cord tissue post-injury. The results of this study show that the injured spinal tissue displays complex viscoelastic behavior, likely indicating changes in tissue permeability and diffusivity.

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.

Spinal subarachnoid space pressure measurements in an in vitro spinal stenosis model: implications on syringomyelia theories

Full explanation for the pathogenesis of syringomyelia (SM), a neuropathology characterized by the formation of a cystic cavity (syrinx) in the spinal cord (SC), has not yet been provided. It has been hypothesized that abnormal cerebrospinal fluid (CSF) pressure, caused by subarachnoid space (SAS) flow blockage (stenosis), is an underlying cause of syrinx formation and subsequent pain in the patient. However, paucity in detailed in vivo pressure data has made theoretical explanations for the syrinx difficult to reconcile. In order to understand the complex pressure environment, four simplified in vitro models were constructed to have anatomical similarities with post-traumatic SM and Chiari malformation related SM. Experimental geometry and properties were based on in vivo data and incorporated pertinent elements such as a realistic CSF flow waveform, spinal stenosis, syrinx, flexible SC, and flexible spinal column. The presence of a spinal stenosis in the SAS caused peak-to-peak cerebrospinal fluid CSF pressure fluctuations to increase rostral to the stenosis. Pressure with both stenosis and syrinx present was complex. Overall, the interaction of the syrinx and stenosis resulted in a diastolic valve mechanism and rostral tensioning of the SC. In all experiments, the blockage was shown to increase and dissociate SAS pressure, while the axial pressure distribution in the syrinx remained uniform. These results highlight the importance of the properties of the SC and spinal SAS, such as compliance and permeability, and provide data for comparison with computational models. Further research examining the influence of stenosis size and location, and the importance of tissue properties, is warranted.

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.