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

Chromatin folding through nonuniform motorization by responsive motor proteins

Chromatin is partially structured through the effects of biological motors. "Swimming motors" such as RNA polymerases and chromatin remodelers are thought to act differentially on the active parts of the genome and the stored inactive part. By systematically expanding the many-body master equation for chromosomes driven by swimming motors, we show that this nonuniform aspect of motorization leads to heterogeneously folded conformations, thereby contributing to chromosome compartmentalization.

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.

Motorized chain models of the ideal chromosome

An array of motor proteins consumes chemical energy in setting up the architectures of chromosomes. Here, we explore how the structure of ideal polymer chains is influenced by two classes of motors. The first class which we call "swimming motors" acts to propel the chromatin fiber through three-dimensional space. They represent a caricature of motors such as RNA polymerases. Previously, they have often been described by adding a persistent flow onto Brownian diffusion of the chain. The second class of motors, which we call "grappling motors" caricatures the loop extrusion processes in which segments of chromatin fibers some distance apart are brought together. We analyze these models using a self-consistent variational phonon approximation to a many-body Master equation incorporating motor activities. We show that whether the swimming motors lead to contraction or expansion depends on the susceptibility of the motors, that is, how their activity depends on the forces they must exert. Grappling motors in contrast to swimming motors lead to long-ranged correlations that resemble those first suggested for fractal globules and that are consistent with the effective interactions inferred by energy landscape analyses of Hi-C data on the interphase chromosome.

Development of a capillary zone electrophoresis method to monitor magnesium ion consumption during in vitro transcription for mRNA production

Messenger RNA (mRNA) vaccines represent a landmark in vaccinology, especially with their success in COVID-19 vaccines, which have shown great promise for future vaccine development and disease prevention. As a platform technology, synthetic mRNA can be produced with high fidelity using in vitro transcription (IVT). Magnesium plays a vital role in the IVT process, facilitating the phosphodiester bond formation between adjacent nucleotides and ensuring accurate transcription to produce high-quality mRNA. The development of the IVT process has prompted key inquiries about in-process characterization of magnesium ion (Mg++) consumption, relating to the RNA polymerase (RNAP) activation, fed-batch mode production yield, and mRNA quality. Hence, it becomes crucial to monitor the free Mg++ concentration throughout the IVT process. However, no free Mg++ analysis method has been reported for complex IVT reactions. Here we report a robust capillary zone electrophoresis (CZE) method with indirect UV detection. The assay allows accurate quantitation of free Mg++ for the complex IVT reaction where it is essential to preserve IVT samples in their native-like state during analysis to avoid dissociation of bound Mg complexes. By applying this CZE method, the relationships between free Mg++ concentration, the mRNA yield, and dsRNA impurity level were investigated. Such mechanistic understanding facilitates informed decisions regarding the quantity and timing of feeding starting materials to increase the yield. Furthermore, this approach can serve as a platform method for analyzing the free Mg++ in complex sample matrices where preserving the native-like state of Mg++ binding is key for accurate quantitation.

Systematic identification and repurposing of FDA-approved drugs as antibacterial agents against Streptococcus pyogenes: In silico and in vitro studies

Streptococcus pyogenes (Group A Streptococcus - GAS) is a human pathogen causing wide range of infections and toxin-mediated diseases in human beings of all age groups with fatality of 10-30 %. The limited success of antibiotics and the non-availability of vaccines makes GAS a global burden. The multi-subunit RNA polymerase (RNAP) is a validated bacterial therapeutic target as it is involved in transcription and can arrest growth. Of the five subunits of this enzyme complex, the β-subunit (RpoC) has attracted specific attention as a drug target, particularly in the switch region. Here we attempt to repurpose non-antimicrobial drugs to act as RpoC inhibitors against S. pyogenes. In this study, 1826 FDA approved drugs have been identified through high-throughput virtual screening. Free Energy Perturbation (FEP) based binding free energy calculations have been performed at the final step of the virtual screening funnel to ensure high accuracy in silico results. Three compounds identified have been tested for susceptibility of S. pyogenes MTCC 442 strain and two antibiotic-resistant clinical isolates of S. pyogenes using microdilution assay. Among the three, two drugs Amlodipine Besylate (Amd) and Ranitidine hydrochloride (Rnt) have shown inhibition against all the tested strains and its mechanism of interaction with RpoC has been studied. The docked complexes were analyzed to understand the binding mode of the drugs to the target. Classical Molecular Dynamics studies for RpoC-Rnt complex and the two stable conformations of RpoC-Amd complex was carried out. Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of Gyration (RoG) and Solvent Accessible Surface Area (SASA) of the complexes were plotted and studied. The thermodynamic parameters of protein-drug were experimentally determined using Isothermal Titration Calorimetry (ITC). Infrared spectroscopic studies and Fluorescence quenching studies provided insights into the secondary structural changes in RpoC on binding to the drugs.

DNA opening during transcription initiation by RNA polymerase II in atomic detail

RNA polymerase II (RNAP II) synthesizes RNA by reading the DNA code. During transcription initiation, RNAP II opens the double-stranded DNA to expose the DNA template to the active site. The molecular interactions driving and controlling DNA opening are not well understood. We used all-atom steered molecular dynamics simulations to derive a continuous pathway of DNA opening in human RNAP II, involving a 55 Å DNA strand displacement and a nearly 360° DNA helix rotation. To drive such large-scale transitions, we used a combination of RMSD-based collective variables, a newly designed rotational coordinate, and a path collective variable. The simulations reveal extensive interactions of the DNA with three conserved protein loops near the active site, namely with the rudder, fork loop 1, and fork loop 2. According to the simulations, DNA-protein interactions support DNA opening by a twofold mechanism; they catalyze DNA opening by attacking Watson-Crick hydrogen bonds, and they stabilize the open DNA bubble by the formation of a wide set of DNA-protein salt bridges.

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.

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 cutting edge of archaeal transcription

The archaeal RNA polymerase (RNAP) is a double-psi β-barrel enzyme closely related to eukaryotic RNAPII in terms of subunit composition and architecture, promoter elements and basal transcription factors required for the initiation and elongation phase of transcription. Understanding archaeal transcription is, therefore, key to delineate the universally conserved fundamental mechanisms of transcription as well as the evolution of the archaeo-eukaryotic transcription machineries. The dynamic interplay between RNAP subunits, transcription factors and nucleic acids dictates the activity of RNAP and ultimately gene expression. This review focusses on recent progress in our understanding of (i) the structure, function and molecular mechanisms of known and less characterized factors including Elf1 (Elongation factor 1), NusA (N-utilization substance A), TFS4, RIP and Eta, and (ii) their evolution and phylogenetic distribution across the expanding tree of Archaea.

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.

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.

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.

The Top1 paradox: Friend and foe of the eukaryotic genome

Topoisomerases manage the torsional stress associated with the separation of DNA strands during transcription and DNA replication. Eukaryotic Topoisomerase I (Top1) is a Type IB enzyme that nicks and rejoins only one strand of duplex DNA, and it is especially important during transcription. By resolving transcription-associated torsional stress, Top1 reduces the accumulation of genome-destabilizing R-loops and non-B DNA structures. The DNA nicking activity of Top1, however, can also initiate genome instability in the form of illegitimate recombination, homologous recombination and mutagenesis. In this review, we focus on the diverse, and often opposing, roles of Top1 in regulating eukaryotic genome stability.

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.

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.