Recognition modes of RNA tetraloops and tetraloop-like motifs by RNA-binding proteins

Investigating the Intramolecular Competition of Different RNA Binding Motifs for Neomycin B by Native Top-Down Mass Spectrometry

The ongoing search for small molecule drugs that target ribonucleic acids (RNA) is complicated by a limited understanding of the principles that govern RNA-small molecule interactions. Here we have used stoichiometry-resolved native top-down mass spectrometry (MS) to study the binding of neomycin B to small model hairpin RNAs, an unstructured RNA, and a viral RNA construct. For 15-22 nt model RNAs with hairpin structure, we found that neomycin B binding to hairpin loops relies on interactions with both the nucleobases and the 2'-OH groups, and that a simple 5' or 3' overhang can introduce an additional binding motif. For a 47 nt RNA construct derived from stem IA of the human immunodeficiency virus 1 (HIV-1) rev response element (RRE) RNA, native top-down MS identified four different binding motifs, of which the purine-rich internal loop showed the highest affinity for neomycin B. Stoichiometry-resolved binding site mapping by native top-down MS allows for a new perspective on binding specificity, and has the potential to reveal unexpected principles of small molecule binding to RNA.

Development of Better Aptamers: Structured Library Approaches, Selection Methods, and Chemical Modifications

Systematic evolution of ligands by exponential enrichment (SELEX) has been used to discover thousands of aptamers since its development in 1990. Aptamers are short single-stranded oligonucleotides capable of binding to targets with high specificity and selectivity through structural recognition. While aptamers offer advantages over other molecular recognition elements such as their ease of production, smaller size, extended shelf-life, and lower immunogenicity, they have yet to show significant success in real-world applications. By analyzing the importance of structured library designs, reviewing different SELEX methodologies, and the effects of chemical modifications, we provide a comprehensive overview on the production of aptamers for applications in drug delivery systems, therapeutics, diagnostics, and molecular imaging.

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

The N6-methyladenosine modification is one of the most abundant post-transcriptional modifications in ribonucleic acid (RNA) molecules. Using molecular dynamics simulations and alchemical free-energy calculations, we studied the structural and energetic implications of incorporating this modification in an adenine mononucleotide and an RNA hairpin structure. At the mononucleotide level, we found that the syn configuration is more favorable than the anti configuration by 2.05 ± 0.15 kcal/mol. The unfavorable effect of methylation was due to the steric overlap between the methyl group and a nitrogen atom in the purine ring. We then probed the effect of methylation in an RNA hairpin structure containing an AUCG tetraloop, which is recognized by a "reader" protein (YTHDC1) to promote transcriptional silencing of long noncoding RNAs. While methylation had no significant conformational effect on the hairpin stem, the methylated tetraloop showed enhanced conformational flexibility compared to the unmethylated tetraloop. The increased flexibility was associated with the outward flipping of two bases (A6 and U7) which formed stacking interactions with each other and with the C8 and G9 bases in the tetraloop, leading to a conformation similar to that in the RNA/reader protein complex. Therefore, methylation-induced conformational flexibility likely facilitates RNA recognition by the reader protein.

Structural basis for cloverleaf RNA-initiated viral genome replication

The genomes of positive-strand RNA viruses serve as a template for both protein translation and genome replication. In enteroviruses, a cloverleaf RNA structure at the 5' end of the genome functions as a switch to transition from viral translation to replication by interacting with host poly(C)-binding protein 2 (PCBP2) and the viral 3CDpro protein. We determined the structures of cloverleaf RNA from coxsackievirus and poliovirus. Cloverleaf RNA folds into an H-type four-way junction and is stabilized by a unique adenosine-cytidine-uridine (A•C-U) base triple involving the conserved pyrimidine mismatch region. The two PCBP2 binding sites are spatially proximal and are located on the opposite end from the 3CDpro binding site on cloverleaf. We determined that the A•C-U base triple restricts the flexibility of the cloverleaf stem-loops resulting in partial occlusion of the PCBP2 binding site, and elimination of the A•C-U base triple increases the binding affinity of PCBP2 to the cloverleaf RNA. Based on the cloverleaf structures and biophysical assays, we propose a new mechanistic model by which enteroviruses use the cloverleaf structure as a molecular switch to transition from viral protein translation to genome replication.

Conserved Sequences in the 5' and 3' Untranslated Regions of Jingmenvirus Group Representatives

The Jingmenvirus group (JVG), with members such as Jingmen tick virus (JMTV), Alongshan virus (ALSV), Yanggou tick virus (YGTV), and Takachi virus (TAKV), is drawing attention due to evidence of it causing disease in humans and its unique genome architecture. In the current work, complete untranslated regions (UTRs) of four strains of ALSV and eight strains of YGTV were obtained. An analysis of these sequences, as well as JVG sequences from GenBank, uncovered several regions within viral UTRs that were highly conserved for all the segments and viruses. Bioinformatics predictions suggested that the UTRs of all the segments of YGTV, ALSV, and JMTV could form similar RNA structures. The most notable feature of these structures was a stable stem-loop with one (5' UTR) or two (3' UTR) AAGU tetraloops on the end of a hairpin.

Characterizing the Conformational Free-Energy Landscape of RNA Stem-Loops Using Single-Molecule Field-Effect Transistors

We have developed and used single-molecule field-effect transistors (smFETs) to characterize the conformational free-energy landscape of RNA stem-loops. Stem-loops are one of the most common RNA structural motifs and serve as building blocks for the formation of complex RNA structures. Given their prevalence and integral role in RNA folding, the kinetics of stem-loop (un)folding has been extensively characterized using both experimental and computational approaches. Interestingly, these studies have reported vastly disparate timescales of (un)folding, which has been interpreted as evidence that (un)folding of even simple stem-loops occurs on a highly rugged conformational energy landscape. Because smFETs do not rely on fluorophore reporters of conformation or mechanical (un)folding forces, they provide a unique approach that has allowed us to directly monitor tens of thousands of (un)folding events of individual stem-loops at a 200 μs time resolution. Our results show that under our experimental conditions, stem-loops (un)fold over a 1-200 ms timescale during which they transition between ensembles of unfolded and folded conformations, the latter of which is composed of at least two sub-populations. The 1-200 ms timescale of (un)folding we observe here indicates that smFETs report on complete (un)folding trajectories in which unfolded conformations of the RNA spend long periods of time wandering the free-energy landscape before sampling one of several misfolded conformations or the natively folded conformation. Our findings highlight the extremely rugged landscape on which even the simplest RNA structural elements fold and demonstrate that smFETs are a unique and powerful approach for characterizing the conformational free-energy of RNA.

Acylation probing of "generic" RNA libraries reveals critical influence of loop constraints on reactivity

The reactivity of RNA 2'-OH acylation is broadly useful both in probing structure and in preparing conjugates. To date, this reactivity has been analyzed in limited sets of biological RNA sequences, leaving open questions of how reactivity varies inherently without regard to sequence in structured contexts. We constructed and probed "generic" structured RNA libraries using homogeneous loop sequences, employing deep sequencing to carry out a systematic survey of reactivity. We find a wide range of RNA reactivities among single-stranded sequences, with nearest neighbors playing substantial roles. Remarkably, certain small loops are found to be far more reactive on average (up to 4,000-fold) than single-stranded RNAs, due to conformational constraints that enhance reactivity. Among loops, we observe large variations in reactivity based on size, type, and position. The results lend insights into RNA designs for achieving high-efficiency local conjugation and provide new opportunities to refine structure analysis.

Decoding the Mechanism of Specific RNA Targeting by Ribosomal Methyltransferases

Methylation of specific nucleotides is integral for ribosomal biogenesis and also serves as a common mechanism to confer antibiotic resistance by pathogenic bacteria. Here, by determining the high-resolution structure of the 30S-KsgA complex by cryo-electron microscopy, a state was captured, where KsgA juxtaposes between helices h44 and h45 of the 30S ribosome, separating them, thereby enabling remodeling of the surrounded rRNA and allowing the cognate site to enter the methylation pocket. With the structure as a guide, several mutant versions of the ribosomes, where interacting bases in the catalytic helix h45 and surrounding helices h44, h24, and h27, were mutated and evaluated for their methylation efficiency revealing factors that direct the enzyme to its cognate site with high fidelity. The biochemical studies show that the three-dimensional environment of the ribosome enables the interaction of select loop regions in KsgA with the ribosome helices paramount to maintain selectivity.

Sequence-specific RNA recognition by an RGG motif connects U1 and U2 snRNP for spliceosome assembly

In mammals, the structural basis for the interaction between U1 and U2 small nuclear ribonucleoproteins (snRNPs) during the early steps of splicing is still elusive. The binding of the ubiquitin-like (UBL) domain of SF3A1 to the stem-loop 4 of U1 snRNP (U1-SL4) contributes to this interaction. Here, we determined the 3D structure of the complex between the UBL of SF3A1 and U1-SL4 RNA. Our crystallography, NMR spectroscopy, and cross-linking mass spectrometry data show that SF3A1-UBL recognizes, sequence specifically, the GCG/CGC RNA stem and the apical UUCG tetraloop of U1-SL4. In vitro and in vivo mutational analyses support the observed intermolecular contacts and demonstrate that the carboxyl-terminal arginine-glycine-glycine-arginine (RGGR) motif of SF3A1-UBL binds sequence specifically by inserting into the RNA major groove. Thus, the characterization of the SF3A1-UBL/U1-SL4 complex expands the repertoire of RNA binding domains and reveals the capacity of RGG/RG motifs to bind RNA in a sequence-specific manner.

Characteristic chemical probing patterns of loop motifs improve prediction accuracy of RNA secondary structures

RNA structures play a fundamental role in nearly every aspect of cellular physiology and pathology. Gaining insights into the functions of RNA molecules requires accurate predictions of RNA secondary structures. However, the existing thermodynamic folding models remain less accurate than desired, even when chemical probing data, such as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reactivities, are used as restraints. Unlike most SHAPE-directed algorithms that only consider SHAPE restraints for base pairing, we extract two-dimensional structural features encoded in SHAPE data and establish robust relationships between characteristic SHAPE patterns and loop motifs of various types (hairpin, internal, and bulge) and lengths (2-11 nucleotides). Such characteristic SHAPE patterns are closely related to the sugar pucker conformations of loop residues. Based on these patterns, we propose a computational method, SHAPELoop, which refines the predicted results of the existing methods, thereby further improving their prediction accuracy. In addition, SHAPELoop can provide information about local or global structural rearrangements (including pseudoknots) and help researchers to easily test their hypothesized secondary structures.

In vivo structural characterization of the SARS-CoV-2 RNA genome identifies host proteins vulnerable to repurposed drugs

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 2019 (COVID-19) pandemic. Understanding of the RNA virus and its interactions with host proteins could improve therapeutic interventions for COVID-19. By using icSHAPE, we determined the structural landscape of SARS-CoV-2 RNA in infected human cells and from refolded RNAs, as well as the regulatory untranslated regions of SARS-CoV-2 and six other coronaviruses. We validated several structural elements predicted in silico and discovered structural features that affect the translation and abundance of subgenomic viral RNAs in cells. The structural data informed a deep-learning tool to predict 42 host proteins that bind to SARS-CoV-2 RNA. Strikingly, antisense oligonucleotides targeting the structural elements and FDA-approved drugs inhibiting the SARS-CoV-2 RNA binding proteins dramatically reduced SARS-CoV-2 infection in cells derived from human liver and lung tumors. Our findings thus shed light on coronavirus and reveal multiple candidate therapeutics for COVID-19 treatment.

Plausible Minimal Substrate for Erm Protein

Erm proteins methylate a specific adenine residue (A2058, Escherichia coli coordinates) conferring macrolide-lincosamide-streptogramin B (MLS_B) antibiotic resistance on a variety of microorganisms, ranging from antibiotic producers to pathogens. To identify the minimal motif required to be recognized and methylated by the Erm protein, various RNA substrates from 23S rRNA were constructed, and the substrate activity of these constructs was studied using three Erm proteins, namely, ErmB from Firmicutes and ErmE and ErmS from Actinobacteria.

Erm proteins methylate a specific adenine residue (A2058, Escherichia coli coordinates) conferring macrolide-lincosamide-streptogramin B (MLS_B) antibiotic resistance on a variety of microorganisms, ranging from antibiotic producers to pathogens. To identify the minimal motif required to be recognized and methylated by the Erm protein, various RNA substrates from 23S rRNA were constructed, and the substrate activity of these constructs was studied using three Erm proteins, namely, ErmB from Firmicutes and ErmE and ErmS from Actinobacteria. The shortest motif of 15 nucleotides (nt) could be recognized and methylated by ErmS, consisting of A2051 to the methylatable adenine (A2058) and its base-pairing counterpart strand, presumably assuming a quite similar structure to that in 23S rRNA, an unpaired target adenine immediately followed by an irregular double-stranded RNA region. This observation confirms the ultimate end of each side in helix 73 for methylation, determined by the approaches described above, and could reveal the mechanism behind the binding, recognition, induced fit, methylation, and conformational change for product release in the minimal context of substrate, presumably with the help of structural determination of the protein-RNA complex. In the course of determining the minimal portion of substrate from domain V, protein-specific features could be observed among the Erm proteins in terms of the methylation of RNA substrate and cooperativity and/or allostery between the region in helix 73 furthest away from the target adenine and the large portion of domain V above the methylatable adenine.

Structural insights into dimethylation of 12S rRNA by TFB1M: indispensable role in translation of mitochondrial genes and mitochondrial function

Mitochondria are essential molecular machinery for the maintenance of cellular energy supply by the oxidative phosphorylation system (OXPHOS). Mitochondrial transcription factor B1 (TFB1M) is a dimethyltransferase that maintains mitochondrial homeostasis by catalyzing dimethylation of two adjacent adenines located in helix45 (h45) of 12S rRNA. This m⁶₂A modification is indispensable for the assembly and maturation of human mitochondrial ribosomes. However, both the mechanism of TFB1M catalysis and the precise function of TFB1M in mitochondrial homeostasis are unknown. Here we report the crystal structures of a ternary complex of human (hs) TFB1M–h45–S-adenosyl-methionine and a binary complex hsTFB1M–h45. The structures revealed a distinct mode of hsTFB1M interaction with its rRNA substrate and with the initial enzymatic state involved in m⁶₂A modification. The suppression of hsTFB1M protein level or the overexpression of inactive hsTFB1M mutants resulted in decreased ATP production and reduced expression of components of the mitochondrial OXPHOS without affecting transcription of the corresponding genes and their localization to the mitochondria. Therefore, hsTFB1M regulated the translation of mitochondrial genes rather than their transcription via m⁶₂A modification in h45.

Developing and Characterization of Chemically Modified RNA Aptamers for Targeting Wild Type and Mutated c-KIT Receptor Tyrosine Kinases

The c-KIT receptor represents an attractive target for cancer therapy. Aptamers are emerging as a new promising class of nucleic acid therapeutics. In this study, a conventional SELEX approach was applied against the kinase domain of a group of c-KIT proteins (c-KIT^WT, c-KIT^D816V, and c-KIT^D816H) to select aptamers from a random RNA pool that can bind to the kinase domain of each target with high affinity and can selectively interfere with their kinase activities. Interestingly, our data indicated that one candidate aptamer, called V15, can specifically inhibit the in vitro kinase activity of mutant c-KIT^D816V with an IC₅₀ value that is 9-fold more potent than the sunitinib drug tested under the same conditions. Another aptamer, named as H5/V36, showed the potential to distinguish between the c-KIT kinases by modulating the phosphorylation activity of each in a distinct mechanism of action and in a different potency.

Identification of receptors for UNCG and GNRA Z-turns and their occurrence in rRNA

In contrast to GNRA tetraloop receptors that are common in RNA, receptors for the more thermostable UNCG loops have remained elusive for almost three decades. An analysis of all RNA structures with resolution ≤3.0 Å from the PDB allowed us to identify three previously unnoticed receptors for UNCG and GNRA tetraloops that adopt a common UNCG fold, named ‘Z-turn’ in agreement with our previously published nomenclature. These receptors recognize the solvent accessible second Z-turn nucleotide in different but specific ways. Two receptors participating in a complex network of tertiary interactions are associated with the rRNA UUCG and GAAA Z-turns capping helices H62 and H35a in rRNA large subunits. Structural comparison of fully assembled ribosomes and comparative sequence analysis of >6500 rRNA sequences helped us recognize that these motifs are almost universally conserved in rRNA, where they may contribute to organize the large subunit around the subdomain-IV four-way junction. The third UCCG receptor was identified in a rRNA/protein construct crystallized at acidic pH. These three non-redundant Z-turn receptors are relevant for our understanding of the assembly of rRNA and other long-non-coding RNAs, as well as for the design of novel folding motifs for synthetic biology.

Ageritin, a Ribotoxin from Poplar Mushroom ( Agrocybe aegerita) with Defensive and Antiproliferative Activities

Ribotoxins make up a group of extracellular rRNA endoribonucleases produced by ascomycetes that display cytotoxicity toward animal cells, having been proposed as insecticidal agents. Recently, the ribotoxin Ageritin has been isolated from the basidiomycetes Agrocybe aegerita (poplar mushroom), suggesting that ribotoxins are widely distributed among fungi. To gain insights into the protective properties of Ageritin against pathogens and its putative biotechnological applications, we have tested several biological activities of Ageritin, comparing them with those of the well-known ribotoxin α-sarcin, and we found that Ageritin displayed, in addition to the already reported activities, (i) antibacterial activity against Micrococcus lysodeikticus, (ii) activity against the tobacco mosaic virus RNA, (iii) endonuclease activity against a supercoiled plasmid, (iv) nuclease activity against genomic DNA, (v) cytotoxicity to COLO 320, HeLa, and Raji cells by promoting apoptosis, and (vi) antifungal activity against the green mold Penicillium digitatum. Therefore, Ageritin and α-sarcin can induce resistance not only to insects but also to viruses, bacteria, and fungi. The multiple biological activities of Ageritin could be exploited to improve resistance to different pathogens by engineering transgenic plants. Furthermore, the induction of cell death by different mechanisms turns these ribotoxins into useful tools for cancer therapy.

RNA structural analysis of the MYC mRNA reveals conserved motifs that affect gene expression

The MYC gene encodes a human transcription factor and proto-oncogene that is dysregulated in over half of all known cancers. To better understand potential post-transcriptional regulatory features affecting MYC expression, we analyzed secondary structures in the MYC mRNA using a program that is optimized for finding small locally-folded motifs with a high propensity for function. This was accomplished by calculating folding metrics across the MYC sequence using a sliding analysis window and generating unique consensus base pairing models weighted by their lower-than-random predicted folding energy. A series of 30 motifs were identified, primarily in the 5' and 3' untranslated regions, which show evidence of structural conservation and compensating mutations across vertebrate MYC homologs. This analysis was able to recapitulate known elements found within an internal ribosomal entry site, as well as discover a novel element in the 3' UTR that is unusually stable and conserved. This novel motif was shown to affect MYC expression, potentially via the modulation of miRNA target accessibility or other trans-regulatory factors. In addition to providing basic insights into mechanisms that regulate MYC expression, this study provides numerous, potentially druggable RNA targets for the MYC gene, which is considered “undruggable” at the protein level.

RNA Structural Differentiation: Opportunities with Pattern Recognition

Our awareness and appreciation of the many regulatory roles of RNA have dramatically increased in the past decade. This understanding, in addition to the impact of RNA in many disease states, has renewed interest in developing selective RNA-targeted small molecule probes. However, the fundamental guiding principles in RNA molecular recognition that could accelerate these efforts remain elusive. While high-resolution structural characterization can provide invaluable insight, examples of well-characterized RNA structures, not to mention small molecule:RNA complexes, remain limited. This Perspective provides an overview of the current techniques used to understand RNA molecular recognition when high-resolution structural information is unavailable. We will place particular emphasis on a new method, pattern recognition of RNA with small molecules (PRRSM), that provides rapid insight into critical components of RNA recognition and differentiation by small molecules as well as into RNA structural features.

Processive Antitermination

Transcription is a discontinuous process, where each nucleotide incorporation cycle offers a decision between elongation, pausing, halting, or termination. Many cis-acting regulatory RNAs, such as riboswitches, exert their influence over transcription elongation. Through such mechanisms, certain RNA elements can couple physiological or environmental signals to transcription attenuation, a process where cis-acting regulatory RNAs directly influence formation of transcription termination signals. However, through another regulatory mechanism called processive antitermination (PA), RNA polymerase can bypass termination sites over much greater distances than transcription attenuation. PA mechanisms are widespread in bacteria, although only a few classes have been discovered overall. Also, although traditional, signal-responsive riboswitches have not yet been discovered to promote PA, it is increasingly clear that small RNA elements are still oftentimes required. In some instances, small RNA elements serve as loading sites for cellular factors that promote PA. In other instances, larger, more complicated RNA elements participate in PA in unknown ways, perhaps even acting alone to trigger PA activity. These discoveries suggest that what is now needed is a systematic exploration of PA in bacteria, to determine how broadly these transcription elongation mechanisms are utilized, to reveal the diversity in their molecular mechanisms, and to understand the general logic behind their cellular applications. This review covers the known examples of PA regulatory mechanisms and speculates that they may be broadly important to bacteria.

A second RNA-binding protein is essential for ethanol tolerance provided by the bacterial OLE ribonucleoprotein complex

OLE (ornate, large, extremophilic) RNAs comprise a class of structured noncoding RNAs (ncRNAs) found in many extremophilic bacteria species. OLE RNAs constitute one of the longest and most widespread bacterial ncRNA classes whose major biochemical function remains unknown. In the Gram-positive alkaliphile Bacillus halodurans, OLE RNA is abundant, and localizes to the cell membrane by association with the transmembrane OLE-associated protein called OapA (formerly OAP). These characteristics, along with the well-conserved sequence and structural features of OLE RNAs, suggest that the OLE ribonucleoprotein (RNP) complex performs important biological functions. B. halodurans strains lacking OLE RNA (∆ole) or OapA (∆oapA) are less tolerant of cold (20 °C) and short-chain alcohols (e.g., ethanol). Here, we describe the effects of a mutant OapA (called PM1) that more strongly inhibits growth under cold or ethanol stress compared with strains lacking the oapA gene, even when wild-type OapA is present. This dominant-negative effect of PM1 is reversed by mutations that render OLE RNA nonfunctional. This finding demonstrates that the deleterious PM1 phenotype requires an intact RNP complex, and suggests that the complex has one or more additional undiscovered components. A genetic screen uncovered PM1 phenotype suppressor mutations in the ybzG gene, which codes for a putative RNA-binding protein of unknown biological function. We observe that YbzG protein (also called OapB) selectively binds OLE RNA in vitro, whereas a mutant version of the protein is not observed to bind OLE RNA. Thus, YbzG/OapB is an important component of the functional OLE RNP complex in B. halodurans.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Effect of GCAA stabilizing loops on three- and four-way intramolecular junctions

Tetraloops are a common way of changing the melting behavior of a DNA or RNA structure without changing the sequence of the stem. Because of the ubiquitous nature of tetraloops, our goal is to understand the effect a GCAA tetraloop, which belongs to the GNRA family of tetraloops, has on the unfolding thermodynamics of intramolecular junctions. Specifically, we have described the melting behavior of intramolecular three-way and four-way junctions where a T₅ loop has been replaced with a GCAA tetraloops in different positions. Their thermodynamic profiles, including Δn_Na⁺ and Δn_W, were analyzed based on the position of the tetraloop. We obtained between -16.7 and -27.5 kcal mol^-1 for all junctions studied. The experimental data indicates the influence of the GCAA tetraloop is primarily dictated by the native unfolding of the junction; if the tetraloop is placed on a stem that unfolds as a single domain when the tetraloop is not present, it will unfold as a single domain when the tetraloop is present but with a higher thermal stability. Conversely, if the tetraloop is placed on a stem which unfolds cooperatively with other stems when the tetraloop is not present, the tetraloop will increase the thermal stability of all the stems in the melting domain. The oligonucleotide structure and not the tetraloop itself affects ion uptake; three-way junctions do not gain an increase in ion uptake, but four-way junctions do. This is not the case for water immobilization, where the position of the tetraloop dictates the amount of water immobilized.

Insights into Structural and Mechanistic Features of Viral IRES Elements

Internal ribosome entry site (IRES) elements are cis-acting RNA regions that promote internal initiation of protein synthesis using cap-independent mechanisms. However, distinct types of IRES elements present in the genome of various RNA viruses perform the same function despite lacking conservation of sequence and secondary RNA structure. Likewise, IRES elements differ in host factor requirement to recruit the ribosomal subunits. In spite of this diversity, evolutionarily conserved motifs in each family of RNA viruses preserve sequences impacting on RNA structure and RNA–protein interactions important for IRES activity. Indeed, IRES elements adopting remarkable different structural organizations contain RNA structural motifs that play an essential role in recruiting ribosomes, initiation factors and/or RNA-binding proteins using different mechanisms. Therefore, given that a universal IRES motif remains elusive, it is critical to understand how diverse structural motifs deliver functions relevant for IRES activity. This will be useful for understanding the molecular mechanisms beyond cap-independent translation, as well as the evolutionary history of these regulatory elements. Moreover, it could improve the accuracy to predict IRES-like motifs hidden in genome sequences. This review summarizes recent advances on the diversity and biological relevance of RNA structural motifs for viral IRES elements.

Independent evolution of tetraloop in enterovirus oriL replicative element and its putative binding partners in virus protein 3C

Background

Enteroviruses are small non-enveloped viruses with a (+) ssRNA genome with one open reading frame. Enterovirus protein 3C (or 3CD for some species) binds the replicative element oriL to initiate replication. The replication of enteroviruses features a low-fidelity process, which allows the virus to adapt to the changing environment on the one hand, and requires additional mechanisms to maintain the genome stability on the other. Structural disturbances in the apical region of oriL domain d can be compensated by amino acid substitutions in positions 154 or 156 of 3C (amino acid numeration corresponds to poliovirus 3C), thus suggesting the co-evolution of these interacting sequences in nature. The aim of this work was to understand co-evolution patterns of two interacting replication machinery elements in enteroviruses, the apical region of oriL domain d and its putative binding partners in the 3C protein.

Methods

To evaluate the variability of the domain d loop sequence we retrieved all available full enterovirus sequences (>6, 400 nucleotides), which were present in the NCBI database on February 2017 and analysed the variety and abundance of sequences in domain d of the replicative element oriL and in the protein 3C.

Results

A total of 2,842 full genome sequences was analysed. The majority of domain d apical loops were tetraloops, which belonged to consensus YNHG (Y = U/C, N = any nucleotide, H = A/C/U). The putative RNA-binding tripeptide 154–156 (Enterovirus C 3C protein numeration) was less diverse than the apical domain d loop region and, in contrast to it, was species-specific.

Discussion

Despite the suggestion that the RNA-binding tripeptide interacts with the apical region of domain d, they evolve independently in nature. Together, our data indicate the plastic evolution of both interplayers of 3C-oriL recognition.

Effect of Loop Length and Sequence on the Stability of DNA Pyrimidine Triplexes with TAT Base Triplets

We report the thermodynamic contributions of loop length and loop sequence to the overall stability of DNA intramolecular pyrimidine triplexes. Two sets of triplexes were designed: in the first set, the C₅ loop closing the triplex stem was replaced with ^5'-CT_nC loops (n = 1-5), whereas in the second set, both the duplex and triplex loops were replaced with a ^5'-GCAA or ^5'-AACG tetraloop. For the triplexes with a ^5'-CT_nC loop, the triplex with five bases in the loop has the highest stability relative to the control. A loop length lower than five compromises the strength of the base-pair stacks without decreasing the thermal stability, leading to a decreased enthalpy, whereas an increase in the loop length leads to a decreased enthalpy and a higher entropic penalty. The incorporation of the GCAA loop yielded more stable triplexes, whereas the incorporation of AACG in the triplex loop yielded a less stable triplex due to an unfavorable enthalpy term. Thus, addition of the GCAA tetraloop can cause an increase in the thermodynamics of the triplex without affecting the sequence or melting behavior and may result in an additional layer of genetic regulation.

Biophysical properties, thermal stability and functional impact of 8-oxo-7,8-dihydroguanine on oligonucleotides of RNA-a study of duplex, hairpins and the aptamer for preQ1 as models

A better understanding of the effects that oxidative lesions have on RNA is of importance to understand their role in the development/progression of disease. 8-oxo-7,8-dihydroguanine was incorporated into RNA to understand its structural and functional impact on RNA:RNA and RNA:DNA duplexes, hairpins and pseudoknots. One to three modifications were incorporated into dodecamers of RNA [AAGAGGGAUGAC] resulting in thermal destabilization (ΔT_m – 10°C per lesion). Hairpins with tetraloops c-UUCG*-g* (8-10), a-ACCG-g* (11-12), c-UUG*G*-g* (13-16) and c-ACG*G*-g* (17-20) were modified and used to determine thermal stabilities, concluding that: (i) modifying the stem leads to destabilization unless adenosine is the opposing basepair of 8-oxoGua; (ii) modification at the loop is position- and sequence-dependent and varies from slight stabilization to large destabilization, in some cases leading to formation of other secondary structures (hairpin→duplex). Functional effects were established using the aptamer for preQ₁ as model. Modification at G5 disrupted the stem P1 and inhibited recognition of the target molecule 7-methylamino-7-deazaguanine (preQ₁). Modifying G11 results in increased thermal stability, albeit with a K_d 4-fold larger than its canonical analog. These studies show the capability of 8-oxoG to affect structure and function of RNA, resulting in distinct outcomes as a function of number and position of the lesion.

Structure and Dynamics of the Tetra-A Loop and (A-A)-U Sequence Motif within the Coliphage GA Replicase RNA Operator

The three-dimensional structure of a RNA hairpin containing the RNA operator binding site for bacteriophage GA coat protein is presented. The phage GA operator contains the asymmetric (A-A)-U sequence motif and is capped by a four-adenine (tetra-A) loop. The uridine of the (A-A)-U motif preferentially pairs with the 5'-proximal cross-strand adenine, and the 3'-proximal adenine stacks into the helix. The tetra-A loop is well-ordered with adenine residues 2-4 forming a 3' stack. This loop conformation stands in contrast to the structure of the 5'-AUUA loop of the related phage MS2 operator in which residues 1 and 2 form a 5' stack. The context dependence of the (A-A)-U sequence motif conformation was examined using structures of 76 unique occurrences from the Protein Data Bank. The motif almost always has one adenine bulged and the other adenine adopting an A-U base pair. In the case in which the (A-A)-U motif is flanked by only one Watson-Crick base pair, the adenine adjacent to the flanking base pair tends to bulge; 80% of motifs with a 3' flanking pair have a 3' bulged adenine, and 84% of motifs with a 5' flanking pair have a 5' bulged adenine. The frequencies of 3'- and 5'-proximal adenines bulging are 33 and 67%, respectively, when the (A-A)-U motif is flanked by base pairs on both sides. Although a 3' flanking cytidine correlates (88%) with bulging of the 5'-proximal adenine, no strict dependence on flanking nucleotide identity was identified for the 5' side.

Engineering Structurally Interacting RNA (sxRNA)

RNA-based three-way junctions (3WJs) are naturally occurring structures found in many functional RNA molecules including rRNA, tRNA, snRNA and ribozymes. 3WJs are typically characterized as resulting from an RNA molecule folding back on itself in cis but could also form in trans when one RNA, for instance a microRNA binds to a second structured RNA, such as a mRNA. Trans-3WJs can influence the final shape of one or both of the RNA molecules and can thus provide a means for modulating the availability of regulatory motifs including potential protein or microRNA binding sites. Regulatory 3WJs generated in trans represent a newly identified regulatory category that we call structurally interacting RNA or sxRNA for convenience. Here we show that they can be rationally designed using familiar cis-3WJ examples as a guide. We demonstrate that an sxRNA “bait” sequence can be designed to interact with a specific microRNA “trigger” sequence, creating a regulatable RNA-binding protein motif that retains its functional activity. Further, we show that when placed downstream of a coding sequence, sxRNA can be used to switch “ON” translation of that sequence in the presence of the trigger microRNA and the amount of translation corresponded with the amount of microRNA present.

Integrated structural biology to unravel molecular mechanisms of protein-RNA recognition

Recent advances in RNA sequencing technologies have greatly expanded our knowledge of the RNA landscape in cells, often with spatiotemporal resolution. These techniques identified many new (often non-coding) RNA molecules. Large-scale studies have also discovered novel RNA binding proteins (RBPs), which exhibit single or multiple RNA binding domains (RBDs) for recognition of specific sequence or structured motifs in RNA. Starting from these large-scale approaches it is crucial to unravel the molecular principles of protein-RNA recognition in ribonucleoprotein complexes (RNPs) to understand the underlying mechanisms of gene regulation. Structural biology and biophysical studies at highest possible resolution are key to elucidate molecular mechanisms of RNA recognition by RBPs and how conformational dynamics, weak interactions and cooperative binding contribute to the formation of specific, context-dependent RNPs. While large compact RNPs can be well studied by X-ray crystallography and cryo-EM, analysis of dynamics and weak interaction necessitates the use of solution methods to capture these properties. Here, we illustrate methods to study the structure and conformational dynamics of protein-RNA complexes in solution starting from the identification of interaction partners in a given RNP. Biophysical and biochemical techniques support the characterization of a protein-RNA complex and identify regions relevant in structural analysis. Nuclear magnetic resonance (NMR) is a powerful tool to gain information on folding, stability and dynamics of RNAs and characterize RNPs in solution. It provides crucial information that is complementary to the static pictures derived from other techniques. NMR can be readily combined with other solution techniques, such as small angle X-ray and/or neutron scattering (SAXS/SANS), electron paramagnetic resonance (EPR), and Förster resonance energy transfer (FRET), which provide information about overall shapes, internal domain arrangements and dynamics. Principles of protein-RNA recognition and current approaches are reviewed and illustrated with recent studies.

LoaP is a broadly conserved antiterminator protein that regulates antibiotic gene clusters in Bacillus amyloliquefaciens

A valuable resource available in the search for new natural products is the diverse microbial life that spans the planet. A large subset of these microorganisms synthesize complex specialized metabolites exhibiting biomedically important activities. A limiting step to the characterization of these compounds is an elucidation of the genetic regulatory mechanisms that oversee their production. Although proteins that control transcription initiation of specialized metabolite gene clusters have been identified, those affecting transcription elongation have not been broadly investigated. In this study, we analysed the phylogenetic distribution of the large, widespread NusG family of transcription elongation proteins and found that it includes a cohesive outgroup of paralogues (herein coined LoaP), which are often positioned adjacent or within gene clusters for specialized metabolites. We established Bacillus amyloliquefaciens LoaP as a paradigm for this protein subgroup and showed that it regulated the transcriptional readthrough of termination sites located within two different antibiotic biosynthesis operons. Both of these antibiotics have been implicated in plant-protective activities, demonstrating that LoaP controls an important regulon of specialized metabolite genes for this microorganism. These data therefore reveal transcription elongation as a point of regulatory control for specialized metabolite pathways and introduce a subgroup of NusG proteins for this purpose.

Specific Recognition of a Single-Stranded RNA Sequence by a Synthetic Antibody Fragment

Antibodies that bind RNA represent an unrealized source of reagents for synthetic biology and for characterizing cellular transcriptomes. However, facile access to RNA-binding antibodies requires the engineering of effective Fab libraries guided by the knowledge of the principles that govern RNA recognition. Here, we describe a Fab identified from a minimalist synthetic library during phage display against a branched RNA target. The Fab (BRG) binds with 20nM dissociation constant to a single-stranded RNA (ssRNA) sequence adjacent to the branch site and can block the action of debranchase enzyme. We report the crystal structure in complex with RNA target at 2.38Å. The Fab traps the RNA in a hairpin conformation that contains a 2-bp duplex capped by a tetraloop. The paratope surface consists of residues located in four complementarity-determining regions including a major contribution from H3, which adopts a helical structure that projects into a deep, wide groove formed by the RNA. The amino acid composition of the paratope reflects the library diversity, consisting mostly of tyrosine and serine residues and a small but significant contribution from a single arginine residue. This structure, involving the recognition of ssRNA via a stem-loop conformation, together with our two previous structures involving the recognition of an RNA hairpin loop and an RNA tertiary structure, reveals the capacity of minimalist libraries biased with tyrosine, serine, glycine, and arginine to form binding surfaces for specific RNA conformations and distinct levels of RNA structural hierarchy.

PCBP2 enables the cadicivirus IRES to exploit the function of a conserved GRNA tetraloop to enhance ribosomal initiation complex formation

The cadicivirus IRES diverges structurally from canonical Type 1 IRESs (e.g. poliovirus) but nevertheless also contains an essential GNRA tetraloop in a subdomain (d10c) that is homologous to poliovirus dIVc. In addition to canonical initiation factors, the canonical Type 1 and divergent cadicivirus IRESs require the same IRES trans-acting factor, poly(C)-binding protein 2 (PCBP2). PCBP2 has three KH domains and binds poliovirus IRES domain dIV in the vicinity of the tetraloop. How PCBP2 binds the cadicivirus IRES, and the roles of PCBP2 and the tetraloop in Type 1 IRES function are unknown. Here, directed hydroxyl radical probing showed that KH1 also binds near the cadicivirus tetraloop. KH2 and KH3 bind adjacently to an IRES subdomain (d10b) that is unrelated to dIV, with KH3 in an inverted orientation. KH3 is critical for PCBP2's binding to this IRES whereas KH1 is essential for PCBP2's function in promoting initiation. PCBP2 enforced the wild-type structure of d10c when it contained minor destabilizing substitutions, exposing the tetraloop. Strikingly, PCBP2 enhanced initiation on mutant IRESs that retained consensus GNRA tetraloops, whereas mutants with divergent sequences did not respond to PCBP2. These studies show that PCBP2 enables the IRES to exploit the GNRA tetraloop to enhance initiation.

The RNA 3D Motif Atlas: Computational methods for extraction, organization and evaluation of RNA motifs

RNA 3D motifs occupy places in structured RNA molecules that correspond to the hairpin, internal and multi-helix junction "loops" of their secondary structure representations. As many as 40% of the nucleotides of an RNA molecule can belong to these structural elements, which are distinct from the regular double helical regions formed by contiguous AU, GC, and GU Watson-Crick basepairs. With the large number of atomic- or near atomic-resolution 3D structures appearing in a steady stream in the PDB/NDB structure databases, the automated identification, extraction, comparison, clustering and visualization of these structural elements presents an opportunity to enhance RNA science. Three broad applications are: (1) identification of modular, autonomous structural units for RNA nanotechnology, nanobiology and synthetic biology applications; (2) bioinformatic analysis to improve RNA 3D structure prediction from sequence; and (3) creation of searchable databases for exploring the binding specificities, structural flexibility, and dynamics of these RNA elements. In this contribution, we review methods developed for computational extraction of hairpin and internal loop motifs from a non-redundant set of high-quality RNA 3D structures. We provide a statistical summary of the extracted hairpin and internal loop motifs in the most recent version of the RNA 3D Motif Atlas. We also explore the reliability and accuracy of the extraction process by examining its performance in clustering recurrent motifs from homologous ribosomal RNA (rRNA) structures. We conclude with a summary of remaining challenges, especially with regard to extraction of multi-helix junction motifs.

MATtrack: A MATLAB-Based Quantitative Image Analysis Platform for Investigating Real-Time Photo-Converted Fluorescent Signals in Live Cells

We introduce here MATtrack, an open source MATLAB-based computational platform developed to process multi-Tiff files produced by a photo-conversion time lapse protocol for live cell fluorescent microscopy. MATtrack automatically performs a series of steps required for image processing, including extraction and import of numerical values from Multi-Tiff files, red/green image classification using gating parameters, noise filtering, background extraction, contrast stretching and temporal smoothing. MATtrack also integrates a series of algorithms for quantitative image analysis enabling the construction of mean and standard deviation images, clustering and classification of subcellular regions and injection point approximation. In addition, MATtrack features a simple user interface, which enables monitoring of Fluorescent Signal Intensity in multiple Regions of Interest, over time. The latter encapsulates a region growing method to automatically delineate the contours of Regions of Interest selected by the user, and performs background and regional Average Fluorescence Tracking, and automatic plotting. Finally, MATtrack computes convenient visualization and exploration tools including a migration map, which provides an overview of the protein intracellular trajectories and accumulation areas. In conclusion, MATtrack is an open source MATLAB-based software package tailored to facilitate the analysis and visualization of large data files derived from real-time live cell fluorescent microscopy using photoconvertible proteins. It is flexible, user friendly, compatible with Windows, Mac, and Linux, and a wide range of data acquisition software. MATtrack is freely available for download at eleceng.dit.ie/courtney/MATtrack.zip.

Modulating RNA Alignment Using Directional Dynamic Kinks: Application in Determining an Atomic-Resolution Ensemble for a Hairpin using NMR Residual Dipolar Couplings

Approaches that combine experimental data and computational molecular dynamics (MD) to determine atomic resolution ensembles of biomolecules require the measurement of abundant experimental data. NMR residual dipolar couplings (RDCs) carry rich dynamics information, however, difficulties in modulating overall alignment of nucleic acids have limited the ability to fully extract this information. We present a strategy for modulating RNA alignment that is based on introducing variable dynamic kinks in terminal helices. With this strategy, we measured seven sets of RDCs in a cUUCGg apical loop and used this rich data set to test the accuracy of an 0.8 μs MD simulation computed using the Amber ff10 force field as well as to determine an atomic resolution ensemble. The MD-generated ensemble quantitatively reproduces the measured RDCs, but selection of a sub-ensemble was required to satisfy the RDCs within error. The largest discrepancies between the RDC-selected and MD-generated ensembles are observed for the most flexible loop residues and backbone angles connecting the loop to the helix, with the RDC-selected ensemble resulting in more uniform dynamics. Comparison of the RDC-selected ensemble with NMR spin relaxation data suggests that the dynamics occurs on the ps-ns time scales as verified by measurements of R(1ρ) relaxation-dispersion data. The RDC-satisfying ensemble samples many conformations adopted by the hairpin in crystal structures indicating that intrinsic plasticity may play important roles in conformational adaptation. The approach presented here can be applied to test nucleic acid force fields and to characterize dynamics in diverse RNA motifs at atomic resolution.

Structural basis for specific recognition of single-stranded RNA by Toll-like receptor 13

Toll-like receptors (TLRs) have crucial roles in innate immunity, functioning as pattern-recognition receptors. TLR13 recognizes a conserved sequence from bacterial 23S rRNA and then triggers an immune response. Here we report the crystal structure of the mouse TLR13 ectodomain bound by a 13-nt single-stranded (ss) RNA derived from 23S rRNA. The ssRNA induces TLR13 dimerization but assumes a stem-loop-like structure that is completely different from that in the bacterial ribosome but nevertheless is crucial for TLR13 recognition. Most of the RNA nucleotides are splayed out to make base-specific contacts with the concave surface of TLR13, and RNA-specific interactions are important to allow TLR13 to distinguish RNA from DNA. Interestingly, a viral-derived 16-nt ssRNA predicted to form a similar stem-loop-like structure also induces TLR13 activation. Together, our results reveal the structural mechanism of TLR13's sequence- and conformation-specific recognition of ssRNA.

Hydroxylated histidine of human ribosomal protein uL2 is involved in maintaining the local structure of 28S rRNA in the ribosomal peptidyl transferase center

Protein uL2 is essential for the catalytic activity of the ribosome and has a conserved shape in ribosomes from all domains of life. However, the sequence of its unstructured C-terminal loop apex that contacts the conserved 23S/28S rRNA helix (H) 93 near the ribosomal peptidyl transferase center differs in bacteria, archaea and eukaryotes. Eukaryote-specific residue His216 located in this loop in mammalian uL2 is hydroxylated in ribosomes. We used a set of chemical probes to explore the structure of an RNA that mimicked a segment of 28S rRNA domain V containing part of the uL2 binding site including H93, complexed with either natural (hydroxylated) or recombinant (unmodified) human uL2. It was found that both protein forms engage H93 during binding, but only natural uL2 (uL2n) protects it from hydroxyl radicals. The association of uL2n with RNA leads to changes in its structure at U4532 adjacent to the universally conserved U4531 (U2585, Escherichia coli numbering) involved in peptidyl transferase center formation, and at the universally conserved C4447 (2501) located in the ribosome near A4397 (2451) and C3909 (2063) belonging to the peptidyl transferase center. As a result, both nucleotides become strongly exposed to hydroxyl radicals. Our data argue that the hydroxyl group at His216 in the C-terminal loop apex of mammalian uL2 contributes to stabilization of a protein conformation that is favorable for binding to H93 of 28S rRNA and that this binding induces structural rearrangement in the regions close to the peptidyl transferase center in the mature ribosome.

Structure-specific nucleic acid recognition by L-motifs and their diverse roles in expression and regulation of the genome

The high-mobility group (HMG) domain containing proteins regulate transcription, DNA replication and recombination. They adopt L-shaped folds and are structure-specific DNA binding motifs. Here, I define the L-motif super-family that consists of DNA-binding HMG-box proteins and the L-motif of the histone mRNA binding domain of stem-loop binding protein (SLBP). The SLBP L-motif and HMG-box domains adopt similar L-shaped folds with three α-helices and two or three small hydrophobic cores that stabilize the overall fold, but have very different and distinct modes of nucleic acid recognition. A comparison of the structure, dynamics, protein-protein and nucleic acid interactions, and regulation by PTMs of the SLBP and the HMG-box L-motifs reveals the versatile and diverse modes by which L-motifs utilize their surfaces for structure-specific recognition of nucleic acids to regulate gene expression.

RNA sociology: group behavioral motifs of RNA consortia

RNA sociology investigates the behavioral motifs of RNA consortia from the social science perspective. Besides the self-folding of RNAs into single stem loop structures, group building of such stem loops results in a variety of essential agents that are highly active in regulatory processes in cellular and non-cellular life. RNA stem loop self-folding and group building do not depend solely on sequence syntax; more important are their contextual (functional) needs. Also, evolutionary processes seem to occur through RNA stem loop consortia that may act as a complement. This means the whole entity functions only if all participating parts are coordinated, although the complementary building parts originally evolved for different functions. If complementary groups, such as rRNAs and tRNAs, are placed together in selective pressure contexts, new evolutionary features may emerge. Evolution initiated by competent agents in natural genome editing clearly contrasts with statistical error replication narratives.

Contribution of protein phosphorylation to binding-induced folding of the SLBP-histone mRNA complex probed by phosphorus-31 NMR

•

SLBP is an intrinsically disordered protein (IDP) in the absence of RNA.

•

A phosphothreonine in the SLBP RNA-binding domain stabilizes the SLBP–histone mRNA complex.

•

This phosphate exhibits torsional strain as revealed by its ³¹P NMR chemical shift.

•

Phosphates can play structural roles in stabilizing tertiary structure in IDPs.

•

³¹P NMR can be a good spectroscopic probe for folding of phosphorylated IDPs.

Phosphorus-31 (³¹P) NMR can be used to characterize the structure and dynamics of phosphorylated proteins. Here, I use ³¹P NMR to report on the chemical nature of a phosphothreonine that lies in the RNA binding domain of SLBP (stem-loop binding protein). SLBP is an intrinsically disordered protein and phosphorylation at this threonine promotes the assembly of the SLBP–RNA complex. The data show that the ³¹P chemical shift can be a good spectroscopic probe for phosphate-coupled folding and binding processes in intrinsically disordered proteins, particularly where the phosphate exhibits torsional strain and is involved in a network of hydrogen-bonding interactions.

Trans-regulation of RNA-binding protein motifs by microRNA

The wide array of vital functions that RNA performs is dependent on its ability to dynamically fold into different structures in response to intracellular and extracellular changes. RNA-binding proteins regulate much of this activity by targeting specific RNA structures or motifs. One of these structures, the 3-way RNA junction, is characteristically found in ribosomal RNA and results from the RNA folding in cis, to produce three separate helices that meet around a central unpaired region. Here we demonstrate that 3-way junctions can also form in trans as a result of the binding of microRNAs in an unconventional manner with mRNA by splinting two non-contiguous regions together. This may be used to reinforce the base of a stem-loop motif being targeted by an RNA-binding protein. Trans interactions between non-coding RNA and mRNA may be used to control the post-transcriptional regulatory code and suggests a possible role for some of the recently described transcripts of unknown function expressed from the human genome.