Intrinsic Cleavage of RNA Polymerase II Adopts a Nucleobase-independent Mechanism Assisted by Transcript Phosphate

Enzymatic-related network of catalysis, polyamine, and tumors for acetylpolyamine oxidase: from calculation to experiment

The regulation of enzymes and development of polyamine analogs capable of controlling the dynamics of endogenous polyamines to achieve anti-tumor effects is one of the biggest challenges in polyamine research. However, the root of the problem remains unsolved. This study represents a significant milestone as it unveils, for the first time, the comprehensive catalytic map of acetylpolyamine oxidase that includes chemical transformation and product release kinetics, by utilizing multiscale simulations with over six million dynamical snapshots. The transportation of acetylspermine is strongly exothermic, and high binding affinity of enzyme and reactant is observed. The transfer of hydride from polyamine to FAD is the rate-limiting step, via an H-shift coupled electron transfer mechanism. The two products are released in a detour stepwise mechanism, which also impacts the enzymatic efficiency. Inspired by these mechanistic insights into enzymatic catalysis, we propose a novel strategy that regulates the polyamine level and catalytic progress through the action of His64. Directly suppressing APAO by mutating His64 further inhibited growth and migration of tumor cells and tumor tissue in vitro and in vivo. Therefore, the network connecting microcosmic and macroscopic scales opens up new avenues for designing polyamine compounds and conducting anti-tumor research in the future.

Chemically Modified Platforms for Better RNA Therapeutics

RNA-based therapies have catalyzed a revolutionary transformation in the biomedical landscape, offering unprecedented potential in disease prevention and treatment. However, despite their remarkable achievements, these therapies encounter substantial challenges including low stability, susceptibility to degradation by nucleases, and a prominent negative charge, thereby hindering further development. Chemically modified platforms have emerged as a strategic innovation, focusing on precise alterations either on the RNA moieties or their associated delivery vectors. This comprehensive review delves into these platforms, underscoring their significance in augmenting the performance and translational prospects of RNA-based therapeutics. It encompasses an in-depth analysis of various chemically modified delivery platforms that have been instrumental in propelling RNA therapeutics toward clinical utility. Moreover, the review scrutinizes the rationale behind diverse chemical modification techniques aiming at optimizing the therapeutic efficacy of RNA molecules, thereby facilitating robust disease management. Recent empirical studies corroborating the efficacy enhancement of RNA therapeutics through chemical modifications are highlighted. Conclusively, we offer profound insights into the transformative impact of chemical modifications on RNA drugs and delineates prospective trajectories for their future development and clinical integration.

Nucleotide addition and cleavage by RNA polymerase II: Coordination of two catalytic reactions using a single active site

RNA polymerase II (Pol II) incorporates complementary ribonucleotides into the growing RNA chain one at a time via the nucleotide addition cycle. The nucleotide addition cycle, however, is prone to misincorporation of noncomplementary nucleotides. Thus, to ensure transcriptional fidelity, Pol II backtracks and then cleaves the misincorporated nucleotides. These two reverse reactions, nucleotide addition and cleavage, are catalyzed in the same active site of Pol II, which is different from DNA polymerases or other endonucleases. Recently, substantial progress has been made to understand how Pol II effectively performs its dual role in the same active site. Our review highlights these recent studies and provides an overall model of the catalytic mechanisms of Pol II. In particular, RNA extension follows the two-metal-ion mechanism, and several Pol II residues play important roles to facilitate the catalysis. In sharp contrast, the cleavage reaction is independent of any Pol II residues. Interestingly, Pol II relies on its residues to recognize the misincorporated nucleotides during the backtracking process, prior to cleavage. In this way, Pol II efficiently compartmentalizes its two distinct catalytic functions using the same active site. Lastly, we also discuss a new perspective on the potential third Mg²⁺ in the nucleotide addition and intrinsic cleavage reactions.

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.

1'-Ribose cyano substitution allows Remdesivir to effectively inhibit nucleotide addition and proofreading during SARS-CoV-2 viral RNA replication

COVID-19 has recently caused a global health crisis and an effective interventional therapy is urgently needed. Remdesivir is one effective inhibitor for SARS-CoV-2 viral RNA replication. It supersedes other NTP analogues because it not only terminates the polymerization activity of RNA-dependent RNA polymerase (RdRp), but also inhibits the proofreading activity of intrinsic exoribonuclease (ExoN). Even though the static structure of Remdesivir binding to RdRp has been solved and biochemical experiments have suggested it to be a "delayed chain terminator", the underlying molecular mechanisms is not fully understood. Here, we performed all-atom molecular dynamics (MD) simulations with an accumulated simulation time of 24 microseconds to elucidate the inhibitory mechanism of Remdesivir on nucleotide addition and proofreading. We found that when Remdesivir locates at an upstream site in RdRp, the 1'-cyano group experiences electrostatic interactions with a salt bridge (Asp865-Lys593), which subsequently halts translocation. Our findings can supplement the current understanding of the delayed chain termination exerted by Remdesivir and provide an alternative molecular explanation about Remdesivir's inhibitory mechanism. Such inhibition also reduces the likelihood of Remdesivir to be cleaved by ExoN acting on 3'-terminal nucleotides. Furthermore, our study also suggests that Remdesivir's 1'-cyano group can disrupt the cleavage site of ExoN via steric interactions, leading to a further reduction in the cleavage efficiency. Our work provides plausible and novel mechanisms at the molecular level of how Remdesivir inhibits viral RNA replication, and our findings may guide rational design for new treatments of COVID-19 targeting viral replication.

Catalytic Mechanism of Non-Target DNA Cleavage in CRISPR-Cas9 Revealed by Ab Initio Molecular Dynamics

CRISPR-Cas9 is a cutting-edge genome editing technology, which uses the endonuclease Cas9 to introduce mutations at desired sites of the genome. This revolutionary tool is promising to treat a myriad of human genetic diseases. Nevertheless, the molecular basis of DNA cleavage, which is a fundamental step for genome editing, has not been established. Here, quantum-classical molecular dynamics (MD) and free energy methods are used to disclose the two-metal-dependent mechanism of phosphodiester bond cleavage in CRISPR-Cas9. Ab initio MD reveals a conformational rearrangement of the Mg2+-bound RuvC active site, which entails the relocation of H983 to act as a general base. Then, the DNA cleavage proceeds through a concerted associative pathway fundamentally assisted by the joint dynamics of the two Mg2+ ions. This clarifies previous controversial experimental evidence, which could not fully establish the catalytic role of the conserved H983 and the metal cluster conformation. The comparison with other two-metal-dependent enzymes supports the identified mechanism and suggests a common catalytic strategy for genome editing and recombination. Overall, the non-target DNA cleavage catalysis described here resolves a fundamental open question in the CRISPR-Cas9 biology and provides valuable insights for improving the catalytic efficiency and the metal-dependent function of the Cas9 enzyme, which are at the basis of the development of genome editing tools.

Elucidating molecular mechanisms of functional conformational changes of proteins via Markov state models

Functional conformational changes of proteins can facilitate numerous biological events in cells. The Markov state model (MSM) built from molecular dynamics simulations provide a powerful approach to study them. We here introduce a protocol that is tailor-made for constructing MSMs to study the functional conformational changes of proteins. In this protocol, one of the important steps is to select proper molecular features that can collectively describe the slowest timescales of conformational changes of interest. We recommend spectral oASIS, the modified version of oASIS, as a promising approach for automatic feature selection. Recently developed deep learning methods could also serve efficient approaches for selecting features and finding collective variables. Using DNA repair enzymes and RNA polymerases as examples, we review recent applications of MSMs to elucidate molecular mechanisms of functional conformational changes. Finally, we discuss remaining challenges and future perspectives for constructing MSMs to study functional conformational changes of proteins.

Identifying Transcription Error-Enriched Genomic Loci Using Nuclear Run-on Circular-Sequencing Coupled with Background Error Modeling

RNA polymerase transcribes certain genomic loci with higher errors rates. These transcription error-enriched genomic loci (TEELs) have implications in disease. Current deep-sequencing methods cannot distinguish TEELs from post-transcriptional modifications, stochastic transcription errors, and technical noise, impeding efforts to elucidate the mechanisms linking TEELs to disease. Here, we describe background error model-coupled precision nuclear run-on circular-sequencing (EmPC-seq) to discern genomic regions enriched for transcription misincorporations. EmPC-seq innovatively combines a nuclear run-on assay for capturing nascent RNA before post-transcriptional modifications, a circular-sequencing step that sequences the same nascent RNA molecules multiple times to improve accuracy, and a statistical model for distinguishing error-enriched regions among stochastic polymerase errors. Applying EmPC-seq to the ribosomal RNA transcriptome, we show that TEELs of RNA polymerase I are not randomly distributed but clustered together, with higher error frequencies at nascent transcript 3' ends. Our study establishes a reliable method of identifying TEELs with nucleotide precision, which can help elucidate their molecular origins.

RNA polymerase II stalls on oxidative DNA damage via a torsion-latch mechanism involving lone pair-π and CH-π interactions

Oxidation of guanine generates several types of DNA lesions, such as 8-oxoguanine (8OG), 5-guanidinohydantoin (Gh), and spiroiminodihydantoin (Sp). These guanine-derived oxidative DNA lesions interfere with both replication and transcription. However, the molecular mechanism of transcription processing of Gh and Sp remains unknown. In this study, by combining biochemical and structural analysis, we revealed distinct transcriptional processing of these chemically related oxidized lesions: 8OG allows both error-free and error-prone bypass, whereas Gh or Sp causes strong stalling and only allows slow error-prone incorporation of purines. Our structural studies provide snapshots of how polymerase II (Pol II) is stalled by a nonbulky Gh lesion in a stepwise manner, including the initial lesion encounter, ATP binding, ATP incorporation, jammed translocation, and arrested states. We show that while Gh can form hydrogen bonds with adenosine monophosphate (AMP) during incorporation, this base pair hydrogen bonding is not sufficient to hold an ATP substrate in the addition site and is not stable during Pol II translocation after the chemistry step. Intriguingly, we reveal a unique structural reconfiguration of the Gh lesion in which the hydantoin ring rotates ∼90° and is perpendicular to the upstream base pair planes. The perpendicular hydantoin ring of Gh is stabilized by noncanonical lone pair-π and CH-π interactions, as well as hydrogen bonds. As a result, the Gh lesion, as a functional mimic of a 1,2-intrastrand crosslink, occupies canonical -1 and +1 template positions and compromises the loading of the downstream template base. Furthermore, we suggest Gh and Sp lesions are potential targets of transcription-coupled repair.

Cleaning Up Mechanistic Debris Generated by Twister Ribozymes Using Computational RNA Enzymology

The catalytic properties of RNA have been a subject of fascination and intense research since their discovery over 30 years ago. Very recently, several classes of nucleolytic ribozymes have emerged and been characterized structurally. Among these, the twister ribozyme has been center-stage, and a topic of debate about its architecture and mechanism owing to conflicting interpretations of different crystal structures, and in some cases conflicting interpretations of the same functional data. In the present work, we attempt to clean up the mechanistic "debris" generated by twister ribozymes using a comprehensive computational RNA enzymology approach aimed to provide a unified interpretation of existing structural and functional data. Simulations in the crystalline environment and in solution provide insight into the origins of observed differences in crystal structures, and coalesce on a common active site architecture, and dynamical ensemble in solution. We use GPU-accelerated free energy methods with enhanced sampling to ascertain microscopic nucleobase pK a values of the implicated general acid and base, from which predicted activity-pH profiles can be compared directly with experiments. Next, ab initio quantum mechanical/molecular mechanical (QM/MM) simulations with full dynamic solvation under periodic boundary conditions are used to determine mechanistic pathways through multi-dimensional free energy landscapes for the reaction. We then characterize the rate-controlling transition state, and make predictions about kinetic isotope effects and linear free energy relations. Computational mutagenesis is performed to explain the origin of rate effects caused by chemical modifications and make experimentally testable predictions. Finally, we provide evidence that helps to resolve conflicting issues related to the role of metal ions in catalysis. Throughout each stage, we highlight how a conserved L-platform structural motif, to- gether with a key L-anchor residue, forms the characteristic active site scaffold enabling each of the catalytic strategies to come together not only for the twister ribozyme, but the majority of the known small nucleolytic ribozyme classes.

A Transient and Flexible Cation-π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases

Metal-dependent DNA and RNA nucleases are enzymes that cleave nucleic acids with great efficiency and precision. These enzyme-mediated hydrolytic reactions are fundamental for the replication, repair, and storage of genetic information within the cell. Here, extensive classical and quantum-based free-energy molecular simulations show that a cation-π interaction is transiently formed in situ at the metal core of Bacteriophage-λ Exonuclease (Exo-λ), during catalysis. This noncovalent interaction (Lys131-Tyr154) triggers nucleophile activation for nucleotide excision. Then, our simulations also show the oscillatory dynamics and swinging of the newly formed cation-π dyad, whose conformational change may favor proton release from the cationic Lys131 to the bulk solution, thus restoring the precatalytic protonation state in Exo-λ. Altogether, we report on the novel mechanistic character of cation-π interactions for catalysis. Structural and bioinformatic analyses support that flexible orientation and transient formation of mobile cation-π interactions may represent a common catalytic strategy to promote nucleic acid hydrolysis in DNA and RNA nucleases.