help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Poberaj, I.
Right arrow Articles by Zorec, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Poberaj, I.
Right arrow Articles by Zorec, R.

Biophys J, January 2002, p. 226-232, Vol. 82, No. 1

Modeling Excess Retrieval in Rat Melanotroph Membrane Capacitance Records

Igor Poberaj,dagger Marjan Rupnik,* Marko Kreft,*§ Sujit K. Sikdar,*Dagger and Robert Zorec*

 *Laboratory of Neuroendocrinology-Molecular Cell Physiology, Institute of Pathophysiology, Medical School, and  dagger Department of Physics, Faculty of Mathematics and Physics, 1001 Ljubljana, Slovenia;  Dagger Molecular Biophysics Unit, Indian Institute of Science, Bangalore 560012, India; and  §Celica, Biomedical Sciences Center, 1000 Ljubljana, Slovenia


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

We have used the patch-clamp technique to monitor changes in membrane capacitance (Cm) elicited by fast and spatially homogeneous rises in cytosolic calcium concentration ([Ca2+]i) using flash photolysis of NP-EGTA. Average peak [Ca2+]i amplitudes of 20-25 µM triggered three different types of responses in Cm: (i) In 42% of cells, a rise in [Ca2+]i activated a monotonic increase in Cm followed by a slow decline to resting values; (ii) In 30% of cells, the rise in Cm was clearly characterized by two dynamic components, consisting of a rapid and a slow exo-endocytosis cycle; (iii) In 28% of cells, after the initial rapid rise in Cm, endocytosis exhibited excess retrieval that was characterized by a decline in Cm below resting Cm. The aim of this work is to develop a unified mathematical model with a minimum number of parameters that would describe all the observed types of responses. Three models were considered: Model A, a model with a single component of exo-endocytosis cycle; model B, a model consisting of a sum of two independent dynamic components; and model C, a model in which, in addition to the two dynamic components as in model B, excess retrieval due to a lipid flow through the reversal closing of the fusion pore during the rapid component of exo-endocytosis cycle was considered. The results show that the latter model describes all the types of responses in Cm recorded in rat melanotrophs. The association of excess retrieval exclusively with the rapid, but not the slow, exocytosis indicates that some fusing vesicles mediate a lipidic flux during the reversal closing of the fusion pore, whereas those entering the slow phase of exocytosis may fuse with the plasma membrane completely and are retrieved by other endocytic machinery, independent of the lipid flow that might have occurred as the fusion pore opened permanently.


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

The application of caged-Ca2+ compounds to study Ca2+-dependent exocytosis by membrane capacitance (Cm) measurements has revealed multiple kinetic components (Kasai, 1999). There are two views about the nature of this kinetic diversity. First, a sequence of intermediates of a homogeneous population of vesicles may result in multiple kinetic components (Heinemann et al., 1994, Thomas et al., 1993, Xu et al., 1998). Second, heterogeneous populations of vesicles engaged in distinct pathways of exocytosis may also result in multiple kinetic components (Voets, 2000). Although the quantitative Cm measurements cannot readily distinguish the contributions of heterogeneous vesicles or vesicle intermediates, a combination with other approaches was used to address this problem. Using serotonin-loaded dense-core vesicles and amperometry, it was shown that, in pancreatic beta -cells, two populations of dense-core granules enter distinct pathways of exocytosis (Takahashi et al., 1997). Using myoballs to monitor the release of acetylcholine from PC12 cells, it was demonstrated that exocytosis of synaptic vesicles and dense-core granules consists of distinct pathways (Ninomiya et al., 1997). By using an antibody to Ca2+-dependent activator protein for secretion (CAPS), a neural/endocrine-specific 145-kD protein, originally characterized as a brain cytosolic factor that reconstitutes Ca2+-dependent secretion in permeable neuroendocrine cells (Walent et al., 1992), it was shown that the rapid component of secretory response in rat melanotrophs was selectively abolished (Rupnik et al., 2000), whereas an antibody selective for the heterotrimeric Galpha i3 attenuated the slow phase of exocytosis (Kreft et al. 1999), suggesting that the two kinetically and biochemically distinct phases of Cm increase represent distinct exocytotic pathways in melanotrophs. To understand the nature of this highly complex process, a suitable model of membrane turnover dynamics would be of great benefit. Therefore, the aim of this work is to establish a model describing flash-induced time-dependent changes in Cm.

Using rat melanotrophs, which secrete pro-opiomelanocortin-derived peptides via dense-core vesicle exocytosis (Mains and Eipper, 1979), we have used Cm measurements (Neher and Marty, 1982) combined with flash photolysis to deliver rapid and spatially homogeneous steps in cytosolic Ca2+ (Neher and Zucker, 1993), which elicit multiple kinetic components in secretory activity. As reported previously (Thomas et al., 1993, 1994; Rupnik et al., 2000), we confirmed that multiple kinetic components characterized these responses in melanotrophs. It was shown previously that distinct kinetic components display distinct biochemical mechanisms (Kreft et al., 1999; Rupnik et al., 2000), therefore we modeled the time course in Cm to be due to one or a sum of two independent endo-exocytosis cycles. Furthermore, excess retrieval was observed after the rapid rise in Cm in ~28% of recordings, as reported by Thomas et al. (1994). We model this process by assuming it is due to a lipid flux during reversible fusion-pore opening (Monck et al., 1990). The results show that the model consisting of two dynamic components of exo-endocytosis cycle, with lipid flow taken into account, fits experimental results significantly better than other tested models. Excess retrieval observed after the rapid exocytosis indicates the existence of a functional state of docked vesicles, which mediates a significant lipid flux during the reversal closing of the fusion pore.


    MATERIALS AND METHODS
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Cell preparation

Melanotrophs were prepared as described by Rupnik and Zorec (1992).

Compensated membrane capacitance measurements

Compensated membrane capacitance measurements were performed as described (Neher and Marty, 1982; Zorec et al. 1991; Rupnik and Zorec; 1995), using a SWAM CELL or SWAM IIB patch-clamp-lock-in amplifier (Celica, Ljubljana, Slovenia), operating at 1.6 kHz. Upon establishment of the whole-cell configuration, Cm and Ga (access conductance) were compensated by Cslow and Ga compensation controls. Sine voltage of 11-111 mVrms was applied. The phase angle setting was determined by applying a 100 fF or 1 pF calibration pulse while monitoring the pulse on the C (signal proportional to Cm) and G outputs of the lock-in amplifier. These two signals were stored unfiltered (C-DAT4, Cygnus Technology, Inc., Delaware Water Gap, PA) for off-line analysis. Simultaneously, we recorded filtered (300 Hz, 4 pole Bessel) C and G signals, the fluorescence intensity from a C660 photon counter (Thorn EMI, Middlesex, UK) and membrane current (0-10 Hz, low pass Bessel). PhoCal program (LSR, Cambridge, UK) was used to acquire signals every 5 ms. For temporally high resolution measurements of Cm, the data recorded on DAT recorder were played back in 10-s sequences and digitized at 50 kHz using a CDR program (J. Dempster, Strachclyde, UK). Signals were digitally low-pass filtered at 1 kHz (2-way 150th order FIR filter, Math Works MATLAB, Math Works, Natik, MA) and resampled at 10 kHz. The pipette solution contained (in mM): KCl 110, TEACl 10, KOH/HEPES (N-2-hydroxyethylpiperazine-N'-2-ethanesulphonic acid) 40, Na2ATP 2, MgCl2 2, K4-NP-EGTA [o-nitrophenyl ethyleneglycol-bis-(beta -aminoethyl-ether)-N,N,N',N'-tetrapotassium salt, Molecular Probes, Eugene, OR] 4, CaCl2 3.6, furaptra 0.5, pH = 7.2. The bath contained (in mM): NaCl 131.8, KCl 5, MgCl2 2, NaH2 PO4 0.5, NaHCO3 5, Na HEPES 10, D-glucose 10, CaCl2 1.8, pH = 7.2. Recordings were made at room temperature. Pipette resistances ranged from 1 to 4 MOmega .

Flash photolysis and [Ca2+]i measurements

We used o-nitrophenyl ethyleneglycol-bis-(beta -aminoethyl-ether)-N,N,N',N'-tetraacetic acid (NP-EGTA) (Molecular Probes) to manipulate [Ca2+]i. A UV flash from a Xe arc flash lamp (Hi-Tech, Salisbury, UK) illuminated cells through a 40× fluor oil immersion objective of a Nikon Diaphot microscope. The same optical pathway was used to illuminate the fluorescent [Ca2+] indicator furaptra (Molecular Probes). Calibration of [Ca2+]i measurements was performed in each cell (Carter and Ogden, 1994; Kreft et al., 1999) by using the autofluorescence in a cell-attached configuration and the fluorescence in a resting whole-cell recording.

Data analysis and modeling

Measured data were analyzed with the Microcal Origin program. First, the time interval of the measurements was truncated to the same length chosen to be five seconds after flash photolysis excitation. Then the data were fitted with functions corresponding to different models. The functions were coded in C and compiled in DLL which was used by the Nonlinear Least Square Fitting (NLSF) procedure built into Origin. Statistics are in the format mean ± SEM.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Classification of flash-induced time-dependent changes in Cm

Resting membrane capacitance (Cm) was stable, averaging 5.3 ± 0.2 pF (mean ± SEM, n = 50), as reported by (Zupancic et al., 1994). Figure 1 shows a typical time course for [Ca2+]i and Cm during flash photolysis of an NP-EGTA-loaded cell. The arrow on the top indicates the flash (150 V/3.9 mF discharge) delivery to the patched cell through the objective, which photolyzed around 10% of the Ca2+-bound NP-EGTA (see Materials and Methods). This flash transiently increased [Ca2+]i to 25 µM (Fig. 1, top trace), after which [Ca2+]i dropped back to the baseline with an exponential time constant of ~4 s, due to extrusion from the cell, and due to diffusion into the pipette and mixing with the unphotolyzed cage. After flash delivery, Cm increased by 805 fF and later slowly decreased due to slow endocytosis toward resting Cm of 4.8 pF (Fig. 1, middle trace). The prominent increase in Cm was preceded by a smaller (160 fF) transient rise in Cm, termed rapid exocytosis (Rupnik et al., 2000), which is clearly resolved in Fig. 2 B at a higher time resolution. The rapid phase is followed by a slower rise in Cm, which was termed slow exocytosis (Rupnik et al., 2000). Conductance trace (Fig. 1, bottom trace), reflecting contributions of access conductance, membrane conductance and membrane capacitance (Lindau and Neher, 1988) was not correlated to changes in Cm. When the NP-EGTA filling the cell contained no Ca2+ (six cells), the flash did not change [Ca2+]i and both Cm and conductance remained constant (data not shown).



View larger version (11K):
[in this window]
[in a new window]
 
FIGURE 1   Time-dependent changes in [Ca2+]i (top trace), membrane capacitance (Cm), and conductance (G, reflecting changes in access conductance, membrane conductance, and Cm) after flash photolysis of Ca2+-loaded NP-EGTA (arrow). Dashed lines in the middle and bottom traces indicate resting capacitance of 4.8 pF and access conductance of 150 nS, respectively. Both parameters were determined using the controls on the SWAM patch-clamp amplifier.



View larger version (9K):
[in this window]
[in a new window]
 
FIGURE 2   Types of responses in flash-induced increases in Cm. Top: Type A response is a monotonic increase in Cm, which resembles a single slow phase response. Middle: Type B response consists of two dynamically distinct phases in Cm. Bottom: Type C response is similar to type B response, but the signal exhibits excess retrieval after the rapid rise in Cm. Curves with dots represent best fits to Eq. 3, where the parameters were as follows: Top: B0F = 0.61 pF, B0S = 3.62 pF, alpha  = 1.78, kF = 1.34 s-1, kexoS = 0.21 s-1, kendoS = 0.2 s-1; middle: B0F = 0.14 pF, B0S = 1.04 pF, alpha  = 5.04, kF = 6.49 s-1, kexoS = 1.27 s-1, kendoS = 0 s-1; bottom: B0F = 0.27 pF, B0S = 0.67 pF, alpha  = 2.98, kF = 4.43 s-1, kexoS = 1.01 s-1, kendoS = 0 s-1.

Membrane capacitance measurements were performed on 54 cells with a transient increase in [Ca2+]i similar to that seen in Fig. 1 (ranging between 10 and 40 µM). Experiments revealed two other types of changes in Cm during the period of 5 s after the flash delivery (Fig. 2, A and C). Therefore, the responses in Cm were grouped into three types (Fig. 2). Type A: In 23 cells (42%), Cm responses consisted of a monotonic increase in Cm followed by a decline in Cm similar to that seen in Fig. 1 (middle trace). The increase in Cm resembled the slow phase of exocytosis (Rupnik et al., 2000). The presence of an initial rapid component in these responses cannot be excluded, however it was not clearly visible, possibly due to a smaller amplitude representing only 10% of the amplitude of the slow component (see Rupnik et al., 2000). Type B: In 16 cells (30%), Cm rise consisted of a superposition of two clearly visible dynamic components of exo-endocytosis cycle, the rapid and the slow components (Rupnik et al., 2000). Type C: In 15 cells (28%), Cm responses exhibited excess retrieval after the initial rise in Cm. Peak values in [Ca2+]i eliciting different types of responses were not significantly different: 26.6 ± 0.8 µM (n = 23) for Type A; 23.5 ± 5.6 µM (n = 15) for type B; and 19.9 ± 5.8 µM (n = 16) for Type C response.

Fitting different models to the time-courses in Cm

To describe the time course of the three types of responses in Cm (Fig. 2), we consider three models (see Fig. 3 and Appendix). Model A is the simplest and takes only one dynamic component of exo-endocytosis cycle into account (Fig. 3, top trace). Note that its applicability is limited only to Cm responses of type A (Fig. 2). To describe responses of type B, Model A is expanded to Model B by adding a second independent dynamic component of exo-endocytic cycle (Fig. 3, middle trace). The sum of two independent components implies that the whole process is parallel rather than sequential (Rupnik et al., 2000).



View larger version (17K):
[in this window]
[in a new window]
 
FIGURE 3   Diagram of models of exo-endocytosis cycle in rat melanotrophs. The models are based on the vesicle evolution scheme, where vesicles are grouped into three different pools A, B, and C, according to their evolutionary state. Pool A contains vesicles that have not reached release-ready state yet and thus cannot directly contribute to the measured signal. The contribution of pool A to the capacitance signals in our models can be neglected. Pool B contains vesicles that can be immediately released upon Ca2+ stimulation. Released vesicles enter pool C. Measured capacitance is proportional to the number of vesicles in the pool C. Model A (top) consists of one exo-endocytosis cycle, kexo, kend represent rate constants. Model B (middle) consists of two parallel exo-endocytosis cycles, where subscripts kS, kF denote slow and fast exo-endocytic pathways, respectively. Model C (bottom) is similar to Model B, but a process of excess retrieval was considered to take place in the fast exo-endocytic cycle, depicted by the dashed curve.

Excess retrieval was observed exclusively in conjunction with the rapid phase of the Cm increase (Fig. 2), consistent with previous reports (Thomas et al., 1994). It has been reported that, during pore opening, a substantial lipid flow between the plasmalemma and the vesicle membrane can take place, which is driven by a tension gradient between the two membranes (Monck et al., 1990; Chizmadhzev et al., 2000). To describe type C responses in Cm, we refine Model B into Model C by taking the lipid flow associated with the reversal opening of a fusion pore during rapid exocytosis into account (Fig. 3, bottom, dashed line).

Functional representations of the three models (derived in Appendix) are presented by the Eqs. 1-3: Model A
C<SUB><UP>m</UP></SUB>(t)=B<SUB><UP>o</UP></SUB> <FR><NU><IT>k</IT><SUB><UP>exo</UP></SUB></NU><DE>k<SUB><UP>exo</UP></SUB>−k<SUB><UP>endo</UP></SUB></DE></FR> (e<SUP><UP>−k<SUB>endo</SUB></UP><UP>t</UP></SUP>−e<SUP><UP>−k<SUB>exo</SUB>t</UP></SUP>) (1)
Cm denotes measured Cm, B0 is the initial vesicle pool size, kexo and kendo are rate constants for exo- and endocytosis (Fig. 3, top trace).
Model B
C<SUB><UP>m</UP></SUB>(<UP>t</UP>)<UP>=B<SUB>0F</SUB></UP>k<SUB><UP>F</UP></SUB>te<SUP><UP>−k<SUB>F</SUB></UP><UP>t</UP></SUP> (2)

<UP>+B<SUB>0S </SUB></UP><FR><NU>k<SUB><UP>exoS</UP></SUB></NU><DE>k<SUB><UP>exoS</UP></SUB>−k<SUB><UP>endoS</UP></SUB></DE></FR> (e<SUP><UP>−k<SUB>endoS</SUB></UP><UP>t</UP></SUP>−e<SUP><UP>−k<SUB>endoS</SUB></UP><UP>t</UP></SUP>)
Cm denotes measured capacitance. B0F and kF are the initial vesicle pool size and rate constant of the rapid phase. B0S, kexoS, and kendoS are the initial vesicle pool size and respective rate constants for the slow phase exo- and endocytosis. It was determined experimentally by fitting that the rate constants for rapid exo- and endocytosis are the same within the experimental error.
Model C
C<SUB><UP>m</UP></SUB>(t)=B<SUB><UP>0F</UP></SUB>(&agr;(k<SUB><UP>F</UP></SUB>t+1)e<SUP><UP>−k<SUB>F</SUB></UP><UP>t</UP></SUP>−e<SUP>−k<SUB>F</SUB><UP>t</UP></SUP>−&agr;+1) (3)

+B<SUB><UP>0S</UP></SUB> <FR><NU>K<SUB><UP>exoS</UP></SUB></NU><DE>k<SUB><UP>exoS</UP></SUB>−k<SUB><UP>endoS</UP></SUB></DE></FR> (e<SUP><UP>−k<SUB>endoS</SUB></UP><UP>t</UP></SUP>−e<SUP><UP>−k<SUB>exoS</SUB></UP><UP>t</UP></SUP>)
Cm denotes measured capacitance. B0F and kF are the initial vesicle pool size and rate constant of the fast phase, respectively. alpha  is the ratio of vesicle membrane area after endocytosis and before exocytosis. B0S, kexoS, and kendoS are the initial vesicle pool size and respective rate constants for slow exo- and endocytosis

The models were tested by fitting the corresponding functions to experimental data and by comparing normalized chi 2 of the fits. In addition, the quality of the fits was examined visually (Fig. 4). Because Models A and B are subsets of Model C, it was not surprising that Model C produced lower chi 2 values than Models A or B. However, only Model C reproduced with high accuracy all measured membrane capacitance records. In contrast, Models A and B completely failed to reproduce excess retrieval (group C) even qualitatively. The only difference between Models B and C is that the latter contains an additional parameter, alpha , describing the lipid flow through a fusion pore, which is created during exocytosis. We would like to point out that, within the framework of the presented modeling, excess retrieval cannot be reproduced by simply increasing the number of dynamic components unless lipid flow through a fusion pore is taken into account. Although a model with a higher number of dynamic exo-endocytic components would consist of a higher number of free-fitting parameters than Model C, it would certainly fail to reproduce excess retrieval. Thus, we conclude that, within the presented scheme, Model C has the minimum number of fitting parameters that can reproduce all the measured flash-induced membrane capacitance responses in rat melanotrophs with high fidelity.



View larger version (16K):
[in this window]
[in a new window]
 
FIGURE 4   Fitting membrane capacitance changes of Type B (top) to equations representing Models A, B, and C (see inset for expanded traces). Residuals of fits to respective models are shown below. Note that the best fit was obtained by fitting the equation representing Model C to the data. Parameters of the curves shown on the top panel were as follows: Model C (dots on line): B0F = 0.14 pF, B0S = 1.04 pF, alpha  = 5.04, kF = 6.49 s-1, kexoS = 1.27 s-1, kendoS = 0; Model B (dashed line): B0F = 0.13 pF, B0S = 0.55 pF, alpha  = 1, kF = 16.3 s-1, kexoS = 0. 57 s-1, kendoS = 0; Model A (thin line): B0F = 0, B0S = 0.55 pF, alpha  = 1, kF = 0, kexoS = 0.57 s-1, kendoS = 0.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

The aim of this study was to develop a model describing time-dependent changes in Cm elicited by photolysis of caged calcium (NP-EGTA) in rat melanotrophs. As reported previously (Thomas et al., 1993; 1994; Rupnik et al., 2000), recorded changes in Cm consist of multiple kinetic components. Here, we observed that the time course of Cm during the first 5 s immediately after the UV flash delivery consists of three types of responses, termed type A (monotonic increase in Cm), type B (biphasic increase in Cm), and type C (biphasic increase in Cm with excess retrieval after the first rise in Cm) (Fig. 2). These different categories of responses in Cm appeared to be independent of the peak [Ca2+]i.

Multiphasic Cm responses in neuroendocrine cells have generally been interpreted by a sequential model in which vesicles staged at varying distances from exocytosis undergo progressive fusion (Heinemann et al., 1994; Thomas et al., 1993, Xu et al., 1998). The initial Cm increase is presumed to reflect triggered exocytosis of those vesicles that are fusion competent at the time of the Ca2+ increase, whereas later phases of the Cm increase are assumed to reflect the progressive fusion of vesicles that require time-dependent recruitment, docking, or priming steps. However, results showing the selective CAPS-antibody inhibition of the rapid rise in Cm increase (Rupnik et al., 2000) are not easily compatible with a strictly sequential model. Moreover, one cannot neglect the excess retrieval, which appears in approximately 30% of all experiments (Thomas et al., 1994; and Fig. 2) and is also difficult to explain within the framework of a sequential model. Thus, we considered a hypothesis in which rapid and slow exocytosis are mediated by two distinct parallel pathways that use a common pool of vesicles (Fig. 3). Mathematical models of vesicle secretion dynamics based on the above hypothesis were developed and used in quantitative data analysis and interpretation. Model C, which consists of two dynamic components of exo-endocytosis cycle, was shown to reproduce with high fidelity all types of measured responses in rat melanotrophs (Figs. 2 and 4), an observation supporting previously reported existence of biochemically dissimilar mechanisms of exocytosis in melanotrophs (Rupnik et al., 2000) and consistent with a report on chromaffin cells from tissue slices (Voets, 2000). The existence of multiple pathways of exocytosis also appears to be present in other cell types, such as adipocytes (Bogan and Lodish, 1999) and neutrophils (Nusse et al., 1998). In the latter cell type, parallel pathways of exocytosis are associated with distinct calcium sensitivity, which may indicate that vesicles undergoing different exocytic pathways are characterized by particular molecular mechanisms. Indeed, two calcium-sensitive pathways of exocytosis have been associated with different biochemical characteristics (Rupnik et al., 2000). It is likely that different calcium-sensitive mechanisms of glutamate secretion from bipolar neurones (Heidelberger et al., 1994) and from calyx-type synapses (Schneggenburger and Neher, 2000; Bollmann et al., 2000) reflect differences in the molecular organization of exocytic apparatus in these types of neurons.

A possibly important aspect of this work is the attempt to include, in modeling of time-dependent changes in membrane capacitance, a flux of lipids through a reversal closing of the fusion pore (Monck et al., 1990). Inclusion of this process in the model with two dynamic components of exo-endocytosis cycles improves the fitting fidelity in such a way that the model describes well not only distinct types of responses in Cm, but can describe all types of responses observed in melanotrophs (Figs. 2, 3, and 4). The flux of lipids is driven by the tension gradient between the vesicle membrane and the plasmalemma (Monck et al., 1990; Solsona et al., 1998; Chizmadzhev et al., 1999; 2000). Excess retrieval was exclusively observed to be associated with the rapid exocytosis in rat melanotrophs (Fig. 2; Thomas et al., 1994). This may indicate that the tension gradient between the vesicle and the plasma membrane is different for vesicles in the rapid and slow exocytotic pools. Alternatively, vesicles entering the slow phase of exocytosis may be fusing with the plasma membrane permanently and are retrieved by other endocytic machinery, which then controls the size of the retrieved membrane area, independent of the lipid flow that might have occurred as the fusion pore opened. The mechanism of such a difference in the properties of vesicles in the rapid and slow exocytotic pools is not known. On one hand, different proteins or lipidic structures between interacting membranes may contribute to the two functional populations of vesicles in melanotrophs. Distinct biochemical characteristics of rapid and slow exocytosis (Rupnik et al., 2000; Kreft et al., 1999) support such a mechanism. On the other hand, an interesting question to be elucidated is whether the rapid component of Cm increase reflects exocytosis of vesicles, of which the fusion pore expands more quickly in comparison to vesicles undergoing slow exocytosis. A recent mathematical analysis showed that fusion-pore growth is not affected by tension-driven lipid flux from one membrane to another (Chizmadzhev et al., 2000). It is therefore likely that vesicles undergoing rapid exocytosis are "frozen" in a functional state that favors lipidic flux, hindering the dissipation of tension difference after establishment of the fusion pore. Hence, in such a relatively stable state (i.e., reversible and repetitive fusion-pore opening and closing) the time for lipidic flux to take place would be longer in comparison to the functional state of vesicles undergoing slow exocytosis where the fusion pore may open completely and permanently.

In summary, we have studied Cm elicited by flash photolysis of caged Ca2+. The different types of responses appeared not to be associated with differences in cytoplasmic calcium concentration, an observation consistent with previous reports (Thomas et al., 1993, 1994). To interpret the complex responses, we considered a new mathematical model consisting of two dynamic cycles of exo-endocytosis with a component representing lipid flow during fusion-pore opening after the onset of the rapid rise in Cm. We have demonstrated that the model describes all measured responses with high fidelity. Our results are consistent with the view that complex kinetic changes in Cm are due to distinct pathways of regulated exocytosis. Unlike vesicles entering slow exocytosis, vesicles in rapid exocytosis are characterized by a functional state during which a significant lipid flux can take place.


    APPENDIX
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

For this paper to be self-contained, we hereby derive functions describing capacitance signals as a function of time for different models. The models are based on the vesicle evolution scheme presented in Fig. 3, where vesicles are grouped into three different pools A, B, and C according to their evolutionary state. Pool A contains vesicles that have not reached release-ready state yet and thus cannot directly contribute to the measured signal. The contribution of pool A to the capacitance signals in our models can be neglected. Pool B contains vesicles that can be immediately released upon Ca2+ stimulation. Released vesicles enter pool C. Measured capacitance is proportional to the number of vesicles in pool C.

We start our modeling with the simplest case of monophasic slow exocytic response of type A. We assume that the vesicle transition dynamics between different pools is governed by the following set of rate equations:
<FR><NU><UP>d</UP>N<SUB><UP>B</UP></SUB></NU><DE><UP>d</UP>t</DE></FR>=−k<SUB><UP>exo</UP></SUB>N<SUB><UP>B</UP></SUB><UP>,</UP> (A1)

<FR><NU><UP>d</UP>N<SUB><UP>C</UP></SUB></NU><DE><UP>d</UP>t</DE></FR>=k<SUB><UP>exo</UP></SUB>N<SUB><UP>B</UP></SUB>−k<SUB><UP>endo</UP></SUB>N<SUB><UP>C</UP></SUB><UP>,</UP> (A2)
where NB, and NC are respective numbers of vesicles in pools B and C. The constants kexo and kendo are rate constants for exo- and endocytosis, respectively. The solution of the above set of equations is given by
N<SUB><UP>C</UP></SUB>(t)=<FR><NU>N<SUB><UP>B<SUB>0</SUB></UP></SUB>k<SUB><UP>exo</UP></SUB></NU><DE>k<SUB><UP>exo</UP></SUB>−k<SUB><UP>endo</UP></SUB></DE></FR> (e<SUP><UP>−k<SUB>endo</SUB></UP><UP>t</UP></SUP>−e<SUP><UP>−k<SUB>exo</SUB></UP><UP>t</UP></SUP>). (A3)
Multiplying Eq. A3 by the capacitance of a single vesicle gives a time dependent capacitance Cm(t) signal, Eq. 1.

A similar set of equations with different coefficients was applied to the rapid phase. However, fitting of the rapid phase data revealed that kexo and kendo for the rapid phase are equal within experimental error. The solution for the rapid-phase capacitance contribution can be found as a limiting value of the Eq. 1 and is given by
C<SUB>m</SUB>(t)=B<SUB>0F</SUB>k<SUB>F</SUB>te<SUP>−k<SUB>F</SUB><UP>t</UP></SUP>, (A4)
where B0F and kF represent initial pool size and rate constant, respectively. Assuming that rapid and slow phase secretory responses evolve in parallel independent pathways, the final solution (Eq. 2) can be readily obtained by summing Eqs. 1 and A4.

Until now, we have not taken into account the possibility that added and retrieved membrane areas per vesicle are not the same. This may occur if the tensions in plasma and vesicle membrane are different, and the fusion pore opens temporarily. In such a case, a significant lipid flow through a fusion pore in a direction of a tension gradient is present. Usually, vesicle-membrane tension is higher than plasma-membrane tension (Monck et al., 1990).

The change in the area of the plasma membrane as a function of time can be written as
&Dgr;S(t)=S<SUB>0</SUB>N<SUB>exo</SUB>(t)−S<SUB>endo</SUB>N<SUB>endo</SUB>(t)=S<SUB>0</SUB>(N<SUB>exo</SUB>(t)−&agr;N<SUB>endo</SUB>(t)), (A5)
where Nexo(t) is the number of vesicles that were fused and Nendo(t), the number of vesicles that were retrieved by time t. S0 and Sendo are the initial and the final vesicle membrane areas. The parameter alpha  is defined by the ratio
&agr;=<FR><NU>S<SUB>endo</SUB></NU><DE>S<SUB>0</SUB></DE></FR>. (A6)
Nexo(t) and Nendo(t) can be calculated from Eqs. A1 and A2. We will do it for the rapid phase, which is associated with excess retrieval. The result is given by
N<SUB>exo</SUB>=N<SUB>B<SUB>0</SUB></SUB>(1−e<SUP>−k<SUB>F</SUB><UP>t</UP></SUP>), (A7)

N<SUB>endo</SUB>=N<SUB>B<SUB>0</SUB></SUB>(1−e<SUP>−k<SUB>F</SUB><UP>t</UP></SUP>−k<SUB>F</SUB>te<SUP>−k<SUB>F</SUB><UP>t</UP></SUP>). (A8)
From Eqs. A6, A7, and A8, an expression for the plasma membrane surface as a function of time,
&Dgr;S<SUB>F</SUB>(t)=S<SUB>0F</SUB>N<SUB>B<SUB>oF</SUB></SUB>(&agr;k<SUB>F</SUB>te<SUP>−k<SUB>F</SUB><UP>t</UP></SUP>+(&agr;−1)(e<SUP>−k<SUB>F</SUB><UP>t</UP></SUP>−1)) (A9)
is obtained, which is, again, proportional to the measured capacitance. By adding the expression for the slow phase, one obtains the complete solution for the capacitance signal as a function of time (Eq. 3), which is the most general one and contains Eqs. 1, A4, and A6, which describe signals of type A and B as special cases, respectively. Eq. A6 is obtained by setting parameter alpha  = 1, corresponding to zero lipid flow through a fusion pore. Similarly, Eq. A4 can be obtained by setting the amplitude of the rapid phase B0F = 0.

    ACKNOWLEDGMENTS

This work was supported by a Ministry of Sciences and Technology of The Republic of Slovenia (#P3 521 381 and #J3 2344-7421-00) awarded to R.Z. and M.K. We thank S. Grilc for cell cultures.

    FOOTNOTES

Received for publication 23 October 2000 and in final form 1 October 2001.

Address reprint requests to Robert Zorec, PAFI, Institute of Pathophysiology, Medical School, Laboratory of Neuro-Endocrinology-Molecular Cell Physiology, P.O.B. 2211, 1001 Ljubljana, Slovenia. Tel.: +386-1-543-70-20; Fax: +386-1-543-70-21; E-mail: robert.zorec{at}mf.uni-lj.si.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Biophys J, January 2002, p. 226-232, Vol. 82, No. 1
© 2002 by the Biophysical Society   0006-3495/02/01/226/07  $2.00



This article has been cited by other articles:


Home page
J. Biol. Chem.Home page
J. H. Kim, J. H. Nam, M.-H. Kim, D.-S. Koh, S.-J. Choi, S. J. Kim, J. E. Lee, K. M. Min, D.-Y. Uhm, and S. J. Kim
Purinergic Receptors Coupled to Intracellular Ca2+ Signals and Exocytosis in Rat Prostate Neuroendocrine Cells
J. Biol. Chem., June 25, 2004; 279(26): 27345 - 27356.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
M. Kreft, D. Krizaj, S. Grilc, and R. Zorec
Properties of Exocytotic Response in Vertebrate Photoreceptors
J Neurophysiol, July 1, 2003; 90(1): 218 - 225.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Cell Physiol.Home page
M. Kreft, V. Kuster, S. Grilc, M. Rupnik, I. Milisav, and R. Zorec
Synaptotagmin I increases the probability of vesicle fusion at low [Ca2+] in pituitary cells
Am J Physiol Cell Physiol, February 1, 2003; 284(2): C547 - C554.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Poberaj, I.
Right arrow Articles by Zorec, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Poberaj, I.
Right arrow Articles by Zorec, R.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2002 by the Biophysical Society.