Unveiling interactions between intervertebral disc morphologies and mechanical behavior through personalized finite element modeling.

Introduction: Intervertebral Disc (IVD) Degeneration (IDD) is a significant health concern, potentially influenced by mechanotransduction. However, the relationship between the IVD phenotypes and mechanical behavior has not been thoroughly explored in local morphologies where IDD originates. This work unveils the interplays among morphological and mechanical features potentially relevant to IDD through Abaqus UMAT simulations. Methods: A groundbreaking automated method is introduced to transform a calibrated, structured IVD finite element (FE) model into 169 patient-personalized (PP) models through a mesh morphing process. Our approach accurately replicates the real shapes of the patient's Annulus Fibrosus (AF) and Nucleus Pulposus (NP) while maintaining the same topology for all models. Using segmented magnetic resonance images from the former project MySpine, 169 models with structured hexahedral meshes were created employing the Bayesian Coherent Point Drift++ technique, generating a unique cohort of PP FE models under the Disc4All initiative. Machine learning methods, including Linear Regression, Support Vector Regression, and eXtreme Gradient Boosting Regression, were used to explore correlations between IVD morphology and mechanics. Results: We achieved PP models with AF and NP similarity scores of 92.06\% and 92.10\% compared to the segmented images. The models maintained good quality and integrity of the mesh. The cartilage endplate (CEP) shape was represented at the IVD-vertebra interfaces, ensuring personalized meshes. Validation of the constitutive model against literature data showed a minor relative error of 5.20%. Discussion: Analysis revealed the influential impact of local morphologies on indirect mechanotransduction responses, highlighting the roles of heights, sagittal areas, and volumes. While the maximum principal stress was influenced by morphologies such as heights, the disc's ellipticity influenced the minimum principal stress. Results suggest the CEPs are not influenced by their local morphologies but by those of the AF and NP. The generated free-access repository of individual disc characteristics is anticipated to be a valuable resource for the scientific community with a broad application spectrum.

Dataset of Finite Element Models of Normal and Deformed Thoracolumbar Spine.

Adult spine deformity (ASD) is prevalent and leads to a sagittal misalignment in the vertebral column. Computational methods, including Finite Element (FE) Models, have emerged as valuable tools for investigating the causes and treatment of ASD through biomechanical simulations. However, the process of generating personalised FE models is often complex and time-consuming. To address this challenge, we present a dataset of FE models with diverse spine morphologies that statistically represent real geometries from a cohort of patients. These models are generated using EOS images, which are utilized to reconstruct 3D surface spine models. Subsequently, a Statistical Shape Model (SSM) is constructed, enabling the adaptation of a FE hexahedral mesh template for both the bone and soft tissues of the spine through mesh morphing. The SSM deformation fields facilitate the personalization of the mean hexahedral FE model based on sagittal balance measurements. Ultimately, this new hexahedral SSM tool offers a means to generate a virtual cohort of 16807 thoracolumbar FE spine models, which are openly shared in a public repository.

Proximal Junction Failure in Spine Surgery: Integrating Geometrical and Biomechanical Global Descriptors Improves GAP Score-Based Assessment.

STUDY DESIGN: Retrospective observational study. OBJECTIVE: Biomechanical and geometrical descriptors are used to improve global alignment and proportion (GAP) prediction accuracy to detect proximal junctional failure (PJF). SUMMARY OF BACKGROUND DATA: PJF is probably the most important complication after sagittal imbalance surgery. The GAP score has been introduced as an effective predictor for PJF, but it fails in certain situations. In this study, 112 patient records were gathered (57 PJF; 55 controls) with biomechanical and geometrical descriptors measured to stratify control and failure cases. PATIENTS AND METHODS: Biplanar EOS radiographs were used to build 3-dimensional full-spine models and determine spinopelvic sagittal parameters. The bending moment (BM) was calculated as the upper body mass times, the effective distance to the body center of mass at the adjacent upper instrumented vertebra +1. Other geometrical descriptors such as full balance index (FBI), spino-sacral angle (SSA), C7 plumb line/sacrofemoral distance ratio (C7/SFD ratio), T1-pelvic angle (TPA), and cervical inclination angle (CIA) were also evaluated. The respective abilities of the GAP, FBI, SSA, C7/SFD, TPA, CIA, body weight, body mass index, and BM to discriminate PJF cases were analyzed through receiver operating characteristic curves and corresponding areas under the curve (AUC). RESULTS: GAP (AUC = 0.8816) and FBI (AUC = 0.8933) were able to discriminate PJF cases but the highest discrimination power (AUC = 0.9371) was achieved with BM at upper instrumented vertebra + 1. Parameter cutoff analyses provided quantitative thresholds to characterize the control and failure groups and led to improved PJF discrimination, with GAP and BM being the most important contributors. SSA (AUC = 0.2857), C7/SFD (AUC = 0.3143), TPA (AUC = 0.5714), CIA (AUC = 0.4571), body weight (AUC = 0.6319), and body mass index (AUC = 0.7716) did not adequately predict PJF. CONCLUSION: BM reflects the quantitative biomechanical effect of external loads and can improve GAP accuracy. Sagittal alignments and mechanical integrated scores could be used to better prognosticate the risk of PJF.