Nonlinear Diffusion and Compression Rim in Biopolymer Condensates with Sources ✛ ✦

Original: "Nonlinear diffusion and compressive rims in source-driven biopolymer condensates"
· Avraham Moriel, Howard A. Stone
arXiv:2606.31786 · 2026-06-30 · CC BY 4.0 · 3 min
Scientists uncovered the mechanism behind the compressed layer at the boundary of growing biopolymer droplets and applied it to nucleolus dynamics.

Context

Many membrane-less subcellular structures function as microscale factories, continuously producing biopolymers. The nucleolus is a prime example: it synthesizes ribosomal RNA and exports it into the nucleoplasm. Recent experiments have shown that at the periphery of the nucleolus, RNA molecules are in a compressed state, and the organelle itself behaves like a viscoelastic gel. Understanding how entropic forces, flows, and mechanical stresses give rise to such architecture remained an open question. This work fills the gap by combining polymer physics with nonequilibrium thermodynamics of active cellular condensates.

Methods

The authors built a theoretical framework based on a two-fluid model, in which long biopolymers (A-mers) diffuse through a medium of short ones (B-mers) under osmotic pressure given by the Flory-Huggins potential. Viscoelastic stresses were described by the Maxwell model in the convective derivative—a cornerstone of rheology laid down by Maxwell. The equations were solved using numerical simulation in one and two dimensions, varying chain length, entanglement degree, and source time dependence. Similarity analysis was used to extract scaling laws, and experimental rRNA distribution profiles in nucleoli were digitized and rescaled.

Results

It was found that regardless of details, a compression rim always forms at the diffusion front of the polymer droplet—a zone where chain conformations are compressed (A−1 < 0). In the one-dimensional case with a constant source, the droplet expands according to ℓ ∼ $t^{3/4}$, indicating a nonlinear diffusion regime. Scaling analysis yields a universal solution φ(x,t) = (Q² t)^{1/4} Φ(ξ) with ξ = x/(Q² t³)^{1/4}, where Q is the source amplitude. Applying this to nucleolus data allowed the radial rRNA distribution curves to collapse with exponents δ₁ ≈ –0.05 and δ₂ ≈ 0.45. This corresponds to a source time-dependence exponent q ≈ 1.3, intermediate between constant flux (q=1) and the regime of maintaining a constant boundary concentration (q=1.5).

Implications

This work demonstrates for the first time that polymer compression at the periphery of active condensates is a direct consequence of osmotically-induced flows, not just specific interactions. This provides a mechanical basis for the previously observed 'thrust' of droplets against their surroundings. The model links the entropic nature of mixing to macroscopic deformations, opening the way to quantitative estimation of biopolymer synthesis rates from imaging data.

Future development

The proposed approach can be extended by including enthalpic, electrostatic, and specific intermolecular interactions, bringing it closer to atomistic descriptions. It is also promising to apply the model to other active biomolecular condensates—from stress granules to transcription clusters—and to test the prediction of the compression rim using super-resolution microscopy methods sensitive to local anisotropy.

Impact

The results are significant for biosignatures of cellular activity, since rRNA synthesis rates can serve as an indicator of the cell's functional state, including cancerous transformation. Moreover, they deepen understanding of nonequilibrium soft matter physics and open new engineering principles for creating synthetic self-sustaining microreactors.

Next steps

Next steps include verification of the compression rim in the nucleolus using polarization microscopy or FRET sensors, as well as adapting the model to three-dimensional geometries and more realistic rheologies.

Key open problems

The study touches on a fundamental problem of nonequilibrium self-organization in active matter: how a continuous influx of energy and mass shapes stable spatial structures in living systems. It also resonates with the question of the role of mechanical stresses in regulating cellular processes, which remains one of the major unsolved problems in biophysics.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The authors developed a two-fluid model with Flory–Huggins thermodynamics and Maxwell viscoelasticity to describe biopolymer condensates with sources. They showed numerically and analytically the formation of a compression rim at the diffusion front and derived scaling laws, which they compared with data on the nucleolus.

Where its strength is

A universal scaling analysis links observed density profiles to the synthesis rate, providing a quantitative tool for interpreting experiments.

Where it could go next

Worth trying

  • Test how periodic modulation of the source affects the width of the compression rim and the droplet growth rate.
  • Assess the contribution of viscoelastic stresses to the collective motion of multiple droplets interacting through a common polymer field.
  • Use the scaling relations to estimate synthesis rates in other condensates from published concentration profiles.

Original abstract

as written by the authors arXiv:2606.31786

Many subcellular condensates continuously produce biopolymers. Coupling Flory-Huggins thermodynamics to two-fluid viscoelasticity, we probe the diffusion of such source-driven polymeric droplets, and identify a universal structural compressive rim at their diffusion front. Integrating analytical scaling laws, numerical simulations, and experimental data, we show that this framework captures key structural and dynamic characteristics of the nucleolus, demonstrating the role of polymer diffusion in non-equilibrium biological transport.

The nucleolus is the largest membrane-less organelle in the nucleus, and despite lacking a shell, it manages to compress synthesized RNA molecules at its boundary, like a miniature piston powered by entropy.

f(\phi) = \frac{\phi}{N_A} \log \phi + \frac{1-\phi}{N_B} \log(1-\phi) + \chi \phi (1-\phi)
Defines the thermodynamic contribution to the osmotic pressure of a polymer mixture
\partial_t \phi \simeq \nabla \cdot \left[ \frac{(1-\phi)^2 \phi}{n_0 \zeta_0 \mu(\phi)} \nabla \Pi \right]
Describes the evolution of the volume fraction of A-mers with concentration-dependent mobility

Key numbers

  • scaling exponent of source time dependence q: ≈1.3
  • amplitude rescaling exponent δ₁: ≈–0.05
  • coordinate rescaling exponent δ₂: ≈0.45
  • entanglement parameter α: up to 100
  • polymer chain length N_A: up to 104
Scientists
Ludwig BoltzmannRudolf ClausiusJacob BekensteinStephen HawkingJames Clerk MaxwellEdward Witten
Concepts
diffusion microorganism swimming dynamics biofilm diffusion mri polymer fluid dynamics in porous media fluid flow microswimmer fluctuation-dissipation theorem surface tension flux constant Fick's law Gauss's theorem rainfall microplastic dispersion electrolyte transport