Biomechanical properties of a novel nonfusion artificial vertebral body for anterior lumbar vertebra resection and internal fixation

The aim of the study was to evaluate the biomechanical properties of a novel nonfused artificial vertebral body in treating lumbar diseases and to compare with those of the fusion artificial vertebral body. An intact finite element model of the L1–L5 lumbar spine was constructed and validated. Then, the finite element models of the fusion group and nonfusion group were constructed by replacing the L3 vertebral body and adjacent intervertebral discs with prostheses. For all finite element models, an axial preload of 500 N and another 10 N m imposed on the superior surface of L1. The range of motion and stress peaks in the adjacent discs, endplates, and facet joints were compared among the three groups. The ranges of motion of the L1–2 and L4–5 discs in flexion, extension, left lateral bending, right lateral bending, left rotation and right rotation were greater in the fusion group than those in the intact group and nonfusion group. The fusion group induced the greatest stress peaks in the adjacent discs and adjacent facet joints compared to the intact group and nonfusion group. The nonfused artificial vertebral body could better retain mobility of the surgical site after implantation (3.6°–8.7°), avoid increased mobility and stress of the adjacent discs and facet joints.

www.nature.com/scientificreports/ and stability of the novel prosthesis through in vitro mechanical tests with human specimens, and the results showed that the novel prosthesis could not only replace vertebral lesions but also reconstruct the stability and mobility of the surgical site 10,11 . Although we found that the new prosthesis could reduce the activity of the adjacent intervertebral space compared with the fused artificial vertebral body, the stress effect on the surrounding tissues needed further study. Therefore, finite element (FE) analysis was used to study the biomechanical properties of this new prosthesis to provide a reference for its long-term biomechanical safety after implantation.
Disc ROM. The ROMs of the L1-2, L2-3, L3-4 and L4-5 discs in flexion, extension, left lateral bending, right lateral bending, left rotation and right rotation among the different groups are shown in Fig. 2. The ROMs of the discs adjacent to the surgical site (L1-2 and L4-5) in all the above directions were obviously greater in the fusion group than in the intact group. However, the ROMs of the discs in the surgical site (L2-3 and L3-4) were significantly decreased in the fusion group compared with those in the intact group. The ROMs of the discs in the nonfusion group showed an opposite trend to those in the fusion group. The ROMs of the L1-2 and L4-5 discs in flexion, extension, left lateral bending, right lateral bending, left rotation and right rotation were obvi-   www.nature.com/scientificreports/

Discussion
Anterior lumbar subtotal resection is widely used to treat lumbar tuberculosis, fracture, tumors and other diseases and often uses titanium cages combined with titanium plates or screw rods. Although this operation can achieve full decompression of the spinal cord and reconstruct the height and stability of the vertebra under direct vision, it often requires the fusion of three or more vertebrae, resulting in loss of mobility in the surgical area. Long-term follow-up results showed that the incidence of complications such as adjacent disc degeneration was significantly increased after this procedure 14 . The appearance of artificial intervertebral discs changed the fusion concept and integrated the concept of mobility into the spinal field. The long-term follow-up results showed that the artificial intervertebral disc had a significant effect on recovery of the height of the intervertebral space and retention of spinal motor function. However, artificial intervertebral discs are difficult to use in patients with vertebrotomy due to their inherent structural limitations. Therefore, there is an urgent need to design a new type of prosthesis that can not only preserve the mobility of the surgical site but also restore the height and stability of the lesion site. In the early stage, a new type of lumbar prosthesis was designed according to the anatomical parameters of lumbar vertebrae -the MALV. We tested the mobility and stability of the MALV with an in vitro mechanics test 11 , and the results showed that the new prosthesis could not only replace the vertebral body of the lesion but also reconstruct the stability and mobility of the surgical site. However, the stress effect of the new prosthesis on the surrounding tissues after implantation still needed to be further studied. The FE method is a kind of discrete method in numerical calculations and is based on the matrix method in structural mechanics, elastic mechanics and other fields. Belytschko et al. 15 first applied FE analysis to the study of spine biomechanics in 1973. Subsequently, many scholars verified the reliability of the stress and strain results of FE analysis through relevant experiments and obtained relatively positive results [16][17][18][19][20] . In this study, a FE model of the whole lumbar spine was established in the intact group using the CT data of a healthy adult and validated using published experimental data. On this basis, FE models of the fusion and nonfusion groups were successfully established. The ROMs were evaluated in the intact model and models with prostheses. The stresses exerted on the adjacent discs, endplates and facet joints were analyzed.
According to our results, the ROMs of the L2-3 and L3-4 intervertebral discs in the fusion group were significantly lower than those in the intact group in the directions of flexion, extension, lateral bending and rotation, while the ROMs of the L2-3 and L3-4 intervertebral discs in the nonfusion group were higher than those in the intact group. In the fusion group, the ROMs of the adjacent discs (L1-2 disc and L4-5 disc) in the above directions were significantly increased compared with those in the intact group and nonfusion group, while the ROMs of the adjacent discs in the nonfusion group in the flexion, extension, lateral bending and rotation directions were decreased compared with those in the intact group. This is consistent with our previous results from in vitro mechanics tests on fresh cadaver specimens 11 . It was suggested that fusion surgery would seriously affect the physiological motor function of the lumbar spine at the surgical site and increase the intervertebral disc ROM in the adjacent segments, which was consistent with the previous research results of Biswas, JK et al 21 . Although the ROMs of the adjacent discs to the surgical site were reduced in the nonfusion group compared with those in the intact group, the difference between the two groups was not obvious. Therefore, after implantation of the new prosthesis, not only can the mobility of the surgical site be better retained but the increased mobility of the adjacent discs caused by fusion can also be avoided to some extent.
The maximum von Mises stresses of the adjacent discs (L1-2 disc and L4-5 disc) and facet joints (J12 and J45) in the fusion group in the directions of flexion, extension, lateral bending and rotation were significantly increased compared with those in the intact group and the nonfusion group, while the maximum von Mises stresses of the adjacent discs (L1-2 disc and L4-5 disc) and facet joints (J12 and J45) in the nonfusion group were decreased compared with those in the intact group. Thus, to some extent, the stress in the adjacent segments will increase after fusion, leading to accelerated degeneration, which is consistent with the increased stress in the The results of the maximum von Mises stresses on the endplate of the surgical site (lower endplate of L2 and upper endplate of L4) showed that the maximum von Mises stresses in the fusion group and nonfusion group in the flexion, extension, lateral bending and rotation directions were significantly higher than those in the intact group. The reason is that the intervertebral discs in the intact group have a certain elastic cushioning, which can reduce the stress to the adjacent endplates to a certain extent. However, the prostheses in the fusion group and nonfusion group were made of metal materials, which lack an elastic buffer, resulting in a significant stress increase in the adjacent endplate. However, the maximum von Mises stresses of the adjacent endplates in the above directions in the nonfusion group were much higher than those in the fusion group. Due to the use of titanium plate or pin rod auxiliary fixation, the load of the titanium cage can be shared to a certain extent, reducing the stress in the contact area between the titanium cage and the end plate. Therefore, to avoid the high sinking rate after traditional titanium cage implantation, the contact area between the new prosthesis and the endplate was carefully designed to not only increase the contact area between the new prosthesis and the endplate but also match the anatomical morphology of the adjacent endplate.
In conclusion, this study confirmed that this new type of lumbar prosthesis could better retain the mobility of the surgical site after implantation, avoid increases in mobility and stress of the adjacent discs and stress of the facet joints, and reduce the incidence of adjacent segment degeneration after long-term implantation. However, the long-term stability, fatigue resistance, wear resistance and other mechanical properties of the new prosthesis after implantation still need to be further studied.

Materials and methods
Development of MALV. The MALV prosthesis consists of three parts: vertebral body part, intervertebral disc part and composite material ball (Fig. 5). The vertebral body part is an irregular cylinder with a depression in the back, grooves on the upper and lower ends and a fixed column in the center. The intervertebral disc part consists of a sunken back plate, curved side plate, and cylindrical protuberances with ends connected to the spherical shell structure. The composite material ball is composed of Ultrahigh molecular weight polyethylene (UHMWPE) that forms a ball and socket joint with a spherical shell structure. The vertebral body part and intervertebral disc part were made of Ti6Al4V.  www.nature.com/scientificreports/ model of the L1-L5 vertebrae (Fig. 6a). The model was saved as an STL file and imported into Geomagic Studio 12.0 (Raindrop Inc., USA) to edit and de-noise the triangular surfaces. The triangular surfaces enveloping each vertebra were ensured to be spatially closed, and finally, the spatially triangular surfaces were fitted into spatially closed NURBS and saved in IGES format (Fig. 6b). Models of the bone, cartilage endplates, annulus fibrosus, nucleus pulposus, and ligaments were constructed using Pro/Engineer5.0 (PTC Inc., Massachusetts, USA) (Fig. 6c). Abaqus (Hibbitt, Karlsson, and Sorensen, Inc., Providence, Rhode Island, USA) was used to define material properties and for FE analysis (Fig. 6d).
The vertebral body consisted of cortical bone and cancellous bone. The thickness of the cortical shell was 1 mm. The thickness of the endplate was set at 0.5 mm 22 . The intervertebral disc was divided into the nucleus pulposus and annulus fibrosus. The nucleus pulposus accounted for 30% ~ 40% of the intervertebral volume 22 and was modeled as a linearly elastic fluid element. The annulus fibrosus consisted of the annulus ground substance and fibers. Six layers of annulus fibers were embedded into the annulus ground substance at an inclination of ± 30°. The elastic strength of the annulus fibers proportionally increased and varied from the innermost layer (360 MPa) to the outermost layer (550 MPa). The articular cartilage of the facet joints was set as 0.5 mm 23 . A total of seven ligaments were modeled, including the anterior longitudinal ligament, posterior longitudinal ligament, ligamentum flavum, capsular ligament, interspinous ligament, supraspinous ligament, and intertransverse ligament. The element types and material properties used in the FE model were defined according to previous reports [22][23][24][25] and are shown in Table 1.

FE modeling after implant implantation.
The FE models after prothesis implantation are shown in Fig. 7. The modeling processes were as follows. First, the L3 vertebral body and adjacent discs (L2-3 disc, L3-4 disc) were removed. The anterior longitudinal ligament and posterior longitudinal ligament of the surgical site were also removed. However, the posterior vertebral body parts, such as the pedicle and facet joints, were retained. Second, a titanium cage filled with cancellous bone was implanted at the surgical site, a titanium plate and screw were used to assist the lateral fixation, and then the FE model for the fusion group was constructed. The inner diameter of the titanium cage was 20 mm, the height was 46 mm, the thickness was 2 mm, the titanium plate length was 100 mm, the width was 23 mm, and the thickness was 3 mm. The screw was 6.5 mm in diameter and 45 mm in length. Third, the MALV was implanted into the surgical site and fixed with the adjacent vertebral body with four screws to construct the FE model of the nonfusion group. The diameter and length of the screws used were 6.5 mm and 45 mm, respectively. For all FE models above, geometric matching at the prosthesis-endplate interface was achieved using the "Boolean calculation" to remove the portion of the cage and MALV that overlapped with the vertebral body. The effect of the teeth on the titanium cage and MALV surface was minimized by assigning a friction coefficient of 0.2 to the prosthesis-endplate interface 26 . A "tie" constraint was assigned to the interfaces between the screw plate and screw bone to simulate rigid fixation. The contact surface between the ball and socket joint and MALV was smooth.
Boundary and loading conditions. For all FE models, the lower surface of the L5 vertebra was constrained. An axial preload of 500 N was imposed on the superior surface of L1 to simulate the corresponding physiological compression. Another 10 Nm was applied on L1 to simulate flexion (FLX), extension (EXT), left lateral bending (LLB), right lateral bending (RLB), left rotation (LT) and right rotation (RT). To validate the L1-L5 FE model of the intact group, the segmental ranges of motion (ROMs) (L1-L2, L2-L3, L3-L4 and L4-L5)    www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.