Quantitative theory for the diffusive dynamics of liquid condensates
Figures
![](https://iiif.elifesciences.org/lax/68620%2Felife-68620-fig1-v2.tif/full/617,/0/default.jpg)
Quantitatively measuring by extracting intensity at condensate interface.
(a) (1) Spatial time course of FRAP recovery after full bleach for PLYS/ATP coacervate droplet. Azimuthal average (orange) is highlighted. (2) To obtain the boundary condition , fluorescence intensity is extracted approximately away from the interface (see methods). Data closer to the droplet boundary can be subject to optical artefacts giving rise to an artificially broad interface. Lower panel shows azimuthal averages at different time points. (3) Cost function for fitting. For each fitting step, a numerical solution to Equation (1a) is calculated with a trial and the experimentally measured boundary condition from step (2). Subsequently simulation and experiment are compared according to the cost function, before choosing the next trial . Scale bar, . (b) Flurorescence recovery inside a PGL-3 condensate, along the radial direction (azimuthal averages, see (a)). is extracted by global fitting of Equation (1b) to the experimental profiles, using the experimentally extracted initial and boundary conditions (see panel (a) (1) and (2)). The gray line indicates an offset that comes about due to incomplete bleaching, and a small but visible, fast and uniform recovery with unknown cause (see Materials and methods). (c) Same as (b) but for ATP/PLYS coacervate droplet. (d) Same as (b) but for a CMD/PLYS coacervate droplet. (e) Cost function (see panel (a)) used to extract for a PGL-3 droplet. Due to the smoothness of the cost function, the minimum can be determined with high accuracy. (f) Comparison between and viscosity η for PGL-3 condensates at different salt concentrations. Viscosity data taken from Jawerth et al., 2018 for untagged PGL-3. Note, GFP-tagged PGL-3 has a higher viscosity than untagged PGL-3, which means the hydrodynamic radius of PGL-3:GFP cannot be directly computed from this panel (Jawerth et al., 2020). Replicates: salt concentration [mM]/#: 50/2, 60/3, 75/3, 90/3, 100/6, 120/4, 150/3, 180/3. (g) Diffusion coefficients for coacervate systems. Each data point represents for a single droplet time course. Note the low spread of the measured values. Replicates: CMD/PLYS: 9, PLYS/ATP:16.
![](https://iiif.elifesciences.org/lax/68620%2Felife-68620-fig2-v2.tif/full/617,/0/default.jpg)
Ternary mixture accounts for the dynamics of bleached and unbleached molecules.
(a) Before bleaching, a droplet that is composed of fluorescently labeled molecules can be described by a binary mixture, namely unbleached molecules and solvent. After bleaching, the system is composed of three components, bleached molecules, unbleached molecules and solvent. If the system was at equilibrium prior to bleaching, the sum of bleached and unbleached molecules forms a stationary, non-uniform profile (see panel b). (b) Snapshot of model dynamics at . Initial conditions are , corresponding to a fully bleached droplet. Note that at any time we have . (c) Time course of spatial recovery. For long times, when nearly all bleached material has been exchanged, approaches . Panels (c,d) use radial symmetry for illustration purposes, however, the theory is general (see Figure 3).
![](https://iiif.elifesciences.org/lax/68620%2Felife-68620-fig3-v2.tif/full/617,/0/default.jpg)
Impact of droplet environment on recovery dynamics.
(a) Sketch of a typical experimental set-up with a droplet above a passivated coverslip, where the droplet center has a distance to the coverslip. (b) Recovery of average unbleached volume fraction for different heights above the coverslip at different partition coefficients . Results were obtained by solving Equation (6) using the finite element method and considering the geometries depicted in (a). For even larger -values (e.g. no coverslip), results are approximately equal to the blue dashed line. (c) Using the method introduced in Figure 1 on the scenario with the largest influence of the coverslip (droplet sessile on coverslip) in (a) results in an excellent fit and can reliably extract the input . (d) Sketch of neighboring droplets next to a bleached droplet. (e) Total recovery curves for finite element simulations of the geometry depicted in (d), for different distances between neighboring droplet centers, at different partition coefficients . Note the strong dependence on the distance of neighboring droplets. For even larger -values (e.g. no neighboring drops), results are approximately equal to the blue dashed line. (f) Same as (c) but for largest influence of neighboring droplets, that is , where there is no distance between bleached droplet and neighboring droplets.
![](https://iiif.elifesciences.org/lax/68620%2Felife-68620-fig4-v2.tif/full/617,/0/default.jpg)
Varying partition coefficient and diffusivity outside simultaneously can lead to similar recovery kinetics.
(a) Schematics of the fitting approach used to extract the relationship between and . Although there is one specific pair (, ) with the lowest cost (see Figure 5), there is a line in (, ) space, which produces good fits for a range of values. (b) Spatial recovery of a single PLYS/ATP droplet (blue) with a fit to the full model (black). Note, data close to the droplet boundary cannot be fit, due to optical artefacts giving rise to an artificially broad interface (see Materials and methods and Figure 1a). (c) Same as (b) for a CMD/PLYS coacervate droplet. (d) Same as (b) for a PGL-3 droplet. (e) Given the partition coefficient , is found by fitting the coacervate data to the model. Note the convergence to a power law, with for large partition coefficients (gray dashed line; for a discussion see Appendix 2). Shaded area around curves: standard deviation. Replicates: same as in Figure 1g. (f) Same as (e) but for PGL-3 with different salt concentrations. Note the order from top to bottom from highest to lowest salt concentration. (g) Based on (f), the change of partition coefficient can be estimated for a given . Confidence intervals not shown for clarity. Similar to panel e. Replicates: same as in Figure 1f.
![](https://iiif.elifesciences.org/lax/68620%2Felife-68620-fig5-v2.tif/full/617,/0/default.jpg)
Given the recovery dynamics inside a condensate, the key parameters , and can in theory be determined uniquely without measuring the outside dynamics, using Equation (6).
(a) Given the partition coefficient , is found by fitting the model to synthetically generated example data. As example systems, indicated by open circles, we consider in silico data, obtained by solving Equation (6) with known parameters and . To mimic the approach of initially determining (see Figure 1) we keep for all in silico datasets. For each example system (open circles), the best fitting is found given a range of (solid lines). Parameters used were as follows: pink circle: , , orange circle: , , green circle: , , blue circle: , . (b) Cost function (colorbar, log-scale) as a function of and . We note that the global minimum coincides with the parameters used to generate the synthetic data (green circle). The valley in parameter space (dashed line) corresponds to the green line in (a). (c) Minimum of cost function for each , corresponding to curves shown in (a). This minimum corresponds to the valley indicated by the dashed line in (b). Note the minimum at the input parameter set, which indicates uniqueness of the outside dynamics for given values of and .
![](https://iiif.elifesciences.org/lax/68620%2Felife-68620-app2-fig1-v2.tif/full/617,/0/default.jpg)
Solution to dynamics inside droplets is unique, but early and late time behavior can be approximated by different scaling regimes.
(Top) Early time behavior of numerical solutions. If a parameter set is chosen (’Reference’), at early times this can be approximately recovered by scaling (orange). A linear scaling does not recover the same inside dynamics at early times (green). (Bottom) Late time behavior of numerical solutions. At late times, the reference dynamics can be approximately recovered by scaling (green). Meanwhile, does not recover the reference inside dynamics (orange).
Videos
FRAP dynamics in a PGL-3 droplet Left: Representative in vitro droplet after full bleach.
Time course starts after a small time lag due to a fast uniform recovery (see methods). Middle: Diffusion Equation (1) fit to azimuthally averaged droplet intensity with a global fit parameter . Right: Full model Equation (6) fit to azimuthally averaged droplet intensity.
FRAP dynamics in a CMD/PLYS coacervate.
For description see Video 1.
FRAP dynamics in a PLYS/ATP coacervate.
For description see Video 1.
Additional files
-
Source data 1
Quantification_Fig1_Fig4.
- https://cdn.elifesciences.org/articles/68620/elife-68620-data1-v2.zip
-
Transparent reporting form
- https://cdn.elifesciences.org/articles/68620/elife-68620-transrepform-v2.docx