Research Article  Open Access
Explicit Theoretical Analysis of How the Rate of Exocytosis Depends on Local Control by Ca^{2+} Channels
Abstract
Hormones and neurotransmitters are released from cells by calciumregulated exocytosis, and local coupling between Ca^{2+} channels (CaVs) and secretory granules is a key factor determining the exocytosis rate. Here, we devise a methodology based on Markov chain models that allows us to obtain analytic results for the expected rate. First, we analyze the property of the secretory complex obtained by coupling a single granule with one CaV. Then, we extend our results to a more general case where the granule is coupled with n CaVs. We investigate how the exocytosis rate is affected by varying the location of granules and CaVs. Moreover, we assume that the single granule can form complexes with inactivating or noninactivating CaVs. We find that increasing the number of CaVs coupled with the granule determines a much higher rise of the exocytosis rate that, in case of inactivating CaVs, is more pronounced when the granule is close to CaVs, while, surprisingly, in case of noninactivating CaVs, the highest relative increase in rate is obtained when the granule is far from the CaVs. Finally, we exploit the devised model to investigate the relation between exocytosis and calcium influx. We find that the quantities are typically linearly related, as observed experimentally. For the case of inactivating CaVs, our simulations show a change of the linear relation due to nearcomplete inactivation of CaVs.
1. Introduction
Molecules, e.g., neurotransmitters and proteins, are released from the cell by exocytosis [1]. In this paper, we focus on regulated exocytosis in the endocrine cells that release different kinds of hormones regulating various physiological processes [2]. When hormone secretion is defectively regulated, several diseases may develop. For example, in diabetes, the two main pancreatic hormones, insulin and glucagon, are not released appropriately for finetuning glucose homeostasis [3, 4]. Therefore, it is crucial to achieve a better understanding of the main mechanisms underlying hormone exocytosis that determines the control of different physiological processes.
In most endocrine cells, the hormones are contained in secretory granules that, in response to a series of cellular mechanisms culminating with an increase in the intracellular Ca^{2+} levels, fuse with the cell membrane and release the hormone molecules. The main mechanisms regulating hormone exocytosis are shared with exocytosis of synaptic vesicles underlying neurotransmitter release in neurons [1, 5]. The granules contain vSNARE proteins that can form the socalled SNARE complexes with tSNAREs inserted in the cell membrane [1]. SNARE complexes interact with other proteins, notably, Ca^{2+}sensing proteins such as synaptotagmins, which trigger exocytosis upon Ca^{2+} binding. Therefore, the local Ca^{2+} concentration at the Ca^{2+} sensor of the exocytotic machinery is a key factor determining the probability rate of exocytosis of the secretory granule [6].
Recently, we have devised a detailed model of Ca^{2+} dynamics and exocytosis for the glucagonsecreting pancreatic alphacells and showed how exocytosis is dependent on calcium dynamics, in particular, on calcium levels surrounding the Ca^{2+} channels (CaVs) [7], the socalled nanodomains [8]. Here, in order to characterize the local interactions between the single granule and the surrounding CaVs, we will exploit a strategy that is similar to the methodology devised in our recent paper to describe the large conductance BK potassium current that is controlled locally by CaVs [9]. We showed that the number and the type of CaVs coupled with the BK channel affect the electrical activity of neurons and other excitable cells, such as pancreatic betacells and pituitary cells. Therefore, we will implement mathematical modelling for characterizing the local interactions between granules and CaVs and, specifically, Markov chain models that could provide important insight into the exocytosis rate. In particular, by using the Markov chain theory [10], we will achieve analytic results for the expected rate and show how coupling different numbers and types of CaVs with the granule determines different responses.
2. Methods
2.1. CaV Channel Model
We model the Ca^{2+} channel by using the 3state Markov chain of Figure 1(a), where C corresponds to the closed state, O to the open state, and B to the inactivated (blocked) state of the calcium channel [11]. Then, the CaV model takes values in the state space and its transition rate or generator matrix is given bywhere and represent the voltagedependent Ca^{2+} channel opening rate and closing rate, respectively, and have the following forms:
(a)
(b)
(c)
(d)
The rate for channel inactivation, , is Ca^{2+}dependent and has the following form:where is the Ca^{2+} concentration at the Ca^{2+} sensor for inactivation and is given using reactiondiffusion theory [8, 12, 13] bywhere is the singlechannel Ca^{2+} current with the singlechannel conductance and the reverse potential, and represents the distance of the sensor for Ca^{2+}dependent inactivation from the channel pore. Finally, is the constant reverse reactivation rate. Table 1 reports the parameter values for the CaV model defined by above equations.

The deterministic description of the 3state Markov chain model for the CaV channel is given by the following ODE system:where the italic lowercase letters represent the corresponding state variables of the ODE model (h represents the fraction of Ca^{2+} channels not inactivated).
Finally, in order to investigate the relationship between exocytosis and Ca^{2+} loading, we compute the total charge entering via the Ca^{2+} channel at a given step voltage with time window, , as
2.2. Exocytosis Model
We assume a single granule, adjacent to the plasma membrane and primed for exocytosis, that can be in one of four different states depending on the number of Ca^{2+} ions bound to the Ca^{2+} sensor on the granule, likely synaptotagmin [14]: in with no bound Ca^{2+} ions, or in with one, or in with two, or in with three bound ions. Once it is in , the granule can fuse with the membrane and release its hormone content, assuming the final state Y [6, 15]. Therefore, we use a fivestate Markov chain model for describing exocytosis as shown in Figure 1(b), where the model takes values in the state space , and its transition rate or generator matrix is given bywhererepresents the Ca^{2+} binding rate, with the Ca^{2+} concentration at the granule sensor given by Equation (4) with being the distance from the CaV to the Ca^{2+} sensor on the granule. In the following, the distance from the CaV to the granule means the distance from the CaV to the Ca^{2+} sensor on the granule, which will be of the order of tens of nm. For comparison, secretory granules have diameters on the order 100–500 nm [16–19]. We assume a constant number of Ca^{2+} sensor molecules, which is therefore included in the binding parameter . The parameter is the unbinding rate, and is the fusion rate. Table 1 reports the parameter values.
The deterministic description of the 5state Markov chain model for exocytosis is given by the following ODE system:
For the above ODE model of Equations (9)–(13), we exploit quasi steadystate approximation for state , since its dynamics are fastest (the value of is much higher than those of the other parameters). Then, by renaming the state variables asby setting Equation (12) equal to zero yieldingand by summing Equations (11) and (12), we achieve a single ODE model for describing the dynamics of state variable and as follows:
The corresponding Markov chain model takes values in the state space (Figure 1(c)) and is described by the following generating matrix, :
Note that state of the Markov chain described by is an absorbing state: the process can never leave after entering it, reflecting that fusion is an irreversible process. Then can be rewritten aswheredescribes only the transitions between the transient states , , and and is a vector containing the transition intensities from the transient states to the absorbing state . The row vector consists entirely of 0’s since no transitions from to the transient states can occur. The remaining element of the matrix is 0 and gives the transition rate out of the absorbing state.
Using phasetype distribution results for Markov chains [10], we obtain an explicit formula for calculating the expected event rate to reach the absorbing state , given the initial probability row vector for the transient states (), aswhere .
2.3. GranuleCaV Complex Model with 1 : 1 and 1 : n Stoichiometries
2.3.1. 1 : 1 Stoichiometry
By coupling the CaV and exocytosis models, we obtain the 12state Markov chain model of Figure 1(d). The model takes values in the state spaceand its transition matrix, , is as follows:where is defined by Equation (1), by Equation (8) (Equation (15)) with , i.e., the concentration at the granule when the associated CaV is closed (or inactivated, i.e., ), and by Equation (8) (Equation (15)) with , i.e., the concentration at the granule when the associated CaV is open, computed by Equation (4). Then, the expected exocytosis rate for the single granule, , can be estimated by using Equation (20), assuming initially the granule in state and the CaV closed, i.e., the complex in the state , aswhere .
We also consider the particular case with noninactivating CaV (i.e., the Ca^{2+} channel can be only in C or in O). In this case, and is defined by Equation (1) with , and then , given by Equation (22), belongs to .
2.3.2. 1 : n Stoichiometry
In the following, we assume the case where the granule is coupled with more than one CaV. In particular, by considering k Ca^{2+} channels, we have a Markov chain model with possible states describing the k CaVs. In particular, the CaVs model takes values in the state space with and , and its generating matrix, , is given bywhereand .
Then, by coupling the CaVs and exocytosis models, we obtain a state Markov chain model. The model takes values in the state space , with , and , and its transition matrix, , can be written aswhere
Then, the expected exocytosis rate for the single granule coupled with k CaVs, , can be estimated by using Equation (20), assuming initially the granule in state and the k CaVs closed, i.e., the complex is initially in state , which yieldswhere .
For the particular case with noninactivating CaVs channels, by Equations (24) and (25) with , and then, , given by Equation (27), belongs to .
In order to compare the rate for a granule coupled with different number k of CaVs, we define the relative rate, , aswith . Moreover, in order to compare the rate at different distances from the granule to CaVs, we define the relative distance rate, , aswhere is the rate computed at a given distance and , the rate computed at nm.
3. Results and Discussion
We analyze the behavior of the devised exocytosis model where the single granule is coupled with k Ca^{2+} channels by using phasetype distribution results for Markov chains [10] (see Methods). First, we assume that a granule is coupled with one CaV and, then, we extend the results to a more general case with k CaVs. Moreover, we consider for both the cases (1 or k CaVs) that the granule forms complexes with inactivating or noninactivating CaVs. This scenario reflects, e.g., what it is observed in pancreatic betacells where the two main high voltageactivated Ca^{2+} channels, the L and P/Qtype Ca^{2+} channels, are examples of inactivating and noninactivating CaVs, respectively [20].
3.1. Granule Coupled with One Inactivating (or NonInactivating) CaV
Figure 2(a) shows the expected exocytosis rate, , computed by Equation (23), for a granule at different distances from an inactivating CaV channel. Independently of the distance to the CaV, the exocytosis rate has a bellshaped relation to voltage, as seen experimentally [20–22]. The same holds true in the case of noninactivating CaV (Figure 2(b)). As the distance between the granule and the Ca^{2+} channel increases, the expected rate decreases substantially and nonlinearly (for instance, in Figure 2(a), compare the red and blue lines for nm and nm, respectively). This is clearer from Figure 2(c), showing the relative distance rate defined by Equation (33) for different values of . Note that increasing the distance by a factor of two corresponds to a more than fivefold reduction of the exocytosis rate (the relative ratio is less than 0.2, see the red plot in Figure 2(c)). This steep dependence of the distance to the channel is because the calcium levels drop rapidly, moving away from the channel [8, 23].
(a)
(b)
(c)
(d)
We perform a similar analysis for the case where a granule is coupled with a noninactivating CaV (Figure 2(b)). We note an increase about of two orders of magnitudes for the exocytosis rate compared to the case with a granule coupled with an inactivating CaV (Figures 2(a) and 2(b)): the exocytosis proceeds more rapidly since the triggering Ca^{2+} signal is increased due to noninactivation of Ca^{2+} currents. Also in this case, the degree of decrease for the rate is much higher than the relative increase for the distance (Figure 2(d)). However, the benefit in terms of by reducing the distance is slightly less than that obtained with inactivating CaV (compare Figures 2(c) and 2(d)): for the case with inactivating CaV, it seems that moving away from the channel, decreases more due to the inactivation of CaV that determines a further drop of calcium levels.
3.2. Granule Coupled with k Inactivating (or NonInactivating) CaVs
Figures 3(a)–3(d) show the expected exocytosis rate computed by Equation (31), for a granule coupled with different numbers of inactivating CaVs and at fixed distances between the granule and the CaVs. It is clear that increasing the number of CaVs coupled with the granule determines a rise of the exocytosis rate. Moreover, as the number of CaVs coupled with the granule increases, the rise in the rate is more pronounced when the distance of the granule from the CaVs is small. This is evident by considering the relative rate defined by Equation (32) (Figure 3(e)). For instance, consider the cyan curves computed for with different types of lines denoting the different distances of the granule from the CaVs. In this case, the number of CaVs decreases by a factor of 2 (from 8 to 4) while the exocytosis rate drops more than threefold for nm (dashed cyan line, , for mV) and more than fivefold for nm (solid cyan line, , for mV).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
As done for the case with one CaV, we performed the same analysis with k noninactivating CaVs coupled with the granule (Figures 3(f)–3(i)). Also in this case, it is clear that increasing the number of CaVs determines a rise of the exocytosis rate for the granule. Surprisingly and in contrast with the case with inactivating CaVs, as the number of noninactivating CaVs increases, the relative rise in exocytosis rate is much higher at larger distances from the CaVs, as shown in Figure 3(j) reporting the relative rate . In case the number of CaVs is reduced from 8 to 4, the exocytosis rate decreases by 2–2.5fold when the granule is near the CaVs (see the solid cyan curve for nm, with mV mV), while it goes down fivefold when the granule is far from CaVs (see the dotted cyan curve for nm, with mV mV). It seems that when the granule is surrounded by more noninactivating CaVs, it is not necessary that the granule is very close to the CaVs for triggering exocytosis.
3.3. Relationship between Ca^{2+} Influx and Exocytosis
To investigate the relationship between exocytosis and Ca^{2+} loading, we consider a set of scenarios where the granule is coupled with different number of noninactivating or inactivating CaVs, placed very close (10 nm) or far (100 nm) from the granule. Figure 4(a) shows the calcium current at mV, for different numbers of noninactivating CaVs, while Figure 4(b) shows the corresponding cases with inactivating CaVs. In the latter, it is evident how the calcium influx drops after few tens of ms due to the inactivation of the CaVs. Figures 4(c) and 4(d) show the probability of exocytosis () vs. the integral of the Ca^{2+} current, , defined by Equation (6), for the granule placed close to the CaV cluster, for different numbers of CaVs ( nm). For the case of noninactivating CaVs (Figure 4(c)), raises linearly with , with slope that increases with the number of CaVs and then saturates due to the depletion of the granule pool as approaches 1 (see also [24]). For inactivating CaVs, we note a change of the slope of the linearity between and that is not only due to depletion (when ) but also to nearcomplete inactivation of CaVs, in particular after 50 ms (Figure 4(d)). Figures 4(e) and 4(f) show vs. when the granule is placed far from CaVs ( nm). Due to the distance to CaVs, the calcium concentration at the granule increases only modestly; hence, a greater calcium influx is needed to allow the granule to move through the Markov chain from to Y and undergoes exocytosis. This causes an evident initial delay for the granule to be released, resulting in an initial convex relation between and . After this initial phase, for the case of noninactivating CaVs (Figure 4(e)), raises linearly with with slope depending on the number of CaVs. For higher , the slope of slightly decreases in the case with CaVs reflecting slight depletion of the granule pool ( at fC). For inactivating CaVs (Figure 4(f)), as for the case with nm, we note a change of the linearity between and that is due to CaV inactivation.
(a)
(b)
(c)
(d)
(e)
(f)
4. Conclusions
In this paper, we devise a strategy that allows us to characterize the local interactions between granules and CaVs. The methodology is similar to our approach for modelling the local effect of CaVs on wholecell BK currents [9]. We develop Markov chain models describing the dynamics of a single granule coupled with one or more inactivating (or noninactivating) Ca^{2+} channels and use phasetype distribution results [10] for estimating the expected exocytosis rate.
We investigate how the release probability of a granule can be affected by varying the number of CaVs and the distance of the (Ca^{2+} sensor of the) granule from CaVs. In particular, from our analysis, we find that the distance between the granule and CaVs is a major factor in determining the exocytosis rate, as we recently demonstrated and quantified explicitly [23]. Further and in agreement with experiments [23], the simulations presented here show that the increase of the number of CaVs coupled with the granule determines a much higher rise of the exocytosis rate, which in the case of inactivating CaVs is more pronounced when the granule is close to CaVs ( nm), whereas for noninactivating CaVs the highest relative increase in rate is obtained when the CaVs are far from CaVs ( nm).
We also study the relationship between Ca^{2+} influx and exocytosis. The results of the devised exocytosis model confirm that the granule secretion is generally linearly related to the integral of Ca^{2+}current, as experimentally observed [25–29] and theoretically justified [24]. Surprisingly, for the case of inactivating CaVs, our analysis shows a change of the linear relation between and due to nearcomplete inactivation of CaVs. This fact is due to the rather complex exocytosis model where the efficacy of Ca^{2+} influx in triggering exocytosis depends on the number of active CaVs, as clearly seen in the case of noninactivating CaVs (Figures 4(c) and 4(e)), because of multiple steps of Ca^{2+} bindings before exocytosis. During inactivation, the effective number of CaVs declines, which has a similar effect as reducing the number of CaVs, and hence the slope of the relation between exocytosis and decreases. This finding reinforces the notion that a concave relation between exocytosis and Ca^{2+} influx does not necessarily reflect pool depletion [24] and provides a new example of such a scenario.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
F.M. was supported by the University of Padova (Research Grant BIRD 2017). M.G.P. was supported by the University of Padova (Research Project SID and Research Project PROACTIVE).
References
 R. D. Burgoyne and A. Morgan, “Secretory granule exocytosis,” Physiological Reviews, vol. 83, no. 2, pp. 581–632, 2003. View at: Publisher Site  Google Scholar
 S. Misler, “Unifying concepts in stimulussecretion coupling in endocrine cells and some implications for therapeutics,” Advances in Physiology Education, vol. 33, no. 3, pp. 175–186, 2009. View at: Publisher Site  Google Scholar
 S. E. Kahn, S. Zraika, K. M. Utzschneider, and R. L. Hull, “The beta cell lesion in type 2 diabetes: there has to be a primary functional abnormality,” Diabetologia, vol. 52, no. 6, pp. 1003–1012, 2009. View at: Publisher Site  Google Scholar
 M. Frances, M. Ashcroft, and P. Rorsman, “Diabetes mellitus and the betacell: the last ten years,” Cell, vol. 148, no. 6, pp. 1160–1171, 2012. View at: Publisher Site  Google Scholar
 S. Barg, “Mechanisms of exocytosis in insulinsecreting betacells and glucagonsecreting alphacells,” Pharmacology and Toxicology, vol. 92, no. 1, pp. 3–13, 2003. View at: Publisher Site  Google Scholar
 M. G. Pedersen, A. Tagliavini, G. Cortese, M. Riz, and F. Montefusco, “Recent advances in mathematical modeling and statistical analysis of exocytosis in endocrine cells,” Mathematical Biosciences, vol. 283, pp. 60–70, 2017. View at: Publisher Site  Google Scholar
 F. Montefusco and M. G. Pedersen, “Mathematical modelling of local calcium and regulated exocytosis during inhibition and stimulation of glucagon secretion from pancreatic alphacells,” Journal of Physiology, vol. 593, no. 20, pp. 4519–4530, 2015. View at: Publisher Site  Google Scholar
 E. Neher, “Vesicle pools and Ca^{2+} microdomains: new tools for understanding their roles in neurotransmitter release,” Neuron, vol. 20, no. 3, pp. 389–399, 1998. View at: Publisher Site  Google Scholar
 F. Montefusco, A. Tagliavini, M. Ferrante, and M. G. Pedersen, “Concise wholecell modeling of BKCaV activity controlled by local coupling and stoichiometry,” Biophysical Journal, vol. 112, no. 11, pp. 2387–2396, 2017. View at: Publisher Site  Google Scholar
 P. Buchholz, K. Jan, and I. Felko, “Phasetype distributions,” in Input Modeling with PhaseType Distributions and Markov Models, pp. 5–28, Springer International Publishing, New York, NY, USA, 2014. View at: Google Scholar
 A. Sherman, J. Keizer, and J. Rinzel, “Domain model for Ca^{2+}inactivation of Ca^{2+} channels at low channel density,” Biophysical Journal, vol. 58, no. 4, pp. 985–995, 1990. View at: Publisher Site  Google Scholar
 E. Neher, “Concentration profiles of intracellular Ca^{2+} in the presence of diffusible chelator,” in Calcium Electrogenesis and Neuronal Functioning, vol. 14, pp. 80–96, SpringerVerlag, Berlin, Germany, 1986. View at: Google Scholar
 D. H. Cox, “Modeling a Ca^{2+} channel/BK channel complex at the singlecomplex level,” Biophysical Journal, vol. 107, no. 12, pp. 2797–2814, 2014. View at: Publisher Site  Google Scholar
 P. S. Pinheiro, S. Houy, and J. B. Sørensen, “C2domain containing calcium sensors in neuroendocrine secretion,” Journal of Neurochemistry, vol. 139, no. 6, pp. 943–958, 2016. View at: Publisher Site  Google Scholar
 T. Voets, “Dissection of three Ca^{2+}dependent steps leading to secretion in chromaffin cells from mouse adrenal slices,” Neuron, vol. 28, no. 2, pp. 537–545, 2000. View at: Publisher Site  Google Scholar
 C. P. Grabner, S. D. Price, A. Lysakowski, and A. P. Fox, “Mouse chromaffin cells have two populations of dense core vesicles,” Journal of Neurophysiology, vol. 94, no. 3, pp. 2093–2104, 2005. View at: Publisher Site  Google Scholar
 S. A. Andersson, M. G. Pedersen, J. Vikman, and L. Eliasson, “Glucosedependent docking and SNARE proteinmediated exocytosis in mouse pancreatic alphacell,” Pflügers Archiv—European Journal of Physiology, vol. 462, no. 3, pp. 443–454, 2011. View at: Publisher Site  Google Scholar
 C. S. Olofsson, S. O. Göpel, S. Barg et al., “Fast insulin secretion reflects exocytosis of docked granules in mouse pancreatic betacells,” Pflügers Archiv, vol. 444, no. 12, pp. 43–51, 2002. View at: Publisher Site  Google Scholar
 P. M. Dean, “Ultrastructural morphometry of the pancreatic betacell,” Diabetologia, vol. 9, no. 2, pp. 115–119, 1973. View at: Publisher Site  Google Scholar
 M. Braun, R. Ramracheya, M. Bengtsson et al., “Voltagegated ion channels in human pancreatic betacells: electrophysiological characterization and role in insulin secretion,” Diabetes, vol. 57, no. 6, pp. 1618–1628, 2008. View at: Publisher Site  Google Scholar
 K. D. Gillis, R. Y. K. Pun, and S. Misler, “Single cell assay of exocytosis from adrenal chromaffin cells using “perforated patch recording”,” Pflugers Archiv European Journal of Physiology, vol. 418, no. 6, pp. 611–613, 1991. View at: Publisher Site  Google Scholar
 K. L. Engisch and M. C. Nowycky, “Calcium dependence of large densecored vesicle exocytosis evoked by calcium influx in bovine adrenal chromaffin cells,” Journal of Neuroscience, vol. 16, no. 4, pp. 1359–1369, 1996. View at: Publisher Site  Google Scholar
 N. R. Gandasi, P. Yin, M. Riz et al., “Ca^{2+} channel clustering with insulincontaining granules is disturbed in type 2 diabetes,” Journal of Clinical Investigation, vol. 127, no. 6, pp. 2353–2364, 2017. View at: Publisher Site  Google Scholar
 M. G. Pedersen, “On depolarizationevoked exocytosis as a function of calcium entry: possibilities and pitfalls,” Biophysical Journal, vol. 101, no. 4, pp. 793–802, 2011. View at: Publisher Site  Google Scholar
 T. Moser and E. Neher, “Rapid exocytosis in single chromaffin cells recorded from mouse adrenal slices,” Journal of Neuroscience, vol. 17, no. 7, pp. 2314–2323, 1997. View at: Publisher Site  Google Scholar
 S. Barg, X. Ma, L. Eliasson et al., “Fast exocytosis with few Ca^{2+} channels in insulinsecreting mouse pancreatic betacells,” Biophysical Journal, vol. 81, no. 6, pp. 3308–3323, 2001. View at: Publisher Site  Google Scholar
 Y. Z. De Marinis, S. Albert, C. E. Ward et al., “GLP1 inhibits and adrenaline stimulates glucagon release by differential modulation of N and Ltype Ca^{2+} channeldependent exocytosis,” Cell Metabolism, vol. 11, no. 6, pp. 543–553, 2010. View at: Publisher Site  Google Scholar
 R. Thiagarajan, J. Wilhelm, T. Tewolde, Y. Li, M. M. Rich, and K. L. Engisch, “Enhancement of asynchronous and trainevoked exocytosis in bovine adrenal chromaffin cells infected with a replication deficient adenovirus,” Journal of Neurophysiology, vol. 94, no. 5, pp. 3278–3291, 2005. View at: Publisher Site  Google Scholar
 M. G. Pedersen, V. A. Salunkhe, E. Svedin, A. Edlund, and L. Eliasson, “Calcium current inactivation rather than pool depletion explains reduced exocytotic rate with prolonged stimulation in insulinsecreting INS1 832/13 cells,” PLoS One, vol. 9, no. 8, Article ID e103874, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Francesco Montefusco and Morten Gram Pedersen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.