Phosphorothioate substitution can substantially alter RNA conformation

Phosphorothioate substitution-interference experiments, routinely used to stereospecifically identify phosphoryl oxygen sites that participate in RNA-ligand binding and RNA-directed catalysis, rest in their interpretation on the untested assumption that substitution does not alter the conformation of the modified molecule from its biologically active state. Using NMR spectroscopy, we have tested this assumption by determining the structural effect of stereospecific phosphorothioate substitution at five positions in an RNA hairpin containing the binding site for bacteriophage MS2 capsid protein. At most sites, substitution has little or no effect, causing minor perturbations in the phosphate backbone and increasing the stacking among nucleotides in the hairpin loop. At one site, however, phosphorothioate substitution causes an unpaired adenine necessary for formation of the capsid protein-RNA complex to loop out of the RNA helix into the major groove. These results indicate that phosphorothioate substitution can substantially alter the conformation of RNA at positions of irregular secondary structure, complicating the use of substitution-interference experiments to study RNA structure and function.     

RNANetMotif: Identifying sequence-structure RNA network motifs in RNA-protein binding sites

RNA molecules can adopt stable secondary and tertiary structures, which are essential in mediating physical interactions with other partners such as RNA binding proteins (RBPs) and in carrying out their cellular functions. In vivo and in vitro experiments such as RNAcompete and eCLIP have revealed in vitro binding preferences of RBPs to RNA oligomers and in vivo binding sites in cells. Analysis of these binding data showed that the structure properties of the RNAs in these binding sites are important determinants of the binding events; however, it has been a challenge to incorporate the structure information into an interpretable model. Here we describe a new approach, RNANetMotif, which takes predicted secondary structure of thousands of RNA sequences bound by an RBP as input and uses a graph theory approach to recognize enriched subgraphs. These enriched subgraphs are in essence shared sequence-structure elements that are important in RBP-RNA binding. To validate our approach, we performed RNA structure modeling via coarse-grained molecular dynamics folding simulations for selected 4 RBPs, and RNA-protein docking for LIN28B. The simulation results, e.g., solvent accessibility and energetics, further support the biological relevance of the discovered network subgraphs.     

Exploring the repertoire of RNA secondary motifs using graph theory; implications for RNA design

Understanding the structural repertoire of RNA is crucial for RNA genomics research. Yet current methods for finding novel RNAs are limited to small or known RNA families. To expand known RNA structural motifs, we develop a two-dimensional graphical representation approach for describing and estimating the size of RNA’s secondary structural repertoire, including naturally occurring and other possible RNA motifs. We employ tree graphs to describe RNA tree motifs and more general (dual) graphs to describe both RNA tree and pseudoknot motifs. Our estimates of RNA’s structural space are vastly smaller than the nucleotide sequence space, suggesting a new avenue for finding novel RNAs. Specifically our survey shows that known RNA trees and pseudoknots represent only a small subset of all possible motifs, implying that some of the ‘missing’ motifs may represent novel RNAs. To help pinpoint RNA-like motifs, we show that the motifs of existing functional RNAs are clustered in a narrow range of topological characteristics. We also illustrate the applications of our approach to the design of novel RNAs and automated comparison of RNA structures; we report several occurrences of RNA motifs within larger RNAs. Thus, our graph theory approach to RNA structures has implications for RNA genomics, structure analysis and design.     

Phylogenetic analysis and evolution of RNase P RNA in proteobacteria.

The secondary structures of the eubacterial RNase P RNAs are being elucidated by a phylogenetic comparative approach. Sequences of genes encoding RNase P RNA from each of the recognized subgroups (alpha, beta, gamma, and delta) of the proteobacteria have now been determined. These sequences allow the refinement, to nearly the base pair level, of the phylogenetic model for RNase P RNA secondary structure. Evolutionary change among the RNase P RNAs was found to occur primarily in four discrete structural domains that are peripheral to a highly conserved core structure. The new sequences were used to examine critically the proposed similarity (C. Guerrier-Takada, N. Lumelsky, and S. Altman, Science 246:1578-1584, 1989) between a portion of RNase P RNA and the "exit site" of the 23S rRNA of Escherichia coli. Phylogenetic comparisons indicate that these sequences are not homologous and that any similarity in the structures is, at best, tenuous.     

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.     

A two unit antisense RNA cassette test system for silencing of target genes.

This communication describes a two unit antisense RNA cassette system for use in gene silencing. Cassettes consist of a recognition unit and an inhibitory unit which are transcribed into a single RNA that carries sequences of non-contiguous complementarity to the chosen target RNA. The recognition unit is designed as a stem-loop for rapid formation of long- lived binding intermediates with target sequences and resembles the major stem-loop of a naturally occurring antisense RNA, CopA. The inhibitory unit consists of either a sequence complementary to a ribosome binding site or of a hairpin ribozyme targeted at a site within the chosen mRNA. The contributions of the individual units to inhibition was assessed using the lacI gene as a target. All possible combinations of recognition and inhibitory units were tested in either orientation. In general, inhibition of lacI expression was relatively low. Fifty per cent inhibition was obtained with the most effective of the constructs, carrying the recognition stem-loop in the antisense orientation and the inhibitory unit with an anti-RBS sequence. Several experiments were performed to assess activities of the RNAs in vitro and in vivo : antisense RNA binding assays, cleavage assays, secondary structure analysis as well as Northern blotting and primer extension analysis of antisense and target RNAs. The problems associated with this antisense RNA approach as well as its potential are discussed with respect to possible optimization strategies.     

Phosphorylation of p68 RNA helicase regulates RNA binding by the C-terminal domain of the protein

We previously reported ATPase, RNA unwinding, and RNA-binding activities of recombinant p68 RNA helicase that was expressed in Escherichia coli. Huang et al. The recombinant protein bound both single-stranded (ss) and double-stranded (ds) RNAs. To further characterize the substrate RNA binding by p68 RNA helicase, we expressed and purified the recombinant N-terminal and C-terminal domains of the protein. RNA-binding property and protein phosphorylation of the recombinant domains of p68 were analyzed. Our data demonstrated that the C-terminal domain of p68 RNA helicase bound ssRNA. More interestingly, the C-terminal domain was a target of protein kinase C (PKC). Phosphorylation of the C-terminal domain of p68 abolished its RNA binding. Based on our observations, we propose that the C-terminal domain is an RNA substrate binding site for p68. The protein phosphorylation by PKC regulates the RNA binding of p68 RNA helicase, which consequently controls the enzymatic activities of the protein.     

Identification of phosphates involved in catalysis by the ribozyme RNase P RNA.

The RNA subunit of ribonuclease P (RNase P RNA) is a catalytic RNA that cleaves precursor tRNAs to generate mature tRNA 5' ends. Little is known concerning the identity and arrangement of functional groups that constitute the active site of this ribozyme. We have used an RNase P RNA-substrate conjugate that undergoes rapid, accurate, and efficient self-cleavage in vitro to probe, by phosphorothioate modification-interference, functional groups required for catalysis. We identify four phosphate oxygens where substitution by sulfur significantly reduces the catalytic rate (50-200-fold). Interference at one site was partially rescued in the presence of manganese, suggesting a direct involvement in binding divalent metal ion cofactors required for catalysis. All sites are located in conserved sequence and secondary structure, and positioned adjacent to the substrate phosphate in a tertiary structure model of the ribozyme-substrate complex. The spatial arrangement of phosphorothioate-sensitive sites in RNase P RNA was found to resemble the distribution of analogous positions in the secondary and potential tertiary structures of other large catalytic RNAs.     

An algorithm for searching RNA motifs in genomic sequences

RNA molecules, which are found in all living cells, fold into characteristic structures that account for their diverse functional activities. Many of these RNA structures consist of a collection of fundamental RNA motifs. The various combinations of RNA basic components form different RNA classes and define their unique structural and functional properties. The availability of many genome sequences makes it possible to search computationally for functional RNAs. Biological experiments indicate that functional RNAs have characteristic RNA structural motifs represented by specific combinations of base pairings and conserved nucleotides in the loop regions. The searching for those well-ordered RNA structures and their homologues in genomic sequences is very helpful for the understanding of RNA-based gene regulation. In this paper, we consider the following problem: given an RNA sequence with a known secondary structure, efficiently determine candidate segments in genomic sequences that can potentially form RNA secondary structures similar to the given RNA secondary structure. Our new bottom-up approach searches all potential stem-loops similar to ones of the given RNA secondary structure first, and then based on located stem-loops, detects potential homologous structural RNAs in genomic sequences.     

A designed RNA selection: establishment of a stable complex between a target and selectant RNA via two coordinated interactions

In this paper, we describe a new method for selecting RNA aptamers that cooperatively bind to two specific sites within a target RNA. We designed a selection system in which two RNAs, a target RNA and a RNA pool, were assembled by employing a pre-organized GAAA tetraloop-11-nt receptor interaction. This allows us to select the binding sequence against a targeted internal loop as well as a linker region optimized for binding of the two binding sites. After the selection, the aptamers bound with dissociation constants in the nanomolar range, thereby forming a stable complex with the target RNA. Thus this method enables identification of aptamers for a specific binding site together with a linker for cooperative binding of the two RNAs. It appears that our new method can be applied generally to select RNAs that adhere tightly to a target RNA via two specific sites. The method can also be applicable for further engineering of both natural and artificial RNAs.     

RBPBind: Quantitative Prediction of Protein-RNA Interactions

There are hundreds of RNA binding proteins in the human genome alone and their interactions with messenger and other RNAs in a cell regulate every step in an RNA's life cycle. To understand this interplay of proteins and RNA it is important to be able to know which protein binds which RNA how strongly and where. Here, we introduce RBPBind, a web-based tool for the quantitative prediction of the interaction of single-stranded RNA binding proteins with target RNAs that fully takes into account the effect of RNA secondary structure on binding affinity. Given a user-specified RNA and a protein selected from a set of several RNA-binding proteins, RBPBind computes their binding curve and effective binding constant. The server also computes the probability that, at a given protein concentration, a protein molecule will bind to any particular nucleotide along the RNA. The sequence specificity of the protein-RNA interaction is parameterized from public RNAcompete experiments and integrated into the recursions of the Vienna RNA package to simultaneously take into account protein binding and RNA secondary structure. We validate our approach by comparison to experimentally determined binding affinities of the HuR protein for several RNAs of different sequence contexts from the literature, showing that integration of raw sequence affinities into RNA secondary structure prediction significantly improves the agreement between computationally predicted and experimentally measured binding affinities. Our resource thus provides a quick and easy way to obtain reliable predicted binding affinities and locations for single-stranded RNA binding proteins based on RNA sequence alone.     

Phylogenetic-comparative analysis of the eukaryal ribonuclease P RNA

Ribonuclease P (RNase P) is the ribonucleoprotein enzyme that cleaves 5'-leader sequences from precursor-tRNAs. Bacterial and eukaryal RNase P RNAs differ fundamentally in that the former, but not the latter, are capable of catalyzing pre-tRNA maturation in vitro in the absence of proteins. An explanation of these functional differences will be assisted by a detailed comparison of bacterial and eukaryal RNase P RNA structures. However, the structures of eukaryal RNase P RNAs remain poorly characterized, compared to their bacterial and archaeal homologs. Hence, we have taken a phylogenetic-comparative approach to refine the secondary structures of eukaryal RNase P RNAs. To this end, 20 new RNase P RNA sequences have been determined from species of ascomycetous fungi representative of the genera Arxiozyma, Clavispora, Kluyveromyces, Pichia, Saccharomyces, Saccharomycopsis, Torulaspora, Wickerhamia, and Zygosaccharomyces. Phylogenetic-comparative analysis of these and other sequences refines previous eukaryal RNase P RNA secondary structure models. Patterns of sequence conservation and length variation refine the minimum-consensus model of the core eukaryal RNA structure. In comparison to bacterial RNase P RNAs, the eukaryal homologs lack RNA structural elements thought to be critical for both substrate binding and catalysis. Nonetheless, the eukaryal RNA retains the main features of the catalytic core of the bacterial RNase P. This indicates that the eukaryal RNA remains intrinsically a ribozyme.     

Messenger RNA isolation using novel PNA analogues

Homo-Thy hetero-oligomer probes composed of trans-4-hydroxy-L-proline based PNA-like (HypNA) monomers and phosphono PNA (pPNA) monomers demonstrated strong binding to complementary poly A+ RNA strands. We used a mixture of chimeric oligomers containing both "linear" and "clamping" PNA-analogues to develop an mRNA isolation procedure and demonstrate the improved recovery of RNA molecules with secondary structure at the 3'end as well as RNAs with short poly A tails.     

Interaction of R17 coat protein with its RNA binding site for translational repression

The interaction between bacteriophage R17 coat protein and its RNA binding site for translational repression was studied as an example of a sequence-specific RNA-protein interaction. A nitrocellulose filter retention assay is used to demonstrate equimolar binding between the coat protein and a synthetic 21 nucleotide RNA fragment. The Kd at 2 degrees C in a buffer containing 0.19 M salt is about 1 nM. The relatively weak ionic strength dependence of Ka and a delta H = -19 kcal/mole indicates that most of the binding free energy is due to non-electrostatic interactions. Since a variety of RNAs failed to compete with the 21 nucleotide fragment for coat protein binding, the interaction appears highly sequence specific. We have synthesized more than 30 different variants of the binding site sequence in order to identify the portions of the RNA molecule which are important for protein binding. Out of the five single stranded residues examined, four were essential for protein binding whereas the fifth could be replaced by any nucleotide. One variant was found to bind better than the wild type sequence. Substitution of nucleotides which disrupted the secondary structure of the binding fragment resulted in very poor binding to the protein. These data indicated that there are several points of contact between the RNA and the protein and the correct hairpin secondary structure of the RNA is essential for protein binding.