Morphogenesis driven by nematic defects in active biological networks

Sintesi Tecnica: Morfogenesi Guidata da Difetti Nematici in Reti Biologiche Attive

Problematica
Il documento affronta la sfida fondamentale di modellare la morfogenesi cellulare, specificamente come i tessuti biologici passino da stati disordinati a strutture ordinate con forme specifiche. Sebbene sia stabilito che i filamenti del citoscheletro possano organizzarsi in fasi ematitiche con un ordine orientazionale parziale, il ruolo meccanico preciso dei difetti topologici all'interno di queste fasi nel guidare la deformazione tissutale e la formazione di pattern rimane poco chiaro. I modelli esistenti spesso si affidano ad approssimazioni lineari o a gradienti di morfogeni predefiniti, che possono non essere in grado di catturare l'interazione complessa e non lineare tra ordine nematico, elasticità e crescita guidata dallo stress osservata in organismi rigeneranti come Hydra. Gli autori mirano a sviluppare un framework continuo che tratti i tessuti viventi come reti biologiche attive (ABN), dove i cambiamenti di forma emergono dal accoppiamento tra ordine nematico, elasticità e rilassamento adattivo, senza prescrivere a priori il pattern di crescita.

Metodologia
Gli autori propongono un framework continuo quasi-statico che modella i tessuti viventi come reti biologiche attive composte da catene di filamenti reticolati capaci di formare una fase nematica. Le assunzioni fondamentali e la formulazione matematica includono:

• Cinematica e Crescita: Il tessuto è modellato come localmente incomprimibile. La morfogenesi è descritta come una sequenza di stati di equilibrio. La crescita non è prescritta, ma è guidata dallo stress meccanico locale. Il gradiente di deformazione $\mathbf{F}$ è scomposto moltiplicativamente in un tensore di lunghezza di passo nematico $\mathbf{L}$, una distorsione elastica $\mathbf{F}_e$, un tensore di crescita $\mathbf{F}_g$ e un tensore di rilassamento $\mathbf{F}_{rel}$.
• Formulazione Energetica: L'energia libera totale $\mathcal{F}$ comprende due contributi principali:
  1. Energia Nematica ($\mathcal{F}_n$): Basata sulla teoria del continuo di Ericksen, questo termine tiene conto delle variazioni spaziali nel campo del direttore $\mathbf{n}$ e del parametro d'ordine scalare $s$. Fondamentalmente, $s$ è permesso di variare spazialmente, svanendo nei difetti topologici per regolarizzare le singolarità e permettere transizioni locali a una fase isotropa.
  2. Energia Elastica ($\mathcal{F}_e$): Una funzione di energia di deformazione neo-classica estesa per includere crescita e rilassamento. Essa dipende dalla parte elastica della deformazione, catturando la risposta della rete di fibre alle distorsioni.
• Meccanismi di Accoppiamento: Il modello integra quattro meccanismi: (a) stress meccanico vicino ai difetti topologici, (b) ordinamento indotto dall'attività (aumento di $s$), (c) crescita guidata dallo stress (dove nuovo materiale viene aggiunto proporzionalmente alla densità di energia elastica), e (d) rilassamento strutturale (aggiornamento della configurazione priva di stress per alleviare gli stress interni).
• Implementazione Numerica: Il problema accoppiato non lineare viene risolto utilizzando una formulazione agli elementi finiti mista implementata in FEniCSx. Per gestire la simmetria testa-coda dei direttori nematici e le relative interruzioni di ramo per i difetti a semintero, il campo del direttore è riformulato utilizzando un campo tensoriale simmetrico di rango 2 $\mathbf{M} = \mathbf{n} \otimes \mathbf{n}$. Le simulazioni impiegano un metodo di Newton con ricerca lineare e carico incrementale per tracciare l'evoluzione da uno sferoide isotropo iniziale a morfologie mature.

Contributi Chiave

1. Framework Continuo per ABN: Il documento introduce un modello matematico che accoppia l'ordine nematico, l'elasticità e la crescita/rilassamento indotti dallo stress, andando oltre le teorie di shell lineari per catturare esiti morfogenetici non lineari.
2. Morfogenesi Guidata dai Difetti: Lo studio dimostra che la topologia dei difetti nematici (specificamente la loro carica e disposizione spaziale) agisce come il driver primario per l'evoluzione della forma. Il modello mostra che i difetti non sono meriamente singolarità, ma centri organizzatori che localizzano lo stress e dirigono la crescita.
3. Riproduzione di Morfologie Complesse: Variando la configurazione iniziale dei difetti, il modello riproduce con successo tre distinti esiti morfogenetici osservati nella rigenerazione di Hydra:
   • Rigenerazione Unassiale (Testa-Piede): Due difetti $+1$ ai poli opposti guidano l'allungamento e la formazione di una bocca (rottura apicale) e di un piede.
   • Sviluppo Bicefalico (Due Teste): Due difetti $+1$ nell'emisfero superiore, bilanciati da difetti $-1/2$, portano alla formazione di due teste distinte.
   • Formazione di Tentacoli: Una configurazione con difetti $+1$ alle punte dei tentacoli e difetti $-1/2 alle loro basi guida la crescita dei tentacoli.
4. Meccanismo di Crescita Guidato dallo Stress: Le simulazioni mostrano esplicitamente che la crescita è localizzata nelle regioni di alto stress elastico (vicino ai difetti), mentre il rilassamento permette l'allargamento delle aperture (bocche) dove lo stress può essere geometricamente alleviato, in contrasto con lo stress persistente al piede o alle punte dei tentacoli.

Risultati
Le simulazioni agli elementi finiti della rigenerazione di Hydra rivelano che:

• Localizzazione dello Stress: Le distorsioni nel direttore nematico vicino ai difetti topologici generano alti stress meccanici. Questi stress guidano l'accumulo di nuovo materiale (crescita) nelle loro vicinanze.
• Esiti Asimmetrici: Nella configurazione unassiale, il rilassamento dello stress permette al polo apicale (bocca) di rompersi e allargarsi, mentre il polo basale (piede) rimane chiuso e strutturalmente stressato.
• Ruoli dei Difetti: I difetti $+1$ guidano la formazione di protrusioni e la rottura del tessuto. I difetti $-1/2 agiscono principalmente come stabilizzatori strutturali o "bilanciatori topologici" nelle fasi iniziali, sebbene nella configurazione tentacolare l'accumulo di stress in questi difetti finisca per innescare una crescita localizzata alle basi dei tentacoli.
• Evoluzione Morfologica: Il modello cattura la transizione da un frammento di tessuto sferico e isotropo a strutture complesse ed allungate con caratteristiche specifiche (pori orali, tentacoli) esclusivamente attraverso l'evoluzione del parametro d'ordine nematico e la conseguente crescita guidata dallo stress.

Significato e Rivendicazioni
Gli autori sostengono che il loro lavoro fornisca una solida base teorica per comprendere come i difetti topologici in reti ematiche attive possano guidare processi morfogenetici complessi senza la necessità di pattern di crescita pre-specificati o gradienti di morfogeni esterni. Trattando il tessuto come una rete biologica attiva dove la crescita è una risposta diretta allo stress meccanico generato dalle distorsioni ematiche, il modello colma il divario tra l'organizzazione microscopica del citoscheletro e la forma macroscopica del tessuto. I risultati si allineano con le osservazioni sperimentali in Hydra, suggerendo che la "memoria" della polarità dell'asse corporeo e l'emergere di strutture specifiche siano codificati nella carica topologica e nella disposizione dei difetti nematici. Il documento posiziona questo framework come un passo verso la comprensione dei driver meccanici intrinseci della formazione di pattern nei tessuti viventi, distinti dai modelli puramente biochimici di segnalazione.