Inhomogeneous cartilage properties enhance superficial interstitial fluid support and frictional properties, but do not provide a homogeneous state of stress

It has been well established that articular cartilage is compositionally and mechanically inhomogenous through its depth. To what extent this structural inhomogeneity is a prerequisite for appropriate cartilage function and integrity is not well understood. The first hypothesis to be tested in this study was that the depth-dependent inhomogeneity of the cartilage acts to maximize the interstitial fluid load support at the articular surface, to provide efficient frictional and wear properties. The second hypothesis was that the inhomogeneity produces a more homogeneous state of elastic stress in the matrix than would be achieved with uniform properties. We have, for the first time, simultaneously determined depth-dependent tensile and compressive properties of human patellofemoral cartilage from unconfined compression stress relaxation tests. The results show that the tensile modulus increases significantly from 4.1 +/- 1.9 MPa in the deep zone to 8.3 +/- 3.7 MPa at the superficial zone, while the compressive modulus decreases from 0.73 +/- 0.26 MPa to 0.28 +/- 0.16 MPa. The experimental measurements were then implemented with the finite-element method to compute the response of an inhomogeneous and homogeneous cartilage layer to loading. The finite-element models demonstrate that structural inhomogeneity acts to increase the interstitial fluid load support at the articular surface. However, the state of stress, strain, or strain energy density in the solid matrix remained inhomogeneous through the depth of the articular layer, whether or not inhomogeneous material properties were employed. We suggest that increased fluid load support at the articular surface enhances the frictional and wear properties of articular cartilage, but that the tissue is not functionally adapted to produce homogeneous stress, strain, or strain energy density distributions. Interstitial fluid pressurization, but not a homogeneous elastic stress distribution, appears thus to be a prerequisite for the functional and morphological integrity of the cartilage.     

Depth-dependent Compressive Equilibrium Properties of Articular Cartilage Explained by its Composition

For this study, we hypothesized that the depth-dependent compressive equilibrium properties of articular cartilage are the inherent consequence of its depth-dependent composition, and not the result of depth-dependent material properties. To test this hypothesis, our recently developed fibril-reinforced poroviscoelastic swelling model was expanded to include the influence of intra- and extra-fibrillar water content, and the influence of the solid fraction on the compressive properties of the tissue. With this model, the depth-dependent compressive equilibrium properties of articular cartilage were determined, and compared with experimental data from the literature. The typical depth-dependent behavior of articular cartilage was predicted by this model. The effective aggregate modulus was highly strain-dependent. It decreased with increasing strain for low strains, and increases with increasing strain for high strains. This effect was more pronounced with increasing distance from the articular surface. The main insight from this study is that the depth-dependent material behavior of articular cartilage can be obtained from its depth-dependent composition only. This eliminates the need for the assumption that the material properties of the different constituents themselves vary with depth. Such insights are important for understanding cartilage mechanical behavior, cartilage damage mechanisms and tissue engineering studies.     

Two-dimensional strain fields on the cross-section of the bovine humeral head under contact loading

The objective of this study was to provide a detailed experimental assessment of the two-dimensional cartilage strain distribution on the cross-section of immature and mature bovine humeral heads subjected to contact loading at a relatively rapid physiological loading rate. Six immature and six mature humeral head specimens were loaded against glass and strains were measured at the end of a 5s loading ramp on the textured articular cross-section using digital image correlation analysis. The primary findings indicate that elevated tensile and compressive strains occur near the articular surface, around the center of the contact region. Few qualitative or quantitative differences were observed between mature and immature joints. Under an average contact stress of approximately 1.7 MPa, the peak compressive strains averaged -0.131+/-0.048, which was significantly less than the relative change in cartilage thickness, -0.104+/-0.032 (p<0.05). The peak tensile strains were significantly smaller in magnitude, at 0.0325+/-0.013. These experimental findings differ from a previous finite element analysis of articular contact, which predicted peak strains at the cartilage-bone interface even when accounting for the porous-hydrated nature of the tissue, its depth-dependent inhomogeneity, and the disparity between its tensile and compressive properties. These experimental results yield new insights into the local mechanical environment of the tissue and cells, and suggest that further refinements are needed in the modeling of contacting articular layers.     

Extraction of mechanical properties of articular cartilage from osmotic swelling behavior monitored using high frequency ultrasound

Articular cartilage is a biological weight-bearing tissue covering the bony ends of articulating joints. Negatively charged proteoglycan (PG) in articular cartilage is one of the main factors that govern its compressive mechanical behavior and swelling phenomenon. PG is nonuniformly distributed throughout the depth direction, and its amount or distribution may change in the degenerated articular cartilage such as osteoarthritis. In this paper, we used a 50 MHz ultrasound system to study the depth-dependent strain of articular cartilage under the osmotic loading induced by the decrease of the bathing saline concentration. The swelling-induced strains under the osmotic loading were used to determine the layered material properties of articular cartilage based on a triphasic model of the free-swelling. Fourteen cylindrical cartilage-bone samples prepared from fresh normal bovine patellae were tested in situ in this study. A layered triphasic model was proposed to describe the depth distribution of the swelling strain for the cartilage and to determine its aggregate modulus H(a) at two different layers, within which H(a) was assumed to be linearly dependent on the depth. The results showed that H(a) was 3.0+/-3.2, 7.0+/-7.4, 24.5+/-11.1 MPa at the cartilage surface, layer interface, and deep region, respectively. They are significantly different (p<0.01). The layer interface located at 70%+/-20% of the overall thickness from the uncalcified-calcified cartilage interface. Parametric analysis demonstrated that the depth-dependent distribution of the water fraction had a significant effect on the modeling results but not the fixed charge density. This study showed that high-frequency ultrasound measurement together with triphasic modeling is practical for quantifying the layered mechanical properties of articular cartilage nondestructively and has the potential for providing useful information for the detection of the early signs of osteoarthritis.     

Ultrasonic measurement of depth-dependent transient behaviors of articular cartilage under compression

We previously reported an ultrasound method for measuring the depth-dependent equilibrium mechanical properties of articular cartilage using quasi-static compression. The objective of this paper was to introduce our recent development for nondestructively measuring the transient depth-dependent strains of full-thickness articular cartilage specimens prepared from bovine patellae. A 50 MHz focused ultrasound transducer was used to collect ultrasound echoes from articular cartilage specimens (n=8) and sponge phantoms with open pores (n=10) during tests of compression and subsequent stress-relaxation. The transient displacements of the tissues at different depths along the compression direction were calculated from the ultrasound echoes using a cross-correlation tracking technique. An LVDT sensor and a load cell were used to measure the overall deformation of the tissue and the applied force, respectively. Results showed that the tissues inside the cartilage layer continued to move during the stress-relaxation phase after the compression was completed. In the equilibrium state, the displacements of the cartilage tissues at the depths of 1/4, 1/2, and 3/4 of the full-thickness reduced by 51%+/-22%, 54%+/-17%, and 50+/-17%, respectively, in comparison with its peak value. However, no similar phenomenon was observed in the sponge phantoms. Our preliminary results demonstrated that this ultrasound method may provide a potential tool for the nondestructive measurement of the transient depth-dependent processes involved in biological and bioengineered soft tissues as well as soft biomaterials under dynamic loading.     

Anisotropic strain-dependent material properties of bovine articular cartilage in the transitional range from tension to compression

Articular cartilage exhibits complex mechanical properties such as anisotropy, inhomogeneity and tension-compression nonlinearity. This study proposes and demonstrates that the application of compressive loading in the presence of osmotic swelling can be used to acquire a spectrum of incremental cartilage moduli (EYi) and Poisson's ratios (upsilon ij) from tension to compression. Furthermore, the anisotropy of the tissue can be characterized in both tension and compression by conducting these experiments along three mutually perpendicular loading directions: parallel to split-line (1-direction), perpendicular to split-line (2-direction) and along the depth direction (3-direction, perpendicular to articular surface), accounting for tissue inhomogeneity between the surface and deep layers in the latter direction. Tensile moduli were found to be strain-dependent while compressive moduli were nearly constant. The peak tensile (+) Young's moduli in 0.15M NaCl were E+Y1=3.1+/-2.3, E+Y2=1.3+/-0.3, E+Y3(Surface)=0.65+/-0.29 and E+Y3(Deep)=2.1+/-1.2 MPa. The corresponding compressive (-) Young's moduli were E-Y1=0.23+/-0.07, E-Y2=0.22+/-0.07, E-Y3(Surface)=0.18+/-0.07 and E-Y3(Deep)=0.35+/-0.11 MPa. Peak tensile Poisson's ratios were upsilon+12=0.22+/-0.06, upsilon+21=0.13+/-0.07, upsilon+31(Surface)=0.10+/-0.03 and upsilon+31(Deep)=0.20+/-0.05 while compressive Poisson's ratios were upsilon-12=0.027+/-0.012, upsilon-21=0.017+/-0.07, upsilon-31(Surface)=0.034+/-0.009 and upsilon-31(Deep)=0.065+/-0.024. Similar measurements were also performed at 0.015 M and 2 M NaCl, showing strong variations with ionic strength. Results indicate that (a) a smooth transition occurs in the stress-strain and modulus-strain responses between the tensile and compressive regimes, and (b) cartilage exhibits orthotropic symmetry within the framework of tension-compression nonlinearity. The strain-softening behavior of cartilage (the initial decrease in EYi with increasing compressive strain) can be interpreted in the context of osmotic swelling and tension-compression nonlinearity.     

An automated approach for direct measurement of two-dimensional strain distributions within articular cartilage under unconfined compression

An automated approachfor measuring in situ two-dimensional strain fields was developed and validated for its application to cartilage mechanics. This approach combines video microscopy, optimized digital image correlation (DIC), thin-plate spline smoothing (TPSS) and generalized cross-validation (GCV) techniques to achieve the desired efficiency and accuracy. Results demonstrate that sub-pixel accuracies can be achieved for measuring tissue displacements with this methodology with a measurement uncertainty ranging from 0.25 to 0.30 pixels. The deformational gradients (from which the strains are determined) can be evaluated directly using the optimized DIC, with a measurement uncertainty of 0.017 to approximately 0.032. In actual measurements of strain in cartilage, TPSS and differentiation can be used to achieve a more accurate measurement of the gradients from the displacement data. Using this automated approach, the two-dimensional strain fields inside immature bovine carpometacarpal joint cartilage specimens under unconfined compression were characterized (n=21). The depth-dependent apparent elastic modulus and Poisson's ratio were also determined and found to be smallest at the articular surface and increasing with depth. The apparent Poisson's ratio is found to decrease with increasing compressive strain, with values as low as 0.01 observed near the articular surface at 25% compression. The variation of the apparent Poisson's ratio with depth is found to be consistent with a theoretical model of cartilage which accounts for the disparity in its tensile and compressive moduli.     

Finite element study of a tissue-engineered cartilage transplant in human tibiofemoral joint

Most tissue-engineered cartilage constructs are more compliant than native articular cartilage (AC) and are poorly integrated to the surrounding tissue. To investigate the effect of an implanted tissue-engineered construct (TEC) with these inferior properties on the mechanical environment of both the engineered and adjacent native tissues, a finite element study was conducted. Biphasic swelling was used to model tibial cartilage and an implanted TEC with the material properties of either native tissue or a decreased elastic modulus and fixed charged density. Creep loading was applied with a rigid impermeable indenter that represented the femur. In comparison with an intact joint, compressive strains in the transplant, surface contact stress in the adjacent native AC and load partitioning between different phases of cartilage were affected by inferior properties of TEC. Results of this study may lead to a better understanding of the complex mechanical environment of an implanted TEC.     

Finite element study of a tissue-engineered cartilage transplant in human tibiofemoral joint

Most tissue-engineered cartilage constructs are more compliant than native articular cartilage (AC) and are poorly integrated to the surrounding tissue. To investigate the effect of an implanted tissue-engineered construct (TEC) with these inferior properties on the mechanical environment of both the engineered and adjacent native tissues, a finite element study was conducted. Biphasic swelling was used to model tibial cartilage and an implanted TEC with the material properties of either native tissue or a decreased elastic modulus and fixed charged density. Creep loading was applied with a rigid impermeable indenter that represented the femur. In comparison with an intact joint, compressive strains in the transplant, surface contact stress in the adjacent native AC and load partitioning between different phases of cartilage were affected by inferior properties of TEC. Results of this study may lead to a better understanding of the complex mechanical environment of an implanted TEC.     

Differential regulation of immature articular cartilage compressive moduli and Poisson’s ratios by in vitro stimulation with IGF-1 and TGF-β1

Mechanisms of articular cartilage growth and maturation have been elucidated by studying composition-function dynamics during in vivo development and in vitro culture with stimuli such as insulin-like growth factor-1 (IGF-1) and transforming growth factor-beta 1 (TGF-β1). This study tested the hypothesis that IGF-1 and TGF-β1 regulate immature cartilage compressive moduli and Poisson’s ratios in a manner consistent with known effects on tensile properties. Bovine calf articular cartilage from superficial-articular (S) and middle-growth (M) regions were analyzed fresh or following culture in medium with IGF-1 or TGF-β1. Mechanical properties in confined (CC) and unconfined (UCC) compression, cartilage matrix composition, and explant size were assessed. Culture with IGF-1 resulted in softening in CC and UCC, increased Poisson’s ratios, substantially increased tissue volume, and accumulation of glycosaminoglycan (GAG) and collagen (COL). Culture with TGF-β1 promoted maturational changes in the S layer, including stiffening in CC and UCC and increased concentrations of GAG, COL, and pyridinoline crosslinks (PYR), but little growth. Culture of M layer explants with TGF-β1 was nearly homeostatic. Across treatment groups, compressive moduli in CC and UCC were positively related to GAG, COL, and PYR concentrations, while Poisson’s ratios were negatively related to concentrations of these matrix components. Thus, IGF-1 and TGF-β1 differentially regulate the compressive mechanical properties and size of immature articular cartilage in vitro. Prescribing tissue growth, maturation, or homeostasis by controlling the in vitro biochemical environment with such growth factors may have applications in cartilage repair and tissue engineering.     

In vitro articular cartilage growth with sequential application of IGF-1 and TGF-β1 enhances volumetric growth and maintains compressive properties

In vitro cultures with insulin-like growth factor-1 (IGF-1) and transforming growth factor-β1 (TGF-β1) have previously been shown to differentially modulate the growth of immature bovine articular cartilage. IGF-1 stimulates expansive growth yet decreases compressive moduli and increases compressive Poisson's ratios, whereas TGF-β1 maintains tissue size, increases compressive moduli, and decreases compressive Poisson's ratios. The current study's hypothesis was that sequential application of IGF-1 and TGF-β1 during in vitro culture produces geometric and compressive mechanical properties that lie between extreme values produced when using either growth factor alone. Immature bovine articular cartilage specimens were harvested and either untreated (D0, i.e., day zero) or cultured in vitro for either 6 days with IGF-1 (D6 IGF), 12 days with IGF-1 (D12 IGF), or 6 days with IGF-1 followed by 6 days with TGF-β1 (D12 SEQ, i.e., sequential). Following treatment, all specimens were tested for geometric, biochemical, and compressive mechanical properties. Relative to D0, D12 SEQ treatment enhanced volumetric growth, but to a lower value than that for D12 IGF. Furthermore, D12 SEQ treatment maintained compressive moduli and Poisson's ratios at values higher and lower, respectively, than those for D12 IGF. Considering the previously described effects of 12 days of treatment with TGF-β1 alone, D12 SEQ induced both growth and mechanical property changes between those produced with either IGF-1 or TGF-β1 alone. The results suggest that it may be possible to vary the durations of select growth factors, including IGF-1 and TGF-β1, to more precisely modulate the geometric, biochemical, and mechanical properties of immature cartilage graft tissue in clinical repair strategies.     

Functional tissue engineering of articular cartilage through dynamic loading of chondrocyte-seeded agarose gels

Due to its avascular nature, articular cartilage exhibits a very limited capacity to regenerate and to repair. Although much of the tissue-engineered cartilage in existence has been successful in mimicking the morphological and biochemical appearance of hyaline cartilage, it is generally mechanically inferior to the natural tissue. In this study, we tested the hypothesis that the application of dynamic deformational loading at physiological strain levels enhances chondrocyte matrix elaboration in cell-seeded agarose scaffolds to produce a more functional engineered tissue construct than in free swelling controls. A custom-designed bioreactor was used to load cell-seeded agarose disks dynamically in unconfined compression with a peak-to-peak compressive strain amplitude of 10 percent, at a frequency of 1 Hz, 3 x (1 hour on, 1 hour off)/day, 5 days/week for 4 weeks. Results demonstrated that dynamically loaded disks yielded a sixfold increase in the equilibrium aggregate modulus over free swelling controls after 28 days of loading (100 +/- 16 kPa versus 15 +/- 8 kPa, p < 0.0001). This represented a 21-fold increase over the equilibrium modulus of day 0 (4.8 +/- 2.3 kPa). Sulfated glycosaminoglycan content and hydroxyproline content was also found to be greater in dynamically loaded disks compared to free swelling controls at day 21 (p < 0.0001 and p = 0.002, respectively).     

Anisotropy, inhomogeneity, and tension-compression nonlinearity of human glenohumeral cartilage in finite deformation

The tensile and compressive properties of human glenohumeral cartilage were determined by testing 120 rectangular strips in uniaxial tension and 70 cylindrical plugs in confined compression, obtained from five human glenohumeral joints. Specimens were harvested from five regions across the articular surface of the humeral head and two regions on the glenoid. Tensile strips were obtained along two orientations, parallel and perpendicular to the split-line directions. Two serial slices through the thickness, corresponding to the superficial and middle zones of the cartilage layers, were prepared from each tensile strip and each compressive plug. The equilibrium tensile modulus and compressive aggregate modulus of cartilage were determined from the uniaxial tensile and confined compression tests, respectively. Significant differences in the tensile moduli were found with depth and orientation relative to the local split-line direction. Articular cartilage of the humeral head was significantly stiffer in tension than that of the glenoid. There were significant differences in the aggregate compressive moduli of articular cartilage between superficial and middle zones in the humeral head. Furthermore, tensile and compressive stress-strain responses exhibited nonlinearity under finite strain, while the tensile modulus differed by up to two orders of magnitude from the compressive aggregate modulus at 0% strain, demonstrating a high degree of tension-compression nonlinearity. The complexity of the mechanical properties of human glenohumeral cartilage was exposed in this study, showing anisotropy, inhomogeneity, and tension-compression nonlinearity within the same joint. The observed differences in the tensile properties of human glenohumeral cartilage suggest that the glenoid may be more susceptible to cartilage degeneration than the humeral head.     

Comparison of nonlinear mechanical properties of bovine articular cartilage and meniscus

Nonlinear, linear and failure properties of articular cartilage and meniscus in opposing contact surfaces are poorly known in tension. Relationships between the tensile properties of articular cartilage and meniscus in contact with each other within knee joints are also not known. In the present study, rectangular samples were prepared from the superficial lateral femoral condyle cartilage and lateral meniscus of bovine knee joints. Tensile tests were carried out with a loading rate of 5 mm/min until the tissue rupture. Nonlinear properties of the toe region, linear properties in larger strains, and failure properties of both tissues were analysed. The strain-dependent tensile modulus of the toe region, Young's modulus of the linear region, ultimate tensile stress and toughness were on average 98.2, 8.3, 4.0 and 1.9 times greater (p<0.05) for meniscus than for articular cartilage. In contrast, the toe region strain, yield strain and failure strain were on average 9.4, 3.1 and 2.3 times greater (p<0.05) for cartilage than for meniscus. There was a significant negative correlation between the strain-dependent tensile moduli of meniscus and articular cartilage samples within the same joints (r=−0.690, p=0.014). In conclusion, the meniscus possesses higher nonlinear and linear elastic stiffness and energy absorption capability before rupture than contacting articular cartilage, while cartilage has longer nonlinear region and can withstand greater strains before failure. These findings point out different load carrying demands that both articular cartilage and meniscus have to fulfil during normal physiological loading activities of knee joints.     

Acoustic properties of articular cartilage under mechanical stress

Mechano-acoustic and elastographic techniques may provide quantitative means for the in vivo diagnostics of articular cartilage. These techniques assume that sound speed does not change during tissue loading. As articular cartilage shows volumetric changes during compression, acoustic properties of cartilage may change affecting the validity of mechano-acoustic measurements. In this study, we examined the ultrasound propagation through human, bovine and porcine articular cartilage during stress-relaxation in unconfined compression. The time of flight (TOF) technique with known cartilage thickness (true sound speed) as well as in situ calibration method [Suh, Youn, Fu, J. Biomech. 34 (2001), 1347-1353] were used for the determination of sound speed. Ultrasound speed and attenuation decreased in articular cartilage during ramp compression, but returned towards the level of original values during relaxation. Variations in ultrasound speed induced an error in strain and compressive moduli provided that constant ultrasound speed and time-of-flight data was used to determine the tissue thickness. Highest errors in strain (-11.8 +/- 12.0%) and dynamic modulus (15.4 +/- 17.9%) were recorded in bovine cartilage. TOF and in situ calibration methods yielded different results for changes in sound speed during compression. We speculate that the variations in acoustic properties in loaded cartilage are related to rearrangement of the interstitial matrix, especially to that of collagen fibers. In human cartilage the changes, are, however relatively small and, according to the numerical simulations, mechano-acoustic techniques that assume constant acoustic properties for the cartilage will not be significantly impaired by this phenomenon.     

Nonlinear tensile properties of bovine articular cartilage and their variation with age and depth

Tensile stiffness of articular cartilage is much greater than its compressive stiffness and plays an essential role even in compressive properties by increasing transient fluid pressures during physiological loading. Recent studies of nonlinear properties of articular cartilage in compression revealed several physiologically pertinent nonlinear behaviors, all of which required that cartilage tensile stiffness increase significantly with stretch. We therefore performed sequences of uniaxial tension tests on fresh bovine articular cartilage slices using a protocol that allowed several hours to attain equilibrium and measured longitudinal and transverse tissue strain. By testing bovine cartilage from different ages (6 months to 6 years) we found that equilibrium and transient tensile modulus increased significantly with maturation and age, from 0 to 15 MPa at equilibrium and from 10 to 28 MPa transiently. Our results indicate that cartilage stiffens with age in a manner similar to other highly hydrated connective tissues, possibly due to age-dependent content of enzymatic and nonenzymatic collagen cross links. The long relaxation period used in our tests (5-10 hours) was necessary in order to attain equilibrium and avoid a very significant overestimation of equilibrium modulus that occurs when much shorter times are used (15-30 minutes). We also found that equilibrium and transient tensile modulus increased nonlinearly when cartilage is stretched from 0 to 10% strain without any previous tare load. Although our results estimate a nonlinear increase in tensile stiffness with stretch that is an order of magnitude lower than that required to predict nonlinear properties in compression, they are in agreement with previous results from other uniaxial tension tests of collagenous materials. We therefore speculate that biaxial tensile moduli may be much higher and thereby more compatible with observed nonlinear compressive properties.     

Cartilage interstitial fluid load support in unconfined compression

Under physiological conditions of loading, articular cartilage is subjected to both compressive strains, normal to the articular surface, and tensile strains, tangential to the articular surface. Previous studies have shown that articular cartilage exhibits a much higher modulus in tension than in compression, and theoretical analyses have suggested that this tension–compression nonlinearity enhances the magnitude of interstitial fluid pressurization during loading in unconfined compression, above a theoretical threshold of 33% of the average applied stress. The first hypothesis of this experimental study is that the peak fluid load support in unconfined compression is significantly greater than the 33% theoretical limit predicted for porous permeable tissues modeled with equal moduli in tension and compression. The second hypothesis is that the peak fluid load support is higher at the articular surface side of the tissue samples than near the deep zone, because the disparity between the tensile and compressive moduli is greater at the surface zone. Ten human cartilage samples from six patellofemoral joints, and 10 bovine cartilage specimens from three calf patellofemoral joints were tested in unconfined compression. The peak fluid load support was measured at 79±11% and 69±15% at the articular surface and deep zone of human cartilage, respectively, and at 94±4% and 71±8% at the articular surface and deep zone of bovine calf cartilage, respectively. Statistical analyses confirmed both hypotheses of this study. These experimental results suggest that the tension–compression nonlinearity of cartilage is an essential functional property of the tissue which makes interstitial fluid pressurization the dominant mechanism of load support in articular cartilage.     

Hyaluronic Acid Enhances the Mechanical Properties of Tissue-Engineered Cartilage Constructs

There is a need for materials that are well suited for cartilage tissue engineering. Hydrogels have emerged as promising biomaterials for cartilage repair, since, like cartilage, they have high water content, and they allow cells to be encapsulated within the material in a genuinely three-dimensional microenvironment. In this study, we investigated the mechanical properties of tissue-engineered cartilage constructs using in vitro culture models incorporating human chondrocytes from osteoarthritis patients. We evaluated hydrogels formed from mixtures of photocrosslinkable gelatin-methacrylamide (Gel-MA) and varying concentrations (0–2%) of hyaluronic acid methacrylate (HA-MA). Initially, only small differences in the stiffness of each hydrogel existed. After 4 weeks of culture, and to a greater extent 8 weeks of culture, HA-MA had striking and concentration dependent impact on the changes in mechanical properties. For example, the initial compressive moduli of cell-laden constructs with 0 and 1% HA-MA were 29 and 41 kPa, respectively. After 8 weeks of culture, the moduli of these constructs had increased to 66 and 147 kPa respectively, representing a net improvement of 69 kPa for gels with 1% HA-MA. Similarly the equilibrium modulus, dynamic modulus, failure strength and failure strain were all improved in constructs containing HA-MA. Differences in mechanical properties did not correlate with glycosaminoglycan content, which did not vary greatly between groups, yet there were clear differences in aggrecan intensity and distribution as assessed using immunostaining. Based on the functional development with time in culture using human chondrocytes, mixtures of Gel-MA and HA-MA are promising candidates for cartilage tissue-engineering applications.     

The importance of superficial collagen fibrils for the function of articular cartilage

It has been proposed that the superficial tangential zone (STZ) of articular cartilage is essential to the tissue’s load-distributing function. However, the exact mechanism by which the STZ fulfills this function has not yet been revealed. Using a channel-indentation experiment, it was recently shown that compared to intact tissue, cartilage without STZ behaves slightly stiffer and deforms significantly different in regions adjacent to mechanically compressed areas (Bevill et al. in Osteoarthr Cartil 18:1310–1318, 2010). We aim to further explore the role of STZ in the load-transfer mechanism of AC by thorough biomechanical analysis of these experiments. Using our previously validated fibril-reinforced swelling model of articular cartilage, which accounts for the depth-dependent collagen structure and biochemical composition of articular cartilage, we simulated the above-mentioned channel-indenter compression experiments for both intact and STZ-removed cartilage. First, we show that the composition of the deep zone in cartilage is most effective in carrying cartilage compression, which explains the apparent tissue stiffening after STZ removal. Second, we show that tangential fibrils in the STZ are responsible for transferring compressive loads from directly loaded regions to adjacent tissue. Cartilage with an intact STZ has superior load-bearing properties compared to cartilage in which the STZ is compromised, because the STZ is able to recruit a larger area of deep zone cartilage to carry compressive loads.     

Mapping the depth dependence of shear properties in articular cartilage

Determining the depth dependence of the shear properties of articular cartilage is essential for understanding the structure-function relation in this tissue. Here, we measured spatial variations in the shear modulus G of bovine articular cartilage using a novel technique that combines shear testing, confocal imaging and force measurement. We found that G varied by up to two orders of magnitude across a single sample, exhibited a global minimum 50-250 microm below the articular surface in a region just below the superficial zone and was roughly constant at depths > 1000 microm (the "plateau region"). For plateau strains gamma(plateau) approximately 0.75% and overall compressive strains epsilon approximately 5%, G(min) and G(plateau) were approximately 70 and approximately 650 kPa, respectively. In addition, we found that the shear modulus profile depended strongly on the applied shear and axial strains. The greatest change in G occurred at the global minimum where the tissue was highly nonlinear, stiffening under increased shear strain, and weakening under increased compressive strain. Our results can be explained through a simple thought model describing the observed nonlinear behavior in terms of localized buckling of collagen fibers and suggest that compression may decrease the vulnerability of articular cartilage to shear-induced damage by lowering the effective strain on individual collagen fibrils.