Thermodynamic characterization of RNA triloops

ABSTRACT: Relatively few thermodynamic parameters are available for RNA triloops. Therefore, 24 stemloop sequences containing naturally occurring triloops were optically melted, and the thermodynamic parameters ΔH°, ΔS°, ΔG°3 7 , and T M for each stem-loop were determined. These new experimental values, on average, are 0.5 kcal/mol different from the values predicted for these triloops using the model proposed by Mathews et al. [Mathews, D. H., Disney, M. D., Childs, J. L., Schroeder, S. J., Zuker, M., and Turner, D. H. (2004) Proc. Natl. Acad. Sci. U.S. A. 101, 7287-7292]. The data for the 24 triloops reported here were then combined with the data for five triloops that were published previously. A new model was derived to predict the free energy contribution of previously unmeasured triloops. The average absolute difference between the measured values and the values predicted using this proposed model is 0.3 kcal/mol. These new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA stemloops containing triloops and, furthermore, should allow for improved prediction of secondary structure from sequence. RNA stem-loops containing three nucleotides in the loop, triloops, are common secondary structure motifs found in naturally occurring RNA. For example, bacterial 16S rRNAs strongly favor tetraloops; however, the UUU triloop is the most common replacement (1). In the 16S-like rRNA variable regions, triloops account for 7% of the loops in bacteria and 16% of the loops in eukaryotes (2). Triloops are also found in large subunit rRNAs (3, 4), 5S rRNAs (5), signal recognition particles (6), RNase P RNAs (7), and group I introns (8, 9). More specifically, triloops are found in Brome mosaic virus (þ) strand RNA (10), human rhinovirus isotype 14 (11), iron responsive element RNA (12), and an RNA aptamer for bacteriophage MS2 coat protein (13), to name a few. Although relatively unstable due to the strain in the loop, triloops may be an important structural feature due to the accessibility of the loop nucleotides for recognition by proteins, other nucleic acids, or small molecules. It has been shown that triloops play a role in various biological processes, including virus replication The current model used by secondary structure prediction algorithms to predict the thermodynamic contribution of RNA triloops to stem-loop stability is sequence independent; all triloops contribute 5.4 kcal/mol to stem-loop stability, with the exception of 5 0 CCC3 0 which contributes 6.9 kcal/mol (21). In addition, there are two unstable triloop sequences (5 0 CAACG3 0 and 5 0 GUUAC3 0 ) for which this predictive model is not used; instead, the ΔG°3 7,loop values (6.8 and 6.9 kcal/mol, respectively) for these two triloops are provided in a lookup table (21). An interesting study by the Bevilacqua laboratory (19) used a combinatorial approach and temperature gradient gel electrophoresis to identify stable and unstable RNA triloops. It was discovered that sequence preferences for exceptionally stable triloops included a U-rich loop and C-G as the closing base pair. Although they used 10 mM NaCl during their melting experiments, they suggested that the rules for predicting triloop stability at 1 M NaCl should be modified; however, this has yet to be done. Here, we report the thermodynamic parameters for 24 previously unmeasured RNA triloops in 1 M NaCl and propose a new algorithm for predicting the contribution of triloops to stem-loop stability, which includes two bonuses for stabilizing sequence features. MATERIALS AND METHODS Compiling and Searching a Database for RNA Triloops. The initial aim of this project was to identify the most frequently occurring RNA triloops in nature and to thermodynamically characterize these hairpin triloop sequences. Therefore, a database of 1349 RNA secondary structures containing 123 small subunit rRNAs (22), 223 large subunit rRNAs (3, 4), 309 5S rRNAs (5), 484 tRNAs (23), 91 signal recognition particles (6), 16 RNase P RNAs (7), 100 group I introns (8, 9), and 3 group II introns (24) was compiled. This database was searched for triloops, and the number of occurrences for each type of triloop was tabulated. In this work, G-U pairs are considered to be canonical base pairs. Design of Sequences for Optical Melting Studies. Since most thermodynamic parameters for RNA secondary structure motifs are reported for RNA solutions containing 1 M NaCl, the melting buffer used in this work also contained 1 M NaCl. A major limitation of a thermodynamic analysis of RNA hairpins using this high salt concentration is the possible bimolecular

Revision of AMBER Torsional Parameters for RNA Improves Free Energy Predictions for Tetramer Duplexes with GC and iGiC Base Pairs

All-atom force fields are important for predicting thermodynamic, structural, and dynamic properties of RNA. In this paper, results are reported for thermodynamic integration calculations of free energy differences of duplex formation when CG pairs in the RNA duplexes r(CCGG)2, r(GGCC)2, r(GCGC)2, and r(CGCG)2 are replaced by isocytidine–isoguanosine (iCiG) pairs. Agreement with experiment was improved when ε/ζ, α/γ, β, and χ torsional parameters in the AMBER99 force field were revised on the basis of quantum mechanical calculations. The revised force field, AMBER99TOR, brings free energy difference predictions to within 1.3, 1.4, 2.3, and 2.6 kcal/mol at 300 K, respectively, compared to experimental results for the thermodynamic cycles of CCGG → iCiCiGiG, GGCC → iGiGiCiC, GCGC → iGiCiGiC, and CGCG → iCiGiCiG. In contrast, unmodified AMBER99 predictions for GGCC → iGiGiCiC and GCGC → iGiCiGiC differ from experiment by 11.7 and 12.6 kcal/mol, respectively. In order to test the dynamic stability of the above duplexes with AMBER99TOR, four individual 50 ns molecular dynamics (MD) simulations in explicit solvent were run. All except r(CCGG)2 retained A-form conformation for ≥82% of the time. This is consistent with NMR spectra of r(iGiGiCiC)2, which reveal an A-form conformation. In MD simulations, r(CCGG)2 retained A-form conformation 52% of the time, suggesting that its terminal base pairs may fray. The results indicate that revised backbone parameters improve predictions of RNA properties and that comparisons to measured sequence dependent thermodynamics provide useful benchmarks for testing force fields and computational methods

Synthesis of Janus compounds for the recognition of G-U mismatched nucleobase pairs

The design and synthesis of two Janus-type heterocycles with the capacity to simultaneously recognize guanine and uracyl in G-U mismatched pairs through complementary hydrogen bond pairing is described. Both compounds were conveniently functionalized with a carboxylic function and efficiently attached to a tripeptide sequence by using solid-phase methodologies. Ligands based on the derivatization of such Janus compounds with a small aminoglycoside, neamine, and its guanidinylated analogue have been synthesized, and their interaction with Tau RNA has been investigated by using several biophysical techniques, including UV-monitored melting curves, fluorescence titration experiments, and 1H NMR. The overall results indicated that Janus-neamine/guanidinoneamine showed some preference for the +3 mutated RNA sequence associated with the development of some tauopathies, although preliminary NMR studies have not confirmed binding to G-U pairs. Moreover, a good correlation has been found between the RNA binding affinity of such Janus-containing ligands and their ability to stabilize this secondary structure upon complexation

Improvement of RNA secondary structure prediction using RNase H cleavage and randomized oligonucleotides

RNA secondary structure prediction using free energy minimization is one method to gain an approximation of structure. Constraints generated by enzymatic mapping or chemical modification can improve the accuracy of secondary structure prediction. We report a facile method that identifies single-stranded regions in RNA using short, randomized DNA oligonucleotides and RNase H cleavage. These regions are then used as constraints in secondary structure prediction. This method was used to improve the secondary structure prediction of Escherichia coli 5S rRNA. The lowest free energy structure without constraints has only 27% of the base pairs present in the phylogenetic structure. The addition of constraints from RNase H cleavage improves the prediction to 100% of base pairs. The same method was used to generate secondary structure constraints for yeast tRNAPhe, which is accurately predicted in the absence of constraints (95%). Although RNase H mapping does not improve secondary structure prediction, it does eliminate all other suboptimal structures predicted within 10% of the lowest free energy structure. The method is advantageous over other single-stranded nucleases since RNase H is functional in physiological conditions. Moreover, it can be used for any RNA to identify accessible binding sites for oligonucleotides or small molecules

Elevated antibody to D-alanyl lipoteichoic acid indicates caries experience associated with fluoride and gingival health

BACKGROUND: Acidogenic, acid-tolerant bacteria induce dental caries and require D-alanyl glycerol lipoteichoic acid (D-alanyl LTA) on their cell surface. Because fluoride inhibits acid-mediated enamel demineralization, an elevated antibody response to D-alanyl LTA may indicate subjects with more acidogenic bacteria and, therefore, an association of DMFT with fluoride exposure and gingival health not apparent in low responders. METHODS: Cluster analysis was used to identify low antibody content. Within low and high responders (control and test subjects), the number of teeth that were decayed missing and filled (DMFT), or decayed only (DT) were regressed against fluoride exposure in the water supply and from dentrifice use. The latter was determined from gingival health: prevalences of plaque (PL) and bleeding on probing (BOP), and mean pocket depth (PD). Age was measured as a possible confounding cofactor. RESULTS: In 35 high responders, DMFT associated with length of exposure to fluoridated water (F score), PL and BOP (R(2) = 0.51, p < 0.001), whereas in 67 low D-ala-IgG responders, DMFT associated with PL, age, and PD (R(2) = 0.26, p < 0.001). BOP correlated strongly with number of 7 7 decayed teeth (DT) in 54 high responders (R(2) = 0.57, p < 0.001), but poorly in 97 low responders (R(2) = 0.12, p < 0.001). The strength of the PD association with DMFT, or of BOP with DT, in high responders significantly differed from that in low responders (p < 0.05). CONCLUSION: Caries associates with gingival health and fluoridated water exposure in high D-alanyl LTA antibody responders

The 3′ Splice Site of Influenza A Segment 7 mRNA Can Exist in Two Conformations: A Pseudoknot and a Hairpin

The 3′ splice site of influenza A segment 7 is used to produce mRNA for the M2 ion-channel protein, which is critical to the formation of viable influenza virions. Native gel analysis, enzymatic/chemical structure probing, and oligonucleotide binding studies of a 63 nt fragment, containing the 3′ splice site, key residues of an SF2/ASF splicing factor binding site, and a polypyrimidine tract, provide evidence for an equilibrium between pseudoknot and hairpin structures. This equilibrium is sensitive to multivalent cations, and can be forced towards the pseudoknot by addition of 5 mM cobalt hexammine. In the two conformations, the splice site and other functional elements exist in very different structural environments. In particular, the splice site is sequestered in the middle of a double helix in the pseudoknot conformation, while in the hairpin it resides in a two-by-two nucleotide internal loop. The results suggest that segment 7 mRNA splicing can be controlled by a conformational switch that exposes or hides the splice site

RNA structural analysis of the MYC mRNA reveals conserved motifs that affect gene expression.

The MYC gene encodes a human transcription factor and proto-oncogene that is dysregulated in over half of all known cancers. To better understand potential post-transcriptional regulatory features affecting MYC expression, we analyzed secondary structures in the MYC mRNA using a program that is optimized for finding small locally-folded motifs with a high propensity for function. This was accomplished by calculating folding metrics across the MYC sequence using a sliding analysis window and generating unique consensus base pairing models weighted by their lower-than-random predicted folding energy. A series of 30 motifs were identified, primarily in the 5' and 3' untranslated regions, which show evidence of structural conservation and compensating mutations across vertebrate MYC homologs. This analysis was able to recapitulate known elements found within an internal ribosomal entry site, as well as discover a novel element in the 3' UTR that is unusually stable and conserved. This novel motif was shown to affect MYC expression, potentially via the modulation of miRNA target accessibility or other trans-regulatory factors. In addition to providing basic insights into mechanisms that regulate MYC expression, this study provides numerous, potentially druggable RNA targets for the MYC gene, which is considered "undruggable" at the protein level