Investigating RNA-RNA interactions through computational and biophysical analysis.

Numerous viruses utilize essential long-range RNA-RNA genome interactions, specifically flaviviruses. Using Japanese encephalitis virus (JEV) as a model system, we computationally predicted and then biophysically validated and characterized its long-range RNA-RNA genomic interaction. Using multiple RNA computation assessment programs, we determine the primary RNA-RNA interacting site among JEV isolates and numerous related viruses. Following in vitro transcription of RNA, we provide, for the first time, characterization of an RNA-RNA interaction using size-exclusion chromatography coupled with multi-angle light scattering and analytical ultracentrifugation. Next, we demonstrate that the 5' and 3' terminal regions of JEV interact with nM affinity using microscale thermophoresis, and this affinity is significantly reduced when the conserved cyclization sequence is not present. Furthermore, we perform computational kinetic analyses validating the cyclization sequence as the primary driver of this RNA-RNA interaction. Finally, we examined the 3D structure of the interaction using small-angle X-ray scattering, revealing a flexible yet stable interaction. This pathway can be adapted and utilized to study various viral and human long-non-coding RNA-RNA interactions and determine their binding affinities, a critical pharmacological property of designing potential therapeutics.

The 3' terminal region of Zika virus RNA contains a conserved G-quadruplex and is unfolded by human DDX17.

Zika virus (ZIKV) infection remains a worldwide concern, and currently no effective treatments or vaccines are available. Novel therapeutics are an avenue of interest that could probe viral RNA-human protein communication to stop viral replication. One specific RNA structure, G-quadruplexes (G4s), possess various roles in viruses and all domains of life, including transcription and translation regulation and genome stability, and serves as nucleation points for RNA liquid-liquid phase separation. Previous G4 studies on ZIKV using a quadruplex forming G-rich sequences Mapper located a potential G-quadruplex sequence in the 3' terminal region (TR) and was validated structurally using a 25-mer oligo. It is currently unknown if this structure is conserved and maintained in a large ZIKV RNA transcript and its specific roles in viral replication. Using bioinformatic analysis and biochemical assays, we demonstrate that the ZIKV 3' TR G4 is conserved across all ZIKV isolates and maintains its structure in a 3' TR full-length transcript. We further established the G4 formation using pyridostatin and the BG4 G4-recognizing antibody binding assays. Our study also demonstrates that the human DEAD-box helicases, DDX3X132-607 and DDX17135-555, bind to the 3' TR and that DDX17135-555 unfolds the G4 present in the 3' TR. These findings provide a path forward in potential therapeutic targeting of DDX3X or DDX17's binding to the 3' TR G4 region for novel treatments against ZIKV.

Xinyang flavivirus, from Haemaphysalis flava ticks in Henan Province, China, defines a basal, likely tick-only Orthoflavivirus clade.

Tick-borne orthoflaviviruses (TBFs) are classified into three conventional groups based on genetics and ecology: mammalian, seabird and probable-TBF group. Recently, a fourth basal group has been identified in Rhipicephalus ticks from Africa: Mpulungu flavivirus (MPFV) in Zambia and Ngoye virus (NGOV) in Senegal. Despite attempts, isolating these viruses in vertebrate and invertebrate cell lines or intracerebral injection of newborn mice with virus-containing homogenates has remained unsuccessful. In this study, we report the discovery of Xinyang flavivirus (XiFV) in Haemaphysalis flava ticks from Xìnyáng, Henan Province, China. Phylogenetic analysis shows that XiFV was most closely related to MPFV and NGOV, marking the first identification of this tick orthoflavivirus group in Asia. We developed a reverse transcriptase quantitative PCR assay to screen wild-collected ticks and egg clutches, with absolute infection rates of 20.75 % in adult females and 15.19 % in egg clutches, suggesting that XiFV could be potentially spread through transovarial transmission. To examine potential host range, dinucleotide composition analyses revealed that XiFV, MPFV and NGOV share a closer composition to classical insect-specific orthoflaviviruses than to vertebrate-infecting TBFs, suggesting that XiFV could be a tick-only orthoflavivirus. Additionally, both XiFV and MPFV lack a furin cleavage site in the prM protein, unlike other TBFs, suggesting these viruses might exist towards a biased immature particle state. To examine this, chimeric Binjari virus with XIFV-prME (bXiFV) was generated, purified and analysed by SDS-PAGE and negative-stain transmission electron microscopy, suggesting prototypical orthoflavivirus size (~50 nm) and bias towards uncleaved prM. In silico structural analyses of the 3'-untranslated regions show that XiFV forms up to five pseudo-knot-containing stem-loops and a prototypical orthoflavivirus dumbbell element, suggesting the potential for multiple exoribonuclease-resistant RNA structures.

Biophysical characterisation of human LincRNA-p21 sense and antisense Alu inverted repeats.

Human Long Intergenic Noncoding RNA-p21 (LincRNA-p21) is a regulatory noncoding RNA that plays an important role in promoting apoptosis. LincRNA-p21 is also critical in down-regulating many p53 target genes through its interaction with a p53 repressive complex. The interaction between LincRNA-p21 and the repressive complex is likely dependent on the RNA tertiary structure. Previous studies have determined the two-dimensional secondary structures of the sense and antisense human LincRNA-p21 AluSx1 IRs using SHAPE. However, there were no insights into its three-dimensional structure. Therefore, we in vitro transcribed the sense and antisense regions of LincRNA-p21 AluSx1 Inverted Repeats (IRs) and performed analytical ultracentrifugation, size exclusion chromatography, light scattering, and small angle X-ray scattering (SAXS) studies. Based on these studies, we determined low-resolution, three-dimensional structures of sense and antisense LincRNA-p21. By adapting previously known two-dimensional information, we calculated their sense and antisense high-resolution models and determined that they agree with the low-resolution structures determined using SAXS. Thus, our integrated approach provides insights into the structure of LincRNA-p21 Alu IRs. Our study also offers a viable pipeline for combining the secondary structure information with biophysical and computational studies to obtain high-resolution atomistic models for long noncoding RNAs.

Strukturierte RNAs in Viren.

Evolutionarily conserved RNAs in untranslated regions are key regulators of the viral life cycle. Exoribonuclease-resistant RNAs (xrRNAs) are particularly interesting examples of structurally conserved elements because they actively dysregulate the messenger RNA (mRNA) degradation machinery of host cells, thereby mediating viral pathogenicity. We review the principles of RNA structure conservation in viruses and discuss potential applications of xrRNAs in synthetic biology and future mRNA vaccines.

Insights into the secondary and tertiary structure of the Bovine Viral Diarrhea Virus Internal Ribosome Entry Site.

The internal ribosome entry site (IRES) RNA of bovine viral diarrhoea virus (BVDV), an economically significant Pestivirus, is required for the cap-independent translation of viral genomic RNA. Thus, it is essential for viral replication and pathogenesis. We applied a combination of high-throughput biochemical RNA structure probing (SHAPE-MaP) and in silico modelling approaches to gain insight into the secondary and tertiary structures of BVDV IRES RNA. Our study demonstrated that BVDV IRES RNA in solution forms a modular architecture composed of three distinct structural domains (I-III). Two regions within domain III are represented in tertiary interactions to form an H-type pseudoknot. Computational modelling of the pseudoknot motif provided a fine-grained picture of the tertiary structure and local arrangement of helices in the BVDV IRES. Furthermore, comparative genomics and consensus structure predictions revealed that the pseudoknot is evolutionarily conserved among many Pestivirus species. These studies provide detailed insight into the structural arrangement of BVDV IRES RNA H-type pseudoknot and encompassing motifs that likely contribute to the optimal functionality of viral cap-independent translation element.

Interplay between the catabolite repression control protein Crc, Hfq and RNA in Hfq-dependent translational regulation in Pseudomonas aeruginosa.

In Pseudomonas aeruginosa the RNA chaperone Hfq and the catabolite repression control protein (Crc) act as post-transcriptional regulators during carbon catabolite repression (CCR). In this regard Crc is required for full-fledged Hfq-mediated translational repression of catabolic genes. RNAseq based transcriptome analyses revealed a significant overlap between the Crc and Hfq regulons, which in conjunction with genetic data supported a concerted action of both proteins. Biochemical and biophysical approaches further suggest that Crc and Hfq form an assembly in the presence of RNAs containing A-rich motifs, and that Crc interacts with both, Hfq and RNA. Through these interactions, Crc enhances the stability of Hfq/Crc/RNA complexes, which can explain its facilitating role in Hfq-mediated translational repression. Hence, these studies revealed for the first time insights into how an interacting protein can modulate Hfq function. Moreover, Crc is shown to interfere with binding of a regulatory RNA to Hfq, which bears implications for riboregulation. These results are discussed in terms of a working model, wherein Crc prioritizes the function of Hfq toward utilization of favored carbon sources.

Indications for a moonlighting function of translation factor aIF5A in the crenarchaeum Sulfolobus solfataricus.

Translation factor a/eIF5A is highly conserved in Eukarya and Archaea. The eukaryal eIF5A protein is required for transit of ribosomes across consecutive proline codons, whereas the function of the archaeal orthologue remains unknown. Here, we provide a first hint for an involvement of Sulfolobus solfataricus (Sso) aIF5A in translation. CRISPR-mediated knock down of the aif5A gene resulted in strong growth retardation, underlining a pivotal function. Moreover, in vitro studies revealed that Sso aIF5A is endowed with endoribonucleolytic activity. Thus, aIF5A appears to be a moonlighting protein that might be involved in protein synthesis as well as in RNA metabolism.

Efficient computation of co-transcriptional RNA-ligand interaction dynamics.

Riboswitches form an abundant class of cis-regulatory RNA elements that mediate gene expression by binding a small metabolite. For synthetic biology applications, they are becoming cheap and accessible systems for selectively triggering transcription or translation of downstream genes. Many riboswitches are kinetically controlled, hence knowledge of their co-transcriptional mechanisms is essential. We present here an efficient implementation for analyzing co-transcriptional RNA-ligand interaction dynamics. This approach allows for the first time to model concentration-dependent metabolite binding/unbinding kinetics. We exemplify this novel approach by means of the recently studied I-A 2'-deoxyguanosine (2'dG)-sensing riboswitch from Mesoplasma florum.

In silico design of ligand triggered RNA switches.

This contribution sketches a work flow to design an RNA switch that is able to adapt two structural conformations in a ligand-dependent way. A well characterized RNA aptamer, i.e., knowing its Kd and adaptive structural features, is an essential ingredient of the described design process. We exemplify the principles using the well-known theophylline aptamer throughout this work. The aptamer in its ligand-binding competent structure represents one structural conformation of the switch while an alternative fold that disrupts the binding-competent structure forms the other conformation. To keep it simple we do not incorporate any regulatory mechanism to control transcription or translation. We elucidate a commonly used design process by explicitly dissecting and explaining the necessary steps in detail. We developed a novel objective function which specifies the mechanistics of this simple, ligand-triggered riboswitch and describe an extensive in silico analysis pipeline to evaluate important kinetic properties of the designed sequences. This protocol and the developed software can be easily extended or adapted to fit novel design scenarios and thus can serve as a template for future needs.

The SmAP1/2 proteins of the crenarchaeon Sulfolobus solfataricus interact with the exosome and stimulate A-rich tailing of transcripts.

The conserved Sm and Sm-like proteins are involved in different aspects of RNA metabolism. Here, we explored the interactome of SmAP1 and SmAP2 of the crenarchaeon Sulfolobus solfataricus (Sso) to shed light on their physiological function(s). Both, SmAP1 and SmAP2 co-purified with several proteins involved in RNA-processing/modification, translation and protein turnover as well as with components of the exosome involved in 3΄ to 5΄ degradation of RNA. In follow-up studies a direct interaction with the poly(A) binding and accessory exosomal subunit DnaG was demonstrated. Moreover, elevated levels of both SmAPs resulted in increased abundance of the soluble exosome fraction, suggesting that they affect the subcellular localization of the exosome in the cell. The increased solubility of the exosome was accompanied by augmented levels of RNAs with A-rich tails that were further characterized using RNASeq. Hence, the observation that the Sso SmAPs impact on the activity of the exosome revealed a hitherto unrecognized function of SmAPs in archaea.

The MazF-regulon: a toolbox for the post-transcriptional stress response in Escherichia coli.

Flexible adaptation to environmental stress is vital for bacteria. An energy-efficient post-transcriptional stress response mechanism in Escherichia coli is governed by the toxin MazF. After stress-induced activation the endoribonuclease MazF processes a distinct subset of transcripts as well as the 16S ribosomal RNA in the context of mature ribosomes. As these 'stress-ribosomes' are specific for the MazF-processed mRNAs, the translational program is changed. To identify this 'MazF-regulon' we employed Poly-seq (polysome fractionation coupled with RNA-seq analysis) and analyzed alterations introduced into the transcriptome and translatome after mazF overexpression. Unexpectedly, our results reveal that the corresponding protein products are involved in all cellular processes and do not particularly contribute to the general stress response. Moreover, our findings suggest that translational reprogramming serves as a fast-track reaction to harsh stress and highlight the so far underestimated significance of selective translation as a global regulatory mechanism in gene expression. Considering the reported implication of toxin-antitoxin (TA) systems in persistence, our results indicate that MazF acts as a prime effector during harsh stress that potentially introduces translational heterogeneity within a bacterial population thereby stimulating persister cell formation.

NMR Structural Profiling of Transcriptional Intermediates Reveals Riboswitch Regulation by Metastable RNA Conformations.

Gene repression induced by the formation of transcriptional terminators represents a prime example for the coupling of RNA synthesis, folding, and regulation. In this context, mapping the changes in available conformational space of transcription intermediates during RNA synthesis is important to understand riboswitch function. A majority of riboswitches, an important class of small metabolite-sensing regulatory RNAs, act as transcriptional regulators, but the dependence of ligand binding and the subsequent allosteric conformational switch on mRNA transcript length has not yet been investigated. We show a strict fine-tuning of binding and sequence-dependent alterations of conformational space by structural analysis of all relevant transcription intermediates at single-nucleotide resolution for the I-A type 2'dG-sensing riboswitch from Mesoplasma florum by NMR spectroscopy. Our results provide a general framework to dissect the coupling of synthesis and folding essential for riboswitch function, revealing the importance of metastable states for RNA-based gene regulation.

SHAPE directed RNA folding.

SUMMARY: Chemical mapping experiments allow for nucleotide resolution assessment of RNA structure. We demonstrate that different strategies of integrating probing data with thermodynamics-based RNA secondary structure prediction algorithms can be implemented by means of soft constraints. This amounts to incorporating suitable pseudo-energies into the standard energy model for RNA secondary structures. As a showcase application for this new feature of the ViennaRNA Package we compare three distinct, previously published strategies to utilize SHAPE reactivities for structure prediction. The new tool is benchmarked on a set of RNAs with known reference structure. AVAILABILITY AND IMPLEMENTATION: The capability for SHAPE directed RNA folding is part of the upcoming release of the ViennaRNA Package 2.2, for which a preliminary release is already freely available at http://www.tbi.univie.ac.at/RNA. CONTACT: michael.wolfinger@univie.ac.at SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Predicting RNA secondary structures from sequence and probing data.

RNA secondary structures have proven essential for understanding the regulatory functions performed by RNA such as microRNAs, bacterial small RNAs, or riboswitches. This success is in part due to the availability of efficient computational methods for predicting RNA secondary structures. Recent advances focus on dealing with the inherent uncertainty of prediction by considering the ensemble of possible structures rather than the single most stable one. Moreover, the advent of high-throughput structural probing has spurred the development of computational methods that incorporate such experimental data as auxiliary information.

Memory-efficient RNA energy landscape exploration.

MOTIVATION: Energy landscapes provide a valuable means for studying the folding dynamics of short RNA molecules in detail by modeling all possible structures and their transitions. Higher abstraction levels based on a macro-state decomposition of the landscape enable the study of larger systems; however, they are still restricted by huge memory requirements of exact approaches. RESULTS: We present a highly parallelizable local enumeration scheme that enables the computation of exact macro-state transition models with highly reduced memory requirements. The approach is evaluated on RNA secondary structure landscapes using a gradient basin definition for macro-states. Furthermore, we demonstrate the need for exact transition models by comparing two barrier-based approaches, and perform a detailed investigation of gradient basins in RNA energy landscapes. AVAILABILITY AND IMPLEMENTATION: Source code is part of the C++ Energy Landscape Library available at http://www.bioinf.uni-freiburg.de/Software/.

TSSAR: TSS annotation regime for dRNA-seq data.

BACKGROUND: Differential RNA sequencing (dRNA-seq) is a high-throughput screening technique designed to examine the architecture of bacterial operons in general and the precise position of transcription start sites (TSS) in particular. Hitherto, dRNA-seq data were analyzed by visualizing the sequencing reads mapped to the reference genome and manually annotating reliable positions. This is very labor intensive and, due to the subjectivity, biased. RESULTS: Here, we present TSSAR, a tool for automated de novo TSS annotation from dRNA-seq data that respects the statistics of dRNA-seq libraries. TSSAR uses the premise that the number of sequencing reads starting at a certain genomic position within a transcriptional active region follows a Poisson distribution with a parameter that depends on the local strength of expression. The differences of two dRNA-seq library counts thus follow a Skellam distribution. This provides a statistical basis to identify significantly enriched primary transcripts.We assessed the performance by analyzing a publicly available dRNA-seq data set using TSSAR and two simple approaches that utilize user-defined score cutoffs. We evaluated the power of reproducing the manual TSS annotation. Furthermore, the same data set was used to reproduce 74 experimentally validated TSS in H. pylori from reliable techniques such as RACE or primer extension. Both analyses showed that TSSAR outperforms the static cutoff-dependent approaches. CONCLUSIONS: Having an automated and efficient tool for analyzing dRNA-seq data facilitates the use of the dRNA-seq technique and promotes its application to more sophisticated analysis. For instance, monitoring the plasticity and dynamics of the transcriptomal architecture triggered by different stimuli and growth conditions becomes possible.The main asset of a novel tool for dRNA-seq analysis that reaches out to a broad user community is usability. As such, we provide TSSAR both as intuitive RESTful Web service ( http://rna.tbi.univie.ac.at/TSSAR) together with a set of post-processing and analysis tools, as well as a stand-alone version for use in high-throughput dRNA-seq data analysis pipelines.

A framework for automated scalable designation of viral pathogen lineages from genomic data.

Pathogen lineage nomenclature systems are a key component of effective communication and collaboration for researchers and public health workers. Since February 2021, the Pango dynamic lineage nomenclature for SARS-CoV-2 has been sustained by crowdsourced lineage proposals as new isolates were sequenced. This approach is vulnerable to time-critical delays as well as regional and personal bias. Here we developed a simple heuristic approach for dividing phylogenetic trees into lineages, including the prioritization of key mutations or genes. Our implementation is efficient on extremely large phylogenetic trees consisting of millions of sequences and produces similar results to existing manually curated lineage designations when applied to SARS-CoV-2 and other viruses including chikungunya virus, Venezuelan equine encephalitis virus complex and Zika virus. This method offers a simple, automated and consistent approach to pathogen nomenclature that can assist researchers in developing and maintaining phylogeny-based classifications in the face of ever-increasing genomic datasets.

BarMap: RNA folding on dynamic energy landscapes.

Dynamical changes of RNA secondary structures play an important role in the function of many regulatory RNAs. Such kinetic effects, especially in time-variable and externally triggered systems, are usually investigated by means of extensive and expensive simulations of large sets of individual folding trajectories. Here we describe the theoretical foundations of a generic approach that not only allows the direct computation of approximate population densities but also reduces the efforts required to analyze the folding energy landscapes to a one-time preprocessing step. The basic idea is to consider the kinetics on individual landscapes and to model external triggers and environmental changes as small but discrete changes in the landscapes. A "barmap" links macrostates of temporally adjacent landscapes and defines the transfer of population densities from one "snapshot" to the next. Implemented in the BarMap software, this approach makes it feasible to study folding processes at the level of basins, saddle points, and barriers for many nonstationary scenarios, including temperature changes, cotranscriptional folding, refolding in consequence to degradation, and mechanically constrained kinetics, as in the case of the translocation of a polymer through a pore.

Evolutionary traits of Tick-borne encephalitis virus: Pervasive non-coding RNA structure conservation and molecular epidemiology.

Tick-borne encephalitis virus (TBEV) is the aetiological agent of tick-borne encephalitis, an infectious disease of the central nervous system that is often associated with severe sequelae in humans. While TBEV is typically classified into three subtypes, recent evidence suggests a more varied range of TBEV subtypes and lineages that differ substantially in the architecture of their 3' untranslated region (3'UTR). Building on comparative genomic approaches and thermodynamic modelling, we characterize the TBEV UTR structureome diversity and propose a unified picture of pervasive non-coding RNA structure conservation. Moreover, we provide an updated phylogeny of TBEV, building on more than 220 publicly available complete genomes, and investigate the molecular epidemiology and phylodynamics with Nextstrain, a web-based visualization framework for real-time pathogen evolution.