The internal mechanics of the intervertebral disc under cyclic loading

A novel in-vitro model of intervertebral disc degeneration using hyperphysiological loading

Intervertebral disc (IVD) degeneration includes changes in tissue biomechanics, physical attributes, biochemical composition, disc microstructure, and cellularity, which can all affect the normal function of the IVD, and ultimately may lead to pain. The purpose of this research was to develop an in-vitro model of degeneration that includes the evaluation of physical, biomechanical, and structural parameters, and that does so over several load/recovery periods. Hyperphysiological loading was used as the degenerative initiator with three experimental groups employed using bovine coccygeal IVD specimens: Control; Single-Overload; and Double-Overload. An equilibrium stage comprising a static load followed by two load/recovery periods was followed by six further load/recovery periods. In the Control group all load/recovery periods were the same, comprising physiological cyclic loading. The overload groups differed in that hyperphysiological loading was applied during the 4th loading period (Single-Overload), or the 4th and 5th loading period (Double-Overload). Overloading led to a significant reduction in disc height compared to the Control group, which was not recovered in subsequent physiological load/recovery periods. However, there were no significant changes in stiffness. Overloading also led to significantly more microstructural damage compared to the Control group. Taking all outcome measures into account, the overload groups were evaluated as replicating clinically relevant aspects of moderate IVD degeneration. Further research into a potential dose-effect, and how more severe degeneration can be replicated would provide a model with the potential to evaluate new treatments and interventions for different stages of IVD degeneration.

Regulating Movement Frequency and Speed: Implications for Lumbar Spine Load Management Strategies Demonstrated Using an In Vitro Porcine Model

The relationship between internal loading dose and low-back injury risk during lifting is well known. However, the implications of movement parameters that influence joint loading rates-movement frequency and speed-on time-dependent spine loading responses remain less documented. This study quantified the effect of loading rate and frequency on the tolerated cumulative loading dose and its relation to joint lifespan. Thirty-two porcine spinal units were exposed to biofidelic compression loading paradigms that differed by joint compression rate (4.2 and 8.3 kN/s) and frequency (30 and 60 cycles per minute). Cyclic compression testing was applied until failure was detected or 10,800 continuous cycles were tolerated. Instantaneous weighting factors were calculated to evaluate the cumulative load and Kaplan-Meier survival probability functions were examined following nonlinear dose normalization of the cyclic lifespan. Significant reductions in cumulative compression were tolerated when spinal units were compressed at 8.3 kN/s (P < .001, 67%) and when loaded at 30 cycles per minute (P = .008, 45%). There was a positive moderate relationship between cumulative load tolerance and normalized cyclic lifespan (R2 = .52), which was supported by joint survivorship functions. The frequency and speed of movement execution should be evaluated in parallel to loading dose for the management of low-back training exposures.

Spine biomechanical testing methodologies: The controversy of consensus vs scientific evidence

Biomechanical testing methodologies for the spine have developed over the past 50 years. During that time, there have been several paradigm shifts with respect to techniques. These techniques evolved by incorporating state-of-the-art engineering principles, in vivo measurements, anatomical structure-function relationships, and the scientific method. Multiple parametric studies have focused on the effects that the experimental technique has on outcomes. As a result, testing methodologies have evolved, but there are no standard testing protocols, which makes the comparison of findings between experiments difficult and conclusions about in vivo performance challenging. In 2019, the international spine research community was surveyed to determine the consensus on spine biomechanical testing and if the consensus opinion was consistent with the scientific evidence. More than 80 responses to the survey were received. The findings of this survey confirmed that while some methods have been commonly adopted, not all are consistent with the scientific evidence. This review summarizes the scientific literature, the current consensus, and the authors' recommendations on best practices based on the compendium of available evidence.

A microstructure-based modeling approach to assess aging-sensitive mechanics of human intervertebral disc

BACKGROUND AND OBJECTIVE

The human body soft tissues are hierarchic structures interacting in a complex manner with the surrounding biochemical environment. The loss of soft tissues functionality with age leads to more vulnerability regarding to the external mechanical loadings and increases the risk of injuries. As a main example of the human body soft tissues, the intervertebral disc mechanical response evolution with age is explored. Although the age-dependence of the intervertebral disc microstructure is a well-known feature, no noticeable age effect on the disc stiffness is evidenced in the in-vitro experimental studies of the literature. So, if the disc intrinsic mechanics remains constant, how to explain the correlation of disc degeneration and disc functionality loss with age.

METHODS

A microstructure-based modeling approach was developed to assess in-silico the aging-sensitive mechanics of human intervertebral disc. The model considers the relationship between stress/volumetric macro-response and microstructure along with effective age effects acting at the lamellar and multi-lamellar scales. The stress-stretch and transversal responses of the different disc regions were computed for various age groups (13-18, 36, 58, 69 and 82 years old) and their evolution with age was studied.

RESULTS

While matching with in-vitro experimental data, the predicted stiffness was found to increase while passing from adolescent young discs to mature older discs and then to remain almost constant for the rest of life. Important age-related changes in the disc transversal behavior were also predicted affecting the flexibility of the disc, changing its volumetric behavior, and modifying its dimensions.

CONCLUSION

The developed approach was found able to bring new conclusions about age-dependent mechanical properties including regional dependency. The disc mechanics in terms of rigidity, radial and axial transversal responses were found to alter going from adolescent to middle age where the disc reaches a certain maturity. After reaching maturity, the mechanical properties undergo very slight changes until becoming almost constant with age.

Interlamellar matrix governs human annulus fibrosus multiaxial behavior

Establishing accurate structure–property relationships for intervertebral disc annulus fibrosus tissue is a fundamental task for a reliable computer simulation of the human spine but needs excessive theoretical-numerical-experimental works. The difficulty emanates from multiaxiality and anisotropy of the tissue response along with regional dependency of a complex hierarchic structure interacting with the surrounding environment. We present a new and simple hybrid microstructure-based experimental/modeling strategy allowing adaptation of animal disc model to human one. The trans-species strategy requires solely the basic knowledge of the uniaxial circumferential response of two different animal disc regions to predict the multiaxial response of any human disc region. This work demonstrates for the first time the determining role of the interlamellar matrix connecting the fibers-reinforced lamellae in the disc multiaxial response. Our approach shows encouraging multiaxial predictive capabilities making it a promising tool for human spine long-term prediction.

Parameter identification for the simulation of the periodontal ligament during the initial phase of orthodontic tooth movement

The paper is concerned with simulation of the periodontal ligament response to force in the initial phase of orthodontic tooth movement. This is based on two previous investigations, a in vitro experiment with specimens of porcine mandibular premolars and a in vivo experiment on human upper first incisors. For the curve fit of the in vitro experiment a model function, assuming viscoelasticity, was introduced. The viscoelastic model function was augmented by a ramp rise time term, to account for observed dependence of the response on actuator velocity, and a previous load history term, to account for the effect of the previous tests on the current test. The correlation coefficient of a curve fit for all tests grouped together was R2=0.98. Next, a curve fit of the in vivo experiment was done. Good correlation was found for a simplified model function, without viscoelastic term (R2=0.96). For both tests, in vitro and in vivo, the ramp rise time term improved correlation. A finite element model of the specimen of the in vitro experiment was created. For the PDL a hyperelastic constitutive model for compressible material was used and model parameters were identified. The present work indicates that the macroscopic response of the periodontal ligament to an external load can be simulated with a poro-visco-hyperelastic model. The simulation showed that poroelastic behaviour will gradually cease when viscoelastic relaxation progresses. This followed also from dimensionless analysis. As a consequence, for slow loading, or if initial response to fast loading is not of interest, a visco-hyperelastic model may suffice. To identify parameters of the finite element model several optimisation problems were solved. A model function, which can be regarded as a reduced order model, allowed a full factorial experiment (analysis) at low cost, to identify initial parameters. The thus found parameters were further refined with an optimum interpolation meta-model. That is, for limited number of parameter combinations the response was simulated with the finite element model and a refined parameter study was conducted by means of optimal interpolation. The thus found optimal parameters were verified by simulation with the finite element model. Optimal interpolation is computationally cheap, which allowed full factorial experiments at low cost.

How Osmoviscoelastic Coupling Affects Recovery of Cyclically Compressed Intervertebral Disc

STUDY DESIGN

Osmoviscoelastic behavior of cyclically loaded cervical intervertebral disc.

OBJECTIVE

The aim of this study was to evaluate in vitro the effects of physiologic compressive cyclic loading on the viscoelastic properties of cervical intervertebral disc and, examine how the osmoviscoelastic coupling affects time-dependent recovery of these properties following a long period of unloading.

SUMMARY OF BACKGROUND DATA

The human neck supports repetitive loadings during daily activities and recovery of disc mechanics is essential for normal mechanical function. However, the response of cervical intervertebral disc to cyclic loading is still not very well defined. Moreover, how loading history conditions could affect the time-dependent recovery is still unclear.

METHODS

Ten thousand cycles of compressive loading, with different magnitudes and saline concentrations of the surrounding fluid bath, are applied to 8 motion segments (composed by 2 adjacent vertebrae and the intervening disc) extracted from the cervical spines of mature sheep. Subsequently, specimens are hydrated during 18 hours of unloading. The viscoelastic disc responses, after cyclic loading and recovery phase, are characterized by relaxation tests.

RESULTS

Viscoelastic behaviors are significantly altered following large number of cyclic loads. Moreover, after 18-hour recovery period in saline solution at reference concentration (0.15 mol/L), relaxation behaviors were fully restored. Nonetheless, full recovery is not obtained whether the concentration of the surrounding fluid, that is, hypo-, iso-, or hyper-osmotic conditions.

CONCLUSION

Cyclic loading effects and full recovery of viscoelastic behavior after hydration at iso-osmotic condition (0.15 mol/L) are governed by osmotic attraction of fluid content in the disc due to imbalance between the external load and the swelling pressure of the disc. After removal of the load, the disc recovers its viscoelastic properties following period of rest. Nevertheless, the viscoelastic recovery is a chemically activated process and its dependency on saline concentration is governed by fluid flow due to imbalance of ions between the disc tissues and the surrounding fluid.

LEVEL OF EVIDENCE

3.

Nonlinear stress-dependent recovery behavior of the intervertebral disc

The intervertebral disc exhibits complex mechanics due to its heterogeneous structure, inherent viscoelasticity, and interstitial fluid-matrix interactions. Sufficient fluid flow into the disc during low loading periods is important for maintaining mechanics and nutrient transport. However, there is a lack of knowledge on the effect of loading magnitude on time-dependent recovery behavior and the relative contribution of multiple recovery mechanisms during recovery. In most studies that have evaluated disc recovery behavior, a single load condition has been considered, making it difficult to compare findings across studies. Hence, the objective of this study was to quantify unloaded disc recovery behavior after compressive creep loading under a wide range of physiologically relevant stresses (0.2-2 MPa). First, the repeatability of disc recovery behavior was assessed. Once repeatable recovery behavior was confirmed, each motion segment was subject to three cycles of creep-recovery loading, where each cycle consisted of a 24-h creep at a pre-assigned load (100, 200, 300, 600, 900, or 1200 N), followed by an 18-h recovery period at a nominal load (10 N). Results showed that disc recovery behavior was strongly influenced by the magnitude of loading. The magnitude of instantaneous and time-dependent recovery deformations increased nonlinearly with an increase in compressive stress during creep. In conclusion, this study highlights that elastic deformation, intrinsic viscoelasticity, and poroelasticity all have substantial contributions to disc height recovery during low loading periods. However, their relative contributions to disc height recovery largely depend on the magnitude of loading. While loading history does not influence the contribution of the short-term recovery, the contribution of long-term recovery is highly sensitive to loading magnitude.

Interlamellar-induced time-dependent response of intervertebral disc annulus: A microstructure-based chemo-viscoelastic model

The annulus fibrosus of the intervertebral disc exhibits an unusual transversal behavior for which a constitutive representation that considers as well regional effect, chemical sensitivity and time-dependency has not yet been developed, and it is hence the aim of the present contribution. A physically-based model is proposed by introducing a free energy function that takes into account the actual disc annulus structure in relation with the surrounding biochemical environment. The response is assumed to be dominated by the viscoelastic contribution of the extracellular matrix, the elastic contribution of the oriented collagen fibers and the osmo-induced volumetric contribution of the internal fluid content variation. The regional dependence of the disc annulus response due to variation in fibers content/orientation allows a micromechanical treatment of the soft tissue. A finite element model of the annulus specimen is designed while taking into consideration the 'interlamellar' ground substance zone between lamellae of the layered soft tissue. The kinetics is designed using full-field strain measurements performed on specimens extracted from two disc annulus regions and tested under different osmotic conditions. The time-dependency of the tissue response is reported on stress-free volumetric changes, on hysteretic stress and transversal strains during quasi-static stretching at different strain-rates and on their temporal changes during an interrupted stretching. Considering the effective contributions of the internal fluid transfer and the extracellular matrix viscosity, the microstructure-based chemo-mechanical model is found able to successfully reproduce the significant features of the macro-response and the unusual transversal behavior including the strong regional dependency from inner to outer parts of the disc: Poisson's ratio lesser than 0 (auxetic) in lamellae plane, higher than 0.5 in fibers plane, and their temporal changes towards usual values (between 0 and 0.5) at chemo-mechanical equilibrium. The underlying time-dependent mechanisms occurring in the tissue are analyzed via the local numerical fields and important insights about the effective role of the interlamellar zone are revealed for the different disc localizations. STATEMENT OF SIGNIFICANCE: The structural complexity of the annulus fibrosus has only been appreciated through recent experimental contributions and a constitutive representation that considers as well regional effect, chemical sensitivity and time-dependency of the unusual transversal behavior has not yet been developed. Here, a microstructure-based chemo-viscoelastic model is developed to highlight the interlamellar-induced time-dependent response by means of a two-scale strategy. The model provides important insights about the origin of the time-dependent phenomena in disc annulus along with regional dependency, essential for understanding disc functionality.

Automatic detection System for Degenerative Disk and simulation for artificial disc replacement surgery in the spine

This paper presents an automatic detection system for the degenerative disc and simulates the artificial disc replacement surgery in the spine. It starts by visualizing the Digital Imaging and Communications in Medicine (DICOM) in three views and then identifies regions for the lumber or sacral, followed by applying an adaptive threshold and modified region growing techniques to separate the spine from the surrounding tissues, organs, and bones. Then, segment each vertebral in the spine using K-means clustering algorithm. Then, the distance of intervertebral space was computed to detect the degenerative disk automatically in the lumbar region. Finally, the replacement surgery simulation occurred to design the disc with exact dimensions and identifies the appropriate place that will be implanted in the spine. The proposed system was tested on ten 3D computed tomography (CT) images downloaded from the Cancer Imaging Archive (TCIA). It was accessible to implant the artificial disk instead of the degenerative disk in the space that was specified automatically by the proposed system. The experimental outcomes positively demonstrate that the intended system is efficacious in the detection of degenerative disc and replacement of artificial disc simulation.

Recovering the mechanical properties of denatured intervertebral discs through Platelet-Rich Plasma therapy

Degenerative disc disease is one of the most common causes of low back pain instigating huge socioeconomic costs and posing an immense burden on healthcare systems worldwide. New therapeutic approaches to damaged intervertebral discs are therefore of great interest. Platelet-Rich Plasma (PRP) has been proposed for the repair and regeneration of degenerated discs, but there remains a knowledge gap regarding its effectiveness and influence on disc material properties. The objective of this study was to investigate and quantify the material properties of intact, denatured, and PRP treated discs. A systematic methodology was established in the process, where ex-vivo experiments were conducted and material properties were extracted using an inverse finite element approach. The results showed that PRP is able to recover the mechanical properties of denatured discs, thereby providing a promising effective therapeutic modality.

Stiffening of the nucleus pulposus upon axial loading of the intervertebral disc: An experimental in situ study

Mechanical loading is inherently related to the function and degeneration of the intervertebral disc. We present a series of experiments aimed at measuring the effect of a loading/unloading cycle of the intervertebral disc on the mechanical properties of the nucleus pulposus. The study relies on our new minimally invasive microindenter, which allows us to quantify the storage and loss moduli of the nucleus pulposus by inserting an optomechanical probe in an intact (resected) intervertebral disk through the annulus fibrosis via a small needle. Our results indicate that, under the influence of compressive loading, the nucleus pulposus exhibits a more solid-like behavior.

Matrisome Profiling During Intervertebral Disc Development And Ageing

Intervertebral disc (IVD) degeneration is often the cause of low back pain. Degeneration occurs with age and is accompanied by extracellular matrix (ECM) depletion, culminating in nucleus pulpous (NP) extrusion and IVD destruction. The changes that occur in the disc with age have been under investigation. However, a thorough study of ECM profiling is needed, to better understand IVD development and age-associated degeneration. As so, iTRAQ LC-MS/MS analysis of foetus, young and old bovine NPs, was performed to define the NP matrisome. The enrichment of Collagen XII and XIV in foetus, Fibronectin and Prolargin in elder NPs and Collagen XI in young ones was independently validated. This study provides the first matrisome database of healthy discs during development and ageing, which is key to determine the pathways and processes that maintain disc homeostasis. The factors identified may help to explain age-associated IVD degeneration or constitute putative effectors for disc regeneration.

Laser-induced micropore formation and modification of cartilage structure in osteoarthritis healing

Pores are vital for functioning of avascular tissues. Laser-induced pores play an important role in the process of cartilage regeneration. The aim of any treatment for osteoarthritis is to repair hyaline-type cartilage. The aims of this study are to answer two questions: (1) How do laser-assisted pores affect the cartilaginous cells to synthesize hyaline cartilage (HC)? and (2) How can the size distribution of pores arising in the course of laser radiation be controlled? We have shown that in cartilage, the pores arise predominately near chondrocytes, which promote nutrition of cells and signal molecular transfer that activates regeneration of cartilage. In vivo laser treatment of damaged cartilage of miniature pig joints provides cellular transformation and formation of HC. We propose a simple model of pore formation in biopolymers that paves the way for going beyond the trial-and-error approach when choosing an optimal laser treatment regime. Our findings support the approach toward laser healing of osteoarthritis.

Static axial overloading primes lumbar caprine intervertebral discs for posterior herniation

Introduction

Lumbar hernias occur mostly in the posterolateral region of IVDs and mechanical loading is an important risk factor. Studies show that dynamic and static overloading affect the nucleus and annulus of the IVD differently. We hypothesize there is also variance in the effect of overloading on the IVD’s anterior, lateral and posterior annulus, which could explain the predilection of herniations in the posterolateral region. We assessed the regional mechanical and cellular responses of lumbar caprine discs to dynamic and static overloading.

Material and methods

IVDs (n = 125) were cultured in a bioreactor and subjected to simulated-physiological loading (SPL), high dynamic (HD), or high static (HS) overloading. The effect of loading was determined in five disc regions: nucleus, inner-annulus and anterior, lateral and posterior outer-annulus. IVD height loss and external pressure transfer during loading were measured, cell viability was mapped and quantified, and matrix integrity was assessed.

Results

During culture, overloaded IVDs lost a significant amount of height, yet the distribution of axial pressure remained unchanged. HD loading caused cell death and disruption of matrix in all IVD regions, whereas HS loading particularly affected cell viability and matrix integrity in the posterior region of the outer annulus.

Conclusion

Axial overloading is detrimental to the lumbar IVD. Static overloading affects the posterior annulus more strongly, while the nucleus is relatively spared. Hence, static overloading predisposes the disc for posterior herniation. These findings could have implications for working conditions, in particular of sedentary occupations, and the design of interventions aimed at prevention and treatment of early intervertebral disc degeneration.

A regenerative approach towards recovering the mechanical properties of degenerated intervertebral discs: Genipin and platelet-rich plasma therapies

Degenerative disc disease, associated with discrete structural changes in the peripheral annulus and vertebral endplate, is one of the most common pathological triggers of acute and chronic low back pain, significantly depreciating an individual's quality of life and instigating huge socioeconomic costs. Novel emerging therapeutic techniques are hence of great interest to both research and clinical communities alike. Exogenous crosslinking, such as Genipin, and platelet-rich plasma therapies have been recently demonstrated encouraging results for the repair and regeneration of degenerated discs, but there remains a knowledge gap regarding the quantitative degree of effectiveness and particular influence on the mechanical properties of the disc. This study aimed to investigate and quantify the material properties of intact (N = 8), trypsin-denatured (N = 8), Genipin-treated (N = 8), and platelet-rich plasma-treated (N = 8) discs in 32 porcine thoracic motion segments. A poroelastic finite element model was used to describe the mechanical properties during different treatments, while a meta-model analytical approach was used in combination with ex vivo experiments to extract the poroelastic material properties. The results revealed that both Genipin and platelet-rich plasma are able to recover the mechanical properties of denatured discs, thereby affording promising therapeutic modalities. However, platelet-rich plasma-treated discs fared slightly, but not significantly, better than Genipin in terms of recovering the glycosaminoglycans content, an essential building block for healthy discs. In addition to investigating these particular degenerative disc disease therapies, this study provides a systematic methodology for quantifying the detailed poroelastic mechanical properties of intervertebral disc.

Biomechanical analysis of press-extension technique on degenerative lumbar with disc herniation and staggered facet joint

This study investigates the effect of a new Chinese massage technique named “press-extension” on degenerative lumbar with disc herniation and facet joint dislocation, and provides a biomechanical explanation of this massage technique. Self-developed biomechanical software was used to establish a normal L1–S1 lumbar 3D FE model, which integrated the spine CT and MRI data-based anatomical structure. Then graphic technique is utilized to build a degenerative lumbar FE model with disc herniation and facet joint dislocation. According to the actual press-extension experiments, mechanic parameters are collected to set boundary condition for FE analysis. The result demonstrated that press-extension techniques bring the annuli fibrosi obvious induction effect, making the central nucleus pulposus forward close, increasing the pressure in front part. Study concludes that finite element modelling for lumbar spine is suitable for the analysis of press-extension technique impact on lumbar intervertebral disc biomechanics, to provide the basis for the disease mechanism of intervertebral disc herniation using press-extension technique.

High Rotation Rate Behavior of Cervical Spine Segments in Flexion and Extension

Numerical finite element (FE) models of the neck have been developed to simulate occupant response and predict injury during motor vehicle collisions. However, there is a paucity of data on the response of young cervical spine segments under dynamic loading in flexion and extension, which is essential for the development or validation of tissue-level FE models. This limitation was identified during the development and validation of the FE model used in this study. The purpose of this study was to measure the high rotation rate loading response of human cervical spine segments in flexion and extension, and to investigate a new tissue-level FE model of the cervical spine with the experimental data to address a limitation in available data. Four test samples at each segment level from C2–C3 to C7–T1 were dissected from eight donors and were tested to 10 deg of rotation at 1 and 500 deg/s in flexion and extension using a custom built test apparatus. There was strong evidence (p < 0.05) of increased stiffness at the higher rotation rate above 4 deg of rotation in flexion and at 8 deg and 10 deg of rotation in extension. Cross-correlation software, Cora, was used to evaluate the fit between the experimental data and model predictions. The average rating was 0.771, which is considered to demonstrate a good correlation to the experimental data.

Intervertebral disc creep behavior assessment through an open source finite element solver

Degenerative Disc Disease (DDD) is one of the largest health problems faced worldwide, based on lost working time and associated costs. By means of this motivation, this work aims to evaluate a biomimetic Finite Element (FE) model of the Intervertebral Disc (IVD). Recent studies have emphasized the importance of an accurate biomechanical modeling of the IVD, as it is a highly complex multiphasic medium. Poroelastic models of the disc are mostly implemented in commercial finite element packages with limited access to the algorithms. Therefore, a novel poroelastic formulation implemented on a home-developed open source FE solver is briefly addressed throughout this paper. The combination of this formulation with biphasic osmotic swelling behavior is also taken into account. Numerical simulations were devoted to the analysis of the non-degenerated human lumbar IVD time-dependent behavior. The results of the tests performed for creep assessment were inside the scope of the experimental data, with a remarkable improvement of the numerical accuracy when compared with previously published results obtained with ABAQUS(®). In brief, this in-development open-source FE solver was validated with literature experimental data and aims to be a valuable tool to study the IVD biomechanics and DDD mechanisms.

Molecular immunotherapy might shed a light on the treatment strategies for disc degeneration and herniation

Despite surgical discectomy is one of the most effective treatments for intervertebral disc degeneration and lumbar disc herniation, a number of patients still complain of reserved low back pain, sciatica and numbness post-operatively with decreased life quality. Sciatica in patients with disc herniation is not only due to mechanical compression from herniated nucleus pulposus, but chemical and immunity agents. The intervertebral disc is composed of annulus fibrosus in the wedge and gelatinous nucleus pulposus in the centre with cartilage endplate sandwiched. Similar to other immune privilege organs, human intervertebral disc is one of the biggest avascular structures with FasL expression. Moreover, FasL-Fas and TRAIL death pathways might play roles in the machinery of immune privilege of the disc. We found that down-regulated miR-155 promotes Fas-mediated apoptosis in disc degeneration. Furthermore, once exposed to human immune system, nucleus pulposus can activate multiple specific and non-specific immune responses with cellular and fluid immune cells and molecules involved. Taken together, we hypothesize that a combined molecular immunotherapy with local and systemic immunity regulators might shed a novel light on the treatment strategies for disc degeneration and herniation.

A meta-model analysis of a finite element simulation for defining poroelastic properties of intervertebral discs

Finite element analysis is an effective tool to evaluate the material properties of living tissue. For an interactive optimization procedure, the finite element analysis usually needs many simulations to reach a reasonable solution. The meta-model analysis of finite element simulation can be used to reduce the computation of a structure with complex geometry or a material with composite constitutive equations. The intervertebral disc is a complex, heterogeneous, and hydrated porous structure. A poroelastic finite element model can be used to observe the fluid transferring, pressure deviation, and other properties within the disc. Defining reasonable poroelastic material properties of the anulus fibrosus and nucleus pulposus is critical for the quality of the simulation. We developed a material property updating protocol, which is basically a fitting algorithm consisted of finite element simulations and a quadratic response surface regression. This protocol was used to find the material properties, such as the hydraulic permeability, elastic modulus, and Poisson's ratio, of intact and degenerated porcine discs. The results showed that the in vitro disc experimental deformations were well fitted with limited finite element simulations and a quadratic response surface regression. The comparison of material properties of intact and degenerated discs showed that the hydraulic permeability significantly decreased but Poisson's ratio significantly increased for the degenerated discs. This study shows that the developed protocol is efficient and effective in defining material properties of a complex structure such as the intervertebral disc.

Flow-induced deformation of poroelastic tissues and gels: a new perspective on equilibrium pressure-flow-thickness relations

Hydrostatic pressure-driven flows through soft tissues and gels cause deformations of the solid network to occur, due to drag from the flowing fluid. This phenomenon occurs in many contexts including physiological flows and infusions through soft tissues, in mechanically stimulated engineered tissues, and in direct permeation measurements of hydraulic permeability. Existing theoretical descriptions are satisfactory in particular cases, but none provide a description which is easy to generalize for the design and interpretation of permeation experiments involving a range of different boundary conditions and gel properties. Here a theoretical description of flow-induced permeation is developed using a relatively simple approximate constitutive law for strain-dependent permeability and an assumed constant elastic modulus, using dimensionless parameters which emerge naturally. Analytical solutions are obtained for relationships between fundamental variables, such as flow rate and pressure drop, which were not previously available. Guidelines are provided for assuring that direct measurements of hydraulic permeability are performed accurately, and suggestions emerge for alternative measurement protocols. Insights obtained may be applied to interpretation of flow-induced deformation and related phenomena in many contexts.

Hydraulic permeability of multilayered collagen gel scaffolds under plastic compression-induced unidirectional fluid flow

Under conditions of free fluid flow, highly hydrated fibrillar collagen gels expel fluid and undergo gravity driven consolidation (self-compression; SC). This process can be accelerated by the application of a compressive stress (plastic compression; PC) in order to generate dense collagen scaffolds for tissue engineering. To define the microstructural evolution of collagen gels under PC, this study applied a two-layer micromechanical model that was previously developed to measure hydraulic permeability (k) under SC. Radially confined PC resulted in unidirectional fluid flow through the gel and the formation of a dense lamella at the fluid expulsion boundary which was confirmed by confocal microscopy of collagen immunoreactivity. Gel mass loss due to PC and subsequent SC were measured and applied to Darcy's law to calculate the thickness of the lamella and hydrated layer, as well as their relative permeabilities. Increasing PC level resulted in a significant increase in mass loss fraction and lamellar thickness, while the thickness of the hydrated layer dramatically decreased. Permeability of lamella also decreased from 1.8×10(-15) to 1.0×10(-15) m(2) in response to an increase in PC level. Ongoing SC, following PC, resulted in a uniform decrease in mass loss and k with increasing PC level and as a function SC time. Experimental k data were in close agreement with those estimated by the Happel model. Calculation of average k values for various two-layer microstructures indicated that they each approached 10(-15)-10(-14) m(2) at equilibrium. In summary, the two-layer micromechanical model can be used to define the microstructure and permeability of multi-layered biomimetic scaffolds generated by PC.

The impact of posture and prolonged cyclic compressive loading on vertebral joint mechanics

STUDY DESIGN

An in vitro biomechanics investigation exposing porcine functional spinal units (FSUs) to submaximal cyclic or static compressive forces while in a flexed, neutral, or extended posture.

OBJECTIVE

To investigate the combined effect of cyclically applied compressive force (e.g., vibration) and postural deviation on intervertebral joint mechanics.

SUMMARY OF BACKGROUND DATA

Independently, prolonged vibration exposure and non-neutral postures are known risk factors for development of low back pain and injury. However, there is limited basic scientific evidence to explain how the risk of low back injury from vibration exposure is modified by other mechanical factors. This work examined the influence of static postural deviation on vertebral joint height loss and compressive stiffness under cyclically applied compressive force.

METHODS

Forty-eight FSUs, consisting of 2 adjacent vertebrae, ligaments, and the intervening intervertebral disc were included in the study. Each specimen was randomized to 1 of 3 experimental posture conditions (neutral, flexed, or extended) and assigned to 1 of 2 loading protocols, consisting of (1) cyclic (1500 ± 1200 N applied at 5 Hz using a sinusoidal waveform, resulting in 0.2 g rms acceleration) or (2) 1500 N of static compressive force. RESULTS.: As expected, FSU height loss followed a typical first-order response in both the static and cyclic loading protocols, with the majority (~50%) of the loss occurring in the first 20 minutes of testing. A significant interaction between posture and loading protocol (P < 0.001) was noted in the magnitude of FSU height loss. Subsequent analysis of simple effects revealed significant differences between cyclic and static loading protocols in both a neutral (P = 0.016) and a flexed posture (P < 0.0001). No significant differences (P = 0.320) were noted between pre/postmeasurements of FSU compressive stiffness.

CONCLUSION

Posture is an important mechanical factor to consider when assessing the risk of injury from cyclic loading to the lumbar spine.

Microstructural stress relaxation mechanics in functionally different tendons

Tendons experience widely varying loading conditions in vivo. They may be categorised by their function as either positional tendons, which are used for intricate movements and experience lower stress, or as energy storage tendons which act as highly stressed springs during locomotion. Structural and compositional differences between tendons are thought to enable an optimisation of their properties to suit their functional environment. However, little is known about structure-function relationships in tendon. This study adopts porcine flexor and extensor tendon fascicles as examples of high stress and low stress tendons, comparing their mechanical behaviour at the micro-level in order to understand their stress relaxation response. Stress-relaxation was shown to occur predominantly through sliding between collagen fibres. However, in the more highly stressed flexor tendon fascicles, more fibre reorganisation was evident when the tissue was exposed to low strains. By contrast, the low load extensor tendon fascicles appears to have less capacity for fibre reorganisation or shearing than the energy storage tendon, relying more heavily on fibril level relaxation. The extensor fascicles were also unable to sustain loads without rapid and complete stress relaxation. These findings highlight the need to optimise tendon repair solutions for specific tendons, and match tendon properties when using grafts in tendon repairs.

Biomechanical disc culture system: feasibility study using rat intervertebral discs

A small-scale biomechanical disc culture system was designed to stimulate intervertebral disc (IVD) 'motion segment' in culture environment with load-controlled compression and combined load (compression+shear). After 7 days of diurnal mechanical loading, cell viability of discs stimulated with static compression load (0.25 MPa) and static combined load (compression (0.25 MPa)+shear (1.5N)) were similar (>90 per cent) to unloaded controls. Mechanically stimulated discs showed decrease in static/dynamic moduli, early stress relaxation, and loss of disc height after 7 days of diurnal loading. Histological data of discs indicated load-induced transformations that were not apparent in controls. The feasibility of studying the mechanobiology of intact IVD as a motion segment was demonstrated. Media conditioning (improve tissue stability in long-term culture) and application of biochemical gene expression assays (differential tissue response to types of mechanical stimulation) are proposed as future improvements. The study suggests that the limitations in studying mechanobiology of IVD pathology in vitro can be overcome and it is possible to understand the physiologically relevant mechanism of IVD pathology.

The effect of degenerative morphological changes of the intervertebral disc on the lumbar spine biomechanics: a poroelastic finite element investigation

Intervertebral disc degeneration involves changes in the spinal anatomical structures. The mechanical relevance of the following changes was investigated: disc height, endplate sclerosis, disc water content, permeability and depressurisation. A poroelastic nonlinear finite element model of the L4-L5 human spine segments was employed. Loads represented a daily cycle (500 N compression combined with flexion-extension motion for 16 h followed by 200 N compression for 8 h). In non-degenerative conditions, the model predicted a diurnal axial displacement of 1.32 mm and a peak intradiscal pressure of 0.47 MPa. Axial displacement, facet force and range of motion in flexion-extension are decreased by decreasing disc height. By decreasing the initial water content, axial displacement, facet force and fluid loss were all reduced. Endplate sclerosis did not have a significant influence on the calculated results. Depressurisation determined an increase of the disc effective stress, possibly inducing failure. Degenerative instability was not calculated in any simulations.

Human disc nucleus properties and vertebral endplate permeability

STUDY DESIGN

Experimental quantification of relationships between vertebral endplate morphology, permeability, disc cell density, glycosaminoglycan (GAG) content, and degeneration in samples harvested from human cadaveric spines.

OBJECTIVE

To test the hypothesis that variation in endplate permeability and porosity contributes to changes in intervertebral disc cell density and overall degeneration.

SUMMARY OF BACKGROUND DATA

Cells within the intervertebral disc are dependent on diffusive exchange with capillaries in the adjacent vertebral bone. Previous findings suggest that blocked routes of transport negatively affect disc quality, yet there are no quantitative relationships between human vertebral endplate permeability, porosity, cell density, and disc degeneration. Such relationships would be valuable for clarifying degeneration risk factors and patient features that may impede efforts at disc tissue engineering.

METHODS

Fifty-one motion segments were harvested from 13 frozen cadaveric human lumbar spines (32-85 years) and classified for degeneration using the magnetic resonance imaging-based Pfirrmann scale. A cylindrical core was harvested from the center of each motion segment that included vertebral bony and cartilage endplates along with adjacent nucleus tissue. The endplate mobility, a type of permeability, was measured directly using a custom-made permeameter before and after the cartilage endplate was removed. Cell density within the nucleus tissue was estimated using the picogreen method, while the nuclear GAG content was quantified using the dimethylmethylene blue technique. Specimens were imaged at 8 μm resolution using microCT; bony porosity was calculated. Analysis of variance, linear regression, and multiple comparison tests were used to analyze the data. RESULTS.: Nucleus cell density increased as the disc height decreased (R² = 0.13; P = 0.01) but was not related to subchondral bone porosity (P > 0.5), total mobility (P > 0.4), or age (P > 0.2). When controlling for disc height, however, a significant, negative effect of age on cell density was observed (P = 0.03). In addition to this, GAG content decreased with age nonlinearly (R² = 0.83, P < 0.0001) and a cell function measurement, GAGs/cell, decreased with degeneration (R² = 0.24; P < 0.0001). Total mobility (R² = 0.14; P < 0.01) and porosity (R² = 0.1, P < 0.01) had a positive correlation with age.

CONCLUSION

Although cell density increased with degeneration, cell function indicated that GAGs/cell decreased. Because permeability and porosity increase with age and degeneration, this implies that cell dysfunction, rather than physical barriers to transport, accelerates disc disease.

Axial creep loading and unloaded recovery of the human intervertebral disc and the effect of degeneration

The intervertebral disc maintains a balance between externally applied loads and internal osmotic pressure. Fluid flow plays a key role in this process, causing fluctuations in disc hydration and height. The objectives of this study were to quantify and model the axial creep and recovery responses of nondegenerate and degenerate human lumbar discs. Two experiments were performed. First, a slow compressive ramp was applied to 2000 N, unloaded to allow recovery for up to 24 h, and re-applied. The linear-region stiffness and disc height were within 5% of the initial condition for recovery times greater than 8 h. In the second experiment, a 1000 N creep load was applied for four hours, unloaded recovery monitored for 24 h, and the creep load repeated. A viscoelastic model comprised of a "fast" and "slow" exponential response was used to describe the creep and recovery, where the fast response is associated with flow in the nucleus pulposus (NP) and endplate, while the slow response is associated with the annulus fibrosus (AF). The study demonstrated that recovery is 3-4X slower than loading. The fast response was correlated with degeneration, suggesting larger changes in the NP with degeneration compared to the AF. However, the fast response comprised only 10%-15% of the total equilibrium displacement, with the AF-dominated slow response comprising 40%-70%. Finally, the physiological loads and deformations and their associated long equilibrium times confirm that diurnal loading does not represent "equilibrium" in the disc, but that over time the disc is in steady-state.

Simulation of inhomogeneous rather than homogeneous poroelastic tissue material properties within disc annulus and nucleus better predicts cervical spine response: a C3-T1 finite element model analysis under compression and moment loadings

STUDY DESIGN

A finite element (FE) modeling of homogeneous and inhomogeneous poroelastic tissue material properties within disc anulus fibrosus (AF) and nucleus pulposus (NP).

OBJECTIVE

To test the hypothesis that simulation of inhomogeneous poroelastic tissue material properties within AF and NP quadrants, rather than homogeneous properties within regions of AF and NP without quadrants, would better predict the cervical spine biomechanics.

SUMMARY OF BACKGROUND DATA

In order to represent tissue swelling and creep deformation behavior more physiologically in FE models, disc poroelastic tissue material properties should be modeled appropriately. Past studies show an existence of inhomogeneous rather than homogeneous nature of the tissue properties in various quadrants of AF and NP, and this has been simulated in a single-segment FE lumbar model with only compression analysis. This article simulated these tissue properties in a multisegmental cervical spine and reported the results of both compression and moment loads.

METHODS

Two three-dimensional FE models of a C3-T1 segment were developed. Model I included homogeneous poroelastic tissue properties in AF and NP, whereas Model II included inhomogeneous poroelastic tissue properties in AF and NP quadrants. Biomechanical responses of the FE models under diurnal compression and moment loads were compared with corresponding in vivo published studies.

RESULTS

Model II with disc quadrant-based inhomogeneous poroelastic tissue properties predicted better, mainly in flexion and extension, than the Model I with homogeneous tissue properties when compared with the corresponding in vivo results, thereby confirming the current study hypothesis. Inhomogeneous tissue properties govern segmental behavior mainly during sagittal plane motions, with a root-mean-square difference of nearly 50% across the motion segments.

CONCLUSION

The current data justify the need to simulate inhomogeneous tissue properties within disc quadrants for any FE model analysis. Model II can be further used to understand the biomechanical effects of quadrant-based degenerative poroelastic tissue properties on cervical spine behavior. Future experiments are necessary to support the current study results.

Relative contributions of strain-dependent permeability and fixed charged density of proteoglycans in predicting cervical disc biomechanics: a poroelastic C5-C6 finite element model study

Disc swelling pressure (P(swell)) facilitated by fixed charged density (FCD) of proteoglycans (P(fcd)) and strain-dependent permeability (P(strain)) are of critical significance in the physiological functioning of discs. FCD of proteoglycans prevents any excessive matrix deformation by tissue stiffening, whereas strain-dependent permeability limits the rate of stress transfer from fluid to solid skeleton. To date, studies involving the modeling of FCD of proteoglycans and strain-dependent permeability have not been reported for the cervical discs. The current study objective is to compare the relative contributions of strain-dependent permeability and FCD of proteoglycans in predicting cervical disc biomechanics. Three-dimensional finite element models of a C5-C6 segment with three different disc compositions were analyzed: an SPFP model (strain-dependent permeability and FCD of proteoglycans), an SP model (strain-dependent permeability alone), and an FP model (FCD of proteoglycans alone). The outcomes of the current study suggest that the relative contributions of strain-dependent permeability and FCD of proteoglycans were almost comparable in predicting the physiological behavior of the cervical discs under moment loads. However, under compression, strain-dependent permeability better predicted the in vivo disc response than that of the FCD of proteoglycans. Unlike the FP model (least stiff) in compression, motion behavior of the three models did not vary much from each other and agreed well within the standard deviations of the corresponding in vivo published data. Flexion was recorded with maximum P(fcd) and P(strain), whereas minimum values were found in extension. The study data enhance the understanding of the roles played by the FCD of proteoglycans and strain-dependent permeability and porosity in determining disc tissue swelling behavior. Degenerative changes involving strain-dependent permeability and/or loss of FCD of proteoglycans can further be studied using an SPFP model. Future experiments are necessary to support the current study results.

The effects of exogenous crosslinking on hydration and fluid flow in the intervertebral disc subjected to compressive creep loading and unloading

STUDY DESIGN

In vitro study of genipin crosslinking effect on disc water content changes under compressive loading and unloading.

OBJECTIVE

To investigate the influence of collagen crosslinking on hydration and fluid flow in different regions of intact discs, and to evaluate the nutritional implications.

SUMMARY OF BACKGROUND DATA

Age-related reductions of nutrient supply and waste product removal are critically important factors in disc pathogenesis. Diffusion and fluid flow are blocked by subchondral bone thickening, cartilaginous endplate calcification, loss of hydrophilic proteoglycans, and clogging of anular pores by degraded matrix molecules. Previous studies demonstrated increased hydraulic permeability and macromolecular transport through crosslinked collagenous matrices. Genipin has also demonstrated the capability to increase retention of proteoglycans.

METHODS

A total of 57 bovine lumbar motion segments were divided randomly into phosphate buffered saline and 0.33% genipin-soaked treatment groups. Water content changes were measured using a mass-loss technique in 3 intervertebral disc regions following successive stages of compressive loading and unloading (post-treatment, after 1 hour 750 N compression, and after a subsequent 24-hour period of nominal loading). Net flow of fluid into or out of a region was determined from the percentage change in mean water content from successive groups.

RESULTS

Fluid flow to and from the nucleus doubled with genipin crosslinking. Relative to the buffer-only controls, overall net fluid flow increased 103% in the nucleus pulposus, 36% in the inner anulus, and was 31% less in the outer anulus of genipin treated discs.

CONCLUSION

The effects of genipin crosslinking on matrix permeability and proteoglycan retention can alter hydration levels and fluid flow in the intervertebral disc. Resulting increases in fluid flow, including a doubling of flow to and from the nucleus, could lead to enhanced nutritional inflow and waste product outflow for the disc, and may have implications for emerging cell-based therapies.

The mechanical response of the lumbar spine to different combinations of disc degenerative changes investigated using randomized poroelastic finite element models

Degeneration of the intervertebral disc is related to progressive changes in the disc tissue composition and morphology, such as water loss, disc height loss, endplate calcification, osteophytosis. These changes may be present separately or, more frequently, in various combinations. This work is aimed to the biomechanical investigation of a wide range of clinical scenarios of disc degeneration, in which the most common degenerative changes are present in various combinations. A poroelastic non-linear finite element model of the healthy L4-L5 human spine segment was employed and randomly scaled to represent ten spine segments from different individuals. Six different degenerative characteristics (water loss in the nucleus pulposus and annulus fibrosus; calcification and thickness reduction of endplate cartilage; disc height loss; osteophyte formation; diffuse sclerosis) were modeled in 30 randomly generated models, 10 for each overall degree of degeneration (mild, moderate, severe). For each model, a daily loading cycle including 8 h of rest, 16 h in the standing position with superimposed two flexion-extension motion cycles was simulated. A tendency to an increase of stiffness with progressing overall degeneration was observed, in compression, flexion and extension. Hence, instability for mild degeneration was not predicted. Facet forces and fluid loss decreased with disc degeneration. Nucleus, annulus and endplate degeneration, disc height loss, osteophytosis and diffuse sclerosis all induced a statistically significant decrease in the total daily disc height variation, facet force and flexibility in flexion-extension. Therefore, grading systems for disc degeneration should include all the degenerative changes considered in this work, since all of them had a significant influence on the spinal biomechanics.

Characterization of a biodegradable electrospun polyurethane nanofiber scaffold: Mechanical properties and cytotoxicity

The current study analyzes the biodegradation of a polycarbonate polyurethane scaffold intended for the growth of a tissue-engineered annulus fibrosus (AF) disc component. Electrospun scaffolds with random and aligned nanofiber configurations were fabricated using a biodegradable polycarbonate urethane with and without an anionic surface modifier (anionic dihydroxyl oligomer), and the mechanical behavior of the scaffolds was examined during a 4 week biodegradation study. Both the tensile strength and initial modulus of aligned scaffolds (sigma=14+/-1 MPa, E=46+/-3 MPa) were found to be higher than those of random fiber scaffolds (sigma=1.9+/-0.4 MPa, E=2.1+/-0.2 MPa) prior to degradation. Following initial wetting of the scaffold, the initial modulus of the aligned samples showed a significant decrease (dry: 46+/-3 MPa; pre-wetted: 9+/-1 MPa, p<0.001). The modulus remained relatively constant during the remainder of the 4 week incubation period (aligned at 4 weeks: 8.0+/-0.3 MPa). The tensile strength for aligned fiber scaffolds was affected in the same manner. Similar changes were not observed for the initial modulus of the random scaffold configuration. Biodegradation of the scaffold in the presence of cholesterol esterase (a monocyte derived enzyme) yielded a 0.5 mg week(-1) weight loss. The soluble and non-soluble degradation products were found to be non-toxic to bovine AF cells grown in vitro. The consistent rate of material degradation along with stable mechanical properties comparable to those of native AF tissue and the absence of cytotoxic effects make this polymer a suitable biomaterial candidate for further investigation into its use for tissue-engineering annulus fibrosus.

Refinement of elastic, poroelastic, and osmotic tissue properties of intervertebral disks to analyze behavior in compression

Intervertebral disks support compressive forces because of their elastic stiffness as well as the fluid pressures resulting from poroelasticity and the osmotic (swelling) effects. Analytical methods can quantify the relative contributions, but only if correct material properties are used. To identify appropriate tissue properties, an experimental study and finite element analytical simulation of poroelastic and osmotic behavior of intervertebral disks were combined to refine published values of disk and endplate properties to optimize model fit to experimental data. Experimentally, nine human intervertebral disks with adjacent hemi-vertebrae were immersed sequentially in saline baths having concentrations of 0.015, 0.15, and 1.5 M and the loss of compressive force at constant height (force relaxation) was recorded over several hours after equilibration to a 300-N compressive force. Amplitude and time constant terms in exponential force-time curve-fits for experimental and finite element analytical simulations were compared. These experiments and finite element analyses provided data dependent on poroelastic and osmotic properties of the disk tissues. The sensitivities of the model to alterations in tissue material properties were used to obtain refined values of five key material parameters. The relaxation of the force in the three bath concentrations was exponential in form, expressed as mean compressive force loss of 48.7, 55.0, and 140 N, respectively, with time constants of 1.73, 2.78, and 3.40 h. This behavior was analytically well represented by a model having poroelastic and osmotic tissue properties with published tissue properties adjusted by multiplying factors between 0.55 and 2.6. Force relaxation and time constants from the analytical simulations were most sensitive to values of fixed charge density and endplate porosity.

Rheology of intervertebral disc: an ex vivo study on the effect of loading history, loading magnitude, fatigue loading, and disc degeneration

STUDY DESIGN

An ex vivo biomechanical study on the rheological properties of healthy porcine and degenerated human intervertebral disc.

OBJECTIVE

To quantify the effect of loading history, loading magnitude, fatigue loading, and degeneration on disc rheology.

SUMMARY OF BACKGROUND DATA

Disc rheological parameters, i.e., the aggregate modulus (HA) and hydraulic permeability (k) regulate the mechanical and biologic function of disc. The knowledge of effects of loading condition and degeneration on disc rheology can be beneficial for the design of new disc/nucleus implants or therapy.

METHODS

The following 4 phases of experiments were conducted to find the changes of disc rheological properties: (1) Effect of loading history during 1-hour creep (640 N) and relaxation (20% strain) test. (2) Effect of loading magnitude (420 N vs. 640 N) during the creep test. (3) Effect of fatigue loading (420 N, 5 Hz for 0.5, 1, and 2 hours) on the creep loading behavior. (4) Difference of healthy porcine and degenerated human discs during creep loading. The experimental data were fitted with linear biphasic model.

RESULTS

The aggregate modulus increased but hydraulic permeability decreased during creep loading. The aggregate modulus decreased but the hydraulic permeability did not change significantly during relaxation loading. The higher creep loading increased the aggregate modulus but decreased the hydraulic permeability. The fatigue loading did not change the aggregate modulus significantly, but decreased hydraulic permeability. Comparing the degenerated human disc to the healthy porcine disc, the aggregate modulus was higher but the hydraulic permeability was lower.

CONCLUSION

The external loading and degeneration induce disc structural changes, e.g., the disc water content and interstitial matrix porosity, hence affect the disc rheological properties. The increase of aggregate modulus may be due to the reduction of disc hydration level, whereas the decrease of hydraulic permeability may be because of the shrinkage of disc matrix pores.

Pre-operative blood test for antibody to nucleus pulposus may distinguish types of lumbar disc herniation

The types of lumbar disc herniation not only play a key role in planning therapeutic procedure, but also appear to be connected with prognosis. However, it is difficult to distinguish the types preoperatively. To the best of our knowledge, the intervertebral disc is the largest avascular structure in the human body. Once the nucleus pulposus extrudes from the annulus fibrosus or combining with posterior longitudinal ligament, the immune system recognizes it as a foreign body, which may induce the autoimmune reaction ending with production of auto-antibody, that is so-called antibody to nucleus pulposus. The blood test for antibody to nucleus pulposus is ultimately accurate because the result is not relied upon the patient's complain but specific to herniated nucleus pulposus tissue. If this hypothesis is proved to be correct, it could be possible to improve the accuracy of distinguishing types of lumbar disc herniation just by a pre-operative blood test for antibody to nucleus pulposus.

The immediate effect of repeated loading on the compressive strength of young porcine lumbar spine

The human spine is exposed to repeated loading during daily activities and more extremely during sports. Despite this, there remains a lack of knowledge regarding the immediate effects on the spine due to this mode of loading. Age-specific spinal injury patterns has been demonstrated and this implies differences in reaction to load mode and load history The purpose of the present study was to investigate the impact of cyclic pre-loading on the biomechanical properties and fracture patterns of the adolescent spine in an experimental model. Eight functional spinal units from four young porcine spines were harvested. The functional spinal units were cyclic loaded with 20,000 cycles and then axially compressed to failure. The compression load at failure, ultimate stress and viscoelastic parameters were calculated. The functional spinal units were examined with plain radiography, computer tomography and MRI before and after the loading, and finally macroscopically and histologically. The median compression load at failure in this study was 8.3 kN (range 5.6-8.7 kN). The median deformation for all cases was 2.24 mm (range 2.30-2.7 mm) and stiffness was 3.45 N/mm (range 3.5-4.5 N/mm). A fracture was seen on radiograph in one case, on CT and macroscopically in seven, and on MRI and histologically in all eight cases. The cyclic loaded functional spinal units in the present study were not more sensitive to axial compression than non-cyclic loaded functional spinal units from young porcine. The endplate and the growth zone were the weakest part in the cyclic loaded functional spinal units. Disc signal reduction and disc height reduction was found on MRI. The E-modulus value found in this study was of the same order of magnitude as found by others using a porcine animal model.

Mechanical behavior of annulus fibrosus: a microstructural model of fibers reorientation

Experimental uniaxial tensile tests have been carried out on annulus tissue samples harvested on pig and lamb lumbar intervertebral discs. When subjecting the samples to loading cycles, the stress-strain curves exhibit strong nonlinearities and hysteresis. This particular behavior results from the anisotropic microstructure of annulus tissue composed of woven oriented collagen fibers embedded in the extracellular matrix. During uniaxial tension, the collagen fibers reorient toward the loading direction increasing its global stiffness. To describe this behavior, we propose a heuristic two-dimensional rheological model based on three mechanical and one geometrical characteristics. The latter one is the fibers orientation angle becoming the key parameter that govern the macroscopic mechanical behavior. The experimental results are used to identify the physical properties associated with the rheological model, leading to an accurate representation of the stress-strain curve over a complete loading cycle. In this framework, the fibers reorientation can solely account for the rigidity increase while the hysteresis is associated with liquid viscous flows through the matrix. Based on this representation, unusual coupling effects between strains and fluid flows can be observed, that would significantly affect the cell nutrients transport mechanisms.

Under-arm partial body weight unloading causes spinal elongation and vibration attenuation during treadmill walking

Whilst exercise is beneficial to those suffering from low back pain (LBP), spinal shrinkage and vibration during walking may aggravate the condition. This study investigates the effects of spinal unloading, by means of body-weight unloading (BWU), on spinal length and vibration response. Under-arm partial BWU (40% of bodyweight) was performed on 8 healthy males whilst walking on a treadmill for 1h, and compared to a control condition in the same participants. Motion analysis was used to track four reflective markers attached to the spine between C7 and the lumbar concavity at 100Hz, in 7s samples at regular intervals during the walk. A quintic polynomial was fitted to the coordinates in the sagittal plane, and sub-sectioned into three regions: the upper thoracic (UT), lower thoracic (LT) and upper lumbar (UL). The lengths of the curves were analysed in the time and frequency domains. The length of the spine increased by 4+/-2% (18mm) during the unloading condition with all regions showing an increase in length. The UL and LT regions lost length in the control condition, thereby exhibiting a significant interaction between unloading and time on region length (both P<0.05). In addition, compared to the control condition, the frequencies of the length changes were attenuated between 3 and 6Hz with unloading. Therefore, under-arm BWU facilitates spinal elongation, attenuates the frequency response of the spine in its resonant frequency regime and thus has potential benefits to the LBP population.

An extended biphasic model for charged hydrated tissues with application to the intervertebral disc

Finite element models for hydrated soft biological tissue are numerous but often exhibit certain essential deficiencies concerning the reproduction of relevant mechanical and electro-chemical responses. As a matter of fact, singlephasic models can never predict the interstitial fluid flow or related effects like osmosis. Quite a few models have more than one constituent, but are often restricted to the small-strain domain, are not capable of capturing the intrinsic viscoelasticity of the solid skeleton, or do not account for a collagen fibre reinforcement. It is the goal of this contribution to overcome these drawbacks and to present a thermodynamically consistent model, which is formulated in a very general way in order to reproduce the behaviour of almost any charged hydrated tissue. Herein, the Theory of Porous Media (TPM) is applied in combination with polyconvex Ogden-type material laws describing the anisotropic and intrinsically viscoelastic behaviour of the solid matrix on the basis of a generalised Maxwell model. Moreover, other features like the deformation-dependent permeability, the possibility to include inhomogeneities like varying fibre alignment and behaviour, or osmotic effects based on the simplifying assumption of Lanir are also included. Finally, the human intervertebral disc is chosen as a representative for complex soft biological tissue behaviour. In this regard, two numerical examples will be presented with focus on the viscoelastic and osmotic capacity of the model.

The strain-dependent osmotic pressure and stiffness of the bovine nucleus pulposus apportioned into ionic and non-ionic contributors

Elucidation of the load-bearing mechanism of the nucleus pulposus (NP) facilitates understanding of the mechanical and metabolic functioning of the intervertebral disc and provides key data for mathematical models. Negatively charged proteoglycans in the NP generate an ionic osmotic pressure, pi(i), which contributes to the tissue's resistance to load and, moreover, is the main mechanism by which the unloaded disc rehydrates. Functionally important, pi(i) has seldom been investigated in situ and, crucially, its variation with strain has not been reported. In a confined compression apparatus, we aimed to apportion the strain-dependent load-bearing mechanism of the NP at equilibrium to the tissue matrix and ionic osmotic pressure; and to determine whether any proteoglycan loss occurs during confined compression testing. Forty-eight confined compression experiments were conducted in isotonic (0.15M NaCl) and hypertonic (3.0 and 6.1M NaCl) external solutions in single and multiple step-strain protocols. The 6.1M NaCl external solution was needed to eliminate as much of the ionic effects as possible. The ionic osmotic pressure was well described by pi(i)=19.1lambda(-1.58) (R(2)=0.992), and was approximately 70% of the applied load at equilibrium, independent of lambda. The effective aggregate modulus, H(A)(eff), also increased with strain: H(A)(eff)=59.0lambda(-2.18). Concentrations of sulphated glycosaminoglycans were obtained for the samples tested in isotonic NaCl with no proteoglycan loss detected from the confined compression tests. These results highlight the non-linearity of the stress-strain response of NP tissue and the necessity to include a non-linear function for osmotic pressure in mathematical models of this tissue.

The effect of relative needle diameter in puncture and sham injection animal models of degeneration

STUDY DESIGN

Biomechanical study and literature review.

OBJECTIVES

To quantify the acute effect of needle diameter on the in vitro mechanical properties of cadaver lumbar discs in the rat and sheep. To review published in vivo animal studies and evaluate disc changes with respect to the relative needle size.

SUMMARY OF BACKGROUND DATA

There are many cases where a disc needle puncture or injection is applied to animal models: puncture injuries to induce degeneration, chemonucleolysis to induce degeneration, and delivery of disc therapies. It is not clear what role the size of the needle may have in the outcome.

METHODS

Mechanics were measured after sham phosphate buffered saline injection with a 27 G or 33 G needle in the rat and with a 27 G needle in the sheep. A literature review was performed to evaluate studies in which animal discs were treated with a needle puncture or a sham injection. For each study, the ratio of the needle diameter to disc height (needle:height) was calculated.

RESULTS

When the rat was injected with a 27 G needle (52% of disc height), the compression, tension, and neutral zone stiffnesses were 20% to 60% below preinjected values and the neutral zone length was 130% higher; when injected with a 33 G needle (26% of disc height), the only affected property was the neutral zone length, which was only 20% greater. When the sheep was injected with a 27 G needle (10% of disc height), none of the axial properties were different from intact, the torsion stiffness was not different, and the torque range was 15% smaller. Twenty-three in vivo studies in the rat, rabbit, dog, or sheep were reviewed. The disc changes depended on the ratio of needle diameter to disc height as follows: significant changes were not observed for needle:height less than 40%, although between 25% and 40% results were variable and some minor nonsignificant effects were observed, disc changes were universal for needle:height over 40%.

CONCLUSION

A needle puncture may directly alter mechanical properties via nucleus pulposus depressurization and/or anulus fibrosus damage, depending on the relative needle size. As more basic science research is aimed at treating disc degeneration via injection of therapeutic factors, these findings provide guidance in design of animal studies. Such studies should consider the relative needle size and include sham control groups to account for the potential effects of the needle injection.