A mechanistic framework for co-transcriptional folding of the HDV genomic ribozyme in the presence of downstream sequence

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.

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.

Hepatitis delta virus-like circular RNAs from diverse metazoans encode conserved hammerhead ribozymes

Human hepatitis delta virus (HDV) is a unique infectious agent whose genome is composed of a small circular RNA. Recent data, however, have reported the existence of highly divergent HDV-like circRNAs in the transcriptomes of diverse vertebrate and invertebrate species. The HDV-like genomes described in amniotes such as birds and reptiles encode self-cleaving RNA motifs or ribozymes similar to the ones present in the human HDV, whereas no catalytic RNA domains have been reported for the HDV-like genomes detected in metagenomic data from some amphibians, fish, and invertebrates. Herein, we describe the self-cleaving motifs of the HDV-like genomes reported in newts and fish, which belong to the characteristic class of HDV ribozymes. Surprisingly, HDV-like genomes from a toad and a termite show conserved type III hammerhead ribozymes, which belong to an unrelated class of catalytic RNAs characteristic of plant genomes and plant subviral circRNAs, such as some viral satellites and viroids. Sequence analyses revealed the presence of similar HDV-like hammerhead ribozymes encoded in two termite genomes, but also in the genomes of several dipteran species. In vitro transcriptions confirmed the cleaving activity for these motifs, with moderate rates of self-cleavage. These data indicate that all described HDV-like agents contain self-cleaving motifs from either the HDV or the hammerhead class. Autocatalytic ribozymes in HDV-like genomes could be regarded as interchangeable domains and may have arisen from cellular transcriptomes, although we still cannot rule out some other evolutionary explanations.

Co-transcriptional Analysis of Self-Cleaving Ribozymes and Their Ligand Dependence

Self-cleaving ribozymes are RNA molecules that catalyze a site-specific self-scission reaction. Analysis of self-cleavage is a crucial aspect of the biochemical study and understanding of these molecules. Here we describe a co-transcriptional assay that allows the analysis of self-cleaving ribozymes in different reaction conditions and in the presence of desired ligands and/or cofactors. Utilizing a standard T7 RNA polymerase in vitro transcription system under limiting Mg²⁺ concentration, followed by a 25-fold dilution of the reaction in desired conditions of self-cleavage (buffer, ions, ligands, pH, temperature, etc.) to halt the synthesis of new RNA molecules, allows the study of self-scission of these molecules without the need for purification or additional preparation steps, such as refolding procedures. Furthermore, because the transcripts are not denatured, this assay likely yields RNAs in conformations relevant to co-transcriptionally folded species in vivo.

Single-pass transcription by T7 RNA polymerase

RNA molecules can be conveniently synthesized in vitro by the T7 RNA polymerase (T7 RNAP). In some experiments, such as cotranscriptional biochemical analyses, continuous synthesis of RNA is not desired. Here, we propose a method for a single-pass transcription that yields a single transcript per template DNA molecule using the T7 RNAP system. We hypothesized that stalling the polymerase downstream from the promoter region and subsequent cleavage of the promoter by a restriction enzyme (to prevent promoter binding by another polymerase) would allow synchronized production of a single transcript per template. The single-pass transcription was verified in two different scenarios: a short self-cleaving ribozyme and a long mRNA. The results show that a controlled single-pass transcription using T7 RNAP allows precise measurement of cotranscriptional ribozyme activity, and this approach will facilitate the study of other kinetic events.

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.

Computational Methods for Modeling Aptamers and Designing Riboswitches

Riboswitches, which are located within certain noncoding RNA region perform functions as genetic "switches", regulating when and where genes are expressed in response to certain ligands. Understanding the numerous functions of riboswitches requires computation models to predict structures and structural changes of the aptamer domains. Although aptamers often form a complex structure, computational approaches, such as RNAComposer and Rosetta, have already been applied to model the tertiary (three-dimensional (3D)) structure for several aptamers. As structural changes in aptamers must be achieved within the certain time window for effective regulation, kinetics is another key point for understanding aptamer function in riboswitch-mediated gene regulation. The coarse-grained self-organized polymer (SOP) model using Langevin dynamics simulation has been successfully developed to investigate folding kinetics of aptamers, while their co-transcriptional folding kinetics can be modeled by the helix-based computational method and BarMap approach. Based on the known aptamers, the web server Riboswitch Calculator and other theoretical methods provide a new tool to design synthetic riboswitches. This review will represent an overview of these computational methods for modeling structure and kinetics of riboswitch aptamers and for designing riboswitches.

Deep Sequencing Analysis of Aptazyme Variants Based on a Pistol Ribozyme

Chemically regulated self-cleaving ribozymes, or aptazymes, are emerging as a promising class of genetic devices that allow dynamic control of gene expression in synthetic biology. However, further expansion of the limited repertoire of ribozymes and aptamers, and development of new strategies to couple the RNA elements to engineer functional aptazymes are highly desirable for synthetic biology applications. Here, we report aptazymes based on the recently identified self-cleaving pistol ribozyme class using a guanine aptamer as the molecular sensing element. Two aptazyme architectures were studied by constructing and assaying 17 728 mutants by deep sequencing. Although one of the architectures did not yield functional aptazymes, a novel aptazyme design in which the aptamer and the ribozyme were placed in tandem yielded a number of guanine-inhibited ribozymes. Detailed analysis of the extensive sequence-function data suggests a mechanism that involves a competition between two mutually exclusive RNA structures reminiscent of natural bacterial riboswitches.

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.

NMR Structural Profiling of Transcriptional Intermediates Reveals Riboswitch Regulation by Metastable RNA Conformations

Gene repression induced by the formation of transcriptional terminators represents a prime example for the coupling of RNA synthesis, folding, and regulation. In this context, mapping the changes in available conformational space of transcription intermediates during RNA synthesis is important to understand riboswitch function. A majority of riboswitches, an important class of small metabolite-sensing regulatory RNAs, act as transcriptional regulators, but the dependence of ligand binding and the subsequent allosteric conformational switch on mRNA transcript length has not yet been investigated. We show a strict fine-tuning of binding and sequence-dependent alterations of conformational space by structural analysis of all relevant transcription intermediates at single-nucleotide resolution for the I-A type 2'dG-sensing riboswitch from Mesoplasma florum by NMR spectroscopy. Our results provide a general framework to dissect the coupling of synthesis and folding essential for riboswitch function, revealing the importance of metastable states for RNA-based gene regulation.

Using RNA as Molecular Code for Programming Cellular Function

RNA is involved in a wide-range of important molecular processes in the cell, serving diverse functions: regulatory, enzymatic, and structural. Together with its ease and predictability of design, these properties can lead RNA to become a useful handle for biological engineers with which to control the cellular machinery. By modifying the many RNA links in cellular processes, it is possible to reprogram cells toward specific design goals. We propose that RNA can be viewed as a molecular programming language that, together with protein-based execution platforms, can be used to rewrite wide ranging aspects of cellular function. In this review, we catalogue developments in the use of RNA parts, methods, and associated computational models that have contributed to the programmability of biology. We discuss how RNA part repertoires have been combined to build complex genetic circuits, and review recent applications of RNA-based parts and circuitry. We explore the future potential of RNA engineering and posit that RNA programmability is an important resource for firmly establishing an era of rationally designed synthetic biology.

Detecting riboSNitches with RNA folding algorithms: a genome-wide benchmark

Ribonucleic acid (RNA) secondary structure prediction continues to be a significant challenge, in particular when attempting to model sequences with less rigidly defined structures, such as messenger and non-coding RNAs. Crucial to interpreting RNA structures as they pertain to individual phenotypes is the ability to detect RNAs with large structural disparities caused by a single nucleotide variant (SNV) or riboSNitches. A recently published human genome-wide parallel analysis of RNA structure (PARS) study identified a large number of riboSNitches as well as non-riboSNitches, providing an unprecedented set of RNA sequences against which to benchmark structure prediction algorithms. Here we evaluate 11 different RNA folding algorithms’ riboSNitch prediction performance on these data. We find that recent algorithms designed specifically to predict the effects of SNVs on RNA structure, in particular remuRNA, RNAsnp and SNPfold, perform best on the most rigorously validated subsets of the benchmark data. In addition, our benchmark indicates that general structure prediction algorithms (e.g. RNAfold and RNAstructure) have overall better performance if base pairing probabilities are considered rather than minimum free energy calculations. Although overall aggregate algorithmic performance on the full set of riboSNitches is relatively low, significant improvement is possible if the highest confidence predictions are evaluated independently.

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 kinetics of ribozyme cleavage: a tool to analyze RNA folding as a function of catalysis

As catalytically active RNAs, ribozymes can be characterized by kinetic measurements similar to classical enzyme kinetics. However, in contrast to standard protein enzymes, for which reactions can usually be started by mixing the enzyme with its substrate, ribozymes are typically self-cleaving. The reaction has to be initiated by folding the RNA into its active conformation. Thus, ribozyme kinetics are influenced by both folding and catalytic components and often enable indirect observation of RNA folding. Here, I describe how to obtain quantitative ribozyme cleavage data via denaturing polyacrylamide gel electrophoresis (PAGE) of radioactively labeled in vitro transcripts and discuss general considerations for subsequent kinetic analysis.

Activation of PKR by RNA misfolding: HDV ribozyme dimers activate PKR

Protein Kinase R (PKR), the double-stranded RNA (dsRNA)-activated protein kinase, plays important roles in innate immunity. Previous studies have shown that PKR is activated by long stretches of dsRNA, RNA pseudoknots, and certain single-stranded RNAs; however, regulation of PKR by RNAs with globular tertiary structure has not been reported. In this study, the HDV ribozyme is used as a model of a mostly globular RNA. In addition to a catalytic core, the ribozyme contains a peripheral 13-bp pairing region (P4), which, upon shortening, affects neither the catalytic activity of the ribozyme nor its ability to crystallize. We report that the HDV ribozyme sequence alone can activate PKR. To elucidate the RNA structural basis for this, we prepared a number of HDV variants, including those with shortened or lengthened P4 pairing regions, with the anticipation that lengthening the P4 extension would yield a more potent activator since it would offer more base pairs of dsRNA. Surprisingly, the variant with a shortened P4 was the most potent activator. Through native gel mobility and enzymatic structure mapping experiments we implicate misfolded HDV ribozyme dimers as the PKR-activating species, and show that the shortened P4 leads to enhanced occupancy of the RNA dimer. These observations have implications for how RNA misfolding relates to innate immune response and human disease.

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.

RNA folding in living cells

RNA folding is the most essential process underlying RNA function. While significant progress has been made in understanding the forces driving RNA folding in vitro, exploring the rules governing intracellular RNA structure formation is still in its infancy. The cellular environment hosts a great diversity of factors that potentially influence RNA folding in vivo. For example, the nature of transcription and translation is known to shape the folding landscape of RNA molecules. Trans-acting factors such as proteins, RNAs and metabolites, among others, are also able to modulate the structure and thus the fate of an RNA. Here we summarize the ongoing efforts to uncover how RNA folds in living cells.

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.

mRNA secondary structures fold sequentially but exchange rapidly in vivo

RNAs adopt defined structures to perform biological activities, and conformational transitions among alternative structures are critical to virtually all RNA-mediated processes ranging from metabolite-activation of bacterial riboswitches to pre-mRNA splicing and viral replication in eukaryotes. Mechanistic analysis of an RNA folding reaction in a biological context is challenging because many steps usually intervene between assembly of a functional RNA structure and execution of a biological function. We developed a system to probe mechanisms of secondary structure folding and exchange directly in vivo using self-cleavage to monitor competition between mutually exclusive structures that promote or inhibit ribozyme assembly. In previous work, upstream structures were more effective than downstream structures in blocking ribozyme assembly during transcription in vitro, consistent with a sequential folding mechanism. However, upstream and downstream structures blocked ribozyme assembly equally well in vivo, suggesting that intracellular folding outcomes reflect thermodynamic equilibration or that annealing of contiguous sequences is favored kinetically. We have extended these studies to learn when, if ever, thermodynamic stability becomes an impediment to rapid equilibration among alternative RNA structures in vivo. We find that a narrow thermodynamic threshold determines whether kinetics or thermodynamics govern RNA folding outcomes in vivo. mRNA secondary structures fold sequentially in vivo, but exchange between adjacent secondary structures is much faster in vivo than it is in vitro. Previous work showed that simple base-paired RNA helices dissociate at similar rates in vivo and in vitro so exchange between adjacent structures must occur through a different mechanism, one that likely involves facilitation of branch migration by proteins associated with nascent transcripts.

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.

RNA folding during transcription: protocols and studies

RNA folds during transcription in the cell. Compared to most in vitro studies where the focus is generally on Mg(2+)-initiated refolding of fully synthesized transcripts, cotranscriptional RNA folding studies better replicate how RNA folds in a cellular environment. Unique aspects of cotranscriptional folding include the 5'- to 3'-polarity of RNA, the transcriptional speed, pausing properties of the RNA polymerase, the effect of the transcriptional complex and associated factors, and the effect of RNA-binding proteins. Identifying strategic pause sites can reveal insights on the folding pathway of the nascent transcript. Structural mapping of the paused transcription complexes identifies important folding intermediates along these pathways. Oligohybridization assays and the appearance of the catalytic activity of a ribozyme either in trans or in cis can be used to monitor cotranscriptional folding under a wide range of conditions. In our laboratory, these methodologies have been applied to study the folding of three highly conserved RNAs (RNase P, SRP, and tmRNA), several circularly permuted forms of a bacterial RNase P RNA, a riboswitch (thiM), and an aptamer-activated ribozyme (glmS).

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

RNA viruses are responsible for a variety of human diseases, and the pathogenicity of RNA viruses is often attributed to a high rate of mutation. Self-cleavage activity of the wild-type hepatitis delta virus (HDV) ribozyme as measured in standard divalent ion renaturation assays is biphasic and mostly slow and can be improved by multiple rational changes to ribozyme sequence or by addition of chemical denaturants. This is unusual in the sense that wild type is the most catalytically active sequence for the majority of protein enzymes, and RNA viruses are highly mutable. To see whether the ribozyme takes advantage of fast-reacting sequence changes in vivo, we performed alignment of 76 genomic and 269 antigenomic HDV isolates. Paradoxically, the sequence for the ribozyme was found to be essentially invariant in nature. We therefore tested whether three ribozyme sequence changes that improve self-cleavage under standard divalent ion renaturation assays also improve self-cleavage during transcription. Remarkably, wild type was as fast, or faster, than these mutants under cotranscriptional conditions. Slowing the rate of transcription or adding the hepatitis delta antigen protein only further stimulated cotranscriptional self-cleavage activity. Thus, the relative activity of HDV ribozyme mutants depends critically on whether the reaction is assayed under in vivo-like conditions. A model is presented for how wild-type ribozyme sequence and flanking sequence work in concert to promote efficient self-cleavage during transcription. Wild type being the optimal ribozyme sequence under in vivo-like conditions parallels the behavior of most protein enzymes.

Self-splicing of a group I intron reveals partitioning of native and misfolded RNA populations in yeast

Stable RNAs must form specific three-dimensional structures, yet many RNAs become kinetically trapped in misfolded conformations. To understand the factors that control the accuracy of RNA folding in the cell, the self-splicing activity of the Tetrahymena group I intron was compared in different genetic contexts in budding yeast. The extent of splicing was 98% when the intron was placed in its natural rDNA context, but only 3% when the intron was expressed in an exogenous pre-mRNA. Further experiments showed that the probability of forming the active intron structure depends on local sequence context and transcription by Pol I. Pre-rRNAs decayed at similar rates, whether the intron was wild type or inactivated by an internal deletion, suggesting that most of the unreacted pre-rRNA is incompetent to splice. Northern blots and complementation assays showed that mutations that destabilize the intron tertiary structure inhibited self-splicing and processing of internal transcribed spacer 2. The data are consistent with partitioning of pre-rRNAs into active and inactive populations. The misfolded RNAs are sequestered and degraded without refolding to a significant extent. Thus, the initial fidelity of folding can dictate the intracellular fate of transcripts containing this group I intron.

The role of a metastable RNA secondary structure in hepatitis delta virus genotype III RNA editing

RNA editing plays a critical role in the life cycle of hepatitis delta virus (HDV). The host editing enzyme ADAR1 recognizes specific RNA secondary structure features around the amber/W site in the HDV antigenome and deaminates the amber/W adenosine. A previous report suggested that a branched secondary structure is necessary for editing in HDV genotype III. This branched structure, which is distinct from the characteristic unbranched rod structure required for HDV replication, was only partially characterized, and knowledge concerning its formation and stability was limited. Here, we examine the secondary structures, conformational dynamics, and amber/W site editing of HDV genotype III RNA using a miniaturized HDV genotype III RNA in vitro. Computational analysis of this RNA using the MPGAfold algorithm indicated that the RNA has a tendency to form both metastable and stable unbranched secondary structures. Moreover, native polyacrylamide gel electrophoresis demonstrated that this RNA forms both branched and unbranched rod structures when transcribed in vitro. As predicted, the branched structure is a metastable structure that converts readily to the unbranched rod structure. Only branched RNA was edited at the amber/W site by ADAR1 in vitro. The structural heterogeneity of HDV genotype III RNA is significant because not only are both conformations of the RNA functionally important for viral replication, but the ratio of the two forms could modulate editing by determining the amount of substrate RNA available for modification.

RNA folding during transcription

The evolution of RNA sequence needs to satisfy three requirements: folding, structure, and function. Studies on folding during transcription are related directly to folding in the cell. Understanding RNA folding during transcription requires the elucidation of structure formation and structural changes of the RNA, and the consideration of intrinsic properties of the RNA polymerase and other proteins that interact with the RNA. This review summarizes the research progress in this area and outlines the enormous challenges facing this field. Significant advancement requires the development of new experimental methods and theoretical considerations in all aspects of transcription and RNA folding.

Computer simulation of chaperone effects of Archaeal C/D box sRNA binding on rRNA folding

Archaeal C/D box small RNAs (sRNAs) are homologues of eukaryotic C/D box small nucleolar RNAs (snoRNAs). Their main function is guiding 2'-O-ribose methylation of nucleotides in rRNAs. The methylation requires the pairing of an sRNA antisense element to an rRNA target site with formation of an RNA-RNA duplex. The temporary formation of such a duplex during rRNA maturation is expected to influence rRNA folding in a chaperone-like way, in particular in thermophilic Archaea, where multiple sRNAs with two binding sites are found. Here we investigate possible mechanisms of chaperone function of Archaeoglobus fulgidus and Pyrococcus abyssi C/D box sRNAs using computer simulations of rRNA secondary structure formation by genetic algorithm. The effects of sRNA binding on rRNA structure are introduced as temporary structural constraints during co-transcriptional folding. Comparisons of the final predictions with simulations without sRNA binding and with phylogenetic structures show that sRNAs with two antisense elements may significantly facilitate the correct formation of long-range interactions in rRNAs, in particular at elevated temperatures. The simulations suggest that the main mechanism of this effect is a transient restriction of folding in rRNA domains where the termini are brought together by binding to double-guide sRNAs.

Kinetics and thermodynamics make different contributions to RNA folding in vitro and in yeast

RNAs somehow adopt specific functional structures despite the capacity to form alternative nonfunctional structures with similar stabilities. We analyzed RNA assembly during transcription in vitro and in yeast using hairpin ribozyme self-cleavage to assess partitioning between functional ribozyme structures and nonfunctional stem loops. Complementary insertions located upstream of the ribozyme inhibited ribozyme assembly more than downstream insertions during transcription in vitro, consistent with a sequential folding model in which the outcome is determined by the structure that forms first. In contrast, both upstream and downstream insertions strongly inhibited assembly of the same ribozyme variants when expressed as chimeric mRNAs in yeast, indicating that inhibitory stem loops can form even after the entire ribozyme sequence has been transcribed. Evidently, some feature unique to the intracellular environment modulates the influence of transcription polarity and enhances the contribution of thermodynamic stability to RNA folding in vivo.

Linkage between proton binding and folding in RNA: implications for RNA catalysis

Small ribozymes use their nucleobases to catalyse phosphodiester bond cleavage. The hepatitis delta virus ribozyme employs C75 as a general acid to protonate the 5′-bridging oxygen leaving group, and to accomplish this task efficiently, it shifts its pKa towards neutrality. Simulations and thermodynamic experiments implicate linkage between folding and protonation in nucleobase pKa shifting. Even small oligonucleotides are shown to fold in a highly co-operative manner, although they do so in a context-specific fashion. Linkage between protonation and co-operativity of folding may drive pKa shifting and provide for enhanced function in RNA.

Mechanistic Insights on the Folding of a Large Ribozyme during Transcription†

RNA folding during transcription resembles folding in a cellular environment. We previously investigated the folding of a large ribozyme derived from a bacterial RNase P RNA during its transcription by the Escherichia coli RNA polymerase and the effect of the elongation factor NusA. We found that transcriptional pausing at a specific site induced by NusA significantly altered the folding pathway. In this work, we compare folding during transcription by the E. coli RNA polymerase of circularly permuted variants and site-specific mutants of the RNase P ribozyme to elucidate the molecular mechanism of transcriptional pausing and RNA folding. The effect of NusA-induced pausing depends on the order of RNA synthesis and only affects local folding of the RNA. Pausing likely prevents a misfolded structure between the 5' strand of a helix and its adjacent junction located in the specificity domain and a region known to bind single-stranded RNA located in the catalytic domain. These results lead to a structural model on how transcriptional pausing affects folding of RNase P RNA. Structural rearrangements of a nascent RNA transcript enhanced by transcriptional pausing may be a general feature of RNA folding during transcription.

A protein folding pathway with multiple folding intermediates at atomic resolution

Using native-state hydrogen-exchange-directed protein engineering and multidimensional NMR, we determined the high-resolution structure (rms deviation, 1.1 angstroms) for an intermediate of the four-helix bundle protein: Rd-apocytochrome b562. The intermediate has the N-terminal helix and a part of the C-terminal helix unfolded. In earlier studies, we also solved the structures of two other folding intermediates for the same protein: one with the N-terminal helix alone unfolded and the other with a reorganized hydrophobic core. Together, these structures provide a description of a protein folding pathway with multiple intermediates at atomic resolution. The two general features for the intermediates are (i) native-like backbone topology and (ii) nonnative side-chain interactions. These results have implications for important issues in protein folding studies, including large-scale conformation search, -value analysis, and computer simulations.

Folding thermodynamics and kinetics of YNMG RNA hairpins: specific incorporation of 8-bromoguanosine leads to stabilization by enhancement of the folding rate

Modified nucleotides allow fundamental energetic and kinetic properties of nucleic acids to be probed. Here, we demonstrate that an RNA hairpin containing the nucleotide analogue 8-bromoguanosine (8BrG or G), gcUUCGgc, has enhanced stability relative to the unmodified hairpin, with DeltaDeltaG(37)(degrees)= -0.69 +/- 0.15 kcal mol(-1) and DeltaT(M) = +6.8 +/- 1.4 degrees C. NMR spectroscopic data suggest that the enhanced stability of gcUUCGgc does not arise from the native state; laser temperature-jump experiments support this notion, as gcUUCGgc and gcUUCGgc have similar unfolding rate constants, but the folding rate constant of gcUUCGgc is 4.1-fold faster at 37.5 degrees C and 2.8-fold faster under isoenergetic conditions. On the basis of these findings, we propose that 8BrG reduces the conformational entropy of the denatured state, resulting in an accelerated conformational search for the native state and enhanced stability.

Design of a highly reactive HDV ribozyme sequence uncovers facilitation of RNA folding by alternative pairings and physiological ionic strength

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA that resides in the HDV genome and regulates its replication. The native fold of the ribozyme is complex, having two pseudoknots. Earlier work implicated four non-native pairings in slowing pseudoknot formation: Alt 1, Alt 2, Alt 3, and Alt P1. The goal of the present work was design of a kinetically simplified and maximally reactive construct for in vitro mechanistic and structural studies. The initial approach chosen was site-directed mutagenesis in which known alternative pairings were destabilized while leaving the catalytic core intact. Based on prior studies, the G11C/U27Delta double mutant was prepared. However, biphasic kinetics and antisense oligonucleotide response trends opposite those of the well-studied G11C mutant were observed suggesting that new alternative pairings with multiple registers, termed Alt X and Alt Y, had been created. Enzymatic structure mapping of oligonucleotide models supported this notion. This led to a model wherein Alt 2 and the phylogenetically conserved Alt 3 act as "folding guides", facilitating folding of the major population of the RNA molecules by hindering formation of the Alt X and Alt Y registers. Attempts to eliminate the strongest of the Alt X pairings by rational design of a quadruple mutant only resulted in more complex kinetic behavior. In an effort to simultaneously destabilize multiple alternative pairings, studies were carried out on G11C/U27Delta in the presence of urea or increased monovalent ion concentration. Inclusion of physiological ionic strength allowed the goal of monophasic, fast-folding (kobs approximately 60 min(-1)) kinetics to be realized. To account for this, a model is developed wherein Na+, which destabilizes secondary and tertiary structures in the presence of Mg2+, facilitates native folding by destabilizing the multiple alternative secondary structures with a higher-order dependence.

RNA folding: models and perspectives

Intrinsic events during RNA folding include conformational search and metal ion binding. Several experimentally testable models have been proposed to explain how large ribozymes accomplish folding. Future challenges include the validation of these models, and the correlation of experimental results and theoretical simulations.

Effect of transcription on folding of the Tetrahymena ribozyme

Sequential formation of RNA interactions during transcription can bias the folding pathway and ultimately determine the functional state of a transcript. The kinetics of cotranscriptional folding of the Tetrahymena L-21 ribozyme was compared with refolding of full-length transcripts under the same conditions. Sequential folding after transcription by phage T7 or Escherichia coli polymerase is only twice as fast as refolding, and the yield of native RNA is the same. By contrast, a greater fraction of circularly permuted variants folded correctly at early times during transcription than during refolding. Hybridization of complementary oligonucleotides suggests that cotranscriptional folding enables a permuted RNA beginning at G303 to escape non-native interactions in P3 and P9. We propose that base pairing of upstream sequences during transcription elongation favors branched secondary structures that increase the probability of forming the native ribozyme structure.

Mechanistic characterization of the HDV genomic ribozyme: classifying the catalytic and structural metal ion sites within a multichannel reaction mechanism

Prior studies of the metal ion dependence of the self-cleavage reaction of the HDV genomic ribozyme led to a mechanistic framework in which the ribozyme can self-cleave by multiple Mg2+ ion-independent and -dependent channels [Nakano et al. (2001) Biochemistry 40, 12022]. In particular, channel 2 involves cleavage in the presence of a structural Mg2+ ion without participation of a catalytic divalent metal ion, while channel 3 involves both structural and catalytic Mg2+ ions. In the present study, experiments were performed to probe the nature of the various divalent ion sites and any specificity for Mg2+. A series of alkaline earth metal ions was tested for the ability to catalyze self-cleavage of the ribozyme under conditions that favor either channel 2 or channel 3. Under conditions that populate primarily channel 3, nearly identical K(d)s were obtained for Mg2+, Ca2+, Ba2+, and Sr2+, with a slight discrimination against Ca2+. In contrast, under conditions that populate primarily channel 2, tighter binding was observed as ion size decreases. Moreover, [Co(NH3)6]3+ was found to be a strong competitive inhibitor of Mg2+ for channel 3 but not for channel 2. The thermal unfolding of the cleaved ribozyme was also examined, and two transitions were found. Urea-dependent studies gave m-values that allowed the lower temperature transition to be assigned to tertiary structure unfolding. The effects of high concentrations of Na+ on the melting temperature for RNA unfolding and the reaction rate revealed ion binding to the folded RNA, with significant competition of Na+ (Hill coefficient of 1.5-1.7) for a structural Mg2+ ion and an unusually high intrinsic affinity of the structural ion for the RNA. Taken together, these data support the existence of two different classes of metal ion sites on the ribozyme: a structural site that is inner sphere with a major electrostatic component and a preference for Mg2+, and a weak catalytic site that is outer sphere with little preference for a particular divalent ion.