Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “RNA structure”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43Linked to original sources

The massively parallel genetic algorithm for RNA folding: MIMD implementation and population variation.

A massively parallel Genetic Algorithm (GA) has been applied to RNA sequence folding on three different computer architectures. The GA, an evolution-like algorithm that is applied to a large population of RNA structures based on a pool of helical stems derived from an RNA sequence, evolves this population in parallel. The algorithm was originally designed and developed for a 16384 processor SIMD (Single Instruction Multiple Data) MasPar MP-2. More recently it has been adapted to a 64 processor MIMD (Multiple Instruction Multiple Data) SGI ORIGIN 2000, and a 512 processor MIMD CRAY T3E. The MIMD version of the algorithm raises issues concerning RNA structure data-layout and processor communication. In addition, the effects of population variation on the predicted results are discussed. Also presented are the scaling properties of the algorithm from the perspective of the number of physical processors utilized and the number of virtual processors (RNA structures) operated upon.

Algorithms↗

RNAdraw: an integrated program for RNA secondary structure calculation and analysis under 32-bit Microsoft Windows.

The use of algorithms for calculation and analysis of RNA secondary structures has largely been limited to mainframe computers, mainly due to the 16-bit memory restrictions imposed by MS-DOS. The program presented here, RNAdraw, moves calculations to the 32-bit Microsoft Windows environments with an intuitive user interface with extensive viewing, editing and printing possibilities. The algorithms for secondary structure/basepair probability matrix/heat curve calculation have been ported directly from a 32-bit Unix environment. RNAdraw also offers novel features such as the options to edit energy parameters, extract structures of different probability levels, create de novo secondary structures interactively, and combine viewing of structures and basepair probabilities.

Algorithms↗

A retroviral RNA secondary structure required for efficient initiation of reverse transcription.

Genetic evidence is presented which suggests the existence of an important structural element in the 5' noncoding region of avian retrovirus RNA. The proposed structure, which we term the U5-leader stem, is composed of sequences in the middle of U5 and in the leader, flanking the primer-binding site. U5 and leader mutations which would disrupt this structure caused a partial replication defect. However, nucleotide substitutions in the leader, which would structurally compensate for a U5 deletion mutation, restored normal replication. Analysis of replication intermediates of viruses with the above mutations suggests that the U5-leader stem is required for efficient DNA synthesis in vivo and for initiation of DNA synthesis from the tRNA(Trp) primer in melittin-activated virions. However, this structure does not appear to be required for binding of the tRNA(Trp) primer to viral RNA. These results support a role for the U5-leader stem structure, independent of its primary sequence, in the initiation of retroviral replication.

Animals↗

Monitoring RNA base structure and dynamics using site-directed spin labeling.

Site-directed spin labeling utilizes site-specific attachment of a stable nitroxide radical to probe the structure and dynamics of macromolecules. In the present study, a 4-thiouridine base is introduced at each of six different positions in a 23-nucleotide RNA molecule. The 4-thiouridine derivatives were subsequently modified with one of three methanethiosulfonate nitroxide reagents to introduce a spin label at specific sites. The electron paramagnetic resonance spectra of the labeled RNAs were analyzed in terms of nitroxide motion and the RNA solution structure. At a base-paired site in the RNA helix, where the nitroxide has weak or no local interactions, motion of the nitroxide is apparently dominated by rotation about bonds within the probe. The motion is similar to that found for a structurally related probe on helical sites in proteins, suggesting a similar mode of motion. At other sites that are hydrogen bonded and stacked within the helix, local interactions within the RNA molecule modulate the nitroxide motion in a manner consistent with expectations based on the known structure. For a base that is not structurally constrained, the mobility is higher than at any other site, presumably due to motion of the base itself. These results demonstrate the general utility of the 4-thiouridine/methanethiosulfonate coupling method to introduce nitroxide spin labels into RNA and the ability of the resulting label to probe local structure and dynamics.

Base Sequence↗

Saccharomyces SRP RNA secondary structures: a conserved S-domain and extended Alu-domain.

The contribution made by the RNA component of signal recognition particle (SRP) to its function in protein targeting is poorly understood. We have generated a complete secondary structure for Saccharomyces cerevisiae SRP RNA, scR1. The structure conforms to that of other eukaryotic SRP RNAs. It is rod-shaped with, at opposite ends, binding sites for proteins required for the SRP functions of signal sequence recognition (S-domain) and translational elongation arrest (Alu-domain). Micrococcal nuclease digestion of purified S. cerevisiae SRP separated the S-domain of the RNA from the Alu-domain as a discrete fragment. The Alu-domain resolved into several stable fragments indicating a compact structure. Comparison of scR1 with SRP RNAs of five yeast species related to S. cerevisiae revealed the S-domain to be the most conserved region of the RNA. Extending data from nuclease digestion with phylogenetic comparison, we built the secondary structure model for scR1. The Alu-domain contains large extensions, including a sequence with hallmarks of an expansion segment. Evolutionarily conserved bases are placed in the Alu- and S-domains as in other SRP RNAs, the exception being an unusual GU(4)A loop closing the helix onto which the signal sequence binding Srp54p assembles (domain IV). Surprisingly, several mutations within the predicted Srp54p binding site failed to disrupt SRP function in vivo. However, the strength of the Srp54p-scR1 and, to a lesser extent, Sec65p-scR1 interaction was decreased in these mutant particles. The availability of a secondary structure for scR1 will facilitate interpretation of data from genetic analysis of the RNA.

Alu Elements↗

Alignment of possible secondary structures in multiple RNA sequences using simulated annealing.

Multiple sequence alignment has been a useful technique for identifying RNA secondary structures. In this paper, an algorithm for aligning multiple RNA sequences to identify possible secondary structure is presented. In this algorithm, dot matrices generated from intra-sequence comparisons are used to obtain possible common secondary structures. A hit probability for dot matrices is calculated and a score function based on this hit probability is defined. Simulated annealing is applied to optimize the score function. The solution set of multiple sequence alignment is introduced, and the effects on the solution set of increasing the number of alignment gaps and the alignment length are analyzed. Several additional strategies to reduce simulated annealing time are applied. A method is applied to reduce the computation time based on the solution set. Also, an optimized transition rule, double shuffle, which moves two positions in a sequence with each iteration, is applied to increase the rate of convergence. This algorithm was used to find possible common secondary structures in RNA sequences.

Algorithms↗

Evidence for control of splicing by alternative RNA secondary structures in Dipteran homothorax pre-mRNA.

In a recent study that identified highly evolutionary conserved sequences in three genomes of Diptera species we described an ultraconserved element found at an internal exon-intron junction of the Drosophila melanogaster homothorax (hth) gene that appeared to be involved in the control of hth pre-mRNA splicing. We also discussed a possible role of RNA secondary structure at this site in the regulation of hth pre-mRNA splicing. In this report we identify a shorter evolutionary conserved intronic element within the hth gene that is located downstream of the first element and has sequence complementarity to it. We demonstrate that intramolecular interactions between these two elements would give rise to alternative RNA secondary structures, which in turn may result in differential control of homothorax pre-mRNA splicing. We also provide additional comparative genomic data from several newly available insect genomes supporting our original conclusion that these conserved elements are important in the post-transcriptional regulation of homothorax gene expression in Diptera.

Alternative Splicing↗

Cleavage of highly structured viral RNA molecules by combinatorial libraries of hairpin ribozymes. The most effective ribozymes are not predicted by substrate selection rules.

Combinatorial libraries of hairpin ribozymes representing all possible cleavage specificities (>10(5)) were used to evaluate all ribozyme cleavage sites within a large (4.2-kilobase) and highly structured viral mRNA, the 26 S subgenomic RNA of Sindbis virus. The combinatorial approach simultaneously accounts for target site structure and dynamics, together with ribozyme folding, and the sequences that result in a ribozyme-substrate complex with maximal activity. Primer extension was used to map and rank the relative activities of the ribozyme pool against individual sites and revealed two striking findings. First, only a small fraction of potential recognition sites are effectively cleaved (activity-selected sites). Second, nearly all of the most effectively cleaved sites deviated substantially from the established consensus selection rules for the hairpin ribozyme and were not predicted by examining the sequence, or through the use of computer-assisted predictions of RNA secondary structure. In vitro selection methods were used to isolate ribozymes with increased activity against substrates that deviate from the GUC consensus sequence. trans-Acting ribozymes targeting nine of the activity-selected sites were synthesized, together with ribozymes targeting four sites with a perfect match to the cleavage site consensus (sequence-selected sites). Activity-selected ribozymes have much higher cleavage activity against the long, structured RNA molecules than do sequence-selected ribozymes, although the latter are effective in cleaving oligoribonucleotides, as predicted. These results imply that, for Sindbis virus 26 S RNA, designing ribozymes based on matches to the consensus sequence may be an ineffective strategy.

Base Sequence↗

Evidence for differential ribonucleoprotein complex assembly in vitro on the 5'-untranslated region of the human IGF-IR transcript.

The type I insulin-like growth factor receptor (IGF-IR) plays a key role in the control of cellular proliferation and survival. The human IGF-IR transcript is characterized by an unusually long 1038 nucleotide 5'-untranslated region (5'-UTR). We hypothesized that the contribution of this complex 5'-untranslated RNA sequence to the post-transcriptional regulation of IGF-IR expression would involve a dynamic interplay between RNA structure and specific RNA-binding proteins. Here we have detected and characterized a diverse series of regulatory proteins binding the IGF-IR 5'-UTR under disparate conditions. One pair of proteins ( approximately 42/38 kDa) binds readily to the intact 5'-UTR, which is predicted to adopt a highly base-paired, highly favorable (dG=-498 kcal/mol) three-domain structure. Another protein(s) (p20*) specifically induces formation of a novel RNA structure from within the initial 209 nucleotides of the nascent IGF-IR transcript, but fails to UV crosslink to this RNA sequence. A third group of proteins recognizes and binds the IGF-IR 5'-UTR under highly stringent conditions, but only after higher-ordered RNA structure has been disrupted. Our in vitro results indicate that the IGF-IR 5'-UTR may exist in at least three distinct states, and we propose that interconversion between these states might take place in vivo and differentially alter IGF-IR transcript utilization.

5' Untranslated Regions↗

Effects of secondary structures in RNA on interlocking probabilities.

In 1967 Wang and Schwartz reported on the formation of interlocked rings between linear or circular DNA molecules by the enzyme topoisomerase. We propose viewing the secondary structured loop in RNA (or single stranded DNA) as analogous to a circular DNA molecule. Formation of a catenane between such an RNA loop with a DNA molecule may constitute a probe of the secondary and general three dimensional structure of the RNA molecule. The experimental results may be compared with the theoretical calculation. We suggest here a method for estimating linkage probabilities and calculate them for several cases for which secondary structures of the RNA have been proposed.

DNA Topoisomerases, Type I↗

Compensatory mutations cause excess of antagonistic epistasis in RNA secondary structure folding.

BACKGROUND: The rate at which fitness declines as an organism's genome accumulates random mutations is an important variable in several evolutionary theories. At an intuitive level, it might seem natural that random mutations should tend to interact synergistically, such that the rate of mean fitness decline accelerates as the number of random mutations is increased. However, in a number of recent studies, a prevalence of antagonistic epistasis (the tendency of multiple mutations to have a mitigating rather than reinforcing effect) has been observed. RESULTS: We studied in silico the net amount and form of epistatic interactions in RNA secondary structure folding by measuring the fraction of neutral mutants as a function of mutational distance d. We found a clear prevalence of antagonistic epistasis in RNA secondary structure folding. By relating the fraction of neutral mutants at distance d to the average neutrality at distance d, we showed that this prevalence derives from the existence of many compensatory mutations at larger mutational distances. CONCLUSIONS: Our findings imply that the average direction of epistasis in simple fitness landscapes is directly related to the density with which fitness peaks are distributed in these landscapes.

Animals↗

RNA tertiary structure determination: NOE pathways construction by tabu search.

MOTIVATION: Liquid state nuclear magnetic resonance (NMR) spectroscopy has now been well established as a method for RNA tertiary structure determination. Most of the steps involved in the determination of RNA molecules are performed using computer programs. They however, do not apply to resonance assignment being the starting point of the whole procedure. We propose a tabu search algorithm as a tool for automating this step. Nuclear overhause effect (NOE) pathway, which determines the assignment, is constructed during an analysis of possible connections between resonances within aromatic/anomeric region of two-dimensional NOESY spectrum resulting from appropriate NMR experiment. RESULTS: Computational tests demonstrate the superior performance of the tabu search algorithm as compared with the exact enumerative approach and genetic procedure applied to the experimental and simulated spectral data for RNA molecules. AVAILABILITY: The software package can be obtained upon request from Marta Szachniuk.

Algorithms↗

Selective recognition and cleavage of RNA loop structures by Ni(II).Xaa-Gly-His metallopeptides.

The recognition and cleavage of tRNAPhe and the TAR RNA of HIV-1 by metallopeptides of the general form Ni(II).Xaa-Gly-His (where Xaa is Gly, Lys, or Arg) were investigated. The results of RNA cleavage analyses suggest that KHSO5- or magnesium monoperoxyphthalate-activated metallopeptides (1) induce nucleobase damage which requires aniline acetate for complete RNA strand scission and (2) selectively target the loops of stem-loop structures of the above-named substrates. In targeting RNA loop regions, the metallopeptides may be sensitive to intraloop structural features, including the overall structural environment of the loop itself and possibly the presence of intraloop hydrogen bonding. Overall, these results suggest that the metallopeptides interact selectively within a loop, in a fashion reminiscent of many RNA binding proteins, instead of targeting RNA single-stranded character alone. These observations further suggest a possible metallopeptide-based strategy for the molecular recognition of native RNA structures and insight with regard to the general features available for ligand binding site discrimination.

Base Sequence↗

[Search for conservative secondary structures of RNA].

We suggest a new algorithm to search a given set of the RNA sequences for conserved secondary structures. The algorithm is based on alignment of the sequences for potential helical strands. This procedure can be used to search for new structured RNAs and new regulatory elements. It is efficient for the genome-scale analysis. The results of various tests run with this algorithm are shown.

Algorithms↗

Kissing of the two predominant hairpin loops in the coxsackie B virus 3' untranslated region is the essential structural feature of the origin of replication required for negative-strand RNA synthesis.

Higher-order RNA structures in the 3' untranslated region (3'UTR) of enteroviruses are thought to play a pivotal role in viral negative-strand RNA synthesis. The structure of the 3'UTR was predicted by thermodynamic calculations using the STAR (structural analysis of RNA) computer program and experimentally verified using chemical and enzymatic probing of in vitro-synthesized RNA. A possible pseudoknot interaction between the 3D polymerase coding sequence and domain Y and a "kissing" interaction between domains X and Y was further studied by mutational analysis, using an infectious coxsackie B3 virus cDNA clone (domain designation as proposed by E. V. Pilipenko, S. V. Maslova, A. N. Sinyakov, and V.I. Agol (Nucleic Acids Res. 20:1739-1745, 1992). The higher-order RNA structure of the 3'UTR appeared to be maintained by an intramolecular kissing interaction between the loops of the two predominant hairpin structures (X and Y) within the 3'UTR. Disturbing this interaction had no effect on viral translation and processing of the polyprotein but exerted a primary effect on viral replication, as was demonstrated in a subgenomic coxsackie B3 viral replicon, in which the capsid P1 region was replaced by the luciferase gene. Mutational analysis did not support the existence of the pseudoknot interaction between hairpin loop Y and the 3D polymerase coding sequence. Based on these experiments, we constructed a three-dimensional model of the 3'UTR of coxsackie B virus that shows the kissing interaction as the essential structural feature of the origin of replication required for its functional competence.

Base Sequence↗

Protein-RNA interactions: structural analysis and functional classes.

A data set of 89 protein-RNA complexes has been extracted from the Protein Data Bank, and the nucleic acid recognition sites characterized through direct contacts, accessible surface area, and secondary structure motifs. The differences between RNA recognition sites that bind to RNAs in functional classes has also been analyzed. Analysis of the complete data set revealed that van der Waals interactions are more numerous than hydrogen bonds and the contacts made to the nucleic acid backbone occur more frequently than specific contacts to nucleotide bases. Of the base-specific contacts that were observed, contacts to guanine and adenine occurred most frequently. The most favored amino acid-nucleotide pairings observed were lysine-phosphate, tyrosine-uracil, arginine-phosphate, phenylalanine-adenine and tryptophan-guanine. The amino acid propensities showed that positively charged and polar residues were favored as expected, but also so were tryptophan and glycine. The propensities calculated for the functional classes showed trends similar to those observed for the complete data set. However, the analysis of hydrogen bond and van der Waal contacts showed that in general proteins complexed with messenger RNA, transfer RNA and viral RNA have more base specific contacts and less backbone contacts than expected, while proteins complexed with ribosomal RNA have less base-specific contacts than the expected. Hence, whilst the types of amino acids involved in the interfaces are similar, the distribution of specific contacts is dependent upon the functional class of the RNA bound.

Amino Acid Motifs↗