Detection of unusual RNA folding regions in HIV and SIV sequences

We have developed a method for detecting more stable and significantfolding regions relative to others in the sequence. The algorithmis based on the calculation of the lowest free energy of RNAsecondary structures and Monte Carlo simulation. For any givenRNA segment, the stability and statistical significance of RNAfolding are assessed by two measures: the stability score andthe significance score. The stability score measures the degreeof thermodynamic stability of the segment between all possiblebiological segments in the RNA sequence. The significance scorecharacterizes the specific arrangement of the nucleotides inthe segment that could imply a structural role for the sequenceinformation. Using these two measures, we are able to detecta series of distinct folding regions where highly stable andstatistically significant secondary structures occur in humanimmunodeficiency virus (HIV) and simian immunodeficiency virus(SIV) sequences. Received on April 4, 1990; accepted on October 2, 1990     

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.     

The Paradox of Dual Roles in the RNA World: Resolving the Conflict Between Stable Folding and Templating Ability

The hypothesized dual roles of RNA as both information carrier and biocatalyst during the earliest stages of life require a combination of features: good templating ability (for replication) and stable folding (for ribozymes). However, this poses the following paradox: well-folded sequences are poor templates for copying, but poorly folded sequences are unlikely to be good ribozymes. Here, we describe a strategy to overcome this dilemma through G:U wobble pairing in RNA. Unlike Watson–Crick base pairs, wobble pairs contribute highly to the energetic stability of the folded structure of their sequence, but only slightly, if at all, to the stability of the folded reverse complement. Sequences in the RNA World might thereby combine stable folding of the ribozyme with an unstructured, reverse-complementary genome, resulting in a “division of labor” between the strands. We demonstrate this strategy using computational simulations of RNA folding and an experimental model of early replication, nonenzymatic template-directed RNA primer extension. Additional study is needed to solve other problems associated with a complete replication cycle, including separation of strands after copying. Interestingly, viroid RNA sequences, which have been suggested to be relics of an RNA World (Diener, Proc Natl Acad Sci USA 86:9370–9374, 1989), also show significant asymmetry in folding energy between the infectious (+) and template (?) strands due to G:U pairing, suggesting that this strategy may even be used by replicators in the present day.     

Secondary structure alone is generally not statistically significant for the detection of noncoding RNAs

MOTIVATION: Several results in the literature suggest that biologically interesting RNAs have secondary structures that are more stable than expected by chance. Based on these observations, we developed a scanning algorithm for detecting noncoding RNA genes in genome sequences, using a fully probabilistic version of the Zuker minimum-energy folding algorithm. RESULTS: Preliminary results were encouraging, but certain anomalies led us to do a carefully controlled investigation of this class of methods. Ultimately, our results argue that for the probabilistic model there is indeed a statistical effect, but it comes mostly from local base-composition bias and not from RNA secondary structure. For the thermodynamic implementation (which evaluates statistical significance by doing Monte Carlo shuffling in fixed-length sequence windows, thus eliminating the base-composition effect) the signals for noncoding RNAs are still usually indistinguishable from noise, especially when certain statistical artifacts resulting from local base-composition inhomogeneity are taken into account. We conclude that although a distinct, stable secondary structure is undoubtedly important in most noncoding RNAs, the stability of most noncoding RNA secondary structures is not sufficiently different from the predicted stability of a random sequence to be useful as a general genefinding approach.     

Secondary structure and hybridization accessibility of hepatitis C virus 3'-terminal sequences

The 3'-terminal sequences of hepatitis C virus (HCV) positive- and negative-strand RNAs contribute cis-acting functions essential for viral replication. The secondary structure and protein-binding properties of these highly conserved regions are of interest not only for the further elucidation of HCV molecular biology, but also for the design of antisense therapeutic constructs. The RNA structure of the positive-strand 3' untranslated region has been shown previously to influence binding by various host and viral proteins and is thus thought to promote HCV RNA synthesis and genome stability. Recent studies have attributed analogous functions to the negative-strand 3' terminus. We evaluated the HCV negative-strand secondary structure by enzymatic probing with single-strand-specific RNases and thermodynamic modeling of RNA folding. The accessibility of both 3'-terminal sequences to hybridization by antisense constructs was evaluated by RNase H cleavage mapping in the presence of combinatorial oligodeoxynucleotide libraries. The mapping results facilitated identification of antisense oligodeoxynucleotides and a 10-23 deoxyribozyme active against the positive-strand 3'-X region RNA in vitro.     

Thermodynamic stability and statistical significance of potential stem-loop structures situated at the frameshift sites of retroviruses.

RNA stem-loop structures situated just 3' to the frameshift sites of the retroviral gag-pol or gag-pro and pro-pol regions may make important contributions to frame-shifting in retroviruses. In this study, the thermodynamic stability and statistical significance of such secondary structural features relative to others in the sequence have been assessed using a newly developed method that combines calculations of the lowest free energy of formation of RNA secondary structures and the Monte Carlo simulations. Our results show that stem-loop structures situated just 3' to the frameshift sites are both highly stable and statistically significant relative to others in the gag-pol or gag-pro and pro-pol junction domains (both 300 nucleotides upstream and downstream from the possible frameshift sites are included) of Rous sarcoma virus (RSV), human immunodeficiency virus (HIV-1), bovine leukemia virus (BLV), human T-cell leukemia virus type II (HTLV-II), and mouse mammary tumor virus (MMTV). No other more stable, or significant folding regions are predicted in these domains.     

Kinetics of double-chain reversals bridging contiguous quartets in tetramolecular quadruplexes

Repetitive 5'GGXGG DNA segments abound in, or near, regulatory regions of the genome and may form unusual structures called G-quadruplexes. Using NMR spectroscopy, we demonstrate that a family of 5'GCGGXGGY sequences adopts a folding topology containing double-chain reversals. The topology is composed of two bistranded quadruplex monomeric units linked by formation of G:C:G:C tetrads. We provide a complete thermodynamic and kinetic analysis of 13 different sequences using absorbance spectroscopy and DSC, and compare their kinetics with a canonical tetrameric parallel-stranded quadruplex formed by TG4T. We demonstrate large differences (up to 10(5)-fold) in the association constants of these quadruplexes depending on primary sequence; the fastest samples exhibiting association rate equal or higher than the canonical TG4T quadruplex. In contrast, all sequences studied here unfold at a lower temperature than this quadruplex. Some sequences have thermodynamic stability comparable to the canonical TG4T tetramolecular quadruplex, but with faster association and dissociation. Sequence effects on the dissociation processes are discussed in light of structural data.     

Beta-turn propensities as paradigms for the analysis of structural motifs to engineer protein stability.

The thermodynamic stability of a protein provides an experimental metric for the relationship of protein sequence and native structure. We have investigated an approach based on an analysis of the structural database for stability engineering of an immunoglobulin variable domain. The most frequently occurring residues in specific positions of beta-turn motifs were predicted to increase the folding stability of mutants that were constructed by site-directed mutagenesis. Even in positions in which different residues are conserved in immunoglobulin sequences, the predictions were confirmed. Frequently, mutants with increased beta-turn propensities display increased folding cooperativities, suggesting pronounced effects on the unfolded state independent of the expected effect on conformational entropy. We conclude that structural motifs with predominantly local interactions can serve as templates with which patterns of sequence preferences can be extracted from the database of protein structures. Such preferences can predict the stability effects of mutations for protein engineering and design.     

Thermodynamic and phylogenetic prediction of RNA secondary structures in the coding region of hepatitis C virus

The existence and functional importance of RNA secondary structure in the replication of positive-stranded RNA viruses is increasingly recognized. We applied several computational methods to detect RNA secondary structure in the coding region of hepatitis C virus (HCV), including thermodynamic prediction, calculation of free energy on folding, and a newly developed method to scan sequences for covariant sites and associated secondary structures using a parsimony-based algorithm. Each of the prediction methods provided evidence for complex RNA folding in the core- and NS5B-encoding regions of the genome. The positioning of covariant sites and associated predicted stem-loop structures coincided with thermodynamic predictions of RNA base pairing, and localized precisely in parts of the genome with marked suppression of variability at synonymous sites. Combined, there was evidence for a total of six evolutionarily conserved stem-loop structures in the NS5B-encoding region and two in the core gene. The virus most closely related to HCV, GB virus-B (GBV-B) also showed evidence for similar internal base pairing in its coding region, although predictions of secondary structures were limited by the absence of comparative sequence data for this virus. While the role(s) of stem-loops in the coding region of HCV and GBV-B are currently unknown, the structure predictions in this study could provide the starting point for functional investigations using recently developed self-replicating clones of HCV.     

Molecular characterization and interstrain variability of pHPS1, a plasmid isolated from the Sydney strain (SS1) of Helicobacter pylori

The 5846-bp circular plasmid pHPS1 of Helicobacter pylori Sydney strain, SS1, was cloned, sequenced, and structurally characterized. The SS1 strain is widely used in animal studies of H. pylori infection. The sequence of pHPS1 revealed three open reading frames (ORFs), all of which are transcribed. Two ORFs encode putative plasmid replication proteins, RepA and RepB, similar to replicases resident on theta plasmids. In contrast, the function of ORF2 remains cryptic due to the absence of sequence similarity with any known protein in sequence databases. In addition, species specificity of these three coding regions was shown using DNA dot blot hybridization in 57 diverse clinical H. pylori isolates and 32 Helicobacter and Campylobacter strains. RepA appears to be the predominant plasmid replication protein of H. pylori and the deduced amino acid sequence was highly conserved (76-96%) in 8 H. pylori isolates, including SS1. RepB was detected in 3 H. pylori isolates examined in this study, 2 of which possess only the repB gene. Analysis of the protein sequences of these two replicases, together with previously characterized H. pylori plasmid replication proteins, supports the formation of a distinct class of H. pylori plasmid proteins. Moreover, comprehensive analysis of the whole genome sequence of H. pylori strain 26695, pHPS1, and other H. pylori plasmid sequences that are available revealed interesting insights as to the occurrence of plasmid-mediated recombination within H. pylori. Common regions between plasmids and chromosome sequences of H. pylori were identified in this study which could only have arisen by genetic recombination, thus providing the first line of evidence, albeit indirectly, of the contribution of H. pylori plasmids in generating an extensive genetic heterogeneity characteristic of this important gastroduodenal pathogen.     

Elements of thermodynamics in RNA evolution.

The paper presents some aspects correlating thermal stability of RNA folding and the occurrence of structural motifs in natural ribonucleic acids. Particularly, the thermodynamic stability of 2'-5' and 3'-5' linked RNA and the contribution of unpaired terminal nucleotides (dangling ends) in secondary (2D) and tertiary (3D) structures of RNA are discussed. Both examples suggest that during evolution nature selected sequences and structures of RNA which are the most thermally stable and efficient for their biological function.     

A highly conserved RNA folding region coincident with the Rev response element of primate immunodeficiency viruses.

A series of unusual folding regions (UFR) immediately 3' to the cleavage site of the outer membrane protein (OMP) and transmembrane protein (TMP) were detected in the envelope gene RNA of the human immunodeficiency virus (HIV-1, HIV-2) and simian immunodeficiency virus (SIV) by an extensive Monte Carlo simulation. These RNA secondary structures were predicted to be both highly stable and statistically significant. In the calculation, twenty-five different sequence isolates of HIV-1, three isolates of HIV-2 and eight sequences of SIV were included. Although significant sequence divergence occurs in the env coding regions of these viruses, a distinct UFR of 234-nt is consistently located ten nucleotides 3' to the cleavage site of the OMP/TMP in HIV-1, and a 216-nt UFR occurs forty-six and forty-nine nucleotides downstream from the OMP/TMP cleavage site of HIV-2 and SIV, respectively. Compensatory base changes in the helical stem regions of these conserved RNA secondary structures are identified. These results support the hypothesis that these special RNA folding regions are functionally important and suggest that the role of this sequence as the Rev response element (RRE) is mediated by secondary structure as well as primary RNA sequence.     

Correlation between sequence conservation and structural thermodynamics of microRNA precursors from human,mouse, and chicken genomes

Background

Previous studies have shown that microRNA precursors (pre-miRNAs) have considerably more stable secondary structures than other native RNAs (tRNA, rRNA, and mRNA) and artificial RNA sequences. However, pre-miRNAs with ultra stable secondary structures have not been investigated. It is not known if there is a tendency in pre-miRNA sequences towards or against ultra stable structures? Furthermore, the relationship between the structural thermodynamic stability of pre-miRNA and their evolution remains unclear.

Results

We investigated the correlation between pre-miRNA sequence conservation and structural stability as measured by adjusted minimum folding free energies in pre-miRNAs isolated from human, mouse, and chicken. The analysis revealed that conserved and non-conserved pre-miRNA sequences had structures with similar average stabilities. However, the relatively ultra stable and unstable pre-miRNAs were more likely to be non-conserved than pre-miRNAs with moderate stability. Non-conserved pre-miRNAs had more G+C than A+U nucleotides, while conserved pre-miRNAs contained more A+U nucleotides. Notably, the U content of conserved pre-miRNAs was especially higher than that of non-conserved pre-miRNAs. Further investigations showed that conserved and non-conserved pre-miRNAs exhibited different structural element features, even though they had comparable levels of stability.

Conclusions

We proposed that there is a correlation between structural thermodynamic stability and sequence conservation for pre-miRNAs from human, mouse, and chicken genomes. Our analyses suggested that pre-miRNAs with relatively ultra stable or unstable structures were less favoured by natural selection than those with moderately stable structures. Comparison of nucleotide compositions between non-conserved and conserved pre-miRNAs indicated the importance of U nucleotides in the pre-miRNA evolutionary process. Several characteristic structural elements were also detected in conserved pre-miRNAs.

     

Consensus folding of unaligned RNA sequences revisited.

As one of the earliest problems in computational biology, RNA secondary structure prediction (sometimes referred to as "RNA folding") problem has attracted attention again, thanks to the recent discoveries of many novel non-coding RNA molecules. The two common approaches to this problem are de novo prediction of RNA secondary structure based on energy minimization and the consensus folding approach (computing the common secondary structure for a set of unaligned RNA sequences). Consensus folding algorithms work well when the correct seed alignment is part of the input to the problem. However, seed alignment itself is a challenging problem for diverged RNA families. In this paper, we propose a novel framework to predict the common secondary structure for unaligned RNA sequences. By matching putative stacks in RNA sequences, we make use of both primary sequence information and thermodynamic stability for prediction at the same time. We show that our method can predict the correct common RNA secondary structures even when we are given only a limited number of unaligned RNA sequences, and it outperforms current algorithms in sensitivity and accuracy.     

Funnel-like organization in sequence space determines the distributions of protein stability and folding rate preferred by evolution

To understand the physical and evolutionary determinants of protein folding, we map out the complete organization of thermodynamic and kinetic properties for protein sequences that share the same fold. The exhaustive nature of our study necessitates using simplified models of protein folding. We obtain a stability map and a folding rate map in sequence space. Comparison of the two maps reveals a common organizational principle: optimality decreases more or less uniformly with distance from the optimal sequence in the sequence space. This gives a funnel-shaped optimality surface. Evolutionary dynamics of a sequence population on these two maps reveal how the simple organization of sequence space affects the distributions of stability and folding rate preferred by evolution.     

Entropy-driven folding of an RNA helical junction: an isothermal titration calorimetric analysis of the hammerhead ribozyme

Helical junctions are extremely common motifs in naturally occurring RNAs, but little is known about the thermodynamics that drive their folding. Studies of junction folding face several challenges: non-two-state folding behavior, superposition of secondary and tertiary structural energetics, and drastically opposing enthalpic and entropic contributions to folding. Here we describe a thermodynamic dissection of the folding of the hammerhead ribozyme, a three-way RNA helical junction, by using isothermal titration calorimetry of bimolecular RNA constructs. By using this method, we show that tertiary folding of the hammerhead core occurs with a highly unfavorable enthalpy change, and is therefore entropically driven. Furthermore, the enthalpies and heat capacities of core folding are the same whether supported by monovalent or divalent ions. These properties appear to be general to the core sequence of bimolecular hammerhead constructs. We present a model for the ion-induced folding of the hammerhead core that is similar to those advanced for the folding of much larger RNAs, involving ion-induced collapse to a structured, non-native state accompanied by rearrangement of core residues to produce the native fold. In agreement with previous enzymological and structural studies, our thermodynamic data suggest that the hammerhead structure is stabilized in vitro predominantly by diffusely bound ions. Our approach addresses several significant challenges that accompany the study of junction folding, and should prove useful in defining the thermodynamic determinants of stability in these important RNA motifs.     

Proline residues within spacer peptide p1 are important for human immunodeficiency virus type 1 infectivity,protein processing,and genomic RNA dimer stability

The full-length human immunodeficiency virus type 1 (HIV-1) mRNA encodes two precursor polyproteins, Gag and GagProPol. An infrequent ribosomal frameshifting event allows these proteins to be synthesized from the same mRNA in a predetermined ratio of 20 Gag proteins for each GagProPol. The RNA frameshift signal consists of a slippery sequence and a hairpin stem-loop whose thermodynamic stability has been shown in in vitro translation systems to be critical to frameshifting efficiency. In this study we examined the frameshift region of HIV-1, investigating the effects of altering stem-loop stability in the context of the complete viral genome and assessing the role of the Gag spacer peptide p1 and the GagProPol transframe (TF) protein that are encoded in this region. By creating a series of frameshift region mutants that systematically altered the stability of the frameshift stem-loop and the protein sequences of the p1 spacer peptide and TF protein, we have demonstrated the importance of stem-loop thermodynamic stability in frameshifting efficiency and viral infectivity. Multiple changes to the amino acid sequence of p1 resulted in altered protein processing, reduced genomic RNA dimer stability, and abolished viral infectivity. The role of the two highly conserved proline residues in p1 (position 7 and 13) was also investigated. Replacement of the two proline residues by leucines resulted in mutants with altered protein processing and reduced genomic RNA dimer stability that were also noninfectious. The unique ability of proline to confer conformational constraints on a peptide suggests that the correct folding of p1 may be important for viral function.     

Selective Pressures on RNA Hairpins In Vivo and In Vitro

Comparison of the most stable potential hairpins in the sequences of natural ribozymes with those in the randomized sequences has revealed that the hairpin loop energies are lower than expected by chance. Although these hairpins are not necessarily parts of functional structures, there is a selective pressure to diminish the destabilizing free energies of the hairpin loops. In contrast, no significant bias is observed in the stacking values of the most stable stems. In the ribozymes isolated in vitro the loops of potential hairpins are closer to random values, which can result in less efficient folding rates. Furthermore, the effects of kinetic traps seem to be more significant in the folding pathways of the in vitro isolates due to a potential to form stable stacks incompatible with the functional folds. Similarly to natural ribozyme sequences, the untranslated regions of viral RNAs also form hairpins with relatively low loop free energies. These evolutionary trends suggest ways for efficient engineering of improved RNA constructs on the basis of analysis of in vitro isolates and approaches for the search of regions coding for functional RNA structures in large genome sequences. Received: 12 January 2001 / Accepted: 21 May 2001