Modeling study of exocytosis in neuroendocrine cells: influence of the geometrical parameters

Close agreement between deterministic versus stochastic modeling of first-passage time to vesicle fusion

Ca2+-dependent cell processes, such as neurotransmitter or endocrine vesicle fusion, are inherently stochastic due to large fluctuations in Ca2+ channel gating, Ca2+ diffusion, and Ca2+ binding to buffers and target sensors. However, previous studies revealed closer-than-expected agreement between deterministic and stochastic simulations of Ca2+ diffusion, buffering, and sensing if Ca2+ channel gating is not Ca2+ dependent. To understand this result more fully, we present a comparative study complementing previous work, focusing on Ca2+ dynamics downstream of Ca2+ channel gating. Specifically, we compare deterministic (mean-field/mass-action) and stochastic simulations of vesicle exocytosis latency, quantified by the probability density of the first-passage time (FPT) to the Ca2+-bound state of a vesicle fusion sensor, following a brief Ca2+ current pulse. We show that under physiological constraints, the discrepancy between FPT densities obtained using the two approaches remains small even if as few as ∼50 Ca2+ ions enter per single channel-vesicle release unit. Using a reduced two-compartment model for ease of analysis, we illustrate how this close agreement arises from the smallness of correlations between fluctuations of the reactant molecule numbers, despite the large magnitude of fluctuation amplitudes. This holds if all relevant reactions are heteroreaction between molecules of different species, as is the case for bimolecular Ca2+ binding to buffers and downstream sensor targets. In this case, diffusion and buffering effectively decorrelate the state of the Ca2+ sensor from local Ca2+ fluctuations. Thus, fluctuations in the Ca2+ sensor's state underlying the FPT distribution are only weakly affected by the fluctuations in the local Ca2+ concentration around its average, deterministically computable value.

α3β4 Acetylcholine Nicotinic Receptors Are Components of the Secretory Machinery Clusters in Chromaffin Cells

The heteromeric assembly of α3 and β4 subunits of acetylcholine nicotinic receptors (nAChRs) seems to mediate the secretory response in bovine chromaffin cells. However, there is no information about the localization of these nAChRs in relationship with the secretory active zones in this cellular model. The present work presents the first evidence that, in fact, a population of these receptors is associated through the F-actin cytoskeleton with exocytotic machinery components, as detected by SNAP-25 labeling. Furthermore, we also prove that, upon stimulation, the probability to find α3β4 nAChRs very close to exocytotic events increases with randomized distributions, thus substantiating the clear dynamic behavior of these receptors during the secretory process. Modeling on secretory dynamics and secretory component distributions supports the idea that α3β4 nAChR cluster mobility could help with improving the efficiency of the secretory response of chromaffin cells. Our study is limited by the use of conventional confocal microscopy; in this sense, a strengthening to our conclusions could come from the use of super-resolution microscopy techniques in the near future.

Influence of T-Bar on Calcium Concentration Impacting Release Probability

The relation of form and function, namely the impact of the synaptic anatomy on calcium dynamics in the presynaptic bouton, is a major challenge of present (computational) neuroscience at a cellular level. The Drosophila larval neuromuscular junction (NMJ) is a simple model system, which allows studying basic effects in a rather simple way. This synapse harbors several special structures. In particular, in opposite to standard vertebrate synapses, the presynaptic boutons are rather large, and they have several presynaptic zones. In these zones, different types of anatomical structures are present. Some of the zones bear a so-called T-bar, a particular anatomical structure. The geometric form of the T-bar resembles the shape of the letter “T” or a table with one leg. When an action potential arises, calcium influx is triggered. The probability of vesicle docking and neurotransmitter release is superlinearly proportional to the concentration of calcium close to the vesicular release site. It is tempting to assume that the T-bar causes some sort of calcium accumulation and hence triggers a higher release probability and thus enhances neurotransmitter exocytosis. In order to study this influence in a quantitative manner, we constructed a typical T-bar geometry and compared the calcium concentration close to the active zones (AZs). We compared the case of synapses with and without T-bars. Indeed, we found a substantial influence of the T-bar structure on the presynaptic calcium concentrations close to the AZs, indicating that this anatomical structure increases vesicle release probability. Therefore, our study reveals how the T-bar zone implies a strong relation between form and function. Our study answers the question of experimental studies (namely “Wichmann and Sigrist, Journal of neurogenetics 2010”) concerning the sense of the anatomical structure of the T-bar.

Modeling the influence of co-localized intracellular calcium stores on the secretory response of bovine chromaffin cells

Catecholamines secretion from chromaffin cells is mediated by a Ca²⁺-dependent process in the submembrane space where the exocytotic machinery is located and high-Ca²⁺ microdomains (HCMDs) are formed by the coordinated activity of a functional triad composed of Ca²⁺ channels, endoplasmic reticulum (ER) and mitochondria. It has been observed experimentally that subpopulations of cortical mitochondria and ER associate to secretory sites in bovine chromaffin cells. Here, we study the effect of the geometrical distribution of the co-localized cortical organelles both in the formation of HCMDs in the vicinity of Ca²⁺ channels and on the secretory activity of bovine chromaffin cells in response to a single voltage pulse. Our simulations indicate that co-localized organelles have a dual role in the formation of HCMDs, having, on the one hand, an amplification effect due to the Ca²⁺-induced Ca²⁺-release mechanism from the ER and, on the other, acting as physical barriers to Ca²⁺ diffusion. In addition, our simulations suggest that the increased levels of Ca²⁺ in the microdomain enhances the secretion of the vesicles co-localized to the Ca²⁺ channels. As a whole, our results support the idea that the functional triads formed by Ca²⁺ channels, subplasmalemma ER and mitochondria have a positive effect on the secretion of catecholamines in bovine chromaffin cells.

The Differential Organization of F-Actin Alters the Distribution of Organelles in Cultured When Compared to Native Chromaffin Cells

Cultured bovine chromaffin cells have been used extensively as a neuroendocrine model to study regulated secretion. In order to extend such experimental findings to the physiological situation, it is necessary to study mayor cellular structures affecting secretion in cultured cells with their counterparts present in the adrenomedullary tissue. F-actin concentrates in a peripheral ring in cultured cells, as witnessed by phalloidin–rodhamine labeling, while extends throughout the cytoplasm in native cells. This result is also confirmed when studying the localization of α-fodrin, a F-actin-associated protein. Furthermore, as a consequence of this redistribution of F-actin, we observed that chromaffin granules and mitochondria located into two different cortical and internal populations in cultured cells, whereas they are homogeneously distributed throughout the cytoplasm in the adrenomedullary tissue. Nevertheless, secretion from isolated cells and adrenal gland pieces is remarkably similar when measured by amperometry. Finally, we generate mathematical models to consider how the distribution of organelles affects the secretory kinetics of intact and cultured cells. Our results imply that we have to consider F-actin structural changes to interpret functional data obtained in cultured neuroendocrine cells.

A theoretical study of factors influencing calcium-secretion coupling in a presynaptic active zone model

A theoretical analysis of some of the relevant factors influencing the calcium time course and the strength and timing of release probabilities of vesicles evoked by an action potential in a calyx-type active zone is presented in this paper. In particular, our study focus on the comparison of cooperative vs non-cooperative calcium binding by the release site and the effect of the number of Ca(2+) binding sites on the calcium sensitivity for release. Regarding the comparison of cooperative and non-cooperative kinetic schemes, our simulations show that quite different results are obtained when considering one or another: a reduction in the release probability of more than a 50% is obtained when considering the cooperative kinetic scheme. Also, a delay in the average time for release appears when using this model for the calcium sensor. Our study also shows that a non-cooperative kinetic binding scheme gives rise to a well defined average calcium level for release assuming that the same kinetic constants are considered for all the sites. Our results also suggest that the central value of the calcium sensitivity for release depends on the number of binding sites N and the dissociation constant KD with a scaling law depending on NKD.

The immediately releasable pool of mouse chromaffin cell vesicles is coupled to P/Q-type calcium channels via the synaptic protein interaction site

It is generally accepted that the immediately releasable pool is a group of readily releasable vesicles that are closely associated with voltage dependent Ca²⁺ channels. We have previously shown that exocytosis of this pool is specifically coupled to P/Q Ca²⁺ current. Accordingly, in the present work we found that the Ca²⁺ current flowing through P/Q-type Ca²⁺ channels is 8 times more effective at inducing exocytosis in response to short stimuli than the current carried by L-type channels. To investigate the mechanism that underlies the coupling between the immediately releasable pool and P/Q-type channels we transiently expressed in mouse chromaffin cells peptides corresponding to the synaptic protein interaction site of Cav2.2 to competitively uncouple P/Q-type channels from the secretory vesicle release complex. This treatment reduced the efficiency of Ca²⁺ current to induce exocytosis to similar values as direct inhibition of P/Q-type channels via ω-agatoxin-IVA. In addition, the same treatment markedly reduced immediately releasable pool exocytosis, but did not affect the exocytosis provoked by sustained electric or high K⁺ stimulation. Together, our results indicate that the synaptic protein interaction site is a crucial factor for the establishment of the functional coupling between immediately releasable pool vesicles and P/Q-type Ca²⁺ channels.

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²⁺ 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²⁺ channels suggests that in human chromaffin cells there is a tight control of those releasable vesicles available for fast secretion.

New insights into the role of the cortical cytoskeleton in exocytosis from neuroendocrine cells

The cortical cytoskeleton is a dense network of filamentous actin (F-actin) that participates in the events associated with secretion from neuroendocrine cells. This filamentous web traps secretory vesicles, acting as a retention system that blocks the access of vesicles to secretory sites during the resting state, and it mediates their active directional transport during stimulation. The changes in the cortical cytoskeleton that drive this functional transformation have been well documented, particularly in cultured chromaffin cells. At the biochemical level, alterations in F-actin are governed by the activity of molecular motors like myosins II and V and by other calcium-dependent proteins that influence the polymerization and cross-linking of F-actin structures. In addition to modulating vesicle transport, the F-actin cortical network and its associated motor proteins also influence the late phases of the secretory process, including membrane fusion and the release of active substances through the exocytotic fusion pore. Here, we discuss the potential interactions between the F-actin cortical web and proteins such as SNAREs during secretion. We also discuss the role of the cytoskeleton in organizing the molecular elements required to sustain regulated exocytosis, forming a molecular structure that foments the efficient release of neurotransmitters and hormones.

Combined computational and experimental approaches to understanding the Ca(2+) regulatory network in neurons

Ca(2+) is a ubiquitous signaling ion that regulates a variety of neuronal functions by binding to and altering the state of effector proteins. Spatial relationships and temporal dynamics of Ca(2+) elevations determine many cellular responses of neurons to chemical and electrical stimulation. There is a wealth of information regarding the properties and distribution of Ca(2+) channels, pumps, exchangers, and buffers that participate in Ca(2+) regulation. At the same time, new imaging techniques permit characterization of evoked Ca(2+) signals with increasing spatial and temporal resolution. However, understanding the mechanistic link between functional properties of Ca(2+) handling proteins and the stimulus-evoked Ca(2+) signals they orchestrate requires consideration of the way Ca(2+) handling mechanisms operate together as a system in native cells. A wide array of biophysical modeling approaches is available for studying this problem and can be used in a variety of ways. Models can be useful to explain the behavior of complex systems, to evaluate the role of individual Ca(2+) handling mechanisms, to extract valuable parameters, and to generate predictions that can be validated experimentally. In this review, we discuss recent advances in understanding the underlying mechanisms of Ca(2+) signaling in neurons via mathematical modeling. We emphasize the value of developing realistic models based on experimentally validated descriptions of Ca(2+) transport and buffering that can be tested and refined through new experiments to develop increasingly accurate biophysical descriptions of Ca(2+) signaling in neurons.

Modeling F-actin cortex influence on the secretory properties of neuroendocrine cells

Chromaffin cells are considered as one of the most valuable models to study regulated exocytosis. In these cells, like in other neuroendocrine systems, an intricate cortical cytoskeleton acts as a retentive network impeding vesicle access to plasma membrane. Therefore, during stimulation this structure suffers a transient reorganization allowing active transport of vesicles toward secretory sites. Interestingly, a combination of confocal microscopy studies and mathematical modeling is showing us new aspects of the influence of cortical cytoskeleton in shaping the secretory properties of excitable cells. In this new vision the F-actin-myosin II cortical cytoskeleton is organized forming polygonal cages with the molecular machinery of exocytosis composed by SNARE proteins and voltage-dependent calcium channels associating with its border. In this way the cytoskeleton not only holds together the essential elements acting during secretion, but we proposed that could also act as a structural factor opposing to the free diffusion of the calcium signal and therefore sustains high levels of the intracellular signal triggering exocytosis.

The F-actin cortical network is a major factor influencing the organization of the secretory machinery in chromaffin cells

We have studied how the F-actin cytoskeleton is involved in establishing the heterogeneous intracellular Ca(2+) levels ([Ca(2+)](i)) and in the organization of the exocytotic machinery in cultured bovine chromaffin cells. Simultaneous confocal visualization of [Ca(2+)](i) and transmitted light studies of the cytoskeleton showed that, following cell stimulation, the maximal signal from the Ca(2+)-sensitive fluorescent dye Fluo-3 was in the empty cytosolic spaces left by cytoskeletal cages. This was mostly due to the accumulation of the dye in spaces devoid of cytoskeletal components, as shown by the use of alternative Ca(2+)-insensitive fluorescent cytosolic markers. In addition to affecting the distribution of such compounds in the cytosol, the cytoskeleton influenced the location of L- and P-Q-type Ca(2+) channel clusters, which were associated with the borders of cytoskeletal cages in resting and stimulated cells. Indeed, syntaxin-1 and synaptotagmin-1, which are components of the secretory machinery, were present in the same location. Furthermore, granule exocytosis took place at these sites, indicating that the organization of the F-actin cytoskeletal cortex shapes the preferential sites for secretion by associating the secretory machinery with preferential sites for Ca(2+) entry. The influence of this cortical organization on the propagation of [Ca(2+)](i) can be modelled, illustrating how it serves to define rapid exocytosis.

Pancreatic islet cells: a model for calcium-dependent peptide release

In mammals the concentration of blood glucose is kept close to 5 mmol∕l. Different cell types in the islet of Langerhans participate in the control of glucose homeostasis. β-cells, the most frequent type in pancreatic islets, are responsible for the synthesis, storage, and release of insulin. Insulin, released with increases in blood glucose promotes glucose uptake into the cells. In response to glucose changes, pancreatic α-, β-, and δ-cells regulate their electrical activity and Ca(2+) signals to release glucagon, insulin, and somatostatin, respectively. While all these signaling steps are stimulated in hypoglycemic conditions in α-cells, the activation of these events require higher glucose concentrations in β and also in δ-cells. The stimulus-secretion coupling process and intracellular Ca(2+) ([Ca(2+)](i)) dynamics that allow β-cells to secrete is well-accepted. Conversely, the mechanisms that regulate α- and δ-cell secretion are still under study. Here, we will consider the glucose-induced signaling mechanisms in each cell type and the mathematical models that explain Ca(2+) dynamics.

The organization of the secretory machinery in chromaffin cells as a major factor in modeling exocytosis

The organization of cytoplasm in excitable cells was a largely ignored factor when mathematical models were developed to understand intracellular calcium and secretory behavior. Here we employed a combination of fluorescent evanescent and transmitted light microscopy to explore the F-actin cytoskeletal organization in the vicinity of secretory sites in cultured bovine chromaffin cells. This technique and confocal fluorescent microscopy show chromaffin granules associated with the borders of cortical cytoskeletal cages forming an intricate tridimensional network. Furthermore, the overexpression of SNAP-25 in these cells also reveals the association of secretory machinery clusters with the borders of these cytoskeletal cages. The importance of these F-actin cage borders is stressed when granules appear to interact and remain associated during exocytosis visualized in acridin orange loaded vesicles. These results will prompt us to propose a model of cytoskeletal cages, where the secretory machinery is associated with its borders. Both the calcium level and the secretory response are enhanced in this geometrical arrangement when compared with a random distribution of the secretory machinery that is not restricted to the borders of the cage.

Modeling of quantal neurotransmitter release kinetics in the presence of fixed and mobile calcium buffers

The local calcium concentration in the active zone of secretion determines the number and kinetics of neurotransmitter quanta released after the arrival of a nerve action potential in chemical synapses. The small size of mammalian neuromuscular junctions does not allow direct measurement of the correlation between calcium influx, the state of endogenous calcium buffers determining the local concentration of calcium and the time course of quanta exocytosis. In this work, we used computer modeling of quanta release kinetics with various levels of calcium influx and in the presence of endogenous calcium buffers with varying mobilities. The results of this modeling revealed the desynchronization of quanta release under low calcium influx in the presence of an endogenous fixed calcium buffer, with a diffusion coefficient much smaller than that of free Ca(2+), and synchronization occurred upon adding a mobile buffer. This corresponds to changes in secretion time course parameters found experimentally (Samigullin et al., Physiol Res 54:129-132, 2005; Bukharaeva et al., J Neurochem 100:939-949, 2007).

Calcium microdomains at presynaptic active zones of vertebrate hair cells unmasked by stochastic deconvolution

Signal transduction by auditory and vestibular hair cells involves an impressive ensemble of finely tuned control mechanisms, strictly dependent on the local intracellular Ca(2+) concentration ([Ca(2+)](i)). The study of Ca(2+) dynamics in hair cells typically combines Ca(2+)-sensitive fluorescent indicators (dyes), patch clamp and optical microscopy to produce images of the patterns of fluorescence of a Ca(2+) indicator following various stimulation protocols. Here we describe a novel method that combines electrophysiological recordings, fluorescence imaging and numerical simulations to effectively deconvolve Ca(2+) signals within cytoplasmic microdomains that would otherwise remain inaccessible to direct observation. The method relies on the comparison of experimental data with virtual signals derived from a Monte Carlo reaction-diffusion model based on a realistic reconstruction of the relevant cell boundaries in three dimensions. The model comprises Ca(2+) entry at individual presynaptic active zones followed by diffusion, buffering, extrusion and release of Ca(2+). Our results indicate that changes of the hair cell [Ca(2+)](i) during synaptic transmission are primarily controlled by the Ca(2+) endogenous buffers both at short (<1mu) and at long (tens of microns) distances from the active zones. We provide quantitative estimates of concentration and kinetics of the hair cell endogenous Ca(2+) buffers and Ca(2+)-ATPases. We finally show that experimental fluorescence data collected during Ca(2+) influx are not interpreted correctly if the [Ca(2+)](i) is estimated by assuming that Ca(2+) equilibrates instantly with its reactants. In our opinion, this approach is of potentially general interest as it can be easily adapted to the study of Ca(2+) dynamics in diverse biological systems.

P/Q Ca2+ channels are functionally coupled to exocytosis of the immediately releasable pool in mouse chromaffin cells

Chromaffin cell exocytosis is triggered by Ca(2+) entry through several voltage-dependent channel subtypes. Because it was postulated that immediately releasable vesicles are closely associated with Ca(2+) channels, we wondered what channel types are specifically coupled to the release of this pool. To study this question, cultured mouse chromaffin cell exocytosis was followed by patch-clamp membrane capacitance measurements. The immediately releasable pool was estimated using paired pulse stimulation, resulting in an upper limit of 31+/-3 fF for control conditions (I(Ca): 25+/-2 pA/pF). The N-type channel blocker omega-conotoxin-GVIA affected neither I(Ca) nor the immediately releasable pool exocytosis; although the L channel blocker nitrendipine decreased current by 50%, it did not reduce this pool significantly; and the R channel inhibitor SNX-482 significantly reduced the current but induced only a moderate decrease in the estimated IRP exocytosis. In contrast, the P/Q channel blocker omega-Agatoxin-IVA decreased I(Ca) by 37% but strongly reduced the immediately releasable pool (upper limit: 6+/-1 fF). We used alpha1A subunit knockout mice to corroborate that P/Q Ca(2+) channels were specifically linked to immediately releasable vesicles, and we found that also in this preparation the exocytosis of this pool was severely decreased (6+/-1 fF). On the other hand, application of a strong stimulus that caused the fusion of most of releasable vesicles (3 min, 50 mM K(+)) induced similar exocytosis for wild type and knockout cells. Finally, whereas application of train stimulation on chromaffin cells derived from wild type mice provoked typical early synchronous and delayed asynchronous exocytosis components, the knockout derived cells presented a strongly depressed early exocytosis but showed a prominent delayed asynchronous component. These results demonstrate that P/Q are the dominant calcium channels associated to the release of immediately releasable pool in mouse chromaffin cells.

Determining calcium concentration in heterogeneous model systems using multiple indicators

Intracellular free calcium concentrations ([Ca2+]i) are assessed by measuring indicator fluorescence in entire cells or subcellular regions using fluorescence microscopy. [Ca2+]i is calculated using equations which link fluorescence intensities (or intensity ratios) to calcium concentrations [G. Grynkiewicz, M. Poenie, R.Y. Tsien, A new generation of Ca2+ indicators with greatly improved fluorescence properties, J. Biol. Chem. 260 (1985) 3440-3450]. However, if calcium ions are heterogeneously distributed within a region of interest, then the observed average fluorescence intensity may not reflect average [Ca2+]i. We assessed potential calcium determination errors in mathematical and experimental models consisting of 'low' and 'high' calcium compartments, using indicators with different affinity for calcium. [Ca2+] calculated using average fluorescence intensity was lower than the actual mean concentrations. Low affinity indicators reported higher (more accurate) values than their high affinity counterparts. To estimate compartment dimensions and respective [Ca2+], we extended the standard approach by using different indicator responses to the same [Ca2+]. While two indicators were sufficient to provide a partial characterization of two-compartment model systems, the use of three or more indicators offered full description of the model provided compartmental [Ca2+] were within the indicator sensitivity ranges. These results show that uneven calcium distribution causes underestimation of actual [Ca2+], and offers novel approaches to estimating calcium heterogeneity.

Computational model of excitable cell indicates ATP free energy dynamics in response to calcium oscillations are undampened by cytosolic ATP buffers

Mitochondria in excitable cells are recurrently exposed to pulsatile calcium gradients that activate cell function. Rapid calcium uptake by the mitochondria has previously been shown to cause uncoupling of oxidative phosphorylation. To test (i) if periodic nerve firing may cause oscillation of the cytosolic thermodynamic potential of ATP hydrolysis and (ii) if cytosolic adenylate (AK) and creatine kinase (CK) ATP buffering reactions dampen such oscillations, a lumped kinetic model of an excitable cell capturing major aspects of the physiology has been developed. Activation of ATP metabolism by low-frequency calcium pulses caused large oscillation of the cytosolic, but not mitochondrial ATP/ADP, ratio. This outcome was independent of net ATP synthesis or hydrolysis during mitochondrial calcium uptake. The AK/CK ATP buffering reactions dampened the amplitude and rate of cytosolic ATP/ADP changes on a timescale of seconds, but not milliseconds. These model predictions suggest that alternative sources of capacitance in neurons and striated muscles should be considered to protect ATP-free energy-driven cell functions.

Software for simulating calcium-triggered exocytotic processes

We describe a software package for the simulation of exocytotic events from readily releasable pools of secretory vesicles in neuroendocrine cells and presynaptic terminals. The visual package Ca3D_Exolab simulates the entry of Ca(2+) through the calcium channels, the kinetic reactions of calcium with buffers, the diffusion of calcium and mobile buffers, and the kinetic reactions of calcium with the secretory vesicles. The location of both channels and secretory vesicles can be set by using a graphical interface. Calcium and buffer concentrations at different depths from the cellular membrane and capacitance time courses are obtained as outputs. The software package also provides a descriptive statistical data analysis of the different output data.

Calcium3D: a visual software package for the simulation of calcium buffered diffusion in neuroendocrine cells

We present Calcium3D, a user-friendly software package for simulating calcium triggered processes in neuroendocrine cells. We use Monte Carlo methods for the simulation of the basic processes involved: entry of calcium into the cytoplasm, the diffusion of ions and mobile intracellular calcium buffers inside the intracellular medium, and the kinetics of the reaction of calcium with these buffers. The outputs of the simulation are calcium and buffer concentrations as a function of time and for different depths from the cellular membrane.

Calcium gradients and exocytosis in bovine adrenal chromaffin cells

The relationship between the localized Ca(2+) concentration and depolarization-induced exocytosis was studied in patch-clamped adrenal chromaffin cells using pulsed-laser Ca(2+) imaging and membrane capacitance measurements. Short depolarizing voltage steps induced Ca(2+) gradients and small "synchronous" increases in capacitance during the pulses. Longer pulses increased the capacitance changes, which saturated at 16 fF, suggesting the presence of a small immediately releasable pool of fusion-ready vesicles. A Hill plot of the capacitance changes versus the estimated Ca(2+) concentration in a thin (100 nm) shell beneath the membrane gave n = 2.3 and K(d) = 1.4 microM. Repetitive stimulation elicited a more complex pattern of exocytosis: early pulses induced synchronous capacitance increases, but after five or more pulses there was facilitation of the synchronous responses and gradual increases in capacitance continued between pulses (asynchronous exocytosis) as the steep submembrane Ca(2+) gradients collapsed. Raising the pipette Ca(2+) concentration led to early facilitation of the synchronous response and early appearance of asynchronous exocytosis. We used this data to develop a kinetic model of depolarization-induced exocytosis, where Ca(2+)-dependent fusion of vesicles occurs from a small immediately releasable pool with an affinity of 1-2 microM and vesicles are mobilized to this pool in a Ca(2+)-dependent manner.

Consequences of molecular-level Ca2+ channel and synaptic vesicle colocalization for the Ca2+ microdomain and neurotransmitter exocytosis: a monte carlo study

Morphological and biochemical studies indicate association between voltage-gated Ca2+ channels and the vesicle docking complex at vertebrate presynaptic active zones, which constrain the separation between some Ca2+ channels and vesicles to 20 nm or less. To address the effect of the precise geometrical relationship among the vesicles, the Ca2+ channel, and the proteins of the release machinery on neurotransmitter release, we developed a Monte Carlo simulation of Ca2+ diffusion and buffering with nanometer resolution. We find that the presence of a vesicle as a diffusion barrier alters the shape of the Ca2+ microdomain of a single Ca2+ channel around the vesicle. This effect is maximal in the vicinity of the vesicle and depends critically on the vesicle's distance from the plasmalemma. Ca2+-sensor(s) for release would be exposed to markedly different [Ca2+], varying by up to 13-fold, depending on their position around the vesicle. As a result, the precise position of Ca2+-sensor(s) with respect to the vesicle and the channel can be critical to determining the release probability. Variation in the position of Ca2+-sensor molecule(s) and their accessibility could be an important source of heterogeneity in vesicle release probability.

Novel players in pancreatic islet signaling: from membrane receptors to nuclear channels

Glucose and other nutrients regulate many aspects of pancreatic islet physiology. This includes not only insulin release, but also insulin synthesis and storage and other aspects of beta-cell biology, including cell proliferation, apoptosis, differentiation, and gene expression. This implies that in addition to the well-described signals for insulin release, other intracellular signaling mechanisms are needed. Here we describe the role of global and local Ca(2+) signals in insulin release, the regulation of these signals by new membrane receptors, and the generation of nuclear Ca(2+) signals involved in gene expression. An integrated view of these pathways should improve the present description of the beta-cell biology and provide new targets for novel drugs.

Secretory granule exocytosis

Regulated exocytosis of secretory granules or dense-core granules has been examined in many well-characterized cell types including neurons, neuroendocrine, endocrine, exocrine, and hemopoietic cells and also in other less well-studied cell types. Secretory granule exocytosis occurs through mechanisms with many aspects in common with synaptic vesicle exocytosis and most likely uses the same basic protein components. Despite the widespread expression and conservation of a core exocytotic machinery, many variations occur in the control of secretory granule exocytosis that are related to the specialized physiological role of particular cell types. In this review we describe the wide range of cell types in which regulated secretory granule exocytosis occurs and assess the evidence for the expression of the conserved fusion machinery in these cells. The signals that trigger and regulate exocytosis are reviewed. Aspects of the control of exocytosis that are specific for secretory granules compared with synaptic vesicles or for particular cell types are described and compared to define the range of accessory control mechanisms that exert their effects on the core exocytotic machinery.

In-vitro differentiation of pancreatic beta-cells

Stem cell biology is a new field that holds promise for in-vitro mass production of pancreatic beta-cells, which are responsible for insulin synthesis, storage, and release. Lack or defect of insulin produces diabetes mellitus, a devastating disease suffered by 150 million people in the world. Transplantation of insulin-producing cells could be a cure for type 1 and some cases of type 2 diabetes, however this procedure is limited by the scarcity of material. Obtaining pancreatic beta-cells from embryonic stem cells would overcome this problem. We have derived insulin-producing cells from mouse embryonic stem cells by a 3-step in-vitro differentiation method consisting of directed differentiation, cell-lineage selection, and maturation. These insulin-producing cells normalize blood glucose when transplanted into streptozotocin-diabetic mice. Strategies to increase islet precursor cells from embryonic stem cells include the expression of relevant transcription factors (Pdx1, Ngn3, Isl-1, etc), together with the use of extracellular factors. Once a high enough proportion of islet precursors has been obtained there is a need for cell-lineage selection in order to purify the desired cell population. For this purpose, we designed a cell-trapping method based on a chimeric gene that fuses the human insulin gene regulatory region with the structural gene that confers resistance to neomycin. When incorporated into embryonic stem cells, this fusion gene will generate neomycin resistance in those cells that initiate the synthesis of insulin. Not only embryonic, but also adult stem cells are potential sources for insulin-containing cells. Duct cells from the adult pancreas are committed to differentiate into the four islet cell types; other possibilities may include nestin-positive cells from islets and adult pluripotent stem cells from other origins. Whilst the former are committed to be islet cells but have a reduced capacity to expand, the latter are more pluripotent and more expandable, but a longer pathway separates them from the insulin-producing stage. The aim of this review is to discuss the different strategies that may be followed to in-vitro differentiate pancreatic beta-cells from stem cells.