Millisecond dynamics of RNA polymerase II translocation at atomic resolution

T7 RNA Polymerase Discriminates Correct and Incorrect Nucleoside Triphosphates by Free Energy

RNA polymerase (RNAP) is the primary machine responsible for transcription. Its ability to distinguish between correct (cognate) and incorrect (noncognate) nucleoside triphosphates (NTPs) is important for fidelity control in transcription. In this work, we investigated the substrate selection mechanism of T7 RNAP from the perspective of energetics. The dissociation free energies were determined for matched and unmatched base pairs in the preinsertion complex using the umbrella sampling method. A clear hydrogen-bond-rupture peak is observed in the potential of mean force curve for a matched base pair, whereas no such peaks are present in the position of mean force profiles for unmatched ones. The free-energy barrier could prevent correct substrates from being separated from the active site. Therefore, when NTPs diffuse into the active site, correct ones will stay for chemistry once they establish effective base pairing contacts with the template nucleotide, whereas incorrect ones will be withdrawn from the active site and rejected back to solution. This result provides an important energy evidence for the substrate selection mechanism of RNAP. Then we elucidated energetics and molecular details for correct NTP binding to the active site of the insertion complex. Our observations reveal that strong interactions act on the triphosphate of NTP to constrain its movement, whereas relatively weak interactions serve to position the base in the correct conformation. Triple interactions, hydrophobic contacts from residues M635 and Y639, base stacking from the 3' RNA terminal nucleotide, and base pairing from the template nucleotide act together to position the NTP base in a catalytically competent conformation. At last, we observed that incorrect NTPs cannot be as well-stabilized as the correct one in the active site when they are misincorporated in the insertion site. It is expected that our work can be helpful for comprehensively understanding details of this basic step in genetic transcription.

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.

Protobiotic Systems Chemistry Analyzed by Molecular Dynamics

Systems chemistry has been a key component of origin of life research, invoking models of life's inception based on evolving molecular networks. One such model is the graded autocatalysis replication domain (GARD) formalism embodied in a lipid world scenario, which offers rigorous computer simulation based on defined chemical kinetics equations. GARD suggests that the first pre-RNA life-like entities could have been homeostatically-growing assemblies of amphiphiles, undergoing compositional replication and mutations, as well as rudimentary selection and evolution. Recent progress in molecular dynamics has provided an experimental tool to study complex biological phenomena such as protein folding, ligand-receptor interactions, and micellar formation, growth, and fission. The detailed molecular definition of GARD and its inter-molecular catalytic interactions make it highly compatible with molecular dynamics analyses. We present a roadmap for simulating GARD's kinetic and thermodynamic behavior using various molecular dynamics methodologies. We review different approaches for testing the validity of the GARD model by following micellar accretion and fission events and examining compositional changes over time. Near-future computational advances could provide empirical delineation for further system complexification, from simple compositional non-covalent assemblies towards more life-like protocellular entities with covalent chemistry that underlies metabolism and genetic encoding.

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

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

Functional assays for transcription mechanisms in high-throughput

Dramatic increases in the scale of programmed synthesis of nucleic acid libraries coupled with deep sequencing have powered advances in understanding nucleic acid and protein biology. Biological systems centering on nucleic acids or encoded proteins greatly benefit from such high-throughput studies, given that large DNA variant pools can be synthesized and DNA, or RNA products of transcription, can be easily analyzed by deep sequencing. Here we review the scope of various high-throughput functional assays for studies of nucleic acids and proteins in general, followed by discussion of how these types of study have yielded insights into the RNA Polymerase II (Pol II) active site as an example. We discuss methodological considerations in the design and execution of these experiments that should be valuable to studies in any system.

Structural and biochemical analysis of DNA lesion-induced RNA polymerase II arrest

Transcription, catalyzed by RNA polymerase II (Pol II) in eukaryotes, is the first step in gene expression. RNA Pol II is a 12-subunit enzyme complex regulated by many different transcription factors during transcription initiation, elongation, and termination. During elongation, Pol II encounters various types of obstacles that can cause transcriptional pausing and arrest. Through decades of research on transcriptional pausing, it is widely known that Pol II can distinguish between different types of obstacles by its active site. A major class of obstacles is DNA lesions. While some DNA lesions can cause transient transcriptional pausing, which can be bypassed by Pol II itself or with the help from other elongation factors, bulky DNA damage can cause prolonged transcriptional pausing and arrest, which signals for transcription coupled repair. Using biochemical and structural biology approaches, the outcomes of many different types of DNA lesions, DNA modifications, and DNA binding molecules to transcription were studied. In this mini review, we will describe the in vitro transcription assays with Pol II to investigate the impacts of various DNA lesions on transcriptional outcomes and the crystallization method of lesion-arrested Pol II complex. These methods can provide a general platform for the structural and biochemical analysis of Pol II transcriptional pausing and bypass mechanisms.

An efficient Bayesian kinetic lumping algorithm to identify metastable conformational states via Gibbs sampling

Markov State Model (MSM) has become a popular approach to study the conformational dynamics of complex biological systems in recent years. Built upon a large number of short molecular dynamics simulation trajectories, MSM is able to predict the long time scale dynamics of complex systems. However, to achieve Markovianity, an MSM often contains hundreds or thousands of states (microstates), hindering human interpretation of the underlying system mechanism. One way to reduce the number of states is to lump kinetically similar states together and thus coarse-grain the microstates into macrostates. In this work, we introduce a probabilistic lumping algorithm, the Gibbs lumping algorithm, to assign a probability to any given kinetic lumping using the Bayesian inference. In our algorithm, the transitions among kinetically distinct macrostates are modeled by Poisson processes, which will well reflect the separation of time scales in the underlying free energy landscape of biomolecules. Furthermore, to facilitate the search for the optimal kinetic lumping (i.e., the lumped model with the highest probability), a Gibbs sampling algorithm is introduced. To demonstrate the power of our new method, we apply it to three systems: a 2D potential, alanine dipeptide, and a WW protein domain. In comparison with six other popular lumping algorithms, we show that our method can persistently produce the lumped macrostate model with the highest probability as well as the largest metastability. We anticipate that our Gibbs lumping algorithm holds great promise to be widely applied to investigate conformational changes in biological macromolecules.

Structural basis of RNA polymerase I stalling at UV light-induced DNA damage

RNA polymerase I (Pol I) transcribes ribosomal DNA (rDNA) to produce the ribosomal RNA (rRNA) precursor, which accounts for up to 60% of the total transcriptional activity in growing cells. Pol I monitors rDNA integrity and influences cell survival, but little is known about how this enzyme processes UV-induced lesions. We report the electron cryomicroscopy structure of Pol I in an elongation complex containing a cyclobutane pyrimidine dimer (CPD) at a resolution of 3.6 Å. The structure shows that the lesion induces an early translocation intermediate exhibiting unique features. The bridge helix residue Arg1015 plays a major role in CPD-induced Pol I stalling, as confirmed by mutational analysis. These results, together with biochemical data presented here, reveal the molecular mechanism of Pol I stalling by CPD lesions, which is distinct from Pol II arrest by CPD lesions. Our findings open the avenue to unravel the molecular mechanisms underlying cell endurance to lesions on rDNA.

Structural basis of DNA lesion recognition for eukaryotic transcription-coupled nucleotide excision repair

Eukaryotic transcription-coupled nucleotide excision repair (TC-NER) is a pathway that removes DNA lesions capable of blocking RNA polymerase II (Pol II) transcription from the template strand. This process is initiated by lesion-arrested Pol II and the recruitment of Cockayne Syndrome B protein (CSB). In this review, we will focus on the lesion recognition steps of eukaryotic TC-NER and summarize the recent research progress toward understanding the structural basis of Pol II-mediated lesion recognition and Pol II-CSB interactions. We will discuss the roles of CSB in both TC-NER initiation and transcription elongation. Finally, we propose an updated model of tripartite lesion recognition and verification for TC-NER in which CSB ensures Pol II-mediated recognition of DNA lesions for TC-NER.

An Introduction to the Structure and Function of the Catalytic Core Enzyme of Escherichia coli RNA Polymerase

RNA polymerase (RNAP) is the essential enzyme responsible for transcribing genetic information stored in DNA to RNA. Understanding the structure and function of RNAP is important for those who study basic principles in gene expression, such as the mechanism of transcription and its regulation, as well as translational sciences such as antibiotic development. With over a half-century of investigations, there is a wealth of information available on the structure and function of Escherichia coli RNAP. This review introduces the structural features of E. coli RNAP, organized by subunit, giving information on the function, location, and conservation of these features to early stage investigators who have just started their research of E. coli RNAP.

RNA-DNA and DNA-DNA base-pairing at the upstream edge of the transcription bubble regulate translocation of RNA polymerase and transcription rate

Translocation of RNA polymerase (RNAP) along DNA may be rate-limiting for transcription elongation. The Brownian ratchet model posits that RNAP rapidly translocates back and forth until the post-translocated state is stabilized by NTP binding. An alternative model suggests that RNAP translocation is slow and poorly reversible. To distinguish between these two models, we take advantage of an observation that pyrophosphorolysis rates directly correlate with the abundance of the pre-translocated fraction. Pyrophosphorolysis by RNAP stabilized in the pre-translocated state by bacteriophage HK022 protein Nun was used as a reference point to determine the pre-translocated fraction in the absence of Nun. The stalled RNAP preferentially occupies the post-translocated state. The forward translocation rate depends, among other factors, on melting of the RNA-DNA base pair at the upstream edge of the transcription bubble. DNA-DNA base pairing immediately upstream from the RNA-DNA hybrid stabilizes the post-translocated state. This mechanism is conserved between E. coli RNAP and S. cerevisiae RNA polymerase II and is partially dependent on the lid domain of the catalytic subunit. Thus, the RNA-DNA hybrid and DNA reannealing at the upstream edge of the transcription bubble emerge as targets for regulation of the transcription elongation rate.

Linking time-series of single-molecule experiments with molecular dynamics simulations by machine learning

Single-molecule experiments and molecular dynamics (MD) simulations are indispensable tools for investigating protein conformational dynamics. The former provide time-series data, such as donor-acceptor distances, whereas the latter give atomistic information, although this information is often biased by model parameters. Here, we devise a machine-learning method to combine the complementary information from the two approaches and construct a consistent model of conformational dynamics. It is applied to the folding dynamics of the formin-binding protein WW domain. MD simulations over 400 μs led to an initial Markov state model (MSM), which was then "refined" using single-molecule Förster resonance energy transfer (FRET) data through hidden Markov modeling. The refined or data-assimilated MSM reproduces the FRET data and features hairpin one in the transition-state ensemble, consistent with mutation experiments. The folding pathway in the data-assimilated MSM suggests interplay between hydrophobic contacts and turn formation. Our method provides a general framework for investigating conformational transitions in other proteins.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Improving the productivity of monodisperse polyhedral cages by the rational design of kinetic self-assembly pathways

Hollow polyhedral cages hold great potential for application in nanotechnological and biomedical fields. Understanding the formation mechanism of these self-assembled structures could provide guidance for the rational design of the desired polyhedral cages. Here, by constructing kinetic network models from extensive coarse-grained molecular dynamics simulations, we elucidated the formation mechanism of the dodecahedral cage, which is formed by the self-assembly of patchy particles. We found that the dodecahedral cage is formed through increasing the aggregate size followed by structure rearrangement. Based on this mechanistic understanding, we improved the productivity of the dodecahedral cage through the rational design of the patch arrangement of patchy particles, which promotes the structural rearrangement process. Our results demonstrate that it should be a feasible strategy to achieve the rational design of the desired nanostructures via the kinetic analysis. We anticipate that this methodology could be extended to other self-assembly systems for the fabrication of functional nanomaterials.

Harnessing complexity in molecular self-assembly using computer simulations

In molecular self-assembly, hundreds of thousands of freely-diffusing molecules associate to form ordered and functional architectures in the absence of an actuator. This intriguing phenomenon plays a critical role in biology and has become a powerful tool for the fabrication of advanced nanomaterials. Due to the limited spatial and temporal resolutions of current experimental techniques, computer simulations offer a complementary strategy to explore self-assembly with atomic resolution. Here, we review recent computational studies focusing on both thermodynamic and kinetic aspects. As we shall see, thermodynamic approaches based on modeling and statistical mechanics offer initial guidelines to design nanostructures with modest computational effort. Computationally more intensive analyses based on molecular dynamics simulations and kinetic network models (KNMs) reach beyond it, opening the door to the rational design of self-assembly pathways. Current limitations of these methodologies are discussed. We anticipate that the synergistic use of thermodynamic and kinetic analyses based on computer simulations will provide an important contribution to the de novo design of self-assembly.

Elucidation of the conformational dynamics of multi-body systems by construction of Markov state models

Constructing Markov State Models (MSMs) based on short molecular dynamics simulations is a powerful computational technique to complement experiments in predicting long-time kinetics of biomolecular processes at atomic resolution. Even though the MSM approach has been widely applied to study one-body processes such as protein folding and enzyme conformational changes, the majority of biological processes, e.g. protein-ligand recognition, signal transduction, and protein aggregation, essentially involve multiple entities. Here we review the attempts at constructing MSMs for multi-body systems, point out the challenges therein and discuss recent algorithmic progresses that alleviate these challenges. In particular, we describe an automatic kinetics based partitioning method that achieves optimal definition of the conformational states in a multi-body system, and discuss a novel maximum-likelihood approach that efficiently estimates the slow uphill kinetics utilizing pre-computed equilibrium populations of all states. We expect that these new algorithms and their combinations may boost investigations of important multi-body biological processes via the efficient construction of MSMs.

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

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

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.

Hinge action versus grip in translocation by RNA polymerase

Based on molecular dynamics simulations and functional studies, a conformational mechanism is posited for forward translocation by RNA polymerase (RNAP). In a simulation of a ternary elongation complex, the clamp and downstream cleft were observed to close. Hinges within the bridge helix and trigger loop supported generation of translocation force against the RNA-DNA hybrid resulting in opening of the furthest upstream i-8 RNA-DNA bp, establishing conditions for RNAP sliding. The β flap tip helix and the most N-terminal β' Zn finger engage the RNA, indicating a path of RNA threading out of the exit channel. Because the β flap tip connects to the RNAP active site through the β subunit double-Ψ-β-barrel and the associated sandwich barrel hybrid motif (also called the flap domain), the RNAP active site is coupled to the RNA exit channel and to the translocation of RNA-DNA. Using an exonuclease III assay to monitor translocation of RNAP elongation complexes, we show that K⁺ and Mg²⁺ and also an RNA 3'-OH or a 3'-H₂ affect RNAP sliding. Because RNAP grip to template suggests a sticky translocation mechanism, and because grip is enhanced by increasing K⁺ and Mg²⁺concentration, biochemical assays are consistent with a conformational change that drives forward translocation as observed in simulations. Mutational analysis of the bridge helix indicates that 778-GARKGL-783 (Escherichia coli numbering) is a homeostatic hinge that undergoes multiple bends to compensate for complex conformational dynamics during phosphodiester bond formation and translocation.

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.

Nucleosome mobility and the regulation of gene expression: Insights from single-molecule studies

Nucleosomes at the promoters of genes regulate the accessibility of the transcription machinery to DNA, and function as a basic layer in the complex regulation of gene expression. Our understanding of the role of the nucleosome's spontaneous, thermally driven position changes in modulating expression is lacking. This is the result of the paucity of experimental data on these dynamics, at high-resolution, and for DNA sequences that belong to real, transcribed genes. We have developed an assay that uses partial, reversible unzipping of nucleosomes with optical tweezers to repeatedly probe a nucleosome's position over time. Using the nucleosomes at the promoters of two model genes, Cga and Lhb, we show that the mobility of nucleosomes is modulated by the sequence of DNA and by the use of alternative histone variants, and describe how the mobility can affect transcription, at the initiation and elongation phases.

Mechanism of transcription-coupled DNA modification recognition

As a key enzyme for gene expression, RNA polymerase II (pol II) reads along the DNA template and catalyzes accurate mRNA synthesis during transcription. On the other hand, genomic DNA is under constant attack by endogenous and environmental stresses. These attack cause many DNA lesions. Pol II functions as a specific sensor that is able to recognize changes in DNA sequences and structures and induces different outcomes. A critical question in the field is how Pol II recognizes and senses these DNA modifications or lesions. Recent studies provided new insights into understanding this critical question. In this mini-review, we would like to focus on three classes of DNA lesions/modifications: (1) Bulky, DNA-distorting lesions that block pol II transcription, (2) small DNA lesions that promote pol II pausing and error-prone transcriptional bypass, and (3) endogenous enzyme-catalyzed DNA modifications that lead to pol II pausing and error-free transcriptional bypass.

Elucidating Mechanisms of Molecular Recognition Between Human Argonaute and miRNA Using Computational Approaches

MicroRNA (miRNA) and Argonaute (AGO) protein together form the RNA-induced silencing complex (RISC) that plays an essential role in the regulation of gene expression. Elucidating the underlying mechanism of AGO-miRNA recognition is thus of great importance not only for the in-depth understanding of miRNA function but also for inspiring new drugs targeting miRNAs. In this chapter we introduce a combined computational approach of molecular dynamics (MD) simulations, Markov state models (MSMs), and protein-RNA docking to investigate AGO-miRNA recognition. Constructed from MD simulations, MSMs can elucidate the conformational dynamics of AGO at biologically relevant timescales. Protein-RNA docking can then efficiently identify the AGO conformations that are geometrically accessible to miRNA. Using our recent work on human AGO2 as an example, we explain the rationale and the workflow of our method in details. This combined approach holds great promise to complement experiments in unraveling the mechanisms of molecular recognition between large, flexible, and complex biomolecules.

High-Resolution Phenotypic Landscape of the RNA Polymerase II Trigger Loop

The active sites of multisubunit RNA polymerases have a “trigger loop” (TL) that multitasks in substrate selection, catalysis, and translocation. To dissect the Saccharomyces cerevisiae RNA polymerase II TL at individual-residue resolution, we quantitatively phenotyped nearly all TL single variants en masse. Three mutant classes, revealed by phenotypes linked to transcription defects or various stresses, have distinct distributions among TL residues. We find that mutations disrupting an intra-TL hydrophobic pocket, proposed to provide a mechanism for substrate-triggered TL folding through destabilization of a catalytically inactive TL state, confer phenotypes consistent with pocket disruption and increased catalysis. Furthermore, allele-specific genetic interactions among TL and TL-proximal domain residues support the contribution of the funnel and bridge helices (BH) to TL dynamics. Our structural genetics approach incorporates structural and phenotypic data for high-resolution dissection of transcription mechanisms and their evolution, and is readily applicable to other essential yeast proteins.

Proper regulation of Pol II transcription, the first step of gene expression, is essential for life. Extensive evidence has revealed a widely conserved and dynamic polymerase active site component, termed the Trigger Loop (TL), in balancing transcription rate and fidelity while possibly allowing control of transcription elongation. Coupling high-throughput sequencing with our previously established genetic system, we are able to assess the in vivo phenotypes for almost all possible single substitution Pol II TL mutants in the budding yeast Saccharomyces cerevisiae. We show that mutants in the TL nucleotide interacting and linker regions widely confer dominant and severe growth defects. Clustering of TL mutants’ transcription-related and general stress phenotypes reveals three main classes of TL mutants, including previously identified fast and slow elongating mutants. Comprehensive analyses of the distribution of fast and slow elongation mutants in light of existing Pol II crystal structures reveal critical regions contributing to proper TL dynamics and function. Evidence is presented linking a previously observed hydrophobic pocket to NTP substrate-induced TL closing, the mechanism critical for correct substrates selection and transcription fidelity. Finally, we assess the functional interplay between TL and its proximal domains, and their presumptive roles in the function and evolution of the TL. Utilizing the Pol II TL as a case study, we present a structural genetics approach that reveals insights into a complex, multi-functional, and essential domain in yeast.

Adaptive partitioning by local density-peaks: An efficient density-based clustering algorithm for analyzing molecular dynamics trajectories

We present an efficient density-based adaptive-resolution clustering method APLoD for analyzing large-scale molecular dynamics (MD) trajectories. APLoD performs the k-nearest-neighbors search to estimate the density of MD conformations in a local fashion, which can group MD conformations in the same high-density region into a cluster. APLoD greatly improves the popular density peaks algorithm by reducing the running time and the memory usage by 2-3 orders of magnitude for systems ranging from alanine dipeptide to a 370-residue Maltose-binding protein. In addition, we demonstrate that APLoD can produce clusters with various sizes that are adaptive to the underlying density (i.e., larger clusters at low-density regions, while smaller clusters at high-density regions), which is a clear advantage over other popular clustering algorithms including k-centers and k-medoids. We anticipate that APLoD can be widely applied to split ultra-large MD datasets containing millions of conformations for subsequent construction of Markov State Models. © 2016 Wiley Periodicals, Inc.

Molecular Structures of Transcribing RNA Polymerase I

RNA polymerase I (Pol I) is a 14-subunit enzyme that solely synthesizes pre-ribosomal RNA. Recently, the crystal structure of apo Pol I gave unprecedented insight into its molecular architecture. Here, we present three cryo-EM structures of elongating Pol I, two at 4.0 Å and one at 4.6 Å resolution, and a Pol I open complex at 3.8 Å resolution. Two modules in Pol I mediate the narrowing of the DNA-binding cleft by closing the clamp domain. The DNA is bound by the clamp head and by the protrusion domain, allowing visualization of the upstream and downstream DNA duplexes in one of the elongation complexes. During formation of the Pol I elongation complex, the bridge helix progressively folds, while the A12.2 C-terminal domain is displaced from the active site. Our results reveal the conformational changes associated with elongation complex formation and provide additional insight into the Pol I transcription cycle.

•

Pol I gradually closes its DNA-binding cleft during elongation complex formation

•

Pol I bridge helix folds into a conformation similar to that in other RNA polymerases

•

During elongation, the Pol I A12.2 C-terminal domain is excluded from the cleft

•

A49 tandem winged helix domain contacts upstream DNA, similar to TFIIE

Understanding the core of RNA interference: The dynamic aspects of Argonaute-mediated processes

At the core of RNA interference, the Argonaute proteins (Ago) load and utilize small guide nucleic acids to silence mRNAs or cleave foreign nucleic acids in a sequence specific manner. In recent years, based on extensive structural studies of Ago and its interaction with the nucleic acids, considerable progress has been made to reveal the dynamic aspects of various Ago-mediated processes. Here we review these novel insights into the guide-strand loading, duplex unwinding, and effects of seed mismatch, with a focus on two representative Agos, the human Ago 2 (hAgo2) and the bacterial Thermus thermophilus Ago (TtAgo). In particular, comprehensive molecular simulation studies revealed that although sharing similar overall structures, the two Agos have vastly different conformational landscapes and guide-strand loading mechanisms because of the distinct rigidity of their L1-PAZ hinge. Given the central role of the PAZ motions in regulating the exposure of the nucleic acid binding channel, these findings exemplify the importance of protein motions in distinguishing the overlapping, yet distinct, mechanisms of Ago-mediated processes in different organisms.

A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases

The enzymatic polymerization of DNA and RNA is the basis for genetic inheritance for all living organisms. It is catalyzed by the DNA/RNA polymerase (Pol) superfamily. Here, bioinformatics analysis reveals that the incoming nucleotide substrate always forms an H-bond between its 3'-OH and β-phosphate moieties upon formation of the Michaelis complex. This previously unrecognized H-bond implies a novel self-activated mechanism (SAM), which synergistically connects the in situ nucleophile formation with subsequent nucleotide addition and, importantly, nucleic acid translocation. Thus, SAM allows an elegant and efficient closed-loop sequence of chemical and physical steps for Pol catalysis. This is markedly different from previous mechanistic hypotheses. Our proposed mechanism is corroborated via ab initio QM/MM simulations on a specific Pol, the human DNA polymerase-η, an enzyme involved in repairing damaged DNA. The structural conservation of DNA and RNA Pols supports the possible extension of SAM to Pol enzymes from the three domains of life.

Structure of the initiation-competent RNA polymerase I and its implication for transcription

Eukaryotic RNA polymerase I (Pol I) is specialized in rRNA gene transcription synthesizing up to 60% of cellular RNA. High level rRNA production relies on efficient binding of initiation factors to the rRNA gene promoter and recruitment of Pol I complexes containing initiation factor Rrn3. Here, we determine the cryo-EM structure of the Pol I-Rrn3 complex at 7.5 Å resolution, and compare it with Rrn3-free monomeric and dimeric Pol I. We observe that Rrn3 contacts the Pol I A43/A14 stalk and subunits A190 and AC40, that association re-organizes the Rrn3 interaction interface, thereby preventing Pol I dimerization; and Rrn3-bound and monomeric Pol I differ from the dimeric enzyme in cleft opening, and localization of the A12.2 C-terminus in the active centre. Our findings thus support a dual role for Rrn3 in transcription initiation to stabilize a monomeric initiation competent Pol I and to drive pre-initiation complex formation.

Eukaryotic RNA polymerase I (Pol I) is responsible for the transcription of rRNA genes. Here the authors determine the cryo-EM structure of the Pol I-Rrn3 complex, providing insight into how Rrn3 stabilizes the monomeric initiation competent Pol I to drive pre-initiation complex formation.

Nanoscale Engineering of Designer Cellulosomes

Biocatalysts showcase the upper limit obtainable for high-speed molecular processing and transformation. Efforts to engineer functionality in synthetic nanostructured materials are guided by the increasing knowledge of evolving architectures, which enable controlled molecular motion and precise molecular recognition. The cellulosome is a biological nanomachine, which, as a fundamental component of the plant-digestion machinery from bacterial cells, has a key potential role in the successful development of environmentally-friendly processes to produce biofuels and fine chemicals from the breakdown of biomass waste. Here, the progress toward so-called "designer cellulosomes", which provide an elegant alternative to enzyme cocktails for lignocellulose breakdown, is reviewed. Particular attention is paid to rational design via computational modeling coupled with nanoscale characterization and engineering tools. Remaining challenges and potential routes to industrial application are put forward.

RNA polymerase II acts as a selective sensor for DNA lesions and endogenous DNA modifications

During transcription elongation, RNA polymerase II (pol II) travels along the DNA template across thousands to millions of nucleotides and accurately synthesizes the complementary RNA transcripts. Apart from its canonical function as a key enzyme for DNA-dependent RNA synthesis, pol II also functions as a selective sensor to recognize DNA lesions or epigenetic modifications.

Multiensemble Markov models of molecular thermodynamics and kinetics

We introduce the general transition-based reweighting analysis method (TRAM), a statistically optimal approach to integrate both unbiased and biased molecular dynamics simulations, such as umbrella sampling or replica exchange. TRAM estimates a multiensemble Markov model (MEMM) with full thermodynamic and kinetic information at all ensembles. The approach combines the benefits of Markov state models-clustering of high-dimensional spaces and modeling of complex many-state systems-with those of the multistate Bennett acceptance ratio of exploiting biased or high-temperature ensembles to accelerate rare-event sampling. TRAM does not depend on any rate model in addition to the widely used Markov state model approximation, but uses only fundamental relations such as detailed balance and binless reweighting of configurations between ensembles. Previous methods, including the multistate Bennett acceptance ratio, discrete TRAM, and Markov state models are special cases and can be derived from the TRAM equations. TRAM is demonstrated by efficiently computing MEMMs in cases where other estimators break down, including the full thermodynamics and rare-event kinetics from high-dimensional simulation data of an all-atom protein-ligand binding model.

Enhancing pairwise state-transition weights: A new weighting scheme in simulated tempering that can minimize transition time between a pair of conformational states

Simulated tempering (ST) is a widely used enhancing sampling method for Molecular Dynamics simulations. As one expanded ensemble method, ST is a combination of canonical ensembles at different temperatures and the acceptance probability of cross-temperature transitions is determined by both the temperature difference and the weights of each temperature. One popular way to obtain the weights is to adopt the free energy of each canonical ensemble, which achieves uniform sampling among temperature space. However, this uniform distribution in temperature space may not be optimal since high temperatures do not always speed up the conformational transitions of interest, as anti-Arrhenius kinetics are prevalent in protein and RNA folding. Here, we propose a new method: Enhancing Pairwise State-transition Weights (EPSW), to obtain the optimal weights by minimizing the round-trip time for transitions among different metastable states at the temperature of interest in ST. The novelty of the EPSW algorithm lies in explicitly considering the kinetics of conformation transitions when optimizing the weights of different temperatures. We further demonstrate the power of EPSW in three different systems: a simple two-temperature model, a two-dimensional model for protein folding with anti-Arrhenius kinetics, and the alanine dipeptide. The results from these three systems showed that the new algorithm can substantially accelerate the transitions between conformational states of interest in the ST expanded ensemble and further facilitate the convergence of thermodynamics compared to the widely used free energy weights. We anticipate that this algorithm is particularly useful for studying functional conformational changes of biological systems where the initial and final states are often known from structural biology experiments.

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.

Elucidation of the Dynamics of Transcription Elongation by RNA Polymerase II using Kinetic Network Models

RNA polymerase II (Pol II) is an essential enzyme that catalyzes transcription with high efficiency and fidelity in eukaryotic cells. During transcription elongation, Pol II catalyzes the nucleotide addition cycle (NAC) to synthesize mRNA using DNA as the template. The transitions between the states of the NAC require conformational changes of both the protein and nucleotides. Although X-ray structures are available for most of these states, the dynamics of the transitions between states are largely unknown. Molecular dynamics (MD) simulations can predict structure-based molecular details and shed light on the mechanisms of these dynamic transitions. However, the employment of MD simulations on a macromolecule (tens to hundreds of nanoseconds) such as Pol II is challenging due to the difficulty of reaching biologically relevant timescales (tens of microseconds or even longer). For this challenge to be overcome, kinetic network models (KNMs), such as Markov State Models (MSMs), have become a popular approach to access long-timescale conformational changes using many short MD simulations. We describe here our application of KNMs to characterize the molecular mechanisms of the NAC of Pol II. First, we introduce the general background of MSMs and further explain procedures for the construction and validation of MSMs by providing some technical details. Next, we review our previous studies in which we applied MSMs to investigate the individual steps of the NAC, including translocation and pyrophosphate ion release. In particular, we describe in detail how we prepared the initial conformations of Pol II elongation complex, performed MD simulations, extracted MD conformations to construct MSMs, and further validated them. We also summarize our major findings on molecular mechanisms of Pol II elongation based on these MSMs. In addition, we have included discussions regarding various key points and challenges for applications of MSMs to systems as large as the Pol II elongation complex. Finally, to study the overall NAC, we combine the individual steps of the NAC into a five-state KNM based on a nonbranched Brownian ratchet scheme to explain the single-molecule optical tweezers experimental data. The studies complement experimental observations and provide molecular mechanisms for the transcription elongation cycle. In the long term, incorporation of sequence-dependent kinetic parameters into KNMs has great potential for identifying error-prone sequences and predicting transcription dynamics in genome-wide transcriptomes.

Unfolding mechanism of thrombin-binding aptamer revealed by molecular dynamics simulation and Markov State Model

Thrombin-binding aptamer (TBA) with the sequence 5′GGTTGGTGTGGTTGG3′ could fold into G-quadruplex, which correlates with functionally important genomic regionsis. However, unfolding mechanism involved in the structural stability of G-quadruplex has not been satisfactorily elucidated on experiments so far. Herein, we studied the unfolding pathway of TBA by a combination of molecular dynamics simulation (MD) and Markov State Model (MSM). Our results revealed that the unfolding of TBA is not a simple two-state process but proceeds along multiple pathways with multistate intermediates. One high flux confirms some observations from NMR experiment. Another high flux exhibits a different and simpler unfolding pathway with less intermediates. Two important intermediate states were identified. One is similar to the G-triplex reported in the folding of G-quadruplex, but lack of H-bonding between guanines in the upper plane. More importantly, another intermediate state acting as a connector to link the folding region and the unfolding one, was the first time identified, which exhibits higher population and stability than the G-triplex-like intermediate. These results will provide valuable information for extending our understanding the folding landscape of G-quadruplex formation.

A Flexible Domain-Domain Hinge Promotes an Induced-fit Dominant Mechanism for the Loading of Guide-DNA into Argonaute Protein in Thermus thermophilus

Argonaute proteins (Ago) are core components of the RNA Induced Silencing Complex (RISC) that load and utilize small guide nucleic acids to silence mRNAs or cleave foreign DNAs. Despite the essential role of Ago in gene regulation and defense against virus, the molecular mechanism of guide-strand loading into Ago remains unclear. We explore such a mechanism in the bacterium Thermus thermophilus Ago (TtAgo), via a computational approach combining molecular dynamics, bias-exchange metadynamics, and protein-DNA docking. We show that apo TtAgo adopts multiple closed states that are unable to accommodate guide-DNA. Conformations able to accommodate the guide are beyond the reach of thermal fluctuations from the closed states. These results suggest an induced-fit dominant mechanism for guide-strand loading in TtAgo, drastically different from the two-step mechanism for human Ago 2 (hAgo2) identified in our previous study. Such a difference between TtAgo and hAgo2 is found to mainly originate from the distinct rigidity of their L1-PAZ hinge. Further comparison among known Ago structures from various species indicates that the L1-PAZ hinge may be flexible in general for prokaryotic Ago's but rigid for eukaryotic Ago's.