Application of Finite Element Analysis for Investigation of Intervertebral Disc Degeneration: from Laboratory to Clinic

Abstract

Due to the ethical concern and inability to detect inner stress distributions of intervertebral disc (IVD), traditional methods for investigation of intervertebral disc degeneration (IVDD) have significant limitations. Many researchers have demonstrated that finite element analysis (FEA) is an effective tool for the research of IVDD. However, the specific application of FEA for investigation of IVDD has not been systematically elucidated before. In the present review, we summarize the current finite element models (FEM) used for the investigation of IVDD, including the poroelastic nonlinear FEM, diffusive-reactive theory model and cell-activity coupled mechano-electrochemical theory model. We further elaborate the use of FEA for the research of IVDD pathogenesis especially for nutrition and biomechanics associated etiology, and the biological, biomechanical and clinical influences of IVDD. In addition, the application of FEA for evaluation and exploration of various treatments for IVDD is also elucidated. We conclude that FEA is an excellent technique for research of IVDD, which could be used to explore the etiology, biology and biomechanics of IVDD. In the future, FEA may help us to achieve the goal of individualized precision therapy.

Keywords

finite element analysis / intervertebral disc degeneration / biomechanics / spine

References

[1] AdamsMA, DolanP. Spine biomechanics. J Biomech, 2005, 38(10): 1972-1983

[2] ChoiH, JohnsonZI, RisbudMV. Understanding nucleus pulposus cell phenotype: a prerequisite for stem cell based therapies to treat intervertebral disc degeneration. Curr Stem Cell Res Ther, 2015, 10(4): 307-316

[3] HoyD, MarchL, BrooksP, et al.. The global burden of low back pain: estimates from the Global Burden of Disease 2010 study. Ann Rheum Dis, 2014, 73(6): 968-974

[4] ChenS, LvX, HuB, et al.. RIPK1/RIPK3/MLKLmediated necroptosis contributes to compressioninduced rat nucleus pulposus cells death. Apoptosis, 2017, 22(5): 626-638

[5] AliniM, EisensteinSM, ItoK, et al.. Are animal models useful for studying human disc disorders/degeneration. Eur Spine J, 2008, 17(1): 2-19

[6] GantenbeinB, Illien-JüngerS, ChanS, et al.. Organ Culture Bioreactors-Platforms to Study Human Intervertebral Disc Degeneration and Regenerative Therapy. Curr Stem Cell Res Ther, 2015, 10(4): 339-352

[7] MoonSM, YoderJH, WrightAC, et al.. Evaluation of intervertebral disc cartilaginous endplate structure using magnetic resonance imaging. Eur Spine J, 2013, 22(8): 1820-1828

[8] MieleVJ, PanjabiMM, BenzelEC. Anatomy and biomechanics of the spinal column and cord. Handb Clin Neurol, 2012, 109: 31-43

[9] ChadderdonRC, ShimerAL, GilbertsonLG, et al.. Advances in gene therapy for intervertebral disc degeneration. Spine J, 2004, 4(6Suppl): 341s-347s

[10] FaganMJ, JulianS, MohsenAM. Finite element analysis in spine research. Proc Inst Mech Eng H, 2002, 216(5): 281-298

[11] BrekelmansWA, PoortHW, SlooffTJ. A new method to analyse the mechanical behaviour of skeletal parts. Acta Orthop Scand, 1972, 43(5): 301-317

[12] NikkhooM, HsuYC, HaghpanahiM, et al.. A metamodel analysis of a finite element simulation for defining poroelastic properties of intervertebral discs. Proc Inst Mech Eng H, 2013, 227(6): 672-682

[13] BelytschkoT, KulakRF, SchultzAB, et al.. Finite element stress analysis of an intervertebral disc. J Biomec, 1974, 7(3): 277-285

[14] MalandrinoA, NoaillyJ, LacroixD. The effect of sustained compression on oxygen metabolic transport in the intervertebral disc decreases with degenerative changes. PLoS Comput Biol, 2011, 7(8): e1002112

[15] SongC, LiXF, LiuZD, et al.. Biomechanical assessment of a novel L4/5 level interspinous implant using three dimensional finite element analysis. Eur Rev Med Pharmacol Sci, 2014, 18(1): 86-94

[16] TangS. Does TLIF aggravate adjacent segmental degeneration more adversely than ALIF? A finite element study. Turk Neurosurg, 2012, 22(3): 324-328

[17] TravascioF, ElmasryS, AsfourS. Modeling the role of IGF-1 on extracellular matrix biosynthesis and cellularity in intervertebral disc. J Biomech, 2014, 47(10): 2269-2276

[18] IbarzE, HerreraA, MasY, et al.. Development and kinematic verification of a finite element model for the lumbar spine: application to disc degeneration. Acta Bioeng Biomech, 2013, 2013: 705185

[19] HussainM, NatarajanRN, AnHS, et al.. Progressive disc degeneration at C5-C6 segment affects the mechanics between disc heights and posterior facets above and below the degenerated segment: A flexion-extension investigation using a poroelastic C3-T1 finite element model. Med Eng Phys, 2012, 34(5): 552-558

[20] ZhuQ, GaoX, LeveneHB, et al.. Influences of Nutrition Supply and Pathways on the Degenerative Patterns in Human Intervertebral Disc. Spine, 2016, 41(7): 568-576

[21] NatarajanRN, WilliamsJR, AnderssonGB. Modeling changes in intervertebral disc mechanics with degeneration. J Bone Joint Surg Am, 2006, 88(Suppl2): 36-40

[22] MasseyCJ, van DonkelaarCC, VresilovicE, et al.. Effects of aging and degeneration on the human intervertebral disc during the diurnal cycle: a finite element study. J Orthop Res, 2012, 30(1): 122-128

[23] QasimM, NatarajanRN, AnHS, et al.. Damage accumulation location under cyclic loading in the lumbar disc shifts from inner annulus lamellae to peripheral annulus with increasing disc degeneration. J Biomech, 2014, 47(1): 24-31

[24] von ForellGA, StephensTK, SamartzisD, et al.. Low Back Pain: A Biomechanical Rationale Based on "Patterns" of Disc Degeneration. Spine, 2015, 40(15): 1165-1172

[25] WuY, WangY, WuJ, et al.. Study of Double-level Degeneration of Lower Lumbar Spines by Finite Element Model. World Neurosurg, 2016, 86: 294-299

[26] TangS, RebholzBJ. Does anterior lumbar interbody fusion promote adjacent degeneration in degenerative disc disease? A finite element study. J Orthop Sci, 2011, 16(2): 221-228

[27] HussainM, NatarajanRN, AnHS, et al.. Reduction in segmental flexibility because of disc degeneration is accompanied by higher changes in facet loads than changes in disc pressure: a poroelastic C5-C6 finite element investigation. Spine J, 2010, 10(12): 1069-1077

[28] EllingsonAM, ShawMN, GiambiniH, et al.. Comparative role of disc degeneration and ligament failure on functional mechanics of the lumbar spine. Comput Methods Biomech Biomed Engin, 2016, 19(9): 1009-1018

[29] SchmidtH, GalbuseraF, RohlmannA, et al.. Effect of multilevel lumbar disc arthroplasty on spine kinematics and facet joint loads in flexion and extension: a finite element analysis. Eur Spine J, 2012, 21(Suppl5): S663-S674

[30] LittleJP, AdamCJ, EvansJH, et al.. Nonlinear finite element analysis of anular lesions in the L4/5 intervertebral disc. J Biomech, 2007, 40(12): 2744-2751

[31] HommingaJ, AquariusR, BulsinkVE, et al.. Can vertebral density changes be explained by intervertebral disc degeneration. Med Eng Phys, 2012, 34(4): 453-458

[32] ElmasryS, AsfourS, de Rivero VaccariJP, et al.. Effects of Tobacco Smoking on the Degeneration of the Intervertebral Disc: A Finite Element Study. PLoS One, 2015, 10(8): e0136137

[33] HanI, RopperAE, KonyaD, et al.. Biological Approaches to Treating Intervertebral Disk Degeneration: Devising Stem Cell Therapies. Cell Transplantation, 2015, 24(11): 2197-2208

[34] AsfourS, TravascioF, ElmasryS, et al.. A computational analysis on the implications of age-related changes in the expression of cellular signals on the role of IGF-1 in intervertebral disc homeostasis. J Biomech, 2015, 48(2): 332-339

[35] ZhuQ, GaoX, GuW. Temporal changes of mechanical signals and extracellular composition in human intervertebral disc during degenerative progression. J Biomech, 2014, 47(15): 3734-3743

[36] ZhuQ, JacksonAR, GuWY. Cell viability in intervertebral disc under various nutritional and dynamic loading conditions: 3d finite element analysis. J Biomech, 2012, 45(16): 2769-2777

[37] JacksonAR, HuangCY, BrownMD, et al.. 3D finite element analysis of nutrient distributions and cell viability in the intervertebral disc: effects of deformation and degeneration. J Biomech Eng, 2011, 133(9): 091006

[38] ZhuQ, GaoX, BrownMD, et al.. Simulation of water content distributions in degenerated human intervertebral discs. J Orthop Res, 2017, 35(1): 147-153

[39] ChagnonA, AubinCE, VillemureI. Biomechanical influence of disk properties on the load transfer of healthy and degenerated disks using a poroelastic finite element model. J Biomech Eng, 2010, 132(11): 111006

[40] GalbuseraF, SchmidtH, Neidlinger-WilkeC, et al.. The mechanical response of the lumbar spine to different combinations of disc degenerative changes investigated using randomized poroelastic finite element models. Eur Spine J, 2011, 20(4): 563-571

[41] HuangCY, TravascioF, GuWY. Quantitative analysis of exogenous IGF-1 administration of intervertebral disc through intradiscal injection. J Biomech, 2012, 45(7): 1149-1155

[42] HussainM, NatarajanRN, AnHS, et al.. Motion changes in adjacent segments due to moderate and severe degeneration in C5-C6 disc: a poroelastic C3-T1 finite element model study. Spine, 2010, 35(9): 939-947

[43] GuW, ZhuQ, GaoX, et al.. Simulation of the progression of intervertebral disc degeneration due to decreased nutritional supply. Spine, 2014, 39(24): E1411-E1417

[44] SmithLJ, NerurkarNL, ChoiKS, et al.. Degeneration and regeneration of the intervertebral disc: lessons from development. Dis Model Mech, 2011, 4(1): 31-41

[45] HuangCY, GuWY. Effects of mechanical compression on metabolism and distribution of oxygen and lactate in intervertebral disc. J Biomech, 2008, 41(6): 1184-1196

[46] SoukaneDM, Shirazi-AdlA, UrbanJP. Computation of coupled diffusion of oxygen, glucose and lactic acid in an intervertebral disc. J Biomech, 2007, 40(12): 2645-2654

[47] WuY, CisewskiS, SachsBL, et al.. Effect of cartilage endplate on cell based disc regeneration: a finite element analysis. Mol Cell Biomech, 2013, 10(2): 159-182

[48] JacksonAR, HuangCY, GuWY. Effect of endplate calcification and mechanical deformation on the distribution of glucose in intervertebral disc: a 3D finite element study. Comput Methods Biomech Biomed Engin, 2011, 14(2): 195-204

[49] Shirazi-AdlA, TaheriM, UrbanJP. Analysis of cell viability in intervertebral disc: Effect of endplate permeability on cell population. J Biomech, 2010, 43(7): 1330-1336

[50] MalandrinoA, NoaillyJ, LacroixD. Numerical exploration of the combined effect of nutrient supply, tissue condition and deformation in the intervertebral disc. J Biomech, 2014, 47(6): 1520-1525

[51] NachemsonA, LewinT, MaroudasA, et al.. In vitro diffusion of dye through the end-plates and the annulus fibrosus of human lumbar inter-vertebral discs. Acta Orthop Scand, 1970, 41(6): 589-607

[52] DeLuccaJF, CortesDH, JacobsNT, et al.. Human cartilage endplate permeability varies with degeneration and intervertebral disc site. J Biomech, 2016, 49(4): 550-557

[53] AyturkUM, GadomskiB, SchuldtD, et al.. Modeling degenerative disk disease in the lumbar spine: a combined experimental, constitutive, and computational approach. J Biomech Eng, 2012, 134(10): 101003

[54] GalbuseraF, SchmidtH, Neidlinger-WilkeC, et al.. The effect of degenerative morphological changes of the intervertebral disc on the lumbar spine biomechanics: a poroelastic finite element investigation. Comput Methods Biomech Biomed Engin, 2011, 14(8): 729-739

[55] MaquerG, SchwiedrzikJ, HuberG, et al.. Compressive strength of elderly vertebrae is reduced by disc degeneration and additional flexion. J Mech Behav Biomed Mater, 2015, 42: 54-66

[56] HussainM, NatarajanRN, AnHS, et al.. Patterns of height changes in anterior and posterior cervical disc regions affects the contact loading at posterior facets during moderate and severe disc degeneration: a poroelastic C5-C6 finite element model study. Spine, 2010, 35(18): E873-E881

[57] KimYE, GoelVK, WeinsteinJN, et al.. Effect of disc degeneration at one level on the adjacent level in axial mode. Spine, 1991, 16(3): 331-335

[58] RuberteLM, NatarajanRN, AnderssonGB. Influence of single-level lumbar degenerative disc disease on the behavior of the adjacent segments—a finite element model study. J Biomech, 2009, 42(3): 341-348

[59] LuYM, HuttonWC, GharpurayVM. Do bending, twisting, and diurnal fluid changes in the disc affect the propensity to prolapse? A viscoelastic finite element model. Spine, 1996, 21(22): 2570-2579

[60] von ForellGA, NelsonTG, SamartzisD, et al.. Changes in vertebral strain energy correlate with increased presence of Schmorl’s nodes in multi-level lumbar disk degeneration. J Biomech Eng, 2014, 136(6): 061002

[61] MaidhofR, AlipuiDO, RafiuddinA, et al.. Emerging trends in biological therapy for intervertebral disc degeneration. Discov Med, 2012, 14(79): 401-411

[62] ZagraA, ScaramuzzoL, GalbuseraF, et al.. Biomechanical and clinical study of single posterior oblique cage POLIF in the treatment of degenerative diseases of the lumbar spine. Eur Spine J, 2015, 24(Suppl7): 924-930

[63] ErbulutDU, KiapourA, OktenogluT, et al.. A computational biomechanical investigation of posterior dynamic instrumentation: combination of dynamic rod and hinged (dynamic) screw. J Biomech Eng, 2014, 136(5): 051007

[64] ChienCY, KuoYJ, LinSC, et al.. Kinematic and mechanical comparisons of lumbar hybrid fixation using Dynesys and Cosmic systems. Spine (Phila Pa 1976), 2014, 39(15): E878-E884

[65] CegoninoJ, Calvo-EcheniqueA, Perez-del PalomarA. Influence of different fusion techniques in lumbar spine over the adjacent segments: A 3D finite element study. J Orthop Res, 2015, 33(7): 993-1000

[66] FaizanA, GoelVK, BiyaniA, et al.. Adjacent level effects of bi level disc replacement, bi level fusion and disc replacement plus fusion in cervical spine—a finite element based study. Clin Biomech, 2012, 27(3): 226-233

[67] ChungTT, HuengDY, LinSC. Hybrid Strategy of Two-Level Cervical Artificial Disc and Intervertebral Cage: Biomechanical Effects on Tissues and Implants. Medicine, 2015, 94(47): e2048

[68] ZhuQ, GaoX, TempleHT, et al.. Simulation of biological therapies for degenerated intervertebral discs. J Orthop Res, 2016, 34(4): 699-708

[69] GanJC, DucheyneP, VresilovicE, et al.. Bioactive glass serves as a substrate for maintenance of phenotype of nucleus pulposus cells of the intervertebral disc. J Biomed Mater Res, 2000, 51(4): 596-604

[70] YaoJ, TurteltaubSR, DucheyneP. A three-dimensional nonlinear finite element analysis of the mechanical behavior of tissue engineered intervertebral discs under complex loads. Biomaterials, 2006, 27(3): 377-387

[71] StrangeDG, FisherST, BoughtonPC, et al.. Restoration of compressive loading properties of lumbar discs with a nucleus implant-a finite element analysis study. Spine J, 2010, 10(7): 602-609