Determination of thermodynamic parameters for HIV DIS type loop-loop kissing complexes

Engineering Light-Up Aptamers for the Detection of RNA Hairpins through Kissing Interaction

Aptasensors are biosensors that include aptamers for detecting a target of interest. We engineered signaling aptasensors for the detection of RNA hairpins from the previously described malachite green (MG) RNA aptamer. The top part of this imperfect hairpin aptamer was modified in such a way that it can engage loop-loop (so-called kissing) interactions with RNA hairpins displaying partly complementary apical loops. These newly derived oligonucleotides named malaswitches bind their cognate fluorogenic ligand (MG) exclusively when RNA-RNA kissing complexes are formed, whereas MG does not bind to malaswitches alone. Consequently, the formation of the ternary target RNA-malaswitch RNA-MG complex results in fluorescence emission, and malaswitches constitute sensors for detecting RNA hairpins. Malaswitches were designed that specifically detect precursors of microRNAs let7b and miR-206.

Structure folding of RNA kissing complexes in salt solutions: predicting 3D structure, stability, and folding pathway

RNA kissing complexes are essential for genomic RNA dimerization and regulation of gene expression, and their structures and stability are critical to their biological functions. In this work, we used our previously developed coarse-grained model with an implicit structure-based electrostatic potential to predict three-dimensional (3D) structures and stability of RNA kissing complexes in salt solutions. For extensive RNA kissing complexes, our model shows great reliability in predicting 3D structures from their sequences, and our additional predictions indicate that the model can capture the dependence of 3D structures of RNA kissing complexes on monovalent/divalent ion concentrations. Moreover, the comparisons with extensive experimental data show that the model can make reliable predictions on the stability for various RNA kissing complexes over wide ranges of monovalent/divalent ion concentrations. Notably, for RNA kissing complexes, our further analyses show the important contribution of coaxial stacking to the 3D structures and stronger stability than the corresponding kissing-interface duplexes at high salts. Furthermore, our comprehensive analyses for RNA kissing complexes reveal that the thermally folding pathway for a complex sequence is mainly determined by the relative stability of two possible folded states of kissing complex and extended duplex, which can be significantly modulated by its sequence.

Thermodynamic investigation of kissing-loop interactions

Kissing loop interactions (KLIs) are a common motif that is critical in retroviral dimerization, viroid replication, mRNA, and riboswitches. In addition, KLIs are currently used in a variety of biotechnology applications, such as in aptamer sensors, RNA scaffolds and to stabilize vaccines for therapeutics. Here we describe the thermodynamics of a basic intramolecular DNA capable of engaging in a KLI, consisting of two hairpins connected by a flexible linker. Each hairpin loop has a five-nucleotide complementary sequence theoretically capable of engaging in a KLI. On either side of each loop is two thymines which will not engage in kissing but are present to provide more flexibility and optimal KLI positioning. Our results suggest that the KLI occurs even at physiological salt levels, and that the KLI does not alter the thermodynamics and stability of the two stem structures. The KLI does not involve all five nucleotides, or at least each base-pair stack is not making full contact. Adding a second strand complementary to the bottom of the kissing complex removes flexibility and causes destabilization of the stems. The KLI of this less flexible complex is maintained but the T_M is reduced, indicating an entopic penalty to its formation.

Structural Insight into the Unbound State of the DNA Analogue of the PreQ1 Riboswitch: A Thermodynamic Approach

The preQ₁ riboswitch aptamer domain is very dynamic in its unbound state with the ability to form multiple structures: a hairpin, kissing hairpins, and pseudoknot-like structure. The aim of this study is to determine whether the DNA analogue (PreQ₁) is able to form structures similar to that of the reported RNA aptamer. Using a thermodynamic approach, we report on structural determination using differential scanning calorimetry under different salt conditions. Further analysis of the primary sequence allowed us to design modified molecules to determine what potential structures are forming in this single-stranded DNA analogue. We found, in a 16 mM Na⁺ solution, PreQ₁ has three transitions with T_M values of 14.8, 19.4, and 26.2 °C and a total ΔH of -44.7 kcal/mol. With the increase in salt concentration to 116 mM, there are T_M values of 22.3, 28.7, and 38.9 °C and a ΔH of -69.1 kcal/mol, while at 216 mM, the three transitions have T_M values of 24.4, 31.6, and 42.9 °C with a total ΔH of -71.5 kcal/mol. Therefore, the increase in enthalpy is due to the formation of additional base-pair stacks. The modified molecules, which would inhibit pseudoknot formation, kissing hairpins, and internal loop interactions, were fully characterized and compared to the native DNA analogue. The analysis of the enthalpy and differential binding of counterions allows us to conclude this single-stranded DNA analogue under physiological conditions is not forming a pseudoknot-like structure. Instead, two potential structures, Compact-Hairpin and Kissing-Complex, are more likely and could be in equilibrium.

Imaging HIV-1 RNA dimerization in cells by multicolor super-resolution and fluctuation microscopies

Dimerization is a unique and vital characteristic of retroviral genomes. It is commonly accepted that genomic RNA (gRNA) must be dimeric at the plasma membrane of the infected cells to be packaged during virus assembly. However, where, when and how HIV-1 gRNA find each other and dimerize in the cell are long-standing questions that cannot be answered using conventional approaches. Here, we combine two state-of-the-art, multicolor RNA labeling strategies with two single-molecule microscopy technologies to address these questions. We used 3D-super-resolution structured illumination microscopy to analyze and quantify the spatial gRNA association throughout the cell and monitored the dynamics of RNA-RNA complexes in living-cells by cross-correlation fluctuation analysis. These sensitive and complementary approaches, combined with trans-complementation experiments, reveal for the first time the presence of interacting gRNA in the cytosol, a challenging observation due to the low frequency of these events and their dilution among the bulk of other RNAs, and allow the determination of the subcellular orchestration of the HIV-1 dimerization process.

Computational Design of a Circular RNA with Prionlike Behavior

RNA molecules engineered to fold into predefined conformations have enabled the design of a multitude of functional RNA devices in the field of synthetic biology and nanotechnology. More complex designs require efficient computational methods, which need to consider not only equilibrium thermodynamics but also the kinetics of structure formation. Here we present a novel type of RNA design that mimics the behavior of prions, that is, sequences capable of interaction-triggered autocatalytic replication of conformations. Our design was computed with the ViennaRNA package and is based on circular RNA that embeds domains amenable to intermolecular kissing interactions.

DNA-DNA kissing complexes as a new tool for the assembly of DNA nanostructures

Kissing-loop annealing of nucleic acids occurs in nature in several viruses and in prokaryotic replication, among other circumstances. Nucleobases of two nucleic acid strands (loops) interact with each other, although the two strands cannot wrap around each other completely because of the adjacent double-stranded regions (stems). In this study, we exploited DNA kissing-loop interaction for nanotechnological application. We functionalized the vertices of DNA tetrahedrons with DNA stem-loop sequences. The complementary loop sequence design allowed the hybridization of different tetrahedrons via kissing-loop interaction, which might be further exploited for nanotechnology applications like cargo transport and logical elements. Importantly, we were able to manipulate the stability of those kissing-loop complexes based on the choice and concentration of cations, the temperature and the number of complementary loops per tetrahedron either at the same or at different vertices. Moreover, variations in loop sequences allowed the characterization of necessary sequences within the loop as well as additional stability control of the kissing complexes. Therefore, the properties of the presented nanostructures make them an important tool for DNA nanotechnology.

RNA folding: structure prediction, folding kinetics and ion electrostatics

Beyond the "traditional" functions such as gene storage, transport and protein synthesis, recent discoveries reveal that RNAs have important "new" biological functions including the RNA silence and gene regulation of riboswitch. Such functions of noncoding RNAs are strongly coupled to the RNA structures and proper structure change, which naturally leads to the RNA folding problem including structure prediction and folding kinetics. Due to the polyanionic nature of RNAs, RNA folding structure, stability and kinetics are strongly coupled to the ion condition of solution. The main focus of this chapter is to review the recent progress in the three major aspects in RNA folding problem: structure prediction, folding kinetics and ion electrostatics. This chapter will introduce both the recent experimental and theoretical progress, while emphasize the theoretical modelling on the three aspects in RNA folding.

Interstrand cross-link and bioconjugate formation in RNA from a modified nucleotide

RNA oligonucleotides containing a phenyl selenide derivative of 5-methyluridine were chemically synthesized by solid-phase synthesis. The phenyl selenide is rapidly converted to an electrophilic, allylic phenyl seleneate under mild oxidative conditions. The phenyl seleneate yields interstrand cross-links when part of a duplex and is useful for synthesizing oligonucleotide conjugates. Formation of the latter is illustrated by reaction of an oligonucleotide containing the phenyl selenide with amino acids in the presence of mild oxidant. The products formed are analogous to those observed in tRNA that are believed to be formed posttranslationally via a biosynthetic intermediate that is chemically homologous to the phenyl seleneate.

A remarkably stable kissing-loop interaction defines substrate recognition by the Neurospora Varkud Satellite ribozyme

Kissing loops are tertiary structure elements that often play key roles in functional RNAs. In the Neurospora VS ribozyme, a kissing-loop interaction between the stem-loop I (SLI) substrate and stem-loop V (SLV) of the catalytic domain is known to play an important role in substrate recognition. In addition, this I/V kissing-loop interaction is associated with a helix shift in SLI that activates the substrate for catalysis. To better understand the role of this kissing-loop interaction in substrate recognition and activation by the VS ribozyme, we performed a thermodynamic characterization by isothermal titration calorimetry using isolated SLI and SLV stem-loops. We demonstrate that preshifted SLI variants have higher affinity for SLV than shiftable SLI variants, with an energetic cost of 1.8-3 kcal/mol for the helix shift in SLI. The affinity of the preshifted SLI for SLV is remarkably high, the interaction being more stable by 7-8 kcal/mol than predicted for a comparable duplex containing three Watson-Crick base pairs. The structural basis of this remarkable stability is discussed in light of previous NMR studies. Comparative thermodynamic studies reveal that kissing-loop complexes containing 6-7 Watson-Crick base pairs are as stable as predicted from comparable RNA duplexes; however, those with 2-3 Watson-Crick base pairs are more stable than predicted. Interestingly, the stability of SLI/ribozyme complexes is similar to that of SLI/SLV complexes. Thus, the I/V kissing loop interaction represents the predominant energetic contribution to substrate recognition by the trans-cleaving VS ribozyme.

Predicting structure and stability for RNA complexes with intermolecular loop-loop base-pairing

RNA loop-loop interactions are essential for genomic RNA dimerization and regulation of gene expression. In this article, a statistical mechanics-based computational method that predicts the structures and thermodynamic stabilities of RNA complexes with loop-loop kissing interactions is described. The method accounts for the entropy changes for the formation of loop-loop interactions, which is a notable advancement that other computational models have neglected. Benchmark tests with several experimentally validated systems show that the inclusion of the entropy parameters can indeed improve predictions for RNA complexes. Furthermore, the method can predict not only the native structures of RNA/RNA complexes but also alternative metastable structures. For instance, the model predicts that the SL1 domain of HIV-1 RNA can form two different dimer structures with similar stabilities. The prediction is consistent with experimental observation. In addition, the model predicts two different binding sites for hTR dimerization: One binding site has been experimentally proposed, and the other structure, which has a higher stability, is structurally feasible and needs further experimental validation.

HIV-1 DIS stem loop forms an obligatory bent kissing intermediate in the dimerization pathway

The HIV-1 dimerization initiation sequence (DIS) is a conserved palindrome in the apical loop of a conserved hairpin motif in the 5′-untranslated region of its RNA genome. DIS hairpin plays an important role in genome dimerization by forming a ‘kissing complex’ between two complementary hairpins. Understanding the kinetics of this interaction is key to exploiting DIS as a possible human immunodeficiency virus (HIV) drug target. Here, we present a single-molecule Förster resonance energy transfer (smFRET) study of the dimerization reaction kinetics. Our data show the real-time formation and dissociation dynamics of individual kissing complexes, as well as the formation of the mature extended duplex complex that is ultimately required for virion packaging. Interestingly, the single-molecule trajectories reveal the presence of a previously unobserved bent intermediate required for extended duplex formation. The universally conserved A272 is essential for the formation of this intermediate, which is stabilized by Mg²⁺, but not by K⁺ cations. We propose a 3D model of a possible bent intermediate and a minimal dimerization pathway consisting of three steps with two obligatory intermediates (kissing complex and bent intermediate) and driven by Mg²⁺ ions.

Analysis of diffuse K+ and Mg2+ ion binding to a two-base-pair kissing complex by single-molecule mechanical unfolding

The folding and stability of RNA tertiary interactions depend critically on cationic conditions. It is usually difficult, however, to isolate such effects on tertiary interactions from those on the entire RNA. By manipulating conformations of single RNA molecules using optical tweezers, we distinguished individual steps of breaking and forming of a two-base-pair kissing interaction from those of secondary folding. The binding of metal ions to the small tertiary structure appeared to be saturable with an apparent Kd of 160 mM for K(+) and 1.5 mM for Mg(2+). The kissing formation was estimated to be associated with binding of ~2-3 diffuse K(+) or Mg(2+) ions. At their saturated binding, Mg(2+) provided ~3 kcal/mol more stabilizing energy to the structure than K(+). Furthermore, the cations change the unkissing forces significantly more than the kissing ones. For example, the presence of Mg(2+) ions increased the average unkissing force from 21 pN to 44 pN, surprisingly high for breaking merely two base pairs; in contrast, the mean kissing force was changed by only 4.5 pN. Interestingly, the differential salt effects on the transition forces were not caused by different changes in the height of the kinetic barriers but were instead attributed to how different molecular structures respond to the applied force. Our results showed the importance of diffuse cation binding to the stability of tertiary interaction and demonstrated the utility of mechanical unfolding in studying tertiary interactions.

Molecular determinants of HIV-1 NCp7 chaperone activity in maturation of the HIV-1 dimerization initiation site

Human immunodeficiency virus genome dimerization is initiated through an RNA-RNA kissing interaction formed via the dimerization initiation site (DIS) loop sequence, which has been proposed to be converted to a more thermodynamically stable linkage by the viral p7 form of the nucleocapsid protein (NC). Here, we systematically probed the role of specific amino acids of NCp7 in its chaperone activity in the DIS conversion using 2-aminopurine (2-AP) fluorescence and nuclear magnetic resonance spectroscopy. Through comparative analysis of NCp7 mutants, the presence of positively charged residues in the N-terminus was found to be essential for both helix destabilization and strand transfer functions. It was also observed that the presence and type of the Zn finger is important for NCp7 chaperone activity, but not the order of the Zn fingers. Swapping single aromatic residues between Zn fingers had a significant effect on NCp7 activity; however, these mutants did not exhibit the same activity as mutants in which the order of the Zn fingers was changed, indicating a functional role for other flanking residues. RNA chaperone activity is further correlated with NCp7 structure and interaction with RNA through comparative analysis of nuclear magnetic resonance spectra of NCp7 variants, and complexes of these proteins with the DIS dimer.

The role of salt concentration and magnesium binding in HIV-1 subtype-A and subtype-B kissing loop monomer structures

The subtype-B monomers of the human immunodeficiency virus type-1 (HIV-1) have experimentally been shown to dimerize at high salt concentration or in the presence of magnesium, while the dimerization of the subtype-A monomers requires magnesium binding at the G273 or G274 phosphate groups regardless of salt concentration. We used explicit solvent molecular dynamics (MD) simulations to investigate the conformational changes in subtype-A and -B monomers in different salt concentrations, and we found that our MD simulation results are consistent with those of experiments. At low salt concentration, hairpin loop structures of both subtypes were deformed and bases in the hairpin loop were turned inside. At high salt concentrations, the subtype-B monomer maintained the hairpin loop shape and most bases in the hairpin loop pointed out, while the subtype-A monomer showed a severe deformation. We also found that the flanking bases in the subtype-B stabilize the hairpin loop, while the flanking base G273 in the subtype-A caused a significant deformation. However, a bound magnesium ion at the G273 or G274 phosphate groups controlled the behavior of the G273 base and prevented the subtype-A monomer from deformation. We also applied restraints to both subtypes to examine the role of high salt concentration or magnesium binding. While restraints were applied, both subtypes at 0 M salt concentration maintained their shapes. However, when restraints were removed, they deformed significantly. Therefore, we suggest that the dimerization of both subtypes requires the proper conformation of the monomers which is induced by the appropriate salt strength and magnesium ion binding.

Mechanism of enhanced mechanical stability of a minimal RNA kissing complex elucidated by nonequilibrium molecular dynamics simulations

An RNA kissing loop from the Moloney murine leukemia virus (MMLV) exhibits unusual mechanical stability despite having only two intermolecular base pairs. Mutations at this junction have been shown to destabilize genome dimerization, with concomitant reductions in viral packaging efficiency and infectivity. Optical tweezers experiments have shown that it requires as much force to break the MMLV kissing-loop complex as is required to unfold an 11-bp RNA hairpin [Li PTX, Bustamante C, Tinoco I (2006) Proc Natl Acad Sci USA 103:15847-15852]. Using nonequilibrium all-atom molecular dynamics simulations, we have developed a detailed model for the kinetic intermediates of the force-induced dissociation of the MMLV dimerization initiation site kissing loop. Two hundred and eight dissociation events were simulated (approximately 16 μs total simulation time) under conditions of constant applied external force, which we use to construct a Markov state model for kissing-loop dissociation. We find that the complex undergoes a conformational rearrangement, which allows for equal distribution of the applied force among all of the intermolecular hydrogen bonds, which is intrinsically more stable than the sequential unzipping of an ordinary hairpin. Stacking interactions with adjacent, unpaired loop adenines further stabilize the complex by increasing the repair rate of partially broken H-bonds. These stacking interactions are prominently featured in the transition state, which requires additional coordinates orthogonal to the end-to-end extension to be uniquely identified. We propose that these stabilizing features explain the unusual stability of other retroviral kissing-loop complexes such as the HIV dimerization site.

Dominant role of the 5' TAR bulge in dimerization of HIV-1 genomic RNA, but no evidence of TAR-TAR kissing during in vivo virus assembly

The 5' untranslated region of HIV-1 genomic RNA (gRNA) contains two stem-loop structures that appear to be equally important for gRNA dimerization: the 57-nucleotide 5' TAR, at the very 5' end, and the 35-nucleotide SL1 (nucleotides 243-277). SL1 is well-known for containing the dimerization initiation site (DIS) in its apical loop. The DIS is a six-nucleotide palindrome. Here, we investigated the mechanism of TAR-directed gRNA dimerization. We found that the trinucleotide bulge (UCU24) of the 5' TAR has dominant impacts on both formation of HIV-1 RNA dimers and maturation of the formed dimers. The ΔUCU trinucleotide deletion strongly inhibited the first process and blocked the other, thus impairing gRNA dimerization as severely as deletion of the entire 5' TAR, and more severely than deletion of the DIS, inactivation of the viral protease, or most severe mutations in the nucleocapsid protein. The apical loop of TAR contains a 10-nucleotide palindrome that has been postulated to stimulate gRNA dimerization by a TAR-TAR kissing mechanism analogous to the one used by SL1 to stimulate dimerization. Using mutations that strongly destabilize formation of the TAR palindrome duplex, as well as compensatory mutations that restore duplex formation to a wild-type-like level, we found no evidence of TAR-TAR kissing, even though mutations nullifying the kissing potential of the TAR palindrome could impair dimerization by a mechanism other than hindering of SL1. However, nullifying the kissing potential of TAR had much less severe effects than ΔUCU. By not uncovering a dimerization mechanism intrinsic to TAR, our data suggest that TAR mutations exert their effect 3' of TAR, yet not on SL1, because TAR and SL1 mutations have synergistic effects on gRNA dimerization.

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.

Structure and stability of RNA/RNA kissing complex: with application to HIV dimerization initiation signal

We develop a statistical mechanical model to predict the structure and folding stability of the RNA/RNA kissing-loop complex. One of the key ingredients of the theory is the conformational entropy for the RNA/RNA kissing complex. We employ the recently developed virtual bond-based RNA folding model (Vfold model) to evaluate the entropy parameters for the different types of kissing loops. A benchmark test against experiments suggests that the entropy calculation is reliable. As an application of the model, we apply the model to investigate the structure and folding thermodynamics for the kissing complex of the HIV-1 dimerization initiation signal. With the physics-based energetic parameters, we compute the free energy landscape for the HIV-1 dimer. From the energy landscape, we identify two minimal free energy structures, which correspond to the kissing-loop dimer and the extended-duplex dimer, respectively. The results support the two-step dimerization process for the HIV-1 replication cycle. Furthermore, based on the Vfold model and energy minimization, the theory can predict the native structure as well as the local minima in the free energy landscape. The root-mean-square deviations (RMSDs) for the predicted kissing-loop dimer and extended-duplex dimer are ~3.0 Å. The method developed here provides a new method to study the RNA/RNA kissing complex.

Salt contribution to RNA tertiary structure folding stability

Accurate quantification of the ionic contribution to RNA folding stability could greatly enhance our ability to understand and predict RNA functions. Recently, motivated by the potential importance of ion correlation and fluctuation in RNA folding, we developed the tightly bound ion (TBI) model. Extensive experimental tests showed that the TBI model can lead to better treatment of multivalent ions than the Poisson-Boltzmann equation. In this study, we use the model to quantify the contribution of salt (Na(+) and Mg(2+)) to the RNA tertiary structure folding free energy. Folding of the RNA tertiary structure often involves intermediates. We focus on the folding transition from an intermediate state to the native state, and compute the electrostatic folding free energy of the RNA. Based on systematic calculations for a variety of RNA molecules, we derive a set of formulas for the electrostatic free energy for tertiary structural folding as a function of the sequence length and compactness of the RNA and the Na(+) and Mg(2+) concentrations. Extensive comparisons with experimental data suggest that our model and the extracted empirical formulas are quite reliable.

Melting of model HIV-1 stem-loop 1 RNA dimers monitored by 2-aminopurine fluorescence

Viral maturation of HIV-1 involves refolding of its genomic RNA, which is believed to include a rearrangement of the SL1 stem-loop from a metastable conformation called kissing loop dimer (KD) to a stable one termed extended dimer (ED). To investigate this rearrangement in vitro we have studied the thermal melting of the RNA dimers formed by slightly modified 23-nucleotide SL1 RNA of HIV-1 Mal. Local structural changes in the RNA dimers during the melting were monitored by changes in the fluorescence of 2-aminopurine (2AP) incorporated in predetermined positions of RNA. We have shown that the stem regions of both preformed KD and ED melt in the temperature interval from 75 ° C to 90 ° C. Kissing loop interface of the KD RNA is found to be disrupted at lower temperatures from 20 ° C to 55 ° C, at which the stem regions remain intact. Conversion of the preformed KD to ED overcoming the kinetic barrier occurs between 55 ° C and 65 ° C. The melting of "loop-loop" regions in both preformed and newly formed EDs takes place around 70 ° C. Our finding that thermoinduced KD-to-ED conversion is preceded by transient dissociation of loop-loop interface disagrees with a common idea of strand exchange without disruption of loop-loop-contact.

A long-distance RNA-RNA interaction plays an important role in programmed -1 ribosomal frameshifting in the translation of p88 replicase protein of Red clover necrotic mosaic virus

Programmed -1 ribosomal frameshifting (-1 PRF) is one viral translation strategy to express overlapping genes in positive-strand RNA viruses. Red clover necrotic mosaic virus (RCNMV) uses this strategy to express its replicase component protein p88. In this study, we used a cell-free translation system to map cis-acting RNA elements required for -1 PRF. Our results show that a small stem-loop structure adjacent to the cap-independent translation element in the 3' untranslated region (UTR) of RCNMV RNA1 is required for -1 PRF. Site-directed mutagenesis experiments suggested that this stem-loop regulates -1 PRF via base-pairing with complementary sequences in a bulged stem-loop adjacent to the shifty site. The existence of RNA elements responsible for -1 PRF and the cap-independent translation of replicase proteins in the 3' UTR of RNA1 might be important for switching translation to replication and for regulating the ratio of p88 to p27.

Importance of diffuse metal ion binding to RNA

RNAs are highly charged polyanionic molecules. RNA structure and function are strongly correlated with the ionic condition of the solution. The primary focus of this article is on the role of diffusive ions in RNA folding. Due to the long-range nature of electrostatic interactions, the diffuse ions can contribute significantly to RNA structural stability and folding kinetics. We present an overview of the experimental findings as well as the theoretical developments on the diffuse ion effects in RNA folding. This review places heavy emphasis on the effect of magnesium ions. Magnesium ions play a highly efficient role in stabilizing RNA tertiary structures and promoting tertiary structural folding. The highly efficient role goes beyond the mean-field effect such as the ionic strength. In addition to the effects of specific ion binding and ion dehydration, ion-ion correlation for the diffuse ions can contribute to the efficient role of the multivalent ions such as the magnesium ions in RNA folding.

Predicting ion binding properties for RNA tertiary structures

Recent experiments pointed to the potential importance of ion correlation for multivalent ions such as Mg(2+) ions in RNA folding. In this study, we develop an all-atom model to predict the ion electrostatics in RNA folding. The model can treat ion correlation effects explicitly by considering an ensemble of discrete ion distributions. In contrast to the previous coarse-grained models that can treat ion correlation, this new model is based on all-atom nucleic acid structures. Thus, unlike the previous coarse-grained models, this new model allows us to treat complex tertiary structures such as HIV-1 DIS type RNA kissing complexes. Theory-experiment comparisons for a variety of tertiary structures indicate that the model gives improved predictions over the Poisson-Boltzmann theory, which underestimates the Mg(2+) binding in the competition with Na(+). Further systematic theory-experiment comparisons for a series of tertiary structures lead to a set of analytical formulas for Mg(2+)/Na(+) ion-binding to various RNA and DNA structures over a wide range of Mg(2+) and Na(+) concentrations.

Compact intermediates in RNA folding

Large noncoding RNAs fold into their biologically functional structures via compact yet disordered intermediates, which couple the stable secondary structure of the RNA with the emerging tertiary fold. The specificity of the collapse transition, which coincides with the assembly of helical domains, depends on RNA sequence and counterions. It determines the specificity of the folding pathways and the magnitude of the free energy barriers to the ensuing search for the native conformation. By coupling helix assembly with nascent tertiary interactions, compact folding intermediates in RNA also play a crucial role in ligand binding and RNA-protein recognition.

Identification of antisense RNA stem-loops that inhibit RNA-protein interactions using a bacterial reporter system

Many well-characterized examples of antisense RNAs from prokaryotic systems involve hybridization of the looped regions of stem–loop RNAs, presumably due to the high thermodynamic stability of the resulting loop–loop and loop–linear interactions. In this study, the identification of RNA stem–loops that inhibit U1A protein binding to the hpII RNA through RNA–RNA interactions was attempted using a bacterial reporter system based on phage λ N-mediated antitermination. As a result, loop sequences possessing 7–8 base complementarity to the 5′ region of the boxA element important for functional antitermination complex formation, but not the U1 hpII loop, were identified. In vitro and in vivo mutational analysis strongly suggested that the selected loop sequences were binding to the boxA region, and that the structure of the antisense stem–loop was important for optimal inhibitory activity. Next, in an attempt to demonstrate the ability to inhibit the interaction between the U1A protein and the hpII RNA, the rational design of an RNA stem–loop that inhibits U1A-binding to a modified hpII was carried out. Moderate inhibitory activity was observed, showing that it is possible to design and select antisense RNA stem–loops that disrupt various types of RNA–protein interactions.

Mechanical unfolding of two DIS RNA kissing complexes from HIV-1

An RNA kissing complex formed by the dimerization initiation site plays a critical role in the survival and infectivity of human immunodeficiency virus. Two dimerization initiation site kissing sequences, Mal and Lai, have been found in most human immunodeficiency virus 1 variants. Formation and stability of these RNA kissing complexes depend crucially on cationic conditions, particularly Mg 2+. Using optical tweezers, we investigated the mechanical unfolding of single RNA molecules with either Mal-type (GUGCAC) or Lai-type (GCGCGC) kissing complexes under various ionic conditions. The force required to disrupt the kissing interaction of the two structures, the rip force, is sensitive to concentrations of KCl and MgCl2; addition of 3 mM MgCl2 to 100 mM KCl changes the rip force of Mal from 21 +/- 4 to 46 +/- 3 pN. From the rip force distribution, the kinetics of breaking the kissing interaction is calculated as a function of force and cation concentration. The two kissing complexes have distinct unfolding transition states, as shown by different values of deltaX(++), which is the distance from the folded structure to the unfolding transition state. The deltaX(++) of Mal is approximately 0.6 nm smaller than that of Lai, suggesting that fewer kissing base pairs are broken at the transition state of the former, consistent with observations that the Lai-type kissing complex is more stable and requires significantly more force to unfold than the Mal type. More importantly, neither K+ nor Mg 2+ significantly changes the position of the transition state along the reaction coordinate. However, increasing concentrations of cations increase the kinetic barrier. We derived a cation-specific parameter, m, to describe how the height of the kinetic barrier depends on the concentration of cations. Our results suggest that Mg 2+ greatly slows down the unfolding of the kissing complex but has moderate effects on the formation kinetics of the structure.

Dissecting structural transitions in the HIV-1 dimerization initiation site RNA using 2-aminopurine fluorescence

A highly conserved 35 nucleotide RNA stem-loop, the dimerization initiation site (DIS), in the 5' untranslated region (UTR) of the human immunodeficiency virus type I (HIV-1) genome has been identified as the sequence primarily responsible for initiation of viral genome dimerization. The DIS initiates viral genome dimerization through a loop-loop 'kissing' interaction and is converted from an intermediate 'kissing' to a more thermodynamically stable extended duplex dimer in a conformational rearrangement that is chaperoned by the HIV-1 nucleocapsid protein (NCp7). Here we describe fluorescence methods designed to probe local RNA dynamics and structural transitions associated with the DIS dimer formation and its NCp7 chaperoned structural conversion. These methods take advantage of the exquisite sensitivity of the quantum yield of the fluorescent nucleotide base analog, 2-aminopurine (2-AP), to its immediate structural and dynamic environment. The 2-AP fluorescence methods described allow a detailed kinetic and thermodynamic examination of this type of RNA-RNA interaction, as well as an analysis of the molecular mechanism of NCp7 chaperone activity.

Effect of Mg(2+) and Na(+) on the nucleic acid chaperone activity of HIV-1 nucleocapsid protein: implications for reverse transcription

The human immunodeficiency virus type 1 (HIV-1) nucleocapsid protein (NC) is an essential protein for retroviral replication. Among its numerous functions, NC is a nucleic acid (NA) chaperone protein that catalyzes NA rearrangements leading to the formation of thermodynamically more stable conformations. In vitro, NC chaperone activity is typically assayed under conditions of low or no Mg(2+), even though reverse transcription requires the presence of divalent cations. Here, the chaperone activity of HIV-1 NC was studied as a function of varying Na(+) and Mg(2+) concentrations by investigating the annealing of complementary DNA and RNA hairpins derived from the trans-activation response domain of the HIV genome. This reaction mimics the annealing step of the minus-strand transfer process in reverse transcription. Gel-shift annealing and sedimentation assays were used to monitor the annealing kinetics and aggregation activity of NC, respectively. In the absence of protein, a limited ability of Na(+) and Mg(2+) cations to facilitate hairpin annealing was observed, whereas NC stimulated the annealing 10(3)- to 10(5)-fold. The major effect of either NC or the cations is on the rate of bimolecular association of the hairpins. This effect is especially strong under conditions wherein NC induces NA aggregation. Titration with NC and NC/Mg(2+) competition studies showed that the annealing kinetics depends only on the level of NA saturation with NC. NC competes with Mg(2+) or Na(+) for sequence-nonspecific NA binding similar to a simple trivalent cation. Upon saturation, NC induces attraction between NA molecules corresponding to approximately 0.3 kcal/mol/nucleotide, in agreement with an electrostatic mechanism of NC-induced NA aggregation. These data provide insights into the variable effects of NC's chaperone activity observed during in vitro studies of divalent metal-dependent reverse transcription reactions and suggest the feasibility of NC-facilitated proviral DNA synthesis within the mature capsid core.

Predicting electrostatic forces in RNA folding

Metal ion-mediated electrostatic interactions are critical to RNA folding. Although considerable progress has been made in mechanistic studies, the problem of accurate predictions for the ion effects in RNA folding remains unsolved, mainly due to the complexity of several potentially important issues such as ion correlation and dehydration effects. In this chapter, after giving a brief overview of the experimental findings and theoretical approaches, we focus on a recently developed new model, the tightly bound ion (TBI) model, for ion electrostatics in RNA folding. The model is unique because it can treat ion correlation and fluctuation effects for realistic RNA 3D structures. For monovalent ion (such as Na(+)) solutions, where ion correlation is weak, TBI and the Poisson-Boltzmann (PB) theory give the same results and the results agree with the experimental data. For multivalent ion (such as Mg(2+)) solutions, where ion correlation can be strong, however, TBI gives much improved predictions than the PB. Moreover, the model suggests an ion correlation-induced mechanism for the unusual efficiency of Mg(2+) ions in the stabilization of RNA tertiary folds. In this chapter, after introducing the theoretical framework of the TBI model, we will describe how to apply the model to predict ion-binding properties and ion-dependent folding stabilities.

How RNA unfolds and refolds

Understanding how RNA folds and what causes it to unfold has become more important as knowledge of the diverse functions of RNA has increased. Here we review the contributions of single-molecule experiments to providing answers to questions such as: How much energy is required to unfold a secondary or tertiary structure? How fast is the process? How do helicases unwind double helices? Are the unwinding activities of RNA-dependent RNA polymerases and of ribosomes different from other helicases? We discuss the use of optical tweezers to monitor the unfolding activities of helicases, polymerases, and ribosomes, and to apply force to unfold RNAs directly. We also review the applications of fluorescence and fluorescence resonance energy transfer to measure RNA dynamics.

Salt dependence of nucleic acid hairpin stability

Single-stranded junctions/loops are frequently occurring structural motifs in nucleic acid structures. Due to the polyanionic nature of the nucleic acid backbone, metal ions play a crucial role in the loop stability. Here we use the tightly bound ion theory, which can account for the possible ion correlation and ensemble (fluctuation) effects, to predict the ion-dependence of loop and stem-loop (hairpin) free energies. The predicted loop free energy is a function of the loop length, the loop end-to-end distance, and the ion (Na(+) and Mg(2+) in this study) concentrations. Based on the statistical mechanical calculations, we derive a set of empirical formulas for the loop thermodynamic parameters as functions of Na(+) and Mg(2+) concentrations. For three specific types of loops, namely, hairpin, bulge, and internal loops, the predicted free energies agree with the experimental data. Further applications of these empirical formulas to RNA and DNA hairpin stability lead to good agreements with the available experimental data. Our results indicate that the ion-dependent loop stability makes significant contribution to the overall ion-dependence of the hairpin stability.

Nucleocapsid protein-mediated maturation of dimer initiation complex of full-length SL1 stemloop of HIV-1: sequence effects and mechanism of RNA refolding

Specific binding of HIV-1 viral protein NCp7 to a unique 35-base RNA stem-loop SL1 is critical for formation and packaging of the genomic RNA dimer found within HIV-1 virions. NCp7 binding stimulates refolding of SL1 from a metastable kissing dimer (KD) into thermodynamically stable linear dimer (LD). Using UV melting, gel electrophoresis and heteronuclear NMR, we investigated effects of various site-specific mutations within the full-length SL1 on temperature- or NCp7-induced refolding in vitro. Refolding involved intramolecular melting of SL1 stems but not dissociation of the intermolecular KD interface. Refolding required only two NCp7 molecules per KD but was limited by the amount of NCp7 present, implying that the protein does not catalytically promote refolding. Efficient refolding depended strictly on the presence and, to a lesser degree, on sequence of a highly conserved G-rich internal loop that normally limits thermal stability of the SL1 stem. Adding two base pairs to the lower stem created a hyperstable SL1 mutant that failed to refold, even when bound by NCp7 at high stoichiometries. NMR analysis of these kinetically trapped mutant RNA–protein complexes indicated that NCp7 initiates refolding by dissociating base pairs in the upper stem of SL1. This study illuminates structural transitions critical for HIV-1 assembly and replication.

RNA helix stability in mixed Na+/Mg2+ solution

A recently developed tightly bound ion model can account for the correlation and fluctuation (i.e., different binding modes) of bound ions. However, the model cannot treat mixed ion solutions, which are physiologically relevant and biologically significant, and the model was based on B-DNA helices and thus cannot directly treat RNA helices. In the present study, we investigate the effects of ion correlation and fluctuation on the thermodynamic stability of finite length RNA helices immersed in a mixed solution of monovalent and divalent ions. Experimental comparisons demonstrate that the model gives improved predictions over the Poisson-Boltzmann theory, which has been found to underestimate the roles of multivalent ions such as Mg2+ in stabilizing DNA and RNA helices. The tightly bound ion model makes quantitative predictions on how the Na+-Mg2+ competition determines helix stability and its helix length-dependence. In addition, the model gives empirical formulas for the thermodynamic parameters as functions of Na+/Mg2+ concentrations and helix length. Such formulas can be quite useful for practical applications.

Fluorescent HIV-1 Dimerization Initiation Site: design, properties, and use for ligand discovery

The HIV-1 Dimerization Initiation Site (DIS) is an intriguing, yet underutilized, viral RNA target for potential antiretroviral therapy. To study the recognition features of this target and to provide a quantitative, rapid, and real-time tool for the discovery of new binders, a fluorescence-based assay has been constructed. It relies on strategic incorporation of 2-aminopurine, an isosteric fluorescent adenosine analogue, into short hairpin RNA constructs. These oligomers self-associate to form a kissing loop that thermally rearranges into a more stable extended duplex, thereby mimicking the association and structural features of the native RNA sequence. We demonstrate the ability of two fluorescent DIS constructs, DIS272(2AP) and DIS273(2AP), to report the binding of known DIS binders via changes in their emission intensity. Binding of aminoglycosides such as paromomycin to DIS272(2AP) results in significant fluorescence enhancement, while ligand binding to DIS273(2AP) results in fluorescence quenching. These observations are rationalized by comparison to the sequence-analogous bacterial A-site, where the relative emission of the fluorescent probe is dependent on the placement of the flexible purine residues inside or outside the helical domain. Analysis of binding isotherms generated using DIS272(2AP) yields submicromolar EC50 values for paromomycin (0.5 +/- 0.2 microM) and neomycin B (0.6 +/- 0.2 microM). Other neomycin-family aminoglycosides are less potent binders with neamine, the core pharmacophore, displaying the lowest affinity of 21 +/- 1 microM. Screening of additional aminoglycosides and their derivatives led to the discovery of new, previously unreported, aminoglycoside binders of the HIV DIS RNA, among them butirosin A (5.5 +/- 0.6 microM) and apramycin (7.6 +/- 1.0 microM). A conformationally constrained neomycin B analogue displays a rather high affinity to the DIS (1.9 +/- 0.2 microM). Among a series of nucleobase aminoglycoside conjugates, only the uracil derivatives display a measurable affinity using this assay with EC50 values in the 2 microM range. In addition, similarity between the solution behavior of HIV-1 DIS and the bacterial decoding A-site has been observed, particularly with respect to the intra- and extra-helical residence of the conformationally flexible A residues within the bulge. Taken together, the observations reported here shed light on the solution behavior of this important RNA target and are likely to facilitate the design of new DIS selective ligands as potential antiretroviral agents.

Influence of cationic molecules on the hairpin to duplex equilibria of self-complementary DNA and RNA oligonucleotides

A self-complementary nucleotide sequence can form both a unimolecular hairpin and a bimolecular duplex. In this study, the secondary structures of the self-complementary DNA and RNA oligonucleotides with different sequences and lengths were investigated under various solution conditions by gel electrophoresis, circular dichroism (CD) and electron paramagnetic resonance (EPR) spectroscopy and a ultraviolet (UV) melting analysis. The DNA sequences tended to adopt a hairpin conformation at low cation concentrations, but a bimolecular duplex was preferentially formed at an elevated cationic strength. On the other hand, fully matched RNA sequences adopted a bimolecular duplex regardless of the cation concentration. The thermal melting experiments indicated a greater change in the melting temperature of the bimolecular duplexes (by approximately 20 degrees C) than that of the hairpin (by approximately 10 degrees C) by increasing the NaCl concentration from 10 mM to 1 M. Hairpin formations were also observed for the palindrome DNA sequences derived from Escherichia coli, but association of the complementary palindrome sequences was observed when spermine, one of the major cationic molecules in a cell, existed at the physiological concentration. The results indicate the role of cations for shifting the structural equilibrium toward a nucleotide assembly and implicate nucleotide structures in cells.

DNA hairpins: fuel for autonomous DNA devices

We present a study of the hybridization of complementary DNA hairpin loops, with particular reference to their use as fuel for autonomous DNA devices. The rate of spontaneous hybridization between complementary hairpins can be reduced by increasing the neck length or decreasing the loop length. Hairpins with larger loops rapidly form long-lived kissed complexes. Hairpin loops may be opened by strand displacement using an opening strand that contains the same sequence as half of the neck and a "toehold" complementary to a single-stranded domain adjacent to the neck. We find loop opening via an external toehold to be 10-100 times faster than via an internal toehold. We measure rates of loop opening by opening strands that are at least 1000 times faster than the spontaneous interaction between hairpins. We discuss suitable choices for loop, neck, and toehold length for hairpin loops to be used as fuel for autonomous DNA devices.

Mechanistic studies of mini-TAR RNA/DNA annealing in the absence and presence of HIV-1 nucleocapsid protein

HIV-1 reverse transcription involves several nucleic acid rearrangements, which are catalyzed by the nucleocapsid protein (NC). Annealing of the trans-activation response element (TAR) DNA hairpin to a complementary TAR RNA hairpin, resulting in the formation of an extended 98-base-pair duplex, is an essential step in the minus-strand transfer step of reverse transcription. To elucidate the TAR RNA/DNA annealing reaction pathway, annealing kinetics were studied systematically by gel-shift assays performed in the presence or absence of HIV-1 NC. Truncated 27 nucleotide mini-TAR RNA and DNA constructs were used in this work. In the absence of NC, the annealing is slow, and involves the fast formation of an unstable extended "kissing" loop intermediate, followed by a slower strand exchange between the terminal stems. This annealing is very sensitive to loop-loop complementarity, as well as to nucleic acid concentration, ionic strength and temperature. NC stimulates the annealing approximately 5000-fold by stabilizing the bimolecular intermediate approximately 100 to 200-fold, and promoting the subsequent strand exchange reaction approximately 10 to 20-fold. NC concentration dependence studies suggest that there is a direct correlation between the amount of NC required to stabilize the intermediate and the amount needed to induce mini-TAR aggregation. Whereas saturating levels of NC are required to efficiently aggregate nucleic acids, sub-saturating NC is sufficient to significantly enhance duplex destabilization. Equilibrium levels of mini-TAR RNA/DNA annealing were also measured under a variety of conditions. Taken together, the results presented here provide a quantitative accounting of HIV-1 NC's aggregation and duplex destabilizing activity, and provide insights into the universal nucleic acid chaperone activity of this essential viral protein.

Partition function and base pairing probabilities of RNA heterodimers

BACKGROUND

RNA has been recognized as a key player in cellular regulation in recent years. In many cases, non-coding RNAs exert their function by binding to other nucleic acids, as in the case of microRNAs and snoRNAs. The specificity of these interactions derives from the stability of inter-molecular base pairing. The accurate computational treatment of RNA-RNA binding therefore lies at the heart of target prediction algorithms.

METHODS

The standard dynamic programming algorithms for computing secondary structures of linear single-stranded RNA molecules are extended to the co-folding of two interacting RNAs.

RESULTS

We present a program, RNAcofold, that computes the hybridization energy and base pairing pattern of a pair of interacting RNA molecules. In contrast to earlier approaches, complex internal structures in both RNAs are fully taken into account. RNAcofold supports the calculation of the minimum energy structure and of a complete set of suboptimal structures in an energy band above the ground state. Furthermore, it provides an extension of McCaskill's partition function algorithm to compute base pairing probabilities, realistic interaction energies, and equilibrium concentrations of duplex structures.

Solution RNA structures of the HIV-1 dimerization initiation site in the kissing-loop and extended-duplex dimers

Dimer formation of HIV-1 genomic RNA through its dimerization initiation site (DIS) is crucial to maintaining infectivity. Two types of dimers, the initially generated kissing-loop dimer and the subsequent product of the extended-duplex dimer, are formed in the stem-bulge-stem region with a loop including a self-complementary sequence. NMR chemical shift analysis of a 39mer RNA corresponding to DIS, DIS39, in the kissing-loop and extended-duplex dimers revealed that the three dimensional structures of the stem-bulge-stem region are extremely similar between the two types of dimers. Therefore, we designed two shorter RNA molecules, loop25 and bulge34, corresponding to the loop-stem region and the stem-bulge-stem region of DIS39, respectively. Based upon the chemical shift analysis, the conformation of the loop region of loop25 is identical to that of DIS39 for each of the two types of dimers. The conformation of bulge34 was also found to be the same as that of the corresponding region of DIS39. Thus, we determined the solution structures of loop25 in each of the two types of dimers as well as that of bulge34. Finally, the solution structures of DIS39 in the kissing-loop and extended-duplex dimers were determined by combining the parts of the structures. The solution structures of the two types of dimers were similar to each other in general with a difference found only in the A16 residue. The elucidation of the structures of DIS39 is important to understanding the molecular mechanism of the conformational dynamics of viral RNA molecules.

Thermodynamics of RNA-RNA binding

BACKGROUND

Reliable prediction of RNA-RNA binding energies is crucial, e.g. for the understanding on RNAi, microRNA-mRNA binding and antisense interactions. The thermodynamics of such RNA-RNA interactions can be understood as the sum of two energy contributions: (1) the energy necessary to 'open' the binding site and (2) the energy gained from hybridization.

METHODS

We present an extension of the standard partition function approach to RNA secondary structures that computes the probabilities Pu[i, j] that a sequence interval [i, j] is unpaired.

RESULTS

Comparison with experimental data shows that Pu[i, j] can be applied as a significant determinant of local target site accessibility for RNA interference (RNAi). Furthermore, these quantities can be used to rigorously determine binding free energies of short oligomers to large mRNA targets. The resource consumption is comparable with a single partition function computation for the large target molecule. We can show that RNAi efficiency correlates well with the binding energies of siRNAs to their respective mRNA target.

AVAILABILITY

RNAup will be distributed as part of the Vienna RNA Package, www.tbi.univie.ac.at/~ivo/RNA/

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.

A yeast RNA-hybrid system for the detection of RNA-RNA interactions in vivo

RNA-RNA interactions play a crucial role at many different levels of the cellular metabolism such as plasmid replication control, viral encapsidation, or transcriptional and translational regulation. Therefore, methods are necessary to investigate the molecular determinants of given interactions, including their stabilities, or to screen for new interacting partners. We designed an RNA-hybrid system in S. cerevisiae, based on the yeast three-hybrid system. In this setup, the activation of a reporter gene is dependent on the interaction of two RNAs. A loop-loop interaction similar to the dimerization initiation site of the HIV genome was used as a model system, demonstrating that in this novel RNA-hybrid system only cognate RNAs promote the activation of the reporter gene. Levels of reporter activation correlate well with interaction stabilities determined in vitro by UV melting analyses, suggesting that conditions used for the analysis of in vitro structural stabilities translate well into the intracellular environment. Furthermore, the system was applicable for a screen against a test library. Nine out of ten selected clones were identified as predicted interaction partners for the bait RNA. In summary, we present a yeast reporter system depending on RNA-RNA interactions, which can be used alternatively for analysis of known interactions or for screening libraries in search for new interaction partners.

Long-distance RNA-RNA interactions in plant virus gene expression and replication

The vast majority of plant and animal viruses have RNA genomes. Viral gene expression and replication are controlled by cis-acting elements in the viral genome, which have been viewed conventionally as localized structures. However, recent research has altered this perception and provided compelling evidence for cooperative activity involving distantly positioned RNA elements. This chapter focuses on viral RNA elements that interact across hundreds or thousands of intervening nucleotides to control translation, genomic RNA synthesis, and subgenomic mRNA transcription. We discuss evidence supporting the existence and function of the interactions, and speculate on the regulatory roles that such long-distance interactions play in the virus life cycle. We emphasize viruses in the Tombusviridae and Luteoviridae families in which long-distance interactions are best characterized, but similar phenomena in other viruses are also discussed. Many more examples likely remain undiscovered.