RNA loop-loop interactions as dynamic functional motifs

RNA loop-loop interactions are frequently used to trigger initial recognition between two RNA molecules. In this review, we present selected well-documented cases that illustrate the diversity of biological processes using RNA loop-loop recognition properties. The first one is related to natural antisense RNAs that play a variety of regulatory functions in bacteria and their extra-chromosomal elements. The second one concerns the dimerization of HIV-1 genomic RNA, which is responsible for the encapsidation of a diploid RNA genome. The third one concerns RNA interactions involving double-loop interactions. These are used by the bicoid mRNA to form dimers, a property that appears to be important for mRNA localization in drosophila embryo, and by bacteriophage phi29 pRNA which forms hexamers that participate in the translocation of the DNA genome through the portal vertex of the capsid. Despite the high diversity of systems and mechanisms, some common features can be highlighted. (1) Efficient recognition requires rapid bi-molecular binding rates, regardless of the RNA pairing scheme. (2) The initial recognition is favored by particular conformations of the loops enabling a proper presentation of nucleotides (generally a restricted number) that initiate the recognition process. (3) The fate of the initial reversible loop-loop complex is dictated by both functional and structural constraints. RNA structures have evolved either to "freeze" the initial complex, or to convert it into a more stable one, which involves propagation of intermolecular interactions along topologically feasible pathways. Stabilization of the initial complex may also be assisted by proteins and/or formation of additional contacts.     

Sequence requirements for self-splicing of the Tetrahymena thermophila pre-ribosomal RNA.

The sequence requirements for splicing of the Tetrahymena pre-rRNA have been examined by altering the rRNA gene to produce versions that contain insertions and deletions within the intervening sequence (IVS). The altered genes were transcribed and the RNA tested for self-splicing in vitro. A number of insertions (8-54 nucleotides) at three locations had no effect on self-splicing activity. Two of these insertions, located at a site 5 nucleotides preceding the 3'-end of the IVS, did not alter the choice of the 3' splice site. Thus the 3' splice site is not chosen by its distance from a fixed point within the IVS. Analysis of deletions constructed at two sites revealed two structures, a hairpin loop and a stem-loop, that are entirely dispensable for IVS excision in vitro. Three other regions were found to be necessary. The regions that are important for self-splicing are not restricted to the conserved sequence elements that define this class of intervening sequences. The requirement for structures within the IVS for pre-rRNA splicing is in sharp contrast to the very limited role of IVS structure in nuclear pre-mRNA splicing.     

Importance of long-range interactions in protein folding

Long-range interactions play an active role in the stability of protein molecules. In this work, we have analyzed the importance of long-range interactions in different structural classes of globular proteins in terms of residue distances. We found that 85% of residues are involved in long-range contacts. The residues occurring in the range of 4-10 residues apart contribute more towards long-range contacts in all-alpha proteins while the range is 11-20 in all-beta proteins. The hydrophobic residues Cys, Ile and Val prefer the 11-20 range and all other residues prefer the 4-10 range. The residues in all-beta proteins have an average of 3-8 long-range contacts whereas the residues in other classes have 1-4 long-range contracts. Furthermore, the preference of residue pairs to the folding and stability will be discussed.     

Effect of short- and long-range interactions on protein folding

The relative importance of short- and long-range interactions is examined using a Monte Carlo simulation of protein folding on bovine pancreatic trypsin inhibitor. The model of the protein and the interaction energies were parametrized using X-ray structures of 30 native proteins. A nearest neighbor Ising model is used to determine the conformational state at each stage of the Monte Carlo procedure. Long-range interactions are simulated by contact free energies which become effective as two residues, separated by four or more residues along the chain, approach each other, and by disulfide-bond energies. Short-range interactions for residues separated by one, two, or three residues along the chain are also modeled by contact free energies and by -helical hydrogen bonds. A hard-sphere model is used to represent repulsive interactions. The ratios of short- to long-range interactions studied are 1:1, 2:1, 1:2, 0:1, and 1:0; e.g., for the 2:1 ratio, short-range interactions are weighted twice as much as long-range interactions, and for the 1:0 ratio, long-range interactions are omitted. For each ratio of short- to long-range interactions, a native conformation is found by a Monte Carlo procedure, a segment of 11 residues (residue numbers 1–11) is then rotated away from the rest of the molecule [breaking the 5–55 native disulfide bond, and moving this segment so that the distance between the sulfur atoms of the 5 and 55 cystine side chains (averaged for all native conformations) increases from 3.9 to 7.3 Å], and the Monte Carlo simulation is carried out (allowing the conformation of the whole molecule to change) until equilibrium is attained. For each ratio, the refolded conformation is compared to the native one using triangular distance maps and differential geometry distance criteria. With ratios of short- to long-range interaction energies of 1:1 and 0:1, the native disulfide bond could be re-formed; with ratios of 2:1 and 1:2 it did not; and with the 1:0 ratio, even a stable native conformation was not achieved. Therefore, long-range interactions (in addition to short-range ones) are required to bring remote parts of the protein together and to stabilize its native conformation.NIH Postdoctoral Fellow, 1977–1978.     

Bimodal loop-loop interactions increase the affinity of RNA aptamers for HIV-1 RNA structures

Multiple loop-loop interactions between adjacent RNA hairpins regulate gene expression in different organisms. To demonstrate that such natural interactions could be mimicked for generating RNA ligands that are able to recognize simultaneously at least two structured RNA targets, a double kissing complex model was designed. The target consisted of two HIV-1 transactivating responsive (TAR) RNA variants, BRU and MAL, connected by a non-nucleotidic linker. The double ligand was generated by combining the corresponding hairpin aptamers, R06BRU and R06MAL, identified previously by in vitro selection [Ducongé, F., and Toulmé, J. J (1999) RNA 5, 1605-1614]. The resulting interaction was analyzed by thermal denaturation monitored by UV spectroscopy, electrophoretic mobility shift assays (EMSAs), and surface plasmon resonance (SPR) experiments. The bimodal complex was characterized by a binding equilibrium constant increased by at least 1 order of magnitude compared to that of the complexes between the individual parent hairpins. This resulted from a slower dissociation rate. We then made use of such a strategy for targeting two structured functional motifs of the folded 5' untranslated region (5'UTR) of HIV-1. Two bivalent RNA ligands were designed that targeted simultaneously the TAR and dimerization initiation site (DIS) hairpins or the TAR and poly(A) ones. The results show that these ligands also displayed enhanced affinity for their target compared to the individual molecules. The work reported here suggests that bimodal structured RNA ligands might provide a way of increasing the affinity of aptamers for folded RNA targets.     

The self-splicing RNA of Tetrahymena is trapped in a less active conformation by gel purification

When the circular form of the self-splicing intervening sequence of Tetrahymena thermophila was purified by denaturing polyacrylamide gel electrophoresis by standard methods, the rate of its reaction with tetrauridylate decreased 150-fold at 30 degrees C and at least 1000-fold at 0 degrees C. The activity of the self-splicing RNA was restored by heating it to high temperature and letting it renature in the presence of Mg2+. The rate of reaction of tetrauridylate with the self-splicing RNA flanked by exons was also greatly decreased by gel purification. The difference in activation energies for the reaction of native and denatured intervening sequences suggests that a substantial conformational rearrangement of the gel-purified RNA occurs prior to reaction.     

Anchoring hairpin ribozymes to long target RNAs by loop-loop RNA interactions

Efficient ribozyme-mediated gene silencing requires the effective binding of a ribozyme to its specific target sequence. Stable stem-loop domains are key elements for efficiency of natural antisense RNAs. This work tests the possibility of using such naturally existing structural motifs for anchoring hairpin ribozymes when targeting long RNAs. Assays were performed with four catalytic antisense RNAs, based on the hairpin ribozyme (HP), that carried a stable stem-loop motif at their 3' end. Extensions consisted of one of the following motifs: the stem-loop II of the natural antisense RNA-CopA, its natural target in CopT, the TAR-RNA motif, or its complementary sequence alphaTAR. Interestingly, the presence of any of these antisense motifs resulted in an enhancement of catalytic performance against the ribozyme's 14-nucleotide-long target RNA (Swt). A series of artificial, long RNA substrates containing the Swt sequence and the natural TAR-RNA stem-loop were constructed and challenged with a catalytic antisense RNA carrying the TAR-complementary stem-loop. This cleaves each of these substrates significantly more efficiently than HP. The deletion of the TAR domain in the substrate, or its substitution by its complementary counterpart alphaTAR, abolishes the positive effect. These results suggest that the enhancement is owed to the interaction of both complementary stem-loop domains. Moreover, they demonstrate that the TAR domain can be used as an anchoring site to facilitate the access of hairpin ribozymes to their specific target sequences within TAR-containing RNAs.     

Origin of the phosphate at the ligation junction produced by self-splicing of Tetrahymena thermophila pre-ribosomal RNA

The exons of the self-splicing pre-ribosomal RNA of Tetrahymena thermophila are joined accurately in vitro, even when only 33 nucleotides of the natural 5' exon and 38 nucleotides of the natural 3' exon remain. RNA fingerprint analysis was used to identify the unique ribonuclease T1 oligonucleotide generated by exon ligation. Secondary digests of the ligation junction oligonucleotide with ribonuclease A confirmed the identity of the fragment and demonstrated that the phosphate group that forms the phosphodiester bond at the ligation junction is derived from the 5' position of a uridine nucleotide in the RNA. This observation supports the prediction that the splice junction phosphate is derived from the 3' splice site. These results emphasize the mechanistic similarities of RNA splicing reactions of the group I introns, group II introns and nuclear pre-mRNA introns.     

Two major tertiary folding transitions of the Tetrahymena catalytic RNA.

The L-21 Tetrahymena ribozyme, an RNA molecule with sequence-specific endoribonuclease activity derived from a self-splicing group I intron, provides a model system for studying the RNA folding problem. A 160 nucleotide, independently folding domain of tertiary structure (the P4-P6 domain) comprises about half of the ribozyme. We now apply Fe(II)-EDTA cleavage to mutants of the ribozyme to explore the role of individual structural elements in tertiary folding of the RNA at equilibrium. Deletion of peripheral elements near the 3' end of the ribozyme destabilizes a region of the catalytic core (P3-P7) without altering the folding of the P4-P6 domain. Three different mutations within the P4-P6 domain that destabilize its folding also shift the folding of the P3-P7 region of the catalytic core to higher MgCl2 concentrations. We conclude that the role of the extended P4-P6 domain and of the 3'-terminal peripheral elements is at least in part to stabilize the catalytic core. The organization of RNA into independently folding domains of tertiary structure may be common in large RNAs, including ribosomal RNAs. Furthermore, the observation of domain-domain interactions in a catalytic RNA supports the feasibility of a primitive spliceosome without any proteins.     

Transient RNA-protein interactions in RNA folding

The RNA folding trajectory features numerous off-pathway folding traps, which represent conformations that are often equally as stable as the native functional ones. Therefore, the conversion between these off-pathway structures and the native correctly folded ones is the critical step in RNA folding. This process, referred to as RNA refolding, is slow, and is represented by a transition state that has a characteristic high free energy. Because this kinetically limiting process occurs in vivo, proteins (called RNA chaperones) have evolved that facilitate the (re)folding of RNA molecules. Here, we present an overview of how proteins interact with RNA molecules in order to achieve properly folded states. In this respect, the discrimination between static and transient interactions is crucial, as different proteins have evolved a multitude of mechanisms for RNA remodeling. For RNA chaperones that act in a sequence-unspecific manner and without the use of external sources of energy, such as ATP, transient RNA-protein interactions represent the basis of the mode of action. By presenting stretches of positively charged amino acids that are positioned in defined spatial configurations, RNA chaperones enable the RNA backbone, via transient electrostatic interactions, to sample a wider conformational space that opens the route for efficient refolding reactions.     

The effect of long-range interactions on the secondary structure formation of proteins

The influence of long-range residue interactions on defining secondary structure in a protein has long been discussed and is often cited as the current limitation to accurate secondary structure prediction. There are several experimental examples where a local sequence alone is not sufficient to determine its secondary structure, but a comprehensive survey on a large data set has not yet been done. Interestingly, some earlier studies denied the negative effect of long-range interactions on secondary structure prediction accuracy. Here, we have introduced the residue contact order (RCO), which directly indicates the separation of contacting residues in terms of the position in the sequence, and examined the relationship between the RCO and the prediction accuracy. A large data set of 2777 nonhomologous proteins was used in our analysis. Unlike previous studies, we do find that prediction accuracy drops as residues have contacts with more distant residues. Moreover, this negative correlation between the RCO and the prediction accuracy was found not only for beta-strands, but also for alpha-helices. The prediction accuracy of beta-strands is lower if residues have a high RCO or a low RCO, which corresponds to the situation that a beta-sheet is formed by beta-strands from different chains in a protein complex. The reason why the current study draws the opposite conclusion from the previous studies is examined. The implication for protein folding is also discussed.     

New pathways in folding of the Tetrahymena group I RNA enzyme.

Previous studies have shown that the earliest detectable step in folding of the Tetrahymena ribozyme is tertiary structure formation of the peripheral element P5abc. This, along with other results, has suggested that P5abc may serve as a scaffold upon which additional tertiary structure is built. Herein we use the onset of oligonucleotide cleavage activity as a readout for native state formation and investigate the effect of P5abc on the rate of folding to the native structure. Despite the early folding of P5abc, its removal to give the E delta P5abc variant decreases the rate of attainment of an active structure less than fivefold (20-100 mM Mg2+, 15-50 degrees C). Furthermore, P5abc added in trans is able to bind the folded E delta P5abc ribozyme and promote oligonucleotide cleavage at least tenfold more rapidly than folding of the wild-type ribozyme, indicating that E delta P5abc does not have to first unfold before productively binding P5abc to form the true native state. This suggests that a state with the overall tertiary structure formed but with P5abc unfolded represents a viable on-pathway intermediate for the wild-type ribozyme. These results provide strong evidence for the existence of two pathways to the native state: in one pathway P5abc forms tertiary structure first, and in another it forms late. The pathway in which P5abc forms first is favored because P5abc can fold quickly and because its tertiary structure is stable in the absence of additional structured elements, not because P5abc formation is required for subsequent folding steps. In the course of these experiments, we also found that most of the ribozyme population does not reach the native state directly under standard conditions in vitro, but instead forms an inactive structure that is stable for hours. Finally, the fraction that does fold to the native state folds with a single rate constant of 1 min-1, suggesting that there are no significantly populated "fast-track" pathways that reach the native state directly by avoiding slow folding steps.     

Stabilities of HIV-1 DIS type RNA loop-loop interactions in vitro and in vivo

RNA loop–loop interactions are a prevalent motif in the formation of tertiary structure and are well suited to trigger molecular recognition between RNA molecules. We determined the stabilities of several loop–loop interactions with a constant 6 bp core sequence and varying unpaired flanking nucleotides and found that the flanking bases have a strong influence on the stability and ion dependence of the kissing complex. In general, the stabilities determined in 1 M Na⁺ are equivalent to those in the presence of near physiological Mg²⁺ concentrations. Therefore we further tested whether the stabilities determined in vitro and within yeast cells correlate, using a recently developed yeast RNA-hybrid system. For the majority of the loop types analyzed here, the melting temperatures determined in vitro are in good agreement with the relative β-galactosidase activity in yeast cells, showing that data derived from in vitro measurements reflect in vivo properties. The most stable interactions are the naturally occurring HIV-1 DIS MAL and LAI derived loops with the motif (5′ A^A/_GN₆A 3′), emphasizing the crucial role of stable kissing complexes in HIV genome dimerization.     

Role of conserved sequence elements 9L and 2 in self-splicing of the Tetrahymena ribosomal RNA precursor

Oligonucleotide-directed mutagenesis has been used to alter highly conserved sequences within the intervening sequence (IVS) of the Tetrahymena large ribosomal RNA precursor. Mutations within either sequence element 9L or element 2 eliminate splicing activity under standard in vitro splicing conditions. A double mutant with compensatory base changes in elements 9L and 2 has accurate splicing activity restored. Thus, the targeted nucleotides of elements 9L and 2 base-pair with one another in the IVS RNA, and pairing is important for self-splicing. Mutant splicing activities are restored by increased magnesium ion concentrations, supporting the conclusion that the role of the targeted bases in splicing is primarily structural. Based on the temperature dependence, we propose that a conformational switch involving pairing and unpairing of elements 9L and 2 is required for splicing.     

Interactions of cations with RNA loop-loop complexes

RNA loop-loop interactions are essential in many biological processes, including initiation of RNA folding into complex tertiary shapes, promotion of dimerization, and viral replication. In this article, we examine interactions of metal ions with five RNA loop-loop complexes of unique biological significance using explicit-solvent molecular-dynamics simulations. These simulations revealed the presence of solvent-accessible tunnels through the major groove of loop-loop interactions that attract and retain cations. Ion dynamics inside these loop-loop complexes were distinctly different from the dynamics of the counterion cloud surrounding RNA and depend on the number of basepairs between loops, purine sequence symmetry, and presence of unpaired nucleotides. The cationic uptake by kissing loops depends on the number of basepairs between loops. It is interesting that loop-loop complexes with similar functionality showed similarities in cation dynamics despite differences in sequence and loop size.     

The effect of long-range interactions between adsorbed ligands on the DNA helix-coil transition

The influence of different types of long-range interaction of ligands adsorbed on DNA on the helix-coil transition was theoretically considered. The contact interaction was shown to differ significantly from the long-rang one. It was shown also that even weak dependence of a long-range potential on a degree of helicity resulted in the strong changes of a DNA melting curve. This result allowed to understand the different experimental data on DNA melting in the presence of different substances which reduced AT-and GC-base pairs thermostability difference.     

Comparison between long-range interactions and contact order in determining the folding rate of two-state proteins: application of long-range order to folding rate prediction.

The contact order is believed to be an important factor for understanding protein folding mechanisms. In our earlier work, we have shown that the long-range interactions play a vital role in protein folding. In this work, we analyzed the contribution of long-range contacts to determine the folding rate of two-state proteins. We found that the residues that are close in space and are separated by at least ten to 15 residues in sequence are important determinants of folding rates, suggesting the presence of a folding nucleus at an interval of approximately 25 residues. A novel parameter "long-range order" has been proposed to predict protein folding rates. This parameter shows as good a relationship with the folding rate of two-state proteins as contact order. Further, we examined the minimum limit of residue separation to determine the long-range contacts for different structural classes. We observed an excellent correlation between long-range order and folding rate for all classes of globular proteins. We suggest that in mixed-class proteins, a larger number of residues can serve as folding nuclei compared to all-alpha and all-beta proteins. A simple statistical method has been developed to predict the folding rates of two-state proteins using the long-range order that produces an agreement with experimental results that is better or comparable to other methods in the literature.     

Characterization of an authentic intermediate in the self-splicing process of ribosomal precursor RNA in macronuclei of Tetrahymena thermophila.

We have characterized a 1.5 kb RNA species in T. thermophila macronuclei previously found in vivo and including intron sequences linked to the 3' exon. This IVS-3' exon RNA could be detected in gels as a discrete molecule only after denaturation of nuclear RNA. After addition of 32P-GTP, as splicing cofactor in a nuclear in vitro system, the IVS-3' exon RNA was labeled at its 5' terminus, as was the by-product of splicing, the excised IVS RNA. The time course of labeling indicates that the IVS-3' exon RNA acts like a reaction intermediate and specifically a kinetic precursor to IVS RNA. Partial nuclease digestions showed that the IVS-3' exon RNA and the IVS RNA have the same 5' terminal sequence. In addition the IVS-3' exon RNA can release the 15-mer oligonucleotide cleaved off during circularization of IVS RNA under conditions of high temperature. Taken together, the structural, functional, and kinetic properties of the IVS-3' exon RNA strongly suggest that it represents a previously postulated in vivo intermediate in the splicing pathway.