Sequence Design for RNA-RNA Interactions.

The design of RNA sequences with desired structural properties presents a challenging computational problem with promising applications in biotechnology and biomedicine. Most regulatory RNAs function by forming RNA-RNA interactions, e.g., in order to regulate mRNA expression. It is therefore natural to consider problems where a sequence is designed to form a desired RNA-RNA interaction and switch between structures upon binding. This contribution demonstrates the use of the Infrared framework to design interacting sequences. Specifically, we consider the regulation of the rpoS mRNA by the sRNA DsrA and design artificial 5 ' UTRs that place a downstream protein coding gene under control of DsrA. The design process is explained step by step in a Jupyter notebook, accompanied by Python code. The text discusses setting up design constraints for sampling sequences in Infrared, computing quality measures, constructing a suitable cost function, as well as the optimization procedure. We show that not only thermodynamic but also kinetic folding features can be relevant. Kinetics of interaction formation can be estimated efficiently using the RRIkinDP tool, and the chapter explains how to include kinetic folding features from RRIkinDP directly in the cost function. The protocol implemented in our Jupyter notebook can easily be extended to consider additional requirements or adapted to novel design scenarios.

Sequence Design Using RNAstructure.

RNA is present in all domains of life. It was once thought to be solely involved in protein expression, but recent advances have revealed its crucial role in catalysis and gene regulation through noncoding RNA. With a growing interest in exploring RNAs with specific structures, there is an increasing focus on designing RNA structures for in vivo and in vitro experimentation and for therapeutics. The development of RNA secondary structure prediction methods has also spurred the growth of RNA design software. However, there are challenges to designing RNA sequences that meet secondary structure requirements. One major challenge is that the secondary structure design problem is likely NP-hard, making it computationally intensive. Another issue is that objective functions need to consider the folding ensemble of RNA molecules to avoid off target structures. In this chapter, we provide protocols for two software tools from the RNAstructure package: "Design" for structured RNA sequence design and "orega" for unstructured RNA sequence design.

SHAPE Probing to Screen Computationally Designed RNA.

RNA design is a major challenge for the future development of synthetic biology and RNA-based therapy. The development of efficient and accurate RNA design pipelines is based on trial and error strategies. The fast progression of such algorithms requires assaying the properties of many RNA sequences in a short time frame. High throughput RNA structure chemical probing technologies such as SHAPE-MaP allow for assaying RNA structure and interaction rapidly and at a very large scale. However, the promiscuity of the designed sequences that may differ only by one nucleotide requires special care. In addition, it necessitates the analysis and evaluation of many experimental results that may reveal to be very tedious. Here we propose an experimental and analytical workflow that eases the screening of thousands of designed RNA sequences at once. In particular, we have developed shapemap tools a customized software suite available at https://github.com/sargueil-citcom/shapemap-tools .

Datasets for Benchmarking RNA Design Algorithms.

RNA molecules play vital roles in many biological processes, such as gene regulation or protein synthesis. The adoption of a specific secondary and tertiary structure by RNA is essential to perform these diverse functions, making RNA a popular tool in bioengineering therapeutics. The field of RNA design responds to the need to develop novel RNA molecules that possess specific functional attributes. In recent years, computational tools for predicting RNA sequences with desired folding characteristics have improved and expanded. However, there is still a lack of well-defined and standardized datasets to assess these programs. Here, we present a large dataset of internal and multibranched loops extracted from PDB-deposited RNA structures that encompass a wide spectrum of design difficulties. Furthermore, we conducted benchmarking tests of widely utilized open-source RNA design algorithms employing this dataset.

The role of structure in regulatory RNA elements.

Regulatory RNA elements fulfill functions such as translational regulation, control of transcript levels, and regulation of viral genome replication. Trans-acting factors (i.e. RNA-binding proteins) bind the so-called cis elements and confer functionality to the complex. The specificity during protein-RNA complex (RNP) formation often exploits the structural plasticity of RNA. Functional integrity of cis-trans pairs depends on the availability of properly folded RNA elements, and RNA conformational transitions can cause diseases. Knowledge of RNA structure and the conformational space is needed for understanding complex formation and deducing functional effects. However, structure determination of RNAs under in vivo conditions remains challenging. This review provides an overview of structured eukaryotic and viral RNA cis elements and discusses the effect of RNA structural equilibria on RNP formation. We showcase implications of RNA structural changes for diseases, outline strategies for RNA structure-based drug targeting, and summarize the methodological toolbox for deciphering RNA structures.

Selective deuteration of an RNA:RNA complex for structural analysis using small-angle scattering.

The structures of RNA:RNA complexes regulate many biological processes. Despite their importance, protein-free RNA:RNA complexes represent a tiny fraction of experimentally-determined structures. Here, we describe a joint small-angle X-ray and neutron scattering (SAXS/SANS) approach to structurally interrogate conformational changes in a model RNA:RNA complex. Using SAXS, we measured the solution structures of the individual RNAs in their free state and of the overall RNA:RNA complex. With SANS, we demonstrate, as a proof-of-principle, that isotope labeling and contrast matching (CM) can be combined to probe the bound state structure of an RNA within a selectively deuterated RNA:RNA complex. Furthermore, we show that experimental scattering data can validate and improve predicted AlphaFold 3 RNA:RNA complex structures to reflect its solution structure. Our work demonstrates that in silico modeling, SAXS, and CM-SANS can be used in concert to directly analyze conformational changes within RNAs when in complex, enhancing our understanding of RNA structure in functional assemblies.

Putting a Kink in HIV-1 Particle Infectivity: Rocaglamide Inhibits HIV-1 Replication by Altering Gag-Genomic RNA Interaction.

Our examination of RNA helicases for effects on HIV-1 protein production and particle assembly identified Rocaglamide (RocA), a known modulator of eIF4A1 function, as an inhibitor of HIV-1 replication in primary CD4+ T cells and three cell systems. HIV-1 attenuation by low-nM RocA doses was associated with reduced viral particle formation without a marked decrease in Gag production. Rather, the co-localization of Gag and HIV-1 genomic RNA (gRNA) assemblies was impaired by RocA treatment in a reversible fashion. Ribonucleoprotein (RNP) immunoprecipitation studies recapitulated the loss of Gag-gRNA assemblies upon RocA treatment. Parallel biophysical studies determined that neither RocA nor eIF4A1 independently affected the ability of Gag to interact with viral RNA, but together, they distorted the structure of the HIV-1 RNP visualized by electron microscopy. Taken together, several lines of evidence indicate that RocA induces stable binding of eIF4A1 onto the viral RNA genome in a manner that interferes with the ordered assembly of Gag along Gag-gRNA assemblies required to generate infectious virions.

Alu insertion-mediated dsRNA structure formation with pre-existing Alu elements as a disease-causing mechanism.

We previously identified a homozygous Alu insertion variant (Alu_Ins) in the 3'-untranslated region (3'-UTR) of SPINK1 as the cause of severe infantile isolated exocrine pancreatic insufficiency. Although we established that Alu_Ins leads to the complete loss of SPINK1 mRNA expression, the precise mechanisms remained elusive. Here, we aimed to elucidate these mechanisms through a hypothesis-driven approach. Initially, we speculated that, owing to its particular location, Alu_Ins could independently disrupt mRNA 3' end formation and/or affect other post-transcriptional processes such as nuclear export and translation. However, employing a 3'-UTR luciferase reporter assay, Alu_Ins was found to result in only an ∼50% reduction in luciferase activity compared to wild type, which is insufficient to account for the severe pancreatic deficiency in the Alu_Ins homozygote. We then postulated that double-stranded RNA (dsRNA) structures formed between Alu elements, an upstream mechanism regulating gene expression, might be responsible. Using RepeatMasker, we identified two Alu elements within SPINK1's third intron, both oriented oppositely to Alu_Ins. Through RNAfold predictions and full-length gene expression assays, we investigated orientation-dependent interactions between these Alu repeats. We provide compelling evidence to link the detrimental effect of Alu_Ins to extensive dsRNA structures formed between Alu_Ins and pre-existing intronic Alu sequences, including the restoration of SPINK1 mRNA expression by aligning all three Alu elements in the same orientation. Given the widespread presence of Alu elements in the human genome and the potential for new Alu insertions at almost any locus, our findings have important implications for detecting and interpreting Alu insertions in disease genes.

Alternative translation initiation produces synaptic organizer proteoforms with distinct localization and functions.

While many mRNAs contain more than one translation initiation site (TIS), the functions of most alternative TISs and their corresponding protein isoforms (proteoforms) remain undetermined. Here, we showed that alternative usage of CUG and AUG TISs in neuronal pentraxin receptor (NPR) mRNA produced two proteoforms, of which the ratio was regulated by RNA secondary structure and neuronal activity. Downstream AUG initiation truncated the N-terminal transmembrane domain and produced a secreted NPR proteoform sufficient in promoting synaptic clustering of AMPA-type glutamate receptors. Mutations that altered the ratio of NPR proteoforms reduced AMPA receptors in parvalbumin-positive interneurons and affected learning behaviors in mice. In addition to NPR, upstream AUU-initiated N-terminal extension of C1q-like synaptic organizers anchored these otherwise secreted factors to the membrane. Together, these results uncovered the plasticity of N-terminal signal sequences regulated by alternative TIS usage as a potentially widespread mechanism in diversifying protein localization and functions.

Translation of HIV-1 unspliced RNA is regulated by 5' untranslated region structure.

HIV-1 must generate infectious virions to spread to new hosts and HIV-1 unspliced RNA (HIV-1 RNA) plays two central roles in this process. HIV-1 RNA serves as an mRNA that is translated to generate proteins essential for particle production and replication, and it is packaged into particles as the viral genome. HIV-1 uses several transcription start sites to generate multiple RNAs that differ by a few nucleotides at the 5' end, including those with one (1G) or three (3G) 5' guanosines. The virus relies on host machinery to translate its RNAs in a cap-dependent manner. Here, we demonstrate that the 5' context of HIV-1 RNA affects the efficiency of translation both in vitro and in cells. Although both RNAs are competent for translation, 3G RNA is translated more efficiently than 1G RNA. The 5' untranslated region (UTR) of 1G and 3G RNAs has previously been shown to fold into distinct structural ensembles. We show that HIV-1 mutants in which the 5' UTR of 1G and 3G RNAs fold into similar structures were translated at similar efficiencies. Thus, the host machinery translates two 99.9% identical HIV-1 RNAs with different efficiencies, and the translation efficiency is regulated by the 5' UTR structure.IMPORTANCEHIV-1 unspliced RNA contains all the viral genetic information and encodes virion structural proteins and enzymes. Thus, the unspliced RNA serves distinct roles as viral genome and translation template, both critical for viral replication. HIV-1 generates two major unspliced RNAs with a 2-nt difference at the 5' end (3G RNA and 1G RNA). The 1G transcript is known to be preferentially packaged over the 3G transcript. Here, we showed that 3G RNA is favorably translated over 1G RNA based on its 5' untranslated region (UTR) RNA structure. In HIV-1 mutants in which the two major transcripts have similar 5' UTR structures, 1G and 3G RNAs are translated similarly. Therefore, HIV-1 generates two 9-kb RNAs with a 2-nt difference, each serving a distinct role dictated by differential 5' UTR structures.

Evaluation of the effect of RNA secondary structure on Cas13d-mediated target RNA cleavage.

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas13d system was adapted as a powerful tool for targeting viral RNA sequences, making it a promising approach for antiviral strategies. Understanding the influence of template RNA structure on Cas13d binding and cleavage efficiency is crucial for optimizing its therapeutic potential. In this study, we investigated the effect of local RNA secondary structure on Cas13d activity. To do so, we varied the stability of a hairpin structure containing the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) target sequence, allowing us to determine the threshold RNA stability at which Cas13d activity is affected. Our results demonstrate that Cas13d possesses the ability to effectively bind and cleave highly stable RNA structures. Notably, we only observed a decrease in Cas13d activity in the case of exceptionally stable RNA hairpins with completely base-paired stems, which are rarely encountered in natural RNA molecules. A comparison of Cas13d and RNA interference (RNAi)-mediated cleavage of the same RNA targets demonstrated that RNAi is more sensitive for local target RNA structures than Cas13d. These results underscore the suitability of the CRISPR-Cas13d system for targeting viruses with highly structured RNA genomes.

Atomistic Simulations Reveal Crucial Role of Metal Ions for Ligand Binding in Guanidine-I Riboswitch.

Riboswitches are structured ribonucleic acid (RNA) segments that act as specific sensors for small molecules in bacterial metabolism. Due to the flexible nature of these highly charged macromolecules, molecular dynamics simulations are instrumental to investigating the mechanistic details of their regulatory function. In the present study, the guanidine-I riboswitch serves as an example of how atomistic simulations can shed light on the effect of ions on the structure and dynamics of RNA and on ligand binding. Relying on two orthologous crystal structures from different bacterial species, it is demonstrated how the ion setup crucially determines whether the simulation yields meaningful insights into the conformational stability of the RNA, functionally relevant residues and RNA-ligand interactions. The ion setup in this context includes diffuse ions in solution and bound ions associated directly with the RNA, in particular a triad of 2 Mg2+ ions and a K+ ion in close proximity to the guanidinium binding site. A detailed investigation of the binding pocket reveals that the K+ from the ion triad plays a decisive role in stabilizing the ligand binding by stabilizing important localized interactions, which in turn contribute to the overall shape of the folded state of the RNA.

Phase separation and viral factories: unveiling the physical processes supporting RNA packaging in dsRNA viruses.

Understanding of the physicochemical properties and functions of biomolecular condensates has rapidly advanced over the past decade. More recently, many RNA viruses have been shown to form cytoplasmic replication factories, or viroplasms, via phase separation of their components, akin to numerous cellular membraneless organelles. Notably, diverse viruses from the Reoviridae family containing 10-12 segmented double-stranded RNA genomes induce the formation of viroplasms in infected cells. Little is known about the inner workings of these membraneless cytoplasmic inclusions and how they may support stoichiometric RNA assembly in viruses with segmented RNA genomes, raising questions about the roles of phase separation in coordinating viral genome packaging. Here, we discuss how the molecular composition of viroplasms determines their properties, highlighting the interplay between RNA structure, RNA remodelling, and condensate self-organisation. Advancements in RNA structural probing and theoretical modelling of condensates can reveal the mechanisms through which these ribonucleoprotein complexes support the selective enrichment and stoichiometric assembly of distinct viral RNAs.

Synthetic antibodies for accelerated RNA crystallography.

RNA structure is crucial to a wide range of cellular processes. The intimate relationship between macromolecular structure and function necessitates the determination of high-resolution structures of functional RNA molecules. X-ray crystallography is the predominant technique used for macromolecular structure determination; however, solving RNA structures has been more challenging than their protein counterparts, as reflected in their poor representation in the Protein Data Bank (<1%). Antibody-assisted RNA crystallography is a relatively new technique that promises to accelerate RNA structure determination by employing synthetic antibodies (Fabs) as crystallization chaperones that are specifically raised against target RNAs. Antibody chaperones facilitate the formation of ordered crystal lattices by minimizing RNA flexibility and replacing unfavorable RNA-RNA contacts with contacts between chaperone molecules. Atomic coordinates of these antibody fragments can also be used as search models to obtain phase information during structure determination. Antibody-assisted RNA crystallography has enabled the structure determination of 15 unique RNA targets, including 11 in the last 6 years. In this review, I cover the historical development of antibody fragments as crystallization chaperones and their application to diverse RNA targets. I discuss how the first structures of antibody-RNA complexes informed the design of second-generation antibodies and led to the development of portable crystallization modules that have greatly reduced the uncertainties associated with RNA crystallography. Finally, I outline unexplored avenues that can increase the impact of this technology in structural biology research and discuss potential applications of antibodies as affinity reagents for interrogating RNA biology outside of their use in crystallography. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Structure of Essential RNA Regulatory Elements in the West Nile Virus 3'-Terminal Stem Loop.

Flaviviruses, such as West Nile and Dengue Virus, pose a significant and growing threat to global health. Central to the flavivirus life cycle are highly structured 5'- and 3'-untranslated regions (UTRs), which harbor conserved cis-acting RNA elements critical for viral replication and host adaptation. Despite their essential roles, detailed molecular insights into these RNA elements have been limited. By employing nuclear magnetic resonance (NMR) spectroscopy in conjunction with SAXS experiments, we determined the three-dimensional structure of the West Nile Virus (WNV) 3'-terminal stem-loop core, a highly conserved element critical for viral genome cyclization and replication. Single nucleotide mutations at several sites within this RNA abolish the ability of the virus to replicate. These critical sites are located within a short 18-nucleotide hairpin stem, a substructure notable for its conformational flexibility, while the adjoining main stem-loop adopts a well-defined extended helix interrupted by three non-Watson-Crick pairs. This study enhances our understanding of several metastable RNA structures that play key roles in regulating the flavivirus lifecycle, and thereby also opens up potential new avenues for the development of antivirals targeting these conserved RNA structures. In particular, the structure we observe suggests that the plastic junction between the small hairpin and the tail of the longer stem-loop could provide a binding pocket for small molecules, for example potentially stabilizing the RNA in a conformation which hinders the conformational rearrangements critical for viral replication.

Introns with branchpoint-distant 3' splice sites: Splicing mechanism and regulatory roles.

The two transesterification reactions of pre-mRNA splicing require highly complex yet well-controlled rearrangements of small nuclear RNAs and proteins (snRNP) in the spliceosome. The efficiency and accuracy of these reactions are critical for gene expression, as almost all human genes pass through pre-mRNA splicing. Key parameters that determine the splicing outcome are the length of the intron, the strengths of its splicing signals and gaps between them, and the presence of splicing controlling elements. In particular, the gap between the branchpoint (BP) and the 3' splice site (ss) of introns is a major determinant of the splicing efficiency. This distance falls within a small range across the introns of an organism. The constraints exist possibly because BP and 3'ss are recognized by BP-binding proteins, U2 snRNP and U2 accessory factors (U2AF) in a coordinated manner. Furthermore, varying distances between the two signals may also affect the second transesterification reaction since the intervening RNA needs to be accurately positioned within the complex RNP machinery. Splicing such pre-mRNAs requires cis-acting elements in the RNA and many trans-acting splicing regulators. Regulated pre-mRNA splicing with BP-distant 3'ss adds another layer of control to gene expression and promotes alternative splicing.

Computational exploration of protein structure dynamics and RNA structural consequences of PKD1 missense variants: implications in ADPKD pathogenesis.

We analyzed the impact of nine previously identified missense PKD1 variants from our studies, including c.6928G > A p.G2310R, c.8809G > A p.E2937K, c.2899 T > C p.W967R, c.6284A > G p.D2095G, c.6644G > A p.R2215Q, c.7810G > A p.D2604N, c.11249G > C p.R3750P, c.1001C > T p.T334M, and c.3101A > G p.N1034S on RNA structures and PC1 protein structure dynamics utilizing computational tools. RNA structure analysis was done using short RNA snippets of 41 nucleotides with the variant position at the 21st nucleotide, ensuring 20 bases on both sides. The secondary structures of these RNA snippets were predicted using RNAstructure. Structural changes of the mutants compared to the wild type were analyzed using the MutaRNA webserver. Molecular dynamics (MD) simulation of PC1 wild-type and mutant protein regions were performed using GROMACS 2018 (GROMOS96 54a7 force field). Findings revealed that five variants including c.8809G > A (p.E2937K), c.11249G > C (p.R3750P), c.3101A > G (p.N1034S), c.6928G > A (p.G2310R), c.6644G > A (p.R2215Q) exhibited major alterations in RNA structures and thereby their interactions with other proteins or RNAs affecting protein structure dynamics. While certain variants have minimal impact on RNA conformations, their observed alterations in MD simulations indicate impact on protein structure dynamics highlighting the importance of evaluating the functional consequences of genetic variants by considering both RNA and protein levels. The study also emphasizes that each missense variant exerts a unique impact on RNA stability, and protein structure dynamics, potentially contributing to the heterogeneous clinical manifestations and progression observed in Autosomal Dominant Polycystic Kidney Disease (ADPKD) patients offering a novel perspective in this direction. Thus, the utility of studying the structure dynamics through computational tools can help in prioritizing the variants for their functional implications, understanding the molecular mechanisms underlying variability in ADPKD presentation and developing targeted therapeutic interventions. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-024-04057-9.

Rho-dependent transcriptional switches regulate the bacterial response to cold shock.

Binding of the bacterial Rho helicase to nascent transcripts triggers Rho-dependent transcription termination (RDTT) in response to cellular signals that modulate mRNA structure and accessibility of Rho utilization (Rut) sites. Despite the impact of temperature on RNA structure, RDTT was never linked to the bacterial response to temperature shifts. We show that Rho is a central player in the cold-shock response (CSR), challenging the current view that CSR is primarily a posttranscriptional program. We identify Rut sites in 5'-untranslated regions of key CSR genes/operons (cspA, cspB, cspG, and nsrR-rnr-yjfHI) that trigger premature RDTT at 37°C but not at 15°C. High concentrations of RNA chaperone CspA or nucleotide changes in the cspA mRNA leader reduce RDTT efficiency, revealing how RNA restructuring directs Rho to activate CSR genes during the cold shock and to silence them during cold acclimation. These findings establish a paradigm for how RNA thermosensors can modulate gene expression.

RNA thermometers are widespread upstream of ABC transporter genes in bacteria.

RNA thermometers are temperature-sensing non-coding RNAs that regulate the expression of downstream genes. A well-characterized RNA thermometer motif discovered in bacteria is the ROSE-like element (repression of heat shock gene expression). ATP-binding cassette (ABC) transporters are a superfamily of transmembrane proteins that harness ATP hydrolysis to facilitate the export and import of substrates across cellular membranes. Through structure-guided bioinformatics, we discovered that ROSE-like RNA thermometers are widespread upstream of ABC transporter genes in bacteria. X-ray crystallography, biochemistry, and cellular assays indicate that these RNA thermometers are functional regulatory elements. This study expands the known biological role of RNA thermometers to these key membrane transporters.

Structural Analysis of the Hepatitis B Virus RNA Encapsidation Signal ε by Selective 2'-Hydroxyl Acylation Analyzed by Primer Extension (SHAPE).

RNA structure is crucial for RNA function, including in viral cis-elements such as the hepatitis B virus (HBV) RNA encapsidation signal ε. Interacting with the viral polymerase ε mediates packaging of the pregenomic (pg) RNA into capsids, initiation of reverse transcription, and it affects the mRNA functions of pgRNA. As free RNA, the 61-nucleotide (nt) ε sequence adopts a bipartite stem-loop structure with a central bulge and an apical loop. Due to stable Watson-Crick base pairing, this was already predicted by early RNA folding programs and confirmed by classical enzymatic and chemical structure probing. A newer, high-resolution probing technique exploits the selective acylation of solvent-accessible 2'-hydroxyls in the RNA backbone by electrophilic compounds such as 2-methylnicotinic acid imidazolide (NAI), followed by mapping of the modified sites by primer extension. This SHAPE principle has meanwhile been extended to numerous applications. Here we provide a basic protocol for NAI-based SHAPE of isolated HBV ε RNA which already provided insights into the impact of mutations, and preliminarily, of polymerase binding on the RNA structural dynamics. While the focus is on NAI modification, we also briefly cover target RNA preparation by in vitro transcription, primer extension using a radiolabeled primer, and analysis of the resulting cDNAs by denaturing polyacrylamide gelelectrophoresis (PAGE). Given the high tolerance of SHAPE chemistry to different conditions, including applicability in live cells, we expect this technique to greatly facilitate deciphering the conformational dynamics underlying the various functions of the ε element, especially in concert with the recently solved three-dimensional structure of the free RNA.