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.

Combining 3D tracking and surgical instrumentation to determine the stiffness of spinal motion segments: a validation study

The spine is a complex structure that provides motion in three directions: flexion and extension, lateral bending and axial rotation. So far, the investigation of the mechanical and kinematic behavior of the basic unit of the spine, a motion segment, is predominantly a domain of in vitro experiments on spinal loading simulators. Most existing approaches to measure spinal stiffness intraoperatively in an in vivo environment use a distractor. However, these concepts usually assume a planar loading and motion. The objective of our study was to develop and validate an apparatus, that allows to perform intraoperative in vivo measurements to determine both the applied force and the resulting motion in three dimensional space. The proposed setup combines force measurement with an instrumented distractor and motion tracking with an optoelectronic system. As the orientation of the applied force and the three dimensional motion is known, not only force-displacement, but also moment-angle relations could be determined. The validation was performed using three cadaveric lumbar ovine spines. The lateral bending stiffness of two motion segments per specimen was determined with the proposed concept and compared with the stiffness acquired on a spinal loading simulator which was considered to be gold standard. The mean values of the stiffness computed with the proposed concept were within a range of ±15% compared to data obtained with the spinal loading simulator under applied loads of less than 5 Nm.

Stochastic amplitude-modulated stretching of Rabbit flexor digitorum profundus tendons reduces stiffness compared to Cyclic Loading but does not affect Tenocyte Metabolism

Background It has been demonstrated that frequency modulation of loading influences cellular response and metabolism in 3D tissues such as cartilage, bone and intervertebral disc. However, the mechano-sensitivity of cells in linear tissues such as tendons or ligaments might be more sensitive to changes in strain amplitude than frequency. Here, we hypothesized that tenocytes in situ are mechano-responsive to random amplitude modulation of strain. Methods We compared stochastic amplitude-modulated versus sinusoidal cyclic stretching. Rabbit tendon were kept in tissue-culture medium for twelve days and were loaded for 1h/day for six of the total twelve culture days. The tendons were randomly subjected to one of three different loading regimes: i) stochastic (2 – 7% random strain amplitudes), ii) cyclic_RMS (2–4.42% strain) and iii) cyclic_high (2 - 7% strain), all at 1 Hz and for 3,600 cycles, and one unloaded control. Results At the end of the culture period, the stiffness of the “stochastic” group was significantly lower than that of the cyclic_RMS and cyclic_high groups (both, p < 0.0001). Gene expression of eleven anabolic, catabolic and inflammatory genes revealed no significant differences between the loading groups. Conclusions We conclude that, despite an equivalent metabolic response, stochastically stretched tendons suffer most likely from increased mechanical microdamage, relative to cyclically loaded ones, which is relevant for tendon regeneration therapies in clinical practice.

Mechanical and radiological assessment of the influence of rhTGFbeta-3 on bone regeneration in a segmental defect in the ovine tibia: pilot study

Limitations in the use of autologous bone graft, which is the gold standard therapy in bone defect healing, drive the search for alternative treatments. In this study the influence of rhTGFbeta-3 on mechanical and radiological parameters of a healing bone defect in the sheep tibia was assessed. In the sheep, an 18-mm long osteoperiosteal defect in the tibia was treated by rhTGFbeta-3 seeded on a poly(L/DL-lactide) carrier (n = 4). In a second group (n = 4), the defect was treated by the carrier only, in a third group (n = 4) by autologous cancellous bone graft, and in a fourth group (n = 2) the defect remained blank. The healing process of the defect was assessed by weekly in vivo stiffness measurements and radiology as well as by quantitative computed tomographic assessment of bone mineral density (BMD) every 4 weeks. The duration of the experiment was 12 weeks under loading conditions. In the bone graft group, a marginally significant higher increase in stiffness was observed than in the PLA/rhTGFbeta-3 group (p = 0.06) and a significantly higher increase than in the PLA-only group (p = 0.03). The radiographic as well as the computed tomographic evaluation yielded significant differences between the groups (p = 0.03), indicating the bone graft treatment (bone/per area, 83%; BMD, 0.57 g/cm(3)) performing better than the PLA/rhTGFbeta-3 (38%; 0.23 g/cm(3)) and the PLA-only treatment (2.5%; 0.09 g/cm(3)), respectively. Regarding the mechanical and radiological parameters assessed in this study, we conclude that rhTGFbeta-3 has a promoting effect on bone regeneration. However, under the conditions of this study, this effect does not reach the potential of autologous cancellous bone graft transplantation.

Variation of the mechanical properties of PMMA to suit osteoporotic cancellous bone

Poly(methyl methacrylate) (PMMA) is by far the most frequently used bone substitute material for vertebroplasty. However, there are serious complications, such as cement leakage and an increased fracture rate of the adjacent vertebral bodies. The latter may be related to the mechanical properties of the augmented segment within the osteoporotic spine. A possible counter-measure is prophylactic augmentation at additional levels, but this aggravates the risk for the patient. Introduction of pores is a possible method to reduce the inherent high stiffness of PMMA. This study investigates the effect of porosity on the mechanical properties of PMMA bone cement. Different fractions of a highly viscous liquid were mixed into the PMMA during preparation. An open-porous material with adjustable mechanical properties resulted after removal of the aqueous phase. Different radiopacifiers were admixed to investigate their suitability for vertebroplasty. The final material was characterized mechanically by compressive testing, microscopically and radiologically. In addition, the monomer release subsequent to hardening was measured by means of gas chromatography. The Young's modulus in compression could be varied between 2800 +/- 70 MPa and 120 +/- 150 MPa, and the compression ultimate strength between 170 +/- 5 MPa and 8 +/- 9 MPa for aqueous fractions ranging between 0 and 50% of volume. Only a slight decrease of the Young's modulus and small changes of ultimate strength were found when the mixing time was increased. An organic hydrophilic and lipophilic radiopacifier led to a higher Young's modulus of the porous material; however, the ultimate strength was not significantly affected by adding different radiopacifiers to the porous cement. The radiopacity was lost after washing the aqueous phase out of the pores. No separation occurred between the aqueous and the PMMA phase during injection into an open porous ceramic material. The monomer released was found to increase for increasing aqueous fractions, but remained comparable in magnitude to standard PMMA. This study demonstrates that a conventional PMMA can be modified to obtain a range of mechanical properties, including those of osteoporotic bone.

Augmentation of mechanical properties in osteoporotic vertebral bones--a biomechanical investigation of vertebroplasty efficacy with different bone cements

Recent clinical trials have reported favorable early results for transpedicular vertebral cement reinforcement of osteoporotic vertebral insufficiencies. There is, however, a lack of basic data on the application, safety and biomechanical efficacy of materials such as polymethyl-methacrylate (PMMA) and calciumphospate (CaP) cements. The present study analyzed 33 vertebral pairs from five human cadaver spines. Thirty-nine vertebrae were osteoporotic (bone mineral density < 0.75 g/cm2), 27 showed nearly normal values. The cranial vertebra of each pair was augmented with either PMMA (Palacos E-Flow) or experimental brushite cement (EBC), with the caudal vertebra as a control. PMMA and EBC were easy to inject, and vertebral fillings of 20-50% were achieved. The maximal possible filling was inversely correlated to the bone mineral density (BMD) values. Cement extrusion into the spinal canal was observed in 12% of cases. All specimens were subjected to axial compression tests in a displacement-controlled mode. From load-displacement curves, the stiffness, S, and the maximal force before failure, Fmax, were determined. Compared with the native control vertebrae, a statistically significant increase in vertebral stiffness and Fmax was observed by the augmentation. With PMMA the stiffness increased by 174% (P = 0.018) and Fmax by 195% (P = 0.001); the corresponding augmentation with EBC was 120% (P = 0.03) and 113% (P = 0.002). The lower the initial BMD, the more pronounced was the augmentation effect. Both PMMA and EBC augmentation reliably and significantly raised the stiffness and maximal tolerable force until failure in osteoporotic vertebral bodies. In non-porotic specimens, no significant increase was achieved.

Fabric-mechanical property relationships of trabecular bone allografts are altered by supercritical CO2 treatment and gamma sterilization

Tissue grafts are implanted in orthopedic surgery every day. In order to minimize infection risk, bone allografts are often delipidated with supercritical CO2 and sterilized prior to implantation. This treatment may, however, impair the mechanical behavior of the bone graft tissue. The goal of this study was to determine clinically relevant mechanical properties of treated/sterilized human trabecular bone grafts, e.g. the apparent modulus, strength, and the ability to absorb energy during compaction. They were compared with results of identical experiments performed previously on untreated/fresh frozen human trabecular bone from the same anatomical site (Charlebois, 2008). We tested the hypothesis that the morphology–mechanical property relationships of treated cancellous allografts are similar to those of fresh untreated bone. The morphology of the allografts was determined by μCT. Subsequently, cylindrical samples were tested in unconfined and confined compression. To account for various morphologies, the experimental data was fitted to phenomenological mechanical models for elasticity, strength, and dissipated energy density based on bone volume fraction (BV/TV) and the fabric tensor determined by MIL. The treatment/sterilization process does not appear to influence bone graft stiffness. However, strength and energy dissipation of the bone grafts were found to be significantly reduced by 36% to 47% and 66% to 81%, respectively, for a broad range of volume fraction (0.14 < BV/TV < 0.39) and degree of anisotropy (1.24 < DA < 2.18). Since the latter properties are strongly dominated by BV/TV, the clinical consequences of this reduction can be compensated by using grafts with lower porosity. The data of this study suggests that an increase of 5–10% in BV/TV is sufficient to compensate for the reduced post-yield mechanical properties of treated/sterilized bone in monotonic compression. In applications where graft stiffness needs to be matched and strength is not a concern, treated allograft with the same BV/TV as an appropriate fresh bone graft may be used.

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.

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.

Patient-specific spinal stiffness in AIS: a preoperative and noninvasive method.

INTRODUCTION The clinical tests currently used to assess spinal biomechanics preoperatively are unable to assess true mechanical spinal stiffness. They rely on spinal displacement without considering the force required to deform a patient's spine. We propose a preoperative method for noninvasively quantifying the three-dimensional patient-specific stiffness of the spines of adolescent idiopathic scoliosis patients. METHODS The technique combines a novel clinical test with numerical optimization of a finite element model of the patient's spine. RESULTS A pilot study conducted on five patients showed that the model was able to provide accurate 3D reconstruction of the spine's midline and predict the spine's stiffness for each patient in flexion, bending, and rotation. Statistically significant variation of spinal stiffness was observed between the patients. CONCLUSION This result confirms that spinal biomechanics is patient-specific, which should be taken into consideration to individualize surgical treatment.

Bone Volume Fraction and Fabric Anisotropy Are Better Determinants of Trabecular Bone Stiffness than Other Morphological Variables

As our population ages, more individuals suffer from osteoporosis. This disease leads to impaired trabecular architecture and increased fracture risk. It is essential to understand how morphological and mechanical properties of the cancellous bone are related. Morphologyelasticity relationships based on bone volume fraction (BV/TV) and fabric anisotropy explain up to 98% of the variation in elastic properties. Yet, other morphological variables such as individual trabeculae segmentation (ITS) and trabecular bone score (TBS) could improve the stiffness predictions. A total of 743 micro-computed tomography reconstructions of cubic trabecular bone samples extracted from femur, radius, vertebrae and iliac crest were analysed. Their morphology was assessed via 25 variables and their stiffness tensor (inline image) was computed from six independent load cases using micro finite element analyses. Variance inflation factors were calculated to evaluate collinearity between morphological variables and decide upon their inclusion in morphology-elasticity relationships. The statistically admissible morphological variables were included in a multi-linear regression modelling the dependent variable inline image. The contribution of each independent variable was evaluated (ANOVA). Our results show that BV/TV is the best determinant of inline image (inline image=0.889), especially in combination with fabric (inline image=0.968). Including the other independent predictors hardly affected the amount of variance explained by the model (inline image=0.975). Across all anatomical sites, BV/TV explained 87% of the variance of the bone elastic properties. Fabric further described 10% of the bone stiffness, but the improvement in variance explanation by adding other independent factors was marginal (<1%). These findings confirm that BV/TV and fabric are the best determinants of trabecular bone stiffness and show, against common belief, that other morphological variables do not bring any further contribution. These overall conclusions remain to be confirmed for specific bone diseases and post-elastic properties.

Measuring the mechanical properties of plant cells by combining micro-indentation with osmotic treatments

Growth in plants results from the interaction between genetic and signalling networks and the mechanical properties of cells and tissues. There has been a recent resurgence in research directed at understanding the mechanical aspects of growth, and their feedback on genetic regulation. This has been driven in part by the development of new micro-indentation techniques to measure the mechanical properties of plant cells in vivo. However, the interpretation of indentation experiments remains a challenge, since the force measures results from a combination of turgor pressure, cell wall stiffness, and cell and indenter geometry. In order to interpret the measurements, an accurate mechanical model of the experiment is required. Here, we used a plant cell system with a simple geometry, Nicotiana tabacum Bright Yellow-2 (BY-2) cells, to examine the sensitivity of micro-indentation to a variety of mechanical and experimental parameters. Using a finite-element mechanical model, we found that, for indentations of a few microns on turgid cells, the measurements were mostly sensitive to turgor pressure and the radius of the cell, and not to the exact indenter shape or elastic properties of the cell wall. By complementing indentation experiments with osmotic experiments to measure the elastic strain in turgid cells, we could fit the model to both turgor pressure and cell wall elasticity. This allowed us to interpret apparent stiffness values in terms of meaningful physical parameters that are relevant for morphogenesis.

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.