RNA polymerase II flexibility during translocation from normal mode analysis

T7 RNA polymerase translocation is facilitated by a helix opening on the fingers domain that may also prevent backtracking

Here, we studied the complete process of a viral T7 RNA polymerase (RNAP) translocation on DNA during transcription elongation by implementing extensive all-atom molecular dynamics (MD) simulations to construct a Markov state model (MSM). Our studies show that translocation proceeds in a Brownian motion, and the RNAP thermally transits among multiple metastable states. We observed non-synchronized backbone movements of the nucleic acid (NA) chains with the RNA translocation accomplished first, while the template DNA lagged. Notably, both the O-helix and Y-helix on the fingers domain play key roles in facilitating NA translocation through the helix opening. The helix opening allows a key residue Tyr639 to become inserted into the active site, which pushes the RNA–DNA hybrid forward. Another key residue, Phe644, coordinates the downstream template DNA motions by stacking and un-stacking with a transition nucleotide (TN) and its adjacent nucleotide. Moreover, the O-helix opening at pre-translocation (pre-trans) likely resists backtracking. To test this hypothesis, we computationally designed mutants of T7 RNAP by replacing the amino acids on the O-helix with counterpart residues from a mitochondrial RNAP that is capable of backtracking. The current experimental results support the hypothesis.

Flexibility-rigidity index for protein-nucleic acid flexibility and fluctuation analysis

Protein-nucleic acid complexes are important for many cellular processes including the most essential functions such as transcription and translation. For many protein-nucleic acid complexes, flexibility of both macromolecules has been shown to be critical for specificity and/or function. The flexibility-rigidity index (FRI) has been proposed as an accurate and efficient approach for protein flexibility analysis. In this article, we introduce FRI for the flexibility analysis of protein-nucleic acid complexes. We demonstrate that a multiscale strategy, which incorporates multiple kernels to capture various length scales in biomolecular collective motions, is able to significantly improve the state of art in the flexibility analysis of protein-nucleic acid complexes. We take the advantage of the high accuracy and O(N) computational complexity of our multiscale FRI method to investigate the flexibility of ribosomal subunits, which are difficult to analyze by alternative approaches. An anisotropic FRI approach, which involves localized Hessian matrices, is utilized to study the translocation dynamics in an RNA polymerase.

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.

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.

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.

Energetic and structural details of the trigger-loop closing transition in RNA polymerase II

An evolutionarily conserved element in RNA polymerase II, the trigger loop (TL), has been suggested to play an important role in the elongation rate, fidelity of selection of the matched nucleoside triphosphate (NTP), catalysis of transcription elongation, and translocation in both eukaryotes and prokaryotes. In response to NTP binding, the TL undergoes large conformational changes to switch between distinct open and closed states to tighten the active site and avail catalysis. A computational strategy for characterizing the conformational transition pathway is presented to bridge the open and closed states of the TL. Information from a large number of independent all-atom molecular dynamics trajectories from Hamiltonian replica exchange and targeted molecular dynamics simulations is gathered together to assemble a connectivity map of the conformational transition. The results show that with a cognate NTP, TL closing should be a spontaneous process. One major intermediate state is identified along the conformational transition pathway, and the key structural features are characterized. The complete pathway from the open TL to the closed TL provides a clear picture of the TL closing.

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.

Efficient prediction of protein conformational pathways based on the hybrid elastic network model

Various computational models have gained immense attention by analyzing the dynamic characteristics of proteins. Several models have achieved recognition by fulfilling either theoretical or experimental predictions. Nonetheless, each method possesses limitations, mostly in computational outlay and physical reality. These limitations remind us that a new model or paradigm should advance theoretical principles to elucidate more precisely the biological functions of a protein and should increase computational efficiency. With these critical caveats, we have developed a new computational tool that satisfies both physical reality and computational efficiency. In the proposed hybrid elastic network model (HENM), a protein structure is represented as a mixture of rigid clusters and point masses that are connected with linear springs. Harmonic analyses based on the HENM have been performed to generate normal modes and conformational pathways. The results of the hybrid normal mode analyses give new physical insight to the 70S ribosome. The feasibility of the conformational pathways of hybrid elastic network interpolation (HENI) was quantitatively evaluated by comparing three different overlap values proposed in this paper. A remarkable observation is that the obtained mode shapes and conformational pathways are consistent with each other. Our timing results show that HENM has some advantage in computational efficiency over a coarse-grained model, especially for large proteins, even though it takes longer to construct the HENM. Consequently, the proposed HENM will be one of the best alternatives to the conventional coarse-grained ENMs and all-atom based methods (such as molecular dynamics) without loss of physical reality.

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.

The RNA polymerase bridge helix YFI motif in catalysis, fidelity and translocation

The bridge α-helix in the β' subunit of RNA polymerase (RNAP) borders the active site and may have roles in catalysis and translocation. In Escherichia coli RNAP, a bulky hydrophobic segment near the N-terminal end of the bridge helix is identified (β' 772-YFI-774; the YFI motif). YFI is located at a distance from the active center and adjacent to a glycine hinge (β' 778-GARKG-782) involved in dynamic bending of the bridge helix. Remarkably, amino acid substitutions in YFI significantly alter intrinsic termination, pausing, fidelity and translocation of RNAP. F773V RNAP largely ignores the λ tR2 terminator at 200μM NTPs and is strongly reduced in λ tR2 recognition at 1μM NTPs. F773V alters RNAP pausing and backtracking and favors misincorporation. By contrast, the adjacent Y772A substitution increases fidelity and exhibits other transcriptional defects generally opposite to those of F773V. All atom molecular dynamics simulation revealed two separate functional connections emanating from YFI explaining the distinct effects of substitutions: Y772 communicates with the active site through the link domain in the β subunit, whereas F773 communicates through the fork domain in the β subunit. I774 interacts with the F-loop, which also contacts the glycine hinge of the bridge helix. These results identified negative and positive circuits coupled at YFI and employed for regulation of catalysis, elongation, termination and translocation.

Basic mechanisms of RNA polymerase II activity and alteration of gene expression in Saccharomyces cerevisiae

Transcription by RNA polymerase II (Pol II), and all RNA polymerases for that matter, may be understood as comprising two cycles. The first cycle relates to the basic mechanism of the transcription process wherein Pol II must select the appropriate nucleoside triphosphate (NTP) substrate complementary to the DNA template, catalyze phosphodiester bond formation, and translocate to the next position on the DNA template. Performing this cycle in an iterative fashion allows the synthesis of RNA chains that can be over one million nucleotides in length in some larger eukaryotes. Overlaid upon this enzymatic cycle, transcription may be divided into another cycle of three phases: initiation, elongation, and termination. Each of these phases has a large number of associated transcription factors that function to promote or regulate the gene expression process. Complicating matters, each phase of the latter transcription cycle are coincident with cotranscriptional RNA processing events. Additionally, transcription takes place within a highly dynamic and regulated chromatin environment. This chromatin environment is radically impacted by active transcription and associated chromatin modifications and remodeling, while also functioning as a major platform for Pol II regulation. This review will focus on our basic knowledge of the Pol II transcription mechanism, and how altered Pol II activity impacts gene expression in vivo in the model eukaryote Saccharomyces cerevisiae. This article is part of a Special Issue entitled: RNA Polymerase II Transcript Elongation.

Molecular dynamics and mutational analysis of the catalytic and translocation cycle of RNA polymerase

BACKGROUND

During elongation, multi-subunit RNA polymerases (RNAPs) cycle between phosphodiester bond formation and nucleic acid translocation. In the conformation associated with catalysis, the mobile "trigger loop" of the catalytic subunit closes on the nucleoside triphosphate (NTP) substrate. Closing of the trigger loop is expected to exclude water from the active site, and dehydration may contribute to catalysis and fidelity. In the absence of a NTP substrate in the active site, the trigger loop opens, which may enable translocation. Another notable structural element of the RNAP catalytic center is the "bridge helix" that separates the active site from downstream DNA. The bridge helix may participate in translocation by bending against the RNA/DNA hybrid to induce RNAP forward movement and to vacate the active site for the next NTP loading. The transition between catalytic and translocation conformations of RNAP is not evident from static crystallographic snapshots in which macromolecular motions may be restrained by crystal packing.

RESULTS

All atom molecular dynamics simulations of Thermus thermophilus (Tt) RNAP reveal flexible hinges, located within the two helices at the base of the trigger loop, and two glycine hinges clustered near the N-terminal end of the bridge helix. As simulation progresses, these hinges adopt distinct conformations in the closed and open trigger loop structures. A number of residues (described as "switch" residues) trade atomic contacts (ion pairs or hydrogen bonds) in response to changes in hinge orientation. In vivo phenotypes and in vitro activities rendered by mutations in the hinge and switch residues in Saccharomyces cerevisiae (Sc) RNAP II support the importance of conformational changes predicted from simulations in catalysis and translocation. During simulation, the elongation complex with an open trigger loop spontaneously translocates forward relative to the elongation complex with a closed trigger loop.

CONCLUSIONS

Switching between catalytic and translocating RNAP forms involves closing and opening of the trigger loop and long-range conformational changes in the atomic contacts of amino acid side chains, some located at a considerable distance from the trigger loop and active site. Trigger loop closing appears to support chemistry and the fidelity of RNA synthesis. Trigger loop opening and limited bridge helix bending appears to promote forward nucleic acid translocation.

Active site opening and closure control translocation of multisubunit RNA polymerase

Multisubunit RNA polymerase (RNAP) is the central information-processing enzyme in all cellular life forms, yet its mechanism of translocation along the DNA molecule remains conjectural. Here, we report direct monitoring of bacterial RNAP translocation following the addition of a single nucleotide. Time-resolved measurements demonstrated that translocation is delayed relative to nucleotide incorporation and occurs shortly after or concurrently with pyrophosphate release. An investigation of translocation equilibrium suggested that the strength of interactions between RNA 3′ nucleotide and nucleophilic and substrate sites determines the translocation state of transcription elongation complexes, whereas active site opening and closure modulate the affinity of the substrate site, thereby favoring the post- and pre-translocated states, respectively. The RNAP translocation mechanism is exploited by the antibiotic tagetitoxin, which mimics pyrophosphate and induces backward translocation by closing the active site.

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.

iMod: multipurpose normal mode analysis in internal coordinates

MOTIVATION

Dynamic simulations of systems with biologically relevant sizes and time scales are critical for understanding macromolecular functioning. Coarse-grained representations combined with normal mode analysis (NMA) have been established as an alternative to atomistic simulations. The versatility and efficiency of current approaches normally based on Cartesian coordinates can be greatly enhanced with internal coordinates (IC).

RESULTS

Here, we present a new versatile tool chest to explore conformational flexibility of both protein and nucleic acid structures using NMA in IC. Consideration of dihedral angles as variables reduces the computational cost and non-physical distortions of classical Cartesian NMA methods. Our proposed framework operates at different coarse-grained levels and offers an efficient framework to conduct NMA-based conformational studies, including standard vibrational analysis, Monte-Carlo simulations or pathway exploration. Examples of these approaches are shown to demonstrate its applicability, robustness and efficiency.

CONTACT

pablo@chaconlab.org

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

RNA polymerase II with open and closed trigger loops: active site dynamics and nucleic acid translocation

RNA polymerase II is the central eukaryotic enzyme in transcription from DNA to RNA. The dynamics of RNA polymerase II is described from molecular-dynamics simulations started from two crystal structures with open and closed trigger loop (TL) forms. Dynamic transitions between neutral and forward translocated states were observed, especially for the downstream DNA duplex. Dynamic rearrangements were also seen in the active site environment, including conformations in which the active site nucleotide assumed a possibly precatalytic conformation in close proximity to the terminal 3'-hydroxyl of the nascent RNA. Because nucleic acid translocation was observed primarily in the simulations with an open TL structure, whereas close approach of the active site nucleotide to the terminal RNA ribose predominantly occurred in the closed TL structure, a modified Brownian ratchet mechanism is proposed whereby thermally driven translocation is only possible with an open TL, and fidelity control and catalysis require TL closing.

The nucleotide addition cycle of RNA polymerase is controlled by two molecular hinges in the Bridge Helix domain

Background

Cellular RNA polymerases (RNAPs) are complex molecular machines that combine catalysis with concerted conformational changes in the active center. Previous work showed that kinking of a hinge region near the C-terminus of the Bridge Helix (BH-H_C) plays a critical role in controlling the catalytic rate.

Results

Here, new evidence for the existence of an additional hinge region in the amino-terminal portion of the Bridge Helix domain (BH-H_N) is presented. The nanomechanical properties of BH-H_N emerge as a direct consequence of the highly conserved primary amino acid sequence. Mutations that are predicted to influence its flexibility cause corresponding changes in the rate of the nucleotide addition cycle (NAC). BH-H_N displays functional properties that are distinct from BH-H_C, suggesting that conformational changes in the Bridge Helix control the NAC via two independent mechanisms.

Conclusions

The properties of two distinct molecular hinges in the Bridge Helix of RNAP determine the functional contribution of this domain to key stages of the NAC by coordinating conformational changes in surrounding domains.

Conformational coupling, bridge helix dynamics and active site dehydration in catalysis by RNA polymerase

Molecular dynamics simulation of Thermus thermophilus (Tt) RNA polymerase (RNAP) in a catalytic conformation demonstrates that the active site dNMP-NTP base pair must be substantially dehydrated to support full active site closing and optimum conditions for phosphodiester bond synthesis. In silico mutant beta R428A RNAP, which was designed based on substitutions at the homologous position (Rpb2 R512) of Saccharomyces cerevisiae (Sc) RNAP II, was used as a reference structure to compare to Tt RNAP in simulations. Long range conformational coupling linking a dynamic segment of the bridge alpha-helix, the extended fork loop, the active site, and the trigger loop-trigger helix is apparent and adversely affected in beta R428A RNAP. Furthermore, bridge helix bending is detected in the catalytic structure, indicating that bridge helix dynamics may regulate phosphodiester bond synthesis as well as translocation. An active site "latch" assembly that includes a key trigger helix residue Tt beta' H1242 and highly conserved active site residues beta E445 and R557 appears to help regulate active site hydration/dehydration. The potential relevance of these observations in understanding RNAP and DNAP induced fit and fidelity is discussed.