A combined numerical and experimental technique for estimation of the forces and moments in the lumbar intervertebral disc

Evaluation of the loads on lumbar intervertebral discs (IVD) is critically important since it is closely related to spine biomechanics, pathology and prosthesis design. Non-invasive estimation of the loads in the discs remains a challenge. In this study, we proposed a new technique to estimate in vivo loads in the IVD using a subject-specific finite element (FE) model of the disc and the kinematics of the disc endplates as input boundary conditions. The technique was validated by comparing the forces and moments in the discs calculated from the FE analyses to the in vitro experiment measurements of three corresponding lumbar discs. The results showed that the forces and moments could be estimated within an average error of 20%. Therefore, this technique can be a promising tool for non-invasive estimation of the loads in the discs and may be extended to be used on living subjects.     

Effect of microgravity on the biomechanical properties of lumbar and caudal intervertebral discs in mice

Prolonged exposure to microgravity has shown to have deleterious effects on the human spine, indicated by low back pain during spaceflight and increased incidence of post-spaceflight herniated nucleus pulposus. We examined the effect of microgravity on biomechanical properties of lumbar and caudal discs from mice having been on 15-day shuttle mission STS-131. Sixteen C57BL/C mice (spaceflight group, n=8; ground-based control group, n=8) were sacrificed immediately after spaceflight. Physiological disc height (PDH) was measured in situ, and compressive creep tests were performed to parameterize biomechanical properties into endplate permeability (k), nuclear swelling pressure strain dependence (D), and annular viscoelasticity (G). For caudal discs, the spaceflight group exhibited 32% lower PDH, 70% lower D and crept more compared to the control mice (p=0.03). For lumbar discs, neither PDH nor D was significantly different between murine groups. Initial modulus, osmotic pressure, k and G for lumbar and caudal discs did not appear influenced by microgravity (p>0.05). Decreases in both PDH and D suggest prolonged microgravity effectively diminished biomechanical properties of caudal discs. By contrast, differences were not noted for lumbar discs. This potentially deleterious interaction between prolonged weightlessness and differential ranges of motion along the spine may underlie the increased cervical versus lumbar disc herniation rates observed among astronauts.     

Nucleus implant parameters significantly change the compressive stiffness of the human lumbar intervertebral disc

Nucleus replacement by a synthetic material is a recent trend for treatment of lower back pain. Hydrogel nucleus implants were prepared with variations in implant modulus, height, and diameter Human lumbar intervertebral discs (IVDs) were tested in compression for intact, denucleated, and implanted condition. Implantation of nucleus implants with different material and geometric parameters into a denucleated IVD significantly altered the IVD compressive stiffness. Variations in the nucleus implant parameters significantly change the compressive stiffness of the human lumbar IVD. Implant geometrical variations were more effective than those of implant modulus variations in the range examined.     

Prostheses designed for vertebral body replacement

Prostheses are proposed to restore the spinal stability of patients suffering from metastatic malignant tumours in their vertebral bodies. They are designed to replace only one vertebral body and two neighbouring intervertebral discs of the spine. Experiments performed on cervical, thoracic and lumbar sections, which were obtained from fresh cadavers, have shown that the reduction in average compressive strengths of these regions due to the placement of prostheses is about 9%. This seems acceptable for those patients in performing their daily activities. The same amount of reduction has also been observed in average compressive strengths of the neighbouring healthy vertebrae due to the placement of prosthesis heads by bone cement. Developed prostheses have a number of advantages over the existing fusion constructs for the cases considered in this work.     

The influence of exogenous cross-linking and compressive creep loading on intradiscal pressure

This study involves a biomechanical evaluation of a prospective injectable treatment for degenerative discs. The high osmolarity of the non-degenerated nucleus pulposus attracts water contributing to the hydrostatic behavior of the tissue. This intradiscal pressure is known to drop as fluid is exuded from the matrix due to compressive loading. The objective of this study was to compare the changes in intradiscal pressure in control and genipin cross-linked intervertebral discs. Thirty bovine lumbar motion segments were randomly divided into a phosphate-buffered saline control group and a 0.33% genipin group and soaked at room temperature for 2 days. A needle pressure sensor was held in the center of the disc while short-term and static creep compressive loads were applied. The control group demonstrated a 25% higher average intradiscal pressure compared to genipin-treated discs under 750 N compressive load (p=0.029). Depressurization during static compressive creep was 56% higher in the control than in the genipin group (p=0.014). These results suggest cross-linking induced changes in the poroelastic properties of the involved tissues affected the mechanics of compressive load support in the disc with lower levels of nucleus pressure, a corresponding decrease in the elastic expansion of the annulus, and an increased axial compressive loading of the inner and outer annulus tissues. It is possible that concurrent changes in hydraulic permeability and proteoglycan retention known to be associated with genipin cross-linking were also contributors to poroelastic changes. Reduction of peak pressures and moderation of pressure fluctuations could be beneficial relative to discogenic pain.     

Annulus Fibrosus Cell Characteristics Are a Potential Source of Intervertebral Disc Pathogenesis

In the end stage of intervertebral disc degeneration, cartilage, bone, endothelial cells, and neurons appear in association with the worsening condition. The origin of the abnormal cells is not clear. This study investigated the properties of progenitor cells in the annulus fibrosus (AF) using one in vitro and two in vivo models. Cultivation of rabbit AF cells with chondrogenic media significantly increased expressions of collagen and aggrecan. Upon exposure to osteogenic conditions, the cultures showed increased mineralization and expression of osteopontin, runx2, and bmp2 genes. Two models were used in the in vivo subcutaneous implantation experiments: 1) rabbit AF tissue in a demineralized bone matrix (DBM) cylinder (DBM/AF), and, 2) rat intact and needle punctured lumbar discs. Bone formation in the AF tissue was detected and hypertrophic chondrocytes and osteoblasts were present 1 month after implantation of the DBM/AF to nude mice. In addition to collagen I and II, immunostaining shows collagen X and osteocalcin expression in DBM/AF specimens 4 months after implantation. Similar changes were detected in the injured discs. Almost the entire needle punctured disc had ossified at 6 months. The results suggest that AF cells have characteristics of progenitor cells and, under appropriate stimuli, are capable of differentiating into chondrocytes and osteoblasts in vitro as well as in vivo. Importantly, these cells may be a target for biological treatment of disc degeneration.     

Role of endplates in contributing to compression behaviors of motion segments and intervertebral discs

The purpose of this study was to gain an improved understanding of the mechanical behavior of the intervertebral disc in the presence and absence of the vertebral endplates. Mechanical behaviors of rat caudal motion segments, vertebrae and isolated disc explants under two different permeability conditions were investigated and viscoelastic behaviors were evaluated using a stretched-exponential function to describe creep and recovery behaviors. The results demonstrated that both vertebrae and discs underwent significant deformations in the motion segment even under relatively low-loading conditions. Secondly, disruption of the collagenous network had minimal impact on equilibrium deformations of disc explants as compared to disc deformations occurring in the motion segments provided that vertebral deformations were accounted for; however, differences in endplate permeability conditions had a significant effect on viscoelastic behaviors. Creep occurred more quickly than recovery for motion segment and explant specimens. In addition, disc explants and motion segments both exhibited non-recoverable deformations under axial compression under low- and high-loading conditions. Results have important implications for interpreting the role of vertebral endplates in contributing to disc mechanical behaviors and direct application to mechanobiology studies involving external loading to rodent tail intervertebral discs.     

Contribution of vertebral [corrected] bodies, endplates, and intervertebral discs to the compression creep of spinal motion segments

Spinal segments show non-linear behavior under axial compression. It is unclear to what extent this behavior is attributable to the different components of the segment. In this study, we quantified the separate contributions of vertebral bodies and intervertebral discs to creep of a segment. Secondly, we investigated the contribution of bone and osteochondral endplate (endplates including cartilage) to the deformation of the vertebral body. From eight porcine spines a motion segment, a disc and a vertebral body were dissected and subjected to mechanical testing. In an additional test, cylindrical samples, machined from the lowest thoracic vertebrae of 11 porcine spines, were used to compare the deformation of vertebral bone and endplate. All specimens were subjected to three loading cycles, each comprising a loading phase (2.0 MPa, 15 min) and a recovery phase (0.001 MPa, 30 min). All specimens displayed substantial time-dependent height changes. Average creep was the largest in motion segments and smallest in vertebral bodies. Bone samples with endplates displayed substantially more creep than samples without. In the early phase, behavior of the vertebra was similar to that of the disc. Visco-elastic deformation of the endplate therefore appeared dominant. In the late creep phase, behavior of the segment was similar to that of isolated discs, suggesting that in this phase the disc dominated creep behavior, possibly by fluid flow from the nucleus. We conclude that creep deformation of vertebral bodies contributes substantially to creep of motion segments and that within a vertebral body endplates play a major role.     

Computational study of the role of fluid content and flow on the lumbar disc response in cyclic compression: Replication of in vitro and in vivo conditions

The intervertebral disc viscoelastic response is governed primarily by its fluid content and flow. In vivo measurements demonstrate that the disc volume, fluid content, height and nucleus pressure completely recover during resting even after diurnal loading with twice longer duration (16 vs. 8 h). In view of much longer periods required for the recovery of disc height and pressure in vitro, concerns have been raised on the fluid inflow through the endplates that might be hampered by clogged blood vessels post mortem. This in silico study aimed to identify fluid-flow dependent response of discs and conditions essential to replicate in vitro and in vivo observations.An osmo-poroelastic finite element model of the human lumbar L4-L5 disc-bone unit was used. Simulating earlier in vitro experiments on bovine discs, the loading protocol started with 8 h preload at 0.06 MPa followed by 30 high/low compression loading cycles each lasting 7.5 min at 0.5/0.06 MPa, respectively. Three different endplate configurations were investigated: free in- and outflow, no inflow and closed endplates with no flow. Additionally, the preload magnitude was increased from 0.06 MPa to 0.28 MPa and 0.50 MPa, or the initial nucleus hydration was reduced from 83% to 50%.For 0.06 MPa preload, the model with no inflow best matched in vitro trends. The model with free inflow increased segment height and nucleus pressure while the model with no fluid inflow resulted in a relatively small recovery in segment height and a rather constant nucleus pressure during unloading periods.Results highlight an excessive mobile fluid content as well as a restricted fluid inflow through endplates as likely causes of the discrepancies between in vivo and in vitro studies. To replicate in vivo conditions in vitro and in silico, disc hydration level should be controlled by adequate selection of preload magnitude/period and/or mobile fluid porosity.     

Optimization of compressive loading parameters to mimic in vivo cervical spine kinematics in vitro

The human cervical spine supports substantial compressive load in vivo. However, the traditional in vitro testing methods rarely include compressive loads, especially in investigations of multi-segment cervical spine constructs. Previously, a systematic comparison was performed between the standard pure moment with no compressive loading and published compressive loading techniques (follower load – FL, axial load – AL, and combined load – CL). The systematic comparison was structured a priori using a statistical design of experiments and the desirability function approach, which was chosen based on the goal of determining the optimal compressive loading parameters necessary to mimic the segmental contribution patterns exhibited in vivo. The optimized set of compressive loading parameters resulted in in vitro segmental rotations that were within one standard deviation and 10% of average percent error of the in vivo mean throughout the entire motion path. As hypothesized, the values for the optimized independent variables of FL and AL varied dynamically throughout the motion path. FL was not necessary at the extremes of the flexion–extension (FE) motion path but peaked through the neutral position, whereas, a large negative value of AL was necessary in extension and increased linearly to a large positive value in flexion. Although further validation is required, the long-term goal is to develop a “physiologic” in vitro testing method, which will be valuable for evaluating adjacent segment effect following spinal fusion surgery, disc arthroplasty instrumentation testing and design, as well as mechanobiology experiments where correct kinematics and arthrokinematics are critical.     

Creation of an in vitro biomechanical model of the trachea using rapid prototyping

Previous in vitro models of the airways are either rigid or, if flexible, have not matched in vivo compliance characteristics. Rapid prototyping provides a quickly evolving approach that can be used to directly produce in vitro airway models using either rigid or flexible polymers. The objective of this study was to use rapid prototyping to directly produce a flexible hollow model that matches the biomechanical compliance of the trachea. The airway model consisted of a previously developed characteristic mouth–throat region, the trachea, and a portion of the main bronchi. Compliance of the tracheal region was known from a previous in vivo imaging study that reported cross-sectional areas over a range of internal pressures. The compliance of the tracheal region was matched to the in vivo data for a specific flexible resin by iteratively selecting the thicknesses and other dimensions of tracheal wall components. Seven iterative models were produced and illustrated highly non-linear expansion consisting of initial rapid size increase, a transition region, and continued slower size increase as pressure was increased. Thickness of the esophageal interface membrane and initial trachea indention were identified as key parameters with the final model correctly predicting all phases of expansion within a value of 5% of the in vivo data. Applications of the current biomechanical model are related to endotracheal intubation and include determination of effective mucus suctioning and evaluation of cuff sealing with respect to gases and secretions.     

Transdetermination of Drosophila imaginal discs cultured in vitro

Imaginal discs of Drosophila melanogaster undergo transdetermination when cultured in vivo in the abdominal cavity of adult female hosts. We report here that leg discs cultured in vitro, in a recently developed system, also undergo transdetermination. Whether cultured in vivo or in vitro, leg discs produce a similar range of specific transdetermined structures. Moreover, in comparison to discs cultured in vivo, the discs cultured in vitro exhibit a similar correlation between the amount of growth and the total frequency of transdetermination.     

Neural arch load-bearing in old and degenerated spines

We validate a technique for measuring neural arch load-bearing in cadaveric spines, and use it to test the hypothesis that such load-bearing rises to high levels in old and degenerated spines. Fifty-nine cadaveric lumbar motion segments, aged 19-92 yr, were subjected to compressive creep loading to reduce intervertebral disc water content and height to in vivo levels. The distribution of compressive "stress" within the disc was then measured by pulling a miniature pressure transducer, side-mounted in a 1.3mm-diameter needle, along its mid-sagittal diameter. During these measurements, the motion segment was subjected to a compressive load of 2 kN, and positioned in 2 degrees of extension to simulate erect standing. Measurements of compressive "stress" were integrated over disc area, and this force subtracted from the applied 2 kN to give the force resisted by the neural arch. An empirical calibration factor was applied to normalise results from each disc to values obtained under conditions when all of the compressive force could be assumed to pass through the disc. Disc degeneration was graded macroscopically on a scale of 1-4. Validation tests showed that calculated values of disc loading were proportional to actual applied load (r(2)>0.96) and predicted it with errors of 2-8%. Neural arch load-bearing in non-degenerated specimens was generally less than 20%, but averaged 49% for specimens aged over 70 yr. Multiple regression showed that neural arch load bearing (%)=14.4 x disc degeneration score+0.46 x age-35. These results indicate a substantial shift in vertebral load-bearing with increasing age and degeneration.     

Preliminary Investigations on Intradiscal Pressures during Daily Activities: An In Vivo Study Using the Merino Sheep

Purpose

Currently, no studies exist, which attest the suitability of the ovine intervertebral disc as a biomechanical in vivo model for preclinical tests of new therapeutic strategies of the human disc. By measuring the intradiscal pressure in vivo, the current study attempts to characterize an essential biomechanical parameter to provide a more comprehensive physiological understanding of the ovine intervertebral disc.

Methods

Intradiscal pressure (IDP) was measured for 24 hours within the discs L2-L3 and L4-L5 via a piezo-resistive pressure sensor in one merino sheep. The data were divided into an activity and a recovery phase and the corresponding average pressures for both phases were determined. Additionally, IDPs for different static and dynamic activities were analyzed and juxtaposed to human data published previously. After sacrificing the sheep, the forces corresponding to the measured IDPs were examined ex vivo in an axial compression test.

Results

The temporal patterns of IDP where pressure decreased during activity and increased during rest were comparable between humans and sheep. However, large differences were observed for different dynamic activities such as standing up or walking. Here, IDPs averaged 3.73 MPa and 1.60 MPa respectively, approximately two to four times higher in the ovine disc compared to human. These IDPs correspond to lower ex vivo derived axial compressive forces for the ovine disc in comparison to the human disc. For activity and rest, average ovine forces were 130 N and 58 N, compared to human forces of 400-600 N and 100 N, respectively.

Conclusions

In vivo IDPs were found to be higher in the ovine than in the human disc. In contrast, axial forces derived ex vivo were markedly lower in comparison to humans. Both should be considered in future preclinical tests of intradiscal therapies using the sheep. The techniques used in the current study may serve as a protocol for measuring IDP in a variety of large animal models.     

A non-optimized follower load path may cause considerable intervertebral rotations

Osseoligamentous spinal specimens buckle under even a small vertical compressive force. To allow higher axial forces, a compressive follower load (FL) was suggested previously that approximates the curvature of the spine without inducing intervertebral rotation in both the frontal and the sagittal planes. In in vitro experiments and finite element analyses, the location of the FL path is subjected to estimation by the investigator. Such non-optimized FLs may induce bending and so far it is still unknown how this affects the results of the study and their comparability.A symmetrical finite element model of the lumbar spine was employed to simulate upright standing while applying a follower load. In analogy to in vitro experiments, the path of this FL was estimated seven times by different members of our institute’s spine group. Additionally, an optimized FL path was determined and additional moments of ±7.5 Nm were applied to simulate flexion and extension.Application of the optimized 500 N compressive FL causes only a marginal alteration of the curvature (cardan angle L1–S1 in sagittal plane <0.25°). An individual estimation of the FL path, however, results in flexions of up to 10.0° or extensions of up to 12.3°. The resulting angles for the different non-optimized FL paths depend on the magnitude of the bending moment applied and whether a differential or an absolute measurement is taken.A preceding optimization of the location of the FL path would increase the comparability of different studies.     

Validation of single-segment and three-segment spinal models used to represent lumbar flexion

Changes in spinal posture between the erect and flexed positions were calculated using angular measurements from lateral photographs and radiographs of ten adult male subjects. For photographic measurements, the thoracolumbar vertebral column was modelled as either a single segment or as three segments. In the three-segment model, there was a non-significant correlation between the decrease in lumbar concavity and intervertebral motion. In addition, there was a non-significant negative correlation between the increase in thoracic convexity and lumbar motion determined radiographically. In the single-segment model, the decrease in angulation between the thoracolumbar spine and pelvis was a good representation of lumbar spine flexion as determined by the mean lumbar intervertebral angular change. Therefore, modelling the thoracolumbar vertebral column as a single segment allowed better estimation of lumbar intervertebral angular change during flexion than a three-segment model. The results indicate that large range dynamic motion of the lumbar vertebral column can be represented using photographic analysis of the positions of three easily identified anatomical landmarks: the anterior superior iliac spine, posterior superior iliac spine and the spinous process of the first thoracic vertebra.     

Potential use of human adipose mesenchymal stromal cells for intervertebral disc regeneration: a preliminary study on biglycan-deficient murine model of chronic disc degeneration

Introduction

Biglycan is an important proteoglycan of the extracellular matrix of intervertebral disc (IVD), and its decrease with aging has been correlated with IVD degeneration. Biglycan deficient (Bgn^−/0) mice lack this protein and undergo spontaneous IVD degeneration with aging, thus representing a valuable in vivo model for preliminary studies on therapies for human progressive IVD degeneration. The purpose of the present study was to assess the possible beneficial effects of adipose-derived stromal cells (ADSCs) implants in the Bgn^−/0 mouse model.

Methods

To evaluate ADSC implant efficacy, Bgn^−/0 mice were intradiscally (L1-L2) injected with 8x10⁴ ADSCs at 16 months old, when mice exhibit severe and complete IVD degeneration, evident on both 7Tesla Magnetic Resonance Imaging (7TMRI) and histology. Placebo and ADSCs treated Bgn^−/0 mice were assessed by 7TMRI analysis up to 12 weeks post-transplantation. Mice were then sacrificed and implanted discs were analyzed by histology and immunohistochemistry for the presence of human cells and for the expression of biglycan and aggrecan in the IVD area.

Results

After in vivo treatment, 7TMRI revealed evident increase in signal intensity within the discs of mice that received ADSCs, while placebo treatment did not show any variation. Ultrastructural analyses demonstrated that human ADSC survival occurred in the injected discs up to 12 weeks after implant. These cells acquired a positive expression for biglycan, and this proteoglycan was specifically localized in human cells. Moreover, ADSC treatment resulted in a significant increase of aggrecan tissue levels.

Conclusion

Overall, this work demonstrates that ADSC implant into degenerated disc of Bgn^−/0 mice ameliorates disc damage, promotes new expression of biglycan and increased levels of aggrecan. This suggests a potential benefit of ADSC implant in the treatment of chronic degenerative disc disease and prompts further studies in this field.     

The role of the nucleus pulposus in neutral zone human lumbar intervertebral disc mechanics

To study the effect of denucleation on the mechanical behavior of the human lumbar intervertebral disc through a 2mm incision, two groups of six human cadaver lumbar spinal units were tested in axial compression, axial rotation, lateral bending and flexion/extension after incremental steps of "partial" denucleation. Neutral zone, range of motion, stiffness, intradiscal pressure and energy dissipation were measured; the results showed that the contribution of the nucleus pulposus to the mechanical behavior of the intervertebral disc was more dominant through the neutral zone than at the farther limits of applied loads and moments.