Role of load history in intervertebral disc mechanics and intradiscal pressure generation

Solid–fluid interactions play an important role in mediating viscoelastic behaviour of biological tissues. In the intervertebral disc, water content is governed by a number of factors, including age, disease and mechanical loads, leading to changes in stiffness characteristics. We hypothesized that zonal stress distributions depend on load history, or the prior stresses experienced by the disc. To investigate these effects, rat caudal motion segments were subjected to compressive creep biomechanical testing in vitro using a protocol that consisted of two phases: a Prestress Phase (varied to represent different histories of load) followed immediately by an Exertion Phase, identical across all Prestress groups. Three analytical models were used to fit the experimental data in order to evaluate load history effects on gross and zonal disc mechanics. Model results indicated that while gross transient response was insensitive to load history, there may be changes in the internal mechanics of the disc. In particular, a fluid transport model suggested that the role of the nucleus pulposus in resisting creep during Exertion depended on Prestress conditions. Separate experiments using similarly defined load history regimens were performed to verify these predictions by measuring intradiscal pressure with a fibre optic sensor. We found that the ability for intradiscal pressure generation was load history-dependent and exhibited even greater sensitivity than predicted by analytical models. A 0.5?MPa Exertion load resulted in 537.2?kPa IDP for low magnitude Prestress compared with 373.7?kPa for high magnitude Prestress. Based on these measurements, we developed a simple model that may describe the pressure-shear environment in the nucleus pulposus. These findings may have important implications on our understanding of how mechanical stress contributes to disc health and disease etiology.     

Modeling the role of IGF-1 on extracellular matrix biosynthesis and cellularity in intervertebral disc

The insulin-like growth factor-1 (IGF-1) is a well-known anabolic agent for intervertebral disc (IVD), promoting both proteoglycan (PG) biosynthesis and cell proliferation. Accordingly, it is believed that IGF-1 may play a central role in IVD homeostasis. Furthermore, the exogenous administration of IGF-1 has been proposed as a possible therapeutic strategy for disc degeneration. The objectives of this study were to develop a new computational framework for describing the mechanisms regulating IGF-mediated homeostasis in IVD, and to apply this numerical tool for investigating the effectiveness of exogenous administration of IGF-1 for curing disc degeneration. A diffusive–reactive model was developed for describing competitive binding of IGF-1 to its binding proteins and cell surface receptors, with the latter reaction initiating the intracellular signaling mechanism leading to PG production and cell proliferation. Because PG production increases cell metabolic rate, and cell proliferation increases nutritional demand, nutrients transport and metabolism were also included into the model, and co-regulated, together with IGF-1, IVD cellularity. The sustainability and the effectiveness of IGF-mediated anabolism were investigated for conditions of pathologically insufficient nutrient supply, and for the case of exogenous administration of IGF-1 to degenerated IVD. Results showed that pathological nutrients deprivation, by decreasing cellularity, caused a reduction of PG biosynthesis. Also, exogenous administration of IGF-1 was only beneficial in well-nourished regions of IVD, and exacerbated cell mortality in malnourished regions. These findings remark the central role of nutrition in IVD health, and suggest that adequate nutritional supply is paramount for achieving a successful IGF-based therapy for disc degeneration.     

Thermal analysis of the human intervertebral disc

Intervertebral disc (IVD) degeneration is one of the most common musculuskeletal disorders affecting western society. Degeneration alters the morphology and the mechanical properties of the discs. According to previous reports, DSC proved to be a suitable method for the demonstration of thermal consequences of local as well as global conformational changes in the structure of the human intervertebral discs. In the present study, a wide spectrum of degenerated IVD was examined by DSC. The results suggest that definitive differences exist between the stages of disc degeneration in calorimetric measures.     

Linking continuous and discrete intervertebral disc models through homogenisation

At present, there are two main numerical approaches that are frequently used to simulate the mechanical behaviour of the human spine. Researchers with a continuum-mechanical background often utilise the finite-element method (FEM), where the involved biological soft and hard tissues are modelled on a macroscopic (continuum) level. In contrast, groups associated with the science of human movement usually apply discrete multi-body systems (MBS). Herein, the bones are modelled as rigid bodies, which are connected by Hill-type muscles and non-linear rheological spring-dashpot models to represent tendons and cartilaginous connective tissue like intervertebral discs (IVD). A possibility to benefit from both numerical methods is to couple them and use each approach, where it is most appropriate. Herein, the basic idea is to utilise MBS in simulations of the overall body and apply the FEM only to selected regions of interest. In turn, the FEM is used as homogenisation tool, which delivers more accurate non-linear relationships describing the behaviour of the IVD in the multi-body dynamics model. The goal of this contribution is to present an approach to couple both numerical methods without the necessity to apply a gluing algorithm in the context of a co-simulation. Instead, several pre-computations of the intervertebral disc are performed offline to generate an approximation of the homogenised finite-element (FE) result. In particular, the discrete degrees of freedom (DOF) of the MBS, that is, three displacements and three rotations, are applied to the FE model of the IVD, and the resulting homogenised forces and moments are recorded. Moreover, a polynomial function is presented with the discrete DOF of the MBS as variables and the discrete forces an moments as function values. For the sake of a simple verification, the coupling method is applied to a simplified motion segment of the spine. Herein, two stiff cylindrical vertebrae with an interjacent homogeneous cylindrical IVD are examined under the restriction of purely elastic deformations in the sagittal plane.     

Quantitative MRI water content mapping of porcine intervertebral disc during uniaxial compression

Background: Intervertebral disc (IVD) diseases are major public health problem in industrialized countries where they affect a large proportion of the population. In particular, IVD degeneration is considered to be one of the leading causes of pain consultation and sick leave. The aim of this study was to develop a new method for assessing the functionality of IVD in order to diagnose IVD degeneration. Methods: For this purpose, we have designed a specific device that enables to mechanically load porcine IVD ex vivo in the 4.7-Tesla horizontal superconducting magnet of a magnetic resonance (MR) scanner. Proton density weighted imaging (ρ_H-MRI) of the samples was acquired. Findings: The post-processing on MR images allowed (1) to reconstruct the 3D deformation under a known mechanical load and (2) to infer the IVD porosity assuming an incompressible poroelastic model. Interpretation: This study demonstrates the ability to follow the change in morphology and hydration of an IVD using MR measurements, thereby providing valued information for a better understanding of IVD function.     

Gene expression profile analysis of human intervertebral disc degeneration

In this study, we used microarray analysis to investigate the biogenesis and progression of intervertebral disc degeneration. The gene expression profiles of 37 disc tissue samples obtained from patients with herniated discs and degenerative disc disease collected by the National Cancer Institute Cooperative Tissue Network were analyzed. Differentially expressed genes between more and less degenerated discs were identified by significant analysis of microarray. A total of 555 genes were significantly overexpressed in more degenerated discs with a false discovery rate of < 3%. Functional annotation showed that these genes were significantly associated with membrane-bound vesicles, calcium ion binding and extracellular matrix. Protein-protein interaction analysis showed that these genes, including previously reported genes such as fibronectin, COL2A1 and β-catenin, may play key roles in disc degeneration. Unsupervised clustering indicated that the widely used morphology-based Thompson grading system was only marginally associated with the molecular classification of intervertebral disc degeneration. These findings indicate that detailed, systematic gene analysis may be a useful way of studying the biology of intervertebral disc degeneration.     

Tissue-engineered intervertebral disc and chondrogenesis using human nucleus pulposus regulated through TGF-beta1 in platelet-rich plasma

Human intervertebral disc (IVD) degeneration often initiated from the human nucleus pulposus (hNP) with aging leading to IVD destruction and extracellular matrix (ECM) depletion. Previously, we have successfully employed transforming growth factor-beta1 (TGF-beta1) to promote chondrogenesis of mesenchymal progenitor cells (MPCs) and immortalized human mesenchymal stem cells. In this study, we examine the role of TGF-beta1 in platelet-rich plasma (PRP) on disc regeneration, including proliferation, redifferentiation, and the reconstitution of tissue-engineered NP. hNP cells were isolated from volunteers with different ages and cultured in the presence of PRP. We found that the most effective concentration for hNP proliferation was 1 ng/ml TGF-beta1 in PRP, which was further applied in the following experiments. hNP cell proliferation in all age groups were increased time-dependently by PRP and cell morphologies showed aggregation. The mRNA of Sox9, type II collagen, and aggrecan were all significantly upregulated by PRP through RT-PCR. Glycosaminoglycan (GAG) accumulation reached the highest value at day 7 and continued to day 9 culture. PRP promoted NP regeneration via the Smad pathway was also determined and highly activated p-Smad2/3 at 30 min and continuously sustained to 120 min. Immunostaining of type II collagen indicates that PRP participates in chondrogenesis of tissue-engineered NP with collagen scaffolds. We concluded that growth factors in PRP can effectively react as a growth factor cocktail to induce hNP proliferation and differentiation, and also promote tissue-engineered NP formation. These findings are the first to demonstrate that PRP might be a therapeutic candidate for prevention of disc degeneration.     

Diurnal variations in intervertebral disc height affect spine flexibility,intradiscal pressure and contact compressive forces in the facet joints

Diurnal changes of intervertebral disc height are caused by high compressive loading during the day, which expulses fluid from the disc, and by osmotic pressure, which imbibes fluid into the disc at low loading. The aim of the present study was to determine the magnitude of diurnal changes in spine flexibility, intradiscal pressures and contact forces in the facet joints. A validated osseoligamentous finite element model of the lumbar spine was used to determine these quantities for morning and evening situations. Disc height varied by 10% for these two situations. Spine flexibility and facet joint forces were markedly higher in the evening than in the morning. Intradiscal pressures were higher in the morning than in the evening. The different spine flexibilities in the morning and evening should be taken into account during kinematical measurements. Predicted facet joint forces may be used for the designing and pre-clinical testing of artificial facet joint replacements.     

Degeneration of the intervertebral disc

The intervertebral disc is a cartilaginous structure that resembles articular cartilage in its biochemistry, but morphologically it is clearly different. It shows degenerative and ageing changes earlier than does any other connective tissue in the body. It is believed to be important clinically because there is an association of disc degeneration with back pain. Current treatments are predominantly conservative or, less commonly, surgical; in many cases there is no clear diagnosis and therapy is considered inadequate. New developments, such as genetic and biological approaches, may allow better diagnosis and treatments in the future.     

Intradiscal injection of simvastatin retards progression of intervertebral disc degeneration induced by stab injury

Introduction  

Earlier work indicates that the cholesterol-lowering drug, simvastatin, is anabolic to chondrogenic expression of rat intervertebral disc (IVD) cells, which suggests a potential role for simvastatin in IVD regeneration. In this study, we expand on our earlier work to test the effectiveness of simvastatin on disc degeneration utilizing a rat tail disc degeneration model.     

Diffusion of antibiotics in intervertebral disc

Delivering charged antibiotics to the intervertebral disc is challenging because of the avascular, negatively charged extracellular matrix (ECM) of the tissue. The purpose of this study was to measure the apparent diffusion coefficient of two clinically relevant, charged antibiotics, vancomycin (positively charged) and oxacillin (negatively charged) in IVD. A one-dimensional steady state diffusion experiment was employed to measure the apparent diffusion coefficient of the two antibiotics in bovine coccygeal annulus fibrosus (AF) tissue. The averaged apparent diffusion coefficient for vancomycin under 20% compressive strain was 7.94 ± 2.00 × 10⁻¹² m²/s (n = 10), while that of oxacillin was 2.26 ± 0.68 × 10⁻¹⁰ m²/s (n = 10). A student’s t-test showed that the diffusivity of vancomycin was significantly lower than that of oxacillin. This finding may be attributed to two factors: solute size and possible binding effects. Vancomycin is approximately 3 times larger in molecular weight than oxacillin, meaning that steric hindrance likely plays a role in the slower transport. Reversible binding between positive vancomycin and the negative ECM could also slow down the rate of diffusion. Therefore, more investigation is necessary to determine the specific relationship between net charge on antibiotic and diffusion coefficients in IVD. This study provides essential quantitative information regarding the transport rates of antibiotics in the IVD, which is critical in using computational modeling to design effective strategies to treat disc infection.     

Pathophysiology of the human intervertebral disc

Intervertebral disc degeneration is a common invalidating disorder that can affect the musculoskeletal apparatus in both younger and older ages. The chief component of the intervertebral disc is the highly organized extracellular matrix; maintenance of its organization is essential for correct spinal mechanics. The matrix components, mainly proteoglycans and collagens, undergo a slow and continuous cell-mediated turnover process that enables disc cells to adapt their environment to external stimuli. Cellular senescence and a history of chronic abnormal loading can upset this balance, leading to progressive tissue failure that results in disc degeneration. Although biological treatment approaches to disc repair are still far to come, advances in our understanding of disc biochemistry and in defining the role of genetic inheritance have provided a starting point for developing new concepts in the diagnosis, therapy and prevention of disc degeneration.     

The influence of exogenous cross-linking and compressive creep loading on intradiscal pressure

This study involves a biomechanical evaluation of a prospective injectable treatment for degenerative discs. The high osmolarity of the non-degenerated nucleus pulposus attracts water contributing to the hydrostatic behavior of the tissue. This intradiscal pressure is known to drop as fluid is exuded from the matrix due to compressive loading. The objective of this study was to compare the changes in intradiscal pressure in control and genipin cross-linked intervertebral discs. Thirty bovine lumbar motion segments were randomly divided into a phosphate-buffered saline control group and a 0.33% genipin group and soaked at room temperature for 2 days. A needle pressure sensor was held in the center of the disc while short-term and static creep compressive loads were applied. The control group demonstrated a 25% higher average intradiscal pressure compared to genipin-treated discs under 750 N compressive load (p=0.029). Depressurization during static compressive creep was 56% higher in the control than in the genipin group (p=0.014). These results suggest cross-linking induced changes in the poroelastic properties of the involved tissues affected the mechanics of compressive load support in the disc with lower levels of nucleus pressure, a corresponding decrease in the elastic expansion of the annulus, and an increased axial compressive loading of the inner and outer annulus tissues. It is possible that concurrent changes in hydraulic permeability and proteoglycan retention known to be associated with genipin cross-linking were also contributors to poroelastic changes. Reduction of peak pressures and moderation of pressure fluctuations could be beneficial relative to discogenic pain.     

The biological basis of degenerative disc disease: proteomic and biomechanical analysis of the canine intervertebral disc

IntroductionIn the present study, we sought to quantify and contrast the secretome and biomechanical properties of the non-chondrodystrophic (NCD) and chondrodystrophic (CD) canine intervertebral disc (IVD) nucleus pulposus (NP).MethodsWe used iTRAQ proteomic methods to quantify the secretome of both CD and NCD NP. Differential levels of proteins detected were further verified using immunohistochemistry, Western blotting, and proteoglycan extraction in order to evaluate the integrity of the small leucine-rich proteoglycans (SLRPs) decorin and biglycan. Additionally, we used robotic biomechanical testing to evaluate the biomechanical properties of spinal motion segments from both CD and NCD canines.ResultsWe detected differential levels of decorin, biglycan, and fibronectin, as well as of other important extracellular matrix (ECM)-related proteins, such as fibromodulin and HAPLN1 in the IVD NP obtained from CD canines compared with NCD canines. The core proteins of the vital SLRPs decorin and biglycan were fragmented in CD NP but were intact in the NP of the NCD animals. CD and NCD vertebral motion segments demonstrated significant differences, with the CD segments having less stiffness and a more varied range of motion.ConclusionsThe CD NP recapitulates key elements of human degenerative disc disease. Our data suggest that at least some of the compromised biomechanical properties of the degenerative disc arise from fibrocartilaginous metaplasia of the NP secondary to fragmentation of SLRP core proteins and associated degenerative changes affecting the ECM. This study demonstrates that the degenerative changes that naturally occur within the CD NP make this animal a valuable animal model with which to study IVD degeneration and potential biological therapeutics.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0733-z) contains supplementary material, which is available to authorized users.     

Morphometric analysis of vertebrae and intervertebral discs as a basis of disc replacement

In this study measurements were obtained from 359 vertebrae and 215 intervertebral discs in an attempt to classify discs by their size. At the cervical and thoracic levels, this attempt was unproductive because of extensive variations. In the lumbar spine, discs were allocated to six size-matched groups and to two height-matched groups. The breadth of marginal rims were also measured with a view to provide surgeons operating on the spine with precise data on disc morphology to facilitate disc replacement.