Derivation of inter-lamellar behaviour of the intervertebral disc annulus

Influence of Radial Variations in Biochemical Concentrations in Collagen Type and Water on Mechanical Stability of Annulus Fibrosus' Collagen-Hyaluronan Interfaces at Nanoscale: A Molecular Dynamics Investigation

Multidirectional load transmission ability by annulus fibrosus (AF) requires substantial mechanical stability. Additionally, AF exhibits a unique biochemical concentration gradient with outer AF (OA) dominated by type I collagen (COL-I) and inner AF dominated by type II collagen (COL-II) with higher water and proteoglycan concentration. This indicates an intricate relationship between biochemistry and mechanical stability, which remains unclear. This study uses molecular dynamics (MD) simulations to investigate the impact of water, COL-I and COL-II, concentration gradients on mechanical stability of AF's collagen-hyaluronan (COL-HYL) nano-interfaces during tensile and compressive deformation. For this, COL-HYL atomistic models are created by increasing COL-II concentrations from 0% to 75% and water from 65% to 75%. Additional tensile and compressive deformation simulations are conducted for COL-I-HYL interface (COL-HYL interfaces with 0% COL-II) by increasing water concentration from 65% to 75% to segregate the effects of increasing water concentration alone. Results show that increasing water concentration alone to 75% results in marginal changes in local hydration indicating increase in bulk water. This enhances HYL and COL segment sliding-leading to reduction in mechanical stability in tension, indicated by drop in stress-strain characteristics. Additionally, increase in bulk water shifts load-bearing characteristics toward water-leading to reduction in modulus from 3.7 GPa to 1.9 GPa. Conversely, increasing COL-II and water concentration facilitates stable water bridge formation which impedes sliding in HYL and COL-enhancing mechanical stability. These water bridges further improve compressive load sustenance leading to lower reduction in compressive modulus from 3.7 GPa to 2.8 GPa.

Role of intra-lamellar collagen and hyaluronan nanostructures in annulus fibrosus on lumbar spine biomechanics: insights from molecular mechanics-finite element-based multiscale analyses

Annulus fibrosus' (AF) ability to transmit multi-directional spinal motion is contributed by a combination of chemical interactions among biomolecular constituents-collagen type I (COL-I), collagen type II (COL-II), and proteoglycans (aggrecan and hyaluronan)-and mechanical interactions at multiple length scales. However, the mechanistic role of such interactions on spinal motion is unclear. The present work employs a molecular mechanics-finite element (FE) multiscale approach to investigate the mechanistic role of molecular-scale collagen and hyaluronan nanostructures in AF, on spinal motion. For this, an FE model of the lumbar segment is developed wherein a multiscale model of AF collagen fiber, developed from COL-I, COL-II, and hyaluronan using the molecular dynamics-cohesive finite element multiscale method, is incorporated. Analyses show AF collagen fibers primarily contribute to axial rotation (AR) motion, owing to angle-ply orientation. Maximum fiber strain values of 2.45% in AR, observed at the outer annulus, are 25% lower than the reported values. This indicates native collagen fibers are softer, attributed to the softer non-fibrillar matrix and higher interfibrillar sliding. Additionally, elastic zone stiffness of 8.61 Nm/° is observed to be 20% higher than the reported range, suggesting native AF lamellae exhibit lower stiffness, resulting from inter-collagen fiber bundle sliding. The presented study has further implications towards the hierarchy-driven designing of AF-substitute materials.

Impact of extracellular matrix and collagen network properties on the cervical intervertebral disc response to physiological loads: A parametric study

Current intervertebral disc finite element models are hard to validate since they describe multi-physical phenomena and contain a huge number of material properties. This work aims to simplify numerical validation/identification studies by prioritizing the sensitivity of intervertebral disc behavior to mechanical properties. A 3D fiber-reinforced hyperelastic model of a C6-C7 intervertebral disc is used to carry out the parametric study. 10 parameters describing the extracellular matrix and the collagen network behaviors are included in the parametric study. The influence of varying these parameters on the disc response is estimated during physiological movements of the head, including compression, lateral bending, flexion, and axial rotation. The obtained results highlight the high sensitivity of the disc behavior to the stiffness of the annulus fibrosus extracellular matrix for all the studied loads with a relative increase in the disc apparent stiffness by 67% for compression and by 57% for axial rotation when the annulus stiffness increases from 0.4 to 2 MPa. It is also shown that varying collagen network orientation, stiffness, and stiffening in the studied configuration range have a noticeable effect on rotational motions with a relative apparent stiffness difference reaching 6.8%, 10%, and 22%, respectively, in lateral bending. However, the collagen orientation does not affect disc response to axial load.

Biomechanical effect of age-related structural changes on cervical intervertebral disc: A finite element study

Previous literature has investigated the biomechanical response of healthy and degenerative discs, but the biomechanical response of suboptimal healthy intervertebral discs received less attention. The purpose was to compare the biomechanical responses and risk of herniation of young healthy, suboptimal healthy, and degenerative intervertebral discs. A cervical spine model was established and validated using the finite element method. Suboptimal healthy, mildly, moderately, and severely degenerative disc models were developed. Disc height deformation, range of motion, intradiscal pressure, and von Mises stress in annulus fibrosus were analyzed by applying a moment of 4 Nm in flexion, extension, lateral bending, and axial rotation with 100 N compressive loads. Disc height deformation in young healthy, suboptimal healthy, mildly, moderately, and severely degenerative discs was 40%, 37%, 21%, 12%, and 8%, respectively. The decreasing order of the range of motion was young healthy spine > suboptimal healthy spine > mildly degenerative spine > moderately degenerative spine > severely degenerative spine. The mean stress of annulus ground substance in the suboptimal healthy disc was higher than in the young healthy disc. The mean stress of inter-lamellar matrix and annulus ground substance in moderately and severely degenerative discs was higher than in other discs. Age-related structural changes and degenerative changes increased the stiffness and reduced the elastic deformation of intervertebral discs. Decreased range of motion due to the effects of aging or degeneration on the intervertebral disc, may cause compensation of adjacent segments and lead to progressive degeneration of multiple segments. The effect of aging on the intervertebral disc increased the risk of annulus fibrosus damage from the biomechanical point of view. Moderately and severely degenerative discs may have a higher risk of herniation due to the higher risk of damage and layers separation of annulus fibrosus caused by increased stress in the annulus ground substance and inter-lamellar matrix.

Elastic Fibers in the Intervertebral Disc: From Form to Function and toward Regeneration

Despite extensive efforts over the past 40 years, there is still a significant gap in knowledge of the characteristics of elastic fibers in the intervertebral disc (IVD). More studies are required to clarify the potential contribution of elastic fibers to the IVD (healthy and diseased) function and recommend critical areas for future investigations. On the other hand, current IVD in-vitro models are not true reflections of the complex biological IVD tissue and the role of elastic fibers has often been ignored in developing relevant tissue-engineered scaffolds and realistic computational models. This has affected the progress of IVD studies (tissue engineering solutions, biomechanics, fundamental biology) and translation into clinical practice. Motivated by the current gap, the current review paper presents a comprehensive study (from the early 1980s to 2022) that explores the current understanding of structural (multi-scale hierarchy), biological (development and aging, elastin content, and cell-fiber interaction), and biomechanical properties of the IVD elastic fibers, and provides new insights into future investigations in this domain.

Using Statistical Shape and Appearance Modelling to characterise the 3D shape and material properties of human lumbar vertebrae: A proof of concept study

Patient variation affects the outcomes of a range of spinal interventions, from disc replacement to vertebral fixation and vertebroplasty. Statistical Shape and Appearance Modelling (SSAM) can be used to describe anatomical variation and pathological differences within the population. To better understand how bone density and shape variation affect load transfer with respect to surgical treatments, Finite Element (FE) models can be generated from a SSAM. The aim for this study is to understand whether geometric and density variation as well as multiple vertebral levels can be incorporated into a single SSAM and whether this can be used to investigate the relationships between, and effects of, the various modes of variation. FE models of 14 human lumbar vertebrae that had been μCT imaged and validated through experimental testing were used as input specimens for a SSAM. The validity of the SSAM was evaluated by using principal component analysis to identify the primary modes of geometric and bone density variation and comparing to those in the input set. FE models were generated from the SSAM to examine the response to loading. The mean error between the input set and generated models for volume, mean density and FE compressive stiffness were 10%, 3% and 10% respectively. Principal Component (PC) 1 captured the majority of the bone density variation. The remaining PCs described specific geometric variation. The FE models generated from the SSAM showed the variations in vertebral stiffness as a result of complex relationships between bone density and shape. The SSAM created has limited data for its input set, however, it acts as a proof of concept for the novel combination of material and shape variation into a single shape model. This approach and the tools developed can be applied to wider patient groups and treatment scenarios to improve patient stratification and to optimise treatments.

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.

Architecture-Promoted Biomechanical Performance-Tuning of Tissue-Engineered Constructs for Biological Intervertebral Disc Replacement

Background: Biological approaches to intervertebral disc (IVD) restoration and/or regeneration have become of increasing interest. However, the IVD comprises a viscoelastic system whose biological replacement remains challenging. The present study sought to design load-sharing two-component model systems of circular, nested, concentric elements reflecting the nucleus pulposus and annulus fibrosus. Specifically, we wanted to investigate the effect of architectural design variations on (1) model system failure loads when testing the individual materials either separately or homogeneously mixed, and (2) also evaluate the potential of modulating other mechanical properties of the model systems. Methods: Two sets of softer and harder biomaterials, 0.5% and 5% agarose vs. 0.5% agarose and gelatin, were used for fabrication. Architectural design variations were realized by varying ring geometries and amounts while keeping the material composition across designs comparable. Results: Variations in the architectural design, such as lamellar width, number, and order, combined with choosing specific biomaterial properties, strongly influenced the biomechanical performance of IVD constructs. Biomechanical characterization revealed that the single most important parameter, in which the model systems vastly exceeded those of the individual materials, was failure load. The model system failure loads were 32.21- and 84.11-fold higher than those of the agarose materials and 55.03- and 2.14-fold higher than those of the agarose and gelatin materials used for system fabrication. The compressive strength, dynamic stiffness, and viscoelasticity of the model systems were always in the range of the individual materials. Conclusions: Relevant architecture-promoted biomechanical performance-tuning of tissue-engineered constructs for biological IVD replacement can be realized by slight modifications in the design of constructs while preserving the materials’ compositions. Minimal variations in the architectural design can be used to precisely control structure–function relations for IVD constructs rather than choosing different materials. These fundamental findings have important implications for efficient tissue-engineering of IVDs and other load-bearing tissues, as potential implants need to withstand high in situ loads.

A microstructure-based modeling approach to assess aging-sensitive mechanics of human intervertebral disc

BACKGROUND AND OBJECTIVE

The human body soft tissues are hierarchic structures interacting in a complex manner with the surrounding biochemical environment. The loss of soft tissues functionality with age leads to more vulnerability regarding to the external mechanical loadings and increases the risk of injuries. As a main example of the human body soft tissues, the intervertebral disc mechanical response evolution with age is explored. Although the age-dependence of the intervertebral disc microstructure is a well-known feature, no noticeable age effect on the disc stiffness is evidenced in the in-vitro experimental studies of the literature. So, if the disc intrinsic mechanics remains constant, how to explain the correlation of disc degeneration and disc functionality loss with age.

METHODS

A microstructure-based modeling approach was developed to assess in-silico the aging-sensitive mechanics of human intervertebral disc. The model considers the relationship between stress/volumetric macro-response and microstructure along with effective age effects acting at the lamellar and multi-lamellar scales. The stress-stretch and transversal responses of the different disc regions were computed for various age groups (13-18, 36, 58, 69 and 82 years old) and their evolution with age was studied.

RESULTS

While matching with in-vitro experimental data, the predicted stiffness was found to increase while passing from adolescent young discs to mature older discs and then to remain almost constant for the rest of life. Important age-related changes in the disc transversal behavior were also predicted affecting the flexibility of the disc, changing its volumetric behavior, and modifying its dimensions.

CONCLUSION

The developed approach was found able to bring new conclusions about age-dependent mechanical properties including regional dependency. The disc mechanics in terms of rigidity, radial and axial transversal responses were found to alter going from adolescent to middle age where the disc reaches a certain maturity. After reaching maturity, the mechanical properties undergo very slight changes until becoming almost constant with age.

Using inverse finite element analysis to identify spinal tissue behaviour in situ

In computational modelling of musculoskeletal applications, one of the critical aspects is ensuring that a model can capture intrinsic population variability and not only representative of a "mean" individual. Developing and calibrating models with this aspect in mind is key for the credibility of a modelling methodology. This often requires calibration of complex models with respect to 3D experiments and measurements on a range of specimens or patients. Most Finite Element (FE) software's do not have such a capacity embedded in their core tools. This paper presents a versatile interface between Finite Element (FE) software and optimisation tools, enabling calibration of a group of FE models on a range of experimental data. It is provided as a Python toolbox which has been fully tested and verified on Windows platforms. The toolbox is tested in three case studies involving in vitro testing of spinal tissues.

Interlamellar matrix governs human annulus fibrosus multiaxial behavior

Establishing accurate structure–property relationships for intervertebral disc annulus fibrosus tissue is a fundamental task for a reliable computer simulation of the human spine but needs excessive theoretical-numerical-experimental works. The difficulty emanates from multiaxiality and anisotropy of the tissue response along with regional dependency of a complex hierarchic structure interacting with the surrounding environment. We present a new and simple hybrid microstructure-based experimental/modeling strategy allowing adaptation of animal disc model to human one. The trans-species strategy requires solely the basic knowledge of the uniaxial circumferential response of two different animal disc regions to predict the multiaxial response of any human disc region. This work demonstrates for the first time the determining role of the interlamellar matrix connecting the fibers-reinforced lamellae in the disc multiaxial response. Our approach shows encouraging multiaxial predictive capabilities making it a promising tool for human spine long-term prediction.

Estimating tissue-level properties of porcine talar subchondral bone

Tissue-level properties of bone play an important role when characterising apparent-level bone biomechanical behaviour and yet little is known about its effect at this hierarchical level. In combination with trabecular morphological data these properties can be used to predict bone strength, which becomes an invaluable tool for clinicians in patient treatment planning. This study developed specimen-specific micro-finite element (μFE) models using validated continuum-level models, containing grayscale-derived material properties, to indirectly establish tissue-level properties of porcine talar subchondral bone. Specimen-specific continuum finite element (hFE) models of subchondral trabecular bone were setup using μCT data of ten cylindrical specimens extracted from juvenile porcine tali. The models were validated using quasi-static uniaxial compression testing. Validated hFE models were used to calibrate the tissue modulus of corresponding μFE models by minimising the difference between the μFE and hFE stiffness values. Key trabecular morphological indices (BV/TV, DA, Conn.D, Tb.Th, EF) were evaluated. Good agreement was observed between hFE models and experiment (CCC = 0.66). Calibrated E_tiss was 504 ± 37.65 MPa. Average BV/TV and DA for μFE specimens were 0.37 ± 0.05 and 0.68 ± 0.11, respectively. BV/TV (r² = 0.667) correlated highly with μFE stiffness. The small intra-specimen variation to tissue-level properties suggests that variations to apparent-level stiffness originate from variations to microarchitecture rather than tissue mechanical properties.

Osmo-inelastic response of the intervertebral disc annulus fibrosus tissue

The aim of this article is to provide some insights on the osmo-inelastic response under stretching of annulus fibrosus of the intervertebral disc. Circumferentially oriented specimens of square cross section, extracted from different regions of bovine cervical discs (ventral-lateral and dorsal-lateral), are tested under different strain-rates and saline concentrations within normal range of strains. An accurate optical strain measuring technique, based upon digital image correlation, is used in order to determine the full-field displacements in the lamellae and fibers planes of the layered soft tissue. Annulus stress-stretch relationships are measured along with full-field transversal strains in the two planes. The mechanical response is found hysteretic, rate-dependent and osmolarity-dependent with a Poisson's ratio higher than 0.5 in the fibers plane and negative (auxeticity) in the lamellae plane. While the stiffness presents a regional-dependency due to variations in collagen fibers content/orientation, the strain-rate sensitivity of the response is found independent on the region. A significant osmotic effect is found on both the auxetic response in the lamellae plane and the stiffness rate-sensitivity. These local experimental observations will result in more accurate chemo-mechanical modeling of the disc annulus and a clearer multi-scale understanding of the disc intervertebral function.

Interlamellar-induced time-dependent response of intervertebral disc annulus: A microstructure-based chemo-viscoelastic model

The annulus fibrosus of the intervertebral disc exhibits an unusual transversal behavior for which a constitutive representation that considers as well regional effect, chemical sensitivity and time-dependency has not yet been developed, and it is hence the aim of the present contribution. A physically-based model is proposed by introducing a free energy function that takes into account the actual disc annulus structure in relation with the surrounding biochemical environment. The response is assumed to be dominated by the viscoelastic contribution of the extracellular matrix, the elastic contribution of the oriented collagen fibers and the osmo-induced volumetric contribution of the internal fluid content variation. The regional dependence of the disc annulus response due to variation in fibers content/orientation allows a micromechanical treatment of the soft tissue. A finite element model of the annulus specimen is designed while taking into consideration the 'interlamellar' ground substance zone between lamellae of the layered soft tissue. The kinetics is designed using full-field strain measurements performed on specimens extracted from two disc annulus regions and tested under different osmotic conditions. The time-dependency of the tissue response is reported on stress-free volumetric changes, on hysteretic stress and transversal strains during quasi-static stretching at different strain-rates and on their temporal changes during an interrupted stretching. Considering the effective contributions of the internal fluid transfer and the extracellular matrix viscosity, the microstructure-based chemo-mechanical model is found able to successfully reproduce the significant features of the macro-response and the unusual transversal behavior including the strong regional dependency from inner to outer parts of the disc: Poisson's ratio lesser than 0 (auxetic) in lamellae plane, higher than 0.5 in fibers plane, and their temporal changes towards usual values (between 0 and 0.5) at chemo-mechanical equilibrium. The underlying time-dependent mechanisms occurring in the tissue are analyzed via the local numerical fields and important insights about the effective role of the interlamellar zone are revealed for the different disc localizations. STATEMENT OF SIGNIFICANCE: The structural complexity of the annulus fibrosus has only been appreciated through recent experimental contributions and a constitutive representation that considers as well regional effect, chemical sensitivity and time-dependency of the unusual transversal behavior has not yet been developed. Here, a microstructure-based chemo-viscoelastic model is developed to highlight the interlamellar-induced time-dependent response by means of a two-scale strategy. The model provides important insights about the origin of the time-dependent phenomena in disc annulus along with regional dependency, essential for understanding disc functionality.

Microstructural characterization of annulus fibrosus by ultrasonography: a feasibility study with an in vivo and in vitro approach

The main function of the intervertebral disc is biomechanical function, since it must resist repetitive high loadings, while giving the spine its flexibility and protecting the spinal cord from over-straining. It partially owes its mechanical characteristics to the lamellar architecture of its outer layer, the annulus fibrosus. Today, no non-invasive means exist to characterize annulus lamellar structure in vivo. The aim of this work was to test the feasibility of imaging annulus fibrosus microstructure in vivo with ultrasonography. Twenty-nine healthy adolescents were included. Ultrasonographies of L3-L4 disc were acquired with a frontal approach. Annulus fibrosus was segmented in the images to measure the thickness of the lamellae. To validate lamellar appearance in ultrasonographies, multimodality images of two cow tail discs were compared: ultrasonography, magnetic resonance and optical microscopy. In vivo average lamellar thickness was 229.7 ± 91.5 μm, and it correlated with patient body mass index and age. Lamellar appearance in the three imaging modalities in vitro was consistent. Lamellar measurement uncertainty was 7%, with good agreement between two operators. Feasibility of ultrasonography for the analysis of lumbar annulus fibrosus structure was confirmed. Further work should aim at validating measurement reliability, and to assess the relevance of the method to characterize annulus alterations, for instance in disc degeneration or scoliosis.

Assessment of intervertebral disc degeneration-related properties using finite element models based on ρ H -weighted MRI data

Quantitative magnetic resonance imaging (MRI) provides useful information about intervertebral disc (IVD) biomechanical properties, especially those in relation to the fluid phase. These properties may improve IVD finite element (FE) models using data closer to physiological reality. The aim of this study is to investigate IVD degeneration-related properties using a coupling between MRI and FE modeling. To this end, proton density (ρ H )-weighted MRI sequences of a porcine lumbar IVD were carried out to develop two biphasic swelling models with hyperelastic extracellular matrix behavior. The first model is isotropic, and the second one is anisotropic and takes into account the role of collagen fibers in the mechanical behavior of the IVD. MRI sequences permitted to determine the geometry and the real porosity mapping within the disc. The differentiation between disc components (nucleus pulposus, annulus fibrosus and cartilaginous end plates) was taken into account using spatial continuous distributions of the mechanical properties. The validation of the FE models was performed through two steps: the identification of the model's mechanical properties using relaxation compressive test and the comparison between the MRI after load porosity distributions and those numerically obtained using the set of identified properties. The results confirmed that the two developed FE models were able to predict the mechanical response of uniaxial time-dependent compressive test and the redistribution of porosity after load. A slight difference between the measured and the numerical local bulges of the disc was found. This study suggests that from the coupling between MRI imaging in different state of load and finite element modeling we can deduce relevant information that can be used in the assessment of the early intervertebral disc degeneration changes.

GEORG SCHMORL PRIZE OF THE GERMAN SPINE SOCIETY (DWG) 2018: combined inflammatory and mechanical stress weakens the annulus fibrosus: evidences from a loaded bovine AF organ culture

PURPOSE

The pathomechanism of annulus fibrosus (AF) failure is still unknown. We hypothesise that mechanical overload and an inflammatory microenvironment contribute to AF structural weakening. Therefore, the objective of this study was to investigate the influence of these factors on the AF, particularly the translamellar bridging network (TLBN) which connects the AF lamellae.

METHODS

A bovine AF organ culture (AF-OC) model of standardised AF rings was used to study the individual and combined effects of cyclic tensile strain (CTS) and IL-1β (1 ng/mL) culture medium supplementation. AF-OCs were analysed for PGE₂ production (ELISA) and deposition of IL-6, COX-2, fibrillin, and MMP3 in the tissue (immunohistochemistry, IHC). The mechanical strength of the TLBN was evaluated using a peel test to measure the strength required to separate an AF segment along a lamellar bound.

RESULTS

The combination of CTS + IL-1β led to a significant increase in PGE₂ production compared to Control (p < 0.01). IHC evaluations showed that the CTS + IL-1β group exhibited higher production of COX-2 and MMP3 within the TLBN regions compared to the adjacent lamellae and a significant increase in IL-6 ratio compared to Control (p < 0.05). A significant decrease in the annular peel strength was observed in the CTS + IL1β group compared to Control (p < 0.05).

CONCLUSION

Our findings suggest that CTS and IL-1β act synergistically to increase pro-inflammatory and catabolic molecules within the AF, particularly the TLBN, leading to a weakening of the tissue. This standardised model enables the investigation of AF/TLBN structure-function relationship and is a platform to test AF-focused therapeutics. These slides can be retrieved under Electronic Supplementary Material.

The Mechanical Role of the Radial Fiber Network Within the Annulus Fibrosus of the Lumbar Intervertebral Disc: A Finite Elements Study

The annulus fibrosus (AF) of the intervertebral disc (IVD) consists of a set of concentric layers composed of a primary circumferential collagen fibers arranged in an alternating oblique orientation. Moreover, there exists an additional secondary set of radial translamellar collagen fibers which connects the concentric layers, creating an interconnected fiber network. The aim of this study was to investigate the mechanical role of the radial fiber network. Toward that goal, a three-dimensional (3D) finite element model of the L3–L4 spinal segment was generated and calibrated to axial compression and pure moment loading. The AF model explicitly recognizes the two heterogeneous networks of fibers. The presence of radial fibers demonstrated a pronounced effect on the local disc responses under lateral bending, flexion, and extension modes. In these modes, the radial fibers were in a tensile state in the disc region that subjected to compression. In addition, the circumferential fibers, on the opposite side of the IVD, were also under tension. The local stress in the matrix was decreased in up to 9% in the radial fibers presence. This implies an active fiber network acting collectively to reduce the stresses and strains in the AF lamellae. Moreover, a reduction of 26.6% in the matrix sideways expansion was seen in the presence of the radial fibers near the neutral bending axis of the disc. The proposed biomechanical model provided a new insight into the mechanical role of the radial collagen fibers in the AF structure. This model can assist in the design of future IVD substitutes.

Biomaterials for intervertebral disc regeneration: Current status and looming challenges

A biomaterial-based strategy is employed to regenerate the degenerated intervertebral disc, which is considered a major generator of neck and back pain. Although encouraging enhancements in the anatomy and kinematics of the degenerative disc have been gained by biomaterials with various formulations in animals, the number of biomaterials tested in humans is rare. At present, most studies that involve the use of newly developed biomaterials focus on regeneration of the degenerative disc, but not pain relief. In this review, we summarise the current state of the art in the field of biomaterial-based regeneration or repair for the nucleus pulposus, annulus fibrosus, and total disc transplantation in animals and humans, and we then provide essential suggestions for the development and clinical translation of biomaterials for disc regeneration. It is important for researchers to consider the commonly neglected issues instead of concentrating solely on biomaterial development and fabrication.

New insights into the viscoelastic and failure mechanical properties of the elastic fiber network of the inter-lamellar matrix in the annulus fibrosus of the disc

The mechanical role of elastic fibers in the inter-lamellar matrix (ILM) is unknown; however, it has been suggested that they play a role in providing structural integrity to the annulus fibrosus (AF). Therefore, the aim of this study was to measure the viscoelastic and failure properties of the elastic fiber network in the ILM of ovine discs under both tension and shear directions of loading. Utilizing a technique, isolated elastic fibers within the ILM from ovine discs 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), followed by a ramp test to failure at 10% s^-1. A significant strain-rate dependent response was found, particularly at the fastest rate for phase angle and normalized stiffness (p < 0.001). The elastic fibers in the ILM demonstrated a significantly higher capability for energy absorption at slow compared to medium and fast strain rates (p < 0.001). These finding suggests that the elastic fiber network of the ILM exhibits nonlinear elastic behavior. When tested to failure, a significantly higher normalized failure force was found in tension compared to shear loading (p = 0.011), which is consistent with the orthotropic structure of elastic fibers in the ILM. The results of this study confirmed the mechanical contribution of the elastic fiber network to the ILM and the structural integrity of the AF. This research serves as a foundation for future studies to investigate the relationship between degeneration and ILM mechanical properties.

STATEMENT OF SIGNIFICANCE

The mechanical role of elastic fibres in the inter-lamellar matrix (ILM) of the disc is unknown. The viscoelastic and failure properties of the elastic fibre network in the ILM in both tension and shear directions of loading was measured for the first time. We found a strain-rate dependent response for the elastic fibres in the ILM. The elastic fibres in the ILM demonstrated a significantly higher capability for energy absorption at slow compared to medium and fast strain rates. When tested to failure, a significantly higher normalized failure force was found in tension compared to shear loading, which is consistent with the orthotropic structure of elastic fibres in the ILM.

Towards engineering distinct multi-lamellated outer and inner annulus fibrosus tissues

The annulus fibrosus (AF) of the intervertebral disc (IVD) has a zonal distribution of phenotypically distinct cells. The outer AF (OAF) cells produce an extracellular matrix (ECM) rich in type I collagen with little proteoglycans, whereas the ECM of the inner AF (IAF) has abundant type II collagen and proteoglycans. The inhomogeneous distribution of the ECM in the AF may reflect the complex mechanical forces that the IVD experiences. A bioengineered AF tissue should recapitulate both the inner and outer zones in order to have proper functionality. The aim of this study is to generate multi-lamellated OAF and IAF tissues with ECM compositions that resemble their zonal origin using polycarbonate urethane (PU) scaffolds. It was observed that supplementation of the media with insulin-transferrin-selenium (ITS) and proline yielded tissues with good cellularity. However, IAF cells accumulated only type I collagen, similar to OAF cells. Addition of dexamethasone and sodium pyruvate induced the accumulation of IAF tissues rich in type II collagen and aggrecan, without altering the accumulation of type I collagen in OAF tissues. Dexamethasone stimulated mitochondrial membrane potential in both tissues in the presence of sodium pyruvate, suggesting a relationship between the mitochondrial aerobic respiratory state and dexamethasone signalling during the in vitro-tissue formation by OAF and IAF cells. Inhibition of the glucocorticoid receptor blocked the stimulation of mitochondrial membrane potentials and type II collagen accumulation. In summary, biologically distinct multi-lamellated OAF and IAF tissues can be generated, which will facilitate advancement towards the goal of engineering a biological IVD replacement. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1346-1355, 2018.

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.

Ultrastructural organization of elastic fibres in the partition boundaries of the annulus fibrosus within the intervertebral disc

The relationship between elastic fibre disorders and disc degeneration, aging and progression of spine deformity have been discussed in a small number of studies. However, the clinical relevance of elastic fibres in the annulus fibrosus (AF) of the disc is poorly understood. Ultrastructural visualization of elastic fibres is an important step towards understanding their structure-function relationship. In our previous studies, a novel technique for visualization of elastic fibres across the AF was presented and their ultrastructural organization in intra- and inter-lamellar regions was compared. Using the same novel technique in the present study, the ultrastructural organization of elastic fibres in the partition boundaries (PBs), which are located between adjacent collagen bundles, is presented for the first time. Visualization of elastic fibres in the PBs in control and partially digested (digested) samples was compared, and their orientation in two different cutting planes (transverse and oblique) were discussed. The ultrastructural analysis revealed that elastic fibres in PBs were a well-organized dense and complex network having different size and shape. Adjacent collagen bundles in a cross section (CS) lamella appear to be connected to each other, where elastic fibres in the PBs were merged in parallel or penetrated into the collagen bundles. There was no significant difference in directional coherency coefficient of elastic fibres between the two different cutting planes (p = .35). The present study revealed that a continuous network of elastic fibres may provide disc integrity by connecting adjacent bundles of CS lamellae together. Compared to our previous studies, the density of the elastic fibre network in PBs was lower, and fibre orientation was similar to the intra-lamellar space and inter-lamellar matrix.

STATEMENT OF SIGNIFICANCE

A detailed ultrastructural study in the partition boundaries of the annulus fibrosus within the disc revealed a well-organized elastic fibre network with a complex ultrastructure. The continuous network of elastic fibres may provide disc integrity by connecting adjacent bundles of cross section lamellae together. The density of the elastic fibre network in PBs was lower, and fibre orientation was similar to the intra-lamellar space and the inter-lamellar matrix.

A study of the damage behaviour of porcine intervertebral discs in a bioreactor environment

The intervertebral discs are cartilaginous, articulating structures that lie between vertebral bodies, allowing flexibility, transmission, modification, and also distribution of the forces to the spinal column. Disc degeneration is characterised by progressive loss of disc height and exaggerated radial bulging. Therefore, the spine becomes shorter, stiffer, and less mobile. In the last several decades, there is a strong need for a tissue engineering strategy that alleviates pain and restores spine function by directly addressing the underlying biological causes of disc degeneration. Numerous studies that are currently showing potential have been conducted on developing regenerative and reparative strategies for treating this condition. In this study, to numerically describe the anisotropic mechanical damage behaviour of discs, the pseudo-elastic damage model was applied. To experimentally picture the biomechanical response of discs and to study the damage mechanisms as well as the spinal disc herniation, a special bioreactor was evolved. The specimens were obtained from pigs aged six months. A total of eight functional spine units were taken from porcine lumbar spines (L1-L2). Firstly, the experiments were performed by using long-term cyclic uniaxial compression tests. Secondly, the mean value of experimental results with consideration of the different shapes and sizes of the samples was calculated. Afterwards, the experimental results were compared with outcomes of numerical simulations.

The ultra-structural organization of the elastic network in the intra- and inter-lamellar matrix of the intervertebral disc

The inter-lamellar matrix (ILM)-located between adjacent lamellae of the annulus fibrosus-consists of a complex structure of elastic fibers, while elastic fibers of the intra-lamellar region are aligned predominantly parallel to the collagen fibers. The organization of elastic fibers under low magnification, in both inter- and intra-lamellar regions, was studied by light microscopic analysis of histologically prepared samples; however, little is known about their ultrastructure. An ultrastructural visualization of elastic fibers in the inter-lamellar matrix is crucial for describing their contribution to structural integrity, as well as mechanical properties of the annulus fibrosus. The aims of this study were twofold: first, to present an ultrastructural analysis of the elastic fiber network in the ILM and intra-lamellar region, including cross section (CS) and in-plane (IP) lamellae, of the AF using Scanning Electron Microscopy (SEM) and second, to -compare the elastic fiber orientation between the ILM and intra-lamellar region. Four samples (lumbar sheep discs) from adjacent sections (30μm thickness) of anterior annulus were partially digested by a developed NaOH-sonication method for visualization of elastic fibers by SEM. Elastic fiber orientation and distribution were quantified relative to the tangential to circumferential reference axis. Visualization of the ILM under high magnification revealed a dense network of elastic fibers that has not been previously described. Within the ILM, elastic fibers form a complex network, consisting of different size and shape fibers, which differed to those located in the intra-lamellar region. For both regions, the majority of fibers were oriented near 0° with respect to tangential to circumferential (TCD) direction and two minor symmetrical orientations of approximately±45°. Statistically, the orientation of elastic fibers between the ILM and intra-lamellar region was not different (p=0.171). The present study used extracellular matrix partial digestion to address significant gaps in understanding of disc microstructure and will contribute to multidisciplinary ultrastructure-function studies.

STATEMENT OF SIGNIFICANCE

Visualization of the intra-lamellar matrix under high magnification revealed a dense network of elastic fibers that has not been previously described. The present study used extracellular matrix partial digestion to address significant gaps in understanding of disc microstructure and will contribute to multidisciplinary ultrastructure-function studies.

Annulus fibrosus functional extrafibrillar and fibrous mechanical behaviour: experimental and computational characterisation

The development of current surgical treatments for intervertebral disc damage could benefit from virtual environment accounting for population variations. For such models to be reliable, a relevant description of the mechanical properties of the different tissues and their role in the functional mechanics of the disc is of major importance. The aims of this work were first to assess the physiological hoop strain in the annulus fibrosus in fresh conditions (n = 5) in order to extract a functional behaviour of the extrafibrillar matrix; then to reverse-engineer the annulus fibrosus fibrillar behaviour (n = 6). This was achieved by performing both direct and global controlled calibration of material parameters, accounting for the whole process of experimental design and in silico model methodology. Direct-controlled models are specimen-specific models representing controlled experimental conditions that can be replicated and directly comparing measurements. Validation was performed on another six specimens and a sensitivity study was performed. Hoop strains were measured as 17 ± 3% after 10 min relaxation and 21 ± 4% after 20-25 min relaxation, with no significant difference between the two measurements. The extrafibrillar matrix functional moduli were measured as 1.5 ± 0.7 MPa. Fibre-related material parameters showed large variability, with a variance above 0.28. Direct-controlled calibration and validation provides confidence that the model development methodology can capture the measurable variation within the population of tested specimens.

Methodology to Produce Specimen-Specific Models of Vertebrae: Application to Different Species

Image-based continuum-level finite element models have been used for bones to evaluate fracture risk and the biomechanical effects of diseases and therapies, capturing both the geometry and tissue mechanical properties. Although models of vertebrae of various species have been developed, an inter-species comparison has not yet been investigated. The purpose of this study was to derive species-specific modelling methods and compare the accuracy of image-based finite element models of vertebrae across species. Vertebral specimens were harvested from porcine (N = 12), ovine (N = 13) and bovine (N = 14) spines. The specimens were experimentally loaded to failure and apparent stiffness values were derived. Image-based finite element models were generated reproducing the experimental protocol. A linear relationship between the element grayscale and elastic modulus was calibrated for each species matching in vitro and in silico stiffness values, and validated on independent sets of models. The accuracy of these relationships were compared across species. Experimental stiffness values were significantly different across species and specimen-specific models required species-specific linear relationship between image grayscale and elastic modulus. A good agreement between in vitro and in silico values was achieved for all species, reinforcing the generality of the developed methodology.

Structure and mechanical function of the inter-lamellar matrix of the annulus fibrosus in the disc

The inter-lamellar matrix (ILM) has an average thickness of less than 30 µm and lies between adjacent lamellae in the annulus fibrosus (AF). The microstructure and composition of the ILM have been studied in various anatomic regions of the disc; however, their contribution to AF mechanical properties and structural integrity is unknown. It was suggested that the ILM components, mainly elastic fibers and cross-bridges, play a role in providing mechanical integrity of the AF. Therefore, the manner in which they respond to different loadings and stabilize adjacent lamellae structure will influence AF tear formation and subsequent herniation. This review paper summarizes the composition, microstructure, and potential role of the ILM in the progression of disc herniation, clarifies the micromechanical properties of the ILM, and proposes critical areas for future studies. There are a number of unknown characteristics of the ILM, such as its mechanical role, impact on AF integrity, and ultrastructure of elastic fibers at the ILM-lamella boundary. Determining these characteristics will provide important information for tissue engineering, repair strategies, and the development of more-physiological computational models to study the initiation and propagation of AF tears that lead to herniation and degeneration. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1307-1315, 2016.

Modelling the failure precursor mechanism of lamellar fibrous tissues, example of the annulus fibrosus

The aims of this study were to assess the damage and failure strengths of lamellar fibrous tissues, such as the anterior annulus fibrosus (AF), and to develop a mathematical model of damage propagation of the lamellae and inter-lamellar connections. This level of modelling is needed to accurately predict the effect of damage and failure induced by trauma or clinical interventions. 26 ovine anterior AF cuboid specimens from 11 lumbar intervertebral discs were tested in radial tension and mechanical parameters defining damage and failure were extracted from the in-vitro data. Equivalent 1D analytical models were developed to represent the specimen strength and the damage propagation, accounting for the specimen dimensions and number of lamellae. Model parameters were calibrated on the in-vitro data. Similar to stiffness values reported for other orientations, the outer annulus was found stronger than the inner annulus in the radial direction, with failure at higher stress values. The inner annulus failed more progressively, showing macroscopic failure at a higher strain value. The 1D analytical model of damage showed that lamellar damage is predominant in the failure mechanism of the AF. The analytical model of the connections between lamellae allowed us to represent separately damage processes in the lamellae and the inter-lamellar connections, which cannot be experimentally tested individually.

Lamellar and fibre bundle mechanics of the annulus fibrosus in bovine intervertebral disc

The intervertebral disc is a multicomposite structure, with an outer fibrous ring, the annulus fibrosus, retaining a gel-like core, the nucleus pulposus. The disc presents complex mechanical behaviour, and it is of high importance for spine biomechanics. Advances in multiscale modelling and disc repair raised a need for new quantitative data on the finest details of annulus fibrosus mechanics. In this work we explored inter-lamella and inter-bundle behaviour of the outer annulus using micromechanical testing and second harmonic generation microscopy. Twenty-one intervertebral discs were dissected from cow tails; the nucleus and inner annulus were excised to leave a ring of outer annulus, which was tested in circumferential loading while imaging the tissue's collagen fibres network with sub-micron resolution. Custom software was developed to determine local tissue strains through image analysis. Inter-bundle linear and shear strains were 5.5 and 2.8 times higher than intra-bundle strains. Bundles tended to remain parallel while rotating under loading, with large slipping between them. Inter-lamella linear strain was almost 3 times the intra-lamella one, but no slipping was observed at the junction between lamellae. This study confirms that outer annulus straining is mainly due to bundles slipping and rotating. Further development of disc multiscale modelling and repair techniques should take into account this modular behaviour of the lamella, rather than considering it as a homogeneous fibre-reinforced matrix.

STATEMENT OF SIGNIFICANCE

The intervertebral disc is an organ tucked between each couple of vertebrae in the spine. It is composed by an outer fibrous layer retaining a gel-like core. This organ undergoes severe and repeated loading during everyday life activities, since it is the compliant component that gives the spine its flexibility. Its properties are affected by pathologies such as disc degeneration, a major cause of back pain. In this article we explored the micromechanical behaviour of the disc's outer layer using second harmonic generation, a technique which allowed us to visualize, with unprecedented detail, how bundles of collagen fibres slide relative to each other when loaded. Our results will help further the development of new multiscale numerical models and repairing techniques.