Multicellular simulations with shape and volume constraints using optimal transport

Technische Samenvatting: Meercellige Simulaties met Vorm- en Volumedwang met Behulp van Optimale Transport

Probleemstelling
Levende systemen, variërend van intracellulaire componenten tot macroscopische weefsels en bacteriële kolonies, vertonen complexe dynamiek gedreven door volumexclusie, vorminteracties en vervormbaarheid. Het begrijpen van hoe deze beperkingen leiden tot macroscopische zelforganisatie is een fundamentele uitdaging in de ontwikkelingsbiologie en biofysica. Bestaande computationele modellen staan voor een afweging tussen fysieke nauwkeurigheid, numerieke kosten en wiskundige hanteerbaarheid. Gangbare benaderingen omvatten:

• Punt-deeltjesystemen: Wiskundig eenvoudig maar missen vaak een realistische vormrepresentatie.
• Faseveld- en level-set-modellen: Capabel om willekeurige vormen weer te geven, maar kunnen computergewijs duur zijn en uitdagend voor contactinteracties.
• Cellulaire automata en roostermmodellen: Discretiseren de ruimte in voxels, wat de geometrische flexibiliteit potentieel beperkt.
• Voronoi- en vertexmodellen: Efficiënt voor weefselachtige structuren, maar vertrouwen vaak op zachte dwang voor volume of missen controle over specifieke vormvervormingen.

Er bestaat een aanzienlijke kloof in het modelleren van systemen met strikte volumedwang, willekeurige dynamische vormen en hoge deeltjesaantallen (tientallen tot tienduizenden) binnen een verenigd raamwerk dat zowel individuele vervormbaarheid als collectieve interacties ondersteunt.

Methodologie
De auteurs introduceren een mesoschaal simulatiekader gebaseerd op theorie van optimale transport, specifiek met gebruikmaking van Laguerre-tessellaties.

1. Statisch Model (Laguerre-tessellatie):
   De ruimte die door elk deeltje $i$ wordt ingenomen, wordt gemodelleerd als een Laguerre-cel $L_i$, gedefinieerd als de verzameling punten $x$ die voldoen aan:
   $$L_i = \{x \in \Omega, c(x, x_i) - w_i \leq c(x, x_j) - w_j \text{ voor alle } j\}$$
   Hierbij is $x_i$ de deeltjespositie, $v_i$ het doelvolume, en $c(x, x_i)$ een kostenfunctie die de vorm- en vervormbaarheidseigenschappen van het deeltje codeert. De gewichten $w_i$ (Kantorovich-potentialen) worden uniek bepaald om aan de strikte volumedwang $|L_i| = v_i$ te voldoen.

   • Deze formulering generaliseert Voronoi-diagrammen (waarbij $w_i=0$ en $c$ de Euclidische afstand is) om gebogen grenzen en niet-sferische vormen toe te staan door willekeurige kostenfuncties te kiezen (bijvoorbeeld power-kosten, anisotrope kosten).
   • Lege ruimte wordt gemodelleerd als een extra deeltje met nul kosten.

2. Dynamisch Model:
   De evolutie van de deeltjesposities $x_i$ volgt een differentiaalvergelijking voor eerste-orde gradiëntafdaal:
   $$\dot{x}_i = -\tau_i \nabla_{x_i} T_c + F_i^{\text{point}} + F_i^{\text{surf}} + \dot{\xi}_i$$

   • Ondrukbaarheidskracht: De term $-\nabla_{x_i} T_c$ werkt als een afstotende kracht afgeleid van de totale transportkosten $T_c$. Deze zorgt ervoor dat deeltjes niet overlappen en dwingt op natuurlijke wijze volumedwang af. In tegenstelling tot parenkrachten hangt deze kracht af van de globale geometrie van alle contactoppervlakken.
   • Externe/Interactiekrachten ($F^{\text{point}}$): Omvat zelfaandrijving, zwaartekracht of aantrekking, vergelijkbaar met punt-deeltjesmodellen.
   • Oppervlaktekrachten ($F^{\text{surf}}$): Modelleert druk- en oppervlaktespanningseffecten langs interfaces, waardoor fenomenen zoals celsortering en adhesie mogelijk zijn.
   • Vormdynamiek: Vormen kunnen evolueren via geparametriseerde kosten (bijvoorbeeld ellipsoïden gedefinieerd door covariantiematrices) of level-set potentiaalfuncties die reageren op externe velden (bijvoorbeeld chemotaxis).

3. Numerieke Implementatie:

   • Het probleem wordt opgelost als een semi-discreet optimale transport-probleem op een discreet rooster (voxels).
   • De auteurs maken gebruik van een verdel-en-verover-meerschaalstrategie en een quasi-Newton-oplosser (L-BFGS-B) om de Kantorovich-potentialen $w_i$ efficiënt te berekenen.
   • De implementatie is massaal geparallelliseerd op GPU's met behulp van de KEOPS- en GEOMLOSS-bibliotheken, wat simulaties van tot 50.000 deeltjes in 3D mogelijk maakt met state-of-the-art prestaties.

Belangrijkste Bijdragen

• Verenigd Kader: De methode verenigt concepten uit Voronoi-tessellaties, level-set-methoden en vertexmodellen onder de paraplu van optimale transport. Het ondersteunt willekeurige vormen, strikte volumedwang en complexe interactiemechanismen binnen één enkel kader.
• Strikte Volumedwang: In tegenstelling tot veel bestaande modellen die zachte dwang of relaxatiekrachten gebruiken, garandeert deze benadering wiskundig op elk moment volumedwang door de optimalisatie van Kantorovich-potentialen.
• Schaalbaarheid: De door GPU's versnelde oplosser maakt simulatie van grootschalige systemen (tot 50.000 deeltjes) in 3D mogelijk, wat computationeel onhaalbaar is voor veel traditionele eindige-elementen- of faseveldmethoden.
• Veelzijdigheid in Modelleren: Het kader modelleert succesvol diverse fenomenen, waaronder:
  • Vloeistof-vast faseovergangen in actieve Brownse deeltjes gedreven door vervormbaarheid.
  • Ontstaande oriëntatieorde in zwermende staafvormige deeltjes uitsluitend door volumexclusie.
  • Chemotactische migratie gedreven door vormvervormingen.
  • Celsortering en inslikking-fenomenen in 3D heterogene aggregaten, waarbij resultaten worden gereproduceerd die consistent zijn met de Hypothese van het Differentiële Interfaciale Spanning (DITH).

Resultaten
De auteurs valideren het kader via verschillende benchmarks:

• Groei en Deling: Een 3D-aggregaat groeit van één cel tot 50.000 cellen via deling, waarbij een domein wordt gevuld in ongeveer één dag rekentijd.
• Faseovergangen: In een systeem van actieve Brownse deeltjes drijft de vervormbaarheidsparameter $\alpha$ een overgang van een vloeibare fase (diffusief gedrag) naar een vaste fase (vastgelopen, hexagonale pakking). Het model vangt het verdwijnen van de effectieve diffusie en de verzadiging van de gemiddelde kwadratische verplaatsing in de vaste fase.
• Staafvormige Uitlijning: Simulaties van zelfaangedreven staafvormige deeltjes tonen het spontane ontstaan van langafstands oriëntatieorde aan door botsing-geïnduceerde vervormingen, consistent met theorieën uit de polymerfysica.
• 3D Celsortering: Het model reproduceert diverse sorteerpatronen (inslikking, scheiding, schaakbord) in 3D-mengsels van twee celtypen met verschillende oppervlaktespanningen en zachtheid, overeenkomend met theoretische voorspellingen en experimentele waarnemingen.
• Menigtedynamiek: Het model simuleert menigte-evacuatie, waarbij blijkt dat vervormbare deeltjes (vloeistof-achtig) veilig evacueren, terwijl harde-sfeer deeltjes (vast-achtig) stabiele bogen vormen en vastlopen bij uitgangen.

Betekenis en Beweringen
Het artikel beweert een nieuwe benadering te introduceren voor het modelleren van meercellige systemen die de kloof overbrugt tussen gedetailleerde fysieke beschrijvingen en computationele efficiëntie.

• Validatie: De primaire bijdrage is de validatie van dit nieuwe kader tegen klassieke computationele modellen, waarbij wordt aangetoond dat het gevestigde resultaten kan reproduceren (bijvoorbeeld celsortering, faseovergangen) terwijl het superieure behandeling van volumedwang en 3D-complexiteit biedt.
• Biologische Relevantie: De auteurs benadrukken het potentieel van de methode voor in silico modellering van complexe ontwikkelingsprocessen, zoals blastocyst-ontwikkeling, waarbij zowel individuele celvervormingen als het gedrag van dichte aggregaten cruciaal zijn.
• Theoretische Uitbreiding: Het werk suggereert een nieuw perspectief van optimale transport voor bestaande methoden (punt-deeltje, level-set, vertex) en biedt een weg naar het afleiden van grofkorrelige continuümmodellen (mean-field benaderingen) voor systemen met volume-interacties, een grotendeels open vraag in de wiskundige biologie.
• Toegankelijkheid: De auteurs bieden een vrij beschikbaar Python-implementatie, waardoor de methode toegankelijk is voor verder onderzoek in computationele fysica en biologie.

De auteurs blijven bescheiden, met de opmerking dat hoewel het kader nieuw is, de onmiddellijke focus lag op validatie tegen bestaande modellen. Zij identificeren toekomstige richtingen, zoals het uitbreiden van de actieve vertextheorie naar Laguerre-tessellaties, het verfijnen van de afleiding van continuümlimieten, en het verbeteren van de behandeling van complexe dynamische kosten en oppervlakte-interacties voor simulatie in real-time.