Understanding the etiopathogenesis of lumbar intervertebral disc herniation: From clinical evidence to basic scientific research

Lumbar intervertebral disc herniation, as a leading cause of low back pain, productivity loss, and disability, is a common musculoskeletal disorder that results in significant socioeconomic burdens. Despite extensive clinical and basic scientific research efforts, herniation etiopathogenesis, particularly its initiation and progression, is not well understood. Understanding herniation etiopathogenesis is essential for developing effective preventive measures and therapeutic interventions. Thus, this review seeks to provide a thorough overview of the advances in herniation-oriented research, with a discussion on ongoing challenges and potential future directions for clinical, translational, and basic scientific investigations to facilitate innovative interdisciplinary research aimed at understanding herniation etiopathogenesis. Specifically, risk factors for herniation are identified and summarized, including familial predisposition, obesity, diabetes mellitus, smoking tobacco, selected cardiovascular diseases, disc degeneration, and occupational risks. Basic scientific experimental and computational research that aims to understand the link between excessive mechanical load, catabolic tissue remodeling due to inflammation or insufficient nutrient supply, and herniation, are also reviewed. Potential future directions to address the current challenges in herniation-oriented research are explored by combining known progressive development in existing research techniques with ongoing technological advances. More research on the relationship between occupational risk factors and herniation, as well as the relationship between degeneration and herniation, is needed to develop preventive measures for working-age individuals. Notably, researchers should explore using or modifying existing degeneration animal models to study herniation etiopathogenesis, as such models may allow for a better understanding of how to prevent mild-to-moderately degenerated discs from herniating.

An Approach to Quantify Anisotropic Multiaxial Failure of the Annulus Fibrosus

Tears in the annulus fibrosus (AF) of the intervertebral disk (IVD) occur due to multiaxial loading on the spine. However, most existing AF failure studies measure uniaxial stress, not the multiaxial stress at failure. Delamination theory, which requires advanced structural knowledge and knowledge about the interactions between the AF fibers and matrix, has historically been used to understand and predict AF failure. Alternatively, a simple method, the Tsai-Hill yield criteria, could describe multiaxial failure of the AF. This yield criteria uses the known tissue fiber orientation and an equation to establish the multiaxial failure stresses that cause failure. This paper presents a method to test the multiaxial failure stress of the AF experimentally and evaluate the potential for the Tsai-Hill model to predict these failure stresses. Porcine AF was cut into a dogbone shape at three distinct angles relative to the primary lamella direction (parallel, transverse, and oblique). Then, each dogbone was pulled to complete rupture. The Cauchy stress in the material's fiber coordinates was calculated. These multiaxial stress parameters were used to optimize the coefficients of the Tsai-Hill yield. The coefficients obtained for the Tsai-Hill model vary by an order of magnitude between the fiber and transverse directions, and these coefficients provide a good description of the AF multiaxial failure stress. These results establish both an experimental approach and the use of the Tsai-Hill model to explain the anisotropic failure behavior of the tissue.

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.

Comparison of biomechanical parameters of two Chinese cervical spine rotation manipulations based on motion capture and finite element analysis

Objective: The purpose of this study was to obtain the stress-strain of the cervical spine structure during the simulated manipulation of the oblique pulling manipulation and the cervical rotation-traction manipulation in order to compare the mechanical mechanism of the two manipulations. Methods: A motion capture system was used to record the key kinematic parameters of operating the two manipulations. At the same time, a three-dimensional finite element model of the C0-T1 full healthy cervical spine was established, and the key kinematic parameters were loaded onto the finite element model in steps to analyze and simulate the detailed process of the operation of the two manipulations. Results: A detailed finite element model of the whole cervical spine including spinal nerve roots was established, and the validity of this 3D finite element model was verified. During the stepwise simulation of the two cervical spine rotation manipulations to the right, the disc (including the annulus fibrosus and nucleus pulposus) and facet joints stresses and displacements were greater in the oblique pulling manipulation group than in the cervical rotation-traction manipulation group, while the spinal cord and nerve root stresses were greater in the cervical rotation-traction manipulation group than in the oblique pulling manipulation group. The spinal cord and nerve root stresses in the cervical rotation-traction manipulation group were mainly concentrated in the C4/5 and C5/6 segments. Conclusion: The oblique pulling manipulation may be more appropriate for the treatment of cervical spondylotic radiculopathy, while cervical rotation-traction manipulation is more appropriate for the treatment of cervical spondylosis of cervical type. Clinicians should select cervical rotation manipulations for different types of cervical spondylosis according to the patient's symptoms and needs.

Non-enzymatic glycation of annulus fibrosus alters tissue-level failure mechanics in tension

Advanced-glycation end products (AGEs) are known to accumulate in biological tissues with age and at an accelerated rate in patients with diabetes and chronic kidney disease. Clinically, diabetes has been linked to increased frequency and severity of back pain, accelerated disc degeneration, and an increased risk of disc herniation. Despite significant clinical evidence suggesting that diabetes-induced AGEs may play a role in intervertebral disc failure and substantial previous work investigating the effects of AGEs on bone, cartilage, and tendon mechanics, the effects of AGEs on annulus fibrosus (AF) failure mechanics have not yet been reported. Thus, the aim of this study was to determine the relationship between physiological levels of AGEs and AF tensile mechanics at two distinct loading rates. In vitro glycation treatments with methylglyoxal were applied to minimize changes in tissue hydration and induce two distinct levels of AGEs based on values measured from human AF tissues. In vitro glycation increased modulus by 48-99% and failure stress by 45-104% versus control and decreased post-failure energy absorption capacity by 15-32% versus control (ANOVA p < 0.0001 on means; range given across two loading rates and glycation levels). AGE content correlated strongly with modulus (R = 0.74, p < 0.0001) and failure stress (R = 0.70, p < 0.0001) and moderately with post-failure energy absorption capacity (R = 0.62, p < 0.0001). Failure strain was reduced by 10-17% at the high-glycation level (ANOVA p = 0.01). Tissue water content remained near or just above fresh-tissue levels for all groups. The alterations in mechanics with glycation reported here are consistent with trends from other connective tissues but do not fully explain the clinical predisposition of diabetics to disc herniation. The results from this study may be used in the development of advanced computational models that aim to study disc disease progression and to provide a deeper understanding of altered structure-function relationships that may lead to tissue dysfunction and failure with aging and disease.

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.

A Robust Multiscale and Multiphasic Structure-Based Modeling Framework for the Intervertebral Disc

A comprehensive understanding of multiscale and multiphasic intervertebral disc mechanics is crucial for designing advanced tissue engineered structures aiming to recapitulate native tissue behavior. The bovine caudal disc is a commonly used human disc analog due to its availability, large disc height and area, and similarities in biochemical and mechanical properties to the human disc. Because of challenges in directly measuring subtissue-level mechanics, such as in situ fiber mechanics, finite element models have been widely employed in spinal biomechanics research. However, many previous models use homogenization theory and describe each model element as a homogenized combination of fibers and the extrafibrillar matrix while ignoring the role of water content or osmotic behavior. Thus, these models are limited in their ability in investigating subtissue-level mechanics and stress-bearing mechanisms through fluid pressure. The objective of this study was to develop and validate a structure-based bovine caudal disc model, and to evaluate multiscale and multiphasic intervertebral disc mechanics under different loading conditions and with degeneration. The structure-based model was developed based on native disc structure, where fibers and matrix in the annulus fibrosus were described as distinct materials occupying separate volumes. Model parameters were directly obtained from experimental studies without calibration. Under the multiscale validation framework, the model was validated across the joint-, tissue-, and subtissue-levels. Our model accurately predicted multiscale disc responses for 15 of 16 cases, emphasizing the accuracy of the model, as well as the effectiveness and robustness of the multiscale structure-based modeling-validation framework. The model also demonstrated the rim as a weak link for disc failure, highlighting the importance of keeping the cartilage endplate intact when evaluating disc failure mechanisms in vitro. Importantly, results from this study elucidated important fluid-based load-bearing mechanisms and fiber-matrix interactions that are important for understanding disease progression and regeneration in intervertebral discs. In conclusion, the methods presented in this study can be used in conjunction with experimental work to simultaneously investigate disc joint-, tissue-, and subtissue-level mechanics with degeneration, disease, and injury.

Influence of crack length and anatomical location on the fracture toughness of annulus fibrosus

Fracture toughness (J_c) of a soft biological tissue is an important mechanical property that characterizes its resistance to crack or tear extension. To date, no information is available on fracture toughness of annulus fibrosus (AF); therefore, its defect tolerance is not known. The present study modified a previously introduced method to determine J_c of ovine AF. Then, the effect of the notch length on the failure pattern and J_c was investigated. Also, the test samples of anterior and lateral regions were collected to determine the effect of the location on J_c. Results showed that for a notch length of less than 45% of total width, no crack extension occurred, but for a notch length above 45% of the width, crack propagation and ultimately the failure of the AF were observed. However, statistical analysis indicated no significant difference on J_c (p = 0.5) for the initial notch length of 50% and 70% of total width. The fracture toughness was significantly higher for the samples extracted from the lateral site than those from the anterior site (p < 0.05). Dissimilar failure patterns were observed for different initial notch lengths. Among the parameters studied, the defect tolerance of AF was dependent on the initial tear size.

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.

Historical Review of Combined Experimental and Computational Approaches for Investigating Annulus Fibrosus Mechanics

Intervertebral disc research has sought to develop a deeper understanding of spine biomechanics, the complex relationship between disc health and back pain, and the mechanisms of spinal injury and repair. To do so, many researchers have focused on characterizing tissue-level properties of the disc, where the roles of tissue subcomponents can be more systematically investigated. Unfortunately, experimental challenges often limit the ability to measure important disc tissue- and subtissue-level behaviors, including fiber-matrix interactions, transient nutrient and electrolyte transport, and damage propagation. Numerous theoretical and numerical modeling frameworks have been introduced to explain, complement, guide, and optimize experimental research efforts. The synergy of experimental and computational work has significantly advanced the field, and these two aspects have continued to develop independently and jointly. Meanwhile, the relationship between experimental and computational work has become increasingly complex and interdependent. This has made it difficult to interpret and compare results between experimental and computational studies, as well as between solely computational studies. This paper seeks to explore issues of model translatability, robustness, and efficient study design, and to propose and motivate potential future directions for experimental, computational, and combined tissue-level investigations of the intervertebral disc.

Comparison of Two Intervertebral Disc Failure Models in a Numerical C4-C5 Trauma Model

The intervertebral disc (IVD) is essential for the mobility and stability of the spine. During flexion-distraction injuries, which are frequent at the cervical spine level, the IVD is often disrupted. Finite element studies have been done to investigate injury mechanisms and patterns at the cervical spine. However, they rarely include IVD failure model. The aim of this paper was to implement and compare two types of IVD failure models and their impact on hyperflexion and hyperflexion-compression injuries simulations. The failure models were tested on a detailed C4-C5 finite elements model. The first failure model consisted in a maximal strain model applied to the elements of the annulus and nucleus. The second failure model consisted in the implementation of a rupture plane in the middle of the IVD with a tied interface created between the two sections. This interface is defined by threshold stress values of detachment in traction and shearing. The two failure models were tested in flexion only and in flexion-compression. The model without inclusion of an IVD failure model was also tested. Loads at failure and injury patterns were reported. Both failure models produce failure loads that were consistent with experimental data. Injury patterns observed were in agreement with experimental and numerical studies. However, in flexion-compression, the rupture plane model simulation reached important energy error due to high deformations in the IVD elements. Also, without inclusion of an IVD failure model, energy error forced the end of the simulation in flexion-compression. Therefore, inclusion of IVD failure model is important since it leads to realistic results, but the maximal strain failure model is recommended.

Multiscale composite model of fiber-reinforced tissues with direct representation of sub-tissue properties

In many fiber-reinforced tissues, collagen fibers are embedded within a glycosaminoglycan-rich extrafibrillar matrix. Knowledge of the structure-function relationship between the sub-tissue properties and bulk tissue mechanics is important for understanding tissue failure mechanics and developing biological repair strategies. Difficulties in directly measuring sub-tissue properties led to a growing interest in employing finite element modeling approaches. However, most models are homogeneous and are therefore not sufficient for investigating multiscale tissue mechanics, such as stress distributions between sub-tissue structures. To address this limitation, we developed a structure-based model informed by the native annulus fibrosus structure, where fibers and the matrix were described as distinct materials occupying separate volumes. A multiscale framework was applied such that the model was calibrated at the sub-tissue scale using single-lamellar uniaxial mechanical test data, while validated at the bulk scale by predicting tissue multiaxial mechanics for uniaxial tension, biaxial tension, and simple shear (13 cases). Structure-based model validation results were compared to experimental observations and homogeneous models. While homogeneous models only accurately predicted bulk tissue mechanics for one case, structure-based models accurately predicted bulk tissue mechanics for 12 of 13 cases, demonstrating accuracy and robustness. Additionally, six of eight structure-based model parameters were directly linked to tissue physical properties, further broadening its future applicability. In conclusion, the structure-based model provides a powerful multiscale modeling approach for simultaneously investigating the structure-function relationship at the sub-tissue and bulk tissue scale, which is important for studying multiscale tissue mechanics with degeneration, disease, or injury.

Mechanics of pulmonary airways: Linking structure to function through constitutive modeling, biochemistry, and histology

Breathing involves fluid-solid interactions in the lung; however, the lack of experimental data inhibits combining the mechanics of air flow to airway deformation, challenging the understanding of how biomaterial constituents contribute to tissue response. As such, lung mechanics research is increasingly focused on exploring the relationship between structure and function. To address these needs, we characterize mechanical properties of porcine airways using uniaxial tensile experiments, accounting for bronchial orientation- and location- dependency. Structurally-reinforced constitutive models are developed to incorporate the role of collagen and elastin fibers embedded within the extrafibrillar matrix. The strain-energy function combines a matrix description (evaluating six models: compressible NeoHookean, unconstrained Ogden, uncoupled Mooney-Rivlin, incompressible Ogden, incompressible Demiray and incompressible NeoHookean), superimposed with non-linear fibers (evaluating two models: exponential and polynomial). The best constitutive formulation representative of all bronchial regions is determined based on curve-fit results to experimental data, accounting for uniqueness and sensitivity. Glycosaminoglycan and collagen composition, alongside tissue architecture, indicate fiber form to be primarily responsible for observed airway anisotropy and heterogeneous mechanical behavior. To the authors' best knowledge, this study is the first to formulate a structurally-motivated constitutive model, augmented with biochemical analysis and microstructural observations, to investigate the mechanical function of proximal and distal bronchi. Our systematic pulmonary tissue characterization provides a necessary foundation for understanding pulmonary mechanics; furthermore, these results enable clinical translation through simulations of airway obstruction in disease, fluid-structure interaction insights during breathing, and potentially, predictive capabilities for medical interventions. STATEMENT OF SIGNIFICANCE: The advancement of pulmonary research relies on investigating the biomechanical response of the bronchial tree. Experiments demonstrating the non-linear, heterogeneous, and anisotropic material behavior of porcine airways are used to develop a structural constitutive model representative of proximal and distal bronchial behavior. Calibrated material parameters exhibit regional variation in biomaterial properties, initially hypothesized to originate from tissue constituents. Further exploration through biochemical and histological analysis indicates mechanical function is primarily governed by microstructural form. The results of this study can be directly used in finite element and fluid-structure interaction models to enable physiologically relevant and more accurate computational simulations aimed to help diagnose and monitor pulmonary disease.