Quantifying the energetic interplay of RNA tertiary and secondary structure interactions

To understand the RNA-folding problem, we must know the extent to which RNA structure formation is hierarchical (tertiary folding of preformed secondary structure). Recently, nuclear magnetic resonance (NMR) spectroscopy was used to show that Mg2+-dependent tertiary interactions force secondary structure rearrangement in the 56-nt tP5abc RNA, a truncated subdomain of the Tetrahymena group I intron. Here we combine mutagenesis with folding computations, nondenaturing gel electrophoresis, high-resolution NMR spectroscopy, and chemical-modification experiments to probe further the energetic interplay of tertiary and secondary interactions in tP5abc. Point mutations predicted to destabilize the secondary structure of folded tP5abc greatly disrupt its Mg2+-dependent folding, as monitored by nondenaturing gels. Imino proton assignments and sequential NOE walks of the two-dimensional NMR spectrum of one of the tP5abc mutants confirm the predicted secondary structure, which does not change in the presence of Mg2+. In contrast to these data on tP5abc, the same point mutations in the context of the P4-P6 domain (of which P5abc is a subdomain) shift the Mg2+ dependence of P4-P6 folding only moderately, and dimethyl sulfate (DMS) modification experiments demonstrate that Mg2+ does cause secondary structure rearrangement of the P4-P6 mutants' P5abc subdomains. Our data provide experimental support for two simple conclusions: (1) Even single point mutations at bases involved only in secondary structure can be enough to tip the balance between RNA tertiary and secondary interactions. (2) Domain context must be considered in evaluating the relative importance of tertiary and secondary contributions. This tertiary/secondary interplay is likely relevant to the folding of many large RNA and to bimolecular snRNA-snRNA and snRNA-intron RNA interactions.     

Topological constraints: using RNA secondary structure to model 3D conformation, folding pathways, and dynamic adaptation

Accompanying recent advances in determining RNA secondary structure is the growing appreciation for the importance of relatively simple topological constraints, encoded at the secondary structure level, in defining the overall architecture, folding pathways, and dynamic adaptability of RNA. A new view is emerging in which tertiary interactions do not define RNA 3D structure, but rather, help select specific conformers from an already narrow, topologically pre-defined conformational distribution. Studies are providing fundamental insights into the nature of these topological constraints, how they are encoded by the RNA secondary structure, and how they interplay with other interactions, breathing new meaning to RNA secondary structure. New approaches have been developed that take advantage of topological constraints in determining RNA backbone conformation based on secondary structure, and a limited set of other, easily accessible constraints. Topological constraints are also providing a much-needed framework for rationalizing and describing RNA dynamics and structural adaptation. Finally, studies suggest that topological constraints may play important roles in steering RNA folding pathways. Here, we review recent advances in our understanding of topological constraints encoded by the RNA secondary structure.     

Effects of osmolytes on RNA secondary and tertiary structure stabilities and RNA-Mg2+ interactions

Osmolytes are small organic molecules accumulated by cells in response to osmotic stress. Although their effects on protein stability have been studied, there has been no systematic documentation of their influence on RNA. Here, the effects of nine osmolytes on the secondary and tertiary structure stabilities of six RNA structures of differing complexity and stability have been surveyed. Using thermal melting analysis, m-values (change in ΔG° of RNA folding per molal concentration of osmolyte) have been measured. All the osmolytes destabilize RNA secondary structure, although to different extents, probably because they favor solubilization of base surfaces. Osmolyte effects on tertiary structure, however, can be either stabilizing or destabilizing. We hypothesize that the stabilizing osmolytes have unfavorable interactions with the RNA backbone, which becomes less accessible to solvent in most tertiary structures. Finally, it was found that as a larger fraction of the negative charge of an RNA tertiary structure is neutralized by hydrated Mg²⁺, the RNA becomes less responsive to stabilizing osmolytes and may even be destabilized. The natural selection of osmolytes as protective agents must have been influenced by their effects on the stabilities of functional RNA structures, though in general, the effects of osmolytes on RNA and protein stabilities do not parallel each other. Our results also suggest that some osmolytes can be useful tools for studying intrinsically unstable RNA folds and assessing the mechanisms of Mg²⁺-induced RNA stabilization.     

RNA structure prediction from evolutionary patterns of nucleotide composition

Structural elements in RNA molecules have a distinct nucleotide composition, which changes gradually over evolutionary time. We discovered certain features of these compositional patterns that are shared between all RNA families. Based on this information, we developed a structure prediction method that evaluates candidate structures for a set of homologous RNAs on their ability to reproduce the patterns exhibited by biological structures. The method is named SPuNC for ‘Structure Prediction using Nucleotide Composition’. In a performance test on a diverse set of RNA families we demonstrate that the SPuNC algorithm succeeds in selecting the most realistic structures in an ensemble. The average accuracy of top-scoring structures is significantly higher than the average accuracy of all ensemble members (improvements of more than 20% observed). In addition, a consensus structure that includes the most reliable base pairs gleaned from a set of top-scoring structures is generally more accurate than a consensus derived from the full structural ensemble. Our method achieves better accuracy than existing methods on several RNA families, including novel riboswitches and ribozymes. The results clearly show that nucleotide composition can be used to reveal the quality of RNA structures and thus the presented technique should be added to the set of prediction tools.     

Insights into the structure and substrate interactions of the P-glycoprotein multidrug transporter from spectroscopic studies

The P-glycoprotein multidrug transporter is a 170-kDa efflux pump which exports a diverse group of natural products, chemotherapeutic drugs, and hydrophobic peptides across the plasma membrane, driven by ATP hydrolysis. The transporter has been proposed to interact with its drug substrates within the membrane environment; however, much remains to be learned about the nature and number of the drug binding site(s). The two nucleotide binding domains are responsible for ATP binding and hydrolysis, which is coupled to drug movement across the membrane. In recent years, P-glycoprotein has been purified and functionally reconstituted in amounts large enough to allow biophysical studies. The use of spectroscopic techniques has led to insights into both its secondary and tertiary structure, and its interaction with nucleotides and drugs. In this review, we will summarise what has been learned by application to purified P-glycoprotein of fluorescence spectroscopy, circular dichroism spectroscopy and infra-red spectroscopy.     

Circles: automating the comparative analysis of RNA secondary structure

SUMMARY: Circles is a program for inferring RNA secondary structure using maximum weight matching. The program can read in an alignment in FASTA, ClustalW, or NEXUS format, compute a maximum weight matching, and export one or more secondary structures in various file formats. AVAILABILITY: The program is available at no cost from http://taxonomy.zoology.gla.ac.uk/rod/circles/ and requires Windows 95/98/NT. CONTACT: r.page@bio.gla.ac.uk     

Predicting protein folding pathways at the mesoscopic level based on native interactions between secondary structure elements

Background  

Since experimental determination of protein folding pathways remains difficult, computational techniques are often used to simulate protein folding. Most current techniques to predict protein folding pathways are computationally intensive and are suitable only for small proteins.     

RNAMotif, an RNA secondary structure definition and search algorithm

RNA molecules fold into characteristic secondary and tertiary structures that account for their diverse functional activities. Many of these RNA structures are assembled from a collection of RNA structural motifs. These basic building blocks are used repeatedly, and in various combinations, to form different RNA types and define their unique structural and functional properties. Identification of recurring RNA structural motifs will therefore enhance our understanding of RNA structure and help associate elements of RNA structure with functional and regulatory elements. Our goal was to develop a computer program that can describe an RNA structural element of any complexity and then search any nucleotide sequence database, including the complete prokaryotic and eukaryotic genomes, for these structural elements. Here we describe in detail a new computational motif search algorithm, RNAMotif, and demonstrate its utility with some motif search examples. RNAMotif differs from other motif search tools in two important aspects: first, the structure definition language is more flexible and can specify any type of base–base interaction; second, RNAMotif provides a user controlled scoring section that can be used to add capabilities that patterns alone cannot provide.     

Solution structure and RNA interactions of the RNA recognition motif from eukaryotic translation initiation factor 4B

Eukaryotic initiation factor 4B (eIF4B) is a multidomain protein with a range of activities that serves primarily to promote association of messenger RNA to the 40S ribosomal subunit during translation initiation. We report here the solution structure of the eIF4B RNA recognition motif (RRM) domain. It adopts a classical RRM fold, with a beta alpha beta beta alpha beta topology. The most striking difference with other RRM structures is in the disposition of loop 3, which connects the beta 2 and beta 3 strands and is implicated in RNA recognition. This loop folds down against the body of the RRM and exhibits restricted motion on a milli- to microsecond time scale. Although it contributes to a large basic patch on the RNA binding surface, it does not protrude out from the domain as observed in other RRM structures, possibly implying a different mode of RNA binding. On its own, the core RRM domain provides only a relative weak interaction with RNA targets and appears to require extensions at the N- and C-terminus for high-affinity binding.     

Measuring the thermodynamics of RNA secondary structure formation

The thermodynamics of RNA secondary structure formation in small model systems provides a database for predicting RNA structure from sequence. Methods for making these measurements are reviewed with emphasis on optical methods and treatment of experimental errors. Analysis of experimental results in terms of simple nearest-neighbor models is presented. Some measured sequence dependences of non-Watson-Crick motifs are discussed. © 1998 John Wiley & Sons, Inc. Biopoly 44: 309–319, 1997     

Nucleotide sequence, secondary structure and evolution of the 5S ribosomal RNA from five bacterial species

The nucleotide sequences of the 5S ribosomal RNAs of the bacteria Agrobacterium tumefaciens, Alcaligenes faecalis, Pseudomonas cepacia, Aquaspirillum serpens and Acinetobacter calcoaceticus have been determined. The sequences fit in a generally accepted model for 5S RNA secondary structure. However, a closer comparative examination of these and other bacterial 5S RNA primary structures reveals the potential of additional base pairing and of multiple equilibria between a set of slightly different alternative secondary structures in one area of the molecule. The phylogenetic position of the examined bacteria is derived from a 5S RNA sequence alignment by a clustering method and compared with the position derived on the basis of 16S ribosomal RNA oligonucleotide catalogs.