Mechanistic insights into temperature-dependent regulation of the simple cyanobacterial hsp17 RNA thermometer at base-pair resolution

The cyanobacterial hsp17 ribonucleicacid thermometer (RNAT) is one of the smallest naturally occurring RNAT. It forms a single hairpin with an internal 1×3-bulge separating the start codon in stem I from the ribosome binding site (RBS) in stem II. We investigated the temperature-dependent regulation of hsp17 by mapping individual base-pair stabilities from solvent exchange nuclear magnetic resonance (NMR) spectroscopy. The wild-type RNAT was found to be stabilized by two critical CG base pairs (C14-G27 and C13-G28). Replacing the internal 1×3 bulge by a stable CG base pair in hsp17^rep significantly increased the global stability and unfolding cooperativity as evidenced by circular dichroism spectroscopy. From the NMR analysis, remote stabilization and non-nearest neighbour effects exist at the base-pair level, in particular for nucleotide G28 (five nucleotides apart from the side of mutation). Individual base-pair stabilities are coupled to the stability of the entire thermometer within both the natural and the stabilized RNATs by enthalpy–entropy compensation presumably mediated by the hydration shell. At the melting point the Gibbs energies of the individual nucleobases are equalized suggesting a consecutive zipper-type unfolding mechanism of the RBS leading to a dimmer-like function of hsp17 and switch-like regulation behaviour of hsp17^rep. The data show how minor changes in the nucleotide sequence not only offset the melting temperature but also alter the mode of temperature sensing. The cyanobacterial thermosensor demonstrates the remarkable adjustment of natural RNATs to execute precise temperature control.     

Direct observation of the temperature-induced melting process of the Salmonella fourU RNA thermometer at base-pair resolution

In prokaryotes, RNA thermometers regulate a number of heat shock and virulence genes. These temperature sensitive RNA elements are usually located in the 5′-untranslated regions of the regulated genes. They repress translation initiation by base pairing to the Shine–Dalgarno sequence at low temperatures. We investigated the thermodynamic stability of the temperature labile hairpin 2 of the Salmonella fourU RNA thermometer over a broad temperature range and determined free energy, enthalpy and entropy values for the base-pair opening of individual nucleobases by measuring the temperature dependence of the imino proton exchange rates via NMR spectroscopy. Exchange rates were analyzed for the wild-type (wt) RNA and the A8C mutant. The wt RNA was found to be stabilized by the extraordinarily stable G14–C25 base pair. The mismatch base pair in the wt RNA thermometer (A8–G31) is responsible for the smaller cooperativity of the unfolding transition in the wt RNA. Enthalpy and entropy values for the base-pair opening events exhibit linear correlation for both RNAs. The slopes of these correlations coincide with the melting points of the RNAs determined by CD spectroscopy. RNA unfolding occurs at a temperature where all nucleobases have equal thermodynamic stabilities. Our results are in agreement with a consecutive zipper-type unfolding mechanism in which the stacking interaction is responsible for the observed cooperativity. Furthermore, remote effects of the A8C mutation affecting the stability of nucleobase G14 could be identified. According to our analysis we deduce that this effect is most probably transduced via the hydration shell of the RNA.     

Context-sensitivity of isosteric substitutions of non-Watson–Crick basepairs in recurrent RNA 3D motifs

Sequence variation in a widespread, recurrent, structured RNA 3D motif, the Sarcin/Ricin (S/R), was studied to address three related questions: First, how do the stabilities of structured RNA 3D motifs, composed of non-Watson–Crick (non-WC) basepairs, compare to WC-paired helices of similar length and sequence? Second, what are the effects on the stabilities of such motifs of isosteric and non-isosteric base substitutions in the non-WC pairs? And third, is there selection for particular base combinations in non-WC basepairs, depending on the temperature regime to which an organism adapts? A survey of large and small subunit rRNAs from organisms adapted to different temperatures revealed the presence of systematic sequence variations at many non-WC paired sites of S/R motifs. UV melting analysis and enzymatic digestion assays of oligonucleotides containing the motif suggest that more stable motifs tend to be more rigid. We further found that the base substitutions at non-Watson–Crick pairing sites can significantly affect the thermodynamic stabilities of S/R motifs and these effects are highly context specific indicating the importance of base-stacking and base-phosphate interactions on motif stability. This study highlights the significance of non-canonical base pairs and their contributions to modulating the stability and flexibility of RNA molecules.     

Pronounced instability of tandem IU base pairs in RNA

Optical melting was used to determine the stabilities of three series of RNA oligomers containing tandem XU base pairs, GGCXUGCC (5′XU3′), GGCUXGCC (5′UX3′) and GGCXXGGC/CCGUUCCG (5′XX3′), where X is either A, G or I (inosine). The helices containing tandem AU base pairs were the most stable in the first two series (5′XU3′ and 5′UX3′), with an average melting temperature ~11°C higher than the helices with tandem 5′GU3′ base pairs and 25°C higher than the helices with tandem 5′IU3′ base pairs. For the third series (5′XX3′), the helix containing tandem GG is the most stable, with an average melting temperature ~2°C higher than the helix with tandem AA base pairs and ~24°C higher than the helix with tandem II base pairs. The thermodynamic stability of the oligomers with tandem IU base pairs was also investigated as a function of magnesium ion concentration. As with normal A–U or G–U tandem duplexes, the data could best be interpreted as non-specific binding of magnesium ions to the inosine-containing RNA oligonucleotides.     

Noncanonical G(syn)-G(anti) base pairs stabilized by sulphate anions in two X-ray structures of the (GUGGUCUGAUGAGGCC) RNA duplex

The structures of two crystal forms of the RNA 16-mer with the sequence GUGGUCUGAUGAGGCC, grown in the presence of a high concentration of sulphate ions, have been determined using synchrotron radiation at 1.4- and 2.0-Å resolution. RNA with this sequence is known as one of the two strands of the noncleavable form of the hammerhead ribozyme. In both crystal structures, two G(syn)–G(anti) noncanonical base pairs are observed in the middle of a 14 base-pair (bp) duplex having 5′-dangling GU residues. Both structures contain sulphate anions interacting with the G–G bp stabilizing G in its syn conformation and bridging the two RNA strands. In both cases the interactions take place in the major groove, although the anions are accommodated within different helix geometries, most pronounced in the changing width of the major groove. In one structure, where a single sulphate spans both G–G pairs, the major groove is closed around the anion, while in the other structure, where each of the two G–G pairs is associated with a separate sulphate, the groove is open. This work provides the first examples of a G–G pair in syn-anti conformation, which minimizes the purine–purine clash in the center of the duplex, while utilizing its residual hydrogen bonding potential in specific interactions with sulphate anions.     

Predicting loop–helix tertiary structural contacts in RNA pseudoknots

Tertiary interactions between loops and helical stems play critical roles in the biological function of many RNA pseudoknots. However, quantitative predictions for RNA tertiary interactions remain elusive. Here we report a statistical mechanical model for the prediction of noncanonical loop–stem base-pairing interactions in RNA pseudoknots. Central to the model is the evaluation of the conformational entropy for the pseudoknotted folds with defined loop–stem tertiary structural contacts. We develop an RNA virtual bond-based conformational model (Vfold model), which permits a rigorous computation of the conformational entropy for a given fold that contains loop–stem tertiary contacts. With the entropy parameters predicted from the Vfold model and the energy parameters for the tertiary contacts as inserted parameters, we can then predict the RNA folding thermodynamics, from which we can extract the tertiary contact thermodynamic parameters from theory–experimental comparisons. These comparisons reveal a contact enthalpy (ΔH) of −14 kcal/mol and a contact entropy (ΔS) of −38 cal/mol/K for a protonated C⁺•(G–C) base triple at pH 7.0, and (ΔH = −7 kcal/mol, ΔS = −19 cal/mol/K) for an unprotonated base triple. Tests of the model for a series of pseudoknots show good theory–experiment agreement. Based on the extracted energy parameters for the tertiary structural contacts, the model enables predictions for the structure, stability, and folding pathways for RNA pseudoknots with known or postulated loop–stem tertiary contacts from the nucleotide sequence alone.     

Nearest neighbor parameters for Watson-Crick complementary heteroduplexes formed between 2'-O-methyl RNA and RNA oligonucleotides

Results from optical melting studies of Watson–Crick complementary heteroduplexes formed between 2′-O-methyl RNA and RNA oligonucleotides are used to determine nearest neighbor thermodynamic parameters for predicting the stabilities of such duplexes. The results are consistent with the physical model assumed by the individual nearest neighbor-hydrogen bonding model, which contains terms for helix initiation, base pair stacking and base pair composition. The sequence dependence is similar to that for Watson–Crick complementary RNA/RNA duplexes, which suggests that the sequence dependence may also be similar to that for other backbones that favor A-form RNA conformations.     

Vfold: A Web Server for RNA Structure and Folding Thermodynamics Prediction

Background

The ever increasing discovery of non-coding RNAs leads to unprecedented demand for the accurate modeling of RNA folding, including the predictions of two-dimensional (base pair) and three-dimensional all-atom structures and folding stabilities. Accurate modeling of RNA structure and stability has far-reaching impact on our understanding of RNA functions in human health and our ability to design RNA-based therapeutic strategies.

Results

The Vfold server offers a web interface to predict (a) RNA two-dimensional structure from the nucleotide sequence, (b) three-dimensional structure from the two-dimensional structure and the sequence, and (c) folding thermodynamics (heat capacity melting curve) from the sequence. To predict the two-dimensional structure (base pairs), the server generates an ensemble of structures, including loop structures with the different intra-loop mismatches, and evaluates the free energies using the experimental parameters for the base stacks and the loop entropy parameters given by a coarse-grained RNA folding model (the Vfold model) for the loops. To predict the three-dimensional structure, the server assembles the motif scaffolds using structure templates extracted from the known PDB structures and refines the structure using all-atom energy minimization.

Conclusions

The Vfold-based web server provides a user friendly tool for the prediction of RNA structure and stability. The web server and the source codes are freely accessible for public use at “http://rna.physics.missouri.edu”.     

Characterization of a novel eukaryal nick-sealing RNA ligase from Naegleria gruberi

The proteome of the amoebo-flagellate protozoan Naegleria gruberi is rich in candidate RNA repair enzymes, including 15 putative RNA ligases, one of which, NgrRnl, is a eukaryal homolog of Deinococcus radiodurans RNA ligase, DraRnl. Here we report that purified recombinant NgrRnl seals nicked 3′-OH/5′-PO₄ duplexes in which the 3′-OH strand is RNA. It does so via the “classic” ligase pathway, entailing reaction with ATP to form a covalent NgrRnl–AMP intermediate, transfer of AMP to the nick 5′-PO₄, and attack of the RNA 3′-OH on the adenylylated nick to form a 3′–5′ phosphodiester. Unlike members of the four known families of ATP-dependent RNA ligases, NgrRnl lacks a carboxy-terminal appendage to its nucleotidyltransferase domain. Instead, it contains a defining amino-terminal domain that we show is important for 3′-OH/5′-PO₄ nick-sealing and ligase adenylylation, but dispensable for phosphodiester synthesis at a preadenylylated nick. We propose that NgrRnl, DraRnl, and their homologs from diverse bacteria, viruses, and unicellular eukarya comprise a new “Rnl5 family” of nick-sealing ligases with a signature domain organization.     

Phenotype Bias Determines How Natural RNA Structures Occupy the Morphospace of All Possible Shapes

Morphospaces—representations of phenotypic characteristics—are often populated unevenly, leaving large parts unoccupied. Such patterns are typically ascribed to contingency, or else to natural selection disfavoring certain parts of the morphospace. The extent to which developmental bias, the tendency of certain phenotypes to preferentially appear as potential variation, also explains these patterns is hotly debated. Here we demonstrate quantitatively that developmental bias is the primary explanation for the occupation of the morphospace of RNA secondary structure (SS) shapes. Upon random mutations, some RNA SS shapes (the frequent ones) are much more likely to appear than others. By using the RNAshapes method to define coarse-grained SS classes, we can directly compare the frequencies that noncoding RNA SS shapes appear in the RNAcentral database to frequencies obtained upon a random sampling of sequences. We show that: 1) only the most frequent structures appear in nature; the vast majority of possible structures in the morphospace have not yet been explored; 2) remarkably small numbers of random sequences are needed to produce all the RNA SS shapes found in nature so far; and 3) perhaps most surprisingly, the natural frequencies are accurately predicted, over several orders of magnitude in variation, by the likelihood that structures appear upon a uniform random sampling of sequences. The ultimate cause of these patterns is not natural selection, but rather a strong phenotype bias in the RNA genotype–phenotype map, a type of developmental bias or “findability constraint,” which limits evolutionary dynamics to a hugely reduced subset of structures that are easy to “find.”     

Thermodynamic and hydration effects for the incorporation of a cationic 3-aminopropyl chain into DNA

The introduction of cationic 5-(ω-aminoalkyl)-2′-deoxypyrimidines into duplex DNA has been shown to induce DNA bending. In order to understand the energetic and hydration contributions for the incorporation of a cationic side chain in DNA a combination of spectroscopy, calorimetry and density techniques were used. Specifically, the temperature unfolding and isothermal formation was studied for a pair of duplexes with sequence d(CGTAGUCG TGC)/d(GCACGACTACG), where U represents 2′-deoxyuridine (‘control’) or 5-(3-aminopropyl)-2′-deoxyuridine (‘modified’). Continuous variation experiments confirmed 1:1 stoichiometries for each duplex and the circular dichroism spectra show that both duplexes adopted the B conformation. UV and differential scanning calorimetry melting experiments reveal that each duplex unfolds in two-state transitions. In low salt buffer, the ‘modified’ duplex is more stable and unfolds with a lower endothermic heat and lower release of counterion and water. This electrostatic stabilization is entropy driven and disappears at higher salt concentrations. Complete thermodynamic profiles at 15°C show that the favorable formation of each duplex results from the compensation of a favorable exothermic heat with an unfavorable entropy contribution. However, the isothermal profiles yielded a differential enthalpy of 8.8 kcal/mol, which is 4.3 kcal/mol higher than the differential enthalpy observed in the unfolding profiles. This indicates that the presence of the aminopropyl chain induces an increase in base stacking interactions in the modified single strand and a decrease in base stacking interactions in the modified duplex. Furthermore, the formation of the ‘control’ duplex releases water while the ‘modified’ duplex takes up water. Relative to the control duplex, formation of the modified duplex at 15°C yielded a marginal differential ΔG° term, positive ΔΔH_ITC–Δ(TΔS) compensation, negative ΔΔV and a net release of counterions. The opposite signs of the differential enthalpy–entropy compensation and differential volume change terms show a net uptake of structural water around polar and non-polar groups. This indicates that incorporation of the aminopropyl chain induces a higher exposure of aromatic bases to the solvent, which may be consistent with a small and local bend in the ‘modified’ duplex.     

Prevalence and Correlates of HIV Testing among Young People Enrolled in Non-Formal Education Centers in Urban Chiang Mai,Thailand: A Cross-Sectional Study

BackgroundHIV testing is the gateway to HIV prevention, treatment, and care. Despite the established vulnerability of young Thai people to HIV infection, studies examining the prevalence and correlates of HIV testing among the general population of Thai youth are still very limited. This study investigates socio-demographic, behavioral, and psychosocial factors associated with HIV testing among young Thai people enrolled in Non-formal Education Centers (NFEC) in urban Chiang Mai, Northern Thailand.MethodsThis was a cross-sectional quantitative study conducted among young unmarried Thai youth—between the ages of 15 and 24—who were enrolled in NFEC in urban Chiang Mai. Multiple logistic regressions were used to identify correlates of “ever tested for HIV” among the sexually active participants.FindingsOf the 295 sexually active participants, 27.3% reported “ever tested for HIV;” 65.4% “did not consistently use condom;” and 61.7% “had at least 2 lifetime partners.” We found that “self-efficacy” (AOR, 4.92; CI, 1.22–19.73); “perception that it is easy to find a location nearby to test for HIV” (AOR, 4.67; CI, 1.21–18.06); “having at least 2 lifetime sexual partners” (AOR, 2.05; CI, 1.09–3.85); and “ever been pregnant or made someone pregnant” (AOR, 4.06; CI, 2.69–9.15); were associated with increased odds of having ever been tested. On the other hand, “fear of HIV test results” (AOR, 0.21; CI, 0.08–0.57) was associated with lower odds of ever having been tested for HIV.ConclusionThe main finding is that a substantially high proportion of Thai youth is engaged in risky sexual behaviors—yet reports low rates of ever having been tested for HIV. This highlights an urgent need to develop appropriate interventions—based on the identified correlates of HIV testing. There is also an urgent need to enhance HIV testing and to promote safer sexual behaviors among young Thai people—particularly those who are out-of-school.     

Hydration of short DNA, RNA and 2'-OMe oligonucleotides determined by osmotic stressing

Studies on hydration are important for better understanding of structure and function of nucleic acids. We compared the hydration of self-complementary DNA, RNA and 2′-O-methyl (2′-OMe) oligonucleotides GCGAAUUCGC, (UA)₆ and (CG)₃ using the osmotic stressing method. The number of water molecules released upon melting of oligonucleotide duplexes, Δn_W, was calculated from the dependence of melting temperature on water activity and the enthalpy, both measured with UV thermal melting experiments. The water activity was changed by addition of ethylene glycol, glycerol and acetamide as small organic co-solutes. The Δn_W was 3–4 for RNA duplexes and 2–3 for DNA and 2′-OMe duplexes. Thus, the RNA duplexes were hydrated more than the DNA and the 2′-OMe oligonucleotide duplexes by approximately one to two water molecules depending on the sequence. Consistent with previous studies, GC base pairs were hydrated more than AU pairs in RNA, whereas in DNA and 2′-OMe oligonucleotides the difference in hydration between these two base pairs was relatively small. Our data suggest that the better hydration of RNA contributes to the increased enthalpic stability of RNA duplexes compared with DNA duplexes.     

tRNA anticodon shifts in eukaryotic genomes

Embedded in the sequence of each transfer RNA are elements that promote specific interactions with its cognate aminoacyl tRNA-synthetase. Although many such “identity elements” are known, their detection is difficult since they rely on unique structural signatures and the combinatorial action of multiple elements spread throughout the tRNA molecule. Since the anticodon is often a major identity determinant itself, it is possible to switch between certain tRNA functional types by means of anticodon substitutions. This has been shown to have occurred during the evolution of some genomes; however, the scale and relevance of “anticodon shifts” to the evolution of the tRNA multigene family is unclear. Using a synteny-conservation–based method, we detected tRNA anticodon shifts in groups of closely related species: five primates, 12 Drosophila, six nematodes, 11 Saccharomycetes, and 61 Enterobacteriaceae. We found a total of 75 anticodon shifts: 31 involving switches of identity (alloacceptor shifts) and 44 between isoacceptors that code for the same amino acid (isoacceptor shifts). The relative numbers of shifts in each taxa suggest that tRNA gene redundancy is likely the driving factor, with greater constraint on changes of identity. Sites that frequently covary with alloacceptor shifts are located at the extreme ends of the molecule, in common with most known identity determinants. Isoacceptor shifts are associated with changes in the midsections of the tRNA sequence. However, the mutation patterns of anticodon shifts involving the same identities are often dissimilar, suggesting that alternate sets of mutation may achieve the same functional compensation.