RNA pseudoknots. Stability and loop size requirements

Estimating RNA Secondary Structure Folding Free Energy Changes with efn2

A number of analyses require estimates of the folding free energy changes of specific RNA secondary structures. These predictions are often based on a set of nearest neighbor parameters that models the folding stability of a RNA secondary structure as the sum of folding stabilities of the structural elements that comprise the secondary structure. In the software suite RNAstructure, the free energy change calculation is implemented in the program efn2. The efn2 program estimates the folding free energy change and the experimental uncertainty in the folding free energy change. It can be run through the graphical user interface for RNAstructure, from the command line, or a web server. This chapter provides detailed protocols for using efn2.

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

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

A Polymer Physics Framework for the Entropy of Arbitrary Pseudoknots

The accurate prediction of RNA secondary structure from primary sequence has had enormous impact on research from the past 40 years. Although many algorithms are available to make these predictions, the inclusion of non-nested loops, termed pseudoknots, still poses challenges arising from two main factors: 1) no physical model exists to estimate the loop entropies of complex intramolecular pseudoknots, and 2) their NP-complete enumeration has impeded their study. Here, we address both challenges. First, we develop a polymer physics model that can address arbitrarily complex pseudoknots using only two parameters corresponding to concrete physical quantities-over an order of magnitude fewer than the sparsest state-of-the-art phenomenological methods. Second, by coupling this model to exhaustive enumeration of the set of possible structures, we compute the entire free energy landscape of secondary structures resulting from a primary RNA sequence. We demonstrate that for RNA structures of ∼80 nucleotides, with minimal heuristics, the complete enumeration of possible secondary structures can be accomplished quickly despite the NP-complete nature of the problem. We further show that despite our loop entropy model's parametric sparsity, it performs better than or on par with previously published methods in predicting both pseudoknotted and non-pseudoknotted structures on a benchmark data set of RNA structures of ≤80 nucleotides. We suggest ways in which the accuracy of the model can be further improved.

Identification of the RNA Pseudoknot within the 3' End of the Porcine Reproductive and Respiratory Syndrome Virus Genome as a Pathogen-Associated Molecular Pattern To Activate Antiviral Signaling via RIG-I and Toll-Like Receptor 3

Once infected by viruses, cells can detect pathogen-associated molecular patterns (PAMPs) on viral nucleic acid by host pattern recognition receptors (PRRs) to initiate the antiviral response. Porcine reproductive and respiratory syndrome virus (PRRSV) is the causative agent of porcine reproductive and respiratory syndrome (PRRS), characterized by reproductive failure in sows and respiratory diseases in pigs of different ages. To date, the sensing mechanism of PRRSV has not been elucidated. Here, we reported that the pseudoknot region residing in the 3' untranslated regions (UTR) of the PRRSV genome, which has been proposed to regulate RNA synthesis and virus replication, was sensed as nonself by retinoic acid-inducible gene I (RIG-I) and Toll-like receptor 3 (TLR3) and strongly induced type I interferons (IFNs) and interferon-stimulated genes (ISGs) in porcine alveolar macrophages (PAMs). The interaction between the two stem-loops inside the pseudoknot structure was sufficient for IFN induction, since disruption of the pseudoknot interaction powerfully dampened the IFN induction. Furthermore, transfection of the 3' UTR pseudoknot transcripts in PAMs inhibited PRRSV replication in vitro Importantly, the predicted similar structures of other arterivirus members, including equine arteritis virus (EAV), lactate dehydrogenase-elevating virus (LDV), and simian hemorrhagic fever virus (SHFV), also displayed strong IFN induction activities. Together, in this work we identified an innate recognition mechanism by which the PRRSV 3' UTR pseudoknot region served as PAMPs of arteriviruses and activated innate immune signaling to produce IFNs that inhibit virus replication. All of these results provide novel insights into innate immune recognition during virus infection.IMPORTANCE PRRS is the most common viral disease in the pork industry. It is caused by PRRSV, a positive single-stranded RNA virus, whose infection often leads to persistent infection. To date, it is not yet clear how PRRSV is recognized by the host and what is the exact mechanism of IFN induction. Here, we investigated the nature of PAMPs on PRRSV and the associated PRRs. We found that the 3' UTR pseudoknot region of PRRSV, which has been proposed to regulate viral RNA synthesis, could act as PAMPs recognized by RIG-I and TLR3 to induce type I IFN production to suppress PRRSV infection. This report is the first detailed description of pattern recognition for PRRSV, which is important in understanding the antiviral response of arteriviruses, especially PRRSV, and extends our knowledge on virus recognition.

RNA Structure: Advances and Assessment of 3D Structure Prediction

Biological functions of RNA molecules are dependent upon sustained specific three-dimensional (3D) structures of RNA, with or without the help of proteins. Understanding of RNA structure is frequently based on 2D structures, which describe only the Watson-Crick (WC) base pairs. Here, we hierarchically review the structural elements of RNA and how they contribute to RNA 3D structure. We focus our analysis on the non-WC base pairs and on RNA modules. Several computer programs have now been designed to predict RNA modules. We describe the RNA-Puzzles initiative, which is a community-wide, blind assessment of RNA 3D structure prediction programs to determine the capabilities and bottlenecks of current predictions. The assessment metrics used in RNA-Puzzles are briefly described. The detection of RNA 3D modules from sequence data and their automatic implementation belong to the current challenges in RNA 3D structure prediction.

The Influenza A PB1-F2 and N40 Start Codons Are Contained within an RNA Pseudoknot

Influenza A is a negative-sense RNA virus with an eight-segment genome. Some segments encode more than one polypeptide product, but how the virus accesses alternate internal open reading frames (ORFs) is not completely understood. In segment 2, ribosomal scanning produces two internal ORFs, PB1-F2 and N40. Here, chemical mapping reveals a Mg(2+)-dependent pseudoknot structure that includes the PB1-F2 and N40 start codons. The results suggest that interactions of the ribosome with the pseudoknot may affect the level of translation for PB1-F2 and N40.

The determination of RNA folding nearest neighbor parameters

The stability of RNA secondary structure can be predicted using a set of nearest neighbor parameters. These parameters are widely used by algorithms that predict secondary structure. This contribution introduces the UV optical melting experiments that are used to determine the folding stability of short RNA strands. It explains how the nearest neighbor parameters are chosen and how the values are fit to the data. A sample nearest neighbor calculation is provided. The contribution concludes with new methods that use the database of sequences with known structures to determine parameter values.

Evolution of coordinated mutagenesis and somatic hypermutation in VH5

The VH5 human antibody gene was analyzed using a computer program (mfg) which simulates transcription, to better understand transcription-driven mutagenesis events that occur during "phase 1" of somatic hypermutation. Results show that the great majority of mutations in the non-transcribed strand occur within loops of two predicted high-stability stem-loop structures, termed SLSs 14.9 and 13.9. In fact, 89% of the 2505 mutations reported are within the encoded complementarity-determining region (CDR) and occur in loops of these high-stability structures. In vitro studies were also done and verified the existence of SLS 14.9. Following the formation of SLSs 14.9 and 13.9, a sustained period of transcriptional activity occurs within a window size of 60-70 nucleotides. During this period, the stability of these two SLSs does not change, and may provide the substrate for base exchanges and mutagenesis. The data suggest that many mutable bases are exposed simultaneously at pause sites, allowing for coordinated mutagenesis.

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.

ProbKnot: fast prediction of RNA secondary structure including pseudoknots

It is a significant challenge to predict RNA secondary structures including pseudoknots. Here, a new algorithm capable of predicting pseudoknots of any topology, ProbKnot, is reported. ProbKnot assembles maximum expected accuracy structures from computed base-pairing probabilities in O(N(2)) time, where N is the length of the sequence. The performance of ProbKnot was measured by comparing predicted structures with known structures for a large database of RNA sequences with fewer than 700 nucleotides. The percentage of known pairs correctly predicted was 69.3%. Additionally, the percentage of predicted pairs in the known structure was 61.3%. This performance is the highest of four tested algorithms that are capable of pseudoknot prediction. The program is available for download at: http://rna.urmc.rochester.edu/RNAstructure.html.

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.

Fluorescence competition assay measurements of free energy changes for RNA pseudoknots

RNA pseudoknots have important functions, and thermodynamic stability is a key to predicting pseudoknots in RNA sequences and to understanding their functions. Traditional methods, such as UV melting and differential scanning calorimetry, for measuring RNA thermodynamics are restricted to temperature ranges around the melting temperature for a pseudoknot. Here, we report RNA pseudoknot free energy changes at 37 °C measured by fluorescence competition assays. Sequence-dependent studies for the loop 1−stem 2 region reveal (1) the individual nearest-neighbor hydrogen bonding (INN-HB) model provides a reasonable estimate for the free energy change when a Watson−Crick base pair in stem 2 is changed, (2) the loop entropy can be estimated by a statistical polymer model, although some penalty for certain loop sequences is necessary, and (3) tertiary interactions can significantly stabilize pseudoknots and extending the length of stem 2 may alter tertiary interactions such that the INN-HB model does not predict the net effect of adding a base pair. The results can inform writing of algorithms for predicting and/or designing RNA secondary structures.

An intermolecular RNA triplex provides insight into structural determinants for the pseudoknot stimulator of -1 ribosomal frameshifting

An efficient -1 programmed ribosomal frameshifting (PRF) signal requires an RNA slippery sequence and a downstream RNA stimulator, and the hairpin-type pseudoknot is the most common stimulator. However, a pseudoknot is not sufficient to promote -1 PRF. hTPK-DU177, a pseudoknot derived from human telomerase RNA, shares structural similarities with several -1 PRF pseudoknots and is used to dissect the roles of distinct structural features in the stimulator of -1 PRF. Structure-based mutagenesis on hTPK-DU177 reveals that the -1 PRF efficiency of this stimulator can be modulated by sequential removal of base-triple interactions surrounding the helical junction. Further analysis of the junction-flanking base triples indicates that specific stem-loop interactions and their relative positions to the helical junction play crucial roles for the -1 PRF activity of this pseudoknot. Intriguingly, a bimolecular pseudoknot approach based on hTPK-DU177 reveals that continuing triplex structure spanning the helical junction, lacking one of the loop-closure features embedded in pseudoknot topology, can stimulate -1 PRF. Therefore, the triplex structure is an essential determinant for the DU177 pseudoknot to stimulate -1 PRF. Furthermore, it suggests that -1 PRF, induced by an in-trans RNA via specific base-triple interactions with messenger RNAs, can be a plausible regulatory function for non-coding RNAs.

Identification of antisense RNA stem-loops that inhibit RNA-protein interactions using a bacterial reporter system

Many well-characterized examples of antisense RNAs from prokaryotic systems involve hybridization of the looped regions of stem–loop RNAs, presumably due to the high thermodynamic stability of the resulting loop–loop and loop–linear interactions. In this study, the identification of RNA stem–loops that inhibit U1A protein binding to the hpII RNA through RNA–RNA interactions was attempted using a bacterial reporter system based on phage λ N-mediated antitermination. As a result, loop sequences possessing 7–8 base complementarity to the 5′ region of the boxA element important for functional antitermination complex formation, but not the U1 hpII loop, were identified. In vitro and in vivo mutational analysis strongly suggested that the selected loop sequences were binding to the boxA region, and that the structure of the antisense stem–loop was important for optimal inhibitory activity. Next, in an attempt to demonstrate the ability to inhibit the interaction between the U1A protein and the hpII RNA, the rational design of an RNA stem–loop that inhibits U1A-binding to a modified hpII was carried out. Moderate inhibitory activity was observed, showing that it is possible to design and select antisense RNA stem–loops that disrupt various types of RNA–protein interactions.

Improved free energy parameters for RNA pseudoknotted secondary structure prediction

Accurate prediction of RNA pseudoknotted secondary structures from the base sequence is a challenging computational problem. Since prediction algorithms rely on thermodynamic energy models to identify low-energy structures, prediction accuracy relies in large part on the quality of free energy change parameters. In this work, we use our earlier constraint generation and Boltzmann likelihood parameter estimation methods to obtain new energy parameters for two energy models for secondary structures with pseudoknots, namely, the Dirks-Pierce (DP) and the Cao-Chen (CC) models. To train our parameters, and also to test their accuracy, we create a large data set of both pseudoknotted and pseudoknot-free secondary structures. In addition to structural data our training data set also includes thermodynamic data, for which experimentally determined free energy changes are available for sequences and their reference structures. When incorporated into the HotKnots prediction algorithm, our new parameters result in significantly improved secondary structure prediction on our test data set. Specifically, the prediction accuracy when using our new parameters improves from 68% to 79% for the DP model, and from 70% to 77% for the CC model.

C2'-endo nucleotides as molecular timers suggested by the folding of an RNA domain

A striking and widespread observation is that higher-order folding for many RNAs is very slow, often requiring minutes. In some cases, slow folding reflects the need to disrupt stable, but incorrect, interactions. However, a molecular explanation for slow folding in most RNAs is unknown. The specificity domain of the Bacillus subtilis RNase P ribozyme undergoes a rate-limiting folding step on the minute time-scale. This RNA also contains a C2'-endo nucleotide at A130 that exhibits extremely slow local conformational dynamics. This nucleotide is evolutionarily conserved and essential for tRNA recognition by RNase P. Here we show that deleting this single nucleotide accelerates folding by an order of magnitude even though this mutation does not change the global fold of the RNA. These results demonstrate that formation of a single stacking interaction at a C2'-endo nucleotide comprises the rate-determining step for folding an entire 154 nucleotide RNA. C2'-endo nucleotides exhibit slow local dynamics in structures spanning isolated helices to complex tertiary interactions. Because the motif is both simple and ubiquitous, C2'-endo nucleotides may function as molecular timers in many RNA folding and ligand recognition reactions.

Translation initiation from the ribosomal A site or the P site, dependent on the conformation of RNA pseudoknot I in dicistrovirus RNAs

Translation initiation of the second ORF of insect dicistrovirus RNA depends on an internal ribosomal entry site (IRES) in its intergenic region (IGR) and is exceptional in using a codon other than AUG and in not using the canonical initiator methionine tRNA. Studies in vitro suggest that pseudoknot I (PKI) immediately preceding the initiation codon occupies the ribosomal P site and that an elongator tRNA initiates translation from the ribosomal A site. Using dicistronic reporters carrying mutations in the initiation codon of the second ORF and mutant elongator or initiator tRNAs capable of reading these codons, we provide direct evidence for initiation from the A site in mammalian cells and, under certain conditions, also from the P site. Initiation from the A but not the P site requires PKI. Thus, PKI structure may be dynamic, and optimal IGR IRES-mediated translation of dicistroviral RNAs may require trans-acting factors to stabilize PKI.

Triplex structures in an RNA pseudoknot enhance mechanical stability and increase efficiency of -1 ribosomal frameshifting

Many viruses use programmed -1 ribosomal frameshifting to express defined ratios of structural and enzymatic proteins. Pseudoknot structures in messenger RNAs stimulate frameshifting in upstream slippery sequences. The detailed molecular determinants of pseudoknot mechanical stability and frameshifting efficiency are not well understood. Here we use single-molecule unfolding studies by optical tweezers, and frameshifting assays to elucidate how mechanical stability of a pseudoknot and its frameshifting efficiency are regulated by tertiary stem-loop interactions. Mechanical unfolding of a model pseudoknot and mutants designed to dissect specific interactions reveals that mechanical stability depends strongly on triplex structures formed by stem-loop interactions. Combining single-molecule and mutational studies facilitates the identification of pseudoknot folding intermediates. Average unfolding forces of the pseudoknot and mutants ranging from 50 to 22 picoNewtons correlated with frameshifting efficiencies ranging from 53% to 0%. Formation of major-groove and minor-groove triplex structures enhances pseudoknot stem stability and torsional resistance, and may thereby stimulate frameshifting. Better understanding of the molecular determinants of frameshifting efficiency may facilitate the development of anti-virus therapeutics targeting frameshifting.

RNA pseudoknots and the regulation of protein synthesis

RNA pseudoknots are structural elements found in almost all classes of RNA. Pseudoknots form when a single-stranded region in the loop of a hairpin base-pairs with a stretch of complementary nucleotides elsewhere in the RNA chain. This simple folding strategy is capable of generating a large number of stable three-dimensional folds that display a diverse range of highly specific functions in a variety of biological processes. The present review focuses on pseudoknots that act in the regulation of protein synthesis using cellular and viral examples to illustrate their versatility. Emphasis is placed on structurally well-defined pseudoknots that play a role in internal ribosome entry, autoregulation of initiation, ribosomal frameshifting during elongation and trans-translation.

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.

RNA structures with pseudo-knots: graph-theoretical, combinatorial, and statistical properties

The secondary structures of nucleic acids form a particularly important class of contact structures. Many important RNA molecules, however, contain pseudo-knots, a structural feature that is excluded explicitly from the conventional definition of secondary structures. We propose here a generalization of secondary structures incorporating ‘non-nested’ pseudo-knots, which we call bi-secondary structures, and discuss measures for the complexity of more general contact structures based on their graph-theoretical properties. Bi-secondary structures are planar trivalent graphs that are characterized by special embedding properties. We derive exact upper bounds on their number (as a function of the chain length n) implying that there are fewer different structures than sequences. Computational results show that the number of bi-secondary structures grows approximately like 2.35ⁿ. Numerical studies based on kinetic folding and a simple extension of the standard energy model show that the global features of the sequence-structure map of RNA do not change when pseudo-knots are introduced into the secondary structure picture. We find a large fraction of neutral mutations and, in particular, networks of sequences that fold into the same shape. These neutral networks percolate through the entire sequence space.

An RNA pseudoknot in the 3' end of the arterivirus genome has a critical role in regulating viral RNA synthesis

In the life cycle of plus-strand RNA viruses, the genome initially serves as the template for both translation of the viral replicase gene and synthesis of minus-strand RNA and is ultimately packaged into progeny virions. These various processes must be properly balanced to ensure efficient viral proliferation. To achieve this, higher-order RNA structures near the termini of a variety of RNA virus genomes are thought to play a key role in regulating the specificity and efficiency of viral RNA synthesis. In this study, we have analyzed the signals for minus-strand RNA synthesis in the prototype of the arterivirus family, equine arteritis virus (EAV). Using site-directed mutagenesis and an EAV reverse genetics system, we have demonstrated that a stem-loop structure near the 3' terminus of the EAV genome is required for RNA synthesis. We have also obtained evidence for an essential pseudoknot interaction between the loop region of this stem-loop structure and an upstream hairpin residing in the gene encoding the nucleocapsid protein. We propose that the formation of this pseudoknot interaction may constitute a molecular switch that could regulate the specificity or timing of viral RNA synthesis. This hypothesis is supported by the fact that phylogenetic analysis predicted the formation of similar pseudoknot interactions near the 3' end of all known arterivirus genomes, suggesting that this interaction has been conserved in evolution.

Viral RNA pseudoknots: versatile motifs in gene expression and replication

RNA pseudoknots are structural motifs in RNA that are increasingly recognized in viral and cellular RNAs. They have been shown to have a various roles in virus and cellular gene expression.

Pseudoknots are formed upon base pairing of a single-stranded region of RNA in the loop of a hairpin to a stretch of complementary nucleotides elsewhere in the RNA chain. This simple folding strategy can generate a large number of stable three-dimensional folds, which display a diverse range of highly specific functions.

Pseudoknot function is frequently associated with interactions with ribosomes. The inclusion of pseudoknots in an mRNA can thus confer unusual translational properties.

Many RNA viruses use pseudoknots in the control of viral RNA translation, replication and the switch between the two processes. Some satellite viruses encode ribozymes with active sites that are folded by a pseudoknot.

In cellular RNAs, pseudoknots are associated with all aspects of mRNA function and also ribosome function, as ribosomal RNAs contain numerous pseudoknots. Other essential cellular pseudoknots have been described in telomerase RNA and transfer messenger RNA.

Future research into pseudoknots will focus on structure–function relationships and bioinformatics identification of pseudoknots in genomes. The use of pseudoknots in antiviral applications could also become more widespread.

RNA pseudoknots have been identified in many different viral and cellular RNAs and are known to have various roles in virus and cellular gene expression. Here, Ian Brierley and colleagues review viral pseudoknots and the role of these structural motifs in virus gene expression and genome replication.

RNA pseudoknots are structural elements found in almost all classes of RNA. First recognized in the genomes of plant viruses, they are now established as a widespread motif with diverse functions in various biological processes. This Review focuses on viral pseudoknots and their role in virus gene expression and genome replication. Although emphasis is placed on those well defined pseudoknots that are involved in unusual mechanisms of viral translational initiation and elongation, the broader roles of pseudoknots are also discussed, including comparisons with relevant cellular counterparts. The relationship between RNA pseudoknot structure and function is also addressed.

Biphasic folding kinetics of RNA pseudoknots and telomerase RNA activity

Using a combined master equation and kinetic cluster approach, we investigate RNA pseudoknot folding and unfolding kinetics. The energetic parameters are computed from a recently developed Vfold model for RNA secondary structure and pseudoknot folding thermodynamics. The folding kinetics theory is based on the complete conformational ensemble, including all the native-like and non-native states. The predicted folding and unfolding pathways, activation barriers, Arrhenius plots, and rate-limiting steps lead to several findings. First, for the PK5 pseudoknot, a misfolded 5' hairpin emerges as a stable kinetic trap in the folding process, and the detrapping from this misfolded state is the rate-limiting step for the overall folding process. The calculated rate constant and activation barrier agree well with the experimental data. Second, as an application of the model, we investigate the kinetic folding pathways for human telomerase RNA (hTR) pseudoknot. The predicted folding and unfolding pathways not only support the proposed role of conformational switch between hairpin and pseudoknot in hTR activity, but also reveal molecular mechanism for the conformational switch. Furthermore, for an experimentally studied hTR mutation, whose hairpin intermediate is destabilized, the model predicts a long-lived transient hairpin structure, and the switch between the transient hairpin intermediate and the native pseudoknot may be responsible for the observed hTR activity. Such finding would help resolve the apparent contradiction between the observed hTR activity and the absence of a stable hairpin.

Structure-function analysis of the ribosomal frameshifting signal of two human immunodeficiency virus type 1 isolates with increased resistance to viral protease inhibitors

Expression of the pol-encoded proteins of human immunodeficiency virus type 1 (HIV-1) requires a programmed -1 ribosomal frameshift at the junction of the gag and pol coding sequences. Frameshifting takes place at a heptanucleotide slippery sequence, UUUUUUA, and is enhanced by a stimulatory RNA structure located immediately downstream. In patients undergoing viral protease (PR) inhibitor therapy, a p1/p6(gag) L449F cleavage site (CS) mutation is often observed in resistant isolates and frequently generates, at the nucleotide sequence level, a homopolymeric and potentially slippery sequence (UUUUCUU to UUUUUUU). The mutation is located within the stimulatory RNA downstream of the authentic slippery sequence and could act to augment levels of pol-encoded enzymes to counteract the PR deficit. Here, RNA secondary structure probing was employed to investigate the structure of a CS-containing frameshift signal, and the effect of this mutation on ribosomal frameshift efficiency in vitro and in tissue culture cells was determined. A second mutation, a GGG insertion in the loop of the stimulatory RNA that could conceivably lead to resistance by enhancing the activity of the structure, was also tested. It was found, however, that the CS and GGG mutations had only a very modest effect on the structure and activity of the HIV-1 frameshift signal. Thus the increased resistance to viral protease inhibitors seen with HIV-1 isolates containing mutations in the frameshifting signal is unlikely to be accounted for solely by enhancement of frameshift efficiency.

Predicting RNA pseudoknot folding thermodynamics

Based on the experimentally determined atomic coordinates for RNA helices and the self-avoiding walks of the P (phosphate) and C₄ (carbon) atoms in the diamond lattice for the polynucleotide loop conformations, we derive a set of conformational entropy parameters for RNA pseudoknots. Based on the entropy parameters, we develop a folding thermodynamics model that enables us to compute the sequence-specific RNA pseudoknot folding free energy landscape and thermodynamics. The model is validated through extensive experimental tests both for the native structures and for the folding thermodynamics. The model predicts strong sequence-dependent helix-loop competitions in the pseudoknot stability and the resultant conformational switches between different hairpin and pseudoknot structures. For instance, for the pseudoknot domain of human telomerase RNA, a native-like and a misfolded hairpin intermediates are found to coexist on the (equilibrium) folding pathways, and the interplay between the stabilities of these intermediates causes the conformational switch that may underlie a human telomerase disease.

Characterization of the frameshift signal of Edr, a mammalian example of programmed -1 ribosomal frameshifting

The ribosomal frameshifting signal of the mouse embryonal carcinoma differentiation regulated (Edr) gene represents the sole documented example of programmed −1 frameshifting in mammalian cellular genes [Shigemoto,K., Brennan,J., Walls,E,. Watson,C.J., Stott,D., Rigby,P.W. and Reith,A.D. (2001), Nucleic Acids Res., 29, 4079–4088]. Here, we have employed site-directed mutagenesis and RNA structure probing to characterize the Edr signal. We began by confirming the functionality and magnitude of the signal and the role of a GGGAAAC motif as the slippery sequence. Subsequently, we derived a model of the Edr stimulatory RNA and assessed its similarity to those stimulatory RNAs found at viral frameshift sites. We found that the structure is an RNA pseudoknot possessing features typical of retroviral frameshifter pseudoknots. From these experiments, we conclude that the Edr signal and by inference, the human orthologue PEG10, do not represent a novel ‘cellular class’ of programmed −1 ribosomal frameshift signal, but rather are similar to viral examples, albeit with some interesting features. The similarity to viral frameshift signals may complicate the design of antiviral therapies that target the frameshift process.

Site-specific modification of functional groups in genomic hepatitis delta virus (HDV) ribozyme

Human hepatitis delta (HDV) ribozyme is one of small ribozymes, such as hammerhead and hairpin ribozymes, etc. Its secondary structure shows pseudoknot structure composed of four stems (I to IV) and three single-stranded regions (SSrA, -B and -C). The 3D structure of 3'-cleaved product of genomic HDV ribozyme provided extensive information about tertiary hydrogen bonding interactions between nucleotide bases, phosphate oxygens and 2'OHs including new stem structure P1.1. To analyze the role of these hydrogen bond networks in the catalytic reaction, site-specific atomic-level modifications (such as deoxynucleotides, deoxyribosyl-2-aminopurine, deoxyribosylpurine, 7-deaza-ribonucleotide and inosine) were incorporated in the smallest trans-acting HDV ribozyme (47-mer). Kinetic analysis of these ribozyme variants demonstrated the importance of the two W-C base pairs of P1.1 for cleavage; in addition, the results suggest that all hydrogen bond interactions detected in the crystal structure involving 2'-OH and N7 atoms are present in the active ribozyme structure. In most of the variants, the relative reduction in kobs caused by substitution of the 2'-OH group correlated with the number of hydrogen bonds affected by the substitution. However G74 and C75 may have more than one hydrogen bond involving the 2'-OH in both the trans- and cis-acting HDV ribozyme. Moreover, in variants in which N7 was deleted, kobs was reduced 5- to 15-fold, it may suggest that N7 assists in coordinating Mg2+ ions or water molecules which bind with weak affinity in the active structure.

Do mRNA and rRNA binding sites of E.coli ribosomal protein S15 share common structural determinants?

Escherichia coli ribosomal protein S15 recognizes two RNA targets: a three-way junction in 16S rRNA and a pseudoknot structure on its own mRNA. Binding to mRNA occurs when S15 is expressed in excess over its rRNA target, resulting in an inhibition of translation start. The sole apparent similarity between the rRNA and mRNA targets is the presence of a G-U/G-C motif that contributes only modestly to rRNA binding but is essential for mRNA. To get more information on the structural determinants used by S15 to bind its mRNA target as compared to its rRNA site, we used site-directed mutagenesis, substitution by nucleotide analogs, footprinting experiments on both RNA and protein, and graphic modeling. The size of the mRNA-binding site could be reduced to 45 nucleotides, without loss of affinity. This short RNA preferentially folds into a pseudoknot, the formation of which depends on magnesium concentration and temperature. The size of the loop L2 that bridges the two stems of the pseudoknot through the minor groove could not be reduced below nine nucleotides. Then we showed that the pseudoknot recognizes the same side of S15 as 16S rRNA, although shielding a smaller surface area. It turned out that the G-U/G-C motif is recognized from the minor groove in both cases, and that the G-C pair is recognized in a very similar manner. However, the wobble G-U pair of the mRNA is not directly contacted by S15, as in rRNA, but is most likely involved in building a precise conformation of the RNA, essential for binding. Otherwise, unique specific features are utilized, such as the three-way junction in the case of 16S rRNA and the looped out A(-46) for the mRNA pseudoknot.

In vitro selection of hairpin ribozymes activated with short oligonucleotides

We have carried out an in vitro selection to obtain an allosteric hairpin ribozyme, which has cleavage activity in the presence of an exogenous short oligonucleotide as a regulator. Random sequences were inserted in a region corresponding to the hairpin loop of the ribozyme. After 12 rounds of selection, DNA templates were cloned. Of a total of 34 clones, 18 contained the same sequence, and the obtained hairpin ribozymes showed the cleavage activity specifically in the presence of the regulator oligonucleotide. All of the clones contained sequences complementary to the regulator oligonucleotide. The ribozymes with high cleavage activities gained characteristic hairpin loops at the random domain, which were similar to each other. In the absence of the oligonucleotide, the loop domain within the allosteric ribozyme probably forms a slipped hairpin loop, and the complementary sequence, with the regulator oligonucleotide located at the single stranded loop, would allow easy access of the oligonucleotide. The binding of the regulator oligonucleotide triggers a structural change of the hairpin loop to form an active conformation. Furthermore, we constructed an allosteric hammerhead ribozyme by introducing the characteristic hairpin loop. The modified hammerhead ribozyme was also changed to an allosteric ribozyme, which was activated by the addition of the regulator oligonucleotide. The characteristic hairpin loop, which was proved to be regulated by an exogenous oligonucleotide in this report, may be used to control RNA functions in various fields.

In vitro structure-function studies of the Bacillus subtilis tyrS mRNA antiterminator: evidence for factor-independent tRNA acceptor stem binding specificity

Expression of many aminoacyl-tRNA synthetase, amino acid biosynthesis and transport genes in Bacillus subtilis is controlled at the level of transcription termination using the T box system and requires the formation of specific secondary structures in the mRNA leader region. One structure functions as a transcriptional terminator, while an alternate form, the antiterminator, is necessary for transcription of the downstream coding regions. We have investigated the interaction of antiterminator model RNAs, based on the B.subtilis tyrS antiterminator with tRNA(Tyr) and tRNA acceptor stem models, using a gel shift assay. Binding of the antiterminator RNA to tRNA(Tyr) was dependent on complimentarity with the acceptor end of the tRNA or microhelix; affinity for the microhelix RNA was reduced relative to the tRNA. Alteration of a conserved position in the non-base pairing region of the bulge greatly reduced tRNA binding, consistent with in vivo studies. Therefore, it appears that some of the antiterminator-tRNA binding specificity is dependent on the structure of the antiterminator bulge alone and the complex it forms with tRNA in the absence of additional trans-acting factors. During the course of these studies we also discovered that the antiterminator can form a 'kissing' bulge complex, a unique RNA motif. The ease of formation of this RNA homodimer illustrates the propensity for the bulge of the antiterminator to bind RNA.

Modulation of RNA function by oligonucleotides recognizing RNA structure

Numerous RNA structures are responsible for regulatory processes either because they constitute a signal, like the hairpins or pseudoknots involved in ribosomal frameshifting, or because they are binding sites for proteins such as the trans-activating responsive RNA element of the human immunodeficiency virus whose binding to the viral protein Tat and cellular proteins allows full-length transcription of the retroviral genome. Selective ligands able to bind with high affinity to such RNA motifs may serve as tools for dissecting the molecular mechanisms in which they are involved. Such ligands might also constitute prototypes of therapeutic agents when RNA structures play a role in the expression of dysfunctional genes or in the multiplication of pathogens. Different classes of ligands (aminoglycosides, interacalating agents, peptides) are of interest to this aim. However, oligonucleotides deserve particular consideration. They have been extensively used in the frame of the antisense strategy. The apparent simplicity of this rational approach is, at first sight, very attractive. Indeed, numerous successful studies have been published describing the efficient inhibition of translation, splicing, or reverse transcription in cell-free systems, in cultured cells, or in vivo by oligomers complementary to an RNA region. However, RNA structures restrict the access of the target site to the antisense sequence: The competition between the intramolecular association of RNA regions weakens or even abolishes the antisense effect. Various possibilities have been developed to circumvent this limitation. This includes both rational and combinatorial strategies. High-affinity oligomers were designed to invade the RNA structure. Alternatively, triplex-forming oligonucleotides (TFO) and aptamers may recognize the folded RNA motif. Whereas the use of TFOs is rather limited owing to the strong sequence constraints for triple-helix formation, in vitro selection offers a way to explore vast oligoribo or oligodeoxyribo libraries to identify strong, selective oligonucleotide binders. The candidates (aptamers) selected against the TAR RNA element of HIV-1, which form stable loop-loop (kissing) complexes with the target, provide interesting examples of oligonucleotides recognizing a functional RNA structure through an important contribution of tertiary interactions.

On the role of magnesium ions in RNA stability

Divalent cations, like magnesium, are crucial for the structural integrity and biological activity of RNA. In this article, we present a picture of how magnesium stabilizes a particular folded form of RNA. The overall stabilization of RNA by Mg2+ is given by the free energy of transferring RNA from a reference univalent salt solution to a mixed salt solution. This term has favorable energetic contributions from two distinct modes of binding: diffuse binding and site binding. In diffuse binding, fully hydrated Mg ions interact with the RNA via nonspecific long-range electrostatic interactions. In site binding, dehydrated Mg2+ interacts with anionic ligands specifically arranged by the RNA fold to act as coordinating ligands for the mental ion. Each of these modes has a strong coulombic contribution to binding; however, site binding is also characterized by substantial changes in ion solvation and other nonelectrostatic contributions. We will show how these energetic differences can be exploited to experimentally distinguish between these two classes of ions using analyses of binding polynomials. We survey a number of specific systems in which Mg(2+)-RNA interactions have been studied. In well-characterized systems such as certain tRNAs and some rRNA fragments these studies show that site-bound ions can play an important role in RNA stability. However, the crucial role of diffusely bound ions is also evident. We emphasize that diffuse binding can only be described rigorously by a model that accounts for long-range electrostatic forces. To fully understand the role of magnesium ions in RNA stability, theoretical models describing electrostatic forces in systems with complicated structures must be developed.

Structure, stability and function of RNA pseudoknots involved in stimulating ribosomal frameshifting

Programmed -1 ribosomal frameshifting has become the subject of increasing interest over the last several years, due in part to the ubiquitous nature of this translational recoding mechanism in pathogenic animal and plant viruses. All cis-acting frameshift signals encoded in mRNAs are minimally composed of two functional elements: a heptanucleotide "slippery sequence" conforming to the general form X XXY YYZ, followed by an RNA structural element, usually an H-type RNA pseudoknot, positioned an optimal number of nucleotides (5 to 9) downstream. The slippery sequence itself promotes a low level ( approximately 1 %) of frameshifting; however, downstream pseudoknots stimulate this process significantly, in some cases up to 30 to 50 %. Although the precise molecular mechanism of stimulation of frameshifting remains poorly understood, significant advances have been made in our knowledge of the three-dimensional structures, thermodynamics of folding, and functional determinants of stimulatory RNA pseudoknots derived from the study of several well-characterized frameshift signals. These studies are summarized here and provide new insights into the structural requirements and mechanism of programmed -1 ribosomal frameshifting.

Energetics of a strongly pH dependent RNA tertiary structure in a frameshifting pseudoknot

Retroviruses employ -1 translational frameshifting to regulate the relative concentrations of structural and non-structural proteins critical to the viral life cycle. The 1.6 A crystal structure of the -1 frameshifting pseudoknot from beet western yellows virus reveals, in addition to Watson-Crick base-pairing, many loop-stem RNA tertiary structural interactions and a bound Na(+). Investigation of the thermodynamics of unfolding of the beet western yellows virus pseudoknot reveals strongly pH-dependent loop-stem tertiary structural interactions which stabilize the molecule, contributing a net of DeltaH approximately -30 kcal mol(-1) and DeltaG degrees (37) of -3.3 kcal mol(-1) to a total DeltaH and DeltaG degrees (37) of -121 and -16 kcal mol(-1), respectively, at pH 6.0, 0.5 M K(+) by DSC. Characterization of mutant RNAs supports the presence of a C8(+).G12-C26 loop 1-stem 2 base-triple (pK(a)=6.8), protonation of which contributes nearly -3.5 kcal mol(-1) in net stability in the presence of a wild-type loop 2. Substitution of the nucleotides in loop 2 with uridine bases, which would eliminate the minor groove triplex, destroys pseudoknot formation. An examination of the dependence of the monovalent ion and type on melting profiles suggests that tertiary structure unfolding occurs in a manner quantitatively consistent with previous studies on the stabilizing effects of K(+), NH(4)(+) and Na(+) on other simple duplex and pseudoknotted RNAs.

Requirement for canonical base pairing in the short pseudoknot structure of genomic hepatitis delta virus ribozyme

The tertiary structure of the 3'-cleaved product of the genomic hepatitis delta virus (HDV) ribozyme was solved by X-ray crystallographic analysis. In this structure, three single-stranded regions (SSrA, -B and -C) interact intricately with one another via hydrogen bonds between nucleotide bases, phosphate oxygens and 2'-OHs to form a nested double pseudoknot structure. Among these interactions, two Watson-Crick (W-C) base pairs, 726G-710C and 727G-709C, that form between SSrA and SSrC (P1.1) seem to be especially important for compact folding. To characterize the importance of these base pairs, ribozymes were subjected to in vitro selection from a pool of RNA molecules randomly substituted at positions 709, 710, 726 and 727. The results establish the importance of the two W-C base pairs for activity, although some mutants are active with one G-C base pair. In addition, the kinetic parameters were analyzed in all 16 combinations with two canonical base pairs. Comparison of variant ribozymes with the wild-type ribozyme reveals that the difference in reaction rates for these variants (DeltaDelta G (double dagger)) is not simply accounted for by the differences in the stability of P1.1 (DeltaDelta G (0)(37)). The role played by Mg(2+)ions in formation of the P1.1 structure is also discussed.

Identification of the minimal essential RNA sequences responsible for site-specific targeting of the Leishmania RNA virus 1-4 capsid endoribonuclease

The Leishmania RNA virus 1-4 capsid protein possesses an endoribonuclease activity responsible for single-site-specific cleavage within the 450-nucleotide 5' untranslated region of its own viral RNA transcript. To characterize the minimal essential RNA determinants required for site-specific cleavage, mutated RNA transcripts were examined for susceptibility to cleavage by the virus capsid protein in an in vitro assay. Deletion analyses revealed that all determinants necessary for accurate cleavage are encoded in viral nucleotides 249 to 342. Nuclease mapping and site-specific mutagenesis of the minimal RNA sequence defined a stem-loop structure that is located 40 nucleotides upstream from the cleavage site (nucleotide 320) and that is essential for accurate RNA cleavage. Abrogation of cleavage by disruption of base pairing within the stem-loop was reversed through the introduction of complementary nucleotide substitutions that reestablished the structure. We also provide evidence that divalent cations, essential components of the cleavage reaction, stabilized the stem-loop structure in solution. That capsid-specific antiserum eliminated specific RNA cleavage provides further evidence that the virus capsid gene encodes the essential endoribonuclease activity.

Stitching together RNA tertiary architectures

The powerful explanatory paradigm of molecular biology requiring form to co-evolve with function has again been proven successful when, over the recent two decades, a wealth of biological functions have been uncovered for RNA. Previously considered as a mere mediator of the genetic code, RNA is now acknowledged as a key player in a wide variety of cellular processes. Along with the discovery of novel biological functions of RNA molecules, a number of RNA three-dimensional structures have been solved which beautifully demonstrate the molecular adaptability which allows RNA to participate as a key player in these functions. A distinct repertoire of molecular motifs provides a basis for the assembly of complex RNA tertiary architectures.

How RNA folds

We describe the RNA folding problem and contrast it with the much more difficult protein folding problem. RNA has four similar monomer units, whereas proteins have 20 very different residues. The folding of RNA is hierarchical in that secondary structure is much more stable than tertiary folding. In RNA the two levels of folding (secondary and tertiary) can be experimentally separated by the presence or absence of Mg2+. Secondary structure can be predicted successfully from experimental thermodynamic data on secondary structure elements: helices, loops, and bulges. Tertiary interactions can then be added without much distortion of the secondary structure. These observations suggest a folding algorithm to predict the structure of an RNA from its sequence. However, to solve the RNA folding problem one needs thermodynamic data on tertiary structure interactions, and identification and characterization of metal-ion binding sites. These data, together with force versus extension measurements on single RNA molecules, should provide the information necessary to test and refine the proposed algorithm.

A phylogenetically conserved hairpin-type 3' untranslated region pseudoknot functions in coronavirus RNA replication

Secondary and tertiary structures in the 3' untranslated region (UTR) of plus-strand RNA viruses have been postulated to function as control elements in RNA replication, transcription, and translation. Here we describe a 54-nucleotide (nt) hairpin-type pseudoknot within the 288-nt 3' UTR of the bovine coronavirus genome and show by mutational analysis of both stems that the pseudoknotted structure is required for the replication of a defective interfering RNA genome. The pseudoknot is phylogenetically conserved among coronaviruses both in location and in shape but only partially in nucleotide sequence, and evolutionary covariation of bases to maintain G. U pairings indicates that it functions in the plus strand. RNase probing of synthetic transcripts provided additional evidence of its tertiary structure and also identified the possible existence of two conformational states. These results indicate that the 3' UTR pseudoknot is involved in coronavirus RNA replication and lead us to postulate that it functions as a regulatory control element.

An NMR and mutational analysis of an RNA pseudoknot of Escherichia coli tmRNA involved in trans-translation

Transfer-messenger RNA (tmRNA) is a unique molecule that combines properties from both tRNA and mRNA, and facilitates a novel translation reaction termed trans -translation. According to phylogenetic sequence analysis among various bacteria and chemical probing analysis, the secondary structure of the 350-400 nt RNA is commonly characterized by a tRNA-like structure, and four pseudoknots with different sizes. A mutational analysis using a number of Escherichia coli tmRNA variants as well as a chemical probing analysis has recently demonstrated not only the presence of the smallest pseudoknot, PK1, upstream of the internal coding region, but also its direct implication in trans -translation. Here, NMR methods were used to investigate the structure of the 31 nt pseudoknot PK1 and its 11 mutants in which nucleotide substitutions are introduced into each of two stems or the linking loops. NMR results provide evidence that the PK1 RNA is folded into a pseudoknot structure in the presence of Mg(2+). Imino proton resonances were observed consistent with formation of two helical stem regions and these stems stacked to each other as often seen in pseudoknot structures, in spite of the existence of three intervening nucleo-tides, loop 3, between the stems. Structural instability of the pseudoknot structure, even in the presence of Mg(2+), was found in the PK1 mutants except in the loop 3 mutants which still maintained the pseudoknot folding. These results together with their biological activities indicate that trans -translation requires the pseudoknot structure stabilized by Mg(2+)and specific residues G61 and G62 in loop 3.

Sequence-dependent conformational differences of small RNAs revealed by native gel electrophoresis

In this study, we use native polyacrylamide gel electrophoresis and one-dimensional NMR spectroscopy to analyze small RNA hairpins containing a UUCG tetraloop. The aggregation state of one RNA 16-mer (5'-CGGCUUCGGUCGACCA-3') in the presence of Mg(2+) was confirmed by laser light scattering. Although it is widely known in the RNA field that some RNAs tend to aggregate, especially when present at high concentrations, the sequence elements responsible for this effect are rarely identified. In this work, we show that Mg(2+)-induced aggregation of the 16-mer RNA hairpin is sensitive to the presence of the 3'-terminal base and a specific 2'-hydroxyl group. Our study highlights the fact that even small changes in a particular RNA sequence can increase its tendency to undergo Mg(2+)-dependent aggregation in an unpredictable manner. Our analysis also shows that native gel electrophoresis is a sensitive probe of RNA conformation with the capability to detect differences apparently caused by subtle base stacking effects at the ends of helices.

Solution structure and thermodynamics of a divalent metal ion binding site in an RNA pseudoknot

Identification and characterization of a metal ion binding site in an RNA pseudoknot was accomplished using cobalt (III) hexammine, Co(NH3)63+, as a probe for magnesium (II) hexahydrate, Mg(H2O)62+, in nuclear magnetic resonance (NMR) structural studies. The pseudoknot causes efficient -1 ribosomal frameshifting in mouse mammary tumor virus. Divalent metal ions, such as Mg2+, are critical for RNA structure and function; Mg2+preferentially stabilizes the pseudoknot relative to its constituent hairpins. The use of Co(NH3)63+as a substitute for Mg2+was investigated by ultraviolet absorbance melting curves, NMR titrations of the imino protons, and analysis of NMR spectra in the presence of Mg2+or Co (NH3)63+. The structure of the pseudoknot-Co(NH3)63+complex reveals an ion-binding pocket formed by a short, two-nucleotide loop and the major groove of a stem. Co(NH3)63+stabilizes the sharp loop-to-stem turn and reduces the electrostatic repulsion of the phosphates in three proximal strands. Hydrogen bonds are identified between the Co(NH3)63+protons and non-bridging phosphate oxygen atoms, 2' hydroxyl groups, and nitrogen and oxygen acceptors on the bases. The binding site is significantly different from that previously characterized in the major groove surface of tandem G.U base-pairs, but is similar to those observed in crystal structures of a fragment of the 5 S rRNA and the P5c helix of the Tetrahymena thermophila group I intron. Changes in chemical shifts occurred at the same pseudoknot protons on addition of Mg2+as on addition of Co(NH3)63+, indicating that both ions bind at the same site. Ion binding dissociation constants of approximately 0.6 mM and 5 mM (in 200 mM Na+and a temperature of 15 degrees C) were obtained for Co(NH3)63+and Mg2+, respectively, from the change in chemical shift as a function of metal ion concentration. An extensive array of non-sequence-specific hydrogen bond acceptors coupled with conserved structural elements within the binding pocket suggest a general mode of divalent metal ion stabilization of this type of frameshifter pseudoknot. These results provide new thermodynamic and structural insights into the role divalent metal ions play in stabilizing RNA tertiary structural motifs such as pseudoknots.

Equilibrium unfolding pathway of an H-type RNA pseudoknot which promotes programmed -1 ribosomal frameshifting

The equilibrium unfolding pathway of a 41-nucleotide frameshifting RNA pseudoknot from the gag-pro junction of mouse intracisternal A-type particles (mIAP), an endogenous retrovirus, has been determined through analysis of dual optical wavelength, equilibrium thermal melting profiles and differential scanning calorimetry. The mIAP pseudoknot is an H-type pseudoknot proposed to have structural features in common with the gag-pro frameshifting pseudoknots from simian retrovirus-1 (SRV-1) and mouse mammary tumor virus (MMTV). In particular, the mIAP pseudoknot is proposed to contain an unpaired adenosine base at the junction of the two helical stems (A15), as well as one in the middle of stem 2 (A35). A mutational analysis of stem 1 hairpins and compensatory base-pair substitutions incorporated into helical stem 2 was used to assign optical melting transitions to molecular unfolding events. The optical melting profile of the wild-type RNA is most simply described by four sequential two-state unfolding transitions. Stem 2 melts first in two closely coupled low-enthalpy transitions at low tmin which the stem 3' to A35, unfolds first, followed by unfolding of the remainder of the helical stem. The third unfolding transition is associated with some type of stacking interactions in the stem 1 hairpin loop not present in the pseudoknot. The fourth transition is assigned to unfolding of stem 1. In all RNAs investigated, DeltaHvH approximately DeltaHcal, suggesting that DeltaCpfor unfolding is small. A35 has the thermodynamic properties expected for an extrahelical, unpaired nucleotide. Deletion of A15 destabilizes the stem 2 unfolding transition in the context of both the wild-type and DeltaA35 mutant RNAs only slightly, by DeltaDeltaG degrees approximately 1 kcal mol-1(at 37 degrees C). The DeltaA15 RNA is considerably more susceptible to thermal denaturation in the presence of moderate urea concentrations than is the wild-type RNA, further evidence of a detectable global destabilization of the molecule. Interestingly, substitution of the nine loop 2 nucleotides with uridine residues induces a more pronounced destabilization of the molecule (DeltaDeltaG degrees approximately 2.0 kcal mol-1), a long-range, non-nearest neighbor effect. These findings provide the thermodynamic basis with which to further refine the relationship between efficient ribosomal frameshifting and pseudoknot structure and stability.

Structural studies of the RNA pseudoknot required for readthrough of the gag-termination codon of murine leukemia virus

Retroviruses, such as murine leukemia virus (MuLV), whose gag and pol genes are in the same reading frame but separated by a UAG stop codon, require that 5-10 % of ribosomes decode the UAG as an amino acid and continue translation to synthesize the Gag-Pol fusion polyprotein. A specific pseudoknot located eight nucleotides 3' of the UAG is required for this redefinition of the UAG stop codon. The structural probing and mutagenic analyses presented here provide evidence that loop I of the pseudoknot is one nucleotide, stem II has seven base-pairs, and the nucleotides 3' of stem II are important for function. Stem II is more resistant to single-strand-specific probes than stem I. Sequences upstream of the UAG codon allow formation of two competing structures, a stem-loop and the pseudoknot.

The role of RNA pseudoknot stem 1 length in the promotion of efficient -1 ribosomal frameshifting

The ribosomal frameshifting signal present in the genomic RNA of the coronavirus infectious bronchitis virus (IBV) contains a classic hairpin-type RNA pseudoknot that is believed to possess coaxially stacked stems of 11 bp (stem 1) and 6 bp (stem 2). We investigated the influence of stem 1 length on the frameshift process by measuring the frameshift efficiency in vitro of a series of IBV-based pseudoknots whose stem 1 length was varied from 4 to 13 bp in single base-pair increments. Efficient frameshifting depended upon the presence of a minimum of 11 bp; pseudoknots with a shorter stem 1 were either non-functional or had reduced frameshift efficiency, despite the fact that a number of them had a stem 1 with a predicted stability equal to or greater than that of the wild-type IBV pseudoknot. An upper limit for stem 1 length was not determined, but pseudoknots containing a 12 or 13 bp stem 1 were fully functional. Structure probing analysis was carried out on RNAs containing either a ten or 11 bp stem 1; these experiments confirmed that both RNAs formed pseudoknots and appeared to be indistinguishable in conformation. Thus the difference in frameshifting efficiency seen with the two structures was not simply due to an inability of the 10 bp stem 1 construct to fold into a pseudoknot. In an attempt to identify other parameters which could account for the poor functionality of the shorter stem 1-containing pseudoknots, we investigated, in the context of the 10 bp stem 1 construct, the influence on frameshifting of altering the slippery sequence-pseudoknot spacing distance, loop 2 length, and the number of G residues at the bottom of the 5'-arm of stem 1. For each parameter, it was possible to find a condition where a modest stimulation of frameshifting was observable (about twofold, from seven to a maximal 17 %), but we were unable to find a situation where frameshifting approached the levels seen with 11 bp stem 1 constructs (48-57 %). Furthermore, in the next smaller construct (9 bp stem 1), changing the bottom four base-pairs to G.C (the optimal base composition) only stimulated frameshifting from 3 to 6 %, an efficiency about tenfold lower than seen with the 11 bp construct. Thus stem 1 length is a major factor in determining the functionality of this class of pseudoknot and this has implications for models of the frameshift process.

Evidence for an RNA pseudoknot loop-helix interaction essential for efficient -1 ribosomal frameshifting

RNA pseudoknots are structural elements that participate in a variety of biological processes. At -1 ribosomal frameshifting sites, several types of pseudoknot have been identified which differ in their organisation and functionality. The pseudoknot found in infectious bronchitis virus (IBV) is typical of those that possess a long stem 1 of 11-12 bp and a long loop 2 (30-164 nt). A second group of pseudoknots are distinguishable that contain stems of only 5 to 7 bp and shorter loops. The NMR structure of one such pseudoknot, that of mouse mammary tumor virus (MMTV), has revealed that it is kinked at the stem 1-stem 2 junction, and that this kinked conformation is essential for efficient frameshifting. We recently investigated the effect on frameshifting of modulating stem 1 length and stability in IBV-based pseudoknots, and found that a stem 1 with at least 11 bp was needed for efficient frameshifting. Here, we describe the sequence manipulations that are necessary to bypass the requirement for an 11 bp stem 1 and to convert a short non-functional IBV-derived pseudoknot into a highly efficient, kinked frameshifter pseudoknot. Simple insertion of an adenine residue at the stem 1-stem 2 junction (an essential feature of a kinked pseudoknot) was not sufficient to create a functional pseudoknot. An additional change was needed: efficient frameshifting was recovered only when the last nucleotide of loop 2 was changed from a G to an A. The requirement for an A at the end of loop 2 is consistent with a loop-helix contact similar to those described in other RNA tertiary structures. A mutational analysis of both partners of the proposed interaction, the loop 2 terminal adenine residue and two G.C pairs near the top of stem 1, revealed that the interaction was essential for efficient frameshifting. The specific requirement for a 3'-terminal A residue was lost when loop 2 was increased from 8 to 14 nt, suggesting that the loop-helix contact may be required only in those pseudoknots with a short loop 2.