Cell Division Changes Fate Decisions in a Genetic Toggle Switch

Resumen Técnico: Los cambios en la división celular determinan las decisiones de destino en un interruptor de conmutación genética

Planteamiento del problema
Las redes de regulación génica (GRN) gobiernan las decisiones de destino celular mediante la dinámica de la multistabilidad, siendo el interruptor de conmutación genética (genetic toggle switch) un modelo canónico para dicho comportamiento. Mientras que los modelos clásicos de expresión génica suelen incorporar los efectos de la división celular mediante un término de "dilución efectiva" añadido a la tasa de degradación, estos modelos tratan la división como un proceso continuo en lugar de un evento explícito. En contraste, los mecanismos de división explícitos eliminan las proteínas de forma discreta al final de un ciclo celular. Aunque estudios previos han señalado que tener plenamente en cuenta la división puede reducir la región de bistabilidad, aún no está claro si, dentro de regímenes de parámetros donde ambos modelos exhiben bistabilidad, condiciones iniciales idénticas conducen al mismo estado celular estable. Este artículo aborda la brecha en la comprensión de cómo el modelado explícito de la división celular impacta la dinámica y las predicciones de destino de las redes reguladoras multistables.

Metodología
Para investigar analíticamente el impacto de la división, los autores proponen una versión simplificada y analíticamente tratable del interruptor de conmutación: el Interruptor de Conmutación Booleano (BTS, por sus siglas en inglés). Este modelo se deriva tomando el límite de coeficientes de Hill infinitos en las ecuaciones diferenciales ordinarias (ODE) acopladas estándar, reemplazando las funciones de Hill por funciones escalón de Heaviside. Esta binarización de la síntesis proteica (encendido/apagado basado en umbrales de concentración) preserva las características dinámicas fundamentales del interruptor de conmutación al tiempo que permite la derivación de expresiones de forma cerrada.

El estudio compara dos versiones distintas del BTS:

1. BTS Estándar: Sigue las convenciones clásicas donde el volumen celular es fijo ($V=1$). Los efectos de la división se incorporan únicamente a través de una tasa de dilución efectiva ($\delta = \ln 2 / T$) añadida a las tasas de degradación activa.
2. BTS con División: Modela explícitamente el crecimiento exponencial celular ($V(t) = V_0 e^{\delta t}$) y la división celular discreta. En esta versión, los recuentos de proteínas se reducen a la mitad al final de cada ciclo celular ($t=T$), y no existe un término de dilución efectiva en las ecuaciones diferenciales.

Los autores derivan expresiones analíticas para las separatrices (las variedades estables de los equilibrios de silla que dividen el plano de fase en cuencas de atracción) para ambos modelos. Asumen una degradación proteica activa nula ($\beta_A = \beta_B = 0$) para centrarse en los efectos de la división, señalando que las vidas medias de las proteínas en células de mamíferos suelen exceder la duración del ciclo celular. Además, se realizan simulaciones estocásticas (500 ejecuciones por condición inicial) para validar el análisis determinista frente al ruido intrínseco, particularmente en regímenes de recuentos bajos.

Contribuciones Clave y Resultados
La contribución principal del artículo es la derivación de las separatrices analíticas tanto para el modelo BTS estándar como para el modelo con división, revelando que ambos modelos poseen fronteras distintas para las decisiones de destino celular.

• Trayectorias Divergentes: El análisis demuestra que la división celular puede redirigir trayectorias con condiciones iniciales idénticas hacia estados estables opuestos. Específicamente, la separatriz para el BTS con división no coincide con la del BTS estándar.
• Región de Desacuerdo (ROD): Los autores definen una "Región de Desacuerdo" (ROD) en el plano de fase donde los dos modelos predicen diferentes estados estables para los mismos recuentos iniciales de proteínas.
  • Cuando los parámetros de síntesis y umbral satisfacen $\alpha_A \lambda_B > \alpha_B \lambda_A$, la ROD se sitúa por debajo de la separatriz estándar.
  • Cuando $\alpha_A \lambda_B < \alpha_B \lambda_A$, la ROD se sitúa a la izquierda de la separatriz estándar.
  • El área de la ROD es estrictamente positiva (excepto cuando las razones entre las tasas de síntesis y los umbrales de represión están perfectamente equilibradas) y se maximiza cuando las razones entre las tasas de síntesis de proteínas o los umbrales de represión son grandes.
• Validación Estocástica: Las simulaciones estocásticas confirman que la ROD determinista es robusta al ruido intrínseco. Las condiciones iniciales que caen dentro de la ROD determinista resultan en una alta proporción de simulaciones estocásticas donde los estados finales de los dos modelos discrepan.

Significancia y Reivindicaciones
El artículo sostiene que la división celular puede remodelar fundamentalmente las fronteras de destino en las redes de regulación génica multistables. Al omitir el mecanismo explícito de división y confiar únicamente en la dilución efectiva, los modelos pueden predecir incorrectamente los destinos celulares para un subconjunto significativo de condiciones iniciales (la ROD). Los autores concluyen que estas discrepancias no son meros artefactos de elecciones de parámetros específicas, sino que representan una diferencia fundamental en cómo los dos enfoques de modelado manejan la dinámica del crecimiento y la división celular. Por consiguiente, los resultados sugieren que discrepancias similares pueden surgir en redes de regulación génica multistables más complejas, lo que implica que el modelado explícito de la división es crucial para predicciones precisas del proceso de toma de decisiones celulares.