Intermolecular Base Pairing Research Articles

Bipartite viral RNA genome heterodimerization influences genome packaging and virion thermostability.

Multi-segmented viruses often multimerize their genomic segments to ensure efficient and stoichiometric packaging of the correct genetic cargo. In the bipartite Nodaviridae family, genome heterodimerization is also observed and conserved among different species. However, the nucleotide composition and biological function for this heterodimer remain unclear. Using Flock House virus as a model system, we developed a next-generation sequencing approach ("XL-ClickSeq") to probe heterodimer site sequences. We identified an intermolecular base-pairing site which contributed to heterodimerization in both wild-type and defective virus particles. Mutagenic disruption of this heterodimer site exhibited significant deficiencies in genome packaging and encapsidation specificity to viral genomic RNAs. Furthermore, the disruption of this intermolecular interaction directly impacts the thermostability of the mature virions. These results demonstrate that the intermolecular RNA-RNA interactions within the encapsidated genome of an RNA virus have an important role on virus particle integrity and thus may impact its transmission to a new host.IMPORTANCEFlock House virus is a member of Nodaviridae family of viruses, which provides a well-studied model virus for non-enveloped RNA virus assembly, cell entry, and replication. The Flock House virus genome consists of two separate RNA molecules, which can form a heterodimer upon heating of virus particles. Although similar RNA dimerization is utilized by other viruses (such as retroviruses) as a packaging mechanism and is conserved among Nodaviruses, the role of heterodimerization in the Nodavirus replication cycle is unclear. In this research, we identified the RNA sequences contributing to Flock House virus genome heterodimerization and discovered that such RNA-RNA interaction plays an essential role in virus packaging efficiency and particle integrity. This provides significant insight into how the interaction of packaged viral RNA may have a broader impact on the structural and functional properties of virus particles.

Intermolecular base-pairing interactions, a unique topology and exoribonuclease-resistant noncoding RNAs drive formation of viral chimeric RNAs in plants.

In plants, exoribonuclease-resistant RNAs (xrRNAs) are produced by many viruses. Whereas xrRNAs contribute to the pathogenicity of these viruses, the role of xrRNAs in the virus infectious cycle remains elusive. Here, we show that xrRNAs produced by a benyvirus (a multipartite RNA virus with four genomic segments) in plants are involved in the formation of monocistronic coat protein (CP)-encoding chimeric RNAs. Naturally occurring chimeric RNAs, we discovered, are composed of 5'-end of RNA 2 and 3'-end of either RNA 3 or RNA 4 bearing conservative exoribonuclease-resistant 'coremin' region. Using computational tools and site-directed mutagenesis, we show that de novo formation of chimeric RNAs requires intermolecular base-pairing interaction between 'coremin' and 3'-proximal part of the CP gene of RNA 2 as well as a stem-loop structure immediately adjacent to the CP gene. Moreover, knockdown of the expression of the XRN4 gene, encoding 5'→3' exoribonuclease, inhibits biogenesis of both xrRNAs and chimeric RNAs. Our findings suggest a novel mechanism involving a unique tropology of the intermolecular base-pairing complex between xrRNAs and RNA2 to promote formation of chimeric RNAs in plants. XrRNAs, essential for chimeric RNA biogenesis, are generated through the action of cytoplasmic Xrn 4 5'→3' exoribonuclease conserved in all plant species.

LinearCoFold and LinearCoPartition: linear-time algorithms for secondary structure prediction of interacting RNA molecules.

Many RNAs function through RNA-RNA interactions. Fast and reliable RNA structure prediction with consideration of RNA-RNA interaction is useful, however, existing tools are either too simplistic or too slow. To address this issue, we present LinearCoFold, which approximates the complete minimum free energy structure of two strands in linear time, and LinearCoPartition, which approximates the cofolding partition function and base pairing probabilities in linear time. LinearCoFold and LinearCoPartition are orders of magnitude faster than RNAcofold. For example, on a sequence pair with combined length of 26,190nt, LinearCoFold is 86.8× faster than RNAcofold MFE mode, and LinearCoPartition is 642.3× faster than RNAcofold partition function mode. Surprisingly, LinearCoFold and LinearCoPartition's predictions have higher PPV and sensitivity of intermolecular base pairs. Furthermore, we apply LinearCoFold to predict the RNA-RNA interaction between SARS-CoV-2 genomic RNA (gRNA) and human U4 small nuclear RNA (snRNA), which has been experimentally studied, and observe that LinearCoFold's prediction correlates better with the wet lab results than RNAcofold's.

RNA structure controls composition and network dynamics of condensates.

Biomolecular condensates can form by phase separation and likely serve myriad functions in cells. Many studies have characterized protein condensates, and sometimes these include model RNAs that are often nonspecific and/or homopolymeric, unlike native RNAs found in biomolecular condensates. How complex RNA sequences and structures control condensation with proteins remains poorly understood. We studied the influence of mRNA structure in phase separation using the Whi3 protein and a native mRNA binding partner, CLN3. We first created an evolutionary algorithm that designs shuffled CLN3 RNA sequences which maximize or minimize predicted free energies of folding while preserving mass, known Whi3 protein binding sites, encoded amino acid sequence, and nucleotide composition. RNAs with minimized free energies of folding contain many stable duplexes, while RNAs with maximized energies have unstable structures and long single-stranded regions. Using mass photometry, we show that RNAs with stable structures are mostly monomeric, while RNAs with unstable structures can multimerize. When mixed with protein, these differential RNA-RNA interactions affect dense phase compositions. Specifically, unstable RNAs form condensates with high RNA and low protein concentrations, while the converse is true for RNAs with stable structures. Using microbead rheology, we observe that unstable RNAs form condensates with long viscoelastic relaxation times that increase with condensate age, indicating elastic properties at many timescales. By contrast, stably-structured RNAs form less elastic condensates with shorter relaxation times that increase to a lesser degree with time. We hypothesize that such aging is controlled by the exchange of intra- for inter-molecular base pairing amongst RNAs whose folding free energies determine the timescale of these exchanges. Thus, RNA structure can have multi-scale consequences for biomolecular condensates, with impacts on condensate composition and material properties. This work underscores how specific RNA sequences can encode physical properties of condensates.

Double-stranded RNA reduction by chaotropic agents during in vitro transcription of messenger RNA

In-vitro-transcribed messenger RNA (mRNA) has recently shown increasing significance in the development of vaccines and therapeutics. Immunogenic double-stranded RNA (dsRNA) is an undesired byproduct formed during in vitro transcription (IVT), and it is challenging to reduce dsRNA byproduct from mRNA due to their similar sizes and intrinsic characteristics. Removal of dsRNA relies heavily on post-IVT chromatography purifications, such as reverse-phase high-pressure liquid chromatography, which increase manufacturing costs, reduce yield, and often decrease integrity, especially for long mRNA. Thus, it would be ideal to reduce and control the level of dsRNA during IVT. We herein present a simple, scalable, and controllable method to reduce the formation of dsRNA byproducts during IVT. Selected chaotropic agents at optimized concentrations are included during IVT to create a mild denaturing environment to prevent the undesired intermolecular or intramolecular base-pairing that is thought to promote RNA-templated dsRNA formation by RNA polymerase. Compared with regular IVT, our improved method produces mRNA with significantly less dsRNA, much lower immuno-stimulation, and more efficient protein expression. Therefore, this method potentially eliminates dsRNA removal purification steps and does not require reduced magnesium concentration, elevated temperature, or custom reagents, enabling a straightforward, high-yield, and cost-effective scale-up approach for mRNA manufacturing.

Stability and conformation of the dimeric HIV-1 genomic RNA 5′UTR

During HIV-1 assembly, the viral Gag polyprotein specifically selects the dimeric RNA genome for packaging into new virions. The 5′ untranslated region (5′UTR) of the dimeric genome may adopt a conformation that is optimal for recognition by Gag. Further conformational rearrangement of the 5′UTR, promoted by the nucleocapsid (NC) domain of Gag, is predicted during virus maturation. Two 5′UTR dimer conformations, the kissing dimer (KD) and the extended dimer (ED), have been identified in vitro, which differ in the extent of intermolecular basepairing. Whether 5′UTRs from different HIV-1 strains with distinct sequences have access to the same dimer conformations has not been determined. Here, we applied fluorescence cross-correlation spectroscopy and single-molecule Förster resonance energy transfer imaging to demonstrate that 5′UTRs from two different HIV-1 subtypes form (KDs) with divergent stabilities. We further show that both 5′UTRs convert to a stable dimer in the presence of the viral NC protein, adopting a conformation consistent with extensive intermolecular contacts. These results support a unified model in which the genomes of diverse HIV-1 strains adopt an ED conformation.

YB-1 unwinds mRNA secondary structures in vitro and negatively regulates stress granule assembly in HeLa cells.

In the absence of the scanning ribosomes that unwind mRNA coding sequences and 5′UTRs, mRNAs are likely to form secondary structures and intermolecular bridges. Intermolecular base pairing of non polysomal mRNAs is involved in stress granule (SG) assembly when the pool of mRNAs freed from ribosomes increases during cellular stress. Here, we unravel the structural mechanisms by which a major partner of dormant mRNAs, YB-1 (YBX1), unwinds mRNA secondary structures without ATP consumption by using its conserved cold-shock domain to destabilize RNA stem/loops and its unstructured C-terminal domain to secure RNA unwinding. At endogenous levels, YB-1 facilitates SG disassembly during arsenite stress recovery. In addition, overexpression of wild-type YB-1 and to a lesser extent unwinding-defective mutants inhibit SG assembly in HeLa cells. Through its mRNA-unwinding activity, YB-1 may thus inhibit SG assembly in cancer cells and package dormant mRNA in an unfolded state, thus preparing mRNAs for translation initiation.

RNAnue: efficient data analysis for RNA-RNA interactomics.

RNA–RNA inter- and intramolecular interactions are fundamental for numerous biological processes. While there are reasonable approaches to map RNA secondary structures genome-wide, understanding how different RNAs interact to carry out their regulatory functions requires mapping of intermolecular base pairs. Recently, different strategies to detect RNA–RNA duplexes in living cells, so called direct duplex detection (DDD) methods, have been developed. Common to all is the Psoralen-mediated in vivo RNA crosslinking followed by RNA Proximity Ligation to join the two interacting RNA strands. Sequencing of the RNA via classical RNA-seq and subsequent specialised bioinformatic analyses the result in the prediction of inter- and intramolecular RNA–RNA interactions. Existing approaches adapt standard RNA-seq analysis pipelines, but often neglect inherent features of RNA–RNA interactions that are useful for filtering and statistical assessment. Here we present RNAnue, a general pipeline for the inference of RNA–RNA interactions from DDD experiments that takes into account hybridisation potential and statistical significance to improve prediction accuracy. We applied RNAnue to data from different DDD studies and compared our results to those of the original methods. This showed that RNAnue performs better in terms of quantity and quality of predictions.

A Web Server for Designing Molecular Switches Composed of Two Interacting RNAs.

The programmability of RNA–RNA interactions through intermolecular base-pairing has been successfully exploited to design a variety of RNA devices that artificially regulate gene expression. An in silico design for interacting structured RNA sequences that satisfies multiple design criteria becomes a complex multi-objective problem. Although multi-objective optimization is a powerful technique that explores a vast solution space without empirical weights between design objectives, to date, no web service for multi-objective design of RNA switches that utilizes RNA–RNA interaction has been proposed. We developed a web server, which is based on a multi-objective design algorithm called MODENA, to design two interacting RNAs that form a complex in silico. By predicting the secondary structures with RactIP during the design process, we can design RNAs that form a joint secondary structure with an external pseudoknot. The energy barrier upon the complex formation is modeled by an interaction seed that is optimized in the design algorithm. We benchmarked the RNA switch design approaches (MODENA+RactIP and MODENA+RNAcofold) for the target structures based on natural RNA-RNA interactions. As a result, MODENA+RactIP showed high design performance for the benchmark datasets.

IntaRNAhelix-composing RNA-RNA interactions from stable inter-molecular helices boosts bacterial sRNA target prediction.

Efficient computational tools for the identification of putative target RNAs regulated by prokaryotic sRNAs rely on thermodynamic models of RNA secondary structures. While they typically predict RNA-RNA interaction complexes accurately, they yield many highly-ranked false positives in target screens. One obvious source of this low specificity appears to be the disability of current secondary-structure-based models to reflect steric constraints, which nevertheless govern the kinetic formation of RNA-RNA interactions. For example, often- even thermodynamically favorable- extensions of short initial kissing hairpin interactions are kinetically prohibited, since this would require unwinding of intra-molecular helices as well as sterically impossible bending of the interaction helix. Another source is the consideration of instable and thus unlikely subinteractions that enable better scoring of longer interactions. In consequence, the efficient prediction methods that do not consider such effects show a high false positive rate. To increase the prediction accuracy we devise IntaRNAhelix, a dynamic programming algorithm that length-restricts the runs of consecutive inter-molecular base pairs (perfect canonical stackings), which we hypothesize to implicitly model the steric and kinetic effects. The novel method is implemented by extending the state-of-the-art tool IntaRNA. Our comprehensive bacterial sRNA target prediction benchmark demonstrates significant improvements of the prediction accuracy and enables more than 40-times faster computations. These results indicate- supporting our hypothesis- that stable helix composition increases the accuracy of interaction prediction models compared to the current state-of-the-art approach.

Retroviral RNA Dimerization: From Structure to Functions.

The genome of the retroviruses is a dimer composed by two homologous copies of genomic RNA (gRNA) molecules of positive polarity. The dimerization process allows two gRNA molecules to be non-covalently linked together through intermolecular base-pairing. This step is critical for the viral life cycle and is highly conserved among retroviruses with the exception of spumaretroviruses. Furthermore, packaging of two gRNA copies into viral particles presents an important evolutionary advantage for immune system evasion and drug resistance. Recent studies reported RNA switches models regulating not only gRNA dimerization, but also translation and packaging, and a spatio-temporal characterization of viral gRNA dimerization within cells are now at hand. This review summarizes our current understanding on the structural features of the dimerization signals for a variety of retroviruses (HIVs, MLV, RSV, BLV, MMTV, MPMV…), the mechanisms of RNA dimer formation and functional implications in the retroviral cycle.

NMR resonance assignments of RNase P protein from Thermotoga maritima.

Ribonuclase P (RNase P) is an essential metallo-endonuclease that catalyzes 5' precursor-tRNA (ptRNA) processing and exists as an RNA-based enzyme in bacteria, archaea, and eukaryotes. In bacteria, a large catalytic RNA and a small protein component assemble to recognize and accurately cleave ptRNA and tRNA-like molecular scaffolds. Substrate recognition of ptRNA by bacterial RNase P requires RNA-RNA shape complementarity, intermolecular base pairing, and a dynamic protein-ptRNA binding interface. To gain insight into the binding specificity and dynamics of the bacterial protein-ptRNA interface, we report the backbone and side chain 1H, 13C, and 15N resonance assignments of the hyperthermophilic Thermatoga maritima RNase P protein in solution at 318K. Our data confirm the formation of a stable RNA recognition motif (RRM) with intrinsic heterogeneity at both the N- and C-terminus of the protein, consistent with available structural information. Comprehensive resonance assignments of the bacterial RNase P protein serve as an important first step in understanding how coupled RNA binding and protein-RNA conformational changes give rise to ribonucleoprotein function.

A homodimer interface without base pairs in an RNA mimic of red fluorescent protein.

Corn, a 28-nucleotide RNA, induces yellow fluorescence of its cognate ligand (3,5-difluoro-4-hydroxybenzylidene-imidazolinone-2-oxime, DFHO) by >1000-fold. It was selected in vitro to overcome limitations of other fluorogenic RNAs, particularly rapid photobleaching. We now report the Corn-DFHO co-crystal structure, discovering that the functional species is a quasisymmetric homodimer. Unusually, the dimer interface, where six unpaired adenosines break overall 2-fold symmetry, lacks any intermolecular base pairs. The homodimer encapsulates one DFHO at its inter-protomer interface, sandwiching it with a G-quadruplex from each protomer. Corn and the green-fluorescent Spinach RNA are structurally unrelated. Their convergent use of G-quadruplexes underscores the usefulness of this motif for RNA-induced small-molecule fluorescence. The asymmetric dimer interface of Corn could form the basis for the development of mutants that only fluoresce as heterodimers. Such variants would be analogous to Split GFP, and may be useful in analyzing RNA co-expression or association, or in designing self-assembling RNA nanostructures.

Mapping RNA-RNA Interactions Globally Using Biotinylated Psoralen

Knowing how RNAs interact with themselves and with others is key to understanding RNA based gene regulation in the cell. While examples of RNA-RNA interactions such as microRNA-mRNA interactions have been shown to regulate gene expression, the full extent to which RNA interactions occur in the cell is still unknown. Previous methods to study RNA interactions have primarily focused on subsets of RNAs that are interacting with a particular protein or RNA species. Here, we detail a method named Sequencing of Psoralen crosslinked, Ligated, and Selected Hybrids (SPLASH) that allows genome-wide capture of RNA interactions in vivo in an unbiased manner. SPLASH utilizes in vivo crosslinking, proximity ligation, and high throughput sequencing to identify intramolecular and intermolecular RNA base-pairing partners globally. SPLASH can be applied to different organisms including bacteria, yeast and human cells, as well as diverse cellular conditions to facilitate the understanding of the dynamics of RNA organization under diverse cellular contexts. The entire experimental SPLASH protocol takes about 5 days to complete and the computational workflow takes about 7 days to complete.

Mapping RNA-RNA Interactions Globally Using Biotinylated Psoralen.

Knowing how RNAs interact with themselves and with others is key to understanding RNA based gene regulation in the cell. While examples of RNA-RNA interactions such as microRNA-mRNA interactions have been shown to regulate gene expression, the full extent to which RNA interactions occur in the cell is still unknown. Previous methods to study RNA interactions have primarily focused on subsets of RNAs that are interacting with a particular protein or RNA species. Here, we detail a method named Sequencing of Psoralen crosslinked, Ligated, and Selected Hybrids (SPLASH) that allows genome-wide capture of RNA interactions in vivo in an unbiased manner. SPLASH utilizes in vivo crosslinking, proximity ligation, and high throughput sequencing to identify intramolecular and intermolecular RNA base-pairing partners globally. SPLASH can be applied to different organisms including bacteria, yeast and human cells, as well as diverse cellular conditions to facilitate the understanding of the dynamics of RNA organization under diverse cellular contexts. The entire experimental SPLASH protocol takes about 5 days to complete and the computational workflow takes about 7 days to complete.

NMR detection of intermolecular interaction sites in the dimeric 5′-leader of the HIV-1 genome

HIV type-1 (HIV-1) contains a pseudodiploid RNA genome that is selected for packaging and maintained in virions as a noncovalently linked dimer. Genome dimerization is mediated by conserved elements within the 5'-leader of the RNA, including a palindromic dimer initiation signal (DIS) that has been proposed to form kissing hairpin and/or extended duplex intermolecular contacts. Here, we have applied a 2H-edited NMR approach to directly probe for intermolecular interactions in the full-length, dimeric HIV-1 5'-leader (688 nucleotides; 230 kDa). The interface is extensive and includes DIS:DIS base pairing in an extended duplex state as well as intermolecular pairing between elements of the upstream Unique-5' (U5) sequence and those near the gag start site (AUG). Other pseudopalindromic regions of the leader, including the transcription activation (TAR), polyadenylation (PolyA), and primer binding (PBS) elements, do not participate in intermolecular base pairing. Using a 2H-edited one-dimensional NMR approach, we also show that the extended interface structure forms on a time scale similar to that of overall RNA dimerization. Our studies indicate that a kissing dimer-mediated structure, if formed, exists only transiently and readily converts to the extended interface structure, even in the absence of the HIV-1 nucleocapsid protein or other RNA chaperones.

Genome-Wide Analysis of RNA Secondary Structure

Single-stranded RNA molecules fold into extraordinarily complicated secondary and tertiary structures as a result of intramolecular base pairing. In vivo, these RNA structures are not static. Instead, they are remodeled in response to changes in the prevailing physicochemical environment of the cell and as a result of intermolecular base pairing and interactions with RNA-binding proteins. Remarkable technical advances now allow us to probe RNA secondary structure at single-nucleotide resolution and genome-wide, both in vitro and in vivo. These data sets provide new glimpses into the RNA universe. Analyses of RNA structuromes in HIV, yeast, Arabidopsis, and mammalian cells and tissues have revealed regulatory effects of RNA structure on messenger RNA (mRNA) polyadenylation, splicing, translation, and turnover. Application of new methods for genome-wide identification of mRNA modifications, particularly methylation and pseudouridylation, has shown that the RNA "epitranscriptome" both influences and is influenced by RNA structure. In this review, we describe newly developed genome-wide RNA structure-probing methods and synthesize the information emerging from their application.

A Method to Predict the Structure and Stability of RNA/RNA Complexes.

RNA/RNA interactions are essential for genomic RNA dimerization and regulation of gene expression. Intermolecular loop-loop base pairing is a widespread and functionally important tertiary structure motif in RNA machinery. However, computational prediction of intermolecular loop-loop base pairing is challenged by the entropy and free energy calculation due to the conformational constraint and the intermolecular interactions. In this chapter, we describe a recently developed statistical mechanics-based method for the prediction of RNA/RNA complex structures and stabilities. The method is based on the virtual bond RNA folding model (Vfold). The main emphasis in the method is placed on the evaluation of the entropy and free energy for the loops, especially tertiary kissing loops. The method also uses recursive partition function calculations and two-step screening algorithm for large, complicated structures of RNA/RNA complexes. As case studies, we use the HIV-1 Mal dimer and the siRNA/HIV-1 mutant (T4) to illustrate the method.

Co-evolution in a putative bundling signal of bluetongue and epizootic hemorrhagic disease viruses.

Bluetongue virus (BTV) and epizootic hemorrhagic disease virus (EHDV) possess a genome of 10 segmented RNAs (S1-S10), one copy of each of which is considered to be packaged in a virion. This selective packaging is thought to be mediated by supramolecular complex formation of the 10 RNAs, through intermolecular base pairing of complementary nucleotide sequences termed the bundling signal. Here, the whole genomic sequences of BTV and EHDV isolates were analyzed to identify co-evolving pairs of complementary nucleotide sequences within and between genomic segments. One co-evolving pair was identified within S5 and another between S5 and S10. The co-evolving pair between S5 and S10, consisting of six bases in each segment, was a candidate for a bundling signal and was identical to one of two putative bundling signals reported in a previous experimental study, validating the effectiveness of the method used in the present study. The six bases in S10 were confirmed to be located in a loop at the end of a stable stem. Although the six bases in S5 were located in a loop at the end of a stem of only four bases long, the complementary nucleotide sequences constituting this stem were, remarkably, the co-evolving pair within S5. These results highlight the importance not only of loops but also of stems in the intermolecular base pairing of bundling signals.

RNA World research-still evolving.

RNA World research-still evolving.