A short motif in Drosophila SECIS Binding Protein 2 provides differential binding affinity to SECIS RNA hairpins

Selenoproteins contain the amino acid selenocysteine which is encoded by a UGA Sec codon. Recoding UGA Sec requires a complex mechanism, comprising the cis-acting SECIS RNA hairpin in the 3′UTR of selenoprotein mRNAs, and trans-acting factors. Among these, the SECIS Binding Protein 2 (SBP2) is central to the mechanism. SBP2 has been so far functionally characterized only in rats and humans. In this work, we report the characterization of the Drosophila melanogaster SBP2 (dSBP2). Despite its shorter length, it retained the same selenoprotein synthesis-promoting capabilities as the mammalian counterpart. However, a major difference resides in the SECIS recognition pattern: while human SBP2 (hSBP2) binds the distinct form 1 and 2 SECIS RNAs with similar affinities, dSBP2 exhibits high affinity toward form 2 only. In addition, we report the identification of a K (lysine)-rich domain in all SBP2s, essential for SECIS and 60S ribosomal subunit binding, differing from the well-characterized L7Ae RNA-binding domain. Swapping only five amino acids between dSBP2 and hSBP2 in the K-rich domain conferred reversed SECIS-binding properties to the proteins, thus unveiling an important sequence for form 1 binding.     

Structural and thermodynamic signatures that define pseudotriloop RNA hairpins

Pseudotriloop (PTL) structures in RNAs have been recognized as essential elements in RNA folding and recognition of proteins. PTL structures are derived from hexaloops by formation of a cross–loop base pair leaving a triloop and 3′ bulged out residue. Despite their common presence and functional importance, insufficient structural and thermodynamic data are available that can be used to predict formation of PTLs from sequence alone. Using NMR spectroscopy and UV-melting data we established factors that contribute to the formation and stability of PTL structures derived from hepatitis B virus and human foamy virus. The NMR data show that, besides the cross–loop base pair, also a 3′ pyrimidine bulge and a G–C loop-closing base pair are primary determinants of PTL formation. By changing the G–C closing base pair into C–G, the PTL switches into a hexaloop. Comparison of these rules with regular triloop hairpins and PTLs from other sources is discussed as well as the conservation of a PTL in human foamy virus and other spumaretroviruses.     

Forced-unfolding and force-quench refolding of RNA hairpins

Nanomanipulation of individual RNA molecules, using laser optical tweezers, has made it possible to infer the major features of their energy landscape. Time-dependent mechanical unfolding trajectories, measured at a constant stretching force (f_S) of simple RNA structures (hairpins and three-helix junctions) sandwiched between RNA/DNA hybrid handles show that they unfold in a reversible all-or-none manner. To provide a molecular interpretation of the experiments we use a general coarse-grained off-lattice Gō-like model, in which each nucleotide is represented using three interaction sites. Using the coarse-grained model we have explored forced-unfolding of RNA hairpin as a function of f_S and the loading rate (r_f). The simulations and theoretical analysis have been done both with and without the handles that are explicitly modeled by semiflexible polymer chains. The mechanisms and timescales for denaturation by temperature jump and mechanical unfolding are vastly different. The directed perturbation of the native state by f_S results in a sequential unfolding of the hairpin starting from their ends, whereas thermal denaturation occurs stochastically. From the dependence of the unfolding rates on r_f and f_S we show that the position of the unfolding transition state is not a constant but moves dramatically as either r_f or f_S is changed. The transition-state movements are interpreted by adopting the Hammond postulate for forced-unfolding. Forced-unfolding simulations of RNA, with handles attached to the two ends, show that the value of the unfolding force increases (especially at high pulling speeds) as the length of the handles increases. The pathways for refolding of RNA from stretched initial conformation, upon quenching f_S to the quench force f_Q, are highly heterogeneous. The refolding times, upon force-quench, are at least an order-of-magnitude greater than those obtained by temperature-quench. The long f_Q-dependent refolding times starting from fully stretched states are analyzed using a model that accounts for the microscopic steps in the rate-limiting step, which involves the trans to gauche transitions of the dihedral angles in the GAAA tetraloop. The simulations with explicit molecular model for the handles show that the dynamics of force-quench refolding is strongly dependent on the interplay of their contour length and persistence length and the RNA persistence length. Using the generality of our results, we also make a number of precise experimentally testable predictions.     

Structural features in TAR RNA of human and simian immunodeficiency viruses: a phylogenetic analysis.

A comparative analysis of TAR RNA structures in human and simian immunodeficiency viruses reveals the conservation of certain structural features despite the divergence in sequence. Both the TAR elements of HIV-1 and SIV-chimpanzee can be folded into relatively simple one-stem hairpin structures. Chemical and RNAase probes were used to analyze the more complex structure of HIV-2 TAR RNA, which folds into a branched hairpin structure. A surprisingly similar RNA conformation can be proposed for SIV-mandrill, despite considerable divergence in nucleotide sequence. A third structural presentation of TAR sequences is seen for SIV-african green monkey. These results are generally consistent with the classification of HIV-SIV viruses in four subgroups based on sequence analyses (both nucleotide- and amino acid-sequences). However, some conserved TAR structures were detected for members of different virus subgroups. It is therefore proposed that RNA structure analysis might provide an additional tool for determining phylogenetic relationships among the HIV-SIV viruses.     

Replacement of RNA hairpins by in vitro selected tetranucleotides.

An in vitro selection method based on the autolytic cleavage of yeast tRNA(Phe) by Pb2+ was applied to obtain tRNA derivatives with the anticodon hairpin replaced by four single-stranded nucleotides. Based on the rates of the site-specific cleavage by Pb2+ and the presence of a specific UV-induced crosslink, certain tetranucleotide sequences allow proper folding of the rest of the tRNA molecule, whereas others do not. One such successful tetramer sequence was also used to replace the acceptor stem of yeast tRNA(Phe) and the anticodon hairpin of E.coli tRNA(Phe) without disrupting folding. These experiments suggest that certain tetramers may be able to replace structurally nonessential hairpins in any RNA.     

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.     

Exploring the complex folding kinetics of RNA hairpins: I. General folding kinetics analysis

Depending on the nucleotide sequence, the temperature, and other conditions, RNA hairpin-folding kinetics can be very complex. The complexity with a wide range of cooperative and noncooperative kinetic behaviors arises from the interplay between the formation of the loops, the disruption of the misfolded states, and the formation of the rate-limiting base stacks. With a rate constant model and a kinetic-cluster theory, we explore the broad landscape for RNA hairpin-folding kinetics. The model is validated through direct tests against several experimental measurements. The general kinetic folding mechanisms and the predicted great variety of folding kinetics are directly applicable and quantitatively testable in experiments. The results from this study suggest that 1), previous experimental findings based on the individual hairpins revealed only a small fraction of much broader and more complex RNA hairpin-folding landscapes; 2), even for structures as simple as hairpins, universal folding timescales and pathways do not exist; and 3), to treat the loop size as the sole factor to determine the hairpin-folding rate is an oversimplification.     

Recurring features of local tertiary structural elements in RNA molecules exemplified by hepatitis D virus RNA

Elements of local tertiary structure in RNA molecules are important in understanding structure-function relationships. The loop E motif, first identified in several eukaryotic RNAs at functional sites which share an exceptional propensity for UV crosslinking between specific bases, was subsequently shown to have a characteristic tertiary structure. Common sequences and secondary structures have allowed other examples of the E-loop motif to be recognized in a number of RNAs at sites of protein binding or other biological function. We would like to know if more elements of local tertiary structure, in addition to the E-loop, can be identified by such common features. The highly structured circular RNA genome of the hepatitis D virus (HDV) provides an ideal test molecule because it has extensive internal structure, a UV-crosslinkable tertiary element, and specific sites for functional interactions with proteins including host PKR. We have now found a UV-crosslinkable element of local tertiary structure in antigenomic HDV RNA which, although differing from the E-loop, has a very similar pattern of sequence and secondary structure to the UV-crosslinkable element found in the genomic strand. Despite the fact that the two structures map close to one another, the sequences comprising them are not the templates for each other. Instead, the template regions for each element are additional sites for potential higher order structure on their respective complementary strands. This wealth of recurring sequences interspersed with base-paired stems provides a context to examine other RNA species for such features and their correlations with biological function.     

Structural features and stability of an RNA triple helix in solution.

A 30 nt RNA with a sequence designed to form an intramolecular triple helix was analyzed by one-and two-dimensional NMR spectroscopy and UV absorption measurements. NMR data show that the RNA contains seven pyrimidine-purine-pyrimidine base triples stabilized by Watson-Crick and Hoogsteen interactions. The temperature dependence of the imino proton resonances, as well as UV absorption data, indicate that the triple helix is highly stable at acidic pH, melting in a single sharp transition centered at 62 degrees C at pH 4.3. The Watson-Crick and Hoogsteen pairings are disrupted simultaneously upon melting. The NMR data are consistent with a structural model where the Watson-Crick paired strands form an A-helix. Results of model building, guided by NMR data, suggest a possible hydrogen bond between the 2' hydroxyl proton of the Hoogsteen strand and a phosphate oxygen of the purine strand. The structural model is discussed in terms of its ability to account for some of the differences in stability reported for RNA and DNA triple helices and provides insight into features that are likely to be important in the design of RNA binding compounds.     

Structural features of small RNA precursors determine Argonaute loading in Caenorhabditis elegans

In C. elegans, DCR-1 is required for the maturation of both short interfering RNAs (siRNAs) and microRNAs (miRNAs), which are subsequently loaded into different Argonaute proteins to mediate silencing via distinct mechanisms. We used in vivo analyses to show that precursors of small RNAs contain structural features that direct the small RNAs into the RNA interference (RNAi) pathway or the miRNA-processing pathway. Nucleotide changes in the pre-let-7 miRNA precursor that make its stem fully complementary cause the resulting small RNA to be recognized as siRNA and induce binding to RDE-1, which leads to RNAi. Mismatches of 1 to 3 nucleotides at various positions in the stem of the precursor restore direction into the miRNA pathway, as the largest portion of such small RNA variants is associated with ALG-1. The Argonaute proteins to which the small RNAs are bound determine the silencing mode, and no functional overlap between RDE-1 and ALG-1 was detected.     

Local base dynamics and local structural features in RNA and DNA duplexes

Local base motion and local structural base information are derived with a simple motional model from site specifically spin-labeled polyribo- and polydeoxyribonucleotides. The model was developed earlier for some nucleic acids and has now been applied to analyze 22 different nucleic acid systems. We conclude that the base motion of the spin-labeled nucleotide in single-stranded RNA, DNA, or non-base-paired bases in duplexes is of the order of 1 ns and that its base mobility decreases by about a factor of 4 upon base pairing. Also, the tether motion of the probe is slower in an RNA than in a DNA duplex.     

Targeting degradation of RNA by RNase H using DNA hairpins

RNase H degradation of two 15 nt RNA target sites was examined in the presence of hairpin DNAs with a 5 nt loop and a 10 bp stem or single-stranded 15 nt DNAs. One target site was a segment of a 79 nt RNA, and the other was part of a 53 nt RNA. Secondary structure predictions indicate that the 53 nt RNA target site is entirely single stranded, while a portion of the 79 nt RNA target site forms an intramolecular duplex. Less RNase H and DNA were needed to cleave the 53 nt RNA target site than the less accessible 79 nt RNA site. The hairpin DNAs had their 5 nt loop and 3' side of the stem fully complementary to the target sites or had sequence changes that produced one to nine mismatched pairs. T(m) values ranged from 57 to 80 degrees C. The stability of the hairpin DNAs relative to the stability of their corresponding RNA-DNA hybrids influenced the extent of RNase H degradation at 37 degrees C. Under the assay conditions employed, the amount of degradation directed by the hairpin DNAs was correlated with their predicted DeltaG(o) (37) of binding to the RNA targets. A DNA hairpin with one mismatch to the target site of the 79 nt RNA did not induce degradation under conditions where fully complementary DNA hairpins produced 50-80% degradation. The in vitro results indicate that DNA hairpins can enhance the stringency of RNase H targeted degradation of the RNA sites.     

Structural constraints identified with covariation analysis in ribosomal RNA

Covariation analysis is used to identify those positions with similar patterns of sequence variation in an alignment of RNA sequences. These constraints on the evolution of two positions are usually associated with a base pair in a helix. While mutual information (MI) has been used to accurately predict an RNA secondary structure and a few of its tertiary interactions, early studies revealed that phylogenetic event counting methods are more sensitive and provide extra confidence in the prediction of base pairs. We developed a novel and powerful phylogenetic events counting method (PEC) for quantifying positional covariation with the Gutell lab's new RNA Comparative Analysis Database (rCAD). The PEC and MI-based methods each identify unique base pairs, and jointly identify many other base pairs. In total, both methods in combination with an N-best and helix-extension strategy identify the maximal number of base pairs. While covariation methods have effectively and accurately predicted RNAs secondary structure, only a few tertiary structure base pairs have been identified. Analysis presented herein and at the Gutell lab's Comparative RNA Web (CRW) Site reveal that the majority of these latter base pairs do not covary with one another. However, covariation analysis does reveal a weaker although significant covariation between sets of nucleotides that are in proximity in the three-dimensional RNA structure. This reveals that covariation analysis identifies other types of structural constraints beyond the two nucleotides that form a base pair.     

A thermodynamic study of unusually stable RNA and DNA hairpins.

About 70% of the RNA tetra-loop sequences identified in ribosomal RNAs from different organisms fall into either (UNCG) or (GNRA) families (where N = A, C, G, or U; and R = A or G). RNA hairpins with these loop sequences form unusually stable tetra-loop structures. We have studied the RNA hairpin GGAC(UUCG)GUCC and several sequence variants to determine the effect of changing the loop sequence and the loop-closing base pair on the thermodynamic stability of (UNCG) tetra-loops. The hairpin GGAG(CUUG)CUCC with the conserved loop G(CUUG)C was also unusually stable. We have determined melting temperatures (Tm), and obtained thermodynamic parameters for DNA hairpins with sequences analogous to stable RNA hairpins with (UNCG), C(GNRA)G, C(GAUA)G, and G(CUUG)C loops. DNA hairpins with (TTCG), (dUdUCG), and related sequences in the loop, unlike their RNA counterparts, did not form unusually stable hairpins. However, DNA hairpins with the consensus loop sequence C(GNRA)G were very stable compared to hairpins with C(TTTT)G or C(AAAA)G loops. The C(GATA)G and G(CTTG)C loops were also extra stable. The relative stabilities of the unusually stable DNA hairpins are similar to those observed for their RNA analogs.     

Muscleblind-like 1 interacts with RNA hairpins in splicing target and pathogenic RNAs

The MBNL and CELF proteins act antagonistically to control the alternative splicing of specific exons during mammalian postnatal development. This process is dysregulated in myotonic dystrophy because MBNL proteins are sequestered by (CUG)_n and (CCUG)_n RNAs expressed from mutant DMPK and ZNF9 genes, respectively. While these observations predict that MBNL proteins have a higher affinity for these pathogenic RNAs versus their normal splicing targets, we demonstrate that MBNL1 possesses comparably high affinities for (CUG)_n and (CAG)_n RNAs as well as a splicing target, Tnnt3. Mapping of a MBNL1-binding site upstream of the Tnnt3 fetal exon indicates that a preferred binding site for this protein is a GC-rich RNA hairpin containing a pyrimidine mismatch. To investigate how pathogenic RNAs sequester MBNL1 in DM1 cells, we used a combination of chemical/enzymatic structure probing and electron microscopy to determine that MBNL1 forms a ring-like structure which binds to the dsCUG helix. While the MBNL1 N-terminal region is required for RNA binding, the C-terminal region mediates homotypic interactions which may stabilize intra- and/or inter-ring interactions. Our results provide a mechanistic basis for dsCUG-induced MBNL1 sequestration and highlight a striking similarity in the binding sites for MBNL proteins on splicing precursor and pathogenic RNAs.     

cDNA cloning,expression pattern and RNA binding analysis of human selenocysteine insertion sequence (SECIS) binding protein 2

Selenocysteine and selenoprotein synthesis require a complex molecular machinery in mammals. Among the key players is the RNA-protein complex formed by the selenocysteine insertion sequence (SECIS) binding protein (SBP2) and the SECIS element, an RNA hairpin in the 3' untranslated regions of selenoprotein messenger RNAs (mRNAs). We have isolated the DNA complementary to mRNA of the human SBP2, enabling us to establish that it differs from a previously reported human SBP2-like protein. Examination of the expression pattern revealed that the human SBP2 protein is encoded by a 4 kb long mRNA that is over-expressed in testis. Compared to the rat SBP2 sequence, the human SBP2 protein displays two highly conserved domains with 92 and 95% amino acid identity, the latter one containing the RNA binding domain. The inter-domain section carries 55% sequence identity, the remainder of the SBP2 sequences showing about 65% identity, values lower than expected for two mammalian proteins. Interestingly, we could show that the binding of human SBP2 to the SECIS RNA is stimulated by the selenoprotein-specialized elongation translation factor mSelB/eEFsec.