An in-vivo measurement and analysis of viscoelastic properties of the spinal cord of cats

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.

Simulation of cancer progression in bone by a virtual thermal flux with a case study on lumbar vertebrae with multiple myeloma

BACKGROUND

Two main problems examining the mechanism of cancer progression in the tissues using the computational models are lack of enough knowledge on the effective factors for such events in vivo environments and lack of specific parameters in the available computational models to simulate such complicated reactions.

METHODS

In this study, it was tried to simulate the progression of cancerous lesions in the bone tissues by an independent parameter from the anatomical and physiological characteristics of the tissues, so to degrade the orthotropic mechanical properties of the bone tissues, a virtual temperature was determined to be used by a well-known framework for simulation of damages in the composite materials. First, the reliability of the FE model to simulate hyperelastic response in the intervertebral discs (IVDs) and progressive failure in the bony components were verified by simulation of some In-Vitro tests, available in the literature. Then, the progression of the osteolytic damage was simulated in a clinical case with multiple myeloma in the lumbar vertebrae.

RESULTS

The FE model could simulate stress-shielding and diffusion of lesion in the posterior elements of the damaged vertebra which led to spinal stenosis. The load carrying shares associated with the anterior half and the posterior half of the examined vertebral body and the posterior elements were estimated equal to 41 %, 47 % and 12 %, respectively for the intact condition, that changed to 14 %, 16 % and 70 %, when lesion occupied one third of the vertebral body.

CONCLUSION

Correlation of the FE results with the deformation shapes, observed in the MRIs for the clinical case study, indicated appropriateness of the procedure, proposed for simulation of the progressive osteolytic damage in the vertebral segments. The future studies may follow simulation of tumor growth for various metastatic tissues using the method, established here.

Stiffness-matched biomaterial implants for cell delivery: clinical, intraoperative ultrasound elastography provides a 'target' stiffness for hydrogel synthesis in spinal cord injury

Safe hydrogel delivery requires stiffness-matching with host tissues to avoid iatrogenic damage and reduce inflammatory reactions. Hydrogel-encapsulated cell delivery is a promising combinatorial approach to spinal cord injury therapy, but a lack of in vivo clinical spinal cord injury stiffness measurements is a barrier to their use in clinics. We demonstrate that ultrasound elastography – a non-invasive, clinically established tool – can be used to measure spinal cord stiffness intraoperatively in canines with spontaneous spinal cord injury. In line with recent experimental reports, our data show that injured spinal cord has lower stiffness than uninjured cord. We show that the stiffness of hydrogels encapsulating a clinically relevant transplant population (olfactory ensheathing cells) can also be measured by ultrasound elastography, enabling synthesis of hydrogels with comparable stiffness to canine spinal cord injury. We therefore demonstrate proof-of-principle of a novel approach to stiffness-matching hydrogel-olfactory ensheathing cell implants to ‘real-life’ spinal cord injury values; an approach applicable to multiple biomaterial implants for regenerative therapies.

Rapid prototyping fabrication of soft and oriented polyester scaffolds for axonal guidance

Biodegradable polyesters have been extensively used for preparation of nerve guidance scaffolds, due to their high biocompatibility and defined degradation periods. However, conventional methods for fabrication of porous polyester scaffolds provide limited control over shape and micro-architecture. Here, a fabrication procedure based on 3D printing was developed to generate highly ordered and anatomically personalized, polyester scaffolds for soft tissue regeneration. Scaffolds composed of Poly-lactic-glycolic acid (PLGA) and poly-L-lactic acid (PLLA) were specifically customized for nerve injuries. This was obtained by using an oriented multi-layer printing pattern which established a linear structure in the fabricated scaffolds to match the aligned topography of nerve tissues. The oriented scaffold was shown to guide regenerating axons to linear conformations and support growth of induced pluripotent stem cell-derived neurons in vitro and in vivo in a model of spinal cord injury. The described scaffolds may advance the field of nerve regeneration. Furthermore, modifications could be integrated to generate soft implants for various types of tissues.

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.

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.

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.

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.

Advances in Nerve Guidance Conduit-Based Therapeutics for Peripheral Nerve Repair

Peripheral nerve injuries have high incidence rates, limited treatment options and poor clinical outcomes, rendering a significant socioeconomic burden. For effective peripheral nerve repair, the gap or site of injury must be structurally bridged to promote correct reinnervation and functional regeneration. However, effective repair becomes progressively more difficult with larger gaps. Autologous nerve grafting remains the best clinical option for the repair of large gaps (20-80 mm) despite being associated with numerous limitations including permanent donor site morbidity, a lack of available tissue and the formation of neuromas. To meet the clinical demand of large gap repair and overcome these limitations, tissue engineering has led to the development of nerve guidance conduit-based therapeutics. This review focuses on the advances of nerve guidance conduit-based therapeutics in terms of their structural properties including biomimetic composition, permeability, architecture, and surface modifications. Associated biochemical properties, pertaining to the incorporation of cells and neurotrophic factors, are also reviewed. After reviewing the progress in the field, we conclude by presenting an outlook on their clinical translatability and the next generation of therapeutics.

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.

Localized and sustained release of brain-derived neurotrophic factor from injectable hydrogel/microparticle composites fosters spinal learning after spinal cord injury

Injectable hydrogel allows for sustained delivery of growth factor resulting in spinal mediated learning after injury.

Laser 3D printing with sub-microscale resolution of porous elastomeric scaffolds for supporting human bone stem cells

A reproducible method is needed to fabricate 3D scaffold constructs that results in periodic and uniform structures with precise control at sub-micrometer and micrometer length scales. In this study, fabrication of scaffolds by two-photon polymerization (2PP) of a biodegradable urethane and acrylate-based photoelastomer is demonstrated. This material supports 2PP processing with sub-micrometer spatial resolution. The high photoreactivity of the biophotoelastomer permits 2PP processing at a scanning speed of 1000 mm s(-1), facilitating rapid fabrication of relatively large structures (>5 mm(3)). These structures are custom printed for in vitro assay screening in 96-well plates and are sufficiently flexible to enable facile handling and transplantation. These results indicate that stable scaffolds with porosities of greater than 60% can be produced using 2PP. Human bone marrow stromal cells grown on 3D scaffolds exhibit increased growth and proliferation compared to smooth 2D scaffold controls. 3D scaffolds adsorb larger amounts of protein than smooth 2D scaffolds due to their larger surface area; the scaffolds also allow cells to attach in multiple planes and to completely infiltrate the porous scaffolds. The flexible photoelastomer material is biocompatible in vitro and is associated with facile handling, making it a viable candidate for further study of complex 3D-printed scaffolds.

3D in vitro modeling of the central nervous system

There are currently more than 600 diseases characterized as affecting the central nervous system (CNS) which inflict neural damage. Unfortunately, few of these conditions have effective treatments available. Although significant efforts have been put into developing new therapeutics, drugs which were promising in the developmental phase have high attrition rates in late stage clinical trials. These failures could be circumvented if current 2D in vitro and in vivo models were improved. 3D, tissue-engineered in vitro systems can address this need and enhance clinical translation through two approaches: (1) bottom-up, and (2) top-down (developmental/regenerative) strategies to reproduce the structure and function of human tissues. Critical challenges remain including biomaterials capable of matching the mechanical properties and extracellular matrix (ECM) composition of neural tissues, compartmentalized scaffolds that support heterogeneous tissue architectures reflective of brain organization and structure, and robust functional assays for in vitro tissue validation. The unique design parameters defined by the complex physiology of the CNS for construction and validation of 3D in vitro neural systems are reviewed here.

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.

Mitigating spinal cord distraction injuries: the effect of durotomy in decreasing cord interstitial pressure in vitro

STUDY DESIGN

The present study involved an in vitro examination of spinal cord interstitial pressure (CIP) during distraction before and after durotomy in three spinal cord segments obtained from five pigs.

OBJECTIVES

To determine whether durotomy can be used to decrease the elevated CIP associated with spinal cord distraction.

SUMMARY OF BACKGROUND DATA

Spinal cord distraction is a known cause of spinal cord injury. Several articles describing the pathophysiology of cord distraction injuries suggest that the underlying mechanism of injury is a microvascular ischemic event. The authors have previously described an increase in CIP with spinal cord distraction, with average pressures of 23 mmHg at loads of 1,000 g. To date, there are no published studies that have evaluated the efficacy of intentional durotomies as a treatment for elevated CIP.

METHODS

A total of 15 spinal cord sections were harvested from pigs and distracted while immersed in saline, using a fixed 1,000 g distraction force. The CIP decay was then measured at 30-s intervals for 10 min. The distraction/relaxation maneuver was performed six times with continuous CIP monitoring and was subsequently followed by durotomy.

RESULTS

The pressure-decay curves were similar for each specimen, but varied according to individual pigs and anatomical levels. CIP decayed over the first 4 min of distraction and remained constant for the final 6 min. Longitudinal durotomy led to a dramatic drop in CIP toward baseline and appeared to be as effective as transverse durotomy with regard to the normalization of pressure.

CONCLUSION

Spinal cord distraction causes elevations in CIP. Durotomy lowers elevated CIP in vitro and may be effective at lowering CIP in vivo. Further study is required to evaluate the usefulness of durotomy in vivo.

Mechanical properties of the lamprey spinal cord: uniaxial loading and physiological strain

During spinal cord injury, nerves suffer a strain beyond their physiological limits which damages and disrupts their structure. Research has been done to measure the modulus of the spinal cord and surrounding tissue; however the relationship between strain and spinal cord fibers is still unclear. In this work, our objective is to measure the stress-strain response of the spinal cord in vivo and in vitro and model this response as a function of the number of fibers. We used the larvae lamprey (Petromyzon Marinus), a model for spinal cord regeneration and animal locomotion. We found that physiologically the spinal cord is pre-stressed to a longitudinal strain of 10% and this strain increases to 15% during swimming. Tensile measurements show that uniaxial, longitudinal loading is independent of the meninges. Stress values for uniaxial strains below 18%, are homogeneous through the length of the body. However, for higher uniaxial strains the Head section shows more resistance to longitudinal loading than the Tail. These data, together with the number of fibers obtained from histological sections were used in a composite-material model to obtain the properties of the spinal cord fibers (2.4 MPa) and matrix (0.017 MPa) to uniaxial longitudinal loading. This model allowed us to approximate the percentage of fibers in the spinal cord, establishing a relationship between uniaxial longitudinal strains and spinal cord composition. We showed that there is a proportional relationship between the number of fibers and the properties of the spinal cord at large uniaxial strains.

Characterization of hyaluronan-methylcellulose hydrogels for cell delivery to the injured spinal cord

No effective clinical treatment currently exists for traumatic spinal cord injury. Cell replacement therapy holds promise for attaining functional repair. Cells may be delivered directly or near the injury site; however, this strategy requires a delivery vehicle to maintain cell viability. We have identified an injectable, biocompatible, and biodegradable hydrogel scaffold composed of hyaluronan (HA) and methylcellulose (MC) that may be an effective scaffold for therapeutic cell delivery. The purpose of the present study was to determine the effects of polymer concentration on HAMC mechanical strength, gelation time, and cell viability. The yield stress of HAMC, a measure of mechanical stiffness, was tunable via manipulation of MC and HA content. Measurement of the elastic and storage moduli as functions of time revealed that HAMC gels in less than 5 min at physiological temperatures. Human umbilical tissue-derived cells encapsulated in HAMC were homogenously and stably distributed over 3 days in culture and extended processes into the scaffold. Cell viability was stable over this period in all but the most concentrated HAMC formulation. Because of its strength-tunability, rapid gelation, and ability to maintain cell viability, HAMC is a promising vehicle for cell delivery and is being tested in ongoing in vivo studies.

Viscoelastic properties of bovine orbital connective tissue and fat: constitutive models

Reported mechanical properties of orbital connective tissue and fat have been too sparse to model strain-stress relationships underlying biomechanical interactions in strabismus. We performed rheological tests to develop a multi-mode upper convected Maxwell (UCM) model of these tissues under shear loading. From 20 fresh bovine orbits, 30 samples of connective tissue were taken from rectus pulley regions and 30 samples of fatty tissues from the posterior orbit. Additional samples were defatted to determine connective tissue weight proportion, which was verified histologically. Mechanical testing in shear employed a triborheometer to perform: strain sweeps at 0.5-2.0 Hz; shear stress relaxation with 1% strain; viscometry at 0.01-0.5 s(-1) strain rate; and shear oscillation at 1% strain. Average connective tissue weight proportion was 98% for predominantly connective tissue and 76% for fatty tissue. Connective tissue specimens reached a long-term relaxation modulus of 668 Pa after 1,500 s, while corresponding values for fatty tissue specimens were 290 Pa and 1,100 s. Shear stress magnitude for connective tissue exceeded that of fatty tissue by five-fold. Based on these data, we developed a multi-mode UCM model with variable viscosities and time constants, and a damped hyperelastic response that accurately described measured properties of both connective and fatty tissues. Model parameters differed significantly between the two tissues. Viscoelastic properties of predominantly connective orbital tissues under shear loading differ markedly from properties of orbital fat, but both are accurately reflected using UCM models. These viscoelastic models will facilitate realistic global modeling of EOM behavior in binocular alignment and strabismus.

Effects of white, grey, and pia mater properties on tissue level stresses and strains in the compressed spinal cord

Recent demographics demonstrate an increase in the number of elderly spinal cord injury patients, motivating the desire for a better understanding of age effects on injury susceptibility. Knowing that age and disease affect neurological tissue, there is a need to better understand the sensitivity of spinal cord injury mechanics to variations in tissue behavior. To address this issue, a plane-strain, geometrically nonlinear, finite element model of a section of a generic human thoracic spinal cord was constructed to model the response to dorsal compression. The material models and stiffness responses for the grey and white matter and pia mater were varied across a range of reported values to observe the sensitivity of model outcomes to the assigned properties. Outcome measures were evaluated for percent change in magnitude and alterations in spatial distribution. In general, principal stresses (114-244% change) and pressure (75-119% change) were the outcomes most sensitive to material variation. Strain outcome measures were less sensitive (7-27% change) than stresses (74-244% change) to variations in material tangent modulus. The pia mater characteristics had limited (<4% change) effects on outcomes. Using linear elastic models to represent non-linear behavior had variable effects on outcome measures, and resulted in highly concentrated areas of elevated stresses and strains. Pressure measurements in both the grey and white matter were particularly sensitive to white matter properties, suggesting that degenerative changes in white matter may influence perfusion in a compressed spinal cord. Our results suggest that the mechanics of spinal cord compression are likely to be affected by changes in tissue resulting from aging and disease, indicating a need to study the biomechanical aspects of spinal cord injury in these specific populations.

Three-dimensional finite element model of the cervical spinal cord: preliminary results of injury mechanism analysis

STUDY DESIGN

A three-dimensional finite element investigation.

OBJECTIVES

To create a three-dimensional finite element model of the cervical spinal cord enlargement and to simulate a hyperextension injury of the cervical cord.

SUMMARY OF BACKGROUND DATA

Experimental studies are difficult to simulate the complex mechanism of spinal cord injuries. The introduction of three-dimensional modeling technique into neurotrauma studies is essential to further understand mechanical behavior of the nerve tissue during traumatic injuries.

METHODS

Geometrical reconstruction of cervical spinal cord enlargement was developed based on the morphologic features of each segment of the fresh human cervical cord. After the validation of the model, the pinching condition in the hyperextension injuries was simulated with compressive and extension forces applied on the cervical enlargement model. The average von Mises stress of the 9 anatomic regions, such as anterior funiculus, lateral part of the lateral funiculus, medial part of the lateral funiculus, lateral part of the posterior funiculus, medial part of the posterior funiculus, anterior horn, the bottom of anterior horn, the apex of posterior horn, the cervix cornu posterioris, and caput cornu posterioris was recorded.

RESULTS

The force-displacement response of the spinal cord under compression and axial tension loading was close to the experimental results reported in the literature. The stress distribution of the spinal cord according to the numerical simulation and the morphologic features of the in vivo experiment were also in close agreement. Hyperextension injury simulation showed high localized stress at the anterior and posterior horn in the gray matter.

CONCLUSION

The finite element method as a three-dimensional modeling technique can improve the understanding of the biomechanical behavior of the spinal cord. The results of hyperextension injury simulation of the cervical spinal cord probably account for the predominance of the hand weakness in patients with central cord injury.

The mechanical properties of neonatal rat spinal cord in vitro, and comparisons with adult

A number of studies have investigated the mechanical properties of adult spinal cord under tension, however it is not known whether age has an effect on these properties. This is of interest to those aiming to understand the clinical differences between adults and children with spinal cord injury (e.g. severity and recovery), and those developing experimental or computational models for paediatric spinal cord injury. Entire spinal cords were freshly harvested from neonatal rats (14 days) and tested in vitro under uniaxial tension at a range of strain rates (0.2, 0.02, 0.002/s) to a range of strains (2%, 3.5%, 5%), with relaxation responses being recorded for 15 min. These mechanical properties were compared to previously reported data from similar experiments on adult rat spinal cords, and the peak stress and the stress after 15 min of relaxation were found to be significantly higher for spinal cords from adults than neonates (p<0.001). A non-linear viscoelastic model was developed and was observed to adequately predict the mechanical behaviour of this tissue. The model developed in this study may be of use in computational models of paediatric spinal cord. The significant differences between adult and neonatal spinal cord properties may explain the higher initial severity of spinal cord injury in children and may have implications for the development of experimental animal models for paediatric spinal cord injury, specifically for those aiming to match the injury severity with adult experimental models.

Quasilinear viscoelastic behavior of bovine extraocular muscle tissue

PURPOSE

Until now, there has been no comprehensive mathematical model of the nonlinear viscoelastic stress-strain behavior of extraocular muscles (EOMs). The present study describes, with the use of a quasilinear viscoelastic (QLV) model, the nonlinear, history-dependent viscoelastic properties and elastic stress-strain relationship of EOMs.

METHODS

Six oculorotary EOMs were obtained fresh from a local abattoir. Longitudinally oriented specimens were taken from different regions of the EOMs and subjected to uniaxial tensile, relaxation, and cyclic loading testing with the use of an automated load cell under temperature and humidity control. Twelve samples were subjected to uniaxial tensile loading with 1.7%/s strain rate until failure. Sixteen specimens were subjected to relaxation studies over 1500 seconds. Cyclic loading was performed to validate predictions of the QLV model characterized from uniaxial tensile loading and relaxation data.

RESULTS

Uniform and highly repeatable stress-strain behavior was observed for 12 specimens extracted from various regions of all EOMs. Results from 16 different relaxation trials illustrated that most stress relaxation occurred during the first 30 to 60 seconds for 30% extension. Elastic and reduced relaxation functions were fit to the data, from which a QLV model was assembled and compared with cyclic loading data. Predictions of the QLV model agreed with observed peak cyclic loading stress values to within 8% for all specimens and conditions.

CONCLUSIONS

Close agreement between the QLV model and the relaxation and cyclic loading data validates model quantification of EOM mechanical properties and will permit the development of accurate overall models of mechanics of ocular motility and strabismus.

The development of an improved physical surrogate model of the human spinal cord--tension and transverse compression

To prevent spinal cord injury, optimize treatments for it, and better understand spinal cord pathologies such as spondylotic myelopathy, the interaction between the spinal column and the spinal cord during injury and pathology must be understood. The spinal cord is a complex and very soft tissue that changes properties rapidly after death and is difficult to model. Our objective was to develop a physical surrogate spinal cord with material properties closely corresponding to the in vivo human spinal cord that would be suitable for studying spinal cord injury under a variety of injurious conditions. Appropriate target material properties were identified from published studies and several candidate surrogate materials were screened, under uniaxial tension, in a materials testing machine. QM Skin 30, a silicone elastomer, was identified as the most appropriate material. Spinal cords manufactured from QM Skin 30 were tested under uniaxial tension and transverse compression. Rectangular specimens of QM Skin 30 were also tested under uniform compression. QM Skin 30 produced surrogate cords with a Young's modulus in tension and compression approximately matching values reported for in vivo animal spinal cords (0.25 and 0.20 MPa, respectively). The tensile and compressive Young's modulus and the behavior of the surrogate cord simulated the nonlinear behavior of the in vivo spinal cord.

An in vitro study on the effects of freezing, spine segment, repeat measurement, and individual cord properties on cord interstitial pressure

STUDY DESIGN

In vitro study of the spinal cord tension and pressure relationships before and after thawing in 6 different spinal cord segment from 2 individual pigs.

OBJECTIVES

To determine if frozen and thawed spinal cord segments had different tension/cord interstitial pressure(CIP) relationships to fresh spinal cord segments. In addition, we will determine if the cord level, individual cord properties, and repeated CIP measurements affect the tension/CIP relationships.

SUMMARY OF BACKGROUND DATA

Spinal cord distraction is a known cause of spinal cord injury. Several articles published on the pathophysiology of the cord distraction injury suggest that the underlying mechanism of injury is a microvascular ischemic event. We have previously described an increase in CIP with spinal cord distraction, pressures average 23 mmHg at 1 kg loads.

METHODS

Six cord segments harvested from 2 pigs contained cervical, thoracic, and lumbar segments, and underwent distraction using a series of 7 calibrated weights from 0 to 1000 g weight. The cords were measured at each level of distraction. The cords were then frozen at -20 degrees C for a period of 2 weeks, and then thawed and retested. Multiple linear regression was then performed.

RESULTS

There was no difference between the fresh and the frozen-thawed cords; there was statistical difference between the 2 pigs (18 mmHg) (P < 0.001). There are differences between the cervical and the thoracic cord segments (P < 0.001), and between cervical and lumbar cord segments (P = 0.056). There is a significant relation between the tension applied and CIP. Repeated trials showed no drift with repeated measures.

CONCLUSION

Freezing and thawing spinal cords has no effect on the CIP/tension curves. Cord interstitial pressure developed is dependant on cord tension, cord level, individual cord properties, but not on the number of repetitions carried out while testing the spinal cord.

The effect of cerebrospinal fluid on the biomechanics of spinal cord: an ex vivo bovine model using bovine and physical surrogate spinal cord

STUDY DESIGN

A biomechanical study using ex vivo bovine spinal cord and dura, and a synthetic surrogate spinal cord with bovine dura.

OBJECTIVE

To investigate the effect of cerebrospinal fluid (CSF) on spinal cord deformation characteristics and to evaluate the biofidelity of a new surrogate spinal cord using an ex vivo bovine model of the burst fracture process.

SUMMARY OF BACKGROUND DATA

Spinal cord injury is associated with significant personal, economic and social costs. The role of CSF during the injury event and its effect on the spinal cord deformation and neurologic injury is not well understood. Such knowledge could inform preventative strategies and clinical interventions and aid the development and validation of experimental and computational models.

METHODS

The transverse impact of a propelled bone fragment analogue with bovine and surrogate cord models was recorded with high speed video and the images analyzed to determine deformation trajectories. Each cord specimen was tested in 3 states: with dura and CSF, with dura only, and without dura. The effect of these states on deformation magnitude, duration, and energy loss parameters was assessed. RESULTS.: The estimated spinal cord deformation was significantly reduced, although not eliminated, in the presence of CSF when compared to the bare state. The duration of deformation was generally increased in the presence of CSF, though this difference was not statistically significant. This may indicate a reduction in the cord-fragment interaction force for a given impulse. The dura was found to have no significant effect on deformation parameters for the bovine spinal cord. The deformation of the surrogate cord gave similar trends for the different states in comparison to the bovine cord, but was significantly less than the bovine spinal cord for all conditions.

CONCLUSION

The results indicate that the protective mechanism of CSF may not eliminate cord deformationunder the high energy transverse impact characteristic of a burst fracture. However, CSF may contribute to a lessening of cord deformation and applied force.

Rheological properties of the tissues of the central nervous system: a review

Knowledge of the biomechanical properties of central nervous system (CNS) tissues is important for understanding mechanisms and thresholds for injury, and aiding development of computer or surrogate models of these tissues. Many investigations have been conducted to estimate the properties of CNS tissues including under shear, compressive and tensile loading, however there is much variability in this body of literature, making it difficult to separate the material properties from effects that result from a given experimental protocol. This review summarises previous studies of brain and spinal cord properties; discussing their main findings and points of difference, and displays the reported data on comparable scales. Additionally, based on the observed effects of methodological choices on reported tissue properties, recommendations for future studies of brain and spinal cord properties are made.

A three-dimensional finite element model of the cervical spine with spinal cord: an investigation of three injury mechanisms

The spinal cord may be injured through various spinal column injury patterns (e.g., burst fracture, fracture dislocation); however, the relationship between column injury pattern and cord damage is not well understood. A three-dimensional finite element model of a human cervical spine and spinal cord segment was developed, verified using published experimental data, and used to investigate differences in cord strain distributions during various column injury patterns. For a transverse contusion injury, as would occur in a burst fracture, a 33% canal occlusion resulted in two peaks of strain between the indentor and opposing vertebral body and intermediate peak strain values. For a distraction injury, relevant to column distortion injuries, a 2.6 mm axial displacement to the cord resulted in more uniform strains throughout the cord and low peak strain values. For a dislocation injury, as would occur in a fracture dislocation, an anterior displacement of C5 corresponding to 30% of the sagittal dimension of the vertebral body resulted in high peak strain values adjacent to the shearing vertebrae and increased strains in the lateral columns compared to contusion. This model includes more anatomical details compared to previous studies and provides a baseline for mechanical comparisons in spinal cord injury.

The biomechanical response of spinal cord tissue to uniaxial loading

The spinal cord is an integral component of the spinal column and is prone to physical injury during trauma or more long-term pathological insults. The development of computational models to simulate the cord-column interaction during trauma is important in developing a proper understanding of the injury mechanism. Such models would be invaluable in seeking both preventive strategies that reduce the propensity for injury and identifying specific treatment regimes. However, these developments are hampered by the limited information available on the structural and mechanical properties of this soft tissue owing to the difficulty in handling this material in a cadaveric situation. The purpose of the present paper is to report the rapid deterioration in the quality of the tissues once excised, which provides a further challenge to the successful elucidation of the structural properties of the tissue. In particular, the tangent modulus of the tissue is seen to increase sharply over a period of 72 h.

The mechanical properties of rat spinal cord in vitro

Freshly excised rat spinal cords were tested in uniaxial tension, in vitro, at strain rates ranging from 0.002 to 0.2 s-1. Stress relaxation tests were performed for a range of strains from 2% to 5%, with the relaxation behaviour being recorded for a period of at least 30 min. Samples exhibited a characteristic "J" shaped non-linear stress-strain response, with stiffness increasing with applied strain. The cords were labelled with rows of small markers and the uniaxial tension tests were recorded via video. Subsequent image analysis enabled the distribution of strain on the cord surface to be determined. Viscoelastic models were developed to model the mechanical behaviour of the specimens and were found to adequately describe the material behaviour.

Guided cell adhesion and outgrowth in peptide-modified channels for neural tissue engineering

A hydrogel scaffold of well-defined geometry was created and modified with laminin-derived peptides in an aqueous solution, thereby maintaining the geometry of the scaffold while introducing bioactive peptides that enhance cell adhesion and neurite outgrowth. By combining a fiber templating technique to create longitudinal channels with peptide modification, we were able to synthesize a scaffold that guided cell adhesion and neurite outgrowth of primary neurons. Scaffolds were designed to have numerous longitudinally oriented channels with an average channel diameter of 196 +/- 6 microm to ultimately promote fasciculation of regenerating cables and a compressive modulus of 192 +/- 8 kPa to match the modulus of the soft nerve tissue. Copolymerization of 2-hydroxylethyl methacrylate (HEMA) with 2-aminoethyl methacrylate (AEMA) scaffolds, provided primary amine groups to which two sulfhydryl terminated, laminin-derived oligopeptides, CDPGYIGSR and CQAASIKVAV, were covalently bound using the sulfo-(N-maleimidomethyl)cyclohexane-1-carboxylate (sulfo-SMCC) crosslinking agent. The concentration of peptides on the scaffolds was measured at 106 +/- 4 micromol/cm(2) using the ninhydrin method and 92 +/- 9 micromol/cm(2) using the BCA protein assay. The peptide modified P(HEMA-co-AEMA) scaffolds were easily fabricated in aqueous conditions, highly reproducible, well-defined, and enhanced neural cell adhesion and guided neurite outgrowth of primary chick dorsal root ganglia neurons relative to non-peptide-modified controls. The copolymerization of AEMA with HEMA can be extended to other radically polymerized monomers and is advantageous as it facilitates scaffold modification in aqueous solutions thereby obviating the use of organic solvents which can be cytotoxic and often disrupt scaffold geometry. The combination of well-defined chemical and physical stimuli described herein provides a means for guided regeneration both in vitro and in vivo.

Synthetic hydrogel guidance channels facilitate regeneration of adult rat brainstem motor axons after complete spinal cord transection

Synthetic guidance channels or tubes have been shown to promote axonal regeneration within the spinal cord from brainstem motor nuclei with the inclusion of agents such as matrices, cells, or growth factors to the tube. We examined the biocompatibility and regenerative capacity of synthetic hydrogel tubular devices that were composed of poly(2-hydroxyethyl methacrylate-co-methyl methacrylate) (PHEMA-MMA). Two PHEMA-MMA channels, having a mean elastic modulus of either 177 or 311 kPa were implanted into T8-transected spinal cords of adult Sprague Dawley rats. The cord stumps were inserted into the channels and fibrin glue was applied to the cord-channel interface. An expanded polytetrafluoroethylene (ePTFE) membrane was used for duraplasty. Controls underwent cord transection alone. Gross and microscopic examination of the spinal cords showed continuity of tissue within the synthetic guidance channels between the cord stumps at 4 and 8 weeks. There was a trend towards an increased area and width of bridging neural tissue in the 311-kPa guidance channels compared to the 177-kPa channels. Neurofilament stained axons were visualized within the bridging tissue, and serotonergic axons were found to enter the 311-kPa channel. Retrograde axonal tracing revealed regeneration of axons from reticular, vestibular, and raphe brainstem motor nuclei. For both channels, there was minimal scarring at the channel-cord interface, and less scarring at the channel-dura interface compared to that observed next to the ePTFE. The present study is the first to show that axons from brainstem motor nuclei regenerated in unfilled synthetic hydrogel guidance channels after complete spinal cord transection.

Manufacture of poly(2-hydroxyethyl methacrylate-co-methyl methacrylate) hydrogel tubes for use as nerve guidance channels

Hydrogel tubes of poly(2-hydroxyethyl methacrylate-co-methyl methacrylate) (p(HEMA-co-MMA)) made by liquid-liquid centrifugal casting are being investigated as potential nerve guidance channels in the central nervous system. An important criterion for the nerve guidance channel is that its mechanical properties are similar to those of the spinal cord, where it will be implanted. The formulated p(HEMA-co-MMA) tubes are soft and flexible, consisting of a gel-like outer layer, and an interconnected macroporous, inner layer. The relative thickness of the gel phase to macroporous phase is controlled by the formulation chemistry, and specifically by the ratio of co-monomers, HEMA and MMA. By varying the surface chemistry of the mold within which the tubes are synthesized, tubes were prepared with either a "cracked" or a smooth outer morphology. Tubes with the cracked outer morphology had periodic channels that traversed the wall of the tube, which resulted in a lower modulus than smooth outer morphology tubes, yet likely greater diffusive permeability. For tubes (and not rods) to be formed, phase separation must precede gelation as is detailed in a formulation phase diagram for HEMA, MMA and water. The tensile elastic modulus of p(HEMA-co-MMA) tubes reflected the formulation chemistry, with greater moduli (up to 400 kPa) recorded for tubes having 10 wt% MMA. The p(HEMA-co-MMA) tubes therefore had similar mechanical properties to those of the spinal cord, which has a reported elastic modulus range between 200 and 600 kPa.

A single integral finite strain viscoelastic model of ligaments and tendons

A general continuum model for the nonlinear viscoelastic behavior of soft biological tissues was formulated. This single integral finite strain (SIFS) model describes finite deformation of a nonlinearly viscoelastic material within the context of a three-dimensional model. The specific form describing uniaxial extension was obtained, and the idea of conversion from one material to another (at a microscopic level) was then introduced to model the nonlinear behavior of ligaments and tendons. Conversion allowed different constitutive equations to be used for describing a single ligament or tendon at different strain levels. The model was applied to data from uniaxial extension of younger and older human patellar tendons and canine medial collateral ligaments. Model parameters were determined from curve-fitting stress-strain and stress-relaxation data and used to predict the time-dependent stress generated by cyclic extensions.

The mechanical properties of the human cervical spinal cord in vitro

The response of spinal cord tissue to mechanical loadings is not well understood. In this study, isolated fresh cervical spinal cord samples were obtained from cadavers at autopsy and tested in uniaxial tension at moderate strain rates. Stress relaxation experiments were performed with an applied strain rate and peak strain in the physiological range, similar to those seen in the spinal cord during voluntary motion. The spinal cord samples exhibited a nonlinear stress-strain response with increasing strain increasing the tangent modulus. In addition, significant relaxation was observed over 1 min. A quasilinear viscoelastic model was developed to describe the behavior of the spinal cord tissue and was found to describe the material behavior adequately. The data also were fitted to both hyperelastic and viscoelastic fluid models for comparison with other data in the literature.