PMCCPMCCPMCC

Search tips
Search criteria 

Advanced


 
Formats
Article sections
Authors
Related links
Logo of bmcbiophysBioMed CentralBiomed Central Web Sitesearchsubmit a manuscriptregisterthis articleBMC BiophysicsJournal Front Page
 
BMC Biophys. 2011; 4: 10.
Published online 2011 May 9. doi:  10.1186/2046-1682-4-10
PMCID: PMC3126788
Copyright ©2011 Lichten and Swain; licensee BioMed Central Ltd.
A Bayesian method for inferring quantitative information from FRET data
Catherine A Lichten1,2 and Peter S Swaincorresponding author2
1Department of Physiology, McGill University, 3655 Promenade Sir William Osler, Montreal, Quebec H3G 1Y6, Canada
2Centre for Systems Biology at Edinburgh, University of Edinburgh, Edinburgh EH9 3JD, UK
corresponding authorCorresponding author.
Catherine A Lichten: catherine.lichten/at/mail.mcgill.ca; Peter S Swain: peter.swain/at/ed.ac.uk
Received February 14, 2011; Accepted May 9, 2011.
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.
Abstract
Background
Understanding biological networks requires identifying their elementary protein interactions and establishing the timing and strength of those interactions. Fluorescence microscopy and Förster resonance energy transfer (FRET) have the potential to reveal such information because they allow molecular interactions to be monitored in living cells, but it is unclear how best to analyze FRET data. Existing techniques differ in assumptions, manipulations of data and the quantities they derive. To address this variation, we have developed a versatile Bayesian analysis based on clear assumptions and systematic statistics.
Results
Our algorithm infers values of the FRET efficiency and dissociation constant, Kd, between a pair of fluorescently tagged proteins. It gives a posterior probability distribution for these parameters, conveying more extensive information than single-value estimates can. The width and shape of the distribution reflects the reliability of the estimate and we used simulated data to determine how measurement noise, data quantity and fluorophore concentrations affect the inference. We are able to show why varying concentrations of donors and acceptors is necessary for estimating Kd. We further demonstrate that the inference improves if additional knowledge is available, for example of the FRET efficiency, which could be obtained from separate fluorescence lifetime measurements.
Conclusions
We present a general, systematic approach for extracting quantitative information on molecular interactions from FRET data. Our method yields both an estimate of the dissociation constant and the uncertainty associated with that estimate. The information produced by our algorithm can help design optimal experiments and is fundamental for developing mathematical models of biochemical networks.
Background
Proteins work together continuously in the cells of all living things, generating cascades of reactions that are vital for life. To fully understand each individual protein's task requires discovering the timing, location, and strength of its interactions. To acquire this detailed information, fluorescence microscopy methods are ideal because they can provide dynamic, single-cell data at high spatial resolution [1,2]. One fluorescence tool that enables researchers to observe protein interactions in living cells is Förster resonance energy transfer (FRET). FRET data has the potential to yield biochemical constants, which are critical for modeling biological systems, but measuring protein interactions from FRET data requires careful quantitative analysis.
FRET is a physical process where a molecule in an excited energetic state (the donor) transfers energy to a nearby ground-state molecule (the acceptor). The chance that an excited donor will transfer its energy to an acceptor (known as the FRET efficiency, Efr) depends on the distance between the donor and acceptor (which must be between 1 and 10 nm for FRET to occur), their relative orientation, and the extent to which the donor's fluorescence emission peak overlaps the excitation peak of the acceptor [3]. Optimal conditions for energy transfer occur when the distance between donor and acceptor is minimal, the molecules' electric dipoles are aligned, and the spectral overlap is significant.
When FRET is used to study protein interactions in living cells, the proteins under investigation are fused to fluorescent tags (often variants of the green fluorescent protein) that act as the donors and acceptors. When the proteins interact, they bring the fluorescent tags together so that FRET may occur. FRET increases the number of photons emitted by acceptors and reduces both the number of photons emitted by donors and the donor's fluorescence lifetime. To observe these effects, the most common techniques for collecting FRET data include fluorescence lifetime imaging (FLIM) and using a fluorescence microscope or spectrofluorometer to record fluorescence intensity after exposing samples to light that mainly excites either donors or acceptors [4-6].
For FRET data to reveal information about the underlying protein interactions, complicating factors must be dealt with. One confounding issue is spectral contamination, which arises from the requirement that the donor and acceptor must have overlapping spectra for FRET to occur. Due to the overlap, excitation light intended to excite one fluorophore may also excite the other (spectral crosstalk), and, conversely, one fluorophore may emit photons in the emission range of the other (spectral bleed-through). Several FRET analysis methods address this issue by calculating a FRET index, which is the FRET signal corrected for contamination from spectral overlap and normalised by donor or acceptor concentrations [7-12]. Although they are straightforward to calculate, FRET indices often bear an indirect and nonlinear relationship to the underlying concentrations and strength of protein interactions [13,14]. Alternatively, the ratio of donors to acceptors and the apparent FRET efficiency, which is the product of the intrinsic FRET efficiency and the fraction of donors (or acceptors) in complex, can be calculated [15-19,14]. Such quantities can be measured, provided calibrations are carried out using constructs consisting of a donor linked to an acceptor [17,19], but again do not relate directly to the dissocation constant (Kd) of an interaction.
Another significant challenge is that both the Kd and the FRET efficiency, Efr, affect the observed signals, and so neither can be determined independently of the other based on data from a single sample. Efr has been estimated from separate acceptor photobleaching [14] and FLIM experiments [15], but these approaches have drawbacks. Acceptor photo-bleaching is slow, irreversible, often fails to bleach all acceptors, and yields the intrinsic FRET efficiency only when all donors are in complex with acceptors [4-6]. FLIM requires specialized equipment, is often slower than standard fluorescence imaging, and analysis of FLIM data is complicated by the multiexponential fluorescence decays of fluorescent proteins [20]. There is however an alternative to directly measuring Efr: Many values of Kd and Efr will be consistent with data from a single cell or sample, but, by taking a set of data from samples that contain varying concentrations of donors and acceptors and analyzing it altogether, it is possible to find the unique values of Kd and Efr that are consistent with the ensemble [16,21].
Given these different approaches, it is not always obvious which one should be applied in different situations and there is no consensus on the statistical analysis, with each method processing the data differently and most giving no procedure to test the reliability of any estimates. A general method for inferring the Kd along with the uncertainty of that inference is necessary for the standardization of quantitative FRET measurements, the design of informative experiments, and for providing in vivo parameters for developing models of protein networks [22-24].
Here we propose a systematic analysis method that explicitly includes models of the photophysics and underlying chemical interactions and of measurement noise. Building on a spectral model for FRET [25], we develop a Bayesian algorithm to infer both the interaction strength (in terms of the in vivo Kd) and the FRET efficiency. Applying our algorithm to simulated data, we gain insight into how both experimental design factors such as measurement noise, number of measurements, fluorophore concentrations and ratios, and prior information impact our estimate and its uncertainty.
Results and Discussion
Algorithm
Overview of the mathematical model
FRET enables the study of molecular interactions in diverse settings. To design a widely applicable analysis technique, we consider a general system containing proteins (or other molecules) that form complexes and are labelled with fluorophores that act as FRET donors and acceptors (Figure ​(Figure1A).1A). We assume that all the molecules of interest are fluorescently tagged, so instead of referring to 'donor-tagged proteins' or 'acceptor-tagged proteins', we simply refer to 'donors' and 'acceptors'.
Figure 1
Figure 1
Overview of interaction and photophysical model. (A) In the underlying interaction, a donor-tagged protein binds an acceptor-tagged protein with dissociation constant Kd. (B) During a FRET experiment, photophysical events occur such as FRET, which occurs (more ...)
Our model relies on a few other assumptions. First, we assume that donors and acceptors may be free or form bimolecular, donor-acceptor complexes (i.e. [D] + [A] [right harpoon over left harpoon] [DA]) with dissociation constant An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i2.gif. Whether free or in complex, the fluorophores undergo fluorescence excitation and emission (Figure ​(Figure1B).1B). In complexes, the donor and acceptor are together and FRET may also occur. We assume that FRET due to random collisions between donors and acceptors occurs seldom enough to be neglected. Finally, we assume that donor-acceptor complexes all have the same FRET efficiency, which is denoted Efr and refers to the fraction of instances where exciting a donor in complex leads to excitation of the acceptor.
For our general mathematical description of the FRET system, we use a previously described spectral mixing model [25] with slightly modified matrices. This model relates concentrations of fluorescent molecules to the fluorescence observed at various excitation wavelengths and ranges of emission wavelengths (spectral channels). It can be written as a matrix equation in the following form [6]:
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i3.gif
(1)
where C is a vector of the concentrations of fluorophores (C = ([D], [A], [DA])T ), Iobs is a vector of the fluorescence intensities observed for each spectral channel, and M is a matrix of the photophysical parameters that relate the concentrations of fluorophore to the observed intensities (Figure ​(Figure1B1B).
The observed fluorescence intensities could represent data obtained from fluorescence microscopy, spectrophotometry or a flow cytometer. Microscopy images would require standard image processing steps for quantification, such as subtracting background fluorescence and defining regions of interest such as areas within the cytosol or nucleus for localized fluorescence [26].
Photophysics in the mathematical model
The matrix M contains all the relevant information about the photophysical processes we model. In addition to the expected direct excitation and emission, we include resonance energy transfer, cross-talk (unintended excitation) and bleed-through (acceptors emitting in the donors' typical emission range).
The spectral mixing framework can represent any number of spectral channels by increasing the dimensions of the matrix, but we focus on the particular case of three, which describes the most common 'three-cube' FRET experiments [12,7-9,14,18,17,21,16]. The first spectral channel, known as the donor channel, consists of an excitation wavelength that primarily excites donors and an emission filter that collects photons from the donor's emission range. The second channel, the acceptor channel, consists of an excitation wavelength that primarily excites acceptors and an emission filter that collects photons from the acceptor's emission range. The third channel, the FRET channel, combines the donor channel's excitation wavelength with the acceptor channel's emission filter, so that donors are excited preferentially and the emissions are filtered to primarily collect photons from acceptors. For three-cube FRET, the detailed form of Eq (1) is
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i4.gif
(2)
where An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i5.gif. The parameters An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif each represent a combination of excitation and emission information [25]:
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i7.gif
(3)
Here the subscript i refers to the spectral channel (for the three-cube case, it is 1 for the donor channel, 2 for the acceptor channel, and 3 for the FRET channel) and the superscript (S) refers to the fluorescent species (D, the donor; A, the acceptor; and DA, the complex). The variable Ii is the illumination intensity for the excitation wavelength of channel i; An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i8.gif is the molar extinction coefficient of species (S) with excitation wavelength from channel i; Q(S) is the quantum yield of species (S); and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i9.gif is the product of the emission of species (S) at the wavelengths of channel i and the sensitivity of the detector to the emitted photons. We describe a calibration procedure for obtaining these values in the Methods.
To illustrate Eq (2), consider the expression it yields for An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i10.gif, the intensity in the donor channel. The intensity will be the sum of contributions from free donors, free acceptors, and donor- acceptor complexes. The contribution from free donors is An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i11.gif and the contribution from free acceptors is An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i12.gif, which results from cross-talk and bleed-through. The complex, DA, can potentially produce contributions from donors not undergoing FRET An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i13.gif, from acceptors An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i14.gif and from FRET An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i15.gif. Adding these contributions together, we obtain the following expression for the predicted intensity in the donor channel, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i10.gif,
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i16.gif
(4)
Continuing the matrix multiplication yields analogous expressions for An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i17.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i18.gif.
Molecular interactions in the mathematical model
To include Kd in our model, we model the chemical equilibrium, [D] + [A] [right harpoon over left harpoon] [DA]. Defining [D0 ] = [D] + [DA] and [A0] = [A] + [DA], and given that An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i19.gif, we can compute the values of [DA], [D], and [A] in terms of [D0], [A0], and Kd.
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i20.gif
(5)
with [D] = [D0] - [DA] and [A] = [A0] - [DA].
Bayesian inference
Our aim is to infer the values of Kd and Efr from the data. To this end, Bayesian inference is an ideal tool. Given a model, it quantifies the knowledge gained from a set of experimental observations about the parameters of interest, including the uncertainty of that knowledge [27]. That is, it allows us to update the prior probability distribution, P (Kd, Efr), which represents what is known about the parameters before data has been acquired, to the posterior probability distribution, P (Kd, Efr|data), which reflects both the prior information and what is learned from the data. According to Bayes' theorem,
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i21.gif
(6)
The likelihood function, P (data|Kd, Efr), can be explicitly expressed through our photophysical model and a model of measurement noise. It quantifies the likelihood of different values for Kd and Efr given the current data. We assume that the data follows a Gaussian distribution centred around the predicted intensity given by Eq (2). For three measurements, one from each of the three spectral channels, the likelihood is:
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i22.gif
(7)
The observed intensities for each channel are ID, IA, and IF. An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i10.gif, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i17.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i18.gif denote the predicted intensities and σD, σA and σF measure the magnitude of the measurement noise in each channel and are approximated by the variances of the data in that channel (Methods). For fitting, we use a Gaussian approximation of Eq (7) based around the values of D0 and A0 that maximize the likelihood (see section on marginalization of D0 and A0 in the Methods). Combining the likelihood and prior probability distribution allows us to calculate the posterior probability distribution and determine the most probable values for our parameters of interest.
Illustration of method
To illustrate our method, we simulated data from three cells (Figure ​(Figure2A,2A, left), although in practice, data could come instead from cellular compartments or other sub-cellular regions of interest. For simplicity, each cell contains equal concentrations of donor ([D0]) and acceptor ([A0]). In the first cell, [D0] = [A0] = 0.2 μM; in the second, [D0] = [A0] = 1 μM; and in the third, [D0] = [A0] = 5 μM. Other parameter values are given in the Methods. We then simulated three-cube FRET measurements, generating data for the donor, acceptor, and FRET channels (Figure ​(Figure2A,2A, right). We made a series of ten measurements for each channel and each cell, to which we added Gaussian noise so that the standard deviation of the measurements was around 5% of the mean signal strength. We define r as the ratio of the standard deviation of the measurements to the mean signal strength; in this case, r ≈ 0.05.
Figure 2
Figure 2
Analysis of typical FRET data recovers true values of Kd and Efr. A typical three-cube FRET experiment is simulated from three virtual cells, each containing the indicated concentrations of donor- and acceptor- tagged proteins (A, left). The data is summarized (more ...)
To obtain the posterior probability distribution that corresponds to this data for two parameters of interest, Kd and Efr, we used two methods, the results of which should agree. First, because we are restricting our analysis to two parameters, we can visualize the distribution by computing the energy, the negative logarithm of the posterior probability, over a grid of values of Kd and Efr. The resulting energy surface, shown in Figure ​Figure2B,2B, was minimal around the true values of the parameters. If we were to extend our analysis to more than two parameters, including uncertainty in the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif constants for instance, we would use a Markov chain Monte Carlo (MCMC) algorithm to sample from the parameter space. To show that the MCMC algorithm gives results consistent with the energy surface, we used the algorithm to sample the posterior probability distribution for Kd and Efr. We ran the algorithm three times, starting each time at different initial values. All three walks converged on the same small region of equally probable values: they are superimposed on the energy surface in Figure ​Figure2B.2B. Both methods therefore recover the true values used to generate the data as optimal parameter values. The approximate posterior probability distributions for Kd and Efr are shown in Figure ​Figure2C.2C. The distributions are estimated as the histograms of the three walks, excluding the burn-in period before convergence [28]. They peak near the parameter values used to simulate the data, indicating that the highly probable values for Kd and Efr predicted by the methods are accurate.
Inferring other quantities from FRET data
To further illustrate the versatility of our method, we show that we can also use our model and Bayesian framework to estimate other parameters. Two instrument-independent parameters that have been a focus of interest are the apparent FRET efficiency, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i23.gif, and the ratio, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i24.gif[15-19,14]. The analogous apparent FRET efficiency for the acceptor, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i25.gif, is the product of Ed and rda. We can estimate these quantities using our method, provided that a calibration has been carried out with cells expressing only donors, only acceptors, and only a donor-acceptor construct. The FRET efficiency of the construct need not be known. Supposing we have calibrated the parameters An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif of Eq 2 relative to An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i26.gif, which is taken to be unity, then we can express Eq 2 just in terms of [A0] with [D0] = rda[A0] and Efr[DA] = Ed rda[A0]. We can analytically maximize the posterior probability of rda and Ed as a function of [A0] by setting the derivative of the probability with respect to [A0] equal to zero and solving for this optimal [A0] in terms of rda and Ed. Our MCMC algorithm can be used to sample values of rda and Ed that are consistent with the data.
Testing
Algorithm reflects data quality
An algorithm for estimating parameter values should report not just an estimate of the most probable parameter values but also the reliability of that estimate. To evaluate this aspect of our algorithm, we measured both the error and the uncertainty of its output. We defined the error as the discrepancy between the mean of the posterior probability distribution and the true value. We calculated error as An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i27.gif, meaning that a perfect estimate would have zero error. Because under experimental conditions we would not know the true value, we also calculated the uncertainty of the estimate as an indicator of the reliability even when true values are unknown. The uncertainty corresponds to the width of the posterior probability distributions for Kd and Efr. To prevent the magnitude of the parameter of interest from skewing the uncertainty, we use the relative uncertainty, calculated as the coefficient of variation (the standard deviation divided by the mean).
We first tested our algorithm in the presence of varying levels of measurement noise. As the examples in Figure ​Figure3A3A suggest, when measurement noise increased, the optimal regions in Kd and Efr-space grew wider and longer and the peak of the posterior probability distribution for Kd often deviated from the true value (see Figure ​Figure3A3A inset). In Figures ​Figures3C3C and ​and3E,3E, the error (upper plot) and uncertainty (lower plot) for estimating Kd are shown and confirm that both the error and uncertainty of the estimate increased with increasing noise. Similar results were obtained for the estimate of Efr (data not shown).
Figure 3
Figure 3
Kd estimates reflect data quality and quantity. We simulated and analyzed data with varying levels of measurement noise and numbers of measurements/cell/channel. The locations visited by MCMC walks (dots) for the three noise levels (A) and the three numbers (more ...)
Next, we verified that gathering more data would improve the quality of the estimates of Kd and Efr. Such data could be obtained by, for example, collecting several images of the same sample or dividing a region of interest into subregions. The examples in Figure ​Figure3B3B show that with more measurements, the optimal region became smaller and the peak of the posterior probability distribution for Kd moved closer to the true value (see Figure ​Figure3B3B inset). Figures ​Figures3D3D and ​and3F3F show, as expected, that increasing the number of measurements decreases the error and uncertainty for the Kd estimate. The effect of data quantity on the error and uncertainty of the Efr estimate was similar (data not shown).
Overall, our algorithm reliably reflects the quality of the data analyzed in the uncertainty it gives of its estimate and accurately infers the values of Kd and Efr for 10 measurements per cell per channel with 10% measurement noise. We also note that in Figures ​Figures3C3C and ​and3D,3D, the mean error in Kd remained close to zero, indicating that neither measurement noise nor number of measurements systematically biased the location of the peak.
The total amount of donor, [D0], and the total amount of acceptor, [A0], affect inference of Kd and Efr
Inferring unique values of Kd and Efr requires variations in [D0] and [A0]
A challenge in analyzing FRET data is that both Efr and Kd affect the extent of FRET that may occur, making it difficult to tease out the contribution of each parameter. In general, for data from a single pair of donor and acceptor concentrations, a weak Efr and strong Kd can give rise to the same observations as a strong Efr and weak Kd. However, because our algorithm can analyze data from multiple cells at once, it can overcome this challenge.
To demonstrate the problem, we analyzed a set of simulated data from a single sample with [A0] = [D0] = 1 μM. As Figure ​Figure4A4A shows, the resulting region of optimal values had an elongated shape, indicating that many equally probable solutions exist for {Kd, Efr} for that dataset. To show that collecting more data from similar cells does not help, we repeated the experiment, simulating data from 3 cells with the same concentrations as the first. As Figure ​Figure4B4B shows, including this data made the region slightly narrower but did not change its elongated shape.
Figure 4
Figure 4
Concentration variations determine shape of region of high probability. Using the Monte Carlo algorithm to analyze data from a single cell containing [D0] = [A0] = 1 μM produced an elongated region of highly probable values (A), indicated by the (more ...)
The difficulty in determining Kd without knowing Efr or vice versa can be overcome by analyzing data from samples containing varying fluorophore concentrations [16,21]. To quantitatively justify how varying concentrations facilitates parameter fitting, we simulated data from three cells containing different concentrations of donors and acceptors. As Figure ​Figure4C4C shows, analyzing data from two different cells decreased the size of the region of equally probable values. Combining data from a cell containing [D0] = [A0] = 1 μM with data from a cell containing 0.5 μM eliminated high values of Kd (Figure ​(Figure4C);4C); combining data from a 1 μM cell and a 5 μM cell eliminated low values of Kd (Figure ​(Figure4C,4C, inset). As Figure ​Figure4D4D indicates, analyzing data from all three cells (0.5 μM, 1 μM, and 5 μM) together eliminated higher and lower Kd values, leaving a small region of equally probable values that included the true values. To show why varying fluorophore concentrations is effective, we analyzed data from each of the three cells in Figure ​Figure4D4D individually. The optimal regions for each cell are plotted in Figure ​Figure4E4E and clearly intersect at the true value. Varying the concentrations more broadly can further shrink the solution region (Figure ​(Figure4F4F).
Inference is limited if [D0] = [A0] [dbl greater-than sign] Kd or [D0] = [A0] [double less-than sign] Kd
Plots of complex formation vs. Kd (Figure ​(Figure5A,5A, insets) show that only Kd values within a limited range significantly affect complex formation: outside this range, further changes in Kd have little effect. For instance, when [D0] = [A0] ≈ Kd, small changes in Kd in either direction affect [DA] (Figure ​(Figure5A5A centre, inset). However, when [D0] = [A0] [double less-than sign] Kd (Figure ​(Figure5A,5A, left inset) or [D0] = [A0] [dbl greater-than sign] Kd (Figure ​(Figure5A,5A, right inset), only changes to Kd in one direction affect the amount of [DA]. When a range of Kd values are consistent with the data, Kd cannot be fit to a single value.
Figure 5
Figure 5
Gaining insight into optimal experimental design. The approximate posterior probability distributions for Kd (A) have different shapes if the data analyzed was simulated from three cells containing equal concentrations of donors and acceptors which are (more ...)
By finding the posterior probability distribution of Kd, our algorithm predicts when a range of Kd values will be consistent with the data. First, we simulated data with equal concentrations of acceptors and donors that were very low or very high relative to Kd. As shown in Figure ​Figure5A,5A, the approximate posterior probability distributions of Kd for data simulated with relatively low or high concentrations form plateaus instead of single peaks. The plateaus indicate that a range of Kd values is equally probable and consistent with the data. The posterior distribution obtained for [D0] = [A0] [double less-than sign] Kd reveals that Kd could be any number greater than approximately 10-6M; the distribution for [D0] = [A0] [dbl greater-than sign] Kd reveals that Kd < 10-6M.
This analysis illustrates how posterior probability distributions for the parameters of interest can be more informative than single values. For example, information from posterior distributions can be useful for improving experimental design. If one were to obtain plateaued distributions like the two shown in Figure ​Figure5A,5A, one could recognize that the fluorophore concentrations were too low or high relative to Kd and devise a more informative experiment.
The ratio [D0]:[A0] affects inference
We also used our algorithm to explore the effects of varying the ratio, [D0] : [A0], on the algorithm's ability to infer Kd (Figure ​(Figure5B)5B) and how these effects arise. While other authors have reported that FRET data becomes increasingly unreliable as the ratio [D0] : [A0] deviates from unity [13], they measured how noise propagates in FRET formulas and did not directly address how the ratio impacts the measurement of interaction strength. As Figure ​Figure5B5B shows, the width of the approximate posterior distributions for Kd broaden as the ratio increases, meaning that the uncertainty in the estimate increases. The peaks of the distributions however remain close to the true values. The same effect occurs when the ratio decreases (data not shown).
The insets in Figure ​Figure55 illustrate why the ratio affects the uncertainty, showing data from the donor and FRET channels in the presence and absence of FRET. Here, FRET contributes only a fraction of the observed signals. This small contribution is particularly evident in the FRET channel, where the change in the signal due to FRET is small compared to the measurement error: measurement error is 3% and the increase in the FRET channel from FRET is only 2.9% for [D0] : [A0] = 1 and drops to 0.7% for [D0] : [A0] = 10 (FRET is included by changing the FRET efficiency from 0 to 0.4). As the ratio increased, the relative shortage of acceptors meant that fewer complexes could form, making the contributions from FRET smaller and causing the uncertainty of the estimate to grow. In the absence of measurement noise, however, our algorithm can estimate Kd even from very extreme ratios. These data also demonstrate that with our method, we can make inferences about Kd even when FRET causes perceptible changes in only one channel (e.g. the donor channel here).
Knowing how spectral overlap and the ratio of donors to acceptors affect the inference of protein interaction strength helps make it possible to design informative experiments. For example, we find that in each channel, the difference generated by FRET depends on spectral overlap. In the donor channel, where FRET decreases the signal, the relative change is An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i29.gif, meaning that the observed change depends on the difference between acceptor bleed-through, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i30.gif, and donor fluorescence, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i31.gif. Analogous changes occur in the other channels. When [D0] = [A0], the change in the donor channel is proportional to An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i32.gif, the change in the acceptor channel to An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i33.gif and the change in the FRET channel to An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i34.gif. In each channel, increasing the difference between An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i35.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i36.gif would make experiments more informative.
Prior information improves estimate
A further advantage of our method is its flexibility. The Bayesian framework makes it possible to incorporate additional details we know about the system as prior information, allowing us to more accurately represent the system being analyzed and potentially improve our estimate of the parameters of interest. To exploit this feature, we tested whether including additional information about the FRET efficiency would improve our inference of Efr and Kd. Information about Efr could be obtained, for example, from lifetime imaging experiments [29].
Figure ​Figure66 shows the approximate posterior probability distributions for Efr and Kd that result from analyzing the same data in the presence and absence of prior information for Efr. We included slightly inaccurate prior information with large uncertainty, supposing that Efr was measured to be 0.45 (instead of the true value, Efr = 0.4) with an uncertainty of 0.15. Even with this limited information, the resulting histograms for both Kd and Efr have taller, narrower peaks which are closer to the true values compared to the histograms produced by analyzing the same data in the absence of the prior information. As expected, this improvement was more substantial when we analyzed a set of data consisting of just 3 measurements per cell per channel (Figure ​(Figure6,6, right column): prior information becomes more important the less relevant information there is in the data. As another form of prior information, we could also include error in the calibration measurements for the constants that relate molecular concentration to fluorescence in Eq (2). Incorporating this information would make it possible to monitor the impact of that error on the uncertainty of the final estimate of Kd.
Figure 6
Figure 6
Prior information on Efr improves estimate. Approximate posterior distributions for Efr (upper panels) and Kd (lower panels) obtained from using the MCMC algorithm to analyze the same datasets in the absence (red) or presence (blue) of prior knowledge (more ...)
Summary of testing
We have shown, using typical, simulated FRET data, that our algorithm accurately recovers the parameters of interest. It responded consistently and intuitively to changes in the amount of measurement noise present in the data and the quantity of data. We also used the posterior probability distributions we obtain for the parameters to gain insight into how the magnitude and variation of the donor and acceptor concentrations affect the ability of the algorithm to infer Kd. The algorithm gave informative results even when given data simulated with very high or low concentrations of donors and acceptors relative to the Kd although, as expected, with this data it could not converge on a unique pair of values for Kd and Efr but gave upper or lower bounds. Gaining an understanding of how experimental parameters such as measurement noise, number of measurements, spectral overlap and fluorophore concentrations affect inference allows experimenters to design optimal experiments that yield information with less uncertainty [24].
Conclusions
Fluorescence microscopy and FRET open a window onto the cell, allowing us to observe protein interactions as the cell functions as a complete system. However, for protein interaction information from FRET data to be integrated into models and improve our understanding of biological systems, it must be reliably quantified, including the uncertainty in the estimates produced [27].
For this purpose, we have presented an algorithm for inferring the most probable values of the absolute or relative in vivo dissociation constant and the FRET efficiency from three-cube FRET data. Our algorithm produces estimates in terms of the posterior probability distribution for the parameters of interest. Posterior probability distributions yield more extensive information than a single value can provide, conveying the reliability of the estimates through the shape and width of the distribution and giving upper or lower bounds on parameters if the data is not more informative. Our method requires only basic three-cube FRET data but has the flexibility to use data from more spectral channels and incorporate other types of data, such as a separate measure of the FRET efficiency. We have focused on using our algorithm to fit Efr and Kd, but it can also infer the apparent FRET efficiency and ratio of donors to acceptors, other quantities sometimes used as measures of FRET [15-19,14].
In the examples described here, we make a few assumptions, but these are not necessarily part of our methodology. First, we use molar extinction coefficients that must be measured separately or taken from the literature and assume that the literature measurements are valid in the cellular environment. Second, we assume Gaussian measurement noise, but our model could be straightforwardly adapted to include log-normal or other types of noise. Finally, we have not taken into consideration photo-bleaching, incomplete labelling or dark states, but our model could readily be extended to include these factors.
We have used simulated data to illustrate principles that apply to real data, showing how, in practice, it is best to infer the dissociation constant from FRET data. Three-cube FRET experiments are a well-established experimental technique [3-6]. For our method, the calibration procedure we describe in the Methods for quantifying the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif constants of Eq (2) is also necessary. The procedure is similar to existing calibration methods. Measuring relative values of An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i36.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i37.gif is a standard part of FRET data collection, necessary for calibrating the extent of spectral bleed-through and cross-talk. We require calibration with a donor-acceptor construct to relate the brightness of the donor to that of the acceptor, which is necessary for obtaining the relative Kd. This calibration is not always a part of standard FRET data collection, but it has been used in several studies [8,18,17,4,14,19,21]. Alternatively, measuring the absolute brightness per molecule, which is necessary for measuring the absolute Kd, is more technically challenging but has been achieved [30,31].
We have demonstrated that to infer the values of Kd and Efr, it is important that the FRET data analyzed come from samples containing varying concentrations of donor- and acceptor-tagged proteins. The precise amount of variation required is difficult to predict because it depends on the measurement noise present in the data, the number of measurements made, and the extent of spectral bleed-through and cross-talk. Naturally occurring variation in protein expression from cell to cell generated sufficient variability for estimating the FRET efficiency and relative dissociation constant [21], at least in transiently transfected mammalian cells.
While a number of ad hoc methods exist for quantifying protein interactions using FRET, our method contributes something new in that it makes plain the spectral information being used and the bio-chemical assumptions made about the system. Our model is general: it is not specific to three-cube data and could be used to analyze data with any number of spectral channels. It can be straight-forwardly adapted to other experimental situations, for instance to measure the dimerization affinity of a homodimeric protein. Our approach is a Bayesian analysis that reveals the uncertainty in estimated parameters and produces informative results for data from a wide-range of experimental set-ups. It focuses on determining what we can learn from experimental observations about the parameters of interest given clear assumptions and a systematic, statistical analysis based only on those assumptions.
Methods
Calibration of fluorescence constants, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif
The An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif constants, which relate concentrations of single fluorophores to observed fluorescence intensity, can be obtained by calibration. Calibration of spectral cross-talk and bleed-through is a standard part of FRET quantification [4-6]. If the fluorophore concentrations can be determined [30,31], then the values for the three An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i36.gif and three An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i37.gif constants can be obtained from samples containing only donors or only acceptors. This calibration makes it possible to measure the absolute Kd. Three-cube data from donor-only and acceptor-only samples correspond to the following equations (from Eq (2)):
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i38.gif
(8)
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i39.gif
(9)
To infer the values of the constants An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i36.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i37.gif and their uncertainty (due to measurement error), we can sample from posterior distributions for An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i36.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i37.gif given the calibration data. We will illustrate this procedure for a more complicated example below.
Determining the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i35.gif, which relate a complex undergoing FRET to its fluorescence in channel i, is more complicated because the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i35.gif depend on properties of both the donor and the acceptor. However, from Eq (3), An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i40.gif, but An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i41.gif. Therefore, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i42.gif, and so An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i35.gif can be obtained through knowledge of the ratio of molar extinction coefficients at the excitation wavelength for channel i. The values of the molar extinction coefficients, ε (D) and ε(A), may not always be available at these wavelengths, but they can be estimated from literature values of the molar extinction coefficients (usually measured at the fluorophore's excitation peak) and the excitation spectra of the donor and acceptor [25]. Interpreting the excitation spectra as the probability of the fluorophore becoming excited and assuming that the extinction coefficient is proportional to the probability of excitation, we can use the excitation spectra to rescale the literature value of the extinction coefficient to the excitation wavelength of channel i. This estimate is valid provided that the molar extinction coefficients obtained from the literature are not significantly different from the molar extinction coefficients in the cellular environment of the experiment (or provided that the change in environment affects both donor and acceptor similarly).
When the relationship between brightness and concentration cannot be determined absolutely, relative values for An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i36.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i37.gif can be obtained from samples containing only donors and only acceptors, along with samples containing linked donor-acceptor constructs. The constructs consist of a donor and an acceptor fluorophore separated by a short linker of 5-10 amino acids and have been used for FRET calibration in several studies [8,18,17,4,14,19,21]. The construct's FRET efficiency, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i43.gif, need not be known and can be determined using the procedure we describe below.
The three-cube measurements on samples containing only donors or only acceptors would correspond to Eqs (8) and (9). Three-cube data obtained from samples containing donor-acceptor constructs would correspond to the following equations:
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i44.gif
(10)
where we have replaced An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i35.gif by An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i45.gif. To infer the values of the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif constants given calibration data consisting of three-cube measurements made on samples with only donors, only acceptors, and donor-acceptor complexes, we first define the general likelihood function of these experiments for a single sample:
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i46.gif
(11)
where An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i47.gif, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i48.gif, and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i49.gif are the ith measurements in the donor, acceptor, and FRET channels, respectively. For a given sample, we assume that an equal number of measurements, n, will be made in each channel and σ quantifies the error of the measurements. The predicted intensities in channel k, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i50.gif, are given by Eqs (8), (9) and (10) and are a function of the constants An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif, the concentrations of the species (donors, acceptors, or donor-acceptor complexes), and also An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i43.gif in the case of the donor-acceptor construct. The full likelihood for all the calibration data together would be the product of three instances of Eq (11), one specified for each three-cube experiment (donor-only, acceptor-only, or donor-acceptor complex).
We are interested in the values of the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif but indifferent to the concentrations of fluorophores and complexes in the samples. For this reason and because the concentrations of fluorophores and complexes will vary for different samples, it is useful to marginalize the calibration likelihood over [A0], [D0], and [DA]. We can also eliminate σ by marginalization [28], assuming the measurement error is the same for all three spectral channels. If this assumption does not hold, one can define σk for the measurement error in each channel and either approximate each σk as equal to the standard deviation of the measurements in the kth channel or include the σk as parameters to be inferred in the procedure we describe below. Assuming a prior probability that only specifies positive values for the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i43.gif, marginalizing Eq (11) over the concentrations and σ yields the posterior probability of An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i43.gif and relative values of An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif:
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i51.gif
(12)
which is valid for the three cube experiments for each sample (donor, acceptor, or construct). For the donor-only sample, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i52.gif; for the acceptor-only sample, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i53.gif; and for the sample with the donor-acceptor construct, the ai are replaced with An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i54.gif. The posterior probability, including all the calibration data together, is obtained by multiplying the forms of Eq (12) for each of the three 3-cube experiments.
For a given data set, we can infer An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i43.gif and five of the six An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif constants relative to the remaining one, for instance An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i55.gif, which is set to unity. The numerator of Eq (12) remains unaltered by such a rescaling of the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif, and the alteration in the denominator is cancelled by the Jacobian required for the change in variables. One can use a Markov chain Monte Carlo method to sample the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i43.gif from Pcalib. Alternatively, one can use a numerical solver to find the most probable values by solving for the roots of the system of equations consisting of the derivatives of Eq (12) with respect to each of the six variables. Using this relative calibration procedure, the final Kd values obtained would relate to the true Kd by the scaling factor, for instance An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i55.gif.
Data simulation
We designed our simulated data to mimic the key features of experimental data, which could come from various sources, such as fluorescence reader measurements of solutions of purified proteins or images of cells from fluorescence microscopy that have been processed and quantified. To simulate data, we wrote a function in Matlab (The Mathworks, Natick, MA) that takes as input Efr, Kd, A0, D0, n, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif, and r where A0, D0 are vectors of length m (representing m cells or regions within cells), n is the number of measurements made on each cell or vesicle, and r is the strength of measurement noise. The function outputs a set of simulated three-cube FRET data from m samples with n measurements per sample.
For each pair of concentrations, [A0] and [D0], we calculate [D], [A], and [DA] using Eq (5). Next, we calculate the simulated experimental intensities ID, IA, and IF using Eq (1). Finally, we simulate measurement noise by adding Gaussian random numbers to the data. The variable r is a scalar between 0 and 1 that refers to the strength of the measurement noise relative to the mean of the observed signal. Noise with r = 10% would have Gaussian measurement noise with standard deviation that is 10% of the mean signal observed in each channel.
In our examples, we use Efr = 0.4 and Kd = 1 μM. Unless otherwise indicated, we simulate ten measurements per cell per channel and use the following set of An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif values, which were chosen to represent a donor-acceptor pair with considerable spectral cross-talk and bleed-through:
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i56.gif
The constants corresponding to the donor, acceptor, and complex fluorescing in their respective channels An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i57.gif are assigned the largest values. As donors and acceptors are often not equally bright, we set An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i58.gif. In the FRET channel, the constants corresponding to bleed-through from donors An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i59.gif and cross-talk from acceptors An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i60.gif have been set to approximately 20% of the constants for donor and acceptor in their respective channels An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i61.gif because spectral overlap contributes significantly in the FRET channel. In the acceptor channel, smaller values are assigned to An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i62.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i63.gif to describe the donor undergoing cross-talk and, subsequently, either bleed-through (An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i62.gif) or FRET (An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i63.gif) in the acceptor channel. It is unlikely that acceptors would emit photons detectable in the donor channel, so the constants corresponding to that process, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i64.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i65.gif, are very small, but they are non-zero to show that any measurable spectral contamination can be included.
Incorporating prior information
The Bayesian framework makes it possible to incorporate prior information about any of the parameters, including uncertainty in An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif. As an example, we include prior knowledge about Efr (see the Results). Such data could be collected in separate experiments under slightly different experimental conditions and in the presence of measurement noise. We define the prior probability of Efr by a Gaussian distribution with mean at the measured value An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i66.gif and variance σE, which reflects the confidence in the measurement:
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i67.gif
(13)
Markov chain Monte Carlo (MCMC) estimation
To sample Kd and Efr from the posterior probability distribution, we use a biased random walk. Although a simpler approach would be sufficient to explore two dimensions, we use this method because it will allow us to efficiently extend our search to sample other, additional parameters we may wish to infer, such as the values of An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif.
We use the Metropolis-Hastings algorithm [32,33]. It begins at a random location in parameter space and takes random steps, moving in up to 3 dimensions at a time with proposed steps drawn from a distribution that is symmetric about the current location. Proposals that increase the posterior probability are always accepted; those that decrease it are accepted with probability An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i68.gif, where P (xj) is the posterior probability of the proposed step and P (xj-1) is the posterior probability of the current location. The walk is run for a sufficiently long time (about 10,000 steps) to generate independent samples from the posterior distribution and the step size chosen to maintain an acceptance rate of 40 - 60% [28].
Once the walk has converged with the energy fluctuating around a minimum value, we record the steps taken and use the histogram of these sampled values as our estimate of the posterior probability distribution. From this estimated distribution, we obtain the mean and standard deviation of each parameter being inferred. The algorithm is summarized below.
Marginalization of D0 and A0
We have defined the likelihood in Eq 7 and assume that the measurement noise in each channel is independent. Because we have a Gaussian model, the values of the measurement noise parameters, σD, σA and σF , that maximise the posterior probability are approximately equivalent to the observed variances of the data from the donor, acceptor, and FRET channels, respectively [28]. However, it would also be possible to fit An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i69.gif, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i70.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i71.gif directly, along with Kd and Efr.
Although Kd, Efr, D0 and A0 are important parameters, we cannot measure them directly. For our purposes, we are interested in the values of Kd and perhaps Efr, but not the values of D0 and A0. Rather than fit D0 and A0, we marginalize or integrate them out:
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i72.gif
(14)
To indicate that we have no knowledge about the values of D0 and A0, we set the prior, P (D0, A0), to a constant for positive D0 and A0 (and 0 otherwise).
As a result,
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i73.gif
(15)
Because this expression is difficult to integrate analytically, we consider the 'energy', E = - log(P(data|Kd, Efr, D0, A0)), and approximate E with a two-dimensional Taylor expansion about An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i74.gif, where An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i74.gif minimizes E (and thus maximizes the likelihood). Ignoring terms above second order,
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i75.gif
(16)
where [nabla][nabla]E is the Hessian, or matrix of second-derivatives of E with respect to D0 and A0.
Therefore,
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i76.gif
(17)
This approximation of the likelihood results in a Gaussian integrand, which we can then integrate analytically [34]. Note that we do not use the limits of integration (0, ∞) which would be appropriate for non-negative concentration values. We instead use the limits (- ∞, ∞) to make the integral simple. It is a valid approximation provided that An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i77.gif and An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i78.gif are sufficiently large and it consistently yields appropriate results in practice.
In summary,
A mathematical equation, expression, or formula. Object name is 2046-1682-4-10-i79.gif
(18)
where |H| is the determinant of the Hessian, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i80.gif. Although Kd and Efr do not appear explicitly in the final form of Eq (18), both An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i81.gif and H depend on Kd and Efr.
Algorithm summary
Our algorithm for sampling from the posterior probability distribution of (Kd, Efr) with a prior probability distribution for Efr consists of the following steps:
1. Perform calibration to obtain values (absolute or relative) for the An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i6.gif constants.
2. Define prior probability distributions for parameters to be inferred based on initial information that is known, if any.
3. Run Markov chain Monte Carlo algorithm:
Choose initial step An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i82.gif and compute the posterior probability, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i83.gif, as outlined in (b) below.
For j = 2 to n,
(a) Choose proposal step An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i84.gif.
(b) Compute the posterior probability, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i85.gif, for that step:
i. Find An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i74.gif, the values of {D0, A0} that minimise the energy of the posterior probability using nonlinear optimisation (the Nelder-Mead simplex algorithm implemented in Matlab's fminsearch function (The Mathworks, Natick, MA)).
ii. Compute An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i80.gif using analytical expressions for the second derivatives of E computed in Mathematica (Wolfram Research, Champaign, IL) and exported to Matlab.
iii. Using Eq (18), compute the likelihood, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i86.gif.
iv. Using Eq (13), compute the prior, An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i87.gif.
v. Compute the posterior probability distribution, Pj, which is the likelihood times the prior probability.
(c) Check whether to accept the move to An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i88.gif
• If Pj >Pj-1, accept.
• Otherwise, accept with probability An external file that holds a picture, illustration, etc. Object name is 2046-1682-4-10-i89.gif.
Steps (a)-(c) are repeated until j = n. After an initial burn-in period where the energy reaches a minimum and the walk achieves a stationary distribution, record samples of Kd and Efr.
4. Repeat step 3, varying initial (Kd, Efr). To ensure convergence, the stationary distributions of all walks should overlap [28].
Availability
We have made our data simulation and analysis software available at http://swainlab.bio.ed.ac.uk/software/FRET.
Authors' contributions
PSS conceived of the study. CL and PSS developed the methodology. CL implemented and tested the algorithm. Both authors wrote the paper and read and approved the final manuscript.
Acknowledgements
We thank Marko Laine (Finnish Meteorological Institute) for providing his MCMC code and helpful feedback. We also wish to acknowledge Stephane Laporte, May Simaan and Jay Nadeau (McGill University) for useful discussions and Andrea Weisse, Bruno Martins, Christos Josephides, Clive Bowsher (University of Edinburgh), Nathan Scales (McGill University), and Vahid Shahrezaei (Imperial College London) for helpful comments on the manuscript. CL and PSS are supported by the Scottish Universities Life Sciences Alliance.
References
  • Verveer P, Bastiaens P. Quantitative microscopy and systems biology: seeing the whole picture. Histochemistry and cell biology. 2008;130(5):833–843. doi: 10.1007/s00418-008-0517-5. [PubMed] [Cross Ref]
  • Megason S, Fraser S. Imaging in systems biology. Cell. 2007;130(5):784–795. doi: 10.1016/j.cell.2007.08.031. [PubMed] [Cross Ref]
  • Lakowicz J. Principles of Fluorescence Spectroscopy. New York, NY: Springer; 2006.
  • Vogel SS, Thaler C, Koushik SV. Fanciful FRET. Sci STKE. 2006;2006:re2. doi: 10.1126/stke.3312006re2. [PubMed] [Cross Ref]
  • Piston DW, Kremers GJ. Fluorescent protein FRET: the good, the bad and the ugly. Trends Biochem Sci. 2007;32:407–14. doi: 10.1016/j.tibs.2007.08.003. [PubMed] [Cross Ref]
  • Hoppe A. In: Imaging Cellular and Molecular Biological Functions. Shorte SL, Frischknecht F, editor. Heidelberg: Springer; 2007. Quantitative FRET microscopy of live cells; pp. 157–180.
  • Gordon G, Berry G, Liang X, Levine B. Quantitative fluorescence resonance energy transfer measurements using fluorescence microscopy. Biophysical Journal. 1998;74:2702–2713. doi: 10.1016/S0006-3495(98)77976-7. [PubMed] [Cross Ref]
  • Xia Z, Zhou Q, Lin J, Liu Y. Stable SNARE complex prior to evoked synaptic vesicle fusion revealed by fluorescence resonance energy transfer. J Biol Chem. 2001;276:1766–71. [PubMed]
  • Nagy P, Vámosi G, Bodnár A, Lockett SJ, Szöllösi J. Intensity-based energy transfer measurements in digital imaging microscopy. Eur Biophys J. 1998;27:377–89. doi: 10.1007/s002490050145. [PubMed] [Cross Ref]
  • Elangovan M, Wallrabe H, Chen Y, Day RN, Barroso M, Periasamy A. Characterization of one- and two-photon excitation fluorescence resonance energy transfer microscopy. Methods. 2003;29:58–73. doi: 10.1016/S1046-2023(02)00283-9. [PubMed] [Cross Ref]
  • Trón L, Szöllósi J, Damjanovich S, Helliwell SH, Arndt-Jovin DJ, Jovin TM. Flow cytometric measurement of fluorescence resonance energy transfer on cell surfaces. Quantitative evaluation of the transfer efficiency on a cell-by-cell basis. Biophysical Journal. 1984;45:939–46. doi: 10.1016/S0006-3495(84)84240-X. [PubMed] [Cross Ref]
  • Youvan D, Silva C, Bylina E, Coleman W, Dilworth M, Yang M. Calibration of fluorescence resonance energy transfer in microscopy using genetically engineered GFP derivatives on nickel chelating beads. Biotechnology et alia. 1997;3:1–18.
  • Berney C, Danuser G. FRET or no FRET: a quantitative comparison. Biophysical Journal. 2003;84:3992–4010. doi: 10.1016/S0006-3495(03)75126-1. [PubMed] [Cross Ref]
  • Zal T, Gascoigne NRJ. Photobleaching-corrected FRET efficiency imaging of live cells. Biophysical Journal. 2004;86:3923–39. doi: 10.1529/biophysj.103.022087. [PubMed] [Cross Ref]
  • Hoppe A, Christensen K, Swanson JA. Fluorescence resonance energy transfer-based stoichiometry in living cells. Biophysical Journal. 2002;83(6):3652–64. doi: 10.1016/S0006-3495(02)75365-4. [PubMed] [Cross Ref]
  • Erickson MG, Alseikhan BA, Peterson BZ, Yue DT. Preassociation of calmodulin with voltage-gated Ca(2+) channels revealed by FRET in single living cells. Neuron. 2001;31(6):973–85. doi: 10.1016/S0896-6273(01)00438-X. [PubMed] [Cross Ref]
  • Chen H, Puhl H, Koushik S, Vogel S, Ikeda S. Measurement of FRET efficiency and ratio of donor to acceptor concentration in living cells. Biophysical Journal. 2006;91:L39–41. doi: 10.1529/biophysj.106.088773. [PubMed] [Cross Ref]
  • Elder A, Domin A, Schierle GK, Lindon C, Pines J, Esposito A, Kaminski C. A quantitative protocol for dynamic measurements of protein interactions by FRET-sensitized fluorescence emission. Royal Society Interface. 2009;6:1–24. doi: 10.1098/rsif.2008.0327. [Cross Ref]
  • Wlodarczyk J, Woehler A, Kobe F, Ponimaskin E, Zeug A, Neher E. Analysis of FRET signals in the presence of free donors and acceptors. Biophysical Journal. 2008;94:986–1000. doi: 10.1529/biophysj.107.111773. [PubMed] [Cross Ref]
  • Gadella T. FRET and FLIM Techniques. Elsevier Science; 2009.
  • Chen H, III HP, Ikeda S. Estimating protein-protein interaction affinity in living cells using quantitative Förster resonance energy transfer measurements. Journal of Biomedical Optics. 2007;12:054011. doi: 10.1117/1.2799171. [PubMed] [Cross Ref]
  • Bruggeman F, Westerhoff H. The nature of systems biology. TRENDS in Microbiology. 2007;15:45–50. doi: 10.1016/j.tim.2006.11.003. [PubMed] [Cross Ref]
  • Piehler J. New methodologies for measuring protein interactions in vivo and in vitro. Curr Opin Struct Biol. 2005;15:4–14. doi: 10.1016/j.sbi.2005.01.008. [PubMed] [Cross Ref]
  • Kreutz C, Timmer J. Systems biology: experimental design. FEBS Journal. 2009;276(4):923–942. doi: 10.1111/j.1742-4658.2008.06843.x. [PubMed] [Cross Ref]
  • Neher RA, Neher E. Applying spectral fingerprinting to the analysis of FRET images. Microsc Res Tech. 2004;64:185–95. doi: 10.1002/jemt.20078. [PubMed] [Cross Ref]
  • Ljosa V, Carpenter AE. Introduction to the Quantitative Analysis of Two-Dimensional Fluorescence Microscopy Images for Cell-Based Screening. PLoS Comput Biol. 2009;5(12):e1000603. doi: 10.1371/journal.pcbi.1000603. [PMC free article] [PubMed] [Cross Ref]
  • Wilkinson D. Stochastic modelling for quantitative description of heterogeneous biological systems. Nat Rev Genet. 2009;10(2):122–133. doi: 10.1038/nrg2509. [PubMed] [Cross Ref]
  • MacKay D. Information Theory, Inference & Learning Algorithms. Cambridge: Cambridge University Press; 2002. [PubMed]
  • Wallrabe H, Periasamy A. Imaging protein molecules using FRET and FLIM microscopy. Current Opinion in Biotechnology. 2005;16:19–27. doi: 10.1016/j.copbio.2004.12.002. [PubMed] [Cross Ref]
  • Rosenfeld N, Perkins TJ, Alon U, Elowitz MB, Swain PS. A fluctuation method to quantify in vivo fluorescence data. Biophysical Journal. 2006;91(2):759–66. doi: 10.1529/biophysj.105.073098. [PubMed] [Cross Ref]
  • Digman MA, Dalal RB, Horwitz AF, Gratton E. Mapping the number of molecules and brightness in the laser scanning microscope. Biophysical Journal. 2008;94(6):2320–2332. doi: 10.1529/biophysj.107.114645. [PMC free article] [PubMed] [Cross Ref]
  • Metropolis N, Rosenbluth A, Rosenbluth M, Teller A, Teller E. Equation of state calculations by fast computing machines. J Chem Phys. 1953;21:1087–1092. doi: 10.1063/1.1699114. [Cross Ref]
  • Hastings W. Monte Carlo sampling methods using Markov chains and their applications. Biometrika. 1970;57:97–109. doi: 10.1093/biomet/57.1.97. [Cross Ref]
  • Sivia D, Skilling J. Data Analysis: a Bayesian Tutorial. Oxford: Oxford University Press; 2006.
Articles from BMC Biophysics are provided here courtesy of
BioMed Central