Structural characterization and reliable biomechanical assessment of integrative cartilage repair

Structural and functional characterization of integrative cartilage repair in controlled model systems can play a key role in the development of innovative strategies to improve the long-term outcome of many cartilage repair procedures. In this work, we first developed a method to reproducibly generate geometrically defined disk/ring cartilage composites and to remove outgrown fibrous layers which can encapsulate cartilaginous tissues during culture. We then used the model system to test the hypothesis that such fibrous layers lead to an overestimation of biomechanical parameters of integration at the disk/ring interface. Transmission electron microscopy images of the composites after 6 weeks of culture indicated that collagen fibrils in the fibrous tissue layer were well integrated into the collagen network of the cartilage disk and ring, whereas molecular bridging between opposing disk/ring cartilage surfaces was less pronounced and restricted to regions with narrow interfacial regions (< 2 microm). Stress-strain profiles generated from mechanical push-out tests for composites with the layers removed displayed a single and distinct peak, whereas profiles for composites with the layers left intact consisted of multiple superimposed peaks. As compared to composites with removed layers, composites with intact layers had significantly higher adhesive strengths (161+/-9 vs. 71+/-11 kPa) and adhesion energies (15.0+/-0.7 vs. 2.7+/-0.4 mJ/mm2). By combining structural and functional analyses, we demonstrated that the outgrowing tissue formed during in vitro culture of cartilaginous specimens should be eliminated in order to reliably quantify biomechanical parameters related to integrative cartilage repair.

Three-dimensional primary stability of cementless femoral stems

OBJECTIVE: This study investigates by means of a new bone-prosthesis interface motion detector whether conceptual design differences of femoral stems are reflected in their primary stability pattern. DESIGN: An in vitro experiment using a biaxial materials testing machine in combination with three-dimensional motion measurement devices was performed. BACKGROUND: Primary stability of uncemented total hip replacements is considered to be a prerequisite for the quality of bony ongrowth to the femoral stem. Dynamic motion as a response to loading as well as total motion of the prosthesis have to be considered under quasi-physiological cyclic loading conditions. METHODS: Seven paired fresh cadaveric femora were used for the testing of two types of uncemented femoral stems with different anchoring concepts: CLS stem (Spotorno) and Cone Prosthesis (Wagner). Under sinusoidal cyclic loading mimicking in vivo hip joint forces a new measurement technique was applied allowing for the analysis of the three-dimensional interface motion. RESULTS: Considerable differences between the two prostheses could be detected both in their dynamic motion and total motion behaviour. Whereas the CLS stem, due to the wedge-shaped concept, provides smaller total motions, the longitudinal ribs of the Cone prostheses result in a substantially smaller dynamic motion. CONCLUSIONS: The measuring technique provided reliable and accurate data illustrating the three-dimensional interface motion of uncemented femoral stems.

Maxillary sinus floor extension and posterior tooth inclination in adolescent patients with Class II Division 1 malocclusion treated with maxillary first molar extractions

INTRODUCTION Our objective was to investigate potential associations between maxillary sinus floor extension and inclination of maxillary second premolars and second molars in patients with Class II Division 1 malocclusion whose orthodontic treatment included maxillary first molar extractions. METHODS The records of 37 patients (18 boys, 19 girls; mean age, 13.2 years; SD, 1.62 years) treated between 1998 and 2004 by 1 orthodontist with full Begg appliances were used in this study. Inclusion criteria were white patients with Class II Division 1 malocclusion, sagittal overjet of ≥4 mm, treatment plan including extraction of the maxillary first permanent molars, no missing teeth, and no agenesis. Maxillary posterior tooth inclination and lower maxillary sinus area in relation to the palatal plane were measured on lateral cephalograms at 3 time points: at the start and end of treatment, and on average 2.5 years posttreatment. Data were analyzed for the second premolar and second molar inclinations by using mixed linear models. RESULTS The analysis showed that the second molar inclination angle decreased by 7° after orthodontic treatment, compared with pretreatment values, and by 11.5° at the latest follow-up, compared with pretreatment. There was evidence that maxillary sinus volume was negatively correlated with second molar inclination angle; the greater the volume, the smaller the inclination angle. For premolars, inclination increased by 15.4° after orthodontic treatment compared with pretreatment, and by 8.1° at the latest follow-up compared with baseline. The volume of the maxillary sinus was not associated with premolar inclination. CONCLUSIONS We found evidence of an association between maxillary second molar inclination and surface area of the lower sinus in patients treated with maxillary first molar extractions. Clinicians who undertake such an extraction scheme in Class II patients should be aware of this potential association and consider appropriate biomechanics to control root uprighting.

The influence of muscle strength on the gait profile score (GPS) across different patients

BACKGROUND Muscle strength greatly influences gait kinematics. The question was whether this association is similar in different diseases. METHODS Data from instrumented gait analysis of 716 patients were retrospectively assessed. The effect of muscle strength on gait deviations, namely the gait profile score (GPS) was evaluated by means of generalised least square models. This was executed for seven different patient groups. The groups were formed according to the type of disease: orthopaedic/neurologic, uni-/bilateral affection, and flaccid/spastic muscles. RESULTS Muscle strength had a negative effect on GPS values, which did not significantly differ amongst the different patient groups. However, an offset of the GPS regression line was found, which was mostly dependent on the basic disease. Surprisingly, spastic patients, who have reduced strength and additionally spasticity in clinical examination, and flaccid neurologic patients showed the same offset. Patients with additional lack of trunk control (Tetraplegia) showed the largest offset. CONCLUSION Gait kinematics grossly depend on muscle strength. This was seen in patients with very different pathologies. Nevertheless, optimal correction of biomechanics and muscle strength may still not lead to a normal gait, especially in that of neurologic patients. The basic disease itself has an additional effect on gait deviations expressed as a GPS-offset of the linear regression line.

Nitric oxide induces cell death in canine cruciate ligament cells by activation of tyrosine kinase and reactive oxygen species

BACKGROUND: There is increasing evidence suggesting that development of progressive canine cranial cruciate ligament (CCL) rupture involves a gradual degeneration of the CCL itself, initiated by a combination of factors, ranging from mechanical to biochemical. To date, knowledge is lacking to what extent cruciate disease results from abnormal biomechanics on a normal ligament or contrary how far preliminary alterations of the ligament due to biochemical factors provoke abnormal biomechanics. This study is focused on nitric oxide (NO), one of the potential biochemical factors. The NO-donor sodium nitroprusside (SNP) has been used to study NO-dependent cell death in canine cranial and caudal cruciate ligament cells and to characterize signaling mechanisms during NO-stimulation. RESULTS: Sodium nitroprusside increased apoptotic cell death dose- and time-dependently in cruciate ligamentocytes. Cells from the CCL were more susceptible to apoptosis than CaCL cells. Caspase-3 processing in response to SNP was not detected. Testing major upstream and signal transducing pathways, NO-induced cruciate ligament cell death seemed to be mediated on different levels. Specific inhibition of tyrosine kinase significantly decreased SNP-induced cell death. Mitogen activated protein kinase ERK1 and 2 are activated upon NO and provide anti-apoptotic signals whereas p38 kinase and protein kinase C are not involved. Moreover, data showed that the inhibition reactive oxygen species (ROS) significantly reduced the level of cruciate ligament cell death. CONCLUSIONS: Our data support the hypothesis that canine cruciate ligamentocytes, independently from their origin (CCL or CaCL) follow crucial signaling pathways involved in NO-induced cell death. However, the difference on susceptibility upon NO-mediated apoptosis seems to be dependent on other pathways than on these tested in the present study. In both, CCL and CaCL, the activation of the tyrosine kinase and the generation of ROS reveal important signaling pathways. In perspective, new efforts to prevent the development and progression of cruciate disease may include strategies aimed at reducing ROS.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

A generalized anisotropic quadric yield criterion and its application to bone tissue at multiple length scales

Nonlinear computational analysis of materials showing elasto-plasticity or damage relies on knowledge of their yield behavior and strengths under complex stress states. In this work, a generalized anisotropic quadric yield criterion is proposed that is homogeneous of degree one and takes a convex quadric shape with a smooth transition from ellipsoidal to cylindrical or conical surfaces. If in the case of material identification, the shape of the yield function is not known a priori, a minimization using the quadric criterion will result in the optimal shape among the convex quadrics. The convexity limits of the criterion and the transition points between the different shapes are identified. Several special cases of the criterion for distinct material symmetries such as isotropy, cubic symmetry, fabric-based orthotropy and general orthotropy are presented and discussed. The generality of the formulation is demonstrated by showing its degeneration to several classical yield surfaces like the von Mises, Drucker–Prager, Tsai–Wu, Liu, generalized Hill and classical Hill criteria under appropriate conditions. Applicability of the formulation for micromechanical analyses was shown by transformation of a criterion for porous cohesive-frictional materials by Maghous et al. In order to demonstrate the advantages of the generalized formulation, bone is chosen as an example material, since it features yield envelopes with different shapes depending on the considered length scale. A fabric- and density-based quadric criterion for the description of homogenized material behavior of trabecular bone is identified from uniaxial, multiaxial and torsional experimental data. Also, a fabric- and density-based Tsai–Wu yield criterion for homogenized trabecular bone from in silico data is converted to an equivalent quadric criterion by introduction of a transformation of the interaction parameters. Finally, a quadric yield criterion for lamellar bone at the microscale is identified from a nanoindentation study reported in the literature, thus demonstrating the applicability of the generalized formulation to the description of the yield envelope of bone at multiple length scales.

The influence of yield surface shape and damage in the depth-dependent response of bone tissue to nanoindentation using spherical and Berkovich indenters

Prevention and treatment of osteoporosis rely on understanding of the micromechanical behaviour of bone and its influence on fracture toughness and cell-mediated adaptation processes. Postyield properties may be assessed by nonlinear finite element simulations of nanoindentation using elastoplastic and damage models. This computational study aims at determining the influence of yield surface shape and damage on the depth-dependent response of bone to nanoindentation using spherical and conical tips. Yield surface shape and damage were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic-to-total work ratio is well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not statistically significant (p<0.0001). For spherical tips, damage was not a significant parameter (p<0.0001). The gained knowledge can be used for developing an inverse method for identification of postelastic properties of bone from nanoindentation.

Embedding of human vertebral bodies leads to higher ultimate load and altered damage localisation under axial compression

Computer tomography (CT)-based finite element (FE) models of vertebral bodies assess fracture load in vitro better than dual energy X-ray absorptiometry, but boundary conditions affect stress distribution under the endplates that may influence ultimate load and damage localisation under post-yield strains. Therefore, HRpQCT-based homogenised FE models of 12 vertebral bodies were subjected to axial compression with two distinct boundary conditions: embedding in polymethylmethalcrylate (PMMA) and bonding to a healthy intervertebral disc (IVD) with distinct hyperelastic properties for nucleus and annulus. Bone volume fraction and fabric assessed from HRpQCT data were used to determine the elastic, plastic and damage behaviour of bone. Ultimate forces obtained with PMMA were 22% higher than with IVD but correlated highly (R2 = 0.99). At ultimate force, distinct fractions of damage were computed in the endplates (PMMA: 6%, IVD: 70%), cortex and trabecular sub-regions, which confirms previous observations that in contrast to PMMA embedding, failure initiated underneath the nuclei in healthy IVDs. In conclusion, axial loading of vertebral bodies via PMMA embedding versus healthy IVD overestimates ultimate load and leads to distinct damage localisation and failure pattern.

Neuromechanical gait adaptations in women with joint hypermobility - an exploratory study

BACKGROUND Joint hypermobility is known to be associated with joint and muscle pain, joint instability and osteoarthritis. Previous work suggested that those individuals present an altered neuromuscular behavior during activities such as level walking. Therefore, the aim of this study was to explore the differences in ground reaction forces, temporal parameters and muscle activation patterns during gait between normomobile and hypermobile women, including symptomatic and asymptomatic hypermobile individuals. METHODS A total of 195 women were included in this cross-sectional study, including 67 normomobile (mean 24.8 [SD 5.4] years) and 128 hypermobile (mean 25.8 [SD 5.4] years), of which 56 were further classified as symptomatic and 47 as asymptomatic. The remaining 25 subjects could not be further classified. Ground reaction forces and muscle activation from six leg muscles were measured while the subjects walked at a self-selected speed on an instrumented walkway. Temporal parameters were derived from ground reaction forces and a foot accelerometer. The normomobile and hypermobile groups were compared using independent samples t-tests, whereas the normomobile, symptomatic and asymptomatic hypermobile groups were compared using one-way ANOVAs with Tukey post-hoc tests (significance level=0.05). FINDINGS Swing phase duration was higher among hypermobile (P=0.005) and symptomatic hypermobile (P=0.018) compared to normomobile women. The vastus medialis (P=0.049) and lateralis (P=0.030) and medial gastrocnemius (P=0.011) muscles showed higher mean activation levels during stance in the hypermobile compared to the normomobile group. INTERPRETATION Hypermobile women might alter their gait pattern in order to stabilize their knee joint.

Patient-specific finite-element simulation of the human cornea: A clinical validation study on cataract surgery

The planning of refractive surgical interventions is a challenging task. Numerical modeling has been proposed as a solution to support surgical intervention and predict the visual acuity, but validation on patient specific intervention is missing. The purpose of this study was to validate the numerical predictions of the post-operative corneal topography induced by the incisions required for cataract surgery. The corneal topography of 13 patients was assessed preoperatively and postoperatively (1-day and 30-day follow-up) with a Pentacam tomography device. The preoperatively acquired geometric corneal topography – anterior, posterior and pachymetry data – was used to build patient-specific finite element models. For each patient, the effects of the cataract incisions were simulated numerically and the resulting corneal surfaces were compared to the clinical postoperative measurements at one day and at 30-days follow up. Results showed that the model was able to reproduce experimental measurements with an error on the surgically induced sphere of 0.38D one day postoperatively and 0.19D 30 days postoperatively. The standard deviation of the surgically induced cylinder was 0.54D at the first postoperative day and 0.38D 30 days postoperatively. The prediction errors in surface elevation and curvature were below the topography measurement device accuracy of ±5μm and ±0.25D after the 30-day follow-up. The results showed that finite element simulations of corneal biomechanics are able to predict post cataract surgery within topography measurement device accuracy. We can conclude that the numerical simulation can become a valuable tool to plan corneal incisions in cataract surgery and other ophthalmosurgical procedures in order to optimize patients' refractive outcome and visual function.

An experimentally validated finite element method for augmented vertebral bodies

Background Finite element models of augmented vertebral bodies require a realistic modelling of the cement infiltrated region. Most methods published so far used idealized cement shapes or oversimplified material models for the augmented region. In this study, an improved, anatomy-specific, homogenized finite element method was developed and validated to predict the apparent as well as the local mechanical behavior of augmented vertebral bodies. Methods Forty-nine human vertebral body sections were prepared by removing the cortical endplates and scanned with high-resolution peripheral quantitative CT before and after injection of a standard and a low-modulus bone cement. Forty-one specimens were tested in compression to measure stiffness, strength and contact pressure distributions between specimens and loading-plates. From the remaining eight, fourteen cylindrical specimens were extracted from the augmented region and tested in compression to obtain material properties. Anatomy-specific finite element models were generated from the CT data. The models featured element-specific, density-fabric-based material properties, damage accumulation, real cement distributions and experimentally determined material properties for the augmented region. Apparent stiffness and strength as well as contact pressure distributions at the loading plates were compared between simulations and experiments. Findings The finite element models were able to predict apparent stiffness (R2 > 0.86) and apparent strength (R2 > 0.92) very well. Also, the numerically obtained pressure distributions were in reasonable quantitative (R2 > 0.48) and qualitative agreement with the experiments. Interpretation The proposed finite element models have proven to be an accurate tool for studying the apparent as well as the local mechanical behavior of augmented vertebral bodies.

Tissue properties of the human vertebral body sub-structures evaluated by means of microindentation

Purpose The better understanding of vertebral mechanical properties can help to improve the diagnosis of vertebral fractures. As the bone mechanical competence depends not only from bone mineral density (BMD) but also from bone quality, the goal of the present study was to investigate the anisotropic indentation moduli of the different sub-structures of the healthy human vertebral body and spondylophytes by means of microindentation. Methods Six human vertebral bodies and five osteophytes (spondylophytes) were collected and prepared for microindentation test. In particular, indentations were performed on bone structural units of the cortical shell (along axial, circumferential and radial directions), of the endplates (along the anterio-posterior and lateral directions), of the trabecular bone (along the axial and transverse directions) and of the spondylophytes (along the axial direction). A total of 3164 indentations down to a maximum depth of 2.5 µm were performed and the indentation modulus was computed for each measurement. Results The cortical shell showed an orthotropic behavior (indentation modulus, Ei, higher if measured along the axial direction, 14.6±2.8 GPa, compared to the circumferential one, 12.3±3.5 GPa, and radial one, 8.3±3.1 GPa). Moreover, the cortical endplates (similar Ei along the antero-posterior, 13.0±2.9 GPa, and along the lateral, 12.0±3.0 GPa, directions) and the trabecular bone (Ei= 13.7±3.4 GPa along the axial direction versus Ei=10.9±3.7 GPa along the transverse one) showed transversal isotropy behavior. Furthermore, the spondylophytes showed the lower mechanical properties measured along the axial direction (Ei=10.5±3.3 GPa). Conclusions The original results presented in this study improve our understanding of vertebral biomechanics and can be helpful to define the material properties of the vertebral substructures in computational models such as FE analysis.

Finite element based estimation of contact areas and pressures of the human scaphoid in various functional positions of the hand

The scaphoid is the most frequently fractured carpal bone. When investigating fixation stability, which may influence healing, knowledge of forces and moments acting on the scaphoid is essential. The aim of this study was to evaluate cartilage contact forces acting on the intact scaphoid in various functional wrist positions using finite element modeling. A novel methodology was utilized as an attempt to overcome some limitations of earlier studies, namely, relatively coarse imaging resolution to assess geometry, assumption of idealized cartilage thicknesses and neglected cartilage pre-stresses in the unloaded joint. Carpal bone positions and articular cartilage geometry were obtained independently by means of high resolution CT imaging and incorporated into finite element (FE) models of the human wrist in eight functional positions. Displacement driven FE analyses were used to resolve inter-penetration of cartilage layers, and provided contact areas, forces and pressure distribution for the scaphoid bone. The results were in the range reported by previous studies. Novel findings of this study were: (i) cartilage thickness was found to be heterogeneous for each bone and vary considerably between carpal bones; (ii) this heterogeneity largely influenced the FE results and (iii) the forces acting on the scaphoid in the unloaded wrist were found to be significant. As major limitations, accuracy of the method was found to be relatively low, and the results could not be compared to independent experiments. The obtained results will be used in a following study to evaluate existing and recently developed screws used to fix scaphoid fractures.

Morphology–elasticity relationships using decreasing fabric information of human trabecular bone from three major anatomical locations

With improving clinical CT scanning technology, the accuracy of CT-based finite element (FE) models of the human skeleton may be ameliorated by an enhanced description of apparent level bone mechanical properties. Micro-finite element (μFE) modeling can be used to study the apparent elastic behavior of human cancellous bone. In this study, samples from the femur, radius and vertebral body were investigated to evaluate the predictive power of morphology–elasticity relationships and to compare them across different anatomical regions. μFE models of 701 trabecular bone cubes with a side length of 5.3 mm were analyzed using kinematic boundary conditions. Based on the FE results, four morphology–elasticity models using bone volume fraction as well as full, limited or no fabric information were calibrated for each anatomical region. The 5 parameter Zysset–Curnier model using full fabric information showed excellent predictive power with coefficients of determination ( r2adj ) of 0.98, 0.95 and 0.94 of the femur, radius and vertebra data, respectively, with mean total norm errors between 14 and 20%. A constant orthotropy model and a constant transverse isotropy model, where the elastic anisotropy is defined by the model parameters, yielded coefficients of determination between 0.90 and 0.98 with total norm errors between 16 and 25%. Neglecting fabric information and using an isotropic model led to r2adj between 0.73 and 0.92 with total norm errors between 38 and 49%. A comparison of the model regressions revealed minor but significant (p<0.01) differences for the fabric–elasticity model parameters calibrated for the different anatomical regions. The proposed models and identified parameters can be used in future studies to compute the apparent elastic properties of human cancellous bone for homogenized FE models.