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.     

Effects of mechanical compression on metabolism and distribution of oxygen and lactate in intervertebral disc

The objective of this study was to examine the effects of mechanical compression on metabolism and distributions of oxygen and lactate in the intervertebral disc (IVD) using a new formulation of the triphasic theory. In this study, the cellular metabolic rates of oxygen and lactate were incorporated into the newly developed formulation of the mechano-electrochemical mixture model [Huang, C.-Y., Gu, W.Y., 2007. Effect of tension-compression nonlinearity on solute transport in charged hydrated fibrosus tissues under dynamic unconfined compression. Journal of Biomechanical Engineering 129, 423-429]. The model was used to numerically analyze metabolism and transport of oxygen and lactate in the IVD under static or dynamic compression. The theoretical analyses demonstrated that compressive loading could affect transport and metabolism of nutrients. Dynamic compression increased oxygen concentration, reduced lactate accumulation, and promoted oxygen consumption and lactate production (i.e., energy conversion) within the IVD. Such effects of dynamic loading were dependent on strain level and loading frequency, and more pronounced in the IVD with less permeable endplate. In contrast, static compression exhibited inverse effects on transport and metabolism of oxygen and lactate. The theoretical predictions in this study are in good agreement with those in the literature. This study established a new theoretical model for analyzing cellular metabolism of nutrients in hydrated, fibrous soft tissues under mechanical compression.     

Numerical exploration of the combined effect of nutrient supply,tissue condition and deformation in the intervertebral disc

Novel strategies to heal discogenic low back pain could highly benefit from comprehensive biophysical studies that consider both mechanical and biological factors involved in intervertebral disc degeneration. A decrease in nutrient availability at the bone–disc interface has been indicated as a relevant risk factor and as a possible initiator of cell death processes. Mechanical behaviour of both healthy and degenerated discs could highly interact with cell death in these compromised situations. In the present study, a mechano-transport finite element model was used to investigate the nature of mechanical effects on cell death processes via load-induced metabolic transport variations. Cycles of static sustained compression were chosen to simulate daily human activity. Healthy and degenerated cases were simulated as well as a reduced supply of solutes and an increase in solute exchange area at the bone–disc interface. Results showed that a reduction in metabolite concentrations at the bone–disc boundaries induced cell death, even when the increased exchange area was simulated. Slight local mechanical enhancements of glucose in the disc centre were capable of decelerating cell death but occurred only with healthy mechanical properties. However, mechanical deformations were responsible for a worsening in terms of cell death in the inner annulus, a disadvantaged zone far from the boundary supply with both an increased cell demand and a strain-dependent decrease of diffusivity. Such adverse mechanical effects were more accentuated when degenerative properties were simulated. Overall, this study paves the way for the use of biophysical models for a more integrated understanding of intervertebral disc pathophysiology.     

Disc mechanics with trans-endplate partial nucleotomy are not fully restored following cyclic compressive loading and unloaded recovery

Mechanical function of the intervertebral disc is maintained through the interaction between the hydrated nucleus pulposus, the surrounding annulus fibrosus, and the superior and inferior endplates. In disc degeneration the normal transfer of load between disc substructures is compromised. The objective of this study was to explore the mechanical role of the nucleus pulposus in support of axial compressive loads over time. This was achieved by measuring the elastic slow ramp and viscoelastic stress-relaxation mechanical behaviors of cadaveric sheep motion segments before and after partial nucleotomy through the endplate (keeping the annulus fibrosus intact). Mechanics were evaluated at five conditions: Intact, intact after 10,000 cycles of compression, acutely after nucleotomy, following nucleotomy and 10,000 cycles of compression, and following unloaded recovery. Radiographs and magnetic resonance images were obtained to examine structure. Only the short time constant of the stress relaxation was altered due to nucleotomy. In contrast, cyclic loading resulted in significant and large changes to both the stiffness and stress relaxation behaviors. Moreover, the nucleotomy had little to no effect on the disc mechanics after cyclic loading, as there were no significant differences comparing mechanics after cyclic loading with or without the nucleotomy. Following unloaded recovery the mechanical changes that had occurred as a consequence of cyclic loading were restored, leaving only a sustained change in the short time constant due to the trans-endplate nucleotomy. Thus the swelling and redistribution of the remaining nucleus pulposus was not able to fully restore mechanical behaviors. This study reveals insights into the role of the nucleus pulposus in disc function, and provides new information toward the potential role of altered nucleus pulpous function in the degenerative cascade.     

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.     

Biological response of the intervertebral disc to dynamic loading

Disc degeneration is a chronic remodeling process that results in alterations of matrix composition and decreased cellularity. This study tested the hypothesis that dynamic mechanical forces are important regulators in vivo of disc cellularity and matrix synthesis. A murine model of dynamic loading was developed that used an external loading device to cyclically compress a single disc in the tail. Loads alternated at a 50% duty cycle between 0MPa and one of two peak stresses (0.9 or 1.3MPa) at one of two frequencies (0.1 or 0.01Hz) for 6h per day for 7 days. An additional group received static compression at 1.3MPa for 3h/day for 7 days. A control group wore the device with no loading. Sections of treated discs were analyzed for morphology, proteoglycan content, apoptosis, cell areal density, and aggrecan and collagen II gene expression. Dynamic loading induced differential effects that depended on frequency and stress. No significant changes to morphology, proteoglycan content or cell death were found after loading at 0.9MPa, 0.1Hz. Loading at lower frequency and/or higher stress increased proteoglycan content, matrix gene expression and cell death. The results have implications in the prevention of intervertebral disc degeneration, suggesting that loading conditions may be optimized to promote maintenance of normal structure and function.     

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.     

Nucleotomy reduces the effects of cyclic compressive loading with unloaded recovery on human intervertebral discs

The first objective of this study was to determine the effects of physiological cyclic loading followed by unloaded recovery on the mechanical response of human intervertebral discs. The second objective was to examine how nucleotomy alters the disc?s mechanical response to cyclic loading. To complete these objectives, 15 human L5-S1 discs were tested while intact and subsequent to nucleotomy. The testing consisted of 10,000 cycles of physiological compressive loads followed by unloaded hydrated recovery. Cyclic loading increased compression modulus (3%) and strain (33%), decreased neutral zone modulus (52%), and increased neutral zone strain (31%). Degeneration was not correlated with the effect of cyclic loading in intact discs, but was correlated with cyclic loading effects after nucleotomy, with more degenerate samples experiencing greater increases in both compressive and neutral zone strain following cyclic loading. Partial removal of the nucleus pulposus decreased the compression and neutral zone modulus while increasing strain. These changes correspond to hypermobility, which will alter overall spinal mechanics and may impact low back pain via altered motion throughout the spinal column. Nucleotomy also reduced the effects of cyclic loading on mechanical properties, likely due to altered fluid flow, which may impact cellular mechanotransduction and transport of disc nutrients and waste. Degeneration was not correlated with the acute changes of nucleotomy. Results of this study provide an ideal protocol and control data for evaluating the effectiveness of a mechanically-based disc degeneration treatment, such as a nucleus replacement.     

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.     

Damage accumulation location under cyclic loading in the lumbar disc shifts from inner annulus lamellae to peripheral annulus with increasing disc degeneration

It is difficult to study the breakdown of lumbar disc tissue over several years of exposure to bending and lifting by experimental methods. In our earlier published study we have shown how a finite element model of a healthy lumbar motion segment was used to predict the damage accumulation location and number of cyclic to failure under different loading conditions. The aim of the current study was to extend the continuum damage mechanics formulation to the degenerated discs and investigate the initiation and progression of mechanical damage. Healthy disc model was modified to represent degenerative discs (Thompson grade III and IV) by incorporating both geometrical and biochemical changes due to degeneration. Analyses predicted decrease in the number of cycles to failure with increasing severity of disc degeneration. The study showed that the damage initiated at the posterior inner annulus adjacent to the endplates and propagated outwards towards its periphery in healthy and grade III degenerated discs. The damage accumulated preferentially in the posterior region of the annulus. However in grade IV degenerated disc damage initiated at the posterior outer periphery of the annulus and propagated circumferentially. The finite element model predictions were consistent with the infrequent occurrence of rim lesions at early age but a much higher incidence in severely degenerated discs.     

Computation of coupled diffusion of oxygen, glucose and lactic acid in an intervertebral disc

The present numerical study aims to investigate the disc nutrition and factors affecting it by evaluating the concentrations of oxygen, glucose and lactic acid in the disc while accounting for the coupling between these species via the pH level in the tissue and the nonlinear concentration-consumption (for glucose and oxygen) and concentration-production (for lactate) relations. The effects of changes in the endplate exchange area (EA) adjacent to the nucleus or the inner annulus for the transport of nutrients and in the disc geometry as well as tissue diffusivities under static compression loading on species concentrations are also studied. Moreover, alterations in solute diffusion following a central endplate fracture are investigated. An axisymmetric geometry with four distinct regions is considered. Supply sources are assumed at the outer annulus periphery and disc endplates. Coupling between different solutes, pH level, endplate disruptions (calcifications and fractures) and mechanical loads substantially influenced the distribution of nutrients throughout the disc as well as the magnitude and location of critical concentrations; maximum for the lactic acid and minimum for oxygen and glucose. In cases with loss of endplate permeability and/or disruptions therein, as well as changes in geometry and fall in diffusivity associated with fluid expression, the nutrient concentrations could fall to levels inadequate to maintain cellular activity or viability, thus initiating or accelerating disc degeneration.     

Region Specific Response of Intervertebral Disc Cells to Complex Dynamic Loading: An Organ Culture Study Using a Dynamic Torsion-Compression Bioreactor

The spine is routinely subjected to repetitive complex loading consisting of axial compression, torsion, flexion and extension. Mechanical loading is one of the important causes of spinal diseases, including disc herniation and disc degeneration. It is known that static and dynamic compression can lead to progressive disc degeneration, but little is known about the mechanobiology of the disc subjected to combined dynamic compression and torsion. Therefore, the purpose of this study was to compare the mechanobiology of the intervertebral disc when subjected to combined dynamic compression and axial torsion or pure dynamic compression or axial torsion using organ culture. We applied four different loading modalities [1. control: no loading (NL), 2. cyclic compression (CC), 3. cyclic torsion (CT), and 4. combined cyclic compression and torsion (CCT)] on bovine caudal disc explants using our custom made dynamic loading bioreactor for disc organ culture. Loads were applied for 8 h/day and continued for 14 days, all at a physiological magnitude and frequency. Our results provided strong evidence that complex loading induced a stronger degree of disc degeneration compared to one degree of freedom loading. In the CCT group, less than 10% nucleus pulposus (NP) cells survived the 14 days of loading, while cell viabilities were maintained above 70% in the NP of all the other three groups and in the annulus fibrosus (AF) of all the groups. Gene expression analysis revealed a strong up-regulation in matrix genes and matrix remodeling genes in the AF of the CCT group. Cell apoptotic activity and glycosaminoglycan content were also quantified but there were no statistically significant differences found. Cell morphology in the NP of the CCT was changed, as shown by histological evaluation. Our results stress the importance of complex loading on the initiation and progression of disc degeneration.     

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.     

Effects of dynamic loading on solute transport through the human cartilage endplate

Nutrient and metabolite transport through the cartilage endplate (CEP) is important for maintaining proper disc nutrition, but the mechanisms of solute transport remain unclear. One unresolved issue is the role of dynamic loading. In comparison to static loading, dynamic loading is thought to enhance transport by increasing convection. However, the CEP has a high resistance to fluid flow, which could limit solute convection. Here we measure solute transport through site-matched cadaveric human lumbar CEP tissues under static vs. dynamic loading, and we determine how the degree of transport enhancement from dynamic loading depends on CEP porosity and solute size. We found that dynamic loading significantly increased small and large solute transport through the CEP: on average, dynamic loading increased the transport of sodium fluorescein (376 Da) by a factor of 1.85 ± 0.64 and the transport of a large dextran (4000 Da) by a factor of 4.97 ± 3.05. Importantly, CEP porosity (0.65 ± 0.07; range: 0.47–0.76) strongly influenced the degree of transport enhancement. Specifically, for both solutes, transport enhancement was greater for CEPs with low porosity than for CEPs with high porosity. This is because the CEPs with low porosity were susceptible to larger improvements in fluid flow under dynamic loading. The CEP becomes less porous and less hydrated with aging and as disc degeneration progresses. Together, these findings suggest that as those changes occur, dynamic loading has a greater effect on solute transport through the CEP compared to static loading, and thus may play a larger role in disc nutrition.     

Temporal changes of mechanical signals and extracellular composition in human intervertebral disc during degenerative progression

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.     

Preservation and analysis of nonequilibrium solute concentration distributions within mechanically compressed cartilage explants

Solute transport within articular cartilage is of central importance to tissue physiology, and may mediate effects of mechanical compression on cell metabolism. We therefore developed and applied a freeze-substitution method for fixation of cartilage explant disks which had been compressed axially during radial solute desorption. Dextrans were used as model solutes. Explant morphology was well preserved and nonequilibrium solute concentration distributions were stable for several hours at room temperature. For desorption from explants compressed statically to 0-46% strain, analysis of laser confocal images and comparison to a theoretical model permitted measurement of effective diffusivities. Results were consistent with previous studies suggesting a role for transport limitations in mediating the decreases of chondrocyte metabolic rates associated with static compression. In explants compressed dynamically (23+/-5% strain at 0.001 Hz), evidence was obtained for the augmentation of effective transport rate of 3 kDa dextrans by oscillatory interstitial fluid flows. This suggests that augmented solute transport may play a role in mediating the increases of chondrocyte metabolic rates associated with dynamic compression. Methods appear suitable for quantitative studies of transport within mechanically compressed cartilage-like tissues, and may be valuable for identification of loading environments which optimize solute transport in tissue engineering applications.     

Dynamic compression augments interstitial transport of a glucose-like solute in articular cartilage

Solute transport through the extracellular matrix is essential for cellular activities in articular cartilage. Increased solute transport via fluid convection may be a mechanism by which dynamic compression stimulates chondrocyte metabolism. However, loading conditions that optimally augment transport likely vary for different solutes. To investigate effects of dynamic loading on transport of a bioactive solute, triangular mechanical loading waveforms were applied to cartilage explants disks while interstitial transport of a fluorescent glucose analog was monitored. Peak-to-peak compression amplitudes varied from 5-50% and frequencies varied from 0.0006-0.1 Hz to alter the spatial distribution and magnitude of oscillatory fluid flow. Solute transport was quantified by monitoring accumulation of fluorescence in a saline bath circulated around the explant. Individual explants were subjected to a series of compression protocols, so that effects of loading on solute desorption could be observed directly. Maximum increases in solute transport were obtained with 10-20% compression amplitudes at 0.1 Hz; similar loading protocols were previously found to stimulate chondrocyte metabolism in vitro. Results therefore support hypotheses relating to increased solute transport as a mediator of the cartilage biological response to dynamic compression, and may have application in mechanical conditioning of cartilage constructs for tissue engineering.     

Is post-contrast MRI a valuable method for the study of the nutrition of the intervertebral disc?

Decreased nutrition has been proposed as a potential mechanism leading to intervertebral disc degeneration. A method to investigate it in vivo is the MRI evaluation of the transport of a paramagnetic contrast agent, which is assumed to diffuse through the endplate to the disc using the same mechanisms as the cell nutrients. However, previous numerical studies questioned the value of this method as a model to investigate disc nutrition. To assess its validity, a parametric osmoporoelastic finite element model of a lumbar intervertebral disc incorporating diffusion and convection of a solute (representing the contrast agent) was developed. A Taguchi sensitivity analysis was performed in order to assess the relevance of various parameters which influence the solute transport. Subsequently, a full-factorial sensitivity analysis was used to investigate specifically the diffusion coefficients of the contrast agent. The most important parameters in determining the results were the disc height, the diffusion coefficients and the pharmacokinetic of the contrast agent. However, diffusion coefficients values as measured in in vitro studies would lead to insubstantial enhancement of the MRI signal. Thus, transport mechanisms other than pure diffusion should be active in in vivo transport of the contrast agent. In conclusion, the study showed that post-contrast MRI may not be suited for a quantitative analysis, but only for a qualitative examination aimed for example to detect endplate lesions. Open questions remain on the use of post-contrast MRI for the investigation of the relevance of reduced nutrition as a trigger to disc degeneration.     

Fluid flow and convective transport of solutes within the intervertebral disc

Previous experimental and analytical studies of solute transport in the intervertebral disc have demonstrated that for small molecules diffusive transport alone fulfils the nutritional needs of disc cells. It has been often suggested that fluid flow into and within the disc may enhance the transport of larger molecules. The goal of the study was to predict the influence of load-induced interstitial fluid flow on mass transport in the intervertebral disc.An iterative procedure was used to predict the convective transport of physiologically relevant molecules within the disc. An axisymmetric, poroelastic finite-element structural model of the disc was developed. The diurnal loading was divided into discrete time steps. At each time step, the fluid flow within the disc due to compression or swelling was calculated. A sequentially coupled diffusion/convection model was then employed to calculate solute transport, with a constant concentration of solute being provided at the vascularised endplates and outer annulus. Loading was simulated for a complete diurnal cycle, and the relative convective and diffusive transport was compared for solutes with molecular weights ranging from 400 Da to 40 kDa.Consistent with previous studies, fluid flow did not enhance the transport of low-weight solutes. During swelling, interstitial fluid flow increased the unidirectional penetration of large solutes by approximately 100%. Due to the bi-directional temporal nature of disc loading, however, the net effect of convective transport over a full diurnal cycle was more limited (30% increase). Further study is required to determine the significance of large solutes and the timing of their delivery for disc physiology.