... | ... | @@ -38,202 +38,14 @@ td+td{ |
|
|
In this WiKi we provide examples on how to analyze and process FCS and imaging data generated on Zeiss LSM microscopes using [microscopy pipeline constructor](#https://git.embl.de/politi/mypic) or
|
|
|
[FCSRunner](#https://git.embl.de/politi/fcsrunner). The processed data can be used to generate a calibration curve for converting pixel fluorescence intensities to concentrations.
|
|
|
|
|
|
[**Generate bleach corrected correlation functions using Fluctuation Analyzer 4G**](#gencorr)
|
|
|
[**Fit correlation data with Fluctuation Analyzer**](#fitfa)
|
|
|
[**Fit correlation data with Matlab**](#fitmatlab)
|
|
|
|
|
|
|
|
|
|
|
|
## <a name = 'gencorr'> </a> 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 offset and crosstalk parameters](#facorrpar)
|
|
|
|
|
|
4. [Compute correction factors and save table](#facorrfa)
|
|
|
|
|
|
See also [Typical values used for computing correlation functions and correction factors](#faval)
|
|
|
|
|
|
> **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.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
### [<img src='./images/up.png'>](#back) <a name=faload></a>1. Loading data into FA
|
|
|
|
|
|
|
|
|
**FCS raw data**: FA requires the raw photon-counting data. For Zeiss LSM with ZEN (black edition) to save raw data use the `Confocor Options` menu in the `Maintain` tab. When saving the fcs recording be sure to save as type `Fcs files with raw data (*.fcs)`
|
|
|
<div align = "center" > <img src = './images/FA/save_raw_data.PNG' ><img src = './images/FA/save_as_type.PNG' > </div>
|
|
|
|
|
|
|
|
|
<a name='sessionname'></a>**Session name** : Alphanumeric identifier used for re-loading processed data. Typically use **1c** when fitting one component model (e.g. fluorescent dye) or **2c** when fitting a two component model (e.g. fluorescent protein)
|
|
|
|
|
|
| | |
|
|
|
:---: | ----
|
|
|
<img src = './images/FA/import1_idx.png' width = "400px" > <br/> Import settings APD/Confocor and GaAsP <br/> <img src = './images/FA/load_APD.png' width = "200px" > <img src = './images/FA/load_GaAsP.png' width = "200px" > | 1. Choose appropriate file format <br/> 2. Click on import settings. The channel with the lowest wavelength should be Ch1. <br/> FA_Ch1/2 = None if only one channel has been acquired <br/> 2a. With APD FA_Ch1 = Ch2 and FA_Ch2 = Ch1.<br/> 2b. With GaAsP FA_Ch1 = ChS1 and FA_Ch2 = ChS2 <br/> 3. Select file path. <br/> 4. Click if data in all subdirectories need to be processed. <br/> 5. Add data to FA
|
|
|
<img src = './images/FA/import2_idx.png' width = "400px" > | 6. Specify a [session name](#sessionname) <br/> 7. Click on `Check files` <br/> 8. If data `Quick check` if data has already been loaded once
|
|
|
|
|
|
|
|
|
|
|
|
### [<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` all <br/> 7. Click on `Calculate all`
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
### [<img src='./images/up.png'>](#back) 3. <a name=facorrpar> </a>Compute offset and crosstalk parameters
|
|
|
[**Generate bleach corrected correlation functions using Fluctuation Analyzer 4G**](./FA_Load_and_Correlate.md)
|
|
|
|
|
|
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 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` 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
|
|
|
|
|
|
|
|
|
|
|
|
* <a name='xtalk'>**Crosstalk (Ch1 to Ch2)**</a>: This parameter is required when two fluorophores are measured, e.g. to quantify the cross-correlation. To quantify the crosstalk cells expressing each one of the two fluorophores are measured.
|
|
|
* Acquire data from cells expressing only one fluorophore (5-10 cells)
|
|
|
* 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 average of columns in *res* file $`\frac{\text{IntervalCh2} - \text{OffsetCh2}}{\text{IntervalCh1} - \text{OffsetCh1}}`$
|
|
|
|
|
|
|
|
|
|
|
|
| |
|
|
|
:---: | ----
|
|
|
<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.
|
|
|
|
|
|
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
|
|
|
|
|
|
### [<img src='./images/up.png'>](#back) 5. <a name=faval> </a> Typical values used for computing correlation functions and correction factors
|
|
|
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.
|
|
|
|
|
|
|
|
|
## Fitting data using FA
|
|
|
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 (**RunA** and **RunB**). Typically one uses different models and initial values when fitting FCS traces from a [fluorescent protein](#fitprot) or a [fluorescent dye](#fitdye).
|
|
|
|
|
|
|
|
|
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.
|
|
|
|
|
|
### <a name=fitprot></a> Fluorescent protein
|
|
|
|
|
|
|
|
|
1. Start fluctuation analyzer, [load data](#faload), and [compute correlations] (#facorr). Set `Weighted fit`, `Initial guess` on. To perform RUNB the `2nd run` needs to be on.This may need to be reset after a fit has been peformed
|
|
|
**RUN A** <br/><img src = './images/FA/protein_runA_idx.png' >
|
|
|
|
|
|
2. Verify that model is *two-component anomalous diffusion with fluorescent protein-like blinking* for all measurements
|
|
|
3. Set the parameter values and fitting options as shown in the table
|
|
|
4. Change kappa according to value appropriate for objective and microscope settings as determined using the fluorescent dye
|
|
|
5. Verify that `Autosave` is on
|
|
|
6. Click `Apply to` all
|
|
|
|
|
|
**RUN B** <br/><img src = './images/FA/protein_runB_idx.png' >
|
|
|
7. Switch to settings for RUNB
|
|
|
8. Set the parameter values and fitting options as shown in the table for RUNB
|
|
|
9. Click `Apply to` all
|
|
|
10. Click `Fit all`
|
|
|
|
|
|
|
|
|
|
|
|
### <a name=fitdye></a> Fluorescent dye
|
|
|
|
|
|
|
|
|
1. Start fluctuation analyzer, [load data](#faload), and [compute correlations](#facorr) using the highest frequency of 1.000.0000 Hz. Set `Weighted fit`, `Initial guess` on. To perform RUNB the `2nd run` needs to be on. This may need to be reset after a fit has been peformed
|
|
|
**RUN A** <br/><img src = './images/FA/dye_runA.png' >
|
|
|
2. For detectors with after-pulsing (e.g. APD's) set the lower fit boundary to 2 μs
|
|
|
3. Set the upper fit boundary to a value suitable for a fast-diffusing molecule (~ 10230 μs )
|
|
|
4. Change the model to *two-component anomalous diffusion with triplet-like blinking*. Press `Apply to` all in the `Change model ...` menu
|
|
|
5. Set the parameters for RoundA as shown
|
|
|
6. Verify that `Autosave` is on
|
|
|
6. Click `Apply to` all
|
|
|
|
|
|
**RUN B** <br/><img src = './images/FA/dye_runB.png' >
|
|
|
7. Switch to settings for RUNB
|
|
|
8. Set the parameter values and fitting options as shown in the table for RUNB
|
|
|
9. Click `Apply to` all
|
|
|
10. Click `Fit all`
|
|
|
|
|
|
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 correlation data with Fluctuation Analyzer**](./FA_Fit_FCS.md)
|
|
|
|
|
|
[**Fit correlation data with Matlab**](#fitmatlab)
|
|
|
|
|
|
|
|
|
|
|
|
## Fitting correlation curves using Matlab
|
|
|
|
|
|
|
|
|
|