Polycaprolactone-based scaffold plus rhBMP-2 in a sheep thoracic spine fusion model

We report the application of a novel scaffold design in a sheep thoracic spine model for spine deformity correction. The combination of the calcium-phosphate coated polycaprolactone scaffolds with recombinant human bone morphogenic protein-2 (rhBMP-2) are intended as a future bone graft substitute in ensuring the stability of bony intervertebral fusion. A solid free-form fabrication process based on melt extrusion has been utilized in the manufacturing of these scaffolds. To date there are no studies examining the use of such biodegradable implants in a sheep thoracic spine model. The success of anterior scoliosis surgery in humans depends on achieving a solid bony fusion between adjacent vertebrae after the intervertebral discs have been surgically cleared and the disc spaces filled with graft material. Due to limited availability of autograft, there is much current interest in the development of synthetic scaffolds in combination with growth factors such as recombinant human bone morphogenetic protein (rhBMP-2) to achieve a solid bony fusion following scoliosis surgery.

Post-operative CT assessment of interbody fusion two years after thoracoscopic scoliosis surgery

The relationship between radiologic union and clinical outcome in thoracoscopic scoliosis surgery is not clear, as apparent non-union of a spinal fusion does not always correspond to a poor clinical result. The aim of this study was to evaluate CT fusion rates 2yrs after thoracoscopic surgery, and to explore the relationship between fusion scores and rod diameter, graft type, fusion level, implant failure, and lateral position in the disc space. This study suggests that moderate fusion scores secure successful clinical outcomes in thoracoscopic scoliosis surgery.

Polycaprolactone-based scaffold plus recombinant human bone morphogenic protein rhBMP-2) in a sheep thoracic spine fusion model

Adolescent idiopathic scoliosis is a complex three dimensional deformity affecting 2-3% of the general population. The resulting spinal deformity consists of coronal curvature, hypokyphosis of the thoracic spine and vertebral rotation in the axial plane with posterior elements turned into the curve concavity. The potential for curve progression is heightened during the adolescent growth spurt. Success of scoliosis deformity correction depends on solid bony fusion between adjacent vertebrae after the intervertebral (IV) discs have been surgically cleared and the disc spaces filled with graft material. Recently a bioactive and resorbable scaffold fabricated from medical grade polycaprolactone has been developed for bone regeneration at load bearing sites. Combined with rhBMP-2, this has been shown to be successful in acting as a bone graft substitute in a porcine lumbar interbody fusion model when compared to autologous bone graft alone. The study aimed to establish a large animal thoracic spine interbody fusion model, develop spine biodegradable scaffolds (PCL) in combination with biologics (rhBMP-2) and to establish a platform for research into spine tissue engineering constructs. Preliminary results demonstrate higher grades of radiologically evident bony fusion across all levels when comparing fusion scores between the 3 and 6 month postop groups at the PCL CaP coated scaffold level, which is observed to be a similar grade to autograft, while no fusion is seen at the scaffold only level. Results to date suggest that the combination of rhBMP-2 and scaffold engineering actively promotes bone formation, laying the basis of a viable tissue engineered constructs.

Biomechanical performance of polycaprolactone (PCL)-based scaffold with rhBMP-2 in a sheep thoracic spine fusion model

Adolescent idiopathic scoliosis is a complex three dimensional deformity affecting 2-3% of the general population. Resulting spine deformities include progressive coronal curvature, hypokyphosis, or frank lordosis in the thoracic spine and vertebral rotation in the axial plane with posterior elements turned into the curve concavity. The potential for curve progression is heightened during the adolescent growth spurt. Success of scoliosis deformity correction depends on solid bony fusion between adjacent vertebrae after the intervertebral discs have been surgically cleared and the disc spaces filled with graft material. Problems with bone graft harvest site morbidity as well as limited bone availability have led to the search for bone graft substitutes. Recently, a bioactive and resorbable scaffold fabricated from medical grade polycaprolactone (PCL) has been developed for bone regeneration at load bearing sites. Combined with recombinant human bone morphogenic protein–2 (rhBMP-2), this has been shown to be successful in acting as a bone graft substitute in acting as a bone graft substitute in a porcine lumbar interbody fusion model when compared to autologous bone graft. This in vivo sheep study intends to evaluate the suitability of a custom designed medical grade PCL scaffold in combination with rhBMP-2 as a bone graft substitute in the setting of mini–thoracotomy surgery as a platform for ongoing research to benefit patients with adolescent idiopathic scoliosis.

An FE investigation simulating intra-operative corrective forces applied to correct scoliosis deformity

Background: Adolescent idiopathic scoliosis (AIS) is a deformity of the spine, which may 34 require surgical correction by attaching a rod to the patient’s spine using screws 35 implanted in the vertebral bodies. Surgeons achieve an intra-operative reduction in the 36 deformity by applying compressive forces across the intervertebral disc spaces while 37 they secure the rod to the vertebra. We were interested to understand how the 38 deformity correction is influenced by increasing magnitudes of surgical corrective forces 39 and what tissue level stresses are predicted at the vertebral endplates due to the 40 surgical correction. 41 Methods: Patient-specific finite element models of the osseoligamentous spine and 42 ribcage of eight AIS patients who underwent single rod anterior scoliosis surgery were 43 created using pre-operative computed tomography (CT) scans. The surgically altered 44 spine, including titanium rod and vertebral screws, was simulated. The models were 45 analysed using data for intra-operatively measured compressive forces – three load 46 profiles representing the mean and upper and lower standard deviation of this data 47 were analysed. Data for the clinically observed deformity correction (Cobb angle) were 48 compared with the model-predicted correction and the model results investigated to 49 better understand the influence of increased compressive forces on the biomechanics of 50 the instrumented joints. 51 Results: The predicted corrected Cobb angle for seven of the eight FE models were 52 within the 5° clinical Cobb measurement variability for at least one of the force profiles. 53 The largest portion of overall correction was predicted at or near the apical 54 intervertebral disc for all load profiles. Model predictions for four of the eight patients 55 showed endplate-to-endplate contact was occurring on adjacent endplates of one or 56 more intervertebral disc spaces in the instrumented curve following the surgical loading 57 steps. 58 Conclusion: This study demonstrated there is a direct relationship between intra-59 operative joint compressive forces and the degree of deformity correction achieved. The 60 majority of the deformity correction will occur at or in adjacent spinal levels to the apex 61 of the deformity. This study highlighted the importance of the intervertebral disc space 62 anatomy in governing the coronal plane deformity correction and the limit of this 63 correction will be when bone-to-bone contact of the opposing vertebral endplates 64 occurs.

Establishment and characterisation of an open mini-thoracotomy surgical approach to an ovine thoracic spine fusion model

Background A large animal model is required for assessment of minimally invasive, tissue engineering based approaches to thoracic spine fusion, with relevance to deformity correction surgery for human adolescent idiopathic scoliosis. Here we develop a novel open mini–thoracotomy approach in an ovine model of thoracic interbody fusion which allows assessment of various fusion constructs, with a focus on novel, tissue engineering based interventions. Methods The open mini-thoracotomy surgical approach was developed through a series of mock surgeries, and then applied in a live sheep study. Customized scaffolds were manufactured to conform with intervertebral disc space clearances required of the study. Twelve male Merino sheep aged 4 to 6 years and weighing 35 – 45 kg underwent the abovementioned procedure and were divided into two groups of six sheep at survival timelines of 6 and 12 months. Each sheep underwent a 3-level discectomy (T6/7, T8/9 and T10/11) with randomly allocated implantation of a different graft substitute at each of the three levels; (i) polycaprolactone (PCL) based scaffold plus 0.54μg rhBMP-2, (ii) PCL-based scaffold alone or (iii) autograft. The sheep were closely monitored post- operatively for signs of pain (i.e. gait abnormalities/ teeth gnawing/ social isolation). Fusion assessments were conducted post-sacrifice using Computed Tomography and hard-tissue histology. All scientific work was undertaken in accordance with the study protocol has been approved by the Institute's committee on animal research. Results. All twelve sheep were successfully operated on and reached the allotted survival timelines, thereby demonstrating the feasibility of the surgical procedure and post-operative care. There were no significant complications and during the post-operative period the animals did not exhibit marked signs of distress according to the described assessment criteria. Computed Tomographic scanning demonstrated higher fusion grades in the rhBMP-2 plus PCL-based scaffold group in comparison to either PCL-based scaffold alone or autograft. These results were supported by histological evaluation of the respective groups. Conclusion. This novel open mini-thoracotomy surgical approach to the ovine thoracic spine represents a safe surgical method which can reproducibly form the platform for research into various spine tissue engineered constructs (TEC) and their fusion promoting properties.

Segmental torso masses and coronal plane joint torques 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 determined. This study used CT scans of AIS patients to measure segmental torso masses and explores how joint moments in the coronal plane are affected by changes in the position of the intervertebral joint’s axis of rotation; particularly at the apex of a scoliotic major curve. Methods. Existing low dose CT data from the Paediatric Spine Research Group was used to calculate vertebral level-by-level torso masses and joint torques occurring in the spine for a group of 20 female AIS patients (mean age 15.0 ± 2.7 years, mean Cobb angle 53 ± 7.1°). Image processing software, ImageJ (v1.45 NIH USA) was used to threshold the T1 to L5 CT images and calculate the segmental 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 (IV) joint torques at each vertebral level were found using principles of static equilibrium together with the segmental body mass data. Summing the torque contributions for each level above the required joint, allowed the cumulative joint torque at a particular level to be found. Since there is some uncertainty in the position of the coronal plane Instantaneous Axis of Rotation (IAR) for scoliosis patients, it was assumed the IAR was located in the centre of the IV disc. A sensitivity analysis was performed to see what effect the IAR had on the joint torques by moving it laterally 10mm in both directions. Results. The magnitude of the torso masses from T1-L5 increased inferiorly, with a 150% increase in mean segmental torso mass from 0.6kg at T1 to 1.5kg at L5. The magnitudes of the calculated coronal plane joint torques during relaxed standing were typically 5-7 Nm at the apex of the curve, with the highest apex joint torque of 7Nm being found in patient 13. Shifting the assumed IAR by 10mm towards the convexity of the spine, increased the joint torque at that level by a mean 9.0%, showing that calculated joint torques were moderately sensitive to the assumed IAR location. When the IAR midline position was moved 10mm away from the convexity of the spine, the joint torque reduced by a mean 8.9%. Conclusion. Coronal plane joint torques as high as 7Nm can occur during relaxed standing in scoliosis patients, which may help to explain the mechanics of AIS progression. This study provides new anthropometric reference data on vertebral level-by-level torso mass in AIS patients which will be useful for biomechanical models of scoliosis progression and treatment. However, the CT scans were performed in supine (no gravitational load on spine) and curve magnitudes are known to be smaller than those measured in standing.

Segmental torso masses and gravity-induced coronal plane joint moments in adolescent idiopathic scoliosis

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. This study used supine CT scans of AIS patients to measure segmental torso masses and explored the joint moments in the coronal plane, particularly at the apex of a scoliotic major curve. Methods: Existing low dose CT data from the Paediatric Spine Research Group was used to calculate vertebral level-by-level torso masses and joint moments occurring in the spine for a group of 20 female AIS patients with right sided thoracic curves. The mean age was 15.0 ± 2.7 years and all curves were classified Lenke Type 1 with a mean Cobb angle 52 ± 5.9°. Image processing software, ImageJ (v1.45 NIH USA) was used to create reformatted coronal plane images, reconstruct vertebral level-by-level torso segments and subsequently measure the torso volume corresponding to each vertebral level. Segment mass was then determined by assuming a tissue density of 1.04x103 kg/m3. Body segment masses for the head, neck and arms were taken from published anthropometric data (Winter 2009). Intervertebral joint moments in the coronal plane at each vertebral level were found from the position of the centroid of the segment masses relative to the joint centres with the segmental body mass data. Results and Discussion: The magnitude of the torso masses from T1-L5 increased inferiorly, with a 150% increase in mean segmental torso mass from 0.6kg at T1 to 1.5kg at L5. The magnitudes of the calculated coronal plane joint moments during relaxed standing were typically 5-7 Nm at the apex of the curve, with the highest apex joint torque of 7Nm. The CT scans were performed in the supine position and curve magnitudes are known to be 7-10° smaller than those measured in standing, due to the absence of gravity acting on the spine. Hence, it can be expected that the moments produced by gravity in the standing individual will be greater than those calculated here.

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.1 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.2 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.

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.

6 versus 12 month performance of polycaprolactone-based scaffold plus recombinant human morphogenetic protein-2 (rhBMP-2) in an ovine thoracic interbody fusion model

Introduction Well-designed biodegradable scaffolds in combination with bone growth factors offer a valuable alternative to the current gold standard autograft in spinal fusion surgery Yong et al. (2013). Here we report on 6- vs 12- month data set evaluating the longitudinal performance of a CaP coated polycaprolactone (PCL) scaffold loaded with recombinant human bone morphogenetic protein-2 (rhBMP-2) as a bone graft substitute within a large preclinical animal model. Methods Twelve sheep underwent a 3-level (T6/7, T8/9 and T10/11) discectomy with randomly allocated implantation of a different graft substitute at each of the three levels; (i) calcium phosphate (CaP) coated polycaprolactone based scaffold plus 0.54µg rhBMP-2, (ii) CaP coated PCL- based scaffold alone or (iii) autograft (mulched rib head). Fusion assessments were performed via high resolution clinical computed tomography and histological evaluation were undertaken at six (n=6) and twelve (n=6) months post-surgery using the Sucato grading system (Sucato et al. 2004). Results The computed tomography fusion grades of the 6- and 12- months in the rhBMP-2 plus PCL- based scaffold group were 1.9 and 2.1 respectively, in the autograft group 1.9 and 1.3 respectively, and in the scaffold alone group 0.9 and 1.17 respectively. There were no statistically significant differences in the fusion scores between 6- and 12- month for the rhBMP plus PCL- based scaffold or PCL – based scaffold alone group however there was a significant reduction in scores in the autograft group. These scores were seen to correlate with histological evaluations of the respective groups. Conclusions The results of this study demonstrate the efficacy of scaffold-based delivery of rhBMP-2 in promoting higher fusion grades at 6- and 12- months in comparison to the scaffold alone or autograft group within the same time frame. Fusion grades achieved at six months using PCL+rhBMP-2 are not significantly increased at twelve months post-surgery.

A biomechanical analysis of growing rods used in the management of early onset scoliosis

INTRODUCTION Managing spinal deformities in young children is challenging, particularly early onset scoliosis (EOS). Surgical intervention is often required if EOS has been unresponsive to conservative treatment particularly with rapidly progressive curves. An emerging treatment option for EOS is fusionless scoliosis surgery. Similar to bracing, this surgical option potentially harnesses growth, motion and function of the spine along with correcting spinal deformity. Dual growing rods are one such fusionless treatment, which aims to modulate growth of the vertebrae. The aim of this study was to ascertain the extent to which semiconstrained growing rods (Medtronic, Sofamor, Danek, Memphis, TN) with a telescopic sleeve component, reduce rotational constraint on the spine compared with standard "constrained / rigid" rods and hence potentially provide a more physiological mechanical environment for the growing spine. METHODS Six 40-60kg English Large White porcine spines served as a model for the paediatric human spine. Each spine was dissected into a 7 level thoracolumbar multi-segment unit (MSU), removing all non-ligamentous soft tissues and leaving 3cm of ribs either side. Pure nondestructive axial rotation moments of ±4Nm at a constant rotation rate of 8deg.s-1 were applied to the mounted MSU spines using a biaxial Instron testing machine. Displacement of each vertebral level was captured using a 3D motion tracking system (Optotrak 3020, Northern Digital Inc, Waterloo, ON). Each spine was tested in an un-instrumented state first and then with appropriately sized semi-constrained growing rods and rigid rods in alternating sequence. The rods were secured by multi-axial pedicle screws (Medtronic CD Horizon) at levels 2 and 6 of the construct. The range of motion (ROM), neutral zone (NZ) size and stiffness (Nm.deg-1) were calculated from the Instron load-displacement data and intervertebral ROM was calculated through a MATLAB algorithm from Optotrak data. RESULTS Irrespective of the order of testing, rigid rods significantly reduced the total ROM compared with semi-constrained rods (p<0.05) with in a significantly stiffer spine for both left and right axial rotation (p<0.05). Analysing the intervertebral motion within the instrumented levels 2-6, rigid rods showed reduced ROM compared with semi-constrained growing rods and compared with un-instrumented motion segments. CONCLUSION Semi-constrained growing rods maintain similar stiffness in axial rotation to un-instrumented spines, while dual rigid rods significantly reduce axial rotation. Clinically the effect of semi-constrained growing rods as observed in this study is that they would be expected to allow growth via the telescopic rod components while maintaining the axial flexibility of the spine, which may reduce occurrence of the crankshaft phenomenon.

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.

A biomechanical analysis of growing rods used in the management of early onset scoliosis

INTRODUCTION Managing spinal deformities in young children is challenging, particularly early-onset scoliosis (EOS). Any progressive spinal deformity particularly in early life presents significant health risks for the child and a challenge for the treating surgeon. Surgical intervention is often required if EOS has been unresponsive to conservative treatment particularly with rapidly progressive curves. An emerging treatment option particularly for EOS is fusionless scoliosis surgery. Similar to bracing this surgical option potentially harnesses growth, motion and function of the spine along with correcting spinal deformity. Dual growing rods is one such fusionless treatment, which aims to modulate growth of the vertebrae. The aim of this study was to ascertain the extent to which semi-constrained growing rods (Medtronic, Memphis, TN) with a telescopic sleeve component, reduce rotational constraint on the spine compared with standard rigid rods and hence potentially provide a more physiological mechanical environment for the growing spine. METHODS Six 40-60kg English Large White porcine spines served as a model for the paediatric human spine. Each spine was dissected into 7 level thoracolumbar multi-segment unit (MSU) spines, removing all non-ligamentous soft tissues. Appropriately sized semi-constrained growing rods and rigid rods were secured by multi-axial screws (Medtronic) prior to testing in alternating sequences for each spine. Pure nondestructive moments of +/4Nm at a constant rotation rate of 8deg/s was applied to the mounted MSU spines. Displacement of each level was captured using an Optotrak (Northern Digital Inc, Waterloo, ON). The range of motion (ROM), neutral zone (NZ) size and stiffness (Nm/deg) were calculated from the Instron load-displacement data and intervertebral ROM was calculated through a MATLAB algorithm from Optotrak data. RESULTS Irrespective of sequence order rigid rods significantly reduced the total ROM (deg) than compared to semi-constrained rods (p<0.05) and resulted in a significantly stiffer (Nm/deg) spine for both left and right axial rotation testing (p<0.05). Analysing the intervertebral motion within the instrumented levels, rigid rods showed reduced ROM (Deg) than compared to semi-constrained growing rods and the un-instrumented (UN-IN) test sequences. CONCLUSION The semi-constrained growing rods maintained rotation similar to UN-IN spines while the rigid rods showed significantly reduced axial rotation across all instrumented levels. Clinically the effect of semi-constrained growing rods evaluated in this study is that they will allow growth via the telescopic rod components while maintaining the axial rotation ability of the spine, which may also reduce the occurrence of the crankshaft phenomenon.

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.