A poroelastic finite element model of the bone-cartilage unit to determine the effects of changes in permeability with osteoarthritis.

The changes experienced in synovial joints with osteoarthritis involve coupled chemical, biological, and mechanical processes. The aim of this study was to investigate the consequences of increasing permeability in articular cartilage (AC), calcified cartilage (CC), subchondral cortical bone (SCB), and subchondral trabecular bone (STB) as observed with osteoarthritis. Two poroelastic finite element models were developed using a depth-dependent anisotropic model of AC with strain-dependent permeability and poroelastic models of calcified tissues (CC, SCB, and STB). The first model simulated a bone-cartilage unit (BCU) in uniaxial unconfined compression, while the second model simulated spherical indentation of the AC surface. Results indicate that the permeability of AC is the primary determinant of the BCU's poromechanical response while the permeability of calcified tissues exerts no appreciable effect on the force-indentation response of the BCU. In spherical indentation simulations with osteoarthritic permeability properties, fluid velocities were larger in magnitude and distributed over a smaller area compared to normal tissues. In vivo, this phenomenon would likely lead to chondrocyte death, tissue remodeling, alterations in joint lubrication, and the progression of osteoarthritis. For osteoarthritic and normal tissue permeability values, fluid flow was predicted to occur across the osteochondral interface. These results help elucidate the consequences of increases in the permeability of the BCU that occur with osteoarthritis. Furthermore, this study may guide future treatments to counteract osteoarthritis.

An evolutionary model of osteoarthritis including articular cartilage damage, and bone remodeling in a computational study.

With osteoarthritis, a complex set of progressive chemical, biological, and mechanical changes occur in both cartilage and bone. The aim of this study is to develop a high-fidelity computational model of the complete bone-cartilage unit to study the evolution of osterarthritis-induced articular cartilage (AC) damage and remodeling of subchondral cortical bone (SCB) and subchondral trabecular bone (STB). A finite element model of spherical indentation was developed with a depth-dependent anisotropic model of degenerating articular cartilage, a calcified cartilage (CC) zone, and SCB and STB remodeling regions. Calcified tissue (CC, SCB, and STB) and AC material regions were integrated to form an evolutionary bone-cartilage unit model. Results indicate that with indentation loading, articular cartilage damage occurs at the articular surface. Furthermore, bone remodeling was predicted to occur with a net stiffening of the subchondral bone plate. Changes in indentation force were minimal (<2%) between initial and final peak indentation loading. However, additional degradation and wear of AC and/or alterations in loading may have more pronounced effects on the mechanical response of the bone-cartilage unit. Bone remodeling and articular cartilage damage predictions are consistent with experimental observations that cartilage damage begins at the articular surface and subchondral bone experiences a thickening (i.e., stiffening) response with osteoarthritis. Our results provide insight into the early-term initiation behavior of osteoarthritis; the potential consequences of evolutions in AC, SCB, and STB with disease progression; and may guide future experimental and computational studies to elucidate mechanisms of osteoarthritis progression.

An Equilibrium Constitutive Model of Anisotropic Cartilage Damage to Elucidate Mechanisms of Damage Initiation and Progression.

Traumatic injuries and gradual wear-and-tear of articular cartilage (AC) that can lead to osteoarthritis (OA) have been hypothesized to result from tissue damage to AC. In this study, a previous equilibrium constitutive model of AC was extended to a constitutive damage articular cartilage (CDAC) model. In particular, anisotropic collagen (COL) fibril damage and isotropic glycosaminoglycan (GAG) damage were considered in a 3D formulation. In the CDAC model, time-dependent effects, such as viscoelasticity and poroelasticity, were neglected, and thus all results represent the equilibrium response after all time-dependent effects have dissipated. The resulting CDAC model was implemented in two different finite-element models. The first simulated uniaxial tensile loading to failure, while the second simulated spherical indentation with a rigid indenter displaced into a bilayer AC sample. Uniaxial tension to failure simulations were performed for three COL fibril Lagrangian failure strain (i.e., the maximum elastic COL fibril strain) values of 15%, 30%, and 45%, while spherical indentation simulations were performed with a COL fibril Lagrangian failure strain of 15%. GAG damage parameters were held constant for all simulations. Our results indicated that the equilibrium postyield tensile response of AC and the macroscopic tissue failure strain are highly dependent on COL fibril Lagrangian failure strain. The uniaxial tensile response consisted of an initial nonlinear ramp region due to the recruitment of intact fibrils followed by a rapid decrease in tissue stress at initial COL fibril failure, as a result of COL fibril damage which continued until ultimate tissue failure. In the spherical indentation simulation, damage to both the COL fibril and GAG constituents was located only in the superficial zone (SZ) and near the articular surface with tissue thickening following unloading. Spherical indentation simulation results are in agreement with published experimental observations. Our results indicate that the proposed CDAC model is capable of simulating both initial small magnitude damage as well as complete failure of AC tissue. The results of this study may help to elucidate the mechanisms of AC tissue damage, which initiate and propagate OA.

A Hertzian contact mechanics based formulation to improve ultrasound elastography assessment of uterine cervical tissue stiffness.

Clinical practice requires improved techniques to assess human cervical tissue properties, especially at the internal os, or orifice, of the uterine cervix. Ultrasound elastography (UE) holds promise for non-invasively monitoring cervical stiffness throughout pregnancy. However, this technique provides qualitative strain images that cannot be linked to a material property (e.g., Young's modulus) without knowledge of the contact pressure under a rounded transvaginal transducer probe and correction for the resulting non-uniform strain dissipation. One technique to standardize elastogram images incorporates a material of known properties and uses one-dimensional, uniaxial Hooke's law to calculate Young's modulus within the compressed material half-space. However, this method does not account for strain dissipation and the strains that evolve in three-dimensional space. We demonstrate that an analytical approach based on 3D Hertzian contact mechanics provides a reasonable first approximation to correct for UE strain dissipation underneath a round transvaginal transducer probe and thus improves UE-derived estimates of tissue modulus. We validate the proposed analytical solution and evaluate sources of error using a finite element model. As compared to 1D uniaxial Hooke's law, the Hertzian contact-based solution yields significantly improved Young's modulus predictions in three homogeneous gelatin tissue phantoms possessing different moduli. We also demonstrate the feasibility of using this technique to image human cervical tissue, where UE-derived moduli estimations for the uterine cervix anterior lip agreed well with published, experimentally obtained values. Overall, UE with an attached reference standard and a Hertzian contact-based correction holds promise for improving quantitative estimates of cervical tissue modulus.

Integrating qPLM and biomechanical test data with an anisotropic fiber distribution model and predictions of TGF-ß1 and IGF-1 regulation of articular cartilage fiber modulus.

A continuum mixture model with distinct collagen (COL) and glycosaminoglycan elastic constituents was developed for the solid matrix of immature bovine articular cartilage. A continuous COL fiber volume fraction distribution function and a true COL fiber elastic modulus ([Formula: see text] were used. Quantitative polarized light microscopy (qPLM) methods were developed to account for the relatively high cell density of immature articular cartilage and used with a novel algorithm that constructs a 3D distribution function from 2D qPLM data. For specimens untreated and cultured in vitro, most model parameters were specified from qPLM analysis and biochemical assay results; consequently, [Formula: see text] was predicted using an optimization to measured mechanical properties in uniaxial tension and unconfined compression. Analysis of qPLM data revealed a highly anisotropic fiber distribution, with principal fiber orientation parallel to the surface layer. For untreated samples, predicted [Formula: see text] values were 175 and 422 MPa for superficial (S) and middle (M) zone layers, respectively. TGF-[Formula: see text]1 treatment was predicted to increase and decrease [Formula: see text] values for the S and M layers to 281 and 309 MPa, respectively. IGF-1 treatment was predicted to decrease [Formula: see text] values for the S and M layers to 22 and 26 MPa, respectively. A novel finding was that distinct native depth-dependent fiber modulus properties were modulated to nearly homogeneous values by TGF-[Formula: see text]1 and IGF-1 treatments, with modulated values strongly dependent on treatment.

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

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