A method for single molecule tracking using a conventional single-focus confocal setup.

One way to achieve spatial resolution using fluorescence imaging-and track single molecules-is to use wide-field illumination and collect measurements over multiple sensors (camera pixels). Here we propose another way that uses confocal measurements and a single sensor. Traditionally, confocal microscopy has been used to achieve high temporal resolution at the expense of spatial resolution. This is because it utilizes very few, and commonly just one, sensors to collect data. Yet confocal data encode spatial information. Here we show that non-uniformities in the shape of the confocal excitation volume can be exploited to achieve spatial resolution. To achieve this, we formulate a specialized hidden Markov model and adapt a forward filtering-backward sampling Markov chain Monte Carlo scheme to efficiently handle molecular motion within a symmetric confocal volume characteristically used in fluorescence correlation spectroscopy. Our method can be used for single confocal volume applications or incorporated into larger computational schemes for specialized, multi-confocal volume, optical setups.

Dynamic 1D search and processive nucleosome translocations by RSC and ISW2 chromatin remodelers.

Eukaryotic gene expression is linked to chromatin structure and nucleosome positioning by ATP-dependent chromatin remodelers that establish and maintain nucleosome-depleted regions (NDRs) near transcription start sites. Conserved yeast RSC and ISW2 remodelers exert antagonistic effects on nucleosomes flanking NDRs, but the temporal dynamics of remodeler search, engagement, and directional nucleosome mobilization for promoter accessibility are unknown. Using optical tweezers and two-color single-particle imaging, we investigated the Brownian diffusion of RSC and ISW2 on free DNA and sparse nucleosome arrays. RSC and ISW2 rapidly scan DNA by one-dimensional hopping and sliding, respectively, with dynamic collisions between remodelers followed by recoil or apparent co-diffusion. Static nucleosomes block remodeler diffusion resulting in remodeler recoil or sequestration. Remarkably, both RSC and ISW2 use ATP hydrolysis to translocate mono-nucleosomes processively at ~30 bp/s on extended linear DNA under tension. Processivity and opposing push-pull directionalities of nucleosome translocation shown by RSC and ISW2 shape the distinctive landscape of promoter chromatin.

Dynamic 1D Search and Processive Nucleosome Translocations by RSC and ISW2 Chromatin Remodelers.

Eukaryotic gene expression is linked to chromatin structure and nucleosome positioning by ATP-dependent chromatin remodelers that establish and maintain nucleosome-depleted regions (NDRs) near transcription start-sites. Conserved yeast RSC and ISW2 remodelers exert antagonistic effects on nucleosomes flanking NDRs, but the temporal dynamics of remodeler search, engagement and directional nucleosome mobilization for promoter accessibility are unknown. Using optical tweezers and 2-color single-particle imaging, we investigated the Brownian diffusion of RSC and ISW2 on free DNA and sparse nucleosome arrays. RSC and ISW2 rapidly scan DNA by one-dimensional hopping and sliding respectively, with dynamic collisions between remodelers followed by recoil or apparent co-diffusion. Static nucleosomes block remodeler diffusion resulting in remodeler recoil or sequestration. Remarkably, both RSC and ISW2 use ATP hydrolysis to translocate mono-nucleosomes processively at ~30 bp/sec on extended linear DNA under tension. Processivity and opposing push-pull directionalities of nucleosome translocation shown by RSC and ISW2 shape the distinctive landscape of promoter chromatin.

Single-Molecule Reaction-Diffusion.

We propose to capture reaction-diffusion on a molecule-by-molecule basis from the fastest acquirable timescale, namely individual photon arrivals. We illustrate our method on intrinsically disordered human proteins, the linker histone H1.0 as well as its chaperone prothymosin α, as these diffuse through an illuminated confocal spot and interact forming larger ternary complexes on millisecond timescales. Most importantly, single-molecule reaction-diffusion, smRD, reveals single molecule properties without trapping or otherwise confining molecules to surfaces. We achieve smRD within a Bayesian paradigm and term our method Bayes-smRD. Bayes-smRD is further free of the average, bulk, results inherent to the analysis of long photon arrival traces by fluorescence correlation spectroscopy. In learning from thousands of photon arrivals continuous spatial positions and discrete conformational and photophysical state changes, Bayes-smRD estimates kinetic parameters on a molecule-by-molecule basis with two to three orders of magnitude less data than tools such as fluorescence correlation spectroscopy thereby also dramatically reducing sample photodamage.

Building Fluorescence Lifetime Maps Photon-by-Photon by Leveraging Spatial Correlations.

Fluorescence lifetime imaging microscopy (FLIM) has become a standard tool in the quantitative characterization of subcellular environments. However, quantitative FLIM analyses face several challenges. First, spatial correlations between pixels are often ignored as signal from individual pixels is analyzed independently thereby limiting spatial resolution. Second, existing methods deduce photon ratios instead of absolute lifetime maps. Next, the number of fluorophore species contributing to the signal is unknown, while excited state lifetimes with <1 ns difference are difficult to discriminate. Finally, existing analyses require high photon budgets and often cannot rigorously propagate experimental uncertainty into values over lifetime maps and number of species involved. To overcome all of these challenges simultaneously and self-consistently at once, we propose the first doubly nonparametric framework. That is, we learn the number of species (using Beta-Bernoulli process priors) and absolute maps of these fluorophore species (using Gaussian process priors) by leveraging information from pulses not leading to observed photon. We benchmark our framework using a broad range of synthetic and experimental data and demonstrate its robustness across a number of scenarios including cases where we recover lifetime differences between species as small as 0.3 ns with merely 1000 photons.

High Resolution Fluorescence Lifetime Maps from Minimal Photon Counts.

Fluorescence lifetime imaging microscopy (FLIM) may reveal subcellular spatial lifetime maps of key molecular species. Yet, such a quantitative picture of life necessarily demands high photon budgets at every pixel under the current analysis paradigm, thereby increasing acquisition time and photodamage to the sample. Motivated by recent developments in computational statistics, we provide a direct means to update our knowledge of the lifetime maps of species of different lifetimes from direct photon arrivals, while accounting for experimental features such as arbitrary forms of the instrument response function (IRF) and exploiting information from empty laser pulses not resulting in photon detection. Our ability to construct lifetime maps holds for arbitrary lifetimes, from short lifetimes (comparable to the IRF) to lifetimes exceeding interpulse times. As our method is highly data efficient, for the same amount of data normally used to determine lifetimes and photon ratios, working within the Bayesian paradigm, we report direct blind unmixing of lifetimes with subnanosecond resolution and subpixel spatial resolution using standard raster scan FLIM images. We demonstrate our method using a wide range of simulated and experimental data.

Computational Proposal for Tracking Multiple Molecules in a Multifocus Confocal Setup.

Tracking single molecules continues to provide new insights into the fundamental rules governing biological function. Despite continued technical advances in fluorescent and non-fluorescent labeling as well as data analysis, direct observations of trajectories and interactions of multiple molecules in dense environments remain aspirational goals. While confocal methods provide a means to deduce dynamical parameters with high temporal resolution, such as diffusion coefficients, they do so at the expense of spatial resolution. Indeed, on account of a confocal volume's symmetry, typically only distances from the center of the confocal spot can be deduced. Motivated by the need for true three dimensional high speed tracking in densely labeled environments, we propose a computational tool for tracking many fluorescent molecules traversing multiple, closely spaced, confocal measurement volumes providing independent observations. Various realizations of this multiple confocal volumes strategy have previously been used for long term, large area, tracking of one fluorescent molecule in three dimensions. What is more, we achieve tracking by directly using single photon arrival times to inform our likelihood and exploit Hamiltonian Monte Carlo to efficiently sample trajectories from our posterior within a Bayesian nonparametric paradigm. A nonparametric paradigm here is warranted as the number of molecules present are, themselves, a priori unknown. Taken together, we provide a computational framework to infer trajectories of multiple molecules at once, below the diffraction limit (the width of a confocal spot), in three dimensions at sub-millisecond or faster time scales.

Direct Photon-by-Photon Analysis of Time-Resolved Pulsed Excitation Data using Bayesian Nonparametrics.

Lifetimes of chemical species are typically estimated by either fitting time-correlated single-photon counting (TCSPC) histograms or phasor analysis from time-resolved photon arrivals. While both methods yield lifetimes in a computationally efficient manner, their performance is limited by choices made on the number of distinct chemical species contributing photons. However, the number of species is encoded in the photon arrival times collected for each illuminated spot and need not be set by hand a priori. Here, we propose a direct photon-by-photon analysis of data drawn from pulsed excitation experiments to infer, simultaneously and self-consistently, the number of species and their associated lifetimes from a few thousand photons. We do so by leveraging new mathematical tools within the Bayesian nonparametric. We benchmark our method for both simulated and experimental data for 1-4 species.

Pitching single-focus confocal data analysis one photon at a time with Bayesian nonparametrics.

Fluorescence time traces are used to report on dynamical properties of molecules. The basic unit of information in these traces is the arrival time of individual photons, which carry instantaneous information from the molecule, from which they are emitted, to the detector on timescales as fast as microseconds. Thus, it is theoretically possible to monitor molecular dynamics at such timescales from traces containing only a sufficient number of photon arrivals. In practice, however, traces are stochastic and in order to deduce dynamical information through traditional means-such as fluorescence correlation spectroscopy (FCS) and related techniques-they are collected and temporally autocorrelated over several minutes. So far, it has been impossible to analyze dynamical properties of molecules on timescales approaching data acquisition without collecting long traces under the strong assumption of stationarity of the process under observation or assumptions required for the analytic derivation of a correlation function. To avoid these assumptions, we would otherwise need to estimate the instantaneous number of molecules emitting photons and their positions within the confocal volume. As the number of molecules in a typical experiment is unknown, this problem demands that we abandon the conventional analysis paradigm. Here, we exploit Bayesian nonparametrics that allow us to obtain, in a principled fashion, estimates of the same quantities as FCS but from the direct analysis of traces of photon arrivals that are significantly smaller in size, or total duration, than those required by FCS.

An alternative framework for fluorescence correlation spectroscopy.

Fluorescence correlation spectroscopy (FCS), is a widely used tool routinely exploited for in vivo and in vitro applications. While FCS provides estimates of dynamical quantities, such as diffusion coefficients, it demands high signal to noise ratios and long time traces, typically in the minute range. In principle, the same information can be extracted from microseconds to seconds long time traces; however, an appropriate analysis method is missing. To overcome these limitations, we adapt novel tools inspired by Bayesian non-parametrics, which starts from the direct analysis of the observed photon counts. With this approach, we are able to analyze time traces, which are too short to be analyzed by existing methods, including FCS. Our new analysis extends the capability of single molecule fluorescence confocal microscopy approaches to probe processes several orders of magnitude faster and permits a reduction of photo-toxic effects on living samples induced by long periods of light exposure.