User Tools

Site Tools


Courses -2016
Courses 2018-
Downloads (-2016)
Downloads (2016-)
RSS aggregates

EMBL BioImage Data Analysis



—- Contact
CMCI Alumni
EMBL Heidelberg
EMBL Intranet

Popularity Ranking

Timeline of @cmci_

cmci_ avatar

RT @J_Cell_Sci: Registration is now open for our 2023 Journal Meeting 'Imaging Cell Dynamics', organised by Erika Holzbaur @ErikaHolzbaur,…
About 7 hours, 9 mins ago by: Kota Miura (@cmci_)

cmci_ avatar

RT @joachimgoedhart: This online resource that introduces data visualization using #rstats & ggplot2 has 12 complete protocols available fo…
About 11 hours, 59 mins ago by: Kota Miura (@cmci_)

cmci_ avatar

The code for image analysis is the "provided upon request" style ... it would be much better to have it openly acce…
About 13 hours, 24 mins ago by: Kota Miura (@cmci_)

cmci_ avatar

RT @SXS_biophysics: Why do cells move faster in higher viscosity environments? Debonil's paper on this question is out:
About 6 days, 3 hours ago by: Kota Miura (@cmci_)

cmci_ avatar

The first day of online “defragmentation” NEUBIAS academy / EOSC-Life ended, lead by @RogerDAntuono . Great job Roc…
About 6 days, 5 hours ago by: Kota Miura (@cmci_)

cmci_ avatar

RT @hpke1980: Our annotator software is coming soon on BiorXiv. Check it out! #napari @napari_imaging #singlecell h…
About 6 days, 7 hours ago by: Kota Miura (@cmci_)

This is an old revision of the document!

Notes: FRAP internal course

April 20, 2010 @EMBL Kota Miura
will be also further added by Sebastian Huet and Christian Tischer

Google Chrome or Firefox (version > 3.6) is recommended for properly viewing math equations.


in vivo protein kinetics could be analyzed in two ways: measuring particular movement or averaged movement. By tracking labeled single protein molecules, it is possible to estimate their diffusion and transport behavior. Single molecule studies of membrane proteins enabled us to analyze how they are organized with their dynamics, through corrals of membrane territory. Motor protein moving along cytoskeletal tracks were analyzed in detail how they convert chemical energy into physical force was only possible by probing their singular movement and steps. While single particle tracking requires high-temporal and spatial resolution setup for analysis, analysis of averaged movement measured through fluorescence intensity dynamics, could be achieved with larger spatial and temporal resolution (typically in micrometer scale).

Here, we focus on one of such averaged movement analysis technique, Fluorescence Recovery After Photobleaching (FRAP). We first start with a simpler case of monitoring averaged movement that does not need to bleach.

Fluorescence intensity and Protein Dynamics

Measurement of VSV-G protein exit dynamics
First-order Chemical Reaction

Increase in intensity at observed area could be measured to know the net increase in the protein at that region. To calculate biochemical kinetics, we can apply traditional biochemical kinetics. Example case: Kinetics of VSVG protein accumulation to ER exit site.

<jsmath>{dI(t)\over dt}=k_{on}[VSVG_{free}]- k_{off}[VSVG_{ERES}]</jsmath>


  • <jsm>k_{on}</jsm> is the binding rate of VSVG protein to ER exit site
  • <jsm>[VSVG_{free}]</jsm> is the concentration of unbound VSVG protein
  • <jsm>k_{off}</jsm> is the dissociation rate of VSVG protein from ER exit site
  • <jsm>[VSVG_{ERES}]</jsm> is the density of VSVG protein bound to the ER exit site

During the initial phase of binding, when there is almost no VSVG protein bound to ER exit site, we can approximate the initial speed of the density increase at ERES site depends only on binding reaction and assume <jsm>k_{off}[VSVG_{ERES}]\simeq0</jsm>. Then

<jsmath> {dI(t)\over{dt}}=k_{on}[VSVG_{free}] </jsmath>

We then are be able to simply calculate the slope of initial increase of intensity, measure the free VSVG intensity and then calculate <jsm>k_{on}</jsm>.

For details, see Runz et al (2006).

FRAP Simple Measures

Half-Max t and Mobile/Immobile fraction from FRAP curve

Unlike the example explained above, dynamics of protein are not observable in many cases. Even though proteins are exchanging in system, the flux of protein constituting the system is not evident if the flux is steady and constant (e.g. liver). In such cases, we need to some how experimentally measure the system. One way is FRAP. In FRAP, we bleach some population of fluorescence-labeled protein and evaluate the mobility of the protein. Typically, we use confocal microscopy and bleach fluorescence of small area of the system by short pulse of strong laser beam and measure changes in the fluorescence intensity over time at that bleached spot over time. Detail on these measurement protocol has been presented in Stefan and Yury's talks (link?). Here, we focus on how to analyze the curve we obtained by such measurements.

From measured temporal changes in the intensity at the bleached Region of Interest (this curve indeed is the Fluorescence Recovery After Photobleaching, FRAP)

Half Max and Mobile-Immobile fraction
Fitting to Exponential curve

FRAP Measurements based on Modelling

FRAP curve reflects the mobility of proteins. In dilute solution of single protein species, mobility of protein could probably be considered as pure-diffusion. But in many cases, this is not the case. The mobility is affected by the system.

  • reaction with other proteins
  • geometry of the system, that constrains the mobility
  • active transport process

By modeling how the mobility is (generate some hypothesis how the protein mobility is affected in the system), we can set up equation/s to that should describe the FRAP curve. This is done by fitting the theoretical curves to the experimental curves. By evaluating the goodness of fit, we can discuss which models would be the most likely hypothesis.

Currently we have more-or-less standardized protocol to analyze FRAP curve. Starting with simple model of diffusion, we test the fit of different curves and proceed to more complex models. See next section for the protocol.

Choosing models based on trial fittings

Maybe Sebastian's flow chart here.

Diffusion Recovery

Diffusion Recovery equation by Soumpasis, 1984

pure diffusion

Theoretical curve of the diffusion mediated fluorescence recovery was proposed by Soumpasis (1984) and has been widely used.
<jsmath> f(t)=e^{- \frac{\tau_D}{2t}}\left(I_{0}(\frac{\tau_D}{2t})+I_{1}(\frac{\tau_D}{2t})\right) </jsmath> This theoretical equation assumes:

  • 2D
  • circular (cylindrical) bleaching

when above equation could be fitted nicely (evaluated by goodness of fit, such as Pearson's coefficient r or gamma-Q value), one could calculate diffusion coeffecient by using the obtained <jsm>\tau_D</jsm> and radius of the circular ROI <jsm>w</jsm>.
<jsmath> D=\frac{w^{2}}{\tau_D} </jsmath>

effective diffusion

(almost diffusion)

special cases: anomalous diffusion

Reaction Dominant Recovery

Diffusion occurs at the very beginning. Recovery curve is dominated by reaction.

Protein cluster and fluorescence recovery
Modeling Reaction Dominant recovery

Reaction Dominant Recovery with Immobile Binding Partner

Illustrated view of Fluorescence recovery at immobile binding partner
Modeling fluorescence recovery at immobile binding partner 01
Modeling fluorescence recovery at immobile binding partner 02

Diffusion and Reaction combined Recovery

Diffusion-Reaction Combined Model

Diffusion and Transport combined Recovery

Hallen and Endow, 2009

Diffusion and Reaction, along with Spatial Context

Example case considering geometrical and structural parameters that are affecting the recovery, by Sbalzarini et al

Tools for FRAP Analysis


  • ImageJ Plugins: Two available ImageJ plugins for FRAP does intensity measurement and normalization. Fitting module is not implemented.
    • Import data output from Zeiss, Leica, Olympus measurements and do FRAP fitting.
    • Does measurement and fitting. Sprague et al. (2004) Reaction-Diffusion Full model is implemented.

Two Advanced strategies: Analytical and Numerical Approaches

Requires your own coding, customization

Analytical Approach

Sprague et. al. (explained above) is an example case of analytically solving the model for the fitting.

ODE Simulation

Image based simulation with ODE solver


  • Numerical analysis based on ODE. Spatial context.


  • Berkley Madonna

Particle Simulation

One of Particle Simulator, Smodldyn
Custom tool for Particle Simualation

Particle Simulation Packages

  • Smoldyn
  • GridCell
  • MCell

Reviews, Protocols and textbooks

  • McNally, J.G. (2008) Quantitative FRAP in analysis of molecular binding dynamics in vivo. Methods in Cell Biology 2008;85:329-51
    • step-by-step protocol for analyzing FRAP. Flow chart for determining recovery model is useful.
  • Sprague, B. L. and McNally, J. G.(2005). FRAP analysis of binding: proper and fitting. Trends Cell Biol 15, 84-91.
  • Phair, R. D., Gorski, S. A. and Misteli, T. (2004). Measurement of dynamic protein binding to chromatin in vivo, using photobleaching microscopy. Methods Enzymol 375, 393-414.
    • Already a bit of out-dated but “Double normalization” introduced in this protocol is widely used.
  • Jacquez, J. A. (1972). Compartmental analysis in biology and medicine: Elsevier
    • If you want to study basics about calculation of kinetics.

Research Papers

documents/100420frapinternal.1271937578.txt.gz · Last modified: 2016/05/24 12:46 (external edit)