The exchange dynamics of biomolecular condensates

  1. Yaojun Zhang  Is a corresponding author
  2. Andrew GT Pyo
  3. Ross Kliegman
  4. Yoyo Jiang
  5. Clifford P Brangwynne
  6. Howard A Stone
  7. Ned S Wingreen  Is a corresponding author
  1. Center for the Physics of Biological Function, Princeton University, United States
  2. Department of Physics and Astronomy, Johns Hopkins University, United States
  3. Department of Biophysics, Johns Hopkins University, United States
  4. Department of Physics, Princeton University, United States
  5. Department of Chemical and Biological Engineering, Princeton University, United States
  6. Howard Hughes Medical Institute, United States
  7. Department of Mechanical and Aerospace Engineering, Princeton University, United States
  8. Department of Molecular Biology, Princeton University, United States
  9. Lewis-Sigler Institute for Integrative Genomics, United States
11 figures and 2 additional files

Figures

FRAP experiments on droplets.

(A) Schematic of a FRAP experiment in which an entire droplet is photobleached and the recovery of fluorescence is recorded. Material exchange between the condensate and the surrounding dilute phase can be limited by the flux of unbleached molecules coming from the dilute phase, the speed of internal mixing in the dense phase, or the flux passing through the interface. (B) The recovery time τ is defined as the time required for fluorescence to return 63% (i.e. 11/e) of the way back to its original level. (C) Experimental data from Taylor et al., 2019 in which a LAF-1 droplet of radius R=1µm recovers from photobleaching. (Left) Images before bleaching, immediately after bleaching of the entire droplet region, and at two subsequent times. (Right) Expected recovery times τdil=cdenR2/(3cdilDdil)=4.2±3.2s and τden=R2/(π2Dden)=60±18s if the slowest recovery process was either the flux from the dilute phase or diffusion within the droplet, respectively, with Dden=0.0017±0.0005µm2/s, Ddil=94±11µm2/s, and cden/cdil=1190±880 taken from Taylor et al., 2019. While the timescale associated with interface resistance τint is unknown, the measured recovery time τ=4570±470s is much longer than τdil and τden, suggesting the recovery is limited by flux through the interface, with an interface conductance of κ=R/3/(ττdilτden)=(7.4±0.8)×105µm/s.

Coarse-grained molecular-dynamics simulations of multivalent phase-separating polymers.

(A) Each polymer is composed of monomers (‘stickers’) of type A (blue) or B (red), and modeled as a linear chain of spherical particles each with a diameter of 2 nm, connected by stretchable bonds with an equilibrium length of 3.9 nm. Stickers of different types have an attractive interaction, while stickers of the same type interact repulsively, ensuring one-to-one binding between the A and B stickers. (B) Snapshot of a simulation of 1000 A6B6 polymers in a 500 nm × 50 nm × 50 nm box with periodic boundary conditions. The system undergoes phase separation into a dense phase (middle region) and a dilute phase (two sides), driven by the one-to-one A-B bonds. (C) Polymer concentration profile for the simulation in (B) with the center of the dense phase aligned at x=0 and averaged over time and over 10 simulation repeats. (D) Average total polymer concentrations in the dense (top) and dilute (bottom) phases from simulations of the five types of polymers shown in (A). (E) Polymer diffusion coefficients in the dense (top) and dilute (bottom) phases. All simulations were performed and snapshots were obtained using LAMMPS Plimpton, 1995. Please refer to Appendix 2 for simulation details.

Determination of interface conductance from simulations.

(A) Illustration of simulation protocol: At t=0 only polymers in the dilute phase are ‘labeled’ (solid balls), any polymer that enters the dense phase (forms an A-B bond lasting >10 times the average bond lifetime of an isolated A-B pair) becomes permanently ‘unlabeled’ (hollow balls). (B) Fraction of labeled polymers in the dilute phase as a function of time for simulations of the five types of polymers shown in Figure 2A. (C) Decay time of labeled polymers from exponential fits to curves in (B) (top), and corresponding calculated values of interface conductance κ (bottom). (D) For all simulated polymers, interface conductance scaled by κ0 (Equation 13) is approximately a linear function of a parameter u which reflects the fraction of unbound stickers in the dense and dilute phases. Inset illustration: Polymers in the dilute phase with few or no unbound stickers may ‘bounce’ off the dense phase, which contributes to the interface resistance. (E) Example of simulated trajectory in which a dilute-phase A6B6 polymer ‘bounces’ multiple times before finally joining the dense phase. (F) Interface conductance κ of A6B6 system as a function of binding strength U0 between A and B stickers.

Molecular-dynamics simulations of in silico ‘FRAP experiment’ on a small droplet of A6B6 polymers.

(A) Snapshot of 2000 A6B6 polymers in a 286 nm × 286 nm × 286 nm box with periodic boundary conditions. The system phase separates into a dense droplet (middle) and a surrounding dilute phase. (B) Polymer concentration profile for the simulation in (A) with the center of the droplet aligned to the origin and averaged over time and over 10 simulation repeats. The dilute- and dense-phase concentrations obtained from simulations in the slab geometry are denoted by green lines. The effective radius R of the droplet is determined through 4πR3cdend/3+(Vbox4πR3/3)cdild=Ntotal, where Vbox is the volume of the simulation box and Ntotal is the total number of polymers. (C) Average concentration profile of bleached population measured at different times. (D) Comparison of the FRAP recovery curves obtained by tracking the simulated fraction of unbleached molecules inside the droplet as a function of time (black dots) and by numerical integration (red curve) of Equation 1 in a sphere of the same volume as Vbox using the measured parameters of the droplet system: cdend=8.1mM, cdild=0.073mM, Dden=0.013μm2/s, Ddil=17μm2/s, R=0.037μm, and κd=0.20μm/s. For details of simulation and theory, refer to Appendix 2.

FRAP recovery patterns for large versus small droplets can be notably different for condensates with a sufficiently large interface resistance.

(A) Expected relaxation time as a function of droplet radius for in silico ‘FRAP experiments’ on the A6B6 system. The interface resistance dominates recovery times for smaller droplets, whereas dense-phase diffusion dominates recovery times for larger droplets. Green circle: FRAP recovery time obtained from direct simulation of an A6B6 droplet of radius 37 nm. Black curve: the recovery time as a function of droplet radius from a single exponential fit of the exact solution of the recovery curve from Equation 1–4. Gray curve: the recovery time predicted by Equation 6. Yellow, blue, and red curves: the recovery time when dense-phase, dilute-phase, and interface flux limit the exchange dynamics, i.e., the first, second, and last term in Equation 6, respectively. Parameters matched to the simulated A6B6 system in the slab geometry: cden=7.7mM, cdil=0.05mM, Dden=0.013μm2/s, Ddil=17μm2/s, and κ=0.14μm/s. (B) Time courses of fluorescence profiles for A6B6 droplets of radius 1μm (top) and 0.1μm (bottom); red is fully bleached, green is fully recovered. These concentration profiles are the numerical solutions of Equations 1–3 using parameters provided in (A).

Appendix 1—figure 1
Representative curves of g(u,t).

Parameters matched to the simulated A6B6 system: cden=7.7mM, cdil=0.05mM, Dden=0.013µm2/s, Ddil=17µm2/s, and κ=0.14µm/s. Droplet radius is R=1µm (top) and 0.1µm (bottom). Inset: same plot with the y-axis in log scale.

Appendix 1—figure 2
Comparison of exact and approximate solutions of Equation 24.

Inset: percentage error of the approximate solution.

Appendix 2—figure 1
Interface conductance inferred from average survival time of particles initially in the dilute phase for different criteria for having ‘joined’ the dense phase (Figure 3).

The interface conductance is shown for three criteria: a polymer is considered to have joined the dense phase if it is in the dense-phase cluster for a continuous duration of τ, 10τ, or 20τ, where τ is the average bond lifetime of an isolated A-B sticker pair.

Appendix 2—figure 2
Comparison of interface conductance κ obtained with the flux method and with Equation 54 as reported in Figure 3C.
Appendix 2—figure 3
Representative trajectories of molecules ‘bouncing’ multiple times at the interface of the A6B6 system (top) and A10B6 system (bottom).
Appendix 2—figure 4
A simulation box of side length 286 nm stabilizes the droplet size.

Number of polymers in the droplet as a function of time averaged over 10 simulation replicates for simulation boxes of side length 275 nm, 286 nm, and 300 nm.

Additional files

MDAR checklist
https://cdn.elifesciences.org/articles/91680/elife-91680-mdarchecklist1-v1.docx
Source code 1

MATLAB scripts used to generate input files for the LAMMPS Molecular Dynamics Simulator, as well as the corresponding LAMMPS input files.

https://cdn.elifesciences.org/articles/91680/elife-91680-code1-v1.zip

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. Yaojun Zhang
  2. Andrew GT Pyo
  3. Ross Kliegman
  4. Yoyo Jiang
  5. Clifford P Brangwynne
  6. Howard A Stone
  7. Ned S Wingreen
(2024)
The exchange dynamics of biomolecular condensates
eLife 12:RP91680.
https://doi.org/10.7554/eLife.91680.3