Tensile behaviour of individual fibre bundles in the human lumbar anulus fibrosus

Non-enzymatic glycation increases the failure risk of annulus fibrosus by predisposing the extrafibrillar matrix to greater stresses

Growing clinical evidence suggests a correlation between diabetes and more frequent and severe intervertebral disc failure, partially attributed to accelerated advanced glycation end-products (AGE) accumulation in the annulus fibrosus (AF) through non-enzymatic glycation. However, in vitro glycation (i.e., crosslinking) reportedly improved AF uniaxial tensile mechanical properties, contradicting clinical observations. Thus, this study used a combined experimental-computational approach to evaluate the effect of AGEs on anisotropic AF tensile mechanics, applying finite element models (FEMs) to complement experimental testing and examine difficult-to-measure subtissue-level mechanics. Methylglyoxal-based treatments were applied to induce three physiologically relevant AGE levels in vitro. Models incorporated crosslinks by adapting our previously validated structure-based FEM framework. Experimental results showed that a threefold increase in AGE content resulted in a ∼55% increase in AF circumferential-radial tensile modulus and failure stress and a 40% increase in radial failure stress. Failure strain was unaffected by non-enzymatic glycation. Adapted FEMs accurately predicted experimental AF mechanics with glycation. Model predictions showed that glycation increased stresses in the extrafibrillar matrix under physiologic deformations, which may increase tissue mechanical failure or trigger catabolic remodeling, providing insight into the relationship between AGE accumulation and increased tissue failure. Our findings also added to the existing literature regarding crosslinking structures, indicating that AGEs had a greater effect along the fiber direction, while interlamellar radial crosslinks were improbable in the AF. In summary, the combined approach presented a powerful tool for examining multiscale structure-function relationships with disease progression in fiber-reinforced soft tissues, which is essential for developing effective therapeutic measures. STATEMENT OF SIGNIFICANCE: Increasing clinical evidence correlates diabetes with premature intervertebral disc failure, likely due to advanced glycation end-products (AGE) accumulation in the annulus fibrosus (AF). However, in vitro glycation reportedly increases AF tensile stiffness and toughness, contradicting clinical observations. Using a combined experimental-computational approach, our work shows that increases in AF bulk tensile mechanical properties with glycation are achieved at the risk of exposing the extrafibrillar matrix to increased stresses under physiologic deformations, which may increase tissue mechanical failure or trigger catabolic remodeling. Computational results indicate that crosslinks along the fiber direction account for 90% of the increased tissue stiffness with glycation, adding to the existing literature. These findings provide insight into the multiscale structure-function relationship between AGE accumulation and tissue failure.

On the identification of the ultra-structural organization of elastic fibers and their effects on the integrity of annulus fibrosus

Due to the complicated structure of the elastic fiber network in annulus fibrosus, existing in-silico studies either simplified or just overlooked its distribution pattern. Nonetheless, experimental and simulation results have proven that elastic fibers are of great importance to maintaining the structural integrity of annulus fibrosus and therefore to ensuring the load-bearing ability of intervertebral discs. Such needs call for a fine model. This work aims at developing a biphasic annulus fibrosus model by incorporating the accurate distribution pattern of collagen and elastic fibers. Both the structural parameters and intrinsic mechanical parameters were successfully identified using single lamella and inter-lamella microscopy anatomy and micromechanical testing data. The proposed model was then used to implement finite element simulations on various anterior and posterolateral multi-lamellae annulus fibrosus specimens. In general, simulation results agree well with available experimental and simulation data. On this basis, the effects of elastic fibers on the integrity of annulus fibrosus were further investigated. It was found that elastic fibers significantly influence the free swelling, radial stretching and circumferential shear performances of annulus fibrosus. Nonetheless, no significant effects were found for the circumferential stretching capability. The proposed biphasic model considers for the first time the distribution characteristics of elastic fibers at two scales, including both the principal orientations of all fiber families and the detailed distribution pattern within each family. Better understandings on the functions of collagen and elastic fibers can therefore be realized. To further enhance its prediction capability, the current model can be extended in the future by taking the fiber-matrix interaction as well as progressive damages into consideration.

A Novel Testing Method to Quantify Mechanical Properties of the Intact Annulus Fibrosus Ring From Rat-Tail Intervertebral Discs

The annulus fibrosus is the ring-like exterior of the intervertebral disc which is composed of concentrically organized layers of collagen fibre bundles. The mechanical properties of the annulus have been studied extensively; however, tests are typically performed on extracted fragments or multilayered samples of the annulus and not on the annulus as a whole. The purpose of this study was two-fold: 1) to develop a novel testing technique to measure the mechanical properties of the intact, isolated annulus; and 2) to perform a preliminary analysis of the rate-dependency of these mechanical properties. Twenty-nine whole annulus ring samples were dissected from 11 skeletally mature Sprague Dawley rat tails and underwent a tensile failure test at either 2%/s (n=16) or 20%/s (n=13). Force and displacement were sampled at 100Hz and were subsequently normalized to stress and strain. Various mechanical properties were derived from the stress-strain curves and statistically compared between the rates. All mechanical variables, with the exception of initial failure stress, were found to be unaffected by rate. Interestingly, initial failure stress was higher for samples tested at the slower rate compared to the higher rate which is atypical for viscoelastic tissues. Although in general rate did not appear to impact the annulus ring response to tensile loading, this novel, intact annular ring testing technique provides an alternative way to quantify mechanical properties of the annulus.

Detailed mechanical characterization of the transition zone: New insight into the integration between the annulus and nucleus of the intervertebral disc

The Nucleus Pulposus (NP) and Annulus Fibrous (AF) are two primary regions of the intervertebral disc (IVD). The interface between the AF and NP, where the gradual transition in structure and type of fibers are observed, is known as the Transition Zone (TZ). Recent structural studies have shown that the TZ contains organized fibers that appear to connect the NP to the AF. However, the mechanical characteristics of the TZ are yet to be explored. The current study aimed to investigate the mechanical properties of the TZ at the anterolateral (AL) and posterolateral (PL) regions in both radial and circumferential directions of loading using ovine IVDs (N = 28). Young's and toe moduli, maximum stress, failure strain, strain at maximum stress, and toughness were calculated mechanical parameters. The findings from this study revealed that the mechanical properties of the TZ, including young's modulus (p = 0.001), failure strain (p < 0.001), strain at maximum stress (p = 0.002), toughness (p = 0.027), and toe modulus (p = 0.005), were significantly lower for the PL compared to the AL region. Maximum stress was not significantly different between the PL and AL regions (p = 0.164). We found that maximum stress (p = 0.002), failure strain (p < 0.001), and toughness (p = 0.001) were significantly different in different loading directions. No significant differences for modulus (young's; p = 0.169 and toe; p = 0.352) and strain at maximum stress (p = 0.727) were found between the radial and circumferential loading directions. STATEMENT OF SIGNIFICANCE: To date there has not been a study that has investigated the mechanical characterization of the annulus (AF)-nucleus (NP) interface (transition zone; TZ) in the intervertebral disc (IVD), nor is it known whether the posterolateral (PL) and anterolateral (AL) regions of the TZ exhibit different mechanical properties. Accordingly, the TZ mechanical properties have been rarely used in the development of computational IVD models and relevant tissue-engineered scaffolds. The current research reported the mechanical properties of the TZ region and revealed that its mechanical properties were significantly lower for the PL compared to the AL region. These new findings enhance our knowledge about the nature of AF-NP integration and may help to develop more realistic tissue-engineered or computational IVD models.

Saline-polyethylene glycol blends preserve in vitro annulus fibrosus hydration and mechanics: An experimental and finite-element analysis

Precise control of tissue water content is essential for ensuring accurate, repeatable, and physiologically relevant measurements of tissue mechanics and biochemical composition. While previous studies have found that saline and polyethylene glycol (PEG) blends were effective at controlling tendon and ligament hydration levels, this work has yet to be extended to the annulus fibrosus (AF). Thus, the first objective of this study was to determine and validate an optimal buffer solution for targeting and maintaining hydration levels of tissue-level AF specimens in vitro. This was accomplished by measuring the transient swelling behavior of bovine AF specimens in phosphate-buffered saline (PBS) and PEG buffers across a wide range of concentrations. Sub-failure, failure, and post-failure mechanics were measured to determine the relationship between changes in tissue hydration and tensile mechanical response. The effect of each buffer solution on tissue composition was also assessed. The second objective of this study was to assess the feasibility and effectiveness of using multi-phasic finite element models to investigate tissue swelling and mechanical responses in different external buffer solutions. A solution containing 6.25%w/v PBS and 6.25%w/v PEG effectively maintained tissue-level AF specimen hydration at fresh-frozen levels after 18 h in solution. Modulus, failure stress, failure strain, and post-failure toughness of specimens soaked in this solution for 18 h closely matched those of fresh-frozen specimens. In contrast, specimens soaked in 0.9%w/v PBS swelled over 100% after 18 h and exhibited significantly diminished sub-failure and failure properties compared to fresh-frozen controls. The increased cross-sectional area with swelling contributed to but was not sufficient to explain the diminished mechanics of PBS-soaked specimens, suggesting additional sub-tissue scale mechanisms. Computational simulations of these specimens generally agreed with experimental results, highlighting the feasibility and importance of including a fluid-phase description when models aim to provide accurate predictions of biological tissue responses. As numerous previous studies suggest that tissue hydration plays a central role in maintaining proper mechanical and biological function, robust methods for controlling hydration levels are essential as the field advances in probing the relationship between tissue hydration, aging, injury, and disease.

Fiber engagement accounts for geometry-dependent annulus fibrosus mechanics: A multiscale, Structure-Based Finite Element Study

A comprehensive understanding of biological tissue mechanics is crucial for designing engineered tissues that aim to recapitulate native tissue behavior. Tensile mechanics of many fiber-reinforced tissues have been shown to depend on specimen geometry, which makes it challenging to compare data between studies. In this study, a validated multiscale, structure-based finite element model was used to evaluate the effect of specimen geometry on multiscale annulus fibrosus tensile mechanics through a fiber engagement analysis. The relationships between specimen geometry and modulus, Poisson's ratio, tissue stress-strain distributions, and fiber reorientation behaviors were investigated at both tissue and sub-tissue levels. It was observed that annulus fibrosus tissue level tensile properties and stress transmission mechanisms were dependent on specimen geometry. The model also demonstrated that the contribution of fiber-matrix interactions to tissue mechanical response was specimen size- and orientation-dependent. The results of this study reinforce the benefits of structure-based finite element modeling in studies investigating multiscale tissue mechanics. This approach also provides guidelines for developing optimal combined computational-experimental study designs for investigating fiber-reinforced biological tissue mechanics. Additionally, findings from this study help explain the geometry dependence of annulus fibrosus tensile mechanics previously reported in the literature, providing a more fundamental and comprehensive understanding of tissue mechanical behavior. In conclusion, the methods presented here can be used in conjunction with experimental tissue level data to simultaneously investigate tissue and sub-tissue scale mechanics, which is important as the field of soft tissue biomechanics advances toward studies that focus on diminishing length scales.

New findings confirm the viscoelastic behaviour of the inter-lamellar matrix of the disc annulus fibrosus in radial and circumferential directions of loading

While few studies have improved our understanding of composition and organization of elastic fibres in the inter-lamellar matrix (ILM), its clinical relevance is not fully understood. Moreover, no studies have measured the direct tensile and shear failure and viscoelastic properties of the ILM. Therefore, the aim of this study was, for the first time, to measure the viscoelastic and failure properties of the ILM in both the tension and shear directions of loading. Using an ovine model, isolated ILM samples were stretched to 40% of their initial length at three strain rates of 0.1%s^-1 (slow), 1%s^-1 (medium) and 10%s^-1 (fast) and a ramp test to failure was performed at a strain rate of 10%s^-1. The findings from this study identified that the stiffness of the ILM was significantly larger at faster strain rates, and energy absorption significantly smaller, compared to slower strain rates, and the viscoelastic and failure properties were not significantly different under tension and shear loading. We found a strain rate dependent response of the ILM during dynamic loading, particularly at the fastest rate. The ILM demonstrated a significantly higher capability for energy absorption at slow strain rates compared to medium and fast strain rates. A significant increase in modulus was found in both loading directions and all strain rates, having a trend of larger modulus in tension and at faster strain rates. The finding of no significant difference in failure properties in both loading directions, was consistent with our previous ultra-structural studies that revealed a well-organized (±45°) elastic fibre orientation in the ILM. The results from this study can be used to develop and validate finite element models of the AF at the tissue scale, as well as providing new strategies for fabricating tissue engineered scaffolds.

STATEMENT OF SIGNIFICANCE

While few studies have improved our understanding of composition and organization of elastic fibres in the inter-lamellar matrix (ILM) of the annulus in the disc no studies have measured the direct mechanical failure and viscoelastic properties of the ILM. The findings from this study identified that the stiffness of the ILM was significantly larger at faster strain rates, and energy absorption significantly smaller, compared to slower strain rates. The failure properties of the ILM were not significantly different under tension and shear.