Intervertebral disc swelling maintains strain homeostasis throughout the annulus fibrosus: A finite element analysis of healthy and degenerated discs

Geometric determinants of the mechanical behavior of image-based finite element models of the intervertebral disc

The intervertebral disc is an important structure for load transfer through the spine. Its injury and degeneration have been linked to pain and spinal fractures. Disc injury and spine fractures are associated with high stresses; however, these stresses cannot be measured, necessitating the use of finite element (FE) models. These models should include the disc's complex structure, as changes in disc geometry have been linked to altered mechanical behavior. However, image-based models using disc-specific structures have yet to be established. This study describes a multiphasic FE modeling approach for noninvasive estimates of subject-specific intervertebral disc mechanical behavior based on medical imaging. The models (n = 22) were used to study the influence of disc geometry on the predicted global mechanical response (moments and forces), internal local disc stresses, and tractions at the interface between the disc and the bone. Disc geometry was found to have a strong influence on the predicted moments and forces on the disc (R2 = 0.69-0.93), while assumptions regarding the side curvature (bulge) of the disc had only a minor effect. Strong variability in the predicted internal disc stresses and tractions was observed between the models (mean absolute differences of 5.1%-27.7%). Disc height had a systematic influence on the internal disc stresses and tractions at the disc-to-bone interface. The influence of disc geometry on mechanics highlights the importance of disc-specific modeling to estimate disc injury risk, loading on the adjacent vertebral bodies, and the mechanical environment present in disc tissues.

Research on the role and mechanism of IL-17 in intervertebral disc degeneration

Intervertebral disc degeneration (IDD) is one of the primary causes of low back pain (LBP), which seriously affects patients' quality of life. In recent years, interleukin (IL)-17 has been shown to be highly expressed in the intervertebral disc (IVD) tissues and serum of patients with IDD, and IL-17A has been shown to promote IDD through multiple pathways. We first searched databases such as PubMed, Cochrane, Embase, and Web of Science using the search terms "IL-17 or interleukin 17″ and "intervertebral discs". The search period ranged from the inception of the databases to December 2023. A total of 24 articles were selected after full-text screening. The main conclusion of the clinical studies was that IL-17A levels are significantly increased in the IVD tissues and serum of IDD patients. The results from the in vitro studies indicated that IL-17A can activate signaling pathways such as the NF-κB and MAPK pathways; promote inflammatory responses, extracellular matrix degradation, and angiogenesis; and inhibit autophagy in nucleus pulposus cells. The main finding of the in vivo experiments was that puncture of animal IVDs resulted in elevated levels of IL-17A within the IVD, thereby inducing IDD. Clinical studies, in vitro experiments, and in vivo experiments confirmed that IL-17A is closely related to IDD. Therefore, drugs that target IL-17A may be novel treatments for IDD, providing a new theoretical basis for IDD therapy.

Biomechanical Analysis of a Newly Proposed Surgical Combination (MIS Screw-Rod System for Indirect Decompression+ Interspinous Fusion System for long Term Spinal Stability) in Treatment of Lumbar Degenerative Diseases

OBJECTIVE

The aim of this study is to analyze the biomechanical stability of a newly proposed surgical combination (minimally invasive surgery of screw-rod system for indirect decompression + interspinous fusion system for long term spinal stability) in treatment of lumbar degenerative diseases.

METHODS

The three-dimensional (3D) computed tomography (CT) image data of an adult healthy male volunteer were selected. An intact model of L4/5 was further established and validated by using Mimic and 3-matic, 3D slicer, abaqus, Python. Four surgical models were constructed. The biomechanical stability among these surgical modes was compared and analyzed using finite element analysis.

RESULTS

The maximum von mises on fixation system in surgical models 2 and 3 exhibited comparable values. This finding suggested that the increase in interspinous fusion did not result in a significant elevation in maximum von mises on fixation system. Compared with the third surgical model, the fourth model, which received less average von mises experienced by the screw in contact with both cancellous and cortical bone. The findings indicated that the inclusion of facet joint fusion in surgical procedures might not be necessary to increase the average von Mises stress experienced by the screw in contact with both cancellous and cortical bone.

CONCLUSIONS

The biomechanical stability of the newly proposed surgical combination (MIS screw-rod for indirect decompression + interspinous fusion for long term spinal stability technique) was not lower than that of the other surgical combination groups, and it might not be necessary to perform facet joint fusion during the surgery.

Stress-Activated Protein Kinases in Intervertebral Disc Degeneration: Unraveling the Impact of JNK and p38 MAPK

Intervertebral disc degeneration (IDD) is a major cause of lower back pain. The pathophysiological development of IDD is closely related to the stimulation of various stressors, including proinflammatory cytokines, abnormal mechanical stress, oxidative stress, metabolic abnormalities, and DNA damage, among others. These factors prevent normal intervertebral disc (IVD) development, reduce the number of IVD cells, and induce senescence and apoptosis. Stress-activated protein kinases (SAPKs), particularly, c-Jun N-terminal kinase (JNK) and p38 mitogen-activated protein kinase (p38 MAPK), control cell signaling in response to cellular stress. Previous studies have shown that these proteins are highly expressed in degenerated IVD tissues and are involved in complex biological signal-regulated processes. Therefore, we summarize the research reports on IDD related to JNK and p38 MAPK. Their structure, function, and signal regulation mechanisms are comprehensively and systematically described and potential therapeutic targets are proposed. This work could provide a reference for future research and help improve molecular therapeutic strategies for IDD.

From hyperglycemia to intervertebral disc damage: exploring diabetic-induced disc degeneration

The incidence of lumbar disc herniation has gradually increased in recent years, and most patients have symptoms of low back pain and nerve compression, which brings a heavy burden to patients and society alike. Although the causes of disc herniation are complex, intervertebral disc degeneration (IDD) is considered to be the most common factor. The intervertebral disc (IVD) is composed of the upper and lower cartilage endplates, nucleus pulposus, and annulus fibrosus. Aging, abnormal mechanical stress load, and metabolic disorders can exacerbate the progression of IDD. Among them, high glucose and high-fat diets (HFD) can lead to fat accumulation, abnormal glucose metabolism, and inflammation, which are considered important factors affecting the homeostasis of IDD. Diabetes and advanced glycation end products (AGEs) accumulation- can lead to various adverse effects on the IVD, including cell senescence, apoptosis, pyroptosis, proliferation, and Extracellular matrix (ECM) degradation. While current research provides a fundamental basis for the treatment of high glucose-induced IDD patients. further exploration into the mechanisms of abnormal glucose metabolism affecting IDD and in the development of targeted drugs will provide the foundation for the effective treatment of these patients. We aimed to systematically review studies regarding the effects of hyperglycemia on the progress of IDD.

Comprehensive modeling of annulus fibrosus: From biphasic refined characterization to damage accumulation under viscous loading

The annulus fibrosus (AF), a permeable, hydrated, and fiber-reinforced soft tissue, exhibits complex responses influenced by fluid pressure, osmotic pressure, and structural mechanics. Existing models struggle to comprehensively represent these intricate interactions and the heterogeneous solid responses within the AF. Additionally, the mechanisms driving differential damage accumulation between non-degenerative and degenerative intervertebral discs remain poorly understood. In this study, we introduce a biphasic-swelling damage model for the AF. We conceptually develop and rigorously validate this model through tissue-level tests employing various loading modes, consistently aligning model predictions with experimental data. Leveraging parametric geometric algorithms and custom Python scripts, we construct models simulating both non-degenerative and degenerative discs. Following calibration, we subject these models to viscous loading protocols. Our findings reveal the posterior AF's susceptibility to damage, contingent upon loading rate and water content. We elucidate the underlying mechanisms by examining the temporal evolution of fluid pressure, osmotic pressure, and the regionally dependent fiber network. This research presents a highly accurate model of the AF, providing valuable insights into disc damage. Future research endeavors should expand this model to incorporate ionic transport and diffusion, enabling a more profound exploration of intervertebral disc mechanobiology. This comprehensive model contributes to a better understanding of AF behavior and may inform therapeutic strategies for disc-related pathologies. STATEMENT OF SIGNIFICANCE: This research presents a comprehensive model of the annulus fibrosus (AF), a crucial component of the intervertebral disc that provides structural support and resists deformation. The study introduces a biphasic-swelling damage model for the AF and validates it through tissue-level tests. The model accounts for fluid pressure, osmotic pressure, and matrix mechanics, providing a more accurate representation of the AF's behavior. The study also investigates the differential damage accumulation between non-degenerative and degenerative discs, shedding light on the mechanisms driving disc degeneration. The findings have significant implications for medical treatments and interventions, as they highlight the posterior AF's susceptibility to damage. This research is of great interest to readers interested in biomechanics, tissue engineering, and medical treatments for disc degeneration.

Impact of Running Exercise on Intervertebral Disc: A Systematic Review

CONTEXT

Running is one of the most popular sports worldwide. However, controversies exist regarding how running affects runner's intervertebral discs (IVD).

OBJECTIVE

The purpose of this study was to systematically review studies that evaluated IVD morphology or composition changes in response to running exercise, to determine the impact of running exercise on IVD.

DATA SOURCES

A systematic literature search was performed for 4 major databases: PubMed, Cochrane, Embase, and Web of Science.

STUDY SELECTION

Inclusion criteria were as follows: (1) healthy people without known IVD disease or major complications such as tuberculosis (IVD degeneration or low back pain are considered as minor complications); (2) subjects performed 1-time or regular running exercises; (3) pre and post comparison of runners or comparison between runners and healthy control subjects; (4) direct or indirect IVD morphology or composition measured; (5) IVD assessed before and after either acute or chronic running exercise, or compared cross-sectionally between runners and controls. Exclusion criteria were as follows: (1) reviews, editorials, letters or abstracts only; (2) animal studies; (3) subjects performed exercise other than running.

STUDY DESIGN

Systematic review.

LEVEL OF EVIDENCE

Level 3.

DATA EXTRACTION

The extracted data included study design and primary outcomes of the included studies. The Newcastle-Ottawa scale (NOS) was used to evaluate study quality and risk of bias.

RESULTS

A total of 13 studies with 632 participants were included in the final analysis; 4 studies measured IVD changes using stature or spinal height, and the other 9 measured IVD changes using magnetic resonance imaging; 6 studies found that running acutely and negatively impacts IVD; 3 out of 5 cross-sectional studies found that IVD parameters are better for runners than controls; 1 longitudinal study found no significant difference in IVD before and after training for marathon in runners; 1 longitudinal study found no significant difference in changes of IVD between runners and controls after 15 years of follow-up.

CONCLUSION

Negative changes in IVD exist for a short period of time after running, which may be due to the temporary compression pushing water content out of the disc. Cross-sectional studies suggest that long-term running exerts a mild positive effect on IVD; however, this inference has not been confirmed by high-quality longitudinal studies.

Mechanical Factors Regulate Annulus Fibrosus (AF) Injury Repair and Remodeling: A Review

Low back pain is a common chronic disease that can severely affect the patient's work and daily life. The breakdown of spinal mechanical homeostasis caused by intervertebral disc (IVD) degeneration is a leading cause of low back pain. Annulus fibrosus (AF), as the outer layer structure of the IVD, is often the first affected part. AF injury caused by consistent stress overload will further accelerate IVD degeneration. Therefore, regulating AF injury repair and remodeling should be the primary goal of the IVD repair strategy. Mechanical stimulation has been shown to promote AF regeneration and repair, but most studies only focus on the effect of single stress on AF, and lack realistic models and methods that can mimic the actual mechanical environment of AF. In this article, we review the effects of different types of stress stimulation on AF injury repair and remodeling, suggest possible beneficial load combinations, and explore the underlying molecular mechanisms. It will provide the theoretical basis for designing better tissue engineering therapy using mechanical factors to regulate AF injury repair and remodeling in the future.

A Novel Fiber-Reinforced Poroviscoelastic Bovine Intervertebral Disc Finite Element Model for Organ Culture Experiment Simulations

Intervertebral disc (IVD) degeneration and methods for repair and regeneration have commonly been studied in organ cultures with animal IVDs under compressive loading. With the recent establishment of a novel multi-axial organ culture system, accurate predictions of the global and local mechanical response of the IVD are needed for control system development and to aid in experiment planning. This study aimed to establish a finite element model of bovine IVD capable of predicting IVD behavior at physiological and detrimental load levels. A finite element model was created based on the dimensions and shape of a typical bovine IVD used in the organ culture. The nucleus pulposus (NP) was modeled as a neo-Hookean poroelastic material and the annulus fibrosus (AF) as a fiber-reinforced poroviscoelastic material. The AF consisted of 10 lamella layers and the material properties were distributed in the radial direction. The model outcome was compared to a bovine IVD in a compressive stress-relaxation experiment. A parametric study was conducted to investigate the effect of different material parameters on the overall IVD response. The model was able to capture the equilibrium response and the relaxation response at physiological and higher strain levels. Permeability and elastic stiffness of the AF fiber network affected the overall response most prominently. The established model can be used to evaluate the response of the bovine IVD at strain levels typical for organ culture experiments, to define relevant boundaries for such studies, and to aid in the development and use of new multi-axial organ culture systems.

Moderate mechanical stimulation antagonizes inflammation of annulus fibrosus cells through YAP-mediated suppression of NF-κB signaling

Intervertebral disc degeneration (IDD) is a leading cause of low back pain. The inflammatory responses caused by aberrant mechanical loading are one of the major factors leading to annulus fibrosus (AF) degeneration and IDD. Previous studies have suggested that moderate cyclic tensile strain (CTS) can regulate anti-inflammatory activities of AF cells (AFCs), and Yes-associated protein (YAP) as a mechanosensitive coactivator senses diverse types of biomechanical stimuli and translates them into biochemical signals controlling cell behaviors. However, it remains poorly understood whether and how YAP mediates the effect of mechanical stimuli on AFCs. In this study, we aimed to investigate the exact effects of different CTS on AFCs as well as the role of YAP signaling involving in it. Our results found that 5% CTS inhibited the inflammatory response and promoted cell growth through inhibiting the phosphorylation of YAP and nuclear localization of NF-κB, while 12% CTS had a significant proinflammatory effect with the inactivation of YAP activity and the activation of NF-κB signaling in AFCs. Furthermore, moderate mechanical stimulation may alleviate the inflammatory reaction of intervertebral discs through YAP-mediated suppression of NF-κB signaling in vivo. Therefore, moderate mechanical stimulation may serve as a promising therapeutic approach for the prevention and treatment of IDD.

Angiopoietin-2 Promotes Mechanical Stress-induced Extracellular Matrix Degradation in Annulus Fibrosus Via the HIF-1α/NF-κB Signaling Pathway

OBJECTIVE

Mechanical stress is an important risk factor for intervertebral disc degeneration (IVDD). Angiopoietin-2 (ANG-2) is regulated by mechanical stress and is widely involved in the regulation of extracellular matrix metabolism. In addition, the signaling cascade between HIF-1α and NF-κB is critical in matrix degradation. This study aims to investigate the role and molecular mechanism of ANG-2 in regulating the degeneration of annulus fibrosus (AF) through the HIF-1α/NF-κB signaling pathway.

METHODS

The bipedal standing mice IVDD model was constructed, and histological experiments were used to evaluate the degree of IVDD and the expression of ANG-2 in the AF. Mouse primary AF cells were extracted in vitro and subjected to mechanical stretching experiments. Western blot assay was used to detect the effect of mechanical stress on ANG-2, and the role of the ANG-2-mediated HIF-1α/NF-κB pathway in matrix degradation. In addition, the effect of inhibiting ANG-2 expression by siRNA or monoclonal antibody on delaying IVDD was investigated at in vitro and in vivo levels. One-way ANOVA with the least significant difference method was used for pairwise comparison of the groups with homogeneous variance, and Dunnett's method was used to compare the groups with heterogeneous variance.

RESULTS

In IVDD, the expressions of catabolic biomarkers (mmp-13, ADAMTS-4) and ANG-2 were significantly increased in AF. In addition, p65 expression was increased while HIF-1α expression was significantly decreased. The results of western blot assay showed mechanical stress significantly up-regulated the expression of ANG-2 in AF cells, and promoted matrix degradation by regulating the activity of HIF-1α/NF-κB pathway. Exogenous addition of Bay117082 and CoCl2 inhibited matrix degradation caused by mechanical stress. Moreover, injection of neutralizing antibody or treatment with siRNA to inhibit the expression of ANG-2 improved the matrix metabolism of AF and inhibited IVDD progression by regulating the HIF-1α/NF-κB signaling pathway.

CONCLUSION

In IVDD, mechanical stress could regulate the HIF-1α/NF-κB signaling pathway and matrix degradation by mediating ANG-2 expression in AF degeneration.

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

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

A swelling-based biphasic analysis on the quasi-static biomechanical behaviors of healthy and degenerative intervertebral discs

BACKGROUND AND OBJECTIVE

The degeneration of intervertebral discs is significantly dependent of the changes in tissue composition ratio and tissue structure. Up to the present, the effects of degeneration on the quasi-static biomechanical responses of discs have not been well understood. The goal of this study is to quantitatively analyze the quasi-static responses of healthy and degenerative discs.

METHODS

Four biphasic swelling-based finite element models are developed and quantitatively validated. Four quasi-static test protocols, including the free-swelling, slow-ramp, creep and stress-relaxation, are implemented. The double Voigt and double Maxwell models are further used to extract the immediate (or residual), short-term and long-term responses of these tests.

RESULTS

Simulation results show that both the swelling-induced pressure in the nucleus pulposus and the initial modulus decrease with degeneration. In the free-swelling test of discs possessing healthy cartilage endplates, simulation results show that over 80% of the total strain is contributed by the short-term response. The long-term response is dominant for discs with degenerated permeability in cartilage endplates. For the creep test, over 50% of the deformation is contributed by the long-term response. In the stress-relaxation test, the long-term stress contribution occupies approximately 31% of total response and is independent of degeneration. Both the residual and short-term responses vary monotonically with degeneration. In addition, both the glycosaminoglycan content and permeability affect the engineering equilibrium time constants of the rheologic models, in which the determining factor is the permeability.

CONCLUSIONS

The content of glycosaminoglycan in intervertebral soft tissues and the permeability of cartilage endplates are two critical factors that affect the fluid-dependent viscoelastic responses of intervertebral discs. The component proportions of the fluid-dependent viscoelastic responses depend also strongly on test protocols. In the slow-ramp test, the glycosaminoglycan content is responsible for the changes of the initial modulus. Since existing computational models simulate disc degenerations only by altering disc height, boundary conditions and material stiffness, the current work highlights the significance of biochemical composition and cartilage endplates permeability in the biomechanical behaviors of degenerated discs.

Determinants of diurnal variation in lumbar intervertebral discs and paraspinal muscles: A prospective quantitative magnetic resonance imaging study

PURPOSE

To prospectively investigate the determinants of diurnal variations in lumbar intervertebral discs and paraspinal muscles.

METHOD

71 females aged 19 ∼ 31 years were examined by morning-evening T2 mapping/diffusion kurtosis imaging (DKI), with weight and lifestyle information (time in night bed-rest [TIB], bed-napping, activity time, and sitting time) assessed by standardized questionnaires. Diurnal shifts in T2, mean diffusivity and mean kurtosis (T2-DS, MD-DS and MK-DS; morning-value minus evening-value) were evaluated for L4-S1 discs (normal, Pfirrmann grade Ⅰ/Ⅱ; degenerative, III/IV). T2 and T2-DS were assessed for L4/5 multifidus and erector spinalis.

RESULTS

For normal discs, bed-napping correlated with MD-DS and MK-DS in disc entirety (p = 0.001 and 0.004); increased activity time suggested higher T2-DS in nucleus pulposus (p = 0.004); prolonged sitting time predicted greater T2-DS in disc entirety and posterior inner annulus fibrosus (PI-AF, p ≤ 0.011); decreased TIB and weight suggested lower T2-DS and higher MK-DS in PI-AF (p = 0.001 ∼ 0.035). For degenerative discs, bed-napping predicted lower T2-DS in nucleus pulposus and PI-AF (p = 0.019); increased TIBsuggested higher T2-DS and lower MK-DS in PI-AF (p = 0.006 and 0.034); longer sitting time predicted higher MK-DS in PI-AF (p = 0.020). Paraspinal muscles exhibited diurnal T2 variation (p < 0.001) which did not correlate with lifestyle factors (p > 0.050).

CONCLUSIONS

Lifestyle and weight have causal effects on the diurnal variation of lumbar discs. Bed-rest may correlate with disc hydration and microstructural stability reserves for subsequent daytime activities. Sitting behavior could induce greater dehydration in normal discs and may alleviate diurnal microstructural rearrangement in degenerative discs. T2 mapping and DKI are promising tools to evaluate disc biomechanics in clinics.

Development of a Computational Model of the Mechanical Behavior of the L4-L5 Lumbar Spine: Application to Disc Degeneration

Degenerative changes in the lumbar spine significantly reduce the quality of life of people. In order to fully understand the biomechanics of the affected spine, it is crucial to consider the biomechanical alterations caused by degeneration of the intervertebral disc (IVD). Therefore, this study is aimed at the development of a discrete element model of the mechanical behavior of the L4-L5 spinal motion segment, which covers all the degeneration grades from healthy IVD to its severe degeneration, and numerical study of the influence of the IVD degeneration on stress state and biomechanics of the spine. In order to analyze the effects of IVD degeneration on spine biomechanics, we simulated physiological loading conditions using compressive forces. The results of modeling showed that at the initial stages of degenerative changes, an increase in the amplitude and area of maximum compressive stresses in the disc is observed. At the late stages of disc degradation, a decrease in the value of intradiscal pressure and a shift in the maximum compressive stresses in the dorsal direction is observed. Such an influence of the degradation of the geometric and mechanical parameters of the tissues of the disc leads to the effect of bulging, which in turn leads to the formation of an intervertebral hernia.

Design of a flexing organ-chip to model in situ loading of the intervertebral disc

The leading cause of disability of all ages worldwide is severe lower back pain. To address this untreated epidemic, further investigation is needed into the leading cause of back pain, intervertebral disc degeneration. In particular, microphysiological systems modeling critical tissues in a degenerative disc, like the annulus fibrosus (AF), are needed to investigate the effects of complex multiaxial strains on AF cells. By replicating these mechanobiological effects unique to the AF that are not yet understood, we can advance therapies for early-stage degeneration at the cellular level. To this end, we designed, fabricated, and collected proof-of-concept data for a novel microphysiological device called the flexing annulus-on-a-chip (AoC). We used computational models and experimental measurements to characterize the device's ability to mimic complex physiologically relevant strains. As a result, these strains proved to be controllable, multi-directional, and uniformly distributed with magnitudes ranging from - 10 % to 12% in the axial, radial, and circumferential directions, which differ greatly from applied strains possible in uniaxial devices. Furthermore, after withstanding accelerated life testing (66 K cycles of 10% strain) and maintaining 2000 bovine AF cells without loading for more than three weeks the AoC proved capable of long-term cell culture. Additionally, after strain (3.5% strain for 75 cycles at 0.5 Hz) was applied to a monolayer of AF cells in the AoC, a population remained adhered to the channel with spread morphology. The AoC can also be tailored for other annular structures in the body such as cardiovascular vessels, lymphatic vessels, and the cervix.

Diurnal Variation in Hydration of the Cervical Intervertebral Disc Assessed Using T2 Mapping of Magnetic Resonance Imaging

Objective

The study aimed to investigate the diurnal variation in cervical disc hydration and its relationship with cervical degeneration.

Materials and Methods

C3–C7 discs of 86 prospectively enrolled participants (37 males, 49 females; mean age ± standard deviation, 23.5 ± 2.5 years) were assessed using T2 mapping in the morning and evening. All discs were stratified by Miyazaki grade or C2–C7 Cobb angle and T2 values (T2). The degree of diurnal T2 variation (T2-DDV), defined as (morning T2 – evening T2)/morning T2 × 100%, was measured for the entire disc, annulus fibrosus (AF), nucleus pulposus (NP), and endplate zones.

Results

T2 of the entire disc decreased significantly after the daytime load (p < 0.001), with a T2-DDV of 13.3% for all discs and 16.0%, 12.2%, and 13.0% for healthy (grade I), mild degenerative (grade II), and advanced degenerative (grade III/IV) discs, respectively. T2 of regional NPs and AFs decreased significantly from morning to evening (p ≤ 0.049) except in the healthy anterior inner AF (p = 0.092). Compared with healthy discs, mild degenerative discs displayed lower T2 and T2-DDV in regional NPs (p < 0.001). Advanced degenerative discs showed higher T2-DDV in the anterior inner AF compared with healthy discs (p = 0.050). Significant diurnal T2 changes in the endplate zones were observed only in healthy discs (p = 0.013). Cervical discs in the low Cobb angle group showed higher T2-DDV in the anterior AFs and anterior NP and lower T2-DDV in the posterior AF than those in the high Cobb angle group (p ≤ 0.041).

Conclusion

This study characterized the diurnal variation in hydration of the cervical discs as assessed using T2 mapping and revealed early chemo-mechanical coupling dysfunction in degenerating discs. Cervical sagittal alignment on MRI can affect the diurnal stress patterns of the cervical discs. T2 mapping is sensitive to disc biomechanical dysfunction and offers translational potential from biomechanical research to clinical application.

Regional variations in discrete collagen fibre mechanics within intact intervertebral disc resolved using synchrotron computed tomography and digital volume correlation

Many soft tissues, such as the intervertebral disc (IVD), have a hierarchical fibrous composite structure which suffers from regional damage. We hypothesise that these tissue regions have distinct, inherent fibre structure and structural response upon loading. Here we used synchrotron computed tomography (sCT) to resolve collagen fibre bundles (∼5μm width) in 3D throughout an intact native rat lumbar IVD under increasing compressive load. Using intact samples meant that tissue boundaries (such as endplate-disc or nucleus-annulus) and residual strain were preserved; this is vital for characterising both the inherent structure and structural changes upon loading in tissue regions functioning in a near-native environment. Nano-scale displacement measurements along >10,000 individual fibres were tracked, and fibre orientation, curvature and strain changes were compared between the posterior-lateral region and the anterior region. These methods can be widely applied to other soft tissues, to identify fibre structures which cause tissue regions to be more susceptible to injury and degeneration. Our results demonstrate for the first time that highly-localised changes in fibre orientation, curvature and strain indicate differences in regional strain transfer and mechanical function (e.g. tissue compliance). This included decreased fibre reorientation at higher loads, specific tissue morphology which reduced capacity for flexibility and high strain at the disc-endplate boundary.

Statement of significance

The analyses presented here are applicable to many collagenous soft tissues which suffer from regional damage. We aimed to investigate regional intervertebral disc (IVD) structural and functional differences by characterising collagen fibre architecture and linking specific fibre- and tissue-level deformation behaviours. Synchrotron CT provided the first demonstration of tracking discrete fibres in 3D within an intact IVD. Detailed analysis of regions was performed using over 200k points, spaced every 8 μm along 10k individual fibres. Such comprehensive structural characterisation is significant in informing future computational models. Morphological indicators of tissue compliance (change in fibre curvature and orientation) and fibre strain measurements revealed localised and regional differences in tissue behaviour.

A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

The Potential Role of Cytokines in Diabetic Intervertebral Disc Degeneration

Intervertebral disc degeneration (IVDD) is a major cause of low back pain. Diabetes mellitus is a chronic inflammatory disease that may cause or aggravate IVDD; however, the mechanism by which diabetes induce IVDD is currently unclear. Compared to non-diabetic individuals, diabetic patients have higher levels of plasma cytokines, especially TNF-α, IL-1β, IL-5, IL-6, IL-7, IL-10, and IL-18. Due to the crucial role of cytokines in the process of intervertebral disc degeneration, we hypothesized that elevation of these cytokines in plasma of diabetic patients may be involved in the process of diabetes-induced IVDD. In this review, changes in plasma cytokine levels in diabetic patients were summarized and the potential role of elevated cytokines in diabetes-induced IVDD was discussed. Results showed that some cytokines such as TNF-α and IL-1β may accelerate the development of IVDD, while others such as IL-10 is supposed to prevent its development. Apoptosis, senescence, and extracellular matrix metabolism were found to be regulated by these cytokines in IVDD. Further studies are required to validate the cytokines targeted strategy for diabetic IVDD therapy.

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.

Tenomodulin and Chondromodulin-1 Are Both Required to Maintain Biomechanical Function and Prevent Intervertebral Disc Degeneration

OBJECTIVE

The underlying mechanisms and molecular factors influencing intervertebral disc (IVD) homeostasis and degeneration remain clinically relevant. Tenomodulin (Tnmd) and chondromodulin (Chm1) are antiangiogenic transmembrane glycoproteins, with cleavable C-terminus, expressed by IVD cells that are implicated in the onset of degenerative processes. We evaluate the organ-level biomechanical impact of knocking out Tnmd alone, and Tnmd and Chm1, simultaneously.

DESIGN

Caudal (c5-8) and lumbar vertebrae (L1-4) of skeletally mature male and female 9-month-old wildtype (WT), Tnmd knockout (Tnmd^-/-), and Tnmd/Chm1 double knockout (Tnmd^-/-/Chm^-/-) mice were used (n = 9-13 per group). Disc height index (DHI), histomorphological changes, and axial, torsional, creep, and failure biomechanical properties were evaluated. Differences were assessed by one-way ANOVA with post hoc Bonferroni-corrected comparisons (P < 0.05).

RESULTS

Tnmd^-/-/Chm1^-/- IVDs displayed increased DHI and histomorphological scores that indicated increased IVD degeneration compared to the WT and Tnmd^-/- groups. Double knockout IVDs required significantly less torque and energy to initiate torsional failure. Creep parameters were comparable between all groups, except for the slow time constant, which indicated faster outward fluid flow. Tnmd^-/- IVDs lost fluid faster than the WT group, and this effect was amplified in the double knockout IVDs.

CONCLUSION

Knocking out Tnmd and Chm1 affects IVD fluid flow and organ-level biomechanical function and therefore may play a role in contributing to IVD degeneration. Larger effects of the Tnmd and Chm1 double knockout mice compared to the Tnmd single mutant suggest that Chm1 may play a compensatory role in the Tnmd single mutant IVDs.

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.

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.

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.

Multiaxial validation of a finite element model of the intervertebral disc with multigenerational fibers to establish residual strain

Finite element models of the intervertebral disc are used to address research questions that cannot be tested through typical experimentation. A disc model requires complex geometry and tissue properties to be accurately defined to mimic the physiological disc. The physiological disc possesses residual strain in the annulus fibrosus (AF) due to osmotic swelling and due to inherently pre-strained fibers. We developed a disc model with residual contributions due to swelling-only, and a multigeneration model with residual contributions due to both swelling and AF fiber pre-strain and validated it against organ-scale uniaxial, quasi-static and multiaxial, dynamic mechanical tests. In addition, we demonstrated the models' ability to mimic the opening angle observed following radial incision of bovine discs. Both models were validated against organ-scale experimental data. While the swelling only model responses were within the experimental 95% confidence interval, the multigeneration model offered outcomes closer to the experimental mean and had a bovine model opening angle within one SD of the experimental mean. The better outcomes for the multigeneration model, which allowed for the inclusion of inherently pre-strained fibers in AF, is likely due to its uniform fiber contribution throughout the AF. We conclude that the residual contribution of pre-strained fibers in the AF should be included to best simulate the physiological disc and its behaviors.

Microgravity and Radiation Effects on Astronaut Intervertebral Disc Health

INTRODUCTION: The effects of spaceflight on the intervertebral disc (IVD) have not been thoroughly studied, despite the knowledge that spaceflight increases the risk of herniation of IVDs in astronauts upon return to Earth. However, as long duration missions become more common, fully characterizing the mechanisms behind space-induced IVD degeneration becomes increasingly imperative for mission success. This review therefore surveys current literature to outline the results of human, animal, and cell-level studies investigating the effect of microgravity and radiation exposure on IVD health. Overall, recurring study findings include increases in IVD height in microgravity conditions, upregulation of catabolic proteases leading to a weakening extracellular matrix (ECM), and both nucleus pulposus (NP) swelling and loss of annulus fibrosus (AF) fiber alignment which are hypothesized to contribute to the increased risk of herniation when reloading is experienced. However, the limitations of current studies are also discussed. For example, human studies do not allow for invasive measures of the underpinning biochemical mechanisms, correlating animal model results to the human condition may be difficult, and cellular studies lack incorporation of ECM and other complexities that mimic the native IVD microarchitecture and environment. Moving forward, the use of three-dimensional organoid culture models that incorporate IVD-specific human cells, ECM, and signals as well as the development of cell- and ECM-level computational models may further improve our understanding of the impacts that spaceflight has on astronaut IVD health.Smith K, Mercuri J. Microgravity and radiation effects on astronaut intervertebral disc health. Aerosp Med Hum Perform. 2021; 92(5):342352.

Intervertebral Disc Mechanics With Nucleotomy: Differences Between Simple and Dual Loading

Painful herniated discs are treated surgically by removing extruded nucleus pulposus (NP) material (nucleotomy). NP removal through enzymatic digestion is also commonly performed to initiate degenerative changes to study potential biological repair strategies. Experimental and computational studies have shown a decrease in disc stiffness with nucleotomy under single loading modalities, such as compression-only or bending-only loading. However, studies that apply more physiologically relevant loading conditions, such as compression in combination with bending or torsion, have shown contradicting results. We used a previously validated bone-disc-bone finite element model (Control) to create a Nucleotomy model to evaluate the effect of dual loading conditions (compression with torsion or bending) on intradiscal deformations. While disc joint stiffness decreased with nucleotomy under single loading conditions, as commonly reported in the literature, dual loading resulted in an increase in bending stiffness. More specifically, dual loading resulted in a 40% increase in bending stiffness under flexion and extension and a 25% increase in stiffness under lateral bending. The increase in bending stiffness was due to an increase and shift in compressive stress, where peak stresses migrated from the NP-annulus interface to the outer annulus. In contrast, the decrease in torsional stiffness was due to greater fiber reorientation during compression. In general, large radial strains were observed with nucleotomy, suggesting an increased risk for delamination or degenerative remodeling. In conclusion, the effect of nucleotomy on disc mechanics depends on the type and complexity of applied loads.

Spine biomechanical testing methodologies: The controversy of consensus vs scientific evidence

Biomechanical testing methodologies for the spine have developed over the past 50 years. During that time, there have been several paradigm shifts with respect to techniques. These techniques evolved by incorporating state-of-the-art engineering principles, in vivo measurements, anatomical structure-function relationships, and the scientific method. Multiple parametric studies have focused on the effects that the experimental technique has on outcomes. As a result, testing methodologies have evolved, but there are no standard testing protocols, which makes the comparison of findings between experiments difficult and conclusions about in vivo performance challenging. In 2019, the international spine research community was surveyed to determine the consensus on spine biomechanical testing and if the consensus opinion was consistent with the scientific evidence. More than 80 responses to the survey were received. The findings of this survey confirmed that while some methods have been commonly adopted, not all are consistent with the scientific evidence. This review summarizes the scientific literature, the current consensus, and the authors' recommendations on best practices based on the compendium of available evidence.

Level-wise differences in in vivo lateral bending moment are associated with microstructural alterations in bovine caudal intervertebral discs

Despite its common use as a laboratory model, little is known about the in vivo forces and moments applied to the bovine caudal intervertebral disc. Such aspects are crucial, as intervertebral disc tissue is known to remodel in response to repeated loading. We hypothesized that the magnitude of loading from muscle contraction during a typical lateral bending motion varies between caudal levels and is accompanied by variations in tissue microstructure. This hypothesis was tested by estimating level-wise forces and bending moments using two independent approaches: a dynamic analytical model of the motion and analysis of muscle cross-sections obtained via computed tomography. Microstructure was assessed by measuring the collagen fiber crimp period in the annulus fibrosus, and composition was assessed via quantitative histology. Both the analytical model and muscle cross-sections indicated peak bending moments of over 3 N m and peak compressive force of over 125 N at the c1c2 level, decreasing distally. There was a significant downward trend from proximal to distal in the outer annulus fibrosus collagen crimp period in the anterior and posterior regions only, suggesting remodeling in response to the highest lateral bending moments. There were no observed trends in composition. Our results suggest that although the proximal discs in the bovine tail are subjected to forces and moments from muscle contraction that are comparable (relative to disc size) to those acting on human lumbar discs, the distal discs are not. The resulting pattern of microstructural alterations suggests that level-wise differences should be considered when using bovine discs as a research model.

What Is Finite Element Analysis?

Finite element analysis is a computational technique to predict how different materials will react when a range of forces are applied. In the field of orthopedics, this technique has predominantly been used for implant design and testing. As the technology improves, increasing clinical applications are being developed, offering promise in the areas of surgical planning and the opportunity to tailor implants to individual patient characteristics. This article introduces the various preclinical mechanical tests available, as well as providing a brief overview of the finite element analysis technology.

Transient swelling behavior of the bovine caudal disc

The intervertebral disc is an avascular composite structure, comprised of the nucleus pulposus (NP) and the annulus fibrosus (AF). Previous tissue-level experiments either examined relative differences in swelling capacity of the two disc regions at a single time point or tested explant structures that did not replicate in situ boundary conditions. Previous joint-level studies that investigated time-dependent fluid flow into the disc provided limited information about swelling-induced intradiscal strains with respect to time and boundary constraints. Therefore, the objective of this study was to investigate time-dependent swelling behavior of the intervertebral disc ex situ. The first study investigated time-dependent free-swelling response of the whole disc and the disc's subcomponents separately (i.e., NP and AF). Findings from this study showed that the swelling rate and swelling capacity of NP explants under free-swelling conditions were greater than AF explants. The second study evaluated the effect of boundary conditions on in-plane strain distributions of intact discs and AF rings. Swelling-induced strain was highly heterogeneous in AF rings, where negative circumferential strains were observed in the inner AF and tensile circumferential strains were observed in the outer AF. Radial strains in AF rings were an order of magnitude greater than circumferential strains. Restricting fluid flow only to the outer AF periphery reduced the swelling of the inner AF. Swelling of intact discs affected both NP and AF swelling behaviors, where NP hydration decreased by 60%. Furthermore, the presence of the NP reduced peak radial strains in the AF and resulted in uniform strain distribution throughout the AF. In conclusion, these studies highlight that tissue hydration and swelling-induced strains largely depend on regional biochemical composition and geometric boundary constraints.

Finite element models can reproduce the effect of nucleotomy on the multi-axial compliance of human intervertebral discs

Finite element (FE) models can unravel the link between intervertebral disc (IVD) degeneration and its mechanical behaviour. Nucleotomy may provide the data required for model verification. Three human IVDs were scanned with MRI and tested in multiple loading scenarios, prior and post nucleotomy. The resulting data was used to generate, calibrate, and verify the FE models. Nucleotomy increased the experimental range of motion by 26%, a result reproduced by the FE simulation within a 5% error. This work demonstrates the ability of FE models to reproduce the mechanical compliance of human IVDs prior and post nucleotomy.

Alterations in the Geometry, Fiber Orientation, and Mechanical Behavior of the Lumbar Intervertebral Disc by Nucleus Swelling

Soft tissues observed in clinical medical images are often prestrained in tension by internal pressure or tissue hydration. For a native disc, nucleus swelling occurs in equilibrium with osmotic pressure induced by the high concentration of proteoglycan in the nucleus. The objective of this computational study was to investigate the effects of nucleus swelling on disc geometry, fiber orientation, and mechanical behavior by comparing those of prestrained and zero-pressure (unswelled) discs. Thermoelastic analysis techniques were repurposed in order to determine the zero-pressure disc geometry which, when pressurized, matches the prestrained disc geometry observed in clinical images. The zero-pressure geometry was then used in simulations to approximately represent a degenerated disc, which loses the ability of nucleus swelling but has not undergone distinct soft tissue remodeling/disruption. Our simulation results demonstrated that the loss of nucleus swelling caused a slight change in the disc geometry and fiber orientation, but a distinct deterioration in the resistance to intervertebral rotations including sagittal bending, lateral bending, and axial torsion. Different from rotational loading, in compression (with a displacement of 0.45 mm applied), a much larger stiffness (3.02 KN/mm) and a greater intradiscal pressure (IDP) (0.61 MPa) were measured in the zero-pressure disc, compared to the prestrained disc (1.41 KN/mm and 0.52 MPa). This computational study could be useful to understand mechanisms of disc degeneration, and guide the future design of disc tissue engineering material and biomimic disc implants.

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.