This Article Abstract Full Text (PDF) All Versions of this Article: biophysj.104.055319v1 88/6/4319    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Kim, S. A. Articles by Schwille, P. Search for Related Content PubMed PubMed Citation Articles by Kim, S. A. Articles by Schwille, P.

Two-Photon Cross-Correlation Analysis of Intracellular Reactions with Variable Stoichiometry

Sally A. Kim ^* , Katrin G. Heinze ^*, Kirsten Bacia ^*, M. Neal Waxham  and Petra Schwille ^*

^* Institute of Biophysics, BioTec, Dresden University of Technology, Dresden, Germany; and Department of Neurobiology and Anatomy, University of Texas Health Science Center at Houston, Houston, Texas

Correspondence: Address reprint requests to Petra Schwille, Institute of Biophysics, BioTec, TU Dresden, Tatzberg 47-51, D-01307 Dresden, Germany. Tel.: 49-351-463-40328; Fax: 49-351-463-40342; E-mail: petra.schwille{at}biotec.tu-dresden.de.

	ABSTRACT

TOP ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
We successfully demonstrate the effectiveness of two-photon fluorescence cross-correlation spectroscopy (TPCCS) to study the complex binding stoichiometry of calmodulin (CaM) and Ca²⁺/CaM-dependent protein kinase II (CaMKII). Practical considerations are made for developing an intracellular cross-correlation assay, including characterization of the fluorescent molecules involved, calibration procedures of the setup, and optimal measurement conditions. Potential pitfalls and artifacts are discussed, and the complex stoichiometry of the molecular system is accounted for by a new experimental and theoretical framework for TPCCS. Our tailored model accommodates up to 12 red-labeled CaMs binding to a single green-labeled dodecameric CaMKII holoenzyme and accounts for the probability distributions of bound ligand as well as the respective changes in fluorescence emission upon binding. The model was experimentally demonstrated both in solution and in living cells by analyzing the binding of Alexa 633(C2)CaM to eGFP-CaMKII under different biochemical conditions known to induce the basal, activated, and autophosphorylated forms of the enzyme. Key binding parameters, such as binding degree, concentrations of reactants, and binding affinities, were determined under varying conditions with certain assumptions. TPCCS thus offers the unique ability to test our biochemical understanding of protein dynamics in the intracellular milieu.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
This study is devoted to performing quantitative analysis of binding interactions between Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) and calmodulin (CaM) using dual-color two-photon cross-correlation analysis (TPCCS), thereby demonstrating, for the first time, the full analytical power of this ultrasensitive technique in living cells. CaMKII is a ubiquitous mediator of Ca²⁺-dependent signaling and regulates cellular functions, such as synaptic plasticity, gene expression, and neurotransmission (Hudmon and Schulman, 2002). CaMKII assembles into multimeric holoenzymes containing 12 subunits (Kolodziej et al., 2000), each of which individually binds to its regulator, the Ca²⁺-receptor, CaM (Yamauchi, et al., 1989).

Each subunit can exist in multiple functional states that differ in their activity. In its basal state, CaMKII is kept inactive by the presence of an intrinsic autoinhibitory domain. Binding Ca²⁺/CaM to the autoinhibitory domain produces a fully active state, which allows the kinase to phosphorylate substrates as well as itself. Ca²⁺/CaM-dependent autophosphorylation of Thr²⁸⁶ (within the autoinhibitory domain of the subunit) produces a new state of the kinase in which CaM becomes essentially trapped, allowing it to remain fully active for many seconds even when Ca²⁺ levels have dropped below those initially required for activation (Meyer et al., 1992; Putkey and Waxham, 1996). Each of these states can be captured under different in vitro biochemical conditions. Since a range of kinase activity can be elicited by varying the CaM occupancy per holoenzyme, quantitating the degree of binding and binding affinities both in vitro and in cells would reveal information about CaMKII's input and output properties due to its oligomeric structure.

Dual-color fluorescence cross-correlation spectroscopy (FCCS) allows us access to these parameters by following molecular interactions between two differently labeled species. For the case of two binding partners, the presence of coincident signal fluctuations in two distinct emission channels indicates concomitant movement, which confirms binding between the two species. FCCS provides a viable alternative to the current standard of Förster resonance energy transfer (FRET) without distances or protein dimension constraints for the fluorescent labeling. Cross-correlation is also more versatile than FRET for issues of signaling, since it can measure interactions between molecules on different sides of the plasma membrane (Larson et al., 2004), co-localization of different cargos in endocytic vesicles (Bacia et al., 2002), and now binding of molecules with complex binding stoichiometry.

Here we apply fluorescence cross-correlation spectroscopy to understanding the binding dynamics of CaMKII and CaM both in solution and in living cells. For simplicity and stability of the setup along with biological benefits (Kim et al., 2004; Berland et al., 1995; Schwille et al., 1999), two-photon excitation (TPE) was utilized. Practical considerations for assay development, including characterization of the molecules involved, calibration of the setup, and optimal measurement conditions, are described along with warnings about potential pitfalls and artifacts. Moreover, we offer an extension of the basic cross-correlation theory to deal with complex binding stoichiometry (Weidemann et al., 2002). This concept is formalized into a new model for quantitative evaluation of cross-correlation experiments to determine key binding parameters, such as binding degree, total and free concentrations of reactants, and binding affinities, and experimentally demonstrated by analyzing the binding of multiple A633(C2)CaM to a multivalent eGFP-CaMKII holoenzyme both in vitro and in living cells.

	THEORY

TOP ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Theoretical concept of cross-correlation
The correlation function employed in fluorescence correlation spectroscopy (Eigen and Rigler, 1994; Elson and Magde, 1974; Magde, et al., 1972) is defined by

(1)

where F = F(t) – F and ... denotes a temporal mean.

Autocorrelation employs only one fluorescence signal (a = b) to analyze, for instance, diffusive behavior. Cross-correlation (a b; G_ab = G_x) looks for common features in two spectrally distinct emissions. By probing coincident attendance of distinct fluorophores in the volume element, it provides more specific analysis of molecular binding and unbinding than conventional autocorrelation analysis (Rarbach et al., 2001; Schwille, 2001).

The experimentally obtained FCS curves are commonly evaluated by nonlinear least squares fitting to analytical functions. For the simple case of a detection profile that can be approximated by a squared Gaussian distribution (for two-photon excitation) and for only one species of freely diffusing particles, the correlation curve is fitted to the following time-dependent correlation decay function M() multiplied by an amplitude G(0),

(2)

where r₀ and z₀ are the characteristic lateral and axial detection volume dimensions, and is defined as the average lateral diffusion time for a molecule with diffusion coefficient D. In the case of only one diffusing species, a molecular brightness parameter in kHz counts per molecule can be experimentally determined by multiplying the average photon count rate F by G(0). When correlation curves are obtained from mixtures of independently diffusing species i with different diffusion coefficients D_i (corresponding to different diffusion times _Di), fitting with linear combinations of M_i() is required.

Correlation curve amplitudes for populations with mixed species
The relative contributions of the M_i() and the overall amplitudes of the auto- and cross-correlation curves in the case of mixed species depend on the concentrations C_i and the molecular brightnesses _i of the species involved. For two detection channels, green g and red r, one must distinguish the brightness of a species i in the green channel, _i,g, from the brightness of that same species detected in the red channel, _i,r. Ideally, single-labeled species have a non-zero brightness only in one channel, and only double-labeled species exhibit non-zero brightnesses in both channels.

Generally, if the species i with concentrations C_i and brightnesses _i,g in the green channel and _i,r in the red channel contribute to the fluorescence fluctuations, the correlation amplitudes arise from the weighted sums of the different species contributions (see Eq. 17.11 in Schwille, 2001),

(3a)

(3b)

(3c)

where V_eff is the effective detection volume size,

(4)

which is determined from a calibration measurement.

The simple 1:1 binding case
For the simplest binding case, molecules of a green species G with molecular brightness _G (detected only in the green channel) and molecules of a red species R with molecular brightness _R (detected only in the red channel) engage in a binding reaction of 1:1 stoichiometry, forming the species GR. It is assumed that there is no spectral cross-talk for an idealized setup and no changes in brightnesses occur upon binding, so that the brightness of GR in the green channel is _G, and in the red channel _R. The autocorrelation amplitudes of the green and the red channel, G_g and G_r, and the cross-correlation amplitude G_x then become independent of molecular brightnesses,

(5a)

(5b)

(5c)

where C_GR is the concentration of bound molecules, C_R,t = C_R + C_GR is the total concentration of red-fluorescent molecules, and C_G,t = C_G + C_GR is the total concentration of green molecules.

In this standard cross-correlation assay, the maximum value that the cross-correlation amplitude can reach in the case of maximum binding equals the lower of the two autocorrelation amplitudes (Kohl, et al., 2002). This can be immediately verified from the expressions in Eq. 5, a–c, by setting the concentration of bound (C_GR) equal to the lower one of the total concentrations. For example, red excess (C_R,t > C_G,t) and full binding (C_GR = C_G,t) result in the red autocorrelation curve being lower than the green, G_r(0) < G_g(0), and the cross-correlation amplitude reaching the red autocorrelation amplitude, G_x(0)/G_r(0) = 1. It follows from the expressions in Eq. 5, a–c, that, in this simple cross-correlation assay, the absolute concentration of bound species is directly determined from the correlation amplitudes,

(6)

The 1:n full binding case
Here, however, we follow a reaction where a uniformly labeled partner A (eGFP-CaMKII oligomer) binds up to n single-labeled ligands B (Alexa 633(C2)CaM), forming complexes ranging between 1:1 and 1:n ligand stoichiometry (AB₁, AB₂, ... , AB_n). The number n of potential binding sites for our case is 12 (Kolodziej et al., 2000). It is assumed that binding is equivalent and independent (no cooperativity), since no physical model exists to support otherwise. We show that the cross-correlation amplitude can well exceed the lower autocorrelation amplitude obtained in this case.

We defer the treatment of all possible stoichiometries (arbitrary degrees of binding) to the next section and now consider only the limiting case of full binding. In this special case, the final reaction product contains two molecular species: 1), fully saturated bound species AB_n and 2), excess free ligand B. The contrary case, where the multi-binding-site molecule A would be in excess, cannot be treated in this simplified manner, since it does not give rise to a single, fully bound species. We have chosen the multi-binding-site partner A to be labeled with green fluorescence (G) and the ligand B to possess a red-fluorescent label (R). The rationale for this choice is that—despite the use of a spectrally well-separated dye pair—cross-talk artifacts could still hamper the analysis if the green ligand was in excess. However, these issues are conveniently minimized when the red-labeled ligand is in excess.

In the simple full binding case, C_ABn represents the concentration of the fully bound species and C_B the concentration of excess free ligand. The molecular brightnesses (in kHz/molecule) of the bound species AB_n is _G in the green channel and _Rn in the red channel. The brightness of the free ligand in the red channel is _R. According to the general expressions in Eq. 3, b and c, the autocorrelation amplitude in the red detection channel is then represented by

(7)

and the cross-correlation amplitude by

(8)

The relative cross-correlation then yields

(9)

When there is only fully bound AB_n with no free ligand B in excess, then and meaning that the relative cross-correlation amplitude equals 1, the same as in the simple 1:1 binding case. However, when the ligand is in large free excess (), then

(10)

Consequently, if the brightness _Rn of AB_n is larger than _R, the relative cross-correlation will actually exceed unity. In the absence of quenching or FRET, _Rn = n · _R; and the relative cross-correlation directly reflects the number of binding sites,

(11)

When we consider that the red fluorescence emission changes upon binding, a correction factor q_R must be introduced, so that _Rn = q_R · n · _R. The relative cross-correlation amplitude then becomes

(12)

If the fluorescence emission is reduced by quenching upon binding, then q_R < 1 and the relative cross-correlation amplitude will be less than n. On the other hand, if the emission is enhanced by Förster resonance energy transfer (FRET), then q_R > 1 and the relative cross-correlation amplitude may even increase beyond n.

The simple limiting full-binding case can be experimentally realized if the concentrations of the reactants can be chosen at will (i.e., as in an in vitro experiment) and if the K_d is low enough to achieve full binding within the limits of detection (Fig. 4 B). Other cases, however, require a more general treatment that accounts for incomplete binding.

View larger version (30K): [in this window] [in a new window]   FIGURE 4  Monitoring complex binding stoichiometry of eGFP-CaMKII and A633(C2)CaM. Steady-state dual-color TPCCS measurements were obtained from the same reaction of A633(C2)CaM and eGFP-CaMKII in buffer solution under elevated Ca2+ (A), elevated Ca2+ and MgATP (B), and Ca2+-free (C) conditions. In the top panels, autocorrelation curves of the higher concentration of ligand, A633CaM (red), and respective cross-correlation curves (black) are shown; in the middle panel the corresponding autocorrelation curves of eGFP-CaMKII are plotted. The bottom panels show the binding distribution for each case, the average degree of binding , and the dissociation constant Kd derived from the amplitudes of the measured correlation curves above assuming high affinity binding (Kd < 200 nM, see text) for Case (B). *Assumption based on literature value Kd <0.02 nM (Meyer et al., 1992).

(13)

where k represents the single site binding constant k = 1/K_d, and C_B is the concentration of free ligand. Considering that the total concentration of all free and complexed A, C_A,t, is given by

(14)

C_A can be eliminated from Eq. 13 and rewritten as

(15)

These expressions for the concentrations of all fluorescent species can be entered into the correlation amplitude terms (the expressions in Eq. 3, a–c) for the general case.

In addition to the concentrations, the correlation amplitude terms require the molecular brightnesses _i,g and _i,r. The brightnesses of the free species in the absence of crosstalk are given by

(16)

If no changes in fluorescence emission occur upon binding, the brightnesses of the bound species are

(17)

Changes in molecular brightnesses upon binding
If changes in fluorescence emission occur upon binding due to quenching or FRET, Eq. 17 is not valid and the model has to be adapted. We express the change in the fluorescence emission as the ratio q = */ , where is the fluorescence of the free molecule that changes to * upon binding. Assuming independent and equivalent binding for all CaMKII subunits and quenching of both CaM and the CaMKII subunit upon binding to be independent from the binding state of neighboring subunits, the brightnesses of the complexes AB_i in the red and green channels become

(18)

(19)

where the unquenched molecular brightnesses are represented by (B,r) = _R and (A,g) = _G.

Application of the concentration distribution and varying brightnesses to the correlation amplitudes
The three correlation amplitudes, defined by the expressions in Eq. 3, a–c, can now be rewritten in terms of the newly defined concentrations and brightnesses. The concentrations c_i in the expressions in Eq. 3, a–c, are replaced by the concentrations of the different bound species C_i (from Eq. 15), the free concentration of A, C₀ (also from Eq. 15), and the concentration of free B, C_B. The brightnesses _i,g and _i,r are replaced by Eqs. 16, 18, and 19. In doing so, the absolute brightnesses _G and _R cancel out, leaving only the quench factors

(20a)

(20b)

(20c)

where the expressions of C_i and Q_i are given by Eqs. 15 and 19, respectively.

With the experimentally determined correlation amplitudes, the resulting system of three equations can be solved numerically, providing distinct solutions for the parameters C_B, C_A,t, and k = 1/K_d if the stoichiometry n and the quench factors q_G and q_R are known. For all experimental data, these equations gave rise to one unique real solution using Mathematica 5.0 (Wolfram Research, Champaign, IL). Additionally, the concentration of total bound ligand can be determined as and used to calculate the binding degrees. The overall binding degree is defined as = C_B,bound/C_A,t and the per-site binding degree as = /n.

Reduced number of available binding sites
One factor that could interfere with the interpretation of the derived results is the possibility of a reduced number of sites ñ among the potential total binding sites n = 12 due to steric hindrance of the fluorescent labels (e.g., only every second site can bind: ñ = 6). In this case, the summations are performed from i = 0 to ñ instead of i = n (the expressions in Eq. 20, a–c). The apparent green quench factor determined from the EGTA and MgATP cases is represented as

then Q_i is redefined as

(21)

and inserted in the expressions for Eq. 20, a and c.

Analysis of the intracellular case: incorporation of unlabeled ligand
The interpretation of intracellular cross-correlation experiments has the added complication of the presence of endogenous, unlabeled proteins that mix with the labeled proteins of interest. For CaMKII, we assume that the effect is negligible for HEK293 cells, since a brain-specific isoform of the enzyme is used, although endogenous - and -isoforms of CaMKII may be present that can form oligomers with the -subunit. On the other hand, CaM is a highly conserved, ubiquitous Ca²⁺-binding protein (van Eldik and Watterson, 1998), and the concentration of endogenous protein cannot be ignored. Therefore, the ratio of labeled to unlabeled CaM must be taken into account. Binding affinities of labeled and unlabeled CaM are assumed to be indistinguishable, for the sake of simplicity.

To accommodate unlabeled ligand B, we consider the concentration of the oligomer A with i ligands B bound, among which j ligands are red-labeled and i–j are unlabeled (i j). This concentration _iC_j is characterized by a binomial distribution of the red ligands among the i ligands bound to A,

(22)

where r is the overall fraction of red-labeled ligand molecules and C_i is taken from Eq. 15. The concentration of A with exactly j red-labeled ligand bound and any number of additional unlabeled ligand bound is obtained by summing up over all A with at least i = j ligands bound,

(23)

Substituting Eq. 22 into Eq. 23 yields

(24)

which accounts for the intracellular scenario where there is a mixture of labeled and unlabeled ligands. The concentrations of A that have j = 0 to n of fluorescent B bound, C_j, are applied to the expressions in Eq. 3, a–c, together with the concentration of free, fluorescent B, which is now r · C_B. The correlation amplitudes obtained in this manner replace the expressions in Eq. 20, a–c, for the in vivo case.

Note that the assumption for the size of the unlabeled fraction influences the exact distribution of the red-labeled ligands C_j or C_i obtained in the calculation (i.e., Eqs. 24 and 15 are not identical unless r = 1). The reason is that the addition of unlabeled CaM influences the reaction.

Effect of unlabeled ligand on the limiting 1:n full binding case
We have now derived how the presence of a fraction r of fluorescently labeled molecules among the ligands enters the general expressions in the binding theory for any degree of binding. For the special limiting case of full binding with excess CaM described previously, it can be directly illustrated how the mixing with unlabeled CaM reduces the cross-correlation. At full binding, Eq. 24 simplifies to a binomial distribution of red-labeled CaM on all green CaMKII oligomers,

(25)

The correlation amplitudes in the expressions in Eq. 20, b and c, are re-derived, this time with the simplified expression from Eq. 25 for C_j and again with r · C_B as the concentration of free, fluorescent B. Considering that, for a binomial distribution,

and the autocorrelation amplitude of the ligand B in the red detection channel is given by

(26)

the cross-correlation amplitude by

(27)

and the ratio by

(28)

where potential quenching of the green labeled A is not considered. The value is as previously defined (Eq. 9). When the ligand is in large excess (), Eq. 28 reduces to

(29)

For this simple limiting case of full binding with excess ligand, the result could have been predicted from Eq. 10, which states that the relative cross-correlation approaches _Rn/_R. The average brightness of the fully bound species in the red channel is now _Rn = r · q_R · n · _R, and the relative cross-correlation amplitude becomes modified not only by the quench factor q_R, but also by the factor of fractional labeling r.

Nonetheless, in the general case, quenching effects and the presence of unlabeled CaM have to be treated separately, since each affects the fluorescence of the bound CaM and the fluorescence of the free CaM differently. However, the results converge to the analogous appearance of the factors q_R and r in Eq. 29 for the limiting case of full binding with large ligand excess.

To reiterate, the effect of unlabeled, endogenous CaM is significant for intracellular measurements. If ignored, the binding affinity will be underestimated from the reduced cross-correlation amplitude.

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Materials
ATP was purchased from Amersham Biosciences (Piscataway, NJ). All chemicals were of reagent grade or higher and were purchased from Sigma-Aldrich (St. Louis, MO) unless otherwise specified. All fluorescent probes were obtained from Molecular Probes (Eugene, OR). The pEGFP-C2 vector was purchased from Clontech (Palo Alto, CA).

Confocal TPCCS instrumentation
The utilized confocal setup is identical to the setup described previously (Kim et al., 2004) and consists of a tunable mode-locked Titanium:Sapphire laser (Tsunami, Spectra Physics, Darmstadt, Germany) coupled to an inverted microscope (IX-71 Olympus, Hamburg, Germany) which was equipped with a water immersion objective (UPlanApo/IR 60x, NA 1.2, Olympus). Various dichroic mirrors and bandpass emission filters (AHF Analysentechnik AG, Tübingen, Germany) were used for coupling the laser beam into the objective as well as for splitting the fluorescence into the green and red detection channels (see Kim et al., 2004, for details). For digital correlation of the detected signal, either of two hardware correlators (ALV-5000 multiple- correlator card, ALV, Langen, Germany) for simultaneous measurement of the dual auto- and cross-correlations or a more advanced cross-correlator (Flex02, Correlator.com, Bridgewater, NJ) simultaneously operating in auto- and cross-correlation mode was used. Evaluation of the curves was carried out using ORIGIN (MicroCal Software, Northampton, MA) with a Levenberg-Marquardt nonlinear least-square fitting routine.

The size of the effective (open) volume element was calculated as 0.26 ± 0.04 femtoliters based on calibration measurements with pure Rhodamine Green (RG) dye excited at 920 nm. Optimal alignment is ensured by a calibration procedure recommended by Schwille et al. (1997) to confirm the same size and position in space of both detection volumes of the dual channel system.

In vitro and intracellular TPCCS experimental conditions
All solutions were prepared in low fluorescence HPLC grade water (Merck Sharp & Dohme, Haar, Germany) and were filter-sterilized. Twenty-four microwell chambers (Evotec OAI, Hamburg, Germany) were pre-incubated with 10% BSA for 10 min followed by two subsequent washes with an enzyme dilution buffer (25 mM HEPES, pH 7.0, 150 mM KCl, 0.1 mg/ml BSA). Reaction conditions for cross-correlation measurements included: A633(C2)CaM and eGFP-CaMKII in 25 mM HEPES, pH 7.0, 100 mM KCl, 10 mM Mg²⁺ with 0.5 mM CaCl₂, and 1 mM MgATP followed by 1 mM EGTA. (Note that A633(C2)CaM was the same construct as previously described in Kim et al., 2004; however, for A633(C3)CaM), the notation has been modified to reflect the conventional amino acid numbering found in the CaM literature.) Each reaction was well mixed before each measurement, and measurements were performed at room temperature. The data recording time for each curve was 100 s.

HEK293 cells were transfected with eGFP-CaM-kinase II using Effectene (Qiagen, Hilden, Germany) as described by the manufacturer. A633(C2)CaM was introduced into adherent cells using electroporation as previously described (Kim et al., 2004). For some experiments cells were incubated in one of the following conditions: HBSS (Kim et al., 2004) with no added Ca²⁺, 15 µg/ml -hemolysin (-toxin; Sigma, St. Louis, MO), and 200 µM EGTA (Ca²⁺-free); or HBSS with 15 µg/ml -hemolysin, 10 mM CaCl₂, and 1 mM MgATP (elevated Ca²⁺). The -hemolysin permeabilizes the cell membrane allowing rapid equilibrium of extracellular ions and nucleotides with the intracellular environment. HEK293 cells were prepared in parallel for spectroscopy measurements and Ca²⁺ imaging. Fluo-4 was used to confirm that Ca²⁺ was being elevated or decreased uniformly in all cells under the different conditions.

All experiments were performed at room temperature. The detailed methods including construction, expression, purification, and labeling strategies for the proteins, cell culture, and transfection are described in Kim et al. (2004).

CaM-kinase II activation assay
The activity of recombinant CaM-kinase II and eGFP-CaM-kinase II was determined in the presence of increasing concentrations of recombinant bacterial-expressed CaM (wtCaM) or A633(C2)CaM. Serial dilutions of the CaM stocks were made in 5 mM MOPS buffer with 0.1 mg/ml BSA. CaM at the final concentrations indicated in the figure were added to reaction mixtures (final volume = 25 µl) containing: 25 mM HEPES pH 7.0, 70 mM KCl, 0.5 mM CaCl₂, 10 mM MgCl₂, 100 µM ATP, 2 µCi of [-32P]ATP, 50 µM syntide (PLRRTLSVAA), and 0.4 mM dithiothreitol. CaM-kinase II was diluted in 10 mM MOPS pH 7.0, 200 mM KCl, 0.1% T-20, and 1.0 mg/ml BSA. Reactions were prewarmed for 1 min at 30°C and were initiated by the addition of enzyme (10 ng). Each reaction was incubated for 30 s at 30°C before 20 µl of the reaction mix was spotted onto Whatman P-81 filter paper. The filters were processed as previously described (extensively washed with phosphoric acid and ethanol and dried), and bound radioactivity was quantified by Cerenkov counting (Kim et al., 2001).

Activity is expressed as µmol/min/mg. For K_act and V_max determinations, activity data were plotted versus CaM concentration and fit using nonlinear least-squares to the hyperbolic Michaelis-Menten equation,

(30)

where K_act is the concentration of CaM which yielded 50% of maximal activity, and V_max is the maximal rate of activity of the enzyme.

Excitation and emission spectra
Fluorescence excitation and emission maxima for all fluorescently labeled CaMs were accomplished at room temperature (22°C) using a PTI QuantaMaster steady-state spectrofluorometer (Model CGM-9/2003; Photon Technology International, Lawrenceville, NJ). The excitation and emission slit widths were 2 and 5 nm, respectively. The standard buffer for all fluorescence measurements was 50 mM HEPES, pH 7.0, 100 mM KCl, and 0.1 mM EGTA. A background excitation and emission scan was measured before addition of the probe. Fluorescently labeled CaM (200–500 nM as assessed by a protein assay) was added to 1 ml of the standard buffer to achieve an appropriate fluorescent signal. To test calcium sensitivity, samples were first scanned in the base buffer and rescanned after the addition of 1 M CaCl₂ for a 5 mM final calcium concentration. At the end of every run 5 µl of 1 M EGTA was added. To determine pH-dependent changes in fluorescence, separate samples of the protein were tested in the standard spectral buffer at pH 6.5 and pH 8.0.

	RESULTS

TOP ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Establishing an intracellular two-photon cross-correlation assay
Optimizing measurement conditions
Our selected dye system consists of enhanced Green Fluorescent Protein (eGFP, green) coupled to CaMKII and A633 (red) labeling CaM. The emission maxima in aqueous solution using TPE are 510 nm for eGFP and 655 nm for A633. As expected, the emission spectra (Kim et al., 2004) are not notably different compared to their one-photon excited counterparts and are very well separable by optical filter systems.

An optimization procedure of the excitation parameters, wavelength, and intensity, was performed to find an initial compromise of optical conditions for which the dyes exhibit similar performance. First, the optimal common two-photon excitation wavelength was found by scanning with the tunable laser between 800 and 980 nm, and the characteristic molecular brightnesses (kHz/molecule) at fluorescence saturation were measured with FCS. During the scan, the laser intensity was adjusted to the respective saturation level for maximal efficiency, while the pulse-width was kept constant at 100 fs. Fig. 1 A shows a plot of versus the TPE wavelength for the current selected fluorophores eGFP (circles) and A633 (squares) measured independently in buffer solution for maximal resolution. All data were recorded under optimal conditions in a standard single color setup with a broadband emission filter and without specific filter separation for red and green. In Fig. 1 A the curves cross at 860 nm, indicating the optimal excitation wavelength. However, since CaMKII is a holoenzyme, multiple eGFPs were incorporated per molecule, and thus, for our intracellular dual-color system, the "green" was found to be 6–7 times brighter than wild-type eGFP. Taking this into consideration along with any losses due to the optical setup, the optimal excitation wavelength would be 830 nm.

View larger version (35K): [in this window] [in a new window]   FIGURE 1  Optimization of measurement conditions. (A) Wavelength dependence of the fluorescence emission yield for eGFP (circles) and Alexa 633 (squares). Data were measured in a standard two-photon FCS setup under optimal conditions with the pulse width of the laser set at 100 fs. The excitation power was adjusted to the respective fluorescence saturation limit for each fluorophore independently. (Bottom) Plot of fluorescence emission dependence on excitation power for eGFP (circle) and A633 (square) at 830 nm (B) and 920 nm (C). The fluorophores were independently analyzed in aqueous solution in a standard FCS setup. At low power levels (solid symbols) the slope is 2. The boxes (dashed line) mark the range of the respective possible excitation power for each individual fluorophore. Contrary to B where these areas show no overlap, C shows "safe" excitation powers where common excitation can be found.

Accordingly, the optimal wavelength for subsequent cross-correlation experiments at which both dyes were well-excited, while the photobleaching and damage within the cells were minimized, was found to be 920 nm. For longer wavelengths, the A633 emission efficiency decreases, whereas for shorter wavelengths, its saturation levels become too low to work in combination with eGFP.

Effects of photobleaching artifacts on cross-correlation experiments
Next, we used the labeled proteins and bound complex to fine-tune measurement conditions for further experimentation. Fig. 2 highlights various aspects of how photobleaching and saturation affects a cross-correlation experiment. For excitation powers up to 1 MW/cm² the cross-correlation amplitude (Fig. 2, C2), as well as the particle number and the diffusion time of the single species eGFP-CaMKII (Fig. 2, A2) and A633(C2)CaM (Fig. 2, B2), remain constant; however, for higher laser power densities, a typical reduction of the diffusion time (Fig. 2, A4 and B4) and increased apparent particle number (Fig. 2, A3 and B3) are observed.

View larger version (30K): [in this window] [in a new window]   FIGURE 2  Controlling photobleaching artifacts in two-photon cross-correlation spectroscopy experiments. To assess the optimal intensity where power-dependent bleaching and saturation artifacts were not introduced, cross-correlation and autocorrelation curves were measured for a high affinity binding reaction (10 mM Ca2+/1 mM MgATP) of eGFP-CaMKII (A) and A633(C2)CaM (B) at power densities between 0.5 and 3.5 MW/cm2. The data acquisition time for each curve was 60 s. Although the count rate per molecule that determines the quality of the measurement (A1, B1) increases with increasing excitation power, a decrease of the cross-correlation amplitude (C2) due to photobleaching is observed for excitation powers >1 MW/cm2, which can be easily mistaken for decreased binding (C1). The corresponding autocorrelation curves of the single species eGFP-CaMKII (A2) and A633(C2)CaM (B2) show an increase of the apparent particle number (A3, B3) and a decrease of apparent diffusion time (A4, B4), which indicate saturation and bleaching artifacts.

Effects of labeling the proteins
It is essential to characterize the consequences of the dye labeling for the biological and biochemical function of the proteins under investigation. Since CaM does not have any resident Cys residues in its sequence, a single Cys residue was introduced for Asp at amino acid 2 for Alexa 633 labeling. A comparison of the crystal and NMR solution structures of free CaM (Babu et al., 1985) and CaM bound to a peptide derived from CaMKII (Meador et al., 1993) indicated that the side chain of Asp2 does not directly participate in peptide binding. Enhanced green fluorescent protein (eGFP) was attached to CaMKII on the N-terminal end of the protein to avoid interfering with the association domain on the C-terminal end, which is critical for holoenzyme formation (Kolb et al., 1998).

To assure that the A633 dye label on CaM did not seriously affect the activation of CaMKII, labeled and unlabeled CaM were tested for their ability to activate recombinant CaMKII and eGFP-CaMKII (Fig. 3). Fig. 3 shows the activity of CaMKII determined in the presence of increasing concentrations of unlabeled CaM (solid squares) and A633(C2)CaM (shaded circles). The lines represent a fit of the experimental data to the Michaelis-Menten equation (Eq. 30). The value K_act is the concentration of CaM that yielded 50% of maximal activity, and V_max is the maximal activation of the enzyme. The K_act and V_max values were obtained for each protein combination by fitting the data using the Levenberg-Marquardt nonlinear least-squares fit. The K_act values for CaM are within the range consistent with those reported under saturating Ca²⁺ conditions (Kuret and Schulman, 1984; Brickey et al., 1990; Katoh and Fujisawa, 1991). No significant statistical differences (p < 0.05) were seen between CaMKII and eGFP-CaMKII in terms of either V_max or K_act (Table 1). However, although no difference in the K_act was apparent, Alexa633-labeled CaM resulted in a statistically significant decrease in V_max for both unlabeled CaMKII (32% reduction; solid triangles) and eGFP-CaMKII (18% reduction; shaded triangles). A full characterization of eGFP-CaMKII showed that the protein was a holoenzyme and possessed minimal basal-independent activity similar to the wild-type enzyme (data not shown).

View larger version (19K): [in this window] [in a new window]   FIGURE 3  Characterizing labeling effects on CaM and CaMKII. The activity of CaMKII and eGFP-CaMKII was determined in the presence of increasing concentrations of wild-type CaM and A633(C2)CaM as described in Materials and Methods. CaM concentrations were varied from 3.13 nM to 3.2 µM. Activity is expressed as µmol/min/mg. Each data point represents a mean of three independent measurements (N = 3). The lines represent a fit of the experimental data to the hyperbolic form of the Michaelis-Menten equation (Eq. 30). The Kact and Vmax values were obtained for each protein combination by fitting the data using the Levenberg-Marquardt nonlinear least-squares fit (see Table 1 for values). No statistically significant differences were seen between CaMKII (solid squares, solid curve) and eGFP-CaMKII (shaded circles, shaded curve) when activated with wild-type CaM in terms of both Vmax and Kact. Although no difference in the Kact was apparent, A633(C2)CaM resulted in a statistically significant decrease in Vmax for both unlabeled CaMKII (solid triangles, dashed solid curve) and eGFP-CaMKII (shaded triangles, dashed shaded curve).

Cross-correlation analysis facing variable binding stoichiometry
Diffusional analysis of CaM and CaMKII
Fig. 4 shows the steady-state TPCCS measurements for a reaction of A633(C2)CaM and eGFP-CaMKII in buffer solution under different reaction conditions. The autocorrelation curves and corresponding fits of the higher concentration of ligand, A633(C2)CaM (in red), and respective cross-correlations (in black), are seen in the upper panel, whereas the corresponding autocorrelation curves and fits of eGFP-CaMKII are plotted in the middle panel. No term for triplet or other photophysical dynamics was needed in the fast timescales for fitting because of using TPE in conjunction with low laser intensities (Dittrich and Schwille, 2001). Table 2 provides a summary of diffusion coefficients for reactions seen in Fig. 4.

In vitro cross-correlation measurements
Since the biochemistry of CaMKII is well characterized, specific in vitro conditions are known to induce different functional states of the kinase with differential binding of CaM. These in vitro measurements were used to determine the efficiency of the assay and stability of the setup. Simultaneous auto- and cross-correlation measurements were performed to investigate the activated, autophosphorylated, and basal states of CaMKII by varying the Ca²⁺ and ATP concentrations while maintaining a roughly constant excess of A633(C2)CaM (800–950 nM) compared to the eGFP-CaMKII concentration (4–5 nM of holoenzymes or 48–60 nM of subunits) within the same reaction.

Fig. 4 shows the correlation curves from one set of such measurements taken in succession from the same sample under different conditions. The activated state was first monitored under the initial conditions of 0.5 mM Ca²⁺ in the absence of MgATP (A). Subsequent addition of MgATP (B; 1 mM final concentration) induced a large increase in the cross-correlation indicative of the binding of multiple CaMs in the autophosphorylated state. Finally, addition of 2 mM EGTA (C) at the end of the reaction chelated the Ca²⁺ (forming apo-CaM), thereby dissociating the CaM and returning the enzyme back to the unbound state. Individual autocorrelation curves varied <5% with a slight decrease in concentration due to a small dilution effect of adding Ca²⁺, MgATP, and EGTA solutions.

The autophosphorylated state of CaMKII and determination of the red quench factor
Incubation of CaMKII with Ca²⁺/CaM in the presence of MgATP (B) leads to a virtually irreversible binding of CaM upon autophosphorylation of Thr²⁸⁶ (Meyer et al, 1992). The binding of Ca²⁺/CaM to CaMKII subunits has been found to be 1:1 (Katoh and Fujisawa, 1991; Sugiura and Yamauchi, 1993). For the concentrations employed here in the in vitro experiments, the limiting case, exemplified in the theory where only saturated CaMKII and a large excess of CaM are present, is applicable. If no changes in red fluorescence occur upon binding, the relative cross-correlation amplitude G_x(0)/G_r(0) (ratio of cross-correlation and red autocorrelation) yields the number of binding sites on the target molecule, which is n = 12 for CaMKII. The fact that a smaller value of G_x(0)/G_r(0) = 2.6 is observed can be attributed to the quenching of the red CaM fluorescence upon binding to CaMKII. According to the theory, the quench factor q_r reduces the maximum cross-correlation amplitude as Based on the autophosphorylated state (B), the red quench factor is therefore estimated to be q_r = 0.2, meaning that a labeled CaM exhibits 80% less fluorescence upon binding. In the case presented here, determination of q_r is independent of knowledge about the quenching of the green fluorescence upon CaM binding. Under these conditions where all CaMKII molecules are fully saturated with CaM, the green brightness can be considered to be essentially uniform and the actual value of the brightness cancels out.

For a more accurate evaluation of the red quench factor, the general binding theory was applied, and solutions were calculated. Inserting the experimental values for the three correlation amplitudes, G_g, G_r, and G_x, the number of potential binding sites (n = 12), and the size of the detection volume (V_eff = 0.26 fl), one obtains a solution for the K_d if the red quench factor q_r is known, or vice versa, a solution for the quench factor q_r if the K_d is known. Since in our experiments, this quench factor is not known a priori, we therefore determine the q_r by using a known value for the K_d in Case B from the well-documented literature. This yields a quench factor, which can then be applied to the other cases (A and C) to obtain experimental values for the K_d in those cases.

With respect to the case of autophosphorylated CaMKII (Case B), several studies have reported that it has one of the highest affinities to date for Ca²⁺/CaM binding with a K_d in the range of 150 nM down to low pM (Sugiura and Yamauchi, 1993; Meyer et al, 1992; Török et al, 2001; Gaertner et al., 2004). Using the middle value in the range (K_d < 0.02 nM; Meyer et al., 1992), we evaluate q_r to be 0.22 in the autophosphorylated state. Since the exact value of the K_d was not independently determined in our experiments, the quenching factor q_r was then calculated for the range of K_d values in the literature noted above using the expressions in Eq. 20, a–c. The quench factor was found to remain essentially constant within the entire range of reported K_d values (data not shown). Furthermore, K_d values for the other cases were also found to be constant for the q_r values obtained within the reported K_d range for Case B (data not shown).

The basal state of CaMKII and determination of the green quench factor
CaMKII does not significantly bind to CaM in the absence of Ca²⁺ due to the presence of an intrinsic autoinhibitory region in its regulatory domain (Hudmon and Schulman, 2002). Addition of excess EGTA (C) at the end of the reaction serves to chelate the Ca²⁺, dissociating the bound CaM. The brightness in kHz per molecule was determined from the product of the green autocorrelation amplitude and the green fluorescence count rate. The green quench factor q_g was quantified as the ratio of the brightnesses of the CaMKII holoenzyme in the autophosphorylated fully bound state (_g*) described above (Fig. 4 B, bottom graph) and the basal unbound state (_g) and was found to be q_g = _g*/_g = 0.48. Analysis of CaMKII and apo-CaM demonstrate that binding in the absence of Ca²⁺ is very weak with a K_d of 15 µM or more.

The activated state of CaMKII and determination of binding parameters
Each subunit of CaMKII becomes activated when CaM in the presence of Ca²⁺ binds to and displaces the autoinhibitory domain from its active site. For this case using the quench factors determined above and the general 1:n binding theory, the correlation data was evaluated to determine binding parameters of the activated state of CaMKII (A). The analysis yielded the following results: the binding distribution as depicted in the bottom graph of Fig. 4 C; a binding degree of = 0.28 CaM per binding site, or = 3.4 CaM per CaMKII holoenzyme; and a single site dissociation constant of K_d = 2.4 µM. The activated state of CaMKII, as expected, is therefore intermediate to the basal (C) and autophosphorylated (B) cases.

Potential deviation from the ideal case: reduced number of subunits available for binding
Although care was taken to minimize interactions due to labeling of the protein, it cannot be excluded that the binding stoichiometry of n = 12 is somehow affected by the eGFP fusion or other steric interference with the labeled CaM. Therefore the effect of assuming a smaller n (i.e., ñ = 9 or ñ = 6) was also analyzed. Although, of course, the overall binding degree is affected by the assumption for ñ, the (single site) dissociation constant, K_d, and the per-site binding degree, , were seen to remain largely unaffected. If we assume, for instance, ñ = 6 (only every second subunit is available for binding), then the results for these parameters of interest do not show any significant difference compared to the ideal case (ñ = n = 12): K_d (Case A) = 2.2 µM (instead of 2.4 µM); (Case A) = 0.29 (instead of 0.28); K_d (Case C) > 14.5 µM (instead of 15.2 µM); and (Case C) < 0.054 (instead of 0.052).

Cross-correlation analysis extended to the intracellular milieu
Assessing intracellular co-localization of CaM and CaMKII
The confocal laser scanning images (Fig. 5, top) show stably transfected eGFP-CaMKII HEK293 cells electroporated with A633(C2)CaM. Due to the variability in eGFP-CaMKII expression and the gradient of concentrations from the electroporation of A633(C2)CaM, different combinations of protein concentrations were found in the same dish. Higher concentrations of both proteins were used for confocal imaging, whereas lower concentrations were used for cross-correlation measurements (nM to low µM range). Permeabilization of the cell membrane was accomplished using -hemolysin, which allowed for equilibrium of extracellular ions and nucleotides with the intracellular environment. The images show no meaningful change in appearance under elevated 10 mM Ca²⁺ and 10 mM Ca²⁺/1 mM MgATP (Fig. 5, top left and center) in the presence of 15 µg/ml -hemolysin followed by Ca²⁺-free conditions (200 µM EGTA) (Fig. 5, top right) in the same dish. In all images Alexa 633(C2)CaM was localized to the cytoplasm with a higher concentration in the nucleus, whereas eGFP-CaMKII was generally cytoplasmic with very little signal in the nucleus. After longer periods of time in elevated Ca²⁺, eGFP-CaMKII often displayed small punctae throughout the cytoplasm (data not shown). These areas were avoided during cross-correlation measurements. Although the molecules were co-localized in the cytoplasm in both conditions, cross-correlation analyses show very different binding dynamics.

View larger version (55K): [in this window] [in a new window]   FIGURE 5  Binding of A633(C2)CaM and eGFP-CaM-kinase II under varying conditions in living cells in culture. (Top) Confocal images of stably transfected eGFP-CaM-kinase II HEK293 cells electroporated with A633(C2)CaM are shown under 10 mM Ca2+ (left) followed by 10 mM Ca2+/1 mM MgATP (middle) and 200 mM EGTA (right) in the presence of 15 mg/ml -hemolysin in the same dish (scale bar = 10 µm). (Middle) Intracellular auto- and cross-correlation curves were measured simultaneously under each of the conditions but at lower intracellular protein concentrations than used for the LSM images above. Plotted are the six subsequent 10-s acquisitions (thin lines) for both the auto- and cross-correlation measurements and their respective average curves with data fits (corresponding thick lines). (Note that the y axis for the green autocorrelation in A, right side in green, has a different scale than the other curves.) The amplitude of the cross-correlation curve (black) increases under elevated Ca2+ conditions (A, B) and almost disappears in the absence of Ca2+ (C). In comparison to the in vitro case, the cross-correlation amplitudes for both elevated Ca2+ conditions (A, B) are reduced due to the additional presence of unlabeled endogenous CaM in the cell. Distribution plots are shown at the bottom as seen in Fig. 4. The hatched shaded bars represent the distribution of CaM irrespective of the label, and the red bars denote the distribution of Alexa 633-labeled CaM molecules only. In A and B, an assumption of a high binding affinity yields virtually full binding (hatched shaded bar) and the indicated values for the labeled fractions r. In C the hatched bars exemplify a CaM distribution for r = 0.3 based on A and B. Distributions for conditions with other r fractions are also denoted (shaded lines).

Deviation from the in vitro case: mixing of Alexa 633(C2)CaM with endogenous unlabeled CaM
Quantitatively, the cross-correlation amplitudes obtained in the cells were always smaller than in vitro for comparable conditions. As a consequence, analysis of the intracellular data using the general formalism in the same way as for the in vitro data yielded the binding distributions depicted as red bars in Fig. 5 (bottom) and unexpectedly weak binding constants (K_d = 14 µM in the calcium case and K_d = 160 nM in the MgATP case, when using the same quench factors as in vitro). Although the binding and photophysical properties of the proteins may differ in living cells from the in vitro case, it seems unlikely that the binding is so much weaker or that quenching of the red fluorescence is so much stronger (q_r = 0.06) causing such a reduced cross-correlation. A more likely explanation is that the mixing of labeled CaM with the ample endogenous unlabeled CaM in the cell is responsible for the reduced cross-correlation amplitudes and an incorrect value for K_d. We therefore devised a theoretical approach that incorporates the presence of unlabeled CaM (see Theory).

Extracting information about intracellular binding interactions between CaM and CaMKII
The intracellular case of elevated Ca²⁺ and MgATP (B) where CaMKII should be autophosphorylated was analyzed with the aim of studying the dependency between the K_d and the fraction r of labeled CaM, following the same routine used for obtaining the relationship between K_d and q_r from the in vitro data. The essential result is that r = 0.31, if binding can be assumed to be sufficiently strong (K_d < 2 nM, data not shown). Under these conditions, binding in terms of CaM would be virtually complete, as indicated by the hatched bar and a per-site binding degree of = 1.0.

For the elevated Ca²⁺ case in the absence of additional ATP (Fig. 5 A, bottom), the binding distribution for the red-labeled CaM only is strikingly similar to the distribution in the MgATP case (Fig. 5 B, bottom). Strictly speaking, we cannot exclude from this data that the red distribution stems from a weaker binding affinity (K_d = 1.5–3 µM) than in Case B. However, this interpretation would automatically imply a high fraction of labeled CaM (r = 0.5 to 1), which is unlikely. It is therefore more reasonable that binding is again very strong (K_d < 20 nM) and that the fraction of labeled CaM is close to its limiting value for small K_d values, which is r = 0.35. The similar distributions of the red CaM obtained in the two cases thus probably both reflect full binding as indicated by the hatched bars.

In the EGTA case (Fig. 5 C, bottom), the distribution of red-labeled CaM on the CaMKII in the plot is strongly shifted to the left. At least more than half of the CaMKII molecules do not have any red-labeled CaM bound (C_i=0). The quantitative interpretation of this low binding case in terms of the distribution of total CaM, the dissociation constant K_d and the binding degree , only provide limits for these values and would again depend on the fractions of red-labeled and unlabeled CaM. The hatched bars in Fig. 5 C (bottom) provide an example of how the total CaM would distribute, if the red-labeled fraction constituted r = 0.3. Modeled CaM distributions for other values of r are also shown denoted by the shaded curves. If the red-labeled CaM is assumed to constitute, at a minimum, a fraction of r = 0.3 of the total CaM, the calculation yields a minimal K_d of 36 µM. Considering that the cross-correlation amplitude may even be closer to zero due to the small nonspecific cross-talk signal that may contribute, and bearing in mind that it is difficult to determine its value accurately with this low signal, binding can only be concluded to be negligible and the K_d to be probably far above 36 µM.

	DISCUSSION

TOP ABSTRACT INTRODUCTION THEORY MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
In this study we demonstrated the effectiveness of two-photon cross-correlation analysis applied to an intracellular assay for studying the complex binding dynamics of CaM and CaMKII. We first devised the appropriate procedures for setting up a TPCCS assay for intracellular use and thoroughly characterized the biochemical system. Next, an extended cross-correlation analysis was developed for cases beyond those with simple 1:1 stoichiometry and constant brightnesses upon binding. The application of this concept, which incorporated a description of the complex binding stoichiometry and a simple ansatz for fluorescence quenching, yielded values for the in vitro binding parameters. TPCCS was then applied directly to live cells, and the Ca²⁺-dependency of the binding reaction known from in vitro experiments was confirmed. The following detailed discussion of our results and conclusions focuses on the accompanyi