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

Technische Samenvatting: Entropie-gestuurde Initiatie en Cellulaire Opname Gemedieerd door Visco-elastisch Cytoskelet

Probleemstelling
Er bestaat een fundamentele kloof in bestaande theoretische modellen van receptor-gemedieerde endocytose: de "initiatieproblematiek". Eerdere kaders gaan doorgaans uit van een reeds bestaand contact tussen een virus (of nanopartikel) en het celmembraan, waarbij het mechanisme dat liganden en receptoren in voldoende nabijheid brengt om te binden, wordt genegeerd. Bovendien wordt receptordiffusie traditioneel beschouwd als de snelheidsbepalende stap, maar dit beeld verklaart moeizaam hoe initiële contact tot stand komt, met name in drukke cellulaire omgevingen waar de receptor-dichtheden aanzienlijk hoger kunnen zijn dan de virale ligand-dichtheden. Daarnaast behandelen bestaande modellen de cel vaak als een puur elastisch medium of negeren zij cytoskeletale deformatie, waardoor ze er niet in slagen de tijdsafhankelijke visco-elastische weerstand te vatten die de kinetiek van omhulling (engulfment) beheerst.

Methodologie
De auteurs stellen een verenigd continuümmodel voor dat geworteld is in het Onsager-variatieprincipe. Het systeem wordt beschreven door een gegeneraliseerde coördinaat, de omhullingsdiepte $h(t)$, en een vrije energie-landschap $E(h)$ dat vier verschillende bijdragen integreert:

1. Entropische (depletie) energie: Voortvloeiend uit de uitsluiting van kleinere biomoleculen (crowding agents) uit de spleet tussen het virus en het membraan.
2. Bindingsenergie: De energiewinst uit de vorming van ligand-receptorcomplexen.
3. Membraandeformatie-energie: Gemodelleerd via het Helfrich-Canham-Hamiltoniaan, rekening houdend met buigstijfheid en oppervlaktespanning.
4. Cytoskeletale visco-elastische energie: Gemodelleerd met een standaard lineaire vaste stof om tijdsafhankelijke kruip en spanningsrelaxatie te vangen.

De dynamica worden afgeleid door de som van de vrije energie-dissipatiesnelheid en de Rayleigh-dissipatiefunctie te minimaliseren. Het proces wordt verdeeld in twee fasen:

• Fase 1 (Initiatie): Het virus nadert het membraan vanaf een scheiding van $2r$ tot contact ($h=0$). Deze fase wordt uitsluitend gedreven door de entropische depletiekracht tegen de fluïde viskeuze weerstand (Stokes-wrijving).
• Fase 2 (Omhulling): Zodra contact is gevestigd, ondersteunt ligand-receptorbinding het proces. De kinetiek wordt beheerst door de visco-elastische respons van het cytoskelet, waarbij gebruik wordt gemaakt van het elastisch-visco-elastisch correspondentieprincipe om Hertziaanse contacttheorie te mappen naar een tijdsafhankelijke kruipcompliantie $\Phi(t)$.

Het model levert een kinetische vergelijking op de vorm $\zeta_{eff}(h)\dot{h} = F(h)$, waarbij $F(h) = -\partial E/\partial h$ de thermodynamische drijvende kracht is en $\zeta_{eff}$ de gegeneraliseerde wrijvingscoëfficiënt is.

Belangrijkste Bijdragen

1. Oplossing van de Initiatieproblematiek: De paper identificeert entropische depletiekrachten, gegenereerd door nanoschaal biomoleculen in drukke omgevingen, als het primaire drijfmechanisme dat het virus in de nabijheid brengt van het membraan voordat specifieke binding plaatsvindt.
2. Verenigd Kinetisch Kader: Het model integreert de hydrodynamische benaderingsfase en de visco-elastische omhullingsfase in één enkel Onsager-variatiekader, waarmee de kloof tussen niet-evenwichtige zachte materie-fysica en cellulaire mechanica wordt overbrugd.
3. Asymptotische Consistentie: De afgeleide krommingsafhankelijke depletie-energie voor een sferisch deeltje reduceert rigoureus tot de klassieke Asakura-Oosawa-resultaten voor vlakke oppervlakken in het limiet van een grote deeltjesstraal, wat de geometrische afleidingen valideert.
4. Kinetische Fase-diagrammen: De auteurs leiden oplossingsvoorwaarden af voor de kinetische vergelijkingen, waarmee de toegestane parameterruimte (deeltjesgrootte, ligand-dichtheid, celstijfheid) voor volledige omhulling wordt gedefinieerd.

Resultaten

• Tijdschaal Compatibiliteit: De berekende depletie-gedreven benaderingstijd ($t_{approach} \approx 7.14 \times 10^{-5}$ s) is van dezelfde orde grootte als de karakteristieke tijd voor stabiele receptor-ligand binding ($\sim 10^{-5}$ s). Dit zorgt ervoor dat de depletiekracht de nabijheid lang genoeg kan handhaven zodat specifieke bindingen kunnen vormen, wat de initiatie-bottleneck oplost.
• Kritische Concentratie: Een kritische concentratie van crowding agents ($c_{crit}$) wordt geïdentificeerd, waaronder de entropische kracht onvoldoende is om de membraandeforlatiebarrière te overwinnen om contact te initiëren.
• Groottevenster en Optimale Straal: Het model voorspelt een eindig groottevenster ($R_{min} < R < R_{max}$) voor volledige omhulling.
  • $R_{min}$ wordt bepaald door de belemmerende membraanbuigenergie voor kleine deeltjes.
  • $R_{max}$ wordt bepaald door de belemmerende cytoskeletale deformatie-energie voor grote deeltjes.
  • Er bestaat een optimale virusstraal ($R_{opt}$) die de volledige omhullings-tijd minimaliseert. Cruciaal is dat $R_{opt}$ monotoon afneemt naarmate de ligand-receptor bindingsenergiedichtheid ($a$) toeneemt.
• HIV-1 Correlatie: Gebruikmakend van fysiologisch realistische parameters voor HIV-1 gp120-CD4 interacties ($a \sim 2.5 \times 10^{-4}$ J/m$^2$), voorspelt het model een optimale straal van $\approx 48.5$ nm, wat nauw overeenkomt met de werkelijke straal van HIV-1 ($\approx 50$ nm).
• Effecten van Celstijfheid: Stijvere cellen (hogere Young-modulus $E_c$) leiden tot langere omhullings-tijden en een nauwer groottevenster. Er bestaat een kritische celstijfheid waarna volledige omhulling onmogelijk is voor een gegeven ligand-dichtheid.
• Fasegrenzen: De fasegrenzen in de $(R, a)$ en $(R, E_c)$ vlakken worden geïdentificeerd als saddle-node bifurcaties waar de Onsager-oplosbaarheidsvoorwaarde marginaal wordt voldaan.

Betekenis en Claims
De paper claimt een computationeel fundament te bieden voor cellulaire opname dat de langdurige initiatieproblematiek adresseert die door diffusie-gecentreerde modellen werd over het hoofd gezien. Door aan te tonen dat entropische krachten de initiële benadering kunnen aandrijven, biedt het model een fysiek consistente verklaring voor hoe virussen en nanopartikels de energiebarrière naar contact overwinnen.

De auteurs benadrukken dat de voorspellingen van het model — speciffisch de optimale grootte-afhankelijkheid van de bindingsenergie en de cruciale rol van de cytoskeletale visco-elasticiteit — consistent zijn met experimentele observaties, waaronder de grootte van HIV-1 en de opname-efficiëntie van nanopartikels. Het kader suggereert dat de "optimale" grootte van virussen een evolutionaire adaptatie kan zijn om de energetische en kinetische kosten van entree te minimaliseren, gegeven hun specifieke ligand-receptor bindingsaffiniteiten. Het werk wordt gepresenteerd als een stap naar een meer uitgebreid begrip van virale infectiemechanismen en het ontwerp van nanopartikel-gebaseerde medicijnleveringssystemen, waarbij het belang van het overwegen van zowel entropische initiatie als visco-elastische kinetiek wordt onderstreept.