Modeling of neutral solute transport in a dynamically loaded porous permeable gel: implications for articular cartilage biosynthesis and tissue engineering

A primary mechanism of solute transport in articular cartilage is believed to occur through passive diffusion across the articular surface, but cyclical loading has been shown experimentally to enhance the transport of large solutes. The objective of this study is to examine the effect of dynamic loading within a theoretical context, and to investigate the circumstances under which convective transport induced by dynamic loading might supplement diffusive transport. The theory of incompressible mixtures was used to model the tissue (gel) as a mixture of a gel solid matrix (extracellular matrix/scaffold), and two fluid phases (interstitial fluid solvent and neutral solute), to solve the problem of solute transport through the lateral surface of a cylindrical sample loaded dynamically in unconfined compression with frictionless impermeable platens in a bathing solution containing an excess of solute. The resulting equations are governed by nondimensional parameters, the most significant of which are the ratio of the diffusive velocity of the interstitial fluid in the gel to the solute diffusivity in the gel (Rg), the ratio of actual to ideal solute diffusive velocities inside the gel (Rd), the ratio of loading frequency to the characteristic frequency of the gel (f), and the compressive strain amplitude (epsilon 0). Results show that when Rg > 1, Rd < 1, and f > 1, dynamic loading can significantly enhance solute transport into the gel, and that this effect is enhanced as epsilon 0 increases. Based on representative material properties of cartilage and agarose gels, and diffusivities of various solutes in these gels, it is found that the ranges Rg > 1, Rd < 1, correspond to large solutes, whereas f > 1 is in the range of physiological loading frequencies. These theoretical predictions are thus in agreement with the limited experimental data available in the literature. The results of this study apply to any porous hydrated tissue or material, and it is therefore plausible to hypothesize that dynamic loading may serve to enhance solute transport in a variety of physiological processes.     

Static compression of articular cartilage can reduce solute diffusivity and partitioning: implications for the chondrocyte biological response.

Chondrocytes depend upon solute transport within the avascular extracellular matrix of adult articular cartilage for many of their biological activities. Alterations to bioactive solute transport may, therefore, represent a mechanism by which cartilage compression is transduced into cellular metabolic responses. We investigated the effects of cartilage static compression on diffusivity and partitioning of a range of model solutes including dextrans of molecular weights 3 and 40 kDa, and tetramethylrhodamine (a 430 Da fluorophore). New fluorescence methods were developed for real-time visualization and measurement of transport within compressed cartilage explants. Experimental design allowed for multiple measurements on individual explants at different compression levels in order to minimize confounding influences of compositional variations. Results demonstrate that physiological levels of static compression may significantly decrease solute diffusivity and partitioning in cartilage. Effects of compression were most dramatic for the relatively high molecular weight solutes. For 40 kDa dextran, diffusivity decreased significantly (p<0.01) between 8% and 23% compression, while partitioning of 3 and 40 kDa dextran decreased significantly (p<0.01) between free-swelling conditions and 8% compression. Since diffusivity and partitioning can influence pericellular concentrations of bioactive solutes, these observations support a role for perturbations to solute transport in mediating the cartilage biological response to compression.     

Effects of tension-compression nonlinearity on solute transport in charged hydrated fibrous tissues under dynamic unconfined compression

Cartilage is a charged hydrated fibrous tissue exhibiting a high degree of tension-compression nonlinearity (i.e., tissue anisotropy). The effect of tension-compression nonlinearity on solute transport has not been investigated in cartilaginous tissue under dynamic loading conditions. In this study, a new model was developed based on the mechano-electrochemical mixture model [Yao and Gu, 2007, J. Biomech. Model Mechanobiol., 6, pp. 63-72, Lai et al., 1991, J. Biomech. Eng., 113, pp. 245-258], and conewise linear elasticity model [Soltz and Ateshian, 2000, J. Biomech. Eng., 122, pp. 576-586; Curnier et al., 1995, J. Elasticity, 37, pp. 1-38]. The solute desorption in cartilage under unconfined dynamic compression was investigated numerically using this new model. Analyses and results demonstrated that a high degree of tissue tension-compression nonlinearity could enhance the transport of large solutes considerably in the cartilage sample under dynamic unconfined compression, whereas it had little effect on the transport of small solutes (at 5% dynamic strain level). The loading-induced convection is an important mechanism for enhancing the transport of large solutes in the cartilage sample with tension-compression nonlinearity. The dynamic compression also promoted diffusion of large solutes in both tissues with and without tension-compression nonlinearity. These findings provide a new insight into the mechanisms of solute transport in hydrated, fibrous soft tissues.     

Cation diffusion in cartilage measured by pulsed field gradient NMR

In this study, the pulsed field gradient (PFG) nuclear magnetic resonance (NMR) technique was used for the investigation of (1) concentration and compression effects on cation self-diffusion, and (2) restricted diffusion of cations in cartilage. Since physiologically relevant cations like Na+ are difficult to investigate owing to their very short relaxation times, the cations tetramethylammonium (TMA) and tetraethylammonium (TEA) were employed for diffusion studies in samples of explanted cartilage. Results indicated that the diffusion of monovalent cations shows strong similarities to observations already made in studies of the diffusion of water in cartilage: with increasing compression, i.e. decreasing water content, the diffusion coefficient of the cation decreases concomitantly. The diffusion coefficients also showed a decrease with increasing cation concentrations, basically reflecting the corresponding decrease in the water content. Both results could be explained by the well-established model of Mackie and Meares. This, together with the similarity of the diffusion coefficient D in cartilage relative to free solution (about 50%) for both cations, is consistent with the view that the water content and not the charge is the most important determinant of the intratissue diffusivity of monovalent cations. Diffusion studies with increasing observation times showed strong evidence of restricted diffusion, allowing the estimation of the geometry of barriers within cartilage.     

Quantitative spatially resolved measurements of mass transfer through laryngeal cartilage.

The scanning electrochemical microscope (SECM) is a scanned probe microscope that uses the response of a mobile ultramicroelectrode (UME) tip to determine the reactivity, topography, and mass transport characteristics of interfaces with high spatial resolution. SECM strategies for measuring the rates of solute diffusion and convection through samples of cartilage, using amperometric UMEs, are outlined. The methods are used to determine the diffusion coefficients of oxygen and ruthenium(III) hexamine [Ru(NH3)6(3+)] in laryngeal cartilage. The diffusion coefficient of oxygen in cartilage is found to be approximately 50% of that in aqueous electrolyte solution, assuming a partition coefficient of unity for oxygen between cartilage and aqueous solution. In contrast, diffusion of Ru(NH3)6(3+) within the cartilage sample cannot be detected on the SECM timescale, suggesting a diffusion coefficient at least two orders of magnitude lower than that in solution, given a measured partition coefficient for Ru(NH3)6(3+) between cartilage and aqueous solution, Kp = [Ru(NH3)6(3+)]cartilage/[RU(NH3)6(3+)]solution = 3.4 +/- 0.1. Rates of Ru(NH3)6(3+) osmotically driven convective transport across cartilage samples are imaged at high spatial resolution by monitoring the current response of a scanning UME, with an osmotic pressure of approximately 0.75 atm across the slice. A model is outlined that enables the current response to be related to the local flux. By determining the topography of the sample from the current response with no applied osmotic pressure, local transport rates can be correlated with topographical features of the sample surface, at much higher spatial resolution than has previously been achieved.     

Single and Multicomponent Gas Phase Diffusion in a Porous Media: Modeling and Laboratory Measurements

Subsurface vapor migration of volatile chemicals may impact ambient and indoor air quality, increasing the importance to investigate the fate and transport of these chemicals. This project involved both modeling and experimental work to study the vapor phase transport behavior of single, binary, and tertiary component systems present in the gas phase. The experimental phase resulted in the development of a diffusion cell for measuring vapor phase transport. Three organic compounds (toluene, cumene, and isooctane) common to petroleum-based products were selected. The objective of this research project was to evaluate how the rate of a component diffusing alone in a stagnant gas mixture compares to the rate of the same component when diffusing in the presence of multiple diffusing species. The equipment was first validated by measuring the unobstructed gas phase diffusion fluxes for each organic compound. The diffusion coefficients were then calculated from the experimentally measured diffusive fluxes using Fick's Law at 20 and 25°C and compared to the respective literature values. The experimental/literature (E/L) ratio was calculated for toluene, cumene, and isooctane. The range of the average E/L ratio for the single component data sets is 0.93 to 1.05. The validation data provided the baseline for extending the research to multicomponent data. The multi-component systems research was characterized as either binary systems or a three-component system. The binary systems were either isooctane/tolu-ene or isooctane/cumene. The three-component system consisted of a mixture of all three compounds. For both temperatures and all compounds the flux rate decreased for any single component due to the dilution effect by incorporation into a mixture. Applying Fick's Law to calculate the effective diffusion coefficient for each compound that corresponded to the resulting concentration gradient by the mixture, an enhancement in the diffusive flux of each individual species was observed. This enhancement can be explained by a compositional coupling of each component to all others which results in a total vapor phase mass flux comprised of both diffusive and pseudo-advective mass transport. This pseudo-advective component is attributed to simultaneous diffusion of other species in the presence of the one of interest. Since this research project incorporated a mixture of toluene and cumene present in a background carrier solvent of isooctane, by calculating the ratios Dexp(3-component)/ Dexp(2-component) and Dexp(2-compo-nent)/Dexp(single component), an estimate is obtained of the enhancement effect due to the advective component of simultaneously diffusing chemicals. The diffusivity ratios for the three-component system compared to the dual component system ranged from 0.8 to 3.7. The diffusivity ratios for individual compounds were for 1.5-3.7 cumene, 0.8-1.2 for toluene, and 1.0-1.2 for isooctane. The diffusivity ratios for the dual component system to the single component systems ranged from 0.8 to 4.0. The range of diffusivity ratios for individual compounds were for 2.0-4.0 for cumene, 0.8-1.6 for toluene, and 1.1-1.4 for isooctane. A ratio greater than 1.0 indicated an enhancement effect on the molecular diffusion rate due to the presence of one or more additional diffusing chemical species present. The majority of fate and transport models are based on single component behavior modeled by Fick's Law using the pure gas phase diffusion coefficient. The enhancement of the individual diffusive flux in a multicomponent mixture observed in this study and accounted for by pseudo-advec-tive mass transport results in an under-prediction of the actual multicomponent diffusive fluxes. It is recommended that a more rigorous diffusion equation such as the Stefan-Maxwell equation be considered for incorporation into vapor phase transport models when modeling multicomponent/ contaminant systems.     

A method for the determination of diffusion coefficients for small molecules in aqueous solution

A method is described for determining the diffusion coefficients of small solutes in limited volumes (approximately equal to 4-9 ml) of fluid. Diffusion is measured in a three-chamber diffusion cell across a central unstirred compartment. Compartments are separated by nitrocellulose membranes. The instantaneous concentration gradient and the instantaneous flux of solute into the dilute end compartment are derived from changes in the concentration of solute in the two stirred end compartments through time. The diffusion coefficient is calculated from the slope of the least-squares regression line relating the magnitude of the instantaneous solute flux to that of the instantaneous concentration gradient. The apparatus is calibrated with a solute of known diffusivity (KCl). Diffusion coefficients thus determined in water at 25 degrees C for CaCl2 (7.54 X 10(-6) cm2.s-1), Na2-ATP (7.01 X 10(-6) cm2.s-1), 2-deoxyglucose (5.31 X 10(-6) cm2.s-1), and D-Na-lactate (5.62 X 10(-6) cm2.s-1) differed by an average of 3.7% from literature values. The method described results in accurate estimates of diffusion coefficients by a simple and relatively rapid procedure.     

Revisiting the Münch pressure-flow hypothesis for long-distance transport of carbohydrates: modelling the dynamics of solute transport inside a semipermeable tube

A mathematical model of the Münch pressure-flow hypothesis for long-distance transport of carbohydrates via sieve tubes is constructed using the Navier-Stokes equation for the motion of a viscous fluid and the van't Hoff equation for osmotic pressure. Assuming spatial dimensions that are appropriate for a sieve tube and ensuring suitable initial profiles of the solute concentration and solution velocity lets the model become mathematically tractable and concise. In the steady-state case, it is shown via an analytical expression that the solute flux is diffusion-like with the apparent diffusivity coefficient being proportional to the local solute concentration and around seven orders of magnitude greater than a diffusivity coefficient for sucrose in water. It is also shown that, in the steady-state case, the hydraulic conductivity over one metre can be calculated explicitly from the tube radius and physical constants and so can be compared with experimentally determined values. In the time-dependent case, it is shown via numerical simulations that the solute (or water) can simultaneously travel in opposite directions at different locations along the tube and, similarly, change direction of travel over time at a particular location along the tube.     

Coupling of Solute and Solvent Flows in Porous Lipid Bilayer Membranes

The present experiments were designed to evaluate coupling of water and nonelectrolyte flows in porous lipid bilayer membranes (i.e., in the presence of amphotericin B) in series with unstirred layers. Alterations in solute flux during osmosis, with respect to the flux in the absence of net water flow, could be related to two factors: first, changes in the diffusional component of solute flux referable to variations in solute concentrations at the membrane interfaces produced by osmotic flow through the unstirred layers; and second, coupling of solute and solvent flows within the membrane phase. Osmotic water flow in the same direction as solute flow increased substantially the net fluxes of glycerol and erythritol through the membranes, while osmotic flow in the opposite direction to glycerol flow reduced the net flux of that solute. The observed effects of osmotic water flow on the fluxes of these solutes were in reasonable agreement with predictions based on a model for coupling of solute and solvent flows within the membrane phase, and considerably in excess of the prediction for a diffusion process alone.     

Pore network modelling of affinity chromatography: determination of the dynamic profiles of the pore diffusivity of beta-galactosidase and its effect on column performance as the loading of beta-galactosidase onto anti-beta-galactosidase varies with time.

A three-dimensional pore network model for diffusion in porous adsorbent particles was employed in a dynamic adsorption model that simulates the adsorption of a solute in porous particles packed in a chromatographic column. The solution of the combined model yielded the dynamic profiles of the pore diffusion coefficient of beta-galactosidase along the radius of porous adsorbent particles and along the length of the column as the loading of beta-galactosidase onto anti-beta-galactosidase immobilized on the surface of the pores of the particles occurred, and, the dynamic adsorptive capacity of the chromatographic column as a function of the design and operational parameters of the chromatographic system. It was found that for a given column length the dynamic profiles of the pore diffusion coefficient were influenced by (a) the superficial fluid velocity in the column, (b) the diameter of the adsorbent particles, and (c) the pore connectivity of the porous structure of the adsorbent particles. The effect of the magnitude of the pore connectivity on the dynamic profiles of the pore diffusion coefficient of beta-galactosidase increased as the diameter of the adsorbent particles and the superficial fluid velocity in the column increased. The dynamic adsorptive capacity of the column increased as (i) the particle diameter and the superficial fluid velocity in the column decreased, and (ii) the column length and the pore connectivity increased. In preparative affinity chromatography, it is desirable to obtain high throughputs within acceptable pressure gradients, and this may require the employment of larger diameter adsorbent particles. In such a case, longer column lengths satisfying acceptable pressure gradients with adsorbent particles having higher pore connectivity values could provide high dynamic adsorptive capacities. An alternative chromatographic system could be comprised of a long column packed with large particles which have fractal pores (fractal particles) that have high pore connectivities and which allow high intraparticle diffusional and convective flow mass transfer rates providing high throughputs and high dynamic adsorptive capacities. If large scale monoliths could be made to be reproducible and operationally stable, they could also offer an alternative mode of operation that could provide high throughputs and high dynamic adsorptive capacities.     

Transport of neutral solute in articular cartilage: effect of microstructure anisotropy

Due to the avascular nature of articular cartilage, solute transport through its extracellular matrix is critical for the maintenance and the functioning of the tissue. What is more, diffusion of macromolecules may be affected by the microstructure of the extracellular matrix in both undeformed and deformed cartilage and experiments demonstrate diffusion anisotropy in the case of large solute. However, these phenomena have not received sufficient theoretical attention to date. We hypothesize here that the diffusion anisotropy of macromolecules is brought about by the particular microstructure of the cartilage network. Based on this hypothesis, we then propose a mathematical model that correlates the diffusion coefficient tensor with the structural orientation tensor of the network. This model is shown to be successful in describing anisotropic diffusion of macromolecules in undeformed tissue and is capable of clarifying the effects of network reorientation as the tissue deforms under mechanical load. Additionally, our model explains the anomaly that at large strain, in a cylindrical plug under unconfined compression, solute diffusion in the radial direction increases with strain. Our results indicate that in cartilage the degree of diffusion anisotropy is site specific, but depends also on the size of the diffusing molecule. Mechanical loading initiates and/or further exacerbates this anisotropy. At small deformation, solute diffusion is near isotropic in a tissue that is isotropic in its unstressed state, becoming anisotropic as loading progresses. Mechanical loading leads to an attenuation of solute diffusion in all directions when deformation is small. However, loading, if it is high enough, enhances solute transport in the direction perpendicular to the load line, instead of inhibiting it.     

Site-specific effects of compression on macromolecular diffusion in articular cartilage

Articular cartilage is the connective tissue that lines joints and provides a smooth surface for joint motion. Because cartilage is avascular, molecular transport occurs primarily via diffusion or convection, and cartilage matrix structure and composition may affect diffusive transport. Because of the inhomogeneous compressive properties of articular cartilage, we hypothesized that compression would decrease macromolecular diffusivity and increase diffusional anisotropy in a site-specific manner that depends on local tissue strain. We used two fluorescence photobleaching methods, scanning microphotolysis and fluorescence imaging of continuous point photobleaching, to measure diffusion coefficients and diffusional anisotropy of 70 kDa dextran in cartilage during compression, and measured local tissue strain using texture correlation. For every 10% increase in normal strain, the fractional change in diffusivity decreased by 0.16 in all zones, and diffusional anisotropy increased 1.1-fold in the surface zone and 1.04-fold in the middle zone, and did not change in the deep zone. These results indicate that inhomogeneity in matrix structure and composition may significantly affect local diffusive transport in cartilage, particularly in response to mechanical loading. Our findings suggest that high strains in the surface zone significantly decrease diffusivity and increase anisotropy, which may decrease transport between cartilage and synovial fluid during compression.     

Simulation of bone tissue formation within a porous scaffold under dynamic compression

A computational model of mechanoregulation is proposed to predict bone tissue formation stimulated mechanically by overall dynamical compression within a porous polymeric scaffold rendered by micro-CT. Dynamic compressions of 0.5–5% at 0.0025–0.025 s⁻¹ were simulated. A force-controlled dynamic compression was also performed by imposing a ramp of force from 1 to 70 N. The model predicts homogeneous mature bone tissue formation under strain levels of 0.5–1% at strain rates of 0.0025–0.005 s⁻¹. Under higher levels of strain and strain rates, the scaffold shows heterogeneous mechanical behaviour which leads to the formation of a heterogeneous tissue with a mixture of mature bone and fibrous tissue. A fibrous tissue layer was also predicted under the force-controlled dynamic compression, although the same force magnitude was found promoting only mature bone during a strain-controlled compression. The model shows that the mechanical stimulation of bone tissue formation within a porous scaffold closely depends on the loading history and on the mechanical behaviour of the scaffold at local and global scales.     

Solute convection in dynamically compressed cartilage

Chondrocytes depend upon solute transport within the avascular extracellular matrix of articular cartilage for many of their biological activities. Alterations to solute transport parameters may therefore mediate the cell response to tissue compression. While interstitial solute transport may be supplemented by convection during dynamic tissue compression, matrix compression is also associated with decreased diffusivities. Such trade-offs between increased convection and decreased diffusivities of solutes in dynamically compressed cartilage remain largely unexplored. We measured diffusion and convection coefficients of a wide range of solutes in mature bovine cartilage explant disks subjected to radially unconfined axial ramp compression and release. Solutes included approximately 500 Da fluorophores bearing positive and negative charges, and 10 kDa dextrans bearing positive, neutral, and negative charges. Significantly positive values of convection coefficients were measured for several different solutes. Findings therefore support a role for solute convection in mediating the cartilage biological response to dynamic compression.     

Static compression is associated with decreased diffusivity of dextrans in cartilage explants

The chondrocytes of adult articular cartilage rely upon transport phenomena within their avascular extracellular matrix for many biological activities. Therefore, changes in matrix structure which influence cytokine transport parameters may be an important mechanism involved in the chondrocyte response to tissue compression. With this hypothesis in mind, partitioning and diffusion of 3-, 10-, and 40-kDa dextrans conjugated to tetramethylrhodamine, and 430-Da tetramethylrhodamine itself, were measured within statically compressed bovine articular cartilage explants using a novel experimental apparatus and desorption fluorescence method. Partitioning and diffusion were examined as functions of solute molecular weight and matrix proteoglycan density, and diffusion was measured versus static compression up to 35% volumetric strain. In general, partition coefficients and diffusivities were found to decrease with increasing solute molecular weight. In addition, for a given solute, diffusivities decreased significantly with increasing static compression. Results therefore suggest a possible role for transport limitations of relatively large molecular weight solutes within the extracellular matrix in mediating the biological response of chondrocytes to cartilage compression.     

Diffusion in kappa-carrageenan gel beads

Using a well-mixed and temperature-led vessel, the diffusion characteristics of various solutes into spherical kappa-carrageenan gel beads were experimentally investigated. The diffusion coefficient of glucose was markedly affected by the glucose concentration and the operating temperature. In all cases the diffusivity obtained was noticeably smaller than that of glucose in pure water. The experimental data also indicated an inverse relationship between the diffusivity and the polymer concentration used in the gel preparation. As well, the glucose diffusivity was affected by the presence of other solutes in the glucose solution. Electrolytes such as ammonium sulfate, KCl, and CaCl(2) were observed to enhance the diffusion coefficient. On the other hand, the addition of arginine or bovine serum albumin had an adverse effect on the diffusivity. No diffusion of albumin into the gel beads was observed, and such a solute created a significant mass transfer resistance during the diffusion process.     

Modeling water flow through arterial tissue

A simple model of, water flow through deformable porous media has been developed with emphasis on application to arterial walls. The model incorporates a strain-dependent permeability function into Darcy's Law which is coupled, to the force balance for the bulk material. A simple analytical expression relating water flux (volume flux) to pressure differential is developed which shows how strain-dependent permeability can lead to a reduction in hydraulic conductivity with increasing differential pressure as observed in experiments with arteries. The variation of permeability with position in the wall, which may influence the convective diffusion of macromolecules, is determined for both cylindrical and planar segments and a marked influence of geometry is noted.     

Temperature affects transport of polysaccharides and proteins in articular cartilage explants

Solute transport phenomena mediate many aspects of the physiology and contrast agent-based clinical imaging of articular cartilage. Temperatures up to 10°C below standard body temperature (37°C) are common in articulating joints during normal activities and clinically (e.g. cold treatment of injuries). Therefore it is of interest to characterize the effects of temperature changes on solute transport parameters in cartilage. A range of fluorescent solutes including fluorescein isothiocyanate, 4 and 40kDa dextrans, myoglobin, insulin and chondroitin sulfate were prepared and used in assays of solute effective partition coefficient and effective diffusivity in bovine intermediate zone articular cartilage explants maintained at 10, 22 or 37°C. Trends for increasing partition coefficient with increasing temperature were evident for all solutes except chondroitin sulfate, with significant changes between 22 and 37°C for 4kDa dextran, insulin and myoglobin. Diffusivities of most solutes tested also tended to increase with increasing temperature, with significant changes between 10 and 22°C for FITC, 40kDa dextran and myoglobin. Oddly, insulin diffusivity decreased significantly as temperature increased from 22 to 37°C while chondroitin sulfate diffusivity exhibited no clear temperature dependence. These results highlight solute-specific temperature dependences of transport phenomena which may depend upon molecular weight, chemical structure, molecular conformation, and solute-matrix and solute-solute interactions. The articular cartilage explants themselves exhibited small but significant changes in water and glycosaminoglycan contents during experiments, underscoring the importance of solute-matrix interactions. Solute transport parameters in cartilage and their temperature dependences are therefore not easily predicted, and case-by-case experimental determination may be essential.     

Transient solute diffusion in articular cartilage

The one-dimensional transient diffusion of glucose, inulin and dextran into adult bovine knee articular cartilage was determined for transport times of 1, 5, 15 and 60 min, and 4, 12, 24 and 48 h. The apparent diffusion coefficient and apparent interface partition coefficient were calculated from the concentration-depth profiles within the tissue using a theoretical model for non-steady state solute diffusion. The diffusion coefficient was found to decrease with both solute size and transport time. The partition coefficient also decreased with solute size but increased with transport time. Neither coefficient was dependent on normal tissue fluid or proteoglycan content variations.     

Diffusion of MRI and CT Contrast Agents in Articular Cartilage under Static Compression

Cartilage has a limited capacity for self-repair and focal damage can eventually lead to complete degradation of the tissue. Early diagnosis of degenerative changes in cartilage is therefore essential. Contrast agent-based computed tomography and magnetic resonance imaging provide promising tools for this purpose. However, the common assumption in clinical applications that contrast agents reach steady-state distributions within the tissue has been of questionable validity. Characterization of nonequilibrium diffusion of contrast agents rather than their equilibrium distributions may therefore be more effective for image-based cartilage assessment. Transport of contrast agent through the extracellular matrix of cartilage can be affected by tissue compression due to matrix structural and compositional changes including reduced pore size and fluid content. We therefore investigate the effects of static compression on diffusion of three common contrast agents: sodium iodide, sodium diatrizoate, and gadolinium diethylenetriamine-pentaacid (Gd-DTPA). Results showed that static compression was associated with significant decreases in diffusivities for sodium iodide and Gd-DTPA, with similar (but not significant) trends for sodium diatrizoate. Molecular mass of contrast agents affected diffusivities as the smallest one tested, sodium iodide, showed higher diffusivity than sodium diatrizoate and Gd-DTPA. Compression-associated cartilage matrix alterations such as glycosaminoglycan and fluid contents were found to correspond with variations in contrast agent diffusivities. Although decreased diffusivity was significantly correlated with increasing glycosaminoglycan content for sodium iodide and Gd-DTPA only, diffusivity significantly increased for all contrast agents by increasing fluid fraction. Because compounds based on iodine and gadolinium are commonly used for computed tomography and magnetic resonance imaging, present findings can be valuable for more accurate image-based assessment of variations in cartilage composition associated with focal injuries.