# Bioimage Analysis Wiki

### Site Tools

documents:100420frapinternal

# Differences

This shows you the differences between two versions of the page.

 documents:100420frapinternal [2010/04/22 15:38]kota documents:100420frapinternal [2020/11/26 08:05] (current)kota [Fluorescence intensity and Protein Dynamics] Both sides previous revision Previous revision 2020/11/26 08:05 kota [Fluorescence intensity and Protein Dynamics] 2016/05/24 15:53 kota 2016/05/24 12:46 external edit2010/05/11 13:12 kota 2010/04/26 09:49 kota 2010/04/26 09:07 christian 2010/04/23 21:49 kota read again to correct terrible writing2010/04/23 09:15 kota 2010/04/23 08:05 kota 2010/04/23 07:24 kota update on pure diffusion2010/04/22 15:38 kota 2010/04/22 15:38 kota 2010/04/22 15:21 kota 2010/04/22 15:14 kota 2010/04/22 15:10 kota 2010/04/22 15:08 kota 2010/04/22 15:01 kota 2010/04/22 14:22 kota 2010/04/22 12:05 kota 2010/04/22 12:00 kota 2010/04/22 11:59 kota 2010/04/22 11:56 kota 2010/04/22 11:46 kota 2010/04/22 11:45 kota 2010/04/22 11:39 kota 2010/04/22 11:36 kota soumpasis2010/04/22 11:24 kota 2010/04/22 11:07 kota 2010/04/22 10:47 kota 2010/04/22 10:42 kota 2010/04/22 08:23 kota 2010/04/22 08:23 kota 2010/04/22 07:39 kota 2010/04/22 07:38 kota 2010/04/21 15:42 kota 2010/04/21 15:41 kota added slide for vsvg exit study2010/04/21 15:05 kota Next revision Previous revision 2020/11/26 08:05 kota [Fluorescence intensity and Protein Dynamics] 2016/05/24 15:53 kota 2016/05/24 12:46 external edit2010/05/11 13:12 kota 2010/04/26 09:49 kota 2010/04/26 09:07 christian 2010/04/23 21:49 kota read again to correct terrible writing2010/04/23 09:15 kota 2010/04/23 08:05 kota 2010/04/23 07:24 kota update on pure diffusion2010/04/22 15:38 kota 2010/04/22 15:38 kota 2010/04/22 15:21 kota 2010/04/22 15:14 kota 2010/04/22 15:10 kota 2010/04/22 15:08 kota 2010/04/22 15:01 kota 2010/04/22 14:22 kota 2010/04/22 12:05 kota 2010/04/22 12:00 kota 2010/04/22 11:59 kota 2010/04/22 11:56 kota 2010/04/22 11:46 kota 2010/04/22 11:45 kota 2010/04/22 11:39 kota 2010/04/22 11:36 kota soumpasis2010/04/22 11:24 kota 2010/04/22 11:07 kota 2010/04/22 10:47 kota 2010/04/22 10:42 kota 2010/04/22 08:23 kota 2010/04/22 08:23 kota 2010/04/22 07:39 kota 2010/04/22 07:38 kota 2010/04/21 15:42 kota 2010/04/21 15:41 kota added slide for vsvg exit study2010/04/21 15:05 kota 2010/04/21 14:55 kota 2010/04/21 14:53 kota 2010/04/21 14:04 kota 2010/04/21 12:56 kota 2010/04/21 12:08 kota 2010/04/21 12:01 kota 2010/04/21 11:50 kota 2010/04/21 09:32 kota 2010/04/21 08:34 kota 2010/04/21 08:25 kota 2010/04/21 08:23 kota 2010/04/21 08:20 kota 2010/04/21 08:11 kota 2010/04/21 08:09 kota Line 1: Line 1: - ====== Notes: FRAP internal course ====== + ====== Lecture Notes: FRAP internal course ====== - April 20, 2010 @EMBL + April 20, 2010 @EMBL\\ Kota Miura\\ Kota Miura\\ + will be also further added by Sebastian Huet and Christian Tischer\\ will be also further added by Sebastian Huet and Christian Tischer\\ Line 8: Line 9: ===== Introduction ===== ===== Introduction ===== - 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). + in vivo protein kinetics could be analyzed in two ways: measuring particular movement or averaged movement. By tracking labeled single protein molecules, we could estimate their diffusion and transport behavior. Such single molecule studies of membrane proteins, for example, enabled us to analyze how they are organized with their dynamics, such as boundary for movement constrained by membrane corrals. Motor protein moving along cytoskeletal tracks were analyzed in detail to know how they convert chemical energy into physical force. This 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 by temporal changes in fluorescence intensity, 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 the explanation with a simpler case of monitoring averaged movement that does not need to bleach. - 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. - ~~NOCACHE~~ ===== Fluorescence intensity and Protein Dynamics ===== ===== Fluorescence intensity and Protein Dynamics ===== - [{{ :documents:vsvg_exit.jpg?150|Measurement of VSV-G protein exit dynamics}}] + [{{ :documents:vsvg_exit.jpg?200|Measurement of VSV-G protein exit dynamics}}] - [{{ :documents:firstorderchemicalreaction.jpg?150| First-order Chemical Reaction}}] + [{{ :documents:firstorderchemicalreaction.jpg?200| 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. + Increase in intensity at observed area could be measured to know the net increase in the protein at that region. To characterize this dynamics, we can apply traditional biochemical kinetics. Example case: Kinetics of VSVG protein accumulation to ER exit site. - {dI(t)\over dt}=k_{on}[VSVG_{free}]- k_{off}[VSVG_{ERES}] + $${dI(t)\over dt}=k_{on}[VSVG_{free}]- k_{off}[VSVG_{ERES}]$$ Here, Here, - * k_{on} is the binding rate of VSVG protein to ER exit site + * $k_{on}$ is the binding rate of VSVG protein to ER exit site - * [VSVG_{free}] is the concentration of unbound VSVG protein + * $[VSVG_{free}]$ is the concentration of unbound VSVG protein - * k_{off} is the dissociation rate of VSVG protein from ER exit site + * $k_{off}$ is the dissociation rate of VSVG protein from ER exit site - * [VSVG_{ERES}] is the density of VSVG protein bound to the ER exit site + * $[VSVG_{ERES}]$ 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 k_{off}[VSVG_{ERES}]\simeq0. Then + 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: $k_{off}[VSVG_{ERES}]\simeq0$. \\Then - + $${dI(t)\over{dt}}=k_{on}[VSVG_{free}] {dI(t)\over{dt}}=k_{on}[VSVG_{free}] - +$$ - Since there are enough free VSVG, we consider that [VSVG_{free}] is constant, we are able to simply calculate the slope of initial increase of intensity, measure the free VSVG intensity and then calculate k_{on}. + Since there are enough free VSVG, we consider that $[VSVG_{free}]$ is constant, we are able to simply calculate the slope of initial increase of intensity, measure the free VSVG intensity and then calculate $k_{on}$. For details, see [[http://www.ncbi.nlm.nih.gov/pubmed/16794576?dopt=Abstract |Runz et al (2006)]]. For details, see [[http://www.ncbi.nlm.nih.gov/pubmed/16794576?dopt=Abstract |Runz et al (2006)]]. Line 38: Line 39: ===== FRAP Simple Measures ===== ===== FRAP Simple Measures ===== [{{ :documents:halfmax_and_mobile.jpg?150| Half-Max t and Mobile/Immobile fraction from FRAP curve}}] [{{ :documents:halfmax_and_mobile.jpg?150| 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. + Unlike the example shown above, dynamics of protein are not observable to eyes (through microscope) in many cases. Even though proteins are exchanging in system, the flux of protein constituting the system is not evident if the in/out flux of protein is steady and constant (e.g. liver). In such cases, we need to some how experimentally treat 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 following changes in the fluorescence intensity 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 through 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) + From measured temporal changes in the intensity at the bleached Region of Interest (this curve indeed is the Fluorescence Recovery After Photobleaching, FRAP) we can measure two parameters which represents speed of recovery, and fraction of molecules that is moving around in the system. - Half Max and Mobile-Immobile fraction\\ + * Half Max - Fitting to Exponential curve + * Mobile-Immobile fraction + + Fitting the curve to exponential equation eases us to calculate these parameters. Half-Max value (time) is rather qualitative value, but is a simple and straightforward index for comparing different systems. ===== FRAP Measurements based on Modelling ===== ===== 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. + FRAP curve reflects the mobility of proteins. In dilute solution of with single protein solute, mobility of protein could probably be considered as pure-diffusion. But in many cases, this is does not hold. The mobility is often affected by the system. * reaction with other proteins * reaction with other proteins * geometry of the system, that constrains the mobility * geometry of the system, that constrains the mobility * active transport process * 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. + 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 hypothesize what is the bases of FRAP curve. To test the hypothesis, we fit the experimental curves with theoretical curve. By evaluating the goodness of fit, we can discuss which models would be the most likely hypothesis. If the fit is good, then we could know the value of biochemical parameters which governs the recovery curve. - 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. + Currently we have more-or-less standardized protocol to analyze FRAP curves. 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. Line 62: Line 66: === pure diffusion === === pure diffusion === Theoretical curve of the diffusion mediated fluorescence recovery was proposed by Soumpasis (1984) and has been widely used. \\ Theoretical curve of the diffusion mediated fluorescence recovery was proposed by Soumpasis (1984) and has been widely used. \\ - + $$f(t)=e^{- \frac{\tau_D}{2t}}\left(I_{0}(\frac{\tau_D}{2t})+I_{1}(\frac{\tau_D}{2t})\right) f(t)=e^{- \frac{\tau_D}{2t}}\left(I_{0}(\frac{\tau_D}{2t})+I_{1}(\frac{\tau_D}{2t})\right) - +$$ This theoretical equation assumes: This theoretical equation assumes: * 2D * 2D * circular (cylindrical) bleaching * 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 \tau_D and radius of the circular ROI w.\\ + 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 $\tau_D$ and radius of the circular ROI $w$.\\ - + $$D=\frac{w^{2}}{\tau_D} D=\frac{w^{2}}{\tau_D} - +$$ + + For strip-ROI bleaching, empirical formula used by Ellenberg et al. (1997) could be used, and is also possible to use Gaussian curve fitting that **Christian Tischer** developed. For Christian's method, [[:documents:100426FRAPgaussfit]]. === effective diffusion === === effective diffusion === Line 77: Line 83: (almost diffusion) (almost diffusion) - === special cases: anomalous diffusion === + === anomalous diffusion === ==== Reaction Dominant Recovery ==== ==== Reaction Dominant Recovery ==== [{{ :documents:reactiondominant01.jpg?150|Protein cluster and fluorescence recovery}}] [{{ :documents:reactiondominant01.jpg?150|Protein cluster and fluorescence recovery}}] [{{ :documents:reactiondominant02.jpg?150|Modeling Reaction Dominant recovery}}] [{{ :documents:reactiondominant02.jpg?150|Modeling Reaction Dominant recovery}}] - If molecule under study is binding/unbinding with other molecular species, FRAP curve is affected by this interaction. There are two cases on how these two events, diffusion and reaction, are combined in the curve. We could think of two cases. + If molecule under study is binding/unbinding with other molecular species, FRAP curve is affected by these interactions. There are two cases on how these two events, diffusion and reaction, are combined in the curve. We could think of two cases. - - **Reaction-dominant recovery** Cases when reaction binding rate (meaning "number of molecules that binds to the binding partner per second") is much lower than the diffusion rate. In such cases, FRAP-bleached field will be become brighter due to diffusion, and then there will be a slower recovery of intensity due to none-bleached fluorescence exchange with the bleached fluorescence. In such cases, we could consider that the recovery curve is dominated by reaction since duration of diffusion-recovery phase and reaction-recovery phase is much shorter for the diffusion-recovery phase, these two phases could be assumed to be separated(call this **reaction-dominant** or **diffusion uncoupled**; this section). + - **Reaction-dominant recovery** Cases when reaction binding rate (meaning "number of molecules that binds to the binding partner per second") is much lower than the diffusion rate. In such cases, FRAP-bleached field will first become brighter due to diffusion, and then there will be a slower recovery of intensity due to none-bleached fluorescence exchange with the bleached fluorescence. In such cases, we could consider that the recovery curve is dominated by reaction since duration of diffusion-recovery phase is much shorter compared to reaction-recovery phase. In such case, these two phases are considered to be separable (call this **reaction-dominant** or **diffusion uncoupled**; this section). We even might not "see" the diffusion recovery phase, which requires high temporal resolution capturing. - - **Reaction-Diffusion Recovery** Other cases would be when durations of diffusion-recovery phase and reaction-recovery phase are comparable. Then recovery curve consists of a combination of fluorescence that came in with diffusion, and also by the binding of fluorescence molecule at the FRAP bleached field (**"reaction-diffusion"** or **"diffusion-coupled"** see the section below). + - **Reaction-Diffusion Recovery** Other cases would be when durations of diffusion-recovery phase and reaction-recovery phase are with comparable duration. Then recovery curve consists of a combination of fluorescence that came in with diffusion, and also by the binding of fluorescence molecule at the FRAP bleached field (**"reaction-diffusion"** or **"diffusion-coupled"** see the section below). - For a simple chemical reaction with singular type of interaction, we could again think of the reaction model that was already explained above, the first-order chemical reaction could be modeled as compartment system (see figure right)\\ + For a simple chemical reaction with singular type of interaction, we could again think of the reaction model that was already explained above, the first-order chemical reaction modeled as a compartment system (see figure right)\\ - + $$\frac {df(t)} {dt} = k_{on}[free] - k_{off}[bound] \frac {df(t)} {dt} = k_{on}[free] - k_{off}[bound] - +$$ where where - * k_{on} Binding constant + * $k_{on}$ Binding constant - * k_{off} Dissociation constant + * $k_{off}$ Dissociation constant - * [free] Density of free molecules + * $[free]$ Density of free molecules - * [bound] Density of bound-molecules + * $[bound]$ Density of bound-molecules We solve the differential equation We solve the differential equation - + $$f(t)=A(1-e^{- \tau t}) f(t)=A(1-e^{- \tau t}) - +$$ where where - * \tau = k_{on} + k_{off} + * $\tau = k_{on} + k_{off}$ - * A = \frac {k_{on}}{k_{on} + k_{off}} + * $A = \frac {k_{on}}{k_{on} + k_{off}}$ === Reaction Dominant Recovery with Immobile Binding Partner=== === Reaction Dominant Recovery with Immobile Binding Partner=== Line 108: Line 115: [{{ :documents:reactiondominant04_immobileentity.jpg?150| Modeling fluorescence recovery at immobile binding partner 01}}] [{{ :documents:reactiondominant04_immobileentity.jpg?150| Modeling fluorescence recovery at immobile binding partner 01}}] [{{ :documents:reactiondominant05_immobileentity.jpg?150| Modeling fluorescence recovery at immobile binding partner 02}}] [{{ :documents:reactiondominant05_immobileentity.jpg?150| Modeling fluorescence recovery at immobile binding partner 02}}] - Next we modify above model to consider a situation a bit more frequently we see in cell biology. The protein we are analyzing is either freely diffusion in cytoplasm or bound to an immobile structure inside cell. We FRAP at this structure, to know the kinetic constants of the protein against the structure (e.g. microtubule binding protein, structure = microtubule) + Next we modify above model to consider a situation a bit more frequently we see in cell biology. The protein we are analyzing is either freely diffusing in cytoplasm or bound to an immobile structure inside cell. We FRAP this structure, to know the kinetic constants of the protein interaction with the structure (e.g. microtubule binding protein, structure = microtubule) - + $$\frac {df(t)} {dt} = k_{on}[free][s] - k_{off}[bound] \frac {df(t)} {dt} = k_{on}[free][s] - k_{off}[bound] - +$$ where where - * k_{on} Binding constant + * $k_{on}$ Binding constant - * k_{off} Dissociation constant + * $k_{off}$ Dissociation constant - * [free] Density of free molecules + * $[free]$ Density of free molecules - * [s] Density of immobile binding partner + * $[s]$ Density of immobile binding partner - * [bound] Density of bound-molecules + * $[bound]$ Density of bound-molecules - Since [s] is immobile and constant during experiment, we define k*_{on} as + Since [s] is immobile and constant during experiment, we define $k*_{on}$ as - + $$k*_{on}=k_{on}[s] k*_{on}=k_{on}[s] - +$$ - in addition, density of free molecule in cytoplasm is almost constant so we assume [free] = F and does not change. We then solve + in addition, density of free molecule in cytoplasm is almost constant so we assume $[free] = F$ and does not change. We then solve - + $$\frac {df(t)} {dt} = k*_{on}F - k_{off}[bound] \frac {df(t)} {dt} = k*_{on}F - k_{off}[bound] - +$$ We get We get - + $$f(t)=1-Ce^{- \tau t} f(t)=1-Ce^{- \tau t} - +$$ where where - * \tau = k_{off} + * $\tau = k_{off}$ - ... note that the recovery now only depends on k_{off} + ... note that the shape of recovery curve now only depends on $k_{off}$ ==== Diffusion and Reaction combined Recovery ==== ==== Diffusion and Reaction combined Recovery ==== [{{ :documents:reactiondiffusion.jpg?150| Diffusion-Reaction Combined Model}}] [{{ :documents:reactiondiffusion.jpg?150| Diffusion-Reaction Combined Model}}] - + $$\frac{\partial [free]}{\partial t} = D_{free} \nabla ^2[free]-k_{on}[free][s]+k_{off}[bound] \frac{\partial [free]}{\partial t} = D_{free} \nabla ^2[free]-k_{on}[free][s]+k_{off}[bound] - +$$ - + $$\frac{\partial [s]}{\partial t} = D_s \nabla ^2[s]-k_{on}[free][s]+k_{off}[bound] \frac{\partial [s]}{\partial t} = D_s \nabla ^2[s]-k_{on}[free][s]+k_{off}[bound] - +$$ - + $$\frac{\partial [bound]}{\partial t} = D_{bound} \nabla ^2[bound]+k_{on}[free][s]-k_{off}[bound] \frac{\partial [bound]}{\partial t} = D_{bound} \nabla ^2[bound]+k_{on}[free][s]-k_{off}[bound] - +$$ Since Since * [s] is constant and immobile * [s] is constant and immobile - * k*_{on} = k_{on}[s] + * $k*_{on} = k_{on}[s]$ - * \frac{\partial [s]}{\partial t}=0  + * $\frac{\partial [s]}{\partial t}=0$ - * bound molecules do not diffuse so D_{bound}=0 + * bound molecules do not diffuse so $D_{bound}=0$ Then we solve only Then we solve only - + $$\frac{\partial [free]}{\partial t} = D_{free} \nabla ^2[free]-k*_{on}[free]+k_{off}[bound] \frac{\partial [free]}{\partial t} = D_{free} \nabla ^2[free]-k*_{on}[free]+k_{off}[bound] - +$$ - + $$\frac{\partial [bound]}{\partial t} = k_{on}[free][s]-k_{off}[bound] \frac{\partial [bound]}{\partial t} = k_{on}[free][s]-k_{off}[bound] - +$$ We could solve this either analytically (Sprague et al, 2004) or numerically (Beaudouin et al, 2006). In the latter paper, calculation involves spatial context (on-rate was spatially varied; see also "geometry" section below). We could solve this either analytically (Sprague et al, 2004) or numerically (Beaudouin et al, 2006). In the latter paper, calculation involves spatial context (on-rate was spatially varied; see also "geometry" section below). Line 166: Line 173: Analytical solution was made in Laplace transformed equation. Analytical solution was made in Laplace transformed equation. + $$\overline{frap(p)} = \frac 1 p - \frac{F_{eq}}{p}\left(1-2K_1(qw)I_1(qw)\right)\times\left(1+\frac{k*_{on}}{p+k_{off}}\right)-\frac C {p+k_{off}} \overline{frap(p)} = \frac 1 p - \frac{F_{eq}}{p}\left(1-2K_1(qw)I_1(qw)\right)\times\left(1+\frac{k*_{on}}{p+k_{off}}\right)-\frac C {p+k_{off}} +$$ === Beaudouin Method === === Beaudouin Method === ==== Diffusion and Transport combined Recovery ==== ==== Diffusion and Transport combined Recovery ==== - We did not talk about this issue in the course, but there is another factor that could interfere with recovery curve in vivo: active transport. There is some trial on including issue by Hallen and Endow (2009). + We did not talk about this issue in the course, but there is another factor that could interfere with recovery curve in vivo: active transport. There is some trial on including this factor by Hallen and Endow (2009). - ==== Diffusion and Reaction, along with Spatial Context ==== + ==== Diffusion and Reaction, along with Spatial Context, Geometry ==== [{{ :documents:frapanalysis_sbalzarini.jpg?150| Example case considering geometrical and structural parameters that are affecting the recovery, by Sbalzarini et al}}] [{{ :documents:frapanalysis_sbalzarini.jpg?150| Example case considering geometrical and structural parameters that are affecting the recovery, by Sbalzarini et al}}] - ===== Tools for FRAP Analysis ===== + Since molecular behavior inside cell is constrained largely by structure and geometry of intracellular architecture, FRAP measured at single point within a cell does not always represent the biochemical characteristic of that molecule. In fact, functionality of protein molecule is determined not only by switching on/off, but is also regulated by the position of that molecule within cell. This means that spatial context should be included when interpreting the FRAP measurement. + Physical parameters such as Diffusion coefficient measured by FRAP is affected largely by geometrical constraint. Even if the geometry is rather simple, there are many obstacles in intracellular space which will cause longer time for molecules to reach from one point to the other. In such cases (which probably is frequently the case), estimated diffusion coefficient would calculated to be smaller than that of the "true diffusion coefficient". Presence of such obstacles should be somehow taken into account. For this reason, estimation of Diffusion coefficient could be more precise if one analyzes the structural geometry of where the molecule is constrained. Such protocol would be especially important if Diffusion and Reaction coupled recovery curve since wrong estimates on Diffusion coefficient would end up in wrong reaction parameters as well. + + Joel Beaudouin who did PhD study in the Ellenberg lab in the EMBL actually encountered such question in his project on nuclear protein study and solved the problem by using initial image frames of the FRAP experiment and let the molecule to diffuse by simulation, then fit the simulation with the experimental FRAP image sequence. Diffusion-reaction model was used and simulation was done (see above) using ODE solver to scan through the parameter space. Excellent idea. For more details, see Beaudouin et al. (2006). We will later add some protocol to this lecture notes on how to fit FRAP image sequence data using ODE solver (with help of Sebastian Huet). Joel's method was actually made into application called "Tropical" by a team in University of Heidleberg, but it seems that there is no direct link for downloading. + + Other papers we could refer to are Sbalzarini et. al. (2005, 2006). In these papers, authors did two things in parallel: 3D reconstruction of ER membrane structure and FRAP of certain molecule moving around along ER membrane. Using the 3D structure they reconstructed, geometrical constraint on protein diffusion could be determined and include this constraint on the estimation of Diffusion coefficient. + + ===== Pitfalls in FRAP Analysis ===== + + Refer to Mueller et. al.(2008). In their paper, they pointed out + * Shape of the FRAP ROI largely affect estimated value of biochemical rate constants(diffusion is important) + * Problems of fitting double exponential curve + * Initial condition (laser intensity profile) is important + * “blinding” of photomultiplier after the FRAP bleaching + + ===== List of Tools for FRAP Analysis ===== ==== Basic ==== ==== Basic ==== Line 188: Line 210: * [[http://cmci.embl.de/downloads/frap_analysis |FRAP analysis]] (EMBL) * [[http://cmci.embl.de/downloads/frap_analysis |FRAP analysis]] (EMBL) * Import data output from Zeiss, Leica, Olympus measurements and do FRAP fitting. * Import data output from Zeiss, Leica, Olympus measurements and do FRAP fitting. + * [[http://actinsim.uni.lu/eng/Downloads | FRAP analyzer]] (University of Luxemburg) + * Similar to above, but stand alone and also incorporated diffusion-reaction model. * [[http://mipav.cit.nih.gov/index.php | MIPAV]] * [[http://mipav.cit.nih.gov/index.php | MIPAV]] * Does measurement and fitting. Sprague et al. (2004) Reaction-Diffusion Full model is implemented. * Does measurement and fitting. Sprague et al. (2004) Reaction-Diffusion Full model is implemented. Line 194: Line 218: Requires your own coding, customization\\ Requires your own coding, customization\\ - === Analytical Approach === + ==== Analytical Approach ==== Sprague et. al. (explained above) is an example case of analytically solving the model for the fitting. Sprague et. al. (explained above) is an example case of analytically solving the model for the fitting. + ==== ODE Simulation ==== ==== ODE Simulation ==== [{{ :documents:tropical.jpg?150| Image based simulation with ODE solver}}] [{{ :documents:tropical.jpg?150| Image based simulation with ODE solver}}] Tropical Tropical - * Numerical analysis based on ODE. Spatial context. + * Numerical analysis based on ODE. Spatial context. **General Solvers** **General Solvers** - * Berkley Madonna + * [[http://www.berkeleymadonna.com/|Berkley Madonna]] * Joel and Sebastian uses this software for fitting ODE. * Joel and Sebastian uses this software for fitting ODE. - * MATLAB + * [[http://www.mathworks.com/|MATLAB]] ==== Particle Simulation ==== ==== Particle Simulation ==== Line 214: Line 239: * GridCell * GridCell * MCell * MCell + ===== References recommended ===== ===== References recommended ===== Line 237: Line 263: * full reaction-diffusion fitting by numerical approach, with spatial context * full reaction-diffusion fitting by numerical approach, with spatial context * [[http://www.ncbi.nlm.nih.gov/pubmed/16940327 | Ulrich M, Kappel C, Beaudouin J, Hezel S, Ulrich J, Eils R.(2006)]]Tropical--parameter estimation and simulation of reaction-diffusion models based on spatio-temporal microscopy images. Bioinformatics. 22(21):2709-10 * [[http://www.ncbi.nlm.nih.gov/pubmed/16940327 | Ulrich M, Kappel C, Beaudouin J, Hezel S, Ulrich J, Eils R.(2006)]]Tropical--parameter estimation and simulation of reaction-diffusion models based on spatio-temporal microscopy images. Bioinformatics. 22(21):2709-10 + * [[http://www.ncbi.nlm.nih.gov/pubmed/18199661|Mueller F, Wach P, McNally JG.(2008)]] Evidence for a common mode of transcription factor interaction with chromatin as revealed by improved quantitative fluorescence recovery after photobleaching. Biophys J. 2008 Apr 15;94(8):3323-39 + * Pitfalls on FRAP analysis. You will be shocked by how problematic it could be...