In this study, a new cell density model was developed and incorporated into the formulation of the mechano-electrochemical mixture theory to investigate the effects of deprivation of nutrition supply at boundary source, degeneration, and dynamic loading on the cell viability of intervertebral disc (IVD) using finite element methods. The deprivation of nutrition supply at boundary source was simulated by reduction in nutrition level at CEP and AF boundaries. Cases with 100%, 75%, 60%, 50% and 30% of normal nutrition level at both CEP and AF boundaries were modeled. Unconfined axial sinusoidal dynamic compressions with different combinations of amplitude (u=10%±2.5%, ±5%) and frequency (f=1, 10, 20 cycle/day) were applied. Degenerated IVD was modeled with altered material properties. Cell density decreased substantially with reduction of nutrition level at boundaries. Cell death was initiated primarily near the NP–AF interface on the mid-plane. Dynamic loading did not result in a change in the cell density in non-degenerated IVD, since glucose levels did not fall below the minimum value for cell survival; in degenerated IVDs, we found that increasing frequency and amplitude both resulted in higher cell density, because dynamic compression facilitates the diffusion of nutrients and thus increases the nutrition level around IVD cells. The novel computational model can be used to quantitatively predict both when and where cells start to die within the IVD under various kinds of nutritional and mechanical conditions.

Study Design Simulate the progression of human disc degeneration. Objective The objective of this study was to quantitatively analyze and simulate the changes in cell density, nutrition level, proteoglycan content, water content, and volume change during human disc degeneration using a numerical method. Summary of Background Data Understanding the etiology and progression of intervertebral disc (IVD) degeneration is crucial for developing effective treatment strategies for IVD-degeneration related diseases. During tissue degeneration, the disc undergoes losses of cell viability and activities, changes in extracellular matrix composition and structure, and compromise of the tissue-level integrity and function, which is significantly influenced by the inter-coupled biological, chemical, electrical, and mechanical signals in the disc. Characterizing these signals in human discs in vivo is difficult. Methods A realistic 3D finite element model of the human IVD was developed based on biomechano-electrochemical continuum mixture theory. The theoretical framework and the constitutive relationships were all biophysics based. All the material properties were obtained from experimental results. The cell-mediated disc degeneration process caused by lowered nutrition levels at disc boundaries was simulated and validated by comparing with experimental results. Results Cell density reached equilibrium state in 30 days after reduced nutrition supply at the disc boundary, while the proteoglycan (PG) and water contents reached a new equilibrium state in 55 years. The simulated results for the distributions of PG and water contents within the disc were consistent with the results measured in the literature, except for the distribution of PG content in the sagittal direction. Conclusions Poor nutrition supply has a long-term effect on disc degeneration.

In this study, a three-dimensional finite element model was used to investigate the changes in tissue composition and mechanical signals within human lumbar intervertebral disc during the degenerative progression. This model was developed based on the cell-activity coupled mechano-electrochemical mixture theory. The disc degeneration was simulated by lowering nutrition levels at disc boundaries, and the temporal and spatial distributions of the fixed charge density, water content, fluid pressure, Von Mises stress, and disc deformation were analyzed. Results showed that fixed charge density, fluid pressure, and water content decreased significantly in the nucleus pulposus (NP) and the inner to middle annulus fibrosus (AF) regions of the degenerative disc. It was found that, with degenerative progression, the Von Mises stress (relative to that at healthy state) increased within the disc, with a larger increase in the outer AF region. Both the disc volume and height decreased with the degenerative progression. The predicted results of fluid pressure change in the NP were consistent with experimental findings in the literature. The knowledge of the variations of temporal and spatial distributions of composition and mechanical signals within the human IVDs provide a better understanding of the progression of disc degeneration.

Study Design Investigation of the effects of the impairment of different nutritional pathways on the intervertebral disc degeneration patterns in terms of spatial distributions of cell density, glycosaminoglycan content, and water content. Objective To test the hypothesis that impairment of different nutritional pathways would result in different degenerative patterns in human discs. Summary of Background Data Impairment of nutritional pathways has been found to affect cell viability in the disc. However, details on how impairment of different nutritional pathways affects the disc degeneration patterns are unknown. Methods A 3D finite element model was used for this study. This finite element method was based on the cell-activity-coupled mechano-electrochemical theory for cartilaginous tissues. Impairment of the nutritional pathways was simulated by lowering the nutrition level at the disc boundaries. Effects of the impartment of cartilaginous endplate-nucleus pulposus (CEP-NP) pathway only (Case 1), annulus fibrosus (AF) pathway only (Case 2), and both pathways (Case 3) on disc degeneration patterns were studied. Results The predicted critical level of nutrition for Case 1, Case 2, and Case 3 were around 30%, 20%, and 50% of the reference values, respectively. Below this critical level, the disc degeneration would occur. Disc degeneration appeared mainly in the NP for Case 1, in the outer AF for Case 2, and in both the NP and inner to middle AF for Case 3. For Cases 1 and 3, the loss of water content was primarily located in the mid-axial plane, which is consistent with the horizontal gray band seen in some T2-weighted MRI images. For the disc geometry used in this study, it was predicted that there existed a High Intensity Zone (for Case 3), as seen in some T2-weighted MRI images. Conclusion Impairment of different nutrition pathways results in different degenerative patterns.