Workflow

In this section you will learn how to correct single molecule intensity-time traces from experimental bias and obtain state trajectories.

Each time molecule data are processed by pressing UPDATE, UPDATE ALL, or by navigating through molecules, data in modules Histogram analysis and Transition analysis are automatically refreshed.

The procedure includes nine steps:

  1. Import trajectories
  2. Adjust single molecule position
  3. Correct intensities
  4. Sort trajectories into subgroups
  5. Correct for photobleaching and blinking
  6. Smooth trajectories
  7. Correct ratio values
  8. Determine state trajectories
  9. Process all data and merge projects

Import trajectories

Initial trajectories come from Simulation’s output in a simulation-based project, Video processing’s output in a video-based project, or from external ASCII files.

In the two first cases, you can skip this step and go directly to the next one.

In the latter case, a new trajectory-based project must be created. This implies to import the trajectory files, register the associated experiment settings and define the data structure in the files. After the project creation is completed, it is recommended to save it to a .mash file that should regularly be overwritten in order to keep traceability and access to the results.

Informing MASH about the particular experiment settings is crucial to adapt the software functionalities to your own experiment setup. In theory, the software is compatible with:

  • an unlimited number of video channels,
  • an unlimited number of alternated lasers,
  • FRET calculations for an unlimited number of FRET pairs.

To create a new trajectory-based project:

  1. Open the experiment settings window by pressing New project in the project management area and selecting import trajectories.

  2. Import a set of trajectory files and define your experiment setup by configuring tabs:

    Import
    Channels
    Lasers
    Calculations
    Divers

    If necessary, modify settings in Calculations and Divers any time after project creation.

  3. Define how data are structured in the files by configuring tab File structure.

  4. Finalize the creation of your project by pressing Save; the experiment settings window now closes and the interface switches to module Trace processing.

  5. Save modifications to a .mash file by pressing Save project in the project management area.


Adjust single molecule position

Single molecule image in each channel and upon illumination by a specific laser are shown in the Visualization area. Single molecule positions are marked with straight crosses and the dimensions of the pixel area used to integrate intensities are delimited by red solid lines. To obtain the most accurate intensities, single molecule positions must be centered on the brightest pixel of the integration zone and the area of two neighbouring molecules must not overlap.

Adjusting molecule positions

Because of imperfect coordinates transformation, it can happen that positions are shifted by one or two pixels away from the brightest pixel. In that case, the positions must be recentered on the brightest pixel.

MASH offers the possibility to automatically recenter single molecule positions; see Remarks for more details.

To recenter single molecule positions:

  1. Select the molecule index in the Molecule list.

  2. If necessary, adjust the brightness and contrast in Single molecule images to render the molecule profile the most apparent.

  3. Press recenter in Single molecule images to adjust molecule positions in all detection channels; intensity-time traces will automatically be recalculated after re-centering.

Note: If no video or molecule coordinates are available, single molecule images can not be shown and molecule positions can not be adjusted.


Correct intensities

Recorded intensities include the contribution of signals that must be subtracted to calculate reliable FRET and stoichiometry values. These unwanted signals are the background intensity and the cross-talks.

Background correction

The background intensity is channel- and laser-specific. It consists mainly of the dark count of the detector and background light, like chamber auto-fluorescence for instance. The background signal is usually distributed over the single molecule image and is therefore more accurately estimated in the local environments of single molecules. MASH offers a set of local background estimators that can be used for such purpose.

The intensity Iemex(n,t) of molecule n detected in emission channel em upon illumination with laser ex at time t is background-corrected such as:

\left I_{em}^{ex}(n,t)^{*} = I_{em}^{ex}(n,t) - bg_{em}^{ex}(n,t)

with Iemex(n,t)* the background-corrected intensity and bgemex(n,t) the estimated background intensity.

To correct intensities from background:

  1. If not already done, select the molecule index in the Molecule list.

  2. For each intensity-time trace, set parameters:

    Background correction settings
    Apply background correction

  3. Update data correction and display by pressing UPDATE in the Control area

  4. If desired, apply the same parameter settings to all molecules by pressing all in panel Background correction


Trajectory re-sampling

When the signal-to-noise ratio of the intensity trajectories is too low, the overlapping between FRET states becomes to large to be resolved and the determination of states and state transition kinetics becomes nearly impossible. One way to remedy to such issue is to re-sample the intensity trajectories using a sampling time larger than the original one.

Trajectory re-sampling

Larger sampling times immediately produce higher SNRs but worsen the time resolution. It is therefore to the user to determine the best compromise between SNR and the resolution of state transition kinetics to improve the separation between the states.

To re-sample trajectories:

  1. Set the Trajectory sampling time

  2. Update intensity calculation and display by pressing UPDATE in the Control area

Note: The trajectory sampling time is used to re-sample all trajectories in the molecule set.


Cross-talk correction

Cross-talks are due to instrumental imperfections and include two phenomena:

  • the detection of an emitter fluorescence into unspecific video channels, called the bleedthrough
  • the detection of an emitter fluorescence into specific video channel after unspecific laser illumination, called the direct excitation

Instrumental imperfections

The bleedthrough and direct excitation coefficient can be determined from control experiments involving single-labelled species. The bleedthrough coefficient btem,em0 of an emitter em into detection channel of an emitter em0 is calculated with intensities measured from species carrying only emitter em and after background correction, such as:

bt_{em,em0} = {}\frac{I_{em0}^{ex}^{*}}{I_{em}^{ex}^{*}}

The direct excitation coefficient dEem0exex0 of an emitter em0 upon unspecific illumination ex is calculated with intensities measured from species carrying only emitter em0 and after background correction, such as:

dE_{em0}^{ex\neq ex0} = {}\frac{I_{em0}^{ex}^{*}}{I_{em0}^{ex0}^{*}}

with ex0 the specific illumination of emitter em0.

The background-corrected intensity Iem0ex(n,t)* of molecule n detected in emission channel em0 upon illumination with laser ex at time t is corrected from cross-talks such as:

I_{em0}^{ex}(n,t)^{**} = I_{em0}^{ex}(n,t)^{*} - \sum_{em\not\equiv em0}\left ( bt_{em,em0}(n) × I_{em}^{ex}(n,t)^{*} \right )
I_{em0}^{ex\neq ex0}(n,t)^{***} = I_{em0}^{ex\neq ex0}(n,t)^{**} - dE_{em0}^{ex\neq ex0}(n) × I_{em0}^{ex0}(n,t)^{**}

with Iem0ex(n,t)** and Iem0ex(n,t)***, the intensities corrected from bleedthrough-only and from both cross-talks respectively.

To correct intensities from cross-talks:

  1. For each emitter in the Emitter list, set parameters:

    Bleedthrough coefficients
    Direct excitation coefficients

  2. Update data correction and display by pressing UPDATE in the Control area

Note: The set of cross-talk coefficients is used to correct all trajectories in the molecule set.


Sort trajectories into subgroups

Single molecules trajectories can be manually or automatically sorted into subgroups to perform category-specific analysis. For instance, molecules with correct dye-labelling and showing dynamics can be saved into a separate project to determine the state configuration and transition rates involved in molecule dynamics.

It can also be convenient to sort molecules with incomplete dye labelling in order to make statistics on the labelling efficiency in the sample.

Beside, it is important to identify and exclude incoherent intensity-time traces from the analysis set.

Molecule sorting

MASH offers a tool called Trace manager that can be used to perform such task. Otherwise, molecules can be sorted manually in the Trace processing interface; see Remarks for more details.

To sort molecules and save a particular subgroup:

  1. Open Trace manager by pressing TM.

  2. Sort single molecule data into subgroups by referring to
    Use Trace manager; it is recommended to save the sorting by pressing Save and overwriting the project file.

  3. Clear out unselected molecules form the Molecule list and from the project by pressing Clear.

  4. Save selected molecules to a new .mash file by pressing Save.


Correct for photobleaching and blinking

It can happen that emitters get photochemically destroyed after absorbing a certain amount of photons. This phenomenon is called photobleaching and results in the permanent loss of signal in the dye-specific emission channel, which in turn translates into a drop of the corresponding intensity-time trace to zero.

Effect of photobleaching of intensity ratio

On the other hand, reversible changes in the structure or environment of the emitters can induce temporary interruption of dye emission, which results in short-lived drops of intensity-time traces to zero.

These zero-intensity portions of the time trace may bias the following histogram and transition analysis by creating irrelevantly fluctuating FRET data. To prevent such bias, it is necessary to detect and ignore intensity interruptions by truncating the time traces.

Photobleaching is corrected by simply truncating the end of the trajectory.

For temporary intensity interruptions, simply excising zero-intensity portions from the trace would cause the continuity of the time sequence to be disrupted, and thus, would falsify state lifetimes and state transition rates. Therefore, intensity interruptions are corrected by dividing the trajectory in two separate traces at the time when the intensity interruption ends, and by subsequently truncating the end of the left-trajectory by applying photobleaching correction.

To truncate temporary intensity interruptions:

  1. Set parameters

    Photobleaching detection method (to Manual)
    Photobleaching cutoff (to the ending of intensity interruption)

  2. Split traces in two by pressing Split.

  3. Correct the original molecule for photobleaching (as described below)

To detect and truncate photobleached data:

  1. If not already done, select the molecule index in the Molecule list

  2. Set parameters

    Photobleaching detection method
    Method parameters
    Truncate trajectories

  3. Update data correction and display by pressing UPDATE in the Control area

  4. If desired, apply the same parameter settings to all molecules by pressing all in panel
    Photobleaching


Smooth trajectories

Intensity-time traces are naturally noisy due to the stochastic nature of photon emission but also of camera detection, signal amplification and signal conversion. The noise magnitude in a background- and cross-talk-corrected intensity-time trace Iemex(n,t)*** is characterized by the signal-to-noise ratio SNRemex(n) calculated as:

SNR_{em}^{ex}(n) = \frac{\mu_{em}^{ex}(n)^{***}}{\sigma_{em}^{ex}(n)^{***}}

with μemex(n)*** and σemex(n)***, the respective mean and standard deviation of the intensities over the observation time.

The lower the intensity, the lower the SNR and the noisier the intensity-time traces.

A low SNR creates a large broadening in intensity histograms and even larger in intensity ratio histograms, e. g. in FRET histograms. Large broadening in FRET histograms implies greater overlaps between FRET populations which make the identification of the state configuration more difficult.

One way to minimize the population overlap in histograms is to artificially reduce the noise magnitude, or to “smooth” intensity-time traces.

To smooth trajectories:

  1. If not already done, select the molecule index in the Molecule list.

  2. Set parameters

    Denoising settings
    Apply denoising

  3. Update correction and display by pressing UPDATE in the Control area

  4. If desired, apply the same parameter settings to all molecules by pressing all in panel Denoising

Note: Smoothing is not recommended when determining state trajectories as it induces a modification in data distribution, while most state finding algorithms rely on a naturally distributed noise to identify states and detect state transitions.


Correct ratio values

To recover the labelling efficiencies and absolute distances between the donor and the acceptor of a FRET pair, it is necessary, but not sufficient, to put donor and acceptor intensities on the same scale. For this, the differences in detection efficiencies, quantum yields, extinction coefficients and excitation intensities between donor D and acceptor A must be corrected.

The differences in detection efficiencies and quantum yields are accounted for in the FRET pair-specific γD,A factor such as:

\gamma _{D,A} = \frac{\eta_{A}\phi_{A}}{\eta_{D}\phi_{D}}

with ΦA and ΦD the acceptor and donor quantum yield respectively, and ηA and ηD the photon detection efficiency in the respective emission channels.

The differences in extinction coefficients and excitation intensities are accounted for in the FRET pair-specific βD,A factor such as:

\beta _{D,A} = \beta_{D,A} = \frac{I_{A,ex}\sigma_{A,ex}^{A}}{I_{D,ex}\sigma_{D,ex}^{D}}

with IA,ex and ID,ex the intensities of acceptor- and donor-specific excitations respectively, and σA,exA and σD,exD the respective absorption cross-sections at specific excitation.

Apparent FRET E*D,A(n,t) of a donor-acceptor pair D-A is calculated according to FRET calculations and is γ-corrected into ED,A(n,t) such as:

E_{D,A}(n,t) = \left \{ 1+\frac{\gamma_{D,A}(n) I_{D,em}^{D,ex}(n,t)^{***}}{ I_{A,em}^{D,ex}(n,t)^{***} + \sum_{i>A}\left [\frac{E_{A,i}}{1-E_{A,i}}I_{A,em}^{D,ex}(n,t)^{***} \right ] -\sum_{D<i<A}\left[ \frac{E_{i,A}}{1-E_{i,A}}\gamma_{i,A}(n)I_{i,em}^{D,ex}(n,t)^{***}\right ]}  \right \}^{-1}

Apparent stoichiometry S*D,A(n,t) of a donor-acceptor pair D-A is calculated according to Stoichiometry calculations and is γ- and β-corrected into SD,A(n,t) such as:

S_{D,A}(n,t) = \left |1+\frac{ I_{A,em}^{A,ex}(n,t)^{***} + \sum_{i>A}\left [ I_{i,em}^{A,ex}(n,t)^{***} \Big/ \gamma_{A,i}(n) \right] }{\beta_{D,A} \left \{ I_{A,em}^{D,ex}(n,t)^{***} + \sum_{D\leq i<A}\left [\gamma_{i,A}(n) I_{i,em}^{D,ex}(n,t)^{***} \right ] + \sum_{i>A}\left [I_{i,em}^{D,ex}(n,t)^{***} \Big/ \gamma_{A,i}(n) \right ] \right \} }  \right |^{-1}

To correct apparent FRET- and stoichiometry-time traces with γ and β factors:

  1. If not already done, select the molecule index in the Molecule list.

  2. For each donor-acceptor FRET pair in the FRET pair list, set parameters:

    Factor estimation method
    Gamma factor
    Beta factor

  3. Update data correction and display by pressing UPDATE in the Control area

  4. If desired, apply the same parameter settings to all molecules by pressing all


Determine state trajectories

The main goal of module Trace processing is to obtain a reliable set of single molecule intensity- and FRET-time traces in order to infer reliable state trajectories.

State trajectories are idealized time traces, i.e., without noise. In other words, state trajectories are a succession of discrete values, and are here called “discretized” time traces.

Discretized time traces allow to characterize the conformational states visited by the molecule during the observation time and to quantify the duration the molecule dwells in each state.

FRET state configuration in trajectory

Further analysis with the module Transition analysis yields the most probable state configuration of single molecules and the rate coefficients that govern transitions between those states.

To infer state trajectories:

  1. If not already done, select the molecule index in the Molecule list.

  2. Set parameters

    Discretization method
    Data to discretize
    Method parameters
    Post-processing parameters

  3. Update calculations and display by pressing UPDATE in the Control area

  4. If desired, apply the same parameter settings to all molecules by pressing all


Process all data and merge projects

Instead of processing one molecule at a time, a general configuration can be applied to all molecules by using the all buttons and all data can be processed at once using the UPDATE ALL button.

Finally, to constitute a more representative sample of molecules to be analyzed in
Transition analysis and Histogram analysis, molecules from different projects with identical experiment settings can be merged into one big data set.

To process all molecules at once:

  1. Update all data corrections and calculations by pressing UPDATE ALL in the Control area

To merge projects:

  1. Select the multiple projects to merge in the Project list

  2. Right-click in the Project list and choose the option Merge projects

  3. Press Yes to start the merging process

    Note: The merging process induces a loss of single molecule videos that were used in individual projects. Therefore, it is recommended to perform all adjustments of molecule positions and background corrections prior merging; see Remarks for more information.


Remarks

If the initial single molecule position is shifted more than 3 pixels from the brightest pixel, spots coordinates are considered as ill-defined and the automatic recentering option will not function. In this case, it is recommended to review the spot detection and/or transformation procedure to obtain a more decent set of single molecule coordinates; see Spotfinder and Coordinates transformation for more information. However, if desired, molecule x- and y-coordinates must be modified manually; see Single molecule coordinates for more information.

Molecules can be given several tags and/or deselected from the set while browsing the Molecule list in the Trace processing interface; see Molecule status for more information.

To merge multiple projects into one, experiment setups must be the same, i.e., the frame rate, the number of channels and laser excitations, the laser wavelengths, the emitter labels, the FRET pairs and the stoichiometry calculations must be identical.