Pervasive Regulatory Functions of mRNA Structure Revealed by High-Resolution SHAPE Probing

Revealing miRNAs patterns by employing matrix representations and energy analysis

MicroRNAs (miRNAs) are small, non-coding RNA molecules that regulate gene expression. Despite their relatively short length (about 21 nucleotides), they can regulate thousands of transcripts within a cell. Due to their low complementarity to targets, studying their activity and binding region preferences (3'UTR, 5'UTR, or CDS) is challenging. In this paper, we analyzed a set of human miRNAs to uncover their general patterns. We began with a sequence logo to verify conservation at specific positions. To discover long-range correlations, we employed chaos game representation (CGR) and genomatrix, methods that enable both graphical and analytical analysis of sequence sets and are well-established in bioinformatics. Our results showed that miRNAs exhibit strongly non-random and characteristic patterns. To incorporate physicochemical properties into the analysis, we applied the electron-ion interaction potential (EIIP) parameter. An important part of our study was to validate the division of miRNAs into two parts-seed and puzzle. The seed region is responsible for target binding, while the puzzle region likely interacts with the RISC complex. We estimated duplex binding energy within the 3'UTR, 5'UTR, and CDS regions using the miRanda tool. Based on the median energy distribution, we divided the miRNAs into two subsets, reflecting different patterns in chaos game representation. Interestingly, one subset displayed significant similarity to conserved and highly confidential miRNAs. Our results confirm the low complementarity of miRNA/mRNA interactions and support the functional division of miRNA structure. Additionally, we present findings related to the localization of transcript target sites, which form the basis for further analyses.

Genome-Wide Identification of Stable RNA Secondary Structures Across Multiple Organisms Using Chemical Probing Data: Insights into Short Structural Motifs and RNA-Targeting Therapeutics

Small molecules targeting specific RNA binding sites, including stable and transient RNA structures, are emerging as effective pharmacological approaches for modulating gene expression. However, little is understood about how stable RNA secondary structures are shared across organisms, an important factor in controlling drug selectivity. In this study, I provide an analytical pipeline named RNA Secondary Structure Finder (R2S-Finder) to discover short, stable RNA structural motifs for humans, Escherichia coli ( E. coli ), SARS-CoV-2, and Zika virus by leveraging existing in vivo and in vitro genome-wide chemical RNA-probing datasets. I found several common features across organisms. For example, apart from the well-documented tetraloops, AU-rich tetraloops are widely present in different organisms. I also found that the 5' untranslated region (UTR) contains a higher proportion of stable structures than the coding sequences in humans, SARS-CoV-2, and Zika virus. In general, stable structures predicted from in vitro (protein-free) and in vivo datasets are consistent in humans, E. coli , and SARS-CoV-2, indicating that most stable structure formation were driven by RNA folding alone, while a larger variation was found between in vitro and in vivo data with certain RNA types, such as human long intergenic non-coding RNAs (lincRNAs). Finally, I predicted stable three- and four-way RNA junctions that exist both in vivo and in vitro conditions, which can potentially serve as drug targets. All results of stable sequences, stem-loops, internal loops, bulges, and three- and four-way junctions have been collated in the R2S-Finder database ( https://github.com/JingxinWangLab/R2S-Finder ), which is coded in hyperlinked HTML pages and tabulated in CSV files.

Beyond ligand binding: Single molecule observation reveals how riboswitches integrate multiple signals to balance bacterial gene regulation

Riboswitches are specialized RNA structures that orchestrate gene expression in response to sensing specific metabolite or ion ligands, mostly in bacteria. Upon ligand binding, these conformationally dynamic RNA motifs undergo structural changes that control critical gene expression processes such as transcription termination and translation initiation, thereby enabling cellular homeostasis and adaptation. Because RNA folds rapidly and co-transcriptionally, riboswitches make use of the low complexity of RNA sequences to adopt alternative, transient conformations on the heels of the transcribing RNA polymerase (RNAP), resulting in kinetic partitioning that defines the regulatory outcome. This review summarizes single molecule microscopy evidence that has begun to unveil a sophisticated network of dynamic, kinetically balanced interactions between riboswitch architecture and the gene expression machinery that, together, integrate diverse cellular signals.

Identification of RNA structures and their roles in RNA functions

The development of high-throughput RNA structure profiling methods in the past decade has greatly facilitated our ability to map and characterize different aspects of RNA structures transcriptome-wide in cell populations, single cells and single molecules. The resulting high-resolution data have provided insights into the static and dynamic nature of RNA structures, revealing their complexity as they perform their respective functions in the cell. In this Review, we discuss recent technical advances in the determination of RNA structures, and the roles of RNA structures in RNA biogenesis and functions, including in transcription, processing, translation, degradation, localization and RNA structure-dependent condensates. We also discuss the current understanding of how RNA structures could guide drug design for treating genetic diseases and battling pathogenic viruses, and highlight existing challenges and future directions in RNA structure research.

Mechanistic analysis of Riboswitch Ligand interactions provides insights into pharmacological control over gene expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Identification of leader-trailer helices of precursor ribosomal RNA in all phyla of bacteria and archaea

Ribosomal RNAs are transcribed as part of larger precursor molecules. In Escherichia coli, complementary RNA segments flank each rRNA and form long leader-trailer (LT) helices, which are crucial for subunit biogenesis in the cell. A previous study of 15 representative species suggested that most but not all prokaryotes contain LT helices. Here, we use a combination of in silico folding and covariation methods to identify and characterize LT helices in 4464 bacterial and 260 archaeal organisms. Our results suggest that LT helices are present in all phyla, including Deinococcota, which had previously been suspected to lack LT helices. In very few organisms, our pipeline failed to detect LT helices for both 16S and 23S rRNA. However, a closer case-by-case look revealed that LT helices are indeed present but escaped initial detection. Over 3600 secondary structure models, many well supported by nucleotide covariation, were generated. These structures show a high degree of diversity. Yet, all exhibit extensive base-pairing between the leader and trailer strands, in line with a common and essential function.

Hfq mediates transcriptome-wide RNA structurome reprogramming under virulence-inducing conditions in a phytopathogen

Although RNA structures play important roles in regulating gene expression, the mechanism and function of mRNA folding in plant bacterial pathogens remain elusive. Therefore, we perform dimethyl sulfate sequencing (DMS-seq) on the Pseudomonas syringae under nutrition-rich and -deficient conditions, revealing that the mRNA structure changes substantially in the minimal medium (MM) that tunes global translation efficiency (TE), thereby inducing virulence. This process is led by the increased expression of hfq, which is directly activated by transcription regulators RpoS and CysB. The co-occurrence of Hfq and RpoS in diverse bacteria and the deep conservation of Hfq Y25 is critical for RNA-mediated regulation and implicates the wider biological importance of mRNA structure and feedback loops in the control of global gene expression.

Characterization of the dual regulation by a c-di-GMP riboswitch Bc1 with a long expression platform from Bacillus thuringiensis

A riboswitch generally regulates the expression of its downstream genes through conformational change in its expression platform (EP) upon ligand binding. The cyclic diguanosine monophosphate (c-di-GMP) class I riboswitch Bc1 is widespread and conserved among Bacillus cereus group species. In this study, we revealed that Bc1 has a long EP with two typical ρ-independent terminator sequences 28 bp apart. The upstream terminator T1 is dominant in vitro, while downstream terminator T2 is more efficient in vivo. Through mutation analysis, we elucidated that Bc1 exerts a rare and incoherent "transcription-translation" dual regulation with T2 playing a crucial role. However, we found that Bc1 did not respond to c-di-GMP under in vitro transcription conditions, and the expressions of downstream genes did not change with fluctuation in intracellular c-di-GMP concentration. To explore this puzzle, we conducted SHAPE-MaP and confirmed the interaction of Bc1 with c-di-GMP. This shows that as c-di-GMP concentration increases, T1 unfolds but T2 remains almost intact and functional. The presence of T2 masks the effect of T1 unwinding, resulting in no response of Bc1 to c-di-GMP. The high Shannon entropy values of EP region imply the potential alternative structures of Bc1. We also found that zinc uptake regulator can specifically bind to the dual terminator coding sequence and slightly trigger the response of Bc1 to c-di-GMP. This work will shed light on the dual-regulation riboswitch and enrich our understanding of the RNA world.IMPORTANCEIn nature, riboswitches are involved in a variety of metabolic regulation, most of which preferentially regulate transcription termination or translation initiation of downstream genes in specific ways. Alternatively, the same or different riboswitches can exist in tandem to enhance regulatory effects or respond to multiple ligands. However, many putative conserved riboswitches have not yet been experimentally validated. Here, we found that the c-di-GMP riboswitch Bc1 with a long EP could form a dual terminator and exhibit non-canonical and incoherent "transcription-translation" dual regulation. Besides, zinc uptake regulator specifically bound to the coding sequence of the Bc1 EP and slightly mediated the action of Bc1. The application of SHAPE-MaP to the dual regulation mechanism of Bc1 may establish the foundation for future studies of such complex untranslated regions in other bacterial genomes.

Massively parallel dissection of RNA in RNA-protein interactions in vivo

Many of the biological functions performed by RNA are mediated by RNA-binding proteins (RBPs), and understanding the molecular basis of these interactions is fundamental to biology. Here, we present massively parallel RNA assay combined with immunoprecipitation (MPRNA-IP) for in vivo high-throughput dissection of RNA-protein interactions and describe statistical models for identifying RNA domains and parsing the structural contributions of RNA. By using custom pools of tens of thousands of RNA sequences containing systematically designed truncations and mutations, MPRNA-IP is able to identify RNA domains, sequences, and secondary structures necessary and sufficient for protein binding in a single experiment. We show that this approach is successful for multiple RNAs of interest, including the long noncoding RNA NORAD, bacteriophage MS2 RNA, and human telomerase RNA, and we use it to interrogate the hitherto unknown sequence or structural RNA-binding preferences of the DNA-looping factor CTCF. By integrating systematic mutation analysis with crosslinking immunoprecipitation, MPRNA-IP provides a novel high-throughput way to elucidate RNA-based mechanisms behind RNA-protein interactions in vivo.

Phylogenetic and Chemical Probing Information as Soft Constraints in RNA Secondary Structure Prediction

Extrinsic, experimental information can be incorporated into thermodynamics-based RNA folding algorithms in the form of pseudo-energies. Evolutionary conservation of RNA secondary structure elements is detectable in alignments of phylogenetically related sequences and provides evidence for the presence of certain base pairs that can also be converted into pseudo-energy contributions. We show that the centroid base pairs computed from a consensus folding model such as RNAalifold result in a substantial improvement of the prediction accuracy for single sequences. Evidence for specific base pairs turns out to be more informative than a position-wise profile for the conservation of the pairing status. A comparison with chemical probing data, furthermore, strongly suggests that phylogenetic base pairing data are more informative than position-specific data on (un)pairedness as obtained from chemical probing experiments. In this context we demonstrate, in addition, that the conversion of signal from probing data into pseudo-energies is possible using thermodynamic structure predictions as a reference instead of known RNA structures.

A new mRNA structure prediction based approach to identifying improved signal peptides for bone morphogenetic protein 2

BACKGROUND

Signal peptide (SP) engineering has proven able to improve production of many proteins yet is a laborious process that still relies on trial and error. mRNA structure around the translational start site is important in translation initiation and has rarely been considered in this context, with recent improvements in in silico mRNA structure potentially rendering it a useful predictive tool for SP selection. Here we attempt to create a method to systematically screen candidate signal peptide sequences in silico based on both their nucleotide and amino acid sequences. Several recently released computational tools were used to predict signal peptide activity (SignalP), localization target (DeepLoc) and predicted mRNA structure (MXFold2). The method was tested with Bone Morphogenetic Protein 2 (BMP2), an osteogenic growth factor used clinically for bone regeneration. It was hoped more effective BMP2 SPs could improve BMP2-based gene therapies and reduce the cost of recombinant BMP2 production.

RESULTS

Amino acid sequence analysis indicated 2,611 SPs from the TGF-β superfamily were predicted to function when attached to BMP2. mRNA structure prediction indicated structures at the translational start site were likely highly variable. The five sequences with the most accessible translational start sites, a codon optimized BMP2 SP variant and the well-established hIL2 SP sequence were taken forward to in vitro testing. The top five candidates showed non-significant improvements in BMP2 secretion in HEK293T cells. All showed reductions in secretion versus the native sequence in C2C12 cells, with several showing large and significant decreases. None of the tested sequences were able to increase alkaline phosphatase activity above background in C2C12s. The codon optimized control sequence and hIL2 SP showed reasonable activity in HEK293T but very poor activity in C2C12.

CONCLUSIONS

These results support the use of peptide sequence based in silico tools for basic predictions around signal peptide activity in a synthetic biology context. However, mRNA structure prediction requires improvement before it can produce reliable predictions for this application. The poor activity of the codon optimized BMP2 SP variant in C2C12 emphasizes the importance of codon choice, mRNA structure, and cellular context for SP activity.

RNA-Based Sensor Systems for Affordable Diagnostics in the Age of Pandemics

In the era of the COVID-19 pandemic, the significance of point-of-care (POC) diagnostic tools has become increasingly vital, driven by the need for quick and precise virus identification. RNA-based sensors, particularly toehold sensors, have emerged as promising candidates for POC detection systems due to their selectivity and sensitivity. Toehold sensors operate by employing an RNA switch that changes the conformation when it binds to a target RNA molecule, resulting in a detectable signal. This review focuses on the development and deployment of RNA-based sensors for POC viral RNA detection with a particular emphasis on toehold sensors. The benefits and limits of toehold sensors are explored, and obstacles and future directions for improving their performance within POC detection systems are presented. The use of RNA-based sensors as a technology for rapid and sensitive detection of viral RNA holds great potential for effectively managing (dealing/coping) with present and future pandemics in resource-constrained settings.

Effects of mRNA conformational switching on translational noise in gene circuits

Intragenic translational heterogeneity describes the variation in translation at the level of transcripts for an individual gene. A factor that contributes to this source of variation is the mRNA structure. Both the composition of the thermodynamic ensemble, i.e., the stationary distribution of mRNA structures, and the switching dynamics between those play a role. The effect of the switching dynamics on intragenic translational heterogeneity remains poorly understood. We present a stochastic translation model that accounts for mRNA structure switching and is derived from a Markov model via approximate stochastic filtering. We assess the approximation on various timescales and provide a method to quantify how mRNA structure dynamics contributes to translational heterogeneity. With our approach, we allow quantitative information on mRNA switching from biophysical experiments or coarse-grain molecular dynamics simulations of mRNA structures to be included in gene regulatory chemical reaction network models without an increase in the number of species. Thereby, our model bridges a gap between mRNA structure kinetics and gene expression models, which we hope will further improve our understanding of gene regulatory networks and facilitate genetic circuit design.

A new look at the fluorescent protein-based approach for identifying optimal coding sequence for recombinant protein expression in E. coli

Due to the degeneracy of the genetic code, most amino acids are encoded by several codons. The choice among synonymous codons at the N-terminus of genes has a profound effect on protein expression in Escherichia coli. This is often explained by the different contributions of synonymous codons to mRNA secondary structure formation. Strong secondary structures at the 5'-end of mRNA interfere with ribosome binding and affect the process of translation initiation. In silico optimization of the gene 5'-end can significantly increase the level of protein expression; however, this method is not always effective due to the uncertainty of the exact mechanism by which synonymous substitutions affect expression; thus, it may produce nonoptimal variants as well as miss some of the best producers. In this paper, an alternative approach is proposed based on screening a partially randomized library of expression constructs comprising hundreds of selected synonymous variants. The effect of such substitutions was evaluated using the gene of interest fused to the reporter gene of the fluorescent protein with subsequent screening for the most promising candidates according to the reporter's signal intensity. The power of the approach is demonstrated by a significant increase in the prokaryotic expression of three proteins: canine cystatin C, human BCL2-associated athanogene 3 and human cardiac troponin I. This simple approach was suggested which may provide an efficient, easy, and inexpensive optimization method for poorly expressed proteins in bacteria.

Extraribosomal Functions of Bacterial Ribosomal Proteins-An Update, 2023

Ribosomal proteins (r-proteins) are abundant, highly conserved, and multifaceted cellular proteins in all domains of life. Most r-proteins have RNA-binding properties and can form protein-protein contacts. Bacterial r-proteins govern the co-transcriptional rRNA folding during ribosome assembly and participate in the formation of the ribosome functional sites, such as the mRNA-binding site, tRNA-binding sites, the peptidyl transferase center, and the protein exit tunnel. In addition to their primary role in a cell as integral components of the protein synthesis machinery, many r-proteins can function beyond the ribosome (the phenomenon known as moonlighting), acting either as individual regulatory proteins or in complexes with various cellular components. The extraribosomal activities of r-proteins have been studied over the decades. In the past decade, our understanding of r-protein functions has advanced significantly due to intensive studies on ribosomes and gene expression mechanisms not only in model bacteria like Escherichia coli or Bacillus subtilis but also in little-explored bacterial species from various phyla. The aim of this review is to update information on the multiple functions of r-proteins in bacteria.

Mechanistic Analysis of Riboswitch Ligand Interactions Provides Insights into Pharmacological Control over Gene Expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Precise quantification of microRNAs based on proximity ligation of AuNPs-immobilized DNA probes

MiRNAs are critical regulators of target gene expression in many biological processes and are considered promising biomarkers for diseases. In this study, we developed a simple, specific, and sensitive miRNA detection method based on proximity ligation reaction, which is easy to operate. The method uses a pair of target-specific DNA probes immobilized on the same gold nanoparticles (AuNPs), which hybridize to the target miRNA. Hybridization brings the probes close together, allowing the formation of a continuous DNA sequence that can be amplified by Quantitative Real-time PCR (qPCR). This method eliminates the need for complex reverse transcription design and achieves high specificity for discriminating single base mismatches between miRNAs through a simple procedure. This method can sensitively measure three different miRNAs with a detection limit of 20 aM, providing high versatility and sensitivity, even distinguishing single-base variations among members of the miR-200 family with high selectivity. Due to its high selectivity and sensitivity, this method has important implications for the investigation of miRNA biological functions and related biomedical research.

2'-OH as a universal handle for studying intracellular RNAs

RNA plays pivotal roles in most cellular processes, serving as both the traditional carrier of genetic information and as a key regulator of cellular functions. The advent of chemical technologies has contributed critically to the analysis of cellular RNA structures, functions, and interactions. Many of these methods and molecules involve the utilization of chemically reactive handles in RNAs, either introduced externally or inherent within the polymer itself. Among these handles, the 2'-hydroxyl (2'-OH) group has emerged as an exceptionally well-suited and general chemical moiety for the modification and profiling of RNAs in intracellular studies. In this review, we provide an overview of the recent advancements in intracellular applications of acylation at the 2'-OH group of RNA. We outline progress made in probing RNA structure and interactomes, controlling RNA function, RNA imaging, and analyzing RNA-small molecule interactions, all achieved in living cells through this simple chemical handle on the biopolymer.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

Estimating RNA Secondary Structure Folding Free Energy Changes with efn2

A number of analyses require estimates of the folding free energy changes of specific RNA secondary structures. These predictions are often based on a set of nearest neighbor parameters that models the folding stability of a RNA secondary structure as the sum of folding stabilities of the structural elements that comprise the secondary structure. In the software suite RNAstructure, the free energy change calculation is implemented in the program efn2. The efn2 program estimates the folding free energy change and the experimental uncertainty in the folding free energy change. It can be run through the graphical user interface for RNAstructure, from the command line, or a web server. This chapter provides detailed protocols for using efn2.

Efficient Prediction Model of mRNA End-to-End Distance and Conformation: Three-Dimensional RNA Illustration Program (TRIP)

The secondary and tertiary structures of RNA play a vital role in the regulation of biological reactions. These structures have been experimentally studied through in vivo and in vitro analyses, and in silico models have become increasingly accurate in predicting them. Recent technologies have diversified RNA structure predictions, from the earliest thermodynamic and molecular dynamic-based RNA structure predictions to deep learning-based conformation predictions in the past decade. While most research on RNA structure prediction has focused on short non-coding RNAs, there has been limited research on predicting the conformation of longer mRNAs. Our study introduces a computer simulation model called the Three-dimensional RNA Illustration Program (TRIP). TRIP is based on single-chain models and angle restriction of each bead component from previously reported single-molecule fluorescence in situ hybridization (smFISH) experiments. TRIP is a fast and efficient application that only requires up to three inputs to acquire outputs. It can also provide a rough visualization of the 3D conformation of RNA, making it a valuable tool for predicting RNA end-to-end distance.

Targeting non-coding RNAs: Perspectives and challenges of in-silico approaches

The growing information currently available on the central role of non-coding RNAs (ncRNAs) including microRNAs (miRNAS) and long non-coding RNAs (lncRNAs) for chronic and degenerative human diseases makes them attractive therapeutic targets. RNAs carry out different functional roles in human biology and are deeply deregulated in several diseases. So far, different attempts to therapeutically target the 3D RNA structures with small molecules have been reported. In this scenario, the development of computational tools suitable for describing RNA structures and their potential interactions with small molecules is gaining more and more interest. Here, we describe the most suitable strategies to study ncRNAs through computational tools. We focus on methods capable of predicting 2D and 3D ncRNA structures. Furthermore, we describe computational tools to identify, design and optimize small molecule ncRNA binders. This review aims to outline the state of the art and perspectives of computational methods for ncRNAs over the past decade.

Exploring the adaptability of the secondary structure of mRNA to temperature in intertidal snails based on SHAPE experiments

RNA-based thermal regulation is an important strategy for organisms to cope with temperature changes. Inhabiting the intertidal rocky shore, a key interface of the ocean, atmosphere and terrestrial environments, intertidal species have developed variable thermal adaptation mechanisms; however, adaptions at the RNA level remain largely uninvestigated. To examine the relationship between mRNA structural stability and species distribution, in the present study, the secondary structure of cytosolic malate dehydrogenase (cMDH) mRNA of Echinolittorina malaccana, Echinolittorina radiata and Littorina brevicula was determined using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE), and the change in folding free energy of formation (ΔGfold) was calculated. The results showed that ΔGfold increased as the temperature increased. The difference in ΔGfold (ΔΔGfold) between two specific temperatures (25 versus 0°C, 37 versus 0°C and 57 versus 0°C) differed among the three species, and the ΔΔGfold value of E. malaccana was significantly lower than those of E. radiata and L. brevicula. The number of stems of cMDH mRNA of the snails decreased with increasing temperature, and the breakpoint temperature of E. malaccana was the highest among these. The number of loops was also reduced with increasing temperature, while the length of the loop structure increased accordingly. Consequently, these structural changes can potentially affect the translational efficiency of mRNA. These results imply that there were interspecific differences in the thermal stability of RNA secondary structures in intertidal snails, and these differences may be related to snail distribution.

Pervasive transcriptome interactions of protein-targeted drugs

The off-target toxicity of drugs targeted to proteins imparts substantial health and economic costs. Proteome interaction studies can reveal off-target effects with unintended proteins; however, little attention has been paid to intracellular RNAs as potential off-targets that may contribute to toxicity. To begin to assess this, we developed a reactivity-based RNA profiling methodology and applied it to uncover transcriptome interactions of a set of Food and Drug Administration-approved small-molecule drugs in vivo. We show that these protein-targeted drugs pervasively interact with the human transcriptome and can exert unintended biological effects on RNA functions. In addition, we show that many off-target interactions occur at RNA loci associated with protein binding and structural changes, allowing us to generate hypotheses to infer the biological consequences of RNA off-target binding. The results suggest that rigorous characterization of drugs' transcriptome interactions may help assess target specificity and potentially avoid toxicity and clinical failures.

Genome-wide probing of eukaryotic nascent RNA structure elucidates cotranscriptional folding and its antimutagenic effect

The transcriptional intermediates of RNAs fold into secondary structures with multiple regulatory roles, yet the details of such cotranscriptional RNA folding are largely unresolved in eukaryotes. Here, we present eSPET-seq (Structural Probing of Elongating Transcripts in eukaryotes), a method to assess the cotranscriptional RNA folding in Saccharomyces cerevisiae. Our study reveals pervasive structural transitions during cotranscriptional folding and overall structural similarities between nascent and mature RNAs. Furthermore, a combined analysis with genome-wide R-loop and mutation rate approximations provides quantitative evidence for the antimutator effect of nascent RNA folding through competitive inhibition of the R-loops, known to facilitate transcription-associated mutagenesis. Taken together, we present an experimental evaluation of cotranscriptional folding in eukaryotes and demonstrate the antimutator effect of nascent RNA folding. These results suggest genome-wide coupling between the processing and transmission of genetic information through RNA folding.

Mutation signature filtering enables high-fidelity RNA structure probing at all four nucleobases with DMS

Chemical probing experiments have transformed RNA structure analysis, enabling high-throughput measurement of base-pairing in living cells. Dimethyl sulfate (DMS) is one of the most widely used structure probing reagents and has played a pivotal role in enabling next-generation single-molecule probing analyses. However, DMS has traditionally only been able to probe adenine and cytosine nucleobases. We previously showed that, using appropriate conditions, DMS can also be used to interrogate base-pairing of uracil and guanines in vitro at reduced accuracy. However, DMS remained unable to informatively probe guanines in cells. Here, we develop an improved DMS mutational profiling (MaP) strategy that leverages the unique mutational signature of N1-methylguanine DMS modifications to enable high-fidelity structure probing at all four nucleotides, including in cells. Using information theory, we show that four-base DMS reactivities convey greater structural information than current two-base DMS and SHAPE probing strategies. Four-base DMS experiments further enable improved direct base-pair detection by single-molecule PAIR analysis, and ultimately support RNA structure modeling at superior accuracy. Four-base DMS probing experiments are straightforward to perform and will broadly facilitate improved RNA structural analysis in living cells.

Methodologies for bacterial ribonuclease characterization using RNA-seq

Bacteria adjust gene expression at the post-transcriptional level through an intricate network of small regulatory RNAs and RNA-binding proteins, including ribonucleases (RNases). RNases play an essential role in RNA metabolism, regulating RNA stability, decay, and activation. These enzymes exhibit species-specific effects on gene expression, bacterial physiology, and different strategies of target recognition. Recent advances in high-throughput RNA sequencing (RNA-seq) approaches have provided a better understanding of the roles and modes of action of bacterial RNases. Global studies aiming to identify direct targets of RNases have highlighted the diversity of RNase activity and RNA-based mechanisms of gene expression regulation. Here, we review recent RNA-seq approaches used to study bacterial RNases, with a focus on the methods for identifying direct RNase targets.

Pervasive downstream RNA hairpins dynamically dictate start-codon selection

Translational reprogramming allows organisms to adapt to changing conditions. Upstream start codons (uAUGs), which are prevalently present in mRNAs, have crucial roles in regulating translation by providing alternative translation start sites1-4. However, what determines this selective initiation of translation between conditions remains unclear. Here, by integrating transcriptome-wide translational and structural analyses during pattern-triggered immunity in Arabidopsis, we found that transcripts with immune-induced translation are enriched with upstream open reading frames (uORFs). Without infection, these uORFs are selectively translated owing to hairpins immediately downstream of uAUGs, presumably by slowing and engaging the scanning preinitiation complex. Modelling using deep learning provides unbiased support for these recognizable double-stranded RNA structures downstream of uAUGs (which we term uAUG-ds) being responsible for the selective translation of uAUGs, and allows the prediction and rational design of translating uAUG-ds. We found that uAUG-ds-mediated regulation can be generalized to human cells. Moreover, uAUG-ds-mediated start-codon selection is dynamically regulated. After immune challenge in plants, induced RNA helicases that are homologous to Ded1p in yeast and DDX3X in humans resolve these structures, allowing ribosomes to bypass uAUGs to translate downstream defence proteins. This study shows that mRNA structures dynamically regulate start-codon selection. The prevalence of this RNA structural feature and the conservation of RNA helicases across kingdoms suggest that mRNA structural remodelling is a general feature of translational reprogramming.

The RNA interactome in the Hallmarks of Cancer

Ribonucleic acid (RNA) molecules are indispensable for cellular homeostasis in healthy and malignant cells. However, the functions of RNA extend well beyond that of a protein-coding template. Rather, both coding and non-coding RNA molecules function through critical interactions with a plethora of cellular molecules, including other RNAs, DNA, and proteins. Deconvoluting this RNA interactome, including the interacting partners, the nature of the interaction, and dynamic changes of these interactions in malignancies has yielded fundamental advances in knowledge and are emerging as a novel therapeutic strategy in cancer. Here, we present an RNA-centric review of recent advances in the field of RNA-RNA, RNA-protein, and RNA-DNA interactomic network analysis and their impact across the Hallmarks of Cancer. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

RNA structure determination: From 2D to 3D

RNA molecules serve a wide range of functions that are closely linked to their structures. The basic structural units of RNA consist of single- and double-stranded regions. In order to carry out advanced functions such as catalysis and ligand binding, certain types of RNAs can adopt higher-order structures. The analysis of RNA structures has progressed alongside advancements in structural biology techniques, but it comes with its own set of challenges and corresponding solutions. In this review, we will discuss recent advances in RNA structure analysis techniques, including structural probing methods, X-ray crystallography, nuclear magnetic resonance, cryo-electron microscopy, and small-angle X-ray scattering. Often, a combination of multiple techniques is employed for the integrated analysis of RNA structures. We also survey important RNA structures that have been recently determined using various techniques.

RNA Structural Dynamics Modulate EGFR-TKI Resistance Through Controlling YRDC Translation in NSCLC Cells

Epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKIs) positively affect the initial control of non-small cell lung cancer (NSCLC). Rapidly acquired resistance to EGFR-TKIs is a major hurdle in successful treatment. However, the mechanisms that control the resistance of EGFR-TKIs remain largely unknown. RNA structures have widespread and crucial functions in many biological regulations; however, the functions of RNA structures in regulating cancer drug resistance remain unclear. Here, the psoralen analysis of RNA interactions and structures (PARIS) method is used to establish the higher-order RNA structure maps of EGFR-TKIs-resistant and -sensitive cells of NSCLC. Our results show that RNA structural regions are enriched in untranslated regions (UTRs) and correlate with translation efficiency (TE). Moreover, yrdC N6-threonylcarbamoyltransferase domain containing (YRDC) promotes resistance to EGFR-TKIs. RNA structure formation in YRDC 3' UTR suppresses embryonic lethal abnormal vision-like 1 (ELAVL1) binding, leading to EGFR-TKI sensitivity by impairing YRDC translation. A potential therapeutic strategy for cancer treatment is provided using antisense oligonucleotide (ASO) to perturb the interaction between RNA and protein. Our study reveals an unprecedented mechanism through which the RNA structure switch modulates EGFR-TKI resistance by controlling YRDC mRNA translation in an ELAVL1-dependent manner.

End-to-end RT-PCR of long RNA and highly structured RNA

RNA molecules play important roles in numerous normal cellular processes and disease states, from protein coding to gene regulation. RT-PCR, applying the power of polymerase chain reaction (PCR) to RNA by coupling reverse transcription with PCR, is one of the most important techniques to characterize RNA transcripts and monitor gene expression. The ability to analyze full-length RNA transcripts and detect their expression is critical to decipher their biological functions. However, due to the low processivity of retroviral reverse transcriptases (RTs), we can only monitor a small fraction of long RNA transcripts, especially those containing stable secondary and tertiary structures. The full-length sequences can only be deduced by computational analysis, which is often misleading. Group II intron-encoded RTs are a new type of RT enzymes. They have evolved specialized structural elements that unwind template structures and maintain close contact with the RNA template. Therefore, group II intron-encoded RTs are more processive than the retroviral RTs. The discovery, optimization and deployment of processive group II intron RTs provide us the opportunity to analyze RNA transcripts with single molecule resolution. MarathonRT, the most processive group II intron RT, has been extensively optimized for processive reverse transcription. In this chapter, we use MarathonRT to deliver a general protocol for long amplicon generation by RT-PCR, and also provide guidance for troubleshooting and further optimization.

Amplifying gene expression with RNA-targeted therapeutics

Many diseases are caused by insufficient expression of mutated genes and would benefit from increased expression of the corresponding protein. However, in drug development, it has been historically easier to develop drugs with inhibitory or antagonistic effects. Protein replacement and gene therapy can achieve the goal of increased protein expression but have limitations. Recent discoveries of the extensive regulatory networks formed by non-coding RNAs offer alternative targets and strategies to amplify the production of a specific protein. In addition to RNA-targeting small molecules, new nucleic acid-based therapeutic modalities that allow highly specific modulation of RNA-based regulatory networks are being developed. Such approaches can directly target the stability of mRNAs or modulate non-coding RNA-mediated regulation of transcription and translation. This Review highlights emerging RNA-targeted therapeutics for gene activation, focusing on opportunities and challenges for translation to the clinic.

Led-Seq: ligation-enhanced double-end sequence-based structure analysis of RNA

Structural analysis of RNA is an important and versatile tool to investigate the function of this type of molecules in the cell as well as in vitro. Several robust and reliable procedures are available, relying on chemical modification inducing RT stops or nucleotide misincorporations during reverse transcription. Others are based on cleavage reactions and RT stop signals. However, these methods address only one side of the RT stop or misincorporation position. Here, we describe Led-Seq, a new approach based on lead-induced cleavage of unpaired RNA positions, where both resulting cleavage products are investigated. The RNA fragments carrying 2', 3'-cyclic phosphate or 5'-OH ends are selectively ligated to oligonucleotide adapters by specific RNA ligases. In a deep sequencing analysis, the cleavage sites are identified as ligation positions, avoiding possible false positive signals based on premature RT stops. With a benchmark set of transcripts in Escherichia coli, we show that Led-Seq is an improved and reliable approach based on metal ion-induced phosphodiester hydrolysis to investigate RNA structures in vivo.

A viral pan-end RNA element and host complex define a SARS-CoV-2 regulon

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, generates multiple protein-coding, subgenomic RNAs (sgRNAs) from a longer genomic RNA, all bearing identical termini with poorly understood roles in regulating viral gene expression. Insulin and interferon-gamma, two host-derived, stress-related agents, and virus spike protein, induce binding of glutamyl-prolyl-tRNA synthetase (EPRS1), within an unconventional, tetra-aminoacyl-tRNA synthetase complex, to the sgRNA 3'-end thereby enhancing sgRNA expression. We identify an EPRS1-binding sarbecoviral pan-end activating RNA (SPEAR) element in the 3'-end of viral RNAs driving agonist-induction. Translation of another co-terminal 3'-end feature, ORF10, is necessary for SPEAR-mediated induction, independent of Orf10 protein expression. The SPEAR element enhances viral programmed ribosomal frameshifting, thereby expanding its functionality. By co-opting noncanonical activities of a family of essential host proteins, the virus establishes a post-transcriptional regulon stimulating global viral RNA translation. A SPEAR-targeting strategy markedly reduces SARS-CoV-2 titer, suggesting a pan-sarbecoviral therapeutic modality.

Roles of the leader-trailer helix and antitermination complex in biogenesis of the 30S ribosomal subunit

Ribosome biogenesis occurs co-transcriptionally and entails rRNA folding, ribosomal protein binding, rRNA processing, and rRNA modification. In most bacteria, the 16S, 23S and 5S rRNAs are co-transcribed, often with one or more tRNAs. Transcription involves a modified RNA polymerase, called the antitermination complex, which forms in response to cis-acting elements (boxB, boxA and boxC) in the nascent pre-rRNA. Sequences flanking the rRNAs are complementary and form long helices known as leader-trailer helices. Here, we employed an orthogonal translation system to interrogate the functional roles of these RNA elements in 30S subunit biogenesis in Escherichia coli. Mutations that disrupt the leader-trailer helix caused complete loss of translation activity, indicating that this helix is absolutely essential for active subunit formation in the cell. Mutations of boxA also reduced translation activity, but by only 2- to 3-fold, suggesting a smaller role for the antitermination complex. Similarly modest drops in activity were seen upon deletion of either or both of two leader helices, termed here hA and hB. Interestingly, subunits formed in the absence of these leader features exhibited defects in translational fidelity. These data suggest that the antitermination complex and precursor RNA elements help to ensure quality control during ribosome biogenesis.

Initiator AUGs Are Discriminated from Elongator AUGs Predominantly through mRNA Accessibility in C. crescentus

The initiation of translation in bacteria is thought to occur upon base pairing between the Shine-Dalgarno (SD) site in the mRNA and the anti-SD site in the rRNA. However, in many bacterial species, such as Caulobacter crescentus, a minority of mRNAs have SD sites. To examine the functional importance of SD sites in C. crescentus, we analyzed the transcriptome and found that more SD sites exist in the coding sequence than in the preceding start codons. To examine the function of SD sites in initiation, we designed a series of mutants with altered ribosome accessibility and SD content in translation initiation regions (TIRs) and in elongator AUG regions (EARs). A lack of mRNA structure content is required for initiation in TIRs, and, when introduced into EARs, can stimulate initiation, thereby suggesting that low mRNA structure content is a major feature that is required for initiation. SD sites appear to stimulate initiation in TIRs, which generally lack structure content, but SD sites only stimulate initiation in EARs if RNA secondary structures are destabilized. Taken together, these results suggest that the difference in secondary structure between TIRs and EARs directs ribosomes to start codons where SD base pairing can tune the efficiency of initiation, but SDs in EARs do not stimulate initiation, as they are blocked by stable secondary structures. This highlights the importance of studying translation initiation mechanisms in diverse bacterial species. IMPORTANCE Start codon selection is an essential process that is thought to occur via the base pairing of the rRNA to the SD site in the mRNA. This model is based on studies in E. coli, yet whole-genome sequencing revealed that SD sites are absent at start codons in many species. By examining the transcriptome of C. crescentus, we found more SD-AUG pairs in the CDS of mRNAs than preceding start codons, yet these internal sites do not initiate. Instead, start codon regions have lower mRNA secondary structure content than do internal SD-AUG regions. Therefore, we find that start codon selection is not controlled by the presence of SD sites, which tune initiation efficiency, but by lower mRNA structure content surrounding the start codon.

Structural analysis of MALAT1 long noncoding RNA in cells and in evolution

Although not canonically polyadenylated, the long noncoding RNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) is stabilized by a highly conserved 76-nt triple helix structure on its 3' end. The entire MALAT1 transcript is over 8000 nt long in humans. The strongest structural conservation signal in MALAT1 (as measured by covariation of base pairs) is in the triple helix structure. Primary sequence analysis of covariation alone does not reveal the degree of structural conservation of the entire full-length transcript, however. Furthermore, RNA structure is often context dependent; RNA binding proteins that are differentially expressed in different cell types may alter structure. We investigate here the in-cell and cell-free structures of the full-length human and green monkey (Chlorocebus sabaeus) MALAT1 transcripts in multiple tissue-derived cell lines using SHAPE chemical probing. Our data reveal levels of uniform structural conservation in different cell lines, in cells and cell-free, and even between species, despite significant differences in primary sequence. The uniformity of the structural conservation across the entire transcript suggests that, despite seeing covariation signals only in the triple helix junction of the lncRNA, the rest of the transcript's structure is remarkably conserved, at least in primates and across multiple cell types and conditions.

Mutation signature filtering enables high-fidelity RNA structure probing at all four nucleobases with DMS

Chemical probing experiments have transformed RNA structure analysis, enabling high-throughput measurement of base-pairing in living cells. Dimethyl sulfate (DMS) is one of the most widely used structure probing reagents and has played a prominent role in enabling next-generation single-molecule probing analyses. However, DMS has traditionally only been able to probe adenine and cytosine nucleobases. We previously showed that, using appropriate conditions, DMS can also be used to interrogate base-pairing of uracil and guanines in vitro at reduced accuracy. However, DMS remained unable to informatively probe guanines in cells. Here, we develop an improved DMS mutational profiling (MaP) strategy that leverages the unique mutational signature of N 1 -methylguanine DMS modifications to enable robust, high-fidelity structure probing at all four nucleotides, including in cells. Using information theory, we show that four-base DMS reactivities convey greater structural information than comparable two-base DMS and SHAPE probing strategies. Four-base DMS experiments further enable improved direct base-pair detection by single-molecule PAIR analysis, and ultimately support RNA structure modeling at superior accuracy. Four-base DMS probing experiments are easily performed and will broadly facilitate improved RNA structural analysis in living cells.

Structure-first identification of RNA elements that regulate dengue virus genome architecture and replication

The genomes of RNA viruses encode the information required for replication in host cells both in their linear sequence and in complex higher-order structures. A subset of these RNA genome structures show clear sequence conservation, and have been extensively described for well-characterized viruses. However, the extent to which viral RNA genomes contain functional structural elements-unable to be detected by sequence alone-that nonetheless are critical to viral fitness is largely unknown. Here, we devise a structure-first experimental strategy and use it to identify 22 structure-similar motifs across the coding sequences of the RNA genomes for the four dengue virus serotypes. At least 10 of these motifs modulate viral fitness, revealing a significant unnoticed extent of RNA structure-mediated regulation within viral coding sequences. These viral RNA structures promote a compact global genome architecture, interact with proteins, and regulate the viral replication cycle. These motifs are also thus constrained at the levels of both RNA structure and protein sequence and are potential resistance-refractory targets for antivirals and live-attenuated vaccines. Structure-first identification of conserved RNA structure enables efficient discovery of pervasive RNA-mediated regulation in viral genomes and, likely, other cellular RNAs.

Sulfonylation of RNA 2'-OH groups

The nucleophilic reactivity of RNA 2'-OH groups in water has proven broadly useful in probing, labeling, and conjugating RNA. To date, reactions selective to ribose 2'-OH have been limited to bond formation with short-lived carbonyl electrophiles. Here we report that many activated small-molecule sulfonyl species can exhibit extended lifetimes in water and retain 2'-OH reactivity. The data establish favorable aqueous solubility for selected reagents and successful RNA-selective reactions at stoichiometric and superstoichiometric yields, particularly for aryl sulfonyltriazole species. We report that the latter are considerably more stable than most prior carbon electrophiles in aqueous environments and tolerate silica chromatography. Furthermore, an azide-substituted sulfonyltriazole reagent is developed to introduce labels into RNA via click chemistry. In addition to high-yield reactions, we find that RNA sulfonylation can also be performed under conditions that give trace yields necessary for structure mapping. Like acylation, the reaction occurs with selectivity for unpaired nucleotides over those in the duplex structure, and a sulfonate adduct causes reverse transcriptase stops, suggesting potential use in RNA structure analysis. Probing of rRNA is demonstrated in human cells, indicating possible cell permeability. The sulfonyl reagent class enables a new level of control, selectivity, versatility, and ease of preparation for RNA applications.

Machine Learning Informs RNA-Binding Chemical Space

Small molecule targeting of RNA has emerged as a new frontier in medicinal chemistry, but compared to the protein targeting literature our understanding of chemical matter that binds to RNA is limited. In this study, we reported Repository Of BInders to Nucleic acids (ROBIN), a new library of nucleic acid binders identified by small molecule microarray (SMM) screening. The complete results of 36 individual nucleic acid SMM screens against a library of 24 572 small molecules were reported (including a total of 1 627 072 interactions assayed). A set of 2 003 RNA-binding small molecules was identified, representing the largest fully public, experimentally derived library of its kind to date. Machine learning was used to develop highly predictive and interpretable models to characterize RNA-binding molecules. This work demonstrates that machine learning algorithms applied to experimentally derived sets of RNA binders are a powerful method to inform RNA-targeted chemical space.

Probing the dynamic RNA structurome and its functions

RNA is a key regulator of almost every cellular process, and the structures adopted by RNA molecules are thought to be central to their functions. The recent fast-paced evolution of high-throughput sequencing-based RNA structure mapping methods has enabled the rapid in vivo structural interrogation of entire cellular transcriptomes. Collectively, these studies are shedding new light on the long underestimated complexity of the structural organization of the transcriptome - the RNA structurome. Moreover, recent analyses are challenging the view that the RNA structurome is a static entity by revealing how RNA molecules establish intricate networks of alternative intramolecular and intermolecular interactions and that these ensembles of RNA structures are dynamically regulated to finely tune RNA functions in living cells. This new understanding of how RNA can shape cell phenotypes has important implications for the development of RNA-targeted therapeutic strategies.

Chemical RNA Cross-Linking: Mechanisms, Computational Analysis, and Biological Applications

In recent years, RNA has emerged as a multifaceted biomolecule that is involved in virtually every function of the cell and is critical for human health. This has led to a substantial increase in research efforts to uncover the many chemical and biological aspects of RNA and target RNA for therapeutic purposes. In particular, analysis of RNA structures and interactions in cells has been critical for understanding their diverse functions and druggability. In the last 5 years, several chemical methods have been developed to achieve this goal, using chemical cross-linking combined with high-throughput sequencing and computational analysis. Applications of these methods resulted in important new insights into RNA functions in a variety of biological contexts. Given the rapid development of new chemical technologies, a thorough perspective on the past and future of this field is provided. In particular, the various RNA cross-linkers and their mechanisms, the computational analysis and challenges, and illustrative examples from recent literature are discussed.

RNA 3D Structure Comparison Using RNA-Puzzles Toolkit

Computational modeling of RNA three-dimensional (3D) structure may help in unrevealing the molecular mechanisms of RNA molecules and in designing molecules with novel functions. An unbiased blind assessment to benchmark the computational modeling is required to understand the achievements and bottlenecks of the prediction, while a standard structure comparison protocol is necessary. RNA-Puzzles is a community-wide effort on the assessment of blind prediction of RNA tertiary structures. And RNA-Puzzles toolkit is a computational resource derived from RNA-Puzzles, which includes (i) decoy sets generated by different RNA 3D structure prediction methods; (ii) 3D structure normalization, analysis, manipulation, and visualization tools; and (iii) 3D structure comparison metric tools. In this chapter, we illustrate a standard RNA 3D structure prediction assessment protocol using the selected tools from RNA-Puzzles toolkit: rna-tools and RNA_assessment.

Natural selection pressure exerted on "Silent" mutations during the evolution of SARS-CoV-2: Evidence from codon usage and RNA structure

From the first emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) till now, multiple mutations that caused synonymous and nonsynonymous substitutions have accumulated. Among them, synonymous substitutions were regarded as "silent" mutations that received less attention than nonsynonymous substitutions that cause amino acid variations. However, the importance of synonymous substitutions can not be neglected. This research focuses on synonymous substitutions on SARS-CoV-2 and proves that synonymous substitutions were under purifying selection in its evolution. The evidence of purifying selection is provided by comparing the mutation number per site in coding and non-coding regions. We then study the two forces of purifying selection: synonymous codon usage and RNA secondary structure. Results show that the codon usage optimization leads to an adapted codon usage towards humans. Furthermore, our results show that the maintenance of RNA secondary structure causes the purifying of synonymous substitutions in the structural region. These results explain the selection pressure on synonymous substitutions during the evolution of SARS-CoV-2.

Chemical Probing of RNA Structure In Vivo Using SHAPE-MaP and DMS-MaP

The functional roles of RNAs are often regulated by their structure. Selective 2' hydroxyl acylation analyzed by primer extension (SHAPE) and dimethyl sulfate (DMS) base reactivity can be employed to investigate the flexibility of nucleotides and correlate it to the constraints imparted by base-pairing and/or protein-binding. In vivo, a multitude of proteins could bind an RNA molecule, regulating its structure and function. Hence, to obtain a more comprehensive view of the RNA structure-function relationship in vivo, it may be required to characterize both the RNA structure and the RNA-protein interaction network. In this chapter, we describe methods for characterizing the in vivo nucleotide flexibility of RNA in cells by SHAPE-MaP (SHAPE by Mutational Profiling) and DMS-MaP. In another chapter, we will discuss the characterization of RNA-protein interaction network by RNP-MaP.

The Influence of the Nucleotide Composition of Genes and Gene Regulatory Elements on the Efficiency of Protein Expression in Escherichia coli

Recombinant proteins expressed in Escherichia coli are widely used in biochemical research and industrial processes. At the same time, achieving higher protein expression levels and correct protein folding still remains the key problem, since optimization of nutrient media, growth conditions, and methods for induction of protein synthesis do not always lead to the desired result. Often, low protein expression is determined by the sequences of the expressed genes and their regulatory regions. The genetic code is degenerated; 18 out of 20 amino acids are encoded by more than one codon. Choosing between synonymous codons in the coding sequence can significantly affect the level of protein expression and protein folding due to the influence of the gene nucleotide composition on the probability of formation of secondary mRNA structures that affect the ribosome binding at the translation initiation phase, as well as the ribosome movement along the mRNA during elongation, which, in turn, influences the mRNA degradation and the folding of the nascent protein. The nucleotide composition of the mRNA untranslated regions, in particular the promoter and Shine-Dalgarno sequences, also affects the efficiency of mRNA transcription, translation, and degradation. In this review, we describe the genetic principles that determine the efficiency of protein production in Escherichia coli.