Physical observables to determine the nature of membrane-less cellular sub-compartments

  1. Mathias L Heltberg
  2. Judith Miné-Hattab
  3. Angela Taddei
  4. Aleksandra M Walczak  Is a corresponding author
  5. Thierry Mora  Is a corresponding author
  1. Laboratoire de physique de l’École normale supérieure, CNRS, PSL University, Sorbonne Université, and Université de Paris, France
  2. Institut Curie, CNRS, PSL University, Sorbonne Université, France
4 figures, 1 table and 1 additional file

Figures

Schematic setup of the models.

(A) In the middle, the observed signal from a fluorescently tagged Rad52 protein inside the nucleus following a double-stand break. Left: Schematic figure showing the Polymer Bridging Model (PBM). Proteins binding specifically to the chromatin stabilize it, effectively trapping the motion of other molecules. Right: Schematic figure showing the Liquid Phase Model (LPM). Liquid-liquid phase separation results in the formation of a droplet foci with a different potential and different effective diffusion properties than outside the droplet. (B) Details of the PBM model. Particles diffuse freely with diffusivity Dn until they hit one of the N spherical binding sites, themselves diffusing with diffusivity Db. The focus is formed due a high concentration of binding sites. The binding sites are only partially absorbing, so that not all collision events result in a binding even. Once bound, the particle stays attached to the binding site, and then unbinds with rate k-.

Figure 2 with 1 supplement
Diffusion properties and effective free energy.

(A) Example of experimental tracking of Rad52 molecules visiting a double-strand break (DSB) focus. Different connected traces correspond to distinct Rad52 molecules. (B) Example trajectory of a particle visiting the focus from simulations in the PBM (left) and LPM (right). The simulated trajectories are visually similar to the data in B. (C) Displacement histogram (jump sizes) for the PBM, LPM and experiments, for an interval δ=20 ms. D. Displacement histogram for the PBM for small values of k-dt (top) and high values (bottom). Here we varied the interval from δt=1 ms (top) to δt=15 ms (bottom). (E) Hypothesis testing using a two sided KS-test, comparing the displacement histogram of a free diffusion process (black line in D) and the displacement histogram of diffusion inside the focus (green line in D). Parameters are the same as in (D) δt was varied from 1 to 25 ms. (F) Effective diffusion coefficient as a function of distance to the focus center r, estimated from simulations of the PBM calculating D~=δr2/(2dδt) in each radial segment. (G) Particle density p(r) as a function of r, estimated from simulations of the PBM. Error bars are standard errors on the mean. (H) Relation between the ratio (D0Db)/(DnDb) versus the ratio of densities inside and outside the focus (From Equation 11), both estimated from simulations of the PBM (green crosses), compared to the identity prediction (Equation 12, black line). Blue cross shows the experimental observation for Rad52 in DSB loci (see text related to Figure 7 and Table S1 in Miné-Hattab et al., 2021). Parameter values as in Table 1 except: rn=1μm for B; A=2.5kBT for B-C; rn=0.3μm, rf=0.15μm and Dn=0.5μm2/s for D-E, κ=300μm/s for D, rn=0.75μm for F-H. In H we varied κ=1400μm/s, k-=11,500 s-1, and ρ=2.44.8104μm. See Figure 2—source data 1.

Figure 2—source data 1

Compressed ZIP file containing all the data plotted in the panels of Figure 2 as CSV and TXT files.

https://cdn.elifesciences.org/articles/69181/elife-69181-fig2-data1-v2.zip
Figure 2—figure supplement 1
Effect of crowding.

Same as 2 H., but with added red points corresponding to replacing 20 % (leftmost), 40%, 60%, and 80 % (rightmost) of the binding sites by inert, totally reflecting spheres to simulate crowding. While crowding should decrease the measured D0, in practice the effect is negligible.

Figure 3 with 1 supplement
Radial and angular dynamics.

(A) Experimental trace of a single Rad52 in a DSB focus (Miné-Hattab et al., 2021). The inset shows the definition of the radial movement δr. Here, concentric circles are shown to define the radius relative to the focus center, where a particle moves a specific distance away from the center of the focus, as well as the angle θ between consecutive displacements. (B) Data extracted from experiments (Miné-Hattab et al., 2021) to estimate the average radial displacement of the tracked particle multiplied by the radius. Here and throughout the x-axis represents the radial position at the beginning of the timestep. Error bars are standard errors on the mean. (C) Simulations showing the radial displacements in the PBM with slowly (top) and rapidly (bottom) moving binding sites. Black lines are predictions based on the measurement of D~ (Equation 12). Error bars are standard deviations on the mean. (D) Radial displacement from simulations of the LPM. Light-green line shows the (wrong) prediction made while assuming the PBM, using the measurement of the effective diffusion coefficient D~(r) (Equation 13). We call the discrepancy between data and the PBM prediction the ‘‘Maximal positive difference.’’ Error bars are standard errors on the mean. (E) Heatmap showing the maximal positive difference in LPM simulations as a function of D0 and A. (F) Distribution of angles (represented radially on the left, and linearly on the right) between the displacements of consecutive steps of length δt, from experiments and simulations. Multiple curves for the PBM correspond to different parameter choices corresponding to the points of Figure 2H. Parameter values as in Table 1, except rn=1μ m for F. For F parameters are varied with κ=100300μ m/s, k-=5001,500 s-1, rf=0.10.14μ m. See Figure 3—source data 1 and Figure 3—source data 2.

Figure 3—source data 1

Compressed ZIP file containing all the data plotted in the panels of Figure 3 as CSV files.

https://cdn.elifesciences.org/articles/69181/elife-69181-fig3-data1-v2.zip
Figure 3—source data 2

Compressed ZIP file containing all the data plotted in Figure 3—figure supplement 1 as CSV files.

https://cdn.elifesciences.org/articles/69181/elife-69181-fig3-data2-v2.zip
Figure 3—figure supplement 1
Asymmetry coefficient for an infinite PBM focus.

Asymmetry coefficient (following the definition of Izeddin et al., 2014), for an infinite focus, simulated with the PBM, as a function of the time step δt. Asymmetry is due to molecules reflecting off binding sites, causing more consecutive displacements to have 180 degree angles. The faster the binding and unbinding relative to the time step δt, the closer to mean-field limit of standard diffusion, and the more symmetric the angle distribution. k+ is changed alongside k- to keep pu constant. The asymmetry coefficient is defined as log2[P(|θ|<π/6)/P(|θ|>5π/6)], where π<θ<π. Standard deviation is obtained from four independent runs of 1000 s.

First passage times to a target site inside the focus.

A. Schematic figure showing the setup of the tracked molecule and the effective target. B. Time to reach the specific target in simulations for the LPM. Black curve shows the predicted result from the analytical derivation (Equation 15). Here we use parameters: A=-5.5kBT, D0=0.05μ2/s, rn=1.0μm, Dn=0.8μ2/s. C. Same as B, but for the PBM. Same parameters as in B, but with D0=DneU(r)/kBT. D. Heatmap showing the expected search time as a function of the droplet size (x-axis) and the focus diffusion coefficient (y-axis) (Equation 15). Green point corresponds to experimental observations. E. Heatmap showing the expected search time as a function of the droplet size (x-axis) and the height of the surface potential (y-axis) (Equation 15). Green point corresponds to experimental observations. Parameter values as in Table 1, except b=2,000μm-1, D0=0.04μm2/s, A=5.5kBT for B, E, and F, κ=50μm2/s, k-=100 s-1, ρ=2.4104μm-3 for C-D. See Figure 4—source data 1.

Figure 4—source data 1

Compressed ZIP file containing all the data plotted in the panels of Figure 4 as CSV files.

https://cdn.elifesciences.org/articles/69181/elife-69181-fig4-data1-v2.zip

Tables

Table 1
Parameters used in this study with their typical values, and the ranges we have considered.

Experimental values are from Miné-Hattab et al., 2021 (see Materials and Methods for details on estimating diffusion coefficients and free energy differences). D0 and A are model parameters in the LPM, but also effective observables in the PBM. The diffusivity of binding sites is taken to be that of Rfa1 molecules in the focus, which bind to single-stranded DNA in repair foci, and are thus believe to follow the diffusion of the chromatin (Miné-Hattab et al., 2021). The number N of binding site is related to their density ρ inside the focus through N=(4/3)πρrf3.

VariableModelDescriptionValueRangeExp. valueUnits
rfbothradius of focus10050–200nm
rnbothradius of nucleus500300–1000nm
DnbothDiffusion coefficient in nucleus1.00.5–2.01.08μm2/s
σbothExperimental noise level303030nm
D0LPMDiffusion coefficient inside droplet0.050.01–0.50.032μm2/s
ALPMSurface potential5.00–105.5kBT
bLPMSteepness in potential1,000500–10000μm-1
ρPBMDensity of binding sites inside focus4.81041103-8.4104μm-3
DbPBMDiffusion coefficient of binding sites0.0050-0.10.005μm2s-1
rbPBMRadius of binding sites105-20nm
k-PBMUnbinding rate50010-10,000s-1
κPBMAbsorption parameter1000-1000μm/s

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Mathias L Heltberg
  2. Judith Miné-Hattab
  3. Angela Taddei
  4. Aleksandra M Walczak
  5. Thierry Mora
(2021)
Physical observables to determine the nature of membrane-less cellular sub-compartments
eLife 10:e69181.
https://doi.org/10.7554/eLife.69181