Bifurcations and multistability in inducible three-gene toggle switch networks

Problema
Los sistemas biológicos dependen de redes de regulación génica capaces de presentar multistabilidad para impulsar distintos fenotipos celulares, tales como la diferenciación y las decisiones de destino del desarrollo. Aunque los interruptores de conmutación (toggle switches) bistables de dos genes están bien caracterizados, extender estos conceptos a redes de orden superior (por ejemplo, sistemas de tres genes) para comprender la tristabilidad y dinámicas más complejas sigue siendo un desafío. Una limitación principal en los modelos teóricos existentes es el tratamiento de los parámetros regulatorios (constantes de disociación, tasas de transcripción) como variables abstractas desconectadas de los mecanismos biológicos específicos que las controlan in vivo. Específicamente, el papel de las moléculas efectoras (inductores) en la modulación de la actividad de los factores de transcripción alostéricos suele omitirse o simplificarse. Este estudio aborda la brecha entre la topología de red abstracta y los mecanismos biológicos sintonizables experimentalmente mediante el examen de cómo las concentraciones de efectores impulsan cambios dinámicos en interruptores de conmutación de tres genes.

Metodología
Los autores construyen un marco matemático para una red de conmutación de tres genes de represión mutua (un interruptor de conmutación) utilizando modelos de la función de Hill derivados de principios de la mecánica estadística.

• Modelo Base: El sistema se define mediante ecuaciones diferenciales que describen la producción y degradación de tres represores ($R_1, R_2, R_3$). El modelo incorpora una interpretación termodinámica de la inducción utilizando el modelo de Monod-Wyman-Changeux (MWC). Aquí, la concentración total de un factor de transcripción se escala mediante una función de probabilidad, $p_{act}(c)$, que depende de la concentración de un efector alostérico ($c$). Esta probabilidad determina la fracción de represores activos.
• Simetría y Perturbación: El estudio asume inicialmente un sistema simétrico donde todos los represores se unen con la misma afinidad. Posteriormente, introduce sistemáticamente perturbaciones para romper esta simetría, analizando cómo los sesgos en las fuerzas de represión afectan los umbrales de bifurcación.
• Extensiones de Autoactivación: El modelo se extiende para incluir la autoactivación, donde una proteína puede actuar tanto como represor de otros genes como activador de sí misma. Se prueban dos interpretaciones biológicas distintas del control por efectores:
  • Caso I: La unión del efector vuelve a la proteína inactiva, independientemente de su función (represora o activadora).
  • Caso II: La unión del efector altera directamente la conformación de la proteína, cambiando su función de un represor activo a un activador activo.
• Modos de Unión: Se analizan tanto la unión competitiva (represor/activador mutuamente excluyentes) como la unión no exclusiva (unión simultánea).
• Análisis: Los autores realizan un análisis de bifurcación, un análisis de estabilidad lineal (vía autovalores de la matriz Jacobiana) y simulaciones numéricas para mapear puntos fijos, puntos de silla y regímenes dinámicos (monoestabilidad, bistabilidad, tristabilidad, cuadristabilidad) en función de la concentración del inductor y los parámetros del modelo.

Contribuciones Clave y Resultados

1. Cambios Dinámicos Impulsados por Inductores: En el interruptor de conmutación de tres genes base, el aumento de la concentración de un único inductor (dirigido a un represor) conduce al sistema a través de distintos regímenes dinámicos. A concentraciones bajas de inductor, el sistema es tristable (tres estados estables, cada uno dominado por un represor). A medida que la concentración del inductor aumenta, el sistema experimenta bifurcaciones de tipo nodo-silla, transitando primero hacia un régimen tristable restringido (menos caminos accesibles entre estados) y finalmente colapsando en un interruptor bistable entre los dos represores no inducidos.
2. Sensibilidad a la Cooperatividad y la Simetría: La concentración de inductor umbral requerida para cambiar de la tristabilidad a la bistabilidad es relativamente insensible al aumento del coeficiente de Hill (cooperatividad) una vez que $n \geq 4$. Sin embargo, el sistema es altamente sensible a la ruptura de la simetría. Pequeñas perturbaciones en las fuerzas de represión ($\alpha$) pueden eliminar el régimen tristable por completo o estrechar la ventana para dinámicas complejas, lo que sugiere que la "proximidad" de un sistema biológico a la simetría puede inferirse de su respuesta dinámica a la inducción.
3. Impacto de la Autoactivación: La introducción de la autoactivación no cambia necesariamente el número de estados estables en el régimen de unión competitiva, pero profundiza los pozos de potencial de los estados estables, haciendo que las transiciones entre ellos sean menos probables. No obstante, con una autoactivación suficientemente fuerte, el sistema puede exhibir cuadristabilidad (cuatro estados estables), incluyendo un estado donde los tres genes se expresan.
4. Dinámicas Distintas según los Mecanismos del Efector: El estudio demuestra que la interpretación biológica de cómo los efectores controlan los factores de transcripción conduce a paisajes dinámicos fundamentalmente diferentes:
   • En el Caso I (el efector inactiva la proteína), la respuesta del sistema es primordialmente una supresión de la actividad.
   • En el Caso II (el efector cambia la función), el sistema exhibe umbrales de bifurcación y configuraciones de estados estables diferentes. Por ejemplo, en escenarios de unión no exclusiva, el Caso II permite estados que expresan el gen dirigido por el efector solo a altas concentraciones de este, mientras que el Caso I permite esto a concentraciones bajas.
5. Complejidad de la Unión No Exclusiva: Cuando los represores y activadores se unen de forma no exclusiva, el espacio de parámetros se expande para permitir estados estables con múltiples genes expresados simultáneamente. La existencia de estos estados depende altamente del equilibrio entre la fuerza de activación y la fuerza de represión, así como del mecanismo específico de control del efector.

Significado
El artículo sostiene que modelar el control por efectores a través de una interpretación termodinámica concreta (modelo MWC) transforma no solo las tendencias fenotípicas del ajuste dinámico, sino también el conjunto de regímenes dinámicos disponibles. Al vincular el análisis de bifurcación directamente con parámetros experimentalmente accesibles (concentraciones de efectores), el estudio ofrece un marco para:

• Descifrar Mecanismos Regulatorios: Distinguir entre diferentes mecanismos biológicos (por ejemplo, si un efector inactiva una proteína o cambia su función) basándose en la respuesta entrada-salida del sistema y el comportamiento de bifurcación.
• Ajuste Experimental: Proporcionar una base teórica para utilizar las concentraciones de inductores como "perillas de ajuste" para sondear la topología subyacente de las redes de regulación génica e inferir parámetros como la cooperatividad y la simetría.
• Comprender la Multiestabilidad: Iluminar cómo cambios simples en la arquitectura regulatoria (añadir autoactivación, cambiar la exclusividad de unión) y los mecanismos de efectores pueden generar comportamientos de multistabilidad complejos (tristabilidad, cuadristabilidad) relevantes para las decisiones de destino celular en sistemas naturales y sintéticos.

Los autores concluyen que su enfoque complementa el trabajo existente al definir rigurosamente cómo funcionan los efectores, permitiendo así la identificación de parámetros físicos que impulsan la expresión en datos experimentales y ofreciendo una estructura ajustable para codificar el comportamiento multiestable inducible.