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.     

Use of ultra stable UNCG tetraloop hairpins to fold RNA structures: thermodynamic and spectroscopic applications.

RNA molecules of > 20 nucleotides have been the focus of numerous recent NMR structural studies. Several investigators have used the UNCG family of hairpins to ensure proper folding. We show that th UUCG hairpin has a minimum requirement of a two base-pair stem. Hairpins with a CG loop closing base pair and an initial 5'CG or 5'GC base pair have a melting temperature approximately 55 degrees C in 10 mM sodium phosphate. The high stability of even such small hairpins suggests that the hairpin can serve as a nucleation site for folding. For high resolution NMR work, the UNCG loop family (UACG in particular) provides excellent spectroscopic markers in one-dimensional exchangeable spectra, in two-dimensional COSY spectra and in NOESY spectra that clearly define it as forming a hairpin. This allows straightforward initiation of chemical shift assignments.     

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.     

Molecular evolution of transfer RNA from two precursor hairpins: Implications for the origin of protein synthesis

In this paper we are going to present a model for the coevolution of major components of the protein synthesis machinery in a primordial RNA world. We propose that the essential prerequisites for RNA-based protein synthesis, i.e., tRNA-like molecules, ribozymic charging catalysts, small-subunit(SSU) rRNA, and large-subunit(LSU) rRNA, evolved from the same ancestral RNA molecule. Several arguments are considered which suggest that tRNA-like molecules were derived by tandem joining of template-flanking hairpin structures involved in replication control. It is further argued that the ancestors of contemporary group I tRNA introns catalyzed such hairpin joining reactions, themselves also giving rise to the ribosomal RNAs. Our model includes a general stereochemical principle for the interaction between ribozymes and hairpin-derived recognition structures, which can be applied to such seemingly different processes as RNA polymerization, aminoacylation, tRNA decoding, and peptidyl transfer, implicating a common origin for these fundamental functions. These and other considerations suggest that generation and evolution of tRNA were coupled to the evolution of synthetases, ribosomal RNAs, and introns from the beginning and have been a consequence arising from the original function of tRNA precursor hairpins as replication and recombination control elements. Correspondence to: T.P. Dick     

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.     

Investigating the structural basis of purine specificity in the structures of MS2 coat protein RNA translational operator hairpins

We have determined the structures of complexes between the phage MS2 coat protein and variants of the replicase translational operator in order to explore the sequence specificity of the RNA–protein interaction. The 19-nt RNA hairpins studied have substitutions at two positions that have been shown to be important for specific binding. At one of these positions, –10, which is a bulged adenosine (A) in the stem of the wild-type operator hairpin, substitutions were made with guanosine (G), cytidine (C) and two non-native bases, 2-aminopurine (2AP) and inosine (I). At the other position, –7 in the hairpin loop, the native adenine was substituted with a cytidine. Of these, only the G-10, C-10 and C-7 variants showed interpretable density for the RNA hairpin. In spite of large differences in binding affinities, the structures of the variant complexes are very similar to the wild-type operator complex. For G-10 substitutions in hairpin variants that can form bulges at alternative places in the stem, the binding affinity is low and a partly disordered conformation is seen in the electron density maps. The affinity is similar to that of wild-type when the base pairs adjacent to the bulged nucleotide are selected to avoid alternative conformations. Both purines bind in a very similar way in a pocket in the protein. In the C-10 variant, which has very low affinity, the cytidine is partly inserted in the protein pocket rather than intercalated in the RNA stem. Substitution of the wild-type adenosine at position –7 by pyrimidines gives strongly reduced affinities, but the structure of the C-7 complex shows that the base occupies the same position as the A-7 in the wild-type RNA. It is stacked in the RNA and makes no direct contact with the protein.     

Thermal unfolding of triosephosphate isomerase from Entamoeba histolytica: dimer dissociation leads to extensive unfolding

In mesophiles, triosephosphate isomerase (TIM) is an obligated homodimer. We have previously shown that monomeric folding intermediates are common in the chemical unfolding of TIM, where dissociation provides 75% of the overall conformational stability of the dimer. However, analysis of the crystallographic structure shows that, during unfolding, intermonomeric contacts contribute to only 5% of the overall increase in accessible surface area. In this work several methodologies were used to characterize the thermal dissociation and unfolding of the TIM from Entamoeba histolytica (EhTIM) and a monomeric variant obtained by chemical derivatization (mEhTIM). During EhTIM unfolding, sequential transitions corresponding to dimer dissociation into a compact monomeric intermediate followed by unfolding and further aggregation of the intermediate occurred. In the case of mEhTIM, a single transition, analogous to the second transition of EhTIM, was observed. Calorimetric, spectroscopic, hydrodynamic, and functional evidence shows that dimer dissociation is not restricted to localized interface reorganization. Dissociation represents 55% (DeltaH(Diss) = 146.8 kcal mol(-1)) of the total enthalpy change (DeltaH(Tot) = 266 kcal mol(-1)), indicating that this process is linked to substantial unfolding. We propose that, rather than a rigid body process, subunit assembly is best represented by a fly-casting mechanism. In TIM, catalysis is restricted to the dimer; therefore, the interface can be viewed as the final nucleation motif that directs assembly, folding, and function.     

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.     

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.     

Structural analysis of new local features in SECIS RNA hairpins

Decoding of the UGA selenocysteine codon for selenoprotein translation requires the SECIS element, a stem-loop motif in the 3'-UTR of the mRNA carrying short or large apical loops. In previous structural studies, we derived a secondary structure model for SECIS RNAs with short apical loops. Work from others proposed that intra-apical loop base pairing can occur in those SECIS that possess large apical loops, yielding form 2 SECIS versus the form 1 with short loops. In this work, SECIS elements arising from eight different selenoprotein mRNAs were assayed by enzymatic and/or chemical probing showing that seven can adopt form 2. Further, database searches led to the discovery in drosophila and zebrafish of SECIS elements in the selenophosphate synthetase 2, type 1 deiodinase and SelW mRNAs. Alignment of SECIS sequences not only highlighted the predominance of form 2 but also made it possible to classify the SECIS elements according to the type of selenoprotein mRNA they belong to. Interestingly, the alignment revealed that an unpaired adenine, previously thought to be invariant, is replaced by a guanine in four SECIS elements. Tested in vivo, neither the A to G nor the A to U changes at this position greatly affected the activity while the most detrimental effect was provided by a C. The putative contribution of the various SECIS motifs to function and ligand binding is discussed.     

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.     

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.     

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.     

Determination of the elastic properties of short ssDNA molecules by mechanically folding and unfolding DNA hairpins

The characterization of elastic properties of biopolymers is crucial to understand many molecular reactions determined by conformational bending fluctuations of the polymer. Direct measurement of such elastic properties using single‐molecule methods is usually hindered by the intrinsic tendency of such biopolymers to form high‐order molecular structures. For example, single‐stranded deoxyribonucleic acids (ssDNA) tend to form secondary structures such as local double helices that prevent the direct measurement of the ideal elastic response of the ssDNA. In this work, we show how to extract the ideal elastic response in the entropic regime of short ssDNA molecules by mechanically pulling two‐state DNA hairpins of different contour lengths. This is achieved by measuring the force dependence of the molecular extension and stiffness on mechanically folding and unfolding the DNA hairpin. Both quantities are fit to the worm‐like chain elastic model giving values for the persistence length and the interphosphate distance. This method can be used to unravel the elastic properties of short ssDNA and RNA sequences and, more generally, any biopolymer that can exhibit a cooperative two‐state transition between mechanically folded and unfolded states (such as proteins). © 2014 Wiley Periodicals, Inc. Biopolymers 101: 1193–1199, 2014.