The internal mechanics of the intervertebral disc under cyclic loading

The mechanics of the intervertebral disc (IVD) under cyclic loading are investigated via a one-dimensional poroelastic model and experiment. The poroelastic model, based on that of Biot (J. Appl. Phys. 12 (1941) 155; J. Appl. Mech. 23 (1956) 91), includes a power-law relation between porosity and permeability, and a linear relation between the osmotic potential and solidity. The model was fitted to experimental data of the unconfined IVD undergoing 5 cyclic loads of 20 min compression by an applied stress of 1MPa, followed by 40 min expansion. To obtain a good agreement between experiment and theory, the initial elastic deformation of the IVD, possibly associated with the bulging of the IVD into the vertebral bodies or laterally, was removed from the experimental data. Many combinations of the permeability-porosity relationship with the initial osmotic potential (pi(i)) were investigated, and the best-fit parameters for the aggregate modulus (H(A)) and initial permeability (k(i)) were determined. The values of H(A) and k(i) were compared to literature values, and agreed well especially in the context of the adopted high-stress testing regime, and the strain related permeability in the model.     

Novel ex-vivo mechanobiological intervertebral disc culture system

Intervertebral disc degeneration, a leading cause of low back pain, poses a significant socioeconomic burden with a broad array of costly treatment options. Mechanical loading is important in disease progression and treatment. Connecting mechanics and biology is critical for determining how loading parameters affect cellular response and matrix homeostasis. A novel ex-vivo experimental platform was developed to facilitate in-situ loading of rabbit functional spinal units (FSUs) with relevant biological outcome measures. The system was designed for motion outside of an incubator and validated for rigid fixation and physiologic environmental conditions. Specimen motion relative to novel fixtures was assessed using a digitizer; fixture stiffness exceeded specimen stiffness by an order of magnitude. Intradiscal pressure (IDP), measured using a fiber-optic pressure transducer, confirmed rigidity and compressive force selection. Surrounding media was controlled at 37 °C, 5% O(2)/CO(2) using a closed flow loop with an hypoxic incubator and was validated with probes in the specimen chamber. FSUs were subjected to cyclic compression (20 cycles) and four-hour creep at 1.0 MPa. Disc tissue was analyzed for cell viability using 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), which showed high viability (>90%) regardless of loading. Conditioned media was assayed for type-II collagen degradation fragments (CTX-II) and an aggrecan epitope (CS-846) associated with new aggrecan synthesis. CTX-II concentrations were not associated with loading, but CS-846 concentrations appeared to be increased with loading. Preservation of the full FSU allows physiologic load transmission and future multi-axis motion and identification of load-responsive proteins, thereby forming a new niche in intervertebral disc organ culture.     

Reduced nucleus pulposus glycosaminoglycan content alters intervertebral disc dynamic viscoelastic mechanics

The intervertebral disc functions over a range of dynamic loading regimes including axial loads applied across a spectrum of frequencies at varying compressive loads. Biochemical changes occurring in early degeneration, including reduced nucleus pulposus glycosaminoglycan content, may alter disc mechanical behavior and thus may contribute to the progression of degeneration. The objective of this study was to determine disc dynamic viscoelastic properties under several equilibrium loads and loading frequencies, and further, to determine how reduced nucleus glycosaminoglycan content alters dynamic mechanics. We hypothesized that (1) dynamic stiffness would be elevated with increasing equilibrium load and increasing frequency, (2) the disc would behave more elastically at higher frequencies, and finally, (3) dynamic stiffness would be reduced at low equilibrium loads under all frequencies due to nucleus glycosaminoglycan loss. We mechanically tested control and chondroitinase ABC injected rat lumbar motion segments at several equilibrium loads using oscillatory loading at frequencies ranging from 0.05 to 5 Hz. The rat lumbar disc behaved non-linearly with higher dynamic stiffness at elevated compressive loads irrespective of frequency. Phase angle was not affected by equilibrium load, although it decreased as frequency was increased. Reduced glycosaminoglycan decreased dynamic stiffness at low loads but not at high equilibrium loads and led to increased phase angle at all loads and frequencies. The findings of this study demonstrate the effect of equilibrium load and loading frequencies on dynamic disc mechanics and indicate possible mechanical mechanisms through which disc degeneration can progress.     

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.     

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.     

Effect of annular incision type on the change in biomechanical properties in a herniated lumbar intervertebral disc

The technique used to incise the disc during discectomy may play a role in the subsequent healing and change in biomechanical stiffness of the disc. Several techniques of lumbar disc annulotomy have been described in clinical reports. The purpose of this paper was to study the influence of annulotomy technique on motion segment stiffness using a finite element model. Four incision methods (square, circular, cross, and slit) were compared. The analyses showed that each of the annular incisions produced increase in motions under axial moment loadings with circular incision producing the largest change in the corresponding rotational motion. Under shear loading mode, cross and slit-type annular incisions produced slightly larger changes in the principal motions of the disc than square and circular incisions. All other incision types considered in the current study produced negligibly small increase in motion under rest of the loading conditions. In addition to annulotomy, when nucleotomy was also included in the analyses, once again cross and slit incisions produced larger change in motion under shear loading mode as compared to the other two incision types. A comparison between the four types of annular incisions showed that cross incision produced an increase in motion larger than those produced by the other three incisions under flexion/extension and lateral moment loading and both shear force loadings. Circular incision produced the largest increase in motion under axial moment load in comparison to those produced by square, cross, and slit incisions. Sagittal plane symmetry was influenced by the incision injury to the motion segment leading to coupled motions as well as increased facet loads. From the study it can be concluded that the increase inflexibility of the disc due to annulotomy depends on the type of annulotomy and the annulotomy also produce asymmetrical deformations leading to increased facet loading.     

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.     

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.     

Neutral proteinases of the human intervertebral disc

Disc tissue consisting of pooled annuli fibrosus and nuclei pulposus from the cadaver of an adolescent aged 19 years was extracted with 4.0 M Gu-HCl. Proteins of low buoyant density (p less than or equal to 1.38 g/ml) containing the disc enzymes and inhibitors were separated from proteoglycans of high buoyant density (p greater than or equal to 1.50 g/ml) by density gradient ultracentrifugation. Sephadex G-75F gel chromatography followed by trypsin affinity chromatography was then used to resolve disc proteolytic and trypsin inhibitory activities. The results obtained were strongly suggestive of the presence of a high molecular weight zymogen which upon activation generated a population of smaller molecular weight proteinases. The disc proteinases obtained by this process showed similar properties in terms of: their pH optima (7.4-7.6); their inhibition patterns by class-specific proteinase inhibitors; their variation of activity as a function of NaCl and lysine concentrations; and the hydrodynamic size of their proteoglycan degradation products. The activated disc neutral proteinase demonstrated many characteristics in common with plasmin; however, unlike the latter, the disc proteinases also showed some calcium dependence.     

The proteoglycans of the canine intervertebral disc

The high-buoyant-density proteoglycans of the nucleus pulposus and annulus fibrosus of the beagle intervertebral disc have been isolated by CsCl density gradient ultracentrifugation. The sulphated proteoglycans were labelled in vivo with 35SO4, 24 h and 60 days prior to killing. The hydrodynamic size and aggregation of the 24 h, 60 day and resident (from hexuronic acid and hexosamine analysis) proteoglycan subunit populations were determined by Sepharose CL-2B chromatography in the presence or absence of excess hyaluronic acid. The hydrodynamic size of the keratan sulphate-proteoglycan core protein complexes were also determined by Sepharose CL-2B chromatography after chondroitinase ABC digestion of proteoglycans. When initially synthesised (24 h) or after 60 days, the percentage aggregation and hydrodynamic size of the proteoglycans derived from the annulus fibrosus were larger than those present in the nucleus pulposus. Hexosamine, hexuronic and protein determination of the high-buoyant-density fractions showed that the proteoglycans of the nucleus pulposus were richer in chondroitin sulphate than those in the annulus. However there was no difference in Mr of the chondroitin sulphate and keratan sulphate attached to the proteoglycans of the two disc regions, nor were differences detected by HPLC between the proportions of chondroitin 4-sulphate and chondroitin 6-sulphate present in these high-density fractions. In contrast, the low-buoyant-density (1.54 greater than p greater than 1.45) proteoglycan fractions and tissue residues remaining after 4 M GuHCl extraction were found to contain dermatan sulphate, suggesting the presence of a third proteoglycan species possibly associated with the collagen of the fibrocartilagenous matrix.     

Proteoglycans of the intervertebral disc. Absence of degradation during the isolation of proteoglycans from the intervertebral disc.

Proteoglycans extracted with 4M-guanidinium chloride from pig intervetebral discs, and purified by equilibrium density-gradient centrifugation in CsCl, were of smaller hydrodynamic size than those extracted and purified in the same way from the laryngeal cartilage of the same animal. Whether this difference in size arose from degradation during the extraction and purification of the proteoglycans of the disc was investigated. Purified proteoglycans labelled either in the chondroitin sulphate chains or in the core protein were obtained from laryngeal cartilage by short-term organ culture. These labelled proteoglycans were added at the beginning of the extraction of the disc proteoglycans, and labelled cartilage and unlabelled disc proteoglycans were isolated and purified together. There was no appreciable loss of radioactivity after density-gradient centrifugation nor decrease in hydrodynamic size of the labelled cartilage proteoglycans on chromatography on Sepharose 2B, when these were present during the extraction of disc proteoglycans. It is concluded that disc proteoglycans are intrinsically of smaller size than cartilage proteoglycans and this difference in size does not arise from degradation during the extraction.     

Aquaporin expression in the human intervertebral disc

The nucleus pulposus (NP) of the human intervertebral disc (IVD) is a hyperosmotic tissue that is subjected to daily dynamic compressive loads. In order to survive within this environment the resident chondrocyte-like cells must be able to control their cell volume, whilst also controlling the anabolism and catabolism of their extra-cellular matrix. Recent studies have demonstrated expression of a range of bi-directional, transmembrane water and solute transporters, named aquaporins (AQPs), within chondrocytes of articular cartilage. The aim of this study was to use immunohistochemsitry to investigate the expression of aquaporins 1, 2 and 3 within the human IVD. Results demonstrated expression of both AQP-1 and -3 by cells within the NP and inner annulus fibrosus (AF), while outer AF cells lacked expression of AQP-1 and showed very low numbers of AQP-3 immunopositive cells. Cells from all regions were negative for AQP-2. Therefore this study demonstrates similarities in the phenotype of NP cells and articular chondrocytes, which may be due to similarities in tissue osmolarity and mechanobiology. The decrease in expression of AQPs from the NP to the outer AF may signify changes in cellular phenotype in response to differences in mechanbiology, osmolarity and hydration between the gelatinous NP and the fibrous AF.     

Substructuring and poroelastic modelling of the intervertebral disc

We proposed a substructure technique to predict the time-dependant response of biological tissue within the framework of a finite element resolution. Theoretical considerations in poroelasticity preceded the calculation of the sub-structured poroelastic matrix. The transient response was obtained using an exponential fitting method. We computed the creep response of an MRI 3D reconstructed L₅–S₁ intervertebral disc of a scoliotic spine. The FE model was reduced from 10,000 degrees of freedom for the full 3D disc to only 40 degrees of freedom for the sub-structured model defined by 10 nodes attached to junction nodes located on both lower and upper surfaces of the disc. Comparisons of displacement fields were made between the full poroelastic FE model and the sub-structured model in three different loading conditions: compression, offset compression and torsion. Discrepancies in displacement were lower than 10% for the first time steps when time-dependant events were significant. The substructuring technique provided an exact solution in quasi-static behavior after pressure relaxation. Couplings between vertical and transversal displacements predicted by the reference FE model were well stored by the sub-structured model despite the drastic reduction of degrees of freedom. Finally, we demonstrated that substructuring was very efficient to reduce the size of numerical models while respecting the time-dependant behavior of the structure. This result highlighted the potential interest of substructure techniques in large-scale models of musculoskeletal structures.