A conformational switch controls hepatitis delta virus ribozyme catalysis

Structure and sequence at an RNA template 5' end influence insertion of transgenes by an R2 retrotransposon protein

R2 non-long terminal repeat retrotransposons insert site-specifically into ribosomal RNA genes (rDNA) in a broad range of multicellular eukaryotes. R2-encoded proteins can be leveraged to mediate transgene insertion at 28S rDNA loci in cultured human cells. This strategy, precise RNA-mediated insertion of transgenes (PRINT), relies on the codelivery of an mRNA encoding R2 protein and an RNA template encoding a transgene cassette of choice. Here, we demonstrate that the PRINT RNA template 5' module, which as a complementary DNA 3' end will generate the transgene 5' junction with rDNA, influences the efficiency and mechanism of gene insertion. Iterative design and testing identified optimal 5' modules consisting of a hepatitis delta virus-like ribozyme fold with high thermodynamic stability, suggesting that RNA template degradation from its 5' end may limit transgene insertion efficiency. We also demonstrate that transgene 5' junction formation can be either precise, formed by annealing the 3' end of first-strand complementary DNA with the upstream target site, or imprecise, by end-joining, but this difference in junction formation mechanism is not a major determinant of insertion efficiency. Sequence characterization of imprecise end-joining events indicates surprisingly minimal reliance on microhomology. Our findings expand the current understanding of the role of R2 retrotransposon transcript sequence and structure, and especially the 5' ribozyme fold, for retrotransposon mobility and RNA-templated gene synthesis in cells.

Assessing RNA atomistic force fields via energy landscape explorations in implicit solvent

Predicting the structure and dynamics of RNA molecules still proves challenging because of the relative scarcity of experimental RNA structures on which to train models and the very sensitive nature of RNA towards its environment. In the last decade, several atomistic force fields specifically designed for RNA have been proposed and are commonly used for simulations. However, it is not necessarily clear which force field is the most suitable for a given RNA molecule. In this contribution, we propose the use of the computational energy landscape framework to explore the energy landscape of RNA systems as it can bring complementary information to the more standard approaches of enhanced sampling simulations based on molecular dynamics. We apply the EL framework to the study of a small RNA pseudoknot, the Aquifex aeolicus tmRNA pseudoknot PK1, and we compare the results of five different RNA force fields currently available in the AMBER simulation software, in implicit solvent. With this computational approach, we can not only compare the predicted 'native' states for the different force fields, but the method enables us to study metastable states as well. As a result, our comparison not only looks at structural features of low energy folded structures, but provides insight into folding pathways and higher energy excited states, opening to the possibility of assessing the validity of force fields also based on kinetics and experiments providing information on metastable and unfolded states.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-024-01202-9.

New Insights into the Dependence of CPEB3 Ribozyme Cleavage on Mn2+ and Mg2

CPEB3 ribozyme is a self-cleaving RNA that occurs naturally in mammals and requires divalent metal ions for efficient activity. Ribozymes exhibit preferences for specific metal ions, but the exact differences in the catalytic mechanisms of various metal ions on the CPEB3 ribozyme remain unclear. Our findings reveal that Mn2+ functions as a more effective cofactor for CPEB3 ribozyme catalysis compared to Mg2+, as confirmed by its stronger binding affinity to CPEB3 by EPR. Cleavage assays of CPEB3 mutants and molecular docking analyses further showed that excessive Mn2+ ions can bind to a second binding site near the catalytic site, hindering CPEB3 catalytic efficiency and contributing to the Mn2+ bell-shaped curve. These results implicate a pivotal role for the local nucleobase-Mn2+ interactions in facilitating RNA folding and modulating the directed attack of nucleophilic reagents. Our study provides new insights and experimental evidence for exploring the divalent cation dependent cleavage mechanism of the CPEB3 ribozyme.

Identification of HDV-like theta ribozymes involved in tRNA-based recoding of gut bacteriophages

Trillions of microorganisms, collectively known as the microbiome, inhabit our bodies with the gut microbiome being of particular interest in biomedical research. Bacteriophages, the dominant virome constituents, can utilize suppressor tRNAs to switch to alternative genetic codes (e.g., the UAG stop-codon is reassigned to glutamine) while infecting hosts with the standard bacterial code. However, what triggers this switch and how the bacteriophage manipulates its host is poorly understood. Here, we report the discovery of a subgroup of minimal hepatitis delta virus (HDV)-like ribozymes - theta ribozymes - potentially involved in the code switch leading to the expression of recoded lysis and structural phage genes. We demonstrate their HDV-like self-scission behavior in vitro and find them in an unreported context often located with their cleavage site adjacent to tRNAs, indicating a role in viral tRNA maturation and/or regulation. Every fifth associated tRNA is a suppressor tRNA, further strengthening our hypothesis. The vast abundance of tRNA-associated theta ribozymes - we provide 1753 unique examples - highlights the importance of small ribozymes as an alternative to large enzymes that usually process tRNA 3'-ends. Our discovery expands the short list of biological functions of small HDV-like ribozymes and introduces a previously unknown player likely involved in the code switch of certain recoded gut bacteriophages.

Inhibition of Cpeb3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs and enhances object location memory

A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element-binding protein 3 (Cpeb3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the Cpeb3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.

Inhibition of CPEB3 ribozyme elevates CPEB3 protein expression and polyadenylation of its target mRNAs, and enhances object location memory

A self-cleaving ribozyme that maps to an intron of the cytoplasmic polyadenylation element binding protein 3 (CPEB3) gene is thought to play a role in human episodic memory, but the underlying mechanisms mediating this effect are not known. We tested the activity of the murine sequence and found that the ribozyme's self-scission half-life matches the time it takes an RNA polymerase to reach the immediate downstream exon, suggesting that the ribozyme-dependent intron cleavage is tuned to co-transcriptional splicing of the CPEB3 mRNA. Our studies also reveal that the murine ribozyme modulates maturation of its harboring mRNA in both cultured cortical neurons and the hippocampus: inhibition of the ribozyme using an antisense oligonucleotide leads to increased CPEB3 protein expression, which enhances polyadenylation and translation of localized plasticity-related target mRNAs, and subsequently strengthens hippocampal-dependent long-term memory. These findings reveal a previously unknown role for self-cleaving ribozyme activity in regulating experience-induced co-transcriptional and local translational processes required for learning and memory.

General Strategies for RNA X-ray Crystallography

An extremely small proportion of the X-ray crystal structures deposited in the Protein Data Bank are of RNA or RNA-protein complexes. This is due to three main obstacles to the successful determination of RNA structure: (1) low yields of pure, properly folded RNA; (2) difficulty creating crystal contacts due to low sequence diversity; and (3) limited methods for phasing. Various approaches have been developed to address these obstacles, such as native RNA purification, engineered crystallization modules, and incorporation of proteins to assist in phasing. In this review, we will discuss these strategies and provide examples of how they are used in practice.

Thermodynamics of unfolding mechanisms of mouse mammary tumor virus pseudoknot from a coarse-grained loop-entropy model

Pseudoknotted RNA molecules play important biological roles that depend on their folded structure. To understand the underlying principles that determine their thermodynamics and folding/unfolding mechanisms, we carried out a study on a variant of the mouse mammary tumor virus pseudoknotted RNA (VPK), a widely studied model system for RNA pseudoknots. Our method is based on a coarse-grained discrete-state model and the algorithm of PK3D (pseudoknot structure predictor in three-dimensional space), with RNA loops explicitly constructed and their conformational entropic effects incorporated. Our loop entropy calculations are validated by accurately capturing previously measured melting temperatures of RNA hairpins with varying loop lengths. For each of the hairpins that constitutes the VPK, we identified alternative conformations that are more stable than the hairpin structures at low temperatures and predicted their populations at different temperatures. Our predictions were validated by thermodynamic experiments on these hairpins. We further computed the heat capacity profiles of VPK, which are in excellent agreement with available experimental data. Notably, our model provides detailed information on the unfolding mechanisms of pseudoknotted RNA. Analysis of the distribution of base-pairing probability of VPK reveals a cooperative unfolding mechanism instead of a simple sequential unfolding of first one stem and then the other. Specifically, we find a simultaneous "loosening" of both stems as the temperature is raised, whereby both stems become partially melted and co-exist during the unfolding process.

Structure and chemistry of enzymatic active sites that play a role in the switch and conformation mechanism

Enzymes, which are biological molecules, are constructed from polypeptide chains, and these molecules are activated through reaction mechanisms. It is the role of enzymes to speed up chemical reactions that are used to build or break down cell structures. Activation energy is reduced by the enzymes' selective binding of substrates in a protected environment. In enzyme tertiary structures, the active sites are commonly situated in a "cleft," which necessitates the diffusion of substrates and products. The amino acid residues of the active site may be far apart in the primary structure owing to the folding required for tertiary structure. Due to their critical role in substrate binding and attraction, changes in amino acid structure at or near the enzyme's active site usually alter enzyme activity. At the enzyme's active site, or where the chemical reactions occur, the substrate is bound. Enzyme substrates are the primary targets of the enzyme's active site, which is designed to assist in the chemical reaction. This chapter elucidates the summary of structure and chemistry of enzymes, their active site features, charges and role of water in the structures to clarify the biochemistry of the enzymes in the depth of atomic features.

Knotify: An Efficient Parallel Platform for RNA Pseudoknot Prediction Using Syntactic Pattern Recognition

Obtaining valuable clues for noncoding RNA (ribonucleic acid) subsequences remains a significant challenge, acknowledging that most of the human genome transcribes into noncoding RNA parts related to unknown biological operations. Capturing these clues relies on accurate “base pairing” prediction, also known as “RNA secondary structure prediction”. As COVID-19 is considered a severe global threat, the single-stranded SARS-CoV-2 virus reveals the importance of establishing an efficient RNA analysis toolkit. This work aimed to contribute to that by introducing a novel system committed to predicting RNA secondary structure patterns (i.e., RNA’s pseudoknots) that leverage syntactic pattern-recognition strategies. Having focused on the pseudoknot predictions, we formalized the secondary structure prediction of the RNA to be primarily a parsing and, secondly, an optimization problem. The proposed methodology addresses the problem of predicting pseudoknots of the first order (H-type). We introduce a context-free grammar (CFG) that affords enough expression power to recognize potential pseudoknot pattern. In addition, an alternative methodology of detecting possible pseudoknots is also implemented as well, using a brute-force algorithm. Any input sequence may highlight multiple potential folding patterns requiring a strict methodology to determine the single biologically realistic one. We conscripted a novel heuristic over the widely accepted notion of free-energy minimization to tackle such ambiguity in a performant way by utilizing each pattern’s context to unveil the most prominent pseudoknot pattern. The overall process features polynomial-time complexity, while its parallel implementation enhances the end performance, as proportional to the deployed hardware. The proposed methodology does succeed in predicting the core stems of any RNA pseudoknot of the test dataset by performing a 76.4% recall ratio. The methodology achieved a F1-score equal to 0.774 and MCC equal 0.543 in discovering all the stems of an RNA sequence, outperforming the particular task. Measurements were taken using a dataset of 262 RNA sequences establishing a performance speed of 1.31, 3.45, and 7.76 compared to three well-known platforms. The implementation source code is publicly available under knotify github repo.

Self-cleaving ribozymes: substrate specificity and synthetic biology applications

Various self-cleaving ribozymes appearing in nature catalyze the sequence-specific intramolecular cleavage of RNA and can be engineered to catalyze cleavage of appropriate substrates in an intermolecular fashion, thus acting as true catalysts. The mechanisms of the small, self-cleaving ribozymes have been extensively studied and reviewed previously. Self-cleaving ribozymes can possess high catalytic activity and high substrate specificity; however, substrate specificity is also engineerable within the constraints of the ribozyme structure. While these ribozymes share a common fundamental catalytic mechanism, each ribozyme family has a unique overall architecture and active site organization, indicating that several distinct structures yield this chemical activity. The multitude of catalytic structures, combined with some flexibility in substrate specificity within each family, suggests that such catalytic RNAs, taken together, could access a wide variety of substrates. Here, we give an overview of 10 classes of self-cleaving ribozymes and capture what is understood about their substrate specificity and synthetic applications. Evolution of these ribozymes in an RNA world might be characterized by the emergence of a new ribozyme family followed by rapid adaptation or diversification for specific substrates.

Self-cleaving ribozymes have become important tools of synthetic biology. Here we summarize the substrate specificity and applications of the main classes of these ribozymes.

Graph Theoretical Methods and Workflows for Searching and Annotation of RNA Tertiary Base Motifs and Substructures

The increasing number and complexity of structures containing RNA chains in the Protein Data Bank (PDB) have led to the need for automated structure annotation methods to replace or complement expert visual curation. This is especially true when searching for tertiary base motifs and substructures. Such base arrangements and motifs have diverse roles that range from contributions to structural stability to more direct involvement in the molecule's functions, such as the sites for ligand binding and catalytic activity. We review the utility of computational approaches in annotating RNA tertiary base motifs in a dataset of PDB structures, particularly the use of graph theoretical algorithms that can search for such base motifs and annotate them or find and annotate clusters of hydrogen-bond-connected bases. We also demonstrate how such graph theoretical algorithms can be integrated into a workflow that allows for functional analysis and comparisons of base arrangements and sub-structures, such as those involved in ligand binding. The capacity to carry out such automatic curations has led to the discovery of novel motifs and can give new context to known motifs as well as enable the rapid compilation of RNA 3D motifs into a database.

Affinity and Structural Analysis of the U1A RNA Recognition Motif with Engineered Methionines to Improve Experimental Phasing

RNA plays a central role in all organisms and can fold into complex structures to orchestrate function. Visualization of such structures often requires crystallization, which can be a bottleneck in the structure-determination process. To promote crystallization, an RNA-recognition motif (RRM) of the U1A spliceosomal protein has been co-opted as a crystallization module. Specifically, the U1-snRNA hairpin II (hpII) single-stranded loop recognized by U1A can be transplanted into an RNA target to promote crystal contacts and to attain phase information via molecular replacement or anomalous diffraction methods using selenomethionine. Herein, we produced the F37M/F77M mutant of U1A to augment the phasing capability of this powerful crystallization module. Selenomethionine-substituted U1A(F37M/F77M) retains high affinity for hpII (K _D of 59.7 ± 11.4 nM). The 2.20 Å resolution crystal structure reveals that the mutated sidechains make new S-π interactions in the hydrophobic core and are useful for single-wavelength anomalous diffraction. Crystals were also attained of U1A(F37M/F77M) in complex with a bacterial preQ₁-II riboswitch. The F34M/F37M/F77M mutant was introduced similarly into a lab-evolved U1A variant (TBP6.9) that recognizes the internal bulged loop of HIV-1 TAR RNA. We envision that this short RNA sequence can be placed into non-essential duplex regions to promote crystallization and phasing of target RNAs. We show that selenomethionine-substituted TBP6.9(F34M/F37M/F77M) binds a TAR variant wherein the apical loop was replaced with a GNRA tetraloop (K _D of 69.8 ± 2.9 nM), laying the groundwork for use of TBP6.9(F34M/F37M/F77M) as a crystallization module. These new tools are available to the research community.

Synthesis of RNA-based gene regulatory devices for redirecting cellular signaling events mediated by p53

Rationale: The p53 gene is a well-known tumor suppressor, and its mutation often contributes to the occurrence and development of tumors. Due to the diversity and complexity of p53 mutations, there is still no effective p53 gene therapy. In this study, we designed and constructed an aptazyme switch that could effectively sense cellular wild-type p53 protein and regulate downstream gene function flexibly. The application of this artificial device in combination with Cre-LoxP and dCas9-VP64 tools achieved a precisely targeted killing effect on tumor cells. Methods: The affinity of the aptamer to p53 protein was verified by SPR. p53 aptazyme and gene circuits were chemically synthesized. The function of the gene circuit was detected by cell proliferation assay, apoptosis assay and Western blot. The nude mouse transplantation tumor experiment was used to evaluate the inhibitory effect of gene circuits on tumor cells in vivo. Results: The results of the SPR experiment showed that the p53 aptamer RNA sequence had a robust binding effect with p53 protein. The p53 aptazyme could efficiently sense wild-type p53 protein and initiate self-cleavage in cells. The Cre-p53 aptazyme gene circuit and dCas9-VP64/sgRNA mediated gene circuit designed based on p53 aptazyme significantly inhibited the growth and promoted the apoptosis of wild-type p53-deficient cancer cells in vitro. In addition, the gene circuits also had a significant inhibitory effect on tumors in vivo. Conclusion: The study developed a novel and efficient ribozyme switch for p53-specific recognition and provided a modular strategy for aptazyme binding to cellular proteins. In addition, the p53 aptazyme successfully inhibited tumor growth through a combined application with other synthetic biological tools, providing a new perspective for cancer therapy.

Fundamental studies of functional nucleic acids: aptamers, riboswitches, ribozymes and DNAzymes

This review aims at juxtaposing common versus distinct structural and functional strategies that are applied by aptamers, riboswitches, and ribozymes/DNAzymes. Focusing on recently discovered systems, we begin our analysis with small-molecule binding aptamers, with emphasis on in vitro-selected fluorogenic RNA aptamers and their different modes of ligand binding and fluorescence activation. Fundamental insights are much needed to advance RNA imaging probes for detection of exo- and endogenous RNA and for RNA process tracking. Secondly, we discuss the latest gene expression-regulating mRNA riboswitches that respond to the alarmone ppGpp, to PRPP, to NAD⁺, to adenosine and cytidine diphosphates, and to precursors of thiamine biosynthesis (HMP-PP), and we outline new subclasses of SAM and tetrahydrofolate-binding RNA regulators. Many riboswitches bind protein enzyme cofactors that, in principle, can catalyse a chemical reaction. For RNA, however, only one system (glmS ribozyme) has been identified in Nature thus far that utilizes a small molecule - glucosamine-6-phosphate - to participate directly in reaction catalysis (phosphodiester cleavage). We wonder why that is the case and what is to be done to reveal such likely existing cellular activities that could be more diverse than currently imagined. Thirdly, this brings us to the four latest small nucleolytic ribozymes termed twister, twister-sister, pistol, and hatchet as well as to in vitro selected DNA and RNA enzymes that promote new chemistry, mainly by exploiting their ability for RNA labelling and nucleoside modification recognition. Enormous progress in understanding the strategies of nucleic acids catalysts has been made by providing thorough structural fundaments (e.g. first structure of a DNAzyme, structures of ribozyme transition state mimics) in combination with functional assays and atomic mutagenesis.

Light-controlled twister ribozyme with single-molecule detection resolves RNA function in time and space

Small ribozymes such as Oryza sativa twister spontaneously cleave their own RNA when the ribozyme folds into its active conformation. The coupling between twister folding and self-cleavage has been difficult to study, however, because the active ribozyme rapidly converts to product. Here, we describe the synthesis of a photocaged nucleotide that releases guanosine within microseconds upon photosolvolysis with blue light. Application of this tool to O. sativa twister achieved the spatial (75 µm) and temporal (≤30 ms) control required to resolve folding and self-cleavage events when combined with single-molecule fluorescence detection of the ribozyme folding pathway. Real-time observation of single ribozymes after photo-deprotection showed that the precleaved folded state is unstable and quickly unfolds if the RNA does not react. Kinetic analysis showed that Mg²⁺ and Mn²⁺ ions increase ribozyme efficiency by making transitions to the high energy active conformation more probable, rather than by stabilizing the folded ground state or the cleaved product. This tool for light-controlled single RNA folding should offer precise and rapid control of other nucleic acid systems.

The roles of structural dynamics in the cellular functions of RNAs

RNAs fold into 3D structures that range from simple helical elements to complex tertiary structures and quaternary ribonucleoprotein assemblies. The functions of many regulatory RNAs depend on how their 3D structure changes in response to a diverse array of cellular conditions. In this Review, we examine how the structural characterization of RNA as dynamic ensembles of conformations, which form with different probabilities and at different timescales, is improving our understanding of RNA function in cells. We discuss the mechanisms of gene regulation by microRNAs, riboswitches, ribozymes, post-transcriptional RNA modifications and RNA-binding proteins, and how the cellular environment and processes such as liquid-liquid phase separation may affect RNA folding and activity. The emerging RNA-ensemble-function paradigm is changing our perspective and understanding of RNA regulation, from in vitro to in vivo and from descriptive to predictive.

VfoldLA: A web server for loop assembly-based prediction of putative 3D RNA structures

RNA three-dimensional (3D) structures are critical for RNA cellular functions. However, structure prediction for large and complex RNAs remains a challenge, which hampers our understanding of RNA structure-function relationship. We here report a new web server, the VfoldLA server (http://rna.physics.missouri.edu/vfoldLA), for the prediction of RNA 3D structures from nucleotide sequences and base-pair information (2D structure). This server is based on the recently developed VfoldLA, a model that classifies the single-stranded loops (junctions) into four different types and according to the loop-helix connections, assembles RNA 3D structures from the loop/junction templates. The VfoldLA web server provides a user-friendly online interface for a fully automated prediction of putative 3D RNA structures using VfoldLA. With a single-RNA or RNA-RNA complex sequence and 2D structure as input, the server generates structure(s) with the JSmol visualization along with a downloadable PDB file. The output result may serve as useful scaffolds for future structure refinement studies.

Classification of the nucleolytic ribozymes based upon catalytic mechanism

The nucleolytic ribozymes carry out site-specific RNA cleavage reactions by nucleophilic attack of the 2'-oxygen atom on the adjacent phosphorus with an acceleration of a million-fold or greater. A major part of this arises from concerted general acid-base catalysis. Recent identification of new ribozymes has expanded the group to a total of nine and this provides a new opportunity to identify sub-groupings according to the nature of the general base and acid. These include nucleobases, hydrated metal ions, and 2'-hydroxyl groups. Evolution has selected a number of different combinations of these elements that lead to efficient catalysis. These differences provide a new mechanistic basis for classifying these ribozymes.

Hatchet ribozyme structure and implications for cleavage mechanism

Self-cleaving ribozymes are RNAs that catalyze position-specific cleavage of their phosphodiester backbone. The cleavage site of the newly discovered hatchet ribozyme is located at the very 5′ end of its consensus secondary structure motif. Here we report on the 2.1-Å crystal structure of the hatchet ribozyme in the product state, which defines its intricate tertiary fold and identifies key residues lining the catalytic pocket. This in turn has allowed us to propose a model of the precatalytic state structure and a role in catalysis for a conserved guanine. This study therefore provides a structure-based platform toward an improved understanding of the catalytic mechanism of hatchet ribozymes.

Small self-cleaving ribozymes catalyze site-specific cleavage of their own phosphodiester backbone with implications for viral genome replication, pre-mRNA processing, and alternative splicing. We report on the 2.1-Å crystal structure of the hatchet ribozyme product, which adopts a compact pseudosymmetric dimeric scaffold, with each monomer stabilized by long-range interactions involving highly conserved nucleotides brought into close proximity of the scissile phosphate. Strikingly, the catalytic pocket contains a cavity capable of accommodating both the modeled scissile phosphate and its flanking 5′ nucleoside. The resulting modeled precatalytic conformation incorporates a splayed-apart alignment at the scissile phosphate, thereby providing structure-based insights into the in-line cleavage mechanism. We identify a guanine lining the catalytic pocket positioned to contribute to cleavage chemistry. The functional relevance of structure-based insights into hatchet ribozyme catalysis is strongly supported by cleavage assays monitoring the impact of selected nucleobase and atom-specific mutations on ribozyme activity.

An algebraic language for RNA pseudoknots comparison

BACKGROUND

RNA secondary structure comparison is a fundamental task for several studies, among which are RNA structure prediction and evolution. The comparison can currently be done efficiently only for pseudoknot-free structures due to their inherent tree representation.

RESULTS

In this work, we introduce an algebraic language to represent RNA secondary structures with arbitrary pseudoknots. Each structure is associated with a unique algebraic RNA tree that is derived from a tree grammar having concatenation, nesting and crossing as operators. From an algebraic RNA tree, an abstraction is defined in which the primary structure is neglected. The resulting structural RNA tree allows us to define a new measure of similarity calculated exploiting classical tree alignment.

CONCLUSIONS

The tree grammar with its operators permit to uniquely represent any RNA secondary structure as a tree. Structural RNA trees allow us to perform comparison of RNA secondary structures with arbitrary pseudoknots without taking into account the primary structure.

Hepatitis Delta Antigen Regulates mRNA and Antigenome RNA Levels during Hepatitis Delta Virus Replication

Hepatitis delta virus (HDV) is a satellite of hepatitis B virus that increases the severity of acute and chronic liver disease. HDV produces three processed RNAs that accumulate in infected cells: the circular genome; the circular antigenome, which serves as a replication intermediate; and lesser amounts of the mRNA, which encodes the sole viral protein, hepatitis delta antigen (HDAg). The HDV genome and antigenome RNAs form ribonucleoprotein complexes with HDAg. Although HDAg is required for HDV replication, it is not known how the relative amounts of HDAg and HDV RNA affect replication, or whether HDAg synthesis is regulated by the virus. Using a novel transfection system in which HDV replication is initiated using in vitro-synthesized circular HDV RNAs, HDV replication was found to depend strongly on the relative amounts of HDV RNA and HDAg. HDV controls these relative amounts via differential effects of HDAg on the production of HDV mRNA and antigenome RNA, both of which are synthesized from the genome RNA template. mRNA synthesis is favored at low HDAg levels but becomes saturated at high HDAg concentrations. Antigenome RNA accumulation increases linearly with HDAg and dominates at high HDAg levels. These results provide a conceptual model for how HDV antigenome RNA production and mRNA transcription are controlled from the earliest stage of infection onward and also demonstrate that, in this control, HDV behaves similarly to other negative-strand RNA viruses, even though there is no genetic similarity between them.IMPORTANCE Hepatitis delta virus (HDV) is a satellite of hepatitis B virus that increases the severity of liver disease; approximately 15 million people are chronically infected worldwide. There are no licensed therapies available. HDV is not related to any known virus, and few details regarding its replication cycle are known. One key question is whether and how HDV regulates the relative amounts of viral RNA and protein in infected cells. Such regulation might be important because the HDV RNA and protein form complexes that are essential for HDV replication, and the proper stoichiometry of these complexes could be critical for their function. Our results show that the relative amounts of HDV RNA and protein in cells are indeed important for HDV replication and that the virus does control them. These observations indicate that further study of these regulatory mechanisms is required to better understand replication of this serious human pathogen.

Evidence That Nucleophile Deprotonation Exceeds Bond Formation in the HDV Ribozyme Transition State

Steric constraints imposed by the active sites of protein and RNA enzymes pose major challenges to the investigation of structure-function relationships within these systems. As a strategy to circumvent such constraints in the HDV ribozyme, we have synthesized phosphoramidites from propanediol derivatives and incorporated them at the 5'-termini of RNA and DNA oligonucleotides to generate a series of novel substrates with nucleophiles perturbed electronically through geminal fluorination. In nonenzymatic, hydroxide-catalyzed intramolecular transphosphorylation of the DNA substrates, pH-rate profiles revealed that fluorine substitution reduces the maximal rate and the kinetic p K_a, consistent with the expected electron-withdrawing effect. In HDV ribozyme reactions, we observed that the RNA substrates undergo transphosphorylation relatively efficiently, suggesting that the conformational constraints imposed by a ribofuranose ring are not strictly required for ribozyme catalysis. In contrast to the nonenzymatic reactions, however, substrate fluorination modestly increases the ribozyme reaction rate, consistent with a mechanism in which (1) the 2'-hydroxyl nucleophile exists predominantly in its neutral, protonated form in the ground state and (2) the 2'-hydroxyl bears some negative charge in the rate-determining step, consistent with a transition state in which the extent of 2'-OH deprotonation exceeds the extent of P-O bond formation.

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.

Nucleic acids: function and potential for abiogenesis

The emergence of functional cooperation between the three main classes of biomolecules - nucleic acids, peptides and lipids - defines life at the molecular level. However, how such mutually interdependent molecular systems emerged from prebiotic chemistry remains a mystery. A key hypothesis, formulated by Crick, Orgel and Woese over 40 year ago, posits that early life must have been simpler. Specifically, it proposed that an early primordial biology lacked proteins and DNA but instead relied on RNA as the key biopolymer responsible not just for genetic information storage and propagation, but also for catalysis, i.e. metabolism. Indeed, there is compelling evidence for such an 'RNA world', notably in the structure of the ribosome as a likely molecular fossil from that time. Nevertheless, one might justifiably ask whether RNA alone would be up to the task. From a purely chemical perspective, RNA is a molecule of rather uniform composition with all four bases comprising organic heterocycles of similar size and comparable polarity and pK a values. Thus, RNA molecules cover a much narrower range of steric, electronic and physicochemical properties than, e.g. the 20 amino acid side-chains of proteins. Herein we will examine the functional potential of RNA (and other nucleic acids) with respect to self-replication, catalysis and assembly into simple protocellular entities.

Kinetic Parameters of trans Scission by Extended HDV-like Ribozymes and the Prospect for the Discovery of Genomic trans-Cleaving RNAs

Hepatitis delta virus (HDV)-like ribozymes are self-cleaving catalytic RNAs with a widespread distribution in nature and biological roles ranging from self-scission during rolling-circle replication in viroids to co-transcriptional processing of eukaryotic retrotransposons, among others. The ribozymes fold into a double pseudoknot with a common catalytic core motif and highly variable peripheral domains. Like other self-cleaving ribozymes, HDV-like ribozymes can be converted into trans-acting catalytic RNAs by bisecting the self-cleaving variants at non-essential loops. Here we explore the trans-cleaving activity of ribozymes derived from the largest examples of the ribozymes (drz-Agam-2 family), which contain an extended domain between the substrate strand and the rest of the RNA. When this peripheral domain is bisected at its distal end, the substrate strand is recognized through two helices, rather than just one 7 bp helix common among the HDV ribozymes, resulting in stronger binding and increased sequence specificity. Kinetic characterization of the extended trans-cleaving ribozyme revealed an efficient trans-cleaving system with a surprisingly high K_M', supporting a model that includes a recently proposed activation barrier related to the assembly of the catalytically competent ribozyme. The ribozymes also exhibit a very long k_off for the products (∼2 weeks), resulting in a trade-off between sequence specificity and turnover. Finally, structure-based searches for the catalytic cores of these ribozymes in the genome of the mosquito Anopheles gambiae, combined with sequence searches for their putative substrates, revealed two potential ribozyme-substrate pairs that may represent examples of natural trans-cleaving ribozymes.

RNA Degradation in Neurodegenerative Disease

Ribonucleic acid (RNA) homeostasis is dynamically modulated in response to changing physiological conditions. Tight regulation of RNA abundance through both transcription and degradation determines the amount, timing, and location of protein translation. This balance is of particular importance in neurons, which are among the most metabolically active and morphologically complex cells in the body. As a result, any disruptions in RNA degradation can have dramatic consequences for neuronal health. In this chapter, we will first discuss mechanisms of RNA stabilization and decay. We will then explore how the disruption of these pathways can lead to neurodegenerative disease.

Crystal structure of an RNA-cleaving DNAzyme

In addition to storage of genetic information, DNA can also catalyze various reactions. RNA-cleaving DNAzymes are the catalytic DNAs discovered the earliest, and they can cleave RNAs in a sequence-specific manner. Owing to their great potential in medical therapeutics, virus control, and gene silencing for disease treatments, RNA-cleaving DNAzymes have been extensively studied; however, the mechanistic understandings of their substrate recognition and catalysis remain elusive. Here, we report three catalytic form 8-17 DNAzyme crystal structures. 8-17 DNAzyme adopts a V-shape fold, and the Pb²⁺ cofactor is bound at the pre-organized pocket. The structures with Pb²⁺ and the modification at the cleavage site captured the pre-catalytic state of the RNA cleavage reaction, illustrating the unexpected Pb²⁺-accelerated catalysis, intrinsic tertiary interactions, and molecular kink at the active site. Our studies reveal that DNA is capable of forming a compacted structure and that the functionality-limited bio-polymer can have a novel solution for a functional need in catalysis.

Structure-based mechanistic insights into catalysis by small self-cleaving ribozymes

Small self-cleaving ribozymes are widely distributed in nature and are essential for rolling-circle-based replication of satellite and pathogenic RNAs. Earlier structure-function studies on the hammerhead, hairpin, glmS, hepatitis delta virus and Varkud satellite ribozymes have provided insights into their overall architecture, their catalytic active site organization, and the role of nearby nucleobases and hydrated divalent cations in facilitating general acid-base and electrostatic-mediated catalysis. This review focuses on recent structure-function research on active site alignments and catalytic mechanisms of the Rzb hammerhead ribozyme, as well as newly-identified pistol, twister and twister-sister ribozymes. In contrast to an agreed upon mechanistic understanding of self-cleavage by Rzb hammerhead and pistol ribozymes, there exists a divergence of views as to the cleavage site alignments and catalytic mechanisms adopted by twister and twister-sister ribozymes. One approach to resolving this conundrum would be to extend the structural studies from currently available pre-catalytic conformations to their transition state mimic vanadate counterparts for both ribozymes.

Allosteric Modulation of the Faecalibacterium prausnitzii Hepatitis Delta Virus-like Ribozyme by Glucosamine 6-Phosphate: The Substrate of the Adjacent Gene Product

Self-cleaving ribozymes were discovered 30 years ago and have been found throughout nature, from bacteria to animals, but little is known about their biological functions and regulation, particularly how cofactors and metabolites alter their activity. A hepatitis delta virus-like self-cleaving ribozyme maps upstream of a phosphoglucosamine mutase (glmM) open reading frame in the genome of the human gut bacterium Faecalibacterium prausnitzii. The presence of a ribozyme in the untranslated region of glmM suggests a regulation mechanism of gene expression. In the bacterial hexosamine biosynthesis pathway, the enzyme glmM catalyzes the isomerization of glucosamine 6-phosphate into glucosamine 1-phosphate. In this study, we investigated the effect of these metabolites on the co-transcriptional self-cleavage rate of the ribozyme. Our results suggest that glucosamine 6-phosphate, but not glucosamine 1-phosphate, is an allosteric ligand that increases the self-cleavage rate of drz-Fpra-1, providing the first known example of allosteric modulation of a self-cleaving ribozyme by the substrate of the adjacent gene product. Given that the ribozyme is activated by the glmM substrate, but not the product, this allosteric modulation may represent a potential feed-forward mechanism of gene expression regulation in bacteria.

Construction and characterization of an improved DNA-launched infectious clone of duck hepatitis a virus type 1

Background

DNA-launched infectious system is a useful tool with high rescue efficiency that allows the introduction of mutations in specific positions to investigate the function of an individual viral element. Rescued virus particles could be harvested by directly transfecting the DNA-launched recombinant plasmid to the host cells, which will reduce labor and experimental cost by skipping the in vitro transcription assay.

Methods

A total of four overlapping fragments covering the entire viral genome were amplified and then were assembled into a transformation vector based on pIRES2-EGFP to establish the DNA-launched infectious system of duck hepatitis A virus type 1 (DHAV-1), named pIR-DHAV-1. Reverse transcription polymerase chain reaction (RT-PCR) detection, quantitative real-time polymerase chain reaction (qRT-PCR), western blotting assay and indirect immunofluorescence (IFA) were conducted for rescued virus identification. A total of 4.0 μg of recombinant plasmid of pIR-DHAV-1 and in vitro transcribed product of 4.0 μg of RNA-launched infectious clone named pR-DHAV-1 were transfected into BHK-21 cells to analyze the rescue efficiency. Following that, tissue tropism of rescued virus (rDHAV-1) and parental virus (pDHAV-1) were assayed for virulence testing in 1-day-old ducklings.

Results

Rescued virus particles carry the designed genetic marker which could be harvested by directly transfecting pIR-DHAV-1 to BHK-21 cells. The qRT-PCR and western blotting results indicated that rDHAV-1 shared similar growth characteristics with pDHAV-1. Furthermore, DNA-launched infectious system possessed much higher rescue efficiency assay compared to RNA-launched infectious system. The mutation at position 3042 from T to C has no impact on viral replication and tissue tropism. From 1 h post infection (hpi) to 48 hpi, the viral RNA copies of rDHAV-1 in liver were the highest among the six tested tissues (with an exception of thymus at 6 hpi), while the viral RNA copy numbers in heart and kidney were alternately the lowest.

Conclusion

We have constructed a genetically stable and highly pathogenic DNA-launched infectious clone, from which the rescued virus could be harvested by direct transfection with recombinant plasmids. rDHAV-1 shared similar growth characteristics and tissue tropism with pDHAV-1. The DNA-launched infectious system of DHAV-1 possessed higher rescue efficiency compared to the traditional RNA-launched infectious system.

How RNA acts as a nuclease: some mechanistic comparisons in the nucleolytic ribozymes

Recent structural and mechanistic studies have shed considerable light on the catalytic mechanisms of nucleolytic ribozymes. The discovery of several new ribozymes in this class has now allowed comparisons to be made, and the beginnings of mechanistic groupings to emerge.

The structure of a nucleolytic ribozyme that employs a catalytic metal ion

The TS ribozyme (originally called "twister sister") is a catalytic RNA. We present a crystal structure of the ribozyme in a pre-reactive conformation. Two co-axial helical stacks are organized by a three-way junction and two tertiary contacts. Five divalent metal ions are directly coordinated to RNA ligands, making important contributions to the RNA architecture. The scissile phosphate lies in a quasihelical loop region that is organized by a network of hydrogen bonding. A divalent metal ion is directly bound to the nucleobase 5' to the scissile phosphate, with an inner-sphere water molecule positioned to interact with the O2' nucleophile. The rate of ribozyme cleavage correlated in a log-linear manner with divalent metal ion pKa, consistent with proton transfer in the transition state, and we propose that the bound metal ion is a likely general base for the cleavage reaction. Our data indicate that the TS ribozyme functions predominantly as a metalloenzyme.

Structural and Biochemical Properties of Novel Self-Cleaving Ribozymes

Fourteen well-defined ribozyme classes have been identified to date, among which nine are site-specific self-cleaving ribozymes. Very recently, small self-cleaving ribozymes have attracted renewed interest in their structure, biochemistry, and biological function since the discovery, during the last three years, of four novel ribozymes, termed twister, twister sister, pistol, and hatchet. In this review, we mainly address the structure, biochemistry, and catalytic mechanism of the novel ribozymes. They are characterized by distinct active site architectures and divergent, but similar, biochemical properties. The cleavage activities of the ribozymes are highly dependent upon divalent cations, pH, and base-specific mutations, which can cause changes in the nucleotide arrangement and/or electrostatic potential around the cleavage site. It is most likely that a guanine and adenine in close proximity of the cleavage site are involved in general acid-base catalysis. In addition, metal ions appear to play a structural rather than catalytic role although some of their crystal structures have shown a direct metal ion coordination to a non-bridging phosphate oxygen at the cleavage site. Collectively, the structural and biochemical data of the four newest ribozymes could contribute to advance our mechanistic understanding of how self-cleaving ribozymes accomplish their efficient site-specific RNA cleavages.

Characterization of Highly Efficient RNA-Cleaving DNAzymes that Function at Acidic pH with No Divalent Metal-Ion Cofactors

Here, we describe the characterization of new RNA‐cleaving DNAzymes that showed the highest catalytic efficiency at pH 4.0 to 4.5, and were completely inactive at pH values higher than 5.0. Importantly, these DNAzymes did not require any divalent metal ion cofactors for catalysis. This clearly suggests that protonated nucleic bases are involved in the folding of the DNAzymes into catalytically active structures and/or in the cleavage mechanism. The trans‐acting DNAzyme variants were also catalytically active. Mutational analysis revealed a conservative character of the DNAzyme catalytic core that underpins the high structural requirements of the cleavage mechanism. A significant advantage of the described DNAzymes is that they are inactive at pH values close to physiological pH and under a wide range of conditions in the presence of monovalent and divalent metal ions. These pH‐dependent DNAzymes could be used as molecular cassettes in biotechnology or nanotechnology, in molecular processes that consist of several steps. The results expand the repertoire of DNAzymes that are active under nonphysiological conditions and shed new light on the possible mechanisms of catalysis.

Unwinding the twister ribozyme: from structure to mechanism

The twister ribozyme motif has been identified by bioinformatic means very recently. Currently, four crystal structures with ordered active sites together with a series of chemical and biochemical data provide insights into how this RNA accomplishes its efficient self‐cleavage. Of particular interest for a mechanistic proposal are structural distinctions observed in the active sites that concern the conformation of the U‐A cleavage site dinucleotide (in‐line alignment of the attacking 2′‐O nucleophile to the to‐be‐cleaved P—O5′ bond versus suboptimal alignments) as well as the presence/absence of Mg²⁺ ions at the scissile phosphate. All structures support the notion that an active site guanine and the conserved adenine at the cleavage site are important contributors to cleavage chemistry, likely being involved in general acid base catalysis. Evidence for innersphere coordination of a Mg²⁺ ion to the pro‐S nonbridging oxygen of the scissile phosphate stems from two of the four crystal structures. Together with the finding of thio/rescue effects for phosphorothioate substrates, this suggests the participation of divalent ions in the overall catalytic strategy employed by twister ribozymes. In this context, it is notable that twister retains wild‐type activity when the phylogenetically conserved stem P1 is deleted, able to cleave a single nucleotide only. WIREs RNA 2017, 8:e1402. doi: 10.1002/wrna.1402

For further resources related to this article, please visit the WIREs website.

Pistol ribozyme adopts a pseudoknot fold facilitating site-specific in-line cleavage

The field of small self-cleaving nucleolytic ribozymes has been invigorated by the recent discovery of the twister, twister-sister, pistol and hatchet ribozymes. We report the crystal structure of a pistol ribozyme termed env25, which adopts a compact tertiary architecture stabilized by an embedded pseudoknot fold. The G-U cleavage site adopts a splayed-apart conformation with in-line alignment of the modeled 2'-O of G for attack on the adjacent to-be-cleaved P-O5' bond. Highly conserved residues G40 (N1 position) and A32 (N3 and 2'-OH positions) are aligned to act as a general base and a general acid, respectively, to accelerate cleavage chemistry, with their roles confirmed by cleavage assays on variants, and an increased pKa of 4.7 for A32. Our structure of the pistol ribozyme defined how the overall and local topologies dictate the in-line alignment at the G-U cleavage site, with cleavage assays on variants revealing key residues that participate in acid-base-catalyzed cleavage chemistry.

Topological constraints of structural elements in regulation of catalytic activity in HDV-like self-cleaving ribozymes

Self-cleaving ribozymes fold into intricate structures, which orient active site groups into catalytically competent conformations. Most ribozyme families have distinct catalytic cores stabilized by tertiary interactions between domains peripheral to those cores. We show that large hepatitis delta virus (HDV)-like ribozymes are activated by peripheral domains that bring two helical segments, P1 and P2, into proximity – a “pinch” that results in rate acceleration by almost three orders of magnitude. Kinetic analysis of ribozymes with systematically altered length and stability of the peripheral domain revealed that about one third of its free energy of formation is used to lower an activation energy barrier, likely related to a rate-limiting conformational change leading to the pre-catalytic state. These findings provide a quantitative view of enzyme regulation by peripheral domains and may shed light on the energetics of allosteric regulation.

Secondary structure confirmation and localization of Mg2+ ions in the mammalian CPEB3 ribozyme

Most of today's knowledge of the CPEB3 ribozyme, one of the few small self-cleaving ribozymes known to occur in humans, is based on comparative studies with the hepatitis delta virus (HDV) ribozyme, which is highly similar in cleavage mechanism and probably also in structure. Here we present detailed NMR studies of the CPEB3 ribozyme in order to verify the formation of the predicted nested double pseudoknot in solution. In particular, the influence of Mg(2+), the ribozyme's crucial cofactor, on the CPEB3 structure is investigated. NMR titrations, Tb(3+)-induced cleavage, as well as stoichiometry determination by hydroxyquinoline sulfonic acid fluorescence and equilibrium dialysis, are used to evaluate the number, location, and binding mode of Mg(2+)ions. Up to eight Mg(2+)ions interact site-specifically with the ribozyme, four of which are bound with high affinity. The global fold of the CPEB3 ribozyme, encompassing 80%-90% of the predicted base pairs, is formed in the presence of monovalent ions alone. Low millimolar concentrations of Mg(2+)promote a more compact fold and lead to the formation of additional structures in the core of the ribozyme, which contains the inner small pseudoknot and the active site. Several Mg(2+)binding sites, which are important for the functional fold, appear to be located in corresponding locations in the HDV and CPEB3 ribozyme, demonstrating the particular relevance of Mg(2+)for the nested double pseudoknot structure.

A Two-Metal-Ion-Mediated Conformational Switching Pathway for HDV Ribozyme Activation

RNA enzymes serve as a potentially powerful platform from which to design catalysts and engineer new biotechnology. A fundamental understanding of these systems provides insight to guide design. The hepatitis delta virus ribozyme (HDVr) is a small, self-cleaving RNA motif widely distributed in nature, that has served as a paradigm for understanding basic principles of RNA catalysis. Nevertheless, questions remain regarding the precise roles of divalent metal ions and key nucleotides in catalysis. In an effort to establish a reaction mechanism model consistent with available experimental data, we utilize molecular dynamics simulations to explore different conformations and metal ion binding modes along the HDVr reaction path. Building upon recent crystallographic data, our results provide a dynamic model of the HDVr reaction mechanism involving a conformational switch between multiple non-canonical G25:U20 base pair conformations in the active site. These local nucleobase dynamics play an important role in catalysis by modulating the metal binding environments of two Mg²⁺ ions that support catalysis at different steps of the reaction pathway. The first ion plays a structural role by inducing a base pair flip necessary to obtain the catalytic fold in which C75 moves towards to the scissile phosphate in the active site. Ejection of this ion then permits a second ion to bind elsewhere in the active site and facilitate nucleophile activation. The simulations collectively describe a mechanistic scenario that is consistent with currently available experimental data from crystallography, phosphorothioate substitutions, and chemical probing studies. Avenues for further experimental verification are suggested.

The Complementarity of the Loop to the Stem in DNA Pseudoknots Gives Rise to Local TAT Base-Triplets

Pseudoknots belong to an RNA structural motif that has significant roles in the biological function of RNA. An example is ribosomal frameshifting; in this mechanism, the formation of a local triplex changes the reading frame that allows for differences in the translation of mRNAs. In this work, we have used a combination of temperature-dependent UV spectroscopy and differential scanning calorimetry (DSC) to determine the unfolding thermodynamics of a set of DNA pseudoknots with the following sequence: d(TCTCTTnAAAAAAAAGAGAT5TTTTTTT), where "Tn" is a thymine loop with n=5 (PsK-5), 7 (PsK-7), 9 (PsK-9), or 11 (PsK-11). All four oligonucleotides form intramolecular pseudoknots, and the increase in the length of this loop yielded more stable pseudoknots due to higher transition temperatures and higher unfolding enthalpies. This indicates formation of one and three TAT/TAT stacks in PsK-9 and PsK-11, respectively. We have flipped one AT for a TA base pair in the core stem of these pseudoknots, preventing in this way the formation of these base-triplet stacks. The DSC curves of these pseudoknots yielded lower unfolding enthalpies, confirming the formation of a local triplex in PsK-9 and PsK-11. Furthermore, we have investigated the reaction of PsK-5 and PsK-9 with their partially complementary strands: directly by isothermal titration calorimetry and indirectly by creating a Hess cycle with the DSC data. Relative to the PsK-5 reaction, PsK-9 reacts with its complementary strand with less favorable free energy and enthalpy contributions; this indicates PsK-9 is more stable and more compact due to the formation of a local triplex.

Nuclear Magnetic Resonance-Assisted Prediction of Secondary Structure for RNA: Incorporation of Direction-Dependent Chemical Shift Constraints

Knowledge of RNA structure is necessary to determine structure–function relationships and to facilitate design of potential therapeutics. RNA secondary structure prediction can be improved by applying constraints from nuclear magnetic resonance (NMR) experiments to a dynamic programming algorithm. Imino proton walks from NOESY spectra reveal double-stranded regions. Chemical shifts of protons in GH1, UH3, and UH5 of GU pairs, UH3, UH5, and AH2 of AU pairs, and GH1 of GC pairs were analyzed to identify constraints for the 5′ to 3′ directionality of base pairs in helices. The 5′ to 3′ directionality constraints were incorporated into an NMR-assisted prediction of secondary structure (NAPSS-CS) program. When it was tested on 18 structures, including nine pseudoknots, the sensitivity and positive predictive value were improved relative to those of three unrestrained programs. The prediction accuracy for the pseudoknots improved the most. The program also facilitates assignment of chemical shifts to individual nucleotides, a necessary step for determining three-dimensional structure.

Wobble pairs of the HDV ribozyme play specific roles in stabilization of active site dynamics

The hepatitis delta virus (HDV) is the only known human pathogen whose genome contains a catalytic RNA motif (ribozyme). The overall architecture of the HDV ribozyme is that of a double-nested pseudoknot, with two GU pairs flanking the active site. Although extensive studies have shown that mutation of either wobble results in decreased catalytic activity, little work has focused on linking these mutations to specific structural effects on catalytic fitness. Here we use molecular dynamics simulations based on an activated structure to probe the active site dynamics as a result of wobble pair mutations. In both wild-type and mutant ribozymes, the in-line fitness of the active site (as a measure of catalytic proficiency) strongly depends on the presence of a C75(N3H3+)N1(O5') hydrogen bond, which positions C75 as the general acid for the reaction. Our mutational analyses show that each GU wobble supports catalytically fit conformations in distinct ways; the reverse G25U20 wobble promotes high in-line fitness, high occupancy of the C75(N3H3+)G1(O5') general-acid hydrogen bond and stabilization of the G1U37 wobble, while the G1U37 wobble acts more locally by stabilizing high in-line fitness and the C75(N3H3+)G1(O5') hydrogen bond. We also find that stable type I A-minor and P1.1 hydrogen bonding above and below the active site, respectively, prevent local structural disorder from spreading and disrupting global conformation. Taken together, our results define specific, often redundant architectural roles for several structural motifs of the HDV ribozyme active site, expanding the known roles of these motifs within all HDV-like ribozymes and other structured RNAs.

A Mini-Twister Variant and Impact of Residues/Cations on the Phosphodiester Cleavage of this Ribozyme Class

Nucleolytic ribozymes catalyze site-specific cleavage of their phosphodiester backbones. A minimal version of the twister ribozyme is reported that lacks the phylogenetically conserved stem P1 while retaining wild-type activity. Atomic mutagenesis revealed that nitrogen atoms N1 and N3 of the adenine-6 at the cleavage site are indispensable for cleavage. By NMR spectroscopy, a pKa value of 5.1 was determined for a (13) C2-labeled adenine at this position in the twister ribozyme, which is significantly shifted compared to the pKa of the same adenine in the substrate alone. This finding pinpoints at a potential role for adenine-6 in the catalytic mechanism besides the previously identified invariant guanine-48 and a Mg(2+) ion, both of which are directly coordinated to the non-bridging oxygen atoms of the scissile phosphate; for the latter, additional evidence stems from the observation that Mn(2+) or Cd(2+) accelerated cleavage of phosphorothioate substrates. The relevance of this metal ion binding site is further emphasized by a new 2.6 Å X-ray structure of a 2'-OCH3 -U5 modified twister ribozyme.

Chemistry and Biology of Self-Cleaving Ribozymes

Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.