Enzymatic digestion of articular cartilage results in viscoelasticity changes that are consistent with polymer dynamics mechanisms

<p>Abstract</p> <p>Background</p> <p>Cartilage degeneration via osteoarthritis affects millions of elderly people worldwide, yet the specific contributions of matrix biopolymers toward cartilage viscoelastic properties remain unknown despite 30 years of research. Polymer dynamics theory may enable such an understanding, and predicts that cartilage stress-relaxation will proceed faster when the average polymer length is shortened.</p> <p>Methods</p> <p>This study tested whether the predictions of polymer dynamics were consistent with changes in cartilage mechanics caused by enzymatic digestion of specific cartilage extracellular matrix molecules. Bovine calf cartilage explants were cultured overnight before being immersed in type IV collagenase, bacterial hyaluronidase, or control solutions. Stress-relaxation and cyclical loading tests were performed after 0, 1, and 2 days of incubation.</p> <p>Results</p> <p>Stress-relaxation proceeded faster following enzymatic digestion by collagenase and bacterial hyaluronidase after 1 day of incubation (both <it>p </it>≤ 0.01). The storage and loss moduli at frequencies of 1 Hz and above were smaller after 1 day of digestion by collagenase and bacterial hyaluronidase (all <it>p </it>≤ 0.02).</p> <p>Conclusion</p> <p>These results demonstrate that enzymatic digestion alters cartilage viscoelastic properties in a manner consistent with polymer dynamics mechanisms. Future studies may expand the use of polymer dynamics as a microstructural model for understanding the contributions of specific matrix molecules toward tissue-level viscoelastic properties.</p

Identification of material parameters based on Mohr-Coulomb failure criterion for bisphosphonate treated canine vertebral cancellous bone

Nanoindentation has been widely used to study bone tissue mechanical properties. The common method and equations for analyzing nanoindentation, developed by Oliver and Pharr, are based on the assumption that the material is linearly elastic. In the present study, we adjusted the constraint of linearly elastic behavior and use nonlinear finite element analysis to determine the change in cancellous bone material properties caused by bisphosphonate treatment, based on an isotropic form of the Mohr–Coulomb failure model. Thirty-three canine lumbar vertebrae were used in this study. The dogs were treated daily for 1 year with oral doses of alendronate, risedronate, or saline vehicle at doses consistent, on a mg/kg basis, to those used clinically for the treatment of post-menopausal osteoporosis. Two sets of elastic modulus and hardness values were calculated for each specimen using the Continuous Stiffness Measurement (CSM) method (ECSM and HCSM) from the loading segment and the Oliver–Pharr method (EO–P and HO–P) from the unloading segment, respectively. Young's modulus (EFE), cohesion (c), and friction angle (ϕ) were identified using a finite element model for each nanoindentation. The bone material properties were compared among groups and between methods for property identification. Bisphosphonate treatment had a significant effect on several of the material parameters. In particular, Oliver–Pharr hardness was larger for both the risedronate- and alendronate-treated groups compared to vehicle and the Mohr–Coulomb cohesion was larger for the risedronate-treated compared to vehicle. This result suggests that bisphosphonate treatment increases the hardness and shear strength of bone tissue. Shear strength was linearly predicted by modulus and hardness measured by the Oliver–Pharr method (r2 = 0.99). These results show that bisphosphonate-induced changes in Mohr–Coulomb material properties, including tissue shear cohesive strength, can be accurately calculated from Oliver–Pharr measurements of Young's modulus and hardness

RANK/RANKL/OPG pathway: genetic associations with stress fracture period prevalence in elite athletes

Context: The RANK/RANKL/OPG signalling pathway is important in the regulation of bone turnover, with single nucleotide polymorphisms (SNPs) in genes within this pathway associated with bone phenotypic adaptations. Objective: To determine whether four SNPs associated with genes in the RANK/RANKL/OPG signalling pathway were associated with stress fracture injury in elite athletes. Design, Participants, and Methods: Radiologically confirmed stress fracture history was reported in 518 elite athletes, forming the Stress Fracture Elite Athlete (SFEA) cohort. Data were analysed for the whole group, and were sub-stratified into male and cases of multiple stress fracture group. Genotypes were determined using proprietary fluorescence-based competitive allele-specific PCR assays. Results: SNPs rs3018362 (RANK) and rs1021188 (RANKL) were associated with stress fracture injury (p<0.05). 8.1% of stress fracture group and 2.8% of the non-stress fracture group were homozygote for the rare allele of rs1021188. Allele frequency, heterozygotes and homozygotes for the rare allele of rs3018362 were associated with stress fracture period prevalence (p<0.05). Analysis of the male only group showed 8.2% of rs1021188 rare allele homozygotes to have suffered a stress fracture while 2.5% of the non-stress fracture group were homozygous. In cases of multiple stress fractures, homozygotes for the rare allele of rs1021188, and individuals possessing at least one copy of the rare allele of rs4355801 (OPG) were shown to be associated with stress fracture injury (p<0.05). Conclusions: The data support an association between SNPs in the RANK/RANKL/OPG signalling pathway and the development of stress fracture injury. The association of rs3018362 (RANK) and rs1021188 (RANKL) with stress fracture injury susceptibility supports their role in the maintenance of bone health, and offers potential targets for therapeutic interventions

A Novel Method for Curvefitting the Stretched Exponential Function to Experimental Data.

The stretched exponential function has many applications in modeling numerous types of experimental relaxation data. However, problems arise when using standard algorithms to fit this function: we have observed that different initializations result in distinct fitted parameters. To avoid this problem, we developed a novel algorithm for fitting the stretched exponential model to relaxation data. This method is advantageous both because it requires only a single adjustable parameter and because it does not require initialization in the solution space. We tested this method on simulated data and experimental stress-relaxation data from bone and cartilage and found favorable results compared to a commonly-used Quasi-Newton method. For the simulated data, strong correlations were found between the simulated and fitted parameters suggesting that this method can accurately determine stretched exponential parameters. When this method was tested on experimental data, high quality fits were observed for both bone and cartilage stress-relaxation data that were significantly better than those determined with the Quasi-Newton algorithm

Training drives turnover rates in racehorse proximal sesamoid bones

Focal bone lesions are often found prior to clinically relevant stress-fractures. Lesions are characterized by low bone volume fraction, low mineral density, and high levels of microdamage and are hypothesized to develop when bone tissue cannot sufficiently respond to damaging loading. It is difficult to determine how exercise drives the formation of these lesions because bone responds to mechanical loading and repairs damage. In this study, we derive steady-state rate constants for a compartment model of bone turnover using morphometric data from fractured and non-fractured racehorse proximal sesamoid bones (PSBs) and relate rate constants to racing-speed exercise data. Fractured PSBs had a subchondral focus of bone turnover and microdamage typical of lesions that develop prior to fracture. We determined steady-state model rate constants at the lesion site and an internal region without microdamage using bone volume fraction, tissue mineral density, and microdamage area fraction measurements. The derived undamaged bone resorption rate, damage formation rate, and osteoid formation rate had significant robust regression relationships to exercise intensity (rate) variables, layup (time out of exercise), and exercise 2-10&nbsp;months before death. However, the direction of these relationships varied between the damaged (lesion) and non-damaged regions, reflecting that the biological response to damaging-loading differs from the response to non-damaging loading

Exercise history predicts focal differences in bone volume fraction, mineral density and microdamage in the proximal sesamoid bones of Thoroughbred racehorses

Medial proximal sesamoid bones (PSBs) from Thoroughbred racehorses that did (Case) or did not (Control) experience unilateral biaxial PSB fracture were evaluated for bone volume fraction (BVF), apparent mineral density (AMD), tissue mineral density (TMD), and microdamage in Case fractured, Case contralateral limb intact, and Control bones. A majority of Case bones had a subchondral lesion with high microdamage density, and low BVF, AMD, and TMD. Lesion microdamage and densitometric measures were associated with training history by robust linear regression. Exercise intensity was negatively related to BVF (0.07 ≤ R2  ≤ 0.12) and positively related to microcrack areal density (0.21 ≤ R2  ≤ 0.29) in the lesion; however, in an undamaged site, the relationships were opposite in direction. Regardless of location, TMD decreased with event frequency for both Case and Control, suggesting increased bone remodeling with exercise. Measures of how often animals were removed from active training (layups) predicted a decrease in TMD, AMD, BVF, and microdamage at regions away from the lesion site. A steady-state compartment model was used to organize the differences in the correlations between variables within the data set. The overall conclusions are that at the osteopenic lesion site, repair of microdamage by remodeling was not successful (e.g., lower bone mass, increased damage, and lower mineralization) but that in regions away from the lesion remodeling successfully controlled damage (e.g., higher bone mass, less microdamage, and lower mineralization)

Cancellous Bone Properties and Matrix Content of TGF-β2 and IGF-I in Human Tibia: A Pilot Study

Transforming and insulin-like growth factors are important in regulating bone mass. Thus, one would anticipate correlations between matrix concentrations of growth factors and functional properties of bone. We therefore investigated the relationships of (1) TGF-β2 and (2) IGF-I matrix concentrations with the trabecular microstructure, stress distribution, and mechanical properties of tibial cancellous bone from six male human cadavers. Trabecular stress amplification (VMExp/σapp) and variability (VMCOV) were calculated using microcomputed tomography (μCT)-based finite element simulations. Bone volume fraction (BV/TV), surface/volume ratio (BS/BV), trabecular thickness (Tb.Th), number (Tb.N) and separation (Tb.Sp), connectivity (Eu.N), and anisotropy (DA) were measured using 3-D morphometry. Bone stiffness and strength were measured by mechanical testing. Matrix concentrations of TGF-β2 and IGF-I were measured by ELISA. We found higher matrix concentrations of TGF-β2 were associated with higher Tb.Sp and VMExp/σapp for pooled data and within subjects. Similarly, a higher matrix concentration of IGF-I was associated with lower stiffness, strength, BV/TV and Tb.Th and with higher BS/BV, Tb.Sp, VMExp/σapp and VMCOV for pooled data and within subjects. IGF-I and Tb.N were negatively associated within subjects. It appears variations of the stress distribution in cancellous bone correlate with the variation of the concentrations of TGF-β2 and IGF-I in bone matrix: increased local matrix concentrations of growth factors are associated with poor biomechanical and architectural properties of tibial cancellous bone