Gene transcription by an RNA polymerase (RNAP) enzyme requires that double-stranded DNA be locally and transiently opened, which results in an increase of DNA supercoiling downstream of the RNAP and a decrease of supercoiling upstream of it. When the DNA is initially torsionally relaxed and the RNAP experiences sufficiently large rotational drag, these variations lead to positively supercoiled plectonemes ahead of the RNAPs and negatively supercoiled ones behind it, a feature known as “twin supercoiled domain” (TSD). This work aims at deciphering into some more detail the torsional dynamics of circular DNA molecules being transcribed by RNAP enzymes. To this end, we performed Brownian dynamics simulations with a specially designed coarse-grained model. Depending on the superhelical density of the DNA molecule and the ratio of RNAP’s twist injection rate and rotational relaxation speed, simulations reveal a rich panel of behaviors, which sometimes differ markedly from the crude TSD picture. In particular, for sufficiently slow rotational relaxation speed, positively supercoiled plectonemes never form ahead of an RNAP that transcribes a DNA molecule with physiological negative supercoiling. Rather, negatively supercoiled plectonemes form almost periodically at the upstream side of the RNAP and grow up to a certain length before detaching from the RNAP and destabilizing rapidly. The extent to which topological barriers hinder the dynamics of TSDs is also discussed.
Read full abstract