Interlamellar matrix governs human annulus fibrosus multiaxial behavior

A data-driven microstructure-based model for predicting circumferential behavior and failure in degenerated human annulus fibrosus

The degeneration of the intervertebral disc annulus fibrosus poses significant challenges in understanding and predicting its mechanical behavior. In this article, we present a novel approach, enriched with detailed insights into microstructure and degeneration progression, to accurately predict the mechanics of the degenerated human annulus. Central to this framework is a fully three-dimensional continuum-based model that integrates hydration state and multiscale structural features, including proteoglycan macromolecules and interpenetrating collagen fibrillar networks across various hierarchical levels within the multi-layered lamellar/inter-lamellar soft tissue, capable of sustaining deformation-induced damage. To ensure accurate and comprehensive predictions of the degenerated annulus mechanical behavior, we establish a data-driven correlation between disc degeneration grade and individual age, which influences the composition and mechanical integrity of annulus constituents while accounting for regional variations. The methodology includes a thorough identification of age- and grade-related evolutions of model inputs, followed by a detailed quantitative evaluation of the model predictive capabilities, with a focus on circumferential behavior and failure. The model successfully replicates experimental data, accurately capturing stiffness, transverse response (Poisson's ratio), and ultimate properties across different annulus regions, while also accommodating the modulation of the age/grade relationship. The reduction rates between normal and severe degeneration align reasonably well with experimental data, with the inner region exhibiting the largest decrease in stiffness (34.63%) and no significant change observed in the outer region. Failure stress drops considerably in both regions (49.86% in the inner and 45.33% in the outer), while failure strain decreases by 36.39% in the outer and 24.74% in the inner. Our findings demonstrate that the proposed framework significantly enhances the predictive accuracy of annulus mechanics across a spectrum of degeneration levels, from normal to severely degenerated states. This approach promises improved predictive accuracy, deeper insights into disc health and injury risk, and a robust foundation for further research on the impact of degeneration on disc integrity. STATEMENT OF SIGNIFICANCE: Understanding and predicting the mechanical behavior of degenerated human annulus fibrosus remains a significant challenge due to the complex interplay of structural, biochemical, and age-related factors. This study presents a microstructure-based approach to address this challenge by integrating hydration state, detailed structural features across hierarchical scales, and deformation-induced damage and failure, alongside age-related changes and degeneration grade factors. This approach enables accurate simulations of annulus mechanics across regions, with model results thoroughly compared to available data, reinforcing its applicability in capturing degeneration effects. By capturing the intricate interactions between microstructure and mechanical behavior in degenerated discs, the model lays a strong foundation for improving clinical assessments and guiding future treatment strategies for disc-related conditions.

Poisson function and volume ratio of soft anisotropic materials under large deformations

Accurate transverse deformation measurements are required for the estimation of the Poisson function and volume ratio. In this study, pure silicone and soft composite specimens were subjected to uniaxial tension, and the digital image correlation method was used to measure longitudinal and in- and out-of-plane transverse stretches. To minimize the effects of measurement errors on parameter estimation, the measured transverse stretches were defined in terms of the longitudinal stretch using a new formulation based on Poisson's ratios and two stretch-dependent parameters. From this formulation, Poisson functions and volume ratio for soft materials under large deformations were obtained. The results showed that pure silicone can be considered isotropic and nearly incompressible under large deformations, as expected. In contrast, Poisson's ratio of silicone reinforced with extensible fabric can exceed classical bounds, including negative value (auxetic behavior). The incompressibility assumption can be employed for describing the stress-stretch curve of pure silicone, while volume ratios are required for soft composites. Data of human skin, aortic wall, and annulus fibrosus from the literature were selected and analyzed. Except for the aortic wall, which can be considered nearly incompressible, the studied soft tissues must be regarded as compressible. All tissues presented anisotropic behavior.

Comprehensive analysis of damage evolution in human annulus fibrosus: Numerical exploration of mechanical sensitivity to biological age-dependent alteration

BACKGROUND AND OBJECTIVE

The annulus fibrosus is an essential part of the intervertebral disc, critical for its structural integrity. Mechanical deterioration in this component can lead to complete disc failure, particularly through tears development, with radial tears being the most common. These tears are often the result of both mechanical and biological factors. This study aims to numerically investigate the mechanisms of radial failure in the annulus tissue, taking into account the mechanical and age-dependent biological damage origins. A newly developed microstructure-based model was upgraded to predict damage evolution in the different annulus regions.

METHODS

The study employs a computational model to predict mechanical failures in various annulus regions, using experimental data for comparison. The model incorporates age-dependent microstructural changes to evaluate the effects of biological aging on the mechanical behavior. It specifically includes a detailed analysis of the temporal changes in circumferential rigidity and failure strain of the annulus.

RESULTS

The model demonstrated a strong ability to replicate the experimental responses of the different annulus regions to failure. It revealed that age-related microstructural changes significantly impact the rigidity and failure response of the annulus, particularly in the posterior regions and as well the anterior inner side. These changes increase susceptibility to rupture with aging. A correlation was also observed between the composition of collagen fibers, water content, and the annulus transversal response in both radial and axial directions.

CONCLUSION

The findings challenge previous assumptions, showing that age-dependent microstructural changes have a notable effect on the annulus mechanical properties. The computational model closely aligns with experimental observations, underscoring the determinant role of oriented collagen fibers in radial failure. This study enhances the understanding of annulus failure and provides a foundation for further research on the impact of aging on disc mechanical integrity and failure.

A Microstructure-Based Mechanistic Approach to Detect Degeneration Effects on Potential Damage Zones and Morphology of Young and Old Human Intervertebral Discs

There is an increasing demand to develop predictive medicine through the creation of predictive models and digital twins of the different body organs. To obtain accurate predictions, real local microstructure, morphology changes and their accompanying physiological degenerative effects must be taken into account. In this article, we present a numerical model to estimate the long-term aging effect on the human intervertebral disc response by means of a microstructure-based mechanistic approach. It allows to monitor in-silico the variations in disc geometry and local mechanical fields induced by age-dependent long-term microstructure changes. Both lamellar and interlamellar zones of the disc annulus fibrosus are constitutively represented by considering the main underlying microstructure features in terms of proteoglycans network viscoelasticity, collagen network elasticity (along with content and orientation) and chemical-induced fluid transfer. With age, a noticeable increase in shear strain is especially observed in the posterior and lateral posterior regions of the annulus which is in correlation with the high vulnerability of elderly people to back problems and posterior disc hernia. Important insights about the relation between age-dependent microstructure features, disc mechanics and disc damage are revealed using the present approach. These numerical observations are hardly obtainable using current experimental technologies which makes our numerical tool useful for patient-specific long-term predictions.

Modeling multiaxial damage regional variation in human annulus fibrosus

In the present article, a fully three-dimensional human annulus fibrosus model is developed by considering the regional variation of the complex structural organization of collagen network at different scales to predict the regional anisotropic multiaxial damage of the intervertebral disc. The model parameters are identified using experimental data considering as elementary structural unit, the single annulus lamellae stretched till failure along the micro-sized collagen fibers. The multi-layered lamellar/inter-lamellar annulus model is constructed by considering the effective interactions between adjacent layers and the chemical-induced volumetric strain. The regional dependent model predictions are analyzed under various loading modes and compared to experimental data when available. The stretching along the circumferential and radial directions till failure serves to check the predictive capacities of the annulus model. Model results under simple shear, biaxial stretching and plane-strain compression are further presented and discussed. Finally, a full disc model is constructed using the regional annulus model and simulations are presented to assess the most likely failed areas under disc axial compression. STATEMENT OF SIGNIFICANCE: The damage in annulus soft tissues is a complex multiscale phenomenon due to a complex structural arrangement of collagen network at different scales of hierarchical organization. A fully three-dimensional constitutive representation that considers the regional variation of the structural complexity to estimate annulus multiaxial mechanics till failure has not yet been developed. Here, a model is developed to predict deformation-induced damage and failure of annulus under multiaxial loading histories considering as time-dependent physical process both chemical-induced volumetric effects and damage accumulation. After model identification using single lamellae extracted from different disc regions, the model predictability is verified for various multiaxial elementary loading modes representative of the spine movement. The heterogeneous mechanics of a full human disc model is finally presented.

Angle-ply scaffold supports annulus fibrosus matrix expression and remodeling by mesenchymal stromal and annulus fibrosus cells

The angle-ply multilaminate structure of the annulus fibrosus is not reestablished following discectomy which leads to reherniation of the intervertebral disc (IVD). Biomimetic scaffolds developed to repair these defects should be evaluated for their ability to support tissue regeneration by endogenous and exogenous cells. Herein a collagen-based, angle-ply multilaminate patch designed to repair the outer annulus fibrosus was assessed for its ability to support mesenchymal stromal and annulus fibrosus cell viability, elongation, alignment, extracellular matrix gene expression, and scaffold remodeling. Results demonstrated that the cells remained viable, elongated, and aligned along the collagen fiber preferred direction of the scaffold, upregulated genes associated with annulus fibrosus matrix and produced collagen on the scaffold yielding biaxial mechanical properties that resembled native annulus fibrosus tissue. In conclusion, these scaffolds have demonstrated their potential to promote a living repair of defects in the annulus fibrosus and thus may be used to prevent recurrent IVD herniations.

A microstructure-based model for a full lamellar-interlamellar displacement and shear strain mapping inside human intervertebral disc core

The determinant role of the annulus fibrosus interlamellar zones in the intervertebral disc transversal and volumetric responses and hence on their corresponding three-dimensional conducts have been only revealed and appreciated recently. Their consideration in disc modeling strategies has been proven to be essential for the reproduction of correct local strain and displacement fields inside the disc especially in the unconstrained directions of the disc. In addition, these zones are known to be the starting areas of annulus fibrosus circumferential tears and disc delamination failure mode, which is often judged as one of the most dangerous disc failure modes that could evolve with time leading to disc hernia. For this latter reason, the main goal of the current contribution is to incorporate physically for the first time, the interlamellar zones, at the scale of a complete human lumbar intervertebral disc, in order to allow a correct local vision and replication of the different lamellar-interlamellar interactions and an identification of the interlamellar critical zones. By means of a fully tridimensional chemo-viscoelastic constitutive model, which we implemented into a finite element code, the physical, mechanical and chemical contribution of the interlamellar zones is added to the disc. The chemical-induced volumetric response is accounted by the model for both the interlamellar zones and the lamellae using experimentally-based fluid kinetics. Computational simulations are performed and critically discussed upon different simple and complex physiological movements. The disc core and the interlamellar zones are numerically accessed, allowing the observation of the displacement and shear strain fields that are compared to direct MRI experiments from the literature. Important conclusions about the correct lamellar-interlamellar-nucleus interactions are provided thanks to the developed model. The critical interlamellar spots with the highest delamination potentials are defined, analyzed and related to the local kinetics and microstructure.

Sensitivity of Intervertebral Disc Finite Element Models to Internal Geometric and Non-geometric Parameters

Finite element models are useful for investigating internal intervertebral disc (IVD) behaviours without using disruptive experimental techniques. Simplified geometries are commonly used to reduce computational time or because internal geometries cannot be acquired from CT scans. This study aimed to (1) investigate the effect of altered geometries both at endplates and the nucleus-anulus boundary on model response, and (2) to investigate model sensitivity to material and geometric inputs, and different modelling approaches (graduated or consistent fibre bundle angles and glued or cohesive inter-lamellar contact). Six models were developed from 9.4 T MRIs of bovine IVDs. Models had two variations of endplate geometry (a simple curved profile from the centre of the disc to the periphery, and precise geometry segmented from MRIs), and three variations of NP-AF boundary (linear, curved, and segmented). Models were subjected to axial compressive loading (to 0.86 mm at a strain rate of 0.1/s) and the effect on stiffness and strain distributions, and the sensitivity to modelling approaches was investigated. The model with the most complex geometry (segmented endplates, curved NP-AF boundary) was 3.1 times stiffer than the model with the simplest geometry (curved endplates, linear NP-AF boundary), although this difference may be exaggerated since segmenting the endplates in the complex geometry models resulted in a shorter average disc height. Peak strains were close to the endplates at locations of high curvature in the segmented endplate models which were not captured in the curved endplate models. Differences were also seen in sensitivity to material properties, graduated fibre angles, cohesive rather than glued inter-lamellar contact, and NP:AF ratios. These results show that FE modellers must take care to ensure geometries are realistic so that load is distributed and passes through IVDs accurately.