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

Technische Samenvatting: Laagfrequente output-fluctuaties in een open exclusieproces met deeltjespauzes

Probleemstelling
Gedreven transport in interagerende deeltjessystemen, zoals het Totally Asymmetric Simple Exclusion Process (TASEP), zet lokale stochastische gebeurtenissen om in collectieve fluctuaties. Terwijl de standaard TASEP het verkeer van ribosomen en polymerasen beschrijft, bevatten biologische systemen vaak trage interne toestanden waarbij deeltjes reversibel schakelen tussen actieve en gepauzeerde toestanden (bijv. transcriptionele pauzes). Vorig werk heeft vastgesteld dat reversibele pauzes de stationaire stroom verminderen en een eindige-grootte-regime creëren waarin het gemiddelde aantal gepauzeerde deeltjes ($N_p$) de relevante controlevariabele wordt, met name wanneer $N_p$ van de orde grootte één is. Echter, de temporele statistieken van de output (voltooide exits) in dit regime blijven onopgelost. Specifiek is het onduidelijk of het eindige-grootte pauze-regime, geïdentificeerd door de gemiddelde stroom, een onderscheidende signatuur bezit in de stationaire statistieken van voltooide exits, en hoe trage interne toestanden de temporele organisatie van de output hervormen voorbij het louter verlagen van de gemiddelde doorvoer.

Methodologie
De studie maakt gebruik van een eindig open TASEP op een rooster van lengte $L$ ($50 \le L \le 500$). Deeltjes springen unidirectioneel met snelheid $\epsilon$, komen binnen met snelheid $\alpha$ en verlaten het systeem met snelheid $\beta$. Cruciaal is dat deeltjes reversibel schakelen tussen een actieve staat ($A$) en een gepauzeerde staat ($P$) met snelheden $k_p$ (pauzeren) en $k_u$ (ontpauzeren). Een gepauzeerd deeltje op de exit-locatie moet eerst weer actief worden voordat het kan vertrekken.

De dynamica wordt gesimuleerd met een event-driven kinetisch Monte Carlo-algoritme. De primaire observabelen zijn:

1. Output Telling Statistiek: De cumulatieve hoeveelheid exits $N(t)$ wordt geanalyseerd over observatievensters van duur $T_w$. De studie richt zich op de venster-afhankelijke Fano-factor, $F(T_w) = \text{Var}(N_w)/\langle N_w \rangle$, om kortstondige tel-ruis te onderscheiden van trage output-modulatie.
2. Structurele Observabelen: Het aantal gepauzeerde deeltjes ($N_p = L\rho_{\text{paused}}$), de grootte van de grootste aaneengesloten deeltjescluster ($C_{\text{max}}$), en verblijftijden van specifie-ke toestanden (bijv. pauze-vrij versus pauze-bevattend, of zwak versus sterk gecongestioneerd).
3. Analytische Benadering: Een minimaal "constant-birth, linear-death" model wordt geconstrueerd om de pauze-aantal dynamica te benaderen als een immigratie-death proces, met als doel de locatie van de ruis-maximum te rationaliseren.

Belangrijkste Bijdragen en Resultaten

• Niet-monotone Laagfrequente Ruis: Terwijl het verhogen van de gepauzeerde populatie ($N_p$) de gemiddelde stationaire stroom ($J$) geleidelijk verlaagt, vertoont de laagfrequente Fano-factor ($F_\infty$) een sterke niet-monotone gedraging. $F_\infty$ bereikt een uitgesproken maximum wanneer het gemiddelde aantal gepauzeerde deeltjes van de orde grootte één is ($N_p \simeq 1.5\text{--}2$). Dit maximum blijft persistent over roosterlengtes van $L=50$ tot $500$.
• Eindige-grootte Schaling: De microscopische pauze-snelheid ($k_p$) die nodig is om dit ruis-maximum te bereiken, schaalt als $k_p^{\text{max}} \propto L^{-1}$. Deze schaling komt voort uit het feit dat een langer rooster meer deeltjes bevat die blootstaan aan pauzes; aldus is een kleinere per-deeltje snelheid voldoende om de kritieke orde-één aantallen gepauzeerde deeltjes ($N_p^\star$) te handhaven.
• Mechanisme van Fluctuatie: Het maximale ruisniveau treedt niet op in het limiet van zeldzame pauzes (regelmatige output) noch in het limiet van aanhoudende congestie (onderdrukte maar constante output). In plaats daarvan ontstaat het in een intermitterend regime waarin het systeem afwisselt tussen productieve intervallen en uitgebreide perioden van zwakke of verdwijnende output. Deze afwisseling wordt gedreven door de nucleatie en dissolutie van grote deeltjesclusters (jams) geïnitieerd door pauzes.
• Rol van de Ontpauzeer-snelheid ($k_u$): Hoewel $N_p$ de positie van de crossover organiseert, bepaalt het niet volledig de amplitude van de ruis. Bij een vaste $N_p$ verhoogt een tragere ontpauzeringssnelheid (kleinere $k_u$) zowel de correlatietijd ($\tau_c$) als de ruisamplitude ($F_\infty$). Dit geeft aan dat de ruisamplitude afhangt van de ratio tussen de pauze-levensduur en de tijdschalen van jam-vorming en relaxatie.
• Structurele Correlaten:
  • Correlatietijd: De verblijfstijd-schaal geassocieerd met het verschijnen en verdwijnen van gepauzeerde deeltjes ($N_{\text{paused}}=0$ versus $\ge 1$) volgt de gefitte output-correlatietijd ($\tau_c$) nauwlettend. In tegen testelling hiermee faalt een drempelgebaseerde definitie van congestie (gebaseerd op $C_{\text{max}}$) om $\tau_c$ te volgen naarmate $N_p$ toeneemt.
  • Ruisamplitude: De amplitude van de laagfrequente ruis correleert sterk met de variantie van de grootste clustergrootte ($\text{Var}(C_{\text{max}})$), in plaats van de gemiddelde clustergrootte. Het maximale ruisniveau treedt op wanneer de dominante cluster een breed bereik aan groottes verkent, wat reflectie is van herhaalde groei, krimp en reorganisatie.
• Beperkingen van het Minimale Model: Een minimaal analytisch model dat uitgaat van onafhankelijke pauzes (Poisson-verdeling) voorspelt een maximum bij $N_p^\star \simeq 1.50$. Simulaties tonen aan dat het werkelijke maximum optreedt nabij $N_p \simeq 1.9$ en dat het pauze-aantal overgedispergeerd is (variantie > gemiddelde). De afwijkingen ontstaan door ruimtelijke correlaties, randeffecten en door uitsluiting geïnduceerde congestie, die het minimale model verwaarloost.

Betekenis
Dit artikel toont aan dat reversibele pauzes in eindige open exclusieprocessen een distinct ruimere regime creëren waarin de timing van moleculaire productie zeer intermitterend is, zelfs als de gemiddelde doorvoer slechts licht wordt gereduceerd. De studie identificeert dat de meest significante output-ruis niet voortkomt uit de gemiddelde dichtheid van defecten, maar uit de collectieve dynamica van trage defecten die interageren met verkeersopstoppingen wanneer het aantal defecten klein is (orde grootte één).

Het werk maakt onderscheid tussen twee niveaus van organisatie:

1. De pauze-aantal dynamica bepaalt de tijdschaal voor output-decorrelatie.
2. De jam-reorganisatie dynamica (fluctuaties in de grootste cluster) bepaalt de amplitude van de ruis-modulatie.

Dit impliceert dat voor processen die gevoelig zijn voor de voltooiingstijden van transcripten (bijv. downstream genexpressie), de temporele levering van producten zeer variabel kan zijn door enkel elongatieverkeer, zonder dat hiervoor promoter-switching of sequentie-specifieke pauzeplaatsen vereist zijn. De bevindingen breiden het begrip van eindige-grootte effecten in gedreven roostergassen uit, door aan te tonen dat het "eindige-grootte-regime" niet alleen wordt gekenmerkt door een reductie in de stroom, maar door een specifieke, niet-monotone versterking van de laagfrequente output-ruis, gedreven door het samenspel van trage interne toestanden en collectieve jamming.