Human annulus fibrosus material properties from biaxial testing and constitutive modeling are altered with degeneration

Biomechanics of annulus fibrosus: Elastic fiber simplification and degenerative impact on damage initiation and propagation

This study addresses three primary objectives related to lumbar intervertebral disc (IVD) biomechanics under ramping quasi-static loading conditions. First, we explore the conditions justifying the simplification of axisymmetric elastic fiber families into single fiber bundles through discretized strain energy functions. Simulations reveal that a concentration factor exceeding 10 allows for a consistent deviation below 10% between simplified and non-simplified responses. Second, we investigate the impact of elastic fibers on the physiological stiffness in IVDs, revealing minimal influence on biological motions but significant effects on degeneration. Lastly, we examine the initiation and progression of annulus fibrosus (AF) damage. Our findings confirm the validity of simplifying elastic fiber families and underscore the necessity of considering elastic fiber damage in biomechanical studies of AF tissues. Elastic fibers contribute to increased biaxial stretch stiffness, and their damage significantly affects the loading capacity of the inner AF. Additionally, degeneration significantly alters the susceptibility to damage in the AF, with specific regions exhibiting higher vulnerability. Damage tends to extend circumferentially and radially, emphasizing the regional variations in collagen and elastic fiber properties. This study offers useful insights for refining biomechanical models, paving the way for a more comprehensive understanding of IVD responses and potential clinical implications.

Comparative FEM study on intervertebral disc modeling: Holzapfel-Gasser-Ogden vs. structural rebars

Introduction: Numerical modeling of the intervertebral disc (IVD) is challenging due to its complex and heterogeneous structure, requiring careful selection of constitutive models and material properties. A critical aspect of such modeling is the representation of annulus fibers, which significantly impact IVD biomechanics. This study presents a comparative analysis of different methods for fiber reinforcement in the annulus fibrosus of a finite element (FE) model of the human IVD. Methods: We utilized a reconstructed L4-L5 IVD geometry to compare three fiber modeling approaches: the anisotropic Holzapfel-Gasser-Ogden (HGO) model (HGO fiber model) and two sets of structural rebar elements with linear-elastic (linear rebar model) and hyperelastic (nonlinear rebar model) material definitions, respectively. Prior to calibration, we conducted a sensitivity analysis to identify the most important model parameters to be calibrated and improve the efficiency of the calibration. Calibration was performed using a genetic algorithm and in vitro range of motion (RoM) data from a published study with eight specimens tested under four loading scenarios. For validation, intradiscal pressure (IDP) measurements from the same study were used, along with additional RoM data from a separate publication involving five specimens subjected to four different loading conditions. Results: The sensitivity analysis revealed that most parameters, except for the Poisson ratio of the annulus fibers and C01 from the nucleus, significantly affected the RoM and IDP outcomes. Upon calibration, the HGO fiber model demonstrated the highest accuracy (R2 = 0.95), followed by the linear (R2 = 0.89) and nonlinear rebar models (R2 = 0.87). During the validation phase, the HGO fiber model maintained its high accuracy (RoM R2 = 0.85; IDP R2 = 0.87), while the linear and nonlinear rebar models had lower validation scores (RoM R2 = 0.71 and 0.69; IDP R2 = 0.86 and 0.8, respectively). Discussion: The results of the study demonstrate a successful calibration process that established good agreement with experimental data. Based on our findings, the HGO fiber model appears to be a more suitable option for accurate IVD FE modeling considering its higher fidelity in simulation results and computational efficiency.

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.

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.

Applied use of biomechanical measurements from human tissues for the development of medical skills trainers: a scoping review

OBJECTIVE

The objective of this review was to identify quantitative biomechanical measurements of human tissues, the methods for obtaining these measurements, and the primary motivations for conducting biomechanical research.

INTRODUCTION

Medical skills trainers are a safe and useful tool for clinicians to use when learning or practicing medical procedures. The haptic fidelity of these devices is often poor, which may be because the synthetic materials chosen for these devices do not have the same mechanical properties as human tissues. This review investigates a heterogeneous body of literature to identify which biomechanical properties are available for human tissues, the methods for obtaining these values, and the primary motivations behind conducting biomechanical tests.

INCLUSION CRITERIA

Studies containing quantitative measurements of the biomechanical properties of human tissues were included. Studies that primarily focused on dynamic and fluid mechanical properties were excluded. Additionally, studies only containing animal, in silico , or synthetic materials were excluded from this review.

METHODS

This scoping review followed the JBI methodology for scoping reviews and the Preferred Reporting Items for Systematic Reviews and Meta-Analyses extension for Scoping Reviews (PRISMA-ScR). Sources of evidence were extracted from CINAHL (EBSCO), IEEE Xplore, MEDLINE (PubMed), Scopus, and engineering conference proceedings. The search was limited to the English language. Two independent reviewers screened titles and abstracts as well as full-text reviews. Any conflicts that arose during screening and full-text review were mediated by a third reviewer. Data extraction was conducted by 2 independent reviewers and discrepancies were mediated through discussion. The results are presented in tabular, figure, and narrative formats.

RESULTS

Data were extracted from a total of 186 full-text publications. All of the studies, except for 1, were experimental. Included studies came from 33 countries, with the majority coming from the United States. Ex vivo methods were the predominant approach for extracting human tissue samples, and the most commonly studied tissue type was musculoskeletal. In this study, nearly 200 unique biomechanical values were reported, and the most commonly reported value was Young's (elastic) modulus. The most common type of mechanical test performed was tensile testing, and the most common reason for testing human tissues was to characterize biomechanical properties. Although the number of published studies on biomechanical properties of human tissues has increased over the past 20 years, there are many gaps in the literature. Of the 186 included studies, only 7 used human tissues for the design or validation of medical skills training devices. Furthermore, in studies where biomechanical values for human tissues have been obtained, a lack of standardization in engineering assumptions, methodologies, and tissue preparation may implicate the usefulness of these values.

CONCLUSIONS

This review is the first of its kind to give a broad overview of the biomechanics of human tissues in the published literature. With respect to high-fidelity haptics, there is a large gap in the published literature. Even in instances where biomechanical values are available, comparing or using these values is difficult. This is likely due to the lack of standardization in engineering assumptions, testing methodology, and reporting of the results. It is recommended that journals and experts in engineering fields conduct further research to investigate the feasibility of implementing reporting standards.

REVIEW REGISTRATION

Open Science Framework https://osf.io/fgb34.

Polymeric Hydrogels for Intervertebral Disc Replacement/Integration: Playing with the Chemical Composition for Tuning Shear Behavior and Hydrophilicity

Damages to the intervertebral disc (IVD) due to improper loading or degeneration result in back pain, which is a common disease affecting an increasing number of patients. Different strategies for IVD remediation have been developed, from surgical treatment to disc replacement, by using both metallic and non-metallic materials. Hydrogels are very attractive materials due to their ability to simulate the properties of many soft tissues; moreover, their chemical composition can be varied in order to assure performances similar to the natural disc. In particular, for the replacement of the IVD outer ring, namely, the anulus fibrosus, the shear properties are of paramount importance. In this work, we produced hydrogels through the photo-induced crosslinking of different mixtures composed of two hydrophilic monofunctional and difunctional polymers, namely, poly(ethyleneglycol) methyl ether methacrylate (PEGMEMA) and poly(ethyleneglycol) dimethacrylate (PEGDMA), together with a hydrophobic molecule, i.e., tert-butyl acrylate (tBA). By changing the ratio among the precursors, we demonstrated the tunability of both the shear properties and hydrophilicity. The structural properties of hydrogels were studied by solid-state nuclear magnetic resonance (NMR). These experiments provided insights on both the structure and molecular dynamics of polymeric networks and, together with information obtained by differential scanning calorimetry (DSC), allowed for correlating the physical properties of the hydrogels with their chemical composition.

Fiber orientation and crimp level might control the auxetic effect of biological tissues

We propose an analytical micromechanical model for studying the lamellar-composite-like structure of fibrous soft tissue. The tissue under consideration is made up of several lamellae, and is designed to resemble the annulus fibrosus (AF) tissue or media layer of arterial tissue, for example. The collagen fibers are arranged in parallel in each lamella and the fiber orientation differs from one lamella to its neighbors. The parallel fibers in each lamella of AF tissue, for example, have been observed to have a crimped microstructure. The proposed model incorporates this quality, considering fiber waviness as a sinusoidal shape and taking into account the fiber dispersion in different layers, where both fiber and matrix are considered as solid phases. We find that collagen-fiber waviness and layer orientation have a significant influence on Poisson's ratio. The effective Poisson's ratio predicted by the proposed model demonstrates that the crimped collagen fiber microstructure might weaken the auxetic effect of fibrous soft tissue, which might explain why, as the literature suggests, the auxetic behavior is more difficult to observe than large Poisson's ratios. As opposed to the many studies that use the well-known hyperelastic fiber-based constitutive model, in which out-of-plane expansion is often observed, the present work explains the auxetic response found in modeling and in experimental data from the perspective of collagen fiber microstructure.

In Situ Covalent Reinforcement of a Benzene-1,3,5-Tricarboxamide Supramolecular Polymer Enables Biomimetic, Tough, and Fibrous Hydrogels and Bioinks

Synthetic hydrogels often lack the load-bearing capacity and mechanical properties of native biopolymers found in tissue, such as cartilage. In natural tissues, toughness is often imparted via the combination of fibrous noncovalent self-assembly with key covalent bond formation. This controlled combination of supramolecular and covalent interactions remains difficult to engineer, yet can provide a clear strategy for advanced biomaterials. Here, a synthetic supramolecular/covalent strategy is investigated for creating a tough hydrogel that embodies the hierarchical fibrous architecture of the extracellular matrix (ECM). A benzene-1,3,5-tricarboxamide (BTA) hydrogelator is developed with synthetically addressable norbornene handles that self-assembles to form a and viscoelastic hydrogel. Inspired by collagen's covalent cross-linking of fibrils, the mechanical properties are reinforced by covalent intra- and interfiber cross-links. At over 90% water, the hydrogels withstand up to 550% tensile strain, 90% compressive strain, and dissipated energy with recoverable hysteresis. The hydrogels are shear-thinning, can be 3D bioprinted with good shape fidelity, and can be toughened via covalent cross-linking. These materials enable the bioprinting of human mesenchymal stromal cell (hMSC) spheroids and subsequent differentiation into chondrogenic tissue. Collectively, these findings highlight the power of covalent reinforcement of supramolecular fibers, offering a strategy for the bottom-up design of dynamic, yet tough, hydrogels and bioinks.

Textile Design of an Intervertebral Disc Replacement Device from Silk Yarn

Low back pain is often due to degeneration of the intervertebral discs (IVD). It is one of the most common age- and work-related problems in today's society. Current treatments are not able to efficiently restore the full function of the IVD. Therefore, the aim of the present work was to reconstruct the two parts of the intervertebral disc-the annulus fibrosus (AF) and the nucleus pulposus (NP)-in such a way that the natural structural features were mimicked by a textile design. Silk was selected as the biomaterial for realization of a textile IVD because of its cytocompatibility, biodegradability, high strength, stiffness, and toughness, both in tension and compression. Therefore, an embroidered structure made of silk yarn was developed that reproduces the alternating fiber structure of +30° and -30° fiber orientation found in the AF and mimics its lamellar structure. The developed embroidered ribbons showed a tensile strength that corresponded to that of the natural AF. Fiber additive manufacturing with 1 mm silk staple fibers was used to replicate the fiber network of the NP and generate an open porous textile 3D structure that may serve as a reinforcement structure for the gel-like NP.

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

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

Anatomically and Biomechanically Relevant Monolithic Total Disc Replacement Made of 3D-Printed Thermoplastic Polyurethane

Various implant treatments, including total disc replacements, have been tried to treat lumbar intervertebral disc (IVD) degeneration, which is claimed to be the main contributor of lower back pain. The treatments, however, come with peripheral issues. This study proposes a novel approach that complies with the anatomical features of IVD, the so-called monolithic total disc replacement (MTDR). As the name suggests, the MTDR is a one-part device that consists of lattice and rigid structures to mimic the nucleus pulposus and annulus fibrosus, respectively. The MTDR can be made of two types of thermoplastic polyurethane (TPU 87A and TPU 95A) and fabricated using a 3D printing approach: fused filament fabrication. The MTDR design involves two configurations—the full lattice (FLC) and anatomy-based (ABC) configurations. The MTDR is evaluated in terms of its physical, mechanical, and cytotoxicity properties. The physical characterization includes the geometrical evaluations, wettability measurements, degradability tests, and swelling tests. The mechanical characterization comprises compressive tests of the materials, an analytical approach using the Voigt model of composite, and a finite element analysis. The cytotoxicity assays include the direct assay using hemocytometry and the indirect assay using a tetrazolium-based colorimetric (MTS) assay. The geometrical evaluation shows that the fabrication results are tolerable, and the two materials have good wettability and low degradation rates. The mechanical characterization shows that the ABC-MTDR has more similar mechanical properties to an IVD than the FLC-MTDR. The cytotoxicity assays prove that the materials are non-cytotoxic, allowing cells to grow on the surfaces of the materials.

Effect of degenerative factors on cervical spinal cord during flexion and extension: a dynamic finite element analysis

Spinal cord injury (SCI) is a global problem that brings a heavy burden to both patients and society. Recent investigations indicated degenerative disease is taking an increasing part in SCI with the growth of the aging population. However, little insight has been gained about the effect of cervical degenerative disease on the spinal cord during dynamic activities. In this work, a dynamic fluid-structure interaction model was developed and validated to investigate the effect of anterior and posterior encroachment caused by degenerative disease on the spinal cord during normal extension and flexion. Maximum von-Mises stress and maximum principal strain were observed at the end of extension and flexion. The abnormal stress distribution caused by degenerative factors was concentrated in the descending tracts of the spinal cord. Our finding indicates that the excessive motion of the cervical spine could potentially exacerbate spinal cord injury and enlarge injury areas. Stress and strain remained low compared to extension during moderate flexion. This suggests that patients with cervical degenerative disease should avoid frequent or excessive flexion and extension which could result in motor function impairment, whereas moderate flexion is safe. Besides, encroachment caused by degenerative factors that are not significant in static imaging could also cause cord compression during normal activities.

Reproducing Morphological Features Of Intervertebral Disc Using Finite Element Modeling To Predict The Course Of Cervical Spine Dorsopathy

Study objective — To evaluate how morphological features of intervertebral disc would affect the outcomes of finite element modeling of axial load in the cervical spine, C3-C5, in order to predict the risk of occurrence and course of dorsopathies. Material and Methods — Three-dimensional models of the cervical spine vertebrae were generated from the computed tomography data of a volunteer (24 years old male without detected pathology of his neck). Intervertebral disc models were developed in two configurations. For each model, we performed a finite element investigation of the stress-strain state with the same loading conditions. The load-displacement curves were compared with the experimental data generated from the results of previously conducted in vitro experiments. Results — The maximum and mean displacement values for the isotropic model were 1.15 mm and 0.73 ± 0.45 mm, respectively. For anisotropic model, maximum and mean displacement values were 0.86 mm and 0.47 ± 0.24 mm, correspondingly. Predicted displacement values for both models matched the experimental data fairly well. Stress profiles of intervertebral discs and stress diagrams of facet joints were calculated. Conclusion — The proposed geometric and constitutive configurations of the intervertebral disc take into account specific morphological features at low computational costs, thereby facilitating the modeling of degenerative disc changes.

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.

Computational modeling of lumbar disc degeneration before and after spinal fusion

BACKGROUND

Advancing age and degeneration frequently lead to low back pain, which is the most prevalent musculoskeletal disorder worldwide. Degenerative changes in intervertebral discs and musculo-ligamentous incapacity to compensate sagittal imbalance are typically amongst the sources of instability, with spinal fusion techniques being the main treatment options to relieve pain. The aims of this work were to: (i) assess the link between ligament degeneration and spinal instability by determining the role of each ligament per movement, (ii) evaluate the impact of disc height reduction in degenerative changes, and (iii) unveil the most advantageous type of posterior fixation in Oblique Lumbar Interbody Fusion to prevent adjacent disc degeneration.

METHODS

Two L3-L5 finite element models were developed, being the first in healthy condition and the second having reduced L4-L5 height. Different degrees of degeneration were tested, combined with different fixation configurations for Oblique Lumbar Interbody Fusion.

FINDINGS

Facet capsular ligament and anterior longitudinal ligament were the most influential ligaments for spinal stability, particularly with increasing degeneration and disc height reduction. Pre-existent degeneration had lower influence than the fusion procedure for the risk of adjacent disc degeneration, being the highest stability and minimal degeneration achieved with bilateral fixation. Right unilateral fixation was more suited to reduce disc stress than left unilateral fixation.

INTERPRETATION

Bilateral fixation is the best option to stabilize the spinal segment, but unilateral right fixation may suffice. This has direct implications for clinical practice, and the extension to a population-based study will allow for more efficient fusion surgeries.

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

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

Biomechanical analysis of the cervical spine segment as a method for studying the functional and dynamic anatomy of the human neck

BACKGROUND

Traditionally, dynamic and functional anatomy, in particular the dynamic anatomy of the neck, is studied on cadaveric material. However, the development of in vivo visualization technologies and in silico modeling has made it possible to expand these possibilities. Despite significant progress in the study of dynamic and functional anatomy of the neck by means of in silico methods, the issues of validating the developed models and taking into account the pronounced nonlinearity of soft tissues as well as local anisotropy remain open. The aim of this study was to develop a virtual dynamic anatomical model of the human neck and reproduce the dynamic processes in the cervical spine from this model using the finite element method.

MATERIALS AND METHODS

Reverse engineering was used to generate a dynamic anatomical model of the neck from CT data (both male, 24 and 22 years old). Two segments of the cervical spine (C3-C5, C2-T1) were isolated from the resulting model for finite element analysis. Finite element mesh generation and contact interactions were performed using the HyperMesh software (Altair Engineering Inc, Troy, Michigan, USA). The anisotropic hyperelastic Holzapfel-Gasser-Ogden model was used to describe the material behavior of the fibrous rings of the disc. Material modeling and finite element analysis were performed using Abaqus CAE 6.14 software (Simulia, Johnston, Rhode Island, USA).

RESULTS

A technique for creating a virtual dynamic anatomical model of the neck was elaborated and implemented. The model includes 79 major anatomical structures of the neck segmented from radiological data. A finite element analysis of the cervical spine was performed. The results of finite element analysis of the C3-C5 segment under axial load were compared with in vitro data. The proposed model shows nonlinear deformation of the disc under static loading; the model predicted displacement values agree well with the experimental ones. The displacement of the С3-С5 central vertebra with an axial load of 800 N reaches a value of 0.65 mm. For the segment C2-T1, data on intradiscal pressure, stress plots and displacements during flexion were obtained. The maximum stress value of 10.036 MPa is observed in the C3-C4 disc.

CONCLUSION

Simulation results using the proposed methodology are in good agreement with experimental data. The generated biomechanical models allow describing dynamic phenomena in the cervical spine and obtaining a wide range of quantitative properties of anatomical objects, which are otherwise inaccessible to classical methods for studying dynamic and functional anatomy.

Reflected cross-polarized light microscopy as a method for measuring collagen fiber crimp in musculoskeletal tissues

Many musculoskeletal tissues are composed primarily of type I collagen, which takes on a periodic crimp morphology that allows large tensile strains in the tissue. The spatial period of collagen fiber crimp may be used to infer internal strains in a tissue and is typically measured using transmitted cross-polarized light imaging of thin slices. However, slicing may induce specimen distortion and precludes mechanical loading of the specimen during imaging. We hypothesized that reflected cross-polarized light imaging of thick tissue explants would yield crimp period measurements comparable to those obtained from transmitted light imaging of thin slices. We further hypothesized that these measurements would be sensitive to applied uniaxial strain in the fiber direction. These hypotheses were tested by imaging both intervertebral disc outer annulus fibrosus and medial collateral ligament tissue specimens. We found that both transmitted and reflected light yielded similar crimp period measurements for intervertebral disc tissue, with an overall average of 43.5 ± 11.5 μm. Reflected light yielded a significantly higher crimp period with lower variance than transmission through thin specimens (54.1 ± 10.6 μm versus 50.4 ± 16.0 μm) in the ligament. Upon application of axial tension, crimp periods in both fibers increased at a rate of approximately three times the applied strain (with 3.17% applied strain yielding a 9.64 ± 4.4% increase in crimp period in the disc and an 11.7 ± 3.7% increase in the ligament), indicating significant fibril sliding. In support of our hypotheses, these findings suggest that reflected cross-polarized light is a suitable method for measuring collagen fiber crimp in musculoskeletal tissues, both statically and under tension.

NONLINEAR DYNAMIC CHARACTERISTICS OF LUMBAR DISC SUBJECTED TO STOCHASTIC EXCITATION

The lumbar disc has complex structure and material properties, and if it develops disease, it cannot heal itself. Therefore, prevention of lumbar disc herniation is very important. In this paper, the model of lumbar disc is built, and the nonlinear dynamic response of the system is researched. A modified Van der Pol model is imported to describe the stress–strain curves of the lumbar disc. The system’s dynamic model is set up, and harmonic balance method is applied to revise the natural frequency of the system. The product of numerical simulation reveals that the lumbar disc has complex dynamic characteristics, including balance point, limit cycle bifurcation and stochastic Hopf bifurcation. By changing the parameters, we can avoid large-scale periodic movement of the lumbar disc which causes lumbar disc herniation. These results contribute to the prevention of lumbar disc herniation.

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

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

Modeling of human intervertebral disc annulus fibrosus with complex multi-fiber networks

Collagen fibers within the annulus fibrosus (AF) lamellae are unidirectionally aligned with alternating orientations between adjacent layers. AF constitutive models often combine two adjacent lamellae into a single equivalent layer containing two fiber networks with a crisscross pattern. Additionally, AF models overlook the inter-lamellar matrix (ILM) as well as elastic fiber networks in between lamellae. We developed a nonhomogenous micromechanical model as well as two coarser homogenous hyperelastic and microplane models of the human AF, and compared their performances against measurements (tissue level uniaxial and biaxial tests as well as whole disc experiments) and seven published hyperelastic models. The micromechanical model had a realistic non-homogenous distribution of collagen fiber networks within each lamella and elastic fiber network in the ILM. For small matrix linear moduli (<0.2 MPa), the ILM showed substantial anisotropy (>10%) due to the elastic fiber network. However, at moduli >0.2 MPa, the effects of the elastic fiber network on differences in stress-strain responses at different directions disappeared (<10%). Variations in sample geometry and boundary conditions (due to uncertainty) markedly affected stress-strain responses of the tissue in uniaxial and biaxial tests (up to 16 times). In tissue level tests, therefore, simulations should represent testing conditions (e.g., boundary conditions, specimen geometry, preloads) as closely as possible. Stress/strain fields estimated from the single equivalent layer approach (conventional method) yielded different results from those predicted by the anatomically more accurate apparoach (i.e., layerwise). In addition, in a disc under a compressive force (symmetric loading), asymmetric stress-strain distributions were computed when using a layerwise simulation. Although all developed and selected published AF models predicted gross compression-displacement responses of the whole disc within the range of measured data, some showed excessively stiff or compliant responses under tissue-level uniaxial/biaxial tests. This study emphasizes, when constructing and validating constitutive models of AF, the importance of the proper simulation of individual lamellae as distinct layers, and testing parameters (sample geometric dimensions/loading/boundary conditions).

Computational Challenges in Tissue Engineering for the Spine

This paper deals with a brief review of the recent developments in computational modelling applied to innovative treatments of spine diseases. Additionally, it provides a perspective on the research directions expected for the forthcoming years. The spine is composed of distinct and complex tissues that require specific modelling approaches. With the advent of additive manufacturing and increasing computational power, patient-specific treatments have moved from being a research trend to a reality in clinical practice, but there are many issues to be addressed before such approaches become universal. Here, it is identified that the major setback resides in validation of these computational techniques prior to approval by regulatory agencies. Nevertheless, there are very promising indicators in terms of optimised scaffold modelling for both disc arthroplasty and vertebroplasty, powered by a decisive contribution from imaging methods.

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.

Interlamellar matrix governs human annulus fibrosus multiaxial behavior

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

Rapid increase in intradiscal pressure in porcine cervical spine units negatively impacts annulus fibrosus strength

Vertebral endplate fracture is generally accepted to occur as a result of excessive, typically rapid, pressure developed within the intervertebral disc. Clinical evidence of disc disorders later in life following endplate fracture suggests that the disc is also impacted by the fracture event. There is an abundance of evidence to indicate that loss of disc pressure at the time of fracture alters annulus fibrosus loading which could impact disc health long-term. However, little is known regarding the isolated impact of the rapid disc pressurization during the fracture event. This study sought to determine if rapid disc pressurization impacts mechanical properties of the annulus. Intervertebral discs of porcine cervical functional spine units (FSU) were pressurized via rapid injection of hydraulic fluid (standard motor oil). FSUs were subsequently grouped according to whether the pressurization resulted in endplate fracture ('fracture'; n = 9) or not ('non-fracture'; n = 13). Annular tensile properties and lamellar adhesion strength were quantified for both fractured and non-fractured samples. A third group was also examined in which FSUs did not undergo the pressurization protocol ('control'; n = 10). Annular tensile properties were not found to be statistically different between the three conditions (p > 0.05). Lamellar adhesion strength, in contrast, was found to be, on average, 53% lower in specimens which were pressurized compared to control, non-pressurized samples (p = 0.03). Further, no statistical differences in adhesion strength was found between fractured and non-fractured pressurized samples (p = 0.990) indicating that rapid pressurization , and not the presence of an endplate fracture, impacted the mechanical properties of the annulus. Decreased lamellar adhesion strength increases the potential of delamination, which can lead to disc disorders.

Osmo-inelastic response of the intervertebral disc annulus fibrosus tissue

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

Effect of aggrecan degradation on the nanomechanics of hyaluronan in extra-fibrillar matrix of annulus fibrosus: A molecular dynamics investigation

Intervertebral Disc (IVD) Degeneration is one of the primary causes of low back pain among the adult population - the most significant cause being the degradation of aggrecan present in the extra-fibrillar matrix (EFM). Aggrecan degradation is closely associated with loss of water content leading to an alteration in the mechanical behaviour of the IVD. The loss in water content has a significant impact on the chemo-mechanical interplay of IVD biochemical constituents at the fundamental level. This work presents a mechanistic understanding of the effect of hydration, closely associated with aggrecan degradation, on the nanoscale mechanical behaviour of the hyaluronan present in the EFM of the Annulus Fibrosus. For this purpose, explicit three-dimensional molecular dynamics analyses of tensile and compressive tests are performed on a representative atomistic model of the hyaluronan present in the EFM. To account for the degradation of aggrecan, hydration levels are varied from 0 to 75% by weight of water. Analyses show that an increase in the hydration levels decreases the elastic modulus of hyaluronan in tension from ~4.6 GPa to ~2.1 GPa. On the other hand, the increase in hydration level increases the elastic moduli in axial compression from ~1.6 GPa in un-hydrated condition to ~6 GPa in 50% hydrated condition. But as the hydration levels increase to 75%, the elastic modulus reduces to ~3.5 GPa signifying a shift in load-bearing characteristic, from the solid hyaluronan component to the fluid component. Furthermore, analyses show a reduction in the intermolecular energy between hyaluronan and water, under axial tensile loading, indicating a nanoscale intermolecular debonding between hyaluronan and water molecules. This is attributed to the ability of hyaluronan to form stabilizing intra-molecular hydrogen bonds between adjacent residues. Compressive loading, on the other hand, causes intensive coiling of hyaluronan molecule, which traps more water through hydrogen bonding and aids in bearing compressive loads. Overall, study shows that hydration level has a strong influence on the atomistic level interactions between hyaluronan molecules and hyaluronan and water molecules in the EFM which influences the nanoscale mechanics of the Annulus Fibrosus.

Numerical simulation of lateral and transforaminal lumbar interbody fusion, two minimally invasive surgical approaches

The present study aims to compare spinal stability after two different minimally invasive techniques, the lateral lumbar interbody fusion (LLIF) and the transforaminal lumbar interbody fusion (TLIF) approaches. Two nonlinear three-dimensional finite element (FE) models of the L4-L5 functional spinal unit (FSU) were subjected to the loads that usually act on the lumbar spine. Findings show that the LLIF approach yields better results for torsion load case, due to the larger surface area of the implant. For extension, flexion and lateral bending loads, the TLIF approach presents smaller displacements probably due to the anterior placement of the cage and to the smaller damaged area of the annulus fibrosus.

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

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

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.

Slow depressurization following intradiscal injection leads to injectate leakage in a large animal model

Needle injection has been indicated as the most practical method of delivering therapeutic agents to the intervertebral disc due to the disc's largely avascular nature. As the disc is characterized by both high stiffness and low permeability, injection requires substantial pressure, which may not relax on practical time scales. Additionally, needle puncture results in a localized disruption to the annulus fibrosus that can provide a leakage pathway for pressurized injectate. We hypothesized that intradiscal injection would result in slow relaxation of injectate pressure, followed by leakage upon needle retraction. This hypothesis was tested via controlled injection of fluorescently labeled saline into bovine caudal discs via a 21 gauge needle. Injections were performed with 10% of total disc volume injected at 3%/s followed by a 4-minute dwell. An analytical poroelastic model was calibrated to the experimental data and used to estimate injectate delivery with time. Experimental results confirmed both pressurization (with a peak of 199 ± 45 kPa) and slow recovery (final pressure of 81 ± 23 kPa). Injectate leakage through the needle puncture was verified following needle retraction in all samples. Histological sections of the discs displayed a clear defect at each disc's injection site with strong fluorescent labeling indicating a leakage pathway. The modeling results suggest that less than one-fourth of the injected volume was absorbed by the tissue in 4 minutes. Taken together these results suggest that needle injection is a feasible, albeit inefficient method for delivery of therapeutic agents into the intervertebral disc. Particular care should be taken to aspirate un-absorbed injectate prior to needle retraction to prevent leakage and exposure of surrounding tissues.

Relative Nucleus Pulposus Area and Position Alters Disc Joint Mechanics

Aging and degeneration of the intervertebral disc are noted by changes in tissue composition and geometry, including a decrease in nucleus pulposus (NP) area. The NP centroid is positioned slightly posterior of the disc's centroid, but the effect of NP size and location on disc joint mechanics is not well understood. We evaluated the effect of NP size and centroid location on disc joint mechanics under dual-loading modalities (i.e., compression in combination with axial rotation or bending). A finite element model was developed to vary the relative NP area (NP:Disc area ratio range = 0.21 - 0.60). We also evaluated the effect of NP position by shifting the NP centroid anteriorly and posteriorly. Our results showed that compressive stiffness and average first principal strains increased with NP size. Under axial compression, stresses are distributed from the NP to the annulus, and stresses were redistributed towards the NP with axial rotation. Moreover, peak stresses were greater for discs with a smaller NP area. NP centroid location had a greater impact on intradiscal pressure during flexion and extension, where peak pressures in the posterior annulus under extension was greater for discs with a more posteriorly situated NP. In conclusion, the findings from this study highlight the importance of closely mimicking NP size and location in computational models that aim to understand stress/strain distribution during complex loading and for developing repair strategies that aim to recapitulate the mechanical behavior of healthy discs.

Mapping of Intervertebral Disk Annulus Fibrosus Compressive Properties Is Sensitive to Specimen Boundary Conditions

Predicting the mechanical behavior of the intervertebral disk (IVD) in health and in disease requires accurate spatial mapping of its compressive mechanical properties. Previous studies confirmed that residual strains in the annulus fibrosus (AF) of the IVD, which result from nonuniform extracellular matrix deposition in response to in vivo loads, vary by anatomical regions (anterior, posterior, and lateral) and zones (inner, middle, and outer). We hypothesized that as the AF is composed of a nonlinear, anisotropic, viscoelastic material, the state of residual strain in the transverse plane would influence the apparent values of axial compressive properties. To test this hypothesis, axial creep indentation tests were performed, using a 1.6 mm spherical probe, at nine different anatomical locations on bovine caudal AFs in both the intact (residual strain present) and strain relieved states. The results showed a shift toward increased spatial homogeneity in all measured parameters, particularly instantaneous strain. This shift was not observed in control AFs, which were tested twice in the intact state. Our results confirm that time-dependent axial compressive properties of the AF are sensitive to the state of residual strain in the transverse plane, to a degree that is likely to affect whole disk behavior.

Validation of an In Vivo Medical Image-Based Young Human Lumbar Spine Finite Element Model

Mathematical models of the human spine can be used to investigate spinal biomechanics without the difficulties, limitations, and ethical concerns associated with physical experimentation. Validation of such models is necessary to ensure that the modeled system behavior accurately represents the physics of the actual system. The goal of this work was to validate a medical image-based nonlinear lumbosacral spine finite element model of a healthy 20-yr-old female subject under physiological moments. Range of motion (ROM), facet joint forces (FJF), and intradiscal pressure (IDP) were compared with experimental values and validated finite element models from the literature. The finite element model presented in this work was in good agreement with published experimental studies and finite element models under pure moments. For applied moments of 7.5 N·m, the ROM in flexion–extension, axial rotation, and lateral bending were 39 deg, 16 deg, and 28 deg, respectively. Excellent agreement was observed between the finite element model and experimental data for IDP under pure compressive loading. The predicted FJFs were lower than those of the experimental results and validated finite element models for extension and torsion, likely due to the nondegenerate properties chosen for the intervertebral disks and morphology of the young female spine. This work is the first to validate a computational lumbar spine model of a young female subject. This model will serve as a valuable tool for predicting orthopedic spinal injuries, studying the effect of intervertebral disk replacements using advanced biomaterials, and investigating soft tissue degeneration.

Materials for the Spine: Anatomy, Problems, and Solutions

Disc degeneration affects 12% to 35% of a given population, based on genetics, age, gender, and other environmental factors, and usually occurs in the lumbar spine due to heavier loads and more strenuous motions. Degeneration of the extracellular matrix (ECM) within reduces mechanical integrity, shock absorption, and swelling capabilities of the intervertebral disc. When severe enough, the disc can bulge and eventually herniate, leading to pressure build up on the spinal cord. This can cause immense lower back pain in individuals, leading to total medical costs exceeding $100 billion. Current treatment options include both invasive and noninvasive methods, with spinal fusion surgery and total disc replacement (TDR) being the most common invasive procedures. Although these treatments cause pain relief for the majority of patients, multiple challenges arise for each. Therefore, newer tissue engineering methods are being researched to solve the ever-growing problem. This review spans the anatomy of the spine, with an emphasis on the functions and biological aspects of the intervertebral discs, as well as the problems, associated solutions, and future research in the field.

Swelling of fiber-reinforced soft tissues is affected by fiber orientation, fiber stiffness, and lamella structure

Native and engineered fiber-reinforced tissues are composites comprised of stiff collagen fibers embedded within an extrafibrillar matrix that is capable of swelling by absorbing water molecules. Tissue swelling is important for understanding stress distributions between collagen fibers and extrafibrillar matrix, as well as for understanding mechanisms of tissue failure. The swelling behavior of fiber-reinforced tissues in the musculoskeletal system has been largely attributed to the glycosaminoglycan content. Recent work demonstrated anisotropy in the swelling response of the annulus fibrosus in the intervertebral disc. It is well known that collagen fiber orientation affects elastic behavior, but the effect of collagen fiber network on tissue swelling behavior is not well understood. In this study, we developed three series of models to evaluate the effect of collagen fiber orientation, fiber network architecture (i.e., single or multi-fiber families within a layer), and fiber stiffness on bulk tissue swelling, which was simulated by describing the extrafibrillar matrix as a triphasic material, as proposed by Lai et al. Model results were within one standard deviation of reported mean values for changes in tissue volume, width, and thickness under free swelling conditions. The predicted swelling response of single-fiber family structures was highly dependent on fiber orientation and the number of lamellae in the bulk tissue. Moreover, matrix swelling resulted in tissue to twist, which reduced fiber deformations, demonstrating a balance between fiber deformation and matrix swelling. Large changes in fiber stiffness (20 × increase) had a relatively small effect on tissue swelling (~ 2% decrease in swelling). In conclusion, fiber angle, fiber architecture (defined as single- versus multiple fiber families in a layer), and the number of layers in a single fiber family structure directly affected tissue swelling behavior, including fiber stretch, fiber reorientation, and tissue deformation. These findings support the need to develop computational models that closely mimic the native architecture in order to understand mechanisms of stress distributions and tissue failure.

A magnetic resonance imaging framework for quantifying intervertebral disc deformation in vivo: Reliability and application to diurnal variations in lumbar disc shape

Low back pain is a significant socioeconomic burden in the United States and lumbar intervertebral disc degeneration is frequently implicated as a cause. The discs play an important mechanical role in the spine, yet the relationship between disc function and back pain is poorly defined. The objective of this work was to develop a technique using magnetic resonance imaging (MRI) and three-dimensional modeling to measure in vivo disc deformations. Using this method, we found that disc geometry was measurable with precision less than the in-plane dimensions of a voxel (≈100 µm, 10% of the MRI pixel size). Furthermore, there was excellent agreement between mean disc height, disc perimeter, disc volume and regional disc height measurements for multiple trials from an individual rater (standard deviation <3.1% across all measurements) and between mean height, perimeter, and volume measurements made by two independent raters (error <1.5% across all measurements). We then used this measurement system to track diurnal deformations in the L5-S1 disc in a young, healthy population (n = 8; age 24.1 ± 3.3 yrs; 2 M/6F). We measured decreases in the mean disc height (-8%) and volume (-9%) with no changes in perimeter over an eight-hour workday. We found that the largest height losses occurred in the posterior (-13%) and posterior-lateral (-14%) regions adjacent to the outer annulus fibrosus. Diurnal annulus fibrosus (AF) strains induced by posterior and posterior-lateral height loss may increase the risk for posterior disc herniation or posterior AF tears. These preliminary findings lay a foundation for determining how deviations from normal deformations may contribute to back pain.

The effects of needle damage on annulus fibrosus micromechanics

Needle puncture of the intervertebral disc can initiate a mechanical and biochemical cascade leading to disc degeneration. Puncture's mechanical effects have been shown near the puncture site, mechanical effects should be observed far, relative to needle size, from the puncture site, given the disc-wide damage induced by the stab. The aim of this work was to quantify these far-field effects, and to observe the local structural damage provoked by the needle. Strips of cow tail annulus fibrosus underwent two consecutive mechanical loadings to 5% tensile strain; fifteen samples were punctured in a radial direction with a randomly assigned needle between the two loadings (needle gauges between 19 and 23). Ten samples (control group) were not punctured. During loading, the tissue strains were imaged using second harmonic generation microscopy in a <600×800µm region about 4.4mm from the puncture site. After mechanical testing, the puncture site was imaged in 3D. Puncture had no significant effect on annulus elastic modulus. Imaging showed a modest change in the shearing between fibre bundles however, the linear strain between bundles, intra-bundle shear and linear strain were not significantly affected. At the puncture site, detached lumps of tissue were present. These results suggest that the mechanical effects observed in intact discs are due to the depressurization of the disc, rather than the local damage to the annulus. Needle profiles could be designed, aiming at separating fibre bundles rather than cutting through them, to avoid leaving dying tissue behind.

STATEMENT OF SIGNIFICANCE

Needle puncture of the intervertebral disc can initiate a mechanical and biochemical cascade leading to disc degeneration, but the link between the local damage of the puncture and the disc-wide effects is not well understood. This work aimed at determining the micro-mechanical effects of the puncture far from its site, and to observe the damage induced by the puncture with high resolution imaging. Results show that the puncture had modest effect far from the puncture, but lumps of tissue were left by the needle, detached from the disc; these could cause further damage through friction and inflammation of the surrounding tissues. This suggests that the cascade leading to degeneration is probably driven by a biochemical response rather than disc-wide mechanical effects.

A Novel Method for Repeatable Failure Testing of Annulus Fibrosus

Tears in the annulus fibrosus (AF) of the intervertebral disk can result in disk herniation and progressive degeneration. Understanding AF failure mechanics is important as research moves toward developing biological repair strategies for herniated disks. Unfortunately, failure mechanics of fiber-reinforced tissues, particularly tissues with fibers oriented off-axis from the applied load, is not well understood, partly due to the high variability in reported mechanical properties and a lack of standard techniques ensuring repeatable failure behavior. Therefore, the objective of this study was to investigate the effectiveness of midlength (ML) notch geometries in producing repeatable and consistent tissue failure within the gauge region of AF mechanical test specimens. Finite element models (FEMs) representing several notch geometries were created to predict the location of bulk tissue failure using a local strain-based criterion. FEM results were validated by experimentally testing a subset of the modeled specimen geometries. Mechanical testing data agreed with model predictions (∼90% agreement), validating the model's predictive power. Two of the modified dog-bone geometries (“half” and “quarter”) effectively ensured tissue failure at the ML for specimens oriented along the circumferential-radial and circumferential-axial directions. The variance of measured mechanical properties was significantly lower for notched samples that failed at the ML, suggesting that ML notch geometries result in more consistent and reliable data. In addition, the approach developed in this study provides a framework for evaluating failure properties of other fiber-reinforced tissues, such as tendons and meniscus.

Angle-ply biomaterial scaffold for annulus fibrosus repair replicates native tissue mechanical properties, restores spinal kinematics, and supports cell viability

Annulus fibrosus (AF) damage commonly occurs due to intervertebral disc (IVD) degeneration/herniation. The dynamic mechanical role of the AF is essential for proper IVD function and thus it is imperative that biomaterials developed to repair the AF withstand the mechanical rigors of the native tissue. Furthermore, these biomaterials must resist accelerated degradation within the proteolytic environment of degenerate IVDs while supporting integration with host tissue. We have previously reported a novel approach for developing collagen-based, multi-laminate AF repair patches (AFRPs) that mimic the angle-ply architecture and basic tensile properties of the human AF. Herein, we further evaluate AFRPs for their: tensile fatigue and impact burst strength, IVD attachment strength, and contribution to functional spinal unit (FSU) kinematics following IVD repair. Additionally, AFRP resistance to collagenase degradation and cytocompatibility were assessed following chemical crosslinking. In summary, AFRPs demonstrated enhanced durability at high applied stress amplitudes compared to human AF and withstood radially-directed biaxial stresses commonly borne by the native tissue prior to failure/detachment from IVDs. Moreover, FSUs repaired with AFRPs and nucleus pulposus (NP) surrogates had their axial kinematic parameters restored to intact levels. Finally, carbodiimide crosslinked AFRPs resisted accelerated collagenase digestion without detrimentally effecting AFRP tensile properties or cytocompatibility. Taken together, AFRPs demonstrate the mechanical robustness and enzymatic stability required for implantation into the damaged/degenerate IVD while supporting AF cell infiltration and viability.

STATEMENT OF SIGNIFICANCE

The quality of life for millions of individuals globally is detrimentally impacted by IVD degeneration and herniation. These pathologies often result in the structural demise of IVD tissue, particularly the annulus fibrosus (AF). Biomaterials developed for AF repair have yet to demonstrate the mechanical strength and durability required for utilization in the spine. Herein, we demonstrate the development of an angle-ply AF repair patch (AFRP) that can resist the application of physiologically relevant stresses without failure and which contributes to the restoration of functional spinal unit axial kinematics following repair. Furthermore, we show that this biomaterial can resist accelerated degradation in a simulated degenerate environment and supports AF cell viability.

Bovine and degenerated human annulus fibrosus: a microstructural and micromechanical comparison

The complex structure of the annulus fibrosus is strongly related to its mechanical properties. Recent work showed that it is possible to observe the relative movement of fibre bundles in loaded cow tail annulus; the aim of this work was to describe and quantify annulus fibrosus micromechanics in degenerated human disc, and compare it with cow tail annulus, an animal model often used in the literature. Second harmonic generation was used to image the collagen matrix in twenty strips of annulus fibrosus harvested from intervertebral disc of seven patients undergoing surgery. Samples were loaded to 6% tensile strain in 1% steps. Elastic modulus was calculated from loading curves, and micromechanical strains were calculated from the images using custom software. The same protocol was applied to twenty strips of annulus harvested from cow tail discs. Significant morphological differences were found between human and cow tail samples, the most striking being the lack of collagen fibre crimp in the former. Fibres were also observed bending and running from one lamella to the other, forming a strong flexible interface. Interdigitation of fibre bundles was also present at this interface. Quantitative results show complex patterns of inter-bundle and inter-lamellar behaviour, with inter-bundle sliding being the main strain mechanism. Elastic modulus was similar between species, and it was not affected by the degree of degeneration. This work gives an insight into the complex structure and mechanical function of the annulus fibrosus, which should be accounted for in disc numerical modelling.

The fabrication and characterization of a multi-laminate, angle-ply collagen patch for annulus fibrosus repair

One major limitation of intervertebral disc (IVD) repair is that no ideal biomaterial has been developed that effectively mimics the angle-ply collagen architecture and mechanical properties of the native annulus fibrosus (AF). Furthermore, it would be beneficial to devise a simple, scalable process by which to manufacture a biomimetic biomaterial that could function as a mechanical repair patch to be secured over a large defect in the outer AF that will support AF tissue regeneration. Such a biomaterial would: (1) enable the employment of early-stage interventional strategies to treat IVD degeneration (i.e. nucleus pulposus arthroplasty); (2) prevent IVD re-herniation in patients with large AF defects; and (3) serve as a platform to develop full-thickness AF and whole IVD tissue engineering strategies. Due to the innate collagen fibre alignment and mechanical strength of pericardium, a procedure was developed to assemble multi-laminate angle-ply AF patches derived from decellularized pericardial tissue. Patches were subsequently assessed histologically to confirm angle-ply microarchitecture, and mechanically assessed for biaxial burst strength and tensile properties. Additionally, patch cytocompatibility was evaluated following seeding with bovine AF cells. This study demonstrated the effective removal of porcine cell remnants from the pericardium, and the ability to reliably produce multi-laminate patches with angle-ply architecture using a simple assembly technique. Resultant patches demonstrated their inherent ability to resist biaxial burst pressures reminiscent of intradiscal pressures commonly borne by the AF, and exhibited tensile strength and modulus values reported for native human AF. Furthermore, the biomaterial supported AF cell viability, infiltration and proliferation. Copyright © 2016 John Wiley & Sons, Ltd.

Design Requirements for Annulus Fibrosus Repair: Review of Forces, Displacements, and Material Properties of the Intervertebral Disk and a Summary of Candidate Hydrogels for Repair

There is currently a lack of clinically available solutions to restore functionality to the intervertebral disk (IVD) following herniation injury to the annulus fibrosus (AF). Microdiscectomy is a commonly performed surgical procedure to alleviate pain caused by herniation; however, AF defects remain and can lead to accelerated degeneration and painful conditions. Currently available AF closure techniques do not restore mechanical functionality or promote tissue regeneration, and have risk of reherniation. This review determined quantitative design requirements for AF repair materials and summarized currently available hydrogels capable of meeting these design requirements by using a series of systematic PubMed database searches to yield 1500+ papers that were screened and analyzed for relevance to human lumbar in vivo measurements, motion segment behaviors, and tissue level properties. We propose a testing paradigm involving screening tests as well as more involved in situ and in vivo validation tests to efficiently identify promising biomaterials for AF repair. We suggest that successful materials must have high adhesion strength (∼0.2 MPa), match as many AF material properties as possible (e.g., approximately 1 MPa, 0. 3 MPa, and 30 MPa for compressive, shear, and tensile moduli, respectively), and have high tensile failure strain (∼65%) to advance to in situ and in vivo validation tests. While many biomaterials exist for AF repair, few undergo extensive mechanical characterization. A few hydrogels show promise for AF repair since they can match at least one material property of the AF while also adhering to AF tissue and are capable of easy implantation during surgical procedures to warrant additional optimization and validation.

Moderately degenerated lumbar motion segments: Are they truly unstable?

The two main load bearing tissues of the intervertebral disc are the nucleus pulposus and the annulus fibrosus. Both tissues are composed of the same basic components, but differ in their organization and relative amounts. With degeneration, the clear distinction between the two tissues disappears. The changes in biochemical content lead to changes in mechanical behaviour of the intervertebral disc. The aim of the current study was to investigate if well-documented moderate degeneration at the biochemical and fibre structure level leads to instability of the lumbar spine. By taking into account biochemical and ultrastructural changes to the extracellular matrix of degenerating discs, a set of constitutive material parameters were determined that described the individual tissue behaviour. These tissue biomechanical models were then used to simulate dynamic behaviour of the degenerated spinal motion segment, which showed instability in axial rotation, while a stabilizing effect in the other two principle bending directions. When a shear load was applied to the degenerated spinal motion segment, no sign of instability was found. This study found that reported changes to the nucleus pulposus and annulus fibrosus matrix during moderate degeneration lead to a more stable spinal motion segment and that such biomechanical considerations should be incorporated into the general pathophysiological understanding of disc degeneration and how its progress could affect low back pain and its treatments thereof.