Thermodynamics of RNA hairpins containing single internal mismatches.

Thermodynamic parameters and circular dichroism spectra are presented for RNA hairpins containing single internal mismatches in the stem regions. Three different sequence contexts for the G*U mismatch and two contexts for C*A, G*A, U*U, A*C and U*G mismatches were examined and compared with Watson-Crick base-pair stabilities. The RNA hairpins employed were a microhelix and tetraloop representing the Escherichia coli tRNAAlaacceptor stem and sequence variants that have been altered at the naturally occurring G*U mismatch site. UV melting studies were carried out under different conditions to evaluate the effects of sodium ion concentration and pH on the stability of mismatch-containing hairpins. Our main findings are that single internal mismatches exhibit a range of effects on hairpin stability. In these studies, the size and sequence of the loop and stem are shown to influence the overall stability of the RNA, and have a minor effect on the relative mismatch stabilities. The relationship of these results to RNA-ligand interactions involving mismatch base-pairs is discussed.     

Architecture of the U5 small nuclear RNA.

We have used comparative sequence analysis and deletion analysis to examine the secondary structure of the U5 small nuclear RNA (snRNA), an essential component of the pre-mRNA splicing apparatus. The secondary structure of Saccharomyces cerevisiae U5 snRNA was studied in detail, while sequences from six other fungal species were included in the phylogenetic analysis. Our results indicate that fungal U5 snRNAs, like their counterparts from other taxa, can be folded into a secondary structure characterized by a highly conserved stem-loop (stem-loop 1) that is flanked by a moderately conserved internal loop (internal loop 1). In addition, several of the fungal U5 snRNAs include a novel stem-loop structure (ca. 30 nucleotides) that is adjacent to stem-loop 1. By deletion analysis of the S. cerevisiae snRNA, we have demonstrated that the minimal U5 snRNA that can complement the lethal phenotype of a U5 gene disruption consists of (i) stem-loop 1, (ii) internal loop 1, (iii) a stem-closing internal loop 1, and (iv) the conserved Sm protein binding site. Remarkably, all essential, U5-specific primary sequence elements are encoded by a 39-nucleotide domain consisting of stem-loop 1 and internal loop 1. This domain must, therefore, contain all U5-specific sequences that are essential for splicing activity, including binding sites for U5-specific proteins.     

Perturbation of DNA hairpins containing the EcoRI recognition site by hairpin loops of varying size and composition: physical (NMR and UV) and enzymatic (EcoRI) studies.

We have investigated loop-induced structural perturbation of the stem structure in hairpins d(GAATTCXnGAATTC) (X = A, T and n = 3, 4, 5 and 6) that contain an EcoRI restriction site in close proximity to the hairpin loop. Oligonucleotides containing either a T3 or a A3 loop were not hydrolyzed by the restriction enzyme and also showed only weak binding to EcoRI in the absence of the cofactor Mg2+. In contrast, hairpins with larger loops are hydrolyzed by the enzyme at the scission site next to the loop although the substrate with a A4 loop is significantly more resistant than the oligonucleotide containing a T4 loop. The hairpin structures with 3 loop residues were found to be thermally most stable while larger hairpin loops resulted in structures with lower melting temperatures. The T-loop hairpins are thermally more stable than the hairpins containing the same number of A residues in the loop. As judged from proton NMR spectroscopy and the thermodynamic data, the base pair closest to the hairpin loop did form in all cases studied. The hairpin loops did, however, affect the conformation of the stem structure of the hairpins. From 31P and 1H NMR spectroscopy we conclude that the perturbation of the stem structure is stronger for smaller hairpin loops and that the extent of the perturbation is limited to 2-3 base pairs for hairpins with T3 or A4 loops. Our results demonstrate that hairpin loops modulate the conformation of the stem residues close to the loop and that this in turn reduces the substrate activity for DNA sequence specific proteins.     

A novel class of nucleolar RNAs from Tetrahymena.

We describe a family of at least four nucleolar RNAs (snoRNAs) from the ciliate, Tetrahymena. The snoRNAs are 120-140 nucleotides long, moderately AU-rich and contain no modified nucleotides. Their 5' ends are blocked by a cap of unknown nature. The snoRNAs can be folded into similar secondary structures consisting of two hairpins separated by a single-stranded AU-rich spacer. The sequences and secondary structures show no extensive sequence or secondary structure resemblance to any other small RNAs in the public databases.     

Tethering of human Ago proteins to mRNA mimics the miRNA-mediated repression of protein synthesis

MicroRNAs (miRNAs) are approximately 21-nt-long RNAs involved in regulating development, differentiation, and other processes in eukaryotes. In metazoa, nearly all miRNAs control gene expression by imperfectly base-pairing with the 3'-untranslated region (3'-UTR) of target mRNAs and repressing protein synthesis by an unknown mechanism. It is also unknown whether miRNA-mRNA duplexes containing mismatches and bulges provide specific features that are recognized by factors mediating the repression. miRNAs form part of ribonucleoprotein complexes, miRNPs, that contain Argonaute (Ago) and other proteins. Here we demonstrate that effects of miRNAs on translation can be mimicked in human HeLa cells by the miRNA-independent tethering of Ago proteins to the 3'-UTR of a reporter mRNA. Inhibition of protein synthesis occurred without a change in the reporter mRNA level and was dependent on the number, but not the position, of the hairpins tethering hAgo2 to the 3'-UTR. These findings indicate that a primary function of miRNAs is to guide their associated proteins to the mRNA.     

RNA pseudoknots. Stability and loop size requirements.

The effects of ionic conditions, loop size and loop sequence on the formation of pseudoknots by RNA oligonucleotides have been investigated using biochemical and biophysical methods. An oligonucleotide with the sequence 5' GCGAUUUCUGACCGCUUUUUUGUCAG 3' and oligonucleotides with variations in the sequences of the two loop regions, denoted by bold face type, were studied. Each sequence with the potential to form a pseudoknot can also form two stable hairpins. The pseudoknot structure is stabilized relative to the hairpins by addition of Mg2+. Even in the presence of Mg2+, the pseudoknots formed by the sequences investigated are only marginally more stable (1.5 to 2 kcal mol-1 in free energy at 37 degrees C) than either of the constituent hairpins. The pseudoknot structure is the stable conformation in the presence of Mg2+ when the first loop region is at least three nucleotides and the second is at least four nucleotides. Further deletion of nucleotides from the loop regions stabilizes possible hairpin structures relative to the pseudoknot and equilibria among secondary and tertiary structures result.     

Homopurine and homopyrimidine strands complementary in parallel orientation form an antiparallel duplex at neutral pH with A-C,G-T,and T-C mismatched base pairs

DNA sequences d-TGAGGAAAGAAGGT (a 14-mer) and d-CTCCTTTCTTCC (a 12-mer) are complementary in parallel orientation forming either Donahue (reverse Watson-Crick) base pairing at neutral pH or Hoogsteen base pairing at slightly acidic pH. The structure of the complex formed by dissolving the two strands in equimolar ratio in water has been investigated by nmr. At neutral pH, the system forms an ordered antiparallel duplex with five A : T and four G : C Watson-Crick base pairs and three mismatches, namely G-T, A-C, and T-C. The nuclear Overhauser effect cross-peak pattern suggests an overall B-DNA conformation with major structural perturbations near the mismatches. The duplex has a low melting point and dissociates directly into single strands with a broad melting profile. The hydrogen-bonding schemes in the mismatched base pairs have been investigated. It has been shown earlier that in acidic pH, the system prefers a triple-stranded structure with two pyrimidine strands and one purine strand. One of the pyrimidine strands has protonated cytosines, forms Hoogsteen base pairing, and is aligned parallel to the purine strand; the other has nonprotonated cytosines and has base-pairing scheme similar to the one discussed in this paper. The parallel duplex is therefore less stable than either the antiparallel duplex or the triplex, in spite of its perfect complementarity. © 1997 John Wiley & Sons, Inc. Biopoly 41: 773–784, 1997     

Diversity of DIS, SD and psi hairpins in HIV-1 isolates of group M: in silico study

The primary sequence and secondary structure of the region encompassing DIS, SD and psi hairpins in HIV-1 genomic RNAs have been analyzed for 731 group M isolates from NCBI database. The secondary structures have been predicted by the m fold program (M. Zuker). Though the primary sequence of the region studied was found to be highly heterogeneous, this region is folded into DIS, SD and psi hairpins (DIS-, SD- and psi-like hairpins) in 96% of the isolates studied. The phylogenetic analysis showed that the most frequent variants of DIS hairpin (DIS(Lai), DIS(Mal) and DIS(C)) tolerate certain base changes. Particularly, base changes at stem position 23 occur 5 and 33 times more frequently in DIS(Lai) than in DIS(Mal) and DIS(C), respectively, while A insertion at the 5'end of apical loop is tolerated in DIS(Mal) and DIS(C) but not in DIS(Lai). We have revealed that the bottom base pair substitution G-C --> A-U in SD hairpin is highly specific for subtype D isolates. All variants of DIS, SD and psi hairpins found in our database are discussed, systematized and presented in schemes of hypothetical transitions between variants via a single base change. Most variants of DIS and psi hairpins were found to adopt several conformations.     

The secondary structure of mRNAs from Escherichia coli: its possible role in increasing the accuracy of translation.

A secondary structure model was proposed for mRNAs during translation (in a polysome) where the secondary structure is described by a set of small unbranched hairpins. Computer simulation experiments reveal that the number of hairpins is much greater (P less than 10(-6) in highly expressed mRNAs from E. coli as compared with the random sequences coding for the same amino acid sequence, i.e. certain synonymous codons are used in definite mRNA positions to increase the number of hairpins. No constraints on the amino acid sequence, which would affect the secondary structure of mRNAs, were found. The codons UGU, UGC (Cys), GCC (Ala), ACA, ACG (Thr), CCU, CCC (Pro), etc. translated by minor tRNAs were found to occur significantly more frequently in the position 5' to the hairpins than the other codons translated by major tRNAs (P less than 5.10(-6). This correlation leads to the hypothesis that the process of hairpin unfolding can increase the time of translocation from the A to P ribosome site of the codon 5' to the hairpin, thus decreasing the probability of translational error (the latter would likely occur more frequently in the codons translated by minor tRNAs).     

Isolation and characterization of thermodynamically stable and unstable RNA hairpins from a triloop combinatorial library

Hairpins are the most common elements of RNA secondary structure, playing important roles in RNA tertiary architecture and forming protein binding sites.Triloops are common in a variety of naturally occurring RNA hairpins, but little is known about their thermodynamic stability. Reported here are the sequences and thermodynamic parameters for a variety of stable and unstable triloop hairpins. Temperature gradient gel electrophoresis (TGGE) can be used to separate a simple RNA combinatorial library based on thermal stability [Bevilacqua, J. M., and Bevilacqua, P. C. (1998) Biochemistry 45, 15877-15884]. Here we introduce the application of TGGE to separating and analyzing a complex RNA combinatorial library based on thermal stability, using an RNA triloop library. Several rounds of in vitro selection of an RNA triloop library were carried out using TGGE, and preferences for exceptionally stable and unstable closing base pairs and loop sequences were identified. For stable hairpins, the most common closing base pair is CG, and U-rich loop sequences are preferred. Closing base pairs of GC and UA result in moderately stable hairpins when combined with a stable loop sequence. For unstable hairpins, the most common closing base pairs are AU and UG, and U-rich loop sequences are no longer preferred. In general, the contributions of the closing base pair and loop sequence to overall hairpin stability appear to be additive. Thermodynamic parameters for individual hairpins determined by UV melting are generally consistent with outcomes from selection experiments, with hairpins containing a CG closing base pair having a DeltaDeltaG degrees (37) 2.1-2.5 kcal/mol more favorable than hairpins with other closing base pairs. Sequences and thermodynamic rules for triloop hairpins should aid in RNA structure prediction and determination of whether naturally occurring triloop hairpins are thermodynamically stable.     

5' and 3' untranslated regions of pestivirus genome: primary and secondary structure analyses.

Within the conserved 5' untranslated region (UTR) of the pestivirus genome three highly variable regions were identified. Preceding the polyprotein start codon, multiple cryptic AUG codons and several small open reading frames are characteristic for all the five pestiviruses. Inspection of the context of AUGs revealed that the polyprotein initiation AUG of pestivirus has a weak context for efficient translation initiation. The most favorable context was found in two of the cryptic AUGs. Two oligopyrimidine-rich tracts upstream to the conserved either cryptic or authentic AUG in the 5'-UTR of pestivirus were identified and 83.3% of their nucleotide sequences are complementary to the consensus sequence at the 3' terminus of eucaryotic 18S rRNA. A secondary structure model for the 5'-UTR of pestivirus was predicted. Nucleotide sequence comparison among five pestiviruses led to the identification of a variable region and a conserved region in the 3'-UTR. A deletion of 41 nucleotides was found within the variable region in Osloss. A secondary structure model for the 3'-UTR was also predicted. The structural similarity of the 5'-UTR between pestiviruses and picornaviruses and hepatitis C viruses was demonstrated and the possible implications of features of the 5' and 3'-UTR of pestiviruses are discussed.     

The mismatched nucleotides encoded in vaccinia virus flip-and-flop hairpin telomeres serve an essential role in virion maturation

Poxvirus genomes consist of a linear duplex DNA that ends in short inverted and complementary hairpin structures. These elements also encode loops and mismatches that likely serve a role in genome packaging and perhaps replication. We constructed mutant vaccinia viruses (VACV) where the native hairpins were replaced by altered forms and tested effects on replication, assembly, and virulence. Our studies showed that structure, not sequence, likely determines function as one can replace an Orthopoxvirus (VACV) hairpin with one copied from a Leporipoxvirus with no effect on growth. Some loops can be deleted from VACV hairpins with little effect, but VACV bearing too few mismatches grew poorly and we couldn’t recover viruses lacking all mismatches. Further studies were conducted using a mutant bearing only one of six mismatches found in wild-type hairpins (SΔ1Δ3–6). This virus grew to ~20-fold lower titers, but neither DNA synthesis nor telomere resolution was affected. However, the mutant exhibited a particle-to-PFU ratio 10-20-fold higher than wild-type viruses and p4b/4b core protein processing was compromised, indicating an assembly defect. Electron microscopy showed that SΔ1Δ3–6 mutant development was blocked at the immature virus (IV) stage, which phenocopies known effects of I1L mutants. Competitive DNA binding assays showed that recombinant I1 protein had less affinity for the SΔ1Δ3–6 hairpin than the wild-type hairpin. The SΔ1Δ3–6 mutant was also attenuated when administered to SCID-NCR mice by tail scarification. Mice inoculated with viruses bearing wild-type hairpins exhibited a median survival of 30–37 days, while mice infected with SΔ1Δ3–6 virus survived >70 days. Persistent infections favor genetic reversion and genome sequencing detected one example where a small duplication near the hairpin tip likely created a new loop. These observations show that mismatches serve a critical role in genome packaging and provide new insights into how VACV “flip and flop” telomeres are arranged.     

Distinct RNA structural domains cooperate to maintain a specific cleavage site in the 3'-UTR of IGF-II mRNAs

The insulin-like growth factor II mRNAs are targets for site-specific endonucleolytic cleavage in the 3'-UTR, which results in a very stable 3' cleavage product of 1.8 kb, consisting of 3'-UTR sequences and a poly(A) tail. The 5' cleavage product contains the coding region and is rapidly degraded. Thus, cleavage is thought to provide an additional way to control IGF-II protein synthesis. We had established that cleavage requires two widely separated sequence elements (I and II) in the 3'-UTR that form a stable duplex of 83 nucleotides. The cleavage-site itself is located in an internal loop preceded by two stable stem-loop structures. Furthermore, in a study which was based on RNA folding algorithms, we have shown that there are specific sequence and structural requirements for the cleavage reaction. Here, the functions of the different structural domains in cleavage were assessed by deletion/mutational analyses, and biochemical structure probing assays were performed to characterize better the RNA structures formed and to verify the computer folding predictions. The data suggest that the stem-loop domain contributes to maintain a highly specific c leavage-site by preventing the formation of alternative structures in the cleavage-site domain. Involvement of the nucleotides in the cleavage-site loop itself in non-Watson-Crick interactions may be important for providing a specific recognition surface for an endoribonuclease activity.     

RNA apatamers for yeast ribosomal protein L32 have a conserved purine-rich internal loop.

Two in vitro selection experiments were conducted to determine the RNA sequence requirements for binding ribosomal protein L32 (RPL32) from Saccharomyces cerevisiae. To preserve the wild-type stem-internal loop-stem fold, only a limited portion of the RNA comprising the internal loop region was randomized. Most of the selected RNAs have secondary structures similar to that of the wild-type, and four purines on both sides of the internal loop are highly conserved. Indeed, a pair of 5'-GA-3' dinucleotides is found in all but one of the stem-loop-stem L32 aptamers and these conserved purines may contact the protein directly or form a necessary RNA secondary or tertiary structure. These aptamers have a potential G:U pair bordering the loop adjacent to the conserved GAs, but a cytidine replaces a phylogenetically conserved adenosine at one loop position in many of the selected RNAs. In model RNAs, the cytidine-bearing variant binds protein slightly more strongly than does the wild-type RNA. That the seven-member, 2 + 5 internal loop is important for protein binding is reinforced by the finding that the position, but not the size, of the loop is variable. A minority of the RNA aptamers has three consecutive uridines and may fold into a similar structure, but with the internal loop inverted.     

猪繁殖与呼吸障碍综合征病毒(PRRSV)3′末端非翻译区中调控序列的研究

以北美株PRRSV感染性克隆pCBC2为平台进行反向遗传操作,将3′UTR中的一级结构进行了系列缺失或插入突变,并改变二级结构中的一个保守的茎环结构,构建全长PRRSV突变体克隆,解析3′UTR突变对病毒感染性的影响,旨在界定调控PRRSV3′UTR的启动子序列及二级结构,即复制过程中的最小调控元件。以空斑和Northern blot来研究拯救后重组病毒的复制、转录和生长特性,发现重组病毒感染动力学与亲本病毒无可见差别。结果表明PRRSV3′UTR的5′端可耐受一定数目的核苷酸的缺失(41nt)与插入(23nt)突变,但进一步9nt缺失造成保守的环结构突变后就使病毒失去了感染性。证明了这是3′UTR中控制PRRSV复制过程的的必需序列及二级结构,为进一步解析PRRSV复制过程的调控元件奠定了基础。     

Eukaryote specific RNA and protein features facilitate assembly and catalysis of H/ACA snoRNPs

H/ACA Box ribonucleoprotein complexes (RNPs) play a major role in modification of rRNA and snRNA, catalyzing the sequence specific pseudouridylation in eukaryotes and archaea. This enzymatic reaction takes place on a substrate RNA recruited via base pairing to an internal loop of the snoRNA. Eukaryotic snoRNPs contain the four proteins Nop10, Cbf5, Gar1 and Nhp2, with Cbf5 as the catalytic subunit. In contrast to archaeal H/ACA RNPs, eukaryotic snoRNPs contain several conserved features in both the snoRNA as well as the protein components. Here, we reconstituted the eukaryotic H/ACA RNP containing snR81 as a guide RNA in vitro and report on the effects of these eukaryote specific features on complex assembly and enzymatic activity. We compare their contribution to pseudouridylation activity for stand-alone hairpins versus the bipartite RNP. Using single molecule FRET spectroscopy, we investigated the role of the different eukaryote-specific proteins and domains on RNA folding and complex assembly, and assessed binding of substrate RNA to the RNP. Interestingly, we found diverging effects for the two hairpins of snR81, suggesting hairpin-specific requirements for folding and RNP formation. Our results for the first time allow assessing interactions between the individual hairpin RNPs in the context of the full, bipartite snoRNP.     

Differential Targeting of Unpaired Bases within Duplex DNA by the Natural Compound Clerocidin: A Valuable Tool to Dissect DNA Secondary Structure

Non-canonical DNA structures have been postulated to mediate protein-nucleic acid interactions and to function as intermediates in the generation of frame-shift mutations when errors in DNA replication occur, which result in a variety of diseases and cancers. Compounds capable of binding to non-canonical DNA conformations may thus have significant diagnostic and therapeutic potential. Clerocidin is a natural diterpenoid which has been shown to selectively react with single-stranded bases without targeting the double helix. Here we performed a comprehensive analysis on several non-canonical DNA secondary structures, namely mismatches, nicks, bulges, hairpins, with sequence variations in both the single-stranded region and the double-stranded flanking segment. By analysis of clerocidin reactivity, we were able to identify the exposed reactive residues which provided information on both the secondary structure and the accessibility of the non-paired sites. Mismatches longer than 1 base were necessary to be reached by clerocidin reactive groups, while 1-base nicks were promptly targeted by clerocidin; in hairpins, clerocidin reactivity increased with the length of the hairpin loop, while, interestingly, reactivity towards bulges reached a maximum in 3-base-long bulges and declined in longer bulges. Electrophoretic mobility shift analysis demonstrated that bulges longer than 3 bases (i.e. 5- and 7-bases) folded or stacked on the duplex region therefore being less accessible by the compound. Clerocidin thus represents a new valuable diagnostic tool to dissect DNA secondary structures.