... | ... | @@ -40,21 +40,34 @@ In this WiKi we provide examples on how to analyze and process FCS and imaging d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
## Generate bleach corrected correlation functions using FluctuationAnalyzer 4G
|
|
|
|
|
|
The [Fluctuation Analyzer 4G]( https://www.embl.de/~wachsmut/downloads.html) (FA) is a software tool for the interactive as well as automated processing of fluorescence auto- and cross-correlation spectroscopy (FCS/FCCS) data. It can read raw data, i.e., one-or two-channel photon streams, from various commercial suppliers of FCS/FCCS data acquisition equipment, organize such data in processing sessions by file-based management, calculate temporal auto- and cross-correlation functions, correct for photobleaching, cross-talk, and background signal and fit the data with appropriate model functions before saving the results. Refer to the manual of FA and the original article [Wachsmuth et al. (2015)]( http://europepmc.org/abstract/MED/25774713) for details.
|
|
|
|
|
|
|
|
|
In this WiKi we provide a detailed explanation on which settings and parameters to use for the FCS calibration protocol and process the data for the fluorescent dye and fluorescent protein. At least the following 4 steps must be executed
|
|
|
<a name=back></a>
|
|
|
1. [Loading data into FA](#faload)
|
|
|
|
|
|
2. [Compute correlation curves with FA](#facorr)
|
|
|
|
|
|
3. [Compute parameters for correction](#facorrpar)
|
|
|
3. [Compute offset and crosstalk parameters](#facorrpar)
|
|
|
|
|
|
4. [Compute correction factors and save table](#facorrfa)
|
|
|
|
|
|
|
|
|
> <p style='color:green'> **Apply to**: Throughout the different steps using fluctuation analyzer it is always required to click on `Apply to` for changes in settings to be active. </p>
|
|
|
|
|
|
4. [Compute correction factors](#facorrfa)
|
|
|
Typical values used for the computation of the correlation curves and correction factors are:
|
|
|
|
|
|
Parameters | Values |
|
|
|
:--- | :---
|
|
|
Base freq. (Hz) | 1.000.000 (dye) or 100.000 (protein)
|
|
|
Composition | Ch1 <> Ch2
|
|
|
Offset (kHz) | 1-5 kHz <br/> Must be measured
|
|
|
Cross-talk (Ch1-Ch2), when using two color FCS | 0.04-0.05 (Ch1: mEGFP, Ch2: mCherry). <br/> Value depends on protein and optical settings. Must be measured.
|
|
|
|
|
|
5. [Save FA result table ](#fares)
|
|
|
|
|
|
### [<img src='./images/up.png'>](#back) <a name=faload></a>Loading data into FA
|
|
|
|
... | ... | @@ -75,26 +88,33 @@ In this WiKi we provide a detailed explanation on which settings and parameters |
|
|
### [<img src='./images/up.png'>](#back) 2. <a name=facorr></a>Compute correlation curves with FA
|
|
|
FA computes auto- and cross-correlation functions given the raw data. FA corrects for slow time varying trends in the photon-counts, such as photobleaching, by calculating the correlation function in small time-windows where the slow trend can be neglected. Then all correlation functions are averaged.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
|
|
|
:--: | ---
|
|
|
<img src = './images/FA/tab_correlate.png' width = "400px" > | 1. Switch to the page `Modify and correlate`
|
|
|
<img src = './images/FA/correlate_idx.png' width = "400px" > | 2. Choose the correlation curves to compute. <br/> 3. Set the base frequency. For the fluorescent dye use 1.000.000 Hz (i.e. 1 us time interval). For a fluorescent protein 100.000 Hz is enough (i.e. 10 us time interval) <br/> 4. Set `Ch1 <> Ch2` for indipendently correct the channels <br/> 5. Check if `Autosave` is on <br/> 6. Click `Apply to` <br/> 7. Click on `Calculate all`
|
|
|
<img src = './images/FA/correlate_idx.png' width = "400px" > | 2. Choose the correlation curves to compute. <br/> 3. Set the base frequency. For the fluorescent dye use 1.000.000 Hz (i.e. 1 us time interval). For a fluorescent protein 100.000 Hz is enough (i.e. 10 us time interval) <br/> 4. Set `Ch1 <> Ch2` for indipendently correct the channels <br/> 5. Check if `Autosave` is on <br/> 6. Click `Apply to` all <br/> 7. Click on `Calculate all`
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
### [<img src='./images/up.png'>](#back) 3. <a name=facorrpar> </a>Compute offset and crosstalk parameters
|
|
|
|
|
|
###[<img src='./images/up.png'>](#back) 3. <a name=facorrpar> </a>Compute parameters for correction
|
|
|
The [**Offset**](#offset) and [**Crosstalk**](#xtalk) parameters are used for computing the [correction factors](#facorrfa) in the next step. The crosstalk value is only needed in case of two color FCS (FCCS).
|
|
|
|
|
|
|
|
|
* <a name='offset'>**Offset**</a>: Average photon-count rate obtained from cells that do not express the fluorescent protein (WT cells). Alternatively the medium can also be measured, however this typically underestimate the background. The laser power and light-path settings must be the same as for the measurement of the fluorescent protein.
|
|
|
* Acquire data from WT cells (5-10 cells)
|
|
|
* Load data into FA and compute intensity corrections
|
|
|
* Save data into result table (**.res* file)
|
|
|
|
|
|
* Save data into a tab delimited result table (**.res* file)
|
|
|
* Load *res* file into a spreadseet and compute average of the column named *Interval Ch1* (for `Offset [Khz] Ch1`) and/or *Interval Ch2* (for `Offset [Khz] Ch2` )
|
|
|
|
|
|
|
|
|
| |
|
|
|
:---: | ----
|
|
|
<img src = './images/FA/WT_correction_idx.png' width = "400px" > | [Import](#faload) FCS measurements from WT cells into FA. The [Modify and correlate](#facorr) step is not required <br/> 1. Change tab to `Intensity corrections` <br/> 2. Check that `Autosave` is on <br/> 3. Click `Apply to` <br/> 4. Click on `Calculate all`
|
|
|
<img src = './images/FA/save_res.png' width = "400px" > | 5. Change tab to `Save, export and report` <br/> 6. Click on `Save all` to save single taces <br/>7. Click on `FA format` to save a summary result table (here 2c.res)
|
|
|
:---: | ----
|
|
|
<img src = './images/FA/WT_correction_idx.png' width = "400px" > | [Import](#faload) FCS measurements from WT cells into FA. The [Modify and correlate](#facorr) step is not required <br/> 1. Change tab to `Intensity corrections` <br/> 2. Check that `Autosave` is on <br/> 3. Click `Apply to` all <br/> 4. Click on `Calculate all`
|
|
|
<img src = './images/FA/save_res.png' width = "400px" > <br/> <img src = './images/FA/res_file_interval.png'> | 5. Change tab to `Save, export and report` <br/> 6. Click on `Save all` to save single taces <br/>7. Click on `FA format` to save a summary result table (here 2c.res)<br/> 8. Compute offset values from averages of columns named *Interval Ch1/2* in the *res* file
|
|
|
|
|
|
|
|
|
|
... | ... | @@ -103,37 +123,78 @@ The [**Offset**](#offset) and [**Crosstalk**](#xtalk) parameters are used for co |
|
|
* Load data into FA
|
|
|
* Change to `Intensity corrections` and enter the offset values for Ch1 and Ch2
|
|
|
* Calculate all corrections and save *res* table
|
|
|
* Compute crosstalk from Ch1 (to Ch2) from $`\frac{\text{IntervalCh2} - \text{OffsetCh2}}{\text{IntervalCh1} - \text{OffsetCh1}}`$
|
|
|
* Compute crosstalk from Ch1 (to Ch2) from average of columns in *res* file $`\frac{\text{IntervalCh2} - \text{OffsetCh2}}{\text{IntervalCh1} - \text{OffsetCh1}}`$
|
|
|
|
|
|
```math
|
|
|
\frac{\text{IntervalCh2} - \text{OffsetCh2}}{\text{IntervalCh1} - \text{OffsetCh1}}
|
|
|
```
|
|
|
|
|
|
|
|
|
| |
|
|
|
:---: | ----
|
|
|
<img src = './images/FA/xtalk_correction_idx.png' width = "400px" > | [Import](#faload) FCS measurements from WT cells into FA. The [Modify and correlate](#facorr) step is not required <br/> 1. Change tab to `Intensity corrections` <br/> 2. Check that `Autosave` is on <br/> 3. Click `Apply to` <br/> 4. Click on `Calculate all`
|
|
|
<img src = './images/FA/save_res.png' width = "400px" > | 5. Change tab to `Save, export and report` <br/> 6. Click on `Save all` to save single taces <br/>7. Click on `FA format` to save a summary result table (here 2c.res)
|
|
|
### [<img src='./images/up.png'>](#back) 4. <a name=facorrfa> </a> Compute correction factors
|
|
|
Correction parameters are used to then compute correction computed protein numbers in the effective volume for bias induced by background, cross-talk, and photobleaching. The corrected values can be computed from the fitted values and correction factors stored in the result table. Correction factors need to be calculated for fluorescent proteins.
|
|
|
Intensity corrections (for fluorescent proteins only)| |
|
|
|
:--: | ---
|
|
|
<img src = './images/FA/xtalk_correction_idx.png' width = "400px" > | [Import](#faload) FCS measurements from cells expressing a single fluorophore. Proceed [Modify and correlate](#facorr) <br/> 1. Change tab to `Intensity corrections` <br/> 2. Enter previously estimated Offset value for Ch1 <br/> 3. For two color FCS enter Offset value for Ch2 <br/> 4. Click on `Apply to` all <br/> 5. Click on `Calculate all`
|
|
|
|
|
|
### [<img src='./images/up.png'>](#back) 4. <a name=facorrfa> </a> Compute correction factors and save table
|
|
|
Correction factors are used to correct the protein numbers for bias induced by background, cross-talk (for 2 color FCS), and photobleaching. The corrected values can be calculated from the fitted values and correction factors stored in the result table.
|
|
|
|
|
|
## Summary of typical parameters used for computing correlation functions and correction factors
|
|
|
The procedure is similar as for the Offset and Crosstalk value computation:
|
|
|
|
|
|
* Load the data into FA
|
|
|
* Switch to the `Intensity corrections` tab
|
|
|
* Enter the `Offset [kHz]` for Ch1 (and Ch2 for two color FCS)
|
|
|
* Enter the `Crosstalk` value in case of two color FCS
|
|
|
* Click `Appky to`
|
|
|
* Click on `Calculate all`
|
|
|
* Save results to *res* file
|
|
|
|
|
|
Parameters | Values |
|
|
|
:--- | :---
|
|
|
Base freq. (Hz) | 1.000.000 (dye) or 100.000 (protein)
|
|
|
Composition | Ch1 <> Ch2
|
|
|
Offset (kHz) | 1-5 kHz <br/> Must be measured
|
|
|
Cross-talk (Ch1-Ch2), when using two color FCS | 0.04-0.05 (Ch1: mEGFP, Ch2: mCherry). <br/> Value depends on protein and optical settings. Must be measured.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
## Fitting data using FA
|
|
|
Using FA the user can freely decide which model to use and which parameters are let free or fixed. Alternatively the user can use the matlab routine that performs several fitting steps at once using some default parameter sets optimal for typical fluorescent proteins and dyes. In particular for the dye it is difficult to estimate
|
|
|
Using FA the user can freely decide which model to use and which parameters are fitted or not for up to two consecutive fitting settings (Run A and Run B). For the fitting of the fluorescent dye and fluorescent protein we use Run A an Run
|
|
|
|
|
|
|
|
|
In the parameter choice section the fate of a parameter can be changed by clicking on the boxes below a parameter.
|
|
|
* ** Free**: Parameter will be fitted
|
|
|
* ** Fixed**: Parameter will is kept constant
|
|
|
* **Link/Change**: This is used for specifying changes in the fitting state of a parameter between RunA and RunB.
|
|
|
|
|
|
Alternatively the user can use the matlab workflow that combines several steps of the fitting. A summary of initial parameter values is given below. Note that the parameters for the fluorescent dye and protein are not the same.
|
|
|
### Fluorescent protein
|
|
|
| |
|
|
|
:---: | ----
|
|
|
**RUN A** <br/><img src = './images/FA/protein_runA.png' width = "400px" >**RUN B** <br/><img src = './images/FA/protein_runB.png' width = "400px" > | 1. Set `Weighted fit`, `Initial guess` on. To perform RUNB the `2nd run` needs to be on. This may need to be reset after a run <br/> 2. Verify that model is *two-component anomalous diffusion with fluorescent protein-like blinking* for all measurements <br/> 3. Set the parameter values and fitting options as shown in the table <br/> 4. Verify that `Autosave` is on <br/> 5. Click `Apply to` all <br/> 6. Switch to settings for RUNB <br/> Change parameter table and values Click on `Fit all`
|
|
|
|
|
|
|
|
|
|
|
|
Set the model *two-component anomalous diffusion with fluorescent protein-like blinking*
|
|
|
|
|
|
Parameters | Round A | Round B
|
|
|
:--- | :--- | :---
|
|
|
Model | two-comp. anomalous diffusion with triplet-like blinking |
|
|
|
N | 10 | 10
|
|
|
thetaT | 0.2 | *0.2*
|
|
|
tauT | *10* | *100*
|
|
|
f1 | *1* | 0.5
|
|
|
tauD1 | 20 | 500
|
|
|
alpha1 | *1* | *1*
|
|
|
tauD2 | NA | 5000
|
|
|
alpha2 | NA | *1*
|
|
|
kappa | 5.5 | 5.5
|
|
|
offset | *0* | *0*
|
|
|
|
|
|
Parameters (round B) | Fluorescent dye (1 component)| Fluorescent protein (2 components)
|
|
|
:--- | :---: | :---:
|
|
|
N | 10 | 10
|
|
|
thetaT | 0.2 | 0.2
|
|
|
tauT | 10 | **100**
|
|
|
f1 | **1** | 0.5
|
|
|
tauD1 | 20 | 500
|
|
|
alpha1 | **1** | 1
|
|
|
tauD2 | NA | 5000
|
|
|
alpha2 | NA | 1
|
|
|
kappa | 5.5 | 5.5
|
|
|
offset | **0** | **0**
|
|
|
### Fitting of the fluorescent dye
|
|
|
|
|
|
|
|
|
### Fit of fluorescent dye data
|
|
|
For the fitting of autocorrelation curves to a physical model
|
... | ... | @@ -144,7 +205,7 @@ Fit dye data | | |
|
|
|
|
|
|
|
|
|
|
|
## Fitting of data using Matlab
|
|
|
## Fitting correlation curves using Matlab
|
|
|
|
|
|
|
|
|
|