Diffusion of nucleoside triphosphates and role of the entry site to the RNA polymerase II active center

Binding Mechanism of the Active Form of Molnupiravir to RdRp of SARS-CoV-2 and Designing Potential Analogues: Insights from Molecular Dynamics Simulations

Molnupiravir, an FDA-approved nucleoside prodrug for treating COVID-19, converts into N4-hydroxycytidine triphosphate (NHC-TP), which integrates into SARS-CoV-2 RNA by its RNA-dependent RNA polymerase (RdRp) causing lethal mutations in viral proteins. Due to the risk of RdRp-mediated drug resistance and potential off-target effects on host polymerases (e.g., human polymerase II/HPolII), it is crucial to understand NHC-TP interactions at polymerase active sites for developing new, resistance-proof treatments. In this study, we used molecular dynamics (MD) simulations to probe key interactions between NHC-TP and SARS-CoV-2 RdRp and designed novel NHC-TP analogues with greater selectivity for SARS-CoV-2 RdRp over HPolII by a virtual screening workflow. We docked NHC-TP to a modified SARS-CoV-2 RdRp-Remdesivir triphosphate structure (PDB ID: 7BV2) and generated 71 NHC-TP analogues with bulky substituents to increase the interaction with RdRP and to reduce HPolII incorporation. MD simulations assessed the stability, binding affinity, and site interactions of these analogues. The top 7 candidates, with favorable ADMET properties, likely inhibit replication via potential dual mechanisms (the replicative stalling and the induction of lethal mutations) while maintaining selectivity for SARS-CoV-2 RdRp.

Transcription elongation mechanisms of RNA polymerases I, II, and III and their therapeutic implications

Transcription is a tightly regulated, complex, and essential cellular process in all living organisms. Transcription is comprised of three steps, transcription initiation, elongation, and termination. The distinct transcription initiation and termination mechanisms of eukaryotic RNA polymerases I, II, and III (Pols I, II, and III) have long been appreciated. Recent methodological advances have empowered high-resolution investigations of the Pols' transcription elongation mechanisms. Here, we review the kinetic similarities and differences in the individual steps of Pol I-, II-, and III-catalyzed transcription elongation, including NTP binding, bond formation, pyrophosphate release, and translocation. This review serves as an important summation of Saccharomyces cerevisiae (yeast) Pol I, II, and III kinetic investigations which reveal that transcription elongation by the Pols is governed by distinct mechanisms. Further, these studies illustrate how basic, biochemical investigations of the Pols can empower the development of chemotherapeutic compounds.

A unified Watson-Crick geometry drives transcription of six-letter expanded DNA alphabets by E. coli RNA polymerase

Artificially Expanded Genetic Information Systems (AEGIS) add independently replicable unnatural nucleotide pairs to the natural G:C and A:T/U pairs found in native DNA, joining the unnatural pairs through alternative modes of hydrogen bonding. Whether and how AEGIS pairs are recognized and processed by multi-subunit cellular RNA polymerases (RNAPs) remains unknown. Here, we show that E. coli RNAP selectively recognizes unnatural nucleobases in a six-letter expanded genetic system. High-resolution cryo-EM structures of three RNAP elongation complexes containing template-substrate UBPs reveal the shared principles behind the recognition of AEGIS and natural base pairs. In these structures, RNAPs are captured in an active state, poised to perform the chemistry step. At this point, the unnatural base pair adopts a Watson-Crick geometry, and the trigger loop is folded into an active conformation, indicating that the mechanistic principles underlying recognition and incorporation of natural base pairs also apply to AEGIS unnatural base pairs. These data validate the design philosophy of AEGIS unnatural basepairs. Further, we provide structural evidence supporting a long-standing hypothesis that pair mismatch during transcription occurs via tautomerization. Together, our work highlights the importance of Watson-Crick complementarity underlying the design principles of AEGIS base pair recognition.

The A12.2 Subunit Plays an Integral Role in Pyrophosphate Release of RNA Polymerase I

RNA polymerase I (Pol I) synthesizes ribosomal RNA (rRNA), which is the first and rate-limiting step in ribosome biosynthesis. A12.2 (A12) is a critical subunit of Pol I that is responsible for activating Pol I's exonuclease activity. We previously reported a kinetic mechanism for single-nucleotide incorporation catalyzed by Pol I lacking the A12 subunit (ΔA12 Pol I) purified from S. cerevisae and revealed that ΔA12 Pol I exhibited much slower incorporation compared to Pol I. However, it is unknown if A12 influences each nucleotide incorporation in the context of transcription elongation. Here, we show that A12 contributes to every repeating cycle of nucleotide addition and that deletion of A12 results in an entirely different kinetic mechanism compared to WT Pol I. We found that instead of one irreversible step between each nucleotide addition cycle, as reported for wild type (WT) Pol I, the ΔA12 variant requires one reversible step to describe each nucleotide addition. Reversibility fundamentally requires slow PPi release. Consistently, we show that Pol I is more pyrophosphate (PPi) concentration dependent than ΔA12 Pol I. This observation supports the model that PPi is retained in the active site of ΔA12 Pol I longer than WT Pol I. These results suggest that A12 promotes PPi release, revealing a larger role for the A12.2 subunit in the nucleotide addition cycle beyond merely activating exonuclease activity.

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.

Transcriptional processing of an unnatural base pair by eukaryotic RNA polymerase II

The development of unnatural base pairs (UBPs) has greatly increased the information storage capacity of DNA, allowing for transcription of unnatural RNA by the heterologously expressed T7 RNA polymerase (RNAP) in Escherichia coli. However, little is known about how UBPs are transcribed by cellular RNA polymerases. Here, we investigated how synthetic unnatural nucleotides, NaM and TPT3, are recognized by eukaryotic RNA polymerase II (Pol II) and found that Pol II is able to selectively recognize UBPs with high fidelity when dTPT3 is in the template strand and rNaMTP acts as the nucleotide substrate. Our structural analysis and molecular dynamics simulation provide structural insights into transcriptional processing of UBPs in a stepwise manner. Intriguingly, we identified a novel 3'-RNA binding site after rNaM addition, termed the swing state. These results may pave the way for future studies in the design of transcription and translation strategies in higher organisms with expanded genetic codes.

Uncovering the release mechanism of nucleotide import by HIV-1 capsid

Efficient nucleotide import is critical to fuel the reverse DNA synthesis that takes place within the HIV-1 capsid. However, the mechanism by which the HIV-1 capsid imports nucleotides is presently unclear. In this work, we carry out a series of Brownian diffusion simulations to elucidate the nucleotide import process through the hexamer pores of the HIV-1 capsid. Our simulations reveal a significant role of the electrostatic field in the import process and the mechanism by which deoxynucleoside triphosphates (dNTPs) diffuse through the arginine ring: specifically, how IP6s and ATPs, though competing with dNTPs for binding at the pore of the arginine ring, end up accelerating the dNTP import rate by thousands of folds so that it is sufficiently high to fuel the encapsidated DNA synthesis.

Nucleotide Loading Modes of Human RNA Polymerase II as Deciphered by Molecular Simulations

Mapping the route of nucleoside triphosphate (NTP) entry into the sequestered active site of RNA polymerase (RNAP) has major implications for elucidating the complete nucleotide addition cycle. Constituting a dichotomy that remains to be resolved, two alternatives, direct NTP delivery via the secondary channel (CH2) or selection to downstream sites in the main channel (CH1) prior to catalysis, have been proposed. In this study, accelerated molecular dynamics simulations of freely diffusing NTPs about RNAPII were applied to refine the CH2 model and uncover atomic details on the CH1 model that previously lacked a persuasive structural framework to illustrate its mechanism of action. Diffusion and binding of NTPs to downstream DNA, and the transfer of a preselected NTP to the active site, are simulated for the first time. All-atom simulations further support that CH1 loading is transcription factor IIF (TFIIF) dependent and impacts catalytic isomerization. Altogether, the alternative nucleotide loading systems may allow distinct transcriptional landscapes to be expressed.

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.

Structural mechanism of transcription inhibition by lasso peptides microcin J25 and capistruin

Many bacteria produce antimicrobial peptides for survival under stressful conditions. Some of these antimicrobial peptides are lasso peptides, which have a unique lasso-like topology and have generated great interest as a result of their stability in harsh conditions and amenability to functional engineering. In this study, we determined crystal structures of two lasso peptides, microcin J25 and capistruin, bound to their natural enzymatic target, the bacterial RNA polymerase (RNAP). The structures define peptide inhibitor–RNAP interactions that are important for inhibition and provide detailed insight into how the peptides inhibit RNAP function. This work provides a structural basis to guide the design of more potent lasso peptide antimicrobial approaches.

We report crystal structures of the antibacterial lasso peptides microcin J25 (MccJ25) and capistruin (Cap) bound to their natural enzymatic target, the bacterial RNA polymerase (RNAP). Both peptides bind within the RNAP secondary channel, through which NTP substrates enter the RNAP active site, and sterically block trigger-loop folding, which is essential for efficient catalysis by the RNAP. MccJ25 binds deep within the secondary channel in a manner expected to interfere with NTP substrate binding, explaining the partial competitive mechanism of inhibition with respect to NTPs found previously [Mukhopadhyay J, Sineva E, Knight J, Levy RM, Ebright RH (2004) Mol Cell 14:739–751]. The Cap binding determinant on RNAP overlaps, but is not identical to, that of MccJ25. Cap binds further from the RNAP active site and does not sterically interfere with NTP binding, and we show that Cap inhibition is partially noncompetitive with respect to NTPs. This work lays the groundwork for structure determination of other lasso peptides that target the bacterial RNAP and provides a structural foundation to guide lasso peptide antimicrobial engineering approaches.

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.

Allosteric Effector ppGpp Potentiates the Inhibition of Transcript Initiation by DksA

DksA and ppGpp are the central players in the stringent response and mediate a complete reprogramming of the transcriptome. A major component of the response is a reduction in ribosome synthesis, which is accomplished by the synergistic action of DksA and ppGpp bound to RNA polymerase (RNAP) inhibiting transcription of rRNAs. Here, we report the X-ray crystal structures of Escherichia coli RNAP in complex with DksA alone and with ppGpp. The structures show that DksA accesses the template strand at the active site and the downstream DNA binding site of RNAP simultaneously and reveal that binding of the allosteric effector ppGpp reshapes the RNAP-DksA complex. The structural data support a model for transcriptional inhibition in which ppGpp potentiates the destabilization of open complexes by DksA. This work establishes a structural basis for understanding the pleiotropic effects of DksA and ppGpp on transcriptional regulation in proteobacteria.

Mechanism of NTP Binding to the Active Site of T7 RNA Polymerase Revealed by Free-Energy Simulation

In genetic transcription, molecular dynamic details and energetics of NTP binding to the active site of RNA polymerase (RNAP) are poorly understood. In this article, we investigated the NTP binding process in T7 RNAP using all-atom MD simulation combined with the umbrella sampling technique. Based on our simulations, a two-step mechanism was proposed to explain NTP binding: first, substrate NTP in aqueous solution, which carries a magnesium ion, diffuses through a secondary channel of RNAP to attain a pore region, where it undergoes conformational changes to give a correct orientation; next, the NTP establishes initial basepairing contacts with the template nucleoside (TN). Our free-energy calculations suggest that both steps are spontaneous. This mechanism can easily explain the problem of NTP binding with different orientations. Moreover, it is found that the nascent NTP:TN basepair is fragile and easily broken by thermal disturbance. Therefore, we speculate that the fingers domain will be triggered to close, so as to create a steady environment for the next chemical step. The observations from the work provide valuable information for comprehensively understanding the mechanism of the basic step in genetic transcription.

Kinetics of nucleotide entry into RNA polymerase active site provides mechanism for efficiency and fidelity

During transcription, RNA polymerase II elongates RNA by adding nucleotide triphosphates (NTPs) complementary to a DNA template. Structural studies have suggested that NTPs enter and exit the active site via the narrow secondary pore but details have remained unclear. A kinetic model is presented that integrates molecular dynamics simulations with experimental data. Previous simulations of trigger loop dynamics and the dynamics of matched and mismatched NTPs in and near the active site were combined with new simulations describing NTP exit from the active site via the secondary pore. Markov state analysis was applied to identify major states and estimate kinetic rates for transitions between those states. The kinetic model predicts elongation and misincorporation rates in close agreement with experiment and provides mechanistic hypotheses for how NTP entry and exit via the secondary pore is feasible and a key feature for achieving high elongation and low misincorporation rates during RNA elongation.

Nucleotide Selectivity at a Preinsertion Checkpoint of T7 RNA Polymerase Transcription Elongation

Nucleotide selection is crucial for transcription fidelity control, in particular, for viral T7 RNA polymerase (RNAP) lack of proofreading activity. It has been recognized that multiple kinetic checkpoints exist prior to full nucleotide incorporation. In this work, we implemented intensive atomistic molecular dynamics (MD) simulations to quantify how strong the nucleotide selection is at the initial checkpoint of an elongation cycle of T7 RNAP. The incoming nucleotides bind into a preinsertion site where a critical tyrosine residue locates nearby to assist the nucleotide selection. We calculated the relative binding free energy between a noncognate nucleotide and a cognate one at a preinsertion configuration via alchemical simulations, showing that a small selection free energy or the binding free energy difference (∼3 k_BT) exists between the two nucleotides. Indeed, another preinsertion configuration favored by the noncognate nucleotides was identified, which appears to be off path for further nucleotide insertion and additionally assists the nucleotide selection. By chemical master equation (CME) approach, we show that the small selection free energy at the preinsertion site along with the off-path noncognate nucleotide filtering can help substantially to reduce the error rate and to maintain the elongation rate high in the T7 RNAP transcription.

Bridge helix bending promotes RNA polymerase II backtracking through a critical and conserved threonine residue

The dynamics of the RNA polymerase II (Pol II) backtracking process is poorly understood. We built a Markov State Model from extensive molecular dynamics simulations to identify metastable intermediate states and the dynamics of backtracking at atomistic detail. Our results reveal that Pol II backtracking occurs in a stepwise mode where two intermediate states are involved. We find that the continuous bending motion of the Bridge helix (BH) serves as a critical checkpoint, using the highly conserved BH residue T831 as a sensing probe for the 3′-terminal base paring of RNA:DNA hybrid. If the base pair is mismatched, BH bending can promote the RNA 3′-end nucleotide into a frayed state that further leads to the backtracked state. These computational observations are validated by site-directed mutagenesis and transcript cleavage assays, and provide insights into the key factors that regulate the preferences of the backward translocation.

RNA Polymerase II can detect and cleave mis-incorporated nucleotides by a proofreading mechanism that requires backtracking of the enzyme. Here the authors show that Pol II backtracking occurs in a stepwise mode that involves two intermediate states where the fraying of the terminal RNA nucleotide is a prerequisite.

In silico prediction of ebolavirus RNA polymerase inhibition by specific combinations of approved nucleotide analogues

BACKGROUND AND OBJECTIVE

The urgency of ebolavirus drug development is obvious in light of the current local epidemic in Western Africa with high morbidity and a risk of wider spread. We present an in silico study as a first step to identify inhibitors of ebolavirus polymerase activity based on approved antiviral nucleotide analogues.

STUDY DESIGN

Since a structure model of the ebolavirus polymerase is lacking, we performed combined homology and ab initio modeling and report a similarity to known polymerases of human enterovirus, bovine diarrhea virus and foot-and-mouth disease virus. This facilitated the localization of a nucleotide binding domain in the ebolavirus polymerase. We next performed molecular docking studies with nucleotides (ATP, CTP, GTP and UTP) and nucleotide analogues, including a variety of approved antiviral drugs.

RESULTS AND CONCLUSIONS

Specific combinations of nucleotide analogues significantly reduce the ligand-protein interaction energies of the ebolavirus polymerase for natural nucleotides. Any nucleotide analogue on its own did not reduce ligand-protein interaction energies. This prediction encourages specific drug testing efforts and guides future strategies to inhibit ebolavirus replication.

RNA polymerase II transcriptional fidelity control and its functional interplay with DNA modifications

Accurate genetic information transfer is essential for life. As a key enzyme involved in the first step of gene expression, RNA polymerase II (Pol II) must maintain high transcriptional fidelity while it reads along DNA template and synthesizes RNA transcript in a stepwise manner during transcription elongation. DNA lesions or modifications may lead to significant changes in transcriptional fidelity or transcription elongation dynamics. In this review, we will summarize recent progress toward understanding the molecular basis of RNA Pol II transcriptional fidelity control and impacts of DNA lesions and modifications on Pol II transcription elongation.

Concerted Protein and Nucleic Acid Conformational Changes Observed Prior to Nucleotide Incorporation in a Bacterial RNA Polymerase: Raman Crystallographic Evidence

Transcription elongation requires the continuous incorporation of ribonucleotide triphosphates into a growing transcript. RNA polymerases (RNAPs) are able to processively synthesize a growing RNA chain via translocation of the RNAP enzyme along its nucleic acid template strand after each nucleotide addition cycle. In this work, a time-resolved Raman spectroscopic analysis of nucleotide addition in single crystals of the Thermus thermophilus elongation complex (TthEC) is reported. When [(13)C,(15)N]GTP (*GTP) is soaked into crystals of the TthEC, large reversible changes in the Raman spectrum that are assigned to protein and nucleic acid conformational events during a single-nucleotide incorporation are observed. The *GTP population in the TthEC crystal reaches a stable population at 37 min, while substantial and reversible protein conformational changes (mainly ascribed to changes in α-helical Raman features) maximize at approximately 50 min. At the same time, changes in nucleic acid bases and phosphodiester backbone Raman marker bands occur. Catalysis begins at approximately 65-70 min, soon after the maximal protein and DNA changes, and is monitored via the decline in a triphosphate vibrational Raman mode from *GTP. The Raman data indicate that approximately 40% of the total triphosphate population, present as *GTP, reacts in the crystal. This may suggest that a second population of noncovalently bound *GTP resides in a site distinct from the catalytic site. The data reported here are an extension of our recent work on the elongation complex (EC) of a bacterial RNAP, Thermus thermophilus (Tth), where Raman spectroscopy and polyacrylamide gel electrophoresis were employed to monitor incorporation and misincorporation in single TthEC crystals [Antonopoulos, I. H., et al. (2015) Biochemistry 54, 652-665]. Therefore, the initial study establishes the groundwork for this study. In contrast to our previous study, in which incorporation takes place very rapidly inside the crystals, the data on this single crystal exhibit a slower time regime, which allows the dissection of the structural dynamics associated with GMP incorporation within the TthEC crystal.

Structural Model of RNA Polymerase II Elongation Complex with Complete Transcription Bubble Reveals NTP Entry Routes

The RNA polymerase II (Pol II) is a eukaryotic enzyme that catalyzes the synthesis of the messenger RNA using a DNA template. Despite numerous biochemical and biophysical studies, it remains elusive whether the “secondary channel” is the only route for NTP to reach the active site of the enzyme or if the “main channel” could be an alternative. On this regard, crystallographic structures of Pol II have been extremely useful to understand the structural basis of transcription, however, the conformation of the unpaired non-template DNA part of the full transcription bubble (TB) is still unknown. Since diffusion routes of the nucleoside triphosphate (NTP) substrate through the main channel might overlap with the TB region, gaining structural information of the full TB is critical for a complete understanding of Pol II transcription process. In this study, we have built a structural model of Pol II with a complete transcription bubble based on multiple sources of existing structural data and used Molecular Dynamics (MD) simulations together with structural analysis to shed light on NTP entry pathways. Interestingly, we found that although both channels have enough space to allow NTP loading, the percentage of MD conformations containing enough space for NTP loading through the secondary channel is twice higher than that of the main channel. Further energetic study based on MD simulations with NTP loaded in the channels has revealed that the diffusion of the NTP through the main channel is greatly disfavored by electrostatic repulsion between the NTP and the highly negatively charged backbones of nucleotides in the non-template DNA strand. Taken together, our results suggest that the secondary channel is the major route for NTP entry during Pol II transcription.

In eukaryotic cells, the RNA polymerase II (Pol II) is a central enzyme that reads the genetic information encoded in the DNA template to synthetize a messenger RNA. To perform its function, Pol II needs to have the substrate nucleoside triphosphate (NTP) diffuse into its deeply buried active site. Despite numerous efforts, the NTP entry routes remain elusive: NTP could diffuse only through the secondary channel, or also via the main channel. The structural information of the transcription bubble is essential to study this process, however, the unpaired non-template DNA of the transcription bubble is absent in the available X-ray crystal structures. In this regard, we have built a structural model of the Pol II elongation complex with reconstructed transcription bubble using existing experimental data. We then performed Molecular Dynamics (MD) simulations and applied structural analysis to study the routes of NTP diffusion. We found that sterically the probability of NTP loading through the secondary channel is more than twice that of the main channel. Further analysis of the non-bonded energetic contributions to NTP diffusion suggests that NTP diffusion through the main channel is greatly disfavored by the electrostatic repulsion between the substrate and negatively charged backbones of nucleotides in the non-template strand of the transcription bubble. Altogether, our findings suggest that the secondary channel is the more favorable NTP diffusion route for Pol II transcription elongation.

Pre-mRNA splicing is facilitated by an optimal RNA polymerase II elongation rate

Fong et al. examined cotranscriptional pre-mRNA splicing using RNA polymerase II mutants that change average elongation rates genome-wide. Slow and fast elongation affected constitutive and alternative splicing and often both increased or both decreased inclusion of a particular exon or retained intron. These results suggest that an optimal rate of transcriptional elongation is required for normal cotranscriptional pre-mRNA splicing.

Alternative splicing modulates expression of most human genes. The kinetic model of cotranscriptional splicing suggests that slow elongation expands and that fast elongation compresses the “window of opportunity” for recognition of upstream splice sites, thereby increasing or decreasing inclusion of alternative exons. We tested the model using RNA polymerase II mutants that change average elongation rates genome-wide. Slow and fast elongation affected constitutive and alternative splicing, frequently altering exon inclusion and intron retention in ways not predicted by the model. Cassette exons included by slow and excluded by fast elongation (type I) have weaker splice sites, shorter flanking introns, and distinct sequence motifs relative to “slow-excluded” and “fast-included” exons (type II). Many rate-sensitive exons are misspliced in tumors. Unexpectedly, slow and fast elongation often both increased or both decreased inclusion of a particular exon or retained intron. These results suggest that an optimal rate of transcriptional elongation is required for normal cotranscriptional pre-mRNA splicing.

Five checkpoints maintaining the fidelity of transcription by RNA polymerases in structural and energetic details

Transcriptional fidelity, which prevents the misincorporation of incorrect nucleoside monophosphates in RNA, is essential for life. Results from molecular dynamics (MD) simulations of eukaryotic RNA polymerase (RNAP) II and bacterial RNAP with experimental data suggest that fidelity may involve as many as five checkpoints. Using MD simulations, the effects of different active site NTPs in both open and closed trigger loop (TL) structures of RNAPs are compared. Unfavorable initial binding of mismatched substrates in the active site with an open TL is proposed to be the first fidelity checkpoint. The leaving of an incorrect substrate is much easier than a correct one energetically from the umbrella sampling simulations. Then, the closing motion of the TL, required for catalysis, is hindered by the presence of mismatched NTPs. Mismatched NTPs also lead to conformational changes in the active site, which perturb the coordination of magnesium ions and likely affect the ability to proceed with catalysis. This step appears to be the most important checkpoint for deoxy-NTP discrimination. Finally, structural perturbations in the template DNA and the nascent RNA in the presence of mismatches likely hinder nucleotide addition and provide the structural foundation for backtracking followed by removing erroneously incorporated nucleotides during proofreading.

RNA polymerase II conserved protein domains as platforms for protein-protein interactions

RNA polymerase II establishes many protein-protein interactions with transcriptional regulators to coordinate gene expression, but little is known about protein domains involved in the contact with them. We use a new approach to look for conserved regions of the RNA pol II of S. cerevisiae located at the surface of the structure of the complex, hypothesizing that they might be involved in the interaction with transcriptional regulators. We defined five different conserved domains and demonstrate that all of them make contact with transcriptional regulators.

Millisecond dynamics of RNA polymerase II translocation at atomic resolution

Transcription is a central step in gene expression, in which the DNA template is processively read by RNA polymerase II (Pol II), synthesizing a complementary messenger RNA transcript. At each cycle, Pol II moves exactly one register along the DNA, a process known as translocation. Although X-ray crystal structures have greatly enhanced our understanding of the transcription process, the underlying molecular mechanisms of translocation remain unclear. Here we use sophisticated simulation techniques to observe Pol II translocation on a millisecond timescale and at atomistic resolution. We observe multiple cycles of forward and backward translocation and identify two previously unidentified intermediate states. We show that the bridge helix (BH) plays a key role accelerating the translocation of both the RNA:DNA hybrid and transition nucleotide by directly interacting with them. The conserved BH residues, Thr831 and Tyr836, mediate these interactions. To date, this study delivers the most detailed picture of the mechanism of Pol II translocation at atomic level.

Molecular basis of transcriptional fidelity and DNA lesion-induced transcriptional mutagenesis

Maintaining high transcriptional fidelity is essential for life. Some DNA lesions lead to significant changes in transcriptional fidelity. In this review, we will summarize recent progress towards understanding the molecular basis of RNA polymerase II (Pol II) transcriptional fidelity and DNA lesion-induced transcriptional mutagenesis. In particular, we will focus on the three key checkpoint steps of controlling Pol II transcriptional fidelity: insertion (specific nucleotide selection and incorporation), extension (differentiation of RNA transcript extension of a matched over mismatched 3'-RNA terminus), and proofreading (preferential removal of misincorporated nucleotides from the 3'-RNA end). We will also discuss some novel insights into the molecular basis and chemical perspectives of controlling Pol II transcriptional fidelity through structural, computational, and chemical biology approaches.

Dissecting the chemical interactions and substrate structural signatures governing RNA polymerase II trigger loop closure by synthetic nucleic acid analogues

The trigger loop (TL) of RNA polymerase II (Pol II) is a conserved structural motif that is crucial for Pol II catalytic activity and transcriptional fidelity. The TL remains in an inactive open conformation when the mismatched substrate is bound. In contrast, TL switches from an inactive open state to a closed active state to facilitate nucleotide addition upon the binding of the cognate substrate to the Pol II active site. However, a comprehensive understanding of the specific chemical interactions and substrate structural signatures that are essential to this TL conformational change remains elusive. Here we employed synthetic nucleotide analogues as ‘chemical mutation’ tools coupling with α-amanitin transcription inhibition assay to systematically dissect the key chemical interactions and structural signatures governing the substrate-coupled TL closure in Saccharomyces cerevisiae Pol II. This study reveals novel insights into understanding the molecular basis of TL conformational transition upon substrate binding during Pol II transcription. This synthetic chemical biology approach may be extended to understand the mechanisms of other RNA polymerases as well as other nucleic acid enzymes in future studies.

Coupling mRNA processing with transcription in time and space

Maturation of mRNA precursors often occurs simultaneously with their synthesis by RNA polymerase II (Pol II). The co-transcriptional nature of mRNA processing has permitted the evolution of coupling mechanisms that coordinate transcription with mRNA capping, splicing, editing and 3' end formation. Recent experiments using sophisticated new methods for analysis of nascent RNA have provided important insights into the relative amount of co-transcriptional and post-transcriptional processing, the relationship between mRNA elongation and processing, and the role of the Pol II carboxy-terminal domain (CTD) in regulating these processes.

THEORETICAL INVESTIGATIONS ON ELUCIDATING FUNDAMENTAL MECHANISMS OF CATALYSIS AND DYNAMICS INVOLVED IN TRANSCRIPTION BY RNA POLYMERASE

RNA polymerase is the enzyme that synthesizes RNA during the transcription process. To understand its mechanism, structural studies have provided us pictures of the series of steps necessary to add a new nucleotide to the nascent RNA chain, the steps altogether known as the nucleotide addition cycle (NAC). However, these static snapshots do not provide dynamic information of these processes involved in NAC, such as the conformational changes of the protein and the atomistic details of the catalysis. Computational studies have made efforts to fill these knowledge gaps. In this review, we provide examples of different computational approaches that have improved our understanding of the transcription elongation process for RNA polymerase, such as normal mode analysis, molecular dynamic (MD) simulations, Markov state models (MSMs). We also point out some unsolved questions that could be addressed using computational tools in the future.

A two-state model for the dynamics of the pyrophosphate ion release in bacterial RNA polymerase

The dynamics of the PPi release during the transcription elongation of bacterial RNA polymerase and its effects on the Trigger Loop (TL) opening motion are still elusive. Here, we built a Markov State Model (MSM) from extensive all-atom molecular dynamics (MD) simulations to investigate the mechanism of the PPi release. Our MSM has identified a simple two-state mechanism for the PPi release instead of a more complex four-state mechanism observed in RNA polymerase II (Pol II). We observed that the PPi release in bacterial RNA polymerase occurs at sub-microsecond timescale, which is ∼3-fold faster than that in Pol II. After escaping from the active site, the (Mg-PPi)(2-) group passes through a single elongated metastable region where several positively charged residues on the secondary channel provide favorable interactions. Surprisingly, we found that the PPi release is not coupled with the TL unfolding but correlates tightly with the side-chain rotation of the TL residue R1239. Our work sheds light on the dynamics underlying the transcription elongation of the bacterial RNA polymerase.

High-throughput purification of affinity-tagged recombinant proteins

X-ray crystallography is the method of choice for obtaining a detailed view of the structure of proteins. Such studies need to be complemented by further biochemical analyses to obtain detailed insights into structure/function relationships. Advances in oligonucleotide- and gene synthesis technology make large-scale mutagenesis strategies increasingly feasible, including the substitution of target residues by all 19 other amino acids. Gain- or loss-of-function phenotypes then allow systematic conclusions to be drawn, such as the contribution of particular residues to catalytic activity, protein stability and/or protein-protein interaction specificity. In order to attribute the different phenotypes to the nature of the mutation--rather than to fluctuating experimental conditions--it is vital to purify and analyse the proteins in a controlled and reproducible manner. High-throughput strategies and the automation of manual protocols on robotic liquid-handling platforms have created opportunities to perform such complex molecular biological procedures with little human intervention and minimal error rates. Here, we present a general method for the purification of His-tagged recombinant proteins in a high-throughput manner. In a recent study, we applied this method to a detailed structure-function investigation of TFIIB, a component of the basal transcription machinery. TFIIB is indispensable for promoter-directed transcription in vitro and is essential for the recruitment of RNA polymerase into a preinitiation complex. TFIIB contains a flexible linker domain that penetrates the active site cleft of RNA polymerase. This linker domain confers two biochemically quantifiable activities on TFIIB, namely (i) the stimulation of the catalytic activity during the 'abortive' stage of transcript initiation, and (ii) an additional contribution to the specific recruitment of RNA polymerase into the preinitiation complex. We exploited the high-throughput purification method to generate single, double and triple substitution and deletions mutations within the TFIIB linker and to subsequently analyse them in functional assays for their stimulation effect on the catalytic activity of RNA polymerase. Altogether, we generated, purified and analysed 381 mutants--a task which would have been time-consuming and laborious to perform manually. We produced and assayed the proteins in multiplicates which allowed us to appreciate any experimental variations and gave us a clear idea of the reproducibility of our results. This method serves as a generic protocol for the purification of His-tagged proteins and has been successfully used to purify other recombinant proteins. It is currently optimised for the purification of 24 proteins but can be adapted to purify up to 96 proteins.

Dynamics of pyrophosphate ion release and its coupled trigger loop motion from closed to open state in RNA polymerase II

Pyrophosphate ion (PP(i)) release after nucleotide incorporation is a necessary step for RNA polymerase II (pol II) to enter the next nucleotide addition cycle during transcription elongation. However, the role of pol II residues in PP(i) release and the mechanistic relationship between PP(i) release and the conformational change of the trigger loop remain unclear. In this study, we constructed a Markov state model (MSM) from extensive all-atom molecular dynamics (MD) simulations in the explicit solvent to simulate the PP(i) release process along the pol II secondary channel. Our results show that the trigger loop has significantly larger intrinsic motion after catalysis and formation of PP(i), which in turn aids PP(i) release mainly through the hydrogen bonding between the trigger loop residue H1085 and the (Mg-PP(i))(2-) group. Once PP(i) leaves the active site, it adopts a hopping model through several highly conserved positively charged residues such as K752 and K619 to release from the pol II pore region of the secondary channel. These positive hopping sites form favorable interactions with PP(i) and generate four kinetically metastable states as identified by our MSM. Furthermore, our single-mutant simulations suggest that H1085 and K752 aid PP(i) exit from the active site after catalysis, whereas K619 facilitates its passage through the secondary channel. Finally, we suggest that PP(i) release could help the opening motion of the trigger loop, even though PP(i) release precedes full opening of the trigger loop due to faster PP(i) dynamics. Our simulations provide predictions to guide future experimental tests.

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.

Quantifying transcription of clinically relevant immobilized DNA within a continuous flow microfluidic reactor

Flow-through reactors are commonly used to control and optimize reagent delivery and product removal. Although recent research suggests that transcription reactions using picogram quantities of cDNA produce RNA efficiently in a flow-through microreactor, there has not been a detailed study on the mass transport and reagent dependence of microfluidic transcription reactions. We present a novel microreactor that contains H5 influenza cDNA immobilized directly onto the reactor walls to study the kinetics and reagent dependence of in vitro transcription reactions on a microfluidic platform. Enzyme and the rNTP substrate continuously flow over the cDNA and create RNA, which flows to a downstream collection well. Using nanogram quantities of cDNA, we found that enzyme limiting conditions caused by the concentration of cDNA in a small-volume microreactor channel may be partially overcome as the enzyme binds and concentrates near the channel wall. Kinetics confirm this phenomenon and show that the timescale for enzyme binding can be approximated by t(f) = cDNA/Q[E]. Surprisingly, on-chip transcription reactions have a strong dependence on the rNTP concentration from 5 to 9 mM despite a low consumption rate of rNTP molecules that is largely independent of the flow rate. Faster flow rates decrease the time it takes to fill DNA promoter sites with enzyme while additionally refreshing rNTP and MgCl(2) to allow for a greater consumption of rNTP. These two effects cause reactions with higher concentrations of cDNA in the reactor channel to have a greater dependence on the flow rate. At high flow rates (>0.37 nL/s), the reaction rate begins to drop, likely because of the release and escape of enzyme molecules from the cDNA layer. This critical flow rate can be predicted by a new modified Peclet number, Pe(m) = L(c)V/D, where L(c) is the full length of the tightly packed cDNA molecules, V is the velocity at the DNA/fluid interface, and D is the diffusivity of the enzyme molecule. Together, these insights can inspire reactor designs for a variety of applications.

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.

Forks, pincers, and triggers: the tools for nucleotide incorporation and translocation in multi-subunit RNA polymerases

The central role of RNA polymerase (RNAP) is to catalyze the processive synthesis of a growing RNA transcript. Recent structural and biophysical data have led to a deeper understanding of the nucleotide addition cycle and insight into the structure-function relationships that govern transcription elongation. In this review, we discuss kinetic data on nucleotide incorporation in the context of crystal structures, which show RNAP in multiple conformations. We present a facilitated Brownian ratchet model of nucleotide incorporation, in which templated NTP binding to a non-catalytic site in the main channel promotes the conformational changes that lead to opening of the catalytic site and translocation.

"Cotranscriptionality": the transcription elongation complex as a nexus for nuclear transactions

Much of the complex process of RNP biogenesis takes place at the gene cotranscriptionally. The target for RNA binding and processing factors is, therefore, not a solitary RNA molecule but, rather, a transcription elongation complex (TEC) comprising the growing nascent RNA and RNA polymerase traversing a chromatin template with associated passenger proteins. RNA maturation factors are not the only nuclear machines whose work is organized cotranscriptionally around the TEC scaffold. Additionally, DNA repair, covalent chromatin modification, "gene gating" at the nuclear pore, Ig gene hypermutation, and sister chromosome cohesion have all been demonstrated or suggested to involve a cotranscriptional component. From this perspective, TECs can be viewed as potent "community organizers" within the nucleus.

RNA polymerase active center: the molecular engine of transcription

RNA polymerase (RNAP) is a complex molecular machine that governs gene expression and its regulation in all cellular organisms. To accomplish its function of accurately producing a full-length RNA copy of a gene, RNAP performs a plethora of chemical reactions and undergoes multiple conformational changes in response to cellular conditions. At the heart of this machine is the active center, the engine, which is composed of distinct fixed and moving parts that serve as the ultimate acceptor of regulatory signals and as the target of inhibitory drugs. Recent advances in the structural and biochemical characterization of RNAP explain the active center at the atomic level and enable new approaches to understanding the entire transcription mechanism, its exceptional fidelity and control.

Steric effects and mass-transfer limitations surrounding amplification reactions on immobilized long and clinically relevant DNA templates

DNA and RNA are commonly captured on solid substrates during purification and isolation, where they can be transferred to downstream amplification and transcription reactions. When compared to the solution phase, however, immobilized DNA- and RNA-directed reactions are less efficient because of a variety of complex factors. Steric inhibition because of the bead surface and neighboring biological polymers, a change in solution chemistry because of the high local concentration of template molecules, and mass transfer to the bead surface could all affect the overall reaction kinetics. Furthermore, these effects may be particularly evident when working with long clinically relevant molecules, such as mRNA, viral RNA, and cDNA. In this paper, we focus on the in vitro transcription reaction (IVT) of both a long and short strand of H5 influenza A RNA (1777 and 465 nt) on both free and immobilized DNA templates to study these phenomena. We found that transcription was less efficient on immobilized beads than in solution, but that it can be dramatically increased with optimal solution chemistry. Using high ribonucleotide concentrations (>6 mM total rNTP), the RNA yield from long immobilized cDNA templates was boosted to 60% of solution control. Surprisingly, we found that steric effects because of surrounding immobilized molecules were only significant when the DNA molecules were short enough to achieve a high density (9x10(-4) microm2/molecule) on the silica substrate, such that the gap between molecules is on the order of the polymerase diameter. Eventually, these findings can be exploited in an automated microreactor, where isolation, purification, amplification, and detection of nucleic acids can be unified into one portable device.

Millisecond phase kinetic analysis of elongation catalyzed by human, yeast, and Escherichia coli RNA polymerase

Strategies for assembly and analysis of human, yeast, and bacterial RNA polymerase elongation complexes are described, and methods are shown for millisecond phase kinetic analyses of elongation using rapid chemical quench flow. Human, yeast, and bacterial RNA polymerases function very similarly in NTP-Mg2+ commitment and phosphodiester bond formation. A "running start, two-bond, double-quench" protocol is described and its advantages discussed. These studies provide information about stable NTP-Mg2+ loading, phosphodiester bond synthesis, the processive transition between bonds, and sequence-specific effects on transcription elongation dynamics.

Exploring the molecular origin of the high selectivity of multisubunit RNA polymerases by stochastic kinetic models

RNA polymerases are molecular machines of great fidelity, which can recognize matched NTPs from unmatched NTPs and 2'-dNTPs. We investigated by a stochastic simulation algorithm the whole nucleotide addition cycle based on an event-driven model. This model allows us to examine possible molecular origins of the high fidelity of RNA polymerases. For unmatched NTP selectivity, the conclusions drawn from simulated elongation rates corroborate those derived from structural analysis. The presence of two conformations (Esite and pre-insertion site) for the incoming nucleotide before the polymerization reaction is sufficient to allow selectivity. Concerning sugar selectivity, our results indicate that selectivity is only achievable if slow chemical reactions occur for 2'-dNTP. These results can be used to understand recent experimental observations.

Promoter Escape by Escherichia coli RNA Polymerase

Promoter escape is the process that an initiated RNA polymerase (RNAP) molecule undergoes to achieve the initiation-elongation transition. Having made this transition, an RNAP molecule would be relinquished from its promoter hold to perform productive (full-length) transcription. Prior to the transition, this process is accompanied by abortive RNA formation-the amount and pattern of which is controlled by the promoter sequence information. Qualitative and quantitative analysis of abortive/productive transcription from several Escherichia coli promoters and their sequence variants led to the understanding that a strong (RNAP-binding) promoter is more likely to be rate limited (during transcription initiation) at the escape step and produce abortive transcripts. Of the two subelements in a promoter, the PRR (the core Promoter Recognition Region) was found to set the initiation frequency and the rate-limiting step, while the ITS (the Initial Transcribed Sequence region) modulated the ratio of abortive versus productive transcription. The highly abortive behavior of E. coli RNAP could be ameliorated by the presence of Gre (transcript cleavage stimulatory) factor(s), linking the first step in abortive RNA formation by the initial transcribing complexes (ITC) to RNAP backtracking. The discovery that translocation during the initiation stage occurs via DNA scrunching provided the source of energy that converts each ITC into a highly unstable "stressed intermediate." Mapping all of the biochemical information onto an X-ray crystallographic structural model of an open complex gave rise to a plausible mechanism of transcription initiation. The chapter concludes with contemplations of the kinetics and thermodynamics of abortive initiation-promoter escape.

Transient reversal of RNA polymerase II active site closing controls fidelity of transcription elongation

To study fidelity of RNA polymerase II (Pol II), we analyzed properties of the 6-azauracil-sensitive and TFIIS-dependent E1103G mutant of rbp1 (rpo21), the gene encoding the catalytic subunit of Pol II in Saccharomyces cerevisiae. Using an in vivo retrotransposition-based transcription fidelity assay, we observed that rpb1-E1103G causes a 3-fold increase in transcription errors. This mutant showed a 10-fold decrease in fidelity of transcription elongation in vitro. The mutation does not appear to significantly affect translocation state equilibrium of Pol II in a stalled elongation complex. Primarily, it promotes NTP sequestration in the polymerase active center. Furthermore, pre-steady-state analyses revealed that the E1103G mutation shifted the equilibrium between the closed and the open active center conformations toward the closed form. Thus, open conformation of the active center emerges as an intermediate essential for preincorporation fidelity control. Similar mechanisms may control fidelity of DNA-dependent DNA polymerases and RNA-dependent RNA polymerases.