Theory of adhesion-driven self-organisation in growing tissues

1. Enunciazione del Problema

L'articolo affronta una lacuna fondamentale nella comprensione dell'auto-organizzazione cellulare: la distinzione tra invasione cellulare (colonizzazione coesa dello spazio, ad esempio nella guarigione delle ferite) e formazione di pattern spaziali (emergere di eterogeneità, ad esempio nello sviluppo embrionale o nelle metastasi tumorali). Tradizionalmente, questi processi sono modellati come distinti.

• L'invasione è spesso modellata tramite sistemi di reazione-diffusione (ad esempio Fisher-KPP), mentre la formazione di pattern è tipicamente spiegata da meccanismi di Turing che richiedono almeno due specie interagenti (attivatore/inibitore).
• La Sfida: Molti sistemi biologici coinvolgono una singola specie cellulare dove l'adesione guida sia l'aggregazione che la migrazione. I meccanismi di Turing classici falliscono qui a causa della mancanza di una seconda specie o di rapporti di diffusività appropriati. Inoltre, i modelli esistenti di adesione a singola specie sono spesso non locali, rendendo difficile la trattabilità analitica.
• Obiettivo: Sviluppare un quadro teorico unificato che spieghi come una popolazione a singola specie possa transitare tra invasione coesa e formazione di pattern frammentati esclusivamente attraverso la regolazione dell'adesione cellula-cellula e della proliferazione.

2. Metodologia

Gli autori impiegano un approccio matematico multiscala, passando da un modello integrale non locale a un'equazione differenziale alle derivate parziali (PDE) locale suscettibile di analisi asintotica.

A. Derivazione del Modello

1. Punto di partenza non locale: Iniziano con un modello continuo che descrive la densità cellulare $\rho(x,t)$ incorporando:
   • Auto-diffusione: Diffusione dipendente dalla densità $D(\rho)$ (rappresentante l'inibizione da contatto/esclusione di volume).
   • Adesione: Un termine attrattivo non locale modellato come una convoluzione $(W * \rho)$ con una mobilità basata sull'adesione $m(\rho)$.
   • Proliferazione: Crescita logistica $r\rho(1-\rho)$.
2. Approssimazione a corto raggio: Per superare l'intrattabilità analitica dei termini non locali, applicano uno sviluppo asintotico accoppiato assumendo che il raggio di interazione dell'adesione sia piccolo rispetto alle dimensioni del tessuto. Questo converte la convoluzione non locale in una PDE locale del quarto ordine:
   $$\partial_t \rho = \nabla \cdot \left[ (d(\rho) - \omega m(\rho)) \nabla \rho - m(\rho) \nabla \Delta \rho \right] + \rho(1-\rho)$$
   Qui, $\omega$ rappresenta la forza efficace dell'adesione, e il termine $-\nabla \cdot (m(\rho)\nabla \Delta \rho)$ agisce come un termine di aggregazione di ordine superiore (simile alle equazioni di Cahn-Hilliard).

B. Tecniche Analitiche

• Analisi di stabilità lineare: Utilizzata per determinare la stabilità degli stati stazionari omogenei ($\rho=0$ e $\rho=1$) e derivare le relazioni di dispersione per la formazione di pattern.
• Analisi delle onde viaggianti: Gli autori analizzano soluzioni di onde viaggianti unidimensionali $\rho(x-ct)$ per caratterizzare le velocità di invasione e i profili del fronte.
• Sviluppi asintotici accoppiati: Utilizzati per derivare approssimazioni esplicite per le velocità delle onde e le forme dei fronti nei limiti di adesione debole ($\mu \to +\infty$) e forte ($\mu \to -\infty$).
• Simulazioni numeriche: Schemi a volumi finiti (separando le parti conservative e non conservative) validano le previsioni analitiche in 1D e 2D.

C. Due Varianti del Modello

1. Modello I (Adesione non saturata): $m(\rho) = \rho$. L'adesione scala linearmente con la densità.
2. Modello II (Adesione saturata/dipendente dalla densità): $m(\rho) = \rho(1-\rho)$. L'adesione raggiunge un picco a densità intermedia e si annulla alla confluenza ($\rho=1$), mimando la regolazione biologica.

3. Contributi e Risultati Chiave

A. Unificazione di Invasione e Pattern

Lo studio dimostra che l'invasione e la formazione di pattern non sono meccanismi distinti ma regimi dinamici diversi dello stesso sistema guidato dall'adesione, controllati dal parametro $\mu = \alpha - \omega$ (l'equilibrio tra diffusione repulsiva e adesione attrattiva).

B. Regime I: Adesione Debole ($\mu \gg 0$)

• Dinamica: Il tessuto mostra onde viaggianti robuste e monotone.
• Meccanismo: La diffusione domina l'adesione. Il fronte di invasione è "spinto" dalla dinamica di massa piuttosto che "tratto" dal bordo di attacco.
• Velocità: La velocità dell'onda $c$ converge alla previsione del modello Porous-Fisher: $c \sim \sqrt{\mu/2}$.
• Morfologia: Fronti planari e stabili in 2D.

C. Regime II: Adesione Intermedia ($-2 < \mu \leq \mu_c$)

• Transizione: All'aumentare dell'adesione, il fronte perde la monotonia.
• Code oscillanti: In 1D, il profilo di densità dietro il fronte sviluppa oscillazioni smorzate.
• Instabilità 2D: In 2D, i fronti planari diventano instabili rispetto a perturbazioni trasversali, portando a instabilità a dita. Piccole protuberanze crescono in dita, assomigliando all'invasione disregolata osservata nel cancro.
• Soglia critica: Un valore critico $\mu_c \approx 1.18$ segna la transizione da fronti monotoni a non monotoni.

D. Regime III: Adesione Forte ($\mu < -2$)

• Formazione di pattern: Lo stato omogeneo $\rho=1$ diventa instabile.
• Dinamica: Il sistema forma aggregati localizzati spazialmente (macchie) dietro il fronte di invasione.
• Struttura: Gli aggregati seguono un profilo di funzione di Bessel ($J_0$) in 2D, con una scala di lunghezza caratteristica $\ell \sim 1/\sqrt{|\mu|}$.
• Significato: Questo dimostra che un sistema a singola specie può generare pattern spaziali robusti (macchie) senza un secondo morfogeno, guidato puramente dalla competizione tra repulsione a corto raggio e adesione a lungo raggio.

E. Il Ruolo della Regolazione Dipendente dalla Densità (Modello II)

• Problema: Nel Modello I, l'adesione forte porta a frammentazione e instabilità.
• Soluzione: Nel Modello II, dove l'adesione si annulla vicino alla saturazione ($m(1)=0$), l'instabilità è soppressa.
• Risultato: Anche con parametri di adesione forti, il tessuto mantiene un'invasione coesa e monotona. Il fronte rallenta ma non si rompe in dita o macchie. Ciò suggerisce che i tessuti biologici regolano la forza dell'adesione in base alla densità locale per garantire una morfogenesi robusta.

4. Significato

• Avanzamento Teorico: L'articolo fornisce un legame matematico rigoroso tra invasione cellulare e formazione di pattern, mostrando che emergono da un singolo meccanismo (adesione) piuttosto che richiedere quadri distinti (come il modello a due specie di Turing).
• Insight Biologico: Identifica la forza dell'adesione e la sua regolazione dipendente dalla densità come determinanti chiave del comportamento tissutale.
  • Adesione forte non regolata $\rightarrow$ Dita e frammentazione simili al cancro.
  • Adesione regolata (saturante ad alta densità) $\rightarrow$ Espansione coesa del tessuto e guarigione delle ferite.
• Spiegazione Meccanicistica: Offre una spiegazione fisica per le "dita" nei tessuti epiteliali e nell'invasione tumorale senza invocare gradienti chimici esterni o multiple tipi cellulari.
• Trattabilità Analitica: Derivando un'approssimazione locale di modelli di adesione non locali, gli autori abilitano l'uso di potenti strumenti asintotici (onde viaggianti, analisi di stabilità) che sono difficili da applicare a equazioni integrali non locali.

Conclusione

Falcó et al. stabiliscono che l'interazione tra adesione cellula-cellula, auto-diffusione e proliferazione crea uno spettro di comportamenti collettivi. Il modello prevede che, mentre l'adesione debole garantisce un'espansione tissutale liscia, l'adesione forte non regolata porta a instabilità e formazione di pattern. Crucialmente, l'introduzione di una regolazione dell'adesione dipendente dalla densità agisce come una "valvola di sicurezza" biologica, sopprimendo queste instabilità per ripristinare la crescita tissutale coesa, fornendo una teoria unificata per comprendere la morfogenesi e la progressione della malattia.