Computer simulation in spine biomechanics : An application in adolescent idiopathic scoliosis

Scoliosis is a spinal deformity, involving a side-to-side curvature of the spine in the coronal plane as well as a rotation of the spinal column in the transverse plane. The coronal curvature is measured using a Cobb angle. If the deformity is severe, treatment for scoliosis may require surgical intervention whereby a rod is attached to the spinal column to correct the abnormal curvature. In order to provide surgeons with an improved ability to predict the likely outcomes following surgery, techniques to create patient-specific finite element models (FEM) of scoliosis patients treated at the Mater Children’s Hospital (MCH) in Brisbane are being developed and validated. This paper presents a comparison of the simulated and clinical data for a scoliosis patient treated at MCH.

Understanding how axial loads on the spine influence segmental biomechanics for idiopathic scoliosis patients: A magnetic resonance imaging study

Background Segmental biomechanics of the scoliotic spine are important since the overall spinal deformity is comprised of the cumulative coronal and axial rotations of individual joints. This study investigates the coronal plane segmental biomechanics for adolescent idiopathic scoliosis patients in response to physiologically relevant axial compression. Methods Individual spinal joint compliance in the coronal plane was measured for a series of 15 idiopathic scoliosis patients using axially loaded magnetic resonance imaging. Each patient was first imaged in the supine position with no axial load, and then again following application of an axial compressive load. Coronal plane disc wedge angles in the unloaded and loaded configurations were measured. Joint moments exerted by the axial compressive load were used to derive estimates of individual joint compliance. Findings The mean standing major Cobb angle for this patient series was 46°. Mean intra-observer measurement error for endplate inclination was 1.6°. Following loading, initially highly wedged discs demonstrated a smaller change in wedge angle, than less wedged discs for certain spinal levels (+ 2,+1,− 2 relative to the apex, (p < 0.05)). Highly wedged discs were observed near the apex of the curve, which corresponded to lower joint compliance in the apical region. Interpretation While individual patients exhibit substantial variability in disc wedge angles and joint compliance, overall there is a pattern of increased disc wedging near the curve apex, and reduced joint compliance in this region. Approaches such as this can provide valuable biomechanical data on in vivo spinal biomechanics of the scoliotic spine, for analysis of deformity progression and surgical planning.

Anterior vertebral stapling for the fusionless correction of scoliosis

Fusionless scoliosis surgery is an emerging treatment for idiopathic scoliosis as it offers theoretical advantages over current forms of treatment. Currently the treatment options for idiopathic scoliosis are observation, bracing and fusion. While brace treatment is non-invasive, and preserves the growth, motion, and function of the spine, it does not correct deformity and is only modestly successful in preventing curve progression. In adolescents who fail brace treatment, surgical treatment with an instrumented spinal fusion usually results in better deformity correction but is associated with substantially greater risk. Furthermore in younger patients requiring surgical treatment, fusion procedures are known to adversely effect the future growth of the chest and spine. Fusionless treatments have been developed to allow effective surgical treatment of patients with idiopathic scoliosis who are too young for fusion procedures. Anterior vertebral stapling is one such fusionless treatment which aims to modulate the growth of vertebra to allow correction of scoliosis whilst maintaining normal spinal motion The Mater Misericordiae Hospital in Brisbane has begun to use anterior vertebral stapling to treat patients with idiopathic scoliosis who are too young for fusion procedures. Currently the only staple approved for clinical use is manufactured by Medtronic Sofamor Danek (Memphis, TN). This thesis explains the biomechanical and anatomical changes that occur following anterior vertebral staple insertion using in vitro experiments performed on an immature bovine model. Currently there is a paucity of published information about anterior vertebral stapling so it is hoped that this project will provide information that will aid in our understanding of the clinical effects of staple insertion. The aims of this experimental study were threefold. The first phase was designed to determine the changes in the bending stiffness of the spine following staple insertion. The second phase was designed to measure the forces experienced by the staple during spinal movements. The third and final phase of testing was designed to describe the structural changes that occur to a vertebra as a consequence of staple insertion. The first phase of testing utilised a displacement controlled testing robot to compare the change in stiffness of a single spinal motion segment following staple insertion for the three basic spinal motions of flexion-extension, lateral bending, and axial rotation. For the second phase of testing strain gauges were attached to staples and used to measure staple forces during spinal movement. In the third and final phase the staples were removed and a testing specimen underwent micro-computed tomography (CT) scanning to describe the anatomical changes that occur following staple insertion. The displacement controlled testing showed that there was a significant decrease in bending stiffness in flexion, extension, lateral bending away from the staple, and axial rotation away from the staple following staple insertion. The strain gauge measurements showed that the greatest staple forces occurred in flexion and the least in extension. In addition, a reduction in the baseline staple compressive force was seen with successive loading cycles. Micro-CT scanning demonstrated that significant damage to the vertebral body and endplate occurred as a consequence of staple insertion. The clinical implications of this study are significant. Based on the findings of this project it is likely that the clinical effect of the anterior vertebral staple evaluated in this project is a consequence of growth plate damage (also called hemiepiphysiodesis) causing a partial growth arrest of the vertebra rather than simply compression of the growth plate. The surgical creation of a unilateral growth arrest is a well established treatment used in the management of congenital scoliosis but has not previously been considered for use in idiopathic scoliosis.

A biomechanical analysis of anterior instrumentation used in the surgical correction of adolescent idiopathic scoliosis

Scoliosis is a deformity of the spine which affects children and adolescents, and remains a challenge to treat. This study measured the forces used during surgery to correct scoliosis and studied changes to spinal mechanics from the implantation of metal rods used to hold the spine straight. The results of this study will help surgeons and engineers understand how to straighten the spine more efficiently to provide patients with better outcomes.

Quantification of prone thoracic manipulation using inertial sensor–derived accelerations

Objective The aim of this study was to determine the linear acceleration, time-to-peak acceleration, and effect of hand position comparing 2 clinicians completing a thoracic manipulation. Methods Thirteen volunteers received a right- and left-“handed” prone thoracic manipulation while accelerations were recorded by an inertial sensor. Peak thrust acceleration and time-to-peak thrust were measured. Results There were differences in thrust acceleration between right- and left-handed techniques for one therapist. The mean peak thrust acceleration was different between therapists, with the more practiced therapist demonstrating greater peak thrust accelerations. Time-to-peak acceleration also revealed between therapist differences, with the more practiced therapist demonstrating shorter time-to-peak acceleration. Cavitation data suggested that manipulations with greater accelerations were more likely to result in cavitation. Conclusion The results of this study suggest that with greater frequency of use, therapists are likely to achieve greater accelerations and shorter time-to-peak accelerations. Furthermore, this study showed that an inertial sensor can be used to quantify important variables during thoracic manipulation and are able to detect intertherapist differences in technique.

Patient-specific computational biomechanics for simulating adolescent scoliosis surgery : predicted vs clinical correction for a preliminary series of six patients

Scoliosis is a spinal deformity that requires surgical correction in progressive cases. In order to optimize surgical outcomes, patient-specific finite element models are being developed by our group. In this paper, a single rod anterior correction procedure is simulated for a group of six scoliosis patients. For each patient, personalised model geometry was derived from low-dose CT scans, and clinically measured intra-operative corrective forces were applied. However, tissue material properties were not patient-specific, being derived from existing literature. Clinically, the patient group had a mean initial Cobb angle of 47.3 degrees, which was corrected to 17.5 degrees after surgery. The mean simulated post-operative Cobb angle for the group was 18.1 degrees. Although this represents good agreement between clinical and simulated corrections, the discrepancy between clinical and simulated Cobb angle for individual patients varied between -10.3 and +8.6 degrees, with only three of the six patients matching the clinical result to within accepted Cobb measurement error of +-5 degrees. The results of this study suggest that spinal tissue material properties play an important role in governing the correction obtained during surgery, and that patient-specific modelling approaches must address the question of how to prescribe patient-specific soft tissue properties for spine surgery simulation.

Effects of surgical joint destabilization on load sharing between ligamentous structures in the thoracic spine : a finite element investigation

Background: In vitro investigations have demonstrated the importance of the ribcage in stabilising the thoracic spine. Surgical alterations of the ribcage may change load-sharing patterns in the thoracic spine. Computer models are used in this study to explore the effect of surgical disruption of the rib-vertebrae connections on ligament load-sharing in the thoracic spine. Methods: A finite element model of a T7-8 motion segment, including the T8 rib, was developed using CT-derived spinal anatomy for the Visible Woman. Both the intact motion segment and the motion segment with four successive stages of destabilization (discectomy and removal of right costovertebral joint, right costotransverse joint and left costovertebral joint) were analysed for a 2000Nmm moment in flexion/extension, lateral bending and axial rotation. Joint rotational moments were compared with existing in vitro data and a detailed investigation of the load sharing between the posterior ligaments carried out. Findings: The simulated motion segment demonstrated acceptable agreement with in vitro data at all stages of destabilization. Under lateral bending and axial rotation, the costovertebral joints were of critical importance in resisting applied moments. In comparison to the intact joint, anterior destabilization increases the total moment contributed by the posterior ligaments. Interpretation: Surgical removal of the costovertebral joints may lead to excessive rotational motion in a spinal joint, increasing the risk of overload and damage to the remaining ligaments. The findings of this study are particularly relevant for surgical procedures involving rib head resection, such as some techniques for scoliosis deformity correction.

Geometric sensitivity of patient-specific finite element models of the spine to variability in user-selected anatomical landmarks

Software to create individualised finite element (FE) models of the osseoligamentous spine using pre-operative computed tomography (CT) data-sets for spinal surgery patients has recently been developed. This study presents a geometric sensitivity analysis of this software to assess the effect of intra-observer variability in user-selected anatomical landmarks. User-selected landmarks on the osseous anatomy were defined from CT data-sets for three scoliosis patients and these landmarks were used to reconstruct patient-specific anatomy of the spine and ribcage using parametric descriptions. The intra-observer errors in landmark co-ordinates for these anatomical landmarks were calculated. FE models of the spine and ribcage were created using the reconstructed anatomy for each patient and these models were analysed for a loadcase simulating clinical flexibility assessment. The intra-observer error in the anatomical measurements was low in comparison to the initial dimensions, with the exception of the angular measurements for disc wedge and zygapophyseal joint (z-joint) orientation and disc height. This variability suggested that CT resolution may influence such angular measurements, particularly for small anatomical features, such as the z-joints, and may also affect disc height. The results of the FE analysis showed low variation in the model predictions for spinal curvature with the mean intra-observer variability substantially less than the accepted error in clinical measurement. These findings demonstrate that intra-observer variability in landmark point selection has minimal effect on the subsequent FE predictions for a clinical loadcase.

Segmental torso masses and joint torques produced by gravity in the adolescent scoliotic spine

Introduction Calculating segmental torso masses in Adolescent Idiopathic Scoliosis (AIS) patients allows the gravitational loading on the scoliotic spine during relaxed standing to be estimated. Methods Low dose CT data was used to calculate vertebral level-by-level torso masses and spinal joint torques for 20 female AIS patients (mean age 15.0 ± 2.7 years, mean Cobb angle 53 ± 7.1°). ImageJ software (v1.45 NIH USA) was used to threshold the T1 to L5 CT images and calculate the segmental torso volume and mass for each vertebral level. Masses for the head, neck and arms were taken from published data. Intervertebral joint torques in the coronal and sagittal planes at each vertebral level were found from the position of the centroid of the segment masses relative to the joint centres (assumed to be at the centre of the intervertebral disc. The joint torque at each level was found by summing torque contributions for all segments above that joint. Results Segmental torso mass increased from 0.6kg at T1 to 1.5kg at L5. The coronal plane joint torques due to gravity were 5-7Nm at the apex of the curve; sagittal torques were 3-5.4Nm. Conclusion CT scans were in the supine position and curve magnitudes are known to be smaller than those in standing. Hence, this study has shown that gravity produces joint torques potentially of higher than 7Nm in the coronal plane and 5Nm in the sagittal plane during relaxed standing in scoliosis patients. The magnitude of these torques may help to explain the mechanics of AIS progression and the mechanics of bracing. This new data on torso segmental mass in AIS patients will assist biomechanical models of scoliosis.

The effect of intervertebral staple insertion on bovine spine segment stiffness

Introduction There is growing interest in the biomechanics of ‘fusionless’ implant constructs used for deformity correction in the thoracic spine. Intervertebral stapling is a leading method of fusionless corrective surgery. Although used for a number of years, there is limited evidence as to the effect these staples have on the stiffness of the functional spinal unit. Materials and Methods Thoracic spines from 6-8 week old calves were dissected and divided into motion segments including levels T4-T11 (n=14). Each segment was potted in polymethylemethacrylate. An Instron Biaxial materials testing machine with a custom made jig was used for testing. The segments were tested in flexion/extension, lateral bending and axial rotation at 37⁰C and 100% humidity, using moment control to a maximum 1.75 Nm with a loading rate of 0.3 Nm per second. This torque was found sufficient to achieve physiologically representative ranges of movement. The segments were initially tested uninstrumented with data collected from the tenth load cycle. Next a left anterolateral Shape Memory Alloy (SMA) staple was inserted (Medtronic Sofamor Danek, USA). Biomechanical testing was repeated as before with data collected from the tenth load cycle. Results In flexion/extension there was an insignificant drop in stiffness of 3% (p=0.478). In lateral bending there was a significant drop in stiffness of 21% (p<0.001). This was mainly in lateral bending away from the staple, where the stiffness reduced by 30% (p<0.001). This was in contrast to lateral bending towards the staple where it dropped by 12% which was still statistically significant (p=0.036). In axial rotation there was an overall near significant drop in stiffness of 11% (p=0.076). However, this was more towards the side of the staple measuring a decrease of 14% as opposed to 8% away from the staple. In both cases it was a statistically insignificant drop (p=0.134 and p=0.352 respectively). Conclusion Insertion of intervertebral SMA staples results in a significant reduction in motion segment stiffness in lateral bending especially in the direction away from the staple. The staple had less effect on axial rotation stiffness and minimal effect on flexion/extension stiffness.

Gravity-induced coronal plane joint moments in the adolescent scoliotic spine

INTRODUCTION Calculating segmental (vertebral level-by-level) torso masses in Adolescent Idiopathic Scoliosis (AIS) patients allows the gravitational loading on the scoliotic spine during relaxed standing to be estimated. METHODS Existing low dose CT scans were used to calculate vertebral level-by-level torso masses and joint moments occurring in the spine for a group of female AIS patients with right-sided thoracic curves. Image processing software, ImageJ (v1.45 NIH USA) was used to reconstruct the torso segments and subsequently measure the torso volume and mass corresponding to each vertebral level. Body segment masses for the head, neck and arms were taken from published anthropometric data. Intervertebral joint moments at each vertebral level were found by summing each of the torso segment masses above the required joint and multiplying it by the perpendicular distance to the centre of the disc. RESULTS AND DISCUSSION Twenty patients were included in this study with a mean age of 15.0±2.7 years and a mean Cobb angle 52±5.9°. The mean total trunk mass, as a percentage of total body mass, was 27.8 (SD 0.5) %. Mean segmental torso mass increased inferiorly from 0.6kg at T1 to 1.5kg at L5. The coronal plane joint moments during relaxed standing were typically 5-7Nm at the apex of the curve (Figure 1), with the highest apex joint of 7Nm. CT scans were performed in the supine position and curve magnitudes are known to be 7-10° smaller than those measured in standing [1]. Therefore joint moments produced by gravity will be greater than those calculated here. CONCLUSIONS Coronal plane joint moments as high as 7Nm can occur during relaxed standing in scoliosis patients, which may help to explain the mechanics of AIS progression. The body mass distributions calculated in this study can be used to estimate joint moments derived using other imaging modalities such as MRI and subsequently determine if a relationship exists between joint moments and progressive vertebral deformity.

Using patient specific finite element models to understand the biomechanics of paediatric spinal deformity surgery

Adolescent idiopathic scoliosis (AIS) is a spinal deformity, which may require surgical correction by attaching rods to the patient’s spine using screws inserted into the vertebrae. Complication rates for deformity correction surgery are unacceptably high. Determining an achievable correction without overloading the adjacent spinal tissues or implants requires an understanding of the mechanical interaction between these components. We have developed novel patient specific modelling software to create individualized finite element models (FEM) representing the thoracolumbar spine and ribcage of scoliosis patients. We are using these models to better understand the biomechanics of spinal deformity correction.

Comparison of eight published static finite element models of the intact lumbar spine : predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values. Eight well-established FE models of the lumbar spine (L1-5) of different research centres around the globe were subjected to pure and combined loading modes and compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges, and their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with published median in vitro values. However, the ranges of predictions were larger and exceeded those reported in vitro, especially for the facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with measured in vivo values. In light of high inter-subject variability, the generalization of results of a single model to a population remains a concern. This study demonstrated that the pooled median of individual model results, similar to a probabilistic approach, can be used as an improved predictive tool in order to estimate the response of the lumbar spine.