Exocytotic dynamics in human chromaffin cells: experiments and modeling

Chromaffin cells have been widely used to study neurosecretion since they exhibit similar calcium dependence of several exocytotic steps as synaptic terminals do, but having the enormous advantage of being neither as small or fast as neurons, nor as slow as endocrine cells. In the present study, secretion associated to experimental measurements of the exocytotic dynamics in human chromaffin cells of the adrenal gland was simulated by using a model that combines stochastic and deterministic approaches for short and longer depolarizing pulses, respectively. Experimental data were recorded from human chromaffin cells, obtained from healthy organ donors, using the perforated patch configuration of the patch-clamp technique. We have found that in human chromaffin cells, secretion would be mainly managed by small pools of non-equally fusion competent vesicles, slowly refilled over time. Fast secretion evoked by brief pulses can be predicted only when 75% of one of these pools (the “ready releasable pool” of vesicles, abbreviated as RRP) are co-localized to Ca^2?+? channels, indicating an immediately releasable pool in the range reported for isolated cells of bovine and rat (Álvarez and Marengo, J Neurochem 116:155–163, 2011). The need for spatial correlation and close proximity of vesicles to Ca^2?+? channels suggests that in human chromaffin cells there is a tight control of those releasable vesicles available for fast secretion.     

Serial sectioning and electron microscopy of large tissue volumes for 3D analysis and reconstruction: a case study of the calyx of Held

Serial section electron microscopy is typically applied to investigation of small tissue volumes encompassing subcellular structures. However, in neurobiology, the need to relate subcellular structure to organization of neural circuits can require investigation of large tissue volumes at ultrastructural resolution. Analysis of ultrastructure and three-dimensional reconstruction of even one to a few cells is time consuming, and still does not generate the necessary numbers of observations to form well-grounded insights into biological principles. We describe an assemblage of existing computer-based methods and strategies for graphical analysis of large photographic montages to accomplish the study of multiple neurons through large tissue volumes. Sample preparation, data collection and subsequent analyses can be completed within 3-4 months. These methods generate extremely large data sets that can be mined in future studies of nervous system organization.     

Fluorescence microscopy of calcium and synaptic vesicle dynamics during synapse formation in tissue culture

The signal transduction process involved in the development of the nerve terminal is an intriguing question in developmental neurobiology. During the formation of the neuromuscular junction, presynaptic development is induced by growth cone's contact with the target muscle cell. Fluorescence microscopy with specific markers has made it possible to follow signalling events during this process. By using fluorescent calcium indicators, such as fura-2 and fluo-3, we found that a rise in intracellular calcium is elicited in the growth cone upon its contact with a target, and this calcium signal can also be elicited by local application of basic fibroblast growth factor. To monitor the clustering of synaptic vesicles in response to target contact, the fluorescent vesicular probe FMl-43 was used. With this probe, we observed that packets of synaptic vesicle are already present along the length of naive neurite, which has not encountered its synaptic target. The activity-dependent loading of FMl-43 indicates that these packets can undergo exocytosis and endocytosis upon depolarization. Time-lapse recording showed that these packets are quite mobile. Upon target contact, synaptic vesicles become clustered and immobilized at the contact site. The methodology and instrumentation used in these studies are described in this article. 1998 © Chapman & Hall     

Lipid a from lipopolysaccharide recognition: Structure,dynamics and cooperativity by molecular dynamics simulations

Molecular dynamics simulations of Lipid A and its natural precursor Lipid IVA from E.coli have been carried out free in solution, bound to the myeliod differentiation protein 2 (MD2) and in the complex of MD2 with the toll like receptor 4 (TLR4). In addition, simulations of the ligand free MD2 and MD2‐TLR4 complex were performed. A structural and energetic characterization of the bound and unbound states of Lipid A/IVA was generated. As the crystal structures depict, the main driving force for MD2‐Lipid A/IVA are the hydrophobic interactions between the aliphatic tails and the MD2 cavity. The charged phosphate groups do strongly interact with positively charged residues, located at the surface of MD2. However, they are not essential for keeping the lipids in the cavity, indicating a more prominent role in binding recognition and ionic interactions with TLR4 at the MD2/TLR4 interface. Interestingly, in the absence of any ligand MD2 rapidly closes, blocking the binding cavity. The presence of TLR4, though changing the dynamics, was not able to impede the aforementioned closing event. We hypothesize that fluctuations of the H1 region are essential for this phenomenon, and it is plausible that an equilibrium between the open and closed states exists, although the lengths of our simulations are not sufficient to encompass the reversible process. The MD2/Lipid A‐TLR4 complex simulations show that the presence of the ligand energetically stabilizes the complex relative to the ligand‐free structures, indicating cooperativity in the binding process. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Released fraction and total size of a pool of immediately available transmitter quanta at a calyx synapse.

The size of a pool of readily releasable vesicles at a giant brainstem synapse, the calyx of Held, was probed with three independent approaches. Using simultaneous pre- and postsynaptic whole-cell recordings, two forms of presynaptic Ca2+ stimuli were applied in rapid succession: uncaging of Ca2+ by flash photolysis and the opening of voltage-gated Ca2+ channels. The ensuing transmitter release showed a nearly complete cross-inhibition between the two stimuli, indicating the depletion of a limited pool of about 700 transmitter quanta. The pool size was confirmed in experiments using enhanced extracellular Ca2+ concentrations, as well as short, high-frequency stimulus trains. The results reveal a surprisingly large pool of functionally available vesicles, of which a fraction of about 0.2 is released by a single presynaptic action potential under physiological conditions.     

Differential segregation in a cell-cell contact interface: the dynamics of the immunological synapse

Receptor-ligand couples in the cell-cell contact interface between a T cell and an antigen-presenting cell form distinct geometric patterns and undergo spatial rearrangement within the contact interface. Spatial segregation of the antigen and adhesion receptors occurs within seconds of contact, central aggregation of the antigen receptor then occurring over 1-5 min. This structure, called the immunological synapse, is becoming a paradigm for localized signaling. However, the mechanisms driving its formation, in particular spatial segregation, are currently not understood. With a reaction diffusion model incorporating thermodynamics, elasticity, and reaction kinetics, we examine the hypothesis that differing bond lengths (extracellular domain size) is the driving force behind molecular segregation. We derive two key conditions necessary for segregation: a thermodynamic criterion on the effective bond elasticity and a requirement for the seeding/nucleation of domains. Domains have a minimum length scale and will only spontaneously coalesce/aggregate if the contact area is small or the membrane relaxation distance large. Otherwise, differential attachment of receptors to the cytoskeleton is required for central aggregation. Our analysis indicates that differential bond lengths have a significant effect on synapse dynamics, i.e., there is a significant contribution to the free energy of the interaction, suggesting that segregation by differential bond length is important in cell-cell contact interfaces and the immunological synapse.     

The dynamics of protons,aluminium, and calcium in the rhizosphere of maize cultivated in tropical acid soils: experimental study and modelling

The goals of this work were to understand the dynamics of H⁺, Al and Ca in the rhizosphere of maize cultivated in tropical acid soils, and to evaluate the contribution of the dissolution kinetics of the Al-hydroxides to Al dynamics. The study of the dissolution kinetics was based on a comparison between experimental and simulated data, using a model of the chemical processes in the rhizosphere. Two Oxisols, pH 5.1 and 4.6, and one Ultisol, pH 5.2, were studied. An Al-tolerant maize variety (Zea mays L.) was grown for 14 days on a 3-mm thick soil layer. The composition of the soil and the soil solution, together with the concentration of Al in the roots, were determined throughout the experiment. The results showed that root growth (i) decreased the soil solution pH, up to one pH unit, (ii) increased Al concentration in the soil solution, (iii) increased exchangeable Al, and (iv) decreased exchangeable Ca. Soil solution pH, exchangeable Al, and exchangeable Ca were closely linked. Exchangeable Al increased 1.5 – 3.0 times, due to the dissolution of easily mobilised Al components. In addition, Al accumulation in roots depended mainly on Al in the soil solution. Modelling the interactions between H⁺, Al, and Ca proved that the main factor determining Al in the soil solution was the kinetic reactivity of the easily mobilised Al components. These components, probably poorly crystallised Al-hydroxides, are key players in the functioning of the rhizosphere in tropical acid soils.     

Physiological modulation of voltage-dependent inactivation in the cardiac muscle L-type calcium channel: a modelling study

The inactivation of the L-type Ca²⁺ current is composed of voltage-dependent and calcium-dependent mechanisms. The relative contribution of these processes is still under dispute and the idea that the voltage-dependent inactivation could be subject to further modulation by other physiological processes had been ignored. This study sought to model physiological modulation of inactivation of the current in cardiac ventricular myocytes, based upon the recent detailed experimental data that separated total and voltage-dependent inactivation (VDI) by replacing extracellular Ca²⁺ with Mg²⁺ and monitoring L-type Ca²⁺ channel behaviour by outward K⁺ current flowing through the channel in the absence of inward current flow. Calcium-dependent inactivation (CDI) was based upon Ca²⁺ influx and formulated from data that was recorded during β-adrenergic stimulation of the myocytes. Ca²⁺ influx and its competition with non-selective monovalent cation permeation were also incorporated into channel permeation in the model. The constructed model could closely reproduce the experimental Ba²⁺ and Ca²⁺ current results under basal condition where no β-stimulation was added after a slight reduction of the development of fast voltage-dependent inactivation with depolarization. The model also predicted that under β-adrenergic stimulation voltage-dependent inactivation is lost and calcium-dependent inactivation largely compensates it. The developed model thus will be useful to estimate the respective roles of VDI and CDI of L-type Ca²⁺ channels in various physiological and pathological conditions of the heart which would otherwise be difficult to show experimentally.     

Calcium-dependent isoforms of protein kinase C mediate posttetanic potentiation at the calyx of Held

High-frequency stimulation leads to a transient increase in the amplitude of evoked synaptic transmission that is known as posttetanic potentiation (PTP). Here we examine the roles of the calcium-dependent protein kinase C isoforms PKCα and PKCβ in PTP at the calyx of Held synapse. In PKCα/β double knockouts, 80% of PTP is eliminated, whereas basal synaptic properties are unaffected. PKCα and PKCβ produce PTP by increasing the size of the readily releasable pool of vesicles evoked by high-frequency stimulation and by increasing the fraction of this pool released by the first stimulus. PKCα and PKCβ do not facilitate presynaptic calcium currents. The small PTP remaining in double knockouts is mediated partly by an increase in miniature excitatory postsynaptic current amplitude and partly by a mechanism involving myosin light chain kinase. These experiments establish that PKCα and PKCβ are crucial for PTP and suggest that long-lasting presynaptic calcium increases produced by tetanic stimulation may activate these isoforms to produce PTP.     

High mortality and enhanced recovery: modelling the countervailing effects of disturbance on population dynamics

Disturbances cause high mortality in populations while simultaneously enhancing population growth by improving habitats. These countervailing effects make it difficult to predict population dynamics following disturbance events. To address this challenge, we derived a novel form of the logistic growth equation that permits time‐varying carrying capacity and growth rate. We combined this equation with concepts drawn from disturbance ecology to create a general model for population dynamics in disturbance‐prone systems. A river flooding example using three insect species (a fast life‐cycle mayfly, a slow life‐cycle dragonfly and an ostracod) found optimal tradeoffs between disturbance frequency vs. magnitude and a close fit to empirical data in 62% of cases. A savanna fire analysis identified fire frequencies of 3–4 years that maximised population size of a perennial grass. The model shows promise for predicting population dynamics after multiple disturbance events and for management of river flows and fire regimes.     

Calcium dynamics at a plastic synapse in APLYSIA

Because certain primitive behavioral responses in the large sea snail Aplysia have recently been linked to neurophysiological events at a synaptic level, special interest attaches to the role played by calcium ions at such synapses. Using an extended version of the model applied earlier to trace the flow of energy and information through a ganglion of the medicinal leech (Triffet & Green, 1980), the authors investigate the electropotential effects of small transient localized changes in the calcium concentration near the inner membrane surface of a neuron in the resting state.When this state is well below the firing threshold, changes in Ca²⁺ concentration less than 10⁻⁸ M are shown to result only in low-level harmonic background oscillations. When the potential of the neurons is closer to threshold, however, and/or the Ca²⁺ concentration is of the order of 10⁻⁸ M, easily recognizable graded potentials appear, and these grow into firing peaks when the calcium concentration is increased still further.Though no attempt is made to deal with the amplification effects dependent on calcium-vesicle interactions and the related release of transmitter molecules, a unified mechanism for the underlying calcium ion dynamics is proposed. Graded potentials of increasing size are associated with a progressive localized thickening of the inner and outer Debye layers. Moreover, the transverse and longitudinal calcium currents set up in such regions prove adequate to account for both the depletion of Ca²⁺ ions necessary to achieve habituation, and the increase in their concentration required for sensitization.     

Wave propagation in cardiac tissue and effects of intracellular calcium dynamics (computer simulation study)

Computer simulation using Luo-Rudy I1 model of ventricular myocyte showed that intracellular calcium dynamics become irregular in case of high rate stimulation. This causes the transition from stationary to nonstationary spiral wave and its breakup in 2D model of cardiac tissue. Obtained results suggest how ventricular fibrillation may occur due to the abnormalities of intracellular calcium dynamics. The short review of existing cardiac cell models with calcium dynamics is presented.     

Live-cell dynamics and the role of costimulation in immunological synapse formation

Using transfected fibroblasts expressing both wild-type I-E(k) and green fluorescent protein-tagged I-E(k) with covalently attached antigenic peptide, we have monitored movement of specific MHC:peptide complexes during CD4(+) T cell-APC interactions by live-cell video microscopy. Ag recognition occurs within 30 s of T cell-APC contact, as shown by a sharp increase in cytoplasmic calcium ion concentration. Within 1 min, small MHC:peptide clusters form in the contact zone that coalesce into an immunological synapse over 3-20 min. When T cells conjugated to APC move across the APC surface, they appear to drag the synapse with them. This system was used to examine the role of costimulation in the formation of the immunological synapse. Blocking CD80/CD28 or ICAM-1/LFA-1 interactions alters synapse morphology and reduces the area and density of accumulated complexes. These reductions correlate with reduced T cell proliferation, while CD69 and CD25 expression and TCR down-modulation remain unaffected. Thus, costimulation is essential for normal mature immunological synapse formation.     

Density-dependent population dynamics in clonal organisms: a modelling approach

1. Density dependence may act at several stages in an organisms life-cycle (e.g. on mortality, fecundity, etc.), but not all density-dependent processes necessarily regulate population size. In this paper I use a density manipulation experiment to determine the effects of density on the transition rates between different size classes of the clonal zoanthid Palythoa caesia Dana 1846. I then formulate a density-dependent matrix model of population dynamics of Palythoa , and perform a series of sensitivity analyses on the model to determine at what stage in the life-cycle regulation acts.
2. Seven of the 16 transition probabilities decreased with density, most of them being shrinkage (due to loss of tissue or fission) and stasis (the self–self transition) of medium and large colonies. The only probability to increase was for the stasis of large colonies. Recruitment was quadratically dependent on density, peaking at intermediate densities.
3. Equilibrium cover in the model was 84% and was reached in ≈40 years. To determine which density-dependent transitions were involved in population regulation, the strength of density dependence was varied in each independently. This sensitivity analysis showed that only changes in the probabilities of large colonies remaining large and producing medium colonies, were regulating.
4. These results suggest that regulation is primarily acting on fission of large colonies to produce intermediate-sized colonies, in combination with size specific growth rates. Fission rates decrease greatly with density, resulting in a greater proportion of large colonies at high densities and large colonies grow more slowly than small. Overall, this behaviour is very similar to that of clonal plants which have a phalanx type life history.     

Mathematical modelling of human P2X-mediated plasma membrane electrophysiology and calcium dynamics in microglia

Regulation of cytosolic calcium (Ca²⁺) dynamics is fundamental to microglial function. Temporal and spatial Ca²⁺ fluxes are induced from a complicated signal transduction pathway linked to brain ionic homeostasis. In this paper, we develop a novel biophysical model of Ca²⁺ and sodium (Na⁺) dynamics in human microglia and evaluate the contribution of purinergic receptors (P2XRs) to both intracellular Ca²⁺ and Na⁺ levels in response to agonist/ATP binding. This is the first comprehensive model that integrates P2XRs to predict intricate Ca²⁺ and Na⁺ transient responses in microglia. Specifically, a novel compact biophysical model is proposed for the capture of whole-cell patch-clamp currents associated with P2X₄ and P2X₇ receptors, which is composed of only four state variables. The entire model shows that intricate intracellular ion dynamics arise from the coupled interaction between P2X₄ and P2X₇ receptors, the Na⁺/Ca²⁺ exchanger (NCX), Ca²⁺ extrusion by the plasma membrane Ca²⁺ ATPase (PMCA), and Ca²⁺ and Na⁺ leak channels. Both P2XRs are modelled as two separate adenosine triphosphate (ATP) gated Ca²⁺ and Na⁺ conductance channels, where the stoichiometry is the removal of one Ca²⁺ for the hydrolysis of one ATP molecule. Two unique sets of model parameters were determined using an evolutionary algorithm to optimise fitting to experimental data for each of the receptors. This allows the proposed model to capture both human P2X₇ and P2X₄ data (hP2X₇ and hP2X₄). The model architecture enables a high degree of simplicity, accuracy and predictability of Ca²⁺ and Na⁺ dynamics thus providing quantitative insights into different behaviours of intracellular Na⁺ and Ca²⁺ which will guide future experimental research. Understanding the interactions between these receptors and other membrane-bound transporters provides a step forward in resolving the qualitative link between purinergic receptors and microglial physiology and their contribution to brain pathology.     

Impact of sarcoplasmic reticulum calcium release on calcium dynamics and action potential morphology in human atrial myocytes: a computational study

Electrophysiological studies of the human heart face the fundamental challenge that experimental data can be acquired only from patients with underlying heart disease. Regarding human atria, there exist sizable gaps in the understanding of the functional role of cellular Ca2+ dynamics, which differ crucially from that of ventricular cells, in the modulation of excitation-contraction coupling. Accordingly, the objective of this study was to develop a mathematical model of the human atrial myocyte that, in addition to the sarcolemmal (SL) ion currents, accounts for the heterogeneity of intracellular Ca2+ dynamics emerging from a structurally detailed sarcoplasmic reticulum (SR). Based on the simulation results, our model convincingly reproduces the principal characteristics of Ca2+ dynamics: 1) the biphasic increment during the upstroke of the Ca2+ transient resulting from the delay between the peripheral and central SR Ca2+ release, and 2) the relative contribution of SL Ca2+ current and SR Ca2+ release to the Ca2+ transient. In line with experimental findings, the model also replicates the strong impact of intracellular Ca2+ dynamics on the shape of the action potential. The simulation results suggest that the peripheral SR Ca2+ release sites define the interface between Ca2+ and AP, whereas the central release sites are important for the fire-diffuse-fire propagation of Ca2+ diffusion. Furthermore, our analysis predicts that the modulation of the action potential duration due to increasing heart rate is largely mediated by changes in the intracellular Na+ concentration. Finally, the results indicate that the SR Ca2+ release is a strong modulator of AP duration and, consequently, myocyte refractoriness/excitability. We conclude that the developed model is robust and reproduces many fundamental aspects of the tight coupling between SL ion currents and intracellular Ca2+ signaling. Thus, the model provides a useful framework for future studies of excitation-contraction coupling in human atrial myocytes.     

The role of calcium in prolonged modification of a GABAergic synapse.

Caudal hair cell impulses cause postsynaptic inhibition of ipsilateral type B photoreceptors in the snail Hermissenda. This inhibition is shown to be GABAergic according to a number of criteria. HPLC, mass spectrophotometric, and immunocytochemical techniques demonstrated the presence of GABA in the hair cells and their axons. GABA agonists and antagonists mimic and block the synaptic effect in a manner consistent with endogenous GABAergic transmission. Other properties, including I-V relations, conductance changes and reversal potentials, are comparable for exogenous GABA responses and endogenous effects of the hair cell impulses. This inhibitory synapse has been found to undergo a long-lasting transformation into an excitatory synapse if GABA release is paired with post-synaptic depolarization. GABA, via GABAA and GABAB receptors in the B cell, causes the opening of calcium sensitive chloride and potassium channels that leads to the post-synaptic hyperpolarization. GABA also induces a long-lasting intracellular calcium elevation at the terminal branches of the B cell that greatly outlasts the voltage responses. Synaptic transformation induced by pairings is caused by a decrease in both GABA induced chloride and potassium conductances in the post-synaptic B cell, as well as a significant prolongation of the intracellular calcium accumulation in the B cell's terminal axonal branches.