Path-preference cellular-automaton model for traffic flow through transit points and its application to the transcription process in human cells

We all use path routing everyday as we take shortcuts to avoid traffic jams, or by using faster traffic means. Previous models of traffic flow of RNA polymerase II (RNAPII) during transcription, however, were restricted to one dimension along the DNA template. Here we report the modeling and application of traffic flow in transcription that allows preferential paths of different dimensions only restricted to visit some transit points, as previously introduced between the 5′ and 3′ end of the gene. According to its position, an RNAPII protein molecule prefers paths obeying two types of time-evolution rules. One is an asymmetric simple exclusion process (ASEP) along DNA, and the other is a three-dimensional jump between transit points in DNA where RNAPIIs are staying. Simulations based on our model, and comparison experimental results, reveal how RNAPII molecules are distributed at the DNA-loop-formation-related protein binding sites as well as CTCF insulator proteins (or exons). As time passes after the stimulation, the RNAPII density at these sites becomes higher. Apparent far-distance jumps in one dimension are realized by short-range three-dimensional jumps between DNA loops. We confirm the above conjecture by applying our model calculation to the SAMD4A gene by comparing the experimental results. Our probabilistic model provides possible scenarios for assembling RNAPII molecules into transcription factories, where RNAPII and related proteins cooperatively transcribe DNA.