Dissecting the energetic architecture within an RNA tertiary structural motif via high-throughput thermodynamic measurements.

Structured RNAs and RNA/protein complexes perform critical cellular functions. They often contain structurally conserved tertiary contact "motifs," whose occurrence simplifies the RNA folding landscape. Prior studies have focused on the conformational and energetic modularity of intact motifs. Here, we turn to the dissection of one common motif, the 11nt receptor (11ntR), using quantitative analysis of RNA on a massively parallel array to measure the binding of all single and double 11ntR mutants to GAAA and GUAA tetraloops, thereby probing the energetic architecture of the motif. While the 11ntR behaves as a motif, its cooperativity is not absolute. Instead, we uncovered a gradient from high cooperativity amongst base-paired and neighboring residues to additivity between distant residues. As expected, substitutions at residues in direct contact with the GAAA tetraloop resulted in the largest decreases to binding, and energetic penalties of mutations were substantially smaller for binding to the alternate GUAA tetraloop, which lacks tertiary contacts present with the canonical GAAA tetraloop. However, we found that the energetic consequences of base partner substitutions are not, in general, simply described by base pair type or isostericity. We also found exceptions to the previously established stability-abundance relationship for 11ntR sequence variants. These findings of "exceptions to the rule" highlight the power of systematic high-throughput approaches to uncover novel variants for future study in addition to providing an energetic map of a functional RNA.

In vivo secondary structural analysis of Influenza A virus genomic RNA.

Influenza A virus (IAV) is a respiratory virus that causes epidemics and pandemics. Knowledge of IAV RNA secondary structure in vivo is crucial for a better understanding of virus biology. Moreover, it is a fundament for the development of new RNA-targeting antivirals. Chemical RNA mapping using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) coupled with Mutational Profiling (MaP) allows for the thorough examination of secondary structures in low-abundance RNAs in their biological context. So far, the method has been used for analyzing the RNA secondary structures of several viruses including SARS-CoV-2 in virio and in cellulo. Here, we used SHAPE-MaP and dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) for genome-wide secondary structure analysis of viral RNA (vRNA) of the pandemic influenza A/California/04/2009 (H1N1) strain in both in virio and in cellulo environments. Experimental data allowed the prediction of the secondary structures of all eight vRNA segments in virio and, for the first time, the structures of vRNA5, 7, and 8 in cellulo. We conducted a comprehensive structural analysis of the proposed vRNA structures to reveal the motifs predicted with the highest accuracy. We also performed a base-pairs conservation analysis of the predicted vRNA structures and revealed many highly conserved vRNA motifs among the IAVs. The structural motifs presented herein are potential candidates for new IAV antiviral strategies.

New Short RNA Motifs Potentially Relevant in the SARS-CoV-2 Genome.

Background: The coronavirus disease has led to an exhaustive exploration of the SARS-CoV-2 genome. Despite the amount of information accumulated, the prediction of short RNA motifs encoding peptides mediating protein-protein or protein-drug interactions has received limited attention. Objective: The study aims to predict short RNA motifs that are interspersed in the SARS-CoV-2 genome. Methods: A method in which 14 trinucleotide families, each characterized by being composed of triplets with identical nucleotides in all possible configurations, was used to find short peptides with biological relevance. The novelty of the approach lies in using these families to search how they are distributed across genomes of different CoV genera and then to compare the distributions of these families with each other. Results: We identified distributions of trinucleotide families in different CoV genera and also how they are related, using a selection criterion that identified short RNA motifs. The motifs were reported to be conserved in SARS-CoVs; in the remaining CoV genomes analysed, motifs contained, exclusively, different configurations of the trinucleotides A, T, G and A, C, G. Eighty-eight short RNA motifs, ranging in length from 12 to 49 nucleotides, were found: 50 motifs in the 1a polyprotein-encoding orf, 27 in the 1b polyprotein-encoding orf, 5 in the spike-encoding orf, and 6 in the nucleocapsid-encoding orf. Although some motifs (~27%) were found to be intercalated or attached to functional peptides, most of them have not yet been associated with any known functions. Conclusion: Some of the trinucleotide family distributions in different CoV genera are not random; they are present in short peptides that, in many cases, are intercalated or attached to functional sites of the proteome.

Multiple Non-Canonical Base-Stacking Interactions as One of the Major Determinants of RNA Tertiary Structure Organization.

Stacking interactions of heterocyclic bases of ribonucleotides are one of the most important factors in the organization of RNA secondary and tertiary structure. Most of these (canonical) interactions are formed between adjacent residues in RNA polynucleotide chains. However, with the accumulation of data on the atomic tertiary structures of various RNAs and their complexes with proteins, it has become clear that nucleotide residues that are not adjacent in the polynucleotide chains and are sometimes separated in the RNA primary structure by tens or hundreds of nucleotides can interact via (non-canonical) base stacking. This paper presents an exhaustive database of such nonadjacent base-stacking elements (NA-BSEs) and their environment in the macromolecules of natural and synthetic RNAs. Analysis of these data showed that NA-BSE-forming nucleotides, on average, account for about a quarter of all nucleotides in a particular RNA and, therefore, should be considered as bona fide motifs of the RNA tertiary structure. We also classified NA-BSEs by their location in RNA macromolecules. It was shown that the structure-forming role of NA-BSEs involves compact folding of single-stranded RNA loops, transformation of double-stranded bulges into imperfect helices, and binding of RNA regions distant in the primary and secondary RNA structure.

RNA-As-Graphs Motif Atlas-Dual Graph Library of RNA Modules and Viral Frameshifting-Element Applications.

RNA motif classification is important for understanding structure/function connections and building phylogenetic relationships. Using our coarse-grained RNA-As-Graphs (RAG) representations, we identify recurrent dual graph motifs in experimentally solved RNA structures based on an improved search algorithm that finds and ranks independent RNA substructures. Our expanded list of 183 existing dual graph motifs reveals five common motifs found in transfer RNA, riboswitch, and ribosomal 5S RNA components. Moreover, we identify three motifs for available viral frameshifting RNA elements, suggesting a correlation between viral structural complexity and frameshifting efficiency. We further partition the RNA substructures into 1844 distinct submotifs, with pseudoknots and junctions retained intact. Common modules are internal loops and three-way junctions, and three submotifs are associated with riboswitches that bind nucleotides, ions, and signaling molecules. Together, our library of existing RNA motifs and submotifs adds to the growing universe of RNA modules, and provides a resource of structures and substructures for novel RNA design.

Role of RNA Motifs in RNA Interaction with Membrane Lipid Rafts: Implications for Therapeutic Applications of Exosomal RNAs.

RNA motifs may promote interactions with exosomes (EXO-motifs) and lipid rafts (RAFT-motifs) that are enriched in exosomal membranes. These interactions can promote selective RNA loading into exosomes. We quantified the affinity between RNA aptamers containing various EXO- and RAFT-motifs and membrane lipid rafts in a liposome model of exosomes by determining the dissociation constants. Analysis of the secondary structure of RNA molecules provided data about the possible location of EXO- and RAFT-motifs within the RNA structure. The affinity of RNAs containing RAFT-motifs (UUGU, UCCC, CUCC, CCCU) and some EXO-motifs (CCCU, UCCU) to rafted liposomes is higher in comparison to aptamers without these motifs, suggesting direct RNA-exosome interaction. We have confirmed these results through the determination of the dissociation constant values of exosome-RNA aptamer complexes. RNAs containing EXO-motifs GGAG or UGAG have substantially lower affinity to lipid rafts, suggesting indirect RNA-exosome interaction via RNA binding proteins. Bioinformatics analysis revealed RNA aptamers containing both raft- and miRNA-binding motifs and involvement of raft-binding motifs UCCCU and CUCCC. A strategy is proposed for using functional RNA aptamers (fRNAa) containing both RAFT-motif and a therapeutic motif (e.g., miRNA inhibitor) to selectively introduce RNAs into exosomes for fRNAa delivery to target cells for personalized therapy.

Dissecting the hybridization of oligonucleotides to structured complementary sequences.

BACKGROUND: When oligonucleotides hybridize to long target molecules, the process is slowed by the secondary structure in the targets. The phenomenon has been analyzed in several previous studies, but many details remain poorly understood. METHODS: I used a spectrofluorometric strategy, focusing on the formation/breaking of individual base pairs, to study the kinetics of association between a DNA hairpin and >20 complementary oligonucleotides ('antisenses'). RESULTS: Hybridization rates differed by over three orders of magnitude. Association was toehold-mediated, both for antisenses binding to the target's ends and for those designed to interact with the loop. Binding of these latter, besides being consistently slower, was affected to variable, non-uniform extents by the asymmetric loop structure. Divalent metal ions accelerated hybridization, more pronouncedly when nucleation occurred at the loop. Incorporation of locked nucleic acid (LNA) residues in the antisenses substantially improved the kinetics only when LNAs participated to the earliest hybridization steps. The effects of individual LNAs placed along the antisense indicated that the reaction transition state occurred after invading at least the first base pair of the stem. CONCLUSIONS: The experimental approach helps dissect hybridization reactions involving structured nucleic acids. Toehold-dependent, nucleation-invasion models appear fully appropriate for describing such reactions. Estimating the stability of nucleation complexes formed at internal toeholds is the major hurdle for the quantitative prediction of hybridization rates. GENERAL SIGNIFICANCE: While analyzing the mechanisms of a fundamental biochemical process (hybridization), this work also provides suggestions for the improvement of technologies that rely on such process.

Ring Catalog: A resource for designing self-assembling RNA nanostructures.

Designing self-assembling RNA ring structures based on known 3D structural elements connected via linker helices is a challenging task due to the immense number of motif combinations, many of which do not lead to ring-closure. We describe an in silico solution to this design problem by combinatorial assembly of RNA 3-way junctions, bulges, and kissing loops, and tabulating the cases that lead to ring formation. The solutions found are made available in the form of a web-accessible Ring Catalog. As an example of a potential use of this resource, we chose a predicted RNA square structure consisting of five RNA strands and demonstrate experimentally that the self-assembly of those five strands leads to the formation of a square-like complex. This is a demonstration of a novel "design by catalog" approach to RNA nano-structure generation. The URL https://rnajunction.ncifcrf.gov/ringdb can be used to access the resource.

Untranslated regions (UTRs) orchestrate translation reprogramming in cellular stress responses.

Stress is the result of an organism's interaction with environmental challenges. Regulations of gene expression including translation modulations are critical for adaptation and survival under stress. Untranslated regions (UTRs) of the transcripts play significant roles in translation regulation and continue to raise many intriguing questions in our understanding of cellular stress physiology. IRES (Internal ribosome entry site) and uORF (upstream open reading frame) mediated alternative translation initiations are emerging as unique mechanisms. Recent studies have revealed novel means of mRNAs stabilization in stress granules and their reversible modifications. Differential regulation of select transcripts is possible by the interplay between the adenine/uridine-rich elements (AREs) in 3'UTR with their binding proteins (AUBP) and by microRNA-mediated effects. Coordination of these various mechanisms control translation and thereby enables appropriate responses to environmental stress. In this review, we focus on the role of sequence signatures both at 5' and 3'UTRs in translation reprogramming during cellular stress responses.

Exploring the functionality and conservation of Alba proteins in Trypanosoma cruzi: A focus on biological diversity and RNA binding ability.

Trypanosoma cruzi is the causative agent of Chagas disease, as well as a trypanosomatid parasite with a complex biological cycle that requires precise mechanisms for regulating gene expression. In Trypanosomatidae, gene regulation occurs mainly at the mRNA level through the recognition of cis elements by RNA-binding proteins (RBPs). Alba family members are ubiquitous DNA/RNA-binding proteins with representatives in trypanosomatid parasites functionally related to gene expression regulation. Although T. cruzi possesses two groups of Alba proteins (Alba1/2 and Alba30/40), their functional role remains poorly understood. Thus, herein, a characterization of T. cruzi Alba (TcAlba) proteins was undertaken. Physicochemical, structural, and phylogenetic analysis of TcAlba showed features compatible with RBPs, such as hydrophilicity, RBP domains/motifs, and evolutionary conservation of the Alba-domain, mainly regarding other trypanosomatid Alba. However, in silico RNA interaction analysis of T. cruzi Alba proteins showed that TcAlba30/40 proteins, but not TcAlba1/2, would directly interact with the assayed RNA molecules, suggesting that these two groups of TcAlba proteins have different targets. Given the marked differences existing between both T. cruzi Alba groups (TcAlba1/2 and TcAlba30/40), regarding sequence divergence, RNA binding potential, and life-cycle expression patterns, we suggest that they would be involved in different biological processes.

CUCUME: An RNA methylation database integrating systemic mRNAs signals, GWAS and QTL genetic regulation and epigenetics in different tissues of Cucurbitaceae.

As an internal modification of transcripts, RNA methylation determines RNA fate by changing RNA-protein binding affinity. In plants, RNA methylation is ubiquitous and is involved in all aspects of RNA post-transcriptional regulation. For instance, long-distance mobile RNAs, strongly influenced by their methylation status, play important roles in plant growth, development and environmental adaptation. Cucumber/pumpkin heterografts are widely used to improve stress tolerance of cucumber and to study mobile RNA signals due to their strong developed vasculature system. Here, we developed the Cucume (Cucurbit RNA methylation, http://cucume.cn/) database for these two important vegetables, cucumber (Cucumis sativus L.) and pumpkin (Cucurbita moschata) with high productivity worldwide. We identified mRNAs harboring 5-methylcytosine (m5C) and N6-methyladenosine (m6A) sites in pumpkin and cucumber at the whole genome level via Methylated RNA Immunoprecipitation sequencing (MeRIP-seq) of different tissues and the vascular exudates. In addition to RNA methylation sites, the Cucume database includes graft-transmissible systemic mRNAs identified in previous studies using cucumber/pumpkin heterografts. The further integration of cucumber genome-wide association analysis (GWAS) and quantitative trait loci (QTL) allows the study of RNA methylation-related genetic and epigenetic regulation in cucurbits. Therefore, the here developed Cucume database will promote understanding the role of cucurbit RNA methylation in RNA mobility and QTL, ultimately benefitting future breeding of agronomic crop germplasms.

Quantifying RNA structures and interactions with a unified reduced chain representation model.

RNA is a pivotal molecule that plays critical roles in various cellular processes. Quantifying RNA structures and interactions is essential to understanding RNA function and developing RNA-based therapeutics. Using a unified five-bead model and a non-redundant database, this paper investigates the structural features and interactions of five commonly occurring RNA motifs, i.e., double-stranded helices, hairpin loops, internal/bulge loops, multi-branched junctions, and single-stranded terminal tails. Analyzing detailed distributions of RNA local structural features and base-base interactions reveals a preference for helical structures in both local backbone structures and base orientations. The interactions between adjacent bases exhibit motif-specific and sequence-dependent characteristics, reflecting the distinct topological constraints imposed by different loop-helix connection modes and the varying pairing and stacking interactions among different sequences. These findings shed light on the stability of RNA helices, emphasizing their significance in providing dominant base pairing and stacking interactions for RNA structures and stability. The four non-helix motifs encompass unpaired nucleotide loops and exhibit diverse base-base interactions, contributing to the structural diversity observed in RNA. Overall, the complexity of RNA structure arises from the intricate interplay of base-base interactions.

Aberrant splicing in human cancer: An RNA structural code point of view.

Alternative splicing represents an essential process that occurs widely in eukaryotes. In humans, most genes undergo alternative splicing to ensure transcriptome and proteome diversity reflecting their functional complexity. Over the last decade, aberrantly spliced transcripts due to mutations in cis- or trans-acting splicing regulators have been tightly associated with cancer development, largely drawing scientific attention. Although a plethora of single proteins, ribonucleoproteins, complexed RNAs, and short RNA sequences have emerged as nodal contributors to the splicing cascade, the role of RNA secondary structures in warranting splicing fidelity has been underestimated. Recent studies have leveraged the establishment of novel high-throughput methodologies and bioinformatic tools to shed light on an additional layer of splicing regulation in the context of RNA structural elements. This short review focuses on the most recent available data on splicing mechanism regulation on the basis of RNA secondary structure, emphasizing the importance of the complex RNA G-quadruplex structures (rG4s), and other specific RNA motifs identified as splicing silencers or enhancers. Moreover, it intends to provide knowledge on newly established techniques that allow the identification of RNA structural elements and highlight the potential to develop new RNA-oriented therapeutic strategies against cancer.

Enhanced molecular dynamic simulation studies unravel long-range effects caused by sequence variations and partner binding in RNA aptamers.

Intrinsic flexibility and structural modularity are two common features of RNA molecules. Although functionally crucial, RNA plasticity often represents a major complication in high-resolution structural studies. To overcome this problem, RNAs may be rigidified through the complexation with high-affinity partners such as Fab molecules. This approach has been previously used to characterize the DIR2-aptamer. However, possible perturbations induced by the insertion of the Fab binding site on the DIR2-aptamer conformational properties were not investigated. Here, using enhanced molecular dynamics simulations, we compared the dynamics of the DIR2 aptamer holding the Fab binding site with that of the parental sequence. Our results suggest that the L2-loop modification for the Fab recognition leads to a significant increase in local flexibility that also affects the mobility of distant regions. The trajectories provide clear indications of the groups and the interactions mediating the dynamics transfer in DIR2. The effectiveness of our approach in addressing RNA flexibility was further corroborated by showing its ability to reproduce the most important events affecting the NF-κB RNA aptamer upon dissociation from the partner. Therefore, REMD analyses, a rarely adopted technique to unravel the structural/dynamical properties of aptamers, could efficiently complement experimental data guiding the rational design of nucleic acid therapeutics.

Low compositions of human toll-like receptor 7/8-stimulating RNA motifs in the MERS-CoV, SARS-CoV and SARS-CoV-2 genomes imply a substantial ability to evade human innate immunity.

BACKGROUND: The innate immune system especially Toll-like receptor (TLR) 7/8 and the interferon pathway, constitutes an important first line of defense against single-stranded RNA viruses. However, large-scale, systematic comparisons of the TLR 7/8-stimulating potential of genomic RNAs of single-stranded RNA viruses are rare. In this study, a computational method to evaluate the human TLR 7/8-stimulating ability of single-stranded RNA virus genomes based on their human TLR 7/8-stimulating trimer compositions was used to analyze 1,002 human coronavirus genomes. RESULTS: The human TLR 7/8-stimulating potential of coronavirus genomic (positive strand) RNAs followed the order of NL63-CoV > HKU1-CoV >229E-CoV ≅ OC63-CoV > SARS-CoV-2 > MERS-CoV > SARS-CoV. These results suggest that among these coronaviruses, MERS-CoV, SARS-CoV and SARS-CoV-2 may have a higher ability to evade the human TLR 7/8-mediated innate immune response. Analysis with a logistic regression equation derived from human coronavirus data revealed that most of the 1,762 coronavirus genomic (positive strand) RNAs isolated from bats, camels, cats, civets, dogs and birds exhibited weak human TLR 7/8-stimulating potential equivalent to that of the MERS-CoV, SARS-CoV and SARS-CoV-2 genomic RNAs. CONCLUSIONS: Prediction of the human TLR 7/8-stimulating potential of viral genomic RNAs may be useful for surveillance of emerging coronaviruses from nonhuman mammalian hosts.

Modeling and Predicting RNA Three-Dimensional Structures.

Modeling the three-dimensional structure of RNAs is a milestone toward better understanding and prediction of nucleic acids molecular functions. Physics-based approaches and molecular dynamics simulations are not tractable on large molecules with all-atom models. To address this issue, coarse-grained models of RNA three-dimensional structures have been developed. In this chapter, we describe a graphical modeling based on the Leontis-Westhof extended base pair classification. This representation of RNA structures enables us to identify highly conserved structural motifs with complex nucleotide interactions in structure databases. We show how to take advantage of this knowledge to quickly predict three-dimensional structures of large RNA molecules and present the RNA-MoIP web server (http://rnamoip.cs.mcgill.ca) that streamlines the computational and visualization processes. Finally, we show recent advances in the prediction of local 3D motifs from sequence data with the BayesPairing software and discuss its impact toward complete 3D structure prediction.

Engineering Crystal Packing in RNA Structures I: Past and Future Strategies for Engineering RNA Packing in Crystals.

X-ray crystallography remains a powerful method to gain atomistic insights into the catalytic and regulatory functions of RNA molecules. However, the technique requires the preparation of diffraction-quality crystals. This is often a resource- and time-consuming venture because RNA crystallization is hindered by the conformational heterogeneity of RNA, as well as the limited opportunities for stereospecific intermolecular interactions between RNA molecules. The limited success at crystallization explains in part the smaller number of RNA-only structures in the Protein Data Bank. Several approaches have been developed to aid the formation of well-ordered RNA crystals. The majority of these are construct-engineering techniques that aim to introduce crystal contacts to favor the formation of well-diffracting crystals. A typical example is the insertion of tetraloop-tetraloop receptor pairs into non-essential RNA segments to promote intermolecular association. Other methods of promoting crystallization involve chaperones and crystallization-friendly molecules that increase RNA stability and improve crystal packing. In this review, we discuss the various techniques that have been successfully used to facilitate crystal packing of RNA molecules, recent advances in construct engineering, and directions for future research in this vital aspect of RNA crystallography.

Organization of the Influenza A Virus Genomic RNA in the Viral Replication Cycle-Structure, Interactions, and Implications for the Emergence of New Strains.

The influenza A virus is a human pathogen causing respiratory infections. The ability of this virus to trigger seasonal epidemics and sporadic pandemics is a result of its high genetic variability, leading to the ineffectiveness of vaccinations and current therapies. The source of this variability is the accumulation of mutations in viral genes and reassortment enabled by its segmented genome. The latter process can induce major changes and the production of new strains with pandemic potential. However, not all genetic combinations are tolerated and lead to the assembly of complete infectious virions. Reports have shown that viral RNA segments co-segregate in particular circumstances. This tendency is a consequence of the complex and selective genome packaging process, which takes place in the final stages of the viral replication cycle. It has been shown that genome packaging is governed by RNA-RNA interactions. Intersegment contacts create a network, characterized by the presence of common and strain-specific interaction sites. Recent studies have revealed certain RNA regions, and conserved secondary structure motifs within them, which may play functional roles in virion assembly. Growing knowledge on RNA structure and interactions facilitates our understanding of the appearance of new genome variants, and may allow for the prediction of potential reassortment outcomes and the emergence of new strains in the future.

Binding site comparisons for target-centered drug discovery.

INTRODUCTION: The success of binding site comparisons in drug discovery is based on the recognized fact that many different proteins have similar binding sites. Indeed, binding site comparisons have found many uses in drug development and have the potential to dramatically cut the cost and shorten the time necessary for the development of new drugs. Areas covered: The authors review recent methods for comparing protein binding sites and their use in drug repurposing and polypharmacology. They examine emerging fields including the use of binding site comparisons in precision medicine, the prediction of structured water molecules, the search for targets of natural compounds, and their application in the development of protein-based drugs by loop modeling and for comparison of RNA binding sites. Expert opinion: Binding site comparisons have produced many interesting results in drug development, but relatively little work has been done on protein-protein interaction sites, which are particularly relevant in view of the success of biological drugs. Growth of protein loop modeling for modulating biological drugs is anticipated. The fusion of currently distinct methods for the comparison of RNA and protein binding sites into a single comprehensive approach could allow the search for new selective ribosomal antibiotics and initiate pharmaceutical research into other nucleoproteins.

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.