On the importance of cotranscriptional RNA structure formation

IoT Platform for COVID-19 Prevention and Control: A Survey

As a result of the worldwide transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), coronavirus disease 2019 (COVID-19) has evolved into an unprecedented pandemic. Currently, with unavailable pharmaceutical treatments and low vaccination rates, this novel coronavirus results in a great impact on public health, human society, and global economy, which is likely to last for many years. One of the lessons learned from the COVID-19 pandemic is that a long-term system with non-pharmaceutical interventions for preventing and controlling new infectious diseases is desirable to be implemented. Internet of things (IoT) platform is preferred to be utilized to achieve this goal, due to its ubiquitous sensing ability and seamless connectivity. IoT technology is changing our lives through smart healthcare, smart home, and smart city, which aims to build a more convenient and intelligent community. This paper presents how the IoT could be incorporated into the epidemic prevention and control system. Specifically, we demonstrate a potential fog-cloud combined IoT platform that can be used in the systematic and intelligent COVID-19 prevention and control, which involves five interventions including COVID-19 Symptom Diagnosis, Quarantine Monitoring, Contact Tracing & Social Distancing, COVID-19 Outbreak Forecasting, and SARS-CoV-2 Mutation Tracking. We investigate and review the state-of-the-art literatures of these five interventions to present the capabilities of IoT in countering against the current COVID-19 pandemic or future infectious disease epidemics.

CoBold: a method for identifying different functional classes of transient RNA structure features that can impact RNA structure formation in vivo

RNA structure formation in vivo happens co-transcriptionally while the transcript is being made. The corresponding co-transcriptional folding pathway typically involves transient RNA structure features that are not part of the final, functional RNA structure. These transient features can play important functional roles of their own and also influence the formation of the final RNA structure in vivo. We here present CoBold, a computational method for identifying different functional classes of transient RNA structure features that can either aid or hinder the formation of a known reference RNA structure. Our method takes as input either a single RNA or a corresponding multiple-sequence alignment as well as a known reference RNA secondary structure and identifies different classes of transient RNA structure features that could aid or prevent the formation of the given RNA structure. We make CoBold available via a web-server which includes dedicated data visualisation.

Protocol for efficient fluorescence 3' end-labeling of native noncoding RNA domains

Noncoding RNAs (ncRNAs) comprise a class of versatile transcripts that are highly involved in the regulation of a wide range of biological processes. Functional long ncRNAs (> 200 nts in length) often adopt secondary structures that arise co-transcriptionally. To maintain the secondary structure elements as well as preparation homogeneity of such transcripts, native-like conditions should be maintained throughout the in vitro synthesis, purification and chemical tagging processes. In this optimized protocol, we describe a simple method for obtaining homogenous samples followed by chemically tagging the 3' termini of natively-purified structured ncRNA domains that are longer than 200 nts. This protocol replaces traditional hazardous radioactive labeling with fluorescence tagging, and eliminates laborious and time consuming RNA purification and concentration steps and replaces them with straightforward recovery of RNA through centrifugal filtration, preserving the homogeneity and mono-dispersion of the preparations. The protocol provides:•An integrative, simple and straightforward approach for synthesis, purification and labeling of structured ncRNAs whilst maintaining their secondary structure intact.•Replacing hazardous, laborious and time-consuming radioactive labeling of RNA with much simpler fluorescence tagging, thereby facilitating potential downstream applications such as electrophoretic mobility shift assay (EMSA).•A versatile protocol that could be applicable to a wide-range of chemical tags and in principle could be used to label DNA or RNA.

STARRPeaker: uniform processing and accurate identification of STARR-seq active regions

STARR-seq technology has employed progressively more complex genomic libraries and increased sequencing depths. An issue with the increased complexity and depth is that the coverage in STARR-seq experiments is non-uniform, overdispersed, and often confounded by sequencing biases, such as GC content. Furthermore, STARR-seq readout is confounded by RNA secondary structure and thermodynamic stability. To address these potential confounders, we developed a negative binomial regression framework for uniformly processing STARR-seq data, called STARRPeaker. Moreover, to aid our effort, we generated whole-genome STARR-seq data from the HepG2 and K562 human cell lines and applied STARRPeaker to comprehensively and unbiasedly call enhancers in them.

The International Virus Bioinformatics Meeting 2020

The International Virus Bioinformatics Meeting 2020 was originally planned to take place in Bern, Switzerland, in March 2020. However, the COVID-19 pandemic put a spoke in the wheel of almost all conferences to be held in 2020. After moving the conference to 8-9 October 2020, we got hit by the second wave and finally decided at short notice to go fully online. On the other hand, the pandemic has made us even more aware of the importance of accelerating research in viral bioinformatics. Advances in bioinformatics have led to improved approaches to investigate viral infections and outbreaks. The International Virus Bioinformatics Meeting 2020 has attracted approximately 120 experts in virology and bioinformatics from all over the world to join the two-day virtual meeting. Despite concerns being raised that virtual meetings lack possibilities for face-to-face discussion, the participants from this small community created a highly interactive scientific environment, engaging in lively and inspiring discussions and suggesting new research directions and questions. The meeting featured five invited and twelve contributed talks, on the four main topics: (1) proteome and RNAome of RNA viruses, (2) viral metagenomics and ecology, (3) virus evolution and classification and (4) viral infections and immunology. Further, the meeting featured 20 oral poster presentations, all of which focused on specific areas of virus bioinformatics. This report summarizes the main research findings and highlights presented at the meeting.

So close, no matter how far: multiple paths connecting transcription to mRNA translation in eukaryotes

Transcription of DNA into mRNA and translation of mRNA into proteins are two major processes underlying gene expression. Due to the distinct molecular mechanisms, timings, and locales of action, these processes are mainly considered to be independent. During the last two decades, however, multiple factors and elements were shown to coordinate transcription and translation, suggesting an intricate level of synchronization. This review discusses the molecular mechanisms that impact both processes in eukaryotic cells of different origins. The emerging global picture suggests evolutionarily conserved regulation and coordination between transcription and mRNA translation, indicating the importance of this phenomenon for the fine-tuning of gene expression and the adjustment to constantly changing conditions.

Challenges and perspectives for structural biology of lncRNAs-the example of the Xist lncRNA A-repeats

Following the discovery of numerous long non-coding RNA (lncRNA) transcripts in the human genome, their important roles in biology and human disease are emerging. Recent progress in experimental methods has enabled the identification of structural features of lncRNAs. However, determining high-resolution structures is challenging as lncRNAs are expected to be dynamic and adopt multiple conformations, which may be modulated by interaction with protein binding partners. The X-inactive specific transcript (Xist) is necessary for X inactivation during dosage compensation in female placental mammals and one of the best-studied lncRNAs. Recent progress has provided new insights into the domain organization, molecular features, and RNA binding proteins that interact with distinct regions of Xist. The A-repeats located at the 5′ end of the transcript are of particular interest as they are essential for mediating silencing of the inactive X chromosome. Here, we discuss recent progress with elucidating structural features of the Xist lncRNA, focusing on the A-repeats. We discuss the experimental and computational approaches employed that have led to distinct structural models, likely reflecting the intrinsic dynamics of this RNA. The presence of multiple dynamic conformations may also play an important role in the formation of the associated RNPs, thus influencing the molecular mechanism underlying the biological function of the Xist A-repeats. We propose that integrative approaches that combine biochemical experiments and high-resolution structural biology in vitro with chemical probing and functional studies in vivo are required to unravel the molecular mechanisms of lncRNAs.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Two sides of the same medal: Noncoding mutations reveal new pathological mechanisms and insights into the regulation of gene expression

Noncoding sequences constitute the major part of the human genome and also of pre-mRNAs. Single nucleotide variants in these regions are often overlooked, but may be responsible for much of the variation of phenotypes observed. Mutations in the noncoding part of pre-mRNAs often reveal new and meaningful insights into the regulation of cellular gene expression. Thus, the mechanistic analysis of the pathological mechanism of such mutations will both foster a deeper understanding of the disease and the underlying cellular pathways. Even synonymous mutations can cause diseases, since the primary mRNA sequence not only encodes amino acids, but also encrypts information on RNA-binding proteins and secondary structure. In fact, the RNA sequence directs assembly of a specific mRNP complex, which in turn dictates the fate of the mRNA or regulates its biogenesis. The accumulation of genomic sequence information is increasing at a rapid pace. However, much of the diversity uncovered may not explain the phenotype of a certain syndrome or disease. For this reason, we also emphasize the value of mechanistic studies on pathological mechanisms being complementary to genome-wide studies and bioinformatic approaches. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease.

Alternate RNA Structures

RNA molecules fold into complex three-dimensional structures that sample alternate conformations ranging from minor differences in tertiary structure dynamics to major differences in secondary structure. This allows them to form entirely different substructures with each population potentially giving rise to a distinct biological outcome. The substructures can be partitioned along an existing energy landscape given a particular static cellular cue or can be shifted in response to dynamic cues such as ligand binding. We review a few key examples of RNA molecules that sample alternate conformations and how these are capitalized on for control of critical regulatory functions.

Base-pair ambiguity and the kinetics of RNA folding

Background

A pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between “bound” molecules, which are active as part of ribonucleoprotein (RNP) complexes, and “unbound,” with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an “ambiguity index”—a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative (“gold standard”) and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules.

Results

We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule—an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure.

Conclusions

A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.

Self-powered RNA nanomachine driven by metastable structure

Many non-coding and regulatory RNA elements have evolved to exploit transient or metastable structures that emerge during transcription to control complex folding pathways or to encode dynamic functions. However, efforts to engineer synthetic RNA devices have mostly focused on the thermodynamically stable structures. Consequently, significant challenges and opportunities exist in engineering functional RNAs that explicitly take advantage of cotranscriptionally generated transient or metastable structures. In this work, we designed a short RNA sequence that adopts a robust metastable structure when transcribed by an RNA polymerase. Although the metastable structure persists for hours at low temperature, it refolds almost completely into the thermodynamically stable structure upon heat denaturation followed by cooling. The synthetic RNA was also equipped with the Broccoli aptamer so that it can bind its ligand and become fluorescent only in the thermodynamically stable structure. We further demonstrated that the relaxation to the thermodynamically stable and fluorescent structure can be catalyzed by a short trigger RNA in a sequence-specific manner. Finally, the RNA architecture was redesigned to sense and respond to microRNA sequences. In summary, we designed RNA nanomachines that can detect an RNA sequence, amplify signal and produce an optical output, all encoded in a single RNA transcript, self-powered by a metastable structure.

A Modified In Vitro Transcription Approach to Improve RNA Synthesis and Ribozyme Cleavage Efficiency

RNA elements such as catalytic RNA, riboswitch, microRNA, and long non-coding RNA perform a major role in cellular processes. A complete understanding of cellular processes is impossible without knowing the structure-function relationship of participating RNA molecules that ultimately requires large quantities of pure RNAs. Thus, structural/functional analyses of emerging RNAs necessitate revised protocols for improved RNA quantity and quality. Here we present a modified in vitro transcription protocol to enhance ribozyme cleaving efficiency and RNA yield by working on two variables, i.e., incubation temperature and limiting GTPs. Following an improved RNA synthesis, the target RNA is purified from transcription mixture components through denaturing size-exclusion chromatography. The protocol confirms that cyclic elevated incubation temperatures during transcription and increased concentrations of GTPs improve the production rate of RNA. Our modified in vitro transcription method improves the ribozyme cleaving efficiency and targets RNA yield by four- to fivefold that can benefit almost any RNA-related study from protein-RNA interaction analysis to crystallography.

A ligand-gated strand displacement mechanism for ZTP riboswitch transcription control

Cotranscriptional folding is an obligate step of RNA biogenesis that can guide RNA structure formation and function through transient intermediate folds. This process is particularly important for transcriptional riboswitches in which the formation of ligand-dependent structures during transcription regulates downstream gene expression. However, the intermediate structures that comprise cotranscriptional RNA folding pathways, and the mechanisms that enable transit between them, remain largely unknown. Here, we determine the series of cotranscriptional folds and rearrangements that mediate antitermination by the Clostridium beijerinckii pfl ZTP riboswitch in response to the purine biosynthetic intermediate ZMP. We uncover sequence and structural determinants that modulate an internal RNA strand displacement process and identify biases within natural ZTP riboswitch sequences that promote on-pathway folding. Our findings establish a mechanism for pfl riboswitch antitermination and suggest general strategies by which nascent RNA molecules navigate cotranscriptional folding pathways.

Synthetic bionanotechnology: synthetic biology finds a toehold in nanotechnology

Enabled by its central role in the molecular networks that govern cell function, RNA has been widely used for constructing components used in biological circuits for synthetic biology. Nucleic acid nanotechnology, which exploits predictable nucleic acid interactions to implement programmable molecular systems, has seen remarkable advances in in vitro nanoscale self-assembly and molecular computation, enabling the production of complex nanostructures and DNA-based neural networks. Living cells genetically engineered to execute nucleic acid nanotechnology programs thus have outstanding potential to significantly extend the current limits of synthetic biology. This perspective discusses the recent developments and future challenges in the field of synthetic bionanotechnology. Thus far, researchers in this emerging area have implemented dozens of programmable RNA nanodevices that provide precise control over gene expression at the transcriptional and translational levels and through CRISPR/Cas effectors. Moreover, they have employed synthetic self-assembling RNA networks in engineered bacteria to carry out computations featuring up to a dozen inputs and to substantially enhance the rate of chemical synthesis. Continued advancement of the field will benefit from improved in vivo strategies for streamlining nucleic acid network synthesis and new approaches for enhancing network function. As the field matures and the complexity gap between in vitro and in vivo systems narrows, synthetic bionanotechnology promises to have diverse potential applications ranging from intracellular circuits that detect and treat disease to synthetic enzymatic pathways that efficiently produce novel drug molecules.

Ruleset Optimization on Isomorphic Oritatami Systems

We study an optimization problem of a computational folding model, proving its hardness and proposing heuristic algorithms. RNA cotranscriptional folding refers to the phenomenon in which an RNA transcript folds upon itself while being synthesized out of a gene. An oritatami model (OM) is a computational model of this phenomenon that lets its sequence of beads (abstract molecules) fold cotranscriptionally by the interactions between beads, according to its ruleset. We study the problem of reducing the ruleset size, while keeping the terminal conformations geometrically the same. We first prove the hardness of finding the smallest ruleset, and then suggest two approaches that reduce the ruleset size efficiently.

A genetic selection reveals functional metastable structures embedded in a toxin-encoding mRNA

Post-transcriptional regulation plays important roles to fine-tune gene expression in bacteria. In particular, regulation of type I toxin-antitoxin (TA) systems is achieved through sophisticated mechanisms involving toxin mRNA folding. Here, we set up a genetic approach to decipher the molecular underpinnings behind the regulation of a type I TA in Helicobacter pylori. We used the lethality induced by chromosomal inactivation of the antitoxin to select mutations that suppress toxicity. We found that single point mutations are sufficient to allow cell survival. Mutations located either in the 5’ untranslated region or within the open reading frame of the toxin hamper its translation by stabilizing stem-loop structures that sequester the Shine-Dalgarno sequence. We propose that these short hairpins correspond to metastable structures that are transiently formed during transcription to avoid premature toxin expression. This work uncovers the co-transcriptional inhibition of translation as an additional layer of TA regulation in bacteria.

Viroid research and its significance for RNA technology and basic biochemistry

Viroids were described 47 years ago as the smallest RNA molecules capable of infecting plants and autonomously self-replicating without an encoded protein. Work on viroids initiated the development of a number of innovative methods. Novel chromatographic and gelelectrophoretic methods were developed for the purification and characterization of viroids; these methods were later used in molecular biology, gene technology and in prion research. Theoretical and experimental studies of RNA folding demonstrated the general biological importance of metastable structures, and nuclear magnetic resonance spectroscopy of viroid RNA showed the partially covalent nature of hydrogen bonds in biological macromolecules. RNA biochemistry and molecular biology profited from viroid research, such as in the detection of RNA as template of DNA-dependent polymerases and in mechanisms of gene silencing. Viroids, the first circular RNA detected in nature, are important for studies on the much wider spectrum of circular RNAs and other non-coding RNAs.

Effects of Refolding on Large-Scale RNA Structure

Understanding RNA-mediated functions requires a detailed characterization of the underlying RNA structure. In many cases, structure probing experiments are performed on RNA that has been "refolded" in some way, which may cause the conformation to differ from that of the native RNA. We used SHAPE-MaP (selective 2'-hydroxyl acylation analyzed by primer extension, read out by mutational profiling), to probe the structure of the Dengue virus (DENV) RNA genome after gentle extraction of the native RNA from intact virions (ex virion) and after heat denaturation and refolding. Comparison of mutiple SHAPE-informed structural features revealed that refolded RNA is more highly structured and samples fewer conformations than the ex virion RNA. Regions with similar structural features are generally those with low SHAPE reactivity and low Shannon entropy (lowSS regions), which correspond to elements with high levels of well-determined structure. This high-structure and low-entropy analysis framework, previously shown to make possible discovery of functional RNA structures, is thus now shown to allow de novo identification of structural elements in a refolded RNA that are likely to recapitulate RNA structures in the ex virion RNA state. Regions with less well-defined structures, which occurred more frequently in the more native-like ex virion RNA and may contain RNA switches, are challenging to recapitulate using refolded RNA.

Conformational Dynamics of the RNA G-Quadruplex and its Effect on Translation Efficiency

During translation, intracellular mRNA folds co-transcriptionally and must refold following the passage of ribosome. The mRNAs can be entrapped in metastable structures during these folding events. In the present study, we evaluated the conformational dynamics of the kinetically favored, metastable, and hairpin-like structure, which disturbs the thermodynamically favored G-quadruplex structure, and its effect on co-transcriptional translation in prokaryotic cells. We found that nascent mRNA forms a metastable hairpin-like structure during co-transcriptional folding instead of the G-quadruplex structure. When the translation progressed co-transcriptionally before the metastable hairpin-like structure transition to the G-quadruplex, function of the G-quadruplex as a roadblock of the ribosome was sequestered. This suggested that kinetically formed RNA structures had a dominant effect on gene expression in prokaryotes. The results of this study indicate that it is critical to consider the conformational dynamics of RNA-folding to understand the contributions of the mRNA structures in controlling gene expression.

Proteins That Chaperone RNA Regulation

RNA-binding proteins chaperone the biological functions of noncoding RNA by reducing RNA misfolding, improving matchmaking between regulatory RNA and targets, and exerting quality control over RNP biogenesis. Recent studies of Escherichia coli CspA, HIV NCp, and E. coli Hfq are beginning to show how RNA-binding proteins remodel RNA structures. These different protein families use common strategies for disrupting or annealing RNA double helices, which can be used to understand the mechanisms by which proteins chaperone RNA-dependent regulation in bacteria.

Kinetics coming into focus: single-molecule microscopy of riboswitch dynamics

Riboswitches are dynamic RNA motifs that are mostly embedded in the 5'-untranslated regions of bacterial mRNAs, where they regulate gene expression transcriptionally or translationally by undergoing conformational changes upon binding of a small metabolite or ion. Due to the small size of typical ligands, relatively little free energy is available from ligand binding to overcome the often high energetic barrier of reshaping RNA structure. Instead, most riboswitches appear to take advantage of the directional and hierarchical folding of RNA by employing the ligand as a structural 'linchpin' to adjust the kinetic partitioning between alternate folds. In this model, even small, local structural and kinetic effects of ligand binding can cascade into global RNA conformational changes affecting gene expression. Single-molecule (SM) microscopy tools are uniquely suited to study such kinetically controlled RNA folding since they avoid the ensemble averaging of bulk techniques that loses sight of unsynchronized, transient, and/or multi-state kinetic behavior. This review summarizes how SM methods have begun to unravel riboswitch-mediated gene regulation.

Motif independent identification of potential RNA G-quadruplexes by G4RNA screener

Motivation

G-quadruplex structures in RNA molecules are known to have regulatory impacts in cells but are difficult to locate in the genome. The minimal requirements for G-quadruplex folding in RNA (G_≥3N_1-7G_≥3N_1-7G_≥3N_1-7G_≥3) is being challenged by observations made on specific examples in recent years. The definition of potential G-quadruplex sequences has major repercussions on the observation of the structure since it introduces a bias. The canonical motif only describes a sub-population of the reported G-quadruplexes. To address these issues, we propose an RNA G-quadruplex prediction strategy that does not rely on a motif definition.

Results

We trained an artificial neural network with sequences of experimentally validated G-quadruplexes from the G4RNA database encoded using an abstract definition of their sequence. This artificial neural network, G4NN, evaluates the similarity of a given sequence to known G-quadruplexes and reports it as a score. G4NN has a predictive power comparable to the reported G richness and G/C skewness evaluations that are the current state-of-the-art for the identification of potential RNA G-quadruplexes. We combined these approaches in the G4RNA screener, a program designed to manage and evaluate the sequences to identify potential G-quadruplexes.

Availability and implementation

G4RNA screener is available for download at http://gitlabscottgroup.med.usherbrooke.ca/J-Michel/g4rna_screener.

Supplementary information

Supplementary data are available at Bioinformatics online.

Mimicking Co-Transcriptional RNA Folding Using a Superhelicase

Vectorial folding of RNA during transcription can produce intermediates with distinct biochemical activities. Here, we design an artificial minimal system to mimic cotranscriptional RNA folding in vitro. In this system, a presynthesized RNA molecule begins to fold from its 5'-end, as it is released from a heteroduplex by an engineered helicase that translocates on the complementary DNA strand in the 3'-to-5' direction. This chemically stabilized "superhelicase" Rep-X processively unwinds thousands of base pairs of DNA. The presynthesized RNA enables us to flexibly position fluorescent labels on the RNA for single-molecule fluorescence resonance energy transfer analysis and allows us to study real-time conformational dynamics during the vectorial folding process. We observed distinct signatures of the maiden secondary and tertiary folding of the Oryza sativa twister ribozyme. The maiden vectorial tertiary folding transitions occurred faster than Mg²⁺-induced refolding, but were also more prone to misfolding, likely due to sequential formation of alternative secondary structures. This novel assay can be applied to studying other kinetically controlled processes, such as riboswitch control and RNA-protein assembly.

Predicting Cotranscriptional Folding Kinetics For Riboswitch

On the basis of a helix-based transition rate model, we developed a new method for sampling cotranscriptional RNA conformational ensemble and the prediction of cotranscriptional folding kinetics. Applications to E. coli. SRP RNA and pbuE riboswitch indicate that the model may provide reliable predictions for the cotranscriptional folding pathways and population kinetics. For E. coli. SRP RNA, the predicted population kinetics and the folding pathway are consistent with the SHAPE profiles in the recent cotranscriptional SHAPE-seq experiments. For the pbuE riboswitch, the model predicts the transcriptional termination efficiency as a function of the force. The theoretical results show (a) a force-induced transition from the aptamer (antiterminator) to the terminator structure and (b) the different folding pathways for the riboswitch with and without the ligand (adenine). More specifically, without adenine, the aptamer structure emerges as a short-lived kinetic transient state instead of a thermodynamically stable intermediate state. Furthermore, from the predicted extension-time curves, the model identifies a series of conformational switches in the pulling process, where the predicted relative residence times for the different structures are in accordance with the experimental data. The model may provide a new tool for quantitative predictions of cotranscriptional folding kinetics, and results can offer useful insights into cotranscriptional folding-related RNA functions such as regulation of gene expression with riboswitches.

Computational approaches for alternative and transient secondary structures of ribonucleic acids

Transient and alternative structures of ribonucleic acids (RNAs) play essential roles in various regulatory processes, such as translation regulation in living cells. Because experimental analyses for RNA structures are difficult and time-consuming, computational approaches based on RNA secondary structures are promising. In this article, we review computational methods for detecting and analyzing transient/alternative secondary structures of RNAs, including static approaches based on probabilistic distributions of RNA secondary structures and dynamic approaches such as kinetic folding and folding pathway predictions.

Towards Long-Range RNA Structure Prediction in Eukaryotic Genes

The ability to form an intramolecular structure plays a fundamental role in eukaryotic RNA biogenesis. Proximate regions in the primary transcripts fold into a local secondary structure, which is then hierarchically assembled into a tertiary structure that is stabilized by RNA-binding proteins and long-range intramolecular base pairings. While the local RNA structure can be predicted reasonably well for short sequences, long-range structure at the scale of eukaryotic genes remains problematic from the computational standpoint. The aim of this review is to list functional examples of long-range RNA structures, to summarize current comparative methods of structure prediction, and to highlight their advances and limitations in the context of long-range RNA structures. Most comparative methods implement the “first-align-then-fold” principle, i.e., they operate on multiple sequence alignments, while functional RNA structures often reside in non-conserved parts of the primary transcripts. The opposite “first-fold-then-align” approach is currently explored to a much lesser extent. Developing novel methods in both directions will improve the performance of comparative RNA structure analysis and help discover novel long-range structures, their higher-order organization, and RNA⁻RNA interactions across the transcriptome.

Bacterial riboswitches and RNA thermometers: Nature and contributions to pathogenesis

Bacterial pathogens are always challenged by fluctuations of chemical and physical parameters that pose serious threats to cellular integrity and metabolic status. Sudden deprivation of nutrients or key metabolites, changes in surrounding pH, and temperature shifts are the most important examples of such parameters. To elicit a proper response to such fluctuations, bacterial cells coordinate the expression of parameter-relevant genes. Although protein-mediated control of gene expression is well appreciated since many decades, RNA-based regulation has been discovered in early 2000s as a parallel level of regulation. Small regulatory RNAs have emerged as one of the most widespread and important gene regulatory systems in bacteria with rare representatives found in Archaea and Eukarya. Riboswitches and thermosensors are cis-encoded RNA regulatory elements that employ different mechanisms to regulate the expression of related genes controlling key metabolic pathways and genes of temperature relevant proteins including virulence factors. The extent of RNA contributions to gene regulation is not completely known even in well-studied models such E. coli and B. subtilis. In depth understanding of riboswitches is promising for opportunity to discover a narrow spectrum antibacterial drugs that target riboswitches of essential metabolic pathways.

Allele-specific SHAPE-MaP assessment of the effects of somatic variation and protein binding on mRNA structure

The impact of inherited and somatic mutations on messenger RNA (mRNA) structure remains poorly understood. Recent technological advances that leverage next-generation sequencing to obtain experimental structure data, such as SHAPE-MaP, can reveal structural effects of mutations, especially when these data are incorporated into structure modeling. Here, we analyze the ability of SHAPE-MaP to detect the relatively subtle structural changes caused by single-nucleotide mutations. We find that allele-specific sorting greatly improved our detection ability. Thus, we used SHAPE-MaP with a novel combination of clone-free robotic mutagenesis and allele-specific sorting to perform a rapid, comprehensive survey of noncoding somatic and inherited riboSNitches in two cancer-associated mRNAs, TPT1 and LCP1 Using rigorous thermodynamic modeling of the Boltzmann suboptimal ensemble, we identified a subset of mutations that change TPT1 and LCP1 RNA structure, with approximately 14% of all variants identified as riboSNitches. To confirm that these in vitro structures were biologically relevant, we tested how dependent TPT1 and LCP1 mRNA structures were on their environments. We performed SHAPE-MaP on TPT1 and LCP1 mRNAs in the presence or absence of cellular proteins and found that both mRNAs have similar overall folds in all conditions. RiboSNitches identified within these mRNAs in vitro likely exist under biological conditions. Overall, these data reveal a robust mRNA structural landscape where differences in environmental conditions and most sequence variants do not significantly alter RNA structural ensembles. Finally, predicting riboSNitches in mRNAs from sequence alone remains particularly challenging; these data will provide the community with benchmarks for further algorithmic development.

Transcription elongation rate affects nascent histone pre-mRNA folding and 3' end processing

In this study, Saldi et al. investigated how transcription elongation rate influences cotranscriptional pre-mRNA maturation. Their findings show that regulation of transcription speed can modulate pre-mRNA processing by changing nascent RNA structure and suggest a mechanism by which alternative processing could be controlled.

Transcription elongation rate influences cotranscriptional pre-mRNA maturation, but how such kinetic coupling works is poorly understood. The formation of nonadenylated histone mRNA 3′ ends requires recognition of an RNA structure by stem–loop-binding protein (SLBP). We report that slow transcription by mutant RNA polymerase II (Pol II) caused accumulation of polyadenylated histone mRNAs that extend past the stem–loop processing site. UV irradiation, which decelerates Pol II elongation, also induced long poly(A)⁺ histone transcripts. Inhibition of 3′ processing by slow Pol II correlates with failure to recruit SLBP to histone genes. Chemical probing of nascent RNA structure showed that the stem–loop fails to fold in transcripts made by slow Pol II, thereby explaining the absence of SLBP and failure to process 3′ ends. These results show that regulation of transcription speed can modulate pre-mRNA processing by changing nascent RNA structure and suggest a mechanism by which alternative processing could be controlled.

Identifying and Characterizing circRNA-Protein Interaction

Circular RNAs have been identified as naturally occurring RNAs that are highly represented in the eukaryotic transcriptome. Although a large number of circRNAs have been reported, circRNA functions remain largely unknown. CircRNAs can function as miRNA sponges, thereby reducing their ability to target mRNAs. We hypothesize that circRNAs may bind, store, sort, and sequester proteins to particular subcellular locations, and act as dynamic scaffolding molecules that modulate protein-protein interactions. Here, we review the biological implication and function of circRNA-protein interaction, and reveal a dynamic model of the interaction in various tissues, development stages and physiological conditions. Improved techniques to identify and characterize the dynamic RNA-protein interactions may elucidate the molecular mechanisms associated with the expression and functional diversity of circRNAs.

Optimization of a novel biophysical model using large scale in vivo antisense hybridization data displays improved prediction capabilities of structurally accessible RNA regions

Current approaches to design efficient antisense RNAs (asRNAs) rely primarily on a thermodynamic understanding of RNA–RNA interactions. However, these approaches depend on structure predictions and have limited accuracy, arguably due to overlooking important cellular environment factors. In this work, we develop a biophysical model to describe asRNA–RNA hybridization that incorporates in vivo factors using large-scale experimental hybridization data for three model RNAs: a group I intron, CsrB and a tRNA. A unique element of our model is the estimation of the availability of the target region to interact with a given asRNA using a differential entropic consideration of suboptimal structures. We showcase the utility of this model by evaluating its prediction capabilities in four additional RNAs: a group II intron, Spinach II, 2-MS2 binding domain and glgC 5΄ UTR. Additionally, we demonstrate the applicability of this approach to other bacterial species by predicting sRNA–mRNA binding regions in two newly discovered, though uncharacterized, regulatory RNAs.

In vivo probing of nascent RNA structures reveals principles of cotranscriptional folding

Defining the in vivo folding pathway of cellular RNAs is essential to understand how they reach their final native conformation. We here introduce a novel method, named Structural Probing of Elongating Transcripts (SPET-seq), that permits single-base resolution analysis of transcription intermediates' secondary structures on a transcriptome-wide scale, enabling base-resolution analysis of the RNA folding events. Our results suggest that cotranscriptional RNA folding in vivo is a mixture of cooperative folding events, in which local RNA secondary structure elements are formed as they get transcribed, and non-cooperative events, in which 5'-halves of long-range helices get sequestered into transient non-native interactions until their 3' counterparts have been transcribed. Together our work provides the first transcriptome-scale overview of RNA cotranscriptional folding in a living organism.

Splicing and transcription touch base: co-transcriptional spliceosome assembly and function

Several macromolecular machines collaborate to produce eukaryotic messenger RNA. RNA polymerase II (Pol II) translocates along genes that are up to millions of base pairs in length and generates a flexible RNA copy of the DNA template. This nascent RNA harbours introns that are removed by the spliceosome, which is a megadalton ribonucleoprotein complex that positions the distant ends of the intron into its catalytic centre. Emerging evidence that the catalytic spliceosome is physically close to Pol II in vivo implies that transcription and splicing occur on similar timescales and that the transcription and splicing machineries may be spatially constrained. In this Review, we discuss aspects of spliceosome assembly, transcription elongation and other co-transcriptional events that allow the temporal coordination of co-transcriptional splicing.

Lariat capping as a tool to manipulate the 5' end of individual yeast mRNA species in vivo

The 5' cap structure of eukaryotic mRNA is critical for its processing, transport, translation, and stability. The many functions of the cap and the fact that most, if not all, mRNA carries the same type of cap makes it difficult to analyze cap function in vivo at individual steps of gene expression. We have used the lariat capping ribozyme (LCrz) from the myxomycete Didymium to replace the mRNA m⁷G cap of a single reporter mRNA species with a tiny lariat in which the first and the third nucleotide are joined by a 2', 5' phosphodiester bond. We show that the ribozyme functions in vivo in the budding yeast Saccharomyces cerevisiae presumably without cofactors and that lariat capping occurs cotranscriptionally. The lariat-capped reporter mRNA is efficiently exported to the cytoplasm where it is found to be oligoadenylated and evenly distributed. Both the oligoadenylated form and a lariat-capped mRNA with a templated poly(A) tail translates poorly, underlining the critical importance of the m⁷G cap in translation. Finally, the lariat-capped RNA exhibits a threefold longer half-life compared to its m⁷G-capped counterpart, consistent with a key role for the m⁷G cap in mRNA turnover. Our study emphasizes important activities of the m⁷G cap and suggests new utilities of lariat capping as a molecular tool in vivo.

In silico methods for co-transcriptional RNA secondary structure prediction and for investigating alternative RNA structure expression

RNA transcripts are the primary products of active genes in any living organism, including many viruses. Their cellular destiny not only depends on primary sequence signals, but can also be determined by RNA structure. Recent experimental evidence shows that many transcripts can be assigned more than a single functional RNA structure throughout their cellular life and that structure formation happens co-transcriptionally, i.e. as the transcript is synthesised in the cell. Moreover, functional RNA structures are not limited to non-coding transcripts, but can also feature in coding transcripts. The picture that now emerges is that RNA structures constitute an additional layer of information that can be encoded in any RNA transcript (and on top of other layers of information such as protein-context) in order to exert a wide range of functional roles. Moreover, different encoded RNA structures can be expressed at different stages of a transcript's life in order to alter the transcript's behaviour depending on its actual cellular context. Similar to the concept of alternative splicing for protein-coding genes, where a single transcript can yield different proteins depending on cellular context, it is thus appropriate to propose the notion of alternative RNA structure expression for any given transcript. This review introduces several computational strategies that my group developed to detect different aspects of RNA structure expression in vivo. Two aspects are of particular interest to us: (1) RNA secondary structure features that emerge during co-transcriptional folding and (2) functional RNA structure features that are expressed at different times of a transcript's life and potentially mutually exclusive.

Capturing RNA Folding Free Energy with Coarse-Grained Molecular Dynamics Simulations

We introduce a coarse-grained RNA model for molecular dynamics simulations, RACER (RnA CoarsE-gRained). RACER achieves accurate native structure prediction for a number of RNAs (average RMSD of 2.93 Å) and the sequence-specific variation of free energy is in excellent agreement with experimentally measured stabilities (R2 = 0.93). Using RACER, we identified hydrogen-bonding (or base pairing), base stacking, and electrostatic interactions as essential driving forces for RNA folding. Also, we found that separating pairing vs. stacking interactions allowed RACER to distinguish folded vs. unfolded states. In RACER, base pairing and stacking interactions each provide an approximate stability of 3-4 kcal/mol for an A-form helix. RACER was developed based on PDB structural statistics and experimental thermodynamic data. In contrast with previous work, RACER implements a novel effective vdW potential energy function, which led us to re-parameterize hydrogen bond and electrostatic potential energy functions. Further, RACER is validated and optimized using a simulated annealing protocol to generate potential energy vs. RMSD landscapes. Finally, RACER is tested using extensive equilibrium pulling simulations (0.86 ms total) on eleven RNA sequences (hairpins and duplexes).

Integration of mRNP formation and export

Expression of protein-coding genes in eukaryotes relies on the coordinated action of many sophisticated molecular machineries. Transcription produces precursor mRNAs (pre-mRNAs) and the active gene provides an environment in which the pre-mRNAs are processed, folded, and assembled into RNA–protein (RNP) complexes. The dynamic pre-mRNPs incorporate the growing transcript, proteins, and the processing machineries, as well as the specific protein marks left after processing that are essential for export and the cytoplasmic fate of the mRNPs. After release from the gene, the mRNPs move by diffusion within the interchromatin compartment, making up pools of mRNPs. Here, splicing and polyadenylation can be completed and the mRNPs recruit the major export receptor NXF1. Export competent mRNPs interact with the nuclear pore complex, leading to export, concomitant with compositional and conformational changes of the mRNPs. We summarize the integrated nuclear processes involved in the formation and export of mRNPs.

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.

The RNA Epistructurome: Uncovering RNA Function by Studying Structure and Post-Transcriptional Modifications

A large fraction of higher metazoan genomes transcribe RNA molecules whose functions extend far beyond carrying instructions for protein synthesis. Although RNA is apparently a simple molecule, the ways in which it performs many of its functions have remained highly elusive for decades. As learned from studying ribosomal and transfer RNAs, two of the key features influencing the function of RNA are its structure and post-transcriptional modifications. A deep understanding of RNA function therefore requires rapid and straightforward approaches to study the complex and intricate landscape of RNA structures and modifications. In this review we summarize and discuss the most recent methods and findings in the field of RNA biology, with an eye toward new frontiers and open questions.

Inverted repeat Alu elements in the human lincRNA-p21 adopt a conserved secondary structure that regulates RNA function

LincRNA-p21 is a long intergenic non-coding RNA (lincRNA) involved in the p53-mediated stress response. We sequenced the human lincRNA-p21 (hLincRNA-p21) and found that it has a single exon that includes inverted repeat Alu elements (IRAlus). Sense and antisense Alu elements fold independently of one another into a secondary structure that is conserved in lincRNA-p21 among primates. Moreover, the structures formed by IRAlus are involved in the localization of hLincRNA-p21 in the nucleus, where hLincRNA-p21 colocalizes with paraspeckles. Our results underscore the importance of IRAlus structures for the function of hLincRNA-p21 during the stress response.

The transcription initiation sites of eggplant latent viroid strands map within distinct motifs in their in vivo RNA conformations

Eggplant latent viroid (ELVd), like other members of family Avsunviroidae, replicates in plastids through a symmetric rolling-circle mechanism in which elongation of RNA strands is most likely catalyzed by a nuclear-encoded polymerase (NEP) translocated to plastids. Here we have addressed where NEP initiates transcription of viroid strands. Because this step is presumably directed by sequence/structural motifs, we have previously determined the conformation of the monomeric linear (+) and (-) RNAs of ELVd resulting from hammerhead-mediated self-cleavage. In silico predictions with 3 softwares led to similar bifurcated conformations for both ELVd strands. In vitro examination by non-denaturing PAGE showed that they migrate as prominent single bands, with the ELVd (+) RNA displaying a more compact conformation as revealed by its faster electrophoretic mobility. In vitro SHAPE analysis corroborated the ELVd conformations derived from thermodynamics-based predictions in silico. Moreover, sequence analysis of 94 full-length natural ELVd variants disclosed co-variations, and mutations converting canonical into wobble pairs or vice versa, which confirmed in vivo most of the stems predicted in silico and in vitro, and additionally helped to introduce minor structural refinements. Therefore, results from the 3 experimental approaches were essentially consistent among themselves. Application to RNA preparations from ELVd-infected tissue of RNA ligase-mediated rapid amplification of cDNA ends, combined with pretreatments to modify the 5' ends of viroid strands, mapped the transcription initiation sites of ELVd (+) and (-) strands in vivo at different sequence/structural motifs, in contrast with the situation previously observed in 2 other members of the family Avsunviroidae.

Recombinant production, purification and characterization of vessel dilator in E. coli

Vessel dilator is a 3.9-KDa potent anticancer peptide and a valuable candidate in the treatment of conditions such as congestive heart failure and acute renal failure amongst others. Here we report the recombinant production of vessel dilator in Escherichia coli. Three different synthetic ORF's dubbed VDI, VDII and VDIII, each encoding a trimmer of the vessel dilator peptide attached to a His tag sequence at their C- terminal, were synthesized and placed in pET21c expression vectors. The highest yield, following expression in E. coli BL21 (DE3), was recorded with VDII that carried the shortest fusion partner. Subsequent to the initial capture of the fusion protein by a Ni affinity column, the vessel dilator monomers were cleaved by trypsin treatment, and further purified to at least 90% homogeneity by anion exchange chromatography. De-novo sequencing and in vivo anticancer activity tests were used to verify the peptide sequence and its biological activity, respectively. The final yield was estimated to be approximately 15 mg of the purified vessel dilator per gram wet weight of the bacterial cells.

Internal guide RNA interactions interfere with Cas9-mediated cleavage

The CRISPR/Cas system uses guide RNAs (gRNAs) to direct sequence-specific DNA cleavage. Not every gRNA elicits cleavage and the mechanisms that govern gRNA activity have not been resolved. Low activity could result from either failure to form a functional Cas9-gRNA complex or inability to recognize targets in vivo. Here we show that both phenomena influence Cas9 activity by comparing mutagenesis rates in zebrafish embryos with in vitro cleavage assays. In vivo, our results suggest that genomic factors such as CTCF inhibit mutagenesis. Comparing near-identical gRNA sequences with different in vitro activities reveals that internal gRNA interactions reduce cleavage. Even though gRNAs containing these structures do not yield cleavage-competent complexes, they can compete with active gRNAs for binding to Cas9. These results reveal that both genomic context and internal gRNA interactions can interfere with Cas9-mediated cleavage and illuminate previously uncharacterized features of Cas9-gRNA complex formation.

The prediction of virus mutation using neural networks and rough set techniques

Viral evolution remains to be a main obstacle in the effectiveness of antiviral treatments. The ability to predict this evolution will help in the early detection of drug-resistant strains and will potentially facilitate the design of more efficient antiviral treatments. Various tools has been utilized in genome studies to achieve this goal. One of these tools is machine learning, which facilitates the study of structure-activity relationships, secondary and tertiary structure evolution prediction, and sequence error correction. This work proposes a novel machine learning technique for the prediction of the possible point mutations that appear on alignments of primary RNA sequence structure. It predicts the genotype of each nucleotide in the RNA sequence, and proves that a nucleotide in an RNA sequence changes based on the other nucleotides in the sequence. Neural networks technique is utilized in order to predict new strains, then a rough set theory based algorithm is introduced to extract these point mutation patterns. This algorithm is applied on a number of aligned RNA isolates time-series species of the Newcastle virus. Two different data sets from two sources are used in the validation of these techniques. The results show that the accuracy of this technique in predicting the nucleotides in the new generation is as high as 75 %. The mutation rules are visualized for the analysis of the correlation between different nucleotides in the same RNA sequence.

Mutational analysis of structural elements in a class-I cyclic di-GMP riboswitch to elucidate its regulatory mechanism

The Vc2 riboswitch possesses an aptamer domain belonging to the class-I c-di-GMP riboswitch family. This domain has been analysed and the molecular mechanism by which it recognizes the c-di-GMP ligand has been elucidated. On the other hand, the regulatory mechanism of the full-length Vc2 riboswitch to control its downstream open reading frame (ORF) remains largely unknown. In this study, we performed in vivo reporter assays and in vitro biochemical analyses of the full-length riboswitch and its aptamer domain. We evaluated the results of in vivo and in vitro analyses to elucidate the regulatory mechanism of the Vc2 riboswitch. The present results suggest that recognition of c-di-GMP ligand by the Vc2 riboswitch aptamer domain downregulates expression of its downstream ORF primarily at the translational level.