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.     

Ex vivo measurement of lumbar intervertebral disc pressure using fibre-Bragg gratings

Methods were developed to measure intervertebral disc pressure using optical fibre-Bragg gratings (FBGs). The FBG sensor was calibrated for hydrostatic pressure in a purpose-built apparatus and the average sensitivity was determined to be -5.7 +/- 0.085 pm/MPa (mean +/- SD). The average coefficient of determination (r(2)) for the calibration data was 0.99, and the average hysteresis of the sensor was 2.13% of full scale. The FBG was used to measure intradiscal pressure response to compressive load in five lumbar functional spine units. The pressure measured by the FBG sensor varied linearly with applied compressive load with coefficients of determination ranging from 0.84 to 0.97. The FBG sensor's sensitivity to compressive load ranged from 0.702 +/- 0.043 kPa/N (mean +/- SD) in a L1-L2 specimen, to 1.07 +/- 0.069 kPa/N in a L4-L5 specimen. These measurements agree with those of previous studies in lumbar spines. Two strain gauge pressure sensors were also used to measure intradiscal pressure response to compressive load. The measured pressure sensitivity to load ranged from 0.251 kPa/N (L4-L5) to 0.850 kPa/N (L2-L3). The average difference in pressure sensitivity to load between Sensors 1 and 2 was 12.9% of the value for Sensor 1, with a range from 1.1% to 20.4%, which suggests that disc pressure was not purely hydrostatic. This may have contributed to the difference between the responses of the FBG and strain gauge sensors.     

High pressures and asymmetrical stresses in the scoliotic disc in the absence of muscle loading

Background

Loads acting on scoliotic spines are thought to be asymmetric and involved in progression of the scoliotic deformity; abnormal loading patterns lead to changes in bone and disc cell activity and hence to vertebral body and disc wedging. At present however there are no direct measurements of intradiscal stresses or pressures in scoliotic spines. The aim of this study was to obtain quantitative measurements of the intradiscal stress environment in scoliotic intervertebral discs and to determine if loads acting across the scoliotic spine are asymmetric. We performed in vivo measurements of stresses across the intervertebral disc in patients with scoliosis, both parallel (termed horizontal) and perpendicular (termed vertical) to the end plate, using a side mounted pressure transducer (stress profilometry)

Methods

Stress profilometry was used to measure horizontal and vertical stresses at 5 mm intervals across 25 intervertebral discs of 7 scoliotic patients during anterior reconstructive surgery. A state of hydrostatic pressure was defined by identical horizontal and vertical stresses for at least two consecutive readings. Results were compared with similar stress profiles measured during surgery across 10 discs of 4 spines with no lateral curvature and with data from the literature.

Results

Profiles across scoliotic discs were very different from those of normal, young, healthy discs of equivalent age previously presented in the literature. Hydrostatic pressure regions were only seen in 14/25 discs, extended only over a short distance. Non-scoliotic discs of equivalent age would be expected to show large centrally placed hydrostatic nuclear regions in all discs. Mean pressures were significantly greater (0.25 MPa) than those measured in other anaesthetised patients (<0.07 MPa). A stress peak was seen in the concave annulus in 13/25 discs. Stresses in the concave annulus were greater than in the convex annulus indicating asymmetric loading in these anaesthetised, recumbent patients.

Conclusion

Intradiscal pressures and stresses in scoliotic discs are abnormal, asymmetrical and high in magnitude even in the absence of significant applied muscle loading. The origin of these abnormal stresses is unclear.     

A simple model for the function of proteoglycans and collagen in the response to compression of the intervertebral disc.

The nucleus pulposus of the intervertebral disc exerts a pressure which enables it to support axial compression when contained by the annulus fibrosus. The disc was modelled as a thick-walled cylindrical pressure vessel in which the nucleus was contained radially by the annulus. As a result, the stress in the annulus had radial (compressive) as well as tangential (tensile) components. The radial stress at a given point in the annulus was considered to be balanced by the internal pressure which is expected to arise from the attraction of water by proteoglycans. There was a reasonable agreement between the calculated radial stress distribution and published results on the distribution of water within the annulus. As the internal pressure is expected to be isotropic, the annulus was expected to contribute to the axial resistance to compression of the disc; this contribution would be equal, in magnitude, to the radial stress. Predictions of the pressure distribution within the annulus were similar to published experimental measurements made in the radial and axial directions. The tangential stress within the annulus was considered to arise from the restoring stress in its strained collagen fibrils.     

Statistical factorial analysis on the material property sensitivity of the mechanical responses of the C4-C6 under compression, anterior and posterior shear

A systematic approach using factorial analysis was conducted on the C4-C6 finite element model to analyse the influence of six spinal components (cortical shell, vertebral body, posterior elements, endplate, disc annulus and disc nucleus) on the internal stresses and external biomechanical responses under compression, anterior and posterior shear. Results indicated that the material properties variation of the disc annulus has a significant influence on both the external biomechanical responses and internal stress of the disc annulus and its neighboring hard bones. The study reveals for the first time, the significant influence of the cancellous bone under compression, while variation in the cortical shell modulus has a high influence under anterior and posterior shear. The study also reveals that the effects of interaction between two main components are insignificant.     

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.     

Variability of a three-dimensional finite element model constructed using magnetic resonance images of a knee for joint contact stress analysis

Magnetic resonance (MR) imaging has been widely used to evaluate the thickness and volume of articular cartilage both in vivo and in vitro. While morphological information on the cartilage can be obtained using MR images, image processing for extracting geometric boundaries of the cartilage may introduce variations in the thickness of the cartilage. To evaluate the variability of using MR images to construct finite element (FE) knee cartilage models, five investigators independently digitized the same set of MR images of a human knee. The topology of cartilage thickness was determined using a minimal distance algorithm. Less than 8 percent variation in cartilage thickness was observed from the digitized data. The effect of changes in cartilage thickness on contact stress analysis was then investigated using five FE models of the knee. One FE model (average FE model) was constructed using the mean values of the digitized contours of the cartilage, and the other four were constructed by varying the thickness of the average FE model by +/- 5 percent and +/- 10 percent, respectively. The results demonstrated that under axial tibial compressive loading (up to 1,400 N), variations of cartilage thickness caused by digitization of MR images may result in a difference of approximately 10 percent in peak contact stresses (surface pressure, von Mises stress, and hydrostatic pressure) in the cartilage. A reduction of cartilage thickness caused increases of contact stresses, while an increase of cartilage thickness reduced contact stresses. Furthermore, the effect of variation of material properties of the cartilage on contact stress analysis was investigated. The peak contact stress increased almost linearly with the Young's modulus of the cartilage. The peak von Mises stress was dramatically reduced when the Poisson,s ratio was increased from 0.05 to 0.49 under an axial compressive load of 1,400 N, while peak hydrostatic pressure was dramatically increased. Peak surface pressure was also increased with the Poisson's ratio, but with a lower magnitude compared to von Mises stress and hydrostatic pressure. In conclusion, the imaging process may cause 10 percent variations in peak contact stress, and the predicted stress distribution is sensitive to the accuracy of the material properties of the cartilage model, especially to the variation of Poisson's ratio.     

A biomechanical model of the lumbosacral joint during lifting activities

A biomechanical model of the lumbosacral region was constructed for the purpose of systematically studying the combined stresses and strains on the local ligaments, muscles and disc tissue during sagittal plane two-handed lifting. The model was validated in two ways. The first validation was a comparison of experimental study results with model predictions. In general predictions compared very reasonably with observed values of several authors with the exception of strain predictions on the articular ligaments. Second, a sensitivity analysis was performed over a wide range of lifting tasks. The predicted stress/strain values followed anticipated patterns and were of reasonable magnitudes. On the basis of the results of the sensitivity analysis it was concluded that typical lifting tasks can lead to excessive disc compressive forces, muscle moment generation requirements, and possibly lumbodorsal fascia strains. Conversely, annulus rupture of a healthy disc due to overstrain appears very unlikely.     

Osmosis and viscoelasticity both contribute to time-dependent behaviour of the intervertebral disc under compressive load: A caprine in vitro study

The mechanical behaviour of the intervertebral disc highly depends on the content and transport of interstitial fluid. It is unknown, however, to what extent the time-dependent behaviour can be attributed to osmosis. Here we investigate the effect of both mechanical and osmotic loading on water content, nucleus pressure and disc height. Eight goat intervertebral discs, immersed in physiological saline, were subjected to a compressive force with a pressure needle inserted in the nucleus. The loading protocol was: 10 N (6 h); 150 N (42 h); 10 N (24 h). Half-way the 150 N-phase (24 h), we eliminated the osmotic gradient by adding 26% poly-ethylene glycol to the surrounding fluid. For 62 additional discs, we determined the water content of both nucleus and annulus after 6, 24, 48, or 72 h. The compressive load was initially counterbalanced by the hydrostatic pressure in the nucleus. The load forced 4.3% of the water out of the nucleus, which reduced nucleus pressure by 44(±6)%. Reduction of the osmotic gradient disturbed the equilibrium disc height, and a significant loss of annulus water content was found. Remarkably, pressure and water content of the nucleus pulposus remained unchanged. This shows that annulus water content is important in the response to axial loading. After unloading, in the absence of an osmotic gradient, there was substantial viscoelastic recovery of 53(±11)% of the disc height, without a change in water content. However, for restoration of the nucleus pressure and for full restoration of disc height, restoration of the osmotic gradient was needed.     

Dynamic unconfined compression of articular cartilage under a cyclic compressive load

To study the effect of dynamic mechanical force on cartilage metabolism, many investigators have applied a cyclic compressive load to cartilage disc explants in vitro. The most frequently used in vitro testing protocol has been the cyclic unconfined compression of articular cartilage in a bath of culture medium. Cyclic compression has been achieved by applying either a prescribed cyclic displacement or a prescribed cyclic force on a loading platen placed on the top surface of a cylindrical cartilage disc. It was found that the separation of the loading platen from the tissue surface was likely when a prescribed cyclic displacement was applied at a high frequency.The purpose of the present study was to simulate mathematically the dynamic behavior of a cylindrical cartilage disc subjected to cyclic unconfined compression under a dynamic force boundary condition protocol, and to provide a parametric analysis of mechanical deformations within the extracellular matrix. The frequency-dependent dynamic characteristics of dilatation, hydrostatic pressure and interstitial fluid velocity were analyzed over a wide range of loading frequencies without the separation of the loading platen. The result predicted that a cyclic compressive force created an oscillating positive-negative hydrostatic pressure together with a forced circulation of interstitial fluid within the tissue matrix. It was also found that the load partitioning mechanism between the solid and fluid phases was a function of loading frequency. At a relatively high loading frequency, a localized dynamic zone was developed near the peripheral free surface of the cartilage disc, where a large dynamic pressure gradient exists, causing vigorous interstitial fluid flow.     

Failure strength of the bovine caudal disc under internal hydrostatic pressure

The structure of the disc is both complex and inhomogeneous, and it functions as a successful load-bearing organ by virtue of the integration of its various structural regions. These same features also render it impossible to assess the failure strength of the disc from isolated tissue samples, which at best can only yield material properties. This study investigated the intrinsic failure strength of the intact bovine caudal disc under a simple mode of internal hydrostatic pressure. Using a hydraulic actuator, coloured hydrogel was injected under monitored pressure into the nucleus through a hollow screw insert which passed longitudinally through one of the attached vertebrae. Failure did not involve vertebra/endplate structures. Rather, failure of the disc annulus was indicated by the simultaneous manifestation of a sudden loss of gel pressure, a flood of gel colouration appearing in the outer annulus and audible fibrous tearing. A mean hydrostatic failure pressure of 18+/-3 MPa was observed which was approximated as a thick-wall hoop stress of 45+/-7 MPa. The experiment provides a measurement of the intrinsic strength of the disc using a method of internal hydrostatic loading which avoids any disruption of the complex architecture of the annular wall. Although the disc in vivo is subjected to a much more complex pattern of loading than is achieved using simple hydrostatic pressurization, this latter mode provides a useful tool for investigating alterations in intrinsic disc strength associated with prior loading history or degeneration.     

In-vivo mitral annuloplasty ring transducer: Implications for implantation and annular downsizing

Mitral annuloplasty has been a keystone to the success of mitral valve repair in functional mitral regurgitation. Understanding the complex interplay between annular-ring stresses and left ventricular function has significant implications for patient-ring selection, repair failure, and patient safety. A step towards assessing these challenges is developing a transducer that can be implanted in the exact method as commercially available rings and can quantify multidirectional ring loading. An annuloplasty ring transducer was developed to measure stresses at eight locations on both the in-plane and out-of-plane surfaces of an annuloplasty ring's titanium core. The transducer was implanted in an ovine subject using 10 sutures at near symmetric locations. At implantation, the ring was observed to undersize the mitral annulus. The flaccid annulus exerted both compressive (−) and tensile stresses (+) on the ring ranging from −3.17 to 5.34 MPa. At baseline hemodynamics, stresses cyclically changed and peaked near mid-systole. Mean changes in cyclic stress from ventricular diastole to mid-systole ranged from −0.61 to 0.46 MPa (in-plane direction) and from −0.49 to 1.13 MPa (out-of-plane direction). Results demonstrate the variability in ring stresses that can be introduced during implantation and the cyclic contraction of the mitral annulus. Ring stresses at implantation were approximately 4 magnitudes larger than the cyclic changes in stress throughout the cardiac cycle. These methods will be extended to ring transducers of differing size and geometry. Upon additional investigation, these data will contribute to improved knowledge of annulus-ring stresses, LV function, and the safer development of mitral repair techniques.     

Numerical analysis of the influence of nucleus pulposus removal on the biomechanical behavior of a lumbar motion segment

Nucleus replacement was deemed to have therapeutic potential for patients with intervertebral disc herniation. However, whether a patient would benefit from nucleus replacement is technically unclear. This study aimed to investigate the influence of nucleus pulposus (NP) removal on the biomechanical behavior of a lumbar motion segment and to further explore a computational method of biomechanical characteristics of NP removal, which can evaluate the mechanical stability of pulposus replacement. We, respectively, reconstructed three types of models for a mildly herniated disc and three types of models for a severely herniated disc based on a L4–L5 segment finite element model with computed tomography image data from a healthy adult. First, the NP was removed from the herniated disc models, and the biomechanical behavior of NP removal was simulated. Second, the NP cavities were filled with an experimental material (Poisson's ratio = 0.3; elastic modulus = 3 MPa), and the biomechanical behavior of pulposus replacement was simulated. The simulations were carried out under the five loadings of axial compression, flexion, lateral bending, extension, and axial rotation. The changes of the four biomechanical characteristics, i.e. the rotation degree, the maximum stress in the annulus fibrosus (AF), joint facet contact forces, and the maximum disc deformation, were computed for all models. Experimental results showed that the rotation range, the maximum AF stress, and joint facet contact forces increased, and the maximum disc deformation decreased after NP removal, while they changed in the opposite way after the nucleus cavities were filled with the experimental material.     

Influence of material properties on the mechanical behaviour of the L5-S1 intervertebral disc in compression: a nonlinear finite element study

A simple finite element model of the L₅-S₁ intervertebral disc body has been constructed; it is circular and symmetrical about the sagittal plane. The annulus fibrosus of the model was idealized as an inhomogeneous composite of an isotropic ground substance, reinforced by helically oriented collagen fibres so that the model has six different structural components namely: cortical bone, cancellous bone, cartilaginous endplates, nucleus pulposus, ground substance and collagen fibres. A sensitivity analysis of the material properties of each structural component was carried out by varying those properties for one structural component at a time and evaluating the changes in the biomechanical response to compressive displacements. Experimentally available relations between the applied compressive force and the vertical displacements, the nucleus pulposus pressure increase and the disc lateral bulge were used to evaluate the biomechanical responses for each set of material properties. Results showed that both the Poisson's ratio and the Young's modulus of the ground substance play an important role in the prediction of the biomechanical response.     

Determination of expander apparatus displacements and contact pressures on the mucosa using FEM modelling considering mandibular asymmetries

This paper presents a method for prediction of forces and displacements in the expansion screw of a modified mandibular Schwarz appliance and the contact pressure distributions on the mucosa during malocclusions treatment. A 3D finite element biomechanical model of the complete mandible–mucosa–apparatus set was built using computerised tomographic images of a patient's mandible and constructive solid geometry by computer software. An iterative procedure was developed to handle a boundary condition that takes into account the mandibular asymmetries. The results showed asymmetries in the contact pressure distributions that indicated with precision the patient's malocclusion diagnosis. In vivo measurements of contact pressure using piezoelectric sensors agreed with the computational results. It was shown that the left and right ends of the expansion screw move differently with respect to the patient mandible, even though the expansion screw has an opening mechanism to ensure equal stretching at both ends. The contact pressures between the apparatus and the mucosa vary linearly with applied forces, which can simplify the analysis of the biomechanical behaviour of the expander mandible apparatus. The biomechanical modelling proposed in this paper can be a useful tool to improve malocclusions treatment, safely avoiding the use of forces acting on live structures beyond the biological tolerance, which could result in traumatic effects.     

Biomechanical effect after Coflex and Coflex rivet implantation for segmental instability at surgical and adjacent segments: a finite element analysis

The Coflex device may provide stability to the surgical segment in extension but does not restore stability in other motion. Recently, a modified version called the Coflex rivet has been developed. The effects of Coflex and Coflex rivet implantation on the adjacent segments are still not clear; therefore, the purpose of this study was to investigate the biomechanical differences between Coflex and Coflex rivet implantation by using finite element analyses. The results show that the Coflex implantation can provide stability in extension, lateral bending, and axial rotation at the surgical segment, and it had no influence at adjacent segments except for extension. The Coflex rivet implantation can provide stability in all motions and reduce disc annulus stress at the surgical segment. Therefore, the higher range of motion and stress induced by the Coflex rivet at both adjacent discs may result in adjacent segment degeneration in flexion and extension.     

Role of load history in intervertebral disc mechanics and intradiscal pressure generation

Solid–fluid interactions play an important role in mediating viscoelastic behaviour of biological tissues. In the intervertebral disc, water content is governed by a number of factors, including age, disease and mechanical loads, leading to changes in stiffness characteristics. We hypothesized that zonal stress distributions depend on load history, or the prior stresses experienced by the disc. To investigate these effects, rat caudal motion segments were subjected to compressive creep biomechanical testing in vitro using a protocol that consisted of two phases: a Prestress Phase (varied to represent different histories of load) followed immediately by an Exertion Phase, identical across all Prestress groups. Three analytical models were used to fit the experimental data in order to evaluate load history effects on gross and zonal disc mechanics. Model results indicated that while gross transient response was insensitive to load history, there may be changes in the internal mechanics of the disc. In particular, a fluid transport model suggested that the role of the nucleus pulposus in resisting creep during Exertion depended on Prestress conditions. Separate experiments using similarly defined load history regimens were performed to verify these predictions by measuring intradiscal pressure with a fibre optic sensor. We found that the ability for intradiscal pressure generation was load history-dependent and exhibited even greater sensitivity than predicted by analytical models. A 0.5?MPa Exertion load resulted in 537.2?kPa IDP for low magnitude Prestress compared with 373.7?kPa for high magnitude Prestress. Based on these measurements, we developed a simple model that may describe the pressure-shear environment in the nucleus pulposus. These findings may have important implications on our understanding of how mechanical stress contributes to disc health and disease etiology.     

Biomechanical effect after Coflex and Coflex rivet implantation for segmental instability at surgical and adjacent segments: a finite element analysis

The Coflex device may provide stability to the surgical segment in extension but does not restore stability in other motion. Recently, a modified version called the Coflex rivet has been developed. The effects of Coflex and Coflex rivet implantation on the adjacent segments are still not clear; therefore, the purpose of this study was to investigate the biomechanical differences between Coflex and Coflex rivet implantation by using finite element analyses. The results show that the Coflex implantation can provide stability in extension, lateral bending, and axial rotation at the surgical segment, and it had no influence at adjacent segments except for extension. The Coflex rivet implantation can provide stability in all motions and reduce disc annulus stress at the surgical segment. Therefore, the higher range of motion and stress induced by the Coflex rivet at both adjacent discs may result in adjacent segment degeneration in flexion and extension.