3D finite element analysis of nutrient distributions and cell viability in the intervertebral disc: effects of deformation and degeneration

The intervertebral disc (IVD) receives important nutrients, such as glucose, from surrounding blood vessels. Poor nutritional supply is believed to play a key role in disc degeneration. Several investigators have presented finite element models of the IVD to investigate disc nutrition; however, none has predicted nutrient levels and cell viability in the disc with a realistic 3D geometry and tissue properties coupled to mechanical deformation. Understanding how degeneration and loading affect nutrition and cell viability is necessary for elucidating the mechanisms of disc degeneration and low back pain. The objective of this study was to analyze the effects of disc degeneration and static deformation on glucose distributions and cell viability in the IVD using finite element analysis. A realistic 3D finite element model of the IVD was developed based on mechano-electrochemical mixture theory. In the model, the cellular metabolic activities and viability were related to nutrient concentrations, and transport properties of nutrients were dependent on tissue deformation. The effects of disc degeneration and mechanical compression on glucose concentrations and cell density distributions in the IVD were investigated. To examine effects of disc degeneration, tissue properties were altered to reflect those of degenerated tissue, including reduced water content, fixed charge density, height, and endplate permeability. Two mechanical loading conditions were also investigated: a reference (undeformed) case and a 10% static deformation case. In general, nutrient levels decreased moving away from the nutritional supply at the disc periphery. Minimum glucose levels were at the interface between the nucleus and annulus regions of the disc. Deformation caused a 6.2% decrease in the minimum glucose concentration in the normal IVD, while degeneration resulted in an 80% decrease. Although cell density was not affected in the undeformed normal disc, there was a decrease in cell viability in the degenerated case, in which averaged cell density fell 11% compared with the normal case. This effect was further exacerbated by deformation of the degenerated IVD. Both deformation and disc degeneration altered the glucose distribution in the IVD. For the degenerated case, glucose levels fell below levels necessary for maintaining cell viability, and cell density decreased. This study provides important insight into nutrition-related mechanisms of disc degeneration. Moreover, our model may serve as a powerful tool in the development of new treatments for low back pain.     

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.     

Analysis of cell viability in intervertebral disc: Effect of endplate permeability on cell population

Responsible for making and maintaining the extracellular matrix, the cells of intervertebral discs are supplied with essential nutrients by diffusion from the blood supply through mainly the cartilaginous endplates (CEPs) and disc tissue. Decrease in transport rate and increase in cellular activity may adversely disturb the intricate supply–demand balance leading ultimately to cell death and disc degeneration. The present numerical study aimed to introduce for the first time cell viability criteria into nonlinear coupled nutrition transport equations thereby evaluating the dynamic nutritional processes governing viable cell population and concentrations of oxygen, glucose and lactic acid in the disc as CEP exchange area dropped from a fully permeable condition to an almost impermeable one. A uniaxial model of an in vitro cell culture analogue of the disc is first employed to examine and validate cell viability criteria. An axisymmetric model of the disc with four distinct regions was subsequently used to investigate the survival of cells at different CEP exchange areas.In agreement with measurements, predictions of the diffusion chamber model demonstrated substantial cell death as essential nutrient concentrations fell to levels too low to support cells. Cells died away from the nutrient supply and at higher cell densities. In the disc model, the nucleus region being farthest away from supply sources was most affected; cell death initiated first as CEP exchange area dropped below ～40% and continued exponentially thereafter to depletion as CEP calcified further. In cases with loss of endplate permeability and/or disruptions therein, as well as changes in geometry and fall in diffusivity associated with fluid outflow, the nutrient concentrations could fall to levels inadequate to maintain cellular activity or viability, resulting in cell death and disc degeneration.     

Analysis of the influence of disc degeneration on the mechanical behaviour of a lumbar motion segment using the finite element method

Compared to a healthy intervertebral disc, the geometry and the material properties of the involved tissues are altered in a degenerated disc. It is not completely understood how this affects the mechanical behaviour of a motion segment. In order to study the influence of disc degeneration on motion segment mechanics a three-dimensional, nonlinear finite element model of the L3/L4 functional unit was used. Different grades of disc degeneration were simulated by varying disc height and bulk modulus of the nucleus pulposus. The model was loaded with pure moments of 10Nm in the three main anatomic planes. The finite element model predicted the same trends for intersegmental rotation and intradiscal pressure as described in the literature for in vitro studies. A comparison between calculated intersegmental rotation and experimental data showed a mean difference of 1.9 degrees while the mean standard deviation was 2.5 degrees . A mildly degenerated disc increases intersegmental rotation for all loading cases. With further increasing disc degeneration intersegmental rotation is decreased. For axial rotation the decrease takes place in the final stage. Intradiscal pressure is lower while facet joint force and maximum von Mises stress in the annulus are higher in a degenerated compared to a healthy disc.     

Effects of Tobacco Smoking on the Degeneration of the Intervertebral Disc: A Finite Element Study

Tobacco smoking is associated with numerous pathological conditions. Compelling experimental evidence associates smoking to the degeneration of the intervertebral disc (IVD). In particular, it has been shown that nicotine down-regulates both the proliferation rate and glycosaminoglycan (GAG) biosynthesis of disc cells. Moreover, tobacco smoking causes the constriction of the vascular network surrounding the IVD, thus reducing the exchange of nutrients and anabolic agents from the blood vessels to the disc. It has been hypothesized that both nicotine presence in the IVD and the reduced solute exchange are responsible for the degeneration of the disc due to tobacco smoking, but their effects on tissue homeostasis have never been quantified. In this study, a previously presented computational model describing the homeostasis of the IVD was deployed to investigate the effects of impaired solute supply and nicotine-mediated down-regulation of cell proliferation and biosynthetic activity on the health of the disc. We found that the nicotine-mediated down-regulation of cell anabolism mostly affected the GAG concentration at the cartilage endplate, reducing it up to 65% of the value attained in normal physiological conditions. In contrast, the reduction of solutes exchange between blood vessels and disc tissue mostly affected the nucleus pulposus, whose cell density and GAG levels were reduced up to 50% of their normal physiological levels. The effectiveness of quitting smoking on the regeneration of a degenerated IVD was also investigated, and showed to have limited benefit on the health of the disc. A cell-based therapy in conjunction with smoke cessation provided significant improvements in disc health, suggesting that, besides quitting smoking, additional treatments should be implemented in the attempt to recover the health of an IVD degenerated by tobacco smoking.     

Matrix-assisted cell transfer for intervertebral disc cell therapy

Cell therapy seems to be a promising way to reconstitute degenerated discs. We elucidate the basic aspects of intervertebral disc (IVD) cell therapy to estimate its potential in disc regeneration. Cell transfer efficiency and survival was quantified by luciferase expression after injection of recombinant cells into healthy, nucleotomized or mechanically degenerated rabbit IVDs in vitro, in situ or in vivo. A two-component fibrin matrix was adapted to allow injection of a fluid cell suspension that quickly polymerizes in IVDs. Thirty-five to fifty percent of matrix injected cells remained in the nucleus and transition zone in contrast to a rapid loss of medium-injected cells. Nucleotomy, which reduces intradiscal pressure, was crucial to the survival of the transferred cells over 3 days and nutritional enrichment of the fibrin matrix with potent biomolecules from serum significantly enhanced cell viability. In conclusion, advanced matrix substitutes are needed for efficient transfer and improved cell survival in the low-nutrient intradiscal environment to further improve disc cell therapy.     

Influence of fixed charge density magnitude and distribution on the intervertebral disc: applications of a poroelastic and chemical electric (PEACE) model

A 3-dimensional formulation for a poroelastic and chemical electric (PEACE) model is presented and applied to an intervertebral disc slice in a 1-dimensional validation problem and a 2-dimensional plane stress problem. The model was used to investigate the influence of fixed charge density magnitude and distribution on this slice of disc material. Results indicated that the mechanical, chemical, and electrical behaviors were all strongly influenced by the amount as well as the distribution of fixed charges in the matrix. Without any other changes in material properties, alterations in the fixed charge density (proteoglycan content) from a healthy to a degenerated distribution will cause an increase in solid matrix stresses and can affect whether the tissue imbibes or exudes fluid under different loading conditions. Disc tissue with a degenerated fixed charge density distribution exhibited greater solid matrix stresses and decreased streaming potential, all of which have implications for disc nutrition, disc biomechanics, and tissue remodeling. It was also seen that application of an electrical potential across the disc can induce fluid transport.     

Effect of Cartilage Endplate on Cell Based Disc Regeneration: A Finite Element Analysis

This study examines the effects of cartilage endplate (CEP) calcification and the injection of intervertebral disc (IVD) cells on the nutrition distributions inside the human IVD under physiological loading conditions using multiphasic finite element modeling. The human disc was modeled as an inhomogeneous mixture consisting of a charged elastic solid, water, ions (Na⁺ and Cl⁻), and nutrient solute(oxygen,glucose and lactate) phases. The effect of the endplate calcification was simulated by a reduction of the tissue porosity (i.e., water volume faction) from 0.60 to 0.48. The effect of cell injection was simulated by increasing the cell density in the nucleus pulposus (NP) region by 50%, 100%, and 150%. Strain-dependent transport properties(e.g., hydraulic permeability and solute diffusivities) were considered to couple the solute transport and the mechanical loading. The simulation results showed that nutrient solute distribution inside the discis maintained at a stable state during the day and night. The physiological diurnal cyclic loading does not change the nutrient environment in the human IVD. The cartilage endplate plays a significant role in the nutrient supply to human IVD. Calcification of the cartilage endplate significantly reduces the nutrient levels in human IVD. Therefore, in cell based therapy for IVD regeneration, theincreased nutrient demand as a result of cell injection needs to be addressed. Excessive numbers of injected cells may cause further deterioration of the nutrient environment in the degenerated disc. This study is important for understanding the pathology of IVD degeneration and providing new insights into cell based therapies for low back pain.     

Intercellular communication via gap junctions affected by mechanical load in the bovine annulus fibrosus

Cells in the intervertebral disc, as in other connective tissues including tendon, ligament and bone, form interconnected cellular networks that are linked via functional gap junctions. These cellular networks may be necessary to affect a coordinated response to mechanical and environmental stimuli. Using confocal microscopy with fluorescence recovery after photobleaching methods, we explored the in situ strain environment of the outer annulus of an intact bovine disc and the effect of high-level flexion on gap junction signalling. The in situ strain environment in the extracellular matrix of the outer annulus under high flexion load was observed to be non-uniform with the extensive cellular processes remaining crimped sometimes at flexion angles greater than 25°. A significant transient disruption of intercellular communication via functional gap junctions was measured after 10 and 20 min under high flexion load. This study illustrates that in healthy annulus fibrosus tissue, high mechanical loads can impede the functioning of the gap junctions. Future studies will explore more complex loading conditions to determine whether losses in intercellular communication can be permanent and whether gap junctions in aged and degenerated tissues become more susceptible to load. The current research suggests that cellular structures such as gap junctions and intercellular networks, as well as other cell–cell and cell–matrix interconnections, need to be considered in computational models in order to fully understand how macroscale mechanical signals are transmitted across scales to the microscale and ultimately into a cellular biosynthetic response in collagenous tissues.     

Effect of intervertebral disc degeneration on disc cell viability: a numerical investigation

Degeneration of the intervertebral disc may be initiated and supported by impairment of the nutrition processes of the disc cells. The effects of degenerative changes on cell nutrition are, however, only partially understood. In this work, a finite volume model was used to investigate the effect of endplate calcification, water loss, reduction of disc height and cyclic mechanical loading on the sustainability of the disc cell population. Oxygen, lactate and glucose diffusion, production and consumption were modelled with non-linear coupled partial differential equations. Oxygen and glucose consumption and lactate production were expressed as a function of local oxygen concentration, pH and cell density. The cell viability criteria were based on local glucose concentration and pH. Considering a disc with normal water content, cell death was initiated in the centre of the nucleus for oxygen, glucose, and lactate diffusivities in the cartilaginous endplate below 20% of the physiological values. The initial cell population could not be sustained even in the non-calcified endplates when a reduction of diffusion inside the disc due to water loss was modelled. Alterations in the disc shape such as height loss, which shortens the transport route between the nutrient sources and the cells, and cyclic mechanical loads, could enhance cell nutrition processes.     

Intervertebral disc viscoelastic parameters and residual mechanics spatially quantified using a hybrid confined/in situ indentation method

With advancing age, injury, musculoskeletal pathology or a combination of these, a degenerative cascade of biomechanical, biochemical, and nutritional alterations diminish the intervertebral discs' ability to maintain its structure and function. While the biomechanics of isolated disc tissues has been investigated across this degenerative spectrum, none have attempted to retain the in situ disc-endplate morphology during compressive tissue characterization. The objective of this study was to spatially quantify the viscoelastic parameters of the intervertebral disc throughout degeneration, including the as yet unreported residual stress/strain. This required the development of a hybrid confined/in situ indentation methodology, which preserves the disc structural morphology. At four locations of the disc (anterior-AF, right and left lateral AF, and NP) stress-relaxation tests were performed using the hybrid confined/in situ indentation method, which utilizes the vertebral endplate as the porous indenter tip. This method allows the endplate to remain interwoven with the disc tissue, retaining its native orientation. Healthy disc tissue exhibited significantly higher residual stress values compared to both moderate and severe degeneration in all locations (p<0.0156). Furthermore, the equilibrium stress at 15% strain (stress relaxation) was significantly diminished with advancing disc degeneration (p<0.0241). The equilibrium viscoelastic parameters show healthy discs encounter higher forces at the same strain level, and are able to maintain this force, where degenerated discs are unable to maintain this force throughout time. This morphology-conserved method provides insight into the spatial compressive mechanical properties of the intervertebral disc across the degeneration spectrum and will aid in modeling these tissue changes.     

Simulated influence of osteoporosis and disc degeneration on the load transfer in a lumbar functional spinal unit

As life expectancy increases, age-related disorders and the search for related medical care will expand. Osteoporosis is the most frequent skeletal disease in this context with the highest fracture risk existing for vertebrae. The aging process is accompanied by systemic changes, with the earliest degeneration occurring in the intervertebral discs. The influence of various degrees of disc degeneration on the load transfer was examined using the finite element method. The effect of different possible alterations of the bone quality due to osteoporosis was simulated by adjusting the corresponding material properties and their distribution and several loadings were applied. An alteration of the load transfer, characterised by changed compression stiffness and strain distributions as well as magnitudes, due to osteoporotic bone and degenerated discs was found. When osteoporosis was simulated, the stiffness was substantially decreased, larger areas of the cancellous bone were subjected to higher strains and strain maxima were increased. Increasing ratios of transverse isotropy in the osteoporotic bone yielded smaller effects than reduced bone properties. Including a degenerated disc mainly altered the strain distribution. Combining osteoporosis and degenerated discs reduced the areas of cancellous bone subjected to substantial strain. Based on these results, it can be concluded that the definition of a healthy disc in osteoporotic spines might be considered as a worst-case scenario. One attempt to evaluate the progress of osteoporosis can be made by introducing increasing degrees of anisotropy. If several parameters in a model are changed to simulate degeneration, it should be pointed out how each individual definition influences the overall result.     

The micromechanical environment of intervertebral disc cells: effect of matrix anisotropy and cell geometry predicted by a linear model

Cells of the intervertebral disc exhibit spatial variations in phenotype and morphology that may be related to differences in their local mechanical environments. In this study, the stresses, strains, and dilatations in and around cells of the intervertebral disc were studied with an analytical model of the cell as a mechanical inclusion embedded in a transversely isotropic matrix. In response to tensile loading of the matrix, the local mechanical environment of the cell differed among the anatomic regions of the disc and was strongly influenced by changes in both matrix anisotropy and parameters of cell geometry. The results of this study suggest that the local cellular mechanical environment may play a role in determining both cell morphology in situ and the inhomogeneous response to mechanical loading observed in cells of the disc.     

Nutrient distribution and metabolism in the intervertebral disc in the unloaded state: A parametric study

A 2-D finite element model for the intervertebral disc in which quadriphasic theory is coupled to the transport of solutes involved in cellular nutrition was developed for investigating the main factors contributing to disc degeneration. Degeneration is generally considered to result from chronic disc cell nutrition insufficiency, which prevents the cells from renewing the extracellular matrix and thus leads to the loss of proteoglycans. Hence, the osmotic power of the disc is decreased, causing osmomechanical impairments. Cellular metabolism depends strongly on the oxygen, lactate and glucose concentrations and on pH in the disc. To study the diffusion of these solutes in a mechanically or osmotically loaded disc, the osmomechanical and diffusive effects have to be coupled. The intervertebral disc is modeled here using a plane strain formulation at the equilibrium state under physiological conditions after a long rest period (called unloaded state). The correlations between solute distribution and various properties of healthy and degenerated discs are investigated. The numerical simulation shows that solute distribution in the disc depends very little on the elastic modulus or the proteoglycan concentration but greatly on the porosity, diffusion coefficient and endplate diffusion area. This coupled model therefore opens new perspectives for investigating intervertebral disc degeneration mechanisms.     

The effect of sustained compression on oxygen metabolic transport in the intervertebral disc decreases with degenerative changes

Intervertebral disc metabolic transport is essential to the functional spine and provides the cells with the nutrients necessary to tissue maintenance. Disc degenerative changes alter the tissue mechanics, but interactions between mechanical loading and disc transport are still an open issue. A poromechanical finite element model of the human disc was coupled with oxygen and lactate transport models. Deformations and fluid flow were linked to transport predictions by including strain-dependent diffusion and advection. The two solute transport models were also coupled to account for cell metabolism. With this approach, the relevance of metabolic and mechano-transport couplings were assessed in the healthy disc under loading-recovery daily compression. Disc height, cell density and material degenerative changes were parametrically simulated to study their influence on the calculated solute concentrations. The effects of load frequency and amplitude were also studied in the healthy disc by considering short periods of cyclic compression. Results indicate that external loads influence the oxygen and lactate regional distributions within the disc when large volume changes modify diffusion distances and diffusivities, especially when healthy disc properties are simulated. Advection was negligible under both sustained and cyclic compression. Simulating degeneration, mechanical changes inhibited the mechanical effect on transport while disc height, fluid content, nucleus pressure and overall cell density reductions affected significantly transport predictions. For the healthy disc, nutrient concentration patterns depended mostly on the time of sustained compression and recovery. The relevant effect of cell density on the metabolic transport indicates the disturbance of cell number as a possible onset for disc degeneration via alteration of the metabolic balance. Results also suggest that healthy disc properties have a positive effect of loading on metabolic transport. Such relation, relevant to the maintenance of the tissue functional composition, would therefore link disc function with disc nutrition.     

MSC response to pH levels found in degenerating intervertebral discs

Painful degenerative disc disease is a major health problem and for successful tissue regeneration, MSCs must endure and thrive in a harsh disc microenvironment that includes matrix acidity as a critical factor. MSCs were isolated from bone marrow of Sprague-Dawley rats from two different age groups (<1 month, n = 6 and 4-5 months, n = 6) and cultured under four different pH conditions representative of the healthy, mildly or severely degenerated intervertebral disc (pH 7.4, 7.1, 6.8, and 6.5) for 5 days. Acidity caused an inhibition of aggrecan, collagen-1, and TIMP-3 expression, as well as a decrease in proliferation and viability and was associated with a change in cell morphology. Ageing had generally minor effects but young MSCs maintained greater mRNA expression levels. As acidic pH levels are typical of increasingly degenerated discs, our findings demonstrate the importance of early interventions and predifferentiation when planning to use MSCs for reparative treatments.     

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.     

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.     

缺氧诱导因子-1α和血管内皮生长因子在突出腰椎间盘中的表达及意义

目的探讨缺氧诱导因子-1α（hypoxia—induciblefactor-1α,HIF-1α）和血管内皮生长因子（vascular endothelial growth factor,VEGF）在突出腰椎间盘组织中的表达及意义。方法采用链霉亲和素-过氧化物酶复合物（SABC）免疫组化方法,测定40例腰椎间盘突出症患者椎间盘组织中HIF-1α和VEGF的表达情况。结果退变椎间盘组织中HIF-1α和VEGF呈高表达,HIF-1α和VEGF在髓核的表达显著高于纤维环;纤维环破裂型显著高于纤维环完整型;各组中HIF-1α和VEGF的表达均高度相关。结论HIF-1α和VEGF共同参与了椎间盘退变;HIF-1α可能通过上调VEGF的表达来促进椎间盘组织中新生血管的形成,进而延缓椎间盘退变的发生。