Photon counting histograms for diffusing fluorophores

Analysis of photon trajectories from diffusing single molecules

In single-molecule free diffusion experiments, molecules spend most of the time outside a laser spot and generate bursts of photons when they diffuse through the focal spot. Only these bursts contain meaningful information and, therefore, are selected using physically reasonable criteria. The analysis of the bursts must take into account the precise way they were chosen. We present new methods that allow one to accurately determine the brightness and diffusivity of individual molecule species from the photon arrival times of selected bursts. We derive analytical expressions for the distribution of inter-photon times (with and without burst selection), the distribution of the number of photons in a burst, and the distribution of photons in a burst with recorded arrival times. The theory accurately treats the bias introduced due to the burst selection criteria. We use a Maximum Likelihood (ML) method to find the molecule's photon count rate and diffusion coefficient from three kinds of data, i.e., the bursts of photons with recorded arrival times (burstML), inter-photon times in bursts (iptML), and the numbers of photon counts in a burst (pcML). The performance of these new methods is tested on simulated photon trajectories and on an experimental system, the fluorophore Atto 488.

Theory and Analysis of Single-Molecule FRET Experiments

Inter-dye distances and conformational dynamics can be studied using single-molecule FRET measurements. We consider two approaches to analyze sequences of photons with recorded photon colors and arrival times. The first approach is based on FRET efficiency histograms obtained from binned photon sequences. The experimental histograms are compared with the theoretical histograms obtained using the joint distribution of acceptor and donor photons or the Gaussian approximation. In the second approach, a photon sequence is analyzed without binning. The parameters of a model describing conformational dynamics are found by maximizing the appropriate likelihood function. The first approach is simpler, while the second one is more accurate, especially when the population of species is small and transition rates are fast. The likelihood-based analysis as well as the recoloring method has the advantage that diffusion of molecules through the laser focus can be rigorously handled.

FRET-based dynamic structural biology: Challenges, perspectives and an appeal for open-science practices

Single-molecule FRET (smFRET) has become a mainstream technique for studying biomolecular structural dynamics. The rapid and wide adoption of smFRET experiments by an ever-increasing number of groups has generated significant progress in sample preparation, measurement procedures, data analysis, algorithms and documentation. Several labs that employ smFRET approaches have joined forces to inform the smFRET community about streamlining how to perform experiments and analyze results for obtaining quantitative information on biomolecular structure and dynamics. The recent efforts include blind tests to assess the accuracy and the precision of smFRET experiments among different labs using various procedures. These multi-lab studies have led to the development of smFRET procedures and documentation, which are important when submitting entries into the archiving system for integrative structure models, PDB-Dev. This position paper describes the current 'state of the art' from different perspectives, points to unresolved methodological issues for quantitative structural studies, provides a set of 'soft recommendations' about which an emerging consensus exists, and lists openly available resources for newcomers and seasoned practitioners. To make further progress, we strongly encourage 'open science' practices.

Resolving fluorescent species by their brightness and diffusion using correlated photon-counting histograms

Fluorescence fluctuation spectroscopy (FFS) refers to techniques that analyze fluctuations in the fluorescence emitted by fluorophores diffusing in a small volume and can be used to distinguish between populations of molecules that exhibit differences in brightness or diffusion. For example, fluorescence correlation spectroscopy (FCS) resolves species through their diffusion by analyzing correlations in the fluorescence over time; photon counting histograms (PCH) and related methods based on moment analysis resolve species through their brightness by analyzing fluctuations in the photon counts. Here we introduce correlated photon counting histograms (cPCH), which uses both types of information to simultaneously resolve fluorescent species by their brightness and diffusion. We define the cPCH distribution by the probability to detect both a particular number of photons at the current time and another number at a later time. FCS and moment analysis are special cases of the moments of the cPCH distribution, and PCH is obtained by summing over the photon counts in either channel. cPCH is inherently a dual channel technique, and the expressions we develop apply to the dual colour case. Using simulations, we demonstrate that two species differing in both their diffusion and brightness can be better resolved with cPCH than with either FCS or PCH. Further, we show that cPCH can be extended both to longer dwell times to improve the signal-to-noise and to the analysis of images. By better exploiting the information available in fluorescence fluctuation spectroscopy, cPCH will be an enabling methodology for quantitative biology.

Global analysis of autocorrelation functions and photon counting distributions in fluorescence fluctuation spectroscopy

In fluorescence correlation spectroscopy (FCS) and photon counting histogram (PCH) analysis, the same experimental fluorescence intensity fluctuations are used, but each analytical method focuses on a different property of the signal. The time-dependent decay of the correlation of fluorescence fluctuations is measured in FCS yielding molecular diffusion coefficients and triplet-state parameters such as fraction and decay time. The amplitude distribution of these fluctuations is calculated by PCH analysis yielding the molecular brightness. Both FCS and PCH give information about the molecular concentration. Here we describe a global analysis protocol that simultaneously recovers relevant and common parameters in model functions of FCS and PCH from a single fluorescence fluctuation trace. Application of a global analysis approach allows increasing the information content available from a single measurement that results in more accurate values of molecular diffusion coefficients and triplet-state parameters and also in robust, time-independent estimates of molecular brightness and number of molecules.

The Stress-response protein prostate-associated gene 4, interacts with c-Jun and potentiates its transactivation

The Cancer/Testis Antigen (CTA), Prostate-associated Gene 4 (PAGE4), is a stress-response protein that is upregulated in prostate cancer (PCa) especially in precursor lesions that result from inflammatory stress. In cells under stress, translocation of PAGE4 to mitochondria increases while production of reactive oxygen species decreases. Furthermore, PAGE4 is also upregulated in human fetal prostate, underscoring its potential role in development. However, the proteins that interact with PAGE4 and the mechanisms underlying its pleiotropic functions in prostatic development and disease remain unknown. Here, we identified c-Jun as a PAGE4 interacting partner. We show that both PAGE4 and c-Jun are overexpressed in the human fetal prostate; and in cell-based assays, PAGE4 robustly potentiates c-Jun transactivation. Single-molecule Förster resonance energy transfer experiments indicate that upon binding to c-Jun, PAGE4 undergoes conformational changes. However, no interaction is observed in presence of BSA or unilamellar vesicles containing the mitochondrial inner membrane diphosphatidylglycerol lipid marker cardiolipin. Together, our data indicate that PAGE4 specifically interacts with c-Jun and that, conformational dynamics may account for its observed pleiotropic functions. To our knowledge, this is the first report demonstrating crosstalk between a CTA and a proto-oncogene. Disrupting PAGE4/c-Jun interactions using small molecules may represent a novel therapeutic strategy for PCa.

Simultaneous diffusion and brightness measurements and brightness profile visualization from single fluorescence fluctuation traces of GFP in living cells

Fluorescence correlation spectroscopy (FCS) and photon-counting histogram (PCH) analysis use the same experimental fluorescence intensity fluctuations, but each analytical method focuses on a different property of the signal. The time-dependent decay of the correlation of fluorescence fluctuations is measured in FCS yielding, for instance, molecular diffusion coefficients. The amplitude distribution of these fluctuations is calculated by PCH analysis yielding information about the molecular brightness of fluorescent species. Analysis of both FCS and PCH results in the molecular concentration of the sample. Using a previously described global analysis procedure we report here precise, simultaneous measurements of diffusion constants and brightness values from single fluorescence fluctuation traces of green-fluorescent protein (GFP, S65T) in the cytoplasm of Dictyostelium cells. The use of a polynomial profile in PCH analysis, describing the detected three-dimensional shape of the confocal volume, enabled us to obtain well fitting results for GFP in cells. We could visualize the polynomial profile and show its deviation from a Gaussian profile.

Classic maximum entropy recovery of the average joint distribution of apparent FRET efficiency and fluorescence photons for single-molecule burst measurements

We describe a method for analysis of single-molecule Förster resonance energy transfer (FRET) burst measurements using classic maximum entropy. Classic maximum entropy determines the Bayesian inference for the joint probability describing the total fluorescence photons and the apparent FRET efficiency. The method was tested with simulated data and then with DNA labeled with fluorescent dyes. The most probable joint distribution can be marginalized to obtain both the overall distribution of fluorescence photons and the apparent FRET efficiency distribution. This method proves to be ideal for determining the distance distribution of FRET-labeled biomolecules, and it successfully predicts the shape of the recovered distributions.

The spectroscopic basis of fluorescence triple correlation spectroscopy

We have developed fluorescence triple correlation spectroscopy (F3CS) as an extension of the widely used fluorescence microscopy technique fluorescence correlation spectroscopy. F3CS correlates three signals at once and provides additional capabilities for the study of systems with complex stoichiometry, kinetic processes, and irreversible reactions. A general theory of F3CS was developed to describe the interplay of molecular dynamics and microscope optics, leading to an analytical function to predict experimental triple correlations of molecules that freely diffuse through the tight focus of the microscope. Experimental correlations were calculated from raw fluorescence data using triple correlation integrals that extend multiple-tau correlation theory to delay times in two dimensions. The quality of experimental data was improved by tuning specific spectroscopic parameters and employing multiple independent detectors to minimize optoelectronic artifacts. Experiments with the reversible system of freely diffusing 16S rRNA revealed that triple correlation functions contain symmetries predicted from time-reversal arguments. Irreversible systems are shown to break these symmetries, and correlation strategies were developed to detect time-reversal asymmetries in a comprehensive way with respect to two delay times, each spanning many orders of magnitude in time. The correlation strategies, experimental approaches, and theory developed here enable studies of the composition and dynamics of complex systems using F3CS.

Using molecular dynamics and quantum mechanics calculations to model fluorescence observables

We provide a critical examination of two different methods for generating a donor-acceptor electronic coupling trajectory from a molecular dynamics (MD) trajectory and three methods for sampling that coupling trajectory, allowing the modeling of experimental observables directly from the MD simulation. In the first coupling method we perform a single quantum-mechanical (QM) calculation to characterize the excited state behavior, specifically the transition dipole moment, of the fluorescent probe, which is then mapped onto the configuration space sampled by MD. We then utilize these transition dipoles within the ideal dipole approximation (IDA) to determine the electronic coupling between the probes that mediates the transfer of energy. In the second method we perform a QM calculation on each snapshot and use the complete transition densities to calculate the electronic coupling without need for the IDA. The resulting coupling trajectories are then sampled using three methods ranging from an independent sampling of each trajectory point (the independent snapshot method) to a Markov chain treatment that accounts for the dynamics of the coupling in determining effective rates. The results show that the IDA significantly overestimates the energy transfer rate (by a factor of 2.6) during the portions of the trajectory in which the probes are close to each other. Comparison of the sampling methods shows that the Markov chain approach yields more realistic observables at both high and low FRET efficiencies. Differences between the three sampling methods are discussed in terms of the different mechanisms for averaging over structural dynamics in the system. Convergence of the Markov chain method is carefully examined. Together, the methods for estimating coupling and for sampling the coupling provide a mechanism for directly connecting the structural dynamics modeled by MD with fluorescence observables determined through FRET experiments.

Quantitative single-molecule imaging by confocal laser scanning microscopy

A new approach to quantitative single-molecule imaging by confocal laser scanning microscopy (CLSM) is presented. It relies on fluorescence intensity distribution to analyze the molecular occurrence statistics captured by digital imaging and enables direct determination of the number of fluorescent molecules and their diffusion rates without resorting to temporal or spatial autocorrelation analyses. Digital images of fluorescent molecules were recorded by using fast scanning and avalanche photodiode detectors. In this way the signal-to-background ratio was significantly improved, enabling direct quantitative imaging by CLSM. The potential of the proposed approach is demonstrated by using standard solutions of fluorescent dyes, fluorescently labeled DNA molecules, quantum dots, and the Enhanced Green Fluorescent Protein in solution and in live cells. The method was verified by using fluorescence correlation spectroscopy. The relevance for biological applications, in particular, for live cell imaging, is discussed.

Fluctuations as a source of information in fluorescence microscopy

Fluctuations in fluorescence spectroscopy and microscopy have traditionally been regarded as noise—they lower the resolution and contrast and do not permit high acquisition rates. However, fluctuations can also be used to gain additional information about a system. This fact has been exploited in single-point microscopic techniques, such as fluorescence correlation spectroscopy and analysis of single molecule trajectories, and also in the imaging field, e.g. in spatio-temporal image correlation spectroscopy. Here, we discuss how fluctuations are used to obtain more quantitative information from the data than that given by average values, while minimizing the effects of noise due to stochastic photon detection.

Concentration effects in "single-molecule" spectroscopy

A concentration as low as 0.1 molecule per observation volume may not be small enough for single-molecule Förster resonance energy transfer (FRET) efficiency measurements of molecules diffusing through a laser spot. This result follows from a rigorous theory that takes many molecules into account. We consider the distributions of the number of photons (photon counting histograms) and show that multiple-molecule effects are pronounced at large photon counts even at low concentrations. FRET efficiency distributions reveal multiple-molecule effects at large threshold values. This might be misinterpreted as multiple conformational states. Multiple-molecule effects strongly depend on the brightness of fluorophores. A simple test is suggested to determine parameters for which the single-molecule description is applicable.

Single-molecule FRET with diffusion and conformational dynamics

Under relatively mild conditions, we show how one can extract information about conformational dynamics from Förster resonance energy transfer (FRET) experiments on diffusing molecules without modeling diffusion. Starting from a rigorous theory that does treat diffusion, we first examine when the single-molecule FRET efficiency distribution can be decomposed into the measured distribution of the total number of photons and the efficiency distribution of an immobilized molecule in the absence of shot noise. If the conformation does not change during the time the molecule spends in the laser spot, this is possible when (I) the efficiency is independent of the location in the laser spot and (II) the total number of photons does not depend on conformation. This decomposition is approximate when the conformation changes during the diffusion time. However, it does provide a simple framework for analyzing data. This is illustrated for a two-state system where the FRET efficiency distribution can be found analytically for all values of the interconversion rates. If the arrival time of each donor and acceptor photon can be monitored, we introduce an alternative procedure that allows one to rigorously extract the rates of conformational changes when the above two conditions hold. In this case, the pattern of colors in the photon trajectory depends solely on conformational dynamics. This can be exploited in the framework of statistical inference because the likelihood function, which must be optimized with respect to the model rate parameters, depends only on how the conformation changes during the interval between photons with specified colors.

A Comparison between Photon Counting Histogram and Fluorescence Intensity Distribution Analysis

Photon counting histogram (PCH) and fluorescence intensity distribution analysis (FIDA) are two methods that were developed independently but reported almost simultaneously. Both of them have been successfully applied to fluorescence fluctuation spectroscopy (FFS). Though publications have indicated that they are theoretically equivalent, they are still commonly considered as different methods, especially in their ways to treat the point spread functions (PSFs). In this paper, the two methods are examined in detail for a direct comparison. After a direct proof of the theoretical equivalence, the authors further point out that PCH and FIDA are completely equivalent in the way of modeling PSFs; that is, any modeling approach developed from one of them can always be applied to the other. It is also demonstrated that simplified FIDA and PCH formulas in the form of power series can be applied for fast and precise numerical calculations. The two methods are also compared for their merits in the calculation efficiency.

Dual-color time-integrated fluorescence cumulant analysis

We introduce dual-color time-integrated fluorescence cumulant analysis (TIFCA) to analyze fluorescence fluctuation spectroscopy data. Dual-color TIFCA utilizes the bivariate cumulants of the integrated fluorescent intensity from two detection channels to extract the brightness in each channel, the occupation number, and the diffusion time of fluorophores simultaneously. Detecting the fluorescence in two detector channels introduces the possibility of differentiating fluorophores based on their fluorescence spectrum. We derive an analytical expression for the bivariate factorial cumulants of photon counts for arbitrary sampling times. The statistical accuracy of each cumulant is described by its variance, which we calculate by the moments-of-moments technique. A method that takes nonideal detector effects such as dead-time and afterpulsing into account is developed and experimentally verified. We perform dual-color TIFCA analysis on simple dye solutions and a mixture of dyes to characterize the performance and accuracy of our theory. We demonstrate the robustness of dual-color TIFCA by measuring fluorescent proteins over a wide concentration range inside cells. Finally we demonstrate the sensitivity of dual-color TIFCA by resolving EGFP/EYFP binary mixtures in living cells with a single measurement.