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.

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.

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.

Lobe specific Ca2+-calmodulin nano-domain in neuronal spines: a single molecule level analysis.

Calmodulin (CaM) is a ubiquitous Ca(2+) buffer and second messenger that affects cellular function as diverse as cardiac excitability, synaptic plasticity, and gene transcription. In CA1 pyramidal neurons, CaM regulates two opposing Ca(2+)-dependent processes that underlie memory formation: long-term potentiation (LTP) and long-term depression (LTD). Induction of LTP and LTD require activation of Ca(2+)-CaM-dependent enzymes: Ca(2+)/CaM-dependent kinase II (CaMKII) and calcineurin, respectively. Yet, it remains unclear as to how Ca(2+) and CaM produce these two opposing effects, LTP and LTD. CaM binds 4 Ca(2+) ions: two in its N-terminal lobe and two in its C-terminal lobe. Experimental studies have shown that the N- and C-terminal lobes of CaM have different binding kinetics toward Ca(2+) and its downstream targets. This may suggest that each lobe of CaM differentially responds to Ca(2+) signal patterns. Here, we use a novel event-driven particle-based Monte Carlo simulation and statistical point pattern analysis to explore the spatial and temporal dynamics of lobe-specific Ca(2+)-CaM interaction at the single molecule level. We show that the N-lobe of CaM, but not the C-lobe, exhibits a nano-scale domain of activation that is highly sensitive to the location of Ca(2+) channels, and to the microscopic injection rate of Ca(2+) ions. We also demonstrate that Ca(2+) saturation takes place via two different pathways depending on the Ca(2+) injection rate, one dominated by the N-terminal lobe, and the other one by the C-terminal lobe. Taken together, these results suggest that the two lobes of CaM function as distinct Ca(2+) sensors that can differentially transduce Ca(2+) influx to downstream targets. We discuss a possible role of the N-terminal lobe-specific Ca(2+)-CaM nano-domain in CaMKII activation required for the induction of synaptic plasticity.

Brain fiber tract plasticity in experimental spinal cord injury: diffusion tensor imaging.

Diffusion tensor imaging (DTI) and immunohistochemistry were performed in spinal cord injured rats to understand the basis for activation of multiple regions in the brain observed in functional magnetic resonance imaging (fMRI) studies. The measured fractional anisotropy (FA), a scalar measure of diffusion anisotropy, along the region encompassing corticospinal tracts (CST) indicates significant differences between control and injured groups in the 3 to 4 mm area posterior to bregma that correspond to internal capsule and cerebral peduncle. Additionally, DTI-based tractography in injured animals showed increased number of fibers that extend towards the cortex terminating in the regions that were activated in fMRI. Both the internal capsule and cerebral peduncle demonstrated an increase in GFAP-immunoreactivity compared to control animals. GAP-43 expression also indicates plasticity in the internal capsule. These studies suggest that the previously observed multiple regions of activation in spinal cord injury are, at least in part, due to the formation of new fibers.

Neurogranin controls the spatiotemporal pattern of postsynaptic Ca2+/CaM signaling.

Neurogranin (Ng) is a postsynaptic IQ-motif containing protein that accelerates Ca(2+) dissociation from calmodulin (CaM), a key regulator of long-term potentiation and long-term depression in CA1 pyramidal neurons. The exact physiological role of Ng, however, remains controversial. Two genetic knockout studies of Ng showed opposite outcomes in terms of the induction of synaptic plasticity. To understand its function, we test the hypothesis that Ng could regulate the spatial range of action of Ca(2+)/CaM based on its ability to accelerate the dissociation of Ca(2+) from CaM. Using a mathematical model constructed on the known biochemistry of Ng, we calculate the cycle time that CaM molecules alternate between the fully Ca(2+) saturated state and the Ca(2+) unbound state. We then use these results and include diffusion of CaM to illustrate the impact that Ng has on modulating the spatial profile of Ca(2+)-saturated CaM within a model spine compartment. Finally, the first-passage time of CaM to transition from the Ca(2+)-free state to the Ca(2+)-saturated state was calculated with or without Ng present. These analyses suggest that Ng regulates the encounter rate between Ca(2+) saturated CaM and its downstream targets during postsynaptic Ca(2+) transients.

Spike timing dependent plasticity: a consequence of more fundamental learning rules.

Spike timing dependent plasticity (STDP) is a phenomenon in which the precise timing of spikes affects the sign and magnitude of changes in synaptic strength. STDP is often interpreted as the comprehensive learning rule for a synapse - the "first law" of synaptic plasticity. This interpretation is made explicit in theoretical models in which the total plasticity produced by complex spike patterns results from a superposition of the effects of all spike pairs. Although such models are appealing for their simplicity, they can fail dramatically. For example, the measured single-spike learning rule between hippocampal CA3 and CA1 pyramidal neurons does not predict the existence of long-term potentiation one of the best-known forms of synaptic plasticity. Layers of complexity have been added to the basic STDP model to repair predictive failures, but they have been outstripped by experimental data. We propose an alternate first law: neural activity triggers changes in key biochemical intermediates, which act as a more direct trigger of plasticity mechanisms. One particularly successful model uses intracellular calcium as the intermediate and can account for many observed properties of bidirectional plasticity. In this formulation, STDP is not itself the basis for explaining other forms of plasticity, but is instead a consequence of changes in the biochemical intermediate, calcium. Eventually a mechanism-based framework for learning rules should include other messengers, discrete change at individual synapses, spread of plasticity among neighboring synapses, and priming of hidden processes that change a synapse's susceptibility to future change. Mechanism-based models provide a rich framework for the computational representation of synaptic plasticity.

Glutamate iontophoresis induces long-term potentiation in the absence of evoked presynaptic activity

$\rm\underline{L}$ong-$\rm\underline{t}$erm $\rm\underline{p}$otentiation (LTP) is a candidate cellular mechanism underlying mammalian learning and memory. Protocols that induce LTP typically involve afferent stimulation. The experiments described in this dissertation tested the hypothesis that LTP induction does not require presynaptic activity. The significance of this hypothesis is underscored by results suggesting that LTP expression may involve activity-dependent presynaptic changes. An induction protocol using glutamate iontophoresis was developed that reliably induces LTP in hippocampal slices without afferent stimulation (ionto-LTP). Ionto-LTP is induced when excitatory postsynaptic potentials are completely blocked with adenosine and $\rm\underline{t}$etrodo$\rm\underline{t}$o$\rm\underline{x}$in (TTX). These results suggest constraints on the involvement of presynaptic mechanisms and putative retrograde messengers in LTP induction and expression; namely, these processes must function without many forms of activity-dependent presynaptic processes.^ In testing the role of pre-and postsynaptic mechanisms in LTP expression whole-cell recordings were used to examine the frequency and amplitude of $\rm\underline{s}$pontaneous $\rm\underline{e}$xcitatory $\rm\underline{p}$o$\rm\underline{s}$ynaptic $\rm\underline{c}$urrents (sEPSCs) in CA1 pyramidal neurons. sEPSCs where comprised of an equal mixture of TTX insensitive miniature EPSCs and sEPSCs that appeared to result from spontaneous action potentials (i.e., TTX sensitive EPSCs). The detection of all sEPSCs was virtually eliminated by CNQX, suggesting that sEPSCs were glutamate mediated synaptic events. Changes in the amplitude and frequency sEPSCs were examined during the expression of ionto-LTP to obtain new information about the cellular location of mechanisms involved in synaptic plasticity. The findings of this dissertation show that ionto-LTP expression results from increased sEPSC amplitude in the absence of lasting increases in sEPSC frequency. Potentiation of sEPSC amplitude without changes in sEPSC frequency has been previously interpreted to be due to postsynaptic mechanisms. Although this interpretation is supported by findings from peripheral synapses, its application to the central nervous system is unclear. Therefore, alternative mechanisms are also considered in this dissertation. Models based on increased release probability for action potential dependent transmitter release appear insufficient to explain our results. The most straightforward interpretation of the results in this dissertation is that LTP induced by glutamate iontophoresis on dendrites of CA1 pyramidal neurons is mediated by postsynaptic mechanisms. ^

Nerve Injury-Induced Neuropathic Pain Causes Disinhibition of the Anterior Cingulate Cortex

Neuropathic pain caused by peripheral nerve injury is a debilitating neurological condition of high clinical relevance. On the cellular level, the elevated pain sensitivity is induced by plasticity of neuronal function along the pain pathway. Changes in cortical areas involved in pain processing contribute to the development of neuropathic pain. Yet, it remains elusive which plasticity mechanisms occur in cortical circuits. We investigated the properties of neural networks in the anterior cingulate cortex (ACC), a brain region mediating affective responses to noxious stimuli. We performed multiple whole-cell recordings from neurons in layer 5 (L5) of the ACC of adult mice after chronic constriction injury of the sciatic nerve of the left hindpaw and observed a striking loss of connections between excitatory and inhibitory neurons in both directions. In contrast, no significant changes in synaptic efficacy in the remaining connected pairs were found. These changes were reflected on the network level by a decrease in the mEPSC and mIPSC frequency. Additionally, nerve injury resulted in a potentiation of the intrinsic excitability of pyramidal neurons, whereas the cellular properties of interneurons were unchanged. Our set of experimental parameters allowed constructing a neuronal network model of L5 in the ACC, revealing that the modification of inhibitory connectivity had the most profound effect on increased network activity. Thus, our combined experimental and modeling approach suggests that cortical disinhibition is a fundamental pathological modification associated with peripheral nerve damage. These changes at the cortical network level might therefore contribute to the neuropathic pain condition.

Segmentation and Recognition of Multi-Food Meal Images for Carbohydrate Counting

In this paper, we propose novel methodologies for the automatic segmentation and recognition of multi-food images. The proposed methods implement the first modules of a carbohydrate counting and insulin advisory system for type 1 diabetic patients. Initially the plate is segmented using pyramidal mean-shift filtering and a region growing algorithm. Then each of the resulted segments is described by both color and texture features and classified by a support vector machine into one of six different major food classes. Finally, a modified version of the Huang and Dom evaluation index was proposed, addressing the particular needs of the food segmentation problem. The experimental results prove the effectiveness of the proposed method achieving a segmentation accuracy of 88.5% and recognition rate equal to 87%

Norepinephrine drives persistent activity in prefrontal cortex via synergistic α1 and α2 adrenoceptors

Optimal norepinephrine levels in the prefrontal cortex (PFC) increase delay-related firing and enhance working memory, whereas stress-related or pathologically high levels of norepinephrine are believed to inhibit working memory via α1 adrenoceptors. However, it has been shown that activation of Gq-coupled and phospholipase C-linked receptors can induce persistent firing, a cellular correlate of working memory, in cortical pyramidal neurons. Therefore, despite its importance in stress and cognition, the exact role of norepinephrine in modulating PFC activity remains elusive. Using electrophysiology and optogenetics, we report here that norepinephrine induces persistent firing in pyramidal neurons of the PFC independent of recurrent fast synaptic excitation. This persistent excitatory effect involves presynaptic α1 adrenoceptors facilitating glutamate release and subsequent activation of postsynaptic mGluR5 receptors, and is enhanced by postsynaptic α2 adrenoceptors inhibiting HCN channel activity. Activation of α2 adrenoceptors or inhibition of HCN channels also enhances cholinergic persistent responses in pyramidal neurons, providing a mechanism of crosstalk between noradrenergic and cholinergic inputs. The present study describes a novel cellular basis for the noradrenergic control of cortical information processing and supports a synergistic combination of intrinsic and network mechanisms for the expression of mnemonic properties in pyramidal neurons.

Reelin immunoreactivity in dissociated cultures of the postnatal hippocampus

Reelin is an extracellular matrix glycoprotein expressed in different nerve cell populations in the developing, early postnatal and adult central nervous system. During histogenesis of the neocortex and hippocampus, reelin is present in Cajal-Retzius cells and other early neurons and contributes to correct layering of these regions. During early postnatal life, pioneer neurons disappear and reelin expression establishes in a subpopulation of cortical and hippocampal GABAergic interneurons, where it is maintained throughout adult life. We studied the developmental distribution pattern of reelin in dissociated cultures obtained from the early postnatal hippocampus to verify whether or not such a maturation phenomenon is maintained in vitro. Reelin is expressed both in Cajal-Retzius cells and multipolar and pyramidal neurons in younger cultures. The density of reelin-positive Cajal-Retzius cells dropped drastically by about 84% in 4-week-old cultures. Multipolar and pyramidal neurons containing reelin represented 12% of the total cell population in younger cultures and decreased by about 25% after 3 to 4 weeks of cultivation. Their density was significantly lower in cultures of the same age treated with glutamate receptor antagonists. These reelin-positive multipolar and pyramidal neurons were heterogeneous, including a larger amount of non-GABAergic, and 30-40% of GABAergic neurons. Cells double labeled for reelin and the GABA synthesizing enzyme glutamic acid decarboxylase represented about 4% of the total neuron population in culture and their density remained constant with age. It is thus possible that the decrease in the total reelin population may selectively be of importance to the larger non-GABAergic fraction of reelin cells. This study shows that reelin-expressing neurons are maintained in dissociated cultures of the neonatal hippocampus and their distribution and age-dependent changes in density resemble those of the early postnatal hippocampus in vivo.

Modulation of orientation-selective neurons by motion: when additive, when multiplicative?

The recurrent interaction among orientation-selective neurons in the primary visual cortex (V1) is suited to enhance contours in a noisy visual scene. Motion is known to have a strong pop-up effect in perceiving contours, but how motion-sensitive neurons in V1 support contour detection remains vastly elusive. Here we suggest how the various types of motion-sensitive neurons observed in V1 should be wired together in a micro-circuitry to optimally extract contours in the visual scene. Motion-sensitive neurons can be selective about the direction of motion occurring at some spot or respond equally to all directions (pandirectional). We show that, in the light of figure-ground segregation, direction-selective motion neurons should additively modulate the corresponding orientation-selective neurons with preferred orientation orthogonal to the motion direction. In turn, to maximally enhance contours, pandirectional motion neurons should multiplicatively modulate all orientation-selective neurons with co-localized receptive fields. This multiplicative modulation amplifies the local V1-circuitry among co-aligned orientation-selective neurons for detecting elongated contours. We suggest that the additive modulation by direction-specific motion neurons is achieved through synaptic projections to the somatic region, and the multiplicative modulation by pandirectional motion neurons through projections to the apical region of orientation-specific pyramidal neurons. For the purpose of contour detection, the V1-intrinsic integration of motion information is advantageous over a downstream integration as it exploits the recurrent V1-circuitry designed for that task.

Experimental validation of a mean field model of mineralized collagen fiber arrays at two levels of hierarchy

In the course of this study, stiffness of a fibril array of mineralized collagen fibrils modeled with a mean field method was validated experimentally at site-matched two levels of tissue hierarchy using mineralized turkey leg tendons (MTLT). The applied modeling approaches allowed to model the properties of this unidirectional tissue from nanoscale (mineralized collagen fibrils) to macroscale (mineralized tendon). At the microlevel, the indentation moduli obtained with a mean field homogenization scheme were compared to the experimental ones obtained with microindentation. At the macrolevel, the macroscopic stiffness predicted with micro finite element (μFE) models was compared to the experimental stiffness measured with uniaxial tensile tests. Elastic properties of the elements in μFE models were injected from the mean field model or two-directional microindentations. Quantitatively, the indentation moduli can be properly predicted with the mean-field models. Local stiffness trends within specific tissue morphologies are very weak, suggesting additional factors responsible for the stiffness variations. At macrolevel, the μFE models underestimate the macroscopic stiffness, as compared to tensile tests, but the correlations are strong.

Surface Roughness and Microhardness of Two Recent CAD/CAM-Materials After Storage

Objectives: To investigate surface roughness and microhardness of two recent resin-ceramic materials for computer-aided design/computer-aided manufacturing (CAD/CAM) after polishing with three polishing systems. Surface roughness and microhardness were measured immediately after polishing and after six months storage including monthly artificial toothbrushing. Methods: Sixty specimens of Lava Ultimate (3M ESPE) and 60 specimens of VITA ENAMIC (VITA Zahnfabrik) were roughened in a standardized manner and polished with one of three polishing systems (n=20/group): Sof-Lex XT discs (SOFLEX; three-step (medium-superfine); 3M ESPE), VITA Polishing Set Clinical (VITA; two-step; VITA Zahnfabrik), or KENDA Unicus (KENDA; one-step; KENDA Dental). Surface roughness (Ra; μm) was measured with a profilometer and microhardness (Vickers; VHN) with a surface hardness indentation device. Ra and VHN were measured immediately after polishing and after six months storage (tap water, 37°C) including monthly artificial toothbrushing (500 cycles/month, toothpaste RDA ~70). Ra- and VHN-values were analysed with nonparametric ANOVA followed by Wilcoxon rank sum tests (α=0.05). Results: For Lava Ultimate, Ra (mean [standard deviation] before/after storage) remained the same when polished with SOFLEX (0.18 [0.09]/0.19 [0.10]; p=0.18), increased significantly with VITA (1.10 [0.44]/1.27 [0.39]; p=0.0001), and decreased significantly with KENDA (0.35 [0.07]/0.33 [0.08]; p=0.03). VHN (mean [standard deviation] before/after storage) decreased significantly regardless of polishing system (SOFLEX: 134.1 [5.6]/116.4 [3.6], VITA: 138.2 [10.5]/115.4 [5.9], KENDA: 135.1 [6.2]/116.7 [6.3]; all p<0.0001). For VITA ENAMIC, Ra (mean [standard deviation] before/after storage) increased significantly when polished with SOFLEX (0.37 [0.18]/0.41 [0.14]; p=0.01) and remained the same with VITA (1.32 [0.37]/1.31 [0.40]; p=0.58) and with KENDA (0.81 [0.35]/0.78 [0.32]; p=0.21). VHN (mean [standard deviation] before/after storage) remained the same regardless of polishing system (SOFLEX: 284.9 [24.6]/282.4 [31.8], VITA: 284.6 [28.5]/276.4 [25.8], KENDA: 292.6 [26.9]/282.9 [24.3]; p=0.42-1.00). Conclusion: Surface roughness and microhardness of Lava Ultimate was more affected by storage and artificial toothbrushing than was VITA ENAMIC.