Vertebral endplate morphology follows bone remodeling principles

Prediction of heterotopic ossification on the cervical spine with offset of the artificial disc - A finite element study

Heterotopic ossification (HO) is a significant complication of cervical total disc replacement (TDR), often leading to fusion and negating the intended benefits of motion preservation. Although clinical factors associated with HO formation are known, the exact biomechanical mechanism remains unclear. This study aims to predict HO formation after Mobi-C disc replacement at the C5-C6 level using a validated finite element model (FEM) of the cervical spine (C2-T1) under physiological loading. The results revealed that the Mobi-C disc increased the range of motion (ROM) at the implanted level by 52 % under flexion and extension, while adjacent levels exhibited a 2-5 % reduction. Following HO formation, ROM at the implanted level decreased by 67-76 % in flexion and extension, respectively, while adjacent levels showed a moderate increase of 5-8 %. Additionally, intradiscal pressure at the adjacent levels increased by up to 60 % in extension, mimicking fusion-like behavior. HO volume was 678 mm³ for the ideal implant position, to 760 mm³ (+12 %) for a 0.5 mm offset and 800 mm³ (+17 %) for a 1 mm offset. This study highlights the importance of Mobi-C placement to minimize HO formation, preserve motion, and mitigate complications, providing insights for clinical practice.

The Impact of Endplate Coverage on Heterotopic Ossification Following Cervical Disc Replacement: A Systematic Review and Meta Analysis

STUDY DESIGN

Systematic review and meta-analysis.

OBJECTIVE

Describe the impact of endplate coverage on HO in cervical disc replacement (CDR).

SUMMARY OF BACKGROUND DATA

CDR is a motion-sparing alternative to anterior cervical discectomy and fusion. However, the high prevalence of heterotopic ossification threatens to diminish range of motion and limit this benefit associated with CDR.

MATERIALS AND METHODS

A systematic review was performed in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. EMBASE and PubMed databases were queried. Results were deduplicated and screened. Relevant studies were included. All metrics that were reported in ≥3 studies were aggregated for analysis. SPSS was used to perform the meta-analysis.

RESULTS

A total of 10 studies were included in the systematic review. Endplate coverage was assessed using a wide variety of measurements, including anteroposterior implant depth (ID), endplate depth (ED), exposed endplate depth (EED), implant depth to endplate depth ratio (ID:ED), EED to ED ratio (EED:ED), implant width (IW) to endplate width (EW) ratio (IW:EW), and the implant area (IA) to endplate area (EA) ratio (IA:EA). No evidence has linked ID (three studies) to HO. Mixed evidence has linked ID:ED (3/5) and IW:ED (1/2) to HO. All available evidence has linked ED (2), EED (4), EED:ED (2), and IA:EA (1) to HO. In our meta-analysis, ID was not found to be a significant risk factor for HO. However, EED and ID:ED were found to be significant risk factors for HO formation.

CONCLUSIONS

Exposed endplate, especially as assessed by EED and ID:ED, is a significant risk factor for HO. Surgeons should focus on preoperative planning and intraoperative implant selection to maximize endplate coverage. While optimizing technique and implant selection is crucial, improved implant design may also be necessary to ensure that appropriate implant-endplate footprint matching is possible across the anatomic spectrum.

Exploring vertebral bone density changes in a trunk with adolescent idiopathic scoliosis: a mechanobiological modeling investigation of intact and unilaterally paralyzed muscles

This study aimed to elucidate the vertebral bone density variations associated with adolescent idiopathic scoliosis (AIS), specifically examining the impact of unilateral muscle paralysis using an integrated approach combining Frost's Mechanostat theory, a three-dimensional subject-specific finite element model and a musculoskeletal model of the L2 vertebra. The findings revealed a spectrum of bone density values ranging from 0.29 to 0.31 g/cm3, along with vertebral micro-strain levels spanning from 300 to 2200, consistent with existing literature. Furthermore, the ratio of maximum von Mises stress between the concave and convex side in the AIS model with intact muscles was approximately 1.08, which decreased by 4% due following unilateral paralysis of longissimus thoracis pars thoracic muscle. Overall, this investigation contributes to a deeper understanding of AIS biomechanics and lays the groundwork for future research endeavors aimed at optimizing clinical management approaches for individuals with this condition.

Correlation between sagittal morphology of lower lumbar end plate and degenerative changes in patients with lumbar disc herniation

Objective

As an important anatomic factor in the process of lumbar disc herniation (LDH), the correlation between end plate sagittal morphology and intervertebral disc degeneration (IDD) is unclear. Moreover, research on imaging data of lumbar end plate in patients with LDH is still insufficient. Our study aimed to observe the morphological change of the lower lumbar end plate (L3-S1) in patients with LDH on magnetic resonance imaging (MRI) and analyze its correlation with the degree of IDD.

Materials and Methods

A total of 116 patients were included in the study. Based on their MRI, we divided end plates into three types (concave, flat, and irregular), assigned intervertebral discs with Grade I-V given 1-5 points successively according to the Pfirrmann system, and determined whether there was Modic change of each end plate. The correlation between the morphology of the end plate and the degree of IDD was analyzed.

Results

There was an excellent interobserver agreement for each item we analyzed (interclass correlation coefficient >0.75). Concave end plate appeared most frequently (187, 53.7%) and was mainly distributed in L3/4 and L4/5, whereas irregular end plate was the least common type (54, 15.5%) and mainly concentrated in L5/S1. The IDD degree of the corresponding disc increased gradually from concave (3.27 ± 0.81) to irregular end plates (4.25 ± 0.79) (P < 0.05). Irregular end plates were more likely to have Modic changes than concave and flat end plates (P < 0.05).

Conclusion

The sagittal morphology of the lower lumbar end plate is related to modic changes and degree of IDD (based on the Pfirrmann grading system) in patients with LDH, and the concave end plate mostly reflects a lower degree of lumbar disc degeneration, which has substantial clinical significance.

The Influence of Endplate Morphology on Cage Subsidence in Patients With Stand-Alone Oblique Lateral Lumbar Interbody Fusion (OLIF)

STUDY DESIGN

A retrospective study of prospectively collected radiographic and clinical data.

OBJECTIVE

This study aims to investigate the relationship between endplate morphology parameters and the incidence of cage subsidence in patients with mini-open single-level oblique lateral lumbar interbody fusion (OLIF).

METHODS

We included 119 inpatients who underwent OLIF from February 2015 to December 2017. A total of 119 patients with single treatment level of OLIF were included. Plain anteroposterior and lateral radiograph were taken preoperatively, postoperatively, and during follow-up. The correlation between disc height, endplate concave angle/depth, cage position and cage subsidence were investigated. Functional rating index (Visual Analogue Scale for pain, and Roland Morris Disability Questionnaire) were employed to assess clinical outcomes.

RESULTS

Cage subsidence was more commonly seen at the superior endplates (42/119, 35.29%) than at the inferior endplates (6/119, 5.04%) (p < 0.01). More importantly, cage subsidence was significantly less in patients with superior endplates that were without concave angle (3/20, 15%) than with concave angle (37/99, 37.37%) (p < 0.05). Cage subsidence correlated negatively with preoperative anterior disc height (r = -0.21, p < 0.05), but positively with disc distraction rate (r = 0.27, p < 0.01). Lastly, the distance of cage to the anterior edges of the vertebral body showed a positive correlation (r = 0.26, p < 0.01).

CONCLUSIONS

This study for the first time demonstrated that endplate morphology correlates with cage subsidence after OLIF. Since relatively flat endplates with smaller concave angle significantly diminish the incidence of subsidence, the morphology of cage surface should be taken into consideration when designing the next generation of cage. In addition, precise measurement of the disc height to avoid over-distraction, and more anteriorly placement of the cage is suggested to reduce subsidence.

TELD with limited foraminoplasty has potential biomechanical advantages over TELD with large annuloplasty: an in-silico study

Background

Facetectomy, an important procedure in the in–out and out–in techniques of transforaminal endoscopic lumbar discectomy (TELD), is related to the deterioration of the postoperative biomechanical environment and poor prognosis. Facetectomy may be avoided in TELD with large annuloplasty, but iatrogenic injury of the annulus and a high grade of nucleotomy have been reported as risk factors influencing poor prognosis. These risk factors may be alleviated in TELD with limited foraminoplasty, and the grade of facetectomy in this surgery can be reduced by using an endoscopic dynamic drill.

Methods

An intact lumbo-sacral finite element (FE) model and the corresponding model with adjacent segment degeneration were constructed and validated to evaluate the risk of biomechanical deterioration and related postoperative complications of TELD with large annuloplasty and TELD with limited foraminoplasty. Changes in various biomechanical indicators were then computed to evaluate the risk of postoperative complications in the surgical segment.

Results

Compared with the intact FE models, the model of TELD with limited foraminoplasty demonstrated slight biomechanical deterioration, whereas the model of TELD with large annuloplasty revealed obvious biomechanical deterioration. Degenerative changes in adjacent segments magnified, rather than altered, the overall trends of biomechanical change.

Conclusions

TELD with limited foraminoplasty presents potential biomechanical advantages over TELD with large annuloplasty. Iatrogenic injury of the annulus and a high grade of nucleotomy are risk factors for postoperative biomechanical deterioration and complications of the surgical segment.

Effect of heterotopic ossification after bryan-cervical disc arthroplasty on adjacent level range of motion: A finite element study

BACKGROUND

Quantitative bone re-modelling theories suggest that bones adapt to mechanical loading conditions. Follow-up studies have shown that total disc replacement (TDR) modifies stress patterns in the bones, leading to heterotopic ossification (HO). Although there are a few studies on HO using finite element models (FEM), its effect on the adjacent levels and change in range of motion (ROM) have not been adequately investigated. This study interfaces the HO using bone re-modelling algorithm with a finite element solution and investigates the subsequent changes in segmental ROM.

METHODS

A FEM of the human cervical spine (C3-C7) was developed for this study, with material properties obtained from literature. The motion of the segments in the sagittal, frontal and transverse planes under combined loading conditions of 1 Nm moment and 73.6 N compression were validated against experimental corridors. The natural disc between the C5-C6 segment was replaced with the Bryan artificial cervical disc, and changes in sagittal ROM were compared before and after HO. The process of HO was simulated using a bone remodelling algorithm using strain energy density (SED) as the mechanical stimuli.

RESULTS AND CONCLUSION

Our study demonstrates the feasibility of using SED calculations from the flexion-extension loading conditions for prediction of HO after ADR. The current findings suggest that the nature of trabecular stresses, and the subsequent rate and location of HO formation could differ based on the geometric design and nature of constraint for different artificial discs. The Bryan disc significantly reduced ROM at the implanted level in flexion. However, in extension, ROM increased at the implanted level and decreased slightly at the adjacent levels. After HO, ROM drastically reduced at the implanted level in both extension and flexion, and showed a minor increase in the adjacent levels, indicating that biomechanical behavior of high-grade HO is similar to a fused segment, thereby reducing the intended and initial motion preservation.

Numerical simulations of bone remodelling and formation following nucleotomy

Nucleotomy is the gold standard treatment for disc herniation and has proven ability to restore stability by creating a bony bridge without any additional fixation. However, the evolution of mineral density in the extant and new bone after nucleotomy and fixation techniques has to date not been investigated in detail. The main goal of this study is to determine possible mechanisms that may trigger the bone remodelling and formation processes. With that purpose, a finite element model of the L4-L5 spinal segment was used. Bone mineral density (BMD), new tissue composition, and endplate deflection were determined as indicators of lumbar fusion. A bone-remodelling algorithm and a tissue-healing algorithm, both mechanically driven, were implemented to predict vertebral bone alterations and fusion patterns after nucleotomy, internal fixation, and anterior plate placement. When considering an intact disc height, neither nucleotomy nor internal fixation were able to provide the necessary stability to promote bony fusion. However, when 75% of the disc height was considered, bone fusion was predicted for both techniques. By contrast, an anterior plate allowed bone fusion at all disc heights. A 50% disc-height reduction led to osteophyte formation in all cases. Changes in the intervertebral disc tissue caused BMD alterations in the endplates. From this observations it can be drawn that fusion may be self-induced by controlling the mechanical stabilisation without the need of additional fixation. The amount of tissue to be removed to achieve this stabilisation remains to be determined.

Effect of mechanical loading on heterotopic ossification in cervical total disc replacement: a three-dimensional finite element analysis

The development of heterotopic ossification (HO) is considered one of the major complications following cervical total disc replacement (TDR). Even though previous studies have identified clinical and biomechanical conditions that may stimulate HO, the mechanism of HO formation has not been fully elucidated. The objective of this study is to investigate whether mechanical loading is a biomechanical condition that plays a substantial role to decide the HO formation. A finite element model of TDR on the C5-C6 was developed, and HO formation was predicted by simulating a bone adaptation process under various physiological mechanical loadings. The distributions of strain energy on vertebrae were assessed after HO formation. For the compressive force, most of the HO formation occurred on the vertebral endplates uncovered by the implant footplate which was similar to the Type 1 HO. For the anteriorly directed shear force, the HO was predominantly formed in the anterior parts of both the upper and lower vertebrae as the Type 2 HO. For both the flexion and extension moments, the HO shapes were similar to those for the shear force. The total strain energy was reduced after HO formation for all loading conditions. Two distinct types of HO were predicted based on mechanically induced bone adaptation processes, and our findings were consistent with those of previous clinical studies. HO formation might have a role in compensating for the non-uniform strain energy distribution which is one of the mechanical parameters related to the bone remodeling after cervical TDR.

Changes in vertebral strain energy correlate with increased presence of Schmorl's nodes in multi-level lumbar disk degeneration

Patients with skipped-level disk degeneration (SLDD) were recently reported as having a higher prevalence of Schmorl's nodes than patients with contiguous multi-level disk degeneration (CMDD). Fourteen versions of a nonlinear finite element model of a lumbar spine, representing different patterns of single and multi-level disk degeneration, were simulated under physiological loading. Results show that vertebral strain energy is a possible predictor in the development of Schmorl's nodes. The analysis also shows evidence that the development of Schmorl's nodes may be highly dependent on the location of the degeneration disk, with a higher prevalence at superior levels of the lumbar spine.

Parameters that effect spine biomechanics following cervical disc replacement

Total disc replacement (TDR) is expected to provide a more physiologic alternative to fusion. However, long-term clinical data proving the efficacy of the implants is lacking. Limited clinical data suggest somewhat of a disagreement between the in vitro biomechanical studies and in vivo assessments. This conceptual paper presents the potential biomechanical challenges affecting the TDR that should be addressed with a hope to improve the clinical outcomes and our understanding of the devices. Appropriate literature and our own research findings comparing the biomechanics of different disc designs are presented to highlight the need for additional investigations. The biomechanical effects of various surgical procedures are analyzed, reiterating the importance of parameters like preserving uncinate processes, disc placement and its orientation within the cervical spine. Moreover, the need for a 360° dynamic system for disc recipients who may experience whiplash injuries is explored. Probabilistic studies as performed already in the lumbar spine may explore high risk combinations of different parameters and explain the differences between "standard" biomechanical investigations and clinical studies. Development of a patient specific optimized finite element model that takes muscle forces into consideration may help resolve the discrepancies between biomechanics of TDR and the clinical studies. Factors affecting long-term performance such as bone remodeling, subsidence, and wear are elaborated. In vivo assessment of segmental spine motion has been, and continues to be, a challenge. In general, clinical studies while reporting the data have placed lesser emphasis on kinematics following intervertebral disc replacements. Evaluation of in vivo kinematics following TDR to analyze the quality and quantity of motion using stereoradiogrammetric technique may be needed.

The restoration of lumbar intervertebral disc load distribution: a comparison of three nucleus replacement technologies

STUDY DESIGN

A validated L3-L4 nonlinear finite element model was used to evaluate strain and pressure in the surrounding structures for 4 nucleus replacement technologies.

OBJECTIVE

The objective of the current study was to compare subsidence and anular damage potential between 4 current nucleus replacement technologies. It was hypothesized that a fully conforming nucleus replacement would minimize the risk of both subsidence and anular damage.

SUMMARY OF BACKGROUND DATA

Nucleus pulposus replacements are emerging as a less invasive alternative to total disc replacement and fusion as a solution to degenerative intervertebral discs. Multiple technologies have been developed and are currently undergoing clinical investigation.

METHODS

The testing conditions were applied by excavating the nucleus of the intact model and virtually implanting models representing the various nucleus replacement technologies. The implants consisted of a conforming injectable polyurethane (E = 4 MPa), soft hydrogel (E = 4 MPa), stiff hydrogel (E = 20 MPa), and polyether-etherketone (PEEK) on PEEK articulating designs. The model was exercised in flexion, extension, lateral bending, axial rotation (7.5 Nm with 450 N preload), and compression (1000 N). Vertebral body strain, anular maximum shear strain, endplate contact pressure, anulus-implant contact pressure, and bone remodeling stimulus were reported.

RESULTS

The PEEK implant induced strain maxima in the vertebral bodies with associated endplate contact pressure concentrations. For the PEEK and hydrogel implants, areas of nonconformity with the endplate indicated adjacent bone resorption. Lack of conformity between the implant and inner anulus for the PEEK and hydrogel implants resulted in inward anular bulging with associated increased maximum shear strain. The conforming polyurethane implant maintained outward bulging of the inner anular wall and indicated no bone resorption or stress shielding adjacent to the implant.

CONCLUSION

A fully conforming nucleus replacement resulted in a decreased propensity for subsidence, anular bulging, and further degeneration of the anulus when compared with nonconforming implants.