RNA tertiary interactions in a riboswitch stabilize the structure of a kink turn

The kink turn is a widespread RNA motif that introduces an acute kink into the axis of duplex RNA, typically comprising a bulge followed by a G?A and A?G pairs. The kinked conformation is stabilized by metal ions, or the binding of proteins including L7Ae. We now demonstrate a third mechanism for the stabilization of k-turn structure, involving tertiary interactions within a larger RNA structure. The SAM-I riboswitch contains an essential standard k-turn sequence that kinks a helix so that its terminal loop can make?a?long-range interaction. We find that some sequence variations in the k-turn within the riboswitch do not prevent SAM binding, despite preventing the folding of the k-turn in isolation. Furthermore, two crystal structures show that the sequence-variant k-turns are conventionally folded within the riboswitch. This study shows that the folded structure of the k-turn can be stabilized by tertiary interactions within a larger RNA structure.     

Second eigenvalue of the Laplacian matrix for predicting RNA conformational switch by mutation

MOTIVATION: Conformational switching in RNAs is thought to be of fundamental importance in several biological processes, including translational regulation, regulation of self-cleavage in viruses, protein biosynthesis and mRNA splicing. Current methods for detecting bi-stable RNAs that can lead to structural switching when triggered by an outside event rely on kinetics, energetics and properties of the combinatorial structure space of RNAs. Based on these properties, tools have been developed to predict whether a given sequence folds to a structure characterized by a bi-stable conformation, or to design multi-stable RNAs by an iterative algorithm. A useful addition is in developing a local procedure to prescribe, given an initial sequence, the least amount of mutations needed to drive the system into an optimal bi-stable conformation. RESULTS: We introduce a local procedure for predicting mutations, by generating and analyzing eigenvalue tables, that are capable of transforming the wild-type sequence into a bi-stable conformation. The method is independent of the folding algorithms but relies on their success. It can be used in conjunction with existing tools, as well as being incorporated into more general RNA prediction packages. We apply this procedure on three well-studied structures. First, the method is validated on the mutation leading to a conformational switch in the spliced leader RNA from Leptomonas collosoma, a mutation that has already been confirmed by an experiment. Second, the method is used to predict a mutation that can lead to a novel conformational switch in the P5abc subdomain of the group I intron ribozyme in Tetrahymena thermophila. Third, the method is applied on Hepatitis delta virus to predict mutations that transform the wild-type into a bi-stable conformation, a configuration assessed by calculating the free energies using folding prediction algorithms. The predictions in the final examples need to be verified experimentally, whereas the mutation predicted in the first example complies with the experiment. This supports the use of our proposed method on other known structures, as well as genetically engineered ones. AVAILABILITY: An eigenvalue application will be available in the near future attached to one of the existing tools.     

RNA thermoswitches modulate Staphylococcus aureus adaptation to ambient temperatures

Thermoregulation of virulence genes in bacterial pathogens is essential for environment-to-host transition. However, the mechanisms governing cold adaptation when outside the host remain poorly understood. Here, we found that the production of cold shock proteins CspB and CspC from Staphylococcus aureus is controlled by two paralogous RNA thermoswitches. Through in silico prediction, enzymatic probing and site-directed mutagenesis, we demonstrated that cspB and cspC 5′UTRs adopt alternative RNA structures that shift from one another upon temperature shifts. The open (O) conformation that facilitates mRNA translation is favoured at ambient temperatures (22°C). Conversely, the alternative locked (L) conformation, where the ribosome binding site (RBS) is sequestered in a double-stranded RNA structure, is folded at host-related temperatures (37°C). These structural rearrangements depend on a long RNA hairpin found in the O conformation that sequesters the anti-RBS sequence. Notably, the remaining S. aureus CSP, CspA, may interact with a UUUGUUU motif located in the loop of this long hairpin and favour the folding of the L conformation. This folding represses CspB and CspC production at 37°C. Simultaneous deletion of the cspB/cspC genes or their RNA thermoswitches significantly decreases S. aureus growth rate at ambient temperatures, highlighting the importance of CspB/CspC thermoregulation when S. aureus transitions from the host to the environment.     

Conformation of beta-hairpins in protein structures. A systematic classification with applications to modelling by homology, electron density fitting and protein engineering

A systematic classification of beta-hairpin structures which takes into account the polypeptide chain length and hydrogen bonding between the two antiparallel beta-strands is described. We have used this classification of beta-hairpin structures and their specific sequence pattern to derive rules which demonstrate its usefulness in assisting modelling beta-hairpins. These rules can be applied to comparative model building, modelling into electron density and in the prediction of conformation of beta-hairpins to aid protein engineering.     

Measurement of the specific and non-specific binding energies of Mg2+ to RNA

Determining the non-specific and specific electrostatic contributions of magnesium binding to RNA is a challenging problem. We introduce a single-molecule method based on measuring the folding energy of a native RNA in magnesium and at its equivalent sodium concentration. The latter is defined so that the folding energy in sodium equals the non-specific electrostatic contribution in magnesium. The sodium equivalent can be estimated according to the empirical 100/1 rule (1 M NaCl is equivalent to 10 mM MgCl₂), which is a good approximation for most RNAs. The method is applied to an RNA three-way junction (3WJ) that contains specific Mg²⁺ binding sites and misfolds into a double hairpin structure without binding sites. We mechanically pull the RNA with optical tweezers and use fluctuation theorems to determine the folding energies of the native and misfolded structures in magnesium (10 mM MgCl₂) and at the equivalent sodium condition (1 M NaCl). While the free energies of the misfolded structure are equal in magnesium and sodium, they are not for the native structure, the difference being due to the specific binding energy of magnesium to the 3WJ, which equals $\text{Δ}G?$ 10 kcal/mol. Besides stabilizing the 3WJ, Mg²⁺ also kinetically rescues it from the misfolded structure over timescales of tens of seconds in a force-dependent manner. The method should generally be applicable to determine the specific binding energies of divalent cations to other tertiary RNAs.     

A NMR strategy to unambiguously distinguish nucleic acid hairpin and duplex conformations applied to a Xist RNA A-repeat

All RNA sequences that fold into hairpins possess the intrinsic potential to form intermolecular duplexes because of their high self-complementarity. The thermodynamically more stable duplex conformation is favored under high salt conditions and at high RNA concentrations, posing a challenging problem for structural studies of small RNA hairpin conformations. We developed and applied a novel approach to unambiguously distinguish RNA hairpin and duplex conformations for the structural analysis of a Xist RNA A-repeat. Using a combination of a quantitative HNN-COSY experiment and an optimized double isotope-filtered NOESY experiment we could define the conformation of the 26-mer A-repeat RNA. In contrast to a previous secondary structure prediction of a double hairpin structure, the NMR data show that only the first predicted hairpin is formed, while the second predicted hairpin mediates dimerization of the A-repeat by duplex formation with a second A-repeat. The strategy employed here will be generally applicable to identify and quantify populations of hairpin and duplex conformations and to define RNA folding topology from inter- and intra-molecular base-pairing patterns.     

Optimal computer folding of large RNA sequences using thermodynamics and auxiliary information.

This paper presents a new computer method for folding an RNA molecule that finds a conformation of minimum free energy using published values of stacking and destabilizing energies. It is based on a dynamic programming algorithm from applied mathematics, and is much more efficient, faster, and can fold larger molecules than procedures which have appeared up to now in the biological literature. Its power is demonstrated in the folding of a 459 nucleotide immunoglobulin gamma 1 heavy chain messenger RNA fragment. We go beyond the basic method to show how to incorporate additional information into the algorithm. This includes data on chemical reactivity and enzyme susceptibility. We illustrate this with the folding of two large fragments from the 16S ribosomal RNA of Escherichia coli.     

The influence of junction conformation on RNA cleavage by the hairpin ribozyme in its natural junction form

In the natural form of the hairpin ribozyme the two loop-carrying duplexes that comprise the majority of essential bases for activity form two adjacent helical arms of a four-way RNA junction. In the present work we have manipulated the sequence around the junction in a way known to perturb the global folding properties. We find that replacement of the junction by a different sequence that has the same conformational properties as the natural sequence gives closely similar reaction rate and Arrhenius activation energy for the substrate cleavage reaction. By comparison, rotation of the natural sequence in order to alter the three-dimensional folding of the ribozyme leads to a tenfold reduction in the kinetics of cleavage. Replacement with the U1 four-way junction that is resistant to rotation into the antiparallel structure required to allow interaction between the loops also gives a tenfold reduction in cleavage rate. The results indicate that the conformation of the junction has a major influence on the catalytic activity of the ribozyme. The results are all consistent with a role for the junction in the provision of a framework by which the loops are presented for interaction in order to create the active form of the ribozyme.     

From structure prediction to genomic screens for novel non-coding RNAs

Non-coding RNAs (ncRNAs) are receiving more and more attention not only as an abundant class of genes, but also as regulatory structural elements (some located in mRNAs). A key feature of RNA function is its structure. Computational methods were developed early for folding and prediction of RNA structure with the aim of assisting in functional analysis. With the discovery of more and more ncRNAs, it has become clear that a large fraction of these are highly structured. Interestingly, a large part of the structure is comprised of regular Watson-Crick and GU wobble base pairs. This and the increased amount of available genomes have made it possible to employ structure-based methods for genomic screens. The field has moved from folding prediction of single sequences to computational screens for ncRNAs in genomic sequence using the RNA structure as the main characteristic feature. Whereas early methods focused on energy-directed folding of single sequences, comparative analysis based on structure preserving changes of base pairs has been efficient in improving accuracy, and today this constitutes a key component in genomic screens. Here, we cover the basic principles of RNA folding and touch upon some of the concepts in current methods that have been applied in genomic screens for de novo RNA structures in searches for novel ncRNA genes and regulatory RNA structure on mRNAs. We discuss the strengths and weaknesses of the different strategies and how they can complement each other.     

A structural database for k-turn motifs in RNA

The kink-turn (k-turn) is a common structural motif in RNA that introduces a tight kink into the helical axis. k-turns play an important architectural role in RNA structures and serve as binding sites for a number of proteins. We have created a database of known and postulated k-turn sequences and three-dimensional (3D) structures, available via the internet. This site provides (1) a database of sequence and structure, as a resource for the RNA community, and (2) a tool to enable the manipulation and comparison of 3D structures where known.     

Estimating the Contributions of Selection and Self-Organization in RNA Secondary Structure

In addition to characteristic structural properties imposed by evolutionary modification, evolved, single-stranded RNAs also display characteristic structural properties imposed by intrinsic physical constraints on RNA polymer folding. The balance of intrinsic and functionally selected characters in the folded conformation of evolved secondary structures was determined by comparing the predicted secondary structures of evolved and unevolved (random) RNA sequences. Though evolved conformations are significantly more ordered than conformations of random-sequence RNA, this analysis demonstrates that the majority of conformational order within evolved structures results not from evolutionary optimization but from constraints imposed by rules intrinsic to RNA polymer folding. Received: 25 November 1998 / Accepted: 12 February 1999     

A computational approach to identify genes for functional RNAs in genomic sequences

Currently there is no successful computational approach for identification of genes encoding novel functional RNAs (fRNAs) in genomic sequences. We have developed a machine learning approach using neural networks and support vector machines to extract common features among known RNAs for prediction of new RNA genes in the unannotated regions of prokaryotic and archaeal genomes. The Escherichia coli genome was used for development, but we have applied this method to several other bacterial and archaeal genomes. Networks based on nucleotide composition were 80–90% accurate in jackknife testing experiments for bacteria and 90–99% for hyperthermophilic archaea. We also achieved a significant improvement in accuracy by combining these predictions with those obtained using a second set of parameters consisting of known RNA sequence motifs and the calculated free energy of folding. Several known fRNAs not included in the training datasets were identified as well as several hundred predicted novel RNAs. These studies indicate that there are many unidentified RNAs in simple genomes that can be predicted computationally as a precursor to experimental study. Public access to our RNA gene predictions and an interface for user predictions is available via the web.     

Understanding the Role of Three-Dimensional Topology in Determining the?Folding Intermediates of Group I Introns

Many RNA molecules exert their biological function only after folding to unique three-dimensional structures. For long, noncoding RNA molecules, the complexity of finding the native topology can be a major impediment to correct folding to the biologically active structure. An RNA molecule may fold to a near-native structure but not be able to continue to the correct structure due to a topological barrier such as crossed strands or incorrectly stacked helices. Achieving the native conformation thus requires unfolding and refolding, resulting in a long-lived intermediate. We investigate the role of topology in the folding of two phylogenetically related catalytic group I introns, the Twort and Azoarcus group I ribozymes. The kinetic models describing the Mg²⁺-mediated folding of these ribozymes were previously determined by time-resolved hydroxyl (⋅OH) radical footprinting. Two intermediates formed by parallel intermediates were resolved for each RNA. These data and analytical ultracentrifugation compaction analyses are used herein to constrain coarse-grained models of these folding intermediates as we investigate the role of nonnative topology in dictating the lifetime of the intermediates. Starting from an ensemble of unfolded conformations, we folded the RNA molecules by progressively adding native constraints to subdomains of the RNA defined by the ⋅OH time-progress curves to simulate folding through the different kinetic pathways. We find that nonnative topologies (arrangement of helices) occur frequently in the folding simulations despite using only native constraints to drive the reaction, and that the initial conformation, rather than the folding pathway, is the major determinant of whether the RNA adopts nonnative topology during folding. From these analyses we conclude that biases in the initial conformation likely determine the relative flux through parallel RNA folding pathways.     

SSCP primer design based on single-strand DNA structure predicted by a DNA folding program.

To predict alterations in single-strand DNA mobility in non-denaturing electrophoretic gels, Zuker's RNA folding program was modified. Energy files utilized by the LRNA RNA folding algorithm were modified to emulate folding of single-strand DNA. Energy files were modified to disallow G-T base pairing. Stacking energies were corrected for DNA thermodynamics. Constraints on loop nucleotide sequences were removed. The LRNA RNA folding algorithm using the DNA fold energy files was applied to predict folding of PCR generated single-strand DNA molecules from polymorphic human ALDH2 and TPH alleles. The DNA-Fold version 1.0 program was used to design primers to create and abolish SSCP mobility shifts. Primers were made that add a 5' tag sequence or alter complementarity to an internal sequence. Differences in DNA secondary structure were assessed by SSCP analysis and compared to single-strand DNA secondary structure predictions. Results demonstrate that alterations in single-strand DNA conformation may be predicted using DNA-Fold 1.0.     

Effect of monovalent cation-induced telomeric DNA structure on the binding of Oxytricha telomeric protein.

Oligonucleotides bearing 4 repeats of telomeric deoxyguanosine-rich sequence undergo a monovalent cation-induced transition to a folded conformation with G-G base pairs, modeled as a 'G-quartet' structure. We have now deduced the rates of folding and unfolding of d(TTTTGGGG)4, which has four repeats of the Oxytricha telomeric DNA sequence. The estimated average values of delta G for the folded form at 37 degrees C are -2.2 kcal/mol and -4.7 kcal/mol in 50 mM na+ and K+, respectively. The fully folded DNA is not recognized by the Oxytricha telomere-binding protein; the substrate for protein binding has properties consistent with its being partly or fully unfolded. In confirmation of this conclusion, prevention of DNA folding by methylation enables the protein to bind as rapidly in the presence of monovalent cations as in their absence. The slow unfolding (t1/2 = 4 hr and 18 hr at 37 degrees C in Na+ and K+, respectively) of the DNA suggests that such structures would be long-lived if they formed in vivo, unless they can be actively unfolded. The inability of the telomere-binding protein to bind the stable, folded form of the 4-repeat telomeric sequence is a problem that may be circumvented in vivo by avoiding four single-stranded repeats.     

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather than Electrostatic Repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg²⁺ concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.     

Structures of helical junctions in nucleic acids

Our knowledge of the architectural principles of nucleic acid junctions has seen significant recent advances. The conformation of DNA junctions is now well understood, and this provides a new basis for the analysis of important structural elements in RNA. The most significant new data have come from X-ray crystallography of four-way DNA junctions; incidentally showing the great importance of serendipity in science, since none of the three groups had deliberately set out to crystallise a junction. Fortunately the results confirm, and of course extend, the earlier conformational studies of DNA junctions in almost every detail. This is important, because it means that these methods can be applied with greater confidence to new systems, especially in RNA. Methods like FRET, chemical probing and even the humble polyacrylamide gel can be rapid and very powerful, allowing the examination of a large number of sequence variants relatively quickly. Molecular modelling in conjunction with experiments is also a very important component of the general approach. Ultimately crystallography provides the gold standard for structural analysis, but the other, simple approaches have considerable value along the way. At the beginning of this review I suggested two simple folding principles for branched nucleic acids, and it is instructive to review these in the light of recent data. In brief, these were the tendency for pairwise coaxial stacking of helical arms, and the importance of metal ion interactions in the induction of folding. We see that both are important in a wide range of systems, both in DNA and RNA. The premier example is the four-way DNA junction, which undergoes metal ion-induced folding into the stacked X-structure that is based on coaxial stacking of arms. As in many systems, there are two alternative ways to achieve this depending on the choice of stacking partners. Recent data reveal that both forms often exist in a dynamic equilibrium, and that the relative stability of the two conformers depends upon base sequence extending a significant distance from the junction. The three-way junction has provided a good test of the folding principles. Perfect three-way (3H) DNA junctions seem to defy these principles in that they appear reluctant to undergo coaxial stacking of arms, and exhibit little change in conformation with addition of metal ions. Modelling suggests that such a junction is stereochemically constrained in an extended conformation. However, upon inclusion of a few additional base pairs at the centre (to create a 3HS2 junction for example) the additional stereochemical flexibility allows two arms to undergo coaxial stacking. Such a junction exhibits all the properties consistent with the general folding principles, with ion-induced folding into a form based on pairwise coaxial stacking of arms in one of two different conformers. The three-way junction is therefore very much the exception that proves the rule. It is instructive to compare the folding of corresponding species in DNA and RNA, where we find both similarities and differences. The RNA four-way junction can adopt a structure that is globally similar to the stacked X-structure (Duckett et al. 1995a), and the crystal structure of the DNAzyme shows that the stacked X-conformation can include one helical pair in the A-conformation (Nowakowski et al. 1999). However, modelling suggests that the juxtaposition of strands and grooves will be less satisfactory in RNA, and the higher magnesium ion concentration required to fold the RNA junction indicates a lower stability of the antiparallel form. Perhaps the biggest difference between the properties of the DNA and RNA four-way junctions is the lack of an unstacked structure at low salt concentrations for the RNA species. This clearly reflects a major difference in the electrostatic interactions in the RNA junction. In general the folding of branched DNA provides some good indications on the likely folding of the corresponding RNA species, but caution is required in making the extrapolation because the two polymers are significantly different. A number of studies point to the flexibility and malleability of branched nucleic acids, and this turns out to have particular significance in their interactions with proteins. Proteins such as the DNA junction-resolving enzymes exhibit considerable selectivity for the structure of their substrates, which is still not understood at a molecular level. Despite this, it appears to be universally true that these proteins distort the global, and in some cases at least the local, structure of the junctions. The somewhat perplexing result is that the proteins appear to distort the very property that they recognise. In general it seems that four-way DNA junctions are opened to one extent or another by interaction with proteins. (ABSTRACT TRUNCATED)     

The importance of G.A hydrogen bonding in the metal ion- and protein-induced folding of a kink turn RNA

The kink turn (K-turn) is a common motif in RNA structure, found in many RNA species important in translation, RNA modification and splicing, and the control of gene expression. In general the K-turn comprises a three nucleotide bulge followed by trans sugar-Hoogsteen G·A pairs. The RNA adopts a tightly kinked conformation, and is a common target for binding proteins, exemplified by the L7Ae family. We have measured the rates of association and dissociation for the binding of L7Ae to the Kt-7 kink turn, from which we calculate an affinity of K_D = 10 pM. This high affinity is consistent with the role of this binding as the first stage in the assembly of key functional nucleoproteins such as box C/D snoRNP. Kink-turn RNA undergoes a two-state transition between the kinked conformation, and a more extended structure, and folding into the kinked form is induced by divalent metal ions, or by binding of proteins of the L7Ae class. The K-turn provides an excellent, simple model for RNA folding, which can be dissected at the atomic level. We have analyzed the contributions of the hydrogen bonds that form the G·A pairs to the ion- and protein-induced folding of the K-turn. We find that all four hydrogen bonds are important to the stability of the kinked form of the RNA, and we can now define all the important hydrogen bonding interactions that stabilize the K-turn. The high affinity of L7Ae binding is coupled to the induced folding of the K-turn, allowing some sub-optimal variants to adopt the kinked geometry. However, in all such cases the affinity is lowered, and the results underline the importance of both G·A pairs to the stability of the K-turn.