Summary

Automated Two-dimensional Spatiotemporal Analysis of Mobile Single-molecule FRET Probes

Published: November 23, 2021
doi:

Summary

This article presents a method for spatiotemporal analysis of mobile, single-molecule Förster resonance energy transfer (smFRET)-based probes using widefield fluorescence microscopy. The newly developed software toolkit allows the determination of smFRET time traces of moving probes, including the correct FRET efficiency and the molecular positions, as functions of time.

Abstract

Single-molecule Förster resonance energy transfer (smFRET) is a versatile technique reporting on distances in the sub-nanometer to nanometer range. It has been used in a wide range of biophysical and molecular biological experiments, including the measurement of molecular forces, characterization of conformational dynamics of biomolecules, observation of intracellular colocalization of proteins, and determination of receptor-ligand interaction times. In a widefield microscopy configuration, experiments are typically performed using surface-immobilized probes. Here, a method combining single-molecule tracking with alternating excitation (ALEX) smFRET experiments is presented, permitting the acquisition of smFRET time traces of surface-bound, yet mobile probes in plasma membranes or glass-supported lipid bilayers. For the analysis of recorded data, an automated, open-source software collection was developed supporting (i) the localization of fluorescent signals, (ii) single-particle tracking, (iii) determination of FRET-related quantities including correction factors, (iv) stringent verification of smFRET traces, and (v) intuitive presentation of the results. The generated data can conveniently be used as input for further exploration via specialized software, e.g., for the assessment of the diffusional behavior of probes or the investigation of FRET transitions.

Introduction

Förster resonance energy transfer (FRET) has been a major driver in molecular biological and biophysical research, as it allows the investigation of processes at sub-nanometer resolution. As the efficiency of energy transfer between donor and acceptor fluorophores strongly depends on the inter-dye distance in the sub-nanometer to nanometer range, it has been effectively used as a spectroscopic ruler to explore static and dynamic conformation of biomolecules1,2,3,4. Additionally, the FRET phenomenon has been widely used for colocalization studies of membrane-associated and intracellular proteins on a bulk level5,6. In the last two decades, the method was adapted for monitoring smFRET events7, which helped to substantially increase temporal and spatial resolution and resolved even rare subpopulations in heterogeneous samples. Equipped with these techniques, unique insights were gained into the dynamics of molecular machinery such as the transcript processing rate of RNA polymerase II8, replication speed of DNA polymerases9,10, nucleosome translocation rate11, transcript splicing and stalling rate of assembled spliceosomes12, the activity of ribosomal subpopulations13, and the walking speed of kinesin motors14, to name a few. Receptor-ligand interaction durations15 and molecular forces16 have been quantified.

Intensity-based smFRET studies typically rely on sensitized emission to measure FRET efficiency: a beam splitter in the emission path spatially separates light originating from donor and acceptor fluorophores upon donor excitation, allowing for the quantification of individual fluorescence intensities. The efficiency can subsequently be calculated as the fraction of photons emitted by the acceptor with respect to the total photon count17. In addition, acceptor excitation following donor excitation (ALEX) permits measurement of the stoichiometry of the FRET events, aiding in the discrimination between true low FRET signals from signals arising, e.g., from probes featuring a photobleached acceptor fluorophore18.

Single-molecule FRET experiments are commonly carried out in one of two ways. First, a small region in the sample volume is illuminated using a confocal microscope. Single probe molecules in solution are excited when they happen to diffuse within the focal volume. With this technique, fast photon-counting detectors can be used, enabling sub-microsecond time resolution. Second, probes are specifically immobilized on surfaces and monitored via widefield microscopy, often using total internal reflection (TIR) configuration to minimize background fluorescence. Probe immobilization allows for much longer recording times than using the first approach. In addition, the larger field of view permits the monitoring of multiple probes in parallel. The need for a camera makes this method slow compared to the one described above. Time resolution is limited to the millisecond to second range.

If long time traces are required, e.g., for studying dynamic processes on a millisecond to second time scale, the first method is not applicable, as the fluorescence bursts are typically too short. The second approach fails whenever immobilization is not feasible, e.g., in live-cell experiments featuring probes diffusing within the cell membrane. Furthermore, it has been observed that biological model systems can vary their response dramatically depending on the mobility of the contacted surface16.

While combined smFRET and single-particle tracking experiments recording mobile FRET probes have been performed in the past19, there is no publicly available software for the evaluation of the data. This prompted the development of a new analysis platform, which allows for the determination of multiple properties of mobile fluorescent probes, including smFRET efficiency and stoichiometry, positions with sub-pixel accuracy, and fluorescence intensities as functions of time. Methods for filtering the resulting traces by examining stepwise bleaching behavior, nearest-neighbor distances, emission intensities, and other traits were established to exclusively choose correctly synthesized and functional single-probe molecules. The software also supports experimental and analytical techniques recently agreed upon in a multilaboratory study to produce reliable, quantitative smFRET data17. In particular, the implementation adheres to the validated procedures for the calculation of FRET efficiency and stoichiometry. Fluorescence intensities upon donor excitation in the donor emission channel IDD and acceptor emission channel IDA are used for the calculation of the apparent FRET efficiency Eapp using Eq (1).

Equation 1 (1)

With the help of the fluorescence intensity in the acceptor emission channel upon acceptor excitation IAA, the apparent stoichiometry is calculated using Eq (2).

Equation 2 (2)

The FRET efficiency E and the stoichiometry S can be derived from Eapp and Sapp by considering four correction factors.

Equation 3

α describes the leakage of donor fluorescence into the acceptor emission channel and can be determined using a sample containing only donor fluorophores or by analyzing parts of trajectories where the acceptor has been bleached. δ corrects for the direct excitation of the acceptor by the donor excitation light source and can be measured using a sample with only acceptor fluorophores or by analyzing parts of trajectories where the donor has been bleached.

Equation 4.

γ scales IDD to rectify diverging detection efficiencies in donor and acceptor emission channels and different quantum efficiencies of the fluorophores. The factor can be computed by analyzing the increase in donor intensity upon acceptor bleaching in trajectories with high FRET efficiencies20 or by studying a sample featuring multiple discrete FRET states.

Equation 5

β scales IAA to correct for disparate efficiencies of donor and acceptor excitation. If γ was determined via acceptor bleaching analysis, β could be calculated from a sample of known donor-to-acceptor ratio21. Otherwise, the multi-state FRET sample also yields β.

Together, the corrections allow the calculation of the corrected FRET efficiency using Eq (3).

Equation 6 (3)

and the corrected stoichiometry using Eq (4).

Equation 7 (4)

Ideally, the corrected stoichiometry for a 1:1 donor-to-acceptor ratio gives S = 0.5. In practice, a reduced signal-to-noise ratio produces a spread of the measured values of S, hampering the discrimination from donor-only signals (S = 1) and acceptor-only signals (S = 0). The resulting time traces can be used as input for a more detailed analysis of the single-molecule trajectories to obtain information such as spatiotemporal force profiles16, the mobility of the single-molecule events22, or transition kinetics between different states1.

The following protocol describes experimental parameters and procedures for smFRET tracking experiments, as well as the working principle behind data analysis using the newly developed software suite. For the acquisition of experimental data, it is recommended to use a microscopy setup meeting the following requirements: i) capability of detecting the emission of single dye molecules; ii) widefield illumination: in particular for live-cell experiments, total internal reflection (TIR23,24,25) configuration is recommended; iii) spatial separation of emission light according to wavelength such that donor and acceptor fluorescence is projected onto different regions of the same camera chip25 or different cameras; iv) modulation of light sources for donor and acceptor excitation with millisecond precision, e.g., using directly modulatable lasers or modulation via acousto-optic modulators. This permits stroboscopic illumination to minimize photobleaching of fluorophores as well as alternating excitation to determine stoichiometries; v) output of one file per recorded image sequence in a format that can be read by the PIMS Python package26. In particular, multipage TIFF files are supported.

Protocol

1. Software prerequisites

  1. Install the miniconda Python distribution27 (minimum required Python version: 3.7).
  2. Open an Anaconda prompt in the Windows Start menu, or open a terminal and execute conda activate if using Linux or macOS.
  3. Enable the community-maintained conda-forge package repository28 by executing the following commands:
    conda config –add channels conda-forge
    conda config –set channel_priority strict
    conda update –all
  4. Install the required Python packages by executing:
    conda install opencv trackpy lmfit ipympl scikit-learn pyqt sdt-python jupyterlab
  5. Become familiar with JupyterLab, the user interface of the analysis software (refer to the software documentation29).
  6. Install the git version control system, which will be used later to download and update the analysis software. If using Linux, use the distribution's package management software to download and update. Otherwise, execute:
    conda install git
  7. Optionally install the sidecar Python package to display datasets after filtering steps during analysis:
    conda install sidecar

2. Measurement of samples

Figure 1
Figure 1: Image Acquisition. (A) Excitation sequence. After recording an optional image of a dye-loaded cell using the 405 nm laser, donor and acceptor are excited alternately and repeatedly for illumination time till using 532 nm and 640 nm lasers, respectively. The time tbetween donor and acceptor excitation must be long enough to allow for image readout by the camera. The delay time tdelay can be used to adjust the acquisition frame rate and, therefore, the observation time span before photobleaching. This panel is modified from 16. (B) Fiducial markers are used for the calculation of the coordinate transforms between the two emission channels. Matching fiducials are indicated by color. Several shifted images should be recorded to ensure that the whole field of view is covered. (C) Laser profiles for flatfield correction are recorded using a densely labeled sample. The acceptor profile is recorded and photobleached, followed by acquisition of the donor profile. Multiple images should be taken at different sample regions, averaged, and smoothed to mitigate the influence of sample imperfections (e.g., the bright spot in the center-top of the image). (D) Flatfield correction map p(x,y) calculated from 20 laser profile recorded as described in C. Abbreviations: FRET = Förster resonance energy transfer; ImDD = donor emission image upon donor excitation; ImDA = acceptor emission image upon donor excitation; ImAA = acceptor emission image upon donor excitation. Scale bars = 5 µm. Please click here to view a larger version of this figure.

  1. When using an electron-multiplying charge-coupled device (EMCCD) camera, enable the EM gain to observe single-molecule signals at high signal-to-noise ratios (refer to the manufacturer's instructions).
  2. Excitation sequence (see Figure 1A for further details).
    1. Optionally record an image for segmentation to restrict data analysis to certain regions in the field of view. For instance, excite Fura-2-loaded cells using a 405 nm laser and capture their emission around 510 nm to evaluate only probes located in interfaces between cells and supported lipid bilayers (SLBs). Consequently, wait for time tr to permit camera readout.
      NOTE: On EMCCD cameras, tr depends on the number of lines in the chosen region of interest (ROI). Therefore, choosing a small ROI can be advantageous because it reduces the delay between frames and the size of recorded data. Additionally, enabling frame transfer mode allows for further reduction of tr.
    2. Alternately excite donor and acceptor fluorophores repeatedly.
      1. Excite the donor for an illumination time till (5-10 ms is typically short enough to avoid motion blur) while also triggering the camera.
      2. Wait for time tr to permit camera readout.
      3. Excite the acceptor for while triggering the camera.
      4. Wait for a time tdelay.
        NOTE: This must be longer than tr to enable readout by the camera but can otherwise be chosen arbitrarily. It shall balance the requirements for time resolution and trace length.
      5. Repeat steps 2.2.2.1-2.2.2.4. Choose the number of repeats to be large enough to ensure photobleaching of at least one fluorophore per probe within the field of view, which permits stepwise photobleaching analysis for discrimination of single-molecule signals from aggregates.
        NOTE: Choosing appropriate till and excitation laser intensities commonly requires some experimentation: The longer the illumination times and the higher the laser intensities, the better is the signal-to-noise ratio in the resulting images, but the shorter are the resulting time traces.
  3. Record a sufficient number of movies for each sample.

3. Additional measurements for the determination of correction factors

  1. Record a series of randomly placed fiducial markers visible in both emission channels for image registration (i.e., finding the transformation that maps coordinates of the donor emission channel onto the acceptor emission channel and vice versa). See Figure 1B.
    NOTE: Image registration is performed by the software; see step 6.1.4.
  2. Measure the intensity profile for both donor and acceptor excitation light sources for flatfield correction (i.e., correcting for inhomogeneous excitation across the field of view). To this end, prepare a sample featuring a high density of FRET probes and first acquire an image upon acceptor excitation, followed by photobleaching of the acceptor and subsequent recording of an image upon donor excitation. For increased stability, repeat several times in different sample regions. See Figure 1C,D. Alternatively, record a sample decorated with only the donor molecule and a second sample decorated with only the acceptor fluorophores.
    NOTE: Flatfield correction is performed by the analysis software; see step 8.1.2.
  3. Record a single-molecule sample (as in section 2) of a probe without an acceptor fluorophore to determine donor emission leaking into the acceptor channel.
    NOTE: Donor leakage can also be computed from the actual probes' time traces after acceptor bleaching. If a sufficient number of such events is recorded, no additional measurement is necessary. Both options are supported by the analysis software; see Supplemental Information, section 3.15.
  4. Acquire recordings of a probe without a donor fluorophore for the quantification of direct acceptor excitation by the donor excitation light source.
    NOTE: Direct acceptor excitation can also be derived from the actual probes' time traces after donor bleaching. If a sufficient number of such events is recorded, no additional measurement is necessary. Both options are supported by the analysis software; see Supplemental Information, section 3.15.
  5. Record a single-molecule sample featuring two distinct FRET efficiencies to correct for differing detection efficiencies of the donor and acceptor emission channels and differing quantum yields of the dyes.
    NOTE: Such samples could be, for instance, Holliday junctions1, which fluctuate between two conformations, or DNA rods that have FRET pairs attached at different, well-defined distances. If probes feature high and sufficiently constant FRET efficiencies, the correction can also be computed from acceptor bleaching events of probes' time traces, in which case no additional measurements are necessary. Both options are supported by the analysis software; see Supplemental Information, section 3.15.

4. Single-molecule localization algorithms

NOTE: Several analysis steps require single-molecule localization. Choose between a Gaussian fitting algorithm30 and center-of-mass computation31, depending on signal density, background, and signal-to-noise ratio.

  1. To perform Gaussian fitting, choose the 3D-DAOSTORM30 algorithm via the respective user interfaces.
    NOTE: 3D-DAOSTORM is designed to distinguish even signals with overlapping point spread functions. While this is generally an advantage, it comes with a caveat: single, bright signals are occasionally identified as two adjacent ones, which can confuse the tracking algorithm and result in the detection of two short trajectories instead of a single long one.
    Set the following parameters (for details, see the documentation of the sdt-python library32, which provides the algorithm's implementation).
    1. radius: Set the initial σ value of the Gaussian fit function in pixels depending on the effective pixel size.
    2. threshold: Set a minimum amplitude (i.e., brightest pixel value, corrected for the estimated local background) for a local intensity maximum to be fit.
      NOTE: The threshold is arguably the most important parameter. If set too low, noise may be considered a fluorescence signal, and bright signals may be fitted with two Gaussians. If set too high, dim signals are not fit.
    3. model: Set to 2d to fit circular Gaussians.
      NOTE: The other models are not applicable to the smFRET data.
    4. find filter: Apply a filter before finding local maxima to reduce noise, which is helpful in low signal-to-noise ratio situations. This can be i) identity: no filter; ii) Crocker-Grier: bandpass filter from Crocker-Grier algorithm31,33; or iii) Gaussian: a gaussian blur with σ set by the sigma parameter.
      NOTE: For Crocker-Grier, the feat. size parameter should be roughly the radius of a point spread function in pixels.
      NOTE: Fitting is performed using unfiltered raw data.
    5. min. distance: Fit two prospective signals separated by fewer than min. distance pixels by a single Gaussian.
      NOTE: This can help in the aforementioned scenario where a bright signal is wrongly detected as two adjacent signals.
    6. size range: Select minimum and maximum σ of the fits to remove detections from spurious signals due to noise.
  2. Choose the Crocker-Grier algorithm via the respective user interfaces to perform center-of-mass computation (a refined algorithm31 based on Crocker's and Grier's idea33).
    NOTE: This algorithm is very robust even in low signal-to-noise scenarios and in dealing with signals featuring a range of intensities but cannot fit molecules with overlapping point spread functions precisely.
    1. radius: Set the radius (in pixels) of a disk large enough to contain the whole point spread function.
    2. signal thresh.: Set the minimum amplitude (brightest pixel above estimated background) for a local intensity maximum to be analyzed.
      NOTE: If set too low, noise may be considered a fluorescence signal. If set too high, dim signals are not fit.
    3. mass thresh.: Set the minimum total intensity (sum of background-corrected pixel values) of a signal to be analyzed.
      NOTE: Same considerations as above apply.

5. Software initialization

  1. Download analysis scripts. In an Anaconda prompt, navigate to a folder to save the analysis (using the cd command) and execute
    git clone https://github.com/schuetzgroup/fret-analysis.git target folder
    1. Replace the target folder with a descriptive name such as 2021-06-14_Force-FRET-experiment.
      NOTE: The analysis software will end up in this folder; ensure that this folder does not exist beforehand. It is recommended to download a copy of the analysis scripts for each experiment. This way, it is possible to revisit the analysis later, recall the parameters used, and make changes.
  2. Copy Jupyter notebooks (01. Tracking.ipynb, 02. Analysis.ipynb, 03. Plots.ipynb) into the newly created folder (henceforth referred to as the root folder). If this is the first time using the software, get them from the notebooks subfolder of the root folder.
    NOTE: If similar datasets have already been analyzed, copying the notebooks from a previous experiment can be a convenient option, as parameters may have changed only slightly.
  3. Start the JupyterLab server by executing the following command in the Anaconda prompt to open a web browser window displaying JupyterLab.
    jupyter lab
    NOTE: The browser is only the interface, while the process running in the Anaconda prompt is doing the actual work. As a consequence, closing the browser window has only minimal effect; the session can be restored by accessing http://localhost:8888. However, interrupting the JupyterLab process in the prompt or closing the prompt will terminate the analysis, leading to the loss of unsaved work.
  4. In the JupyterLab browser window, use the left pane to navigate to the root folder. Double-click on 01. Tracking.ipynb to launch the first notebook. After launching, look for a new tab to appear, which displays boxes, so-called cells, of Python code.
    NOTE: All Jupyter notebooks feature comments describing the functionality of each code cell. Additionally, documentation for every method call can be displayed by placing the text cursor immediately before the opening ( and press Shift+Tab.
  5. See Figure 2 for an overview of the data analysis process.

Figure 2
Figure 2: Overview of a typical analysis pipeline. Note that filtering steps are subject to adaptation according to the experimental design. This figure is modified from 16. Abbreviation: FRET = Förster resonance energy transfer. Please click here to view a larger version of this figure.

NOTE: Sample data to try out the software can be downloaded from https://github.com/schuetzgroup/fret-analysis/releases/tag/example_files

6. Localization, tracking, and fluorescence intensity analysis of single molecules (01. Tracking.ipynb).

  1. Use the 01. Tracking.ipynb Jupyter notebook for the reliable analysis of fluorescence intensity values of single-molecule signals, which is tailored for the precise quantification, particularly of faint signals often occurring in FRET measurements (e.g., due to low donor signals at high FRET events and vice versa).
    NOTE: To this end, direct integration of the pixel intensities in the raw data with correction for the local background is implemented. For screenshots of each analysis step and description of function call parameters, refer to the Supplemental Information.
    1. Specify the illumination sequence to permit the selection of donor and acceptor excitation frames, as well as frames for image segmentation from recorded image sequences.
      NOTE: As the software allows the processing of data recorded with arbitrary illumination protocols, it is necessary to indicate which frame in an image sequence was acquired while exciting what type of fluorophore; see Supplemental Information, section 1.2, step 3. The frame numbers of the original image sequence are preserved.
    2. Describe and load datasets. Analyze multiple datasets at once, provided they were recorded using the same illumination settings. Assign an identifier and a pattern that matches the respective image sequence file names to each dataset. Additionally, define specific datasets for special purposes, such as recordings of fiducial markers for image registration, excitation light profiles for flatfield correction, and optionally donor-only and acceptor-only samples to determine correction factors.
    3. Select emission channels in raw images if both channels were recorded using a single camera. For this, use the appropriate graphical widget to select the appropriate regions for donor and acceptor emission.
    4. Localize fiducial markers in both emission channels and perform image registration. Use the provided user interface to find the appropriate parameters for the localization algorithm for both the donor and the acceptor emission channels. See section 4 for information about supported localization algorithms.
      NOTE: Randomly distributed fiducial markers can be identified across emission channels by the spatial distribution of their nearest neighbors (Figure 1B). A custom implementation of the algorithm proposed for selective plane illumination microscopy34 in the sdt-python library automatically matches the position of each marker in the donor emission channel with the position in the acceptor emission channel. A transformation T mapping the donor emission channel's coordinates to the acceptor emission channel's coordinates is found via a linear least-squares fit of an affine transformation to the markers' positions35. RANSAC is used to account for outliers, such as wrongly matched positions from the previous step.
    5. Localize FRET probes independently upon donor and acceptor excitation in all frames and merge results into one table that contains the original frame number, 2-dimensional coordinates, and an identifier referring to the source image file.
      NOTE: To this end, the software provides user interfaces to find appropriate options for the localization algorithm.
      1. Localize FRET probes upon donor excitation in the sum of the images obtained from donor emission ImDD and acceptor emission ImDA, which hardly depend on the FRET efficiency. For information concerning the options for the localization algorithm, see section 4.
        NOTE: Each sum image is calculated by transforming ImDD using the transformation T previously obtained from image registration and added pixelwise to ImDA.
      2. Localize probes upon acceptor excitation in the acceptor emission channel ImAA (see section 4 for details about the localization algorithms).
    6. Perform tracking and fluorescence intensity measurement.

Figure 3
Figure 3: Single-molecule intensity measurement. (A) For a fluorophore located at the orange pixel, its uncorrected intensity Iuncorr is determined by summing up all pixels' intensities within a disk (yellow and orange pixels) large enough to cover all pixels affected by the signal: Equation 9. The local background is computed as the mean of the pixels in a ring (blue pixels) around the disk: Equation 10, where nring is the number of pixels in the ring. The fluorescence intensity I is the result of subtracting the background from the uncorrected intensity, I = Iuncorrb × ndisk, where ndisk is the number of pixels in the disk. The circle radius is specified via the feat_radius parameter of the tracking method. The width of the ring is given by the bg_frame parameter. If the point spread function of one signal overlaps with the background ring of another (bottom panel), the affected pixels (red) are excluded from local background analysis. If two point spread functions overlap, fluorescence intensities cannot be calculated reliably and are therefore discarded. (B, C) Simulations show that applying a Gaussian blur with a standard deviation of 1 pixel improves the signal-to-noise ratio up to a factor of close to 2 at low fluorescence intensities (B) and introduces hardly any error (slight underestimation of less than 1%, (C))16. Moreover, the relative error (i.e., (ImeasItruth)/Itruth, where Itruth is the ground truth and Imeas is the outcome of the analysis) is constant over the whole intensity range and therefore cancels out for ratiometric quantities such as FRET efficiencies and stoichiometries. All plots are based on previously published work16. Abbreviations: SNR = signal-to-noise ratio; FRET = Förster resonance energy transfer. Please click here to view a larger version of this figure.

  1. Choose appropriate options for the trackpy36 algorithm used to link FRET probe localizations into trajectories. In particular, set the maximum search distance from one frame to the next and the number of consecutive frames for which a signal may go undetected, which can occur due to bleaching or missed localizations.
    NOTE: These gaps are filled via interpolation between preceding and ensuing positions. These interpolated positions are marked and only used later for reading out intensity values but not for diffusion analysis. Tracks are analyzed in the coordinate system of the acceptor emission channel. For fluorescence intensity analysis (step 6.1.6.2), tracks are additionally transformed to the coordinate system of the donor emission channel using the inverse transformation T-1 obtained via image registration (see step 6.1.4).
  2. Select options for the fluorescence intensity calculation algorithm (see Figure 3A for details). Specify i) the radius of a disk that, when centered on a signal's position, contains all pixels affected by that signal and ii) the width of a ring around each disk used to determine the local background.
    NOTE: To reduce noise in the obtained intensity measurements, a Gaussian blur with a standard deviation of 1 pixel is applied to the images (Figure 3B,C).
  1. Use the analysis software functionality to process auxiliary image data from image sequences.
    1. Extract additional images recorded to facilitate segmentation (see step 2.2.1, marked by s in the excitation sequence (see Supplemental Information, section 1.2, step 3).
    2. Determine donor and acceptor excitation light profiles across the field of view from images recorded on a densely labeled sample (see step 3.2).
      NOTE: The pixel-by-pixel mean is calculated from the images to compute the light profiles. The camera baseline is subtracted. Images are blurred using a Gaussian filter to reduce effects due to sample impurities. Finally, the resulting images are divided pixelwise by their maximum value to get the profile p(x,y) mapping coordinates onto the interval [0,1].

7. Visualization of FRET trajectories (optional)

  1. Use the inspector application to display single-molecule tracks in raw image data and the corresponding fluorescence intensities and apparent FRET efficiencies and stoichiometries.
    ​NOTE: This is a valuable tool to assess the validity of chosen parameters and manually accept or reject individual time traces. See Supplemental Information for a screenshot and detailed usage information.

8. Analysis and filtering of single-molecule data (02. Analysis.ipynb)

  1. Use the 02. Analysis.ipynb Jupyter notebook for analysis and filtering of the single-molecule data obtained via the 01. Tracking.ipynb notebook. See the steps below for a typical analysis pipeline.
    NOTE: Different scientific questions and experimental designs may require adjustments of the settings. The use of Jupyter notebooks permits easy adaption by omitting, rearranging, and amending analysis steps. For screenshots of each analysis step and description of function call parameters, refer to Supplemental Information.
    1. Perform initial filtering steps.
      1. Discard signals with overlapping point spread functions as it is difficult to determine their fluorescence intensities reliably.
      2. In the case of inhomogeneous illumination, accept only signals located in well-illuminated regions within the field of view to ensure a good signal-to-noise ratio.
      3. If studying intramolecular FRET, restrict the analysis to those trajectories present from the beginning of the image sequence to ensure that all bleaching steps are recorded and can be properly evaluated later during the stepwise photobleaching analysis.
        NOTE: When performing experiments with intermolecular FRET probes, wherein donor and acceptor fluorophores are not part of a preformed complex, it may not be feasible to restrict the analysis to initially present trajectories.
    2. Execute flatfield correction, which uses the excitation light source profiles obtained in step 6.1.7.2 to reverse the position-dependent fluorescence intensity variations caused by inhomogeneous illumination.
      NOTE: The fluorescence intensity I(x,y) of a probe at the position (x,y) is corrected via Equation 8 ; see Figure 1C,D.
    3. Compute the apparent FRET efficiency Eapp (i.e., the fraction of energy transmitted from the donor fluorophore to the acceptor fluorophore) and the apparent stoichiometry Sapp (i.e., the number of donor fluorophores divided by the total number of fluorophores within a diffraction-limited spot).
      NOTE: By plotting E vs. S for each data point, it is possible to distinguish alterations in the measured FRET efficiencies due to change in donor-acceptor distance from alterations due to changes in the stoichiometry18. This allows for the differentiation between E = 0 due to dye separation from E = 0 due to the absence of an active acceptor. E-S plots are used throughout the analysis as a tool for quality assessment; see Figure 4 as an example.
    4. Perform stepwise analysis of photobleaching for the discrimination between single-molecular probes and aggregates. Choose to accept one of the following options.
      NOTE: To this end, the analysis software applies a custom implementation32 of the changepoint detection algorithm PELT37 separately to the fluorescence intensity upon donor excitation (IDD + IDA) and acceptor excitation (IAA).
      1. Choose option 1, wherein the acceptor fluorophore bleaches in a single step while the donor shows no partial bleaching (i.e., there is no bleaching step to non-zero intensity).
        NOTE: This option further rejects trajectories where the donor bleaches before the acceptor in a single step. Option 1 is the preferred choice in case of high acceptor photobleaching rates.
      2. Choose option 2, wherein the donor bleaches in a single step while there is no partial acceptor bleaching.
        NOTE: This option further rejects trajectories where the donor bleaches after the acceptor in a single step. Option 2 is the preferred choice in case of high donor photobleaching rates.
      3. Choose option 3, wherein either fluorophore bleaches in a single step while the other does not partially bleach.
        NOTE: Option 3 gives higher flexibility than options 1 and 2 and would be the suggested preference for data analysis.
      4. Choose option 4, wherein donor and acceptor fluorophores show single-step photobleaching or no photobleaching at all.
        NOTE: Option 4 is preferred in case of low photobleaching rates.
    5. Calculate the correction factors for donor emission leakage into the acceptor channel α, direct acceptor excitation δ, detection efficiencies γ, and excitation efficiencies β17.
    6. Use the correction factors to calculate the FRET efficiency E from the apparent efficiency Eapp and the stoichiometry S from the apparent stoichiometry Sapp.
    7. Perform further filtering steps. Select only data points from before the first bleaching event in each trajectory. Additionally, accept only trajectories with at least 75% of data points satisfying 0.35 < S < 0.6 to restrict the analysis to single-molecule probes (numbers are adjustable).
      NOTE: The upper and lower bounds for should be chosen according to the spread of the population of interest versus the populations to be excluded from the analysis (e.g., donor-only and acceptor-only populations). Based on experience, 0.35 < S < 0.6 turned out to be a good choice for many experimental situations.
    8. Perform image segmentation via global or adaptive thresholding methods35 on the appropriate auxiliary images (see steps 2.2.1 and 6.1.7) to restrict the analysis to distinct regions within the field of view.
      ​NOTE: This permits, for instance, exclusive evaluation of probes located in a cell-SLB interface or on a patterned structure.

9. Plotting of results and further analysis (03. Plot.ipynb)

NOTE: Refer to Supplemental Information for screenshots of the Jupyter notebook and description of function call parameters.

  1. Create E-S plots to verify that signals of incorrect stoichiometry have been correctly identified and removed.
  2. Plot histograms of FRET efficiencies to provide a well-established overview of the FRET efficiency distributions. Group the histograms for convenient comparison of results from different experiments.
  3. Evaluate the data further (e.g., diffusion analysis, conversion of FRET efficiencies to forces in experiments using molecular force sensors or transition analysis) within the notebook taking advantage of scientific Python libraries.
    NOTE: Data can also be exported in many file formats as input to other analysis software.

Representative Results

A variety of low- and high-level information can be extracted from smFRET tracks depending on the scientific question of the experiment. Here, examples of analysis pipelines with analog and digital probes are presented: a peptide-based molecular force sensor16 and a DNA probe with stochastic switching of its conformation38, respectively. Refer to Figure 5 for the design and working principle of these probes.

After the localization and tracking algorithms have been executed as described in the protocol, the package offers multiple data visualization tools to optimize the initial parameters and subsequent filter steps: (i) visualization of individual smFRET events, (ii) optional image segmentation to analyze data in certain regions of interests, (iii) monitoring of filter steps via FRET efficiency vs. stoichiometry (E-S) plots. The visualization of the single-molecule data is presented in Figure 6.

Finally, the filtered FRET events are represented by an E-S plot and a FRET efficiency histogram (Figure 4). The E-S plot is a useful tool for optimizing the aforementioned filtering steps and investigating the final result. Partially bleached or incompletely labeled FRET sensors can be excluded by their stoichiometry value. Mobility parameters can be investigated by plotting an individual trajectory path in an x-y plot (Figure 6) or a mean square displacement (MSD) plot (Figure 4). The first method is especially useful for discriminating mobile from immobilized events, while the latter is used to calculate the diffusion coefficient.

Figure 4
Figure 4: Exemplary output. (A) The FRET efficiency is plotted versus stoichiometry (E-S plot) for a population of the molecular force sensor (left panel) decorating a glass-supported lipid bilayer and strained by a T cell. Only one population cloud is visible. The respective histogram of FRET efficiencies exemplifies the difference between a force sensor population in presence and absence of cells (middle panel). No shift to lower FRET efficiencies of the sensor population in presence of T cells can be observed, indicating little to no force-dependent stretching of the sensor module. The MSD plot of these experimental conditions confirms that the force sensor population beneath a T cell moves considerably slower than their unbound counterparts (right panel). (B) The same analysis was performed with Holliday junction DNA sensor decorating a glass-supported fluid lipid bilayer. The ES plot clearly shows two populations, which are also apparent in the FRET efficiency histogram. The MSD plot indicates the presence of one fast-moving sensor population. Abbreviations: FRET = Förster resonance energy transfer; MSD = mean square displacement. Please click here to view a larger version of this figure.

Figure 5
Figure 5: Design and working principle of intramolecular FRET probes. (A) Analog peptide sensor for quantification of mechanical molecular forces. The donor and acceptor fluorophores are covalently attached to either end of the peptide backbone. The sensor module is site-specifically attached to a specific ligand, which in turn binds a cell-resident surface receptor of interest (here, an antibody fragment specifically recognizing the beta chain of the T cell receptor). Upon receptor-ligand binding, force is exerted, and the sensor module extends and eventually recoils after bond cleavage. This panel is modified from 16. (B) Digital DNA sensor for quantification of FRET transitions. The FRET sensor is composed of four DNA strands forming a Holliday junction. The donor and acceptor fluorophore are covalently attached to two strands. Holliday junctions frequently switch their conformation depending on the surrounding buffer conditions. The stochastic switching of these conformations can be monitored by quantifying the FRET efficiency of individual probes. Abbreviations: TCR = T cell receptor; FRET = Förster resonance energy transfer. Please click here to view a larger version of this figure.

Figure 6
Figure 6: Examples of localization and tracking of FRET probes. (A) The FRET efficiency and stoichiometry of individual events are calculated by quantifying the intensity of the donor fluorophore upon donor excitation (D → D), the acceptor fluorophore upon donor excitation (D → A), and the acceptor fluorophore upon acceptor excitation (A → A). Nearest neighbor filtering prevents bias by overlapping point spread functions of close emitters. Image segmentation allows the user to choose certain smFRET events localized within an area of interest (e.g., a cell or a micropattern). As an example of image segmentation, T cells were stained with Fura-2 (displayed on the left) and subjected to adaptive thresholding to identify the cell edges (orange dotted line). Scale bars = 5 µm. (B) smFRET trajectories using the molecular force sensor. Individual trajectories can be plotted in the x-y plane, visualizing their diffusion behavior and localization (left panel). Furthermore, each trajectory's intensities can be plotted over time to identify FRET transitions or bleaching steps (middle panel shows the red trajectory from the left panel). The resulting FRET efficiency and stoichiometry can be visualized similarly (right panel). (C) smFRET trajectories using the Holliday junction DNA sensor. HBSS + 12 mM MgCl2 was used as a buffer during the measurements. Apart from the apparent acceptor bleaching step near the sequence end of these examples, the frequency of FRET transitions for each sensor can be determined. The Holliday junctions switch their conformation with a high frequency, whereas the molecular force sensor does not exhibit FRET transitions. This information makes it possible to adjust the experimental conditions, such as the delay between the frames, to increase or reduce the number of observed transitions. Abbreviations: FRET = Förster resonance energy transfer; smFRET = single-molecule FRET; HBSS = Hank's balanced salt solution. Please click here to view a larger version of this figure.

Supplemental Information: Localization and tracking of single molecules (01. Tracking.ipynb). Please click here to download this File.

Discussion

This article details a pipeline for the automated recordings and quantitative analysis of smFRET data originating from mobile yet surface-tethered probe molecules. It complements the two predominant approaches to smFRET experiments, involving either surface-immobilized probes or probes diffusing in solution into and out of a confocal excitation volume17. It provides the correct FRET efficiency and the molecular positions as a function of time. It can therefore be used as input for specialized analysis programs, e.g., to quantify transition kinetics1, FRET histograms39, or two-dimensional diffusion22.

The software is released under a free and open-source license approved by the Open Source Initiative that grants the user the perpetual right to free usage, modification, and redistribution. Github was chosen as a development and distribution platform to make it as easy as possible to obtain the software and participate in the development process by reporting bugs or contributing code40. Written in Python, the software does not depend on proprietary components. The choice of Jupyter notebooks as user interfaces facilitates the inspection of data at every analysis step and allows for tailoring and extending the pipeline specifically for the experimental system at hand. The sdt-python library32 serves as the foundation and implements functionality to evaluate fluorescence microscopy data, such as single-molecule localization, diffusion analysis, fluorescence intensity analysis, color channel registration, colocalization analysis, and ROI handling.

In principle, single-particle tracking can be performed in one-, two- or three-dimensional systems. Here, the single-molecule analysis pipeline was tailored to the study of 2D mobile systems. This choice mirrors the availability of simple systems, such as planar-supported lipid bilayers (SLBs), to present mobile fluorescent probes. Such lipid bilayer systems are typically composed of two or more phospholipids moieties, where the bulk fraction determines the key physicochemical parameters of the SLB (such as phase and viscosity), and the minor fraction provides attachment sites for biomolecules. These attachment sites can be biotinylated phospholipids for avidin- or streptavidin-based protein platforms or nickel-NTA conjugated phospholipids for protein platforms with histidine tags41. The choice of the appropriate platform for linking proteins to the SLB depends on the scientific question. Readers can refer to the literature16,38,42 for examples of successfully employed strategies. The density of probes in the sample should be sufficiently low to avoid overlapping point spread functions; typically, less than 0.1 molecules per µm2 are recommended. See the representative results section (in particular, Figure 6) for an example showing a suitable probe density. The analysis method is also applicable to single fluorescently labeled protein molecules diffusing in the plasma membrane of live cells.

One critical aspect of smFRET experiments is the production and characterization of the FRET probes themselves. When choosing fluorophores for a FRET pair, their Förster radius should match the expected inter-dye distances43. Dyes resistant to photobleaching are preferred as they yield long time traces. However, for elevated bleaching rates, one fluorophore species can be utilized to recognize multiemitter events originating from colocalized molecules via stepwise photobleaching analysis; see step 8.1.4 in the protocol section. Fluorophore pairs should be site-specifically and covalently attached to the molecules of interest, forming intra- or intermolecular FRET pairs.

Combining smFRET with other readily available techniques can increase its spatial resolution beyond the diffraction limit (via STED44). The smFRET tracking algorithm presented here widens the approach's applicability to new experimental settings and model systems. This includes studies of (i) kinetic changes in the stoichiometry of mobile biomolecules, (ii) dynamic association of mobile biomolecules, (iii) the rate of enzymatic reactions of freely diffusing reactants, and (iv) the kinetics of conformational changes of mobile biomolecules. The first two examples require model systems that show intermolecular FRET, i.e., donor and acceptor are conjugated to separate biomolecular entities of interest. The latter examples may make use of biosensors carrying donor and acceptor within the same molecular entity (intramolecular FRET).

Intramolecular FRET-based sensors can provide insight into intrinsic conformational changes of biomolecules1,2,3,4, conformational changes caused by endogenous or external force load (molecular force sensors16), or ion concentrations in the nano-environment such as calcium45 and pH46. Depending on the model system and the preferred anchoring platform, such smFRET events can either be tracked in 2D or 3D: (i) planar tracking of smFRET events can be employed for the quantification of receptor-ligand interaction times within a plasma membrane, the association of membrane-anchored signal amplification cascades, and the stoichiometry changes of surface receptors; (ii) volume tracking of smFRET events can be used for any intra- or intermolecular FRET probes in living cells or in in vitro reconstituted systems.

The smFRET tracking method was developed mainly with intramolecular FRET probes in mind. These probes feature a fixed and well-known number of fluorescent labels, a fact that was exploited to reject data from agglomerated and incorrectly synthesized (e.g., incompletely labeled) molecules, as well as from probes where one of the fluorophores has been photobleached. However, by adjusting the filtering steps, the method can also be applied to intermolecular FRET probes. For instance, instead of accepting only molecules featuring a single donor and a single acceptor fluorophore, one could examine the spatial trajectories of donor and acceptor dyes and select, for example, for co-diffusing donor-acceptor trajectories.

As the 3D-DAOSTORM algorithm has support for determining a signal's position along the optical axis via the astigmatism due to a cylindrical lens in the emission beam path, 3D experiments could be easily integrated into the analysis pipeline. In this case, the acceptor signal upon acceptor excitation would serve to determine the stoichiometry and the axial position. The analysis software can also be employed to evaluate data from experiments featuring immobilized probes by utilizing its large degree of automation and filtering schemes. In fact, smFRET efficiency datasets from Holliday junctions immobilized on gel-phase bilayers38 were analyzed using an early version of the software.

Disclosures

The authors have nothing to disclose.

Acknowledgements

This work was supported by the Austrian Science Fund (FWF) projects P30214-N36, P32307-B, and by the Vienna Science and Technology Fund (WWTF) LS13-030.

Materials

1,2-dioleoyl-sn-glycero-3-[(N-(5-amino-1-carboxypentyl)iminodiacetic acid)succinyl] (nickel salt) (Ni-NTA-DOGS) Avanti Polar Lipids 790404P
1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) Avanti Polar Lipids 850375P
1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) Avanti Polar Lipids 850457P
α Plan-FLUAR 100x/1.45 oil objective Zeiss 000000-1084-514
Axio Observer microscope body Zeiss
Bandpass filter Chroma Technology Corp ET570/60m donor emission filter
Bandpass filter Chroma Technology Corp ET675/50m acceptor emission filter
conda-forge conda-forge community community-maintaned Python package repository for Anaconda/miniconda
Coverslips 60 mm x 24 mm #1.5 MENZEL
Dichroic mirror Semrock Inc FF640-FDi01-25×36 separation of donor and acceptor emission
Dichroic mirror (quad band) Semrock Inc Di01-R405/488/532/635-25×36 separation of excitation and emission light
DPBS Sigma-Aldrich D8537
FCS Sigma-Aldrich F7524 for imaging buffer
fret-analysis Schütz group at TU Wien Python package for smFRET data analysis; version 3
Fura-2 AM Thermo Fisher Scientific 11524766
HBSS Sigma-Aldrich H8264 for imaging buffer
iBeam Smart 405-S 405 nm laser Toptica Photonics AG
iXon Ultra 897 EMCCD camera Andor Technology Ltd
Lab-Tek chambers (8 wells) Thermo Fisher Scientific 177402PK for sample preparation and imaging
Millenia Prime 532 nm laser Spectra Physics
miniconda Anaconda Inc. Python 3 distribution. Min. version: 3.7
Monovalent streptavidin (plasmids for bacterial expression) Addgene 20860 & 20859
OBIS 640 nm laser Coherent Inc 1185055
Optosplit II Cairn Research
Ovalbumin Sigma-Aldrich A5253 for imaging buffer
Plasma cleaner Harrick Plasma PDC-002
sdt-python Schütz group at TU Wien Python library for data analysis; version 17
TetraSpek bead size kit Thermo Fisher Scientific T14792 Randomly distributed, immobilized fiducial markers for image registration
USC500TH Ultrasound bath VWR for SUV formation

References

  1. McKinney, S. A., Déclais, A. -. C., Lilley, D. M. J., Ha, T. Structural dynamics of individual holliday junctions. Nature Structural Biology. 10 (2), 93-97 (2002).
  2. Wang, S., Vafabakhsh, R., Borschel, W. F., Ha, T., Nichols, C. G. Structural dynamics of potassium-channel gating revealed by single-molecule FRET. Nature Structural & Molecular Biology. 23 (1), 31-36 (2015).
  3. Hellenkamp, B., Wortmann, P., Kandzia, F., Zacharias, M., Hugel, T. Multidomain structure and correlated dynamics determined by self-consistent FRET networks. Nature Methods. 14 (2), 174-180 (2016).
  4. Kilic, S., et al. Single-molecule FRET reveals multiscale chromatin dynamics modulated by HP1α. Nature Communications. 9 (1), 235 (2018).
  5. Stryer, L. Fluorescence energy transfer as a spectroscopic ruler. Annual Review of Biochemistry. 47 (1), 819-846 (1978).
  6. Wu, P. G., Brand, L. Resonance energy transfer: Methods and applications. Analytical Biochemistry. 218 (1), 1-13 (1994).
  7. Qiao, Y., Luo, Y., Long, N., Xing, Y., Tu, J. Single-molecular förster resonance energy transfer measurement on structures and interactions of biomolecules. Micromachines. 12 (5), 492 (2021).
  8. Malkusch, N., Dörfler, T., Nagy, J., Eilert, T., Michaelis, J. smFRET experiments of the RNA polymerase II transcription initiation complex. Methods. 120, 115-124 (2017).
  9. Lee, J. -. B., et al. Single-molecule views of MutS on mismatched DNA. DNA repair. 20, 82-93 (2014).
  10. Phelps, C., Israels, B., Jose, D., Marsh, M. C., von Hippel, P. H., Marcus, A. H. Using microsecond single-molecule FRET to determine the assembly pathways of T4 ssDNA binding protein onto model DNA replication forks. Proceedings of the National Academy of Sciences of the United States of America. 114 (18), E3612-E3621 (2017).
  11. Deindl, S., Zhuang, X. Monitoring conformational dynamics with single-molecule fluorescence energy transfer: Applications in nucleosome remodeling. Methods in Enzymology. 513, 59-86 (2012).
  12. Crawford, D. J., Hoskins, A. A., Friedman, L. J., Gelles, J., Moore, M. J. Single-molecule colocalization FRET evidence that spliceosome activation precedes stable approach of 5′ splice site and branch site. Proceedings of the National Academy of Sciences of the United States of America. 110 (17), 6783-6788 (2013).
  13. Wang, Y., Xiao, M., Li, Y. Heterogeneity of single molecule FRET signals reveals multiple active ribosome subpopulations. Proteins. 82 (1), 1-9 (2014).
  14. Mori, T., Vale, R. D., Tomishige, M. How kinesin waits between steps. Nature. 450 (7170), 750-754 (2007).
  15. Huppa, J. B., et al. TCR-peptide-MHC interactions in situ show accelerated kinetics and increased affinity. Nature. 463 (7283), 963-967 (2010).
  16. Göhring, J., et al. Temporal analysis of T-cell receptor-imposed forces via quantitative single molecule FRET measurements. Nature Communications. 12 (1), 2502 (2021).
  17. Hellenkamp, B., et al. Precision and accuracy of single-molecule FRET measurements-a multi-laboratory benchmark study. Nature Methods. 15 (9), 669 (2018).
  18. Kapanidis, A. N., et al. Fluorescence-aided molecule sorting: Analysis of structure and interactions by alternating-laser excitation of single molecules. Proceedings of the National Academy of Sciences of the United States of America. 101 (24), 8936-8941 (2004).
  19. Sakon, J. J., Weninger, K. R. Detecting the conformation of individual proteins in live cells. Nature Methods. 7 (3), 203-205 (2010).
  20. McCann, J. J., Choi, U. B., Zheng, L., Weninger, K., Bowen, M. E. Optimizing methods to recover absolute FRET efficiency from immobilized single molecules. Biophysical Journal. 99 (3), 961-970 (2010).
  21. Lee, N. K., et al. Accurate FRET measurements within single diffusing biomolecules using alternating-laser excitation. Biophysical Journal. 88 (4), 2939-2953 (2005).
  22. Asher, W. B., et al. Single-molecule FRET imaging of GPCR dimers in living cells. Nature Methods. 18 (4), 397-405 (2021).
  23. Joo, C., Ha, T. . Single-molecule FRET with total internal reflection microscopy. (12), 1223-1237 (2012).
  24. Joo, C., Ha, T. Objective-type total internal reflection microscopy (excitation) for single-molecule FRET. Cold Spring Harbor Protocols. 11, 1189-1191 (2012).
  25. Joo, C., Ha, T. Objective-type total internal reflection microscopy (emission) for single-molecule FRET. Cold Spring Harbor Protocols. 11, 1192-1194 (2012).
  26. Allan, D. B., Caswell, T., van der Wel, C. M., Dimiduk, T. . Soft-matter/pims: PIMS v0.5. , (2020).
  27. Anaconda Inc. . Miniconda. , (2021).
  28. conda-forge community. . The conda-forge project: community-based software distribution built on the conda package format and ecosystem. , (2015).
  29. . . JupyterLab Contributors Notebooks – JupyterLab documentation. , (2021).
  30. Babcock, H., Sigal, Y. M., Zhuang, X. A high-density 3D localization algorithm for stochastic optical reconstruction microscopy. Optical Nanoscopy. 1 (6), (2012).
  31. Gao, Y., Kilfoil, M. L. Accurate detection and complete tracking of large populations of features in three dimensions. Optics Express. 17 (6), 4685 (2009).
  32. Schrangl, L. . sdt-python: Python library for fluorescence microscopy data analysis (v17.1). , (2021).
  33. Crocker, J. C., Grier, D. G. Methods of digital video microscopy for colloidal studies. Journal of Colloid and Interface Science. 179 (1), 298-310 (1996).
  34. Preibisch, S., Saalfeld, S., Schindelin, J., Tomancak, P. Software for bead-based registration of selective plane illumination microscopy data. Nature Methods. 7 (6), 418-419 (2010).
  35. Bradski, G. The OpenCV library. Dr. Dobb’s Journal: Software Tools for the Professional Programmer. 25 (11), 120-123 (2000).
  36. Allan, D. B., Caswell, T., Keim, N. C., van der Wel, C. M., Verweij, R. W. Soft-matter/trackpy: Trackpy v0.5.0. Zenodo. , 4682814 (2021).
  37. Killick, R., Fearnhead, P., Eckley, I. A. Optimal detection of changepoints with a linear computational cost. Journal of the American Statistical Association. 107 (500), 1590-1598 (2012).
  38. Schrangl, L., Göhring, J., Schütz, G. J. Kinetic analysis of single molecule FRET transitions without trajectories. The Journal of Chemical Physics. 148 (12), 123328 (2018).
  39. Santoso, Y., Torella, J. P., Kapanidis, A. N. Characterizing single-molecule FRET dynamics with probability distribution analysis. ChemPhysChem. 11 (10), 2209-2219 (2010).
  40. Schrangl, L. Single-molecule FRET analysis software (3.0). Zenodo. , (2021).
  41. Nye, J. A., Groves, J. T. Kinetic control of histidine-tagged protein surface density on supported lipid bilayers. Langmuir. 24 (8), 4145-4149 (2008).
  42. Platzer, R., et al. Unscrambling fluorophore blinking for comprehensive cluster detection via photoactivated localization microscopy. Nature Communications. 11 (1), 4993 (2020).
  43. Johnson, I., Spence, M. . The molecular probes handbook: A guide to fluorescent probes and labeling technologies. , (2010).
  44. Szalai, A. M., et al. Super-resolution imaging of energy transfer by intensity-based STED-FRET. Nano Letters. 21 (5), 2296-2303 (2021).
  45. Miyawaki, A., et al. Fluorescent indicators for Ca2+ based on green fluorescent proteins and calmodulin. Nature. 388 (6645), 882-887 (1997).
  46. Zhai, B., Zhai, S., Hao, R., Xu, J., Liu, Z. A FRET-based two-photon probe for in vivo tracking of pH during a traumatic brain injury process. New Journal of Chemistry. 43 (43), 17018-17022 (2019).

Play Video

Cite This Article
Schrangl, L., Göhring, J., Kellner, F., Huppa, J. B., Schütz, G. J. Automated Two-dimensional Spatiotemporal Analysis of Mobile Single-molecule FRET Probes. J. Vis. Exp. (177), e63124, doi:10.3791/63124 (2021).

View Video