Low-frequency output fluctuations in an open exclusion process with particle pausing

Resumen Técnico: Fluctuaciones de baja frecuencia en la salida de un proceso de exclusión abierto con pausa de partículas

Planteamiento del problema
El transporte dirigido en sistemas de partículas interactuantes, como el TASEP (Proceso de Exclusión Totalmente Asimétrico), convierte eventos estocásticos locales en fluctuaciones colectivas. Mientras que el TASEP estándar describe el tráfico de ribosomas y polimerasas, los sistemas biológicos a menudo involucran estados internos lentos donde las partículas cambian reversiblemente entre estados activos y de pausa (por ejemplo, la pausa transcripcional). Trabajos previos establecieron que la pausa reversible reduce la corriente estacionaria y crea un régimen de tamaño finito donde el número medio de partículas en pausa ($N_p$) se convierte en la variable de control relevante, particularmente cuando $N_p$ es del orden de la unidad. Sin embargo, la estadística temporal de la salida (salidas completadas) en este régimen permanece sin resolver. Específicamente, no está claro si el régimen de pausa de tamaño finito identificado por la corriente media posee una firma distintiva en la estadística estacionaria de las salidas completadas, y cómo los estados internos lentos remodelan la organización temporal de la salida más allá de simplemente reducir el rendimiento medio.

Metodología
El estudio emplea un TASEP abierto finito en una red de longitud $L$ ($50 \le L \le 500$). Las partículas saltan unidireccionalmente con una tasa $\epsilon$, entran con una tasa $\alpha$ y salen con una tasa $\beta$. Crucialmente, las partículas cambian reversiblemente entre un estado activo ($A$) y un estado de pausa ($P$) con tasas $k_p$ (pausa) y $k_u$ (reanudación). Una partícula en pausa en el sitio de salida debe reanudar su actividad antes de salir.

La dinámica se simula mediante un algoritmo de Monte Carlo cinético basado en eventos. Los observables principales son:

1. Estadísticas de conteo de la salida: Se analiza el número acumulado de salidas $N(t)$ sobre ventanas de observación de duración $T_w$. El estudio se centra en el factor de Fano dependiente de la ventana, $F(T_w) = \text{Var}(N_w)/\langle N_w \rangle$, para distinguir el ruido de conteo de corto tiempo de la modulación lenta de la salida.
2. Observables estructurales: El número de partículas en pausa ($N_p = L\rho_{\text{pausada}}$), el tamaño del mayor clúster de partículas contiguas ($C_{\text{max}}$) y los tiempos de residencia de estados específicos (por ejemplo, sin pausa frente a con pausa, o con congestión débil frente a fuerte).
3. Aproximación analítica: Se construye un modelo mínimo de "nacimiento constante, muerte lineal" para aproximar la dinámica del número de pausas como un proceso de inmigración-muerte, con el objetivo de racionalizar la ubicación del máximo de ruido.

Contribuciones clave y resultados

• Ruido de baja frecuencia no monotónico: Mientras que el aumento de la población en pausa ($N_p$) reduce suavemente la corriente estacionaria media ($J$), el factor de Fano de baja frecuencia ($F_\infty$) exhibe un comportamiento fuertemente no monotónico. $F_\infty$ alcanza un máximo pronunciado cuando el número medio de partículas en pausa es del orden de la unidad ($N_p \simeq 1.5\text{--}2$). Este máximo persiste a través de longitudes de red $L=50$ a $500$.
• Escalamiento de tamaño finito: La tasa de pausa microscópica ($k_p$) requerida para alcanzar este máximo de ruido escala como $k_p^{\text{max}} \propto L^{-1}$. Este escalamiento surge porque una red más larga contiene más partículas expuestas a la pausa; por lo tanto, una tasa por partícula más pequeña es suficiente para mantener el número crítico de orden unidad de partículas en pausa ($N_p^\star$).
• Mecanismo de fluctuación: El máximo de ruido no ocurre en el límite de pausas raras (salida regular) ni en el límite de congestión persistente (salida suprimida pero constante). En cambio, emerge en un régimen intermitente donde el sistema alterna entre intervalos productivos y periodos extendidos de salida débil o nula. Esta alternancia es impulsada por la nucleación y disolución de grandes clústeres de partículas (atascos) iniciados por las pausas.
• Papel de la tasa de reanudación ($k_u$): Aunque $N_p$ organiza la posición del cruce, no determina completamente la amplitud del ruido. Para un $N_p$ fijo, una reanudación más lenta (menor $k_u$) aumenta tanto el tiempo de correlación ($\tau_c$) como la amplitud del ruido ($F_\infty$). Esto indica que la amplitud del ruido depende de la relación entre el tiempo de vida de la pausa y las escalas de tiempo de la formación y relajación de los atascos.
• Correlatos estructurales:
  • Tiempo de correlación: La escala de tiempo de residencia asociada con la aparición y desaparición de partículas en pausa ($N_{\text{pausadas}}=0$ frente a $\ge 1$) sigue de cerca al tiempo de correlación de salida ajustado ($\tau_c$). En contraste, una definición basada en umbrales de congestión (basada en $C_{\text{max}}$) no logra seguir a $\tau_c$ a medida que $N_p$ aumenta.
  • Amplitud del ruido: La amplitud del ruido de baja frecuencia se correlaciona fuertemente con la varianza del tamaño del mayor clúster ($\text{Var}(C_{\text{max}}$), en lugar del tamaño medio del clúster. El máximo de ruido ocurre cuando el clúster dominante explora un amplio rango de tamaños, reflejando el crecimiento, la reducción y la reorganización repetidos.
• Limitaciones de la teoría mínima: Un modelo analítico mínimo que asume pausas independientes (distribución de Poisson) predice un máximo en $N_p^\star \simeq 1.50$. Las simulaciones muestran que el máximo real ocurre cerca de $N_p \simeq 1.9$ y que el número de pausas está sobredispersado (varianza > media). Estas desviaciones surgen de las correlaciones espaciales, los efectos de frontera y la congestión inducida por la exclusión, que la teoría mínima ignora.

Significancia
El artículo demuestra que la pausa reversible en procesos de exclusión abiertos finitos crea un régimen de fluctuación distintivo donde la temporización de la producción molecular es altamente intermitente, incluso si el rendimiento medio es meramente reducido. El estudio identifica que el ruido de salida más significativo no proviene de la densidad media de defectos, sino de la dinámica colectiva de defectos lentos interactuando con atascos de tráfico cuando el número de defectos es pequeño (orden de la unidad).

El trabajo distingue entre dos niveles de organización:

1. La dinámica del número de pausas establece la escala de tiempo para la descorrelación de la salida.
2. La dinámica de reorganización de atascos (fluctuaciones en el mayor clúster) determina la amplitud de la modulación del ruido.

Esto implica que para procesos sensibles a los tiempos de completitud de la transcripción (por ejemplo, la expresión génica aguas abajo), la entrega temporal de productos puede ser altamente variable debido únicamente al tráfico de elongación, sin requerir el cambio de promotores o sitios de pausa específicos de la secuencia. Los hallazgos extienden la comprensión de los efectos de tamaño finito en gases de red impulsados, mostrando que el "régimen de tamaño finito" se caracteriza no solo por una reducción en la corriente, sino por un aumento específico y no monotónico del ruido de salida de baja frecuencia impulsado por la interacción entre estados internos lentos y el atasco colectivo.