DNA Replication under Thermal, Chemical, and Genotoxic Stress

Resumo Técnico: Replicação de DNA sob Estresse Térmico, Químico e Genotóxico

Definição do Problema
A duplicação fiel do genoma é essencial para a integridade genômica, contudo, a replicação do DNA eucariótico deve permanecer robusta apesar das flutuações causadas por estresse térmico, químico e genotóxico. Embora o Saccharomyces cerevisiae (levedura de brotamento) seja um sistema modelo para o estudo da replicação, as visões convencionais frequentemente tratam seu programa de replicação como predominantemente determinístico, com tempos de disparo de origem fixos e velocidades de forquilha constantes. Evidências experimentais indicam uma estocasticidade substancial na eficiência do disparo de origens e uma heterogeneidade significativa nas velocidades das forquilhas de replicação. Modelos estocásticos existentes operam frequentemente em resoluções de grão grosso (ex: ~1 kbp) ou assumem velocidades médias de forquilha constantes, falhando em capturar a dinâmica de escala fina, eventos de replicação prolongados raros ou os mecanismos específicos pelos quais as células se adaptam a diversas condições de estresse. Há uma necessidade de um arcabouço quantitativo que integre o disparo estocástico de origens, distribuições de velocidade de forquilha e limitações de recursos em resolução de par de bases único para entender como a replicação robusta emerge sob estresse.

Metodologia
Os autores desenvolveram um arcabouço de Monte Carlo estocástico baseado em rede para a replicação de genoma completo em S. cerevisiae em resolução de par de bases único (1 bp). O modelo representa o genoma como uma fita unidimensional onde pares de bases transitam de um estado não replicado (0) para um estado replicado (1). O arcabouço é governado por dois parâmetros primários, motivados experimentalmente:

1. Disparo Estocástico de Origens: O modelo utiliza dados experimentalmente determinados para 459 potenciais origens, caracterizados por competência ($P$), tempo mediano de disparo ($t_{1/2}$) e largura do tempo de disparo ($t_w$). O disparo é modelado como uma distribuição de probabilidade Gaussiana dependente do tempo. Para evitar grandes lacunas entre origens, o genoma é particionado e as origens são amostradas estocasticamente com uma taxa constante.
2. Fator Limitante ($L(t)$): Um fator limitante dependente do tempo representa um pool finito de recursos celulares necessários para iniciar o disparo de origens. Uma origem dispara apenas se dois fatores limitantes estiverem disponíveis. Esses fatores são recrutados para formar forquilhas de replicação ativas e são liberados após a convergência (fusão) das forquilhas, permitindo que sejam reciclados para disparos subsequentes de origens. A dinâmica temporal de $L(t)$ segue uma função sigmoide restrita pelas estatísticas experimentais de forquilhas ativas.
3. Distribuição de Velocidade de Forquilha: Diferente de modelos que assumem uma velocidade constante, este arcabouço atribui velocidades de forquilha extraídas de uma distribuição Gaussiana assimétrica ($f(v)$), parametrizada por localização ($\xi$), largura ($\omega$) e assimetria ($\alpha$). Esses parâmetros são derivados de dados experimentais ou extrapolações lineares baseadas em condições de temperatura e estresse.

O modelo é testado contra perfis de replicação experimentais (contagem de forquilhas ativas, fração replicada, distâncias entre origens e perfis de replicação espaço-temporais) antes de ser aplicado a condições de estresse térmico, químico (ureia) e genotóxico.

Principais Contribuições

• Resolução de Par de Bases Único: O arcabouço alcança a resolução mais alta fisicamente significativa (1 bp), rastreando explicitamente o tempo de replicação e a propagação da forquilha, permitindo a análise de dinâmicas de escala fina e padrões de mutação localizados.
• Parametrização Mínima: O modelo reproduz diversas respostas de estresse usando apenas dois parâmetros efetivos: o fator limitante dependente do tempo e a distribuição de velocidade da forquilha, sem invocar mecanismos regulatórios específicos de cada condição.
• Integração da Heterogeneidade: Ele incorpora explicitamente a heterogeneidade da velocidade da forquilha e o disparo estocástico de origens, indo além das suposições determinísticas ou de velocidade constante.

Resultados

• Benchmarking: O modelo reproduz quantitativamente dados experimentais para números de forquilhas ativas, frações de replicação do genoma e perfis de replicação espaço-temporais (STRP) do cromossomo 14 a 23°C.
• Durações de Fase-S com Distribuição Erlang: A inclusão da heterogeneidade da velocidade da forquilha gera uma distribuição de tempos totais de replicação ($T_{rep}$) que segue uma distribuição Erlang com uma cauda pronunciada à direita. Isso explica a emergência de eventos de replicação anormalmente prolongados observados experimentalmente em linhagens de E. coli e células humanas, prevendo comportamento semelhante em S. cerevisiae. Simulações com velocidades de forquilha constantes falham em reproduzir esses eventos de longa duração.
• Estresse Térmico: O modelo prevê uma dependência não monotônica do número total de forquilhas de replicação em relação à temperatura. Os números de forquilhas atingem o pico em torno de 23°C; em temperaturas mais baixas, velocidades de forquilha mais lentas levam ao atraso na fusão e sequestro de fatores limitantes, suprimindo novos disparos de origens. Em temperaturas mais altas, velocidades de forquilha rápidas fazem com que muitas origens sejam replicadas passivamente antes do disparo. O modelo recapitula com precisão os valores experimentais de $T_{rep}$ através de temperaturas usando a extrapolação linear dos parâmetros da distribuição de velocidade da forquilha.
• Estresse Químico (Hidroxiureia): Sob estresse de hidroxiureia (HU), o modelo demonstra que a redução da velocidade da forquilha isoladamente é insuficiente para corresponder aos perfis de replicação experimentais. A reprodução precisa dos dados requer uma redução simultânea na disponibilidade do fator limitante (para ~30% de seu valor original). O modelo prevê uma relação de escala de lei de potência entre a velocidade da forquilha e a concentração de HU.
• Estresse Genotóxico: O arcabouço modela com sucesso a dinâmica de replicação em várias linhagens mutantes de S. cerevisiae (que afetam a acetilação de histonas, pools de dNTP e complexos de proteção de forquilha). Ele captura como mutações específicas alteram as distribuições de velocidade da forquilha e, em alguns casos (ex: rtt109$\Delta$), modulam a função do fator limitante, levando a mudanças no tempo total de replicação.

Significância
O artigo afirma que a replicação robusta do genoma através de diversas condições emerge de um arcabouço mínimo impulsionado pela competição por recursos limitantes de replicação e pela heterogeneidade das velocidades das forquilhas. O estudo fornece uma interpretação física de como a dinâmica global da replicação, incluindo a aceleração acentuada da progressão da fase-S e a distribuição das durações da fase-S, surge desses mecanismos fundamentais. Ao operar em resolução de par de bases único, o arcabouço oferece uma plataforma versátil para investigar dinâmicas de replicação em regiões genômicas locais associadas a mutações pontuais e doenças. Os autores concluem que este trabalho estabelece uma base quantitativa para entender como a replicação do genoma se adapta ao estresse, unificando o controle da replicação sob condições térmicas, químicas e genotóxicas através de um conjunto compartilhado de parâmetros efetivos.