Finite element analysis for prediction of bone strength

Article preview View full access options BoneKEy Reports | Review Print Email Share/bookmark Finite element analysis for prediction of bone strength Philippe K Zysset, Enrico Dall'Ara, Peter Varga & Dieter H Pahr Affiliations Corresponding author BoneKEy Reports (2013) 2, Article number: 386 (2013) doi:10.1038/bonekey.2013.120 Received 03 January 2013 Accepted 25 June 2013 Published online 07 August 2013 Article tools Citation Reprints Rights & permissions Abstract Abstract• References• Author information Finite element (FE) analysis has been applied for the past 40 years to simulate the mechanical behavior of bone. Although several validation studies have been performed on specific anatomical sites and load cases, this study aims to review the predictability of human bone strength at the three major osteoporotic fracture sites quantified in recently completed in vitro studies at our former institute. Specifically, the performance of FE analysis based on clinical computer tomography (QCT) is compared with the ones of the current densitometric standards, bone mineral content, bone mineral density (BMD) and areal BMD (aBMD). Clinical fractures were produced in monotonic axial compression of the distal radii, vertebral sections and in side loading of the proximal femora. QCT-based FE models of the three bones were developed to simulate as closely as possible the boundary conditions of each experiment. For all sites, the FE methodology exhibited the lowest errors and the highest correlations in predicting the experimental bone strength. Likely due to the improved CT image resolution, the quality of the FE prediction in the peripheral skeleton using high-resolution peripheral CT was superior to that in the axial skeleton with whole-body QCT. Because of its projective and scalar nature, the performance of aBMD in predicting bone strength depended on loading mode and was significantly inferior to FE in axial compression of radial or vertebral sections but not significantly inferior to FE in side loading of the femur. Considering the cumulated evidence from the published validation studies, it is concluded that FE models provide the most reliable surrogates of bone strength at any of the three fracture sites.

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.

DXA predictions of human femoral mechanical properties depend on the load configuration

The aim of this study was to evaluate the ability of dual energy X-rays absorptiometry (DXA) areal bone mineral density (aBMD) measured in different regions of the proximal part of the human femur for predicting the mechanical properties of matched proximal femora tested in two different loading configurations. 36 pairs of fresh frozen femora were DXA scanned and tested until failure in two loading configurations: a fall on the side or a one-legged standing. The ability of the DXA output from four different regions of the proximal femur in predicting the femoral mechanical properties was measured and compared for the two loading scenarios. The femoral neck DXA BMD was best correlated to the femoral ultimate force for both configurations and predicted significantly better femoral failure load (R2=0.80 vs. R2=0.66, P<0.05) when simulating a side than when simulating a standing configuration. Conversely, the work to failure was predicted similarly for both loading configurations (R2=0.54 vs. R2=0.53, P>0.05). Therefore, neck BMD should be considered as one of the key factors for discriminating femoral fracture risk in vivo. Moreover, the better predictive ability of neck BMD for femoral strength if tested in a fall compared to a one-legged stance configuration suggests that DXA's clinical relevance may not be as high for spontaneous femoral fractures than for fractures associated to a fall.

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.

Metabolism and action of 9-$\beta$-D-arabinofuranosyl-2-fluoroadenine

9-$\beta$-D-arabinofuranosyl-2-fluoroadenine (F-ara-A) is an analogue of adenosine and 2$\sp\prime$-deoxyadenosine with potent antitumor activity both in vitro and in vivo. The mechanism of action of F-ara-A was evaluated both in whole cells and in experimental systems with purified enzymes. F-ara-A was converted to its 5$\sp\prime$-triphosphate F-ara-ATP in cells and then incorporated into DNA in a self-limiting manner. About 98% of the incorporated F-ara-AMP residues were located at the 3$\sp\prime$-termini of DNA strands, suggesting a chain termination property of this compound. DNA synthesis in CEM cells was inhibited by F-ara-A treatment with an IC$\sb{50}$ value of 1 $\mu$M. Cells were not able to restore the normal level of DNA synthesis even after being cultured in drug-free medium for 40 h. A DNA primer extension assay with M13mp18(+) single-stranded DNA template using purified human DNA polymerases $\alpha$ and further revealed that F-ara-ATP competed with dATP for incorporation into the A sites of the elongating DNA strands. The incorporation of F-ara-AMP into DNA resulted in a termination of DNA synthesis at the incorporated A sites. Pol $\alpha$ and $\delta$ were not able to efficiently extend the DNA primer with F-ara-AMP at its 3$\sp\prime$-end. Furthermore, the presence of F-ara-AMP at the 3$\sp\prime$-end of an oligodeoxyribonucleotide impaired its ligation with an adjacent DNA fragment by human and T4 ligases. Human DNA polymerase $\alpha$ incorporated more F-ara-AMP into DNA than polymerase $\delta$ and was more sensitive to the inhibition by F-ara-ATP, suggesting that polymerase $\alpha$ may be a preferred target for this analogue. On the other hand, DNA-dependent nucleotide turnover experiments and sequencing gel analysis demonstrated that DNA polymerase $\delta$ was able to remove the incorporated F-ara-AMP residue from the 3$\sp\prime$-end of the DNA strand with its 3$\sp\prime$-5$\sp\prime$ exonuclease activity in vitro, subsequently permitting further elongation of the DNA strand.^ The incorporation of F-ara-AMP into DNA was linearly correlated both with the inhibition of DNA synthesis and with the loss of clonogenicity. Termination of DNA synthesis and deletion of genetic material resulted from F-ara-AMP incorporation may be the mechanism responsible for cytotoxicity of F-ara-A. (Abstract shortened with permission of author.) ^

A PHARMACOLOGICAL EVALUATION OF THE MECHANISM OF CYTOTOXICITY OF 9-BETA-D- XYLOFURANOSYLADENINE IN CHINESE HAMSTER OVARY CELLS

The biochemical determinants of cytotoxicity of the purine nucleoside analog, 9-(beta)-D-xylofuranosyladenine (xyl-A) were studied in wild-type Chinese hamster ovary cells and in nucleoside kinase deficient mutants. It was found that {('3)H}xyl-A was readily phosphorylated to the triphosphate level in both the wild-type and deoxycytidine kinase deficient mutant, but not by the adenosine kinase deficient cells. Values for the apparent Km and Vmax of this uptake process were 43.9 (mu)M and 118.7 nmol/min/10('9) cells, respectively. Cloning procedures indicated that the viability of CHO cells was decreased 90 per cent by a 5-hr incubation with 10 (mu)M xyl-A. However, the toxicity of xyl-A was increased 100-fold by the addition of a nontoxic concentration (10 (mu)M) of the adenosine deaminase inhibitor erythro-9-(2-hydroxy-3-nonyl)adenine (EHNA) to the medium. High-pressure liquid chromatographic analysis indicated that after 5 hr, the concentration of 9-(beta)-D-xylofuranosyladenine 5'-triphosphate (xyl-ATP) in cells incubated with xyl-A plus EHNA was 2.0 mM, four times greater than in those cells incubated with xyl-A alone. Incubation with xyl-A plus EHNA had no significant effect on the cellular concentrations of 5-phosphoribosyl-1-pyrophosphate after 1 hr whereas, treatment with 3'-dexoyadenosine (cordycepin) decreased the concentration of this metabolite. Determinations of the cellular nucleoside triphosphates indicated that under conditions that resulted in an intracellular accumulation of 500 (mu)M xyl-ATP, the endogenous concentrations of neither the ribonucleoside triphosphates nor deoxyribonucleoside triphosphates were significantly different from those of control cells. The ID(,50) for {('3)H}thymidine incorporation into DNA, 105 (mu)M xyl-ATP, was four-fold less than the ID(,50) for {('3)H}uridine incorporation into RNA suggesting that the process of DNA synthesis is more sensitive to the presence of xyl-ATP. When removed from exogenous xyl-A, CHO cells failed to recover their ability to synthesize RNA and DNA, although the intracellular xyl-ATP concentration decreased to less than 35 (mu)M. The selective inhibition of RNA synthesis by 6-azauridine did not prevent the expression of toxicity by xyl-ATP. However, the selective inhibition of DNA synthesis by ara-C significantly spared toxicity in cells that had accumulated an otherwise lethal concentration of xyl-ATP. It is shown that in cells which had accumulated 1.27 mM {('3)H}xyl-ATP, {('3)H}xyl-A was found to terminate cellular RNA chains at a frequency of 1.42 (mu)mol of {('3)H}xyl-A 3' termini per mol of mononucleotide. These results indicate that a general mechanism for the toxicity of xyl-A to CHO cells includes the cellular accumulation of xyl-ATP, which serves as a substrate for RNA synthesizing enzymes and subsequently is incorporated into nascent RNA transcripts as a chain terminator. A specific mechanism involving the premature termination of RNA primers required for the initiation of DNA synthesis is proposed to account for the inhibitory action of xyl-ATP on DNA synthesis. ^

BIOCHEMICAL DETERMINANTS IN THE CYTOTOXIC ACTION OF 9-BETA- D- ARABINOFURANOSYLADENINE

Studies with 9-(beta)-D-arabinfuranosyladenine (ara-A) have illustrated that there are numerous cell functions which are affected by the nucleoside analog and its metabolites. Although the precise mechanism responsible for the cytotoxicity of this drug is not known, it is presently thought tha

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.

Local and general anaesthesia do not influence outcome of transfemoral aortic valve implantation.

BACKGROUND There is great variability for the type of anaesthesia used during TAVI, with no clear consensus coming from comparative studies or guidelines. We sought to detect regional differences in the anaesthetic management of patients undergoing transcatheter aortic valve implantation (TAVI) in Europe and to evaluate the relationship between type of anaesthesia and in-hospital and 1year outcome. METHODS Between January 2011 and May 2012 the Sentinel European TAVI Pilot Registry enrolled 2807 patients treated via a transfemoral approach using either local (LA-group, 1095 patients, 39%) or general anaesthesia (GA-group, 1712 patients, 61%). RESULTS A wide variation in LA use was evident amongst the 10 participating countries. The use of LA has increased over time (from a mean of 37.5% of procedures in the first year, to 57% in last 6months, p<0.01). MI, major stroke as well as in-hospital death rate (7.0% LA vs 5.3% GA, p=0.053) had a similar incidence between groups, confirmed in multivariate regression analysis after adjusting for confounders. Dividing our population in tertiles according to the Log-EuroSCORE we found similar mortality under LA, whilst mortality was higher in the highest risk tertile under GA. Survival at 1year, compared by Kaplan-Meier analysis, was similar between groups (log-rank: p=0.1505). CONCLUSIONS Selection of anaesthesia appears to be more influenced by national practice and operator preference than patient characteristics. In the absence of an observed difference in outcomes for either approach, there is no compelling argument to suggest that operators and centres should change their anaesthetic practice.

Percutaneous mitral valve edge-to-edge repair: in-hospital results and 1-year follow-up of 628 patients of the 2011-2012 Pilot European Sentinel Registry.

BACKGROUND The use of transcatheter mitral valve repair (TMVR) has gained widespread acceptance in Europe, but data on immediate success, safety, and long-term echocardiographic follow-up in real-world patients are still limited. OBJECTIVES The aim of this multinational registry is to present a real-world overview of TMVR use in Europe. METHODS The Transcatheter Valve Treatment Sentinel Pilot Registry is a prospective, independent, consecutive collection of individual patient data. RESULTS A total of 628 patients (mean age 74.2 ± 9.7 years, 63.1% men) underwent TMVR between January 2011 and December 2012 in 25 centers in 8 European countries. The prevalent pathogenesis was functional mitral regurgitation (FMR) (n = 452 [72.0%]). The majority of patients (85.5%) were highly symptomatic (New York Heart Association functional class III or higher), with a high logistic EuroSCORE (European System for Cardiac Operative Risk Evaluation) (20.4 ± 16.7%). Acute procedural success was high (95.4%) and similar in FMR and degenerative mitral regurgitation (p = 0.662). One clip was implanted in 61.4% of patients. In-hospital mortality was low (2.9%), without significant differences between groups. The estimated 1-year mortality was 15.3%, which was similar for FMR and degenerative mitral regurgitation. The estimated 1-year rate of rehospitalization because of heart failure was 22.8%, significantly higher in the FMR group (25.8% vs. 12.0%, p[log-rank] = 0.009). Paired echocardiographic data from the 1-year follow-up, available for 368 consecutive patients in 15 centers, showed a persistent reduction in the degree of mitral regurgitation at 1 year (6.0% of patients with severe mitral regurgitation). CONCLUSIONS This independent, contemporary registry shows that TMVR is associated with high immediate success, low complication rates, and sustained 1-year reduction of the severity of mitral regurgitation and improvement of clinical symptoms.

Rabbi Tanchum Jeruschalmi arabischer Commentar zum Buche Josua

hrsg. vom Th. Haarbrücker

The Initial Slope of the Variogram, Foundation of the Trabecular Bone Score, Is Not or Poorly Associated With Vertebral Strength

Trabecular bone score (TBS) rests on the textural analysis of DXA to reflect the decay in trabecular structure characterising osteoporosis. Yet, its discriminative power in fracture studies remains incomprehensible as prior biomechanical tests found no correlation with vertebral strength. To verify this result possibly due to an unrealistic set-up and to cover a wide range of loading scenarios, the data from three previous biomechanical studies using different experimental settings was used. They involved the compressive failure of 62 human lumbar vertebrae loaded 1) via intervertebral discs to mimic the in vivo situation (“full vertebra”), 2) via the classical endplate embedding (“vertebral body”) or 3) via a ball joint to induce anterior wedge failure (“vertebral section”). HR-pQCT scans acquired prior testing were used to simulate anterior-posterior DXA from which areal bone mineral density (aBMD) and the initial slope of the variogram (ISV), the early definition of TBS, were evaluated. Finally, the relation of aBMD and ISV with failure load (Fexp) and apparent failure stress (σexp) was assessed and their relative contribution to a multi-linear model was quantified via ANOVA. We found that, unlike aBMD, ISV did not significantly correlate with Fexp and σexp, except for the “vertebral body” case (r2 = 0.396, p = 0.028). Aside from the “vertebra section” set-up where it explained only 6.4% of σexp (p = 0.037), it brought no significant improvement to aBMD. These results indicate that ISV, a replica of TBS, is a poor surrogate for vertebral strength no matter the testing set-up, which supports the prior observations and raises a fortiori the question of the deterministic factors underlying the statistical relationship between TBS and vertebral fracture risk.

INCREASED EXPRESSION OF ONE OF TWO ADENOSINE DEAMINASE ALLELES IN A HUMAN CHORIOCARCINOMA CELL LINE FOLLOWING SELECTION WITH ADENINE NUCLEOSIDES

The human choriocarcinoma cell line JEG-3 is heterozygous at the adenosine deaminase (ADA) gene locus. Both allelic genes are under strong but incomplete repression causing a very low level expression of the gene locus. Because cytotoxic adenosine analogues such as 9-(beta)-D arabinofuranosyladenine (ara-A) and 9-(beta)-D xylofuranosyladenine (xyl-A) can be specifically detoxified by the action of ADA, these analogues were used to select for JEG-3 derived cells which had increased ADA expression. When JEG-3 cells were subjected to a multi-step, successively increasing dosage of either ara-A or xyl-A, resistant cells with increased ADA expression were generated. This increased ADA expression in the resistant cells was unstable, so that when the selective pressure was removed, cellular ADA expression would decrease. Subclone analysis of xyl-A resistant cells revealed that compared to parental JEG-3 cells, individual resistant cells had either elevated ADA levels or decreased adenosine kinase (ADK) levels or both. This altered ADA and ADK expression in the resistant cells were found to be independent events. Because of high endogenous tissue conversion factor (TCF) expression in the JEG-3 cells, the allelic nature of the increased ADA expression in most of the resistant cells could not be determined. However, several resistant subcloned cells were found to have lost TCF expression. These TCF('-) cells expressed only the ADA*2 allelic gene product. Cell fusion experiments demonstrated that the ADA*1 allelic gene was intact and functional in the A3-1A7 cell line. Chromosomal analysis of the A3-1A7 cells showed that they had no double-minutes or homogeneously staining chromosomal regions, although a pair of new chromosomes were found in these cells. Segregation analysis of the hybrid cells indicated that an ADA*2 allelic gene was probably located on this new chromosome. The analysis of the A3-1A7 cell line suggested that the expression of only ADA 2 in these cells was the result of possibly a cis-deregulation of the ADA gene locus or more probably an amplification of the ADA*2 allelic gene. Two effective positive selection systems for ADA('+) cells were also developed and tested. These selection systems should eventually lead to the isolation of the ADA gene.^

Arabinosylguanine: Metabolism, action, and lineage-specific cytotoxicity in human leukemia cells

9-β-D-arabinosylguanine (ara-G), an analogue of deoxyguanosine, has demonstrated T-lymphoblast selective anti-leukemia activity both in vitro and in vivo in cell lines and primary cells and in phase I investigations. The present work was initiated to identify factors that result in this selectivity. ^ The cytotoxicity of ara-G is manifest only after its phosphorylation. Experiments using cell lines transfected to overexpress specific nucleoside kinases demonstrated that the phosphorylation of ara-G to its monophosphate is by both cytoplasmic deoxycytidine kinase and mitochondria) deoxyguanosine kinase. Ara-G monophosphate is converted to its 5′-triphosphate (ara-GTP) in cells by these kinases and then incorporated into DNA. Mechanistic studies demonstrated that incorporation of ara-GTP into DNA was a necessary event for the induction of cell death. ^ Pharmacokinetic and pharmacodynamic studies utilizing three human acute leukemia cell lines, CEM (T-lymphoblastic), Raji (B-lymphoblastic), and ML-1 (myeloid) were performed. CEM cells were most sensitive to ara-G-induced inhibition of colony formation, accumulated ara-GTP at a faster rate and to a greater degree than either Raji or ML-1, but incorporated the lowest number of ara-G molecules into DNA. The position of incorporation was internal and similar in all cell lines. The terminal elimination phase of ara-GTP was >24 h and similar in these cells. Comparisons between inhibition of colony formation and ara-GTP incorporation into DNA demonstrated that while within a cell line there was correlation among these parameters, between cell lines there was no relationship between number of incorporated ara-G molecules and ara-G(TP)-mediated toxicity suggesting that there were additional factors. ^ The expression of membrane bound Fas and Fast was unchanged in all cell lines. In contrast, there was a 2-fold increase in soluble Fast, which was found exclusively in CEM cells. Ara-G-mediated apoptosis in CEM occurred from all phases of the cell cycle and was abrogated partially by Fas antagonist antibodies. These data suggest that Fas-mediated cell death due to the liberation of sFasL may be responsible for the hypersensitivity to ara-G manifested by immature T-cells such as CEM. The role of Fas in ara-G induced death of acute T-lymphoblastic leukemia cells during therapy needs to be tested. ^