Sequence-dependent kinetic model for transcription elongation by RNA polymerase

Transcription elongation of the plant RNA polymerase IV is prone to backtracking

RNA polymerase IV (Pol IV) forms a complex with RNA-directed RNA polymerase 2 (RDR2) to produce double-stranded RNA (dsRNA) precursors essential for plant gene silencing. In the "backtracking-triggered RNA channeling" model, Pol IV backtracks and delivers its transcript's 3' terminus to RDR2, which synthesizes dsRNA. However, the mechanisms underlying Pol IV backtracking and RNA protection from cleavage are unclear. Here, we determined cryo-electron microscopy structures of Pol IV elongation complexes at four states of its nucleotide addition cycle (NAC): posttranslocation, guanosine triphosphate-bound, pretranslocation, and backtracked states. The structures reveal that Pol IV maintains an open DNA cleft and kinked bridge helix in all NAC states, loosely interacts with the nucleoside triphosphate substrate, and barely contacts proximal backtracked nucleotides. Biochemical data indicate that Pol IV is inefficient in forward translocation and RNA cleavage. These findings suggest that Pol IV transcription elongation is prone to backtracking and incapable of RNA hydrolysis, ensuring efficient dsRNA production by Pol IV-RDR2.

Nanopore tweezers show fractional-nucleotide translocation in sequence-dependent pausing by RNA polymerase

RNA polymerases (RNAPs) carry out the first step in the central dogma of molecular biology by transcribing DNA into RNA. Despite their importance, much about how RNAPs work remains unclear, in part because the small (3.4 Angstrom) and fast (~40 ms/nt) steps during transcription were difficult to resolve. Here, we used high-resolution nanopore tweezers to observe the motion of single Escherichia coli RNAP molecules as it transcribes DNA ~1,000 times improved temporal resolution, resolving single-nucleotide and fractional-nucleotide steps of individual RNAPs at saturating nucleoside triphosphate concentrations. We analyzed RNAP during processive transcription elongation and sequence-dependent pausing at the yrbL elemental pause sequence. Each time RNAP encounters the yrbL elemental pause sequence, it rapidly interconverts between five translocational states, residing predominantly in a half-translocated state. The kinetics and force-dependence of this half-translocated state indicate it is a functional intermediate between pre- and post-translocated states. Using structural and kinetics data, we show that, in the half-translocated and post-translocated states, sequence-specific protein-DNA interaction occurs between RNAP and a guanine base at the downstream end of the transcription bubble (core recognition element). Kinetic data show that this interaction stabilizes the half-translocated and post-translocated states relative to the pre-translocated state. We develop a kinetic model for RNAP at the yrbL pause and discuss this in the context of key structural features.

Analytic delay distributions for a family of gene transcription models

Models intended to describe the time evolution of a gene network must somehow include transcription, the DNA-templated synthesis of RNA, and translation, the RNA-templated synthesis of proteins. In eukaryotes, the DNA template for transcription can be very long, often consisting of tens of thousands of nucleotides, and lengthy pauses may punctuate this process. Accordingly, transcription can last for many minutes, in some cases hours. There is a long history of introducing delays in gene expression models to take the transcription and translation times into account. Here we study a family of detailed transcription models that includes initiation, elongation, and termination reactions. We establish a framework for computing the distribution of transcription times, and work out these distributions for some typical cases. For elongation, a fixed delay is a good model provided elongation is fast compared to initiation and termination, and there are no sites where long pauses occur. The initiation and termination phases of the model then generate a nontrivial delay distribution, and elongation shifts this distribution by an amount corresponding to the elongation delay. When initiation and termination are relatively fast, the distribution of elongation times can be approximated by a Gaussian. A convolution of this Gaussian with the initiation and termination time distributions gives another analytic approximation to the transcription time distribution. If there are long pauses during elongation, because of the modularity of the family of models considered, the elongation phase can be partitioned into reactions generating a simple delay (elongation through regions where there are no long pauses), and reactions whose distribution of waiting times must be considered explicitly (initiation, termination, and motion through regions where long pauses are likely). In these cases, the distribution of transcription times again involves a nontrivial part and a shift due to fast elongation processes.

Synergy among Pausing, Intrinsic Proofreading, and Accessory Proteins Results in Optimal Transcription Speed and Tolerable Accuracy

Cleavage of dinucleotides after the misincorporational pauses serves as a proofreading mechanism that increases transcriptional elongation accuracy. The accuracy is further improved by accessory proteins such as GreA and TFIIS. However, it is not clear why RNAP pauses and why cleavage-factor-assisted proofreading is necessary despite transcriptional errors in vitro being of the same order as those in downstream translation. Here, we developed a chemical-kinetic model that incorporates most relevant features of transcriptional proofreading and uncovers how the balance between speed and accuracy is achieved. We found that long pauses are essential for high accuracy, whereas cleavage-factor-stimulated proofreading optimizes speed. Moreover, in comparison to the cleavage of a single nucleotide or three nucleotides, RNAP backtracking and dinucleotide cleavage improve both speed and accuracy. Our results thereby show how the molecular mechanism and the kinetic parameters of the transcriptional process were evolutionarily optimized to achieve maximal speed and tolerable accuracy.

An NTP-driven mechanism for the nucleotide addition cycle of Escherichia coli RNA polymerase during transcription

The elementary steps of transcription as catalyzed by E. coli RNA polymerase during one and two rounds of the nucleotide addition cycle (NAC) were resolved in rapid kinetic studies. Modelling of stopped-flow kinetic data of pyrophosphate release in a coupled enzyme assay during one round of the NAC indicates that the rate of pyrophosphate release is significantly less than that for nucleotide incorporation. Upon modelling of the stopped-flow kinetic data for pyrophosphate release during two rounds of the NAC, it was observed that the presence of the next nucleotide for incorporation increases the rate of release of the first pyrophosphate equivalent; incorrect nucleotides for incorporation had no effect on the rate of pyrophosphate release. Although the next nucleotide for incorporation increases the rate of pyrophosphate release, it is still significantly less than the rate of incorporation of the first nucleotide. The results from the stopped-flow kinetic studies were confirmed by using quench-flow followed by thin-layer chromatography (QF-TLC) with only the first nucleotide for incorporation labeled on the gamma phosphate with ³²P to monitor pyrophosphate release. Collectively, the results are consistent with an NTP-driven model for the NAC in which the binding of the next cognate nucleotide for incorporation causes a synergistic conformational change in the enzyme that triggers the more rapid release of pyrophosphate, translocation of the enzyme along the DNA template strand and nucleotide incorporation.

Thermodynamic model of bacterial transcription

Transcriptional pausing is highly regulated by the template DNA and nascent transcript sequences. Here, we propose a thermodynamic model of transcriptional pausing, based on the thermal energy of transcription bubbles and nascent RNA structures, to describe the kinetics of the reaction pathways between active translocation, elemental, backtracked, and hairpin-stabilized pauses. The model readily predicts experimentally detected pauses in high-resolution optical-tweezer measurements of transcription. Unlike other models, it also predicts the effect of tension and the GreA transcription factor on pausing.

Theoretical study of RNA-polymerase behavior considering the backtracking state

The dynamical behavior of the RNA polymerase in the transcription process is vital to gene expression. During the transcription process, the 3' end of the transcribed RNA can be dislocated from the active site of the enzyme and as a result, the RNA polymerase goes to the backtracked state. Here, we develop a theoretical model to study the transcription process considering the backtracking state. We aim at describing the behavior of the enzyme in the backtracking state in the presence of an external force, which leads to two possibilities: (i) rescuing from the backtracking state and, (ii) the arresting of the enzyme. We study the probability and the rate of the mentioned processes. In addition, we find that entering the backtracking state behaves like the Brownian ratchet mechanism. This model could shed some light on the modeling of the transcription process and further studies on the energy landscape of the backtracking channel and the gene regulation.

Accuracy and speed of elongation in a minimal model of DNA replication

Being a dual purpose enzyme, the DNA polymerase is responsible for elongation of the newly formed DNA strand as well as cleaving the erroneous growth in case of a misincorporation. The efficiency of replication depends on the coordination of the polymerization and exonuclease activity of DNA polymerase. Here, we propose and analyze a minimal kinetic model of DNA replication and determine exact expressions for the velocity of elongation and the accuracy of replication. We first analyze the case without exonuclease activity. In that case, accuracy is determined by a kinetic competition between stepping and unbinding, with discrimination between correct and incorrect nucleotides in both transitions. We then include exonuclease activity and ask how different modes of additional discrimination in the exonuclease pathway can improve the accuracy while limiting the detrimental effect of exonuclease on the speed of replication. In this way, we ask how the kinetic parameters of the model have to be set to coordinate the two activities of the enzyme for high accuracy and high speed. The analysis also shows that the design of a replication system does not universally have to follow the speed-accuracy trade-off rule, although it does in the biologically realized parameter range. The accuracy of the process is mainly controlled by the crucial role of stepping after erroneous incorporation, which has impact on both polymerase and exonuclease activities of DNA polymerase.

Transcriptional Pausing as a Mediator of Bacterial Gene Regulation

Cellular life depends on transcription of DNA by RNA polymerase to express genetic information. RNA polymerase has evolved not just to read information from DNA and write it to RNA but also to sense and process information from the cellular and extracellular environments. Much of this information processing occurs during transcript elongation, when transcriptional pausing enables regulatory decisions. Transcriptional pauses halt RNA polymerase in response to DNA and RNA sequences and structures at locations and times that help coordinate interactions with small molecules and transcription factors important for regulation. Four classes of transcriptional pause signals are now evident after decades of study: elemental pauses, backtrack pauses, hairpin-stabilized pauses, and regulator-stabilized pauses. In this review, I describe current understanding of the molecular mechanisms of these four classes of pause signals, remaining questions about how RNA polymerase responds to pause signals, and the many exciting directions now open to understand pausing and the regulation of transcript elongation on a genome-wide scale. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Residence time analysis of RNA polymerase transcription dynamics: A Bayesian sticky HMM approach

The time spent by a single RNA polymerase (RNAP) at specific locations along the DNA, termed "residence time," reports on the initiation, elongation, and termination stages of transcription. At the single-molecule level, this information can be obtained from dual ultrastable optical trapping experiments, revealing a transcriptional elongation of RNAP interspersed with residence times of variable duration. Successfully discriminating between long and short residence times was used by previous approaches to learn about RNAP's transcription elongation dynamics. Here, we propose an approach based on the Bayesian sticky hidden Markov model that treats all residence times for an Escherichia coli RNAP on an equal footing without a priori discriminating between long and short residence times. Furthermore, our method has two additional advantages: we provide full distributions around key point statistics and directly treat the sequence dependence of RNAP's elongation rate. By applying our approach to experimental data, we find assigned relative probabilities on long versus short residence times, force-dependent average residence time transcription elongation dynamics, ∼10% drop in the average backtracking durations in the presence of GreB, and ∼20% drop in the average residence time as a function of applied force in the presence of RNaseA.

Basic mechanisms and kinetics of pause-interspersed transcript elongation

RNA polymerase pausing during elongation is an important mechanism in the regulation of gene expression. Pausing along DNA templates is thought to be induced by distinct signals encoded in the nucleic acid sequence and halt elongation complexes to allow time for necessary co-transcriptional events. Pausing signals have been classified as those producing short-lived elemental, long-lived backtracked, or hairpin-stabilized pauses. In recent years, structural microbiology and single-molecule studies have significantly advanced our understanding of the paused states, but the dynamics of these states are still uncertain, although several models have been proposed to explain the experimentally observed pausing behaviors. This review summarizes present knowledge about the paused states, discusses key discrepancies among the kinetic models and their basic assumptions, and highlights the importance and challenges in constructing theoretical models that may further our biochemical understanding of transcriptional pausing.

GC-Content Dependence of Elastic and Overstretching Properties of DNA:RNA Hybrid Duplexes

DNA:RNA hybrid duplex plays important roles in various biological processes. Both its structural stability and interactions with proteins are highly sequence dependent. In this study, we utilize homebuilt optical tweezers to investigate how GC contents in the sequence influence the structural and mechanical properties of DNA:RNA hybrid by measuring its contour length, elasticities, and overstretching dynamics. Our results support that the DNA:RNA hybrid adopts a conformation between the A- and B-form helix, and the GC content does not affect its structural and elastic parameters obviously when varying from 40 to 60% before the overstretching transition of DNA:RNA hybrid occurs. In the overstretching transition, however, our study unravels significant heterogeneity and strong sequence dependence, suggesting the presence of a highly dynamic competition between the two processes, namely the S-form duplex formation (nonhysteretic) and the unpeeling (hysteretic). Analyzing the components left in DNA:RNA hybrid after the overstretching transition suggests that the RNA strand is more easily unpeeled than the DNA strand, whereas an increase in the GC content from 40 to 60% can significantly reduce unpeeling. Large hysteresis is observed between the stretching and relaxation processes, which is also quantitatively correlated with the percentage of unpeeling in the DNA:RNA duplex. Increasing in both the salt concentration and GC content can effectively reduce the hysteresis with the latter being more significant. Together, our study reveals that the mechanical properties of DNA:RNA hybrid duplexes are significantly different from double-stranded DNA and double-stranded RNA, and its overstretching behavior is highly sequence dependent. These results should be taken into account in the future studies on DNA:RNA-hybrid-related functional structures and motor proteins.

Template-directed growth of copolymers

The theory of multistate template-directed reversible copolymerization is developed by extending the method based on iterated function systems to matrices, taking into account the possibility of multiple activation states instead of a single one for the growth process. In this extended theory, the mean growth velocity is obtained with an iterated matrix function system and the probabilities of copolymer sequences are given by matrix products defined along the template. The theory allows us to understand the effects of template heterogeneity, which include a fractal distribution of local growth velocities far enough from equilibrium, and a regime of sublinear growth in time close to equilibrium.

Cooperative and sequence-dependent model for RNAP dynamics: Application to ribosomal gene transcription

Escherichia coli ribosomal genes are a well-established experimental model used to investigate the transcription process. These genes are essential to cell physiology and are therefore strongly expressed. Multiple transcription units collaborate in rrn expression. Experiments involving electron microscopy have shown the non-uniform density of the RNA polymerases transcribing these ribosomal operons. Here, we investigate RNAP collaborative transcription in E. coli ribosomal genes using a stochastic sequence-dependent model that included interactions among the RNAPs. We achieved results consistent with experimental data using a model with variable parametrization for genic and intergenic regions, compared with previous attempts that used uniform parameters for genic and intergenic regions. Our model also showed that cooperative behaviour reduced the dwell times in pause sites predicted by the single-round approach but induced a new pausing event at an upstream position. This work may stimulate new experimental research and provide other scenarios to test our model predictions.

Bayesian inference and comparison of stochastic transcription elongation models

Transcription elongation can be modelled as a three step process, involving polymerase translocation, NTP binding, and nucleotide incorporation into the nascent mRNA. This cycle of events can be simulated at the single-molecule level as a continuous-time Markov process using parameters derived from single-molecule experiments. Previously developed models differ in the way they are parameterised, and in their incorporation of partial equilibrium approximations. We have formulated a hierarchical network comprised of 12 sequence-dependent transcription elongation models. The simplest model has two parameters and assumes that both translocation and NTP binding can be modelled as equilibrium processes. The most complex model has six parameters makes no partial equilibrium assumptions. We systematically compared the ability of these models to explain published force-velocity data, using approximate Bayesian computation. This analysis was performed using data for the RNA polymerase complexes of E. coli, S. cerevisiae and Bacteriophage T7. Our analysis indicates that the polymerases differ significantly in their translocation rates, with the rates in T7 pol being fast compared to E. coli RNAP and S. cerevisiae pol II. Different models are applicable in different cases. We also show that all three RNA polymerases have an energetic preference for the posttranslocated state over the pretranslocated state. A Bayesian inference and model selection framework, like the one presented in this publication, should be routinely applicable to the interrogation of single-molecule datasets.

Transcription is a critical biological process which occurs in all living organisms. It involves copying the organism’s genetic material into messenger RNA (mRNA) which directs protein synthesis on the ribosome. Transcription is performed by RNA polymerases which have been extensively studied using both ensemble and single-molecule techniques. Single-molecule data provides unique insights into the molecular behaviour of RNA polymerases. Transcription at the single-molecule level can be computationally simulated as a continuous-time Markov process and the model outputs compared with experimental data. In this study we use Bayesian techniques to perform a systematic comparison of 12 stochastic models of transcriptional elongation. We demonstrate how equilibrium approximations can strengthen or weaken the model, and show how Bayesian techniques can identify necessary or unnecessary model parameters. We describe a framework to a) simulate, b) perform inference on, and c) compare models of transcription elongation.

Oxazinomycin arrests RNA polymerase at the polythymidine sequences

Oxazinomycin is a C-nucleoside antibiotic that is produced by Streptomyces hygroscopicus and closely resembles uridine. Here, we show that the oxazinomycin triphosphate is a good substrate for bacterial and eukaryotic RNA polymerases (RNAPs) and that a single incorporated oxazinomycin is rapidly extended by the next nucleotide. However, the incorporation of several successive oxazinomycins or a single oxazinomycin in a certain sequence context arrested a fraction of the transcribing RNAP. The addition of Gre RNA cleavage factors eliminated the transcriptional arrest at a single oxazinomycin and shortened the nascent RNAs arrested at the polythymidine sequences suggesting that the transcriptional arrest was caused by backtracking of RNAP along the DNA template. We further demonstrate that the ubiquitous C-nucleoside pseudouridine is also a good substrate for RNA polymerases in a triphosphorylated form but does not inhibit transcription of the polythymidine sequences. Our results collectively suggest that oxazinomycin functions as a Trojan horse substrate and its inhibitory effect is attributable to the oxygen atom in the position corresponding to carbon five of the uracil ring.

Mechanisms of Transcriptional Pausing in Bacteria

Pausing by RNA polymerase (RNAP) during transcription regulates gene expression in all domains of life. In this review, we recap the history of transcriptional pausing discovery, summarize advances in our understanding of the underlying causes of pausing since then, and describe new insights into the pausing mechanisms and pause modulation by transcription factors gained from structural and biochemical experiments. The accumulated evidence to date suggests that upon encountering a pause signal in the nucleic-acid sequence being transcribed, RNAP rearranges into an elemental, catalytically inactive conformer unable to load NTP substrate. The conformation, and as a consequence lifetime, of an elemental paused RNAP is modulated by backtracking, nascent RNA structure, binding of transcription regulators, or a combination of these mechanisms. We conclude the review by outlining open questions and directions for future research in the field of transcriptional pausing.

The Mechanisms of Substrate Selection, Catalysis, and Translocation by the Elongating RNA Polymerase

Multi-subunit DNA-dependent RNA polymerases synthesize all classes of cellular RNAs, ranging from short regulatory transcripts to gigantic messenger RNAs. RNA polymerase has to make each RNA product in just one try, even if it takes millions of successive nucleotide addition steps. During each step, RNA polymerase selects a correct substrate, adds it to a growing chain, and moves one nucleotide forward before repeating the cycle. However, RNA synthesis is anything but monotonous: RNA polymerase frequently pauses upon encountering mechanical, chemical and torsional barriers, sometimes stepping back and cleaving off nucleotides from the growing RNA chain. A picture in which these intermittent dynamics enable processive, accurate, and controllable RNA synthesis is emerging from complementary structural, biochemical, computational, and single-molecule studies. Here, we summarize our current understanding of the mechanism and regulation of the on-pathway transcription elongation. We review the details of substrate selection, catalysis, proofreading, and translocation, focusing on rate-limiting steps, structural elements that modulate them, and accessory proteins that appear to control RNA polymerase translocation.

C, Da LT, Yu J. A Viral T7 RNA Polymerase Ratcheting Along DNA With Fidelity Control

RNA polymerase (RNAP) from bacteriophage T7 is a representative single-subunit viral RNAP that can transcribe with high promoter activities without assistances from transcription factors. We accordingly studied this small transcription machine computationally as a model system to understand underlying mechanisms of mechano-chemical coupling and fidelity control in the RNAP transcription elongation. Here we summarize our computational work from several recent publications to demonstrate first how T7 RNAP translocates via Brownian alike motions along DNA right after the catalytic product release. Then we show how the backward translocation motions are prevented at post-translocation upon successful nucleotide incorporation, which is also subject to stepwise nucleotide selection and acts as a pawl for "selective ratcheting". The structural dynamics and energetics features revealed from our atomistic molecular dynamics (MD) simulations and related analyses on the single-subunit T7 RNAP thus provided detailed and quantitative characterizations on the Brownian-ratchet working scenario of a prototypical transcription machine with sophisticated nucleotide selectivity for fidelity control. The presented mechanisms can be more or less general for structurally similar viral or mitochondrial RNAPs and some of DNA polymerases, or even for the RNAP engine of the more complicated transcription machinery in higher organisms.

Transcription factor regulation of RNA polymerase's torque generation capacity

During transcription, RNA polymerase (RNAP) supercoils DNA as it translocates. The resulting torsional stress in DNA can accumulate and, in the absence of regulatory mechanisms, becomes a barrier to RNAP elongation, causing RNAP stalling, backtracking, and transcriptional arrest. Here we investigate whether and how a transcription factor may regulate both torque-induced Escherichia coli RNAP stalling and the torque generation capacity of RNAP. Using a unique real-time angular optical trapping assay, we found that RNAP working against a resisting torque was highly prone to extensive backtracking. We then investigated transcription in the presence of GreB, a transcription factor known to rescue RNAP from the backtracked state. We found that GreB greatly suppressed RNAP backtracking and remarkably increased the torque that RNAP was able to generate by 65%, from 11.2 pN⋅nm to 18.5 pN·nm. Variance analysis of the real-time positional trajectories of RNAP after a stall revealed the kinetic parameters of backtracking and GreB rescue. These results demonstrate that backtracking is the primary mechanism by which torsional stress limits transcription and that the transcription factor GreB effectively enhances the torsional capacity of RNAP. These findings suggest a broader role for transcription factors in regulating RNAP functionality and elongation.

The elemental mechanism of transcriptional pausing

Transcriptional pausing underlies regulation of cellular RNA biogenesis. A consensus pause sequence that acts on RNA polymerases (RNAPs) from bacteria to mammals halts RNAP in an elemental paused state from which longer-lived pauses can arise. Although the structural foundations of pauses prolonged by backtracking or nascent RNA hairpins are recognized, the fundamental mechanism of the elemental pause is less well-defined. Here we report a mechanistic dissection that establishes the elemental pause signal (i) is multipartite; (ii) causes a modest conformational shift that puts γ-proteobacterial RNAP in an off-pathway state in which template base loading but not RNA translocation is inhibited; and (iii) allows RNAP to enter pretranslocated and one-base-pair backtracked states easily even though the half-translocated state observed in paused cryo-EM structures rate-limits pause escape. Our findings provide a mechanistic basis for the elemental pause and a framework to understand how pausing is modulated by sequence, cellular conditions, and regulators.

Stationary RNA polymerase fluctuations during transcription elongation

We study fluctuation effects of nonsteric molecular interactions between RNA polymerase (RNAP) motors that move simultaneously on the same DNA track during transcription elongation. Based on a stochastic model that allows for the exact analytical computation of the stationary distribution of RNAPs as a function of their density, interaction strength, nucleoside triphosphate concentration, and rate of pyrophosphate release we predict an almost geometric headway distribution of subsequent RNAP transcribing on the same DNA segment. The localization length which characterizes the decay of the headway distribution depends directly only the average density of RNAP and the interaction strength, but not on specific single-RNAP properties. Density correlations are predicted to decay exponentially with the distance (in units of DNA base pairs), with a correlation length that is significantly shorter than the localization length.

Active site closure stabilizes the backtracked state of RNA polymerase

All cellular RNA polymerases (RNAP) occasionally backtrack along the template DNA as part of transcriptional proofreading and regulation. Here, we studied the mechanism of RNAP backtracking by one nucleotide using two complementary approaches that allowed us to precisely measure the occupancy and lifetime of the backtracked state. Our data show that the stability of the backtracked state is critically dependent on the closure of the RNAP active site by a mobile domain, the trigger loop (TL). The lifetime and occupancy of the backtracked state measurably decreased by substitutions of the TL residues that interact with the nucleoside triphosphate (NTP) substrate, whereas amino acid substitutions that stabilized the closed active site increased the lifetime and occupancy. These results suggest that the same conformer of the TL closes the active site during catalysis of nucleotide incorporation into the nascent RNA and backtracking by one nucleotide. In support of this hypothesis, we construct a model of the 1-nt backtracked complex with the closed active site and the backtracked nucleotide in the entry pore area known as the E-site. We further propose that 1-nt backtracking mimics the reversal of the NTP substrate loading into the RNAP active site during on-pathway elongation.

Molecular mechanisms of RNA polymerase II transcription elongation elucidated by kinetic network models

Transcription elongation cycle (TEC) of RNA polymerase II (Pol II) is a process of adding a nucleoside triphosphate to the growing messenger RNA chain. Due to the long timescale events in Pol II TEC, an advanced computational technique, such as Markov State Model (MSM), is needed to provide atomistic mechanism and reaction rates. The combination of MSM and experimental results can be used to build a kinetic network model (KNM) of the whole TEC. This review provides a brief protocol to build MSM and KNM of the whole TEC, along with the latest findings of MSM and other computational studies of Pol II TEC. Lastly, we offer a perspective on potentially using a sequence dependent KNM to predict genome-wide transcription error.

Mfd Dynamically Regulates Transcription via a Release and Catch-Up Mechanism

The bacterial Mfd ATPase is increasingly recognized as a general transcription factor that participates in the resolution of transcription conflicts with other processes/roadblocks. This function stems from Mfd's ability to preferentially act on stalled RNA polymerases (RNAPs). However, the mechanism underlying this preference and the subsequent coordination between Mfd and RNAP have remained elusive. Here, using a novel real-time translocase assay, we unexpectedly discovered that Mfd translocates autonomously on DNA. The speed and processivity of Mfd dictate a "release and catch-up" mechanism to efficiently patrol DNA for frequently stalled RNAPs. Furthermore, we showed that Mfd prevents RNAP backtracking or rescues a severely backtracked RNAP, allowing RNAP to overcome stronger obstacles. However, if an obstacle's resistance is excessive, Mfd dissociates the RNAP, clearing the DNA for other processes. These findings demonstrate a remarkably delicate coordination between Mfd and RNAP, allowing efficient targeting and recycling of Mfd and expedient conflict resolution.

Transcriptional accuracy modeling suggests two-step proofreading by RNA polymerase

We suggest a novel two-step proofreading mechanism with two sequential rounds of proofreading selection in mRNA transcription. It is based on the previous experimental observations that the proofreading RNA polymerase cleaves off transcript fragments of at least 2 nt and that transcript elongation after a nucleotide misincorporation is anomalously slow. Taking these results into account, we extend the description of the accuracy of template guided nucleotide selection beyond previous models of RNA polymerase-dependent DNA transcription. The model derives the accuracy of initial and proofreading base selection from experimentally estimated nearest-neighbor parameters. It is also used to estimate the small accuracy enhancement of polymerase revisiting of previous positions following transcript cleavage.

Iterated function systems for DNA replication

The kinetic equations of DNA replication are shown to be exactly solved in terms of iterated function systems, running along the template sequence and giving the statistical properties of the copy sequences, as well as the kinetic and thermodynamic properties of the replication process. With this method, different effects due to sequence heterogeneity can be studied, in particular, a transition between linear and sublinear growths in time of the copies, and a transition between continuous and fractal distributions of the local velocities of the DNA polymerase along the template. The method is applied to the human mitochondrial DNA polymerase γ without and with exonuclease proofreading.

Transcription elongation

This review is focused on recent progress in understanding how Escherichia coli RNAP polymerase translocates along the DNA template and the factors that affect this movement. We discuss the fundamental aspects of RNAP translocation, template signals that influence forward or backward movement, and host or phage factors that modulate translocation.

Sequence-specific thermodynamic properties of nucleic acids influence both transcriptional pausing and backtracking in yeast

RNA Polymerase II pauses and backtracks during transcription, with many consequences for gene expression and cellular physiology. Here, we show that the energy required to melt double-stranded nucleic acids in the transcription bubble predicts pausing in Saccharomyces cerevisiae far more accurately than nucleosome roadblocks do. In addition, the same energy difference also determines when the RNA polymerase backtracks instead of continuing to move forward. This data-driven model corroborates—in a genome wide and quantitative manner—previous evidence that sequence-dependent thermodynamic features of nucleic acids influence both transcriptional pausing and backtracking.

Template-Directed Copolymerization, Random Walks along Disordered Tracks, and Fractals

In biology, template-directed copolymerization is the fundamental mechanism responsible for the synthesis of DNA, RNA, and proteins. More than 50 years have passed since the discovery of DNA structure and its role in coding genetic information. Yet, the kinetics and thermodynamics of information processing in DNA replication, transcription, and translation remain poorly understood. Challenging issues are the facts that DNA or RNA sequences constitute disordered media for the motion of polymerases or ribosomes while errors occur in copying the template. Here, it is shown that these issues can be addressed and sequence heterogeneity effects can be quantitatively understood within a framework revealing universal aspects of information processing at the molecular scale. In steady growth regimes, the local velocities of polymerases or ribosomes along the template are distributed as the continuous or fractal invariant set of a so-called iterated function system, which determines the copying error probabilities. The growth may become sublinear in time with a scaling exponent that can also be deduced from the iterated function system.

Control of transcriptional pausing by biased thermal fluctuations on repetitive genomic sequences

In the process of transcription elongation, RNA polymerase (RNAP) pauses at highly nonrandom positions across genomic DNA, broadly regulating transcription; however, molecular mechanisms responsible for the recognition of such pausing positions remain poorly understood. Here, using a combination of statistical mechanical modeling and high-throughput sequencing and biochemical data, we evaluate the effect of thermal fluctuations on the regulation of RNAP pausing. We demonstrate that diffusive backtracking of RNAP, which is biased by repetitive DNA sequence elements, causes transcriptional pausing. This effect stems from the increased microscopic heterogeneity of an elongation complex, and thus is entropy-dominated. This report shows a linkage between repetitive sequence elements encoded in the genome and regulation of RNAP pausing driven by thermal fluctuations.

Stochastic Kinetics of Nascent RNA

The stochastic kinetics of transcription is typically inferred from the distribution of RNA numbers in individual cells. However, cellular RNA reflects additional processes downstream of transcription, hampering this analysis. In contrast, nascent (actively transcribed) RNA closely reflects the kinetics of transcription. We present a theoretical model for the stochastic kinetics of nascent RNA, which we solve to obtain the probability distribution of nascent RNA per gene. The model allows us to evaluate the kinetic parameters of transcription from single-cell measurements of nascent RNA. The model also predicts surprising discontinuities in the distribution of nascent RNA, a feature which we verify experimentally.

A Two-Way Street: Regulatory Interplay between RNA Polymerase and Nascent RNA Structure

The vectorial (5'-to-3' at varying velocity) synthesis of RNA by cellular RNA polymerases (RNAPs) creates a rugged kinetic landscape, demarcated by frequent, sometimes long-lived, pauses. In addition to myriad gene-regulatory roles, these pauses temporally and spatially program the co-transcriptional, hierarchical folding of biologically active RNAs. Conversely, these RNA structures, which form inside or near the RNA exit channel, interact with the polymerase and adjacent protein factors to influence RNA synthesis by modulating pausing, termination, antitermination, and slippage. Here, we review the evolutionary origin, mechanistic underpinnings, and regulatory consequences of this interplay between RNAP and nascent RNA structure. We categorize and rationalize the extensive linkage between the transcriptional machinery and its product, and provide a framework for future studies.

Cis-Antisense Transcription Gives Rise to Tunable Genetic Switch Behavior: A Mathematical Modeling Approach

Antisense transcription has been extensively recognized as a regulatory mechanism for gene expression across all kingdoms of life. Despite the broad importance and extensive experimental determination of cis-antisense transcription, relatively little is known about its role in controlling cellular switching responses. Growing evidence suggests the presence of non-coding cis-antisense RNAs that regulate gene expression via antisense interaction. Recent studies also indicate the role of transcriptional interference in regulating expression of neighboring genes due to traffic of RNA polymerases from adjacent promoter regions. Previous models investigate these mechanisms independently, however, little is understood about how cells utilize coupling of these mechanisms in advantageous ways that could also be used to design novel synthetic genetic devices. Here, we present a mathematical modeling framework for antisense transcription that combines the effects of both transcriptional interference and cis-antisense regulation. We demonstrate the tunability of transcriptional interference through various parameters, and that coupling of transcriptional interference with cis-antisense RNA interaction gives rise to hypersensitive switches in expression of both antisense genes. When implementing additional positive and negative feed-back loops from proteins encoded by these genes, the system response acquires a bistable behavior. Our model shows that combining these multiple-levels of regulation allows fine-tuning of system parameters to give rise to a highly tunable output, ranging from a simple-first order response to biologically complex higher-order response such as tunable bistable switch. We identify important parameters affecting the cellular switch response in order to provide the design principles for tunable gene expression using antisense transcription. This presents an important insight into functional role of antisense transcription and its importance towards design of synthetic biological switches.

Sequence-Dependent Promoter Escape Efficiency Is Strongly Influenced by Bias for the Pretranslocated State during Initial Transcription

Abortive transcription initiation can be rate-limiting for promoter escape and therefore represents a barrier to productive gene expression. The mechanism for abortive initiation is unknown, but the amount of abortive transcript is known to vary with the composition of the initial transcribed sequence (ITS). Here, we used a thermodynamic model of translocation combined with experimental validation to investigate the relationship between ITS and promoter escape on a set of phage T5 N25 promoters. We found a strong, negative correlation between RNAP's propensity to occupy the pretranslocated state during initial transcription and the efficiency of promoter escape (r = -0.67; p < 10(-6)). This correlation was almost entirely caused by free energy changes due to variation in the RNA 3' dinucleotide sequence at each step, implying that this sequence element controls the disposition of initial transcribing complexes. We tested our model experimentally by constructing a set of novel N25-ITS promoter variants; quantitative transcription analysis again showed a strong correlation (r = -0.81; p < 10(-6)). Our results support a model in which sequence-directed bias for the pretranslocated state during scrunching results in increased backtracking, which limits the efficiency of promoter escape. This provides an answer to the long-standing issue of how sequence composition of the ITS affects promoter escape efficiency.

Visualizing translocation dynamics and nascent transcript errors in paused RNA polymerases in vivo

BACKGROUND

Transcription elongation is frequently interrupted by pausing signals in DNA, with downstream effects on gene expression. Transcription errors also induce prolonged pausing, which can lead to a destabilized genome by interfering with DNA replication. Mechanisms of pausing associated with translocation blocks and misincorporation have been characterized in vitro, but not in vivo.

RESULTS

We investigate the pausing pattern of RNA polymerase (RNAP) in Escherichia coli by a novel approach, combining native elongating transcript sequencing (NET-seq) with RNase footprinting of the transcripts (RNET-seq). We reveal that the G-dC base pair at the 5' end of the RNA-DNA hybrid interferes with RNAP translocation. The distance between the 5' G-dC base pair and the 3' end of RNA fluctuates over a three-nucleotide width. Thus, the G-dC base pair can induce pausing in post-translocated, pre-translocated, and backtracked states of RNAP. Additionally, a CpG sequence of the template DNA strand spanning the active site of RNAP inhibits elongation and induces G-to-A errors, which leads to backtracking of RNAP. Gre factors efficiently proofread the errors and rescue the backtracked complexes. We also find that pausing events are enriched in the 5' untranslated region and antisense transcription of mRNA genes and are reduced in rRNA genes.

CONCLUSIONS

In E. coli, robust transcriptional pausing involves RNAP interaction with G-dC at the upstream end of the RNA-DNA hybrid, which interferes with translocation. CpG DNA sequences induce transcriptional pausing and G-to-A errors.

DNA template dependent accuracy variation of nucleotide selection in transcription

It has been commonly assumed that the effect of erroneous transcription of DNA genes into messenger RNAs on peptide sequence errors are masked by much more frequent errors of mRNA translation to protein. We present a theoretical model of transcriptional accuracy. It uses experimentally estimated standard free energies of double-stranded DNA and RNA/DNA hybrids and predicts a DNA template dependent transcriptional accuracy variation spanning several orders of magnitude. The model also identifies high-error as well a high-accuracy transcription motifs. The source of the large accuracy span is the context dependent variation of the stacking free energy of pairs of correct and incorrect base pairs in the ever moving transcription bubble. Our model predictions have direct experimental support from recent single molecule based identifications of transcriptional errors in the C. elegans transcriptome. Our conclusions challenge the general view that amino acid substitution errors in proteins are mainly caused by translational errors. It suggests instead that transcriptional error hotspots are the dominating source of peptide sequence errors in some DNA template contexts, while mRNA translation is the major cause of protein errors in other contexts.

Trigger loop folding determines transcription rate of Escherichia coli's RNA polymerase

Two components of the RNA polymerase (RNAP) catalytic center, the bridge helix and the trigger loop (TL), have been linked with changes in elongation rate and pausing. Here, single molecule experiments with the WT and two TL-tip mutants of the Escherichia coli enzyme reveal that tip mutations modulate RNAP's pause-free velocity, identifying TL conformational changes as one of two rate-determining steps in elongation. Consistent with this observation, we find a direct correlation between helix propensity of the modified amino acid and pause-free velocity. Moreover, nucleotide analogs affect transcription rate, suggesting that their binding energy also influences TL folding. A kinetic model in which elongation occurs in two steps, TL folding on nucleoside triphosphate (NTP) binding followed by NTP incorporation/pyrophosphate release, quantitatively accounts for these results. The TL plays no role in pause recovery remaining unfolded during a pause. This model suggests a finely tuned mechanism that balances transcription speed and fidelity.

Thermodynamic modeling of variations in the rate of RNA chain elongation of E. coli rrn operons

Previous electron-microscopic imaging has shown high RNA polymerase occupation densities in the 16S and 23S encoding regions and low occupation densities in the noncoding leader, spacer, and trailer regions of the rRNA (rrn) operons in E. coli. This indicates slower transcript elongation within the coding regions and faster elongation within the noncoding regions of the operon. Inactivation of four of the seven rrn operons increases the transcript initiation frequency at the promoters of the three intact operons and reduces the time for RNA polymerase to traverse the operon. We have used the DNA sequence-dependent standard free energy variation of the transcription complex to model the experimentally observed changes in the elongation rate along the rrnB operon. We also model the stimulation of the average transcription rate over the whole operon by increasing rate of transcript initiation. Monte Carlo simulations, taking into account initiation of transcription, translocation, and backward and forward tracking of RNA polymerase, partially reproduce the observed transcript elongation rate variations along the rrn operon and fully account for the increased average rate in response to increased frequency of transcript initiation.

Interactions between RNA polymerase and the "core recognition element" counteract pausing

Transcription elongation is interrupted by sequences that inhibit nucleotide addition and cause RNA polymerase (RNAP) to pause. Here, by use of native elongating transcript sequencing (NET-seq) and a variant of NET-seq that enables analysis of mutant RNAP derivatives in merodiploid cells (mNET-seq), we analyze transcriptional pausing genome-wide in vivo in Escherichia coli. We identify a consensus pause-inducing sequence element, G₋₁₀Y₋₁G(+1) (where -1 corresponds to the position of the RNA 3' end). We demonstrate that sequence-specific interactions between RNAP core enzyme and a core recognition element (CRE) that stabilize transcription initiation complexes also occur in transcription elongation complexes and facilitate pause read-through by stabilizing RNAP in a posttranslocated register. Our findings identify key sequence determinants of transcriptional pausing and establish that RNAP-CRE interactions modulate pausing.

Tagetitoxin inhibits transcription by stabilizing pre-translocated state of the elongation complex

Transcription elongation consists of repetition of the nucleotide addition cycle: phosphodiester bond formation, translocation and binding of the next nucleotide. Inhibitor of multi-subunit RNA polymerase tagetitoxin (TGT) enigmatically slows down addition of nucleotides in a sequence-dependent manner, only at certain positions of the template. Here, we show that TGT neither affects chemistry of RNA synthesis nor induces backward translocation, nor competes with the nucleoside triphosphate (NTP) in the active center. Instead, TGT increases the stability of the pre-translocated state of elongation complex, thus slowing down addition of the following nucleotide. We show that the extent of inhibition directly depends on the intrinsic stability of the pre-translocated state. The dependence of translocation equilibrium on the transcribed sequence results in a wide distribution (~1-10(3)-fold) of inhibitory effects of TGT at different positions of the template, thus explaining sequence-specificity of TGT action. We provide biochemical evidence that, in pre-translocated state, TGT stabilizes folded conformation of the Trigger Loop, which inhibits forward and backward translocation of the complex. The results suggest that Trigger Loop folding in the pre-translocated state may serve to reduce back-tracking of the elongation complex. Overall, we propose that translocation may be a limiting and highly regulated step of RNA synthesis.

External conditions inversely change the RNA polymerase II elongation rate and density in yeast

Elongation speed is a key parameter in RNA polymerase II (RNA pol II) activity. It affects the transcription rate, while it is conditioned by the physicochemical environment it works in at the same time. For instance, it is well-known that temperature affects the biochemical reactions rates. Therefore in free-living organisms that are able to grow at various environmental temperatures, such as the yeast Saccharomyces cerevisiae, evolution should have not only shaped the structural and functional properties of this key enzyme, but should have also provided mechanisms and pathways to adapt its activity to the optimal performance required. We studied the changes in RNA pol II elongation speed caused by alternations in growth temperature in yeast to find that they strictly follow the Arrhenius equation, and that they also provoke an almost inverse proportional change in RNA pol II density within the optimal growth temperature range (26-37 °C). Moreover, we discovered that yeast cells control the transcription initiation rate by changing the total amount of available RNA pol II.

On the move

Single-molecule experiments have shed new light on the mechanisms responsible for the movement of RNA polymerase along DNA during transcription.

Complete dissection of transcription elongation reveals slow translocation of RNA polymerase II in a linear ratchet mechanism

During transcription elongation, RNA polymerase has been assumed to attain equilibrium between pre- and post-translocated states rapidly relative to the subsequent catalysis. Under this assumption, recent single-molecule studies proposed a branched Brownian ratchet mechanism that necessitates a putative secondary nucleotide binding site on the enzyme. By challenging individual yeast RNA polymerase II with a nucleosomal barrier, we separately measured the forward and reverse translocation rates. Surprisingly, we found that the forward translocation rate is comparable to the catalysis rate. This finding reveals a linear, non-branched ratchet mechanism for the nucleotide addition cycle in which translocation is one of the rate-limiting steps. We further determined all the major on- and off-pathway kinetic parameters in the elongation cycle. The resulting translocation energy landscape shows that the off-pathway states are favored thermodynamically but not kinetically over the on-pathway states, conferring the enzyme its propensity to pause and furnishing the physical basis for transcriptional regulation. DOI:http://dx.doi.org/10.7554/eLife.00971.001.