Experimental, theoretical and numerical investigation of the nonlinear micromechanical properties of bone

Aging societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the bone hierarchical organization. A good understanding has been reached for elastic properties on several length scales, but up to now there is a lack of reliable postyield data on the lower length scales. In order to be able to describe the behavior of bone at the microscale, an anisotropic elastic-viscoplastic damage model was developed using an eccentric generalized Hill criterion and nonlinear isotropic hardening. The model was implemented as a user subroutine in Abaqus and verified using single element tests. A FE simulation of microindentation in lamellar bone was finally performed show-ing that the new constitutive model can capture the main characteristics of the indentation response of bone. As the generalized Hill criterion is limited to elliptical and cylindrical yield surfaces and the correct shape for bone is not known, a new yield surface was developed that takes any convex quadratic shape. The main advantage is that in the case of material identification the shape of the yield surface does not have to be anticipated but a minimization results in the optimal shape among all convex quadrics. The generality of the formulation was demonstrated by showing its degeneration to classical yield surfaces. Also, existing yield criteria for bone at multiple length scales were converted to the quadric formulation. Then, a computational study to determine the influence of yield surface shape and damage on the in-dentation response of bone using spherical and conical tips was performed. The constitutive model was adapted to the quadric criterion and yield surface shape and critical damage were varied. They 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 were found to be very well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not a significant fac-tor, while for spherical tips damage was insignificant. All inverse methods based on microindentation suffer from a lack of uniqueness of the found material properties in the case of nonlinear material behavior. Therefore, monotonic and cyclic micropillar com-pression tests in a scanning electron microscope allowing a straightforward interpretation comple-mented by microindentation and macroscopic uniaxial compression tests were performed on dry ovine bone to identify modulus, yield stress, plastic deformation, damage accumulation and failure mecha-nisms. While the elastic properties were highly consistent, the postyield deformation and failure mech-anisms differed between the two length scales. A majority of the micropillars showed a ductile behavior with strain hardening until failure by localization in a slip plane, while the macroscopic samples failed in a quasi-brittle fashion with microcracks coalescing into macroscopic failure surfaces. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behavior of bone at the microscale to a quasi-brittle response driven by the growth of preexisting cracks along interfaces or in the vicinity of pores at the macroscale. Subsequently, a study was undertaken to quantify the topological variability of indentations in bone and examine its relationship with mechanical properties. Indentations were performed in dry human and ovine bone in axial and transverse directions and their topography measured by AFM. Statistical shape modeling of the residual imprint allowed to define a mean shape and describe the variability with 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was highly consistent and free of any pile up. A few of the topological parameters, in particular depth, showed significant correlations to variations in mechanical properties, but the cor-relations were not very strong or consistent. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small de-viations from the ideal case. As the uniaxial properties measured by micropillar compression are in conflict with the current literature on bone indentation, another dissipative mechanism has to be present. The elastic-viscoplastic damage model was therefore extended to viscoelasticity. The viscoelastic properties were identified from macroscopic experiments, while the quasistatic postelastic properties were extracted from micropillar data. It was found that viscoelasticity governed by macroscale properties has very little influence on the indentation curve and results in a clear underestimation of the creep deformation. Adding viscoplasticity leads to increased creep, but hardness is still highly overestimated. It was possible to obtain a reasonable fit with experimental indentation curves for both Berkovich and spherical indenta-tion when abandoning the assumption of shear strength being governed by an isotropy condition. These results remain to be verified by independent tests probing the micromechanical strength prop-erties in tension and shear. In conclusion, in this thesis several tools were developed to describe the complex behavior of bone on the microscale and experiments were performed to identify its material properties. Micropillar com-pression highlighted a size effect in bone due to the presence of preexisting cracks and pores or inter-faces like cement lines. It was possible to get a reasonable fit between experimental indentation curves using different tips and simulations using the constitutive model and uniaxial properties measured by micropillar compression. Additional experimental work is necessary to identify the exact nature of the size effect and the mechanical role of interfaces in bone. Deciphering the micromechanical behavior of lamellar bone and its evolution with age, disease and treatment and its failure mechanisms on several length scales will help preventing fractures in the elderly in the future.

Block of the hERG channel by bupivacaine: Electrophysiological and modeling insights towards stereochemical optimization

The hERG voltage-gated potassium channel mediates the cardiac I(Kr) current, which is crucial for the duration of the cardiac action potential. Undesired block of the channel by certain drugs may prolong the QT interval and increase the risk of malignant ventricular arrhythmias. Although the molecular determinants of hERG block have been intensively studied, not much is known about its stereoselectivity. Levo-(S)-bupivacaine was the first drug reported to have a higher affinity to block hERG than its enantiomer. This study strives to understand the principles underlying the stereoselectivity of bupivacaine block with the help of mutagenesis analyses and molecular modeling simulations. Electrophysiological measurements of mutated hERG channels allowed for the identification of residues involved in bupivacaine binding and stereoselectivity. Docking and molecular mechanics simulations for both enantiomers of bupivacaine and terfenadine (a non-stereoselective blocker) were performed inside an open-state model of the hERG channel. The predicted binding modes enabled a clear depiction of ligand-protein interactions. Estimated binding affinities for both enantiomers were consistent with electrophysiological measurements. A similar computational procedure was applied to bupivacaine enantiomers towards two mutated hERG channels (Tyr652Ala and Phe656Ala). This study confirmed, at the molecular level, that bupivacaine stereoselectively binds the hERG channel. These results help to lay the foundation for structural guidelines to optimize the cardiotoxic profile of drug candidates in silico.

Three-dimensional surface reconstruction within noncontact diffuse optical tomography using structured light

A main field in biomedical optics research is diffuse optical tomography, where intensity variations of the transmitted light traversing through tissue are detected. Mathematical models and reconstruction algorithms based on finite element methods and Monte Carlo simulations describe the light transport inside the tissue and determine differences in absorption and scattering coefficients. Precise knowledge of the sample's surface shape and orientation is required to provide boundary conditions for these techniques. We propose an integrated method based on structured light three-dimensional (3-D) scanning that provides detailed surface information of the object, which is usable for volume mesh creation and allows the normalization of the intensity dispersion between surface and camera. The experimental setup is complemented by polarization difference imaging to avoid overlaying byproducts caused by inter-reflections and multiple scattering in semitransparent tissue.

Differential roles of Rho-kinase and myosin light chain kinase in regulating shape, adhesion, and migration of HT1080 fibrosarcoma cells

We present evidence for differential roles of Rho-kinase and myosin light chain kinase (MLCK) in regulating shape, adhesion, migration, and chemotaxis of human fibrosarcoma HT1080 cells on laminin-coated surfaces. Pharmacological inhibition of Rho-kinase by Y-27632 or inhibition of MLCK by W-7 or ML-7 resulted in significant attenuation of constitutive myosin light chain phosphorylation. Rho-kinase inhibition resulted in sickle-shaped cells featuring long, thin F-actin-rich protrusions. These cells adhered more strongly to laminin and migrated faster. Inhibition of MLCK in contrast resulted in spherical cells and marked impairment of adhesion and migration. Inhibition of myosin II activation with blebbistatin resulted in a morphology similar to that induced by Y-27632 and enhanced migration and adhesion. Cells treated first with blebbistatin and then with ML-7 also rounded up, suggesting that effects of MLCK inhibition on HT1080 cell shape and motility are independent of inhibition of myosin activity.

Direct combination of nanoparticle fabrication and exposure to lung cell cultures in a closed setup as a method to simulate accidental nanoparticle exposure of humans

The tremendous application potential of nanosized materials stays in sharp contrast to a growing number of critical reports of their potential toxicity. Applications of in vitro methods to assess nanoparticles are severely limited through difficulties in exposing cells of the respiratory tract directly to airborne engineered nanoparticles. We present a completely new approach to expose lung cells to particles generated in situ by flame spray synthesis. Cerium oxide nanoparticles from a single run were produced and simultaneously exposed to the surface of cultured lung cells inside a glovebox. Separately collected samples were used to measure hydrodynamic particle size distribution, shape, and agglomerate morphology. Cell viability was not impaired by the conditions of the glovebox exposure. The tightness of the lung cell monolayer, the mean total lamellar body volume, and the generation of oxidative DNA damage revealed a dose-dependent cellular response to the airborne engineered nanoparticles. The direct combination of production and exposure allows studying particle toxicity in a simple and reproducible way under environmental conditions.

Presence of UV filters in surface water and the effects of phenylbenzimidazole sulfonic acid on rainbow trout (Oncorhynchus mykiss) following a chronic toxicity test

UV filters belong to a group of compounds that are used by humans and are present in municipal waste-waters, effluents from sewage treatment plants and surface waters. Current information regarding UV filters and their effects on fish is limited. In this study, the occurrence of three commonly used UV filters - 2-phenylbenzimidazole-5-sulfonic acid (PBSA), 2-hydroxy-4-methoxybenzophenone (benzophenone-3, BP-3) and 5-benzoyl-4-hydroxy-2-methoxy-benzenesulfonic acid (benzophenone-4, BP-4) - in South Bohemia (Czech Republic) surface waters is presented. PBSA concentrations (up to 13μgL(-1)) were significantly greater than BP-3 or BP-4 concentrations (up to 620 and 390ngL(-1), respectively). On the basis of these results, PBSA was selected for use in a toxicity test utilizing the common model organism rainbow trout (Oncorhynchus mykiss). Fish were exposed to three concentrations of PBSA (1, 10 and 1000µgL(-1)) for 21 and 42 days. The PBSA concentrations in the fish plasma, liver and kidneys were elevated after 21 and 42 days of exposure. PBSA increased activity of certain P450 cytochromes. Exposure to PBSA also changed various biochemical parameters and enzyme activities in the fish plasma. However, no pathological changes were obvious in the liver or gonads.

The influence of bedrock orientation on the landscape evolution, surface morphology and denudation (10Be) at the Niesen, Switzerland

Landscape evolution and surface morphology in mountainous settings are a function of the relative importance between sediment transport processes acting on hillslopes and in channels, modulated by climate variables. The Niesen nappe in the Swiss Penninic Prealps presents a unique setting in which opposite facing flanks host basins underlain by identical lithologies, but contrasting litho-tectonic architectures where lithologies either dip parallel to the topographic slope or in the opposite direction (i.e. dip slope and non-dip slope). The north-western facing Diemtigen flank represents such a dip slope situation and is characterized by a gentle topography, low hillslope gradients, poorly dissected channels, and it hosts large landslides. In contrast, the south-eastern facing Frutigen side can be described as non-dip slope flank with deeply incised bedrock channels, high mean hillslope gradients and high relief topography. Results from morphometric analysis reveal that noticeable differences in morphometric parameters can be related to the contrasts in the relative importance of the internal hillslope-channel system between both valley flanks. While the contrasting dip-orientations of the underlying flysch bedrock has promoted hillslope and channelized processes to contrasting extents and particularly the occurrence of large landslides on the dip slope flank, the flank averaged beryllium-10 (10Be)-derived denudation rates are very similar and range between 0.20 and 0.26 mm yr−1. In addition, our denudation rates offer no direct relationship to basin's slope, area, steepness or concavity index, but reveal a positive correlation to mean basin elevation that we interpret as having been controlled by climatically driven factors such as frost-induced processes and orographic precipitation. Our findings illustrate that while the landscape properties in this part of the northern Alpine border can mainly be related to the tectonic architecture of the underlying bedrock, the denudation rates have a strong orographic control through elevation dependent mean annual temperature and precipitation.

Modeling of ore-forming and geoenvironmental systems: Roles of fluid flow and chemical reaction processes

Ore-forming and geoenviromental systems commonly involve coupled fluid flowand chemical reaction processes. The advanced numerical methods and computational modeling have become indispensable tools for simulating such processes in recent years. This enables many hitherto unsolvable geoscience problems to be addressed using numerical methods and computational modeling approaches. For example, computational modeling has been successfully used to solve ore-forming and mine site contamination/remediation problems, in which fluid flow and geochemical processes play important roles in the controlling dynamic mechanisms. The main purpose of this paper is to present a generalized overview of: (1) the various classes and models associated with fluid flow/chemically reacting systems in order to highlight possible opportunities and developments for the future; (2) some more general issues that need attention in the development of computational models and codes for simulating ore-forming and geoenviromental systems; (3) the related progresses achieved on the geochemical modeling over the past 50 years or so; (4) the general methodology for modeling of oreforming and geoenvironmental systems; and (5) the future development directions associated with modeling of ore-forming and geoenviromental systems.

A Biofilm Pocket Model to Evaluate Different Non-Surgical Periodontal Treatment Modalities in Terms of Biofilm Removal and Reformation, Surface Alterations and Attachment of Periodontal Ligament Fibroblasts.

BACKGROUND AND AIM There is a lack of suitable in vitro models to evaluate various treatment modalities intending to remove subgingival bacterial biofilm. Consequently, the aims of this in vitro-study were: a) to establish a pocket model enabling mechanical removal of biofilm and b) to evaluate repeated non-surgical periodontal treatment with respect to biofilm removal and reformation, surface alterations, tooth hard-substance-loss, and attachment of periodontal ligament (PDL) fibroblasts. MATERIAL AND METHODS Standardized human dentin specimens were colonized by multi-species biofilms for 3.5 days and subsequently placed into artificially created pockets. Non-surgical periodontal treatment was performed as follows: a) hand-instrumentation with curettes (CUR), b) ultrasonication (US), c) subgingival air-polishing using erythritol (EAP) and d) subgingival air-polishing using erythritol combined with chlorhexidine digluconate (EAP-CHX). The reduction and recolonization of bacterial counts, surface roughness (Ra and Rz), the caused tooth substance-loss (thickness) as well as the attachment of PDL fibroblasts were evaluated and statistically analyzed by means of ANOVA with Post-Hoc LSD. RESULTS After 5 treatments, bacterial reduction in biofilms was highest when applying EAP-CHX (4 log10). The lowest reduction was found after CUR (2 log10). Additionally, substance-loss was the highest when using CUR (128±40 µm) in comparison with US (14±12 µm), EAP (6±7 µm) and EAP-CHX (11±10) µm). Surface was roughened when using CUR and US. Surfaces exposed to US and to EAP attracted the highest numbers of PDL fibroblasts. CONCLUSION The established biofilm model simulating a periodontal pocket combined with interchangeable placements of test specimens with multi-species biofilms enables the evaluation of different non-surgical treatment modalities on biofilm removal and surface alterations. Compared to hand instrumentation the application of ultrasonication and of air-polishing with erythritol prevents from substance-loss and results in a smooth surface with nearly no residual biofilm that promotes the reattachment of PDL fibroblasts.

Muscle velocity recovery cycles: Comparison between surface and needle recordings.

INTRODUCTION Recording of muscle velocity recovery cycles (MVRCs) has been developed as a technique to investigate the pathophysiology of muscle diseases. MVRCs have been measured by direct muscle stimulation and concentric electromyographic needle recording. This study was undertaken to determine whether recordings can be made with surface electrodes. METHODS MVRCs with 1 and 2 conditioning stimuli were recorded simultaneously with concentric needle and surface electrodes from the brachioradialis muscle in 12 healthy volunteers. Muscle relative refractory period, early and late supernormality, and extra-late supernormality were compared between the recording techniques. RESULTS Surface recordings were possible in all subjects. The multifiber action potentials recorded with surface electrodes were smaller than those recorded with needles, but there was no significant difference between any of their MVRC properties. CONCLUSIONS MVRCs can be recorded with surface electrodes in healthy subjects. The use of surface electrodes may facilitate the technique of recording MVRCs. Muscle Nerve 53: 205-208, 2016.

How to be in a good shape? The influence of clone morphology on cell competition

Cell competition is a conserved mechanism where slow proliferating cells (so called losers) are eliminated by faster proliferating neighbors (so called winners) through apoptosis.(1) It is an important process which prevents developmental malformations and maintains tissue fitness in aging adults.(2) Recently, we have shown that the probability of elimination of loser cells correlates with the surface of contact between losers and winners in Myc-induced competition.(3) Moreover, we have characterized an active mechanism that increases the surface of contact between losers and winners, hence accelerating the elimination of loser cells. This is the first indication that cell shape and mechanics can influence cell competition. Here, we will discuss the consequence of the relationship between shape and competition, as well as the relevance of this model for other modes of competition.