Relative Nucleus Pulposus Area and Position Alters Disc Joint Mechanics

Morphology and Composition of Lumbar Intervertebral Discs: Comparative Analyses of Manual Measurement and Computer-Assisted Algorithms

BACKGROUND

The morphology and internal composition, particularly the nucleus-to-cross sectional area (NP-to-CSA) ratio of the lumbar intervertebral discs (IVDs), is important information for finite element models (FEMs) of spinal loadings and biomechanical behaviors, and, yet, this has not been well investigated and reported.

METHODS

Anonymized MRI scans were retrieved from a previously established database, including a total of 400 lumbar IVDs from 123 subjects (58 F and 65 M). Measurements were conducted manually by a spine surgeon and using two computer-assisted segmentation algorithms, i.e., fuzzy C-means (FCM) and region growing (RG). The respective results were compared. The influence of gender and spinal level was also investigated.

RESULTS

Ratios derived from manual measurements and the two computer-assisted algorithms (FCM and RG) were 46%, 39%, and 38%, respectively. Ratios derived manually were significantly larger.

CONCLUSIONS

Computer-assisted methods provide reliable outcomes that are traditionally difficult for the manual measurement of internal composition. FEMs should consider the variability of NP-to-CSA ratios when studying the biomechanical behavior of the spine.

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

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

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.

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.

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).

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

Tissues in the intervertebral disc have a large capacity to absorb water, partially due to the high glycosaminoglycan (GAG) content, which decreases linearly from the nucleus pulposus (NP) in the center to the outer annulus. Our recent work showed that fiber network and GAG distribution contributes to development of residual stresses and strains that were compressive in the inner annulus to tensile in the outer annulus. GAG loss in the inner annulus, as observed with early to moderate degeneration, reduced swelling capacity and circumferential-direction stress by over 50%. However, our previous model was not capable of evaluating interactions between the NP and annulus fibrosus (AF) during swelling. In this study, we evaluated the effect of degeneration (GAG content or swelling capacity) on residual stress development throughout the disc. Simulations of moderate to severe degeneration showed a 40% decrease in NP swelling capacity, with a 25% decrease in AF and cartilaginous endplate swelling. Together, these changes in tissue swelling resulted in a decrease in NP pressure (healthy = 0.21 MPa; severe degeneration = 0.03 MPa) that was comparable to observations in human discs. There was a 60% decrease in circumferential-direction residual deformations with early degeneration. Radial-direction stretch switched from compressive to tensile with degeneration, which may increase the risk for tears or delamination. Degeneration had a significant impact on residual stress/stretch and fiber stretch in the posterior AF, which is important for understanding herniation risk. In conclusion, degenerative changes in disc geometry and intradiscal deformations was recreated by only altering NP and AF GAG composition. Since most computational models simulate degeneration by altering material stiffness, this work highlights the importance of directly simulating biochemical composition and distribution to study disc biomechanics with degeneration. STATEMENT OF SIGNIFICANCE: Tissues in the intervertebral disc have a large swelling capacity, due to its high glycosaminoglycan content. Our recent work demonstrated the importance of fiber network and glycosaminoglycan distribution residual stresses and strains development. In this study, we evaluated the effect of swelling on intradiscal deformations between the nucleus pulposus and annulus fibrosus. We also investigated the effect of degenerative glycosaminoglycan loss on swelling-based intradiscal deformations of the intact disc and its subcomponents. Decreases in nucleus glycosaminoglycan content resulted in morphological changes observed with degenerated discs and may help to explain mechanisms behind the increases in annular tears and mechanical dysfunction with degeneration.