Cell Division Changes Fate Decisions in a Genetic Toggle Switch

Resumo Técnico: Mudanças na Divisão Celular Alteram Decisões de Destino em um Interruptor Genético de Alternância

Problema
Redes de regulação gênica (GRNs) governam decisões de destino celular através de dinâmicas de multistabilidade, com o interruptor de alternância genética (genetic toggle switch) servindo como um modelo canônico para tal comportamento. Embora modelos clássicos de expressão gênica frequentemente incorporem os efeitos da divisão celular via um termo de "diluição efetiva" adicionado à taxa de degradação, esses modelos tratam a divisão como um processo contínuo em vez de um evento explícito. Em contraste, mecanismos de divisão explícita removem proteínas discretamente ao final de um ciclo celular. Embora estudos anteriores tenham observado que contabilizar totalmente a divisão pode reduzir a região de bistabilidade, permanece incerto se, dentro de regimes de parâmetros onde ambos os modelos exibem bistabilidade, condições iniciais idênticas levam ao mesmo estado celular estável. Este artigo aborda a lacuna na compreensão de como a modelagem explícita da divisão celular impacta a dinâmica e as previsões de destino de redes regulatórias multistáveis.

Metodologia
Para investigar analiticamente o impacto da divisão, os autores propõem uma versão simplificada e analiticamente tratável do interruptor de alternância: o Interruptor de Alternância Booleano (BTS). Este modelo é derivado ao tomar o limite de coeficientes de Hill infinitos nas equações diferenciais ordinárias (ODEs) acopladas padrão, substituindo as funções de Hill por funções degrau de Heaviside. Essa binarização da síntese proteica (LIGADO/DESLIGADO baseado em limiares de concentração) preserva as características dinâmicas fundamentais do interruptor de alternância, permitindo a derivação de expressões de forma fechada.

O estudo compara duas versões distintas do BTS:

1. BTS Padrão: Segue as convenções clássicas onde o volume celular é fixo ($V=1$). Os efeitos da divisão são incorporados apenas através de uma taxa de diluição efetiva ($\delta = \ln 2 / T$) adicionada às taxas de degradação ativa.
2. BTS com Divisão: Modela explicitamente o crescimento celular exponencial ($V(t) = V_0 e^{\delta t}$) e a divisão celular discreta. Nesta versão, as contagens de proteínas são reduzidas à metade ao final de cada ciclo celular ($t=T$), e não há termo de diluição efetiva nas equações diferenciais.

Os autores derivam expressões analíticas para as separatrizes (as variedades estáveis de equilíbrios de sela que dividem o plano de fase em bacias de atração) para ambos os modelos. Eles assumem degradação proteica ativa zero ($\beta_A = \beta_B = 0$) para focar nos efeitos da divisão, observando que as meias-vidas das proteínas em células de mamíferos frequentemente excedem a duração do ciclo celular. Adicionalmente, simulações estocásticas (500 execuções por condição inicial) são realizadas para validar a análise determinística contra o ruído intrínseco, particularmente em regimes de contagem baixa.

Contribuições Principais e Resultados
A principal contribuição do artigo é a derivação de separatrizes analíticas tanto para o modelo padrão quanto para o modelo BTS com inclusão de divisão, revelando que os dois modelos possuem fronteiras distintas para decisões de destino celular.

• Trajetórias Divergentes: A análise demonstra que a divisão celular pode redirecionar trajetórias com condições iniciais idênticas para estados estáveis opostos. Especificamente, a separatriz para o BTS com divisão não coincide com a do BTS padrão.
• Região de Desacordo (ROD): Os autores definem uma "Região de Desacordo" (ROD) no plano de fase onde os dois modelos preveem estados estáveis diferentes para a mesma contagem inicial de proteínas.
  • Quando os parâmetros de síntese e limiar satisfazem $\alpha_A \lambda_B > \alpha_B \lambda_A$, a ROD situa-se abaixo da separatriz padrão.
  • Quando $\alpha_A \lambda_B < \alpha_B \lambda_A$, a ROD situa-se à esquerda da separatriz padrão.
  • A área da ROD é estritamente positiva (exceto quando as razões entre síntese e limiar estão perfeitamente equilibradas) e é maximizada quando as razões entre as taxas de síntese de proteínas ou limiares de repressão são grandes.
• Validação Estocástica: Simulações estocásticas confirmam que a ROD determinística é robusta ao ruído intrínseco. Condições iniciais que caem dentro da ROD determinística resultam em uma alta proporção de simulações estocásticas onde os estados finais dos dois modelos discordam.

Significância e Alegações
O artigo afirma que a divisão celular pode remodelar fundamentalmente as fronteiras de destino em redes regulatórias multistáveis. Ao negligenciar o mecanismo explícito de divisão e confiar apenas na diluição efetiva, os modelos podem prever incorretamente os destinos celulares para um subconjunto significativo de condições iniciais (a ROD). Os autores concluem que essas discrepâncias não são meros artefatos de escolhas específicas de parâmetros, mas representam uma diferença fundamental em como as duas abordagens de modelagem lidam com a dinâmica de crescimento e divisão celular. Consequentemente, os resultados sugerem que discrepâncias semelhantes podem surgir em redes de regulação gênica multistáveis mais complexas, implicando que a modelagem explícita da divisão é crucial para previsões precisas da tomada de decisão celular.