Original Research ARTICLE
Estimation of presynaptic calcium currents and endogenous calcium buffers at the frog neuromuscular junction with two different calcium fluorescent dyes
- 1Laboratory of the Biophysics of Synaptic Processes, Kazan Scientific Centre, Kazan Institute of Biochemistry and Biophysics, Russian Academy of Sciences, Kazan, Russia
- 2Open Laboratory of Neuropharmacology, Kazan Federal University, Kazan, Russia
- 3Department of Radiophotonics and Microwave Technologies, Kazan National Research Technical University named after A. N. Tupolev, Kazan, Russia
- 4Department of Neurobiology and Radioelectronics, Kazan Federal University, Kazan, Russia
- 5Department of Medical and Biological Physics, Kazan State Medical University, Kazan, Russia
At the frog neuromuscular junction, under physiological conditions, the direct measurement of calcium currents and of the concentration of intracellular calcium buffers—which determine the kinetics of calcium concentration and neurotransmitter release from the nerve terminal—has hitherto been technically impossible. With the aim of quantifying both Ca2+ currents and the intracellular calcium buffers, we measured fluorescence signals from nerve terminals loaded with the low-affinity calcium dye Magnesium Green or the high-affinity dye Oregon Green BAPTA-1, simultaneously with microelectrode recordings of nerve-action potentials and end-plate currents. The action-potential-induced fluorescence signals in the nerve terminals developed much more slowly than the postsynaptic response. To clarify the reasons for this observation and to define a spatiotemporal profile of intracellular calcium and of the concentration of mobile and fixed calcium buffers, mathematical modeling was employed. The best approximations of the experimental calcium transients for both calcium dyes were obtained when the calcium current had an amplitude of 1.6 ± 0.08 pA and a half-decay time of 1.2 ± 0.06 ms, and when the concentrations of mobile and fixed calcium buffers were 250 ± 13 μM and 8 ± 0.4 mM, respectively. High concentrations of endogenous buffers define the time course of calcium transients after an action potential in the axoplasm, and may modify synaptic plasticity.
Calcium ions (Ca2+) play a leading role in the initiation, maintenance, and plasticity of neurotransmitter release at chemical synapses in both the central and peripheral nervous systems (Llinás et al., 1992; Chow et al., 1994; Augustine, 2001; Burnashev and Rozov, 2005; Schneggenburger and Neher, 2005; Pang and Südhof, 2010). The influx of Ca2+ into the nerve terminal following the action potential initiates the release of neurotransmitter and the subsequent postsynaptic potential. The properties of this potential reflect the magnitude of neurotransmitter release and the sensitivity of the postsynaptic membrane to the neurotransmitter.
In some synapses, the Ca2+ current can be measured directly (Borst and Sakmann, 1996, 1998; Yazejian et al., 1997), but such recordings from the frog neuromuscular junction are problematic because the small diameter (1–2 μm) and substantial length (100–500 μm) of the nerve terminal prohibits the use of the patch clamp and two-electrode voltage clamp techniques. Recordings with a loose macropatch clamp formed under the nerve perineurial sheath (Anderson et al., 1988), or with a macropatch electrode attached to the membrane as close as possible to a release site (Brigant and Mallart, 1982), have been accomplished only under non-physiological conditions, such as in the presence of K+-channel blockers.
Alternatively, optical methods may be employed to estimate the action-potential-induced Ca2+ influx and the subsequent distribution of Ca2+ ions in the axoplasm. The basis of these methods is the recordings of fluorescence transients which arise as a result of changes in the fluorescence intensity of specific calcium-sensitive dyes that interact with Ca2+ ions (DiGregorio and Vergara, 1997; Sabatini and Regehr, 1998; Suzuki et al., 2000; Sabatini et al., 2002; Luo et al., 2011).
The profile of intracellular Ca2+ distribution, and consequently the parameters of the fluorescence signal, is defined by the calcium influx and by the concentration of endogenous calcium-binding proteins (Ahmed and Connor, 1988; Burnashev and Rozov, 2005; Collin et al., 2005; Lin et al., 2005; Kreiner and Lee, 2006; Muller et al., 2007). After entering through voltage-gated channels, 95% of the Ca2+ ions bind to binding sites located within 10–50 nm of the site of entry (Neher, 1995). Mobile calcium buffers are represented by soluble Ca2+-binding proteins of the EF-hand family, including parvalbumins, calbindin-D9k, calbindin-D28k, calretinin, calmodulin, and proteins of the S100 family (Schwaller, 2010), or else by small molecules such as ATP and GTP which affect the diffusion of Ca2+ and define their operational distance (Feher et al., 1989; Zhou and Neher, 1993). Fixed endogenous Ca2+ buffers are associated with cellular organelles or elements of the cytoskeleton (Gabso et al., 1997; Muller et al., 2005). Fixed buffers slow the diffusion of Ca2+ and prolong the lifetime of high Ca2+ levels near the release zone (Zhou and Neher, 1993; Kits et al., 1999; Matveev et al., 2004). But for many peripheral synapses, including the frog neuromuscular junction, many of the details of the regulation of axoplasmic Ca2+ and of the properties of endogenous calcium buffers in terminals remain unknown.
In previous work, we have compared experimental and model-based approaches to show that the presence of endogenous fixed calcium buffers in the motor nerve terminal of frogs can explain the desynchronization of neurotransmitter release under conditions of low extracellular calcium and synchronization of quantal release in the presence of a mobile buffer (Samigullin et al., 2005; Bukharaeva et al., 2007; Gilmanov et al., 2008). In the present work, we measure fluorescence transients activated by a single nerve impulse in the motor nerve terminal of a frog using two calcium-sensitive fluorescent dyes with different affinities to calcium. We compare these fluorescence transients with the simultaneous measurement of the postsynaptic response of a muscular fiber. We then use mathematical modeling to estimate the parameters of the calcium current and the contribution of the endogenous calcium buffers of different mobilities to the spatial distribution of Ca2+ in the axoplasm.
Materials and Methods
Preparation and Solutions
The animals (frogs, Rana ridibunda) were anesthetized with ether prior to being decapitated in accordance with the European Communities Council Directive (November 24, 1986; 86/609/EEC). The protocol of the experiments was approved by the Animal Care and Use Committee of Kazan State Medical University. An isolated neuromuscular preparation of the musculus cutaneous pectoris was continuously perfused with Ringer solution (in mM): NaCl 113.0, KCl 2.5, CaCl2 1.8, and NaHCO3 3.0 (pH 7.2–7.4). The temperature was set at 20.0 ± 0.3°C. D-tubocurarine (5 μM) was added to the solution to prevent muscle contractions in response to nerve stimulation.
Calcium Dye Loading
Two fluorescent calcium dyes of different affinities were used: the low-affinity Magnesium Green (MG), which has dissociation constant Kd = 19 μM, and the high-affinity Oregon Green BAPTA-1 (OGB-1), whose dissociation constant is Kd = 0.17 μM (Maravall et al., 2000). The loading concentration of the fluorescent dyes was 50 mM. All dyes were obtained from Molecular Probes (Eugene, OR, USA).
The loading of the dyes was carried out through the nerve stump, as described by Peng and Zucker (1993) and Wu and Betz (1996). After this procedure, all terminals in the proximal part of the nerve trunk had sufficient basal levels of fluorescence. It has been established that the approximate intra-terminal concentration of dyes loaded into a nerve by this method is 40–150 μ M (Suzuki et al., 2000). The loading procedure does not have any appreciable influence on the amplitude of the postsynaptic response or on the frequency of the miniature end-plate potentials (Wu and Betz, 1996).
Recording of the Fluorescent Signals of the Single Action Potentials
Fluorescence signals were monitored using a photometric system based on an Olympus BX-51 microscope with water immersion objective (40×, NA = 0.8) and a high-sensitivity photodiode (S1087, Hamamatsu, Japan) (Samigullin et al., 2010). To choose the area for recording the fluorescence, we used the optical viewfinder (Till Photonics, Munich, Germany). Light excitation at a wavelength of 488 nm was generated by Polychrome V (Till Photonics, Munich, Germany). To minimize bleaching of the dye and to the lower background fluorescence, the recording area of the nerve terminal was restricted using an iris diaphragm. Illumination was gated by a shutter with a typical exposure time of 400 μs and a delivery rate of 0.5 Hz. The signals from the photodiode were digitized with a Digidata 1200B analog-to-digital converter (Axon Instruments, California, USA). The recording of fluorescence and the synchronization of the illumination and stimulation was controlled by WinWCP software (Strathclyde University, Glasgow, UK). The change in fluorescence was represented as ΔF/F0 (the change in fluorescence intensity relative to the background fluorescence as a percentage). For each experiment, about 60 fluorescent responses were recorded and averaged.
The time parameters of the fluorescence signals depend on the interaction of the dye with the calcium—specifically, on the speed of the forward and reverse reactions of the calcium–dye complex formation. The speeds of the forward reaction of the low-affinity (MG) and high-affinity (OGB-1) dyes differ (see Table 1). The dissociation rate of the calcium–dye complex is almost 30 times lower for OGB-1 than for MG, which explains the longer rise time and decay time observed with OGB-1.
Electrophysiological Recording of Nerve-Action Potentials and End-Plate Currents
Suprathreshold stimuli of 0.2 ms duration were applied to the nerve at 0.5 Hz through a suction electrode. The nerve-action potentials and extracellular end-plate currents were recorded using heat-polished Ringer-filled extracellular pipettes with tip diameters of 20 μm and resistances of 1–3 MΩ. For optimal contact with the preparation, the pipette tip was beveled at 35° (Vyshedskiy and Lin, 2000). The calcium transients were recorded in the same region where the extracellular pipette was placed. The electrical and optical signals were sampled simultaneously at 3-μs intervals by a 12-bit analog–digital converter (Digidata 1200B, Axon Instruments, California, USA), stored on a computer, and processed in AxoScope (Axon Instruments, California, USA).
Computer Modeling of Intracellular Calcium Concentrations
Experimental and theoretical studies have shown that, in the frog nerve terminal, a presynaptic action potential only opens a few calcium channels, and that the interaction of calcium influxes through different channels is low (Wachman et al., 2004; Luo et al., 2011). Experiments in non-mammalian synapses indicate that calcium channels are tightly coupled to the sites where synaptic vesicles are docked; the opening of a single Ca2+ channel (producing a Ca2+ nanodomain; Yoshikami et al., 1989; Augustine, 1990; Stanley, 1991) is sufficient to initiate vesicular fusion. We therefore neglected non-linear interactions of currents through the channels and considered only one channel in the model of one active zone, with the assumption that the total Ca2+ current is the linear summation of currents through the individual channels. Such a presentation of the Ca2+ channel as a point source on the presynaptic membrane is common for modeling exocytotic release (Smith, 1996; Bauer, 2001; Gilmanov et al., 2008).
The terminal was modeled as a cylinder with height of 2 μm and a diameter of 2 μm, with the calcium channel as a point source in its center. The height of frog's neuromuscular terminal (the height of the cylinder in the model—Figure S1 in the Supplementary Materials) was taken from Birks et al. (1960) and Pawson et al. (1998). In addition to fluorescent dye molecules in the cytoplasm, the presence of mobile and fixed buffers was also assumed (Klingauf and Neher, 1997). The interaction of Ca2+ in the terminal with the fluorescent dye and endogenous Ca2+ buffers can be described by the equation:
where konBn and kofBn are the association and dissociation rate constants and Bn (where n = d, fix, mob, referring to the fluorescent dye, fixed endogenous buffer, and mobile endogenous buffer, respectively). Assuming mass action kinetics and Fickian diffusion, we established a system of seven partial differential equations (PDE):
where DCa, DBmob, and Dd are the diffusion constants for free Ca2+, the mobile buffer, and the fluorescent dye, respectively, δ (x − x0, y − y0, z − z0) is the Dirac delta function denoting the calcium channel, where x0, y0, z0 are coordinates of the center of cylinder base. Three-dimensional geometry with Cartesian coordinates was used in our modeling.
The action-potential-induced calcium influx iCa was modeled as a Gaussian-shaped calcium current, as measured in the presynaptic terminals (Llinás et al., 1976):
where A is current amplitude and tpeak is the time of the peak of the Gaussian calcium current. In this model, we have used the maximum of the presynaptic action potential as the starting point of the calcium current, since it is precisely the depolarization caused by the action potential that opens the presynaptic calcium channels. The parameter σ defines the width of the calcium current, for which we have also calculated the parameter t1/2 (the width of the current at half-amplitude). The parameters t1/2 and σ are connected by the equation:
We have taken into account two additional calcium currents that pass through the cylinder base: the current representing the work of the membrane pumps and a leak current:
where Vmax is the maximum pump velocity, Kpm is the pump dissociation constant, Kleak is the membrane conductivity for calcium, and [Caextra] is the extracellular concentration of calcium.
For all species, the flux through all boundaries was null, i.e.,
where n is normal to the surface cylinder, A is the designation of some species in the model. However, the current J through the cylinder base containing the calcium channel was present in the model throughout all the simulation. So for this base only:
The initial conditions for all species were the equilibrium concentrations listed in Table 1.
The system of PDEs was solved using the finite element method, taking 0.1 μm as the size of elements in the volume of a cylinder and 0.1 nm near the point source, using Comsol Multiphysics 3.5 and Matlab.
The model parameters are specified in Table 1. To model the Ca2+ transients, we varied the concentrations of the mobile and fixed buffers, as well as the amplitude and duration of the inward Ca2+ current. We employed a human-based approach here, but the variations of all free parameters were made in the full physiological range in increments of 5% of the range. We then compared the resulting fluorescence signal of the model with the signal observed in experiments. The theoretical fluorescence signal was used to derive the total number of calcium–dye complexes in the model at each moment in time, as the fluorescence of these molecules gives rise to the fluorescence signal registered in the experiment. The concentration range of the mobile buffer was 50–500 μM, while the concentration range of the fixed buffer was 2000–12000 μ M. These parameters correspond to the values established for the nerve terminals in synapses of a crayfish (Vyshedskiy and Lin, 2000) and in synapses of the central nervous system (Wu and Saggau, 1994). The selection of parameters that provides the best agreement of the modeling and experimental curves with the calcium transient for the differing concentrations of fixed and mobile buffers are shown on Figure S2 in the Supplementary Materials.
Action Potential Induced Responses of Ca2+ Dyes with Differing Affinities for Ca2+
The peak of the fluorescence signal in nerve terminals loaded with the low-affinity dye MG was about 2.6 times smaller than in terminals loaded with the high-affinity dye OGB-1 (see Table 2 for details). The kinetics of the response of the low-affinity dye was significantly faster than the response of the high-affinity dye (Table 2 and Figure 1). Similar differences have been observed in synapses of the central nervous system (Regehr and Atluri, 1995) and in developing neuromuscular junctions of Xenopus (DiGregorio and Vergara, 1997). As described in the Materials and Methods section, the time parameters of fluorescence signals depend on the speed of the forward and reverse reactions of the formation of the calcium–dye complex. Thus, identical inward calcium currents evoked in nerve terminals loaded with dyes with different properties are accompanied by fluorescence signals with different amplitude and temporal properties. It is therefore impossible to accurately determine the behavior of the intracellular free calcium based on the response of only one fluorescent dye.
Figure 1. The time course of calcium transients obtained in the nerve terminal from one electrical stimulus with two different dyes. Magnesium Green (magenta line) and Oregon Green BAPTA-1 (green line). (A) Averaged unnormalized calcium transients represented as ΔF/F0. (B) Averaged calcium transients normalized to the maximum value to compare the time parameters. The averages were calculated from 60 responses. The beginning of the short (200 μs) electrical stimulus of the nerve terminal is indicated with arrows.
Simultaneous Recording of Presynaptic Fluorescence Signal and Postsynaptic Electrical Response
To reveal the correlation between the amplitude of the presynaptic calcium influx and neurotransmitter release, we simultaneously recorded the presynaptic fluorescence signal and the end-plate potential. Figure 2 shows the calcium signals and the extracellulary recorded action potentials of the nerve terminal and end-plate currents. It is clearly visible in Figure 2 that the postsynaptic current developed much faster than the fluorescence signal; in essence, the postsynaptic current ended before the fluorescence signal even reached its maximum. This means that, although the fluorescence signal occurs as a consequence of the interactions of calcium ions with the indicator in the nerve terminal of the frog neuromuscular junction, it does not directly reflect the change in Ca2+ concentration in the active release zones.
Figure 2. Simultaneous recording of end-plate currents and fluorescence calcium transients caused by a single action potential in a terminal loaded with two different dyes. (A) Fluorescence calcium transient for Oregon Green BAPTA-1(green line), end-plate current (black line); (B) fluorescence calcium transient for Magnesium Green (magenta line), end-plate current (black line). The figure shows typical experimental recordings for each dye obtained from the average of 60 responses per experiment. Traces are normalized to the maximum values. End-plate currents were inverted to allow comparison with fluorescence signals.
Because electrophysiological techniques cannot be used to directly record calcium currents in the frog terminal, we used the mathematical model described in the Materials and Methods section to deduce the parameters of the inward calcium current that reproduces the observed fluorescence transient, and also in order to define a profile of calcium concentration changes at different points of the terminal, including the near-membrane space.
Model Parameters and Quantification of the Endogenous Calcium Buffer Concentration
In order to increase the reliability of the model, we compared the theoretical and experimental fluorescence transients obtained with the different dyes and different regimes of nerve stimulation. As specified in Table 1, which shows the model parameters, the amplitude and duration of the calcium current and the concentrations of the fixed and mobile buffers were the variable parameters in the model. These were selected so that the model best reproduced the kinetics of the fluorescence transient for both dyes. We found that a presynaptic calcium current with amplitude 1.60 ± 0.08 pA and half-width 1.20 ± 0.06 ms resulted in the best agreement between the model and the experimental fluorescence transients. The presence in the axoplasm of a mobile calcium buffer at a concentration of 250 ± 13 μM and of fixed buffer at 8.0 ± 0.4 mM could explain the observed temporal differences in the measured fluorescence signal, compared with the postsynaptic response. We have checked additionally that, when all four model-derived parameters are within the intervals shown in Table 1, the shape of simulated fluorescence reproduces the shape of the experimental fluorescence (taking into account the dispersion of the experimental data; Figure 3). Figure S2 in the Supplementary Materials demonstrate the selection of parameters that provides the best agreement of the modeling and experimental curves with the calcium transient for the differing concentrations of fixed and mobile buffers.
Figure 3. Superposition of experimental and model calcium transients in terminals loaded with two different dyes. (A) Experimental curve for Oregon Green BAPTA-1(green line), model calcium transient (blue line), confidence interval for experimental curve (thin black line); (B) experimental curve for Magnesium Green (magenta line), model calcium transient (blue line), and confidence interval for experimental curve (thin black line). Traces are normalized to the maximum value.
Fluorescence Transients During Rhythmic Stimulation of the Motor Nerve
The neuromuscular junction in vivo operates in a mode of rhythmic activity, so it is important to estimate the intracellular calcium kinetics not only in response to individual stimuli, but also during repetitive stimulation. We estimated the properties of the calcium transients for the three signals occurring in response to a 50-Hz train of action potentials in the presence of dyes of both types.
The parameters of the simulated fluorescence response were in a good agreement with those of the experimental response obtained using the high- and low-affinity dyes (Figure 4), thus confirming the validity of the model.
Figure 4. Superposition of experimental and model calcium transients in the terminal loaded with two different dyes during a stimulus train consisting of three signals at 50 Hz. (A) Experimental curve for Oregon Green BAPTA-1(green line), model calcium transient (thick blue line), model inward calcium current (thin blue line); (B) experimental curve for Magnesium Green (magenta line), model calcium transient (thick blue line), model inward calcium current (thin blue line). Traces are normalized to the maximum value.
The Estimation of Calcium Concentration in the Near-Membrane Space
Our model allowed us to calculate the calcium distribution in different parts of the terminal. It is particularly important to estimate Ca2+ concentration in the near-membrane region, since this determines the number of released quanta and the kinetics of neurotransmitter release. Figure 5 shows the calculated dynamics of the calcium concentration at a distance of 15 nm from the calcium channel, in comparison with the experimental postsynaptic current and fluorescent signal. The beginning of the rise in the calcium concentration in the near-membrane space preceded the postsynaptic response, indicating that the real time between the beginning of the calcium response in the active zone and the beginning of the postsynaptic current (about 600 μs) was sufficient for the activation of the exocytosis machinery, the release of neurotransmitter, the activation of the postsynaptic receptors, and the opening of the postsynaptic ion channels.
Figure 5. Superposition of the calculated calcium concentration dynamics at a distance of 15 nm from the calcium channel, with simultaneous recording of end-plate current and calcium transients induced by a single action potential in terminals loaded with two different dyes. (A) Fluorescence calcium transient for Oregon Green BAPTA-1(green line), experimental end-plate current (black line), model calcium concentration (blue line); (B) fluorescence calcium transient for Magnesium Green (magenta line), experimental end-plate current (black line), model calcium concentration (blue line). Traces are normalized to the maximum value. Experimental data are the same as in Figure 2.
In a nerve terminal loaded with calcium-sensitive dye, calcium ions entering through Ca2+ channels bind to the dye, producing a change in its fluorescence intensity. Analysis of Ca2+-related fluorescent signals reveals the dynamics of calcium concentration in the nerve terminal which demonstrate, for example, that synaptic facilitation caused by the accumulation of intracellular calcium outlasts the increase in Ca2+ in the terminal (Regehr and Atluri, 1995). Analysis of Ca2+-dependent fluorescence changes has also been used to identify the types of voltage-gated calcium channels that participate in neurotransmitter release, and to clarify the degree of cooperativity of Ca2+ channels (Wu and Saggau, 1994; Mintz et al., 1995; Tsentsevitsky et al., 2014). While electrophysiological recording of Ca2+ currents provides a more direct measure (Vyshedskiy and Lin, 2000), the morphological features of some synapses—such as the neuromuscular junctions of mice, frogs, and other vertebrates—render the use of voltage clamp methods difficult. Therefore, the use of calcium-sensitive dyes is an alternative for estimating calcium influx into the nerve terminal.
The fluorescent signal had a very slow rise compared with the postsynaptic response, and presumably the calcium current. Similar observations have been made in the neuromuscular synapses of lizards and frogs (David et al., 1997; DiGregorio and Vergara, 1997; Bélair et al., 2005). However, this temporal difference between the beginning of the fluorescence signal and the postsynaptic potential has not always been observed, and it has generally been assumed that the fluorescence transient reflects the change in calcium concentration in the terminal, while its first derivative reflects the calcium current (Sabatini and Regehr, 1998). The lag observed in our study could be interpreted, for example, as a consequence of the presence of a high concentration of fixed buffer, resulting in slow diffusion of calcium into the terminal (Lin et al., 2005). In this case, calcium entering the terminal interacts directly with proteins of the exocytosis machinery to trigger neurotransmitter release, and only subsequently does the slow diffusion of calcium cause the formation of calcium–dye complexes and an appreciable fluorescence signal. Because we recorded the fluorescence from the entire volume of the terminal, the slow kinetics of the signal are likely to reflect a process of slow calcium distribution. Another reason why the postsynaptic current in Figure 2 develops earlier than the fluorescence signal is that the frog neuromuscular junction is a fast synapse (as opposed to slow NMDA, for example) with relatively short synaptic delays and quickly rising and falling postsynaptic currents. For the frog neuromuscular junction, the time between neurotransmitter release and the start of the postsynaptic current is very short (about 300–500 μs; Bukharaeva et al., 2002). In some other synapses, the slow activation of postsynaptic channels can mask the paradox shown in Figure 2. The current through these channels grows at a similar speed to the fluorescence signal (for example, Vyshedskiy and Lin, 2000) due to the slow activation.
Our mathematical model permits us to clarify the observation that the postsynaptic current develops much faster than the fluorescence signal. After the opening of the voltage-gated channels, Ca2+ enters the terminal and, in the process of diffusing away from the entry site, binds to the dye molecules that are distributed throughout the volume of the terminal. The rising phase of the fluorescent signal of the calcium–dye complex in the area recorded by the photodiode depends on the speed of calcium diffusion, which in turn depends on the concentration and properties of endogenous buffers. Endogenous calcium buffers play the leading role in the spatiotemporal distribution of Ca2+ (Sala and Hernandez-Cruz, 1990; Nowycky and Pinter, 1993). After the action-potential-triggered Ca2+ influx through the voltage-dependent channels, Ca2+ ions bind to the buffers within 10–50 nm of the entry site (Augustine and Neher, 1992; Neher, 1995). Further diffusion of Ca2+ depends mainly on the degree of mobility of Ca2+-binding sites. Immobile buffers, which are localized on cellular organelles or cytoskeleton elements, have a slower speed of diffusion than free calcium. Such fixed buffers delay diffusion, because they take up and bind calcium ions for some time, thus reducing the speed of calcium distribution within the space of the terminal (Zhou and Neher, 1993; Kits et al., 1999). This slowing of diffusion has important consequences for any intracellular calcium-dependent process.
As Ca2+ slowly diffuses through the terminal, dissociation of the calcium–dye complexes begins, and the duration of the rising phase of the fluorescent signal also depends on the speed of this process. Figure 6 shows the changes in the concentration of a calcium–dye complex at various points of a terminal for the two different dyes. Dissociation of the calcium–dye complex is expected to occur faster in the case of a low-affinity dye, as the rate of the reverse reaction is greater than for the high-affinity dye. The phase of the signal rise ends when the rate of formation of calcium–dye complexes becomes slower than the rate of dissociation. For the low-affinity dye, dissociation occurs rapidly, and thus the rising phase of the signal is shorter. Conversely, the slower speed of dissociation of the high-affinity dye leads to a longer rising phase.
Figure 6. Change in concentration of a calcium–dye complex at different distances from the point of calcium entry (15, 50, 100, 200, 300 nm) in the terminal loaded with Oregon Green BAPTA-1 (A) or Magnesium Green (B).
Our mathematical model shows that the kinetics of Ca2+ transients observed in the presence of dyes with different calcium affinity and different calcium–dye complex dissociation rates could be described using a fixed-buffer concentration of 8.0 ± 0.4 mM and a mobile-buffer concentration of 250 ± 13 μM. Similar results (7.5 mM for slow buffers) have been reported by Bennett et al. (2007), who compared experimental and model-derived presynaptic fluorescence signals during rhythmic stimulation of toad motor nerves, considering the mitochondria and endoplasmic reticulum as fixed buffers. Suzuki et al. (2000) have assumed the existence in the frog nerve terminal cytoplasm of a fixed low-affinity Ca2+ buffer at 10 mM concentration. Such fixed buffers can determine the early decay of Ca2+ transients after an action potential in crayfish neuromuscular terminals (Lin et al., 2005). Parameters determined from our modeling thus coincide with the values proposed by other authors (Aharon et al., 1994; Suzuki et al., 2000; Bennett et al., 2007).
The idea that the inclusion of fluorescent dyes itself increases the local Ca2+ buffering, thus altering the Ca2+ dynamics, is often encountered in the literature (e.g., Sabatini et al., 2002). But in our model the interaction of the dye with calcium and diffusion of the dye and the dye-calcium complex in the terminal are simulated separately (Equation 2). Thus, our estimates of buffers concentration are made already taking into account the dye influence on the calcium diffusion. In principle it was shown experimentally (Wu and Betz, 1996) that the calcium dye loading haven't any appreciable influence on the calcium dependent quantal release (parameters of the postsynaptic response or of the frequency of the miniature end-plate potentials do not change). It suggests that endogenous calcium buffers have the buffering capacity much stronger than fluorescent dyes. But our model allows even to check the changes in intracellular calcium distribution at fluorescent dye loading.
The presence of mobile and fixed buffers in the axoplasm explains the time discrepancy between the fluorescent signal and the postsynaptic response (Figures 2, 5). Estimation of the inward calcium current shows that the best agreement between model and experimental fluorescence transients was achieved at a current amplitude of 1.6 ± 0.08 pA and a half-width of 1.2 ± 0.06 ms (Figures 3, 4). Assuming that the conductance of presynaptic N-type calcium channels (which are dominant in frog nerve terminals) is 0.33 pA (Weber et al., 2010), we may conclude that about five channels are activated by a single action potential. These data correspond to the results of Shahrezaei et al. (2006), who showed that from two to six channels open per active zone per action potential; and that fusion of a single vesicle is triggered by the opening of one or two channels.
In summary, the analysis of fluorescence transients in the nerve terminal of the frog neuromuscular junction and the comparison of the kinetics of experimental and model responses allowed us to estimate the intracellular concentrations of fixed and mobile endogenous buffers, the properties of the calcium current, and the kinetics of the changes in calcium concentrations in the nerve terminal.
The high concentration of intracellular calcium buffers, which slow diffusion by binding calcium to cytoplasmic elements, has important consequences for any intracellular Ca2+ -dependent process. Diffusion cannot be ignored if the time course of the process studied is on the order of milliseconds, and may contribute to delayed asynchronous release of neurotransmitters (Rahamimoff and Yaari, 1973), calcium-induced calcium release from intracellular calcium stores (Akita and Kuba, 2008), facilitation and augmentation (Miller, 1991), and modulation of Ca2+-dependent deactivation of Ca2+ channels (Kreiner and Lee, 2006). High concentrations of endogenous buffers modify transmitter release at synapses (Caillard et al., 2000; Matveev et al., 2004).
In conclusion, the concentration of a calcium–dye complex is proportional to the concentration of Ca2+ at each point of the terminal volume. When studying small terminals with fast diffusion and a short time to equilibrium, the concentration of Ca2+ can be accurately estimated by the first derivative of the smoothed fluorescent signal (Sabatini and Regehr, 1998). Because of the slow “delayed” diffusion and the resulting longer time before equilibrium in the frog motor nerve terminal, the fluorescence signal does not correspond directly to Ca2+ concentration in the terminal, particularly near the active zone. However, mathematical modeling allows determination of Ca2+ concentration and mobile and fixed buffers concentrations near the membrane of the nerve terminal where the calcium-sensitive exocytosis proteins are placed.
This research has been performed under the Russian Government's Program for Competitive Growth of Kazan Federal University and was also supported by a grant of President of Russian Federation “Leading Scientific School” and a grants from the Russian Foundation for Basic Research (13-04-00886; 14-04-00790; 15-04-02983).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are grateful to Prof. Lev Magazanik, Prof. Alexei Verkhratsky, and Dr. Najl Valeev for their many helpful discussions.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fnsyn.2014.00029/abstract
Anderson, A. J., Harvey, A. L., Rowan, E. G., and Strong, P. N. (1988). Effects of charybdotoxin, a blocker of Ca2+-activated K+ channels, on motor nerve terminals. Br. J. Pharmacol. 95, 1329–1335. doi: 10.1111/j.1476-5381.1988.tb11772.x
Bélair, E. L., Vallée, J., and Robitaille, R. (2005). Long-term in vivo modulation of synaptic efficacy at the neuromuscular junction of Rana pipiens frogs. J. Physiol.569, 163–178. doi: 10.1113/jphysiol.2005.094805
Bennett, M. R., Farnell, L., Gibson, W. G., and Dickens, P. (2007). Mechanisms of calcium sequestration during facilitation at active zones of an amphibian neuromuscular junction. J. Theor. Biol. 247, 230–241. doi: 10.1016/j.jtbi.2007.03.022
Borst, J. G., and Sakmann, B. (1998). Calcium current during a single action potential in a large presynaptic terminal of the rat brainstem. J. Physiol. 506, 143–157. doi: 10.1111/j.1469-7793.1998.143bx.x
Bukharaeva, E. A., Samigullin, D. V., Nikolsky, E. E., and Magazanik, L. G. (2007). Modulation of the kinetics of evoked quantal release at mouse neuromuscular junctions by calcium and strontium. J. Neurochem. 100, 939–949. doi: 10.1111/j.1471-4159.2006.04282.x
Bukharaeva, E. A., Samigullin, D. V., Nikolsky, E. E., and Vyskocil, F. (2002). Protein kinase A cascade regulates quantal release dispersion at frog muscle endplate. J. Physiol. 538, 837–848. doi: 10.1113/jphysiol.2001.012752
Caillard, O., Moreno, H., Schwaller, B., Llano, I., Celio, M. R., and Marty, A. (2000). Role of the calcium-binding protein parvalbumin in short-term synaptic plasticity. Proc. Natl. Acad. Sci. U.S.A. 97, 13372–13377. doi: 10.1073/pnas.230362997
Chow, R., Klingauf, J., and Neher, E. (1994). Time course of Ca2+ concentration triggering exocytosis in neuroendocrine cells. Proc. Natl. Acad. Sci. U.S.A. 91, 12765–12769. doi: 10.1073/pnas.91.26.12765
Collin, T., Chat, M., Lucas, M. G., Moreno, H., Racay, P., Schwaller, B., et al. (2005). Developmental changes in parvalbumin regulate presynaptic Ca2+ signaling. J. Neurosci. 25, 96–107. doi: 10.1523/JNEUROSCI.3748-04.2005
DiGregorio, D. A., and Vergara, J. L. (1997). Localized detection of action potential-induced presynaptic calcium transients at a Xenopus neuromuscular junction. J. Physiol. 505, 585–592. doi: 10.1111/j.1469-7793.1997.585ba.x
Feher, J. J., Fulmer, C. S., and Fritzsch, G. K. (1989). Comparison of the enhanced steady-state diffusion of calcium by calbindin-D9K and calmodulin, possible importance in intestinal calcium absorption. Cell Calcium 10, 189–203. doi: 10.1016/0143-4160(89)90002-X
Gilmanov, I. R., Samigullin, D. V., Vyskocil, F., Nikolsky, E. E., and Bukharaeva, E. A. (2008). Modeling of quantal neurotransmitter release kinetics in the presence of fixed and mobile calcium buffers. J. Comput. Neurosci. 25, 296–307. doi: 10.1007/s10827-008-0079-5
Helmchen, F., Imoto, K., and Sakmann, B. (1996). Ca2+ buffering and action potential-evoked Ca2+ signaling in dendrites of pyramidal neurons. Biophys. J. 70, 1069–1081. doi: 10.1016/S0006-3495(96)79653-4
Kits, K. S., de Vlieger, T. A., Kooi, B. W., and Mansvelder, H. D. (1999). Diffusion barriers limit the effect of mobile calcium buffers on exocytosis of large dense cored vesicles. Biophys. J. 76, 1693–1705. doi: 10.1016/S0006-3495(99)77328-5
Kreiner, L., and Lee, A. (2006). Endogenous and exogenous Ca2+ buffers differentially modulate Ca2+-dependent inactivation of Cav2.1 Ca2+channels. J. Biol. Chem. 281, 4691–4698. doi: 10.1074/jbc.M511971200
Lin, J. W., Fu, Q., and Allana, T. (2005). Probing the endogenous Ca2+ buffers at the presynaptic terminals of the crayfish neuromuscular junction. J. Neurophysiol. 94, 377–386. doi: 10.1152/jn.00617.2004
Llinás, R. I., Steinberg, Z., and Walton, K. (1976). Presynaptic calcium currents and their relation to synaptic transmission: voltage clamp study in squid giant synapse and theoretical model for the calcium gate. Proc. Natl. Acad. Sci. U.S.A. 73, 2918–2922. doi: 10.1073/pnas.73.8.2918
Luo, F., Dittrich, M., Stiles, J. R., and Meriney, S. D. (2011). Single-pixel optical fluctuation analysis of calcium channel function in active zones of motor nerve terminals. J. Neurosci. 31, 11268–11281. doi: 10.1523/JNEUROSCI.1394-11.2011
Maravall, M., Mainen, Z. F., Sabatini, B. L., and Svoboda, K. (2000). Estimating intracellular calcium concentrations and buffering without wavelength ratioing. Biophys. J. 78, 2655–2667. doi: 10.1016/S0006-3495(00)76809-3
Muller, A., Kukley, M., Stausberg, P., Beck, H., Muller, W., and Dietrich, D. (2005). Endogenous Ca2+ buffer concentration and Ca2+ microdomains in hippocampal neurons. J. Neurosci. 25, 558–565. doi: 10.1523/JNEUROSCI.3799-04.2005
Muller, M., Felmy, F., Schwaller, B., and Schneggenburger, R. (2007). Parvalbumin is a mobile presynaptic Ca2+ buffer in the Calyx of Held that accelerates the decay of Ca2+ and short-term facilitation. J. Neurosci. 27, 2261–2271. doi: 10.1523/JNEUROSCI.5582-06.2007
Pawson, P. A., Grinnell, A. D., and Wolowske, B. (1998). Quantitative freeze-fracture analysis of the frog neuromuscular junction synapse I. Naturally occurring variability in active zone structure. J. Neurocytol. 27, 361–377. doi: 10.1023/A:1006942909544
Peng, Y. Y., and Zucker, R. S. (1993). Release of LHRH is linearly related to the time integral of presynaptic Ca2+ elevation above a threshold level in bullfrog sympathetic ganglia. Neuron 10, 465–473. doi: 10.1016/0896-6273(93)90334-N
Samigullin, D. V., Bukharaeva, E. A., Vyskočil, F., and Nikolsky, E. E. (2005). Calcium dependence of uni-quantal release latencies and quantal content at mouse neuromuscular junction. Physiol. Res. 54, 129–132.
Samigullin, D. V., Vasin, A. L., Bukharaeva, E. A., and Nikolsky, E. E. (2010). Characteristics of calcium transient in different parts of frog nerve terminal in response to nerve impulse. Dokl. Biol. Sci. 431, 83–85. doi: 10.1134/S0012496610020043
Schmidt, H., Brown, E. B., Schwaller, B., and Eilers, J. (2003). Diffusional mobility of parvalbumin in spiny dendrites of cerebellar Purkinje neurons quantified by fluorescence recovery after photobleaching. Biophys. J. 84, 2599–2608. doi: 10.1016/S0006-3495(03)75065-6
Shahrezaei, V., Cao, A., and Delaney, K. R. (2006). Ca2+from one or two channels controls fusion of a single vesicle at the frog neuromuscular junction. J Neurosci. 26, 13240–132499. doi: 10.1523/JNEUROSCI.1418-06.2006
Suzuki, S., Osanai, M., Murase, M., Suzuki, N., Ito, K., Shirasaki, T. K., et al. (2000). Ca2+ dynamics at the frog motor nerve terminal. Pflug. Arch. Eur. J. Phisiol. 440, 351–365. doi: 10.1007/s004240000278
Tsentsevitsky, A. N., Samigullin, D. V., Nurullin, L. F., Khaziev, E. F., Nikolsky, E. E., and Bukharaeva, E. A. (2014). “Presynaptic voltage-dependent calcium channels at the frog neuromuscular junction,” in Frogs: Genetic diversity, Neural Development and Ecological Implication, ed H. Lambert (New York, NY: Nova Science Publishers), 179–195.
Wachman, E. S., Poage, R. E., Stiles, J. R., Farkas, D. L., and Meriney, S. D. (2004). Spatial distribution of calcium entry evoked by single action potentials within the presynaptic active zone. J. Neurosci. 24, 2877–2885. doi: 10.1523/JNEUROSCI.1660-03.2004
Weber, A., Wong, F., Tufford, A., Schlichter, L., Matveev, V., and Stanley, E. (2010). N-type Ca2+ channels carry the largest current: implications for nanodomains and transmitter release. Nature 13, 1348–1350. doi: 10.1038/nn.2657
Wu, L. G., and Betz, W. J. (1996). Nerve activity but not intracellular calcium determines the time course of endocytosis at the frog neuromuscular junction. Neuron 17, 769–779. doi: 10.1016/S0896-6273(00)80208-1
Wu, L. G., and Saggau, P. (1994). Presynaptic calcium is increased during normal synaptic transmission and paired-pulse facilitation, but not in long-term potentiation in area CA1 of hippocampus. J. Neurosci. 14, 645–654.
Xu, T., Naraghi, M., Kang, H., and Neher, E. (1997). Kinetic studies of Ca2+ binding and Ca2+ clearance in the cytosol of adrenal chromaffin cells. Biophys. J. 73, 532–545. doi: 10.1016/S0006-3495(97)78091-3
Yazejian, B., DiGregorio, D. A., Vergara, J. L., Poage, R. E., Meriney, S. D., and Grinnell, A. D. (1997). Direct measurements of presynaptic calcium and calcium-activated potassium currents regulating neurotransmitter release at cultured Xenopus nerve-muscle synapses. J. Neurosci. 17, 2990–3001.
Yoshikami, D., Bagabaldo, Z., and Olivera, B. M. (1989). The inhibitory effects of omega-conotoxins on Ca channels and synapses. Ann. NY Acad. Sci. 560, 230–248. doi: 10.1111/j.1749-6632.1989.tb24100.x
Keywords: neuromuscular junction, presynaptic calcium transient, inward calcium current, endogenous calcium buffers, mathematical modeling
Citation: Samigullin D, Fatikhov N, Khaziev E, Skorinkin A, Nikolsky E and Bukharaeva E (2015) Estimation of presynaptic calcium currents and endogenous calcium buffers at the frog neuromuscular junction with two different calcium fluorescent dyes. Front. Synaptic Neurosci. 6:29. doi: 10.3389/fnsyn.2014.00029
Received: 06 October 2014; Accepted: 12 December 2014;
Published online: 07 January 2015.
Edited by:Lisa Topolnik, Laval University, Canada
Copyright © 2015 Samigullin, Fatikhov, Khaziev, Skorinkin, Nikolsky and Bukharaeva. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dmitry Samigullin, Laboratory of the Biophysics of Synaptic Processes, Kazan Institute of Biochemistry and Biophysics, Russian Academy of Sciences, PO Box 30, 2/31, Lobachevskogo str., Kazan 420111, Russia e-mail: email@example.com