Cell Division Changes Fate Decisions in a Genetic Toggle Switch

Sintesi Tecnica: Le variazioni della divisione cellulare determinano le decisioni del destino in un interruttore genetico a toggle

Problema
Le reti di regolazione genica (GRN) governano le decisioni del destino cellulare attraverso dinamiche multistabili, con l'interruttore a toggle genetico che funge da modello canonico per tale comportamento. Sebbene i modelli classici di espressione genica incorporino spesso gli effetti della divisione cellulare tramite un termine di "diluizione efficace" aggiunto al tasso di degradazione, questi modelli trattano la divisione come un processo continuo piuttosto che come un evento esplicito. Al contrario, i meccanismi di divisione esplicita rimuovono le proteine in modo discreto alla fine di ogni ciclo cellulare. Sebbene studi precedenti abbiano osservato che l'integrazione completa della divisione possa ridurre la regione di bistabilità, rimane poco chiaro se, all'interno di regimi parametrici in cui entrambi i modelli esibiscono bistabilità, condizioni iniziali identiche conducano allo stesso stato cellulare stabile. Questo articolo affronta il vuoto nella comprensione di come la modellazione esplicita della divisione cellulare impatti la dinamica e le previsioni del destino nelle reti regolatorie multistabili.

Metodologia
Per investigare analiticamente l'impatto della divisione, gli autori propongono una versione semplificata e analiticamente trattabile dell'interruttore a toggle: il Boolean Toggle Switch (BTS). Questo modello è derivato prendendo il limite di coefficienti di Hill infiniti nelle classiche equazioni differenziali ordinarie (ODE) accoppiate, sostituendo le funzioni di Hill con funzioni a gradino di Heaviside. Questa binarizzazione della sintesi proteica (ON/OFF basata su soglie di concentrazione) preserva le caratteristiche dinamiche fondamentali del toggle switch pur permettendo la derivazione di espressioni in forma chiusa.

Lo studio confronta due versioni distinte del BTS:

1. Standard BTS: Segue le convenzioni classiche in cui il volume cellulare è fisso ($V=1$). Gli effetti della divisione sono incorporati esclusivamente attraverso un tasso di diluizione efficace ($\delta = \ln 2 / T$) aggiunto ai tassi di degradazione attiva.
2. BTS con Divisione: Modella esplicitamente la crescita cellulare esponenziale ($V(t) = V_0 e^{\delta t}$) e la divisione cellulare discreta. In questa versione, il numero di proteine viene dimezzato alla fine di ogni ciclo cellulare ($t=T$), e non vi è alcun termine di diluizione efficace nelle equazioni differenziali.

Gli autori derivano espressioni analitiche per le separatrici (le varietà stabili degli equilibri di sella che dividono il piano di fase nei bacini di attrazione) per entrambi i modelli. Assumono una degradazione proteica attiva nulla ($\beta_A = \beta_B = 0$) per concentrarsi sugli effetti della divisione, notando che l'emivita delle proteine nelle cellule mammaliani spesso supera la durata del ciclo cellulare. Inoltre, vengono eseguite simulazioni stocastiche (500 run per condizione iniziale) per validare l'analisi deterministica rispetto al rumore intrinseco, particolarmente nei regimi a basso numero di molecole.

Contributi Chiave e Risultati
Il contributo primario del lavoro è la derivazione delle separatrici analitiche per entrambi i modelli BTS (standard e con inclusione della divisione), rivelando che i due modelli possiedono confini distinti per le decisioni del destino cellulare.

• Traiettorie Divergenti: L'analisi dimostra che la divisione cellulare può reindirizzare traiettorie con identiche condizioni iniziali verso stati stabili opposti. Nello specifico, la separatrice per il BTS con divisione non coincide con quella del BTS standard.
• Regione di Disaccordo (ROD): Gli autori definiscono una "Regzione di Disaccordo" (ROD) nel piano di fase in cui i due modelli prevedono stati stabili differenti per lo stesso numero iniziale di proteine.
  • Quando i parametri di sintesi e soglia soddisfano $\alpha_A \lambda_B > \alpha_B \lambda_A$, la ROD si trova al di sotto della separatrice standard.
  • Quando $\alpha_A \lambda_B < \alpha_B \lambda_A$, la ROD si trova a sinistra della separatrice standard.
  • L'area della ROD è strettamente positiva (eccetto quando i rapporti tra sintesi e soglia sono perfettamente bilanciati) ed è massimizzata quando i rapporti tra i tassi di sintesi proteica o le soglie di repressione sono elevati.
• Validazione Stocastica: Le simulazioni stocastiche confermano che la ROD deterministica è robusta al rumore intrinseco. Le condizioni iniziali che cadono all'interno della ROD deterministica risultano in un'alta proporzione di simulazioni stocastiche in cui gli stati finali dei due modelli divergono.

Significato e Rivendicazioni
L'articolo sostiene che la divisione cellulare possa rimodellare fondamentalmente i confini del destino nelle reti regolatorie multistabili. Trascurando il meccanismo esplicito della divisione e affidandosi esclusivamente alla diluizione efficace, i modelli possono prevedere erroneamente i destini cellulari per un sottoinsieme significativo di condizioni iniziali (la ROD). Gli autori concludono che queste discrepanze non sono meri artefatti di specifiche scelte parametriche, ma rappresentano una differenza fondamentale nel modo in cui i due approcci di modellazione gestiscono la dinamica della crescita e della divisione cellulare. Di conseguenza, i risultati suggeriscono che discrepanze simili possano sorgere in reti di regolazione genica multistabili più complesse, implicando che la modellazione esplicita della divisione sia cruciale per previsioni accurate del processo decisionale cellulare.