🧬 biology

Entropy-Driven Initiation and Cellular Uptake Mediated by Viscoelastic Cytoskeleton: A Kinetic Phase Diagram from Onsager Variational Principle

This paper proposes a unified continuum model based on the Onsager variational principle to explain receptor-mediated endocytosis, demonstrating that entropic forces from crowded biomolecules initiate ligand-receptor contact while cytoskeletal viscoelasticity and binding energies govern a kinetic phase diagram that predicts an optimal virus size matching HIV-1 dimensions.

Original authors: Jinjie Liu, Zhong-Can Ou-Yang, Hao Wu

Published 2026-07-15

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Jinjie Liu, Zhong-Can Ou-Yang, Hao Wu

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ⚕️ This is an AI-generated explanation of a preprint that has not been peer-reviewed. It is not medical advice. Do not make health decisions based on this content. Read full disclaimer

Imagine a virus trying to break into a cell. For decades, scientists have been puzzled by a specific "chicken-and-egg" mystery: How does the virus even get close enough to the cell's surface to grab onto it? The old story was that the virus just waits there, hoping a receptor (a tiny hook on the cell) swims over to it. But this paper suggests that story is missing a crucial first step.

The authors, using a sophisticated mathematical toolkit called the Onsager variational principle, propose a new, energetic way to start the process. They suggest that the virus doesn't need to wait for a hook; instead, it gets a massive shove from the crowded environment inside the cell.

The "Mosh Pit" Push

Think of the space around a cell as a super-crowded dance floor, packed with tiny, bouncy molecules like globular proteins. The virus is a large, smooth ball trying to get near the cell membrane (the dance floor's edge).

In this crowded room, the tiny molecules can't squeeze into the tiny gap between the virus and the membrane. This creates a "depletion zone"—a no-go area for the small molecules. When the virus gets close to the membrane, these two no-go zones merge, freeing up a huge amount of space for the tiny molecules to bounce around in. Nature loves freedom; this increase in "wiggle room" (entropy) creates a powerful, invisible push that shoves the virus right against the cell wall.

The paper argues that this entropic force is the real "initiation" mechanism. It's the reason the virus gets close enough to start the handshake. Without this push, the virus might just float away before it ever touches the cell.

The Two-Step Dance

Once the virus is shoved close enough, the process splits into two distinct phases:

The Approach (Phase 1): The virus is driven by that "crowd push" through the fluid. The authors calculate that this approach happens incredibly fast—about 7.14 × 10⁻⁵ seconds. This is just the right amount of time: fast enough to get the virus there, but slow enough that the virus doesn't just bounce off before the cell's hooks (receptors) can grab on. It's a perfect timing match for the receptors to lock in.
The Wrap-Up (Phase 2): Once the hooks grab the virus, the cell starts to wrap its membrane around it like a blanket. But here's the twist: the cell isn't just a soft, floppy bag. Underneath the membrane is a cytoskeleton, a meshwork that acts like a viscoelastic sponge. It's stiff but also stretchy and slow to react.

The paper suggests that this "sponge" resists the virus. The cell has to slowly creep and deform to swallow the intruder. The speed of this swallowing depends on how stiff the cell is. If the cell is stiff (like a firm gel), the process takes longer. If it's soft, it's faster.

The "Goldilocks" Size

One of the most striking findings is about size. The paper predicts that there is a "sweet spot" for how big a virus should be to get eaten most efficiently.

Too small? The cell membrane is too hard to bend around a tiny speck. The energy cost is too high.
Too big? The cell's internal sponge (cytoskeleton) fights back too hard. It's like trying to swallow a beach ball; the resistance is too great.

The authors found that for a typical virus with standard binding strength, the perfect size is around 50 nm (nanometers). Interestingly, this matches the size of the HIV-1 virus almost perfectly. The paper suggests this isn't a coincidence; the virus might have evolved to be this size because it's the most energy-efficient way to get inside a cell.

If the virus has stronger "hooks" (higher binding energy), the sweet spot actually shifts to be even smaller, around 30–37 nm. This explains why different viruses and nanoparticles have different ideal sizes for entering cells.

The Rules of the Game

The paper uses math to draw a "phase diagram"—a map showing exactly when a virus can get in and when it can't.

The Stiffness Limit: If the cell is too stiff (specifically, if the Young's modulus is too high, around 8.5 kPa for the parameters used), the virus can't get in, no matter how big or small it is. The cell simply refuses to deform.
The Hook Limit: If the virus doesn't have enough hooks (ligand density), it can't get in either. There is a minimum number of hooks required to overcome the cell's resistance.

What This Paper Rules Out

The authors are very clear about what their model doesn't rely on. They explicitly argue against the idea that the process starts with the virus already touching the cell and waiting for receptors to swim over. They say that "diffusion-centric" views (where the receptor does all the moving) miss the fundamental problem of how the virus gets close in the first place. Their model says the crowd push happens before any specific binding occurs.

They also note that while their model fits the data beautifully, it treats the virus as a perfectly rigid ball. In reality, some viruses might squish or change shape, which could change the details. They also assume the hooks are spread out evenly, whereas in real life, they might clump together.

How Sure Are They?

The authors are quite confident in their theoretical framework. They didn't just guess; they built a unified mathematical model that connects the "crowd push" to the "sponge resistance."

The Math: They proved that their complex equations for a round virus naturally simplify to a famous, textbook result (the Asakura-Oosawa result) when the virus is huge and flat. This acts as a strong sanity check, suggesting their math is solid.
The Match: Their prediction of a 50 nm optimal size matches the real-world size of HIV-1 and other viruses, which gives them strong reason to believe their theory is on the right track.
The Caveat: While the numbers line up with experiments, this is a theoretical model. The authors suggest that future experiments could test these predictions by changing the stiffness of cells or the density of hooks to see if the "size window" for entry changes exactly as their map predicts.

In short, this paper suggests that viruses don't just luck into a cell; they get a nudge from the crowded molecular world, and then they have to dance perfectly with the cell's internal springs to get inside. And if they are the right size—around 50 nm—they win the dance.

Technical Summary: Entropy-Driven Initiation and Cellular Uptake Mediated by Viscoelastic Cytoskeleton

Problem Statement
A fundamental gap in existing theoretical models of receptor-mediated endocytosis is the "initiation problem": previous frameworks typically assume pre-existing contact between a virus (or nanoparticle) and the cell membrane, overlooking the mechanism that brings ligands and receptors into sufficient proximity to bind. Furthermore, while receptor diffusion has traditionally been viewed as the rate-limiting step, this view struggles to explain how initial contact occurs, particularly in crowded cellular environments where receptor densities can significantly exceed viral ligand densities. Additionally, existing models often treat the cell as a purely elastic medium or ignore cytoskeletal deformation, failing to capture the time-dependent viscoelastic resistance that governs the kinetics of engulfment.

Methodology
The authors propose a unified continuum model rooted in the Onsager variational principle. The system is described by a generalized coordinate, the engulfment depth $h(t)$ , and a free energy landscape $E(h)$ that integrates four distinct contributions:

Entropic (Depletion) Energy: Arising from the exclusion of smaller biomolecules (crowding agents) from the gap between the virus and the membrane.
Binding Energy: The energy gain from ligand-receptor complex formation.
Membrane Deformation Energy: Modeled via the Helfrich-Canham Hamiltonian, accounting for bending rigidity and surface tension.
Cytoskeleton Viscoelastic Energy: Modeled using a standard linear solid to capture time-dependent creep and stress relaxation.

The dynamics are derived by minimizing the sum of the free energy dissipation rate and the Rayleigh dissipation function. The process is divided into two phases:

Phase 1 (Initiation): The virus approaches the membrane from a separation of $2r$ to contact ( $h=0$ ). This phase is driven solely by the entropic depletion force against fluid viscous drag (Stokes friction).
Phase 2 (Engulfment): Once contact is established, ligand-receptor binding sustains the process. The kinetics are governed by the viscoelastic response of the cytoskeleton, utilizing the elastic-viscoelastic correspondence principle to map Hertzian contact theory to a time-dependent creep compliance $\Phi(t)$ .

The model yields a kinetic equation of the form $\zeta_{eff}(h)\dot{h} = F(h)$ , where $F(h) = -\partial E/\partial h$ is the thermodynamic driving force and $\zeta_{eff}$ is the generalized friction coefficient.

Key Contributions

Resolution of the Initiation Problem: The paper identifies entropic depletion forces generated by nanoscale biomolecules in crowded environments as the primary driving mechanism that brings the virus into proximity with the membrane before specific binding occurs.
Unified Kinetic Framework: The model integrates the hydrodynamic approach phase and the viscoelastic wrapping phase into a single Onsager variational framework, bridging the gap between non-equilibrium soft matter physics and cellular mechanics.
Asymptotic Consistency: The derived curvature-dependent depletion energy for a spherical particle rigorously reduces to the classic Asakura-Oosawa result for flat surfaces in the limit of large particle radius, validating the geometric derivations.
Kinetic Phase Diagrams: The authors derive solubility conditions for the kinetic equations, defining the permissible parameter space (particle size, ligand density, cell stiffness) for complete engulfment.

Results

Timescale Compatibility: The calculated depletion-driven approach time ( $t_{approach} \approx 7.14 \times 10^{-5}$ s) is of the same order of magnitude as the characteristic time for stable receptor-ligand binding ( $\sim 10^{-5}$ s). This ensures the depletion force can maintain proximity long enough for specific bonds to form, solving the initiation bottleneck.
Critical Concentration: A critical concentration of crowding agents ( $c_{crit}$ ) is identified, below which the entropic force is insufficient to overcome the membrane deformation barrier to initiate contact.
Size Window and Optimal Radius: The model predicts a finite size window ( $R_{min} < R < R_{max}$ $R_{min} < R < R_{ma x}$ ) for complete engulfment.
- $R_{min}$ is determined by the prohibitive membrane bending energy for small particles.
- $R_{max}$ is determined by the prohibitive cytoskeletal deformation energy for large particles.
- An optimal virus radius ( $R_{opt}$ ) exists that minimizes the complete engulfment time. Crucially, $R_{opt}$ decreases monotonically as the ligand-receptor binding energy density ( $a$ ) increases.
HIV-1 Correlation: Using physiologically realistic parameters for HIV-1 gp120-CD4 interactions ( $a \sim 2.5 \times 10^{-4}$ J/m $^2$ ), the model predicts an optimal radius of $\approx 48.5$ nm, which closely matches the actual radius of HIV-1 ( $\approx 50$ nm).
Cell Stiffness Effects: Stiffer cells (higher Young's modulus $E_c$ ) lead to longer engulfment times and a narrower size window. A critical cell stiffness exists beyond which complete engulfment is impossible for a given ligand density.
Phase Boundaries: The phase boundaries in the $(R, a)$ and $(R, E_c)$ planes are identified as saddle-node bifurcations where the Onsager solubility condition is marginally satisfied.

Significance and Claims
The paper claims to provide a variational foundation for cellular uptake that addresses the long-standing initiation problem overlooked by diffusion-centric models. By demonstrating that entropic forces can drive the initial approach, the model offers a physically consistent explanation for how viruses and nanoparticles overcome the energy barrier to contact the membrane.

The authors emphasize that the model's predictions—specifically the optimal size dependence on binding energy and the critical role of cytoskeletal viscoelasticity—are consistent with experimental observations, including the size of HIV-1 and the uptake efficiency of nanoparticles. The framework suggests that the "optimal" size of viruses may be an evolutionary adaptation to minimize the energetic and kinetic costs of entry given their specific ligand-receptor binding affinities. The work is presented as a step toward a more comprehensive understanding of viral infection mechanisms and the design of nanoparticle-based drug delivery systems, highlighting the necessity of considering both entropic initiation and viscoelastic kinetics.