362
OPTICS LETTERS / Vol. 34, No. 3 / February 1, 2009
Fast-fluorescence dynamics in nonratiometric
calcium indicators
M. Gersbach,1 D. L. Boiko,2,* C. Niclass,1 C. C. H. Petersen,1 and E. Charbon1
1
Ecole Polytechnique Fédérale de Lausanne, Quantum Architecture Group, 1015, Lausanne, Switzerland
2
Centre Suisse d’Electronique et de Microtechnique SA, 2002, Neuchâtel, Switzerland
*Corresponding author: dmitri.boiko@csem.ch
Received September 4, 2008; revised November 10, 2008; accepted November 25, 2008;
posted December 4, 2008 (Doc. ID 101111); published January 30, 2009
The fluorescence decay of high-affinity nonratiometric Ca2+ indicator Oregon Green BAPTA-1 (OGB-1) is
analyzed with unprecedented temporal resolution in the two-photon excitation regime. A triple exponential
decay is shown to best fit the fluorescence dynamics of OGB-1. We provide a model for accurate measurements of the free Ca2+ concentration and dissociation constants of nonratiometric calcium indicators.
© 2009 Optical Society of America
OCIS codes: 170.6280, 160.2540, 260.2510, 170.2520, 300.6280, 180.2520.
Fluorescence lifetime imaging microscopy (FLIM) is
used to locally probe the chemical environment of
fluorophores, e.g., ion concentration, pH, or oxygen
content [1,2]. To acquire time-resolved fluorescence
images, the technique of time-correlated singlephoton counting (TCSPC), in combination with detectors exhibiting single-photon sensitivity, is commonly
used. This technique enables the measurement of
photon time-of-arrival distributions with very high
accuracies, independently of instabilities in the excitation beam intensity. So far, temporal resolutions of
a few hundred picoseconds were considered sufficient
in biomedical FLIM applications.
In this Letter, we analyze the fluorescence decay of
the high-affinity Ca2+ indicator Oregon Green
BAPTA-1 (OGB-1) under two-photon excitation conditions, using a TCSPC system of measured resolution 79 ps based on complementary metal-oxide
semiconductor single-photon avalanche diode (SPAD)
detector technology. Our measurements reveal a
triple-exponential decay of OGB-1 fluorescence,
which we show enables accurate measurements of
Ca2+ concentration and dissociation constants of nonratiometric fluorescent probes. We provide a comparison with previously reported data [1,3,4].
The SPAD-based TCSPC system is depicted in Fig.
1. Our SPADs [5] are integrated in a 32⫻ 32 array
and incorporate on-chip high-bandwidth input/
output circuitry. The active region of a SPAD pixel
[Fig. 1(b)] consists of a p+-n junction operating in Geiger mode. Owing to a small diameter 共7 m兲, SPADs
show an extremely low dark-count rate (⬍10 Hz at
room temperature). The photon-detection probability
is 25% at 500 nm wavelength; the dead time is 25 ns
with negligible afterpulsing 共⬍0.1% 兲 [5].
Fluorescent molecules are excited in the twophoton regime [6] using a mode-locked Ti:sapphire
laser (MaiTai, Spectra Physics) emitting 100 fs optical pulses at 800 nm wavelength [Fig. 1(a)]. The attenuated beam with an average power of 9 mW is focused on a sample using a 20⫻ microscope objective
(XLUMPlanFL, Olympus), which also serves to collect the fluorescent emission. Fluorescence from
0146-9592/09/030362-3/$15.00
OGB-1 samples is directed toward the SPAD by a dichroic beam splitter (DBS) and a filter (E650SP,
Chroma Technology) for suppression of backscattered
excitation pulses. Another 20⫻ objective images the
emission spot onto the SPAD. The low fluorescence
(c)
(d)
Fig. 1. (Color online) (a) Schematic of the experimental
setup for fluorescence lifetime measurements. (b) Scanning
electron microscope image of a SPAD. (c) Measured fluorescence decay of OGB-1 at various calcium concentrations
(gray curves) and hyper-Rayleigh scattering from colloidal
gold particles (open circles). For a 38 nM free Ca2+ concentration buffer, the numerical fits of the double-exponential
(dotted curve, 2r = 1.80) and triple-exponential (dashed
curve, 2r = 1.09) decay models are indicated. The inset
shows a close-up of the initial interval of 600 ps width, and
(d) shows the residues for the two models.
© 2009 Optical Society of America
February 1, 2009 / Vol. 34, No. 3 / OPTICS LETTERS
yield of the two-photon excitation results in a count
rate at the detector of 10 kHz for an 80 MHz repetition rate of excitation pulses. Time discrimination is
performed in a “reversed start–stop” configuration
using a 6 GHz bandwidth oscilloscope (WaveMaster,
LeCroy) incorporating a time-to-digital converter.
Figure 1(c) shows the instrument response function (IRF) of the entire system at 400 nm wavelength, recorded by measuring the hyper-Rayleigh
scattering of 800 nm wavelength pulses in a solution
of colloidal gold particles (G1652, Sigma-Aldrich) [7].
At the incident power of 90 mW, the average count
rate of the detector is just 600 Hz. The measured response time jitter is 79 ps (FWHM). The IRF is only
slightly asymmetrical and assumes a Gaussian-curve
2
兲 (inset).
approximation IRF= exp共−t2 / 2IRF
Fluorescent samples were composed of 2 l OGB-1
dye (O6806, Molecular probes) and 20 l Ca:EGTA
buffer solutions from the calibration kit (C3008MP,
Molecular Probes) with quoted free Ca2+ concentrations in the range of 17 nM to 39 M. Several fluorescence decay curves measured in the two-photon
excitation regime are shown in Fig. 1(c). The fluorescence lifetimes were obtained from the numerical
analysis of these data.
Thorough numerical fit must take into account the
IRF of the system [4]. As opposed to conventional systems with strongly asymmetric IRF, the Gaussianlike IRF of our system assumes a simple analytical
expression for the measured fluorescence decay. For a
train of excitation pulses of period T, each term of a
multiexponential decay process reads
Ik =
冋
1 1 + e−T/k
2 1 − e−T/k
冉
⫻exp −
t
k
冉冑
冊
− erf
+
2
IRF
2k2
IRF
,
2k
−
t
冑2IRF
k = 兵f,i,s其,
363
Wilms et al. [5] made the same observation for
OGB-1 fluorescence in the absence of Ca2+. However,
it was attributed to contaminating dye derivatives.
Previously reported data [1,3,4] for OGB-1 fluorescence lifetimes thus assume double-exponential decay. We argue that at high Ca2+ concentration, large
amplitude of long-leaving component make it difficult
to resolve short-lifetime components with small amplitudes, and the FLIM systems used in [1,3,4] exhibit insufficient temporal resolution.
Using triple-exponential fluorescence decay [Eq.
(1)] with the lifetimes independent of the calcium
buffer [1,4], we applied a global numerical analysis to
our data, yielding the partial intensities Ak in function of free Ca2+ concentration and the lifetimes
k (Fig. 2). Temporal resolution of our system allows
the short-lifetime component 共f ⬃ 188± 6 ps兲 to be
unambiguously resolved in the background of
the intermediate- and long-living fluorescence
(i ⬃ 768± 16 ps and s ⬃ 4.18± 0.01 ns). The partial
intensities Af and Ai decrease, while the slowcomponent intensity As = 1 − Af − Ai increases with
concentration 关Ca2+兴.
We attribute the short- and long-lifetime components to, respectively, unbound and Ca2+-bound
OGB-1 with a simple 1:1 complex stoichiometry. As
冊册
共1兲
where a triple-exponential decay is assumed and index k indicates fast 共f兲, intermediate 共i兲 and slow 共s兲
temporal components; k is the fluorescence emission
lifetime 共f ⬍ i ⬍ s兲; and erf共z兲 = 2 / 冑兰0z exp共−2兲d is
the error function. Equation (1) takes the periodic
train of excitation pulses and response time jitter of
the SPAD into account, such that our data in Fig. 1
do not require deconvolution processing.
Analysis of the OGB-1 fluorescence reveals the best
agreement with a triple-exponential decay approximation. The data are modeled using the function
AfIf + AiIi + 共1 − Af − Ai兲Is, with the fast If共t兲, intermediate Ii共t兲, and slow Is共t兲 decaying components [Eq. (1)]
of normalized partial intensities Af, Ai, and As = 1
− Af − Ai, respectively. Figure 1(c) details a comparison between the double-exponential (dotted curve)
and triple-exponential (dashed curve) models applied
to fluorescence from a 38 nM Ca2+ concentration
sample. The residues [Fig. 1(d)] do not show any bias
that might be caused by the Gaussian curve approximation of the IRF, while the quality of the numerical
fit is improved in case of the triple-exponential
model, as also confirmed by the reduced 2 values.
Fig. 2. (Color online) (a) Partial intensities As, Ai, and
Af of the slow 共s = 4.18± 0.01 ns兲, intermediate 共i
= 768± 16 ps兲, and fast 共f = 188± 6 ps兲 decay components
and (b) the mean lifetime ave = 兺kAkk of OGB-1 fluorescence as a function of the free Ca2+ concentration in comparison with the data from [3] and [4], 2r = 0.96. The numerical fit to the model [Eq. (2)] (curves) yields KD
= 195± 9 nM, KDn = 4 ± 6 M, and n = 9 ± 3. (c) Ratio
Af关Ca2+兴 / As (points) and its variance-weighted mean, yielding KD = 180± 11 nM (dashed line). The inset shows the
structure of OGB-1 molecule.
364
OPTICS LETTERS / Vol. 34, No. 3 / February 1, 2009
such, and because the intensities are defined by complex concentrations, we obtain that Af ⬀ 关D兴 and As
2+
⬀ 关CaD兴 = K−1
D 关Ca 兴关D兴, with [D] and [CaD] being the
unbound and Ca-bound OGB-1 concentrations and
KD the dissociation constant of the 1:1 complex. The
ratio Af关Ca2+兴 / As thus yields an estimate of the dissociation constant. In Fig. 2(c), a weighted-average
ratio 180± 11 nM is in good agreement with the
quoted (by manufacturer) KD of 170 nM.
In Fig. 2, the long-lifetime component As is in good
agreement with the measurements of Agronskaia
et al. [3] and Wilms et al. [4], but in [4], As ⬍ 1 in
Ca2+-saturated buffer and KD is too high (⬃300 nM).
These discrepancies are attributed to the dye impurities and to the difference in two-photon absorption
cross sections of Ca-bound and unbound OGB-1,
which are said to be too difficult to be quantified [4].
For the Ca-bound OGB-1, Lakowicz [1] reported
s of 4 ns, which agrees well with our data. In [3], s
varies in the range 2.6– 3.7 ns, while in [4],
s = 3.63 ns. For the Ca-free OGB-1, a single-lifetime
component f has been measured of 700 [1], 290–420
[3], and 346 ps [4], respectively. These results correspond to the combined effect of the decay processes
with the lifetimes f ⬃ 190 ps and i ⬃ 770 ps in our
measurements, which may not have been resolved in
previous studies. The experimental setup in [4] relies
on a commercially available photomultiplying tube
with a quoted timing jitter of 200 ps. Other data [1,3]
are reported without the timing resolution.
Owing to the Ca-binding features of the octadentate chelator BAPTA [8], the 1:1 stoichiometry is usually attributed to the Ca-bound OGB-1 as well. However, the molecule of OGB-1 has an asymmetric
structure with respect to its BAPTA moiety (Fig. 2,
inset). The asymmetric arrangement of the carboxyl
functional groups, benzol rings, fluorine, and nitrogen atoms as well as the variations of electronic density from low (at hydrogen in carboxyl groups) to high
(at benzol rings, F and O atoms) assume a polarity of
the molecule and, as a consequence, weak dipole–
dipole intermolecular forces. As a result of such interaction forces, several OGB-1 molecules might be coordinated to a Ca-bound OGB-1 thus forming a
polymolecular association with a calcium:indicator
molar ratio 1:n. The intermediate lifetime component
i in Fig. 2 then might be attributed to such polymo−1/n
lecular association, yielding Ai ⬀ 关Ca1/nD兴 = KDn
关D兴
2+ 1/n
⫻关Ca 兴 , with KDn being the corresponding dissociation constant. Note that formation of a polymolecular structure is a very complex processes and requires thorough investigations, but a non-1:1
stoichiometry of Ca:OGB-1 has been reported in [9];
here it allows us to build an accurate model:
关Ca2+兴
KD
=
As
Af
=
As
1 − As
冋 册
1+
Ai
Af
,
关Ca2+兴
KDn
=
冉冊
Ai
Af
n
.
共2兲
In the limit KDn → ⬁ (double-exponential decay), it
agrees with the Hill equation, as opposed to the
model in [4]. The numerical fit of data in Fig. 2 reports small binding affinity 共KDn ⬃ 4 M兲 and n = 9,
indicating the most probable form of polymolecular
association with eight Ca-free molecules coordinated
to the Ca-bound OGB-1. The dissociation constant KD
reported by the fit is 195 nM, in agreement with the
average of Af关Ca2+兴 / As (bottom panel).
The triple-exponential fluorescence decay in nonratiometric Ca probes as a result of dye contaminations
was suggested by Lakowicz et al. [10] for the Calcium
Green (CG-1). The CG-1 molecule has the same
structure as OGB-1, but the two F atoms are replaced
by Cl, yielding [10] f = 50 ps, i = 450 ps, and s
= 3.7 ns. Interestingly, the signature of non-1:1 stoichiometry has been noticed for Ca-bound CG as well
[11]. We find that at low Ca2+ concentration, the ratio
Af关Ca2+兴 / As in [10] is a constant of 170± 20 nM, yielding KD close to the quoted value of 190 nM. At high
Ca concentration, the bandpass of the system (2 GHz
modulation using a frequency-domain FLIM technique) was insufficient to accurately measure the
small component Af. The data in [10] can thus be accurately interpreted in the framework of our model
[Eq. (2)] without appealing to dye impurities. (KD obtained in [10] using a conventional model [12] is
128 nM).
In summary, we have shown that high-temporal
resolution measurements of the triple-exponential
fluorescence decay of nonratiometric Ca2+ indicators
allow the free Ca2+ concentration and dye dissociation constant to be measured precisely.
D. L. B. is grateful to Leonid Zekel and Edwin Constable for discussions on the model. This research
was supported by a grant of the Swiss National Science Foundation and by Centre des Systèms Intégrés
of Ecole Polytechnique Fédérale de Lausanne.
References
1. J. Lakowicz, Principles of Fluorescence Spectroscopy
(Plenum, 1999).
2. K. Suhling, P. M. W. French, and D. Phillips,
Photochem. Photobiol. Sci. 4, 13 (2005).
3. A. V. Agronskaia, L. Tertoolen, and H. C. Gerritsen, J.
Biomed. Opt. 9, 1230 (2004).
4. C. D. Wilms, H. Schnidt, and J. Eilers, Cell Calcium
40, 73 (2006).
5. C. Niclass, A. Rochas, P. A. Besse, and E. Charbon,
IEEE J. Solid-State Circuits 40, 1847 (2005).
6. W. Denk, J. H. Strickler, and W. W. Webb, Science 248,
73 (1990).
7. A. Habenicht, J. Hjelm, E. Mukhtar, F. Bergstrom, and
L. B. A. Johansson, Chem. Phys. Lett. 354, 367 (2002).
8. R. Y. Tsien, Biochemistry 19, 2396 (1980).
9. D. Thomas, S. C. Tovey, T. J. Collins, M. D. Bootman,
M. J. Berridge, and P. Lipp, Cell Calcium 28, 213
(2000).
10. J. R. Lakowicz, H. Szmacinski, and M. L. Johnson, J.
Fluoresc. 2, 47 (1992).
11. M. Eberhard and P. Erne, Biochem. Biophys. Res.
Commun. 180, 209 (1991).
12. G. Grynkiewicz, M. Poenie, and R. Y. Tsien, J. Biol.
Chem. 260, 3440 (1985).