The folding pathway of the genomic hepatitis delta virus ribozyme is dominated by slow folding of the pseudoknots

Optimized periphery-core interface increases fitness of the Bacillus subtilisglmS ribozyme

Like other functional RNAs, ribozymes encode a conserved catalytic center supported by peripheral domains that vary among ribozyme sub-families. To understand how core-periphery interactions contribute to ribozyme fitness, we compared the cleavage kinetics of all single base substitutions at 152 sites across the Bacillus subtilis glmS ribozyme by high-throughput sequencing (k-seq). The in vitro activity map mirrored phylogenetic sequence conservation in glmS ribozymes, indicating that biological fitness reports all biochemically important positions. The k-seq results and folding assays showed that most deleterious mutations lower activity by impairing ribozyme self-assembly. All-atom molecular dynamics simulations of the complete ribozyme revealed how individual mutations in the core or the IL4 peripheral loop introduce a non-native tertiary interface that rewires the catalytic center, eliminating activity. We conclude that the need to avoid non-native helix packing powerfully constrains the evolution of tertiary structure motifs in RNA.

The Functional RNA Identification (FRID) Pipeline: Identification of Potential Pseudoknot-Containing RNA Elements as Therapeutic Targets for SARS-CoV-2

The COVID-19 pandemic persists despite the development of effective vaccines. As such, it remains crucial to identify new targets for antiviral therapies. The causative virus of COVID-19, SARS-CoV-2, is a positive-sense RNA virus with RNA structures that could serve as therapeutic targets. One such RNA with established function is the frameshift stimulatory element (FSE), which promotes programmed ribosomal frameshifting. To accelerate identification of additional functional RNA elements, we introduce a novel computational approach termed the Functional RNA Identification (FRID) pipeline. The guiding principle of our pipeline, which uses established component programs as well as customized component programs, is that functional RNA elements have conserved secondary and pseudoknot structures that facilitate function. To assess the presence and conservation of putative functional RNA elements in SARS-CoV-2, we compared over 6,000 SARS-CoV-2 genomic isolates. We identified 22 functional RNA elements from the SARS-CoV-2 genome, 14 of which have conserved pseudoknots and serve as potential targets for small molecule or antisense oligonucleotide therapeutics. The FRID pipeline is general and can be applied to identify pseudoknotted RNAs for targeted therapeutics in genomes or transcriptomes from any virus or organism.

Regulation of the ThiM riboswitch is facilitated by the trapped structure formed during transcription of the wild-type sequence

The ThiM riboswitch from Escherichia coli is a typical mRNA device that modulates downstream gene expression by sensing TPP. The helix-based RNA folding theory is used to investigate its detailed regulatory behaviors in cells. This RNA molecule is transcriptionally trapped in a state with the unstructured SD sequence in the absence of TPP, which induces downstream gene expression. As a key step to turn on gene expression, formation of this trapped state (the genetic ON state) highly depends on the co-transcriptional folding of its wild-type sequence. Instead of stabilities of the genetic ON and OFF states, the transcription rate, pause, and ligand levels are combined to affect the ThiM riboswitch-mediated gene regulation, which is consistent with a kinetic control model.

Rhinovirus 2A is the key protease responsible for instigating the early block to gene expression in infected cells

Human rhinoviruses (HRVs) express 2 cysteine proteases, 2A and 3C, that are responsible for viral polyprotein processing. Both proteases also suppress host gene expression by inhibiting mRNA transcription, nuclear export and cap-dependent translation. However, the relative contribution that each makes in achieving this goal remains unclear. In this study, we have compared both the combined and individual ability of the two proteases to shut down cellular gene expression using a novel dynamic reporter system. Our findings show that 2A inhibits host gene expression much more rapidly than 3C. By comparing the activities of a representative set of proteases from the three different HRV species, we also find variation in the speed at which host gene expression is suppressed. Our work highlights the key role that 2A plays in early suppression of the infected host cell response and shows that this can be influenced by natural variation in the activity of this enzyme.

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.

Predicting 3D structure and stability of RNA pseudoknots in monovalent and divalent ion solutions

RNA pseudoknots are a kind of minimal RNA tertiary structural motifs, and their three-dimensional (3D) structures and stability play essential roles in a variety of biological functions. Therefore, to predict 3D structures and stability of RNA pseudoknots is essential for understanding their functions. In the work, we employed our previously developed coarse-grained model with implicit salt to make extensive predictions and comprehensive analyses on the 3D structures and stability for RNA pseudoknots in monovalent/divalent ion solutions. The comparisons with available experimental data show that our model can successfully predict the 3D structures of RNA pseudoknots from their sequences, and can also make reliable predictions for the stability of RNA pseudoknots with different lengths and sequences over a wide range of monovalent/divalent ion concentrations. Furthermore, we made comprehensive analyses on the unfolding pathway for various RNA pseudoknots in ion solutions. Our analyses for extensive pseudokonts and the wide range of monovalent/divalent ion concentrations verify that the unfolding pathway of RNA pseudoknots is mainly dependent on the relative stability of unfolded intermediate states, and show that the unfolding pathway of RNA pseudoknots can be significantly modulated by their sequences and solution ion conditions.

Ensemble simulations: folding, unfolding and misfolding of a high-efficiency frameshifting RNA pseudoknot

Massive all-atom molecular dynamics simulations were conducted across a distributed computing network to study the folding, unfolding, misfolding and conformational plasticity of the high-efficiency frameshifting double mutant of the 26 nt potato leaf roll virus RNA pseudoknot. Our robust sampling, which included over 40 starting structures spanning the spectrum from the extended unfolded state to the native fold, yielded nearly 120 μs of cumulative sampling time. Conformational microstate transitions on the 1.0 ns to 10.0 μs timescales were observed, with post-equilibration sampling providing detailed representations of the conformational free energy landscape and the complex folding mechanism inherent to the pseudoknot motif. Herein, we identify and characterize two alternative native structures, three intermediate states, and numerous misfolded states, the latter of which have not previously been characterized via atomistic simulation techniques. While in line with previous thermodynamics-based models of a general RNA folding mechanism, our observations indicate that stem-strand-sequence-separation may serve as an alternative predictor of the order of stem formation during pseudoknot folding. Our results contradict a model of frameshifting based on structural rigidity and resistance to mechanical unfolding, and instead strongly support more recent studies in which conformational plasticity is identified as a determining factor in frameshifting efficiency.

Cooperative RNA Folding under Cellular Conditions Arises From Both Tertiary Structure Stabilization and Secondary Structure Destabilization

RNA folding has been studied extensively in vitro, typically under dilute solution conditions and abiologically high salt concentrations of 1 M Na⁺ or 10 mM Mg²⁺. The cellular environment is very different, with 20-40% crowding and only 10-40 mM Na⁺, 140 mM K⁺, and 0.5-2.0 mM Mg²⁺. As such, RNA structures and functions can be radically altered under cellular conditions. We previously reported that tRNA^phe secondary and tertiary structures unfold together in a cooperative two-state fashion under crowded in vivo-like ionic conditions, but in a noncooperative multistate fashion under dilute in vitro ionic conditions unless in nonphysiologically high concentrations of Mg²⁺. The mechanistic basis behind these effects remains unclear, however. To address the mechanism that drives RNA folding cooperativity, we probe effects of cellular conditions on structures and stabilities of individual secondary structure fragments comprising the full-length RNA. We elucidate effects of a diverse set of crowders on tRNA secondary structural fragments and full-length tRNA at three levels: at the nucleotide level by temperature-dependent in-line probing, at the tertiary structure level by small-angle X-ray scattering, and at the global level by thermal denaturation. We conclude that cooperative RNA folding is induced by two overlapping mechanisms: increased stability and compaction of tertiary structure through effects of Mg²⁺, and decreased stability of certain secondary structure elements through the effects of molecular crowders. These findings reveal that despite having very different chemical makeups RNA and protein can both have weak secondary structures in vivo leading to cooperative folding.

The GlcN6P cofactor plays multiple catalytic roles in the glmS ribozyme

RNA enzymes (ribozymes) have remarkably diverse biological roles despite having limited chemical diversity. Protein enzymes enhance their reactivity through recruitment of cofactors; likewise, the naturally occurring glmS ribozyme uses the glucosamine-6-phosphate (GlcN6P) organic cofactor for phosphodiester bond cleavage. Prior structural and biochemical studies have implicated GlcN6P as the general acid. Here we describe new catalytic roles of GlcN6P through experiments and calculations. Large stereospecific normal thio effects and a lack of metal-ion rescue in the holoribozyme indicate that nucleobases and the cofactor play direct chemical roles and align the active site for self-cleavage. Large stereospecific inverse thio effects in the aporibozyme suggest that the GlcN6P cofactor disrupts an inhibitory interaction of the nucleophile. Strong metal-ion rescue in the aporibozyme reveals that this cofactor also provides electrostatic stabilization. Ribozyme organic cofactors thus perform myriad catalytic roles, thereby allowing RNA to compensate for its limited functional diversity.

Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway

The past decade has seen a wealth of 3D structural information about complex structured RNAs and identification of functional intermediates. Nevertheless, developing a complete and predictive understanding of the folding and function of these RNAs in biology will require connection of individual rate and equilibrium constants to structural changes that occur in individual folding steps and further relating these steps to the properties and behavior of isolated, simplified systems. To accomplish these goals we used the considerable structural knowledge of the folded, unfolded, and intermediate states of P4-P6 RNA. We enumerated structural states and possible folding transitions and determined rate and equilibrium constants for the transitions between these states using single-molecule FRET with a series of mutant P4-P6 variants. Comparisons with simplified constructs containing an isolated tertiary contact suggest that a given tertiary interaction has a stereotyped rate for breaking that may help identify structural transitions within complex RNAs and simplify the prediction of folding kinetics and thermodynamics for structured RNAs from their parts. The preferred folding pathway involves initial formation of the proximal tertiary contact. However, this preference was only ∼10 fold and could be reversed by a single point mutation, indicating that a model akin to a protein-folding contact order model will not suffice to describe RNA folding. Instead, our results suggest a strong analogy with a modified RNA diffusion-collision model in which tertiary elements within preformed secondary structures collide, with the success of these collisions dependent on whether the tertiary elements are in their rare binding-competent conformations.

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.

Four RNA families with functional transient structures

Protein-coding and non-coding RNA transcripts perform a wide variety of cellular functions in diverse organisms. Several of their functional roles are expressed and modulated via RNA structure. A given transcript, however, can have more than a single functional RNA structure throughout its life, a fact which has been previously overlooked. Transient RNA structures, for example, are only present during specific time intervals and cellular conditions. We here introduce four RNA families with transient RNA structures that play distinct and diverse functional roles. Moreover, we show that these transient RNA structures are structurally well-defined and evolutionarily conserved. Since Rfam annotates one structure for each family, there is either no annotation for these transient structures or no such family. Thus, our alignments either significantly update and extend the existing Rfam families or introduce a new RNA family to Rfam. For each of the four RNA families, we compile a multiple-sequence alignment based on experimentally verified transient and dominant (dominant in terms of either the thermodynamic stability and/or attention received so far) RNA secondary structures using a combination of automated search via covariance model and manual curation. The first alignment is the Trp operon leader which regulates the operon transcription in response to tryptophan abundance through alternative structures. The second alignment is the HDV ribozyme which we extend to the 5' flanking sequence. This flanking sequence is involved in the regulation of the transcript's self-cleavage activity. The third alignment is the 5' UTR of the maturation protein from Levivirus which contains a transient structure that temporarily postpones the formation of the final inhibitory structure to allow translation of maturation protein. The fourth and last alignment is the SAM riboswitch which regulates the downstream gene expression by assuming alternative structures upon binding of SAM. All transient and dominant structures are mapped to our new alignments introduced here.

Folding pathways of the Tetrahymena ribozyme

Like many structured RNAs, the Tetrahymena group I intron ribozyme folds through multiple pathways and intermediates. Under standard conditions in vitro, a small fraction reaches the native state (N) with kobs ≈ 0.6 min(-1), while the remainder forms a long-lived misfolded conformation (M) thought to differ in topology. These alternative outcomes reflect a pathway that branches late in folding, after disruption of a trapped intermediate (Itrap). Here we use catalytic activity to probe the folding transitions from Itrap to the native and misfolded states. We show that mutations predicted to weaken the core helix P3 do not increase the rate of folding from Itrap but they increase the fraction that reaches the native state rather than forming the misfolded state. Thus, P3 is disrupted during folding to the native state but not to the misfolded state, and P3 disruption occurs after the rate-limiting step. Interestingly, P3-strengthening mutants also increase native folding. Additional experiments show that these mutants are rapidly committed to folding to the native state, although they reach the native state with approximately the same rate constant as the wild-type ribozyme (~1 min(-1)). Thus, the P3-strengthening mutants populate a distinct pathway that includes at least one intermediate but avoids the M state, most likely because P3 and the correct topology are formed early. Our results highlight multiple pathways in RNA folding and illustrate how kinetic competitions between rapid events can have long-lasting effects because the "choice" is enforced by energy barriers that grow larger as folding progresses.

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.

Molecular crowders and cosolutes promote folding cooperativity of RNA under physiological ionic conditions

Folding mechanisms of functional RNAs under idealized in vitro conditions of dilute solution and high ionic strength have been well studied. Comparatively little is known, however, about mechanisms for folding of RNA in vivo where Mg(2+) ion concentrations are low, K(+) concentrations are modest, and concentrations of macromolecular crowders and low-molecular-weight cosolutes are high. Herein, we apply a combination of biophysical and structure mapping techniques to tRNA to elucidate thermodynamic and functional principles that govern RNA folding under in vivo-like conditions. We show by thermal denaturation and SHAPE studies that tRNA folding cooperativity increases in physiologically low concentrations of Mg(2+) (0.5-2 mM) and K(+) (140 mM) if the solution is supplemented with physiological amounts (∼ 20%) of a water-soluble neutral macromolecular crowding agent such as PEG or dextran. Low-molecular-weight cosolutes show varying effects on tRNA folding cooperativity, increasing or decreasing it based on the identity of the cosolute. For those additives that increase folding cooperativity, the gain is manifested in sharpened two-state-like folding transitions for full-length tRNA over its secondary structural elements. Temperature-dependent SHAPE experiments in the absence and presence of crowders and cosolutes reveal extent of cooperative folding of tRNA on a nucleotide basis and are consistent with the melting studies. Mechanistically, crowding agents appear to promote cooperativity by stabilizing tertiary structure, while those low molecular cosolutes that promote cooperativity stabilize tertiary structure and/or destabilize secondary structure. Cooperative folding of functional RNA under physiological-like conditions parallels the behavior of many proteins and has implications for cellular RNA folding kinetics and evolution.

Toward a molecular understanding of RNA remodeling by DEAD-box proteins

DEAD-box proteins are superfamily 2 helicases that function in all aspects of RNA metabolism. They employ ATP binding and hydrolysis to generate tight, yet regulated RNA binding, which is used to unwind short RNA helices non-processively and promote structural transitions of RNA and RNA-protein substrates. In the last few years, substantial progress has been made toward a detailed, quantitative understanding of the structural and biochemical properties of DEAD-box proteins. Concurrently, progress has been made toward a physical understanding of the RNA rearrangements and folding steps that are accelerated by DEAD-box proteins in model systems. Here, we review the recent progress on both of these fronts, focusing on the mitochondrial DEAD-box proteins Mss116 and CYT-19 and their mechanisms in promoting the splicing of group I and group II introns.

Kinetics of tRNA folding monitored by aminoacylation

We describe a strategy for tracking Mg²⁺-initiated folding of ³²P-labeled tRNA molecules to their native structures based on the capacity for aminoacylation by the cognate aminoacyl-tRNA synthetase enzyme. The approach directly links folding to function, paralleling a common strategy used to study the folding of catalytic RNAs. Incubation of unfolded tRNA with magnesium ions, followed by the addition of aminoacyl-tRNA synthetase and further incubation, yields a rapid burst of aminoacyl-tRNA formation corresponding to the prefolded tRNA fraction. A subsequent slower increase in product formation monitors continued folding in the presence of the enzyme. Further analysis reveals the presence of a parallel fraction of tRNA that folds more rapidly than the majority of the population. The application of the approach to study the influence of post-transcriptional modifications in folding of Escherichia coli tRNA₁(Gln) reveals that the modified bases increase the folding rate but do not affect either the equilibrium between properly folded and misfolded states or the folding pathway. This assay allows the use of ³²P-labeled tRNA in integrated studies combining folding, post-transcriptional processing, and aminoacylation reactions.

Fast folding of RNA pseudoknots initiated by laser temperature-jump

RNA pseudoknots are examples of minimal structural motifs in RNA with tertiary interactions that stabilize the structures of many ribozymes. They also play an essential role in a variety of biological functions that are modulated by their structure, stability, and dynamics. Therefore, understanding the global principles that determine the thermodynamics and folding pathways of RNA pseudoknots is an important problem in biology, both for elucidating the folding mechanisms of larger ribozymes as well as addressing issues of possible kinetic control of the biological functions of pseudoknots. We report on the folding/unfolding kinetics of a hairpin-type pseudoknot obtained with microsecond time-resolution in response to a laser temperature-jump perturbation. The kinetics are monitored using UV absorbance as well as fluorescence of extrinsically attached labels as spectroscopic probes of the transiently populated RNA conformations. We measure folding times of 1-6 ms at 37 °C, which are at least 100-fold faster than previous observations of very slow folding pseudoknots that were trapped in misfolded conformations. The measured relaxation times are remarkably similar to predictions of a computational study by Thirumalai and co-workers (Cho, S. S.; Pincus, D.L.; Thirumalai, D. Proc. Natl. Acad. Sci. U. S. A. 2009, 106, 17349-17354). Thus, these studies provide the first observation of a fast-folding pseudoknot and present a benchmark against which computational models can be refined.

Prevalence of syn nucleobases in the active sites of functional RNAs

Biological RNAs, like their DNA counterparts, contain helical stretches, which have standard Watson-Crick base pairs in the anti conformation. Most functional RNAs also adopt geometries with far greater complexity such as bulges, loops, and multihelical junctions. Occasionally, nucleobases in these regions populate the syn conformation wherein the base resides close to or over the ribose sugar, which leads to a more compact state. The importance of the syn conformation to RNA function is largely unknown. In this study, we analyze 51 RNAs with tertiary structure, including aptamers, riboswitches, ribozymes, and ribosomal RNAs, for number, location, and properties of syn nucleobases. These RNAs represent the set of nonoverlapping, moderate- to high-resolution structures available at present. We find that syn nucleobases are much more common among purines than pyrimidines, and that they favor C2'-endo-like conformations especially among those nucleobases in the intermediate syn conformation. Strikingly, most syn nucleobases participate in tertiary stacking and base-pairing interactions: Inspection of RNA structures revealed that the majority of the syn nucleobases are in regions assigned to function, with many syn nucleobases interacting directly with a ligand or ribozyme active site. These observations suggest that judicious placement of conformationally restricted nucleotides biased into the syn conformation could enhance RNA folding and catalysis. Such changes could also be useful for locking RNAs into functionally competent folds for use in X-ray crystallography and NMR.

HDV-like self-cleaving ribozymes

HDV ribozymes catalyze their own scission from the transcript during rolling circle replication of the hepatitis delta virus. In vitro selection of self-cleaving ribozymes from a human genomic library revealed an HDV-like ribozyme in the second intron of the human CPEB3 gene and recent results suggest that this RNA affects episodic memory performance. Bioinformatic searches based on the secondary structure of the HDV/CPEB3 fold yielded numerous functional ribozymes in a wide variety of organisms. Genomic mapping of these RNAs suggested several biological roles, one of which is the 5' processing of non-LTR retrotransposons. The family of HDV-like ribozymes thus continues to grow in numbers and biological importance.

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.

Probing mechanism and transition state of RNA refolding

Kinetics and the atomic detail of RNA refolding are only poorly understood. It has been proposed that conformations with transient base pairing interaction are populated during RNA refolding, but a detailed description of those states is lacking. By NMR and CD spectroscopy, we examined the refolding of a bistable RNA and the influence of urea, Mg(2+), and spermidine on its refolding kinetics. The bistable RNA serves as a model system and exhibits two almost equally stable ground-state conformations. We designed a photolabile caged RNA to selectively stabilize one of the two ground-state conformations and trigger RNA refolding by in situ light irradiation in the NMR spectrometer. We can show that the refolding kinetics of the bistable RNA is modulated by urea, Mg(2+), and spermidine by different mechanisms. From a statistical analysis based on elementary rate constants, we deduce the required number of base pairs that need to be destabilized during the refolding transition and propose a model for the transition state of the folding reaction.

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.

Functional analysis of the SRV-1 RNA frameshifting pseudoknot

Simian retrovirus type-1 uses programmed ribosomal frameshifting to control expression of the Gag-Pol polyprotein from overlapping gag and pol open-reading frames. The frameshifting signal consists of a heptanucleotide slippery sequence and a downstream-located 12-base pair pseudoknot. The solution structure of this pseudoknot, previously solved by NMR [Michiels,P.J., Versleijen,A.A., Verlaan,P.W., Pleij,C.W., Hilbers,C.W. and Heus,H.A. (2001) Solution structure of the pseudoknot of SRV-1 RNA, involved in ribosomal frameshifting. J. Mol. Biol., 310, 1109-1123] has a classical H-type fold and forms an extended triple helix by interactions between loop 2 and the minor groove of stem 1 involving base-base and base-sugar contacts. A mutational analysis was performed to test the functional importance of the triple helix for -1 frameshifting in vitro. Changing bases in L2 or base pairs in S1 involved in a base triple resulted in a 2- to 5-fold decrease in frameshifting efficiency. Alterations in the length of L2 had adverse effects on frameshifting. The in vitro effects were well reproduced in vivo, although the effect of enlarging L2 was more dramatic in vivo. The putative role of refolding kinetics of frameshifter pseudoknots is discussed. Overall, the data emphasize the role of the triple helix in -1 frameshifting.

Catalytic activity as a probe of native RNA folding

As RNAs fold to functional structures, they traverse complex energy landscapes that include many partially folded and misfolded intermediates. For structured RNAs that possess catalytic activity, this activity can provide a powerful means of monitoring folding that is complementary to biophysical approaches. RNA catalysis can be used to track accumulation of the native RNA specifically and quantitatively, readily distinguishing the native structure from intermediates that resemble it and may not be differentiated by other approaches. Here, we outline how to design and interpret experiments using catalytic activity to monitor RNA folding, and we summarize adaptations of the method that have been used to probe aspects of folding well beyond determination of the folding rates.

Slow formation of a pseudoknot structure is rate limiting in the productive co-transcriptional folding of the self-splicing Candida intron

Pseudoknots play critical roles in packing the active structure of various functional RNAs. The importance of the P3-P7 pseudoknot in refolding of group I intron ribozymes has been recently appreciated, while little is known about the pseudoknot function in co-transcriptional folding. Here we used the Candida group I intron as a model to address the question. We show that co-transcriptional folding of the active self-splicing intron is twice as fast as refolding. The P3-P7 pseudoknot folds slowly during co-transcriptional folding at a rate constant similar to the folding of the active ribozyme, and folding of both P3-P7 and P1-P10 pseudoknots are inhibited by antisense oligonucleotides. We conclude that when RNA folding is coupled with transcription, formation of pseudoknot structures dominates the productive folding pathway and serves as a rate-limiting step in producing the self-splicing competent Candida intron.

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.

Thinking inside the box: designing, implementing, and interpreting thermodynamic cycles to dissect cooperativity in RNA and DNA folding

Double and triple mutant thermodynamic cycles provide a means to dissect the cooperativity of RNA and DNA folding at both the secondary and tertiary structural levels through use of the thermodynamic box or cube. In this article, we describe three steps for applying thermodynamic cycles to nucleic acid folding, with considerations of both conceptual and experimental features. The first step is design of an appropriate system and development of hypotheses regarding which residues might interact. Next is implementing this design in terms of a tractable experimental strategy, with an emphasis on UV melting. The final step, and the one we emphasize the most, is interpreting mutant cycles in terms of coupling between specific residues in the RNA or DNA. Coupling free energy in the absence and presence of changes elsewhere in the molecule is discussed in terms of specific folding models, including stepwise folding and concerted changes. Last, we provide a practical section on the use of commercially available software (KaleidaGraph) to fit melting data, along with a consideration of error propagation. Along the way, specific examples are chosen from the literature to illustrate the methods. This article is intended to be accessible to the biochemist or biologist without extensive thermodynamics background.

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: solvent isotope effects and proton inventories in the absence of divalent metal ions support C75 as the general acid

The hepatitis delta virus (HDV) ribozyme uses the nucleobase C75 and a hydrated Mg(2+) ion as the general acid-base catalysts in phosphodiester bond cleavage at physiological salt. A mechanistic framework has been advanced that involves one Mg(2+)-independent and two Mg(2+)-dependent channels. The rate-pH profile for wild-type (WT) ribozyme in the Mg(2+)-free channel is inverted relative to the fully Mg(2+)-dependent channel, with each having a near-neutral pKa. Inversion of the rate-pH profile was used as the crux of a mechanistic argument that C75 serves as general acid both in the presence and absence of Mg(2+). However, subsequent studies on a double mutant (DM) ribozyme suggested that the pKa observed for WT in the absence of Mg(2+) arises from ionization of C41, a structural nucleobase. To investigate this further, we acquired rate-pH/pD profiles and proton inventories for WT and DM in the absence of Mg(2+). Corrections were made for effects of ionic strength on hydrogen ion activity and pH meter readings. Results are accommodated by a model wherein the Mg(2+)-free pKa observed for WT arises from ionization of C75, and DM reactivity is compromised by protonation of C41. The Brønsted base appears to be water or hydroxide ion depending on pH. The observed pKa's are related to salt-dependent pH titrations of a model oligonucleotide, as well as electrostatic calculations, which support the local environment for C75 in the absence of Mg(2+) being similar to that in the presence of Mg(2+) and impervious to bulk ions. Accordingly, the catalytic role of C75 as the general acid does not appear to depend on divalent ions or the identity of the Brønsted base.

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.

The poly(A) site sequence in HDV RNA alters both extent and rate of self-cleavage of the antigenomic ribozyme

The ribozyme self-cleavage site in the antigenomic sequence of hepatitis delta virus (HDV) RNA is 33-nt downstream of the poly(A) site for the delta antigen mRNA. An HDV antigenomic ribozyme precursor RNA that included the upstream poly(A) processing site was used to test the hypothesis that nonribozyme sequence near the poly(A) site could affect ribozyme activity. Relative to ribozyme precursor without the extra upstream sequences, the kinetic profile for self-cleavage of the longer precursor was altered in two ways. First, only half of the precursor RNA self-cleaved. The cleaved fraction could be increased or decreased with mutations in the upstream sequence. These mutations, which were predicted to alter the relative stability of competing secondary structures within the precursor, changed the distribution of alternative RNA structures that are resolved in native-gel electrophoresis. Second, the active fraction cleaved with an observed rate constant that was higher than that of the ribozyme without the upstream sequences. Moreover, the higher rate constants occurred at lower, near-physiological, divalent metal ion concentrations (1–2 mM). Modulation of ribozyme activity, through competing alternative structures, could be part of a mechanism that allows a biologically significant choice between maturation of the mRNA and processing of replication intermediates.

DMS footprinting of structured RNAs and RNA-protein complexes

We describe a protocol in which dimethyl sulfate (DMS) modification of the base-pairing faces of unpaired adenosine and cytidine nucleotides is used for structural analysis of RNAs and RNA-protein complexes (RNPs). The protocol is optimized for RNAs of small to moderate size (< or = 500 nt). The RNA or RNP is first exposed to DMS under conditions that promote formation of the folded structure or complex, as well as 'control' conditions that do not allow folding or complex formation. The positions and extents of modification are then determined by primer extension, polyacrylamide gel electrophoresis and quantitative analysis. From changes in the extent of modification upon folding or protein binding (appearance of a 'footprint'), it is possible to detect local changes in the secondary and tertiary structure of RNA, as well as the formation of RNA-protein contacts. This protocol takes 1.5-3 d to complete, depending on the type of analysis used.

RNA misfolding and the action of chaperones

RNA folds to a myriad of three-dimensional structures and performs an equally diverse set of functions. The ability of RNA to fold and function in vivo is all the more remarkable because, in vitro, RNA has been shown to have a strong propensity to adopt misfolded, non-functional conformations. A principal factor underlying the dominance of RNA misfolding is that local RNA structure can be quite stable even in the absence of enforcing global tertiary structure. This property allows non-native structure to persist, and it also allows native structure to form and stabilize non-native contacts or non-native topology. In recent years it has become clear that one of the central reasons for the apparent disconnect between the capabilities of RNA in vivo and its in vitro folding properties is the presence of RNA chaperones, which facilitate conformational transitions of RNA and therefore mitigate the deleterious effects of RNA misfolding. Over the past two decades, it has been demonstrated that several classes of non-specific RNA binding proteins possess profound RNA chaperone activity in vitro and when overexpressed in vivo, and at least some of these proteins appear to function as chaperones in vivo. More recently, it has been shown that certain DExD/H-box proteins function as general chaperones to facilitate folding of group I and group II introns. These proteins are RNA-dependent ATPases and have RNA helicase activity, and are proposed to function by using energy from ATP binding and hydrolysis to disrupt RNA structure and/or to displace proteins from RNA-protein complexes. This review outlines experimental studies that have led to our current understanding of the range of misfolded RNA structures, the physical origins of RNA misfolding, and the functions and mechanisms of putative RNA chaperone proteins.

Single-molecule mechanical unfolding and folding of a pseudoknot in human telomerase RNA

RNA unfolding and folding reactions in physiological conditions can be facilitated by mechanical force one molecule at a time. By using force-measuring optical tweezers, we studied the mechanical unfolding and folding of a hairpin-type pseudoknot in human telomerase RNA in a near-physiological solution, and at room temperature. Discrete two-state folding transitions of the pseudoknot are seen at approximately 10 and approximately 5 piconewtons (pN), with ensemble rate constants of approximately 0.1 sec(-1), by stepwise force-drop experiments. Folding studies of the isolated 5'-hairpin construct suggested that the 5'-hairpin within the pseudoknot forms first, followed by formation of the 3'-stem. Stepwise formation of the pseudoknot structure at low forces are in contrast with the one-step unfolding at high forces of approximately 46 pN, at an average rate of approximately 0.05 sec(-1). In the constant-force folding trajectories at approximately 10 pN and approximately 5 pN, transient formation of nonnative structures were observed, which is direct experimental evidence that folding of both the hairpin and pseudoknot takes complex pathways. Possible nonnative structures and folding pathways are discussed.

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.

Direct measurement of a pK(a) near neutrality for the catalytic cytosine in the genomic HDV ribozyme using Raman crystallography

The hepatitis delta virus (HDV) ribozyme uses a cytosine to facilitate general acid-base catalysis. Biochemical studies suggest that C75 has a pKa perturbed to near neutrality. To measure this pKa directly, Raman spectra were recorded on single ribozyme crystals using a Raman microscope. A spectral feature arising from a single neutral cytosine was identified at 1528 cm(-1). At low pH, this mode was replaced with a new spectral feature. Monitoring these features as a function of pH revealed pKa values for the cytosine that couple anticooperatively with Mg2+ binding, with values of 6.15 and 6.40 in the presence of 20 and 2 mM Mg2+, respectively. These pKa values agree well with those obtained from ribozyme activity experiments in solution. To correlate the observed pKa with a specific nucleotide, crystals of C75U, which is catalytically inactive, were examined. The Raman difference spectra show that this mutation does not affect the conformation of the ribozyme. However, crystals of C75U did not produce a signal from a protonatable cytosine, providing strong evidence that protonation of C75 is being monitored in the wild-type ribozyme. These studies provide the first direct physical measurement of a pKa near neutrality for a catalytic residue in a ribozyme and show that ribozymes, like their protein enzyme counterparts, can optimize the pKa of their side chains for proton transfer.

Mechanistic characterization of the HDV genomic ribozyme: a mutant of the C41 motif provides insight into the positioning and thermodynamic linkage of metal ions and protons

Binding of two Mg2+ and two H+ ions influences the self-cleavage activity of the genomic HDV ribozyme. The positioning of these four ligands and their thermodynamic linkage are not fully resolved. Protonated C41 engages in a base triple, whereas protonated C75 has been implicated as an acid-base catalyst in bond cleavage. Prior studies led to the identification of one structural inner-sphere ion and one catalytic outer-sphere ion. In the present study, the contributions of the C41 base triple to the metal ion- and pH-dependence of the reaction are examined. Experiments were conducted on a CG to UA double mutant (DM), which changes the base triple to one involving an unprotonated C41. Below pH 6, the DM has a steeper dependence on pH than the wild-type (WT), consistent with a single protonation misfolding the core; this conclusion is also supported by thermal denaturation studies. Between pH 6 and 8, the WT and DM display nearly identical catalytic metal ion and H+ binding profiles. In contrast, over the same pH range, the WT and DM have distinct structural ion binding profiles; for the WT, binding is favored at lower pH, whereas the DM shows no pH dependence. These data localize the structural ion to the vicinity of the C41 motif. An overall model is presented that accommodates binding affinity, coupling, and positioning of the two metal ions and the two protons within the ribozyme. The data suggest that a protonated base triple allows the WT ribozyme to maintain appreciable activity at acidic pH, which could play an important role in the life cycle of the virus.

Dynamics-based amplification of RNA function and its characterization by using NMR spectroscopy

The ever-increasing cellular roles ascribed to RNA raise fundamental questions regarding how a biopolymer composed of only four chemically similar building-block nucleotides achieves such functional diversity. Here, I discuss how RNA achieves added mechanistic and chemical complexity by undergoing highly controlled conformational changes in response to a variety of cellular signals. I examine pathways for achieving selectivity in these conformational changes that rely to different extents on the structure and dynamics of RNA. Finally, I review solution-state NMR techniques that can be used to characterize RNA structural dynamics and its relationship to function.

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.

Mesoscopic modeling for nucleic acid chain dynamics

To gain a deeper insight into cellular processes such as transcription and translation, one needs to uncover the mechanisms controlling the configurational changes of nucleic acids. As a step toward this aim, we present here a mesoscopic-level computational model that provides a new window into nucleic acid dynamics. We model a single-stranded nucleic as a polymer chain whose monomers are the nucleosides. Each monomer comprises a bead representing the sugar molecule and a pin representing the base. The bead-pin complex can rotate about the backbone of the chain. We consider pairwise stacking and hydrogen-bonding interactions. We use a modified Monte Carlo dynamics that splits the dynamics into translational bead motion and rotational pin motion. By performing a number of tests, we first show that our model is physically sound. We then focus on a study of the kinetics of a DNA hairpin--a single-stranded molecule comprising two complementary segments joined by a noncomplementary loop--studied experimentally. We find that results from our simulations agree with experimental observations, demonstrating that our model is a suitable tool for the investigation of the hybridization of single strands.

Deprotonation stimulates productive folding in allosteric TRAP hammerhead ribozymes

Hammerhead ribozymes in crystals change conformation in response to deprotonation of the nucleophilic 2' OH, thereby aligning the hydroxyl for in-line displacement at the scissile phosphate. Published data do not address whether deprotonation affects folding in solution. Allosteric hammerhead "TRAPs," when activated by the appropriate oligonucleotide, show the expected log-linear relation between initial cleavage rate and pH. In contrast, attenuated TRAPs shows biphasic kinetics in which a rapid burst is followed by slow cleavage that is nearly independent of pH. Attenuated ribozymes are stimulated by urea at both low and high pH, confirming that rearrangement of secondary structure is rate-limiting for the attenuated ribozymes once they have folded. Plots of burst magnitude versus pH in the absence of urea show a sharp transition around pH 8.3, which is near the kinetic pKa for the cleavage reaction in Mg2+. Raising the pH after folding at pH 7.5 did not activate attenuated ribozymes even when the RNA was incubated at the elevated pH for extended periods prior to addition of Mg2+. In contrast, lowering the pH after folding at pH 9.5 rapidly re-established attenuation. Deprotonation of the ribozyme-substrate complex thus appears to alter the folding landscape such that a metastable "pre-activated" complex forms before the thermodynamically more stable attenuated state can be attained. From the initial partition into active and inactive conformers, we estimate that this deprotonation contributes approximately 1.2 kcal/mol toward stabilization of the active fold at a crucial step during folding of the TRAP. Assuming that the nucleophilic 2' OH is the relevant acid, its deprotonation would thus serve a dual role of favoring productive fold and enhancing the nucleophilicity of this oxygen.

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.

Cross-linking experiments reveal the presence of novel structural features between a hepatitis delta virus ribozyme and its substrate

The kinetic pathway of a trans-acting delta ribozyme includes an essential structural rearrangement involving the P1 stem, a stem that is formed between the substrate and the ribozyme. We performed cross-linking experiments to determine the substrate position within the catalytic center of an antigenomic, trans-acting, delta ribozyme. Substrates that included a 4-thiouridine either in position -1, +4, or +8 (i.e., adjacent to the cleavage site, or located either in the middle of or at the 3'-end of the P1 stem, respectively) were synthesized and shown to be efficiently cleaved. Examination of the cross-linking conditions, the use of various mutated ribozymes, as well as the probing and characterization of the resulting ribozyme-substrate complexes, revealed several new features of the molecular mechanism: (1) the close proximity of several bases between nucleotides of the substrate and ribozyme; (2) the active ribozyme-substrate complex folds in a manner that docks the middle of the P1 stem on the P3 stem, while concomitantly the scissile phosphate is in close proximity to the catalytic cytosine; and, (3) some complexes appear to be compatible with being active intermediates along the folding pathway, while others seem to correspond to misfolded structures. To provide a model representation of these data, a three-dimensional structure of the delta ribozyme was developed using several RNA bioinformatic software packages.

Inhibition of Escherichia coli RNase P by oligonucleotide directed misfolding of RNA

Oligonucleotide directed misfolding of RNA (ODMiR) uses short oligonucleotides to inhibit RNA function by exploiting the ability of RNA to fold into different structures with similar free energies. It is shown that the 2'-O-methyl oligonucleotide, m(CAGCCUACCCGG), can trap Escherichia coli RNase P RNA (M1 RNA) in a nonfunctional structure in a transcription mixture containing RNase P protein (C5 protein). At about 200 nM, the 12-mer thus inhibits 50% of pre-tRNA processing by RNase P. Roughly 10-fold more 12-mer is required to inhibit RNase P containing full-length, renatured RNase P RNA. Diethyl pyrocarbonate modification in the presence of 12-mer reveals increased modification of sites in and interacting with P4, suggesting a structural rearrangement of a large pseudoknot important for catalytic activity. Thus, the ODMiR method can be applied to RNAs even when folding is facilitated by a cognate protein.

RNA-tethered phenyl azide photocrosslinking via a short-lived indiscriminant electrophile

Arylazide mediated photocrosslinking has been widely used to obtain structural constraints in biological systems, even though the reactive species generated upon photolysis in aqueous solution have not been well characterized. We establish a mechanistic framework for formation of adducts between photoactivated 3-hydroxyphenyl azide and RNA. Tethered to an internal site in an RNA duplex via a 2'-amido linkage, photolysis of the aryl azide yields a cross-strand cross-link. Analysis of the ability of reagents with diagnostic reactivities to intercept formation of this cross-strand cross-link supports the assignment that the photoactivated intermediate is the ketenimine or a ketenimine-derived ring expansion product. Neither the initially produced singlet nitrene nor the subsequently formed triplet nitrene contribute to cross-link formation. Argon matrix and time-resolved solution experiments show that photolysis of free 3-hydroxyphenyl azide releases (in <or=20 ns) either a ketenimine or azepinone intermediate that reacts with nucleophiles. Adenosine, uridine, and guanosine monophosphate nucleotides have approximately equivalent abilities to quench the cross-strand cross-link, indicating that the photoactivated intermediate reacts broadly with functional groups in RNA. The reactive intermediate forms an adduct with adenosine monophosphate when tethered to both an RNA duplex or unstructured single strand; thus, cross-link formation is independent of the local RNA environment. The lifetime of the reactive intermediate generated upon photolysis of free 3-hydroxyphenyl azide in 50 mM Hepes diamine buffer is found to be 60-160 micros or significantly shorter than large scale RNA folding events. RNA-tethered 3-hydroxyphenyl azide cross-linking in aqueous buffer can thus be used with confidence to map structural neighbors, including most dynamic interactions, in RNA.