Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Fitting models

The software includes several models that describe the autocorrelation curves obtained by fluorescence correlation spectroscopy. It also allows the addition of user-defined models. Below, all default models are described.

One-component simple diffusion

This model describes the simple diffusion of fluorescent molecules. The model assumes there is only one type of fluorescent molecule, meaning the diffusing molecules’ sizes are pretty much the same. Additionally, the molecules exhibit no so-called triplet states. The model function is given by Equation (1).

G(τ)=1Np11+ττd111+τκ2τd1 G\left(\tau\right) = \frac{1}{N_p}\frac{1}{1+\frac{\tau}{\tau_{d1}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^{2}\tau_{d1}}}}

where τ\tau is a lag-time, NpN_p stands for the number of molecules inside the focal volume, τd1\tau_{d1} is a diffusion time - average time the fluorescent probe spends inside the focal volume. κ\kappa is the structure parameter defined as the ratio of semi-axes of the ellipsoidal focal volume.

One-component diffusion with triplet states

This model describes the simple diffusion of fluorescent molecules and also accounts for molecules in a triplet state. The model assumes that there is only one type (size) of fluorescent molecule. The model function is given by Equation (2).

G(τ)=1Np(1+T1T)exp(ττf)11+ττd111+τκ2τd1 G\left(\tau\right) = \frac{1}{N_p} \left(1+\frac{T}{1-T}\right)\exp\left(-\frac{\tau}{\tau_f}\right) \frac{1}{1+\frac{\tau}{\tau_{d1}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^{2}\tau_{d1}}}}

Where τ\tau is a lag-time, NpN_p stands for the number of molecules inside the focal volume, τd1\tau_{d1} is a diffusion time, and the average time the fluorescent probe spends inside the focal volume. κ\kappa is the structure parameter defined as the ratio of semi-axes of the ellipsoidal focal volume. TT stands for the fraction of molecules being in a triplet (dark) state, and τf\tau_f is the average lifetime for triplet states.

Two-component simple diffusion

This model describes the simple diffusion of two types of fluorescent molecules. The model neglects photo-physical processes, such as triplet states. The difference in size of molecules results in two distinct diffusion times τd1\tau_{d1} and τd2\tau_{d2}.

G(τ)=1Np((A111+ττd111+τκ2τd1)+((1A1)11+ττd211+τκ2τd2)) G\left(\tau\right) = \frac{1}{N_p}\left(\left(A_1\frac{1}{1+\frac{\tau}{\tau_{d1}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d1}}}}\right)+\left(\left(1-A_1\right)\frac{1}{1+\frac{\tau}{\tau_{d2}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d2}}}}\right)\right)

Where τ\tau is a lag-time, NpN_p stands for the number of molecules inside the focal volume, τd1\tau_{d1} is a diffusion time of component 1, τd2\tau_{d2} is a diffusion time of component 2, A1A_1 stands for the fraction of component 1, and κ\kappa is the structure parameter defined as the ratio of semi-axes of the ellipsoidal focal volume.

Two-component diffusion with triplet states

This model describes the diffusion of fluorescent molecules with two different diffusion coefficients (different sizes). The model also includes information about molecules in a triplet state. The model function is given by Equation (4).

G(τ)=1Np(1+T1T)exp(ττf)((A111+ττd111+τκ2τd1)+((1A1)11+ττd211+τκ2τd2)) G\left(\tau\right) = \frac{1}{N_p} \left(1+\frac{T}{1-T}\right)\exp\left(-\frac{\tau}{\tau_f}\right) \left(\left(A_1\frac{1}{1+\frac{\tau}{\tau_{d1}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d1}}}}\right)+\left(\left(1-A_1\right)\frac{1}{1+\frac{\tau}{\tau_{d2}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d2}}}}\right)\right)

Where τ\tau is a lag-time, NpN_p stands for the number of molecules inside the focal volume, τd1\tau_{d1} is a diffusion time of component 1, τd2\tau_{d2} is a diffusion time of component 2, A1A_1 stands for the fraction of component 1, and κ\kappa is the structure parameter defined as the ratio of semi-axes of the ellipsoidal focal volume. TT stands for the fraction of molecules being in a triplet (dark) state, and τf\tau_f is the average lifetime for triplet states.

Anomalous diffusion

This model describes the diffusion of fluorescent probes that cannot be described by simple diffusion. In this model, the diffusion of molecules through the focal volume is described by an algebraic distribution of diffusion times, namely ω2=4Dτdα\omega^2 = 4D\tau_d^\alpha, where α<1\alpha < 1. This model is also convenient to describe the motion of polydisperse molecules. i.e. Fluorescently labelled polymers characterised by a large polydispersity index, PDI. This model is recommended for the analysis of polymer diffusion in the cytoplasm of living cells Kalwarczyk & others (2017). The model is described by equation (5).

G(τ)=1Np11+(ττd1)α11+τακ2τd1α G\left(\tau\right) = \frac{1}{N_p}\frac{1}{1+\left(\frac{\tau}{\tau_{d1}}\right)^\alpha}\frac{1}{\sqrt{1+\frac{\tau^\alpha}{\kappa^{2}\tau_{d1}^\alpha}}}

In equation (5) τ\tau is a lag-time, NpN_p stands for the number of molecules inside the focal volume, τd1\tau_{d1} is a diffusion time - average time the fluorescent probe spends inside the focal volume. κ\kappa is the structure parameter defined as the ratio of semi-axes of the ellipsoidal focal volume, and α\alpha is an anomality exponent.

Three-component simple diffusion

This model describes the simple diffusion of three types of fluorescent molecules. The model neglects photo-physical processes, such as triplet states. The difference in size of molecules results in three distinct diffusion times τd1\tau_{d1}, τd2\tau_{d2}, and τd3\tau_{d3}.

G(τ)=1Np(G1+G2+G3)G1=(A111+ττd111+τκ2τd1)G2=((A2)11+ττd211+τκ2τd2)G3=((1A1A2)11+ττd311+τκ2τd3)\begin{align*} G\left(\tau\right) =& \frac{1}{N_p} \left(G_1+G_2+G_3\right)\\ G_1 =& \left(A_1\frac{1}{1+\frac{\tau}{\tau_{d1}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d1}}}}\right)\\ G_2 =& \left(\left(A_2\right)\frac{1}{1+\frac{\tau}{\tau_{d2}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d2}}}}\right)\\ G_3 =& \left(\left(1-A_1-A_2\right)\frac{1}{1+\frac{\tau}{\tau_{d3}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d3}}}}\right) \end{align*}

Here, τ\tau is the lag time, NpN_p is the number of molecules inside the focal volume, and A1A_1 and A2A_2 are the amplitudes of components 1 and 2.

Three-component simple diffusion with triplets

This model describes the simple diffusion of three types of fluorescent molecules differing in size. This model includes triplet states.

G(τ)=1Np(1+T1T)exp(ττf)(G1+G2+G3)G1=(A111+ττd111+τκ2τd1)G2=((A2)11+ττd211+τκ2τd2)G3=((1A1A2)11+ττd311+τκ2τd3)\begin{align*} G\left(\tau\right) =& \frac{1}{N_p} \left(1+\frac{T}{1-T}\right)\exp\left(-\frac{\tau}{\tau_f}\right) \left(G_1+G_2+G_3\right)\\ G_1 =& \left(A_1\frac{1}{1+\frac{\tau}{\tau_{d1}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d1}}}}\right)\\ G_2 =& \left(\left(A_2\right)\frac{1}{1+\frac{\tau}{\tau_{d2}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d2}}}}\right)\\ G_3 =& \left(\left(1-A_1-A_2\right)\frac{1}{1+\frac{\tau}{\tau_{d3}}}\frac{1}{\sqrt{1+\frac{\tau}{\kappa^2\tau_{d3}}}}\right) \end{align*}

Here, τ\tau is the lag time, NpN_p is the number of molecules inside the focal volume, τd1\tau_{d1}, τd2\tau_{d2}, and τd3\tau_{d3} are diffusion times of the first, second and third component, and A1A_1 and A2A_2 are the amplitudes of components 1 and 2. TT stands for the fraction of molecules being in a triplet (dark) state, and τf\tau_f is the average lifetime for triplet states.

Two-component diffusion: 1st translational diffusion only, 2nd Translational and rotational diffusion

The model describes diffusion in a system composed of two diffusing species, where one species undergoes rotations on a timescale comparable to that of translational diffusion across the focal volume. For details see Michalski et al. (2024)
Parameters description:
GinfG_\mathrm{inf} - offset
NpN_p - Number of fluorescent molecules
qq - amplitude of diffusing species
A2A_2 -Second, non-zero amplitude for rotation
A4A_4 - fourth, non-zero amplitude for rotation
ff - fraction relating the rotational to translational diffusion for a given diffusing species.
A=πkBA = \pi k_\mathrm{B}
TT - Absolute temperature
ηR1\eta_{R1} - Rotational viscosity experienced by the first diffusing species
ηT1\eta_{T1} - Translational viscosity experienced by the first diffusing species
ηT2\eta_{T2} - Translational viscosity experienced by the second diffusing species
r1r_1 - Hydrodynamic radius of the first diffusing species
r2r_2 - Hydrodynamic radius of the second diffusing species
ω0\omega_0 - Width of the focal volume
κ\kappa - Structure parameter

G(τ)=Ginf+1Np(qGR1GT1+(1q)GT2)GR1=1+A2fexp(3Tτ4AηR1r13)+A4fexp(5Tτ2AηR1r13)GT1=11+2Tτ3AηT1ω02r111+2Tτ3AηT1κ2ω02r1GT2=11+2Tτ3AηT2ω02r211+2Tτ3AηT2κ2ω02r2\begin{align*} G\left(\tau\right)=&G_{\mathrm{inf}}+\frac{1}{N_p}\left(qG_{R1}G_{T1}+\left(1-q\right)G_{T2}\right)\\ G_{R1} =& 1+A_2f \exp\left(-\frac{3T\tau}{4A\eta_{R1}r_1^3}\right)+A_4f\exp\left(-\frac{5T\tau}{2A\eta_{R1}r_1^3}\right)\\ G_{T1} =& \frac{1}{1+\frac{2T\tau}{3A\eta_{T1}\omega_0^2r_1}}\frac{1}{\sqrt{1+\frac{2T\tau}{3A\eta_{T1}\kappa^2\omega_0^2r_1}}}\\ G_{T2} =& \frac{1}{1+\frac{2T\tau}{3A\eta_{T2}\omega_0^2r_2}}\frac{1}{\sqrt{1+\frac{2T\tau}{3A\eta_{T2}\kappa^2\omega_0^2r_2}}} \end{align*}

One-component diffusion with rotation

The model describes diffusion in a system composed of a single diffusing species with coupled translational and rotational modes. The model also includes triplet states and an anomalous diffusion exponent, which can be set to 1 if normal diffusion is assumed. For details see Michalski et al. (2024)
Parameters description:
GinfG_\mathrm{inf} – offset (baseline of the correlation function for long lag times)
𝑁p𝑁_\mathrm{p} – number of fluorescent molecules in the observation volume
fTf_\mathrm{T} – fraction of molecules in the triplet state (triplet amplitude)
tTt_\mathrm{T} – triplet relaxation time
fR1f_\mathrm{R1} – first non-zero amplitude describing rotational diffusion contribution
fR2f_\mathrm{R2} – second non-zero amplitude describing rotational diffusion contribution
AA – constant defined as A=πkBA=\pi k_\mathrm{B}
TT – absolute temperature
ηR\eta_\mathrm{R} – rotational viscosity experienced by the diffusing species
ηT\eta_\mathrm{T} – translational viscosity experienced by the diffusing species
rr – hydrodynamic radius of the diffusing species
ω0\omega_0 – lateral width (radius) of the focal volume
κ\kappa – structure parameter (ratio of axial to lateral dimensions of the focal volume)
α\alpha – anomalous diffusion exponent describing deviation from normal diffusion

G(τ)=Ginf+1NpGT(τ)Gr(τ)GD(τ)GT(τ)=(1+fTeτ/tT)Gr(τ)=(1+fR1eτ4A(ηR/T)r3+fR2e3τ4A(ηR/T)r3)GD(τ)=(1+(2τ3A(ηT/T)w02r)α)1(1+(2τ)α(3A(ηT/T)w02)ακ2rα)1/2\begin{align*} G(\tau)=&G_{\mathrm{inf}}+\frac{1}{N_p} G_\mathrm{T}(\tau) G_\mathrm{r}(\tau) G_\mathrm{D}(\tau) \\ G_\mathrm{T}(\tau)=&\left(1+f_T e^{-\tau/t_T}\right)\\ G_\mathrm{r}(\tau)=&\left( 1+f_{R1}e^{-\frac{\tau}{4A(\eta_R/T)r^3}} +f_{R2}e^{-\frac{3\tau}{4A(\eta_R/T)r^3}} \right)\\ G_\mathrm{D}(\tau)=&\left( 1+\left(\frac{2\tau}{3A(\eta_T/T)w_0^2 r}\right)^{\alpha} \right)^{-1} \left( 1+\frac{(2\tau)^{\alpha}} {(3A(\eta_T/T)w_0^2)^{\alpha}\kappa^2 r^{\alpha}} \right)^{-1/2} \end{align*}
References
  1. Kalwarczyk, T., & others. (2017). Apparent anomalous diffusion in the cytoplasm of human cells: the effect of probes’ polydispersity. Journal of Physical Chemistry B, 121, 9831–9837. 10.1021/acs.jpcb.7b07158
  2. Michalski, J., Kalwarczyk, T., Kwapiszewska, K., Enderlein, J., Poniewierski, A., Karpińska, A., Kucharska, K., & Hołyst, R. (2024). Rotational and translational diffusion of biomolecules in complex liquids and HeLa cells. Soft Matter, 20(29), 5810–5821. 10.1039/D4SM00422A