Bifurcations and multistability in inducible three-gene toggle switch networks

Problema
I sistemi biologici si affidano a reti di regolazione genica capaci di multistabilità per guidare distinti fenotipi cellulari, come il differenziamento e le decisioni del destino dello sviluppo. Sebbene gli interruttori a due geni (toggle switch) bistabili siano ben caratterizzati, estendere questi concetti a reti di ordine superiore (ad esempio, sistemi a tre geni) per comprendere la tristabilità e dinamiche più complesse rimane una sfida. Un limite primario nei modelli teorici esistenti è il trattamento dei parametri regolatori (costanti di dissociazione, tassi di trascrizione) come variabili astratte scollegate dai meccanismi biologici specifici che le controllano in vivo. Nello specifico, il ruolo delle molecole effettrici (induttori) nel modulare l'attività dei fattori di trascrizione allosterici è spesso omesso o semplificato. Questo studio affronta il divario tra la topologia astratta della rete e i meccanismi biologici sperimentalmente sintonizzabili esaminando come le concentrazioni di effettore guidino gli spostamenti dinamici negli interruttori a tre geni.

Metodologia
Gli autori costruiscono un quadro matematico per una rete di tre geni a repressione mutua (un interruttore a scatto) utilizzando modelli a funzione di Hill derivati da principi di meccanica statistica.

• Modello di base: Il sistema è definito da equazioni differenziali che descrivono la produzione e la degradazione di tre repressori ($R_1, R_2, R_3$). Il modello incorpora un'interpretazione termodinamica dell'induzione utilizzando il modello di Monod-Wyman-Changeux (MWC). Qui, la concentrazione totale di un fattore di trascrizione è scalata da una funzione di probabilità, $p_{act}(c)$, che dipende dalla concentrazione di un effettore allosterico ($c$). Questa probabilità determina la frazione di repressori attivi.
• Simmetria e perturbazione: Lo studio assume inizialmente un sistema simmetrico in cui tutti i repressori legano con uguale affinità. Successivamente, introduce sistematicamente perturbazioni per rompere questa simmetria, analizzando come i bias nelle forze di repressione influenzino le soglie di biforcazione.
• Estensioni dell'auto-attivazione: Il modello viene esteso per includere l'auto-attivazione, in cui una proteina può agire sia come repressore di altri geni sia come attivatore del proprio gene. Vengono testate due distinte interpretazioni biologiche del controllo dell'effettore:
  • Caso I: Il legame dell'effettore rende la proteina inattiva indipendentemente dalla sua funzione (repressore o attivatore).
  • Caso II: Il legame dell'effettore altera direttamente la conformazione della proteina, commutando la sua funzione da repressore attivo ad attivatore attivo.
• Modalità di legame: Vengono analizzati sia il legame competitivo (legame repressore/attivatore mutuamente esclusivo) che il legame non esclusivo (legame simultaneo).
• Analisi: Gli autori eseguono un'analisi di biforcazione, un'analisi di stabilità lineare (tramite gli autovalori della Jacobiana) e simulazioni numeriche per mappare i punti fissi, i punti di sella e i regimi dinamici (monostabilità, bistabilità, tristabilità, quadristabilità) in funzione della concentrazione dell'induttore e dei parametri del modello.

Contributi chiave e risultati

1. Spostamenti dinamici guidati dall'induttore: Nel modello di base dell'interruttore a tre geni, l'aumento della concentrazione di un singolo induttore (che bersaglia un repressore) guida il sistema attraverso regimi dinamici distinti. A basse concentrazioni di induttore, il sistema è tristabile (tre stati stabili, ciascuno dominato da un repressore). All'aumentare della concentrazione dell'induttore, il sistema subisce biforcazioni nodo-sella, passando prima a un regime tristabile ristretto (meno percorsi accessibili tra gli stati) e infine collassando in un interruttore bistabile tra i due repressori non indotti.
2. Sensibilità alla cooperatività e alla simmetria: La concentrazione soglia dell'induttore necessaria per passare dalla tristabilità alla bistabilità è relativamente insensibile all'aumento del coefficiente di Hill (cooperatività) una volta che $n \geq 4$. Tuttavia, il sistema è altamente sensibile alla rottura della simmetria. Piccole perturbazioni nelle forze di repressione ($\alpha$) possono eliminare interamente il regime tristabile o restringere la finestra per dinamiche complesse, suggerendo che la "vicinanza" di un sistema biologico alla simmetria può essere inferita dalla sua risposta dinamica all'induzione.
3. Impatto dell'auto-attivazione: L'introduzione dell'auto-attivazione non cambia necessariamente il numero di stati stabili nel regime di legame competitivo, ma approfondisce i pozzi di potenziale degli stati stabili, rendendo le transizioni tra di essi meno probabili. Tuttavia, con un'auto-attivazione sufficientemente forte, il sistema può esibire la quadristabilità (quattro stati stabili), incluso uno stato in cui tutti e tre i geni sono espressi.
4. Dinamiche distinte dalle meccanismi dell'effettore: Lo studio dimostra che l'interpretazione biologica di come gli effettori controllano i fattori di trascrizione porta a paesaggi dinamici fondamentalmente diversi:
   • Nel Caso I (l'effettore inattiva la proteina), la risposta del sistema è principalmente una soppressione dell'attività.
   • Nel Caso II (l'effettore commuta la funzione), il sistema esibisce diverse soglie di biforcazione e configurazioni di stati stabili. Ad esempio, negli scenari di legame non esclusivo, il Caso II permette stati stabili che esprimono il gene bersaglio dell'effettore solo ad alte concentrazioni di effettore, mentre il Caso I lo permette a basse concentrazioni.
5. Complessità del legame non esclusivo: Quando repressori e attivatori legano in modo non esclusivo, lo spazio dei parametri si espande per consentire stati stabili con più geni espressi simultaneamente. L'esistenza di questi stati dipende fortemente dall'equilibrio tra forza di attivazione e forza di repressione, nonché dal meccanismo specifico di controllo dell'effettore.

Significatività
L'articolo sostiene che modellare il controllo dell'effettore attraverso un'interpretazione termodinamica concreta (modello MWC) trasforma non solo le tendenze fenotipiche della sintonizzazione dinamica, ma anche l'insieme dei regimi dinamici disponibili. Collegando l'analisi di biforcazione direttamente ai parametri sperimentalmente accessibili (concentrazioni di effettore), lo studio offre un quadro per:

• Decifrare i meccanismi regolatori: Distinguere tra diversi meccanismi biologici (ad esempio, se un effettore inattiva una proteina o ne commuta la funzione) basandosi sulla risposta input-output del sistema e sul comportamento di biforcazione.
• Sintonizzazione sperimentale: Fornire una base teorica per l'uso delle concentrazioni di induttore come "manopole di sintonizzazione" per sondare la topologia sottostante delle reti di regolazione genica e inferire parametri come la cooperatività e la simmetria.
• Comprendere la multistabilità: Illuminare come semplici cambiamenti nell'architettura regolatoria (aggiunta di auto-attivazione, cambiamento della esclusività del legame) e nei meccanismi dell'effettore possano generare comportamenti multistabili complessi (tristabilità, quadristabilità) rilevanti per le decisioni del destino cellulare in sistemi naturali e sintetici.

Gli autori concludono che il loro approccio integra il lavoro esistente definendo rigorosamente come funzionano gli effettori, consentendo così l'identificazione dei parametri fisici che guidano l'espressione nei dati sperimentali e offrendo una struttura sintonizzabile per codificare il comportamento multistabile inducibile.