Search

A buffer solution is often used to maintain tissue hydration during mechanical testing. The most commonly used buffer solution is a physiological concentration of phosphate buffered saline (PBS); however, PBS increases the tissue’s water content and decreases its tensile stiffness. In addition, solutes from the buffer can diffuse into the tissue and interact with its structure and mechanics. These bathing solution effects can confound the outcome and interpretation of mechanical tests. Potential bathing solution artifacts, including solute diffusion and the effect on mechanical properties, are not well understood. The objective of this study was to measure the effects of long-term exposure of rat tail tendon fascicles to several concentrations (0.9% to 25%) of NaCl, sucrose, polyethylene glycol (PEG), and SPEG (NaCl + PEG) solutions on water content, solute diffusion, and mechanical properties. We found that with an increase in solute concentration the apparent water content decreased for all solution types. Solutes diffused into the tissue for NaCl and sucrose, however, no solute diffusion was observed for PEG or SPEG. The mechanical properties changed for both of NaCl solutions, in particular after long-term (8 hr) incubation the modulus and equilibrium stress decreased compared to short-term (15 min) for 25% NaCl, and the cross sectional area increased for 0.9% NaCl. However, the mechanical properties were unchanged for both PEG and SPEG except for minor alterations in stress relaxation parameters. This study shows that NaCl and sucrose buffer solutions are not suitable for long-term mechanical tests. We therefore propose using PEG or SPEG as alternative buffer solutions that after long-term incubation can maintain tissue hydration without solute diffusion and produce a consistent mechanical response.

Tendon’s viscoelastic behaviors are important to the mechanical function and mechanobiology. When loaded in longitudinal tension, tendons often have a very large Poisson’s ratio (ν > 2) that exceeds the limit of incompressibility for isotropic material (ν = 0.5), indicating that tendon experiences volume loss, inducing poroelastic fluid exudation in the transverse direction. Therefore, transverse poroelasticity is an important contributor to tendon material behavior. Tendon hydraulic permeability which is required to evaluate the fluid flow contribution to viscoelasticity, is mostly unavailable, and where available, varies by several orders of magnitude. In this manuscript, we quantified the transverse poroelastic material parameters of rat tail tendon fascicles by conducting transverse osmotic loading experiments, in both tension and compression. We used a multi-start optimization method to evaluate the parameters based on biphasic finite element modeling. Our tendon samples had a transverse hydraulic permeability of 10(−4) to 10(−5) mm(4) (Ns)(−1) and showed a significant tension-compression nonlinearity in the transverse direction. Further, using these results, we predict hydraulic permeability during longitudinal (fiber-aligned) tensile loading, and the spatial distribution of fluid flow during osmotic loading. These results reveal novel aspects of tendon mechanics and can be used to study the physiomechanical response of tendon in response to mechanical loading.

Intervertebral disc mechanics are affected by both disc shape and disc degeneration, which in turn each affect the other; disc mechanics additionally have a role in the etiology of disc degeneration. Finite element analysis (FEA) is a favored tool to investigate these relationships, but limited data for intervertebral disc 3D shape has forced the use of simplified or single-subject geometries, with the effect of inter-individual shape variation investigated only in specialized studies. Similarly, most data on disc shape variation with degeneration is based on 2D mid-sagittal images, which incompletely define 3D shape changes. Therefore, the objective of this study was to quantify inter-individual disc shape variation in 3D, classify this variation into independently-occurring modes using a statistical shape model, and identify correlations between disc shape and degeneration. Three-dimensional disc shapes were obtained from MRI of 13 human male cadaver L3L4 discs. An average disc shape and four major modes of shape variation (representing 90% of the variance) were identified. The first mode represented disc axial area and was significantly correlated to degeneration (R(2)=0.44), indicating larger axial area in degenerate discs. Disc height variation occurred in three distinct modes, each also involving non-height variation. The statistical shape model provides an average L3L4 disc shape for FEA that is fully defined in 3D, and makes it convenient to generate a set of shapes with which to represent aggregate inter-individual variation. Degeneration grade-specific shapes can also be generated. To facilitate application, the model is included in this paper׳s supplemental content.

The aim of functional tissue engineering is to repair and replace tissues that have a biomechanical function, i.e., connective orthopaedic tissues. To do this, it is necessary to have accurate benchmarks for the elastic, permeability, and swelling (i.e., biphasic-swelling) properties of native tissues. However, in the case of the intervertebral disc, the biphasic-swelling properties of individual tissues reported in the literature exhibit great variation and even span several orders of magnitude. This variation is probably caused by differences in the testing protocols and the constitutive models used to analyze the data. Therefore, the objective of this study was to measure the human lumbar disc annulus fibrosus (AF), nucleus pulposus (NP), and cartilaginous endplates (CEP) biphasic-swelling properties using a consistent experimental protocol and analyses. The testing protocol was composed of a swelling period followed by multiple confined compression ramps. To analyze the confined compression data, the tissues were modeled using a biphasic-swelling model, which augments the standard biphasic model through the addition of a deformation-dependent osmotic pressure term. This model allows considering the swelling deformations and the contribution of osmotic pressure in the analysis of the experimental data. The swelling stretch was not different between the disc regions (AF: 1.28±0.16; NP: 1.73±0.74; CEP: 1.29±0.26), with a total average of 1.42. The aggregate modulus (Ha) of the extra-fibrillar matrix was higher in the CEP (390kPa) compared to the NP (100kPa) or AF (30kPa). The permeability was very different across tissue regions, with the AF permeability (64 E(-16)m(4)/Ns) higher than the NP and CEP (~5.5 E(-16)m(4)/Ns). Additionally, a normalized time-constant (3000s) for the stress relaxation was similar for all the disc tissues. The properties measured in this study are important as benchmarks for tissue engineering and for modeling the disc's mechanical behavior and transport.

Due to the inability of current clinical practices to restore function to degenerated intervertebral discs, the arena of disc tissue engineering has received substantial attention in recent years. Despite tremendous growth and progress in this field, translation to clinical implementation has been hindered by a lack of well-defined functional benchmarks. Because successful replacement of the disc is contingent upon replication of some or all of its complex mechanical behaviour, it is critically important that disc mechanics be well characterized in order to establish discrete functional goals for tissue engineering. In this review, the key functional signatures of the intervertebral disc are discussed and used to propose a series of native tissue benchmarks to guide the development of engineered replacement tissues. These benchmarks include measures of mechanical function under tensile, compressive and shear deformations for the disc and its substructures. In some cases, important functional measures are identified that have yet to be measured in the native tissue. Ultimately, native tissue benchmark values are compared to measurements that have been made on engineered disc tissues, identifying measures where functional equivalence was achieved, and others where there remain opportunities for advancement. Several excellent reviews exist regarding disc composition and structure, as well as recent tissue engineering strategies; therefore this review will remain focused on the functional aspects of disc tissue engineering.

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.

A transversely isotropic biphasic mixture model was applied to tendon in uniaxial tension. Parametric analyses were performed and the sensitivity in predicting material parameters was evaluated. Our results provide quantitative evidence for fluid flow as a mechanism that contributes to tendon viscoelasticity. Transversely isotropic material properties were calculated for mouse tail tendon fascicles. The average transverse modulus (E1) was 0.046 MPa, the fiber-aligned Poisson's ratio (ν31) was 2.73, and the transverse Poisson's ratio (ν21) was 0.96; these properties were not strain-dependent. The fiber-aligned modulus (E3) was strain-dependent and was 20.7 MPa in the toe region and 86.1 MPa in the linear region. These solid matrix properties were consistent with previously published tendon tissue and fascicle data. The fascicle permeability was strain-dependent and was 5.5×10−18 m4/N s in the toe region and 0.32×10−18 m4/N s in the linear region, similar to previously reported meniscus permeability in tension. The similar permeabilities of both fascicle and tissue-level samples suggest that fluid flow from individual fascicles, not the packing of multiple fascicles together, may be the primary barrier to fluid flow in tendon and thus the primary mechanism for viscoelasticity.

The first objective of this study was to determine the effects of physiological cyclic loading followed by unloaded recovery on the mechanical response of human intervertebral discs. The second objective was to examine how nucleotomy alters the disc׳s mechanical response to cyclic loading. To complete these objectives, 15 human L5-S1 discs were tested while intact and subsequent to nucleotomy. The testing consisted of 10,000 cycles of physiological compressive loads followed by unloaded hydrated recovery. Cyclic loading increased compression modulus (3%) and strain (33%), decreased neutral zone modulus (52%), and increased neutral zone strain (31%). Degeneration was not correlated with the effect of cyclic loading in intact discs, but was correlated with cyclic loading effects after nucleotomy, with more degenerate samples experiencing greater increases in both compressive and neutral zone strain following cyclic loading. Partial removal of the nucleus pulposus decreased the compression and neutral zone modulus while increasing strain. These changes correspond to hypermobility, which will alter overall spinal mechanics and may impact low back pain via altered motion throughout the spinal column. Nucleotomy also reduced the effects of cyclic loading on mechanical properties, likely due to altered fluid flow, which may impact cellular mechanotransduction and transport of disc nutrients and waste. Degeneration was not correlated with the acute changes of nucleotomy. Results of this study provide an ideal protocol and control data for evaluating the effectiveness of a mechanically-based disc degeneration treatment, such as a nucleus replacement.

A recent study (Lake et al. 2009) reported the properties of human supraspinatus tendon (SST) tested along the predominant fiber direction. The SST was found to have a relatively disperse distribution of collagen fibers, which may represent an adaptation to multiaxial loads imposed by the complex loading environment of the rotator cuff. However, the multiaxial mechanical properties of human SST remain unknown. The objective of this study, therefore, was to evaluate the mechanical properties, fiber alignment, change in alignment with applied load, and structure-function relationships of SST in transverse testing. Samples from six SST locations were tested in uniaxial tension with samples oriented transverse to the tendon long-axis. Polarized light imaging was used to quantify collagen fiber alignment and change in alignment under applied load. The mechanical properties of samples taken near the tendon-bone insertion were much greater on the bursal surface compared to the joint surface (e.g., bursal moduli 15–30 times greater than joint; p

The angled, lamellar structure of the annulus fibrosus is integral to its load-bearing function. Reorientation of this fiber structure with applied load may contribute to nonlinear mechanical behavior and to large increases in tensile modulus. Fiber reorientation has not yet been quantified for loaded non-degenerated and degenerated annulus fibrosus tissue. The objective of this study was to measure fiber reorientation and mechanical properties (toe- and linear-region modulus, transition strain, and Poisson's ratio) of loaded outer annulus fibrosus tissue using a new application of FFT image processing techniques. This method was validated for quantification of annulus fiber reorientation during loading in this study. We hypothesized that annulus fibrosus fibers would reorient under circumferential tensile load, and that fiber reorientation would be affine. Additionally, we hypothesized that degeneration would affect fiber reorientation, toe-region modulus and Poisson's ratio. Annulus fibrosus fibers were found to reorient toward the loading direction, and degeneration significantly decreased fiber reorientation (the fiber reorientation parameter, m(FFT)=-1.70 degrees /% strain for non-degenerated and -0.95 degrees /% strain for degenerated tissue). Toe-region modulus was significantly correlated with age (r=0.6). Paired t-tests showed no significant difference in the fiber reorientation parameter calculated experimentally with that calculated using an affine prediction. Thus, an affine prediction is a good approximation of fiber reorientation. The findings of this study add to the understanding of overall disc mechanical behavior and degeneration.

The objective of this study was to use a homogenization model of the anisotropic mechanical behavior of annulus fibrosus (AF) to address some of the issues raised in structural finite element and fiber-reinforced strain energy models. Homogenization theory describes the effect of microstructure on macroscopic material properties by assuming the material is composed of repeating representative volume elements. We first developed the general homogenization model and then specifically prescribed the model to in-plane single lamella and multi-lamellae AF properties. We compared model predictions to experimentally measured AF properties and performed parametric studies. The predicted tensile moduli ( E θ and Ez) and their dependence on fiber volume fraction and fiber angle were consistent with measured values. However, the model prediction for shear modulus ( G θ z ) was two orders of magnitude larger than directly measured values. The values of E θ and Ez were strongly dependent on the model input for matrix modulus, much more so than the fiber modulus. These parametric analyses demonstrated the contribution of the matrix in AF load support, which may play a role when protoeglycans are decreased in disc degeneration, and will also be an important design factor in tissue engineering. We next compared the homogenization model to a 3-D structural finite element model and fiber-reinforced energy models. Similarities between the three model types provided confidence in the ability of these models to predict AF tissue mechanics. This study provides a direct comparison between the several types of AF models and will be useful for interpreting previous studies and elucidating AF structure–function relationships in disc degeneration and for functional tissue engineering.