 Research
 Open access
 Published:
Relating a calcium indicator signal to the unperturbed calcium concentration timecourse
Theoretical Biology and Medical Modelling volumeÂ 4, ArticleÂ number:Â 7 (2007)
Abstract
Background
Optical indicators of cytosolic calcium levels have become important experimental tools in systems and cellular neuroscience. Indicators are known to interfere with intracellular calcium levels by acting as additional buffers, and this may strongly alter the timecourse of various dynamical variables to be measured.
Results
By investigating the underlying reaction kinetics, we show that in some ranges of kinetic parameters one can explicitly link the time dependent indicator signal to the timecourse of the calcium influx, and thus, to the unperturbed calcium level had there been no indicator in the cell.
Background
The use of a fluorescent calcium indicator is a familiar technique for detecting dynamical changes in intracellular calcium levels [1â€“8]. However, introduction of the indicator into the cytosol inevitably perturbs the timecourse of free cytosolic calcium by acting as a buffer, thus altering the quantity to be measured. To address this, traditional approaches to quantifying free cytosolic calcium have often restricted the use of the indicators to minimal concentrations with minimal affinity. While this minimizes the perturbation of the free calcium signal, it leads to the problem of small signaltonoise ratios.
As an alternative approach, we examine here the dynamical equations for this process in various parameter ranges in order to identify the conditions under which approximate solutions can be obtained, allowing calcium influx to be calculated directly from the fluorescence time course measurements. Knowing the calcium influx, the free cytosolic calcium can then be calculated as if there had been no indicator in the cytosol.
In the following, we will denote the temporal derivative of a variable dx(t)/dt by the symbol x'(t). Furthermore, we will use the following symbols with the units shown in Table 1.
Results
Upon activation of a neuron, calcium influx Î±(t) leads to an increase of the cytosolic calcium concentration. Once inside the cell, calcium goes one of two ways: either it is cleared from the cell, in proportion to its concentration, at a rate Î³, or it binds to an indicator with forward and backward binding rates k_{f} and k_{b}, respectively. Processes such as diffusion, internal buffering, and release from internal calcium stores or extrusion by calciumsodium exchangers are not considered here. We call the free calcium concentration x(t) with an intracellular calcium indicator, and the concentration of indicator with bound calcium y(t). We call the total indicator concentration (free and calcium bound) y_{max}. Thus, the free indicator concentration becomes y_{max}y(t). The concentration y_{max} is the initial level of the free indicator dye immediately after injection in the case of synthetic dyes, or the total amount of indicator protein (calcium bound and free) in the case of genetically encoded indicators.
The following system of two coupled nonlinear ordinary differential equations describes the dynamics of the system:
(1) x'(t) = Î±(t)  Î³Â·x(t)  y'(t)
(2) y'(t) = k_{ f }Â·x(t)Â·[y_{max}  y(t)]  k_{ b }Â·y(t)
The first equation state that the rate of change of free calcium, x'(t), is driven by the calcium flux Î±(t) and depleted by the pump Î³x(t) as well as the rate of change of indicator bound to calcium, y'(t). The second equation states that calcium is bound to the indicator at a rate k_{f} and is proportional to the concentration of free calcium, x(t), as well as to the concentration of the free indicator, y_{max}  y(t). Calcium disassociates from the indicator at a rate k_{b} and this dissociation process is proportional to the concentration of calcium bound indicator y(t).
For a constant calcium influx Î±(t) = Î±_{C}, the steadystate solutions are
(3) x_{âˆž} = Î±_{ C }/Î³, and
(4) {y}_{\mathrm{\xe2\u02c6\u017e}}={y}_{\mathrm{max}}\xe2\u2039\dots \frac{{\mathrm{\xce\pm}}_{C}}{{\mathrm{\xce\pm}}_{C}+{K}_{D}\xe2\u2039\dots \mathrm{\xce\xb3}}, or in terms of x: {y}_{\mathrm{\xe2\u02c6\u017e}}={y}_{\mathrm{max}}\xe2\u2039\dots \frac{{x}_{\mathrm{\xe2\u02c6\u017e}}}{{x}_{\mathrm{\xe2\u02c6\u017e}}+{K}_{D}}.
In general, eqs. (1) and (2) can only be solved numerically. However, if the indicator concentration is negligible compared to the calcium concentrations, eq. (1) turns into a simple differential equation describing a 1^{st} order lowpass filter with timeconstant 1/Î³:
(5) x'(t) = Î±(t)  Î³Â·x(t)
In other words: if there is no indicator present, and the pump rate is known, the unperturbed calcium concentration can be calculated as the lowpass filtered response to the calcium influx Î±(t).
Our approach will be to use eqs. (1) and (2) to determine the calcium influx Î±(t) in the presence of the indicator. This tells us how much calcium flows into the neuron as a result of activation, and allows us to remove the action of the indicator mathematically. With Î±(t) known, we may use eq. (5) to determine the timecourse of the unperturbed calcium concentration. As we will show in the following, this approach is feasible only within certain parameter regimes, but is not restricted to the linear regime. Nevertheless, we will start our considerations with an analysis of the linear regime.
The linear regime
To investigate the linear regime, we rewrite eq. (2) as
(6) y'(t)/k_{ f }= y_{max}Â·x(t)  y(t)Â·[x(t)  K_{ D }];
When x(t) is much smaller than the K_{ D }value of the indicator, eq. (6) becomes:
(7) y'(t) = y_{max}Â·k_{ f }Â·x(t)  k_{ b }Â·y(t)
Combining the derivative of this with eq. (1) gives us
(8) y''(t) + y'(t)Â·(k_{ b }+ Î³ + k_{ f }y_{max}) + y(t)Â·Î³Â·k_{ b } Î±(t)Â·k_{ f }Â·y_{max} = 0
This is a linear ordinary differential equation with constant coefficients. The solution of the homogeneous equation is of the form y(t) = cÂ·e^{Î»Â·t}, where Î» satisfies the characteristic equation:
(9) Î»^{2} + Î»A + Î³Â·k_{ b }= 0; A = k_{ b }+ Î³ + k_{ f }y_{max}.
This has solutions Î»_{1,2} with the negative inverses Ï„_{1,2} = 1/Î»_{1,2}, which are timeconstants given by
\begin{array}{cc}\left(10\right)& {\mathrm{\xcf\u201e}}_{1,2}=\frac{1}{2\mathrm{\xce\xb3}{k}_{b}}\left[{k}_{b}+\mathrm{\xce\xb3}+{k}_{f}{y}_{\mathrm{max}}\xc2\pm \sqrt{{\left({k}_{b}+\mathrm{\xce\xb3}+{k}_{f}{y}_{\mathrm{max}}\right)}^{2}\xe2\u02c6\u20194\mathrm{\xce\xb3}{k}_{b}}\right]\end{array}
Since (k_{ b }+ Î³ + k_{ f }y_{max})^{2}  4Î³Â·k_{ b }â‰¥ 0 and {k}_{b}+\mathrm{\xce\xb3}+{k}_{f}{y}_{\mathrm{max}}\xe2\u2030\yen \sqrt{{\left({k}_{b}+\mathrm{\xce\xb3}+{k}_{f}{y}_{\mathrm{max}}\right)}^{2}\xe2\u02c6\u20194\mathrm{\xce\xb3}\xe2\u2039\dots {k}_{b}}, both timeconstants are always real and positive. For small values of k_{f} y_{max}, as well as for large values of k_{b}, these become:
\begin{array}{cc}\left(11\right)& \underset{{k}_{f}{y}_{\mathrm{max}}\xe2\u2020\u20190}{\mathrm{lim}}{\mathrm{\xcf\u201e}}_{1/2}=\underset{{k}_{b}\xe2\u2020\u2019\mathrm{\xe2\u02c6\u017e}}{\mathrm{lim}}{\mathrm{\xcf\u201e}}_{1/2}=1/\mathrm{\xce\xb3}.\end{array}
The dependence of the timeconstant with the larger absolute value, Ï„_{1}, on the dimensionless parameters R_{b} and R_{f},
\begin{array}{cc}\left(12\right)& {\mathrm{\xcf\u201e}}_{1}=\frac{1}{2\mathrm{\xce\xb3}{R}_{b}}\left[{R}_{b}+{R}_{f}+1+\sqrt{{\left({R}_{b}+{R}_{f}+1\right)}^{2}\xe2\u02c6\u20194{R}_{b}}\right]\end{array}
is shown in Figure 1a. The axes are logarithmic. As one can see, the larger we set R_{b}, at fixed R_{f}, the smaller is Ï„_{1}, that is, the faster calcium is released from the bound indicator. Ï„_{1} is larger, at fixed R_{b}, for larger R_{f}; in other words, the faster calcium is bound to the indicator (k_{f}) and the larger the initial indicator concentration (y_{max}).
For the case of a pulse of injected calcium current of sufficient length, we can obtain particular solutions for y(t). For that, we insert y(t)={c}_{1}\xe2\u2039\dots {e}^{{\mathrm{\xce\xbb}}_{1}\xe2\u2039\dots t}+{c}_{2}\xe2\u2039\dots {e}^{{\mathrm{\xce\xbb}}_{2}\xe2\u2039\dots t}+k into eq. (8) and note the following initial conditions: y(0) = 0 and y'(0) = 0 for the rise of y(t) after the pulse is initiated, and y(0) = y_{max}Î±_{ c }/(Î³Â·K_{ D }) and y'(0) = 0 for the decay phase after the pulse is completed.
The initial increase of bound indicator from y(t) = 0, using A = k_{ b }+ Î³ + k_{ f }Â·y_{max} again, is
\begin{array}{cc}\left(13\right)& y(t)=\frac{{y}_{\mathrm{max}}\xe2\u2039\dots {\mathrm{\xce\pm}}_{C}}{2\xe2\u2039\dots \mathrm{\xce\xb3}\xe2\u2039\dots {K}_{D}}\xe2\u2039\dots \left[1\xe2\u02c6\u2019\left(1+\frac{A}{\sqrt{{A}^{2}\xe2\u02c6\u20194{k}_{b}\mathrm{\xce\xb3}}}\right)\xe2\u2039\dots \mathrm{exp}\left(\xe2\u02c6\u2019t/{\mathrm{\xcf\u201e}}_{1}\right)\xe2\u02c6\u2019\left(1\xe2\u02c6\u2019\frac{A}{\sqrt{{A}^{2}\xe2\u02c6\u20194{k}_{b}\mathrm{\xce\xb3}}}\right)\xe2\u2039\dots \mathrm{exp}\left(\xe2\u02c6\u2019t/{\mathrm{\xcf\u201e}}_{2}\right)\right]\end{array}
After calcium influx has stopped, when Î±(t) = 0, the bound indicator decays to zero as
\begin{array}{cc}\left(14\right)& y(t)=\frac{{y}_{\mathrm{max}}\xe2\u2039\dots {\mathrm{\xce\pm}}_{C}}{2\xe2\u2039\dots \mathrm{\xce\xb3}\xe2\u2039\dots {K}_{D}}\xe2\u2039\dots \left[\left(1+\frac{A}{\sqrt{{A}^{2}\xe2\u02c6\u20194{k}_{b}\mathrm{\xce\xb3}}}\right)\xe2\u2039\dots \mathrm{exp}\left(\xe2\u02c6\u2019t/{\mathrm{\xcf\u201e}}_{1}\right)+\left(1\xe2\u02c6\u2019\frac{A}{\sqrt{{A}^{2}\xe2\u02c6\u20194{k}_{b}\mathrm{\xce\xb3}}}\right)\xe2\u2039\dots \mathrm{exp}\left(\xe2\u02c6\u2019t/{\mathrm{\xcf\u201e}}_{2}\right)\right]\end{array}
The time course of the indicator signal under these conditions is shown in Fig. 1b.
The goal of this paper is to use the dynamical equations to determine the Î±(t) associated with an observed indicator signal y(t), and then relate that to the free calcium concentration that would be associated with this Î±(t) when the indicator is absent. In the linear regime under consideration, we need to solve eq. (8) for Î±(t) to obtain
\begin{array}{cc}\left(15\right)& \mathrm{\xce\pm}(t)=\frac{1}{{k}_{f}{y}_{\mathrm{max}}}\left[{y}^{\xe2\u20ac\xb3}(t)+{y}^{\xe2\u20ac\xb2}(t)\left({k}_{b}+\mathrm{\xce\xb3}+{k}_{f}{y}_{\mathrm{max}}\right)+y(t)\xe2\u2039\dots \mathrm{\xce\xb3}\xe2\u2039\dots {k}_{b}\right]\end{array}
From this Î±(t) the unperturbed timecourse of the calcium concentration x*(t) can be calculated from (1). It is the response of a 1^{st} order lowpass filter with timeconstant 1/Î³ to the driving input Î±(t):
\begin{array}{cc}\left(16\right)& {x}^{\xe2\u02c6\u2014}(t)={\displaystyle \underset{0}{\overset{t}{\xe2\u02c6\xab}}d{t}^{\xe2\u20ac\xb2}\mathrm{\xce\pm}(t\xe2\u02c6\u2019{t}^{\xe2\u20ac\xb2}){e}^{\xe2\u02c6\u2019\mathrm{\xce\xb3}\xe2\u2039\dots {t}^{\xe2\u20ac\xb2}}}\end{array}
From eq. (15) it also follows that \underset{{y}_{\mathrm{max}}\xe2\u2020\u2019\mathrm{\xe2\u02c6\u017e}}{\mathrm{lim}}\mathrm{\xce\pm}(t)={y}^{\xe2\u20ac\xb2}(t). This is expected as infinite indicator promptly binds all the available free calcium. So when the cell is overloaded, the indicator signal directly integrates the calcium influx: the influx can conversely be recovered by simply differentiating the indicator signal. This completes our discussion of the linear regime and we turn to the nonlinear equations again.
Approximate solution in the nonlinear regime
If we examine the nonlinear eqs. (1) and (2) we see that an approximate solution with small rate of change in the calcium bound to the indicator y'(t) is given by
\begin{array}{cc}\left(17\right)& y(t)=\frac{{y}_{\mathrm{max}}\xe2\u2039\dots x(t)}{x(t)+{K}_{D}}.\end{array}
This is an exact solution when Î±(t) is constant, and x(t) and y(t) are at the fixed point discussed earlier. So, this might well be a good guess for an approximate solution of the overall equations. We discuss this in the appendix, and argue that as long as x(t) is bounded, perturbations to this solution decay back to it at a rate to be established there. Also, the variations in x(t) are required to be slow compared to the variations in the perturbations. This means the frequency of the low pass filter giving x(t) from the calcium flux should be smaller than the decay frequencies of the perturbation. The time constant for the lowpass filter is 1/Î³.
If we use this solution, i.e. eq. (17), we have
\begin{array}{cccc}\left(18\right)& x(t)=\frac{y(t){K}_{D}}{{y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)}& \text{and}& {x}^{\xe2\u20ac\xb2}(t)=\frac{{y}^{\xe2\u20ac\xb2}(t){K}_{D}{y}_{\mathrm{max}}}{{\left({y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)\right)}^{2}}.\end{array}
Substituting these terms in eq (1), we determine Î±(t) from the observed values of y(t) and y'(t):
\begin{array}{cc}\left(19\right)& \mathrm{\xce\pm}(t)=\frac{\mathrm{\xce\xb3}\xe2\u2039\dots {K}_{D}\xe2\u2039\dots y(t)}{{y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)}+{y}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots \left[1+\frac{{K}_{D}{y}_{\mathrm{max}}}{{\left({y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)\right)}^{2}}\right],\end{array}
once again allowing us to determine the effective calcium flux from observations of the indicator signal, related to y(t) as discussed below. The time course of the equivalent unperturbed calcium signal is determined as in eq. (16). Note again that \underset{{y}_{\mathrm{max}}\xe2\u2020\u2019\mathrm{\xe2\u02c6\u017e}}{\mathrm{lim}}\mathrm{\xce\pm}(t)={y}^{\xe2\u20ac\xb2}(t).
The critical question, of course, is under what circumstances this approximation is good. This requires the perturbation analysis in the appendix where we give the decay time constants (in dimensionless units) for small perturbations from the assumed solution (eq. (17)):
\begin{array}{l}{\mathrm{\xce\xbb}}_{1,2}=\xe2\u02c6\u2019C\xc2\pm D,\hfill \\ C=\frac{1}{2}\left(1+{R}_{b}+{R}_{f}{X}_{0}+\frac{{R}_{b}{R}_{f}}{{R}_{b}+{R}_{f}{X}_{0}}\right),\hfill \\ D=\sqrt{{C}^{2}\xe2\u02c6\u2019\left({R}_{b}+{R}_{f}{X}_{0}\right)}<C.\hfill \end{array}
Here, X_{0} is a positive constant. Both time constants are negative, indicating decay of a perturbation back to the assumed solution. These inverse time constants, in dimensional form, must be greater than the low pass filter time constant 1/Î³ for the free calcium concentration. This is true in the regime of large dimensionless forward and backward rates.
A numerical evaluation of the system of differential equations (eqs. (1) and (2)) is shown in Fig. 2. As calcium influx Î±(t) we used a whitenoise signal with a standard deviation of 5 Î¼ Mol/sec that was subsequently filtered by a 1^{st}order lowpass with 1 sec timeconstant and finally rectified (Fig. 2a). This signal was then fed into eqs. (1) and (2), using the following parameters: pump rate Î³ = 10 Hz, initial free indicator concentration y_{max} = 1 Î¼ Mol, indicator backward rate k_{b} = 10 Hz and indicator forward rate k_{f} = 10 Hz/Î¼ Mol. With these parameters, the resulting timecourse of the indicatorbound calcium is shown in Fig. 2b. As a comparison, we also show in Fig. 2b the indicatorbound calcium approximated by eq. (17). Both curves closely agree. In Fig. 2c, the indicatorbound calcium is shown as a function of the free cytosolic calcium, once (in black) as obtained from numerical integration of eqs. (1) and (2), once (in red) using the approximation using eq. (17). In this plot, certain deviations of the real signal from the approximate one can be observed. We subsequently quantified these deviations by calculating the rootmeansquare of the difference between the real and approximated signals. We did that for a total of 10,000 pairs of the two kinetic parameters R_{b} and R_{f} as defined above. Note that the parameters used in the above examples correspond to the values R_{b} = 1.0 and R_{f} = 1.0. The result is shown in Fig. 2d. The contour plot indicates that the rms values are smaller, i.e. the approximation is better, for larger R_{b} and R_{f} values. This is in close agreement with the result of our perturbation analysis.
Including internal buffering
Our mathematical analysis, for the sake of simplicity, has so far excluded the existence of internal buffers. In the following, we introduce an additional variable z(t), denoting the calcium bound internal buffer. We also give a superscript to the rate constants with 'y' referring to the calcium bound indicator, and 'z' referring to the calcium bound internal buffer. Consequently, we call the total (free and calcium bound) buffer concentration z_{max}. Writing down the basic dynamic equations gives:
(21) x'(t) = Î±(t)  Î³Â·x(t)  y'(t)  z'(t)
\begin{array}{cc}\left(22\right)& y(t)={k}_{f}^{y}\xe2\u2039\dots x(t)\xe2\u2039\dots \left[{y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)\right]\end{array}\xe2\u02c6\u2019{k}_{b}^{y}\xe2\u2039\dots y(t)
\begin{array}{cc}\left(23\right)& {z}^{\xe2\u20ac\xb2}(t)={k}_{f}^{z}\xe2\u2039\dots x(t)\xe2\u2039\dots \left[{z}_{\mathrm{max}}\xe2\u02c6\u2019z(t)\right]\xe2\u02c6\u2019{k}_{b}^{z}\xe2\u2039\dots z(t)\end{array}
Comparing these equations with our initial set (eqs. (1) and (2)), one realizes that an additional loss term has entered in eq. (21) to account for the calcium binding to the internal buffer. Eq. (22), which describes the binding to the indicator, is identical to eq. (2), and eq. (23) is a replication of eq. (2) with the buffer z substituting for the indicator y.
The steadystate solutions are:
\begin{array}{cccc}\left(24\right)& {x}_{\mathrm{\xe2\u02c6\u017e}}={\mathrm{\xce\pm}}_{C}/\mathrm{\xce\xb3};& {y}_{\mathrm{\xe2\u02c6\u017e}}={y}_{\mathrm{max}}\xe2\u2039\dots \frac{{\mathrm{\xce\pm}}_{C}}{{\mathrm{\xce\pm}}_{C}+{K}_{D}^{y}\xe2\u2039\dots \mathrm{\xce\xb3}};& {z}_{\mathrm{\xe2\u02c6\u017e}}={z}_{\mathrm{max}}\xe2\u2039\dots \frac{{\mathrm{\xce\pm}}_{C}}{{\mathrm{\xce\pm}}_{C}+{K}_{D}^{z}\xe2\u2039\dots \mathrm{\xce\xb3}};\end{array}
Thus, the steadystate solutions for free calcium and calciumbound indicator remain the same, no matter whether there is a buffer or not.
In the linear regime, the above equations reduce to the following system, now written in matrix notation for the sake of clarity:
\begin{array}{cc}\left(25\right)& \left(\begin{array}{c}{x}^{\xe2\u20ac\xb2}(t)\\ {y}^{\xe2\u20ac\xb2}(t)\\ {z}^{\xe2\u20ac\xb2}(t)\end{array}\right)=\left[\begin{array}{ccc}\xe2\u02c6\u2019\mathrm{\xce\xb3}\xe2\u02c6\u2019{k}_{f}^{y}{y}_{\mathrm{max}}\xe2\u02c6\u2019{k}_{f}^{z}{z}_{\mathrm{max}}& +{k}_{b}^{y}& +{k}_{b}^{z}\\ +{k}_{f}^{y}{y}_{\mathrm{max}}& \xe2\u02c6\u2019{k}_{b}^{y}& 0\\ +{k}_{f}^{z}{z}_{\mathrm{max}}& 0& \xe2\u02c6\u2019{k}_{b}^{z}\end{array}\right]\xe2\u2039\dots \left(\begin{array}{c}x(t)\\ y(t)\\ z(t)\end{array}\right)+\left(\begin{array}{c}\mathrm{\xce\pm}(t)\\ 0\\ 0\end{array}\right)\end{array}
The homogeneous part of this equation has the solutions kÂ·e^{Î»t}, where Î» is an Eigenvalue of the matrix, and k the respective Eigenvector. The timeconstants can, again, be obtained analytically from the characteristic (cubic) equation of the above matrix. The resulting expressions, however, are extremely lengthy and do not give any insight into the solution.
As a further step, we can also study the approximate solution in the nonlinear regime including an internal buffer. We use again the relationship from eq. (17):
and
\begin{array}{cc}\left(27\right)& z(t)={z}_{\mathrm{max}}\xe2\u2039\dots \frac{x(t)}{{K}_{D}^{z}+x(t)}\end{array}
From eq. (27), obtain the derivative z'(t):
\begin{array}{cc}\left(28\right)& {z}^{\xe2\u20ac\xb2}(t)={z}_{\mathrm{max}}\xe2\u2039\dots \frac{{x}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots {K}_{D}^{z}}{{\left({K}_{D}^{z}+x(t)\right)}^{2}}\end{array}
We rearrange eq. (26) to obtain
, and from that, calculate x'(t):
\begin{array}{cc}\left(30\right)& {x}^{\xe2\u20ac\xb2}(t)={K}_{D}^{y}\xe2\u2039\dots {y}_{\mathrm{max}}\frac{{y}^{\xe2\u20ac\xb2}(t)}{{\left({y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)\right)}^{2}}.\end{array}
Using eqs. (29) and (30), we now substitute x(t) and x'(t) in eq. (28) and obtain:
\begin{array}{cc}\left(31\right)& {z}^{\xe2\u20ac\xb2}(t)={z}_{\mathrm{max}}\xe2\u2039\dots {K}_{D}^{y}\xe2\u2039\dots {y}_{\mathrm{max}}\frac{{y}^{\xe2\u20ac\xb2}(t)}{{\left({y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)\right)}^{2}}\xe2\u2039\dots {K}_{D}^{z}\xe2\u2039\dots \frac{1}{{\left({K}_{D}^{z}+{K}_{D}^{y}\xe2\u2039\dots \frac{y(t)}{{y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)}\right)}^{2}}\end{array}
Now, we use eqs. (29), (30) and (31) and substitute in eq. (21). Rearranging for Î±(t) gives:
\begin{array}{cc}\left(32\right)& \mathrm{\xce\pm}(t)=\frac{\mathrm{\xce\xb3}\xe2\u2039\dots {K}_{D}^{y}\xe2\u2039\dots y(t)}{{y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)}+{y}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots \left[1+\frac{{K}_{D}^{y}\xe2\u2039\dots {y}_{\mathrm{max}}}{{\left({y}_{\mathrm{max}}\xe2\u02c6\u2019y(t)\right)}^{2}}+\frac{{K}_{D}^{y}\xe2\u2039\dots {y}_{\mathrm{max}}\xe2\u2039\dots {K}_{D}^{z}\xe2\u2039\dots {z}_{\mathrm{max}}}{{\left(y(t)\xe2\u2039\dots \left({K}_{D}^{z}\xe2\u02c6\u2019{K}_{D}^{y}\right)+{K}_{D}^{z}\xe2\u2039\dots {y}_{\mathrm{max}}\right)}^{2}}\right]\end{array}
Thus, the calcium influx can be determined in a manner similar to the situation without such a buffer (note how eq. (32) reduces to eq. (19) when z_{max} becomes zero). In order to do so, one also has to know the total amount of calciumbound and free internal buffer plus its binding constant.
Removing the indicator y(t) from eq. (21) and inserting eq. (28), the unperturbed calcium concentration x*(t) is the solution of the following nonlinear differential equation:
\begin{array}{cc}\left(33\right)& {x}^{\xe2\u02c6\u2014\xe2\u20ac\xb2}(t)=\left(\mathrm{\xce\pm}(t)\xe2\u02c6\u2019\mathrm{\xce\xb3}\xe2\u2039\dots {x}^{\xe2\u02c6\u2014}(t)\right)\xe2\u2039\dots {\left[1+\frac{{z}_{\mathrm{max}}\xe2\u2039\dots {K}_{D}^{z}}{{\left({K}_{D}^{z}+{x}^{\xe2\u02c6\u2014}(t)\right)}^{2}}\right]}^{\xe2\u02c6\u20191}\end{array}
This equation can be solved by numerical integration. Note again that eq. (33) reduces to eq. (16) when z_{max} becomes zero.
Discussion
In the work presented above we have derived, from first principles, the dependence of the timecourse of the indicator signal on the calcium influx and the relevant properties of the indicator and the cell under investigation. In order to do so, we assumed that the system approximately follows its steadystate at every point in time (eq. (17). Under these conditions, we were able to calculate the calcium influx from the indicator timecourse, no matter whether the free calcium concentration is in the linear or nonlinear range with respect to the binding constant of the indicator. Ignoring a cellinternal buffer system, this solution is represented by our eq. (19), from which the timecourse of the unperturbed calcium concentration can be derived by a simple convolution with a 1^{st} order lowpass filter, the timeconstant of which is given by the inverse of the pump rate, i.e. 1/Î³. Importantly, by using perturbation analysis, we were also able to indicate the parameter regime within which this solution is valid.
We also included an additional cell internal buffer in our model. Using the same approximation as above, i.e. eq. (17), we could calculate the calcium influx from the indicator timecourse (eq. (32)) and the timecourse of the unperturbed calcium concentration under these conditions (eq. (33)). In contrast to the situation without internal buffer, the unperturbed calcium concentration does not follow the calcium influx as fed through a linear, 1^{st} order lowpass filter but, instead, is altered by the dynamic interaction to and from the cellinternal buffer. In this case, however, we could not indicate the parameter range within which our solution is valid.
It is straightforward to see how the above approach can be extended to include several buffer systems. Nevertheless, our current analysis ignores some of the complexity that real nerve cells exhibit, such as feedback of the intracellular calcium level on to the membrane currents via calciumdependent Ca and Kconductances. While these can be included in numerical simulations of calcium dynamics, analytical treatment of the resulting equations are beyond the scope of the present paper and have to await future investigation.
Feasibility of the approach
To apply the approach outlined above to an experimental situation, one has to realize, first of all, that the indicator bound calcium (y(t) in our terminology) is not a parameter immediately being measured. Instead, what is immediately measured is a fluorescence signal. This is, of course, related to the indicator bound calcium, and the quantification of this relationship is given in Appendix II. Nevertheless, the application of our approach to an experimental situation, in particular in the nonlinear regime, has some shortcomings. First of all, application of eqs. (15) or (19) requires knowledge of parameters such as extrusion rate, initial indicator concentration etc. If these are not known, the calcium influx cannot be calculated. But even if all these parameters are known, the application of eqs. (15) or (19) is problematic since the indicator signal will be subject to noise. In this event, taking the first or second order derivatives of a measured signal will boost the noise, and, dividing by small values of (y_{max}  y(t))^{2} (when the bound indicator is saturating, i.e. approaching the initial free indicator concentration) will further lead to unstable solutions. Therefore, alternative approaches should be considered.
Alternative approach I: linear regime
One such alternative approach is applicable when the relationship between the membrane voltage and the calcium influx and the indicator signal is linear through all stages. While the first relationship, i.e. the one between membrane voltage and calcium influx, is in general not linear, one can either work with small membrane deviations around a potential where calcium channels are already activated, or use the number of action potentials of the actual membrane potential as the signal V(t). The method outlined below requires measuring the voltage signal V(t) and indicator signal y(t) simultaneously. Then we can determine the relationship between the voltage and the bound indicator time course, and from the latter determine V(t). If we know V(t), we can use an equation for calcium dynamics to predict the calcium influx Î±(t).
In the linear regime we can do this by assuming that y(t) is given by a first order kernel g(t) in terms of V(t)

(34)
y(t) = âˆ«dt' g(t  t')V (t').
From several such example recordings, the optimal reverse filter g_{rev}(t) can be calculated in the Fourier domain using the WienerKolmogorov formalism if the calcium concentrations are small compared to the K_{D} value of the indicator, i.e. when the system is in the linear regime. Under these conditions, the bound indicator concentration can be calculated from the calcium influx as a convolution with the following socalled 'forward' filter g_{forw}(t) (see eq. (13)):
\begin{array}{cc}\left(35\right)& {g}_{forw}(t)=\frac{{y}_{\mathrm{max}}}{2\xe2\u2039\dots \mathrm{\xce\xb3}\xe2\u2039\dots {K}_{D}}\xe2\u2039\dots \left[\left(1+\frac{A}{\sqrt{{A}^{2}\xe2\u02c6\u20194{k}_{b}\mathrm{\xce\xb3}}}\right)\xe2\u2039\dots \mathrm{exp}\left(\xe2\u02c6\u2019t/{\mathrm{\xcf\u201e}}_{1}\right)+\left(1\xe2\u02c6\u2019\frac{A}{\sqrt{{A}^{2}\xe2\u02c6\u20194{k}_{b}\mathrm{\xce\xb3}}}\right)\xe2\u2039\dots \mathrm{exp}\left(\xe2\u02c6\u2019t/{\mathrm{\xcf\u201e}}_{2}\right)\right]\end{array}
Given that there is a linear relationship between membrane voltage and calcium influx, the problem of recovering the membrane voltage from indicator measurements is to find the optimal reverse filter, which can then be applied to all those situations where only the optical signal from the calciumbound indicator y(t) has been measured. As can be shown, the optimal reverse filter g_{rev}(t) is not the inverse of the forward filter that turns V(t) into y(t) (as done by Yaksi and Friedrich, [9]), but rather the average crosscorrelation between V(t) and y(t), divided by the power spectrum of y(t) [10, 11]. Denoting the inverse Fourier Transform by F^{1}, y*(f) the complex conjugate of y(f) and the average across n trials by âŸ¨...âŸ©, the optimal reverse filter g_{rev}(t) becomes:
\begin{array}{cc}\left(36\right)& {g}_{rev}(t)={F}^{\xe2\u02c6\u20191}\{g(f)\}={F}^{\xe2\u02c6\u20191}\left\{\frac{\xe3\u20ac\u02c6{V}_{i}(f)\xe2\u2039\dots {y}_{i}^{\xe2\u02c6\u2014}(f)\xe3\u20ac\u2030}{\xe3\u20ac\u02c6{y}_{i}(f)\xe2\u2039\dots {y}_{i}(f)\xe3\u20ac\u2030}\right\}\end{array}
Convolving each new optical signal y(t) with g_{rev}(t) then results in the optimal estimate of the voltage signal, leading to a calcium influx and consequently to the optical signal of bound indicator. Clearly, the advantage of this method is that no parameters need to be known; the disadvantage is that enough dual measurements of membrane voltage and indicator need to be at hand to calculate the optimal reverse filter g_{rev}(t). As another caveat, this method only works as long as calcium concentrations are in the linear regime with respect to the K_{D} of the indicator and to membrane voltage. An example of a reverse reconstruction in the linear regime is shown in Fig. 3. Here, the signal was created by Gaussian noise with an autocorrelation timeconstant of 100 ms that was subsequently rectified. From this influx, the calcium bound indicator concentration y(t) was numerically determined using eqs. (1) and (2) and the following parameters: pump rate Î³ = 20 Hz, k_{f} = 0.01 1/(nMol sec), k_{b} = 10 Hz, resulting in a K_{D} value of 1000 nMol, and an initial indicator concentration y_{max} = 100 nMol. This led to the average time course of calciumbound indicator y(t) shown as a black line in Fig. 3c. Through 100 trials, a Gaussian noise signal was added with an autocorrelation timeconstant of 10 ms, which had an average amplitude of 5% of y(t). Twelve such trials are shown as grey lines superimposed on Fig. 3c. From these trials, optimal forward g_{forw}(t) and reverse filter g_{rev}(t) were calculated according to {g}_{forw}(t)={F}^{\xe2\u02c6\u20191}\left\{\xe3\u20ac\u02c6{y}_{i}(f)\xe2\u2039\dots {\mathrm{\xce\pm}}_{i}^{\xe2\u02c6\u2014}(f)\xe3\u20ac\u2030/\xe3\u20ac\u02c6{\mathrm{\xce\pm}}_{i}(f)\xe2\u2039\dots {\mathrm{\xce\pm}}_{i}^{\xe2\u02c6\u2014}(f)\xe3\u20ac\u2030\right\} and {g}_{rev}(t)={F}^{\xe2\u02c6\u20191}\left\{\xe3\u20ac\u02c6{\mathrm{\xce\pm}}_{i}(f)\xe2\u2039\dots {y}_{i}^{\xe2\u02c6\u2014}(f)\xe3\u20ac\u2030/\xe3\u20ac\u02c6{y}_{i}(f)\xe2\u2039\dots {y}_{i}^{\xe2\u02c6\u2014}(f)\xe3\u20ac\u2030\right\}. These filters are shown in b and d, respectively. Applying the forward filter to Î±(t), the signal shown in red in Fig. 3c was obtained. Applying the reverse filter to y(t), the signal shown in red in Fig. 3a was obtained. Note that while the forward filter leads to an output that is almost indistinguishable from y(t), the reverse filter can only reconstruct the lowfrequency components of Î±(t), since high frequency components are covered by noise in the individual response trials.
Alternative approach II: nonlinear regime
If either of the two relationships, i.e. the one between membrane voltage V(t) and the calcium influx Î±(t) or the one between Î±(t) and bound indicator y(t) (due to a high calcium level with respect to the K_{D} value of the indicator) is nonlinear, the above method will lead to erroneous results. In such a nonlinear regime we must use a different approach. Here, too, we must measure the indicator signal y(t) and the membrane voltage V(t) simultaneously. In the reconstructed state space [12] of the voltage measurement we can fully describe the state of the system (neuron plus indicator) using the voltage and its time lags, or we can use the indicator signal and its time lags. If we do the latter, we create data vectors
(37) U(t) = [y(t), y(t  Ï„), y(t  2Ï„), ..., y(t  (D  1)Ï„)],
where the number of lags D and the time lag Ï„ are respectively determined by the method of false nearest neighbors and by average mutual information. For every U(t) there is an associated indicator signal V(t), and since we have already totally characterized the state of the neuron by U(t) there must be a nonlinear relationship V(t) = f(U(t)). We can discover this nonlinear relation from the simultaneous measurements of y(t) and V(t), then, just as in the linear case, map new measurements of y(t) to allow us to predict V(t).
The method requires determining f(U(t)). To accomplish this, we represent f(U) in terms of some basis functions chosen by the user: Ï†_{ m }(U), and write
\begin{array}{cc}\left(38\right)& f(U)={\displaystyle \underset{m=1}{\overset{M}{\xe2\u02c6\u2018}}{c}_{m}{\mathrm{\xcf\u2022}}_{m}(U)}.\end{array}
In the state space of the U's, each state vector has many neighbors U^{(l)}(t); l = 0,1...,N_{ B }; U^{(0)}(t) = U(t). Each of these neighbors maps into a voltage V^{(l)}(t) = f(U^{(l)}(t)). At any given time, corresponding to a location in U space, we can determine the coefficients c_{ m }by minimizing the least squares form
\begin{array}{cc}\left(39\right)& {\displaystyle \underset{l=0}{\overset{{N}_{B}}{\xe2\u02c6\u2018}}{\left({V}^{(l)}(t)\xe2\u02c6\u2019{\displaystyle \underset{m=1}{\overset{M}{\xe2\u02c6\u2018}}{c}_{m}(t){\mathrm{\xcf\u2022}}_{m}({U}^{(l)}(t))}\right)}^{2}.}\end{array}
This establishes the map V(t) = f(U(t)) locally in U space. Now we make a new measurement of y_{ new }(t). Use this to create a new Ddimensional data vector U_{ new }(t) = [y_{ new }(t), y_{ new }(t  Ï„), y_{ new }(t  2Ï„), ..., y_{ new }(t  (D  1)Ï„)]. Search among all the data vectors in the initial training set and find the one is closest to U_{ new }(t); suppose it is U(t'). Then using the local map attached to U(t') we predict
\begin{array}{cc}\left(40\right)& {V}_{new}(t)={\displaystyle \underset{m=1}{\overset{M}{\xe2\u02c6\u2018}}{c}_{m}({t}^{\xe2\u20ac\xb2}){\mathrm{\xcf\u2022}}_{m}({U}_{new}(t)).}\end{array}
From the time course of new measurements y_{ new }(t), we are thus able to use the learned map to predict the time course of the new membrane voltage V_{ new }(t), which was our goal.
Relationship to previous studies
Previous studies on calcium binding mainly considered steadystate situations or the linear case, i.e. that calcium concentrations are small compared to the dissociation constant K_{D} of the indicator [13;14]. In particular, a number of studies investigated the diffusion of Calcium ions in the presence of buffers [15â€“19]. However, none of these interesting papers focused on the time dependent nonlinear kinetics without diffusion or addressed the temporal stability of our eqs. (1) and (2) for the approximate linear or the approximate nonlinear solutions that we derived above.
Our study can be related to these previous investigations when we combine our approximation about the dynamics of the system without internal buffer (eq. (17)) with the condition of small free calcium concentration. Thus, eq. (17) becomes:
, and
\begin{array}{cc}\left(42\right)& {y}^{\xe2\u20ac\xb2}(t)\end{array}\xe2\u2030\u02c6{x}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots \frac{{y}_{\mathrm{max}}}{{K}_{D}}
Inserting eq. (41) into eq. (1) leads to:
\begin{array}{cc}\left(43\right)& {x}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots \frac{1}{\mathrm{\xce\xb3}}\xe2\u2039\dots \left(1+\frac{{y}_{\mathrm{max}}}{{K}_{D}}\right)=\frac{\mathrm{\xce\pm}(t)}{\mathrm{\xce\xb3}}\xe2\u02c6\u2019x(t)\end{array}
Eq. (43) describes a 1^{st} order lowpass filter with a timeconstant equal to
\begin{array}{cc}\left(44\right)& \mathrm{\xcf\u201e}=\frac{1}{\mathrm{\xce\xb3}}\left(1+\frac{{y}_{\mathrm{max}}}{{K}_{D}}\right)\end{array}
With the indicator concentration being small, the timeconstant becomes 1/Î³. Large indicator concentrations, therefore, increase the timeconstant from 1/Î³ to the value indicated by eq. (44).
Repeating the above for the situation with an internal buffer, eq. (42) remains unaltered. In a similar way, we derive from eq. (27)
, and from that
\begin{array}{cc}\left(46\right)& {z}^{\xe2\u20ac\xb2}(t)\xe2\u2030\u02c6{x}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots \frac{{z}_{\mathrm{max}}}{{K}_{D}^{z}},\end{array}
Substituting eqs. (42) and (46) into eq. (21) gives
\begin{array}{cc}\left(47\right)& {x}^{\xe2\u20ac\xb2}(t)=\mathrm{\xce\pm}(t)\xe2\u02c6\u2019\mathrm{\xce\xb3}\xe2\u2039\dots x(t)\xe2\u02c6\u2019{x}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots \frac{{y}_{\mathrm{max}}}{{K}_{D}^{y}}\xe2\u02c6\u2019{x}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots \frac{{z}_{\mathrm{max}}}{{K}_{D}^{z}}\end{array}
Rearranging leads to
\begin{array}{cc}\left(48\right)& {x}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots \frac{1}{\mathrm{\xce\xb3}}\left(1+\frac{{y}_{\mathrm{max}}}{{K}_{D}^{y}}+\frac{{z}_{\mathrm{max}}}{{K}_{D}^{z}}\right)=\frac{\mathrm{\xce\pm}(t)}{\mathrm{\xce\xb3}}\xe2\u02c6\u2019x(t)\end{array}
This, again, describes a 1^{st} order lowpass filter with a timeconstant equal to:
\begin{array}{cc}\left(49\right)& \mathrm{\xcf\u201e}=\frac{1}{\mathrm{\xce\xb3}}\left(1+\frac{{y}_{\mathrm{max}}}{{K}_{D}^{y}}+\frac{{z}_{\mathrm{max}}}{{K}_{D}^{z}}\right).\end{array}
Comparing this result to eq. (44), one can see that internal buffering enlarges the time constant by an additive term, equivalent to the one introduced by the indicator. Eq. (49) is identical to eq. (2) in [20].
Neher and Augustine [21] defined the calcium binding capacity as the ratio of the change in bound indicator concentration over the change in free Calcium:
\begin{array}{cc}\left(50\right)& \mathrm{\xce\xba}=\frac{\mathrm{\xce\u201d}y(t)}{\mathrm{\xce\u201d}x(t)}\end{array}
For the linear case, i.e. when the calcium concentrations are small compared to the dissociation constant K_{D}, this quantity is identical to y_{max}/K_{D}, as can be derived from eq. (42). For the nonlinear case, i.e. when the calcium concentrations are large compared to the dissociation constant K_{D}, we first derive from eq. (17)
\begin{array}{cc}\left(51\right)& {y}^{\xe2\u20ac\xb2}(t)={y}_{\mathrm{max}}\xe2\u2039\dots \frac{{x}^{\xe2\u20ac\xb2}(t)\xe2\u2039\dots {K}_{D}}{{\left({K}_{D}+x(t)\right)}^{2}}\end{array}
Using eqs. (17) and (50), the calcium binding capacity becomes:
\begin{array}{cc}\left(52\right)& \mathrm{\xce\xba}=\frac{{y}_{\mathrm{max}}\xe2\u2039\dots {K}_{D}}{{\left({K}_{D}+x(t)\right)}^{2}}\end{array}
This, again, is identical to eq. (3) in [20].
In their study on Calcium diffusion, Naraghi and Neher [17] investigated a linearized mathematical model of diffusion and kinetics. For one buffer their results are contained in their eqs (AII.9) and (AII.10). This corresponds to our analysis when one sets the calcium source, our Î±(t), and our kinetic loss terms for free calcium, our Î³ x(t), both to zero, accounting for the zero eigenvalues they report. Majewska et al. [19] focused on determining time constants for intracellular calcium kinetics experimentally. Their results show time scales in the range of 100 s of ms. This is the order of magnitude we have used in selecting our pump rate Î³ whose value we chose as 1/Î³ = 100 ms in our numerical simulations. Wagner and Keizer [16] again focused on diffusion of calcium. Their notation identifies a free calcium concentration [Ca^{2+}], which is just our x(t), and a concentration of calcium bound to a mobile buffer [CaB_{m}], which is precisely our y(t), and finally the concentration of the mobile buffer itself [B_{m}], which is our y_{max}y(t). They do not have source terms for the calcium influx, our Î±(t), or kinetic loss terms for free calcium, our Î³ x(t). Ignoring diffusion and these sources and sinks of [Ca^{2+}], their eqs. (2), (3), and (4), are precisely our eqs. (1) and (2) above. It is important they do not analyze the nonlinear ordinary (kinetic) differential equations that result when diffusion is not important. Using the estimates of Zhou and Neher for the diffusion constants to be about 300 Î¼ m^{2}/s this translates to a time for diffusion over a cellular scale to be about 3 ms which is much shorter than the kinetic time constants we consider or are discussed by Majewska et al. [19]. This gives our rationale for ignoring diffusion and focusing on properties of the nonlinear kinetics.
Appendix I: perturbation analysis of the nonlinear solution
We begin by making eqs. (1) and (2) dimensionless. There are three quantities with the dimensions of (time)^{1}: k_{b, Î³}, and k_{f}y_{max}. We express our indicator kinetic equations in terms of the two dimensionless variables which can be made from these
(A1) R_{ b }= k_{ b }/Î³; R_{ f }= k_{ f }y_{max}/Î³.
We also scale x(t) and y(t) with the initial indicator concentration y_{max} and the time by the pump rate Î³:
(A2) x(t) â†’ y_{max} X(t); y(t) â†’ y_{max} Y(t); t â†’ t/Î³;
Thus, in these new dimensionless variables, free calcium and calciumbound indicator concentrations are given as fractions of the initial free indicator concentration, and the forward and backward rates are given relative to the pump rate.
The kinetic eqs. (1) and (2) now become:
(A3) X'(t) = Î±(t)/(Î³Â·y_{max})  X(t)  Y'(t)
(A4) Y'(t) = R_{ f }X(t)Â·(1  Y(t))  R_{ b }Y(t).
We chose as an approximate solution of these equations functions
\begin{array}{cc}\left(\text{A}5\right)& {Y}_{0}(t)=\frac{{R}_{f}{X}_{0}(t)}{{R}_{b}+{R}_{f}{X}_{0}(t)},\end{array}
for which
\begin{array}{cc}\left(\text{A}6\right)& {{Y}^{\xe2\u20ac\xb2}}_{0}(t)=\frac{{R}_{f}{R}_{b}}{{\left[{R}_{b}+{R}_{f}{X}_{0}(t)\right]}^{2}}\xe2\u2039\dots {{X}^{\xe2\u20ac\xb2}}_{0}(t).\end{array}
This solution is suggested by the vanishing of the right hand side of eq. (A4) as well as by the fixed point solution, true when X(t) is time independent. Another motivation for this approximate solution is that when both R_{b} and R_{f} are large, the right hand side of eq. (A4) would make the rate of change of the calcium bound indicator vary quite rapidly unless the balance indicated by eq. (A5) were maintained.
To determine when this solution is accurate, we make perturbations
(A7) X(t) = X_{0}(t) + Î”_{ X }(t)
(A8) Y(t) = Y_{0}(t) + Î”_{ Y }(t),
and linearize the equations in Î”_{ X }(t) and Î”_{ Y }(t). From eqs. (A3) and (A4) we obtain, to first order in the perturbations,
\begin{array}{cc}\left(\text{A}9\right)& \frac{d{\mathrm{\xce\u201d}}_{X}(t)}{dt}=\xe2\u02c6\u2019{\mathrm{\xce\u201d}}_{X}(t)\left[1+\frac{{R}_{b}{R}_{f}}{\mathrm{\xce\xb7}}\right]+\mathrm{\xce\xb7}{\mathrm{\xce\u201d}}_{Y}(t)\end{array}
\begin{array}{cc}\left(\text{A}10\right)& \frac{d{\mathrm{\xce\u201d}}_{Y}(t)}{dt}={\mathrm{\xce\u201d}}_{X}(t)\frac{{R}_{f}{R}_{b}}{\mathrm{\xce\xb7}}\xe2\u02c6\u2019\mathrm{\xce\xb7}{\mathrm{\xce\u201d}}_{Y}(t)\xe2\u02c6\u2019\frac{d{Y}_{0}(t)}{dt},\end{array}
where
(A11) Î· = R_{ b }+ R_{ f }X_{0}(t),
and X_{0}(t) satisfies
\begin{array}{cc}\left(\text{A}12\right)& {{X}^{\xe2\u20ac\xb2}}_{0}(t)=\frac{\mathrm{\xce\pm}(t)}{\mathrm{\xce\xb3}\xe2\u2039\dots {y}_{\mathrm{max}}}\xe2\u02c6\u2019{X}_{0}(t).\end{array}
These are very similar to those for the linearized problem discussed in the text. The key differences are that Î· = R_{ b }+ R_{ f }X_{0}(t) is time dependent and there is an inhomogeneous term in the equations for Î”_{ Y }(t). Since the solution for the unperturbed X_{0}(t) is a low pass filtered version of the calcium influx, we take it as a positive constant, slowly varying, in the perturbation equations.
The equation for Î”(t) = (Î”_{ X }(t), Î”_{ Y }(t)) written in matrix form is
\begin{array}{cc}\left(\text{A}13\right)& \frac{d\mathrm{\xce\u201d}(t)}{dt}=M\mathrm{\xce\u201d}(t)+F(t),\end{array}
with F(t)=(0,\xe2\u02c6\u2019\frac{d{Y}_{0}(t)}{dt}). The eigenvalues of the matrix M are
(A14) Î»_{1,2} = C Â± D
With C=\frac{1}{2}\left[1+\mathrm{\xce\xb7}+\frac{{R}_{b}{R}_{f}}{\mathrm{\xce\xb7}}\right]
and D=\frac{1}{2}\sqrt{{\left[1+\mathrm{\xce\xb7}+\frac{{R}_{b}{R}_{f}}{\mathrm{\xce\xb7}}\right]}^{2}\xe2\u02c6\u20194\mathrm{\xce\xb7}}=\frac{1}{2}\sqrt{{\left[1\xe2\u02c6\u2019\mathrm{\xce\xb7}+\frac{{R}_{b}{R}_{f}}{\mathrm{\xce\xb7}}\right]}^{2}+4{R}_{b}{R}_{f}},,
so C, D >0 and D<C, so both eigenvalues are negative. This means the solution is trying to drive Î”(t) = (Î”_{ X }(t), Î”_{ Y }(t)) to zero, with "bumps" from the forcing term. If the forcing term is bounded above, that is the derivative (A6) remains below some maximum value while the calcium current is flowing, the solutions Î”(t) = (Î”_{ X }(t), Î”_{ Y }(t)) go to zero. In these weak conditions, the assumed solutions (X_{0}(t),Y_{0}(t)) are stable. The larger timeconstant is shown in Fig. 4 as a function of R_{b} and R_{f}. It agrees with the rms values shown in Fig. 2d.
Appendix II: relating the fluorescence signal to calciumbound indicator concentration
From eqs. (19) and (32), it is important to note that Î±(t) does not scale with y(t). Therefore, the indicator concentration enters these equations as an absolute concentration. Otherwise, the calculated timecourse of the calcium influx will be incorrect (and not just by a factor!). Usually, however, the indicator concentration is not available directly, but rather as fluorescence values, in most cases as Î”F/F, i.e. fluorescence changes relative to a reference fluorescence F(0) obtained just before the start of an experiment. The fluorescence value F(t) is the sum of the fluorescence of the indicator with bound calcium. i.e. y(t), and free indicator concentrations, i.e. z_{0}y(t), each one contributing to the total fluorescence by a factor f_{b} (bound) and f_{f} (free), respectively:
(A15) F(t) = f_{ b }Â·y(t) + f_{ f }Â·(y_{max}  y(t))
The factors f_{b} and f_{f} both can be determined experimentally from the fluorescence of a calciumfree and a calciumsaturated indicator solution. Using the maximum fluorescence change (Î”F/F)_{max} of the indicator as (f_{b}  f_{f})/f_{f}, the following relation then holds between y(t) and Î”F/F:
\begin{array}{cc}\left(\text{A}16\right)& \frac{\mathrm{\xce\u201d}F}{F}(t)=\frac{F(t)\xe2\u02c6\u2019F(0)}{F(0)}=\frac{\left[y(t)\xe2\u02c6\u2019y(0)\right]\xe2\u2039\dots {\left(\mathrm{\xce\u201d}F/F\right)}_{\mathrm{max}}}{y(0)\xe2\u2039\dots {\left(\mathrm{\xce\u201d}F/F\right)}_{\mathrm{max}}+{y}_{\mathrm{max}}}\end{array}
Solving eq (A16) for y(t) yields:
\begin{array}{cc}\left(\text{A}17\right)& y(t)=\frac{\mathrm{\xce\u201d}F}{F}(t)\xe2\u2039\dots \left(y(0)+\frac{{y}_{\mathrm{max}}}{\left(\mathrm{\xce\u201d}F/{F}_{\mathrm{max}}\right)}\right)+y(0).\end{array}
When using indicators based on fluorescence resonance energy transfer ('FRET'), results are usually expressed in the relative change of the fluorescence ratio obtained at two different wavelengths, one from the donor fluorophore F_{1}, and the other from the acceptor fluorophore F_{2}, respectively:
\begin{array}{cc}\left(\text{A}18\right)& \frac{\mathrm{\xce\u201d}R}{R}=\frac{{F}_{1}(t)/{F}_{2}(t)\xe2\u02c6\u2019{F}_{1}(0)/{F}_{2}(0)}{{F}_{1}(0)/{F}_{2}(0)}=\frac{{F}_{1}(t){F}_{2}(0)}{{F}_{2}(t){F}_{1}(0)}\xe2\u02c6\u20191\end{array}
Inserting eq. (A15) for each wavelength, the relation to the indicator concentration becomes:
\begin{array}{cc}\left(\text{A}19\right)& \frac{\mathrm{\xce\u201d}R}{R}=\frac{{y}_{\mathrm{max}}\xe2\u2039\dots \left(y(t)\xe2\u02c6\u2019y(0)\right)\xe2\u2039\dots \left({f}_{b1}{f}_{f2}\xe2\u02c6\u2019{f}_{b2}{f}_{f1}\right)}{\left[y(t)\xe2\u2039\dots \left({f}_{b2}\xe2\u02c6\u2019{f}_{f2}\right)+{y}_{\mathrm{max}}{f}_{f2}\right]\xe2\u2039\dots \left[y(t)\xe2\u2039\dots \left({f}_{b1}\xe2\u02c6\u2019{f}_{f1}\right)+{y}_{\mathrm{max}}{f}_{f1}\right]}\end{array}
References
Grynkiewicz G, Poenie M, Tsien RY: A new generation of Ca indicators with greatly improved flourescence properties. J Biol Chem. 1985, 260: 34403450.
Tank DW, Sugimori M, Connor JA, LlinÃ¡s RR: Spatially resolved calcium dynamics of mammalian purkinje cells in cerebellar slice. Science. 1988, 242: 773777. 10.1126/science.2847315.
Borst A, Egelhaaf M: In vivo imaging of calcium accumulation in fly interneurons as elicited by visual motion stimulation. PNAS. 1992, 89: 41394143. 10.1073/pnas.89.9.4139.
Sobel EC, Tank DW: In vivo Ca^{2+} dynamics in a cricket auditory neuron: an example of chemical computation. Science. 1994, 263: 823826. 10.1126/science.263.5148.823.
Miyawaki A, Llopis J, Heim R, McCaffery JM, Adams JA, Ikura M, Tsien RY: Fluorescent indicators for Ca^{2+} based on green fluorescent proteins and calmodulin. Nature. 1997, 388: 882887. 10.1038/42264.
Single S, Borst A: Dendritic integration and its role in computing image velocity. Science. 1998, 281: 18481850. 10.1126/science.281.5384.1848.
Pologruto TA, Yasuda R, Svoboda K: Monitoring neural activity and [Ca2+] with genetically encoded Calcium indicators. J Neuroscience. 2004, 24: 95729579. 10.1523/JNEUROSCI.285404.2004.
Reiff DF, Ihring A, Guerrero G, Isacoff EY, Joesch M, Nakai J, Borst A: In vivo comparison of 11 GFPbased indicators of neural activity in genetically engineered flies. J Neuroscience. 2005, 25: 47664778. 10.1523/JNEUROSCI.490004.2005.
Yaksi E, Friedrich RW: Reconstruction of firing rate changes across neuronal populations by temporally deconvolved Ca^{2+} imaging. Nature Methods. 2006, 3: 377383. 10.1038/nmeth874.
Bialek W, Rieke F, de Ruyter van Steveninck RR, Warland D: Reading a neural code. Science. 1991, 252: 18541857. 10.1126/science.2063199.
Borst A, Theunissen FE: Information theory and neural coding. Nature Neuroscience. 1999, 2: 947957. 10.1038/14731.
Abarbanel HDI: Analysis of Observed Chaotic Data. 1996, New York: SpringerVerlag
Neher E: The use of Fure2 for estimating Ca buffers and Ca fluxes. Neuropharmacology. 1995, 34: 14231442. 10.1016/00283908(95)00144U.
Neher E: Usefulness and limitations of linear approximations to the understanding of Ca^{2+} signals. Cell Calcium. 1998, 24: 345357. 10.1016/S01434160(98)900586.
Zhou Z, Neher E: Mobile and immobile calcium buffers in bovine adrenal chromafin cells. J Physiol. 1993, 469: 245273.
Wagner Keizer J: Effects of rapid buffers on Ca^{2+} diffusion and Ca^{2+} oscillations. Biophys J. 1994, 67: 447456.
Naraghi M, Neher E: Linearized buffered Ca^{2+} diffusion in microdomains and its implications for calculation of Ca^{2+} at the mouth of a Calcium channel. J Neurosci. 1997, 17: 69616973.
Pape P, DeShien J, Chandler WK: Effects of partial sarcoplasmic reticulum Calcium depletion on Calcium release in frog cut muscle fibers equilibrated with 20 mM EGTA. J Gen Physiol. 1998, 112: 263295. 10.1085/jgp.112.3.263.
Majewska A, Brown E, Ross J, Yuste R: Mechanisms of Calcium decay kinetics in hippocampal spines: Role of spine Calcium pumps and Calcium diffusion through the spine neck in biochemical compartmentalization. J Neurosci. 2000, 20: 17221734.
Jackson MB, Redman SJ: Calcium dynamics, buffering, and buffer saturation in the boutons of dentate granulecell axons in the Hilus. J Neurosci. 2003, 23: 16121621.
Neher E, Augustine GJ: Calcium gradients and buffers in bovine chromaffin cells. J Physiology. 1992, 450: 273301.
Acknowledgements
We are grateful to E. Neher, J. Mueller and F. Theunissen for fruitful discussions and to D. Spavieri and three anonymous referees for carefully reading previous versions of the ms.
Author information
Authors and Affiliations
Corresponding author
Authorsâ€™ original submitted files for images
Below are the links to the authorsâ€™ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Borst, A., Abarbanel, H.D. Relating a calcium indicator signal to the unperturbed calcium concentration timecourse. Theor Biol Med Model 4, 7 (2007). https://doi.org/10.1186/1742468247
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1742468247