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.

Comparison of eight published lumbar spine finite element models

Due to its ability to represent intricate systems with material nonlinearities as well as irregular loading, boundary, geometrical and material domains, the finite element (FE) method has been recognized as an important computational tool in spinal biomechanics. Current FE models generally account for a single distinct spinal geometry with one set of material properties despite inherently large inter-subject variability. The uncertainty and high variability in tissue material properties, geometry, loading and boundary conditions has cast doubt on the reliability of their predictions and comparability with reported in vitro and in vivo values. A multicenter study was undertaken to compare the results of eight well-established models of the lumbar spine that have been developed, validated and applied for many years. Models were subjected to pure and combined loading modes and their predictions were 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; 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 previously published median in vitro values. However, the ranges of predictions were larger and exceeded the in vitro ranges, especially for facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with in vivo values. The simulations yielded median facet joint forces of 0 N in flexion, 38 N in extension, 14 N in lateral bending and 60 N in axial rotation that could not be validated due to the paucity of in vivo facet joint forces. In light of high inter-subject variability, one must be cautious when generalizing predictions obtained from one deterministic model. This study demonstrates however that the predictive power increases when FE models are combined together. The median of individual numerical results can hence be used as an improved tool in order to estimate the response of the lumbar spine.

Medical imaging and Biomechanical analysis of scoliosis progression in the growing adolescent spine

Progression of spinal deformity in children was studied with Computed Tomography (CT) and Magnetic Resonance Imaging (MRI) to identify how gravity affects the deformity and to determine the full three-dimensional character of the deformity. The CT study showed that gravity is significant in deformity progression in some patients which has implications for clinical patient management. The world first MRI study showed that the standard clinical measure used to define the extent of the deformity is inadequate and further use of three-dimensional MRI should be considered by spinal surgeons.

Development of a new ultrasound-based system for tracking motion of the human lumbar spine: Reliability, stability and repeatability during forward bending movement trials

The aim of this study was to develop a new method for quantifying intersegmental motion of the spine in an instrumented motion segment L4–L5 model using ultrasound image post-processing combined with an electromagnetic device. A prospective test–retest design was employed, combined with an evaluation of stability and within- and between-day intra-tester reliability during forward bending by 15 healthy male patients. The accuracy of the measurement system using the model was calculated to be ± 0.9° (standard deviation = 0.43) over a 40° range and ± 0.4 cm (standard deviation = 0.28) over 1.5 cm. The mean composite range of forward bending was 15.5 ± 2.04° during a single trial (standard error of the mean = 0.54, coefficient of variation = 4.18). Reliability (intra-class correlation coefficient = 2.1) was found to be excellent for both within-day measures (0.995–0.999) and between-day measures (0.996–0.999). Further work is necessary to explore the use of this approach in the evaluation of biomechanics, clinical assessments and interventions.

Using axially loaded MRI to investigate the biomechanics of adolescent idiopathic scoliosis

Introduction. Spinal flexibility measurement is an important aspect of pre-operative clinical assessment in the treatment of Adolescent Idiopathic Scoliosis (AIS). Clinically, curve flexibility is a combined measure for all vertebral levels. We propose that in vivo flexibility for individual spinal joints could provide valuable additional information in planning treatment for scoliosis. Methods. Individual spinal joint flexibility in the coronal plane was measured for a series of AIS patients using axially loaded magnetic resonance imaging. Each patient underwent magnetic resonance imaging in the supine position, with no axial load, and then following application of an axial compressive load equal to half the patient’s bodyweight. 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. Results. Fifteen AIS patients were included in the study (mean clinical Cobb angle 46 degrees, mean age 15.3 years). Mean intra-observer measurement error for endplate inclination was 1.6˚. The mean increase in measured major Cobb angle between unloaded and loaded scans was 7.6˚. For certain spinal levels (+2,+1,-2 relative to the apex) there was a statistically significant relationship between change in wedge angle under load and initial wedge angle, such that initially highly wedged discs demonstrated a smaller change in wedge angle than less wedged discs. Highly wedged discs were observed near the apex of the curve, which corresponded to lower joint compliance in the apical region. Conclusion. Approaches such as this can provide valuable biomechanical data on in vivo spinal biomechanics in AIS. Knowledge of individual joint flexibility may assist surgeons to determine which spinal procedure is most appropriate for a patient, which levels should be included in a spinal fusion and the relative mobility of individual joints in the deformed region of the spine.

Sequential MRI reveals individual level deformities in the growing scoliotic spine that are not seen clinically by the cobb angle

Clinically, the Cobb angle method measures the overall scoliotic curve in the coronal plane but does not measure individual vertebra and disc wedging. The contributions of the vertebrae and discs in the growing scoliotic spine were measured to investigate coronal plane deformity progression with growth. Sequential MRI data in this project showed complex patterns of deformity progression. Changes to the wedging of individual vertebrae and discs may occur in patients who have no increase in Cobb angle measure; the Cobb method alone may be insufficient to capture the complex mechanisms of deformity progression.

Influence of posture on the range of axial rotation and coupled lateral flexion of the thoracic spine

Objective
This study examines the influence of posture on the range of axial rotation of the thorax and the range and direction of the coupled lateral flexion.

Methods
The ranges of mid thoracic axial rotation and coupled lateral flexion were measured in 52 asymptomatic subjects (aged 18-43 years) using an optical motion analysis system. To examine the influence of posture on primary and coupled motion, we initiated axial rotation from a neutral sitting posture and from end-range thoracic flexion and extension.

Results
There was a significant decrease in the range of thoracic rotation in flexion compared with the neutral and extended postures (P < .001). The mean range of coupled lateral flexion was 8.9% of the axial rotation range in the neutral posture and increased to 14.3% and 23.2% in the extended and flexed postures, respectively. Patterns of coupled motion varied between subjects, but an ipsilateral pattern was more common in the flexed posture, whereas a contralateral pattern was more common in the neutral and extended postures.

Conclusions
The ranges and patterns of coupled motion of the thorax appear to be strongly influenced by the posture from which the movement is initiated. This has important implications in relation to the interpretation of clinical tests of thoracic motion and in consideration of mechanisms of development of thoracic pain disorders.

Virtual lumbar spine of multi-body model based on simbody

The spine is an important and complex skeletal structure in the human body. It is a vulnerable part of our skeleton that is open to many medical problems. Hence it is necessary to establish a virtual spine model to assist surgeons to understand biomechanics of the spine. In this study, we aim to propose a virtual spine multi-body model. The computational biomechanical modeling of the spine is based on the theory of multi-body dynamics and implemented with SimBody open-source SDK. Simbody was then used to solve the kinetic equations and simulate the movement of spine. The spine model was validated by comparing its simulation results with experimental results from literature. The spine model will be helpful to understand biomechanics of the spine under loading.