Search

Mosquito-borne flaviviruses (MBFVs) including dengue, West Nile, yellow fever, and Zika viruses have an RNA genome encoding one open reading frame flanked by 5′ and 3′ untranslated regions (UTRs). The 3′ UTRs of MBFVs contain regions of high sequence conservation in structured RNA elements known as dumbbells (DBs). DBs regulate translation and replication of the viral RNA genome, functions proposed to depend on the formation of an RNA pseudoknot. To understand how DB structure provides this function, we solved the x-ray crystal structure of the Donggang virus DB to 2.1Å resolution and used structural modeling to reveal the details of its three-dimensional fold. The structure confirmed the predicted pseudoknot and molecular modeling revealed how conserved sequences form a four-way junction that appears to stabilize the pseudoknot. Single-molecule FRET suggests that the DB pseudoknot is a stable element that can regulate the switch between translation and replication during the viral lifecycle by modulating long-range RNA conformational changes.

The rapid spread of COVID-19 is motivating development of antivirals targeting conserved SARS-CoV-2 molecular machinery. The SARS-CoV-2 genome includes conserved RNA elements that offer potential small-molecule drug targets, but most of their 3D structures have not been experimentally characterized. Here, we provide a compilation of chemical mapping data from our and other labs, secondary structure models, and 3D model ensembles based on Rosetta's FARFAR2 algorithm for SARS-CoV-2 RNA regions including the individual stems SL1-8 in the extended 5′ UTR; the reverse complement of the 5′ UTR SL1-4; the frameshift stimulating element (FSE); and the extended pseudoknot, hypervariable region, and s2m of the 3′ UTR. For eleven of these elements (the stems in SL1–8, reverse complement of SL1–4, FSE, s2m and 3′ UTR pseudoknot), modeling convergence supports the accuracy of predicted low energy states; subsequent cryo-EM characterization of the FSE confirms modeling accuracy. To aid efforts to discover small molecule RNA binders guided by computational models, we provide a second set of similarly prepared models for RNA riboswitches that bind small molecules. Both datasets (‘FARFAR2-SARS-CoV-2’, https://github.com/DasLab/FARFAR2-SARS-CoV-2; and ‘FARFAR2-Apo-Riboswitch’, at https://github.com/DasLab/FARFAR2-Apo-Riboswitch’) include up to 400 models for each RNA element, which may facilitate drug discovery approaches targeting dynamic ensembles of RNA molecules., Graphical Abstract Graphical Abstract De novo 3D models of SARS-CoV-2 RNA elements and small-molecule-binding RNAs to aid drug discovery.

Recently, prokaryotic riboswitches have been identified that regulate transcription in response to change of the concentration of secondary messengers. The ZMP (5-Aminoimidazole-4-carboxamide ribonucleotide (AICAR))-sensing riboswitch from Thermosinus carboxydivorans is a transcriptional ON-switch that is involved in purine and carbon-1 metabolic cycles. Its aptamer domain includes the pfl motif, which features a pseudoknot, impeding rho-independent terminator formation upon stabilization by ZMP interaction. We herein investigate the conformational landscape of transcriptional intermediates including the expression platform of this riboswitch and characterize the formation and unfolding of the important pseudoknot structure in the context of increasing length of RNA transcripts. NMR spectroscopic data show that even surprisingly short pre-terminator stems are able to disrupt ligand binding and thus metabolite sensing. We further show that the pseudoknot structure, a prerequisite for ligand binding, is preformed in transcription intermediates up to a certain length. Our results describe the conformational changes of 13 transcription intermediates of increasing length to delineate the change in structure as mRNA is elongated during transcription. We thus determine the length of the key transcription intermediate to which addition of a single nucleotide leads to a drastic drop in ZMP affinity.

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats., Graphical Abstract Graphical AbstractWe present a smart strategy of telomerase RNA (TR) identification based on its conserved type-3 RNA Pol III promoter and TR template elements. We characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. TRs are validated experimentally and show conservation of core TR structural domains. These results shed light on the evolution of a key eukaryotic non-coding RNA across more than a billion years.

As a class of parasitic, non-coding RNAs, satellite RNAs (satRNAs) have to compete with their helper virus for limited amounts of viral and/or host resources for efficient replication, by which they usually reduce viral accumulation and symptom expression. Here, we report a cucumber mosaic virus (CMV)-associated satRNA (sat-T1) that ameliorated CMV-induced symptoms, accompanied with a significant reduction in the accumulation of viral genomic RNAs 1 and 2, which encode components of the viral replicase. Intrans replication assays suggest that the reduced accumulation is the outcome of replication competition. The structural basis of sat-T1 responsible for the inhibition of viral RNA accumulation was determined to be a three-way branched secondary structure that contains two biologically important hairpins. One is indispensable for the helper virus inhibition, and the other engages in formation of a tertiary pseudoknot structure that is essential for sat-T1 survival. The secondary structure containing the pseudoknot is the first RNA element with a biological phenotype experimentally identified in CMV satRNAs, and it is structurally conserved in most CMV satRNAs. Thus, this may be a generic method for CMV satRNAs to inhibit the accumulation of the helper virus via the newly-identified RNA structure.

We have determined the structure of the glutamine-II riboswitch ligand binding domain using X-ray crystallography. The structure was solved using a novel combination of homology modeling and molecular replacement. The structure comprises three coaxial helical domains, the central one of which is a pseudoknot with partial triplex character. The major groove of this helix provides the binding site for L-glutamine, which is extensively hydrogen bonded to the RNA. Atomic mutation of the RNA at the ligand binding site leads to loss of binding shown by isothermal titration calorimetry, explaining the specificity of the riboswitch. A metal ion also plays an important role in ligand binding. This is directly bonded to a glutamine carboxylate oxygen atom, and its remaining inner-sphere water molecules make hydrogen bonding interactions with the RNA.

In response to DNA damage, Escherichia coli cells activate the expression of the toxin gene tisB of the toxin–antitoxin system tisB-istR1. Of three isoforms, only the processed, highly structured +42 tisB mRNA is active. Translation requires a standby site, composed of two essential elements: a single-stranded region located 100 nucleotides upstream of the sequestered RBS, and a structure near the 5′-end of the active mRNA. Here, we propose that this 5′-structure is an RNA pseudoknot which is required for 30S and protein S1-alone binding to the mRNA. Point mutations that prevent formation of this pseudoknot inhibit formation of translation initiation complexes, impair S1 and 30S binding to the mRNA, and render the tisB mRNA non-toxic in vivo. A set of mutations created in either the left or right arm of stem 2 of the pseudoknot entailed loss of toxicity upon overexpression of the corresponding mRNA variants. Combining the matching right-left arm mutations entirely restored toxicity levels to that of the wild-type, active mRNA. Finally, since many pseudoknots have high affinity for S1, we predicted similar pseudoknots in non-homologous type I toxin–antitoxin systems that exhibit features similar to that of tisB-IstR1, suggesting a shared requirement for standby acting at great distances.

In vitro selections are the only known methods to generate catalytic RNAs (ribozymes) that do not exist in nature. Such new ribozymes are used as biochemical tools, or to address questions on early stages of life. In both cases, it is helpful to identify the shortest possible ribozymes since they are easier to deploy as a tool, and because they are more likely to have emerged in a prebiotic environment. One of our previous selection experiments led to a library containing hundreds of different ribozyme clusters that catalyze the triphosphorylation of their 5′-terminus. This selection showed that RNA systems can use the prebiotically plausible molecule cyclic trimetaphosphate as an energy source. From this selected ribozyme library, the shortest ribozyme that was previously identified had a length of 67 nucleotides. Here we describe a combinatorial method to identify short ribozymes from libraries containing many ribozymes. Using this protocol on the library of triphosphorylation ribozymes, we identified a 17-nucleotide sequence motif embedded in a 44-nucleotide pseudoknot structure. The described combinatorial approach can be used to analyze libraries obtained by different in vitro selection experiments.

Ribosomal protein S1 plays important roles in the translation initiation step of many Escherichia coli mRNAs, particularly those with weak Shine-Dalgarno sequences or structured 5′ UTRs, in addition to a variety of cellular processes beyond the ribosome. In all cases, the RNA-binding activity of S1 is a central feature of its function. While sequence determinants of S1 affinity and many elements of the interactions of S1 with simple secondary structures are known, mechanistic details of the protein's interactions with RNAs of more complex secondary and tertiary structure are less understood. Here, we investigate the interaction of S1 with the well-characterized H-type pseudoknot of a class-I translational preQ1 riboswitch as a highly structured RNA model whose conformation and structural dynamics can be tuned by the addition of ligands of varying binding affinity, particularly preQ1, guanine, and 2,6-diaminopurine. Combining biochemical and single molecule fluorescence approaches, we show that S1 preferentially interacts with the less folded form of the pseudoknot and promotes a dynamic, partially unfolded conformation. The ability of S1 to unfold the RNA is inversely correlated with the structural stability of the pseudoknot. These mechanistic insights delineate the scope and limitations of S1-chaperoned unfolding of structured RNAs.

Downstream stable mRNA secondary structures can stall elongating ribosomes by impeding the concerted movements of tRNAs and mRNA on the ribosome during translocation. The addition of a downstream mRNA structure, such as a stem-loop or a pseudoknot, is essential to induce -1 programmed ribosomal frameshifting (-1 PRF). Interestingly, previous studies revealed that -1 PRF efficiencies correlate with conformational plasticity of pseudoknots, defined as their propensity to form incompletely folded structures, rather than with the mechanical properties of pseudoknots. To elucidate the detailed molecular mechanisms of translocation and -1 PRF, we applied several smFRET assays to systematically examine how translocation rates and conformational dynamics of ribosomes were affected by different pseudoknots. Our results show that initial pseudoknot-unwinding significantly inhibits late-stage translocation and modulates conformational dynamics of ribosomal post-translocation complexes. The effects of pseudoknots on the structural dynamics of ribosomes strongly correlate with their abilities to induce -1 PRF. Our results lead us to propose a kinetic scheme for translocation which includes an initial power-stroke step and a following thermal-ratcheting step. This scheme provides mechanistic insights on how selective modulation of late-stage translocation by pseudoknots affects -1 PRF. Overall our findings advance current understanding of translocation and ribosome-induced mRNA structure unwinding.

© The Author(s) 2018. Published by Oxford University Press on behalf of Nucleic Acids Research. Telomerase adds telomeric repeats to chromosome ends by processive copying of a template within the telomerase RNA bound to telomerase reverse transcriptase. Telomerase RNAs have single-stranded regions that separate the template from a 5 stem and 3 pseudoknot, and mammals gained additional stem P2a.1 separating the template from the pseudoknot. Using human telomerase, we show that the length of template 3-flanking single-stranded RNA is a determinant of repeat addition processivity whereas template 5-flanking single-stranded RNA and P2a.1 are critical for activity but not processivity. In comparison, requirements for the template sequence itself are confounding: different substitutions of the same position have strikingly different consequences, from improved processivity and activity to complete inactivation. We discovered that some altered-template sequences stabilize an alternative RNA conformation that precludes the pseudoknot by base-pairing of one pseudoknot strand to the template 3 end. Using mutations to reduce over-stability of the alternative conformation, we restore high activity and processivity to otherwise inactive altered-template telomerase ribonucleoproteins. In cells, over-stabilization or destabilization of the alternative state severely inhibited biogenesis of active telomerase. Our findings delineate roles for human telomerase RNA template-flanking regions, establish a biologically relevant pseudoknot-alternative RNA conformation, and expand the repertoire of human telomerase repeat synthesis.

SAM-V is one of the class of riboswitches that bind S-adenosylmethione, regulating gene expression by controlling translation. We have solved the crystal structure of the metY SAM-V riboswitch bound to its SAM ligand at 2.5 Å resolution. The RNA folds as an H-type pseudoknot, with a major-groove triple helix in which resides the SAM ligand binding site. The bound SAM adopts an elongated conformation aligned with the axis of the triple helix, and is held at either end by hydrogen bonding to the adenine and the amino acid moieties. The central sulfonium cation makes electrostatic interactions with an U:A.U base triple, so conferring specificity. We propose a model in which SAM binding leads to association of the triplex third strand that stabilizes a short helix and occludes the ribosome binding site. Thus the new structure explains both ligand specificity and the mechanism of genetic control.

RNA molecules play important and diverse regulatory roles in the cell. Inspired by this natural versatility, RNA devices are increasingly important for many synthetic biology applications, e.g. optimizing engineered metabolic pathways, gene therapeutics or building up complex logical units. A major advantage of RNA is the possibility of de novo design of RNA-based sensing domains via an in vitro selection process (SELEX). Here, we describe development of a novel ciprofloxacin-responsive riboswitch by in vitro selection and next-generation sequencing-guided cellular screening. The riboswitch recognizes the small molecule drug ciprofloxacin with a KD in the low nanomolar range and adopts a pseudoknot fold stabilized by ligand binding. It efficiently interferes with gene expression both in lower and higher eukaryotes. By controlling an auxotrophy marker and a resistance gene, respectively, we demonstrate efficient, scalable and programmable control of cellular survival in yeast. The applied strategy for the development of the ciprofloxacin riboswitch is easily transferrable to any small molecule target of choice and will thus broaden the spectrum of RNA regulators considerably.

We introduce the SPlit-and-conQueR (SPQR) model, a coarse-grained representation of RNA designed for structure prediction and refinement. In our approach, the representation of a nucleotide consists of a point particle for the phosphate group and an anisotropic particle for the nucleoside. The interactions are, in principle, knowledge-based potentials inspired by the ESCORE function, a base-centered scoring function. However, a special treatment is given to base-pairing interactions and certain geometrical conformations which are lost in a raw knowledge-base model. This results in a representation able to describe planar canonical and non-canonical base pairs and base-phosphate interactions and to distinguish sugar puckers and glycosidic torsion conformations. The model is applied to the folding of several structures, including duplexes with internal loops of non-canonical base pairs, tetraloops, junctions and a pseudoknot. For the majority of these systems, experimental structures are correctly predicted at the level of individual contacts. We also propose a method for efficiently reintroducing atomistic detail from the coarse-grained representation., Accepted for publication on Nucleic Acids Research

Knottins, or inhibitor cystine knots (ICKs), are ultra-stable miniproteins with multiple applications in drug design and medical imaging. These widespread and functionally diverse proteins are characterized by the presence of three interwoven disulfide bridges in their structure, which form a unique pseudoknot. Since 2004, the KNOTTIN database (www.dsimb.inserm.fr/KNOTTIN/) has been gathering standardized information about knottin sequences, structures, functions and evolution. The website also provides access to bibliographic data and to computational tools that have been specifically developed for ICKs. Here, we present a major upgrade of our database, both in terms of data content and user interface. In addition to the new features, this article describes how KNOTTIN has seen its size multiplied over the past ten years (since its last publication), notably with the recent inclusion of predicted ICKs structures. Finally, we report how our web resource has proved usefulness for the researchers working on ICKs, and how the new version of the KNOTTIN website will continue to serve this active community.

Frameshifting is an essential process that regulates protein synthesis in many viruses. The ribosome may slip backward when encountering a frameshift motif on the messenger RNA, which usually contains a pseudoknot structure involving tertiary base pair interactions. Due to the lack of detailed molecular explanations, previous studies investigating which features of the pseudoknot are important to stimulate frameshifting have presented diverse conclusions. Here we constructed a bimolecular pseudoknot to dissect the interior tertiary base pairs and used single-molecule approaches to assess the structure targeted by ribosomes. We found that the first ribosome target stem was resistant to unwinding when the neighboring loop was confined along the stem; such constrained conformation was dependent on the presence of consecutive adenosines in this loop. Mutations that disrupted the distal base triples abolished all remaining tertiary base pairs. Changes in frameshifting efficiency correlated with the stem unwinding resistance. Our results demonstrate that various tertiary base pairs are coordinated inside a highly efficient frameshift-stimulating RNA pseudoknot and suggest a mechanism by which mechanical resistance of the pseudoknot may persistently act on translocating ribosomes.

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.

Telomerase maintains telomere length at the ends of linear chromosomes using an integral telomerase RNA (TER) and telomerase reverse transcriptase (TERT). An essential part of TER is the template/pseudoknot domain (t/PK) which includes the template, for adding telomeric repeats, template boundary element (TBE), and pseudoknot, enclosed in a circle by stem 1. The Tetrahymena telomerase holoenzyme catalytic core (p65-TER-TERT) was recently modeled in our 9 Å resolution cryo-electron microscopy map by fitting protein and TER domains, including a solution NMR structure of the Tetrahymena pseudoknot. Here, we describe in detail the structure and folding of the isolated pseudoknot, which forms a compact structure with major groove U•A-U and novel C•G-A+ base triples. Base substitutions that disrupt the base triples reduce telomerase activity in vitro. NMR studies also reveal that the pseudoknot does not form in the context of full-length TER in the absence of TERT, due to formation of a competing structure that sequesters pseudoknot residues. The residues around the TBE remain unpaired, potentially providing access by TERT to this high affinity binding site during an early step in TERT-TER assembly. A model for the assembly pathway of the catalytic core is proposed.

Design of RNA sequences that adopt functional folds establishes principles of RNA folding and applications in biotechnology. Inverse folding for RNAs, which allows computational design of sequences that adopt specific structures, can be utilized for unveiling RNA functions and developing genetic tools in synthetic biology. Although many algorithms for inverse RNA folding have been developed, the pseudoknot, which plays a key role in folding of ribozymes and riboswitches, is not addressed in most algorithms. For the few algorithms that attempt to predict pseudoknot-containing ribozymes, self-cleavage activity has not been tested. Herein, we design double-pseudoknot HDV ribozymes using an inverse RNA folding algorithm and test their kinetic mechanisms experimentally. More than 90% of the positively designed ribozymes possess self-cleaving activity, whereas more than 70% of negative control ribozymes, which are predicted to fold to the necessary structure but with low fidelity, do not possess it. Kinetic and mutation analyses reveal that these RNAs cleave site-specifically and with the same mechanism as the WT ribozyme. Most ribozymes react just 50- to 80-fold slower than the WT ribozyme, and this rate can be improved to near WT by modification of a junction. Thus, fast-cleaving functional ribozymes with multiple pseudoknots can be designed computationally.

Metabolite-responsive RNA pseudoknots derived from prokaryotic riboswitches have been shown to stimulate -1 programmed ribosomal frameshifting (PRF), suggesting -1 PRF as a promising gene expression platform to extend riboswitch applications in higher eukaryotes. However, its general application has been hampered by difficulty in identifying a specific ligand-responsive pseudoknot that also functions as a ligand-dependent -1 PRF stimulator. We addressed this problem by using the -1 PRF stimulation pseudoknot of SARS-CoV (SARS-PK) to build a ligand-dependent -1 PRF stimulator. In particular, the extra stem of SARS-PK was replaced by an RNA aptamer of theophylline and designed to couple theophylline binding with the stimulation of -1 PRF. Conformational and functional analyses indicate that the engineered theophylline-responsive RNA functions as a mammalian riboswitch with robust theophylline-dependent -1 PRF stimulation activity in a stable human 293T cell-line. Thus, RNA-ligand interaction repertoire provided by in vitro selection becomes accessible to ligand-specific -1 PRF stimulator engineering using SARS-PK as the scaffold for synthetic biology application.

Programmed -1 ribosomal frameshifting (-1PRF) is a recoding mechanism to make alternative proteins from a single mRNA transcript. -1PRF is stimulated by cis-acting signals in mRNA, a seven-nucleotide slippery sequence and a downstream secondary structure element, which is often a pseudoknot. In this study we engineered the frameshifting pseudoknot from the mouse mammary tumor virus to respond to a rationally designed small molecule naphthyridine carbamate tetramer (NCTn). We demonstrate that NCTn can stabilize the pseudoknot structure in mRNA and activate -1PRF both in vitro and in human cells. The results illustrate how NCTn-inducible -1PRF may serve as an important component of the synthetic biology toolbox for the precise control of gene expression using small synthetic molecules.

The dicistrovirus intergenic internal ribosome entry site (IGR IRES) directly recruits the ribosome and initiates translation using a non-AUG codon. A subset of IGR IRESs initiates translation in either of two overlapping open reading frames (ORFs), resulting in expression of the 0 frame viral structural polyprotein and an overlapping +1 frame ORFx. A U-G base pair adjacent to the anticodon-like pseudoknot of the IRES directs +1 frame translation. Here, we show that the U-G base pair is not absolutely required for +1 frame translation. Extensive mutagenesis demonstrates that 0 and +1 frame translation can be uncoupled. Ribonucleic acid (RNA) structural probing analyses reveal that the mutant IRESs adopt distinct conformations. Toeprinting analysis suggests that the reading frame is selected at a step downstream of ribosome assembly. We propose a model whereby the IRES adopts conformations to occlude the 0 frame aminoacyl-tRNA thereby allowing delivery of the +1 frame aminoacyl-tRNA to the A site to initiate translation of ORFx. This study provides a new paradigm for programmed recoding mechanisms that increase the coding capacity of a viral genome.

Folding messenger RNA into specific structures is a common regulatory mechanism involved in translation. In Escherichia coli, the operator of the rpsO gene transcript folds into a pseudoknot or double-hairpin conformation. S15, the gene product, binds only to the pseudoknot, thereby repressing its own synthesis when it is present in excess in the cell. The two RNA conformations have been proposed to exist in equilibrium. However, it remained unclear how structural changes can be achieved between these two topologically distinct conformations. We used optical tweezers to study the structural dynamics and rearrangements of the rpsO operator RNA at the single-molecule level. We discovered that the two RNA structures can be interchanged spontaneously and the pseudoknot can exist in conformations that exhibit various levels of stability. Conversion from the double hairpin to a pseudoknot through potential hairpin–hairpin interactions favoured the high-stability conformation. By contrast, mutations that blocked the formation of a hairpin typically resulted in alternative low-stability pseudoknots. These results demonstrate that specific tertiary interactions of RNA can be established and modulated based on the interactions and rearrangements between secondary structural components. Our findings provide new insight into the RNA folding pathway that leads to a regulatory conformation for target protein binding.

Telomerase is a unique reverse transcriptase that maintains the 3′ ends of eukaryotic chromosomes by adding tandem telomeric repeats. The RNA subunit (TR) of vertebrate telomerase provides a template for reverse transcription, contained within the conserved template/pseudoknot domain, and a conserved regions 4 and 5 (CR4/5) domain, all essential for catalytic activity. We report the nuclear magnetic resonance (NMR) solution structure of the full-length CR4/5 domain from the teleost fish medaka (Oryzias latipes). Three helices emanate from a structured internal loop, forming a Y-shaped structure, where helix P6 stacks on P5 and helix P6.1 points away from P6. The relative orientations of the three helices are Mg2+ dependent and dynamic. Although the three-way junction is structured and has unexpected base pairs, telomerase activity assays with nucleotide substitutions and deletions in CR4/5 indicate that none of these are essential for activity. The results suggest that the junction is likely to change conformation in complex with telomerase reverse transcriptase and that it provides a flexible scaffold that allows P6 and P6.1 to correctly fold and interact with telomerase reverse transcriptase.

Riboswitches are structural elements in the 5' untranslated regions of many bacterial messenger RNAs that regulate gene expression in response to changing metabolite concentrations by inhibition of either transcription or translation initiation. The preQ1 (7-aminomethyl-7-deazaguanine) riboswitch family comprises some of the smallest metabolite sensing RNAs found in nature. Once ligand-bound, the transcriptional Bacillus subtilis and translational Thermoanaerobacter tengcongensis preQ1 riboswitch aptamers are structurally similar RNA pseudoknots; yet, prior structural studies have characterized their ligand-free conformations as largely unfolded and folded, respectively. In contrast, through single molecule observation, we now show that, at near-physiological Mg(2+) concentration and pH, both ligand-free aptamers adopt similar pre-folded state ensembles that differ in their ligand-mediated folding. Structure-based Gō-model simulations of the two aptamers suggest that the ligand binds late (Bacillus subtilis) and early (Thermoanaerobacter tengcongensis) relative to pseudoknot folding, leading to the proposal that the principal distinction between the two riboswitches lies in their relative tendencies to fold via mechanisms of conformational selection and induced fit, respectively. These mechanistic insights are put to the test by rationally designing a single nucleotide swap distal from the ligand binding pocket that we find to predictably control the aptamers' pre-folded states and their ligand binding affinities.

Ribosomal subunit biogenesis in eukaryotes is a complex multistep process. Mrd1 is an essential and conserved small (40S) ribosomal subunit synthesis factor that is required for early cleavages in the 35S pre-ribosomal RNA (rRNA). Yeast Mrd1 contains five RNA-binding domains (RBDs), all of which are necessary for optimal function of the protein. Proteomic data showed that Mrd1 is part of the early pre-ribosomal complexes, and deletion of individual RBDs perturbs the pre-ribosomal structure. In vivo ultraviolet cross-linking showed that Mrd1 binds to the pre-rRNA at two sites within the 18S region, in helix 27 (h27) and helix 28. The major binding site lies in h27, and mutational analyses shows that this interaction requires the RBD1-3 region of Mrd1. RBD2 plays the dominant role in h27 binding, but other RBDs also contribute directly. h27 and helix 28 are located close to the sequences that form the central pseudoknot, a key structural feature of the mature 40S subunit. We speculate that the modular structure of Mrd1 coordinates pseudoknot formation with pre-rRNA processing and subunit assembly.

Telomerase emerged during evolution as a prominent solution to the eukaryotic linear chromosome end-replication problem. Telomerase minimally comprises the catalytic telomerase reverse transcriptase (TERT) and telomerase RNA (TR) that provides the template for telomeric DNA synthesis. While the TERT protein is well-conserved across taxa, TR is highly divergent amongst distinct groups of species. Herein, we have identified the essential functional domains of TR from the basal eukaryotic species Trypanosoma brucei, revealing the ancestry of TR comprising two distinct structural core domains that can assemble in trans with TERT and reconstitute active telomerase enzyme in vitro. The upstream essential domain of T. brucei TR, termed the template core, constitutes three short helices in addition to the 11-nt template. Interestingly, the trypanosome template core domain lacks the ubiquitous pseudoknot found in all known TRs, suggesting later evolution of this critical structural element. The template-distal domain is a short stem-loop, termed equivalent CR4/5 (eCR4/5). While functionally similar to vertebrate and fungal CR4/5, trypanosome eCR4/5 is structurally distinctive, lacking the essential P6.1 stem-loop. Our functional study of trypanosome TR core domains suggests that the functional requirement of two discrete structural domains is a common feature of TRs and emerged early in telomerase evolution.

Telomerase is a ribonucleoprotein with an intrinsic telomerase RNA (TER) component. Within yeasts, TER is remarkably large and presents little similarity in secondary structure to vertebrate or ciliate TERs. To better understand the evolution of fungal telomerase, we identified 74 TERs from Pezizomycotina and Taphrinomycotina subphyla, sister clades to budding yeasts. We initially identified TER from Neurospora crassa using a novel deep-sequencing–based approach, and homologous TER sequences from available fungal genome databases by computational searches. Remarkably, TERs from these non-yeast fungi have many attributes in common with vertebrate TERs. Comparative phylogenetic analysis of highly conserved regions within Pezizomycotina TERs revealed two core domains nearly identical in secondary structure to the pseudoknot and CR4/5 within vertebrate TERs. We then analyzed N. crassa and Schizosaccharomyces pombe telomerase reconstituted in vitro, and showed that the two RNA core domains in both systems can reconstitute activity in trans as two separate RNA fragments. Furthermore, the primer-extension pulse-chase analysis affirmed that the reconstituted N. crassa telomerase synthesizes TTAGGG repeats with high processivity, a common attribute of vertebrate telomerase. Overall, this study reveals the common ancestral cores of vertebrate and fungal TERs, and provides insights into the molecular evolution of fungal TER structure and function.

Phosphoryl transfer onto backbone hydroxyls is a recognized catalytic activity of nucleic acids. We find that kinase ribozyme K28 possesses an unusually complex active site that promotes (thio)phosphorylation of two residues widely separated in primary sequence. After allowing the ribozyme to radiolabel itself by phosphoryl transfer from [γ-(32)P]GTP, DNAzyme-mediated cleavage yielded two radiolabeled cleavage fragments, indicating phosphorylation sites within each of the two cleavage fragments. These sites were mapped by alkaline digestion and primer extension pausing. Enzymatic digestion and mutational analysis identified nucleotides important for activity and established the active structure as being a constrained pseudoknot with unusual connectivity that may juxtapose the two reactive sites. Nuclease sensitivities for nucleotides near the pseudoknot core were altered in the presence of GTPγS, indicating donor-induced folding. The 5' target site was more strongly favored in full-length ribozyme K28 (128 nt) than in truncated RNAs (58 nt). Electrophoretic mobilities of self-thiophosphorylated products on organomercurial gels are distinct from the 5' mono-thiophosphorylated product produced by reaction with polynucleotide kinase, potentially indicating simultaneous labeling of both sites within individual RNA strands. Our evidence supports a single, compact structure with local dynamics, rather than global rearrangement, as being responsible for dual-site phosphorylation.

RNA tetraloops can recognize receptors to mediate long-range interactions in stable natural RNAs. In vitro selected GNRA tetraloop/receptor interactions are usually more 'G/C-rich' than their 'A/U-rich' natural counterparts. They are not as widespread in nature despite comparable biophysical and chemical properties. Moreover, while AA, AC and GU dinucleotide platforms occur in natural GAAA/11 nt receptors, the AA platform is somewhat preferred to the others. The apparent preference for 'A/U-rich' GNRA/receptor interactions in nature might stem from an evolutionary adaptation to avoid folding traps at the level of the larger molecular context. To provide evidences in favor of this hypothesis, several riboswitches based on natural and artificial GNRA receptors were investigated in vitro for their ability to prevent inter-molecular GNRA/receptor interactions by trapping the receptor sequence into an alternative intra-molecular pseudoknot. Extent of attenuation determined by native gel-shift assays and co-transcriptional assembly is correlated to the G/C content of the GNRA receptor. Our results shed light on the structural evolution of natural long-range interactions and provide design principles for RNA-based attenuator devices to be used in synthetic biology and RNA nanobiotechnology.

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.

Computational RNA secondary structure prediction approaches differ by the way RNA pseudoknot interactions are handled. For reasons of computational efficiency, most approaches only allow a limited class of pseudoknot interactions or are not considering them at all. Here we present a computational method for RNA secondary structure prediction that is not restricted in terms of pseudoknot complexity. The approach is based on simulating a folding process in a coarse-grained manner by choosing helices based on established energy rules. The steric feasibility of the chosen set of helices is checked during the folding process using a highly coarse-grained 3D model of the RNA structures. Using two data sets of 26 and 241 RNA sequences we find that this approach is competitive compared to the existing RNA secondary structure prediction programs pknotsRG, HotKnots and UnaFold. The key advantages of the new method are that there is no algorithmic restriction in terms of pseudoknot complexity and a test is made for steric feasibility. Availability: The program is available as web server at the site: http://cylofold.abcc.ncifcrf.gov.

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.

The translational activity of the hepatitis C virus (HCV) internal ribosome entry site (IRES) and other HCV-like IRES RNAs depends on structured RNA elements in domains II and III, which serve to recruit the ribosomal 40S subunit, eukaryotic initiation factor (eIF) 3 and the ternary eIF2/Met-tRNA(i)(Met)/GTP complex and subsequently domain II assists subunit joining. Porcine teschovirus-1 talfan (PTV-1) is a member of the Picornaviridae family, with a predicted HCV-like secondary structure, but only stem-loops IIId and IIIe in the 40S-binding domain display significant sequence conservation with the HCV IRES. Here, we use chemical probing to show that interaction sites with the 40S subunit and eIF3 are conserved between HCV and HCV-like IRESs. In addition, we reveal the functional role of a strictly conserved co-variation between a purine-purine mismatch near the pseudoknot (A-A/G) and the loop sequence of domain IIIe (GAU/CA). These nucleotides are involved in a tertiary interaction, which serves to stabilize the pseudoknot structure and correlates with translational efficiency in both the PTV-1 and HCV IRES. Our data demonstrate conservation of functional domains in HCV and HCV-like IRESs including a more complex structure surrounding the pseudoknot than previously assumed.

A bioinformatic covariation analysis of a collection of 119 novel variants of the antigenomic, self-cleaving hepatitis delta virus (HDV) RNA motif supported the formation of all of the Watson-Crick base pairs (bp) of the catalytic centre except the C19-G81 pair located at the bottom of the P2 stem. In fact, a novel Watson-Crick bp between C19 and G80 is suggested by the data. Both chemical and enzymatic probing demonstrated that initially the C19-G81 pair is formed in the ribozyme (Rz), but upon substrate (S) binding and the formation of the P1.1 pseudoknot C19 switches its base-pairing partner from G81 to G80. As a result of this finding, the secondary structure of this ribozyme has been redrawn. The formation of the C19-G80 bp results in a J4/2 junction composed of four nucleotides, similar to that seen in the genomic counterpart, thereby increasing the similarities between these two catalytic RNAs. Additional mutagenesis, cleavage activity and probing experiments yield an original characterization of the structural features involving the residues of the J4/2 junction.

Telomerase maintains the integrity of telomeres, the ends of linear chromosomes, by adding G-rich repeats to their 3'-ends. Telomerase RNA is an integral component of telomerase. It contains a template for the synthesis of the telomeric repeats by the telomerase reverse transcriptase. Although telomerase RNAs of different organisms are very diverse in their sequences, a functional non-template element, a pseudoknot, was predicted in all of them. Pseudoknot elements in human and the budding yeast Kluyveromyces lactis telomerase RNAs contain unusual triple-helical segments with AUU base triples, which are critical for telomerase function. Such base triples in ciliates have not been previously reported. We analyzed the pseudoknot sequences in 28 ciliate species and classified them in six different groups based on the lengths of the stems and loops composing the pseudoknot. Using miniCarlo, a helical parameter-based modeling program, we calculated 3D models for a representative of each morphological group. In all cases, the predicted structure contains at least one AUU base triple in stem 2, except for that of Colpidium colpoda, which contains unconventional GCG and AUA triples. These results suggest that base triples in a pseudoknot element are a conserved feature of all telomerases.

Most ab initio pseudoknot predicting methods provide very few folding scenarios for a given RNA sequence and have low sensitivities. RNA researchers, in many cases, would rather sacrifice the specificity for a much higher sensitivity for pseudoknot detection. In this study, we introduce the Pseudoknot Local Motif Model and Dynamic Partner Sequence Stacking (PLMM_DPSS) algorithm which predicts all PLM model pseudoknots within an RNA sequence in a neighboring-region-interference-free fashion. The PLM model is derived from the existing Pseudobase entries. The innovative DPSS approach calculates the optimally lowest stacking energy between two partner sequences. Combined with the Mfold, PLMM_DPSS can also be used in predicting complicated pseudoknots. The test results of PLMM_DPSS, PKNOTS, iterated loop matching, pknotsRG and HotKnots with Pseudobase sequences have shown that PLMM_DPSS is the most sensitive among the five methods. PLMM_DPSS also provides manageable pseudoknot folding scenarios for further structure determination.

Visualizing RNA secondary structures and pseudoknot structures is essential to bioinformatics systems that deal with RNA structures. However, many bioinformatics systems use heterogeneous data structures and incompatible software components, so integration of software components (including a visualization component) into a system can be hindered by incompatibilities between the components of the system. This paper presents an XML web service and web application program for visualizing RNA secondary structures with pseudoknots. Experimental results show that the PseudoViewer web service and web application are useful for resolving many problems with incompatible software components as well as for visualizing large-scale RNA secondary structures with pseudoknots of any type. The web service and web application are available at http://pseudoviewer.inha.ac.kr/.

Based on the experimentally determined atomic coordinates for RNA helices and the self-avoiding walks of the P (phosphate) and C4 (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.

The glmS ribozyme is a self-cleaving RNA catalyst that resides in the 5'-untranslated region of glmS mRNA in certain bacteria. The ribozyme is specifically activated by glucosamine-6-phosphate (GlcN6P), the metabolic product of the GlmS protein, and is thus proposed to provide a feedback mechanism of riboswitch regulation. Both phylogenetic and biochemical analyses of the glmS ribozyme have established a highly conserved core sequence and secondary structure required for GlcN6P-dependent self-cleavage. However, the high degree of nucleotide conservation offers few clues regarding the higher-order structural organization of the catalytic core. To further investigate core ribozyme structure, minimal 'consensus-type' glmS ribozymes that retain GlcN6P-dependent activity were produced. Mutational analyses of consensus-type glmS ribozymes support a model for core ribozyme folding through a pseudoknot structure formed by the interaction of two highly conserved sequence segments. Moreover, GlcN6P-dependent function is demonstrated for bimolecular constructs in which substrate interaction with the ribozyme is minimally comprised of sequence representing that involved in putative pseudoknot formation. These studies suggest that the glmS ribozyme adopts an intricate multi-strand catalytic core through the formation of a pseudoknot structure, and provide a refined model for further considering GlcN6P interaction and GlcN6P-dependent ribozyme function.

The -1 ribosomal frameshifting requires the existence of an in cis RNA slippery sequence and is promoted by a downstream stimulator RNA. An atypical RNA pseudoknot with an extra stem formed by complementary sequences within loop 2 of an H-type pseudoknot is characterized in the severe acute respiratory syndrome coronavirus (SARS CoV) genome. This pseudoknot can serve as an efficient stimulator for -1 frameshifting in vitro. Mutational analysis of the extra stem suggests frameshift efficiency can be modulated via manipulation of the secondary structure within the loop 2 of an infectious bronchitis virus-type pseudoknot. More importantly, an upstream RNA sequence separated by a linker 5' to the slippery site is also identified to be capable of modulating the -1 frameshift efficiency. RNA sequence containing this attenuation element can downregulate -1 frameshifting promoted by an atypical pseudoknot of SARS CoV and two other pseudoknot stimulators. Furthermore, frameshift efficiency can be reduced to half in the presence of the attenuation signal in vivo. Therefore, this in cis RNA attenuator represents a novel negative determinant of general importance for the regulation of -1 frameshift efficiency, and is thus a potential antiviral target.

mRNA pseudoknots have a stimulatory function in programmed � 1 ribosomal frameshifting (� 1 PRF). Though we previously presented a model for how mRNA pseudoknots might activate the mechanism for � 1 PRF, it did not address the question of the role that they may play in positioning the mRNA relative to the ribosome in this process [E. P. Plant, K. L. M. Jacobs, J. W. Harger, A. Meskauskas, J. L. Jacobs, J. L. Baxter, A. N. Petrov and J. D. Dinman (2003) RNA, 9, 168–174]. A separate ‘torsional restraint’ model suggests that mRNA pseudoknots act to increase the fraction of ribosomes directed to pause with the upstream heptameric slippery site positioned at the ribosome’s A- and P-decoding sites [J. D. Dinman (1995) Yeast, 11, 1115–1127]. Here, experiments using a series of ‘pseudo-pseudoknots’ having different degrees of rotational freedom were used to test this model. The results of this study support the mechanistic hypothesis that � 1 ribosomal frameshifting is enhanced by torsional resistance of the mRNA pseudoknot.

Modular architecture is a hallmark of RNA structures, implying structural, and possibly functional, similarity among existing RNAs. To systematically delineate the existence of smaller topologies within larger structures, we develop and apply an efficient RNA secondary structure comparison algorithm using a newly developed two-dimensional RNA graphical representation. Our survey of similarity among 14 pseudoknots and subtopologies within ribosomal RNAs (rRNAs) uncovers eight pairs of structurally related pseudoknots with non-random sequence matches and reveals modular units in rRNAs. Significantly, three structurally related pseudoknot pairs have functional similarities not previously known: one pair involves the 3′ end of brome mosaic virus genomic RNA (PKB134) and the alternative hammerhead ribozyme pseudoknot (PKB173), both of which are replicase templates for viral RNA replication; the second pair involves structural elements for translation initiation and ribosome recruitment found in the viral internal ribosome entry site (PKB223) and the V4 domain of 18S rRNA (PKB205); the third pair involves 18S rRNA (PKB205) and viral tRNA-like pseudoknot (PKB134), which probably recruits ribosomes via structural mimicry and base complementarity. Additionally, we quantify the modularity of 16S and 23S rRNAs by showing that RNA motifs can be constructed from at least 210 building blocks. Interestingly, we find that the 5S rRNA and two tree modules within 16S and 23S rRNAs have similar topologies and tertiary shapes. These modules can be applied to design novel RNA motifs via build-up-like procedures for constructing sequences and folds.

Oligoribonucleotides that corresponded to the X regions of the (+) and (−) polarity strands of HCV RNA, as well as several shorter oligomers comprising defined stem-loop motifs of their predicted secondary structure models, were analyzed by Pb2+-induced cleavage, partial digestion with specific nucleases and chemical modification. Patterns characteristic of the motifs were compared with those obtained for the full-length molecules and on the basis of such ‘structural fingerprinting’ conclusions concerning folding of regions X were formulated. It turned out that the secondary structure model of X(+) RNA proposed earlier, the three-stem-loop model composed of hairpins SL1, SL2 and SL3, was only partially consistent with our experimental data. We confirmed the presence of SL1 and SL3 motifs and showed that the single-stranded stretch adjacent to the earlier proposed hairpin SL2 contributed to the folding of that region. It seemed to be arranged into two hairpins, which might form a hypothetical pseudoknot by changing their base-pairing systems. These data were discussed in terms of their possible biological significance. On the other hand, analysis of the X(−) RNA and its sub-fragments supported a three-stem-loop secondary structure model for this RNA.

tmRNA combines tRNA- and mRNA-like properties and ameliorates problems arising from stalled ribosomes. Research on the mechanism, structure and biology of tmRNA is served by the tmRNA website (http://www.indiana.edu/~ tmrna), a collection of sequences, alignments, secondary structures and other information. Because many of these sequences are not in GenBank, a BLAST server has been added; another new feature is an abbreviated alignment for the tRNA-like domain only. Many tmRNA sequences from plastids have been added, five found in public sequence data and another 10 generated by direct sequencing; detection in early-branching members of the green plastid lineage brings coverage to all three primary plastid lineages. The new sequences include the shortest known tmRNA sequence. While bacterial tmRNAs usually have a lone pseudoknot upstream of the mRNA segment and a string of three or four pseudoknots downstream, plastid tmRNAs collectively show loss of pseudoknots at both postions. The pseudoknot-string region is also too short to contain the usual pseudoknot number in another new entry, the tmRNA sequence from a bacterial endosymbiont of insect cells, Tremblaya princeps. Pseudoknots may optimize tmRNA function in free-living bacteria, yet become dispensible when the endosymbiotic lifestyle relaxes selective pressure for fast growth.

The abbreviated name, 'mfold web server', describes a number of closely related software applications available on the World Wide Web (WWW) for the prediction of the secondary structure of single stranded nucleic acids. The objective of this web server is to provide easy access to RNA and DNA folding and hybridization software to the scientific community at large. By making use of universally available web GUIs (Graphical User Interfaces), the server circumvents the problem of portability of this software. Detailed output, in the form of structure plots with or without reliability information, single strand frequency plots and 'energy dot plots', are available for the folding of single sequences. A variety of 'bulk' servers give less information, but in a shorter time and for up to hundreds of sequences at once. The portal for the mfold web server is http://www.bioinfo.rpi.edu/applications/mfold. This URL will be referred to as 'MFOLDROOT'.

Given a set of homologous or functionally related RNA sequences, the consensus motifs may represent the binding sites of RNA regulatory proteins. Unlike DNA motifs, RNA motifs are more conserved in structures than in sequences. Knowing the structural motifs can help us gain a deeper insight of the regulation activities. There have been various studies of RNA secondary structure prediction, but most of them are not focused on finding motifs from sets of functionally related sequences. Although recent research shows some new approaches to RNA motif finding, they are limited to finding relatively simple structures, e.g. stem–loops. In this paper, we propose a novel genetic programming approach to RNA secondary structure prediction. It is capable of finding more complex structures than stem–loops. To demonstrate the performance of our new approach as well as to keep the consistency of our comparative study, we first tested it on the same data sets previously used to verify the current prediction systems. To show the flexibility of our new approach, we also tested it on a data set that contains pseudoknot motifs which most current systems cannot identify. A web-based user interface of the prediction system is set up at http://bioinfo. cis.nctu.edu.tw/service/gprm/.

Minimized trans-acting HDV ribozyme systems consisting of three (Rz-3) and two (Rz-2) RNA strands were prepared and their folding conformations were analyzed by NMR spectroscopy. The guanosine residues in one of the enzyme components of Rz-3 were labeled with 13C and 15N. Imino proton signals were assigned by analysis of NOESY and HSQC spectra. The results are consistent with the nested double pseudoknot model, which contains novel base pairs (P1.1), as observed in the crystal structure of a genomic HDV ribozyme. The NOE connectivities suggest an additional G:G pair at the bottom of P1.1 and at the top of P4. The effects of temperature and Mg2+ ions on base pairs for Rz-3 were examined. The temperature variation experiment on Rz-3 showed that P3 is the most stable and that P1.1 is as stable as P1 and P2. The imino proton signals of the G:U pair at the bottom of P1 and the top of P1.1, which are close to the cleavage site, showed the largest changes upon Mg2+ titration of Rz-3. The results suggest that the catalytic Mg2+ ion binds to the pocket formed by P1 and L3.

The assembly of eukaryotic ribosomes requires numerous factors that transiently associate with evolving pre-ribosomal particles. The Pumilio repeat-containing protein Nop9 briefly associates with the 90S pre-ribosome during its co-transcriptional assembly. Here, we show that Nop9 specifically binds an 11-nucleotide sequence of 18S rRNA that forms the 3΄ side of the central pseudoknot and helix 28 in the mature subunit. Crystal structures of Nop9 in the free and RNA-bound states reveal a new type of Pumilio repeat protein with a distinct structure, target sequence and RNA-binding mode. Nop9 contains 10 Pumilio repeats arranged into a U-shaped scaffold. The target RNA is recognized by two stretches of repeats in a bipartite manner, and three central bases are unrecognized as a result of the degeneracy of repeats 6 and 7. Our data suggest that Nop9 regulates the folding of 18S rRNA at early assembly stages of 90S.

The Alu domain of the signal recognition particle (SRP) arrests protein biosynthesis by competition with elongation factor binding on the ribosome. The mammalian Alu domain is a protein–RNA complex, while prokaryotic Alu domains are protein-free with significant extensions of the RNA. Here we report the crystal structure of the complete Alu domain of Bacillus subtilis SRP RNA at 2.5 Å resolution. The bacterial Alu RNA reveals a compact fold, which is stabilized by prokaryote-specific extensions and interactions. In this ‘closed’ conformation, the 5′ and 3′ regions are clamped together by the additional helix 1, the connecting 3-way junction and a novel minor groove interaction, which we term the ‘minor-saddle motif’ (MSM). The 5′ region includes an extended loop–loop pseudoknot made of five consecutive Watson–Crick base pairs. Homology modeling with the human Alu domain in context of the ribosome shows that an additional lobe in the pseudoknot approaches the large subunit, while the absence of protein results in the detachment from the small subunit. Our findings provide the structural basis for purely RNA-driven elongation arrest in prokaryotes, and give insights into the structural adaption of SRP RNA during evolution.