Comparison of patient-specific computational models vs. clinical follow-up, for adjacent segment disc degeneration and bone remodelling after spinal fusion

Automatic generation of subject-specific finite element models of the spine from magnetic resonance images

The generation of subject-specific finite element models of the spine is generally a time-consuming process based on computed tomography (CT) images, where scanning exposes subjects to harmful radiation. In this study, a method is presented for the automatic generation of spine finite element models using images from a single magnetic resonance (MR) sequence. The thoracic and lumbar spine of eight adult volunteers was imaged using a 3D multi-echo-gradient-echo sagittal MR sequence. A deep-learning method was used to generate synthetic CT images from the MR images. A pre-trained deep-learning network was used for the automatic segmentation of vertebrae from the synthetic CT images. Another deep-learning network was trained for the automatic segmentation of intervertebral discs from the MR images. The automatic segmentations were validated against manual segmentations for two subjects, one with scoliosis, and another with a spine implant. A template mesh of the spine was registered to the segmentations in three steps using a Bayesian coherent point drift algorithm. First, rigid registration was applied on the complete spine. Second, non-rigid registration was used for the individual discs and vertebrae. Third, the complete spine was non-rigidly registered to the individually registered discs and vertebrae. Comparison of the automatic and manual segmentations led to dice-scores of 0.93-0.96 for all vertebrae and discs. The lowest dice-score was in the disc at the height of the implant where artifacts led to under-segmentation. The mean distance between the morphed meshes and the segmentations was below 1 mm. In conclusion, the presented method can be used to automatically generate accurate subject-specific spine models.

Biomechanical and clinical studies on lumbar spine fusion surgery: a review

Low back pain is associated with degenerative disc diseases of the spine. Surgical treatment includes fusion and non-fusion types. The gold standard is fusion surgery, wherein the affected vertebral segment is fused. The common complication of fusion surgery is adjacent segment degeneration (ASD). The ASD often leads to revision surgery, calling for a further fusion of adjacent segments. The existing designs of nonfusion type implants are associated with clinical problems such as subsidence, difficulty in implantation, and the requirement of revision surgeries. Various surgical approaches have been adopted by the surgeons to insert the spinal implants into the affected segment. Over the years, extensive biomechanical investigations have been reported on various surgical approaches and prostheses to predict the outcomes of lumbar spine implantations. Computer models have been proven to be very effective in identifying the best prosthesis and surgical procedure. The objective of the study was to review the literature on biomechanical studies for the treatment of lumbar spinal degenerative diseases. A critical review of the clinical and biomechanical studies on fusion spine surgeries was undertaken. The important modeling parameters, challenges, and limitations of the current studies were identified, showing the future research directions.

Network-based modelling of mechano-inflammatory chondrocyte regulation in early osteoarthritis

Osteoarthritis (OA) is a debilitating joint disease characterized by articular cartilage degradation, inflammation and pain. An extensive range of in vivo and in vitro studies evidences that mechanical loads induce changes in chondrocyte gene expression, through a process known as mechanotransduction. It involves cascades of complex molecular interactions that convert physical signals into cellular response(s) that favor either chondroprotection or cartilage destruction. Systematic representations of those interactions can positively inform early strategies for OA management, and dynamic modelling allows semi-quantitative representations of the steady states of complex biological system according to imposed initial conditions. Yet, mechanotransduction is rarely integrated. Hence, a novel mechano-sensitive network-based model is proposed, in the form of a continuous dynamical system: an interactome of a set of 118 nodes, i.e., mechano-sensitive cellular receptors, second messengers, transcription factors and proteins, related among each other through a specific topology of 358 directed edges is developed. Results show that under physio-osmotic initial conditions, an anabolic state is reached, whereas initial perturbations caused by pro-inflammatory and injurious mechanical loads leads to a catabolic profile of node expression. More specifically, healthy chondrocyte markers (Sox9 and CITED2) are fully expressed under physio-osmotic conditions, and reduced under inflammation, or injurious loadings. In contrast, NF-κB and Runx₂, characteristic of an osteoarthritic chondrocyte, become activated under inflammation or excessive loading regimes. A literature-based evaluation shows that the model can replicate 94% of the experiments tested. Sensitivity analysis based on a factorial design of a treatment shows that inflammation has the strongest influence on chondrocyte metabolism, along with a significant deleterious effect of static compressive loads. At the same time, anti-inflammatory therapies appear as the most promising ones, though the restoration of structural protein production seems to remain a major challenge even in beneficial mechanical environments. The newly developed mechano-sensitive network model for chondrocyte activity reveals a unique potential to reflect load-induced chondroprotection or articular cartilage degradation in different mechano-chemical-environments.

Accelerated Degradation of Poly-ε-caprolactone Composite Scaffolds for Large Bone Defects

This research investigates the accelerated hydrolytic degradation process of both anatomically designed bone scaffolds with a pore size gradient and a rectangular shape (biomimetically designed scaffolds or bone bricks). The effect of material composition is investigated considering poly-ε-caprolactone (PCL) as the main scaffold material, reinforced with ceramics such as hydroxyapatite (HA), β-tricalcium phosphate (TCP) and bioglass at a concentration of 20 wt%. In the case of rectangular scaffolds, the effect of pore size (200 μm, 300 μm and 500 μm) is also investigated. The degradation process (accelerated degradation) was investigated during a period of 5 days in a sodium hydroxide (NaOH) medium. Degraded bone bricks and rectangular scaffolds were measured each day to evaluate the weight loss of the samples, which were also morphologically, thermally, chemically and mechanically assessed. The results show that the PCL/bioglass bone brick scaffolds exhibited faster degradation kinetics in comparison with the PCL, PCL/HA and PCL/TCP bone bricks. Furthermore, the degradation kinetics of rectangular scaffolds increased by increasing the pore size from 500 μm to 200 μm. The results also indicate that, for the same material composition, bone bricks degrade slower compared with rectangular scaffolds. The scanning electron microscopy (SEM) images show that the degradation process was faster on the external regions of the bone brick scaffolds (600 μm pore size) compared with the internal regions (200 μm pore size). The thermal gravimetric analysis (TGA) results show that the ceramic concentration remained constant throughout the degradation process, while differential scanning calorimetry (DSC) results show that all scaffolds exhibited a reduction in crystallinity (Xc), enthalpy (Δm) and melting temperature (Tm) throughout the degradation process, while the glass transition temperature (Tg) slightly increased. Finally, the compression results show that the mechanical properties decreased during the degradation process, with PCL/bioglass bone bricks and rectangular scaffolds presenting higher mechanical properties with the same design in comparison with the other materials.

Recent advancement in finite element analysis of spinal interbody cages: A review

Finite element analysis (FEA) is a widely used tool in a variety of industries and research endeavors. With its application to spine biomechanics, FEA has contributed to a better understanding of the spine, its components, and its behavior in physiological and pathological conditions, as well as assisting in the design and application of spinal instrumentation, particularly spinal interbody cages (ICs). IC is a highly effective instrumentation for achieving spinal fusion that has been used to treat a variety of spinal disorders, including degenerative disc disease, trauma, tumor reconstruction, and scoliosis. The application of FEA lets new designs be thoroughly "tested" before a cage is even manufactured, allowing bio-mechanical responses and spinal fusion processes that cannot easily be experimented upon in vivo to be examined and "diagnosis" to be performed, which is an important addition to clinical and in vitro experimental studies. This paper reviews the recent progress of FEA in spinal ICs over the last six years. It demonstrates how modeling can aid in evaluating the biomechanical response of cage materials, cage design, and fixation devices, understanding bone formation mechanisms, comparing the benefits of various fusion techniques, and investigating the impact of pathological structures. It also summarizes the various limitations brought about by modeling simplification and looks forward to the significant advancement of spine FEA research as computing efficiency and software capabilities increase. In conclusion, in such a fast-paced field, the FEA is critical for spinal IC studies. It helps in quantitatively and visually demonstrating the cage characteristics after implanting, lowering surgeons' learning costs for new cage products, and probably assisting them in determining the best IC for patients.

Immuno-Modulatory Effects of Intervertebral Disc Cells

Low back pain is a highly prevalent, chronic, and costly medical condition predominantly triggered by intervertebral disc degeneration (IDD). IDD is often caused by structural and biochemical changes in intervertebral discs (IVD) that prompt a pathologic shift from an anabolic to catabolic state, affecting extracellular matrix (ECM) production, enzyme generation, cytokine and chemokine production, neurotrophic and angiogenic factor production. The IVD is an immune-privileged organ. However, during degeneration immune cells and inflammatory factors can infiltrate through defects in the cartilage endplate and annulus fibrosus fissures, further accelerating the catabolic environment. Remarkably, though, catabolic ECM disruption also occurs in the absence of immune cell infiltration, largely due to native disc cell production of catabolic enzymes and cytokines. An unbalanced metabolism could be induced by many different factors, including a harsh microenvironment, biomechanical cues, genetics, and infection. The complex, multifactorial nature of IDD brings the challenge of identifying key factors which initiate the degenerative cascade, eventually leading to back pain. These factors are often investigated through methods including animal models, 3D cell culture, bioreactors, and computational models. However, the crosstalk between the IVD, immune system, and shifted metabolism is frequently misconstrued, often with the assumption that the presence of cytokines and chemokines is synonymous to inflammation or an immune response, which is not true for the intact disc. Therefore, this review will tackle immunomodulatory and IVD cell roles in IDD, clarifying the differences between cellular involvements and implications for therapeutic development and assessing models used to explore inflammatory or catabolic IVD environments.

Bone Remodelling Around Solid and Porous Interbody Cages in the Lumbar Spine

Spinal fusion is an effective surgical treatment for intervertebral disc degeneration. However, the consequences of implantation with interbody cages on load transfer and bone remodelling in the vertebral bodies has scarcely been investigated. Using detailed 3D models of an intact and implanted lumbar spine and the strain energy density based bone remodelling algorithm, this study investigated the evolutionary changes in bone density distributions around porous and solid interbody cages. Follower load technique and submodelling approach were employed to simulate applied loading conditions on the lumbar spine models. The study determined the relationship between mechanical properties and parametrical characteristics of porous Body-centered-cubic (BCC) models, which corroborated well with Gibson-Ashby and exponential regression models. Variations in porosity affected the peri-prosthetic stress distributions and bone remodelling around the cages. In comparison to the solid cage, stresses and strains in the cancellous bone decreased with an increase in cage porosity; whereas the range of motion increased. For the solid cage, increase in bone density of 20-28% was predicted in the L4 inferior and L5 superior regions; whereas the model with 78% porosity exhibited a small 3-5% change in bone density. An overall increase of 9-14% bone density was predicted in the L4 and L5 vertebrae after remodelling for solid interbody cages, which may influence disc degeneration in the adjacent segment. In comparison to the solid cage, an interbody cage with 65-78% porosity could be a viable and promising alternative, provided sufficient mechanical strength is offered.

Computational Modeling Intervertebral Disc Pathophysiology: A Review

Lower back pain is a medical condition of epidemic proportion, and the degeneration of the intervertebral disc has been identified as a major contributor. The etiology of intervertebral disc (IVD) degeneration is multifactorial, depending on age, cell-mediated molecular degradation processes and genetics, which is accelerated by traumatic or gradual mechanical factors. The complexity of such intertwined biochemical and mechanical processes leading to degeneration makes it difficult to quantitatively identify cause–effect relationships through experiments. Computational modeling of the IVD is a powerful investigative tool since it offers the opportunity to vary, observe and isolate the effects of a wide range of phenomena involved in the degenerative process of discs. This review aims at discussing the main findings of finite element models of IVD pathophysiology with a special focus on the different factors contributing to physical changes typical of degenerative phenomena. Models presented are subdivided into those addressing role of nutritional supply, progressive biochemical alterations stemming from an imbalance between anabolic and catabolic processes, aging and those considering mechanical factors as the primary source that induces morphological change within the disc. Limitations of the current models, as well as opportunities for future computational modeling work are also discussed.

Patient-Specific Variations in Local Strain Patterns on the Surface of a Trussed Titanium Interbody Cage

Introduction: 3D printed trussed titanium interbody cages may deliver bone stimulating mechanobiological strains to cells attached at their surface. The exact size and distribution of these strains may depend on patient-specific factors, but the influence of these factors remains unknown. Therefore, this study aimed to determine patient-specific variations in local strain patterns on the surface of a trussed titanium interbody fusion cage.

Materials and Methods: Four patients eligible for spinal fusion surgery with the same cage size were selected from a larger database. For these cases, patient-specific finite element models of the lumbar spine including the same trussed titanium cage were made. Functional dynamics of the non-operated lumbar spinal segments, as well as local cage strains and caudal endplate stresses at the operated segment, were evaluated under physiological extension/flexion movement of the lumbar spine.

Results: All patient-specific models revealed physiologically realistic functional dynamics of the operated spine. In all patients, approximately 30% of the total cage surface experienced strain values relevant for preserving bone homeostasis and stimulating bone formation. Mean caudal endplate contact pressures varied up to 10 MPa. Both surface strains and endplate contact pressures varied more between loading conditions than between patients.

Conclusions: This study demonstrates the applicability of patient-specific finite element models to quantify the impact of patient-specific factors such as bone density, degenerative state of the spine, and spinal curvature on interbody cage loading. In the future, the same framework might be further developed in order to establish a pipeline for interbody cage design optimizations.

Stability Evaluation of Different Oblique Lumbar Interbody Fusion Constructs in Normal and Osteoporotic Condition - A Finite Element Based Study

Introduction: In developed countries, the age structure of the population is currently undergoing an upward shift, resulting a decrease in general bone quality and surgical durability. Over the past decade, oblique lumbar interbody fusion (OLIF) has been globally accepted as a minimally invasive surgical technique. There are several stabilization options available for OLIF cage fixation such as self-anchored stand-alone (SSA), lateral plate-screw (LPS), and bilateral pedicle screw (BPS) systems. The constructs’ stability are crucial for the immediate and long-term success of the surgery. The aim of this study is to investigate the biomechanical effect of the aforementioned constructs, using finite element analysis with different bone qualities (osteoporotic and normal).

Method: A bi-segmental (L2–L4) finite element (FE) model was created, using a CT scan of a 24-year-old healthy male. After the FE model validation, CAD geometries of the implants were inserted into the L3–L4 motion segment during a virtual surgery. For the simulations, a 150 N follower load was applied on the models, then 10 Nm of torque was used in six general directions (flexion, extension, right/left bending, and right/left rotation), with different bone material properties.

Results: The smallest segmental (L3–L4) ROM (range of motion) was observed in the BPS system, except for right bending. Osteoporosis increased ROMs in all constructs, especially in the LPS system (right bending increase: 140.26%). Osteoporosis also increased the caudal displacement of the implanted cage in all models (healthy bone: 0.06 ± 0.03 mm, osteoporosis: 0.106 ± 0.07 mm), particularly with right bending, where the displacement doubled in SSA and LPS constructs. The displacement of the screws inside the L4 vertebra increased by 59% on average (59.33 ± 21.53%) due to osteoporosis (100% in LPS, rotation). BPS-L4 screw displacements were the least affected by osteoporosis.

Conclusions: The investigated constructs provide different levels of stability to the spine depending on the quality of the bone, which can affect the outcome of the surgery. In our model, the BPS system was found to be the most stable construct in osteoporosis. The presented model, after further development, has the potential to help the surgeon in planning a particular spinal surgery by adjusting the stabilization type to the patient’s bone quality.

Evidence-Based Network Modelling to Simulate Nucleus Pulposus Multicellular Activity in Different Nutritional and Pro-Inflammatory Environments

Initiation of intervertebral disc degeneration is thought to be biologically driven. This reflects a process, where biochemical and mechanical stimuli affect cell activity (CA) that compromise the tissue strength over time. Experimental research enhanced our understanding about the effect of such stimuli on different CA, such as protein synthesis or mRNA expression. However, it is still unclear how cells respond to their native environment that consists of a "cocktail" of different stimuli that might locally vary. This work presents an interdisciplinary approach of experimental and in silico research to approximate Nucleus Pulposus CA within multifactorial biochemical environments. Thereby, the biochemical key stimuli glucose, pH, and the proinflammatory cytokines TNF-α and IL1β were considered that were experimentally shown to critically affect CA. To this end, a Nucleus Pulposus multicellular system was modelled. It integrated experimental findings from in vitro studies of human or bovine Nucleus Pulposus cells, to relate the individual effects of targeted stimuli to alterations in CA. Unknown stimulus-CA relationships were obtained through own experimental 3D cultures of bovine Nucleus Pulposus cells in alginate beads. Translation of experimental findings into suitable parameters for network modelling approaches was achieved thanks to a new numerical approach to estimate the individual sensitivity of a CA to each stimulus type. Hence, the effect of each stimulus type on a specific CA was assessed and integrated to approximate a multifactorial stimulus environment. Tackled CA were the mRNA expressions of Aggrecan, Collagen types I & II, MMP3, and ADAMTS4. CA was assessed for four different proinflammatory cell states; non-inflamed and inflamed for IL1β, TNF-α or both IL1β&TNF-α. Inflamed cell clusters were eventually predicted in a multicellular 3D agent-based model. Experimental results showed that glucose had no significant impact on proinflammatory cytokine or ADAMTS4 mRNA expression, whereas TNF-α caused a significant catabolic shift in most explored CA. In silico results showed that the presented methodology to estimate the sensitivity of a CA to a stimulus type importantly improved qualitative model predictions. However, more stimuli and/or further experimental knowledge need to be integrated, especially regarding predictions about the possible progression of inflammatory environments under adverse nutritional conditions. Tackling the multicellular level is a new and promising approach to estimate manifold responses of intervertebral disc cells. Such a top-down high-level network modelling approach allows to obtain information about relevant stimulus environments for a specific CA and could be shown to be suitable to tackle complex biological systems, including different proinflammatory cell states. The development of this methodology required a close interaction with experimental research. Thereby, specific experimental needs were derived from systematic in silico approaches and obtained results were directly used to enhance model predictions, which reflects a novelty in this research field. Eventually, the presented methodology provides modelling solutions suitable for multiscale approaches to contribute to a better understanding about dynamics over multiple spatial scales. Future work should focus on an amplification of the stimulus environment by integrating more key relevant stimuli, such as mechanical loading parameters, in order to better approximate native physiological environments.

Biomechanical Evaluation of the Effect of Minimally Invasive Spine Surgery Compared with Traditional Approaches in Lifting Tasks

Fusion of spinal vertebrae can be accomplished by different surgical approaches. We investigated Traditional Open Spine Surgery (TOSS) versus Minimally Invasive Spine Surgery (MISS). While TOSS sacrifices spine muscles originating or inserting on the affected vertebrae, MISS seeks to minimize the approach-related morbidity and preserve the tendon attachments of the muscles in the area. We captured 3-D motions of the full body of one healthy subject performing a variety of 10 kg box lifting operations representing activities-of-daily-living that are likely to challenge the spine biomechanically. The motion data were transferred to a full-body biomechanical model with a detailed representation of the biomechanics of the spine, and simulations of the internal spine loads and muscle forces were performed under a baseline configuration and muscle configurations typical for TOSS respectively MISS for the cases of L3/L4, L4/L5, L5/S1, L4/S1 and L3/L5 fusions. The computational model was then used to investigate the biomechanical differences between surgeries. The simulations revealed that joint reaction forces are more affected by both surgical approaches for lateral lifting motions than for sagittal plane motions, and there are indications that individuals with fused joints, regardless of the approach, should be particularly careful with asymmetrical lifts. The MISS and TOSS approaches shift the average loads of different muscle groups in different ways. TOSS generally leads to higher post-operative muscle loads than MISS in the investigated cases, but the differences are smaller than could be expected, given the differences of surgical technique.

Saline-polyethylene glycol blends preserve in vitro annulus fibrosus hydration and mechanics: An experimental and finite-element analysis

Precise control of tissue water content is essential for ensuring accurate, repeatable, and physiologically relevant measurements of tissue mechanics and biochemical composition. While previous studies have found that saline and polyethylene glycol (PEG) blends were effective at controlling tendon and ligament hydration levels, this work has yet to be extended to the annulus fibrosus (AF). Thus, the first objective of this study was to determine and validate an optimal buffer solution for targeting and maintaining hydration levels of tissue-level AF specimens in vitro. This was accomplished by measuring the transient swelling behavior of bovine AF specimens in phosphate-buffered saline (PBS) and PEG buffers across a wide range of concentrations. Sub-failure, failure, and post-failure mechanics were measured to determine the relationship between changes in tissue hydration and tensile mechanical response. The effect of each buffer solution on tissue composition was also assessed. The second objective of this study was to assess the feasibility and effectiveness of using multi-phasic finite element models to investigate tissue swelling and mechanical responses in different external buffer solutions. A solution containing 6.25%w/v PBS and 6.25%w/v PEG effectively maintained tissue-level AF specimen hydration at fresh-frozen levels after 18 h in solution. Modulus, failure stress, failure strain, and post-failure toughness of specimens soaked in this solution for 18 h closely matched those of fresh-frozen specimens. In contrast, specimens soaked in 0.9%w/v PBS swelled over 100% after 18 h and exhibited significantly diminished sub-failure and failure properties compared to fresh-frozen controls. The increased cross-sectional area with swelling contributed to but was not sufficient to explain the diminished mechanics of PBS-soaked specimens, suggesting additional sub-tissue scale mechanisms. Computational simulations of these specimens generally agreed with experimental results, highlighting the feasibility and importance of including a fluid-phase description when models aim to provide accurate predictions of biological tissue responses. As numerous previous studies suggest that tissue hydration plays a central role in maintaining proper mechanical and biological function, robust methods for controlling hydration levels are essential as the field advances in probing the relationship between tissue hydration, aging, injury, and disease.

Simulating intervertebral disc cell behaviour within 3D multifactorial environments

MOTIVATION

Low back pain is responsible for more global disability than any other condition. Its incidence is closely related to intervertebral disc (IVD) failure, which is likely caused by an accumulation of microtrauma within the IVD. Crucial factors in microtrauma development are not entirely known yet, probably because their exploration in vivo or in vitro remains tremendously challenging. In silico modelling is, therefore, definitively appealing, and shall include approaches to integrate influences of multiple cell stimuli at the microscale. Accordingly, this study introduces a hybrid Agent-based (AB) model in IVD research and exploits network modelling solutions in systems biology to mimic the cellular behaviour of Nucleus Pulposus cells exposed to a 3D multifactorial biochemical environment, based on mathematical integrations of existing experimental knowledge. Cellular activity reflected by mRNA expression of Aggrecan, Collagen type I, Collagen type II, MMP-3 and ADAMTS were calculated for inflamed and non-inflamed cells. mRNA expression over long periods of time is additionally determined including cell viability estimations. Model predictions were eventually validated with independent experimental data.

RESULTS

As it combines experimental data to simulate cell behaviour exposed to a multifactorial environment, the present methodology was able to reproduce cell death within 3 days under glucose deprivation and a 50% decrease in cell viability after 7 days in an acidic environment. Cellular mRNA expression under non-inflamed conditions simulated a quantifiable catabolic shift under an adverse cell environment, and model predictions of mRNA expression of inflamed cells provide new explanation possibilities for unexpected results achieved in experimental research.

AVAILABILITYAND IMPLEMENTATION

The AB model as well as used mathematical functions were built with open source software. Final functions implemented in the AB model and complete AB model parameters are provided as Supplementary Material. Experimental input and validation data were provided through referenced, published papers. The code corresponding to the model can be shared upon request and shall be reused after proper training.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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

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

Clinical-Instrumental Results and Analysis of Functional Activity Restoration in Professional Athletes After Lumbar Total Disk Replacement

OBJECTIVE

To analyze the clinical and biomechanical outcome of professional athletes after lumbar total disk replacement, with a focus on restoration of the functional activity.

METHODS

This nonrandomized retrospective single-center study included 11 professional athletes who underwent lumbar disc replacement surgery using the prosthesis M6-L (Spinal Kinetics, Schaffhausen, Switzerland). The average postoperative follow-up was 3.18 ± 1.14 years. The following outcomes were evaluated: intensity of pain in the lumbar spine and lower limbs, Oswestry Disability Index, Short-Form 36, complications, time of return to previous sports activity, range of movement, degree of lumbar lordosis, degenerative changes of the adjacent levels, and degree of heterotopic ossification.

RESULTS

The operated patients reported significant decrease of pain on visual analog scale (P < 0.001) as well as significant improvement of Oswestry Disability Index (P = 0.001) and Short-Form 36 (P < 0.001). For the duration of follow-up, the patients maintained segmental range of motion at L4-L5 (P = 0.04) and L5-S1 (P = 0.03) levels. There was also some statistically insignificant increase of global lumbar lordosis (P = 0.84). We did not identify any significant degeneration of the adjacent intervertebral disks (P > 0.05) or progression of the facet joint degenerative changes at the implantation level and in the adjacent segments (P > 0.05). One patient (9.1%) developed grade I heterotopic ossification 5 years after surgery and in 1 patient (9.1%), a lesion of superior hypogastric plexus was recorded. The average time of return to previous sports activity was 9.72 ± 3.03 weeks.

CONCLUSIONS

Total lumbar disc replacement using M6-L prosthesis in professional athletes made it possible to achieve statistically significant reduction of pain and facilitated early return to normal sports activities. In our opinion, preservation of movement of the operated lumbar segment can help to reduce the mechanical stress with beneficial impact on the rate of degeneration of the adjacent level.

Maintaining Bone Health in the Lumbar Spine: Routine Activities Alone Are Not Enough

Public health organisations typically recommend a minimum amount of moderate intensity activities such as walking or cycling for two and a half hours a week, combined with some more demanding physical activity on at least 2 days a week to maintain a healthy musculoskeletal condition. For populations at risk of bone loss in the lumbar spine, these guidelines are particularly relevant. However, an understanding of how these different activities are influential in maintaining vertebral bone health is lacking. A predictive structural finite element modelling approach using a strain-driven algorithm was developed to study mechanical stimulus and bone adaptation in the lumbar spine under various physiological loading conditions. These loading conditions were obtained with a previously developed full-body musculoskeletal model for a range of daily living activities representative of a healthy lifestyle. Activities of interest for the simulations include moderate intensity activities involving limited spine movements in all directions such as, walking, stair ascent and descent, sitting down and standing up, and more demanding activities with large spine movements during reaching and lifting tasks. For a combination of moderate and more demanding activities, the finite element model predicted a trabecular and cortical bone architecture representative of a healthy vertebra. When more demanding activities were removed from the simulations, areas at risk of bone degradation were observed at all lumbar levels in the anterior part of the vertebral body, the transverse processes and the spinous process. Moderate intensity activities alone were found to be insufficient in providing a mechanical stimulus to prevent bone degradation. More demanding physical activities are essential to maintain bone health in the lumbar spine.

A novel three-dimensional volumetric method to measure indirect decompression after percutaneous cement discoplasty

PURPOSE

Percutaneous cement discoplasty (PCD) is a minimally invasive surgical option to treat patients who suffer from the consequences of advanced disc degeneration. As the current two-dimensional methods can inappropriately measure the difference in the complex 3D anatomy of the spinal segment, our aim was to develop and apply a volumetric method to measure the geometrical change in the surgically treated segments.

METHODS

Prospective clinical and radiological data of 10 patients who underwent single- or multilevel PCD was collected. Pre- and postoperative CT scan-based 3D reconstructions were performed. The injected PMMA (Polymethylmethacrylate) induced lifting of the cranial vertebra and the following volumetric change was measured by subtraction of the geometry of the spinal canal from a pre- and postoperatively predefined cylinder. The associations of the PMMA geometry and the volumetric change of the spinal canal with clinical outcome were determined.

RESULTS

Change in the spinal canal volume (ΔV) due to the surgery proved to be significant (mean ΔV = 2266.5 ± 1172.2 mm3, n = 16; p = 0.0004). A significant, positive correlation was found between ΔV, the volume and the surface of the injected PMMA. A strong, significant association between pain intensity (low back and leg pain) and the magnitude of the volumetric increase of the spinal canal was shown (ρ = 0.772, p = 0.009 for LBP and ρ = 0.693, p = 0.026 for LP).

CONCLUSION

The developed method is accurate, reproducible and applicable for the analysis of any other spinal surgical method. The volume and surface area of the injected PMMA have a predictive power on the extent of the indirect spinal canal decompression. The larger the ΔV the higher clinical benefit was achieved with the PCD procedure.

THE TRANSLATIONAL POTENTIAL OF THIS ARTICLE

The developed method has the potential to be integrated into clinical software's to evaluate the efficacy of different surgical procedures based on indirect decompression effect such as PCD, anterior lumbar interbody fusion (ALIF), lateral lumbar interbody fusion (LLIF), oblique lumbar interbody fusion (OLIF), extreme lateral interbody fusion (XLIF). The intraoperative use of the method will allow the surgeon to respond if the decompression does not reach the desired level.

Biomechanical Investigation Between Rigid and Semirigid Posterolateral Fixation During Daily Activities: Geometrically Parametric Poroelastic Finite Element Analyses

While spinal fusion using rigid rods remains the gold standard treatment modality for various lumbar degenerative conditions, its adverse effects, including accelerated adjacent segment disease (ASD), are well known. In order to better understand the performance of semirigid constructs using polyetheretherketone (PEEK) in fixation surgeries, the objective of this study was to analyze the biomechanical performance of PEEK versus Ti rods using a geometrically patient-specific poroelastic finite element (FE) analyses. Ten subject-specific preoperative models were developed, and the validity of the models was evaluated with previous studies. Furthermore, FE models of those lumbar spines were regenerated based on postoperation images for posterolateral fixation at the L4–L5 level. Biomechanical responses for instrumented and adjacent intervertebral discs (IVDs) were analyzed and compared subjected to static and cyclic loading. The preoperative model results were well comparable with previous FE studies. The PEEK construct demonstrated a slightly increased range of motion (ROM) at the instrumented level, but decreased ROM at adjacent levels, as compared with the Ti. However, no significant changes were detected during axial rotation. During cyclic loading, disc height loss, fluid loss, axial stress, and collagen fiber strain in the adjacent IVDs were higher for the Ti construct when compared with the intact and PEEK models. Increased ROM, experienced stress in AF, and fiber strain at adjacent levels were observed for the Ti rod group compared with the intact and PEEK rod group, which can indicate the risk of ASD for rigid fixation. Similar to the aforementioned pattern, disc height loss and fluid loss were significantly higher at adjacent levels in the Ti rod group after cycling loading which alter the fluid–solid interaction of the adjacent IVDs. This phenomenon debilitates the damping quality, which results in disc disability in absorbing stress. Such finding may suggest the advantage of using a semirigid fixation system to decrease the chance of ASD.

Misaligned spinal rods can induce high internal forces consistent with those observed to cause screw pullout and disc degeneration

BACKGROUND CONTEXT

Manual contouring of spinal rods is often required intraoperatively for proper alignment of the rods within the pedicle screw heads. Residual misalignments are frequently reduced by using dedicated reduction devices. The forces exerted by these devices, however, are uncontrolled and may lead to excessive reaction forces. As a consequence, screw pullout might be provoked and surrounding tissue may experience unfavorable biomechanical loads. The corresponding loads and induced tissue deformations are however not well identified. Additionally, whether the forced reduction alters the biomechanical behavior of the lumbar spine during physiological movements postoperatively, remains unexplored.

PURPOSE

To predict whether the reduction of misaligned posterior instrumentation might result in clinical complications directly after reduction and during a subsequent physiological flexion movement.

STUDY DESIGN

Finite element analysis.

METHODS

A patient-specific, total lumbar (L1-S1) spine finite element model was available from previous research. The model consists of poro-elastic intervertebral discs with Pfirrmann grade-dependent material parameters, with linear elastic bone tissue with stiffness values related to the local bone density, and with the seven major ligaments per spinal motion segment described as nonlinear materials. Titanium instrumentation was implemented in this model to simulate a L4, L5, and S1 posterolateral fusion. Next, coronal and sagittal misalignments of 6 mm each were introduced between the rod and the screw head at L4. These misalignments were computationally reduced and a physiological flexion movement of 15° was prescribed. Non-instrumented and well-aligned instrumented models were used as control groups.

RESULTS

Pulling forces up to 1.0 kN were required to correct the induced misalignments of 6 mm. These forces affected the posture of the total lumbar spine, as motion segments were predicted to rotate up to 3 degrees and rotations propagated proximally to and even affect the L1-2 level. The facet contact pressures in the corrected misaligned models were asymmetrical suggesting non-physiological joint loading in the misaligned models. In addition, the discs and vertebrae experienced abnormally high forces as a result of the correction procedure. These effects were more pronounced after a 15° flexion movement following forced reduction.

CONCLUSIONS

The results of this study indicate that the correction of misaligned posterior instrumentation can result in high forces at the screws consistent with those reported to cause screw pullout, and may cause high-tissue strains in adjacent and downstream spinal segments.

CLINICAL SIGNIFICANCE

Proper alignment of spinal posterior instrumentation may reduce clinical complications secondary to unfavorable biomechanics.

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.

Automated Pipeline to Generate Anatomically Accurate Patient-Specific Biomechanical Models of Healthy and Pathological FSUs

State-of-the-art preoperative biomechanical analysis for the planning of spinal surgery not only requires the generation of three-dimensional patient-specific models but also the accurate biomechanical representation of vertebral joints. The benefits offered by computational models suitable for such purposes are still outweighed by the time and effort required for their generation, thus compromising their applicability in a clinical environment. In this work, we aim to ease the integration of computerized methods into patient-specific planning of spinal surgery. We present the first pipeline combining deep learning and finite element methods that allows a completely automated model generation of functional spine units (FSUs) of the lumbar spine for patient-specific FE simulations (FEBio). The pipeline consists of three steps: (a) multiclass segmentation of cropped 3D CT images containing lumbar vertebrae using the DenseVNet network, (b) automatic landmark-based mesh fitting of statistical shape models onto 3D semantic segmented meshes of the vertebral models, and (c) automatic generation of patient-specific FE models of lumbar segments for the simulation of flexion-extension, lateral bending, and axial rotation movements. The automatic segmentation of FSUs was evaluated against the gold standard (manual segmentation) using 10-fold cross-validation. The obtained Dice coefficient was 93.7% on average, with a mean surface distance of 0.88 mm and a mean Hausdorff distance of 11.16 mm (N = 150). Automatic generation of finite element models to simulate the range of motion (ROM) was successfully performed for five healthy and five pathological FSUs. The results of the simulations were evaluated against the literature and showed comparable ROMs in both healthy and pathological cases, including the alteration of ROM typically observed in severely degenerated FSUs. The major intent of this work is to automate the creation of anatomically accurate patient-specific models by a single pipeline allowing functional modeling of spinal motion in healthy and pathological FSUs. Our approach reduces manual efforts to a minimum and the execution of the entire pipeline including simulations takes approximately 2 h. The automation, time-efficiency and robustness level of the pipeline represents a first step toward its clinical integration.

Development of a Computer-Aided Design and Finite Element Analysis Combined Method for Affordable Spine Surgical Navigation With 3D-Printed Customized Template

Introduction: Revision surgery of a previous lumbosacral non-union is highly challenging, especially in case of complications, such as a broken screw at the first sacral level (S1). Here, we propose the implementation of a new method based on the CT scan of a clinical case using 3D reconstruction, combined with finite element analysis (FEA), computer-assisted design (CAD), and 3D-printing technology to provide accurate surgical navigation to aid the surgeon in performing the optimal surgical technique by inserting a pedicle screw at the S1 level.

Materials and Methods: A step-by-step approach was developed and performed as follows: (1) Quantitative CT based patient-specific FE model of the sacrum was created. (2) The CAD model of the pedicle screw was inserted into the sacrum model in a bicortical convergent and a monocortical divergent position, by overcoming the geometrical difficulty caused by the broken screw. (3) Static FEAs (Abaqus, Dassault Systemes) were performed using 500 N tensile load applied to the screw head. (4) A template with two screw guiding structures for the sacrum was designed and manufactured using CAD design and 3D-printing technologies, and investment casting. (5) The proposed surgical technique was performed on the patient-specific physical model created with the FDM printing technology. The patient-specific model was CT scanned and a comparison with the virtual plan was performed to evaluate the template accuracy

Results: FEA results proved that the modified bicortical convergent insertion is stiffer (6,617.23 N/mm) compared to monocortical divergent placement (2,989.07 N/mm). The final template was created via investment casting from cobalt-chrome. The template design concept was shown to be accurate (grade A, Gertzbein-Robbins scale) based on the comparison of the simulated surgery using the patient-specific physical model and the 3D virtual surgical plan.

Conclusion: Compared to the conventional surgical navigation techniques, the presented method allows the consideration of the patient-specific biomechanical parameters; is more affordable, and the intraoperative X-ray exposure can be reduced. This new patient- and condition-specific approach may be widely used in revision spine surgeries or in challenging primary cases after its further clinical validation.

Multiscale Regulation of the Intervertebral Disc: Achievements in Experimental, In Silico, and Regenerative Research

Intervertebral disc (IVD) degeneration is a major risk factor of low back pain. It is defined by a progressive loss of the IVD structure and functionality, leading to severe impairments with restricted treatment options due to the highly demanding mechanical exposure of the IVD. Degenerative changes in the IVD usually increase with age but at an accelerated rate in some individuals. To understand the initiation and progression of this disease, it is crucial to identify key top-down and bottom-up regulations' processes, across the cell, tissue, and organ levels, in health and disease. Owing to unremitting investigation of experimental research, the comprehension of detailed cell signaling pathways and their effect on matrix turnover significantly rose. Likewise, in silico research substantially contributed to a holistic understanding of spatiotemporal effects and complex, multifactorial interactions within the IVD. Together with important achievements in the research of biomaterials, manifold promising approaches for regenerative treatment options were presented over the last years. This review provides an integrative analysis of the current knowledge about (1) the multiscale function and regulation of the IVD in health and disease, (2) the possible regenerative strategies, and (3) the in silico models that shall eventually support the development of advanced therapies.

Conventional finite element models estimate the strength of metastatic human vertebrae despite alterations of the bone's tissue and structure

INTRODUCTION

Pathologic vertebral fractures are a major clinical concern in the management of cancer patients with metastatic spine disease. These fractures are a direct consequence of the effect of bone metastases on the anatomy and structure of the vertebral bone. The goals of this study were twofold. First, we evaluated the effect of lytic, blastic and mixed (both lytic and blastic) metastases on the bone structure, on its material properties, and on the overall vertebral strength. Second, we tested the ability of bone mineral content (BMC) measurements and standard FE methodologies to predict the strength of real metastatic vertebral bodies.

METHODS

Fifty-seven vertebral bodies from eleven cadaver spines containing lytic, blastic, and mixed metastatic lesions from donors with breast, esophageal, kidney, lung, or prostate cancer were scanned using micro-computed tomography (μCT). Based on radiographic review, twelve vertebrae were selected for nanoindentation testing, while the remaining forty-five vertebrae were used for assessing their compressive strength. The μCT reconstruction was exploited to measure the vertebral BMC and to establish two finite element models. 1) a micro finite element (μFE) model derived at an image resolution of 24.5 μm and 2) homogenized FE (hFE) model derived at a resolution of 0.98 mm. Statistical analyses were conducted to measure the effect of the bone metastases on BV/TV, indentation modulus (E_it), ratio of plastic/total work (W_Pl/W_tot), and in vitro vertebral strength (F_exp). The predictive value of BMC, μFE stiffness, and hFE strength were evaluated against the in vitro measurements.

RESULTS

Blastic vertebral bodies exhibit significantly higher BV/TV compared to the mixed (p = 0.0205) and lytic (p = 0.0216) vertebral bodies. No significant differences were found between lytic and mixed vertebrae (p = 0.7584). Blastic bone tissue exhibited a 5.8% lower median E_it (p< 0.001) and a 3.3% lower median W_pl/W_tot (p<0.001) compared to non-involved bone tissue. No significant differences were measured between lytic and non-involved bone tissues. F_exp ranged from 1.9 to 13.8 kN, was strongly associated with hFE strength (R²=0.78, p< 0.001) and moderately associated with BMC (R²=0.66, p< 0.001) and μFE stiffness (R²=0.66, p< 0.001), independently of the lesion type.

DISCUSSION

Our findings show that tumour-induced osteoblastic metastases lead to slightly, but significantly lower bone tissue properties compared to controls, while osteolytic lesions appear to have a negligible impact. These effects may be attributed to the lower mineralization and woven nature of bone forming in blastic lesions whilst the material properties of bone in osteolytic vertebrae appeared little changed. The moderate association between BMC- and FE-based predictions to fracture strength suggest that vertebral strength is affected by the changes of bone mass induced by the metastatic lesions, rather than altered tissue properties. In a broader context, standard hFE approaches generated from CTs at clinical resolution are robust to the lesion type when predicting vertebral strength. These findings open the door for the development of FE-based prediction tools that overcomes the limitations of BMC in accounting for shape and size of the metastatic lesions. Such tools may help clinicians to decide whether a patient needs the prophylactic fixation of an impending fracture.

Finite element modelling of hybrid stabilization systems for the human lumbar spine

Intersomatic fusion is a very popular treatment for spinal diseases associated with intervertebral disc degeneration. The effects of three different hybrid stabilization systems on both range of motion and intradiscal pressure were investigated, as there is no consensus in the literature about the efficiency of these systems. Finite element simulations were designed to predict the variations of range of motion and intradiscal pressure from intact to implanted situations. After hybrid stabilization system implantation, L4-L5 level did not lose its motion completely, while L5-S1 had no mobility as a consequence of disc removal and fusion process. BalanC hybrid stabilization system represented higher mobility at the index level, reduced intradiscal pressure of adjacent level, but caused to increment in range of motion by 20% under axial rotation. Higher tendency by 93% to the failure was also detected under axial rotation. Dynesys hybrid stabilization system represented more restricted motion than BalanC, and negligible effects to the adjacent level. B-DYN hybrid stabilization system was the most rigid one among all three systems. It reduced intradiscal pressure and range of motion at the adjacent level except from motion under axial rotation being increased by 13%. Fracture risk of B-DYN and Dynesys Transition Optima components was low when compared with BalanC. Mobility of the adjacent level around axial direction should be taken into account in case of implantation with BalanC and B-DYN systems, as well as on the development of new designs. Having these findings in mind, it is clear that hybrid systems need to be further tested, both clinically and numerically, before being considered for common use.

Influence of cement-augmented pedicle screw instrumentation in an osteoporotic lumbosacral spine over the adjacent segments: a 3D finite element study

Abstract

Purpose

To compare the effect of conventional pedicle screw (CPS) and cement-augmented pedicle screw instrumentation (CAPSI) on adjacent segment degeneration (ASD).

Methods

A normal male volunteer without a history of spinal disease was selected, lumbar CT data was collected, an intact L3-S1 three-dimensional finite element model was created by software including Mimics, Geomagic, and SolidWorks, and the fixation methods were performed accordingly. A common pedicle screw model and a cement-augmented pedicle screw model of L4–L5 with fusion and internal fixation were constructed. With ANSYS Workbench 17.0, a 500 N load was applied to the upper surface of L3 to simulate the weight of a human body, and a 7.5 N m moment was applied at the neutral point to simulate flexion, extension, left/right bending, left/right rotation of the spine. The peak von Mises stress of intervertebral disc and the range of motion (ROM) on the adjacent segments (L3–4 and L5–S1) were compared.

Results

The validity of the intact model shows that the ROM of the model is similar to that of a cadaveric study. Compared with the intact model, CPS model and CAPSI model in all motion patterns increased the ROM of adjacent segments. The intervertebral disc stress and the ROM of adjacent segments were found to be higher in the CAPSI model than in the CPS model, especially in L3–4.

Conclusion

In general, the biomechanical analysis of an osteoporotic lumbar spine showed that both CPS and CAPSI can increase the ROM and disc stresses of osteoporotic lumbar models, and compared with CPS, CAPSI is more likely to increase the potential risk of adjacent segment degeneration.