Morphogenesis driven by nematic defects in active biological networks

Technische Samenvatting: Morfogenese Gedreven door Nematiche Defecten in Actieve Biologische Netwerken

Probleemstelling
Het artikel behandelt de fundamentele uitdaging van het modelleren van cellulaire morfogenese, specif으로 hoe biologische weefsels overgaan van ongeordende toestanden naar geordende structuren met specifieke vormen. Hoewel is vastgesteld dat cytoskeletfilamenten zich kunnen organiseren in nematische fasen met een gedeeltelijke oriëntatieorde, blijft de exacte mechanische rol van topologische defecten binnen deze fasen bij het aansturen van weefseldeformatie en patroonvorming onduidelijk. Bestaande modellen vertrouwen vaak op lineaire benaderingen of vooraf gedefinieerde morfogeen-gradiënten, die mogelijk niet in staat zijn om de complexe, niet-lineaire interactie tussen nematische orde, elasticiteit en door spanning gedreven groei te vangen, zoals waargenomen bij regenererende organismen zoals Hydra. De auteurs streven ernaar een continuümframework te ontwikkelen waarbij levende weefsels worden behandeld als actieve biologische netwerken (ABN's), waarbij vormveranderingen voortkomen uit de koppeling van nematische orde, elasticiteit en adaptieve relaxatie, zonder de groeipatronen vooraf voor te schrijven.

Methodologie
De auteurs stellen een quasi-statisch continuümframework voor waarbij levende weefsels worden gemodelleerd als actieve biologische netwerken bestaande uit gekruiste filamentketens die in staat zijn een nematische fase te vormen. De kernaannames en wiskundige formulering omvatten:

• Kinematica en Groei: Het weefsel wordt gemodelleerd als lokaal incompressibel. Morfogenese wordt beschreven als een sequentie van evenwichtstoestanden. Groei wordt niet voorgeschreven, maar wordt gedreven door lokale mechanische spanning. De vervormingsgradiënt $\mathbf{F}$ wordt multiplicatief gedecomposeerd in een nematische staplengte-tensor $\mathbf{L}$, een elastische distortie $\mathbf{F}_e$, een groeitensor $\mathbf{F}_g$ en een relaxatietensor $\mathbf{F}_{rel}$.
• Energetische Formulering: De totale vrije energie $\mathcal{F}$ bestaat uit twee hoofdcomponenten:
  1. Nematische Energie ($\mathcal{F}_n$): Gebaseerd op de Ericksen-continuümtheorie, houdt deze term rekening met ruimtelijke variaties in het directorveld $\mathbf{n}$ en de scalaire ordeparameter $s$. Cruciaal is dat $s$ ruimtelijk mag variëren, waardoor singulariteiten bij topologische defecten verdwijnen en lokale overgangen naar een isotrope fase mogelijk worden.
  2. Elastische Energie ($\mathcal{F}_e$): Een neo-klassieke rekenergie-functie die is uitgebreid om groei en relaxatie te bevatten. Deze is afhankelijk van het elastische deel van de vervorming en vangt de respons van het vezelnetwerk op distorties op.
• Koppelingsmechanismen: Het model integreert vier mechanismen: (a) mechanische spanning nabij topologische defecten, (b) activiteit-geïnduceerde ordening (toename van $s$), (c) door spanning gedreven groei (waarbij nieuw materiaal wordt toegevoegd proportioneel aan de elastische energiedichtheid), en (d) structurele relaxatie (het bijwerken van de spanningvrije configuratie om interne spanningen te verlichten).
• Numerieke Implementatie: Het gekoppelde niet-lineaire probleem wordt opgelost met behulp van een gemengde eindelementformulering, geïmplementeerd in FEniCSx. Om de kop-staart symmetrie van nematische directors en de bijbehorende taksplitsingen (branch cuts) voor half-integer defecten te hanteren, is het directorveld geherformuleerd met behulp van een symmetrische rang-2 tensorveld $\mathbf{M} = \mathbf{n} \otimes \mathbf{n}$. De simulaties maken gebruik van een Newton-methode met line search en incrementele belasting om de evolutie van een initiële isotrope sferoïde naar volwassen morfologieën te volgen.

Belangrijkste Bijdragen

1. Continuümframework voor ABN's: Het artikel introduceert een wiskundig model dat nematische orde, elasticiteit en door spanning gedreven groei/relaxatie koppelt, waarmee voorbij wordt gegaan aan lineaire schaaltheorieën om niet-lineaire morfogenetische uitkomsten te vangen.
2. Door Defecten Gedreven Morfogenese: De studie toont aan dat de topologie van nematische defecten (specifiek hun lading en ruimtelijke arrangement) fungeert als de primaire drijfveer voor vormevolutie. Het model laat zien dat defecten niet louter singulariteiten zijn, maar organiserende centra die spanning lokaliseren en groei sturen.
3. Reproductie van Complexe Morfologieën: Door de initiële defectconfiguratie te variëren, reproduceert het model succesvol drie verschillende morfogenetische uitkomsten waargenomen bij Hydra-regeneratie:
   • Uniaxiale (Hoofd-Voet) Regeneratie: Twee $+1$ defecten aan tegenovergestelde polen drijven verlenging aan en de vorming van een mond (apicale ruptuur) en een voet.
   • Bicefale (Twee-Hoofden) Ontwikkeling: Twee $+1$ defecten in de bovenste hemisfeer, gebalanceerd door $-1/2$ defecten, leiden tot de vorming van twee afzonderlijke koppen.
   • Tentakelvorming: Een configuratie met $+1$ defecten bij tentakelpunten en $-1/2$ defecten bij hun basis drijft de uitgroei van tentakels aan.
4. Mechanisme van Door Spanning Gedreven Groei: De simulaties tonen expliciet aan dat groei gelokaliseerd is in regio's met hoge elastische spanning (nabij defecten), terwijl relaxatie zorgt voor de vergroting van openingen (monden) waar spanning geometrisch kan worden verlicht, in contrast met de persistente spanning bij de voet of tentakelpunten.

Resultaten
Eindelementsimulaties van Hydra-regeneratie onthullen dat:

• Spanningslokalisatie: Distorties in het nematische directorveld nabij topologische defecten genereren hoge mechanische spanningen. Deze spanningen drijven de accumulatie van nieuw materiaal (groei) in hun nabijheid aan.
• Asymmetrische Uitkomsten: In de uniaxiale configuratie zorgt spanningsrelaxatie ervoor dat de apicale pool (mond) openbreekt en vergroot, terwijl de basale pool (voet) gesloten en structureel gespannen blijft.
• Rol van Defecten: $+1$ defecten drijven de vorming van uitstulpingen (protrusions) en weefselruptuur aan. $-1/2$ defecten fungeren primair als structurele stabilisatoren of "topologische balansers" in vroege stadia, hoewel in de tentakelconfiguratie de spanningsaccumulatie bij deze defecten uiteindelijk lokale groei aan de basis van de tentakels triggert.
• Morfologische Evolutie: Het model vangt de transitie van een sferisch, isotroop weefselfragment naar complexe, verlengde structuren met specifieke kenmerken (orale poriën, tentakels) uitsluitend door de evolutie van de nematische ordeparameter en de daaruit voortvloeiende door spanning gedreven groei.

Betekenis en Claims
De auteurs claimen dat hun werk een robuust theoretisch fundament biedt voor het begrijpen van hoe topologische defecten in actieve nematische netwerken complexe morfogenetische processen kunnen aansturen zonder de noodzaak van vooraf gespecificeerde groeipatronen of externe morfogeen-gradiënten. Door het weefsel te behanden als een actief biologisch netwerk waarbij groei een directe reactie is op de mechanische spanning gegenereerd door nematische distorties, overbrugt het model de kloof tussen microscopische cytoskeletorganisatie en macroscopische weefselvorm. De resultaten komen overeen met experimentele observaties bij Hydra, wat suggereert dat de "geheugen" van de lichaamsas-polariteit en de opkomst van specifieke structuren gecodeerd zijn in de topologische lading en de arrangement van nematische defecten. Het artikel positioneert dit framework als een stap naar het begrijpen van de intrinsieke mechanische drijfveren van patroonvorming in levende weefsels, onderscheidend van puur biochemische signaleringsmodellen.