Sequence-dependent kinetic model for transcription elongation by RNA polymerase

Backtracking dynamics of RNA polymerase: pausing and error correction

Transcription by RNA polymerases is frequently interrupted by pauses. One mechanism of such pauses is backtracking, where the RNA polymerase translocates backward with respect to both the DNA template and the RNA transcript, without shortening the transcript. Backtracked RNA polymerases move in a diffusive fashion and can return to active transcription either by diffusive return to the position where backtracking was initiated or by cleaving the transcript. The latter process also provides a mechanism for proofreading. Here we present some exact results for a kinetic model of backtracking and analyse its impact on the speed and the accuracy of transcription. We show that proofreading through backtracking is different from the classical (Hopfield-Ninio) scheme of kinetic proofreading. Our analysis also suggests that, in addition to contributing to the accuracy of transcription, backtracking may have a second effect: it attenuates the slow down of transcription that arises as a side effect of discriminating between correct and incorrect nucleotides based on the stepping rates.

Control of transcriptional elongation

Elongation is becoming increasingly recognized as a critical step in eukaryotic transcriptional regulation. Although traditional genetic and biochemical studies have identified major players of transcriptional elongation, our understanding of the importance and roles of these factors is evolving rapidly through the recent advances in genome-wide and single-molecule technologies. Here, we focus on how elongation can modulate the transcriptional outcome through the rate-liming step of RNA polymerase II (Pol II) pausing near promoters and how the participating factors were identified. Among the factors we describe are the pausing factors--NELF (negative elongation factor) and DSIF (DRB sensitivity-inducing factor)--and P-TEFb (positive elongation factor b), which is the key player in pause release. We also describe the high-resolution view of Pol II pausing and propose nonexclusive models for how pausing is achieved. We then discuss Pol II elongation through the bodies of genes and the roles of FACT and SPT6, factors that allow Pol II to move through nucleosomes.

Kinetic competition between elongation rate and binding of NELF controls promoter-proximal pausing

Pausing of RNA polymerase II (Pol II) 20-60 bp downstream of transcription start sites is a major checkpoint during transcription in animal cells. Mechanisms that control pausing are largely unknown. We developed permanganate-ChIP-seq to evaluate the state of Pol II at promoters throughout the Drosophila genome, and a biochemical system that reconstitutes promoter-proximal pausing to define pausing mechanisms. Stable open complexes of Pol II are largely absent from the transcription start sites of most mRNA genes but are present at snRNA genes and the highly transcribed heat shock genes following their induction. The location of the pause is influenced by the timing between when NELF loads onto Pol II and how fast Pol II escapes the promoter region. Our biochemical analysis reveals that the sequence-specific transcription factor, GAF, orchestrates efficient pausing by recruiting NELF to promoters before transcription initiation and by assisting in loading NELF onto Pol II after initiation.

Transcription under torsion

In cells, RNA polymerase (RNAP) must transcribe supercoiled DNA, whose torsional state is constantly changing, but how RNAP deals with DNA supercoiling remains elusive. We report direct measurements of individual Escherichia coli RNAPs as they transcribed supercoiled DNA. We found that a resisting torque slowed RNAP and increased its pause frequency and duration. RNAP was able to generate 11 ± 4 piconewton-nanometers (mean ± standard deviation) of torque before stalling, an amount sufficient to melt DNA of arbitrary sequence and establish RNAP as a more potent torsional motor than previously known. A stalled RNAP was able to resume transcription upon torque relaxation, and transcribing RNAP was resilient to transient torque fluctuations. These results provide a quantitative framework for understanding how dynamic modification of DNA supercoiling regulates transcription.

Cooperative RNA polymerase molecules behavior on a stochastic sequence-dependent model for transcription elongation

The transcription process is crucial to life and the enzyme RNA polymerase (RNAP) is the major component of the transcription machinery. The development of single-molecule techniques, such as magnetic and optical tweezers, atomic-force microscopy and single-molecule fluorescence, increased our understanding of the transcription process and complements traditional biochemical studies. Based on these studies, theoretical models have been proposed to explain and predict the kinetics of the RNAP during the polymerization, highlighting the results achieved by models based on the thermodynamic stability of the transcription elongation complex. However, experiments showed that if more than one RNAP initiates from the same promoter, the transcription behavior slightly changes and new phenomenona are observed. We proposed and implemented a theoretical model that considers collisions between RNAPs and predicts their cooperative behavior during multi-round transcription generalizing the Bai et al. stochastic sequence-dependent model. In our approach, collisions between elongating enzymes modify their transcription rate values. We performed the simulations in Mathematica® and compared the results of the single and the multiple-molecule transcription with experimental results and other theoretical models. Our multi-round approach can recover several expected behaviors, showing that the transcription process for the studied sequences can be accelerated up to 48% when collisions are allowed: the dwell times on pause sites are reduced as well as the distance that the RNAPs backtracked from backtracking sites.

A traffic flow model for bio-polymerization processes

Bio-polymerization processes like transcription and translation are central to proper function of a cell. The speed at which the bio-polymer grows is affected both by the number of pauses of elongation machinery, as well the number of bio-polymers due to crowding effects. In order to quantify these effects in fast transcribing ribosome genes, we rigorously show that a classical traffic flow model is the limit of a mean occupancy ODE model. We compare the simulation of this model to a stochastic model and evaluate the combined effect of the polymerase density and the existence of pauses on the instantaneous transcription rate of ribosomal genes.

TEMPLATE-DIRECTED BIOPOLYMERIZATION: TAPE-COPYING TURING MACHINES

DNA, RNA and proteins are among the most important macromolecules in a living cell. These molecules are polymerized by molecular machines. These natural nano-machines polymerize such macromolecules, adding one monomer at a time, using another linear polymer as the corresponding template. The machine utilizes input chemical energy to move along the template which also serves as a track for the movements of the machine. In the Alan Turing year 2012, it is worth pointing out that these machines are "tape-copying Turing machines". We review the operational mechanisms of the polymerizer machines and their collective behavior from the perspective of statistical physics, emphasizing their common features in spite of the crucial differences in their biological functions. We also draw the attention of the physics community to another class of modular machines that carry out a different type of template-directed polymerization. We hope this review will inspire new kinetic models for these modular machines.

Nucleosomal elements that control the topography of the barrier to transcription

The nucleosome represents a mechanical barrier to transcription that operates as a general regulator of gene expression. We investigate how each nucleosomal component-the histone tails, the specific histone-DNA contacts, and the DNA sequence-contributes to the strength of the barrier. Removal of the tails favors progression of RNA polymerase II into the entry region of the nucleosome by locally increasing the wrapping-unwrapping rates of the DNA around histones. In contrast, point mutations that affect histone-DNA contacts at the dyad abolish the barrier to transcription in the central region by decreasing the local wrapping rate. Moreover, we show that the nucleosome amplifies sequence-dependent transcriptional pausing, an effect mediated through the structure of the nascent RNA. Each of these nucleosomal elements controls transcription elongation by affecting distinctly the density and duration of polymerase pauses, thus providing multiple and alternative mechanisms for control of gene expression by chromatin remodeling and transcription factors.

Intrinsic translocation barrier as an initial step in pausing by RNA polymerase II

Pausing of RNA polymerase II (RNAP II) by backtracking on DNA is a major regulatory mechanism in control of eukaryotic transcription. Backtracking occurs by extrusion of the 3' end of the RNA from the active center after bond formation and before translocation of RNAP II on DNA. In several documented cases, backtracking requires a special signal such as A/T-rich sequences forming an unstable RNA-DNA hybrid in the elongation complex. However, other sequence-dependent backtracking signals and conformations of RNAP II leading to backtracking remain unknown. Here, we demonstrate with S. cerevisiae RNAP II that a cleavage-deficient elongation factor TFIIS (TFIIS(AA)) enhances backtracked pauses during regular transcription. This is due to increased efficiency of formation of an intermediate that leads to backtracking. This intermediate may involve misalignment at the 3' end of the nascent RNA in the active center of the yeast RNAP II, and TFIIS(AA) promotes formation of this intermediate at the DNA sequences, presenting a high-energy barrier to translocation. We proposed a three-step mechanism for RNAP II pausing in which a prolonged dwell time in the pre-translocated state increases the likelihood of the 3' RNA end misalignment facilitating a backtrack pausing. These results demonstrate an important role of the intrinsic blocks to forward translocation in pausing by RNAP II.

Efficient reconstitution of transcription elongation complexes for single-molecule studies of eukaryotic RNA polymerase II

Single-molecule studies of RNA polymerase II (RNAP II) require high yields of transcription elongation complexes (TECs) with long DNA tethers upstream and downstream of the TEC. Here we report on a robust system to reconstitute both yeast and mammalian RNAP II with an efficiency of ~80% into TECs that elongate with an efficiency of ~90%, followed by rapid, high-efficiency tripartite ligation of long DNA fragments upstream and downstream of the reconstituted TECs. Single mammalian and yeast TECs reconstituted with this method have been successfully used in an optical-trapping transcription assay capable of applying forces that either assist or hinder transcript elongation.

Single-molecule studies of RNAPII elongation

Elongation, the transcriptional phase in which RNA polymerase (RNAP) moves processively along a DNA template, occurs via a fundamental enzymatic mechanism that is thought to be universally conserved among multi-subunit polymerases in all kingdoms of life. Beyond this basic mechanism, a multitude of processes are integrated into transcript elongation, among them fidelity control, gene regulatory interactions involving elongation factors, RNA splicing or processing factors, and regulatory mechanisms associated with chromatin structure. Many kinetic and molecular details of the mechanism of the nucleotide addition cycle and its regulation, however, remain elusive and generate continued interest and even controversy. Recently, single-molecule approaches have emerged as powerful tools for the study of transcription in eukaryotic organisms. Here, we review recent progress and discuss some of the unresolved questions and ongoing debates, while anticipating future developments in the field. This article is part of a Special Issue entitled: RNA Polymerase II Transcript Elongation.

A dynamic model for processive transcription elongation and backtracking long pauses by multisubunit RNA polymerases

RNA polymerases are enzymes that transcribe genes from DNA onto strands of RNA. The transcription elongation by multisubunit RNA polymerases is processive but nonuniform: one enzyme can translocate along the DNA template for thousands of nucleotide addition steps but, sometimes, it can enter backtracking long pauses. Here, we present a Brownian ratchet model for the processive transcription elongation and the backtracking long pauses, which is developed based on the available structural and biochemical studies. Using the model, we analytically study the dynamics of the transcription elongation, such as the effects of external force and NTP concentration on the transcription velocity free of pauses, and the dynamics of backtracking long pauses, such as the probabilities of entering and returning from the backtracking pauses, with the analytical results in good agreement with the available single molecule experimental data. Values of several parameters for both Escherichia coli and Saccharomyces cerevisiae RNA polymerases such as their affinities for the DNA/RNA substrate during transcription elongation are determined. Moreover, some testable predictions are presented.

Proofreading of misincorporated nucleotides in DNA transcription

The accuracy of DNA transcription is crucial for the proper functioning of the cell. Although RNA polymerases demonstrate selectivity for correct nucleotides, additional active mechanisms of transcriptional error correction are required to achieve observed levels of fidelity. Recent experimental findings have shed light on a particular mechanism of transcriptional error correction involving: (i) diffusive translocation of the RNA polymerase along the DNA (backtracking) and (ii) irreversible RNA cleavage. This mechanism achieves preferential cleavage of misincorporated nucleotides by biasing the local rates of translocation. Here, we study how misincorporated nucleotides affect backtracking dynamics and how this effect determines the level of transcriptional fidelity. We consider backtracking as a diffusive process in a periodic, one-dimensional energy landscape, which at a coarse-grained level gives rise to a hopping process between neighbouring local minima. We propose a model for how misincorporated nucleotides deform this energy landscape and hence affect the hopping rates. In particular, we show that this model can be used to derive both the theoretical limit on the fidelity (i.e. the minimum fraction of misincorporated nucleotides) and the actual fidelity relative to this optimum, achieved for specific combinations of the cleavage and polymerization rates. Finally, we study how external factors influencing backtracking dynamics affect transcriptional fidelity. We show that biologically relevant loads, similar to those exerted by nucleosomes or other transcriptional barriers, increase error correction.

Trigger loop dynamics mediate the balance between the transcriptional fidelity and speed of RNA polymerase II

During transcription, RNA polymerase II (RNAPII) must select the correct nucleotide, catalyze its addition to the growing RNA transcript, and move stepwise along the DNA until a gene is fully transcribed. In all kingdoms of life, transcription must be finely tuned to ensure an appropriate balance between fidelity and speed. Here, we used an optical-trapping assay with high spatiotemporal resolution to probe directly the motion of individual RNAPII molecules as they pass through each of the enzymatic steps of transcript elongation. We report direct evidence that the RNAPII trigger loop, an evolutionarily conserved protein subdomain, serves as a master regulator of transcription, affecting each of the three main phases of elongation, namely: substrate selection, translocation, and catalysis. Global fits to the force-velocity relationships of RNAPII and its trigger loop mutants support a Brownian ratchet model for elongation, where the incoming NTP is able to bind in either the pre- or posttranslocated state, and movement between these two states is governed by the trigger loop. Comparison of the kinetics of pausing by WT and mutant RNAPII under conditions that promote base misincorporation indicate that the trigger loop governs fidelity in substrate selection and mismatch recognition, and thereby controls aspects of both transcriptional accuracy and rate.

Complexity of RNA polymerase II elongation dynamics

Transcription of protein-coding genes by RNA polymerase II can be regulated at multiple points during the process of RNA synthesis, including initiation, elongation, and termination. In vivo data suggests that elongating polymerases exhibit heterogeneity throughout the gene body, suggestive of changes in elongation rate and/or pausing. Here, we review evidence from a variety of different experimental approaches for understanding regulation of transcription elongation. We compare steady-state measurements of nascent RNA density and polymerase occupancy to time-resolved measurements and point out areas of disagreement. Finally, we discuss future avenues of investigation for understanding this critically important step in gene regulation. This article is part of a Special Issue entitled: Chromatin in time and space.

A small post-translocation energy bias aids nucleotide selection in T7 RNA polymerase transcription

The RNA polymerase (RNAP) of bacteriophage T7 is a single subunit enzyme that can transcribe DNA to RNA in the absence of additional protein factors. In this work, we present a model of T7 RNAP translocation during elongation. Based on structural information and experimental data from single-molecule force measurements, we show that a small component of facilitated translocation or power stroke coexists with the Brownian-ratchet-driven motions, and plays a crucial role in nucleotide selection at pre-insertion. The facilitated translocation is carried out by the conserved Tyr(639) that moves its side chain into the active site, pushing aside the 3'-end of the RNA, and forming a locally stabilized post-translocation intermediate. Pre-insertion of an incoming nucleotide into this stabilized intermediate state ensures that Tyr(639) closely participates in selecting correct nucleotides. A similar translocation mechanism has been suggested for multi-subunit RNAPs involving the bridge-helix bending. Nevertheless, the bent bridge-helix sterically prohibits nucleotide binding in the post-transolocation intermediate analog; moreover, the analog is not stabilized unless an inhibitory protein factor binds to the enzyme. Using our scheme, we also compared the efficiencies of different strategies for nucleotide selection, and examined effects of facilitated translocation on forward tracking.

Applied force provides insight into transcriptional pausing and its modulation by transcription factor NusA

Transcriptional pausing by RNA polymerase (RNAP) plays an essential role in gene regulation. Pausing is modified by various elongation factors, including prokaryotic NusA, but the mechanisms underlying pausing and NusA function remain unclear. Alternative models for pausing invoke blockade events that precede translocation (on-pathway), enzyme backtracking (off-pathway), or isomerization to a nonbacktracked, elemental pause state (off-pathway). We employed an optical trapping assay to probe the motions of individual RNAP molecules transcribing a DNA template carrying tandem repeats encoding the his pause, subjecting these enzymes to controlled forces. NusA significantly decreased the pause-free elongation rate of RNAP while increasing the probability of entry into short- and long-lifetime pauses, in a manner equivalent to exerting a ~19 pN force opposing transcription. The effects of force and NusA on pause probabilities and lifetimes support a reaction scheme where nonbacktracked, elemental pauses branch off the elongation pathway from the pretranslocated state of RNAP.

Fitting experimental transcription data with a comprehensive template-dependent modular kinetic model

In the companion article, we developed a modular scheme for representing the kinetics of transcription elongation by RNA polymerase. As an example of how to use these approaches, in this article we use a comprehensive modular model of this sort to fit experimental transcript elongation results obtained on the canonical tR2 template of phage λ by means of complementary bulk gel electrophoresis and surface plasmon resonance assays. The gel electrophoresis results, obtained in experiments quenched at various times after initiation of transcription, provide distributions of RNA lengths as a function of time. The surface plasmon resonance methods were used to monitor increases and decreases in the total mass of transcription elongation complexes in the same experiments. The different measures of transcription dynamics that these methods provide allow us to use them in combination to obtain a set of largely robust and well-defined kinetic parameters. The results show that our modular approach can be used to develop and test predictive kinetic schemes that can be fit to real transcription elongation data. They also suggest that these approaches can be extended to simulate the kinetics of other processes that involve the processive extension or shortening of nucleic acid chains and related systems of sequential branching reaction events.

Development of a "modular" scheme to describe the kinetics of transcript elongation by RNA polymerase

Transcript elongation by RNA polymerase involves the sequential appearance of several alternative and off-pathway states of the transcript elongation complex (TEC), and this complicates modeling of the kinetics of the transcription elongation process. Based on solutions of the chemical master equation for such transcription systems as a function of time, we here develop a modular scheme for simulating such kinetic transcription data. This scheme deals explicitly with the problem of TEC desynchronization as transcript synthesis proceeds, and develops kinetic modules to permit the various alternative states of the TECs (paused states, backtracked states, arrested states, and terminated states) to be introduced one-by-one as needed. In this way, we can set up a comprehensive kinetic model of appropriate complexity to fit the known transcriptional properties of any given DNA template and set of experimental conditions, including regulatory cofactors. In the companion article, this modular scheme is successfully used to model kinetic transcription elongation data obtained by bulk-gel electrophoresis quenching procedures and real-time surface plasmon resonance methods from a template of known sequence that contains defined pause, stall, and termination sites.

Factor-independent transcription pausing caused by recognition of the RNA-DNA hybrid sequence

RNA polymerase pausing during transcription is implicated in controlling gene expression. This study identifies a new type of pausing mechanism, by which the RNAP core recognizes the shape of base pairs of the RNA–DNA hybrid, which determines the rate of translocation and the nucleotide addition cycle. The expression of a number of viral and bacterial genes is shown to be subject to this mechanism.

Pausing of transcription is an important step of regulation of gene expression in bacteria and eukaryotes. Here we uncover a factor-independent mechanism of transcription pausing, which is determined by the ability of the elongating RNA polymerase to recognize the sequence of the RNA–DNA hybrid. We show that, independently of thermodynamic stability of the elongation complex, RNA polymerase directly ‘senses' the shape and/or identity of base pairs of the RNA–DNA hybrid. Recognition of the RNA–DNA hybrid sequence delays translocation by RNA polymerase, and thus slows down the nucleotide addition cycle through ‘in pathway' mechanism. We show that this phenomenon is conserved among bacterial and eukaryotic RNA polymerases, and is involved in regulatory pauses, such as a pause regulating the production of virulence factors in some bacteria and a pause regulating transcription/replication of HIV-1. The results indicate that recognition of RNA–DNA hybrid sequence by multi-subunit RNA polymerases is involved in transcription regulation and may determine the overall rate of transcription elongation.

A systems view of the protein expression process

Many biological processes are regulated by changing the concentration and activity of proteins. The presence of a protein at a given subcellular location at a given time with a certain conformation is the result of an apparently sequential process. The rate of protein formation is influenced by chromatin state, and the rates of transcription, translation, and degradation. There is an exquisite control system where each stage of the process is controlled both by seemingly unregulated proteins as well as through feedbacks mediated by RNA and protein products. Here we review the biological facts and mathematical models for each stage of the protein production process. We conclude that advances in experimental techniques leading to a detailed description of the process have not been matched by mathematical models that represent the details of the process and facilitate analysis. Such an exercise is the first step towards development of a framework for a systems biology analysis of the protein production process.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1007/s11693-011-9088-1) contains supplementary material, which is available to authorized users.

RNA transcript 3'-proximal sequence affects translocation bias of RNA polymerase

Translocation of RNA polymerase on DNA is thought to involve oscillations between pretranslocated and posttranslocated states that are rectified by nucleotide addition or pyrophosphorolysis. The pretranslocated register is also a precursor to transcriptional pause states that mediate regulation of transcript elongation. However, the determinants of bias between the pretranslocated and posttranslocated states are incompletely understood. To investigate translocation bias in multisubunit RNA polymerases, we measured rates of pyrophosphorolysis, which occurs in the pretranslocated register, in minimal elongation complexes containing T. thermophilus or E. coli RNA polymerase. Our results suggest that the identity of RNA:DNA nucleotides in the active site are strong determinants of susceptibility to pyrophosphorolysis, and thus translocation bias, with the 3' RNA nucleotide favoring the pretranslocated state in the order U > C > A > G. The preference of 3' U vs G for the pretranslocated register appeared to be universal among both bacterial and eukaryotic RNA polymerases and was confirmed by exonuclease III footprinting of defined elongation complexes. However, the relationship of pyrophosphate concentration to the rate of pyrophosphorolysis of 3' U-containing versus 3' G-containing elongation complexes did not match predictions of a simple mechanism in which 3'-RNA seqeunce affects only translocation bias and pyrophosphate (PPi) binds only to the pretranslocated state.

A unified model of transcription elongation: what have we learned from single-molecule experiments?

The transcription of the genetic information encoded in DNA into RNA is performed by RNA polymerase (RNAP), a complex molecular motor, highly conserved across species. Despite remarkable progress in single-molecule techniques revealing important mechanistic details of transcription elongation (TE) with up to base-pair resolution, some of the results and interpretations of these studies are difficult to reconcile, and have not yet led to a minimal unified picture of transcription. We propose a simple model that accounts quantitatively for many of the experimental observations. This model belongs to the class of isothermal ratchet models of TE involving the thermally driven stochastic backward and forward motion (backtracking and forward tracking) of RNAP along DNA between single-nucleotide incorporation events. We uncover two essential features for the success of the model. The first is an intermediate state separating the productive elongation pathway from nonelongating backtracked states. The rates of entering and exiting this intermediate state modulate pausing by RNAP. The second crucial ingredient of the model is the cotranscriptional folding of the RNA transcript, sterically inhibiting the extent of backtracking. This model resolves several apparent differences between single-molecule studies and provides a framework for future work on TE.

Rapid pyrophosphate release from transcriptional elongation complexes appears to be coupled to a nucleotide-induced conformational change in E. coli core polymerase

In the nucleotide addition cycle, pyrophosphate is generated upon incorporation of each nucleotide. Rapid release of pyrophosphate is essential for facile transcription elongation. Stopped-flow kinetic studies involving alterations in the intrinsic protein fluorescence of the core polymerase upon the binding of pyrophosphate to well-defined elongation complexes (ECs) indicate that the intrinsic off-rate of pyrophosphate (k=5.7-8.1 s(-1)) is too slow to account for the rapid rate of nucleotide incorporation that occurs during processive transcription elongation. Stopped-flow kinetic studies on UTP binding followed by UMP incorporation into an EC as monitored by alterations in the intrinsic protein fluorescence of the core polymerase resulted in a set of first-order rate constants that varied in a hyperbolic manner as a function of UTP concentration. This is consistent with a binding step (K(UTP)=17±6 μM) followed by a conformational change (k=623±54 s(-1)) in the core polymerase. In comparable studies on ATP binding and AMP incorporation into an EC, the data were also consistent with a binding step (K(ATP)=44±6 μM) followed by a conformational change (k=411±51 s(-1)) in the core polymerase. In stopped-flow kinetic studies with α,β-methyleneadenosine 5' triphosphate, which can bind to the EC but cannot lead to nucleotide incorporation, the analysis of the hyperbolic dependence of the observed first-order rate constant on α,β-methyleneadenosine 5' triphosphate concentration yielded a value of 20±13 μM for the apparent dissociation constant and a value of 221±36 s(-1) for the first-order rate constant for the associated conformational change in the core polymerase. This indicates that the conformational change in the core polymerase precedes chemistry. In conjunction with previously reported results on the increase in the rate of pyrophosphate release in the presence of the next cognate nucleotide for incorporation, the data are consistent with a model in which rapid pyrophosphate release is coupled to a conformational change in the core polymerase that precedes chemistry and that occurs upon the binding of the next cognate nucleotide for incorporation.

Templated nucleoside triphosphate binding to a noncatalytic site on RNA polymerase regulates transcription

The regulation of RNA synthesis by RNA polymerase (RNAP) is essential for proper gene expression. Crystal structures of RNAP reveal two channels: the main channel that contains the downstream DNA and a secondary channel that leads directly to the catalytic site. Although nucleoside triphosphates (NTPs) have been seen only in the catalytic site and the secondary channel in these structures, several models of transcription elongation, based on biochemical studies, propose that template-dependent binding of NTPs in the main channel regulates RNA synthesis. These models, however, remain controversial. We used transient state kinetics and a mutant of RNAP to investigate the role of the main channel in regulating nucleotide incorporation. Our data indicate that a NTP specific for the i + 2 template position can bind to a noncatalytic site and increase the rate of RNA synthesis and that the NTP bound to this site can be shuttled directly into the catalytic site. We also identify fork loop 2, which lies across from the downstream DNA, as a functional component of this site. Taken together, our data support the existence of a noncatalytic template-specific NTP binding site in the main channel that is involved in the regulation of nucleotide incorporation. NTP binding to this site could promote high-fidelity processive synthesis under a variety of environmental conditions and allow DNA sequence-mediated regulatory signals to be communicated to the active site.

Comparison of Pause Predictions of Two Sequence-Dependent Transcription Models

Two recent theoretical models, Bai et al. (2004, 2007) and Tadigotla et al. (2006), formulated thermodynamic explanations of sequence-dependent transcription pausing by RNA polymerase (RNAP). The two models differ in some basic assumptions and therefore make different yet overlapping predictions for pause locations, and different predictions on pause kinetics and mechanisms. Here we present a comprehensive comparison of the two models. We show that while they have comparable predictive power of pause locations at low NTP concentrations, the Bai et al. model is more accurate than Tadigotla et al. at higher NTP concentrations. Pausing kinetics predicted by Bai et al. is also consistent with time-course transcription reactions, while Tadigotla et al. is unsuited for this type of kinetic prediction. More importantly, the two models in general predict different pausing mechanisms even for the same pausing sites, and the Bai et al. model provides an explanation more consistent with recent single molecule observations.

Synergistic action of RNA polymerases in overcoming the nucleosomal barrier

During gene expression, RNA polymerase (RNAP) encounters a major barrier at a nucleosome and yet must access the nucleosomal DNA. Previous in vivo evidence has suggested that multiple RNAPs might increase transcription efficiency through nucleosomes. Here we have quantitatively investigated this hypothesis using Escherichia coli RNAP as a model system by directly monitoring its location on the DNA via a single-molecule DNA-unzipping technique. When an RNAP encountered a nucleosome, it paused with a distinctive 10-base pair periodicity and backtracked by approximately 10-15 base pairs. When two RNAPs elongate in close proximity, the trailing RNAP apparently assists in the leading RNAP's elongation, reducing its backtracking and enhancing its transcription through a nucleosome by a factor of 5. Taken together, our data indicate that histone-DNA interactions dictate RNAP pausing behavior, and alleviation of nucleosome-induced backtracking by multiple polymerases may prove to be a mechanism for overcoming the nucleosomal barrier in vivo.

The origins of time-delay in template biopolymerization processes

Time-delays are common in many physical and biological systems and they give rise to complex dynamic phenomena. The elementary processes involved in template biopolymerization, such as mRNA and protein synthesis, introduce significant time delays. However, there is not currently a systematic mapping between the individual mechanistic parameters and the time delays in these networks. We present here the development of mathematical, time-delay models for protein translation, based on PDE models, which in turn are derived through systematic approximations of first-principles mechanistic models. Theoretical analysis suggests that the key features that determine the time-delays and the agreement between the time-delay and the mechanistic models are ribosome density and distribution, i.e., the number of ribosomes on the mRNA chain relative to their maximum and their distribution along the mRNA chain. Based on analytical considerations and on computational studies, we show that the steady-state and dynamic responses of the time-delay models are in excellent agreement with the detailed mechanistic models, under physiological conditions that correspond to uniform ribosome distribution and for ribosome density up to 70%. The methodology presented here can be used for the development of reduced time-delay models of mRNA synthesis and large genetic networks. The good agreement between the time-delay and the mechanistic models will allow us to use the reduced model and advanced computational methods from nonlinear dynamics in order to perform studies that are not practical using the large-scale mechanistic models.

Translocation by multi-subunit RNA polymerases

DNA template and RNA/DNA hybrid movement through RNA polymerase (RNAP) is referred to as "translocation". Because nucleic acid movement is coupled to NTP loading, pyrophosphate release, and conformational changes, the precise ordering of events during bond addition is consequential. Moreover, based on several lines of experimental evidence, translocation, pyrophosphate release or an associated conformational change may determine the transcription elongation rate. In this review we discuss various models of translocation, the data supporting the hypothesis that translocation rate determines transcription elongation rate and also data that may be inconsistent with this point of view. A model of the nucleotide addition cycle accommodating available experimental data is proposed. On the basis of this model, the molecular mechanisms regulating translocation and potential routes for NTP entry are discussed.

Structure and dynamics of a processive Brownian motor: the translating ribosome

There is mounting evidence indicating that protein synthesis is driven and regulated by mechanisms that direct stochastic, large-scale conformational fluctuations of the translational apparatus. This mechanistic paradigm implies that a free-energy landscape governs the conformational states that are accessible to and sampled by the translating ribosome. This scenario presents interdependent opportunities and challenges for structural and dynamic studies of protein synthesis. Indeed, the synergism between cryogenic electron microscopic and X-ray crystallographic structural studies, on the one hand, and single-molecule fluorescence resonance energy transfer (smFRET) dynamic studies, on the other, is emerging as a powerful means for investigating the complex free-energy landscape of the translating ribosome and uncovering the mechanisms that direct the stochastic conformational fluctuations of the translational machinery. In this review, we highlight the principal insights obtained from cryogenic electron microscopic, X-ray crystallographic, and smFRET studies of the elongation stage of protein synthesis and outline the emerging themes, questions, and challenges that lie ahead in mechanistic studies of translation.

Maximum rrn promoter activity in Escherichia coli at saturating concentrations of free RNA polymerase

During fast growth, the rrn P1 promoters of Escherichia coli operate at their maximum strength, but below their maximum activity (V(max)), since they are not saturated with RNA polymerase. Since higher concentrations of free RNA polymerase are expected to be found in strains carrying rrn deletions, we have analyzed reported electron micrographs of rrn operons from rrn deletion strains growing at maximal rates (at 37 degrees C) in LB medium [1]. We conclude that, in a strain with four of the seven rrn operons inactivated by partial deletions, transcripts are initiated at rrn P1 promoters 1.6-fold more rapidly than in the wild-type strain and the entirety of the rrn operon is transcribed at a 1.5-fold higher average elongation rate due to shortened pauses in the 16S and 23S regions. Under this condition, traffic congestion occurs in front of a pause site in the 5' leader region of the rrn operon near the beginning of the 16S gene; the congestion extends all the way back to the promoter, impedes promoter clearance and limits the promoter activity to one initiation per 0.56 s. This corresponds to a promoter activity of 107 transcripts/min and is assumed to be close to the V(max) of rrn P1 promoters.

From sequence to dynamics: the effects of transcription factor and polymerase concentration changes on activated and repressed promoters

Background

The fine tuning of two features of the bacterial regulatory machinery have been known to contribute to the diversity of gene expression within the same regulon: the sequence of Transcription Factor (TF) binding sites, and their location with respect to promoters. While variations of binding sequences modulate the strength of the interaction between the TF and its binding sites, the distance between binding sites and promoters alter the interaction between the TF and the RNA polymerase (RNAP).

Results

In this paper we estimated the dissociation constants (K_d) of several E. coli TFs in their interaction with variants of their binding sequences from the scores resulting from aligning them to Positional Weight Matrices. A correlation coefficient of 0.78 was obtained when pooling together sites for different TFs. The theoretically estimated K_d values were then used, together with the dissociation constants of the RNAP-promoter interaction to analyze activated and repressed promoters. The strength of repressor sites -- i.e., the strength of the interaction between TFs and their binding sites -- is slightly higher than that of activated sites. We explored how different factors such as the variation of binding sequences, the occurrence of more than one binding site, or different RNAP concentrations may influence the promoters' response to the variations of TF concentrations. We found that the occurrence of several regulatory sites bound by the same TF close to a promoter -- if they are bound by the TF in an independent manner -- changes the effect of TF concentrations on promoter occupancy, with respect to individual sites. We also found that the occupancy of a promoter will never be more than half if the RNAP concentration-to-K_p ratio is 1 and the promoter is subject to repression; or less than half if the promoter is subject to activation. If the ratio falls to 0.1, the upper limit of occupancy probability for repressed drops below 10%; a descent of the limits occurs also for activated promoters.

Conclusion

The number of regulatory sites may thus act as a versatility-producing device, in addition to serving as a source of robustness of the transcription machinery. Furthermore, our results show that the effects of TF concentration fluctuations on promoter occupancy are constrained by RNAP concentrations.

A look-ahead model for the elongation dynamics of transcription

This article introduces a chemical kinetic model of the transcriptional elongation dynamics of RNA polymerase. The model's novel concept is a look-ahead feature, in which nucleotides bind reversibly to the DNA before being incorporated covalently into the nascent RNA chain. Analytical and computational methods for studying the behavior of the look-ahead model are introduced, and several approaches to parameter estimation are tested on synthetic and also on actual experimental data. Two types of experimental data are considered: 1), the mean velocity of RNA polymerase as a function of the ambient concentrations of the ribonucleoside triphosphates; and 2), the distribution of time intervals between the forward steps of RNA polymerase. By separately fitting the look-ahead model to these two types of data, we obtain estimates of the model parameters. The most difficult parameter to estimate is the width of the look-ahead window. Both types of data suggest a small window size, but the second type does a better job of distinguishing the different window sizes. These latter data rule out a window size of 1, and they strongly suggest a look-ahead window that is approximately four bases in width. Additional experiments to determine the window size are proposed.

Active nucleosome displacement: a theoretical approach

Three-quarters of eukaryotic DNA are wrapped around protein cylinders forming so-called nucleosomes that block the access to the genetic information. Nucleosomes need therefore to be repositioned, either passively (by thermal fluctuations) or actively (by molecular motors). Here we introduce a theoretical model that allows us to study the interplay between a motor protein that moves along DNA (e.g., an RNA polymerase) and a nucleosome that it encounters on its way. We aim at describing the displacement mechanisms of the nucleosome and the motor protein on a microscopic level to understand better the intricate interplay between the active step of the motor and the nucleosome-repositioning step. Different motor types (Brownian ratchet versus power-stroke mechanism) that perform very similarly under a constant load are shown to have very different nucleosome repositioning capacities.

Varying rate of RNA chain elongation during rrn transcription in Escherichia coli

The value of the rRNA chain elongation rate in bacteria is an important physiological parameter, as it affects not only the rRNA promoter activity but also the free-RNA polymerase concentration and thereby the transcription of all genes. On average, rRNA chains elongate at a rate of 80 to 90 nucleotides (nt) per s, and the transcription of an entire rrn operon takes about 60 s (at 37 degrees C). Here we have analyzed a reported distribution obtained from electron micrographs of RNA polymerase molecules along rrn operons in E. coli growing at 2.5 doublings per hour (S. Quan, N. Zhang, S. French, and C. L. Squires, J. Bacteriol. 187:1632-1638, 2005). The distribution exhibits two peaks of higher polymerase density centered within the 16S and 23S rRNA genes. An evaluation of this distribution indicates that RNA polymerase transcribes the 5' leader region at speeds up to or greater than 250 nt/s. Once past the leader, transcription slows down to about 65 nt/s within the 16S gene, speeds up in the spacer region between the 16S and 23S genes, slows again to about 65 nt/s in the 23S region, and finally speeds up to a rate greater than 400 nt/s near the end of the operon. We suggest that the slowing of transcript elongation in the 16S and 23S sections is the result of transcriptional pauses, possibly caused by temporary interactions of the RNA polymerase with secondary structures in the nascent rRNA.

The origin of short transcriptional pauses

RNA polymerases are protein molecular machines that transcribe genetic information from DNA into RNA. The elongation of the RNA molecule is frequently interrupted by pauses, the detailed nature of which remains controversial. Here we ask whether backtracking, the central mechanism behind long pauses, could also be responsible for short pauses normally attributed to the ubiquitous pause state. To this end, we model backtracking as a force-biased random walk, giving rise to a broad distribution of pause durations as observed in experiments. Importantly, we find that this single mechanism naturally generates two populations of pauses that are distinct both in duration and trajectory: long-time pauses with the expected behavior of diffusive backtracks, and a new class of short-time backtracks with characteristics similar to those of the ubiquitous pause. These characteristics include an apparent force insensitivity and immobility of the polymerase. Based on these results and a quantitative comparison to published pause trajectories measured with optical tweezers, we suggest that a significant fraction of short pauses are simply due to backtracking.

Mechanism of sequence-specific pausing of bacterial RNA polymerase

Sequence-specific pausing of multisubunit RNA polymerases (RNAPs) represents a rate-limiting step during transcription elongation. Pausing occurs on average every 100 bases of DNA. Several models have been proposed to explain pausing, including backtracking of the ternary elongation complex, delay of translocation of the enzyme along DNA, or a conformational change in the active site preventing formation of the phosphodiester bond. Here, we performed biochemical characterization of previously-reported pauses of Escherichia coli RNAP and found that they are not associated with backtracking or a translocation delay. Instead, the paused complex contains the 3' end of the transcript in the active center and is capable of binding the next cognate NTP. However, bond formation occurs much slower in the paused complex compared with its fully-active counterpart. The pausing is dramatically decreased by a substitution of the base encoding the next incoming NTP and the base encoding the 3' end of the nascent RNA, suggesting that (mis)-alignment of the 3' end of the RNA and the incoming NTP in the active site is crucial for pausing. These pause sites are conserved between E. coli and Thermus thermophilus RNAPs, but are not recognized by Saccharomyces cerevisiae RNAP II, indicating that prokaryotic RNAPs might be more sensitive to the changes in the alignment of the nascent transcript and the substrate NTP in the active site.

Molecular dynamics studies of the energetics of translocation in model T7 RNA polymerase elongation complexes

Translocation in the single subunit T7 RNA polymerase elongation complex was studied by molecular dynamics simulations using the posttranslocated crystal structure with the fingers domain open, an intermediate stable in the absence of pyrophosphate, magnesium ions, and nucleotide substrate. Unconstrained and umbrella sampling simulations were performed to examine the energetics of translocations. The extent of translocation was quantified using reaction coordinates representing the average and individual displacements of the RNA-DNA hybrid base pairs with respect to a reference structure. In addition, an unconstrained simulation was also performed for the product complex with the fingers domain closed, but with the pyrophosphate and magnesium removed, in order to examine the local stability of the pretranslocated closed state after the pyrophosphate release. The average spatial movement of the entire hybrid was found to be energetically costly in the post- to pretranslocated direction in the open state, while the pretranslocated state was stable in the closed complex, supporting the notion that the conformational state dictates the global stability of translocation states. However, spatial fluctuations of the RNA 3'-end in the open conformation were extensive, with the typical range reaching 3-4 A. Our results suggest that thermal fluctuations play more important roles in the translocation of individual nucleotides than in the movement of large sections of nucleotide strands: RNA 3'-end can move into and out of the active site within a single conformational state, while a global movement of the hybrid may be thermodynamically unfavorable without the conformational change.

Interacting RNA polymerase motors on a DNA track: effects of traffic congestion and intrinsic noise on RNA synthesis

RNA polymerase (RNAP) is an enzyme that synthesizes a messenger RNA (mRNA) strand which is complementary to a single-stranded DNA template. From the perspective of physicists, an RNAP is a molecular motor that utilizes chemical energy input to move along the track formed by DNA. In many circumstances, which are described in this paper, a large number of RNAPs move simultaneously along the same track; we refer to such collective movements of the RNAPs as RNAP traffic. Here we develop a theoretical model for RNAP traffic by incorporating the steric interactions between RNAPs as well as the mechanochemical cycle of individual RNAPs during the elongation of the mRNA. By a combination of analytical and numerical techniques, we calculate the rates of mRNA synthesis and the average density profile of the RNAPs on the DNA track. We also introduce, and compute, two different measures of fluctuations in the synthesis of RNA. Analyzing these fluctuations, we show how the level of intrinsic noise in mRNA synthesis depends on the concentrations of the RNAPs as well as on those of some of the reactants and the products of the enzymatic reactions catalyzed by RNAP. We suggest appropriate experimental systems and techniques for testing our theoretical predictions.

Single-molecule studies of RNA polymerase: motoring along

Single-molecule techniques have advanced our understanding of transcription by RNA polymerase (RNAP). A new arsenal of approaches, including single-molecule fluorescence, atomic-force microscopy, magnetic tweezers, and optical traps (OTs) have been employed to probe the many facets of the transcription cycle. These approaches supply fresh insights into the means by which RNAP identifies a promoter, initiates transcription, translocates and pauses along the DNA template, proofreads errors, and ultimately terminates transcription. Results from single-molecule experiments complement the knowledge gained from biochemical and genetic assays by facilitating the observation of states that are otherwise obscured by ensemble averaging, such as those resulting from heterogeneity in molecular structure, elongation rate, or pause propensity. Most studies to date have been performed with bacterial RNAP, but work is also being carried out with eukaryotic polymerase (Pol II) and single-subunit polymerases from bacteriophages. We discuss recent progress achieved by single-molecule studies, highlighting some of the unresolved questions and ongoing debates.

Fluctuations, pauses, and backtracking in DNA transcription

Transcription is a vital stage in the process of gene expression and a major contributor to fluctuations in gene expression levels for which it is typically modeled as a single-step process with Poisson statistics. However, recent single molecule experiments raise questions about the validity of such a simple single-step picture. We present a molecular multistep model of transcription elongation that demonstrates that transcription times are in general non-Poisson-distributed. In particular, we model transcriptional pauses due to backtracking of the RNA polymerase as a first passage process. By including such pauses, we obtain a broad, heavy-tailed distribution of transcription elongation times, which can be significantly longer than would be otherwise. When transcriptional pauses result in long transcription times, we demonstrate that this naturally leads to bursts of mRNA production and non-Poisson statistics of mRNA levels. These results suggest that transcriptional pauses may be a significant contributor to the variability in transcription rates with direct implications for noise in cellular processes as well as variability between cells.

A Central Role of the RNA Polymerase Trigger Loop in Active-Site Rearrangement during Transcriptional Pausing

Transcriptional pausing by RNA polymerase is an underlying event in the regulation of transcript elongation, yet the physical changes in the transcribing complex that create the initially paused conformation remain poorly understood. We report that this nonbacktracked elemental pause results from an active-site rearrangement whose signature includes a trigger-loop conformation positioned near the RNA 3' nucleotide and a conformation of betaDloopII that allows fraying of the RNA 3' nucleotide away from the DNA template. During nucleotide addition, trigger-loop movements or folding appears to assist NTP-stimulated translocation and to be crucial for catalysis. At a pause, the trigger loop directly contributes to the paused conformation, apparently by restriction of its movement or folding, whereas a previously postulated unfolding of the bridge helix does not. This trigger-loop-centric model can explain many properties of transcriptional pausing.

Mechanochemical kinetics of transcription elongation

The mechanochemical kinetics of transcription elongation was examined with a combination of theoretical and experimental approaches. The predictive power of a sequence-dependent thermal ratchet model for transcription elongation was tested by establishing model parameters based solely on measurements under chemical perturbations and then directly predicting responses under mechanical perturbations without additional model parameters. Agreement between predicted and measured force-velocity curves provides strong support for a simple mechanochemical coupling mechanism.

Traffic of interacting ribosomes: effects of single-machine mechanochemistry on protein synthesis

Many ribosomes simultaneously move on the same messenger RNA (mRNA), each synthesizing separately a copy of the same protein. In contrast to the earlier models, here we develop a "unified" theoretical model that not only incorporates the mutual exclusions of the interacting ribosomes, but also describes explicitly the mechanochemistry of each of these macromolecular machines during protein synthesis. Using analytical and numerical techniques of nonequilibrium statistical mechanics, we analyze the rates of protein synthesis and the spatiotemporal organization of the ribosomes in this model. We also predict how these properties would change with the changes in the rates of the various chemomechanical processes in each ribosome. Finally, we illustrate the power of this model by making experimentally testable predictions on the rates of protein synthesis and the density profiles of the ribosomes on some mRNAs in E-coli.

Sequence-resolved detection of pausing by single RNA polymerase molecules

Transcriptional pausing by RNA polymerase (RNAP) plays an important role in the regulation of gene expression. Defined, sequence-specific pause sites have been identified biochemically. Single-molecule studies have also shown that bacterial RNAP pauses frequently during transcriptional elongation, but the relationship of these "ubiquitous" pauses to the underlying DNA sequence has been uncertain. We employed an ultrastable optical-trapping assay to follow the motion of individual molecules of RNAP transcribing templates engineered with repeated sequences carrying imbedded, sequence-specific pause sites of known regulatory function. Both the known and ubiquitous pauses appeared at reproducible locations, identified with base-pair accuracy. Ubiquitous pauses were associated with DNA sequences that show similarities to regulatory pause sequences. Data obtained for the lifetimes and efficiencies of pauses support a model where the transition to pausing branches off of the normal elongation pathway and is mediated by a common elemental state, which corresponds to the ubiquitous pause.