Wild-type is the optimal sequence of the HDV ribozyme under cotranscriptional conditions

A comparative survey of the influence of small self-cleaving ribozymes on gene expression in human cell culture

Self-cleaving ribozymes are versatile tools for synthetic biologists when it comes to controlling gene expression. Up to date, 12 different classes are known, and over the past decades more and more details about their structure, cleavage mechanisms and natural environments have been uncovered. However, when these motifs are applied to mammalian gene expression constructs, the outcome can often be unexpected. A variety of factors, such as surrounding sequences and positioning of the ribozyme influences the activity and hence performance of catalytic RNAs. While some information about the efficiency of individual ribozymes (each tested in specific contexts) is known, general trends obtained from standardized, comparable experiments are lacking, complicating decisions such as which ribozyme to choose and where to insert it into the target mRNA. In many cases, application-specific optimization is required, which can be very laborious. Here, we systematically compared different classes of ribozymes within the 3'-UTR of a given reporter gene. We then examined position-dependent effects of the best-performing ribozymes. Moreover, we tested additional variants of already widely used hammerhead ribozymes originating from various organisms. We were able to identify functional structures suited for aptazyme design and generated highly efficient hammerhead ribozyme variants originating from the human genome. The present dataset will aide decisions about how to apply ribozymes for affecting gene expression as well as for developing ribozyme-based switches for controlling gene expression in human cells.

RNA sequence to structure analysis from comprehensive pairwise mutagenesis of multiple self-cleaving ribozymes

Self-cleaving ribozymes are RNA molecules that catalyze the cleavage of their own phosphodiester backbones. These ribozymes are found in all domains of life and are also a tool for biotechnical and synthetic biology applications. Self-cleaving ribozymes are also an important model of sequence-to-function relationships for RNA because their small size simplifies synthesis of genetic variants and self-cleaving activity is an accessible readout of the functional consequence of the mutation. Here, we used a high-throughput experimental approach to determine the relative activity for every possible single and double mutant of five self-cleaving ribozymes. From this data, we comprehensively identified non-additive effects between pairs of mutations (epistasis) for all five ribozymes. We analyzed how changes in activity and trends in epistasis map to the ribozyme structures. The variety of structures studied provided opportunities to observe several examples of common structural elements, and the data was collected under identical experimental conditions to enable direct comparison. Heatmap-based visualization of the data revealed patterns indicating structural features of the ribozymes including paired regions, unpaired loops, non-canonical structures, and tertiary structural contacts. The data also revealed signatures of functionally critical nucleotides involved in catalysis. The results demonstrate that the data sets provide structural information similar to chemical or enzymatic probing experiments, but with additional quantitative functional information. The large-scale data sets can be used for models predicting structure and function and for efforts to engineer self-cleaving ribozymes.

Concentrating all helper protein functions on a single entity allows rescue of recombinant measles virus by transfection of just two plasmids

The generation of recombinant measles virus (MeV) from manipulated genomes on plasmid DNA is quite a complex and inefficient process. As a member of the order Mononegavirales its single-stranded ssRNA genome in negative sense orientation is not infectious, but requires co-availability of the viral RNA-dependent RNA polymerase L, the polymerase co-factor phosphoprotein P, and the nucleocapsid protein N in defined relative amounts to establish infectious centres in transfected cell cultures that release replication-competent recombinant MeV particles. For this so-called rescue, different rescue systems were developed that rely on at least four different components. In this work, we establish a functional MeV rescue system just being composed of two components: the plasmid encoding the (modified) viral genome, and a one-helper-plasmid bundling all helper functions. In contrast to a rescue-system for Newcastle Disease Virus, another paramyxovirus, co-expression of all helper proteins by the same promoter failed. Instead, adaptation of the strength of the respective promoters to drive each helper gene´s expression to the relative expression found in MeV-infected cells or other rescue systems, which indeed adjusted respective mRNA and protein expression, yielded success, albeit not yet to the same efficacy as the four-component system. Thereby, our study paves the way for the development of easier and, after further optimization, more efficient rescue systems to generate recombinant MeV for e.g. the application as a vaccine platform or oncolytic virus, for example.

Upstream Flanking Sequence Assists Folding of an RNA Thermometer

Many heat shock genes in bacteria are regulated through a class of temperature-sensitive stem-loop (SL) RNAs called RNA thermometers (RNATs). One of the most widely studied RNATs is the Repression Of heat Shock Expression (ROSE) element associated with expression of heat shock proteins. Located in the 5'UTR, the RNAT contains one to three auxiliary hairpins upstream of it. Herein, we address roles of these upstream SLs in the folding and function of an RNAT. Bradyrhizobium japonicum is a nitrogen-fixing bacterium that experiences a wide range of temperatures in the soil and contains ROSE elements, each having multiple upstream SLs. The 5'UTR of the messenger (mRNA) for heat shock protein A (hspA) in B. japonicum has an intricate secondary structure containing three SLs upstream of the RNAT SL. While structure-function studies of the hspA RNAT itself have been reported, it has been unclear if these auxiliary SLs contribute to the temperature-sensing function of the ROSE elements. Herein, we show that the full length (FL) sequence has several melting transitions indicating that the ROSE element unfolds in a non-two-state manner. The upstream SLs are more stable than the RNAT itself, and a variant with disrupted base pairing in the SL immediately upstream of the RNAT has little influence on the melting of the RNAT. On the basis of these results and modeling of the co-transcriptional folding of the ROSE element, we propose that the upstream SLs function to stabilize the transcript and aid proper folding and dynamics of the RNAT.

Inspecting the Ribozyme Region of Hepatitis Delta Virus Genotype 1: Conservation and Variability

The hepatitis delta virus (HDV) genome has an autocatalytic region called the ribozyme, which is essential for viral replication. The aim of this study was to use next-generation sequencing (NGS) to analyze the ribozyme quasispecies (QS) in order to study its evolution and identify highly conserved regions potentially suitable for a gene-silencing strategy. HDV RNA was extracted from 2 longitudinal samples of chronic HDV patients and the ribozyme (nucleotide, nt 688-771) was analyzed using NGS. QS conservation, variability and genetic distance were analyzed. Mutations were identified by aligning sequences with their specific genotype consensus. The main relevant mutations were tested in vitro. The ribozyme was conserved overall, with a hyper-conserved region between nt 715-745. No difference in QS was observed over time. The most variable region was between nt 739-769. Thirteen mutations were observed, with three showing a higher frequency: T23C, T69C and C64 deletion. This last strongly reduced HDV replication by more than 1 log in vitro. HDV Ribozyme QS was generally highly conserved and was maintained during follow-up. The most conserved portion may be a valuable target for a gene-silencing strategy. The presence of the C64 deletion may strongly impair viral replication, as it is a potential mechanism of viral persistence.

A scenario for the emergence of protoviroids in the RNA world and for their further evolution into viroids and viroid-like RNAs by modular recombinations and mutations

Viroids are tiny, circular, and noncoding RNAs that are able to replicate and systemically infect plants. The smallest known pathogens, viroids have been proposed to represent survivors from the RNA world that likely preceded the cellular world currently dominating life on the earth. Although the small, circular, and compact nature of viroid genomes, some of which are also endowed with catalytic activity mediated by hammerhead ribozymes, support this proposal, the lack of feasible evolutionary routes and the identification of hammerhead ribozymes in a large number of DNA genomes of organisms along the tree of life have led some to question such a proposal. Here, we reassess the origin and subsequent evolution of viroids by complementing phylogenetic reconstructions with molecular data, including the primary and higher-order structure of the genomic RNAs, their replication, and recombination mechanisms and selected biological information. Features of some viroid-like RNAs found in plants, animals, and possibly fungi are also considered. The resulting evolutionary scenario supports the emergence of protoviroids in the RNA world, mainly as replicative modules, followed by a further increase in genome complexity based on module/domain shuffling and combination and mutation. Such a modular evolutionary scenario would have facilitated the inclusion in the protoviroid genomes of complex RNA structures (or coding sequences, as in the case of hepatitis delta virus and delta-like agents), likely needed for their adaptation from the RNA world to a life based on cells, thus generating the ancestors of current infectious viroids and viroid-like RNAs. Other noninfectious viroid-like RNAs, such as retroviroid-like RNA elements and retrozymes, could also be derived from protoviroids if their reverse transcription and integration into viral or eukaryotic DNA, respectively, are considered as a possible key step in their evolution. Comparison of evidence supporting a general and modular evolutionary model for viroids and viroid-like RNAs with that favoring alternative scenarios provides reasonable reasons to keep alive the hypothesis that these small RNA pathogens may be relics of a precellular world.

Rescue of an Infectious cDNA Clone of Barley Yellow Dwarf Virus-GAV

Barley yellow dwarf virus-GAV (BYDV-GAV) is one of the most prevalent viruses causing yellow dwarf disease in wheat in China. The biology and pathology of BYDV-GAV are well studied; however, gene functions and molecular mechanisms of BYDV-GAV disease development are unclear because of the lack of a reverse genetics system. In this study, a full-length complementary DNA (cDNA) clone of BYDV-GAV was constructed and expressed via Agrobacterium-mediated inoculation of Nicotiana benthamiana. Virions produced by BYDV-GAV in N. benthamiana were transmitted to wheat by an aphid vector after acquisition via a sandwich feeding method. Infectivity of the cDNA clone in wheat was verified via reverse transcription PCR and western blot assays, and the recombinant virus elicited typical reddening symptoms in oats and was transmitted between wheat plants. These results confirm the production of biologically active transmissible virions. Using the BYDV-GAV infectious clone, we demonstrate that viral protein P4 was involved in cell-to-cell movement and stunting symptoms in wheat. This is the first report describing the development of an infectious full-length cDNA clone of BYDV-GAV and provides a useful tool for virus-host-vector interaction studies.

The Borna Disease Virus 2 (BoDV-2) Nucleoprotein Is a Conspecific Protein That Enhances BoDV-1 RNA-Dependent RNA Polymerase Activity

An RNA virus-based episomal vector (REVec) based on Borna disease virus 1 (BoDV-1) is a promising viral vector that achieves stable and long-term gene expression in transduced cells. However, the onerous procedure of reverse genetics used to generate an REVec is one of the challenges that must be overcome to make REVec technologies practical for use. In this study, to resolve the problems posed by reverse genetics, we focused on BoDV-2, a conspecific virus of BoDV-1 in the Mammalian 1 orthobornavirus. We synthesized the BoDV-2 nucleoprotein (N) and phosphoprotein (P) according to the reference sequences and evaluated their effects on the RNA polymerase activity of the BoDV-1 large protein (L) and viral replication. In the minireplicon assay, we found that BoDV-2 N significantly enhanced BoDV-1 polymerase activity and that BoDV-2 P supported further enhancement of this activity by N. A single amino acid substitution assay identified serine at position 30 of BoDV-2 N and alanine at position 24 of BoDV-2 P as critical amino acid residues for the enhancement of BoDV-1 polymerase activity. In reverse genetics, conversely, BoDV-2 N alone was sufficient to increase the rescue efficiency of the REVec. We showed that the REVec can be rescued directly from transfected 293T cells by using BoDV-2 N as a helper plasmid without cocultivation with Vero cells and following several weeks of passage. In addition, a chimeric REVec harboring the BoDV-2 N produced much higher levels of transgene mRNA and genomic RNA than the wild-type REVec in transduced cells. Our results contribute to not only improvements to the REVec system but also to understanding of the molecular regulation of orthobornavirus polymerase activity. IMPORTANCE Borna disease virus 1 (BoDV-1), a prototype virus of the species Mammalian 1 orthobornavirus, is a nonsegmented negative-strand RNA virus that persists in the host nucleus. The nucleoprotein (N) of BoDV-1 encapsidates genomic and antigenomic viral RNA, playing important roles in viral transcription and replication. In this study, we demonstrated that the N of BoDV-2, another genotype in the species Mammalian 1 orthobornavirus, can participate in the viral ribonucleoprotein complex of BoDV-1 and enhance the activity of BoDV-1 polymerase (L) in both the BoDV-1 minireplicon assay and reverse genetics system. Chimeric recombinant BoDV-1 expressing BoDV-2 N but not BoDV-1 N showed higher transcription and replication levels, whereas the propagation and infectious particle production of the chimeric virus were comparable to those of wild-type BoDV-1, suggesting that the level of viral replication in the nucleus is not directly involved in the progeny virion production of BoDVs. Our results demonstrate a molecular mechanism of bornaviral polymerase activity, which will contribute to further development of vector systems using orthobornaviruses.

Experimental Resurrection of Ancestral Mammalian CPEB3 Ribozymes Reveals Deep Functional Conservation

Self-cleaving ribozymes are genetic elements found in all domains of life, but their evolution remains poorly understood. A ribozyme located in the second intron of the cytoplasmic polyadenylation binding protein 3 gene (CPEB3) shows high sequence conservation in mammals, but little is known about the functional conservation of self-cleaving ribozyme activity across the mammalian tree of life or during the course of mammalian evolution. Here, we use a phylogenetic approach to design a mutational library and a deep sequencing assay to evaluate the in vitro self-cleavage activity of numerous extant and resurrected CPEB3 ribozymes that span over 100 My of mammalian evolution. We found that the predicted sequence at the divergence of placentals and marsupials is highly active, and this activity has been conserved in most lineages. A reduction in ribozyme activity appears to have occurred multiple different times throughout the mammalian tree of life. The in vitro activity data allow an evaluation of the predicted mutational pathways leading to extant ribozyme as well as the mutational landscape surrounding these ribozymes. The results demonstrate that in addition to sequence conservation, the self-cleavage activity of the CPEB3 ribozyme has persisted over millions of years of mammalian evolution.

Reversion to ancestral Zika virus NS1 residues increases competence of Aedes albopictus

Both mosquito species-specific differences and virus strain -specific differences impact vector competence. Previous results in our laboratory with individual populations of N. American mosquitoes support studies suggesting Aedes aegypti are more competent than Ae. albopictus for American Zika virus (ZIKV) strains and demonstrate that U.S. Ae. albopictus have higher competence for an ancestral Asian ZIKV strain. A982V, an amino acid substitution in the NS1 gene acquired prior to the American outbreak, has been shown to increase competence in Ae. aegypti. We hypothesized that variability in the NS1 could therefore contribute to species-specific differences and developed a reverse genetics system based on a 2016 ZIKV isolate from Honduras (ZIKV-WTic) to evaluate the phenotypic correlates of individual amino acid substitutions. In addition to A982V, we evaluated G894A, which was acquired during circulation in the Americas. Reversion of 982 and 894 to ancestral residues increased infectivity, transmissibility and viral loads in Ae. albopictus but had no effect on competence or replication in Ae. aegypti. In addition, while host cell-specific differences in NS1 secretion were measured, with significantly higher secretion in mammalian cells relative to mosquito cells, strain-specific differences in secretion were not detected, despite previous reports. These results demonstrate that individual mutations in NS1 can influence competence in a species-specific manner independent of differences in NS1 secretion and further indicate that ancestral NS1 residues confer increased competence in Ae. albopictus. Lastly, experimental infections of Ifnar1-/- mice demonstrated that these NS1 substitutions can influence viral replication in the host and, specifically, that G894A could represent a compensatory change following a fitness loss from A982V with some viral genetic backgrounds. Together these data suggest a possible role for epistatic interactions in ZIKV fitness in invertebrate and vertebrate hosts and demonstrate that strains with increased transmission potential in U.S. Ae. albopictus could emerge.

Single-nucleotide control of tRNA folding cooperativity under near-cellular conditions

RNA folding is often studied by renaturing full-length RNA in vitro and tracking folding transitions. However, the intracellular transcript folds as it emerges from the RNA polymerase. Here, we investigate the folding pathways and stability of numerous late-transcriptional intermediates of yeast and Escherichia coli transfer RNAs (tRNAs). Transfer RNA is a highly regulated functional RNA that undergoes multiple steps of posttranscriptional processing and is found in very different lengths during its lifetime in the cell. The precursor transcript is extended on both the 5' and 3' ends of the cloverleaf core, and these extensions get trimmed before addition of the 3'-CCA and aminoacylation. We studied the thermodynamics and structures of the precursor tRNA and of late-transcriptional intermediates of the cloverleaf structure. We examined RNA folding at both the secondary and tertiary structural levels using multiple biochemical and biophysical approaches. Our findings suggest that perhaps nature has selected for a single-base addition to control folding to the functional 3D structure. In near-cellular conditions, yeast tRNA^Phe and E. coli tRNA^Ala transcripts fold in a single, cooperative transition only when nearly all of the nucleotides in the cloverleaf are transcribed by indirectly enhancing folding cooperativity. Furthermore, native extensions on the 5' and 3' ends do not interfere with cooperative core folding. This highly controlled cooperative folding has implications for recognition of tRNA by processing and modification enzymes and quality control of tRNA in cells.

Novel ribozymes: discovery, catalytic mechanisms, and the quest to understand biological function

Small endonucleolytic ribozymes promote the self-cleavage of their own phosphodiester backbone at a specific linkage. The structures of and the reactions catalysed by members of individual families have been studied in great detail in the past decades. In recent years, bioinformatics studies have uncovered a considerable number of new examples of known catalytic RNA motifs. Importantly, entirely novel ribozyme classes were also discovered, for most of which both structural and biochemical information became rapidly available. However, for the majority of the new ribozymes, which are found in the genomes of a variety of species, a biological function remains elusive. Here, we concentrate on the different approaches to find catalytic RNA motifs in sequence databases. We summarize the emerging principles of RNA catalysis as observed for small endonucleolytic ribozymes. Finally, we address the biological functions of those ribozymes, where relevant information is available and common themes on their cellular activities are emerging. We conclude by speculating on the possibility that the identification and characterization of proteins that we hypothesize to be endogenously associated with catalytic RNA might help in answering the ever-present question of the biological function of the growing number of genomically encoded, small endonucleolytic ribozymes.

Effects of flanking regions on HDV cotranscriptional folding kinetics

Hepatitis delta virus (HDV) ribozyme performs the self-cleavage activity through folding to a double pseudoknot structure. The folding of functional RNA structures is often coupled with the transcription process. In this work, we developed a new approach for predicting the cotranscriptional folding kinetics of RNA secondary structures with pseudoknots. We theoretically studied the cotranscriptional folding behavior of the 99-nucleotide (nt) HDV sequence, two upstream flanking sequences, and one downstream flanking sequence. During transcription, the 99-nt HDV can effectively avoid the trap intermediates and quickly fold to the cleavage-active state. It is different from its refolding kinetics, which folds into an intermediate trap state. For all the sequences, the ribozyme regions (from 1 to 73) all fold to the same structure during transcription. However, the existence of the 30-nt upstream flanking sequence can inhibit the ribozyme region folding into the active native state through forming an alternative helix Alt1 with the segments 70-90. The longer upstream flanking sequence of 54 nt itself forms a stable hairpin structure, which sequesters the formation of the Alt1 helix and leads to rapid formation of the cleavage-active structure. Although the 55-nt downstream flanking sequence could invade the already folded active structure during transcription by forming a more stable helix with the ribozyme region, the slow transition rate could keep the structure in the cleavage-active structure to perform the activity.

Analyzing and Tuning Ribozyme Activity by Deep Sequencing To Modulate Gene Expression Level in Mammalian Cells

Self-cleaving ribozymes, in combination with aptamers and various classes of RNAs, have been heavily engineered to create RNA devices to control gene expression. Although understanding of sequence-function relationships of ribozymes is critical for such efforts, our current knowledge of self-cleaving ribozymes is mostly limited to the results from small scale mutational studies performed under different conditions, or qualitative results of mutate-and-select experiments that may contain experimental biases. Here, we applied our strategy based on deep sequencing to comprehensively assay a large number of mutants to systematically examine the effect of the P4 stem sequence on the activity of an HDV-like ribozyme. We discovered that the ribozyme activity is highly sensitive to the sequence and the apparent stability of the varied positions. Furthermore, we demonstrated that the collection of the ribozyme variants with different activities can be used as a convenient device to fine-tune the level of gene expression in mammalian cells.

Co-Transcriptional Folding and Regulation Mechanisms of Riboswitches

Riboswitches are genetic control elements within non-coding regions of mRNA. These self-regulatory elements have been found to sense a range of small metabolites, ions, and other physical signals to exert regulatory control of transcription, translation, and splicing. To date, more than a dozen riboswitch classes have been characterized that vary widely in size and secondary structure. Extensive experiments and theoretical studies have made great strides in understanding the general structures, genetic mechanisms, and regulatory activities of individual riboswitches. As the ligand-dependent co-transcriptional folding and unfolding dynamics of riboswitches are the key determinant of gene expression, it is important to investigate the thermodynamics and kinetics of riboswitches both in the presence and absence of metabolites under the transcription. This review will provide a brief summary of the studies about the regulation mechanisms of the pbuE, S_MK, yitJ, and metF riboswitches based on the ligand-dependent co-transcriptional folding of the riboswitches.

Iluminated by foreign letters - Strategies for site-specific cyclopropene modification of large functional RNAs via in vitro transcription

The synthesis of sequence-specifically modified long RNA molecules, which cannot entirely be prepared via solid phase synthesis methods is experimentally challenging. We are using a new approach based on an expanded genetic alphabet preparing site-specifically modified RNA molecules via standard in vitro transcription. In this report, the site-specific labeling of functional RNAs, in particular ribozymes and a long non-coding RNA with cyclopropene moieties, is presented. We provide detailed instructions for RNA labeling via in vitro transcription and include required analytical methods to verify production and identity of the transcript. We further present post-transcriptional inverse electron demand Diels-Alder cycloaddition reactions on the cyclopropene-modified sequences and discuss applications of the genetic alphabet expansion transcription for in vitro preparation of labeled functional RNAs with complex foldings. In detail, the glmS and CPEB3 ribozymes were site-specifically decorated with methyl cyclopropene moieties using the unnatural TPT3^CP triphosphate and were proven to be still functional. In addition, the structurally complex A region of the Xist lncRNA (401nt) was site-specifically modified with methyl cyclopropene and detected by fluorescence after cycloaddition reaction with a tetrazine-BODIPY conjugate.

Probing fast ribozyme reactions under biological conditions with rapid quench-flow kinetics

Reaction kinetics on the millisecond timescale pervade the protein and RNA fields. To study such reactions, investigators often perturb the system with abiological solution conditions or substrates in order to slow the rate to timescales accessible by hand mixing; however, such perturbations can change the rate-limiting step and obscure key folding and chemical steps that are found under biological conditions. Mechanical methods for collecting data on the millisecond timescale, which allow these perturbations to be avoided, have been developed over the last few decades. These methods are relatively simple and can be conducted on affordable and commercially available instruments. Here, we focus on using the rapid quench-flow technique to study the fast reaction kinetics of RNA enzymes, or ribozymes, which often react on the millisecond timescale under biological conditions. Rapid quench of ribozymes is completely parallel to the familiar hand-mixing approach, including the use of radiolabeled RNAs and fractionation of reactions on polyacrylamide gels. We provide tips on addressing and preventing common problems that can arise with the rapid-quench technique. Guidance is also offered on ensuring the ribozyme is properly folded and fast-reacting. We hope that this article will facilitate the broader use of rapid-quench instrumentation to study fast-reacting ribozymes under biological reaction conditions.

Genome-Wide Analysis of RNA Secondary Structure

Single-stranded RNA molecules fold into extraordinarily complicated secondary and tertiary structures as a result of intramolecular base pairing. In vivo, these RNA structures are not static. Instead, they are remodeled in response to changes in the prevailing physicochemical environment of the cell and as a result of intermolecular base pairing and interactions with RNA-binding proteins. Remarkable technical advances now allow us to probe RNA secondary structure at single-nucleotide resolution and genome-wide, both in vitro and in vivo. These data sets provide new glimpses into the RNA universe. Analyses of RNA structuromes in HIV, yeast, Arabidopsis, and mammalian cells and tissues have revealed regulatory effects of RNA structure on messenger RNA (mRNA) polyadenylation, splicing, translation, and turnover. Application of new methods for genome-wide identification of mRNA modifications, particularly methylation and pseudouridylation, has shown that the RNA "epitranscriptome" both influences and is influenced by RNA structure. In this review, we describe newly developed genome-wide RNA structure-probing methods and synthesize the information emerging from their application.

Bridging the gap between in vitro and in vivo RNA folding

Deciphering the folding pathways and predicting the structures of complex three-dimensional biomolecules is central to elucidating biological function. RNA is single-stranded, which gives it the freedom to fold into complex secondary and tertiary structures. These structures endow RNA with the ability to perform complex chemistries and functions ranging from enzymatic activity to gene regulation. Given that RNA is involved in many essential cellular processes, it is critical to understand how it folds and functions in vivo. Within the last few years, methods have been developed to probe RNA structures in vivo and genome-wide. These studies reveal that RNA often adopts very different structures in vivo and in vitro, and provide profound insights into RNA biology. Nonetheless, both in vitro and in vivo approaches have limitations: studies in the complex and uncontrolled cellular environment make it difficult to obtain insight into RNA folding pathways and thermodynamics, and studies in vitro often lack direct cellular relevance, leaving a gap in our knowledge of RNA folding in vivo. This gap is being bridged by biophysical and mechanistic studies of RNA structure and function under conditions that mimic the cellular environment. To date, most artificial cytoplasms have used various polymers as molecular crowding agents and a series of small molecules as cosolutes. Studies under such in vivo-like conditions are yielding fresh insights, such as cooperative folding of functional RNAs and increased activity of ribozymes. These observations are accounted for in part by molecular crowding effects and interactions with other molecules. In this review, we report milestones in RNA folding in vitro and in vivo and discuss ongoing experimental and computational efforts to bridge the gap between these two conditions in order to understand how RNA folds in the cell.

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.

Comprehensive analysis of mutations in the hepatitis delta virus genome based on full-length sequencing in a nationwide cohort study and evolutionary pattern during disease progression

Delta hepatitis, caused by co-infection or super-infection of hepatitis D virus (HDV) in hepatitis B virus (HBV) -infected patients, is the most severe form of chronic hepatitis, often progressing to liver cirrhosis and liver failure. Although 15 million individuals are affected worldwide, molecular data on the HDV genome and its proteins, small and large delta antigen (S-/L-HDAg), are limited. We therefore conducted a nationwide study in HBV-HDV-infected patients from Iran and successfully amplified 38 HDV full genomes and 44 L-HDAg sequences from 34 individuals. Phylogenetic analyses of full-length HDV and L-HDAg isolates revealed that all strains clustered with genotype 1 and showed high genotypic distances to HDV genotypes 2 to 8, with a maximal distance to genotype 3. Longitudinal analyses in individual patients indicated a reverse evolutionary trend, especially in L-HDAg amino acid composition, over time. Besides multiple sequence variations in the hypervariable region of HDV, nucleotide substitutions preferentially occurred in the stabilizing P4 domain of the HDV ribozyme. A high rate of single amino acid changes was detected in structural parts of L-HDAg, whereas its post-translational modification sites were highly conserved. Interestingly, several non-synonymous mutations were positively selected that affected immunogenic epitopes of L-HDAg towards CD8 T-cell- and B-cell-driven immune responses. Hence, our comprehensive molecular analysis comprising a nationwide cohort revealed phylogenetic relationships and provided insight into viral evolution within individual hosts. Moreover, preferential areas of frequent mutations in the HDV ribozyme and antigen protein were determined in this study.

Disparate HDV ribozyme crystal structures represent intermediates on a rugged free-energy landscape

The hepatitis delta virus (HDV) ribozyme is a member of the class of small, self-cleaving catalytic RNAs found in a wide range of genomes from HDV to human. Both pre- and post-catalysis (precursor and product) crystal structures of the cis-acting genomic HDV ribozyme have been determined. These structures, together with extensive solution probing, have suggested that a significant conformational change accompanies catalysis. A recent crystal structure of a trans-acting precursor, obtained at low pH and by molecular replacement from the previous product conformation, conforms to the product, raising the possibility that it represents an activated conformer past the conformational change. Here, using fluorescence resonance energy transfer (FRET), we discovered that cleavage of this ribozyme at physiological pH is accompanied by a structural lengthening in magnitude comparable to previous trans-acting HDV ribozymes. Conformational heterogeneity observed by FRET in solution appears to have been removed upon crystallization. Analysis of a total of 1.8 µsec of molecular dynamics (MD) simulations showed that the crystallographically unresolved cleavage site conformation is likely correctly modeled after the hammerhead ribozyme, but that crystal contacts and the removal of several 2'-oxygens near the scissile phosphate compromise catalytic in-line fitness. A cis-acting version of the ribozyme exhibits a more dynamic active site, while a G-1 residue upstream of the scissile phosphate favors poor fitness, allowing us to rationalize corresponding changes in catalytic activity. Based on these data, we propose that the available crystal structures of the HDV ribozyme represent intermediates on an overall rugged RNA folding free-energy landscape.

The wide expansion of hepatitis delta virus-like ribozymes throughout trypanosomatid genomes is linked to the spreading of L1Tc/ingi clade mobile elements

Background

Hepatitis Delta Virus (HDV)-like ribozymes have recently been found in many mobile elements in which they take part in a mechanism that releases intermediate RNAs from cellular co-transcripts. L1Tc in Trypanosoma cruzi is one of the elements in which such a ribozyme is located. It lies in the so-called Pr77-hallmark, a conserved region shared by retrotransposons belonging to the trypanosomatid L1Tc/ingi clade. The wide distribution of the Pr77-hallmark detected in trypanosomatid retrotransposons renders the potential catalytic activity of these elements worthy of study: their distribution might contribute to host genetic regulation at the mRNA level. Indeed, in Leishmania spp, the pervasive presence of these HDV-like ribozyme-containing mobile elements in certain 3′-untranslated regions of protein-coding genes has been linked to mRNA downregulation.

Results

Intensive screening of publicly available trypanosomatid genomes, combined with manual folding analyses, allowed the isolation of putatively Pr77-hallmarks with HDV-like ribozyme activity. This work describes the conservation of an HDV-like ribozyme structure in the Pr77 sequence of retrotransposons in a wide range of trypanosomatids, the catalytic function of which is maintained in the majority.

These results are consistent with the previously suggested common phylogenetic origin of the elements that belong to this clade, although in some cases loss of functionality appears to have occurred and/or perhaps molecular domestication by the host.

Conclusions

These HDV-like ribozymes are widely distributed within retrotransposons across trypanosomatid genomes. This type of ribozyme was once thought to be rare in nature, but in fact it would seem to be abundant in trypanosomatid transcripts. It can even form part of the pool of mRNA 3′-untranslated regions, particularly in Leishmania spp. Its putative regulatory role in host genetic expression is discussed.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-340) contains supplementary material, which is available to authorized users.

HDV family of self-cleaving ribozymes

The hepatitis delta virus (HDV) ribozymes are catalytic RNAs capable of cleaving their own sugar-phosphate backbone. The HDV virus possesses the ribozymes in both sense and antisense genomic transcripts, where they are essential for processing during replication. These ribozymes have been the subject of intense biochemical scrutiny and have yielded a wealth of mechanistic insights. In recent years, many HDV-like ribozymes have been identified in nearly all branches of life. The ribozymes are implicated in a variety of biological events, including episodic memory in mammals and retrotransposition in many eukaryotes. Detailed analysis of additional HDV-like ribozyme isolates will likely reveal many more biological functions and provide information about the evolution of this unique RNA.

The hok mRNA family

The hok/sok toxin-antitoxin system of Escherichia coli plasmid R1 increases plasmid maintenance by killing plasmid-free daughter cells. The hok/sok locus specifies two RNAs: hok mRNA, which encodes a toxic transmembrane protein, and sok antisense RNA, which binds a complementary region in the hok mRNA and induces transcript degradation. During cell growth, the cis-encoded sok RNA inhibits expression of the Hok toxin. In plasmid-free segregants, the rapid decay of sok RNA relative to hok mRNA permits Hok translation, leading to cell death. This post-segregational killing mechanism relies upon the ability of the hok mRNA to adopt alternative structural configurations, which affect ease of translation and the susceptibility of the molecule to degradation. The full-length hok transcript is stable, highly structured and immune to ribosome and antisense RNA binding. Gradual 3' end processing produces dramatic structural rearrangements in the mRNA, which render the molecule translationally active and expose the sok RNA binding site. During transcription, premature ribosome and sok binding are prevented through the formation of transient metastable hairpins in the 5' end of the nascent transcript. Several hok mRNA paralogs have been identified in the genome of E. coli, and Hok protein orthologs found in the genomes of Enterobacteria. Using a combination of automated search and extensive manual editing, we compiled a multiple sequence alignment for the hok mRNA. All three experimentally validated hok mRNA structures are mapped onto this alignment, which has been submitted to the Rfam database for RNA families.

Mechanistic analysis of the hepatitis delta virus (HDV) ribozyme: methods for RNA preparation, structure mapping, solvent isotope effects, and co-transcriptional cleavage

Small ribozymes such as the hairpin, hammerhead, VS, glm S, and hepatitis delta virus (HDV) are self-cleaving RNAs that are typically characterized by kinetics and structural methods. Working with these RNAs requires attention to numerous experimental details. In this chapter we focus on four different experimental aspects of ribozyme studies: preparing the RNA, mapping its structure with reverse transcription and end-labeled techniques, solvent isotope experiments, and co-transcriptional cleavage assays. Although the focus of these methods is the HDV ribozyme, the methods should be applicable to other ribozymes.

Identification of an hepatitis delta virus-like ribozyme at the mRNA 5'-end of the L1Tc retrotransposon from Trypanosoma cruzi

L1Tc is a non-LTR LINE element from Trypanosoma cruzi that encodes its transposition machinery and bears an internal promoter. Herewith, we report the identification of an in vitro active hepatitis delta virus-like ribozyme located in the first 77 nt at the 5′-end of the L1Tc mRNA (L1TcRz). The data presented show that L1TcRz has a co-transcriptional function. Using gel-purified uncleaved RNA transcripts, the data presented indicate that the kinetics of the self-cleaving, in a magnesium-dependent reaction, fits to a two-phase decay curve. The cleavage point identified by primer extension takes place at +1 position of the element. The hydroxyl nature of the 5′-end of the 3′-fragment generated by the cleavage activity of L1TcRz was confirmed. Since we have previously described that the 77-nt long fragment located at the 5′-end of L1Tc has promoter activity, the existence of a ribozyme in L1Tc makes this element to be the first described non-LTR retroelement that has an internal promoter–ribozyme dual function. The L1Tc nucleotides located downstream of the ribozyme catalytic motif appear to inhibit its activity. This inhibition may be influenced by the existence of a specific L1Tc RNA conformation that is recognized by RNase P.

Mechanistic strategies in the HDV ribozyme: chelated and diffuse metal ion interactions and active site protonation

The crystal structure of the precleaved form of the hepatitis delta virus (HDV) ribozyme reveals two G•U wobbles near the active site: a rare reverse G•U wobble involving a syn G base, and a standard G•U wobble at the cleavage site. The catalytic mechanism for this ribozyme has been proposed to involve a Mg(2+) ion bound to the reverse G•U wobble, as well as a protonated C75 base. We carried out molecular dynamics simulations to analyze metal ion interaction with the reverse and standard G•U wobbles and to investigate the impact of C75 protonation on the structure and motions of the ribozyme. We identified two types of Mg(2+) ions associated with the ribozyme, chelated and diffuse, at the reverse and standard G•U wobbles, respectively, which appear to contribute to catalysis and stability, respectively. These two metal ion sites exhibit relatively independent behavior. Protonation of C75 was observed to locally organize the active site in a manner that facilitates the catalytic mechanism, in which C75(+) acts as a general acid and Mg(2+) as a Lewis acid. The simulations also indicated that the overall structure and thermal motions of the ribozyme are not significantly influenced by the catalytic Mg(2+) interaction or C75 protonation. This analysis suggests that the reaction pathway of the ribozyme is dominated by small local motions at the active site rather than large-scale global conformational changes. These results are consistent with a wealth of experimental data.

Roles of DEAD-box proteins in RNA and RNP Folding

RNAs and RNA-protein complexes (RNPs) traverse rugged energy landscapes as they fold to their native structures, and many continue to undergo conformational rearrangements as they function. Due to the inherent stability of local RNA structure, proteins are required to assist with RNA conformational transitions during initial folding and in exchange between functional structures. DEAD-box proteins are superfamily 2 RNA helicases that are ubiquitously involved in RNA-mediated processes. Some of these proteins use an ATP-dependent cycle of conformational changes to disrupt RNA structure nonprocessively, accelerating structural transitions of RNAs and RNPs in a manner that bears a strong resemblance to the activities of certain groups of protein chaperones. This review summarizes recent work using model substrates and tractable self-splicing intron RNAs, which has given new insights into how DEAD-box proteins promote RNA folding steps and conformational transitions, and it summarizes recent progress in identifying sites and mechanisms of DEAD-box protein activity within more complex cellular targets.

The human HDV-like CPEB3 ribozyme is intrinsically fast-reacting

Self-cleaving RNAs have recently been identified in mammalian genomes. A small ribozyme related in structure to the hepatitis delta virus (HDV) ribozyme occurs in a number of mammals, including chimpanzees and humans, within an intron of the CPEB3 gene. The catalytic mechanisms for the CPEB3 and HDV ribozymes appear to be similar, generating cleavage products with 5'-hydroxyl and 2',3'-cyclic phosphate termini; nonetheless, the cleavage rate reported for the CPEB3 ribozyme is more than 6000-fold slower than for the fastest HDV ribozyme. Herein, we use full-length RNA and cotranscriptional self-cleavage assays to compare reaction rates among human CPEB3, chimp CPEB3, and HDV ribozymes. Our data reveal that a single base change of the upstream flanking sequence, which sequesters an intrinsically weak P1.1 pairing in a misfold, increases the rate of the wild-type human CPEB3 ribozyme by approximately 250-fold; thus, the human ribozyme is intrinsically fast-reacting. Secondary structure determination and native gel analyses reveal that the cleaved population of the CPEB3 ribozyme has a single, secondary structure that closely resembles the HDV ribozyme. In contrast, the precleavage population of the CPEB3 ribozyme appears to have a more diverse secondary structure, possibly reflecting misfolding with the upstream sequence and dynamics intrinsic to the ribozyme. Prior identification of expressed sequence tags (ESTs) in human cells indicated that cleavage activity of the human ribozyme is tissue-specific. It is therefore possible that cellular factors interact with regions upstream of the CPEB3 ribozyme to unmask its high intrinsic reactivity.

Long-distance communication in the HDV ribozyme: insights from molecular dynamics and experiments

The hepatitis delta virus ribozyme is a small, self-cleaving RNA with a compact tertiary structure and buried active site that is important in the life cycle of the virus. The ribozyme's function in nature is to cleave an internal phosphodiester bond and linearize concatemers during rolling circle replication. Crystal structures of the ribozyme have been solved in both pre-cleaved and post-cleaved (product) forms and reveal an intricate network of interactions that conspire to catalyze bond cleavage. In addition, extensive biochemical studies have been performed to work out a mechanism for bond cleavage in which C75 and a magnesium ion catalyze the reaction by general acid-base chemistry. One issue that has remained unclear in this ribozyme and in other ribozymes is the nature of long-distance communication between peripheral regions of the RNA and the buried active site. We performed molecular dynamics simulations on the hepatitis delta virus ribozyme in the product form and assessed communication between a distal structural portion of the ribozyme-the protonated C41 base triple-and the active site containing the critical C75. We varied the ionization state of C41 in both the wild type and a C41 double mutant variant and determined the impact on the active site. In all four cases, effects at the active site observed in the simulations agree with experimental studies on ribozyme activity. Overall, these studies indicate that small functional RNAs have the potential to communicate interactions over long distances and that wild-type RNAs may have evolved ways to prevent such interactions from interfering with catalysis.

Efficient inhibition of hepatitis B virus replication by hepatitis delta virus ribozymes delivered by targeting retrovirus

Background

Hepatitis delta virus (HDV) ribozyme is an attractive molecular tool that can specifically recognize and catalyze the self-cleavage of the viral RNA phosphodiester backbone. However, a major obstacle in the medical application of the HDV ribozyme is the lack of specificity in the delivery of the ribozyme to defined target cells.

Results

The objective of this study was to determine whether retroviral vectors can deliver the HDV ribozyme into the target cells and to elucidate whether HDV ribozyme plays a role in hepatitis B virus (HBV) replication. In our study, the transduction of helper-free pseudotyped retrovirus, which showed a broad host range, in human hepatoma cells was performed under 2 conditions, that is, in the presence of polymerized human serum albumin (pHSA) and in the absence of pHSA. The transduction ability in the presence of pHSA was higher than in the absence of pHSA. Moreover, HBsAg and HBeAg levels after transductions with pHSA were significantly lower than those in the absence of pHSA, thus indicating that the recombinant retrovirus had HBV-specific cleavage activity and targeted HepG2215 cells.

Conclusions

These data suggest that this system provides a new approach for targeting hepatocytes and has a great potential in gene therapy for HBV infection.

Raman crystallography of RNA

Raman crystallography is the application of Raman spectroscopy to single crystals. This technique has been applied to a variety of protein molecules where it has provided unique information about biopolymer folding, substrate binding, and catalysis. Here, we describe the application of Raman crystallography to functional RNA molecules. RNA represents unique opportunities and challenges for Raman crystallography. One issue that confounds studies of RNA is its tendency to adopt multiple non-functional folds. Raman crystallography has the advantage that it isolates a single state of the RNA within the crystal and can evaluate its fold, metal ion binding properties (ligand identity, stoichiometry, and affinity), proton binding properties (identity, stoichiometry, and affinity), and catalytic potential. In particular, base-specific stretches can be identified and then associated with the binding of metal ions and protons. Because measurements are carried out in the hanging drop at ambient, rather than cryo, conditions and because RNA crystals tend to be approximately 70% solvent, RNA dynamics and conformational changes become experimentally accessible. This review focuses on experimental setup and procedures, acquisition and interpretation of Raman data, and determination of physicochemical properties of the RNA. Raman crystallographic and solution biochemical experiments on the HDV RNA enzyme are summarized and found to be in excellent agreement. Remarkably, characterization of the crystalline state has proven to help rather than hinder functional characterization of functional RNA, most likely because the tendency of RNA to fold heterogeneously is limited in a crystalline environment. Future applications of Raman crystallography to RNA are briefly discussed.

Comparative enzymology and structural biology of RNA self-cleavage

Self-cleaving hammerhead, hairpin, hepatitis delta virus, and glmS ribozymes comprise a family of small catalytic RNA motifs that catalyze the same reversible phosphodiester cleavage reaction, but each motif adopts a unique structure and displays a unique array of biochemical properties. Recent structural, biochemical, and biophysical studies of these self-cleaving RNAs have begun to reveal how active site nucleotides exploit general acid-base catalysis, electrostatic stabilization, substrate destabilization, and positioning and orientation to reduce the free energy barrier to catalysis. Insights into the variety of catalytic strategies available to these model RNA enzymes are likely to have important implications for understanding more complex RNA-catalyzed reactions fundamental to RNA processing and protein synthesis.

The fraction of RNA that folds into the correct branched secondary structure determines hepatitis delta virus type 3 RNA editing levels

RNA editing by the host RNA adenosine deaminase ADAR1 at the amber/W site of hepatitis delta virus RNA plays a central role in the viral replication cycle by affecting the balance between viral RNA synthesis and packaging. Previously, we found that HDV genotype III (HDV-3) RNA can form two secondary structures following transcription: an unbranched rod structure, which is characteristic of HDV, and a metastable branched structure that serves as the substrate for editing. The unstable nature of the branched editing substrate structure raised the possibility that structural dynamics of the RNA following transcription could determine the rate at which editing occurs. Here, editing and its control are examined in two HDV-3 isolates, from Peru and Ecuador. Analysis of editing in vitro by ADAR1 indicated that the branched structure formed by RNA derived from the Peruvian isolate is edited more efficiently than that from the Ecuadorian isolate. In contrast, in the context of replication, Peruvian RNA is edited less efficiently than RNA containing Ecuadorian sequences. Computational analyses of RNA folding using the massively parallel genetic algorithm (MPGAfold) indicated that the Peruvian RNA is less likely to form the branched structure required for editing than the Ecuadorian isolate. This difference was confirmed by in vitro transcription of these RNAs. Overall, our data indicate that HDV-3 controls RNA editing levels via (1) the fraction of the RNA that folds, during transcription, into the metastable branched structure required for editing and (2) the efficiency with which ADAR1 edits this branched substrate RNA.

A catalytic metal ion interacts with the cleavage Site G.U wobble in the HDV ribozyme

The HDV ribozyme self-cleaves by a chemical mechanism involving general acid-base catalysis to generate 2',3'-cyclic phosphate and 5'-hydroxyl termini. Biochemical studies from several laboratories have implicated C75 as the general acid and hydrated magnesium as the general base. We have previously shown that C75 has a pK(a) shifted >2 pH units toward neutrality [Gong, B., Chen, J. H., Chase, E., Chadalavada, D. M., Yajima, R., Golden, B. L., Bevilacqua, P. C., and Carey, P. R. (2007) J. Am. Chem. Soc. 129, 13335-13342], while in crystal structures, it is well-positioned for proton transfer. However, no evidence for a hydrated magnesium poised to serve as a general base in the reaction has been observed in high-resolution crystal structures of various reaction states and mutants. Herein, we use solution kinetic experiments and parallel Raman crystallographic studies to examine the effects of pH on the rate and Mg(2+) binding properties of wild-type and 7-deazaguanosine mutants of the HDV ribozyme. These data suggest that a previously unobserved hydrated magnesium ion interacts with N7 of the cleavage site G.U wobble base pair. Integrating this metal ion binding site with the available crystal structures provides a new three-dimensional model for the active site of the ribozyme that accommodates all available biochemical data and appears competent for catalysis. The position of this metal is consistent with a role of a magnesium-bound hydroxide as a general base as dictated by biochemical data.

Mechanistic characterization of the HDV genomic ribozyme: the cleavage site base pair plays a structural role in facilitating catalysis

The hepatitis delta virus (HDV) ribozyme occurs in the genomic and antigenomic strands of the HDV RNA and within mammalian transcriptomes. Previous kinetic studies suggested that a wobble pair (G*U or A(+)*C) is preferred at the cleavage site; however, the reasons for this are unclear. We conducted sequence comparisons, which indicated that while G*U is the most prevalent combination at the cleavage site, G-C occurs to a significant extent in genomic HDV isolates, and G*U, G-C, and A-U pairs are present in mammalian ribozymes. We analyzed the folding of genomic HDV ribozymes by free energy minimization and found that variants with purine-pyrimidine combinations at the cleavage site are predicted to form native structures while pyrimidine-purine combinations misfold, consistent with earlier kinetic data and sequence comparisons. To test whether the cleavage site base pair contributes to catalysis, we characterized the pH and Mg(2+)-dependence of reaction kinetics of fast-folding genomic HDV ribozymes with cleavage site base pair purine-pyrimidine combinations: G*U, A-U, G-C, and A(+)*C. Rates for these native-folding ribozymes displayed highly similar pH and Mg(2+) concentration dependencies, with the exception of the A(+)*C ribozyme, which deviated at high pH. None of the four ribozymes underwent miscleavage. These observations support the A(+)*C ribozyme as being more active with a wobble pair at the cleavage site than with no base pair at all. Overall, the data support a model in which the cleavage site base pair provides a structural role in catalysis and does not need to be a wobble pair.