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

Depth-Dependent Strain Model (1D) for Anisotropic Fibrils in Articular Cartilage

The mechanical response of articular cartilage (AC) under compression is anisotropic and depth-dependent. AC is osmotically active, and its intrinsic osmotic swelling pressure is balanced by its collagen fibril network. This mechanism requires the collagen fibers to be under a state of tensile pre-strain. A simple mathematical model is used to explain the depth-dependent strain calculations observed in articular cartilage under 1D axial compression (perpendicular to the articular surface). The collagen fibers are under pre-strain, influenced by proteoglycan concentration (fixed charged density, FCD) and collagen stiffness against swelling stress. The stiffness is introduced in our model as an anisotropic modulus that varies with fibril orientation through tissue depth. The collagen fibers are stiffer to stretching parallel to their length than perpendicular to it; when combined with depth-varying FCD, the model successfully predicts how tissue strains decrease with depth during compression. In summary, this model highlights that the mechanical properties of cartilage depend not only on proteoglycan concentration but also on the intrinsic properties of the pre-strained collagen network. These properties are essential for the proper functioning of articular cartilage.

Poroelastic behavior and water permeability of human skin at the nanoscale

Topical skin care products and hydrating compositions (moisturizers or injectable fillers) have been used for years to improve the appearance of, for example facial wrinkles, or to increase "plumpness". Most of the studies have addressed these changes based on the overall mechanical changes associated with an increase in hydration state. However, little is known about the water mobility contribution to these changes as well as the consequences to the specific skin layers. This is important as the biophysical properties and the biochemical composition of normal stratum corneum, epithelium, and dermis vary tremendously from one another. Our current studies and results reported here have focused on a novel approach (dynamic atomic force microscopy-based nanoindentation) to quantify biophysical characteristics of individual layers of ex vivo human skin. We have discovered that our new methods are highly sensitive to the mechanical properties of individual skin layers, as well as their hydration properties. Furthermore, our methods can assess the ability of these individual layers to respond to both compressive and shear deformations. In addition, since human skin is mechanically loaded over a wide range of deformation rates (frequencies), we studied the biophysical properties of skin over a wide frequency range. The poroelasticity model used helps to quantify the hydraulic permeability of the skin layers, providing an innovative method to evaluate and interpret the impact of hydrating compositions on water mobility of these different skin layers.

Computational assessment on the impact of collagen fiber orientation in cartilages on healthy and arthritic knee kinetics/kinematics

BACKGROUND

The inhomogeneous distribution of collagen fiber in cartilage can substantially influence the knee kinematics. This becomes vital for understanding the mechanical response of soft tissues, and cartilage deterioration including osteoarthritis (OA). Though the conventional computational models consider geometrical heterogeneity along with fiber reinforcements in the cartilage model as material heterogeneity, the influence of fiber orientation on knee kinetics and kinematics is not fully explored. This work examines how the collagen fiber orientation in the cartilage affects the healthy (intact knee) and arthritic knee response over multiple gait activities like running and walking.

METHODS

A 3D finite element knee joint model is used to compute the articular cartilage response during the gait cycle. A fiber-reinforced porous hyper elastic (FRPHE) material is used to model the soft tissue. A split-line pattern is used to implement the fiber orientation in femoral and tibial cartilage. Four distinct intact cartilage models and three OA models are simulated to assess the impact of the orientation of collagen fibers in a depth wise direction. The cartilage models with fibers oriented in parallel, perpendicular, and inclined to the articular surface are investigated for multiple knee kinematics and kinetics.

FINDINGS

The comparison of models with fiber orientation parallel to articulating surface for walking and running gait has the highest elastic stress and fluid pressure compared with inclined and perpendicular fiber-oriented models. Also, the maximum contact pressure is observed to be higher in the case of intact models during the walking cycle than for OA models. In contrast, the maximum contact pressure is higher during running in OA models than in intact models. Additionally, parallel-oriented models produce higher maximum stresses and fluid pressure for walking and running gait than proximal-distal-oriented models. Interestingly, during the walking cycle, the maximum contact pressure with intact models is approximately three times higher than on OA models. In contrast, the OA models exhibit higher contact pressure during the running cycle.

INTERPRETATION

Overall, the study indicates that collagen orientation is crucial for tissue responsiveness. This investigation provides insights into the development of tailored implants.

Load-induced fluid pressurisation in hydrogel systems before and after reinforcement by melt-electrowritten fibrous meshes

Fluid pressure develops transiently within mechanically-loaded, cell-embedding hydrogels, but its magnitude depends on the intrinsic material properties of the hydrogel and cannot be easily altered. The recently developed melt-electrowriting (MEW) technique enables three-dimensional printing of structured fibrous mesh with small fibre diameter (20 μm). The MEW mesh with 20 μm fibre diameter can synergistically increase the instantaneous mechanical stiffness of soft hydrogels. However, the reinforcing mechanism of the MEW meshes is not well understood, and may involve load-induced fluid pressurisation. Here, we examined the reinforcing effect of MEW meshes in three hydrogels: gelatin methacryloyl (GelMA), agarose and alginate, and the role of load-induced fluid pressurisation in the MEW reinforcement. We tested the hydrogels with and without MEW mesh (i.e., hydrogel alone, and MEW-hydrogel composite) using micro-indentation and unconfined compression, and analysed the mechanical data using biphasic Hertz and mixture models. We found that the MEW mesh altered the tension-to-compression modulus ratio differently for hydrogels that are cross-linked differently, which led to a variable change to their load-induced fluid pressurisation. MEW meshes only enhanced the fluid pressurisation for GelMA, but not for agarose or alginate. We speculate that only covalently cross-linked hydrogels (GelMA) can effectively tense the MEW meshes, thereby enhancing the fluid pressure developed during compressive loading. In conclusion, load-induced fluid pressurisation in selected hydrogels was enhanced by MEW fibrous mesh, and may be controlled by MEW mesh of different designs in the future, thereby making fluid pressure a tunable cell growth stimulus for tissue engineering involving mechanical stimulation.

Control of hydrostatic pressure and osmotic stress in 3D cell culture for mechanobiological studies

Hydrostatic pressure (HP) and osmotic stress (OS) play an important role in various biological processes, such as cell proliferation and differentiation. In contrast to canonical mechanical signals transmitted through the anchoring points of the cells with the extracellular matrix, the physical and molecular mechanisms that transduce HP and OS into cellular functions remain elusive. Three-dimensional cell cultures show great promise to replicate physiologically relevant signals in well-defined host bioreactors with the goal of shedding light on hidden aspects of the mechanobiology of HP and OS. This review starts by introducing prevalent mechanisms for the generation of HP and OS signals in biological tissues that are subject to pathophysiological mechanical loading. We then revisit various mechanisms in the mechanotransduction of HP and OS, and describe the current state of the art in bioreactors and biomaterials for the control of the corresponding physical signals.

A Numerical Scheme for Anisotropic Reactive Nonlinear Viscoelasticity

Reactive viscoelasticity is a theoretical framework based on the theory of reactive constrained mixtures that encompasses nonlinear viscoelastic responses. It models a viscoelastic solid as a mixture of strong and weak bonds that maintain the cohesiveness of the molecular constituents of the solid matter. Strong bonds impart the elastic response while weak bonds break and reform into a stress-free state in response to loading. The process of bonds breaking and reforming is modeled as a reaction where loaded bonds are the reactants and bonds reformed into a stress-free state are the products of a reaction. The reaction is triggered by the evolving state of loading. The state of stress in strong bonds is a function of the total strain in the material, whereas the state of stress in weak bonds is based on the state of strain relative to the time that these bonds were reformed. This study introduces two important practical contributions to the reactive nonlinear viscoelasticity framework: (1) normally, the evaluation of the stress tensor involves taking a summation over a continually increasing number of weak bond generations, which is poorly suited for a computational scheme. Therefore, this study presents an effective numerical scheme for evaluating the strain energy density, the Cauchy stress, and spatial elasticity tensors of reactive viscoelastic materials. (2) We provide the conditions for satisfying frame indifference for anisotropic nonlinear viscoelasticity, including for tension-bearing fiber models. Code verifications and model validations against experimental data provide evidence in support of this updated formulation.

Minimum design requirements for a poroelastic mimic of articular cartilage

The exceptional functional performance of articular cartilage (load-bearing and lubrication) is attributed to its poroelastic structure and resulting interstitial fluid pressure. Despite this, there remains no engineered cartilage repair material capable of achieving physiologically relevant poroelasticity. In this work we develop in silico models to guide the design approach for poroelastic mimics of articular cartilage. We implement the constitutive models in FEBio, a PDE solver for multiphasic mechanics problems in biological and soft materials. We investigate the influence of strain rate, boundary conditions at the contact interface, and fiber modulus on the reaction force and load sharing between the solid and fluid phases. The results agree with the existing literature that when fibers are incorporated the fraction of load supported by fluid pressure is greatly amplified and increases with the fiber modulus. This result demonstrates that a stiff fibrous phase is a primary design requirement for poroelastic mimics of articular cartilage. The poroelastic model is fit to experimental stress-relaxation data from bovine and porcine cartilage to determine if sufficient design constraints have been identified. In addition, we fit experimental data from FiHy™, an engineered material which is claimed to be poroelastic. The fiber-reinforced poroelastic model was able to capture the primary physics of these materials and demonstrates that FiHy™ is beginning to approach a cartilage-like poroelastic response. We also develop a fiber-reinforced poroelastic model with a bonded interface (rigid contact) to fit stress relaxation data from an osteochondral explant and FiHy™ + bone substitute. The model fit quality is similar for both the chondral and osteochondral configurations and clearly captures the first order physics. Based on this, we propose that physiological poroelastic mimics of articular cartilage should be developed under a fiber-reinforced poroelastic framework.

The intrinsic quality of proteoglycans, but not collagen fibres, degrades in osteoarthritic cartilage

The function of articular cartilage as a load-bearing connective tissue is derived primarily from a balanced interaction between the swelling proteoglycan (PG) matrix and tension-resistant collagen fibrous network. Such balance is compromised during joint disease such as osteoarthritis (OA) due to degradation to PGs and/or collagens. While the PG degradation is generally thought to be related to a loss of protein abundance, the collagenous degradation is more complex as it can be caused independently by a decrease of collagen content, disorganisation of fibrous structure and softening of individual collagen fibrils. A comprehensive understanding of the initial trajectories of degradation of PGs and collagen network can improve our chance of finding potential therapeutic solutions for OA. Here, we developed geometrically, structurally, and compositionally realistic and sample-specific Finite Element (FE) models under the framework of multiphasic mixture theory, from which the elastic moduli of collagen fibres and the PG load-bearing quality in healthy and diseased cartilages were estimated by numerical optimisation of the multi-step indentation stress relaxation force-time curves. We found the intrinsic quality of collagen fibres, measured by their elastic moduli, to stay constant for healthy and diseased cartilages. Combining with previous findings which show unaltered collagen content during early stages of OA, our results suggest the disorganisation of collagen fibrous network as the first form of collagenous degradation in osteoarthritic cartilage. We also found that PG degradation involves not only a loss of protein abundance, but also the quality of the remaining PGs in generating sufficient osmotic pressure for load bearing. This study sheds light on the mechanism of OA pathogenesis and highlights the restoration of collageneous organisation in cartilage as key medical intervention for OA. STATEMENT OF SIGNIFICANCE: Collagen network in articular cartilage consists of individual fibres that are organised into depth-dependent structure specialised for joint load-bearing and lubrication. During osteoarthritis, the collagen network undergoes mechanical degradation, but it is unclear if a loss of content, disorganisation of fibrous structure, or softening of individual fibres causes this degeneration. Using mechanical indentation, Finite Element modelling, and numerical optimisation methods, we determined that individual fibres did not soften in early disease stage. Together with previous findings showing unaltered collagen content, our results pinpoint the disorganisation of collagen structure as the main culprit for early collagenous degradation in osteoarthritic cartilage. Thus, early restoration in cartilage of collagen organisation, instead of individual fibre quality, may be key to slow osteoarthritis development.

Deformation behaviors and mechanical impairments of tissue cracks in immature and mature cartilages

Degeneration of articular cartilage is often triggered by a small tissue crack. As cartilage structure and composition change with age, the mechanics of cracked cartilage may depend on the tissue age, but this relationship is poorly understood. Here, we investigated cartilage mechanics and crack deformation in immature and mature cartilage exposed to a full-thickness tissue crack using indentation testing and histology, respectively. When a cut was introduced, tissue cracks opened wider in the mature cartilage compared to the immature cartilage. However, the opposite occurred upon mechanical indentation over the cracked region. Functionally, the immature-cracked cartilages stress-relaxed faster, experienced increased tissue strain, and had reduced instantaneous stiffness, compared to the mature-cracked cartilages. Taken together, mature cartilage appears to withstand surface cracks and maintains its mechanical properties better than immature cartilage and these superior properties can be explained by the structure of their collagen fibrous network.

Pericellular heparan sulfate proteoglycans: Role in regulating the biosynthetic response of nucleus pulposus cells to osmotic loading

Background

Daily physiologic loading causes fluctuations in hydration of the intervertebral disc (IVD); thus, the embedded cells experience cyclic alterations to their osmotic environment. These osmotic fluctuations have been described as a mechanism linking mechanics and biology, and have previously been shown to promote biosynthesis in chondrocytes. However, this phenomenon has yet to be fully interrogated in the IVD. Additionally, the specialized extracellular matrix surrounding the cells, the pericellular matrix (PCM), transduces the biophysical signals that cells ultimately experience. While it is known that the PCM is altered in disc degeneration, whether it disrupts normal osmotic mechanotransduction has yet to be determined. Thus, our objectives were to assess: (1) whether dynamic osmotic conditions stimulate biosynthesis in nucleus pulposus cells, and (2) whether pericellular heparan sulfate proteoglycans (HSPGs) modulate the biosynthetic response to osmotic loading.

Methods

Bovine nucleus pulposus cells isolated with retained PCM were encapsulated in 1.5% alginate beads and treated with or without heparinase III, an enzyme that degrades the pericellular HSPGs. Beads were subjected to 1 h of daily iso-osmotic, hyper-osmotic, or hypo-osmotic loading for 1, 2, or 4 weeks. At each timepoint the total amount of extracellular and pericellular sGAG/DNA were quantified. Additionally, whether osmotic loading triggered alterations to HSPG sulfation was assessed via immunohistochemistry for the heparan sulfate 6-O-sulfertransferase 1 (HS6ST1) enzyme.

Results

Osmotic loading significantly influenced sGAG/DNA accumulation with a hyper-osmotic change promoting the greatest sGAG/DNA accumulation in the pericellular region compared with iso-osmotic conditions. Heparanase-III treatment significantly reduced extracellular sGAG/DNA but pericellular sGAG was not affected. HS6ST1 expression was not affected by osmotic loading.

Conclusion

Results suggest that hyper-osmotic loading promotes matrix synthesis and that modifications to HSPGs directly influence the metabolic responses of cells to osmotic fluctuations. Collectively, results suggest degeneration-associated modifications to pericellular HSPGs may contribute to the altered mechanobiology observed in disease.

A Finite Element Algorithm for Large Deformation Biphasic Frictional Contact Between Porous-Permeable Hydrated Soft Tissues

The frictional response of porous and permeable hydrated biological tissues such as articular cartilage is significantly dependent on interstitial fluid pressurization. To model this response, it is common to represent such tissues as biphasic materials, consisting of a binary mixture of a porous solid matrix and an interstitial fluid. However, no computational algorithms currently exist in either commercial or open-source software that can model frictional contact between such materials. Therefore, this study formulates and implements a finite element algorithm for large deformation biphasic frictional contact in the open-source finite element software FEBio. This algorithm relies on a local form of a biphasic friction model that has been previously validated against experiments, and implements the model into our recently-developed surface-to-surface (STS) contact algorithm. Contact constraints, including those specific to pressurized porous media, are enforced with the penalty method regularized with an active-passive augmented Lagrangian scheme. Numerical difficulties specific to challenging finite deformation biphasic contact problems are overcome with novel smoothing schemes for fluid pressures and Lagrange multipliers. Implementation accuracy is verified against semi-analytical solutions for biphasic frictional contact, with extensive validation performed using canonical cartilage friction experiments from prior literature. Essential details of the formulation are provided in this paper, and the source code of this biphasic frictional contact algorithm is made available to the general public.

Effects of solvent osmolarity and viscosity on cartilage energy dissipation under high-frequency loading

Articular cartilage is a spatially heterogeneous, dissipative biological hydrogel with a high fluid volume fraction. Although energy dissipation is important in the context of delaying cartilage damage, the dynamic behavior of articular cartilage equilibrated in media of varied osmolarity and viscosity is not widely understood. This study investigated the mechanical behaviors of cartilage when equilibrated to media of varying osmolarity and viscosity. Dynamic moduli and phase shift were measured at both low (1 Hz) and high (75-300 Hz) frequency, with cartilage samples compressed to varied offset strain levels. Increasing solution osmolarity and viscosity both independently resulted in larger energy dissipation and decreased dynamic modulus of cartilage at both low and high frequency. Mechanical property alterations induced by varying osmolarity are likely due to the change in permeability and fluid volume fraction within the tissue. The effects of solution viscosity are likely due to frictional interactions at the solid-fluid interface, affecting energy dissipation. These findings highlight the significance of interstitial fluid on the energy dissipation capabilities of the tissue, which can influence the onset of cartilage damage.

Through-thickness regional variation in the mechanical characteristics of the lumbar facet capsular ligament

The human lumbar facet capsule, with the facet capsular ligament (FCL) that forms its primary constituent, is a common source of lower back pain. Prior studies on the FCL were limited to in-plane tissue behavior, but due to the presence of two distinct yet mechanically different regions, a novel out-of-plane study was conducted to further characterize the roles of the collagen and elastin regions. An experimental technique, called stretch-and-bend, was developed to study the tension-compression asymmetry of the FCL due to varying collagen fiber density throughout the thickness of the tissue. Each healthy excised cadaveric FCL sample was tested in four conditions depending on primary collagen fiber alignment and regional loading. Our results indicate that the FCL is stiffest when the collagen fibers (1) are aligned in the direction of loading, (2) are in tension, and (3) are stretched - 16% from its off-the-bone, undeformed state. An optimization routine was used to fit a four-parameter anisotropic, hyperplastic model to the experimental data. The average elastin modulus, E, and the average collagen fiber modulus, ξ, were 13.15 ± 3.59 kPa and 18.68 ± 13.71 MPa (95% CI), respectively.

Reproducing the Biomechanical Environment of the Chondrocyte for Cartilage Tissue Engineering

It is well known that the biomechanical and tribological performance of articular cartilage is inextricably linked to its extracellular matrix (ECM) structure and zonal heterogeneity. Furthermore, it is understood that the presence of native ECM components, such as collagen II and aggrecan, promote healthy homeostasis in the resident chondrocytes. What is less frequently discussed is how chondrocyte metabolism is related to the extracellular mechanical environment, at both the macro and microscale. The chondrocyte is in immediate contact with the pericellular matrix of the chondron, which acts as a mechanocoupler, transmitting external applied loads from the ECM to the chondrocyte. Therefore, components of the pericellular matrix also play essential roles in chondrocyte mechanotransduction and metabolism. Recreating the biomechanical environment through tuning material properties of a scaffold and/or the use of external cyclic loading can induce biosynthetic responses in chondrocytes. Decellularized scaffolds, which retain the native tissue macro- and microstructure also represent an effective means of recapitulating such an environment. The use of such techniques in tissue engineering applications can ensure the regeneration of skeletally mature articular cartilage with appropriate biomechanical and tribological properties to restore joint function. Despite the pivotal role in graft maturation and performance, biomechanical and tribological properties of such interventions is often underrepresented. This review outlines the role of biomechanics in relation to native cartilage performance and chondrocyte metabolism, and how application of this theory can enhance the future development and successful translation of biomechanically relevant tissue engineering interventions. Impact statement Physiological cartilage function is a key criterion in the success of a cartilage tissue engineering solution. The in situ performance is dependent on the initial scaffold design as well as extracellular matrix deposition by endogenous or exogenous cells. Both biological and biomechanical stimuli serve as key regulators of cartilage homeostasis and maturation of the resulting tissue-engineered graft. An improved understanding of the influence of biomechanics on cellular function and consideration of the final biomechanical and tribological performance will help in the successful development and translation of tissue-engineered grafts to restore natural joint function postcartilage trauma or osteoarthritic degeneration, delaying the requirement for prosthetic intervention.

Evolution of Meniscal Biomechanical Properties with Growth: An Experimental and Numerical Study

The menisci of the knee are complex fibro-cartilaginous tissues that play important roles in load bearing, shock absorption, joint lubrication, and stabilization. The objective of this study was to evaluate the interaction between the different meniscal tissue components (i.e., the solid matrix constituents and the fluid phase) and the mechanical response according to the developmental stage of the tissue. Menisci derived from partially and fully developed pigs were analyzed. We carried out biochemical analyses to quantify glycosaminoglycan (GAG) and DNA content according to the developmental stage. These values were related to tissue mechanical properties that were measured in vitro by performing compression and tension tests on meniscal specimens. Both compression and tension protocols consisted of multi-ramp stress-relaxation tests comprised of increasing strains followed by stress-relaxation to equilibrium. To better understand the mechanical response to different directions of mechanical stimulus and to relate it to the tissue structural composition and development, we performed numerical simulations that implemented different constitutive models (poro-elasticity, viscoelasticity, transversal isotropy, or combinations of the above) using the commercial software COMSOL Multiphysics. The numerical models also allowed us to determine several mechanical parameters that cannot be directly measured by experimental tests. The results of our investigation showed that the meniscus is a non-linear, anisotropic, non-homogeneous material: mechanical parameters increase with strain, depend on the direction of load, and vary among regions (anterior, central, and posterior). Preliminary numerical results showed the predominant role of the different tissue components depending on the mechanical stimulus. The outcomes of biochemical analyses related to mechanical properties confirmed the findings of the numerical models, suggesting a specific response of meniscal cells to the regional mechanical stimuli in the knee joint. During maturation, the increase in compressive moduli could be explained by cell differentiation from fibroblasts to metabolically active chondrocytes, as indicated by the found increase in GAG/DNA ratio. The changes of tensile mechanical response during development could be related to collagen II accumulation during growth. This study provides new information on the changes of tissue structural components during maturation and the relationship between tissue composition and mechanical response.

Collagen fibres determine the crack morphology in articular cartilage

Cracks in articular cartilage compromise tissue integrity and mechanical properties and lead to chondral lesions if untreated. An understanding of the mechanics of cracked cartilage may help in the prevention of cartilage deterioration and the development of tissue-engineered substitutes. The degeneration of cartilage in the presence of cracks may depend on the ultrastructure and composition of the tissue, which changes with aging, disease and habitual loading. It is unknown if the structural and compositional differences between immature and mature cartilage affect the mechanics of cartilage cracks, possibly predisposing one to a greater risk of degeneration than the other. We used a fibre-reinforced poro-viscoelastic swelling material model that accounts for large deformations and tension-compression non-linearity, and the finite element method to investigate the role of cartilage structure and composition on crack morphology and tissue mechanics. We demonstrate that the crack morphology predicted by our theoretical model agrees well with the histo-morphometric images of young and mature cracked cartilages under indentation loading. We also determined that the crack morphology was primarily dependent on collagen fibre orientation which differs as a function of cartilage depth and tissue maturity. The arcade-like collagen fibre orientation, first discussed by Benninghoff in his classical 1925 paper, appears to be beneficial for slowing the progression of tissue cracks by 'sealing' the crack and partially preserving fluid pressure during loading. Preservation of the natural load distribution between solid and fluid constituents of cartilage may be a key factor in slowing or preventing the propagation of tissue cracks and associated tissue matrix damage. STATEMENT OF SIGNIFICANCE: Cracks in articular cartilage can be detrimental to joint health if not treated, but it is not clear how they propagate and lead to tissue degradation. We used an advanced numerical model to determine the role of cartilage structure and composition on crack morphology under loading. Based on the structure and composition found in immature and mature cartilages, our model successfully predicts the crack morphology in these cartilages and determines that collagen fibre as the major determinant of crack morphology. The arcade-like Benninghoff collagen fibre orientation appears to be crucial in 'sealing' the tissue crack and preserves normal fluid-solid load distribution in cartilage. Inclusion of the arcade-like fibre orientation in tissue-engineered construct may help improve its integration within the host tissue.

Direct Osmotic Pressure Measurements in Articular Cartilage Demonstrate Nonideal and Concentration-Dependent Phenomena

The osmotic pressure in articular cartilage serves an important mechanical function in healthy tissue. Its magnitude is thought to play a role in advancing osteoarthritis. The aims of this study were to: (1) isolate and quantify the magnitude of cartilage swelling pressure in situ; and (2) identify the effect of salt concentration on material parameters. Confined compression stress-relaxation testing was performed on 18 immature bovine and six mature human cartilage samples in solutions of varying osmolarities. Direct measurements of osmotic pressure revealed nonideal and concentration-dependent osmotic behavior, with magnitudes approximately 1/3 those predicted by ideal Donnan law. A modified Donnan constitutive behavior was able to capture the aggregate behavior of all samples with a single adjustable parameter. Results of curve-fitting transient stress-relaxation data with triphasic theory in febio demonstrated concentration-dependent material properties. The aggregate modulus HA increased threefold as the external concentration decreased from hypertonic 2 M to hypotonic 0.001 M NaCl (bovine: HA=0.420±0.109 MPa to 1.266±0.438 MPa; human: HA=0.499±0.208 MPa to 1.597±0.455 MPa), within a triphasic theory inclusive of osmotic effects. This study provides a novel and simple analytical model for cartilage osmotic pressure which may be used in computational simulations, validated with direct in situ measurements. A key finding is the simultaneous existence of Donnan osmotic and Poisson-Boltzmann electrostatic interactions within cartilage.

Human articular cartilage is orthotropic where microstructure, micromechanics, and chemistry vary with depth and split-line orientation

OBJECTIVE

Quantitative, micrometer length scale assessment of human articular cartilage is essential to enable progress toward new functional tissue engineering approaches, including utilization of emerging 3D bioprinting technologies, and for improved computational modeling of the osteochondral unit. Thus the objective of this study was to characterize the structural organization, material properties, and chemical composition of human skeletally mature articular cartilage with respect to depth and defined morphological features: normal to the articulating surface, parallel to the split-line, and transverse to the split-line.

METHOD

Three samples from the lateral femoral condyles of 4 healthy adult donors (55-61 years old) were evaluated via histology, second harmonic generation, microindentation, and Raman spectroscopy. All metrics were evaluated as a function of depth and direction relative to the split-line.

RESULTS

All donors presented with intact and healthy tissue. Collagen fiber orientation varied significantly between testing directions and with increasing depth from the articular surface. Both compressive and tensile modulus increased significantly with depth and differed across the middle and deep zones and depended on orthogonal direction relative to the split-line. Similarly, matrix components varied with both depth and direction, where chondroitin sulfate steadily increased with depth while collagen prevalence was highest in the surface layer.

CONCLUSIONS

Microscale measurements of human articular cartilage demonstrate that properties are both depth-dependent and orthotropic and depend on the underlying tissue structure and composition. These findings improve upon existing knowledge establishing more accurate measurements, with greater degree of depth and spatial specificity, as inputs for tissue engineering and computational modeling.

Transient stiffening of cartilage during joint articulation: A microindentation study

As a mechanoactive tissue, articular cartilage undergoes compression and shear on a daily basis. With the advent of high resolution and sensitive mechanical testing methods, such as micro- and nanoindentation, it has become possible to assess changes in small-scale mechanical properties due to compression and shear of the tissue. However, investigations on the changes of these properties before and after joint articulation have been limited. To simulate articular loading of cartilage in the context of human gait, a previously developed bioreactor system was used. Immediately after bioreactor testing, the stiffness was measured using microindentation. Specifically, we investigated whether the mechanical response of the tissue was transient or permanent, dependent on counterface material, and an effect limited to the superficial zone of cartilage. We found that cartilage surface stiffness increases immediately after articular loading and returns to baseline values within 3 hr. Cartilage-on-cartilage stiffening was found to be higher compared to both alumina- and cobalt chromium-on-cartilage stiffening, which were not significantly different from each other. This stiffening response was found to be unique to the superficial zone, as articular loading on cartilage with the superficial zone removed showed no changes in stiffness. The findings of this study suggest that the cartilage superficial zone may adapt its stiffness as a response to articular loading. As the superficial zone is often compromised during the course of osteoarthritic disease, this finding is of clinical relevance, suggesting that the load-bearing function deteriorates over time.

Apparatus and Method for Rapid Detection of Acoustic Anisotropy in Cartilage

Purpose

Articular cartilage is known to be mechanically anisotropic. In this paper, the acoustic anisotropy of bovine articular cartilage and the effects of freeze-thaw cycling on acoustic anisotropy were investigated.

Methods

We developed apparatus and methods that use a magnetic L-shaped sample holder, which allowed minimal handling of a tissue, reduced the number of measurements compared to previous studies, and produced highly reproducible results.

Results

SOS was greater in the direction perpendicular to the articular surface compared to the direction parallel to the articular surface (N=17, P = 0.00001). Average SOS was 1,758 ± 107 m/s perpendicular to the surface, and 1,617 ± 55 m/s parallel to it. The average percentage difference in SOS between the perpendicular and parallel directions was 8.2% (95% CI: 5.4% to 11%). Freeze-thaw cycling did not have a significant effect on SOS (P>0.4).

Conclusion

Acoustic measurement of tissue properties is particularly attractive for work in our laboratory since it has the potential for nondestructive characterization of the properties of developing engineered cartilage. Our approach allowed us to observe acoustic anisotropy of articular cartilage rapidly and reproducibly. This property was not significantly affected by freeze-thawing of the tissue samples, making cryopreservation practical for these assays.

Evidence that reduction in volume protects in situ articular chondrocytes from mechanical impact

Chondrocytes, the resident cells in articular cartilage, carry the burden of producing and maintaining the extracellular matrix (ECM). However, as these cells have a low proliferative capacity and are not readily replaced, chondrocyte death due to extreme forces may contribute to the pathogenesis of osteoarthritis (OA) after injury or may inhibit healing after osteochondral transplantation, a restorative procedure for damaged cartilage that requires a series of mechanical impacts to insert the graft. Consequently, there is a need to understand what factors influence the vulnerability of in situ chondrocytes to mechanical trauma. To this end, the objective of this study was to investigate how altering cell volume by different means (hydrostatic pressure, uniaxial load, and osmotic challenge with and without inhibition of regulatory volume decrease) affects the vulnerability of in situ chondrocytes to extreme mechanical forces. Using a custom experimental platform enabling testing of viable and intact murine cartilage-on-bone explants, we established a strong correlation between chondrocyte volume and vulnerability to impact injury wherein reduced volume was protective. Moreover, we found that the volume-perturbing interventions did not affect cartilage ECM mechanical properties, suggesting that their effects on chondrocyte vulnerability occurred at the cellular level. The findings of this study offer new avenues for novel strategies aimed at preventing chondrocyte loss during osteochondral grafting or to halting the progression of cell death after a joint destabilizing injury.

Experimental Study on Creep Characteristics of Microdefect Articular Cartilages in the Damaged Early Stage

Traumatic joint injury is known to cause cartilage deterioration and osteoarthritis. In order to study the mechanical mechanism of damage evolution on articular cartilage, taking the fresh porcine articular cartilage as the experimental samples, the creep experiments of the intact cartilages and the cartilages with different depth defect were carried out by using the noncontact digital image correlation technology. And then, the creep constitutive equations of cartilages were established. The results showed that the creep curves of different layers changed exponentially and were not coincident for the cartilage sample. The defect affected the strain values of the creep curves. The creep behavior of cartilage was dependent on defect depth. The deeper the defect was, the larger the strain value was. The built three-parameter viscoelastic constitutive equation had a good correlation with the experimental results and could predict the creep performance of the articular cartilage. The creep values of the microdefective cartilage in the damaged early stage were different from the diseased articular cartilage. These findings pointed out that defect could accelerate the damage of cartilage. It was helpful to study the mechanical mechanism of damage evolution.

Raman Spectroscopy: Guiding Light for the Extracellular Matrix

The extracellular matrix (ECM) consists of a complex mesh of proteins, glycoproteins, and glycosaminoglycans, and is essential for maintaining the integrity and function of biological tissues. Imaging and biomolecular characterization of the ECM is critical for understanding disease onset and for the development of novel, disease-modifying therapeutics. Recently, there has been a growing interest in the use of Raman spectroscopy to characterize the ECM. Raman spectroscopy is a label-free vibrational technique that offers unique insights into the structure and composition of tissues and cells at the molecular level. This technique can be applied across a broad range of ECM imaging applications, which encompass in vitro, ex vivo, and in vivo analysis. State-of-the-art confocal Raman microscopy imaging now enables label-free assessments of the ECM structure and composition in tissue sections with a remarkably high degree of biomolecular specificity. Further, novel fiber-optic instrumentation has opened up for clinical in vivo ECM diagnostic measurements across a range of tissue systems. A palette of advanced computational methods based on multivariate statistics, spectral unmixing, and machine learning can be applied to Raman data, allowing for the extraction of specific biochemical information of the ECM. Here, we review Raman spectroscopy techniques for ECM characterizations over a variety of exciting applications and tissue systems, including native tissue assessments (bone, cartilage, cardiovascular), regenerative medicine quality assessments, and diagnostics of disease states. We further discuss the challenges in the widespread adoption of Raman spectroscopy in biomedicine. The results of the latest discovery-driven Raman studies are summarized, illustrating the current and potential future applications of Raman spectroscopy in biomedicine.

A Systematic Review and Guide to Mechanical Testing for Articular Cartilage Tissue Engineering

Articular cartilage is integral to the mechanical function of many joints in the body. When injured, cartilage lacks the capacity to self-heal, and thus, therapies and replacements have been developed in recent decades to treat damaged cartilage. Given that the primary function of articular cartilage is mechanical in nature, rigorous physical evaluation of cartilage tissues undergoing treatment and cartilage constructs intended for replacement is an absolute necessity. With the large number of groups developing cartilage tissue engineering strategies, however, a variety of mechanical testing protocols have been reported in the literature. This lack of consensus in testing methods makes comparison between studies difficult at times, and can lead to misinterpretation of data relative to native tissue. Therefore, the purpose of this study was to systematically review mechanical testing of articular cartilage and cartilage repair constructs over the past 10 years (January 2009-December 2018), to highlight the most common testing configurations, and to identify key testing parameters. For the most common tests, key parameters identified in this systematic review were validated by characterizing both cartilage tissue and hydrogels commonly used in cartilage tissue engineering. Our findings show that compression testing was the most common test performed (80.2%; 158/197), followed by evaluation of frictional properties (18.8%; 37/197). Upon further review of those studies performing compression testing, the various modes (ramp, stress relaxation, creep, dynamic) and testing configurations (unconfined, confined, in situ) are described and systematically reviewed for parameters, including strain rate, equilibrium time, and maximum strain. This systematic analysis revealed considerable variability in testing methods. Our validation testing studies showed that such variations in testing criteria could have large implications on reported outcome parameters (e.g., modulus) and the interpretation of findings from these studies. This analysis is carried out for all common testing methods, followed by a discussion of less common trends and directions in the mechanical evaluation of cartilage tissues and constructs. Overall, this work may serve as a guide for cartilage tissue engineers seeking to rigorously evaluate the physical properties of their novel treatment strategies. Impact Statement Articular cartilage tissue engineering has made significant strides with regard to treatments and replacements for injured tissue. The evaluation of these approaches typically involves mechanical testing, yet the plethora of testing techniques makes comparisons between studies difficult, and often leads to misinterpretation of data compared with native tissue. This study serves as a guide for the mechanical testing of cartilage tissues and constructs, highlighting recent trends in test conditions and validating these common procedures. Cartilage tissue engineers, especially those unfamiliar with mechanical testing protocols, will benefit from this study in their quest to physically evaluate novel treatment and regeneration approaches.

Simultaneous tractography and elastography imaging of the zone-specific structural and mechanical responses in articular cartilage under compressive loading

We quantified the precise zonal cartilage structural and mechanical responses to unconfined compressive loading by using simultaneous PSOCT based optical tractography and elastography imaging. Twelve bovine knee articular cartilage samples from six animals were imaged under bulk compression from 4% to 20%. The results revealed strong evidence that the conventional radial zone could be divided into two sub-zones with distinct mechanical properties. The "upper" part of the radial zone played a critical role in "absorbing" the mechanical compression. The study also showed that the zonal fiber organization greatly affected the cartilage structural and mechanical responses. A strong correlation was observed between the optical birefringence and logarithm of the Young's modulus. These new results provide useful information for improving mechanical modeling of articular cartilage and developing better cartilage-mimetic biomaterials.

Influences of sodium and glycosaminoglycans on skin oedema and the potential for ulceration: a finite-element approach

Venous ulcers are chronic transcutaneous wounds common in the lower legs. They are resistant to healing and have a 78% chance of recurrence within 2 years. It is commonly accepted that venous ulcers are caused by the insufficiency of the calf muscle pump, leading to blood pooling in the lower legs, resulting in inflammation, skin oedema, tissue necrosis and eventually skin ulceration. However, the detailed physiological events by which inflammation contributes to wound formation are poorly understood. We therefore sought to develop a model that simulated the inflammation, using it to determine the internal stresses and pressure on the skin that contribute to venous ulcer formation. A three-layer finite-element skin model (epidermis, dermis and hypodermis) was developed to explore the roles in wound formation of two inflammation identifiers: glycosaminoglycans (GAG) and sodium. A series of parametric studies showed that increased GAG and sodium content led to oedema and increased tissue stresses of 1.5 MPa, which was within the reported range of skin tissue ultimate tensile stress (0.1-40 MPa). These results suggested that both the oedema and increased fluid pressure could reach a threshold for tissue damage and eventual ulcer formation. The models presented here provide insights to the pathological events associated with venous insufficiency, including inflammation, oedema and skin ulceration.

Biological connective tissues exhibit viscoelastic and poroelastic behavior at different frequency regimes: Application to tendon and skin biophysics

In this study, a poroviscoelastic finite element model (FEM) was developed and used in conjunction with an AFM-based wide-bandwidth nanorheology system to predict the frequency-dependent mechanical behavior of tendon and dermis subjected to compression via nanoindentation. The aim was to distinguish between loading rates that are dominated by either poroelasticity, viscoelasticity, or the superposition of these processes. Using spherical probe tips having different radii, the force and tip displacement were measured and the magnitude, E^∗, and phase angle, ϕ, of the dynamic complex modulus were evaluated for mouse supraspinatus tendon and mouse dermis. The peak frequencies of the phase angle were associated with the characteristic time constants of poroelastic and viscoelastic material behavior. The developed FE model could predict the separate poroelastic and viscoelastic responses of these soft tissues over a 4 decade frequency range, showing good agreement with experimental results. We observed that poroelasticity was the dominant energy dissipation mechanism for mouse dermis and supraspinatus tendon at higher indentation frequencies (10² to 10⁴ Hz) whereas viscoelasticity was typically dominant at lower frequencies (<10² Hz). These findings show the underlying mechanical behavior of biological connective tissues and give insight into the role played by these different energy dissipation mechanisms in governing the function of these tissues at nanoscale.

STATEMENT OF SIGNIFICANCE

Soft biological tissues exhibit complex, load- and time-dependent mechanical behavior. Evaluating their mechanical behavior requires sophisticated experimental tools and numerical models that can capture the fundamental mechanisms governing tissue function. Using an Atomic-force-microscopy-based rheology system and finite element models, the roles of the two most dominant time-dependent mechanisms (poroelasticity and viscoelasticity) that govern the dynamic loading behavior of mouse skin and tendon have been investigated. FE models were able to predict and quantify the contribution of each mechanism to the overall dynamic response and confirming the presence of these two distinct mechanisms in the mechanical response. Overall, these results provide novel insight into the viscoelastic and poroelastic properties of mouse skin and tendon and promote better understanding of the underlying origins of each mechanism.

Indentation mapping revealed poroelastic, but not viscoelastic, properties spanning native zonal articular cartilage

Osteoarthrosis is a debilitating disease affecting millions, yet engineering materials for cartilage regeneration has proven difficult because of the complex microstructure of this tissue. Articular cartilage, like many biological tissues, produces a time-dependent response to mechanical load that is critical to cell's physiological function in part due to solid and fluid phase interactions and property variations across multiple length scales. Recreating the time-dependent strain and fluid flow may be critical for successfully engineering replacement tissues but thus far has largely been neglected. Here, microindentation is used to accomplish three objectives: (1) quantify a material's time-dependent mechanical response, (2) map material properties at a cellular relevant length scale throughout zonal articular cartilage and (3) elucidate the underlying viscoelastic, poroelastic, and nonlinear poroelastic causes of deformation in articular cartilage. Untreated and trypsin-treated cartilage was sectioned perpendicular to the articular surface and indentation was used to evaluate properties throughout zonal cartilage on the cut surface. The experimental results demonstrated that within all cartilage zones, the mechanical response was well represented by a model assuming nonlinear biphasic behavior and did not follow conventional viscoelastic or linear poroelastic models. Additionally, 10% (w/w) agarose was tested and, as anticipated, behaved as a linear poroelastic material. The approach outlined here provides a method, applicable to many tissues and biomaterials, which reveals and quantifies the underlying causes of time-dependent deformation, elucidates key aspects of material structure and function, and that can be used to provide important inputs for computational models and targets for tissue engineering.

STATEMENT OF SIGNIFICANCE

Elucidating the time-dependent mechanical behavior of cartilage, and other biological materials, is critical to adequately recapitulate native mechanosensory cues for cells. We used microindentation to map the time-dependent properties of untreated and trypsin treated cartilage throughout each cartilage zone. Unlike conventional approaches that combine viscoelastic and poroelastic behaviors into a single framework, we deconvoluted the mechanical response into separate contributions to time-dependent behavior. Poroelastic effects in all cartilage zones dominated the time-dependent behavior of articular cartilage, and a model that incorporates tension-compression nonlinearity best represented cartilage mechanical behavior. These results can be used to assess the success of regeneration and repair approaches, as design targets for tissue engineering, and for development of accurate computational models.

Reactive Constrained Mixtures for Modeling the Solid Matrix of Biological Tissues

This article illustrates our approach for modeling the solid matrix of biological tissues using reactive constrained mixtures. Several examples are presented to highlight the potential benefits of this approach, showing that seemingly disparate fields of mechanics and chemical kinetics are actually closely interrelated and may be elegantly expressed in a unified framework. Thus, constrained mixture models recover classical theories for fibrous materials with bundles oriented in different directions or having different reference configurations, that produce characteristic fiber recruitment patterns under loading. Reactions that exchange mass among various constituents of a mixture may be used to describe tissue growth and remodeling, which may also alter the material's anisotropy. Similarly, reactions that describe the breaking and reforming of bonds may be used to model free energy dissipation in a viscoelastic material. Therefore, this framework is particularly well suited for modeling biological tissues.

How a decreased fibrillar interconnectivity influences stiffness and swelling properties during early cartilage degeneration

OBJECTIVE

The functional coupling between the fibrillar network and the high-swelling proteoglycans largely determines the mechanical properties of the articular cartilage matrix. The objective of this new study was to show specifically how changes in fibrillar interconnectivity arising from early cartilage degeneration influence transverse stiffness and swelling properties at the tissue level.

DESIGN

Radial zone transverse layers of cartilage matrix were obtained from intact and mildly degenerate bovine patellae. Each layer was then subdivided to assess tensile stiffness, free-swelling response, glycosaminoglycan (GAG) content, and micro- and ultra-structural features.

RESULTS

The tensile modulus was significantly lower and the degree of swelling significantly higher for the degenerate matrix compared to the intact. Scanning electron microscopy revealed a homogeneous response to transverse strain in the intact cartilage, whereas large non-fibrillar spaces between fibril aggregates were visible in the degenerate matrix. Although there were no significant differences in GAG content it did correlate significantly with stiffness and swelling in the intact samples but not in the degenerate.

CONCLUSIONS

The lower degree of fibril network interconnectivity in the degenerate matrix led to both a decreased transverse stiffness and reduced resistance to osmotic swelling. This network 'de-structuring' also resulted in a reduced functional interaction between the fibrillar network and the proteoglycans. The study provides new insights into the role of the fibrillar network and how changes in the network arising from the degenerative cascade will influence tissue level behaviour.

Finite Element Modeling of Avascular Tumor Growth Using a Stress-Driven Model

Tumor growth being a multistage process has been investigated from different aspects. In the present study, an attempt is made to represent a constitutive-structure-based model of avascular tumor growth in which the effects of tensile stresses caused by collagen fibers are considered. Collagen fibers as a source of anisotropy in the structure of tissue are taken into account using a continuous fiber distribution formulation. To this end, a finite element modeling is implemented in which a neo-Hookean hyperelastic material is assigned to the tumor and its surrounding host. The tumor is supplied with a growth term. The growth term includes the effect of parameters such as nutrient concentration on the tumor growth and the tumor's solid phase content in the formulation. Results of the study revealed that decrease of solid phase is indicative of decrease in growth rate and the final steady-state value of tumor's radius. Moreover, fiber distribution affects the final shape of the tumor, and it could be used to control the shape and geometry of the tumor in complex morphologies. Finally, the findings demonstrated that the exerted stresses on the tumor increase as time passes. Compression of tumor cells leads to the reduction of tumor growth rate until it gradually reaches an equilibrium radius. This finding is in accordance with experimental data. Hence, this formulation can be deployed to evaluate both the residual stresses induced by growth and the mechanical interactions with the host tissue.

Reinforcement of articular cartilage with a tissue-interpenetrating polymer network reduces friction and modulates interstitial fluid load support

OBJECTIVE

Osteoarthritis (OA) is associated with increased articular cartilage hydraulic permeability and decreased maintenance of high interstitial fluid load support (IFLS) during articulation, resulting in increased friction on the cartilage solid matrix. This study assesses frictional response following in situ synthesis of an interpenetrating polymer network (IPN) designed to mimic glycosaminoglycans (GAGs) depleted during OA.

METHODS

Cylindrical osteochondral explants containing various interpenetrating polymer concentrations were subjected to a torsional friction test under unconfined creep compression. Time-varying coefficient of friction, compressive engineering strain, and normalized strain values (ε/εeq) were calculated and analyzed.

RESULTS

The polymer network reduced friction coefficient over the duration of the friction test, with statistically significantly reduced friction coefficients (95% confidence interval 14-34% reduced) at equilibrium compressive strain upon completion of the test (P = 0.015). A positive trend was observed relating polymer network concentration with magnitude of friction reduction compared to non-treated tissue.

CONCLUSION

The cartilage-interpenetrating polymer treatment improves lubrication by augmenting the biphasic tissue's interstitial fluid phase, and additionally improves the friction dissipation of the tissue's solid matrix. This technique demonstrates potential as a therapy to augment tribological function of articular cartilage.

Ultrasound palpation for fast in-situ quantification of articular cartilage stiffness, thickness and relaxation capacity

Most current cartilage testing devices require the preparation of excised samples and therefore do not allow intra-operative application for diagnostic purposes. The gold standard during open or arthroscopic surgery is still the subjective perception of manual palpation. This work presents a new diagnostic method of ultrasound palpation (USP) to acquire applied stress and strain data during manual palpation of articular cartilage. With the proposed method, we obtain cartilage thickness and stiffness. Moreover, repeated palpations allow the quantification of relaxation effects. USP measurements on elastomer phantoms demonstrated very good repeatability for both, stage-guided (97.2%) and handheld (96.0%) applications. The USP measurements were compared with conventional indentation experiments and revealed very good agreement on elastomer phantoms ([Formula: see text]) and good agreement on porcine cartilage samples ([Formula: see text]). Artificially degenerated cartilage samples showed reduced stiffness, weak capacity to relax after palpation and an increase of stiffness of approximately 50% with each single palpation. Intact cartilage was measured by USP directly at the patella (in situ) and after excision and removal of the subchondral bone (ex situ), leading to stiffness values of [Formula: see text] and [Formula: see text] ([Formula: see text]), respectively. The results demonstrate the potential of the USP system for cartilage testing, its sensitivity to degenerative changes and as a method for quantifying relaxation processes by means of repeated palpations. Furthermore, the differences in the results of in-situ and ex-situ measurements are of general interest, since such comparison has not been reported previously. We point out the limited comparability of ex-situ cartilage with its in-situ biomechanical behavior.

Effect of Hydration on Healthy Intervertebral Disk Mechanical Stiffness

The intervertebral disk has an excellent swelling capacity to absorb water, which is thought to be largely due to the high proteoglycan composition. Injury, aging, degeneration, and diurnal loading are all noted by a significant decrease in water content and tissue hydration. The objective of this study was to evaluate the effect of hydration, through osmotic loading, on tissue swelling and compressive stiffness of healthy intervertebral disks. The wet weight of nucleus pulposus (NP) and annulus fibrosus (AF) explants following swelling was 50% or greater, demonstrating significant ability to absorb water under all osmotic loading conditions (0.015 M-3.0 M phosphate buffered saline (PBS)). Estimated NP residual strains, calculated from the swelling ratio, were approximately 1.5 × greater than AF residual strains. Compressive stiffness increased with hyperosmotic loading, which is thought to be due to material compaction from osmotic-loading and the nonlinear mechanical behavior. Importantly, this study demonstrated that residual strains and material properties are greatly dependent on osmotic loading. The findings of this study support the notion that swelling properties from osmotic loading will be important for accurately describing the effect of degeneration and injury on disk mechanics. Furthermore, the tissue swelling will be an important consideration for developing biological repair strategies aimed at restoring mechanical behavior toward a healthy disk.

Microstructural heterogeneity directs micromechanics and mechanobiology in native and engineered fibrocartilage

Treatment strategies to address pathologies of fibrocartilaginous tissue are in part limited by an incomplete understanding of structure-function relationships in these load-bearing tissues. There is therefore a pressing need to develop micro-engineered tissue platforms that can recreate the highly inhomogeneous tissue microstructures that are known to influence mechanotransductive processes in normal and diseased tissue. Here, we report the quantification of proteoglycan-rich microdomains in developing, ageing and diseased fibrocartilaginous tissues, and the impact of these microdomains on endogenous cell responses to physiologic deformation within a native-tissue context. We also developed a method to generate heterogeneous tissue-engineered constructs (hetTECs) with non-fibrous proteoglycan-rich microdomains engineered into the fibrous structure, and show that these hetTECs match the microstructural, micromechanical and mechanobiological benchmarks of native tissue. Our tissue-engineered platform should facilitate the study of the mechanobiology of developing, homeostatic, degenerating and regenerating fibrous tissues.

Dexamethasone Release from Within Engineered Cartilage as a Chondroprotective Strategy Against Interleukin-1α

While significant progress has been made toward engineering functional cartilage constructs with mechanical properties suitable for in vivo loading, the impact on these grafts of inflammatory cytokines, chemical factors that are elevated with trauma or osteoarthritis, is poorly understood. Previous work has shown dexamethasone to be a critical compound for cultivating cartilage with functional properties, while also providing chondroprotection from proinflammatory cytokines. This study tested the hypothesis that the incorporation of poly(lactic-co-glycolic acid) (PLGA) (75:25) microspheres that release dexamethasone from within chondrocyte-seeded agarose hydrogel constructs would promote development of constructs with functional properties and protect constructs from the deleterious effects of interleukin-1α (IL-1α). After 28 days of growth culture, experimental groups were treated with IL-1α (10 ng/mL) for 7 days. Reaching native equilibrium moduli and proteoglycan levels, dexamethasone-loaded microsphere constructs exhibited tissue properties similar to microsphere-free control constructs cultured in dexamethasone-supplemented culture media and were insensitive to IL-1α exposure. These findings are in stark contrast to constructs containing dexamethasone-free microspheres or no microspheres, cultured without dexamethasone, where IL-1α exposure led to significant tissue degradation. These results support the use of dexamethasone delivery from within engineered cartilage, through biodegradable microspheres, as a strategy to produce mechanically functional tissues that can also combat the deleterious effects of local proinflammatory cytokine exposure.

Continuum theory of fibrous tissue damage mechanics using bond kinetics: application to cartilage tissue engineering

This study presents a damage mechanics framework that employs observable state variables to describe damage in isotropic or anisotropic fibrous tissues. In this mixture theory framework, damage is tracked by the mass fraction of bonds that have broken. Anisotropic damage is subsumed in the assumption that multiple bond species may coexist in a material, each having its own damage behaviour. This approach recovers the classical damage mechanics formulation for isotropic materials, but does not appeal to a tensorial damage measure for anisotropic materials. In contrast with the classical approach, the use of observable state variables for damage allows direct comparison of model predictions to experimental damage measures, such as biochemical assays or Raman spectroscopy. Investigations of damage in discrete fibre distributions demonstrate that the resilience to damage increases with the number of fibre bundles; idealizing fibrous tissues using continuous fibre distribution models precludes the modelling of damage. This damage framework was used to test and validate the hypothesis that growth of cartilage constructs can lead to damage of the synthesized collagen matrix due to excessive swelling caused by synthesized glycosaminoglycans. Therefore, alternative strategies must be implemented in tissue engineering studies to prevent collagen damage during the growth process.

Human cartilage endplate permeability varies with degeneration and intervertebral disc site

Despite the critical functions the human cartilage endplate (CEP) plays in the intervertebral disc, little is known about its structural and mechanical properties and their changes with degeneration. Quantifying these changes with degeneration is important for understanding how the CEP contributes to the function and pathology of the disc. Therefore the objectives of this study were to quantify the effect of disc degeneration on human CEP mechanical properties, determine the influence of superior and inferior disc site on mechanics and composition, and simulate the role of collagen fibers in CEP and disc mechanics using a validated finite element model. Confined compression data and biochemical composition data were used in a biphasic-swelling model to calculate compressive extrafibrillar elastic and permeability properties. Tensile properties were obtained by applying published tensile test data to an ellipsoidal fiber distribution. Results showed that with degeneration CEP permeability decreased 50-60% suggesting that transport is inhibited in the degenerate disc. CEP fibers are organized parallel to the vertebrae and nucleus pulposus and may contribute to large shear strains (0.1-0.2) and delamination failure of the CEP commonly seen in herniated disc tissue. Fiber-reinforcement also reduces CEP axial strains thereby enhancing fluid flux by a factor of 1.8. Collectively, these results suggest that the structure and mechanics of the CEP may play critical roles in the solute transport and disc mechanics.

On the stiffness of demineralized dentin matrices

Resin bonding to dentin requires the use of self-etching primers or acid etching to decalcify the surface and expose a layer of collagen fibrils of the dentin matrix. Acid-etching reduces the stiffness of demineralized dentin from approximately 19 GPa-1 MPa, requiring that it floats in water to prevent it from collapsing during bonding procedures. Several publications show that crosslinking agents like gluteraladehyde, carbodiimide or grape seed extract can stiffen collagen and improve resin-dentin bond strength.

OBJECTIVE

The objective was to assess a new approach for evaluating the changes in stiffness of decalcified dentin by polar solvents and a collagen cross-linker.

METHODS

Fully demineralized dentin beams and sections of etched coronal dentin were subjected to indentation loading using a cylindrical flat indenter in water, and after treatment with ethanol or ethyl-3-(3-dimethylaminopropyl) carbodiimide (EDC). The stiffness was measured as a function of strain and as a function of loading rate from 1 to 50 μm/s.

RESULTS

At a strain of 0.25% the elastic modulus of the fully demineralized dentin was approximately 0.20 MPa. It increased to over 0.90 MPa at strains of 1%. Exposure to ethanol caused an increase in elastic modulus of up to four times. Increasing the loading rate from 1 to 50 μm/s caused an increase in the apparent modulus of up to three times in both water and ethanol. EDC treatment caused increases in the stiffness in fully demineralized samples and in acid-etched demineralized dentin surfaces in situ.

SIGNIFICANCE

Changes in the mechanical behavior of demineralized collagen matrices can be measured effectively under hydration via indentation with cylindrical flat indenters. This approach can be used for quantifying the effects of bonding treatments on the properties of decalcified dentin after acid etching, as well as to follow the loss of stiffness over time due to enzymatic degradation.

Determining Tension-Compression Nonlinear Mechanical Properties of Articular Cartilage from Indentation Testing

The indentation test is widely used to determine the in situ biomechanical properties of articular cartilage. The mechanical parameters estimated from the test depend on the constitutive model adopted to analyze the data. Similar to most connective tissues, the solid matrix of cartilage displays different mechanical properties under tension and compression, termed tension-compression nonlinearity (TCN). In this study, cartilage was modeled as a porous elastic material with either a conewise linear elastic matrix with cubic symmetry or a solid matrix reinforced by a continuous fiber distribution. Both models are commonly used to describe the TCN of cartilage. The roles of each mechanical property in determining the indentation response of cartilage were identified by finite element simulation. Under constant loading, the equilibrium deformation of cartilage is mainly dependent on the compressive modulus, while the initial transient creep behavior is largely regulated by the tensile stiffness. More importantly, altering the permeability does not change the shape of the indentation creep curves, but introduces a parallel shift along the horizontal direction on a logarithmic time scale. Based on these findings, a highly efficient curve-fitting algorithm was designed, which can uniquely determine the three major mechanical properties of cartilage (compressive modulus, tensile modulus, and permeability) from a single indentation test. The new technique was tested on adult bovine knee cartilage and compared with results from the classic biphasic linear elastic curve-fitting program.

Measuring microscale strain fields in articular cartilage during rapid impact reveals thresholds for chondrocyte death and a protective role for the superficial layer

Articular cartilage is a heterogeneous soft tissue that dissipates and distributes loads in mammalian joints. Though robust, cartilage is susceptible to damage from loading at high rates or magnitudes. Such injurious loads have been implicated in degenerative changes, including chronic osteoarthritis (OA), which remains a leading cause of disability in developed nations. Despite decades of research, mechanisms of OA initiation after trauma remain poorly understood. Indeed, although bulk cartilage mechanics are measurable during impact, current techniques cannot access microscale mechanics at those rapid time scales. We aimed to address this knowledge gap by imaging the microscale mechanics and corresponding acute biological changes of cartilage in response to rapid loading. In this study, we utilized fast-camera and confocal microscopy to achieve roughly 85 µm spatial resolution of both the cartilage deformation during a rapid (~3 ms), localized impact and the chondrocyte death following impact. Our results showed that, at these high rates, strain and chondrocyte death were highly correlated (p<0.001) with a threshold of 8% microscale strain norm before any cell death occurred. Additionally, chondrocyte death had developed by two hours after impact, suggesting a time frame for clinical therapeutics. Moreover, when the superficial layer was removed, strain - and subsequently chondrocyte death - penetrated deeper into the samples (p<0.001), suggesting a protective role for the superficial layer of articular cartilage. Combined, these results provide insight regarding the detailed biomechanics that drive early chondrocyte damage after trauma and emphasize the importance of understanding cartilage and its mechanics on the microscale.

Wear and damage of articular cartilage with friction against orthopedic implant materials

The objective of this study was to measure the wear response of immature bovine articular cartilage tested against glass or alloys used in hemiarthroplasties. Two cobalt chromium alloys and a stainless steel alloy were selected for these investigations. The surface roughness of one of the cobalt chromium alloys was also varied within the range considered acceptable by regulatory agencies. Cartilage disks were tested in a configuration that promoted loss of interstitial fluid pressurization to accelerate conditions believed to occur in hemiarthroplasties. Results showed that considerably more damage occurred in cartilage samples tested against stainless steel (10 nm roughness) and low carbon cobalt chromium alloy (27 nm roughness) compared to glass (10 nm) and smoother low or high carbon cobalt chromium (10 nm). The two materials producing the greatest damage also exhibited higher equilibrium friction coefficients. Cartilage damage occurred primarily in the form of delamination at the interface between the superficial tangential zone and the transitional middle zone, with much less evidence of abrasive wear at the articular surface. These results suggest that cartilage damage from frictional loading occurs as a result of subsurface fatigue failure leading to the delamination. Surface chemistry and surface roughness of implant materials can have a significant influence on tissue damage, even when using materials and roughness values that satisfy regulatory requirements.

Micro-poromechanics model of fluid-saturated chemically active fibrous media

We have developed a micromechanics based model for chemically active saturated fibrous media that incorporates fiber network microstructure, chemical potential driven fluid flow, and micro-poromechanics. The stress-strain relationship of the dry fibrous media is first obtained by considering the fiber behavior. The constitutive relationships applicable to saturated media are then derived in the poromechanics framework using Hill's volume averaging. The advantage of this approach is that the resultant continuum model accounts for the discrete nature of the individual fibers while retaining a form suitable for porous materials. As a result, the model is able to predict the influence of micro-scale phenomena, such as the fiber pre-strain caused by osmotic effects and evolution of fiber network structure with loading, on the overall behavior and in particular, on the poromechanics parameters. Additionally, the model can describe fluid-flow related rate-dependent behavior under confined and unconfined conditions and varying chemical environments. The significance of the approach is demonstrated by simulating unconfined drained monotonic uniaxial compression under different surrounding fluid bath molarity, and fluid-flow related creep and relaxation at different loading-levels and different surrounding fluid bath molarity. The model predictions conform to the experimental observations for saturated soft fibrous materials. The method can potentially be extended to other porous materials such as bone, clays, foams and concrete.

Finite element prediction of transchondral stress and strain in the human hip

Cartilage fissures, surface fibrillation, and delamination represent early signs of hip osteoarthritis (OA). This damage may be caused by elevated first principal (most tensile) strain and maximum shear stress. The objectives of this study were to use a population of validated finite element (FE) models of normal human hips to evaluate the required mesh for converged predictions of cartilage tensile strain and shear stress, to assess the sensitivity to cartilage constitutive assumptions, and to determine the patterns of transchondral stress and strain that occur during activities of daily living. Five specimen-specific FE models were evaluated using three constitutive models for articular cartilage: quasilinear neo-Hookean, nonlinear Veronda Westmann, and tension-compression nonlinear ellipsoidal fiber distribution (EFD). Transchondral predictions of maximum shear stress and first principal strain were determined. Mesh convergence analysis demonstrated that five trilinear elements were adequate through the depth of the cartilage for precise predictions. The EFD model had the stiffest response with increasing strains, predicting the largest peak stresses and smallest peak strains. Conversely, the neo-Hookean model predicted the smallest peak stresses and largest peak strains. Models with neo-Hookean cartilage predicted smaller transchondral gradients of maximum shear stress than those with Veronda Westmann and EFD models. For FE models with EFD cartilage, the anterolateral region of the acetabulum had larger peak maximum shear stress and first principal strain than all other anatomical regions, consistent with observations of cartilage damage in disease. Results demonstrate that tension-compression nonlinearity of a continuous fiber distribution exhibiting strain induced anisotropy incorporates important features that have large effects on predictions of transchondral stress and strain. This population of normal hips provides baseline data for future comparisons to pathomorphologic hips. This approach can be used to evaluate these and other mechanical variables in the human hip and their potential role in the pathogenesis of osteoarthritis (OA).