Bulge loops used to measure the helical twist of RNA in solution.

Bulge loops are commonly found in helical segments of cellular RNAs. When incorporated into long double-stranded RNAs, they may introduce points of flexibility or permanent bend that can be detected by the altered electrophoretic gel mobility of the RNA. We find that a single An or Un bulge loop near the middle of a long RNA helix significantly retards the RNA during polyacrylamide gel electrophoresis if n greater than or equal to 2. The mobility of an RNA containing two A2 bulges various periodically with the number of base pairs between the bulges. We interpret this to mean that A2 bulges varies periodically with the number of base pairs between the bulges. We interpret this to mean that Z2 bulges form torsionally stiff bends in the helix; the gel mobility reaches a minimum when the total helical twist between the bulges rotates the arms of the molecule into a cis conformation. The gel mobilities are proportional to the predicted end-to-end distance of the RNA if the average RNA helical repeat is 11.8 +/- 0.2 bp/turn and there is no helical twist (3 +/- 9 degrees) associated with the bulge (data obtained in 0.15 M Na+). Other sizes and sequences of bulges have very different effects on RNA helix conformation and flexibility. U2 bulges bend the helix to a much smaller degree than A2 bulges, while longer A or U bulge sequences probably allow bends of 90 degrees or more; all of these may be fairly flexible joints.(ABSTRACT TRUNCATED AT 250 WORDS)     

Investigation of potential RNA bulge stabilizing elements

As a part of our interest in recognition and cleavage of RNA we carried out thermal melting studies with the aim of screening a number of simple oligonucleotide modifications for their potential as modifying elements for RNA bulge stabilizing oligonucleotides. A specific model system from our studies on oligonucleotide-based artificial nuclease (OBAN) systems was chosen and the bulge size was varied from one to five nucleotides. Introduction of single 2'-modified nucleoside moieties (2'-O-methyl, 2'-deoxy and 2'-deoxy-2'-amino) with different conformational preferences adjacent to the bulge revealed that a higher preference for the north conformers gave more stable bulges across the whole range of bulge sizes. Changing a bulge closing a G-U wobble base pair to a G-C pair resulted in the interesting observation that, although the fully complementary complex and small bulges were highly stabilized, there was little difference in the stability of the larger bulges. The wobble base pair even gave a slight stabilization of the 5 nt bulge system. Introduction of a uridine C-5 linker with a single ammonium group was clearly bulge stabilizing (DeltaT(m) + 4.6 to + 5.4 degrees C for the three most stabilized bulges), although with limited selectivity for different bulge sizes since the fully complementary duplex was also stabilized. Introduction of a naphthoyl group on a 2'-aminolinker mostly gave a destabilizing effect, while introduction of a 5-aminoneocuproine moiety on the same linker resulted in stabilization of all bulges, in particular those with two or four unpaired nucleotides (DeltaT(m) + 3.6 and + 2.9 degrees C respectively). The aromatic groups destabilize the fully complementary duplex, resulting in higher selectivity towards stabilization of bulges. A combination of the studied partial element should be suitable for future designs of modified oligonucleotides that, apart from standard base pairing, can also provide additional non-Watson-Crick recognition of RNA.     

Effect of single-residue bulges on RNA double-helical structures: crystallographic database analysis and molecular dynamics simulation studies

Asymmetric bulge loop motifs are widely dispersed in all types of functional RNAs. They are frequently occurring structural motifs in folded RNA structures and appear commonly in pre-microRNA and ribosomes, where they are involved in specific RNA–RNA and RNA–protein interactions. It is therefore necessary to understand such motifs from a structural point of view. We analyzed all available RNA structures and identified quite a few fragments of double helices that contain bulges. We found that these discontinuities often introduce kinks into the double helices, which also affects the stacking overlap between the base pairs across the irregularity. In order to understand the influence of these bulges on stability and flexibility, we carried out molecular dynamics simulations of three different single-residue bulge-containing RNA helices using the CHARMM36 force field. The structural variability at the junctions of RNA bulges is expected to differ from that in continuous double-helical stretches. The structural features of the junction region were observed to vary noticeably depending on the orientation of the bulge residue. When the base of the bulge residue is looped out, the RNA stretch behaves like a standard long A-form RNA double helix, whereas the entire RNA behaves differently when the base of the bulge residue is intercalated between base pairs inside the RNA stem. Such single-base intercalation was found to introduce a permanent kink into the composite double helix, which could be a recognition element for Dicer during the maturation of miRNA.     

A conserved RNA pseudoknot in a putative molecular switch domain of the 3'-untranslated region of coronaviruses is only marginally stable

The 3'-untranslated region (UTR) of the group 2 coronavirus mouse hepatitis virus (MHV) genome contains a predicted bulged stem-loop (designated P0ab), a conserved cis-acting pseudoknot (PK), and a more distal stem-loop (designated P2). Base-pairing to create the pseudoknot-forming stem (P1(pk)) is mutually exclusive with formation of stem P0a at the base of the bulged stem-loop; as a result, the two structures cannot be present simultaneously. Herein, we use thermodynamic methods to evaluate the ability of individual subdomains of the 3' UTR to adopt a pseudoknotted conformation. We find that an RNA capable of forming only the predicted PK (58 nt; 3' nucleotides 241-185) adopts the P2 stem-loop with little evidence for P1(pk) pairing in 0.1 M KCl and the absence of Mg(2+); as Mg(2+) or 1 M KCl is added, a new thermal unfolding transition is induced and assignable to P1(pk) pairing. The P1(pk) helix is only marginally stable, ΔG(25) ≈ 1.2 ± 0.3 kcal/mol (5.0 mM Mg(2+), 100 mM K(+)), and unfolded at 37°C. Similar findings characterize an RNA 5' extended through the P0b helix only (89 nt; 294-185). In contrast, an RNA capable of forming either the P0a helix or the pseudoknot (97 nt; 301-185) forms no P1(pk) helix. Thermal unfolding simulations are fully consistent with these experimental findings. These data reveal that the PK forms weakly and only when the competing double-hairpin structure cannot form; in the UTR RNA, the double hairpin is the predominant conformer under all solution conditions.     

Fingerprinting the junctions of RNA structure by an open-paddlewheel diruthenium compound

RNA function is determined by its structural organization. The RNA structure consists of the combination of distinct secondary structure motifs connected by junctions that play an essential role in RNA folding. Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) probing is an established methodology to analyze the secondary structure of long RNA molecules in solution, which provides accurate data about unpaired nucleotides. However, the residues located at the junctions of RNA structures usually remain undetected. Here we report an RNA probing method based on the use of a novel open-paddlewheel diruthenium (OPW-Ru) compound [Ru₂Cl₂(µ-DPhF)₃(DMSO)] (DPhF = N,N′-diphenylformamidinate). This compound has four potential coordination sites in a singular disposition to establish covalent bonds with substrates. As a proof of concept, we have analyzed the reactivity of OPW-Ru toward RNA using two viral internal ribosome entry site (IRES) elements whose function depends on the structural organization of the molecule. Our study suggests that the compound OPW-Ru preferentially attacks at positions located one or two nucleotides away from junctions or bulges of the RNA structure. The OPW-Ru fingerprinting data differ from that obtained by other chemical reagents and provides new information about RNA structure features.     

Thermodynamics of the pseudo‐knot in helix 18 of 16S ribosomal RNA

A fragment of E. coli 16S rRNA formed by nucleotides 500 to 545 is termed helix 18. Nucleotides 505‐507 and 524‐526 form a pseudo‐knot and its distortion affects ribosome function. Helix 18 isolated from the ribosome context is thus an interesting fragment to investigate the structural properties and folding of RNA with pseudo‐knots. With all‐atom molecular dynamics simulations, spectroscopic and gel electrophoresis experiments, we investigated thermodynamics of helix 18, with a focus on its pseudo‐knot. In solution studies at ambient conditions we observed dimerization of helix 18. We proposed that the loop, containing nucleotides forming the pseudo‐knot, interacts with another monomer of helix 18. The native dimer is difficult to break but introducing mutations in the pseudo‐knot indeed assured a monomeric form of helix 18. Molecular dynamics simulations at 310 K confirmed the stability of the pseudo‐knot but at elevated temperatures this pseudo‐knot was the first part of helix 18 to lose the hydrogen bond pattern. To further determine helix 18 stability, we analyzed the interactions of helix 18 with short oligomers complementary to a nucleotide stretch containing the pseudo‐knot. The formation of higher‐order structures by helix 18 impacts hybridization efficiency of peptide nucleic acid and 2'‐O methyl RNA oligomers.     

DbpA is a region‐specific RNA helicase

DbpA is a DEAD‐box RNA helicase implicated in RNA structural rearrangements in the peptidyl transferase center. DbpA contains an RNA binding domain, responsible for tight binding of DbpA to hairpin 92 of 23S ribosomal RNA, and a RecA‐like catalytic core responsible for double‐helix unwinding. It is not known if DbpA unwinds only the RNA helices that are part of a specific RNA structure, or if DbpA unwinds any RNA helices within the catalytic core's grasp. In other words, it is not known if DbpA is a site‐specific enzyme or region‐specific enzyme. In this study, we used protein and RNA engineering to investigate if DbpA is a region‐specific or a site‐specific enzyme. Our data suggest that DbpA is a region‐specific enzyme. This conclusion has an important implication for the physiological role of DbpA. It suggests that during ribosome assembly, DbpA could bind with its C‐terminal RNA binding domain to hairpin 92, while its catalytic core may unwind any double‐helices in its vicinity. The only requirement for a double‐helix to serve as a DbpA substrate is for the double‐helix to be positioned within the catalytic core's grasp.     

Determinants of a protein-induced RNA switch in the large domain of signal recognition particle identified by systematic-site directed mutagenesis.

Signal recognition particle (SRP) is a ribonucleoprotein complex that associates with ribosomes to promote the co-translational translocation of proteins across biological membranes. Human SRP RNA molecules exist in two distinct conformations, SR-A and SR-B, which may exchange during the assembly of the particle or could play a functional role in the SRP cycle. We have used systematic site-directed mutagenesis of the SRP RNA to determine the electrophoretic mobilities of altered RNA molecules, and we have identified the nucleotides that avert the formation of the two conformers. The conformer behavior of the various RNAs was addressed quantitatively by calculating a value zeta as an indicator of conformational variability. Single loose A-like forms were induced by changes in helix 5 [nucleotides (nt) at positions 111-128 or 222-231], helix 6 (nt at positions 141-150) and helix 7 (nt at position 169 and 170), whereas other mutations in helix 6 and helix 8 preserved the conformational variability of the mutant RNA molecules. The more compact B-like form was induced only when a small region (129-CAAUAU-134), located in the 5'-proximal portion of helix 6, was altered. Since this region is part of the binding site for SRP19, we suggest that protein SRP19 uses nucleotides at 129-134 to trigger the formation of the favored SR-B-form, and thus directs an early step in the hierarchical assembly of the large SRP domain.     

Effects of single-base bulges on intercalator binding to small RNA and DNA hairpins and a ribosomal RNA fragment

The way in which a single-base bulge might affect the structure of an RNA helix has been examined by preparing a series of six RNA hairpins, all with seven base pairs and a four-nucleotide loop. Five of the hairpins have single-base bulges at different positions. The intercalating cleavage reagent (methidiumpropyl)-EDTA-Fe(II) [MPE-Fe(II)] binds preferentially at a CpG sequence in the helix lacking a bulge and in four of the five hairpins with bulges. Hairpins with a bulge one or two bases to the 3' side of the CpG sequence bind ethidium 4-5-fold more strongly than the others. V1 RNase, which is sensitive to RNA backbone conformation in helices, detects a conformational change in all of the helices when ethidium binds; the most dramatic changes, involving the entire hairpin stem, are in one of the two hairpins with enhanced ethidium affinity. Only a slight conformational change is detected in the hairpin lacking a bulge. A bulge adjacent to a CpG sequence in a 100-nucleotide ribosomal RNA fragment enhances MPE-Fe(II) binding by an order of magnitude. These results extend our previous observations of bulges at a single position in an RNA hairpin [White, S. A., & Draper, D.E. (1987) Nucleic Acids Res. 15, 4049] and show that (1) a structural change in an RNA helix may be propagated for several base pairs, (2) bulges tend to increase the number of conformations available to a helix, and (3) the effects observed in small RNA hairpins are relevant to larger RNAs with more extensive structure. A bulge in a DNA hairpin identical in sequence with the RNA hairpins does not enhance MPE-Fe(II) binding affinity, relative to a control DNA hairpin. The effects of bulges on ethidium intercalation are evidently modulated by helix structure.     

Selection on the structural stability of a ribosomal RNA expansion segment in Daphnia obtusa

The high rate of sequence divergence in nuclear ribosomal RNA (rRNA) expansion segments offers a unique opportunity to study the importance of natural selection in their evolution. To this end, we polymerase chain reaction amplified and cloned a 589-nt fragment of the 18S rRNA gene containing expansion segments 43/e1 and 43/e4 from six individual Daphnia obtusa from four populations. We screened 2,588 clones using single-stranded conformation polymorphism analysis and identified 103 unique haplotype sequences. We detected two pairs of indel sites in segment 43/e4 that complement each other when the secondary structure of the linear sequence is formed. Seven of the 12 observed combinations of length variants at these four sites (haplotypes) are shared between individuals from different populations, which may suggest that some of the length variation was present in their common ancestor. Haplotypes with uncompensated indels were only observed at low frequencies, while compensated indel haplotypes were found at a wide range of frequencies, supporting the hypothesis that the energetic stability of expansion segments is a trait under natural selection. In addition, there was strong linkage disequilibrium between the four complementary indel sites, particularly those that pair with one another in the secondary structure. Despite selection against unpaired bulges at these four indel sites, some nucleotides that form unpaired bulges are highly conserved in segment 43/e4, indicating that they are under a different selective constraint, possibly due to their role in higher level structural interactions.     

Bulged-out nucleotides protect an antisense RNA from RNase III cleavage.

Bulged-out nucleotides or internal loops are present in the stem-loop structures of several antisense RNAs. We have used the antisense/target RNA system (CopA/CopT) that controls the copy number of plasmid R1 to examine the possible biological function of bulged-out nucleotides. Two regions within the major stem-loop of the antisense RNA, CopA, carry bulged-out nucleotides. Base pairing in either one or both of these regions of the stem was restored by site-specific mutagenesis and in one case a new internal loop was introduced. The set of mutant and wild-type CopA variants was characterized structurally in vitro. The results reported here indicate a possible function of the bulges: their presence protects CopA RNA from being a substrate for the double-strand-specific enzyme RNase III. In vitro cleavage rates were drastically increased when either the lower or both bulges were absent. This is paralleled by a similar, but not identical, effect of the bulges on metabolic stability of the CopA RNAs in vivo. The degradation pathways of wild-type and mutant CopA in various strain backgrounds are discussed. In the accompanying paper, we address the significance of bulges in CopA for binding to the target RNA in vitro and for its inhibitory efficiency in vivo.     

Different conformational forms of Escherichia coli and rat liver 5S rRNA revealed by Pb(II)-induced hydrolysis.

Different stable forms of Escherichia coli and rat liver 5S rRNA have been probed by Pb(II)-induced hydrolysis. In the native A forms of 5S rRNA, Pb2+ reveal single-stranded RNA stretches and regions of increased conformational flexibility or distorted by the presence of bulged nucleotides. Hydrolysis of urea/EDTA-treated E. coli 5S rRNA (B form) shows the presence of two strong helical domains; helix A retained from the A form and a helix composed of RNA regions G33-C42 and G79-C88. Other RNA regions resistant to hydrolysis may be involved in alternative base pairing, causing conformational heterogeneity of that form. Pb(II)-induced hydrolysis distinguishes two different forms of rat liver 5S rRNA; the native A form and the form obtained by renaturation of 5S rRNA in the presence of EDTA. Pb(II)-hydrolysis data suggest that both forms are highly structured. In the latter form, the orientation of the bulged C66 is changed with respect to helix B. At the same time, a new helical segment is possibly formed, composed of nucleotides from helix C and loop c on one side and from helix E and loop d' on the other.     

Three-dimensional structures of bulge-containing DNA fragments.

The three-dimensional structure of a DNA tridecamer d(CGCAGAATTCGCG)2 containing bulged adenine bases was determined by single crystal X-ray diffraction methods, at 120 K, to 2.6 A resolution. The structure is a B-DNA type double helix with a single duplex in the asymmetric unit. One of the bulged adenine bases loops out from the double helix, while the other stacks in to it. This is in contrast to our preliminary finding, which indicated that both adenine bases were looped out. This revised model was confirmed by the use of a covalently bound heavy-atom derivative. The conformation of the looped-out bulge hardly disrupts base stacking interactions of the bases flanking it. This is achieved by the backbone making a "loop-the-loop" curve with the extra adenine flipping over with respect to the other nucleotides in the strand. The looped-out base intercalates into the stacked-in bulge site of a symmetrically related duplex. The looped-out and stacked-in bases form an A.A reversed Hoogsteen base-pair that stacks between the surrounding base-pairs, thus stabilizing both bulges. The double helix is frayed at one end with the two "melted" bases participating in intermolecular interactions. A related structure, of the same tridecamer, after soaking the crystals with proflavin, was determined to 3.2 A resolution. The main features of this B-DNA duplex are basically similar to the native tridecamer but differ in detail especially in the conformation of the bulged-out base. Accommodation of a large perturbation such as that described here with minimal disruption of the double helix shows both the flexibility and resiliency of the DNA molecule.     

RNA editing sites in plant mitochondria can share cis-elements

RNA editing in flowering plant mitochondria alters numerous C nucleotides in a given mRNA molecule to U residues. To investigate whether neighbouring editing sites can influence each other we analyzed in vitro RNA editing of two sites spaced 30 nt apart. Deletion and competition experiments show that these two sites carry independent essential specificity determinants in the respective upstream 20-30 nucleotides. However, deletion of a an upstream sequence region promoting editing of the upstream site concomitantly decreases RNA editing of the second site 50-70 nucleotides downstream. This result suggests that supporting cis-/trans-interactions can be effective over larger distances and can affect more than one editing event.     

G-ribo: a new structural motif in ribosomal RNA

Analysis of the available crystal structures of the ribosome and of its subunits has revealed a new RNA motif that we call G-ribo. The motif consists of two double helices positioned side-by-side and connected by an unpaired region. The juxtaposition of the two helices is kept by a complex system of tertiary interactions spread over several layers of stacked nucleotides. In the center of this arrangement, the ribose of a nucleotide from one helix is specifically packed with the ribose and the minor-groove edge of a guanosine from the other helix. In total, we found eight G-ribo motifs in both ribosomal subunits. The location of these motifs suggests that at least some of them play an important role in the formation of the ribosome structure and/or in its function.     

On loop folding in nucleic acid hairpin-type structures

In a series of studies, combining NMR, optical melting and T-jump experiments, it was found that DNA hairpins display a maximum stability when the loop part of the molecule comprises four or five nucleotide residues. This is in contrast with the current notion based on RNA hairpin studies, from which it had been established that a maximum hairpin stability is obtained for six or seven residues in the loop. Here we present a structural model to rationalize these observations. This model is based on the notion that to a major extent base stacking interactions determine the stability of nucleic acid conformations. The model predicts that loop folding in RNA is characterized by an extension of the base stacking at the 5'-side of the double helix by five or six bases; the remaining gap can then easily be closed by two nucleotides. Conversely, loop folding in DNA is characterized by extending base stacking at the 3'-side of the double helical stem by two or three residues; again bridging of the remaining gap can then be achieved by one or two nucleotides. As an example of loop folding in RNA the anticodon loop of yeast tRNAPhe is discussed. For the DNA hairpin formed by d(ATCCTAT4TAGGAT) it is shown that the loop structure obtained from molecular mechanics calculations obeys the above worded loop folding principles.     

基于堆积协变信息与最小自由能预测含伪结的RNA二级结构

RNA伪结预测是RNA研究的一个难点问题。文中提出一种基于堆积协变信息与最小自由能的RNA伪结预测方法。该方法使用已知结构的RNA比对序列(ClustalW比对和结构比对)测试此方法, 侧重考虑相邻碱基对之间相互作用形成的堆积协变信息, 并结合最小自由能方法对碱基配对综合评分, 通过逐步迭代求得含伪结的RNA二级结构。结果表明, 此方法能正确预测伪结, 其平均敏感性和特异性优于参考算法, 并且结构比对的预测性能比ClustalW比对的预测性能更加稳定。文中同时讨论了不同协变信息权重因子对预测性能的影响, 发现权重因子比值在l1: l2=5:1时, 预测性能达到最优。     

Nonenzymatic Template-directed Synthesis of RNA from Monomers

Reviewed are the latest achievements in studying the information transfer mechanisms and the evolutionary significance of prebiotic RNA synthesis, the double helix structures most preferred in this respect, and the possible reasons for the prevalence of particular enantiomeric forms of nucleotides in template-directed synthesis.     

Characterization of an RNA bulge structure by Fourier transform infrared spectroscopy

There may be several advantages associated with an antisense oligonucleotide that induces a bulged structure into its RNA target molecule. Many structures of RNA bulges are elucidated from single-stranded RNA models. However, a two-component system is the minimum requirement for a realistic antisense model. We have used Fourier transform infrared spectroscopy to investigate a single-stranded RNA oligonucleotide with known NMR solution structure, constructed to model a five nucleotide bulge, and its two-component oligonucleotide counterpart. The infrared spectra show A-helical base-paired stems and non-base-paired loops in both systems. The nucleosides are mainly in an anti-conformation. Both N-type and S-type of sugar puckers can be inferred from the infrared region sensitive to sugar conformations. The S-type of sugar pucker is likely to be associated with the nucleotides in the bulge. The FTIR results display an overall structural similarity between the two model systems.