Mechanical restoration and failure analyses of a hydrogel and scaffold composite strategy for annulus fibrosus repair

Unrepaired defects in the annulus fibrosus of intervertebral disks are associated with degeneration and persistent back pain. A clinical need exists for a disk repair strategy that can seal annular defects, be easily delivered during surgical procedures, and restore biomechanics with low risk of herniation. Multiple annulus repair strategies were developed using poly(trimethylene carbonate) scaffolds optimized for cell delivery, polyurethane membranes designed to prevent herniation, and fibrin-genipin adhesive tuned to annulus fibrosus shear properties. This three-part study evaluated repair strategies for biomechanical restoration, herniation risk and failure mode in torsion, bending and compression at physiological and hyper-physiological loads using a bovine injury model. Fibrin-genipin hydrogel restored some torsional stiffness, bending ROM and disk height loss, with negligible herniation risk and failure was observed histologically at the fibrin-genipin mid-substance following rigorous loading. Scaffold-based repairs partially restored biomechanics, but had high herniation risk even when stabilized with sutured membranes and failure was observed histologically at the interface between scaffold and fibrin-genipin adhesive. Fibrin-genipin was the simplest annulus fibrosus repair solution evaluated that involved an easily deliverable adhesive that filled irregularly-shaped annular defects and partially restored disk biomechanics with low herniation risk, suggesting further evaluation for disk repair may be warranted. Statement of significance Lower back pain is the leading cause of global disability and commonly caused by defects and failure of intervertebral disk tissues resulting in herniation and compression of adjacent nerves. Annulus fibrosus repair materials and techniques have not been successful due to the challenging mechanical and chemical microenvironment and the needs to restore biomechanical behaviors and promote healing with negligible herniation risk while being delivered during surgical procedures. This work addressed this challenging biomaterial and clinical problem using novel materials including an adhesive hydrogel, a scaffold capable of cell delivery, and a membrane to prevent herniation. Composite repair strategies were evaluated and optimized in quantitative three-part study that rigorously evaluated disk repair and provided a framework for evaluating alternate repair techniques.

Biomechanical analysis of torsion and shear forces in lumbar and lumbosacral spine segments of nonchondrodystrophic dogs

OBJECTIVE: To determine stiffness and load-displacement curves as a biomechanical response to applied torsion and shear forces in cadaveric canine lumbar and lumbosacral specimens. STUDY DESIGN: Biomechanical study. ANIMALS: Caudal lumbar and lumbosacral functional spine units (FSU) of nonchondrodystrophic large-breed dogs (n=31) with radiographically normal spines. METHODS: FSU from dogs without musculoskeletal disease were tested in torsion in a custom-built spine loading simulator with 6 degrees of freedom, which uses orthogonally mounted electric motors to apply pure axial rotation. For shear tests, specimens were mounted to a custom-made shear-testing device, driven by a servo hydraulic testing machine. Load-displacement curves were recorded for torsion and shear. RESULTS: Left and right torsion stiffness was not different within each FSU level; however, torsional stiffness of L7-S1 was significantly smaller compared with lumbar FSU (L4-5-L6-7). Ventral/dorsal stiffness was significantly different from lateral stiffness within an individual FSU level for L5-6, L6-7, and L7-S1 but not for L4-5. When the data from 4 tested shear directions from the same specimen were pooled, level L5-6 was significantly stiffer than L7-S1. CONCLUSIONS: Increased range of motion of the lumbosacral joint is reflected by an overall decreased shear and rotational stiffness at the lumbosacral FSU. CLINICAL RELEVANCE: Data from dogs with disc degeneration have to be collected, analyzed, and compared with results from our chondrodystrophic large-breed dogs with radiographically normal spines.

Papain-induced in vitro disc degeneration model for the study of injectable nucleus pulposus therapy

BACKGROUND CONTEXT Proteolytic enzyme digestion of the intervertebral disc (IVD) offers a method to simulate a condition of disc degeneration for the study of cell-scaffold constructs in the degenerated disc. PURPOSE To characterize an in vitro disc degeneration model (DDM) of different severities of glycosaminoglycans (GAG) and water loss by using papain, and to determine the initial response of the human mesenchymal stem cells (MSCs) introduced into this DDM. STUDY DESIGN Disc degeneration model of a bovine disc explant with an end plate was induced by the injection of papain at various concentrations. Labeled MSCs were later introduced in this model. METHODS Phosphate-buffered saline (PBS control) or papain in various concentrations (3, 15, 30, 60, and 150 U/mL) were injected into the bovine caudal IVD explants. Ten days after the injection, GAG content of the discs was evaluated by dimethylmethylene blue assay and cell viability was determined by live/dead staining together with confocal microscopy. Overall matrix composition was evaluated by histology, and water content was visualized by magnetic resonance imaging. Compressive and torsional stiffness of the DDM were also recorded. In the second part, MSCs were labeled with a fluorescence cell membrane tracker and injected into the nucleus of the DDM or a PBS control. Mesenchymal stem cell viability and distribution were evaluated by confocal microscopy. RESULTS A large drop of GAG and water content of the bovine disc were obtained by injecting >30 U/mL papain. Magnetic resonance imaging showed Grade II, III, and IV disc degeneration by injecting 30, 60, and 150 U/mL papain. A cavity in the center of the disc could facilitate later injection of the nucleus pulposus tissue engineering construct while retaining an intact annulus fibrosus. The remaining disc cell viability was not affected. Mesenchymal stem cells injected into the protease-treated DDM disc showed significantly higher cell viability than when injected into the PBS-injected control disc. CONCLUSIONS By varying the concentration of papain for injection, an increasing amount of GAG and water loss could be induced to simulate the different severities of disc degeneration. MSC suspension introduced into the disc has a very low short-term survival. However, it should be clear that this bovine IVD DDM does not reflect a clinical situation but offers exciting possibilities to test novel tissue engineering protocols.

Clinical versus pre-clinical FE models for vertebral body strength predictions

The finite element analysis is an accepted method to predict vertebral body compressive strength. This study compares measurements obtained from in vitro tests with the ones from two different simulation models: clinical quantitative computer tomography (QCT) based homogenized finite element (hFE) models and pre-clinical high-resolution peripheral QCT-based (HR-pQCT) hFE models. About 37 vertebral body sections were prepared by removing end-plates and posterior elements, scanned with QCT (390/450μm voxel size) as well as HR-pQCT (82μm voxel size), and tested in compression up to failure. Non-linear viscous damage hFE models were created from QCT/HT-pQCT images and compared to experimental results based on stiffness and ultimate load. As expected, the predictability of QCT/HR-pQCT-based hFE models for both apparent stiffness (r2=0.685/0.801r2=0.685/0.801) and strength (r2=0.774/0.924r2=0.774/0.924) increased if a better image resolution was used. An analysis of the damage distribution showed similar damage locations for all cases. In conclusion, HR-pQCT-based hFE models increased the predictability considerably and do not need any tuning of input parameters. In contrast, QCT-based hFE models usually need some tuning but are clinically the only possible choice at the moment.

High resolution quantitative computed tomography-based assessment of trabecular microstructure and strength estimates by finite-element analysis of the spine, but not DXA, reflects vertebral fracture status in men with glucocorticoid-induced osteoporosis

High-resolution quantitative computed tomography (HRQCT)-based analysis of spinal bone density and microstructure, finite element analysis (FEA), and DXA were used to investigate the vertebral bone status of men with glucocorticoid-induced osteoporosis (GIO). DXA of L1–L3 and total hip, QCT of L1–L3, and HRQCT of T12 were available for 73 men (54.6±14.0years) with GIO. Prevalent vertebral fracture status was evaluated on radiographs using a semi-quantitative (SQ) score (normal=0 to severe fracture=3), and the spinal deformity index (SDI) score (sum of SQ scores of T4 to L4 vertebrae). Thirty-one (42.4%) subjects had prevalent vertebral fractures. Cortical BMD (Ct.BMD) and thickness (Ct.Th), trabecular BMD (Tb.BMD), apparent trabecular bone volume fraction (app.BV/TV), and apparent trabecular separation (app.Tb.Sp) were analyzed by HRQCT. Stiffness and strength of T12 were computed by HRQCT-based nonlinear FEA for axial compression, anterior bending and axial torsion. In logistic regressions adjusted for age, glucocorticoid dose and osteoporosis treatment, Tb.BMD was most closely associated with vertebral fracture status (standardized odds ratio [sOR]: Tb.BMD T12: 4.05 [95% CI: 1.8–9.0], Tb.BMD L1–L3: 3.95 [1.8–8.9]). Strength divided by cross-sectional area for axial compression showed the most significant association with spine fracture status among FEA variables (2.56 [1.29–5.07]). SDI was best predicted by a microstructural model using Ct.Th and app.Tb.Sp (r2=0.57, p<0.001). Spinal or hip DXA measurements did not show significant associations with fracture status or severity. In this cross-sectional study of males with GIO, QCT, HRQCT-based measurements and FEA variables were superior to DXA in discriminating between patients of differing prevalent vertebral fracture status. A microstructural model combining aspects of cortical and trabecular bone reflected fracture severity most accurately.

Finite Element Based Nonlinear Normalization of Human Lumbar Intervertebral Disc Stiffness to Account for Its Morphology

Disc degeneration, usually associated with low back pain and changes of intervertebral stiffness, represents a major health issue. As the intervertebral disc (IVD) morphology influences its stiffness, the link between mechanical properties and degenerative grade is partially lost without an efficient normalization of the stiffness with respect to the morphology. Moreover, although the behavior of soft tissues is highly nonlinear, only linear normalization protocols have been defined so far for the disc stiffness. Thus, the aim of this work is to propose a nonlinear normalization based on finite elements (FE) simulations and evaluate its impact on the stiffness of human anatomical specimens of lumbar IVD. First, a parameter study involving simulations of biomechanical tests (compression, flexion/extension, bilateral torsion and bending) on 20 FE models of IVDs with various dimensions was carried out to evaluate the effect of the disc's geometry on its compliance and establish stiffness/morphology relations necessary to the nonlinear normalization. The computed stiffness was then normalized by height (H), cross-sectional area (CSA), polar moment of inertia (J) or moments of inertia (Ixx, Iyy) to quantify the effect of both linear and nonlinear normalizations. In the second part of the study, T1-weighted MRI images were acquired to determine H, CSA, J, Ixx and Iyy of 14 human lumbar IVDs. Based on the measured morphology and pre-established relation with stiffness, linear and nonlinear normalization routines were then applied to the compliance of the specimens for each quasi-static biomechanical test. The variability of the stiffness prior to and after normalization was assessed via coefficient of variation (CV). The FE study confirmed that larger and thinner IVDs were stiffer while the normalization strongly attenuated the effect of the disc geometry on its stiffness. Yet, notwithstanding the results of the FE study, the experimental stiffness showed consistently higher CV after normalization. Assuming that geometry and material properties affect the mechanical response, they can also compensate for one another. Therefore, the larger CV after normalization can be interpreted as a strong variability of the material properties, previously hidden by the geometry's own influence. In conclusion, a new normalization protocol for the intervertebral disc stiffness in compression, flexion, extension, bilateral torsion and bending was proposed, with the possible use of MRI and FE to acquire the discs' anatomy and determine the nonlinear relations between stiffness and morphology. Such protocol may be useful to relate the disc's mechanical properties to its degree of degeneration.

Inactive Matrix Gla-Protein Is Associated With Arterial Stiffness in an Adult Population-Based Study

Increased pulse wave velocity (PWV) is a marker of aortic stiffness and an independent predictor of mortality. Matrix Gla-protein (MGP) is a vascular calcification inhibitor that needs vitamin K to be activated. Inactive MGP, known as desphospho-uncarboxylated MGP (dp-ucMGP), can be measured in plasma and has been associated with various cardiovascular markers, cardiovascular outcomes, and mortality. In this study, we hypothesized that high levels of dp-ucMGP are associated with increased PWV. We recruited participants via a multicenter family-based cross-sectional study in Switzerland. Dp-ucMGP was quantified in plasma by sandwich ELISA. Aortic PWV was determined by applanation tonometry using carotid and femoral pulse waveforms. Multiple regression analysis was performed to estimate associations between PWV and dp-ucMGP adjusting for age, renal function, and other cardiovascular risk factors. We included 1001 participants in our analyses (475 men and 526 women). Mean values were 7.87±2.10 m/s for PWV and 0.43±0.20 nmol/L for dp-ucMGP. PWV was positively associated with dp-ucMGP both before and after adjustment for sex, age, body mass index, height, systolic and diastolic blood pressure (BP), heart rate, renal function, low- and high-density lipoprotein, glucose, smoking status, diabetes mellitus, BP and cholesterol lowering drugs, and history of cardiovascular disease (P≤0.01). In conclusion, high levels of dp-ucMGP are independently and positively associated with arterial stiffness after adjustment for common cardiovascular risk factors, renal function, and age. Experimental studies are needed to determine whether vitamin K supplementation slows arterial stiffening by increasing MGP carboxylation.

Biomechanical model of human cornea based on stromal microstructure

The optical characteristics of the human cornea depends on the mechanical balance between the intra-ocular pressure and intrinsic tissue stiffness. A wide range of ophthalmic surgical procedures alter corneal biomechanics to induce local or global curvature changes for the correction of visual acuity. Due to the large number of surgical interventions performed every day, a deeper understanding of corneal biomechanics is needed to improve the safety of these procedures and medical devices. The aim of this study is to propose a biomechanical model of the human cornea, based on stromal microstructure. The constitutive mechanical law includes collagen fiber distribution based on X-ray scattering analysis, collagen cross-linking, and fiber uncrimping. Our results showed that the proposed model reproduced inflation and extensiometry experimental data [Elsheikh et al., Curr. Eye Res., 2007; Elsheikh et al., Exp. Eye Res., 2008] successfully. The mechanical properties obtained for different age groups demonstrated an increase in collagen cross-linking for older specimens. In future work such a model could be used to simulate non-symmetric interventions, and provide better surgical planning.

Numerical model of the human cornea based on stromal microstructure: identification of mechanical properties for two age groups

The optical quality of the human eye mainly depends on the refractive performance of the cornea. The shape of the cornea is a mechanical balance between intraocular pressure and tissue intrinsic stiffness. Several surgical procedures in ophthalmology alter the biomechanics of the cornea to provoke local or global curvature changes for vision correction. Legitimated by the large number of surgical interventions performed every day, the demand for a deeper understanding of corneal biomechanics is rising to improve the safety of procedures and medical devices. The aim of our work is to propose a numerical model of corneal biomechanics, based on the stromal microstructure. Our novel anisotropic constitutive material law features a probabilistic weighting approach to model collagen fiber distribution as observed on human cornea by Xray scattering analysis (Aghamohammadzadeh et. al., Structure, February 2004). Furthermore, collagen cross-linking was explicitly included in the strain energy function. Results showed that the proposed model is able to successfully reproduce both inflation and extensiometry experimental data (Elsheikh et. al., Curr Eye Res, 2007; Elsheikh et. al., Exp Eye Res, May 2008). In addition, the mechanical properties calculated for patients of different age groups (Group A: 65-79 years; Group B: 80-95 years) demonstrate an increased collagen cross-linking, and a decrease in collagen fiber elasticity from younger to older specimen. These findings correspond to what is known about maturing fibrous biological tissue. Since the presented model can handle different loading situations and includes the anisotropic distribution of collagen fibers, it has the potential to simulate clinical procedures involving nonsymmetrical tissue interventions. In the future, such mechanical model can be used to improve surgical planning and the design of next generation ophthalmic devices.

Influence of Smoothing on Voxel-Based Mesh Accuracy in Micro-Finite Element

The interest in automatic volume meshing for finite element analysis (FEA) has grown more since the appearance of microfocus CT (μCT), due to its high resolution, which allows for the assessment of mechanical behaviour at a high precision. Nevertheless, the basic meshing approach of generating one hexahedron per voxel produces jagged edges. To prevent this effect, smoothing algorithms have been introduced to enhance the topology of the mesh. However, whether smoothing also improves the accuracy of voxel-based meshes in clinical applications is still under question. There is a trade-off between smoothing and quality of elements in the mesh. Distorted elements may be produced by excessive smoothing and reduce accuracy of the mesh. In the present work, influence of smoothing on the accuracy of voxel-based meshes in micro-FE was assessed. An accurate 3D model of a trabecular structure with known apparent mechanical properties was used as a reference model. Virtual CT scans of this reference model (with resolutions of 16, 32 and 64 μm) were then created and used to build voxel-based meshes of the microarchitecture. Effects of smoothing on the apparent mechanical properties of the voxel-based meshes as compared to the reference model were evaluated. Apparent Young’s moduli of the smooth voxel-based mesh were significantly closer to those of the reference model for the 16 and 32 μm resolutions. Improvements were not significant for the 64 μm, due to loss of trabecular connectivity in the model. This study shows that smoothing offers a real benefit to voxel-based meshes used in micro-FE. It might also broaden voxel-based meshing to other biomechanical domains where it was not used previously due to lack of accuracy. As an example, this work will be used in the framework of the European project ContraCancrum, which aims at providing a patient-specific simulation of tumour development in brain and lungs for oncologists. For this type of clinical application, such a fast, automatic, and accurate generation of the mesh is of great benefit.

Pulse contour analysis: a valid assessment of central arterial stiffness in children?

In adults the contour analysis of peripheral pressure waves in the upper limb reflects central aortic stiffness. Here, we wanted to demonstrate the appropriateness of pulse contour analysis to assess large artery stiffness in children. Digital volume pulse analysis, with the computation of the stiffness index and pulse wave velocity between carotid and femoral artery, were simultaneously determined in 79 healthy children between 8 years and 15 years (mean age 11.4 years, 32 girls). The stiffness index of 42 healthy adults (mean age 45.6 years, 26 women) served as control. Pulse wave velocity between carotid and femoral artery was directly correlated with systolic pressure and mean blood pressure, as well as with pulse pressure. The results from the stiffness index of children revealed the expected values extrapolated from the linear regression of adulthood stiffness index vs. age. Childhood stiffness index positively correlated with pulse wave velocity (r(2) = 0.07, P = 0.02) but not with blood pressure parameters. The exclusion of individuals with an increased vascular tone, as indicated by a reflexion index > 90%, improved the correlation between stiffness index and pulse wave velocity (r(2) = 0.13, P = 0.001). Our data indicate that digital volume pulse-based analysis has limitations if compared with pulse wave velocity to measure arterial stiffness, mostly in patients with a high vascular tone.

Cement distribution, volume, and compliance in vertebroplasty: some answers from an anatomy-based nonlinear finite element study

STUDY DESIGN: The biomechanics of vertebral bodies augmented with real distributions of cement were investigated using nonlinear finite element (FE) analysis. OBJECTIVES: To compare stiffness, strength, and stress transfer of augmented versus nonaugmented osteoporotic vertebral bodies under compressive loading. Specifically, to examine how cement distribution, volume, and compliance affect these biomechanical variables. SUMMARY OF BACKGROUND DATA: Previous FE studies suggested that vertebroplasty might alter vertebral stress transfer, leading to adjacent vertebral failure. However, no FE study so far accounted for real cement distributions and bone damage accumulation. METHODS: Twelve vertebral bodies scanned with high-resolution pQCT and tested in compression were augmented with various volumes of cements and scanned again. Nonaugmented and augmented pQCT datasets were converted to FE models, with bone properties modeled with an elastic, plastic and damage constitutive law that was previously calibrated for the nonaugmented models. The cement-bone composite was modeled with a rule of mixture. The nonaugmented and augmented FE models were subjected to compression and their stiffness, strength, and stress map calculated for different cement compliances. RESULTS: Cement distribution dominated the stiffening and strengthening effects of augmentation. Models with cement connecting either the superior or inferior endplate (S/I fillings) were only up to 2 times stiffer than the nonaugmented models with minimal strengthening, whereas those with cement connecting both endplates (S + I fillings) were 1 to 8 times stiffer and 1 to 12 times stronger. Stress increases above and below the cement, which was higher for the S + I cases and was significantly reduced by increasing cement compliance. CONCLUSION: The developed FE approach, which accounts for real cement distributions and bone damage accumulation, provides a refined insight into the mechanics of augmented vertebral bodies. In particular, augmentation with compliant cement bridging both endplates would reduce stress transfer while providing sufficient strengthening.

Human intervertebral disc stiffness correlates better with the Otsu threshold computed from axial T2 map of its posterior annulus fibrosus than with clinical classifications

Degeneration of the intervertebral disc, sometimes associated with low back pain and abnormal spinal motions, represents a major health issue with high costs. A non-invasive degeneration assessment via qualitative or quantitative MRI (magnetic resonance imaging) is possible, yet, no relation between mechanical properties and T2 maps of the intervertebral disc (IVD) has been considered, albeit T2 relaxation time values quantify the degree of degeneration. Therefore, MRI scans and mechanical tests were performed on 14 human lumbar intervertebral segments freed from posterior elements and all soft tissues excluding the IVD. Degeneration was evaluated in each specimen using morphological criteria, qualitative T2 weighted images and quantitative axial T2 map data and stiffness was calculated from the load-deflection curves of in vitro compression, torsion, lateral bending and flexion/extension tests. In addition to mean T2, the OTSU threshold of T2 (TOTSU), a robust and automatic histogram-based method that computes the optimal threshold maximizing the distinction of two classes of values, was calculated for anterior, posterior, left and right regions of each annulus fibrosus (AF). While mean T2 and degeneration schemes were not related to the IVDs' mechanical properties, TOTSU computed in the posterior AF correlated significantly with those classifications as well as with all stiffness values. TOTSU should therefore be included in future degeneration grading schemes.

A nonlinear QCT-based finite element model validation study for the human femur tested in two configurations in vitro

Purpose Femoral fracture is a common medical problem in osteoporotic individuals. Bone mineral density (BMD) is the gold standard measure to evaluate fracture risk in vivo. Quantitative computed tomography (QCT)-based homogenized voxel finite element (hvFE) models have been proved to be more accurate predictors of femoral strength than BMD by adding geometrical and material properties. The aim of this study was to evaluate the ability of hvFE models in predicting femoral stiffness, strength and failure location for a large number of pairs of human femora tested in two different loading scenarios. Methods Thirty-six pairs of femora were scanned with QCT and total proximal BMD and BMC were evaluated. For each pair, one femur was positioned in one-legged stance configuration (STANCE) and the other in a sideways configuration (SIDE). Nonlinear hvFE models were generated from QCT images by reproducing the same loading configurations imposed in the experiments. For experiments and models, the structural properties (stiffness and ultimate load), the failure location and the motion of the femoral head were computed and compared. Results In both configurations, hvFE models predicted both stiffness (R2=0.82 for STANCE and R2=0.74 for SIDE) and femoral ultimate load (R2=0.80 for STANCE and R2=0.85 for SIDE) better than BMD and BMC. Moreover, the models predicted qualitatively well the failure location (66% of cases) and the motion of the femoral head. Conclusions The subject specific QCT-based nonlinear hvFE model cannot only predict femoral apparent mechanical properties better than densitometric measures, but can additionally provide useful qualitative information about failure location.