Microstructural Modeling of Collagen Network Mechanics and Interactions with the Proteoglycan Gel in Articular Cartilage

Cartilage matrix mechanical function is largely determined by interactions between the collagen fibrillar network and the proteoglycan gel. Although the molecular physics of these matrix constituents have been characterized and modern imaging methods are capable of localized measurement of molecular densities and orientation distributions, theoretical tools for using this information for prediction of cartilage mechanical behavior are lacking. We introduce a means to model collagen network contributions to cartilage mechanics based upon accessible microstructural information (fibril density and orientation distributions) and which self-consistently follows changes in microstructural geometry with matrix deformations. The interplay between the molecular physics of the collagen network and the proteoglycan gel is scaled up to determine matrix material properties, with features such as collagen fibril pre-stress in free-swelling cartilage emerging naturally and without introduction of ad hoc parameters. Methods are developed for theoretical treatment of the collagen network as a continuum-like distribution of fibrils, such that mechanical analysis of the network may be simplified by consideration of the spherical harmonic components of functions of the fibril orientation, strain, and stress distributions. Expressions for the collagen network contributions to matrix stress and stiffness tensors are derived, illustrating that only spherical harmonic components of orders 0 and 2 contribute to the stress, while orders 0, 2, and 4 contribute to the stiffness. Depth- and compression-dependent equilibrium mechanical properties of cartilage matrix are modeled, and advantages of the approach are illustrated by exploration of orientation and strain distributions of collagen fibrils in compressed cartilage. Results highlight collagen-proteoglycan interactions, especially for very small physiological strains where experimental data are relatively sparse. These methods for determining matrix mechanical properties from measurable quantities at the microscale (composition, structure, and molecular physics) may be useful for investigating cartilage structure-function relationships relevant to load-bearing, injury, and repair.     

How changes in interconnectivity affect the bulk properties of articular cartilage: a fibre network study

The remarkable compressive strength of articular cartilage arises from the mechanical interactions between the tension-resisting collagen fibrils and swelling proteoglycan proteins within the tissue. These interactions are facilitated by a significant level of interconnectivity between neighbouring collagen fibrils within the extracellular matrix. A reduction in interconnectivity is suspected to occur during the early stages of osteoarthritic degeneration. However, the relative contribution of these interconnections towards the bulk mechanical properties of articular cartilage has remained an open question. In this study, we present a simple 2D fibre network model which explicitly represents the microstructure of articular cartilage as collection of discrete nodes and linear springs. The transverse stiffness and swelling properties of this fibre network are studied, and a semi-analytic relationship which relates these two macroscopic properties via microscopic interconnectivity is derived. By comparing this derived expression to previously published experimental data, we show that although a reduction in network interconnectivity accounts for some of the observed changes in the mechanical properties of articular cartilage as degeneration occurs, a decrease in matrix interconnectivity alone do not provide a full account of this process.     

Importance of collagen orientation and depth-dependent fixed charge densities of cartilage on mechanical behavior of chondrocytes

The collagen network and proteoglycan matrix of articular cartilage are thought to play an important role in controlling the stresses and strains in and around chondrocytes, in regulating the biosynthesis of the solid matrix, and consequently in maintaining the health of diarthrodial joints. Understanding the detailed effects of the mechanical environment of chondrocytes on cell behavior is therefore essential for the study of the development, adaptation, and degeneration of articular cartilage. Recent progress in macroscopic models has improved our understanding of depth-dependent properties of cartilage. However, none of the previous works considered the effect of realistic collagen orientation or depth-dependent negative charges in microscopic models of chondrocyte mechanics. The aim of this study was to investigate the effects of the collagen network and fixed charge densities of cartilage on the mechanical environment of the chondrocytes in a depth-dependent manner. We developed an anisotropic, inhomogeneous, microstructural fibril-reinforced finite element model of articular cartilage for application in unconfined compression. The model consisted of the extracellular matrix and chondrocytes located in the superficial, middle, and deep zones. Chondrocytes were surrounded by a pericellular matrix and were assumed spherical prior to tissue swelling and load application. Material properties of the chondrocytes, pericellular matrix, and extracellular matrix were obtained from the literature. The loading protocol included a free swelling step followed by a stress-relaxation step. Results from traditional isotropic and transversely isotropic biphasic models were used for comparison with predictions from the current model. In the superficial zone, cell shapes changed from rounded to elliptic after free swelling. The stresses and strains as well as fluid flow in cells were greatly affected by the modulus of the collagen network. The fixed charge density of the chondrocytes, pericellular matrix, and extracellular matrix primarily affected the aspect ratios (height/width) and the solid matrix stresses of cells. The mechanical responses of the cells were strongly location and time dependent. The current model highlights that the collagen orientation and the depth-dependent negative fixed charge densities of articular cartilage have a great effect in modulating the mechanical environment in the vicinity of chondrocytes, and it provides an important improvement over earlier models in describing the possible pathways from loading of articular cartilage to the mechanical and biological responses of chondrocytes.     

On the anisotropy and inhomogeneity of permeability in articular cartilage

Articular cartilage is known to be anisotropic and inhomogeneous because of its microstructure. In particular, its elastic properties are influenced by the arrangement of the collagen fibres, which are orthogonal to the bone-cartilage interface in the deep zone, randomly oriented in the middle zone, and parallel to the surface in the superficial zone. In past studies, cartilage permeability has been related directly to the orientation of the glycosaminoglycan chains attached to the proteoglycans which constitute the tissue matrix. These studies predicted permeability to be isotropic in the undeformed configuration, and anisotropic under compression. They neglected tissue anisotropy caused by the collagen network. However, magnetic resonance studies suggest that fluid flow is "directed" by collagen fibres in biological tissues. Therefore, the aim of this study was to express the permeability of cartilage accounting for the microstructural anisotropy and inhomogeneity caused by the collagen fibres. Permeability is predicted to be anisotropic and inhomogeneous, independent of the state of strain, which is consistent with the morphology of the tissue. Looking at the local anisotropy of permeability, we may infer that the arrangement of the collagen fibre network plays an important role in directing fluid flow to optimise tissue functioning.     

The correspondence between equilibrium biphasic and triphasic material properties in mixture models of articular cartilage

Mixture models have been successfully used to describe the response of articular cartilage to various loading conditions. Mow et al. (J. Biomech. Eng. 102 (1980) 73) formulated a biphasic mixture model of articular cartilage where the collagen-proteoglycan matrix is modeled as an intrinsically incompressible porous-permeable solid matrix, and the interstitial fluid is modeled as an incompressible fluid. Lai et al. (J. Biomech. Eng. 113 (1991) 245) proposed a triphasic model of articular cartilage as an extension of their biphasic theory, where negatively charged proteoglycans are modeled to be fixed to the solid matrix, and monovalent ions in the interstitial fluid are modeled as additional fluid phases. Since both models co-exist in the cartilage literature, it is useful to show how the measured properties of articular cartilage (the confined and unconfined compressive and tensile moduli, the compressive and tensile Poisson's ratios, and the shear modulus) relate to both theories. In this study, closed-form expressions are presented that relate biphasic and triphasic material properties in tension, compression and shear. These expressions are then compared to experimental findings in the literature to provide greater insight into the measured properties of articular cartilage as a function of bathing solutions salt concentrations and proteoglycan fixed-charge density.     

Collagen in tissue-engineered cartilage: Types,structure, and crosslinks

The function of articular cartilage as a weight-bearing tissue depends on the specific arrangement of collagen types II and IX into a three-dimensional organized collagen network that can balance the swelling pressure of the proteoglycan/ water gel. To determine whether cartilage engineered in vitro contains a functional collagen network, chondrocyte-polymer constructs were cultured for up to 6 weeks and analyzed with respect to the composition and ultrastructure of collagen by using biochemical and immunochemical methods and scanning electron microscopy. Total collagen content and the concentration of pyridinium crosslinks were significantly (57% and 70%, respectively) lower in tissue-engineered cartilage that in bovine calf articular cartilage. However, the fractions of collagen types II, IX, and X and the collagen network organization, density, and fibril diameter in engineered cartilage were not significantly different from those in natural articular cartilage. The implications of these findings for the field of tissue engineering are that differentiated chondrocytes are capable of forming a complex structure of collagen matrix in vitro, producing a tissue similar to natural articular cartilage on an ultrastructural scale. J. Cell. Biochem. 71:313–327, 1998. © 1998 Wiley-Liss, Inc.     

Engineering superficial zone features in tissue engineered cartilage

A major challenge in cartilage tissue engineering is the need to recreate the native tissue's anisotropic extracellular matrix structure. This anisotropy has important mechanical and biological consequences and could be crucial for integrative repair. Here, we report that hydrodynamic conditions that mimic the motion‐induced flow fields in between the articular surfaces in the synovial joint induce the formation of a distinct superficial layer in tissue engineered cartilage hydrogels, with enhanced production of cartilage matrix proteoglycan and Type II collagen. Moreover, the flow stimulation at the surface induces the production of the surface zone protein Proteoglycan 4 (aka PRG4 or lubricin). Analysis of second harmonic generation signature of collagen in this superficial layer reveals a highly aligned fibrillar matrix that resembles the alignment pattern in native tissue's surface zone, suggesting that mimicking synovial fluid flow at the cartilage surface in hydrodynamic bioreactors could be key to creating engineered cartilage with superficial zone features. Biotechnol. Bioeng. 2013; 110: 1476–1486. © 2012 Wiley Periodicals, Inc.     

力学环境对软骨基质代谢的影响

正常关节软骨所受压力是由动态压力与静态压力交替完成。压力引起软骨一系列生理变化包括细胞及细胞外基质成分变形、组织内液体流动、水流电位和生理生化变化。这些变化直接调控细胞外基质代谢。体外构建有良好功能的组织工程化软骨是目前软骨病变、缺损理想的修复方法。研究力学环境对软骨基质代谢的影响,对构建组织工程化软骨有深远意义。     

The role of viscoelasticity of collagen fibers in articular cartilage: theory and numerical formulation

The relative importance of fluid-dependent and fluid-independent transient mechanical behavior in articular cartilage was examined for tensile and unconfined compression testing using a fibril reinforced model. The collagen matrix of articular cartilage was modeled as viscoelastic using a quasi-linear viscoelastic formulation with strain-dependent elastic modulus, while the proteoglycan matrix was considered as linearly elastic. The collagen viscoelastic properties were obtained by fitting experimental data from a tensile test. These properties were used to investigate unconfined compression testing, and the sensitivity of the properties was also explored. It was predicted that the stress relaxation observed in tensile tests was not caused by fluid pressurization at the macroscopic level. A multi-step tensile stress relaxation test could be approximated using a hereditary integral in which the elastic fibrillar modulus was taken to be a linear function of the fibrillar strain. Applying the same formulation to the radial fibers in unconfined compression, stress relaxation could not be simulated if fluid pressurization were absent. Collagen viscoelasticity was found to slightly weaken fluid pressurization in unconfined compression, and this effect was relatively more significant at moderate strain rates. Therefore, collagen viscoelasticity appears to play an import role in articular cartilage in tensile testing, while fluid pressurization dominates the transient mechanical behavior in compression. Collagen viscoelasticity plays a minor role in the mechanical response of cartilage in unconfined compression if significant fluid flow is present.     

The effects of selective matrix degradation on the short-term compressive properties of adult human articular cartilage.

The effects of proteoglycan and collagen digestion on the transient response of human articular cartilage when tested in unconfined compression were determined. Small cylindrical specimens of cartilage, isolated from the femoral head of the hip joint and from the femoral condyles of the knee joint, were subjected to a suddenly applied compressive load using a test apparatus designed to yield a transient oscillatory response. From this response values of the elastic stiffness and the viscous damping coefficient were determined. Cathepsin D and cathepsin B1 were used to digest the proteoglycan in some specimens, while in other specimens leukocyte elastase was used to attack the non-helical terminal regions of the Type II tropocollagen molecules and possibly the Type IX collagen molecule and thereby disturb the integrity of the collagen mesh. The results showed that proteoglycan digestion alone reduced the viscous damping coefficient but it did not significantly alter the elastic stiffness as determined from the oscillatory response. In contrast, the action of elastase reduced both the damping coefficient and the elastic stiffness of the cartilage. The results demonstrated the role of proteoglycans in regulating fluid transport in cartilage and hence controlling the time-dependent viscous properties. The elastic stiffness was shown to be dependent on the integrity of the collagen fibre network and not on the proteoglycans.     

A microstructurally based continuum model of cartilage viscoelasticity and permeability incorporating measured statistical fiber orientations

The remarkable mechanical properties of cartilage derive from an interplay of isotropically distributed, densely packed and negatively charged proteoglycans; a highly anisotropic and inhomogeneously oriented fiber network of collagens; and an interstitial electrolytic fluid. We propose a new 3D finite strain constitutive model capable of simultaneously addressing both solid (reinforcement) and fluid (permeability) dependence of the tissue’s mechanical response on the patient-specific collagen fiber network. To represent fiber reinforcement, we integrate the strain energies of single collagen fibers—weighted by an orientation distribution function (ODF) defined over a unit sphere—over the distributed fiber orientations in 3D. We define the anisotropic intrinsic permeability of the tissue with a structure tensor based again on the integration of the local ODF over all spatial fiber orientations. By design, our modeling formulation accepts structural data on patient-specific collagen fiber networks as determined via diffusion tensor MRI. We implement our new model in 3D large strain finite elements and study the distributions of interstitial fluid pressure, fluid pressure load support and shear stress within a cartilage sample under indentation. Results show that the fiber network dramatically increases interstitial fluid pressure and focuses it near the surface. Inhomogeneity in the tissue’s composition also increases fluid pressure and reduces shear stress in the solid. Finally, a biphasic neo-Hookean material model, as is available in commercial finite element codes, does not capture important features of the intra-tissue response, e.g., distributions of interstitial fluid pressure and principal shear stress.     

Collagen Network of Articular Cartilage Modulates Fluid Flow and Mechanical Stresses in Chondrocyte

The extracellular matrix of articular cartilage modulates the mechanical signals sensed by the chondrocytes. In the present study, a finite element model (FEM) of the chondrocyte and its microenvironment was reconstructed using the information from fourier transform infrared imaging spectroscopy. This environment consisted of pericellular, territorial (mainly proteoglycans), and inter-territorial (mainly collagen) matrices. The chondrocyte, pericellular, and territorial matrix were assumedto be mechanically isotropic and poroelastic, whereas the inter-territorial matrix, due to its high collagen content, was assumed to be transversely isotropic and poroelastic. Under instantaneous strain-controlled compression, the FEM indicated that the fluid pressure within the chondrocyte increased nonlinearly as a function of the in-plane Young’s modulus of the collagen network. Under instantaneous force-controlled compression, the chondrocyte experienced the highest fluid pressure when the in-plane Young’s modulus of the collagen network was ~4 MPa. Based on the present results, the mechanical characteristics of the collagen network of articular cartilage can modify fluid flow and stresses in chondrocytes. Therefore, the integrity of the collagen network may be an important determinant in cell stimulation and in the control of the matrix maintenance.     

Elastic anisotropy of articular cartilage is associated with the microstructures of collagen fibers and chondrocytes

Chondrocyte shape and volumetric concentration change as a function of depth in articular cartilage. A given chondrocyte shape produces different effects on the global material properties depending on the structure of the collagen fiber network. The shape and volumetric concentration of chondrocytes in articular cartilage appear to be related to the mechanical stability of the matrix. The present study was aimed to investigate, theoretically, the effects of the structural arrangement of the collagen fiber network, and the shape and distribution of chondrocytes, on the global material behavior of articular cartilage. Articular cartilage was assumed to be a four-phasic composite comprised of a matrix (associated with the properties of the proteoglycan structure), vertically and horizontally distributed collagen fibers, and spheroidal inclusions representing chondrocytes. A solution for composite materials was used to estimate the global, effective material properties of cartilage. Only the elasticity of the solid phase was investigated in the present study. Our simulations suggest that a soft, spheroidal cell inclusion in a fiber-reinforced proteoglycan matrix affects the material properties differently depending on the shape of the spheroidal inclusions. If the long axis of the inclusions is parallel to the collagen fibers, as in the deep zone, the soft inclusions increase the stiffness of the composite in the fiber direction, and reduce the stiffness of the composite in the direction normal to the fibers. Furthermore, we found that Young's modulus normal to the contact surface increases from the superficial to the deep zone in articular cartilage by a factor of 10-50, a finding that agrees well with experimental observations. Our analysis suggests that the combination of proteoglycan matrix, fiber orientation, and shape of chondrocytes are intimately related and are likely adapted to optimize the mechanical stability and load carrying capacity of the structure.     

Interaction of cartilage proteoglycans with collagen-substituted agarose gels.

1. Rat tail-tendon collagen was coupled to activated Sepharose 4B at 2.5 mg of collagen/ml of gel. Chromatographic columns of this gel were calibrated with T2 virus (Vo) and Dnp-alanine (Vt). 2. The chromatographic behaviour of cartilage proteoglycans on the collagen-substituted gel was studied under conditions of varying ionic strength. Proteoglycan subunit obtained from bovine nasal cartilage, the proteoglycan obtained after digestion with chondroitnase ABC and purified chondriotin sulphate were all retarded on the collagen gel by an interaction that abolished at I0.17. Purified keratan sulphate and hyaluronic acid were not retarded. 3. A strong ionic interaction between cartilage proteoglycan and collagen was demonstrated to depend on the structure of the protein core of the proteoglycan.     

Fibril reinforced poroelastic model predicts specifically mechanical behavior of normal,proteoglycan depleted and collagen degraded articular cartilage

Degradation of collagen network and proteoglycan (PG) macromolecules are signs of articular cartilage degeneration. These changes impair cartilage mechanical function. Effects of collagen degradation and PG depletion on the time-dependent mechanical behavior of cartilage are different. In this study, numerical analyses, which take the compression-tension nonlinearity of the tissue into account, were carried out using a fibril reinforced poroelastic finite element model. The study aimed at improving our understanding of the stress-relaxation behavior of normal and degenerated cartilage in unconfined compression. PG and collagen degradations were simulated by decreasing the Young's modulus of the drained porous (nonfibrillar) matrix and the fibril network, respectively. Numerical analyses were compared to results from experimental tests with chondroitinase ABC (PG depletion) or collagenase (collagen degradation) digested samples. Fibril reinforced poroelastic model predicted the experimental behavior of cartilage after chondroitinase ABC digestion by a major decrease of the drained porous matrix modulus (-64+/-28%) and a minor decrease of the fibril network modulus (-11+/-9%). After collagenase digestion, in contrast, the numerical analyses predicted the experimental behavior of cartilage by a major decrease of the fibril network modulus (-69+/-5%) and a decrease of the drained porous matrix modulus (-44+/-18%). The reduction of the drained porous matrix modulus after collagenase digestion was consistent with the microscopically observed secondary PG loss from the tissue. The present results indicate that the fibril reinforced poroelastic model is able to predict specifically characteristic alterations in the stress-relaxation behavior of cartilage after enzymatic modifications of the tissue. We conclude that the compression-tension nonlinearity of the tissue is needed to capture realistically the mechanical behavior of normal and degenerated articular cartilage.     

Relative contribution of articular cartilage’s constitutive components to load support depending on strain rate

Cartilage is considered a biphasic material in which the solid is composed of proteoglycans and collagen. In biphasic tissue, the hydraulic pressure is believed to bear most of the load under higher strain rates and its dissipation due to fluid flow determines creep and relaxation behavior. In equilibrium, hydraulic pressure is zero and load bearing is transferred to the solid matrix. The viscoelasticity of the collagen network also contributes to its time-dependent behavior, and the osmotic pressure to load bearing in equilibrium. The aim of the present study was to determine the relative contributions of hydraulic pressure, viscoelastic collagen stress, solid matrix stiffness and osmotic pressure to load carriage in cartilage under transient and equilibrium conditions. Unconfined compression experiments were simulated using a fibril-reinforced poroviscoelastic model of articular cartilage, including water, fibrillar viscoelastic collagen and non-fibrillar charged glycosaminoglycans. The relative contributions of hydraulic and osmotic pressures and stresses in the fibrillar and non-fibrillar network were evaluated in the superficial, middle and deep zone of cartilage under five different strain rates and after relaxation. Initially upon loading, the hydraulic pressure carried most of the load in all three zones. The osmotic swelling pressure carried most of the equilibrium load. In the surface zone, where the fibers were loaded in tension, the collagen network carried 20 % of the load for all strain rates. The importance of these fibers was illustrated by artificially modifying the fiber architecture, which reduced the overall stiffness of cartilage in all conditions. In conclusion, although hydraulic pressure dominates the transient behavior during cartilage loading, due to its viscoelastic nature the superficial zone collagen fibers carry a substantial part of the load under transient conditions. This becomes increasingly important with higher strain rates. The interesting and striking new insight from this study suggests that under equilibrium conditions, the swelling pressure generated by the combination of proteoglycans and collagen reinforcement accounts cartilage stiffness for more than 90 % of the loads carried by articular cartilage. This finding is different from the common thought that load is transferred from fluid to solid and is carried by the aggregate modulus of the solid. Rather, it is transformed from hydraulic to osmotic swelling pressure. These results show the importance of considering both (viscoelastic) collagen fibers as well as swelling pressure in studies of the (transient) mechanical behavior of cartilage.     

Characterization of articular cartilage by combining microscopic analysis with a fibril-reinforced finite-element model

Load-bearing characteristics of articular cartilage are impaired during tissue degeneration. Quantitative microscopy enables in vitro investigation of cartilage structure but determination of tissue functional properties necessitates experimental mechanical testing. The fibril-reinforced poroviscoelastic (FRPVE) model has been used successfully for estimation of cartilage mechanical properties. The model includes realistic collagen network architecture, as shown by microscopic imaging techniques. The aim of the present study was to investigate the relationships between the cartilage proteoglycan (PG) and collagen content as assessed by quantitative microscopic findings, and model-based mechanical parameters of the tissue. Site-specific variation of the collagen network moduli, PG matrix modulus and permeability was analyzed. Cylindrical cartilage samples (n=22) were harvested from various sites of the bovine knee and shoulder joints. Collagen orientation, as quantitated by polarized light microscopy, was incorporated into the finite-element model. Stepwise stress-relaxation experiments in unconfined compression were conducted for the samples, and sample-specific models were fitted to the experimental data in order to determine values of the model parameters. For comparison, Fourier transform infrared imaging and digital densitometry were used for the determination of collagen and PG content in the same samples, respectively. The initial and strain-dependent fibril network moduli as well as the initial permeability correlated significantly with the tissue collagen content. The equilibrium Young's modulus of the nonfibrillar matrix and the strain dependency of permeability were significantly associated with the tissue PG content. The present study demonstrates that modern quantitative microscopic methods in combination with the FRPVE model are feasible methods to characterize the structure-function relationships of articular cartilage.     

Hip chondrolabral mechanics during activities of daily living: Role of the labrum and interstitial fluid pressurization

Osteoarthritis of the hip can result from mechanical factors, which can be studied using finite element (FE) analysis. FE studies of the hip often assume there is no significant loss of fluid pressurization in the articular cartilage during simulated activities and approximate the material as incompressible and elastic. This study examined the conditions under which interstitial fluid load support remains sustained during physiological motions, as well as the role of the labrum in maintaining fluid load support and the effect of its presence on the solid phase of the surrounding cartilage. We found that dynamic motions of gait and squatting maintained consistent fluid load support between cycles, while static single-leg stance experienced slight fluid depressurization with significant reduction of solid phase stress and strain. Presence of the labrum did not significantly influence fluid load support within the articular cartilage, but prevented deformation at the cartilage edge, leading to lower stress and strain conditions in the cartilage. A morphologically accurate representation of collagen fibril orientation through the thickness of the articular cartilage was not necessary to predict fluid load support. However, comparison with simplified fibril reinforcement underscored the physiological importance. The results of this study demonstrate that an elastic incompressible material approximation is reasonable for modeling a limited number of cyclic motions of gait and squatting without significant loss of accuracy, but is not appropriate for static motions or numerous repeated motions. Additionally, effects seen from removal of the labrum motivate evaluation of labral reattachment strategies in the context of labral repair.     

Role of pericellular matrix in development of a mechanically functional neocartilage

The role of the chondrocyte pericellular matrix (PCM) was examined in a three-dimensional chondrocyte culture system to determine whether retention of the native pericellular matrix could stimulate collagen and proteoglycan accumulation and also promote the formation of a mechanically functional hyaline-like neocartilage. Porcine chondrocytes and chondrons, consisting of the chondrocyte with its intact pericellular matrix, were maintained in pellet culture for up to 12 weeks. Sulfated glycosaminoclycans and type II collagen were measured biochemically. Immunocytochemistry was used to examine collagen localization as well as cell distribution within the pellets. In addition, the equilibrium compressive moduli of developing pellets were measured to determine whether matrix deposition contributed to the mechanical stiffness of the cartilage constructs. Pellets increased in size and weight over a 6-week period without apparent cell proliferation. Although chondrocytes quickly rebuilt a PCM rich in type VI collagen, chondron pellets accumulated significantly more proteoglycan and type II collagen than did chondrocyte pellets, indicating a greater positive effect of the native PCM. After 5 weeks in chondron pellets, matrix remodeling was evident by microscopy. Cells that had been uniformly distributed throughout the pellets began to cluster between large areas of interterritorial matrix rich in type II collagen. After 12 weeks, clusters were stacked in columns. A rapid increase in compressive strength was observed between 1 and 3 weeks in culture for both chondron and chondrocyte pellets and, by 6 weeks, both had achieved 25% of the equilibrium compressive stiffness of cartilage explants. Retention of the in vivo PCM during chondrocyte isolation promotes the formation of a mechanically functional neocartilage construct, suitable for modeling the responses of articular cartilage to chemical stimuli or mechanical compression.