Interaction between cell membranes and protein inclusions in the large-deformation regime

Technische Samenvatting: Interactie tussen Celmembranen en Eiwitinsluitingen in het Regime van Grote Deformaties

Probleemstelling
Biologische membranen zijn dynamische oppervlakken waarvan de morfologie en functie kritisch worden bepaald door eiwitinsluitingen (protein inclusions, PIs). Hoewel de interactie tussen PIs en lipidemembranen uitgebreid is gemodelleerd in het regime van kleine deformaties (SD) met behulp van perturbatietheorie rond vlakke configuraties, omvatten veel biologisch relevante processen sterke membraandeformaties. Dit omvat scenario's waarin PIs met hoge snelheden bewegen, wat significante invaginaties induceert, of waar meerdere PIs dicht bij elkaar gepositioneerd zijn, wat steile gradiënten creëert die niet beschreven kunnen worden door lineaire theorieën. Dit werk adresseert de kloof in het begrip van membraan-eiwitinteracties in het regime van grote deformaties (LD), met een specifieke focus op het falen van SD-benaderingen en het ontstaan van complexe elastische gedragingen.

Methodologie
De auteurs maken gebruik van een hybride aanpak waarbij eindelement-numerieke simulaties (FE) worden gecombineerd met benaderende analytische afleidingen.

• Numeriek Kader: De studie maakt gebruik van de eindelementsoftware IRENE om de volledige niet-lineaire vormvergelijking op te lossen, afgeleid van de Helfrich vrije energie. Het membraan wordt gemodelleerd als een tweedimensionaal manifold ingebed in drie dimensies met behulp van de Monge-gauge ($z(x)$). De leidende vergelijking bevat termen voor de gemiddelde kromming ($H$), Gaussische kromming ($K$), buigstijfheid ($\kappa$) en oppervlaktespanning ($\sigma$).
• Randvoorwaarden: Simulaties leggen specifieke beperkingen op, waaronder gepinde membraanranden op een straal $R$, vaste contacthoeken ($\alpha$) bij de eiwitinterface (gemodelleerd als oneindig rigide), en no-slip condities voor membraanstromingen.
• Analytische Benadering: Om fysiek inzicht te verkrijgen, leiden de auteurs een benaderende analytische oplossing af voor het LD-regime door een aanname van nul gemiddelde kromming ($H=0$) te doen. Dit leidt tot een catenoïde profiel, wat expliciete uitdrukkingen voor membraanhoogte en kracht mogelijk maakt als functie van de contacthoek en verplaatsing.
• Stromingsdynamica: In aanwezigheid van membraanstromingen lost de studie gekoppelde partiële differentiaalvergelijkingen (PDE's) op die massa-behoud, momentumbalans (inclusclusief convectieve termen en viskeuze spanningen) en de vormvergelijking beschrijven, terwijl de externe vloeistof wordt genegeerd om de Stokes-paradox binnen het computationele domein te vermijden.

Belangrijkste Bijdragen en Resultaten

1. Beperkingen van Lineaire Theorie:
   Door FE-oplossingen te vergelijken met de lineaire vergelijking (geldig voor kleine gradiënten), demonstreren de auteurs dat de lineaire benadering faalt wanneer de contacthoek groter is dan $\tan \alpha \approx 0,4$ of wanneer de verticale verplaatsingen groot zijn. In deze LD-regimes moet de volledige niet-lineaire vormvergelijking worden opgelost om de correcte membraangeometrie te vatten.

2. Benaderende Analytische Oplossing (Nul Gemiddelde Kromming):
   De auteurs stellen een analytische oplossing voor op basis van de aanname van nul gemiddelde kromming, wat een catenoïde profiel oplevert. Hoewel dit een benadering is voor het algemene geval, legt het de kwalitatieve kenmerken van het LD-regime succesvol vast. Het voorspelt dat de spontane verticale verplaatsing ($h$) van een eiwit met een vaste contacthoek logaritmisch schaalt met de domeingrootte ($h \sim \ln R$) voor grote $R$.

3. Kracht-Verplaatsingsrelatie:
   De studie berekent de kracht die het membraan uitoefent op een PI. Een belangrijke bevinding is dat deze kracht een niet-monotoon gedrag vertoont ten opzichte van de verticale verplaatsing van het eiwit. De kracht neemt toe tot een maximum en neemt daarna af, een kenmerk dat onafhankelijk lijkt te zijn van de specifieke eiwitgeometrie. Deze niet-monotoniciteit suggereert het potentieel voor het ontstaan van dynamische instabiliteiten die relevant zijn voor tubulusvorming en -retractie. De grootte van deze krachten wordt geschat in de orde van 1–10 pN, wat consistent is met de krachten uitgeoefend door moleculaire motoren.

4. Membraan-Gemedieerde Interactie Tussen Twee PIs:
   Voor een systeem van twee interagerende PIs volgt de membraan-gemedieerde potentiële energie geen eenvoudige machtswetmatige afname zoals in SD-regimes. In plaats daarvan vertoont het potentieel een sub-machtswetmatige afname, die langzamer afneemt dan een machtswet met de afstand tussen de eiwitten. Dit weerspiegelt de complexe aard van het elastische medium. Bovendien is de interactiepotentiaal een oneven functie van de eiitorientatie:

• PIs met identieke oriëntaties (dezelfde "lading") stoten elkaar af.
• PIs met tegengestelde oriëntaties (tegenovergestelde "ladingen") trekken elkaar aan.
  Dit bevestigt de analogie tussen PI-oriëntatie en elektrische lading in het LD-regime. Het systeem energetisch gunstig een antisymmetrische configuratie, wat potentieel leidt tot antiferromagnetische-achtige patronen van afwisselende eiwitten.

5. Door Stroming Geïnduceerde Deformaties en Karakteristieke Snelheid:
   In aanwezigheid van membraanstromingen identificeren de auteurs een karakteristieke snelheidschaal, $v^* = \kappa / (L\eta)$, die overeenkomt met een Scriven-Love getal van één. Deze snelheid scheidt twee verschillende regimes:

• $v < v^*$: Stroming heeft een verwaarloosbaar effect op de membraanvorm; buigstijfheid domineert.
• $v > v^*$: Door stroming geïnduceerde deformaties treden op, met een karakteristieke golflengte $\lambda \sim \kappa/(v\eta)$.
  De studie suggereert dat boven deze drempel de competitie tussen viskeuze en buigkrachten leidt tot significante vormveranderingen.

Significantie en Claims
Het artikel claimt kwantitatieve voorspellingen te bieden voor membraan-eiwitsystemen in biologisch relevante scenario's die grote deformaties omvatten, waarmee zij voorbij de beperkingen van perturbatieve SD-theorieën gaat. De resultaten hebben directe implicaties voor het begrijpen van:

• EiwitSortering en Clustering: De niet-monotone krachten en de specifieke aantrekkings-/afstotingsregels op basis van oriëntatie suggereren mechanismen voor hoe eiwitten zich kunnen organiseren in clusters of patronen op het membraan.
• Membraantransport (Trafficking): De identificatie van een kritische snelheid voor door stroming geïnduceerde deformaties impliceert dat eiwitbeweging actief membraanvormveranderingen kan aansturen, wat potentieel de vorming en het transport van vesikels beïnvloedt.
• Biologische Voorbeelden: De auteurs illustreren de relevantie van hun bevindingen met behulp van bacteriorhodopsine (BR) op gigantische unilamellaire vesikels (GUVs). Zij schatten dat de laterale diffusiesnelheid van BR-moleculen voldoende is om door stroming gedreven deformaties met golflengten vergelijkbaar met de afstand tussen de eiwitten te induceren, wat suggereert dat stromingsbemiddelde interacties de BR-patronering aanzienlijk kunnen beïnvloeden.

In algemene zin vestigt dit werk een kader voor het analyseren van de complexe, niet-lineaire wisselwerking tussen membraanelasticiteit, eiwitgeometrie en hydrodynamische stromingen, wat een realistischerere beschrijving biedt van de membraanmechanica in dynamische cellulaire omgevingen.