RNA Secondary Structure Research Articles

The ALS/FTD-related C9orf72 hexanucleotide repeat expansion forms RNA condensates through multimolecular G-quadruplexes

Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are neurodegenerative diseases that exist on a clinico-pathogenetic spectrum, designated ALS/FTD. The most common genetic cause of ALS/FTD is expansion of the intronic hexanucleotide repeat (GGGGCC)n in C9orf72. Here, we investigate the formation of nucleic acid secondary structures in these expansion repeats, and their role in generating condensates characteristic of ALS/FTD. We observe significant aggregation of the hexanucleotide sequence (GGGGCC)n, which we associate to the formation of multimolecular G-quadruplexes (mG4s) by using a range of biophysical techniques. Exposing the condensates to G4-unfolding conditions leads to prompt disassembly, highlighting the key role of mG4-formation in the condensation process. We further validate the biological relevance of our findings by detecting an increased prevalence of G4-structures in C9orf72 mutant human motor neurons when compared to healthy motor neurons by staining with a G4-selective fluorescent probe, revealing signal in putative condensates. Our findings strongly suggest that RNA G-rich repetitive sequences can form protein-free condensates sustained by multimolecular G-quadruplexes, highlighting their potential relevance as therapeutic targets for C9orf72 mutation-related ALS/FTD.

The distinct translational landscapes of gram-negative Salmonella and gram-positive Listeria

Translational control in pathogenic bacteria is fundamental to gene expression and affects virulence and other infection phenotypes. We used an enhanced ribosome profiling protocol coupled with parallel transcriptomics to capture accurately the global translatome of two evolutionarily distant pathogenic bacteria—the Gram-negative bacterium Salmonella and the Gram-positive bacterium Listeria. We find that the two bacteria use different mechanisms to translationally regulate protein synthesis. In Salmonella, in addition to the expected correlation between translational efficiency and cis-regulatory features such as Shine–Dalgarno (SD) strength and RNA secondary structure around the initiation codon, our data reveal an effect of the 2nd and 3rd codons, where the presence of tandem lysine codons (AAA-AAA) enhances translation in both Salmonella and E. coli. Strikingly, none of these features are seen in efficiently translated Listeria transcripts. Instead, approximately 20% of efficiently translated Listeria genes exhibit 70 S footprints seven nt upstream of the authentic start codon, suggesting that these genes may be subject to a novel translational initiation mechanism. Our results show that SD strength is not a direct hallmark of translational efficiency in all bacteria. Instead, Listeria has evolved additional mechanisms to control gene expression level that are distinct from those utilised by Salmonella and E. coli.

Dissecting the role of GGGGCC repeat RNA secondary structure in cellular toxicity associated with C9 FTD/ALS

AbstractBackgroundAberrant expansions of simple nucleotide repeat sequences can cause more than fifty human diseases. For example, the most common genetic cause of Frontotemporal dementia (FTD) and Amyotrophic Lateral Sclerosis (ALS) is the expansion of an intronic GGGGCC hexanucleotide repeat sequence in C9orf72 gene (C9 FTD/ALS). Transcribed GGGGCC repeat RNAs can form complex secondary structures, which facilitate unusual interactions with cellular RNA‐binding proteins (RBPs). Sequestration of crucial RBPs into nuclear foci can inhibit the normal functions of those proteins. In addition, expanded GGGGCC repeat RNAs can support the translation of highly toxic repetitive peptides through a non‐canonical process known as repeat associated non‐AUG translation (RAN translation), which may require some specialized factors. Here we sought to understand how RNA secondary structure influence GGGGCC RAN translation and repeat elicited toxicity in C9 FTD/ALS.MethodWe used chemical modification coupled with traditional high‐throughput sequencing and nano‐pore sequencing approaches to probe GGGGCC repeat RNA secondary structures inside cells. We further use in vitro and in vivo RAN translation assays to determine whether repeat RNA secondary structures influence RAN translation efficiency.ResultOur RNA structural studies revealed that GGGGCC repeat RNAs form primarily form hairpin structures in cells. With nanopore sequencing, we were also able to detect the possible presence of unfolded and/or g‐quadruplex folded GGGGCC repeat RNAs in cells. Folding GGGGCC RNAs into hairpins supports RAN translation whereas folding repeat RNAs into g‐quadruplex structures strongly inhibits RAN translation in both in vitro and in vivo assays. Our results suggest that RNA secondary structure can affect the RAN translation elongation efficiency and g‐quadruplex structures may lead to the termination of elongating ribosomes.ConclusionWe present a deep sequencing‐based chemical probing technique to evaluate the roles of repeat RNA secondary structures in RNA toxicity underlying dementia‐associated repeat expansion disorders. Repeat RNA structural information can be used to identify chemical compounds that can bind to the RNA and prevent cytotoxicity in cells.

Automated design of dynamic programming schemes for RNA folding with pseudoknots

Although RNA secondary structure prediction is a textbook application of dynamic programming (DP) and routine task in RNA structure analysis, it remains challenging whenever pseudoknots come into play. Since the prediction of pseudoknotted structures by minimizing (realistically modelled) energy is NP-hard, specialized algorithms have been proposed for restricted conformation classes that capture the most frequently observed configurations. To achieve good performance, these methods rely on specific and carefully hand-crafted DP schemes. In contrast, we generalize and fully automatize the design of DP pseudoknot prediction algorithms. For this purpose, we formalize the problem of designing DP algorithms for an (infinite) class of conformations, modeled by (a finite number of) fatgraphs, and automatically build DP schemes minimizing their algorithmic complexity. We propose an algorithm for the problem, based on the tree-decomposition of a well-chosen representative structure, which we simplify and reinterpret as a DP scheme. The algorithm is fixed-parameter tractable for the treewidth tw of the fatgraph, and its output represents a O(ntw+1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathcal {O}}(n^{tw+1})$$\\end{document} algorithm (and even possibly O(ntw)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathcal {O}}(n^{tw})$$\\end{document} in simple energy models) for predicting the MFE folding of an RNA of length n. We demonstrate, for the most common pseudoknot classes, that our automatically generated algorithms achieve the same complexities as reported in the literature for hand-crafted schemes. Our framework supports general energy models, partition function computations, recursive substructures and partial folding, and could pave the way for algebraic dynamic programming beyond the context-free case.

RNAknot2S: a New RNA Secondary Structure Prediction Method that Combines the Comparative Approach with Genetic Algorithms

RNAknot2S: a New RNA Secondary Structure Prediction Method that Combines the Comparative Approach with Genetic Algorithms

Bioactive nutraceuticals as G4 stabilizers: potential cancer prevention and therapy-a critical review.

G-quadruplexes (G4) are non-canonical, four-stranded, nucleic acid secondary structures formed in the guanine-rich sequences, where guanine nucleotides associate with each other via Hoogsteen hydrogen bonding. These structures are widely found near the functional regions of the mammalian genome, such as telomeres, oncogenic promoters, and replication origins, and play crucial regulatory roles in replication and transcription. Destabilization of G4 by various carcinogenic agents allows oncogene overexpression and extension of telomeric ends resulting in dysregulation of cellular growth-promoting oncogenesis. Therefore, targeting and stabilizing these G4 structures with potential ligands could aid cancer prevention and therapy. The field of G-quadruplex targeting is relatively nascent, although many articles have demonstrated the effect of G4 stabilization on oncogenic expressions; however, no previous study has provided a comprehensive analysis about the potency of a wide variety of nutraceuticals and some of their derivatives in targeting G4 and the lattice of oncogenic cell signaling cascade affected by them. In this review, we have discussed bioactive G4-stabilizing nutraceuticals, their sources, mode of action, and their influence on cellular signaling, and we believe our insight would bring new light to the current status of the field and motivate researchers to explore this relatively poorly studied arena.

Role of DEAD-box RNA helicases in low-temperature adapted growth of Antarctic Pseudomonas syringae Lz4W.

RNA metabolism is important as RNA acts as a link between genomic information and functional biomolecules, thereby playing a critical role in cellular response to environment. We investigated the role of DEAD-box RNA helicases in low-temperature adapted growth of P. syringae, as this group of enzymes play an essential role in modulation of RNA secondary structures. This is the first report on the assessment of all major DEAD-box RNA helicases in any Antarctic bacterium. Of the five RNA helicases, three (srmB, csdA, and dbpA) are important for the growth of the Antarctic P. syringae at low temperature. However, the requisite role of dbpA and the indispensable requirement of csdA for low-temperature adapted growth are a novel finding of this study. Growth analysis of combinatorial deletion strains was performed to understand the functional interaction among helicase genes. Similarly, genetic complementation of RNA helicase mutants was conducted for identification of gene redundancy in P. syringae.

The Secondary Structure of Potato Spindle Tuber Viroid Determines Its Infectivity in Nicotiana benthamiana.

The function of RNAs is determined by their structure. However, studying the relationship between RNA structure and function often requires altering RNA sequences to modify the structures, which leads to the neglect of the importance of RNA sequences themselves. In our research, we utilized potato spindle tuber viroid (PSTVd), a circular-form non-coding infectious RNA, as a model with which to investigate the role of a specific rod-like structure in RNA function. By generating linear RNA transcripts with different start sites, we established 12 PSTVd forms with different secondary structures while maintaining the same sequence. The RNA secondary structures were predicted using the mfold tool and validated through native PAGE gel electrophoresis after in vitro RNA folding. Analysis using plant infection assays revealed that the formation of a correct rod-like structure is crucial for the successful infection of PSTVd. Interestingly, the inability of PSTVd forms with non-rod-like structures to infect plants could be partially compensated by increasing the amount of linear viroid RNA transcripts, suggesting the existence of additional RNA secondary structures, such as the correct rod-like structure, alongside the dominant structure in the RNA inoculum of these forms. Our study demonstrates the critical role of RNA secondary structures in determining the function of infectious RNAs.

N-methyl mesoporphyrin IX (NMM) as electrochemical probe for detection of guanine quadruplexes

G-quadruplexes (G4) are stable alternative secondary structures of nucleic acids. With increasing understanding of their roles in biological processes and their application in bio- and nanotechnology, the exploration of novel methods for the analysis of these structures is becoming important. In this work, N-methyl mesoporphyrin IX (NMM) was used as a voltammetric probe for an easy electrochemical detection of G4s. Cyclic voltammetry on a hanging mercury drop electrode (HMDE) was used to detect NMM with a limit of detection (LOD) of 40 nM. Characteristic reduction signal of NMM was found to be substantially higher in the presence of G4 oligodeoxynucleotides (ODNs) than in the presence of single- or double-stranded ODNs and even ODNs susceptible to form G4s but in their unfolded, single-stranded forms. Gradual transition from unstructured single strand to G4, induced by increasing concentrations of the G4 stabilizing K+ ions, was detected by an electrochemical method for the first time. All obtained results were supported by circular dichroism spectroscopy. This work expands on the concept of electrochemical probes utilization in DNA secondary structure recognition and offers a proof of principle that can be potentially employed in the development of novel electroanalytical methods for nucleic acid structure studies.

Association between leptin receptor polymorphisms and polycystic ovary syndrome risk: a meta-analysis based on 11 studies

Objective Published evidence indicated that the leptin receptor (LEPR) gene polymorphisms are associated with polycystic ovary syndrome (PCOS) risk. However, studies on the association between the polymorphisms of LEPR gene are inconsistent or even controversial. Material and Methods We conducted this meta-analysis to explore the more precise relationship between LEPR polymorphisms and PCOS risk. Relevant articles were searched with five online databases up to March 1 2023. Odds ratios (OR) with 95% confidence intervals (CI) were selected to examine the statistical strength of each genetic model. Moreover, RNA secondary structure and variant effects of these loci were examined with in silico analysis. Results Overall, 11 publications were analyzed, and the pooled results did not present any significant association between rs1137101 A/G polymorphism and PCOS risk in general population and some subgroup analysis. But the significant association were observed in Asian population (AG vs. AA: OR = 0.51, 95%CI = 0.32–0.81, p = .01, I2=0%; AG + GG vs. AA: OR = 0.41, 95%CI = 0.26–0.65, p < .01, I2=25.9%). Moreover, similar positive associations were also observed in rs1805096 polymorphism with PCOS risk. Conclusion In summary, our meta-analysis suggested that the LEPR gene polymorphisms might be associated with PCOS susceptibility. Owing to the limited studies and small sample size in our meta-analysis, more well-designed studies from different races were needed to be conducted to verify the current results.

Abstract 16901: A Computational Framework to Estimate the Consequences of Cardiac RNA Editing Events at Single-Cell Resolution

Introduction: RNA is a highly dynamic molecule subject to various epigenetic modifications. One such process is RNA editing, a tightly regulated mechanism that introduces alterations to RNA nucleotides without affecting the DNA template. RNA editing serves as a mechanism for expanding and diversifying the functional repertoire of proteins. In the context of the heart, the RNA editors assume crucial roles in preserving cardiac function. Despite the importance of RNA editing, there are no framework to evaluate consequence of synonymous RNA editing events. Synonymous RNA editing events can affect mRNA stability, translation efficiency, protein folding, and even protein-protein interactions. These subtle changes in the coding region can influence gene expression, protein synthesis, and ultimately, cellular function. Research question: What are the functional consequences of RNA editing events in different types of cardiac cells? Methods: To estimate the RNA structural metrics, we employed the Vienna RNA package software. This comprehensive package provides a range of metrics for quantifying RNA secondary structures. Results: We have developed a framework to profile RNA editing and analyse synonymous events using new and existing single-cell datasets. Our framework enables users to explore the distribution of RNA editing events across different cell types and states, as well as analyse the functional consequences of these events. Here, to demonstrate its efficacy, we have applied our framework on healthy hearts. We identified highly edited cell types such as cardiomyocytes, fibroblasts, endothelial cells, macrophages, and pericytes, with over 1,000 edited sites. Interestingly, more than 50% of these edited sites enhanced secondary RNA structures, while 25% of the events destabilized the secondary structure. Conclusions: We have developed a comprehensive framework to estimate consequence of RNA editing events at single-cell resolution. Our results demonstrate the power of our framework in elucidating the impact of RNA editing on cellular processes. Synonymous RNA editing expands our understanding of the complexity of genetic regulation and its implications for cellular processes and disease.

An exploratory in silico comparison of open-source codon harmonization tools

BackgroundNot changing the native constitution of genes prior to their expression by a heterologous host can affect the amount of proteins synthesized as well as their folding, hampering their activity and even cell viability. Over the past decades, several strategies have been developed to optimize the translation of heterologous genes by accommodating the difference in codon usage between species. While there have been a handful of studies assessing various codon optimization strategies, to the best of our knowledge, no research has been performed towards the evaluation and comparison of codon harmonization algorithms. To highlight their importance and encourage meaningful discussion, we compared different open-source codon harmonization tools pertaining to their in silico performance, and we investigated the influence of different gene-specific factors.ResultsIn total, 27 genes were harmonized with four tools toward two different heterologous hosts. The difference in %MinMax values between the harmonized and the original sequences was calculated (ΔMinMax), and statistical analysis of the obtained results was carried out. It became clear that not all tools perform similarly, and the choice of tool should depend on the intended application. Almost all biological factors under investigation (GC content, RNA secondary structures and choice of heterologous host) had a significant influence on the harmonization results and thus must be taken into account. These findings were substantiated using a validation dataset consisting of 8 strategically chosen genes.ConclusionsDue to the size of the dataset, no complex models could be developed. However, this initial study showcases significant differences between the results of various codon harmonization tools. Although more elaborate investigation is needed, it is clear that biological factors such as GC content, RNA secondary structures and heterologous hosts must be taken into account when selecting the codon harmonization tool.

G-Quadruplexes in the Viral Genome: Unlocking Targets for Therapeutic Interventions and Antiviral Strategies.

G-quadruplexes (G4s) are unique non-canonical four-stranded nucleic acid secondary structures formed by guanine-rich DNA or RNA sequences. Sequences with the potential to form quadruplex motifs (pG4s) are prevalent throughout the genomes of all organisms, spanning from prokaryotes to eukaryotes, and are enriched within regions of biological significance. In the past few years, the identification of pG4s within most of the Baltimore group viruses has attracted increasing attention due to their occurrence in regulatory regions of the genome and the subsequent implications for regulating critical stages of viral life cycles. In this context, the employment of specific G4 ligands has aided in comprehending the intricate G4-mediated regulatory mechanisms in the viral life cycle, showcasing the potential of targeting viral G4s as a novel antiviral strategy. This review offers a thorough update on the literature concerning G4s in viruses, including their identification and functional significance across most of the human-infecting viruses. Furthermore, it delves into potential therapeutic avenues targeting G4s, encompassing various G4-binding ligands, G4-interacting proteins, and oligonucleotide-based strategies. Finally, the article highlights both progress and challenges in the field, providing valuable insights into leveraging this unusual nucleic acid structure for therapeutic purposes.

Effect of montmorillonite K10 clay on RNA structure and function

One of the earliest living systems was likely based on RNA (“the RNA world”). Mineral surfaces have been postulated to be an important environment for the prebiotic chemistry of RNA. In addition to adsorbing RNA and thus potentially reducing the chance of parasitic takeover through limited diffusion, minerals have been shown to promote a range of processes related to the emergence of life, including RNA polymerization, peptide bond formation, and self-assembly of vesicles. In addition, self-cleaving ribozymes have been shown to retain activity when adsorbed to the clay mineral montmorillonite. However, simulation studies suggest that adsorption to minerals is likely to interfere with RNA folding and, thus, function. To further evaluate the plausibility of a mineral-adsorbed RNA world, here we studied the effect of the synthetic clay montmorillonite K10 on the malachite green RNA aptamer, including binding of the clay to malachite green and RNA, as well as on the formation of secondary structures in model RNA and DNA oligonucleotides. We evaluated the fluorescence of the aptamer complex, adsorption to the mineral, melting curves, Förster resonance energy transfer interactions, and 1H-NMR signals to study the folding and functionality of these nucleic acids. Our results indicate that while some base pairings are unperturbed, the overall folding and binding of the malachite green aptamer are substantially disrupted by montmorillonite. These findings suggest that minerals would constrain the structures, and possibly the functions, available to an adsorbed RNA world.

AliNA - a deep learning program for RNA secondary structure prediction.

Nowadays there are numerous discovered natural RNA variations participating in different cellular processes and artificial RNA, e.g., aptamers, riboswitches. One of the required tasks in the investigation of their functions and mechanism of influence on cells and interaction with targets is the prediction of RNA secondary structures. The classic thermodynamic-based prediction algorithms do not consider the specificity of biological folding and deep learning methods that were designed to resolve this issue suffer from homology-based methods problems. Herein, we present a method for RNA secondary structure prediction based on deep learning - AliNA (ALIgned Nucleic Acids). Our method successfully predicts secondary structures for non-homologous to train-data RNA families thanks to usage of the data augmentation techniques. Augmentation extends existing datasets with easily-accessible simulated data. The proposed method shows a high quality of prediction across different benchmarks including pseudoknots. The method is available on GitHub for free (https://github.com/Arty40m/AliNA).

Single-molecule quantitation of RNA-binding protein occupancy and stoichiometry defines a role for Yra1 (Aly/REF) in nuclear mRNP organization

RNA-binding proteins (RBPs) interact with mRNA to form supramolecular complexes called messenger ribonucleoprotein (mRNP) particles. These dynamic assemblies direct and regulate individual steps of gene expression; however, their composition and functional importance remain largely unknown. Here, we develop a total internal reflection fluorescence-based single-molecule imaging assay to investigate stoichiometry and co-occupancy of 15 RBPs within mRNPs from Saccharomyces cerevisiae. We show compositional heterogeneity of single mRNPs and plasticity across different growth conditions, with major co-occupants of mRNPs containing the nuclear cap-binding complex identified as Yra1 (1-10 copies), Nab2 (1-6 copies), and Npl3 (1-6 copies). Multicopy Yra1-bound mRNPs are specifically co-occupied by the THO complex and assembled on mRNAs biased by transcript length and RNA secondary structure. Yra1 depletion results in decreased compaction of nuclear mRNPs demonstrating a packaging function. Together, we provide a quantitative framework for gene- and condition-dependent RBP occupancy and stoichiometry in individual nuclear mRNPs.

Global genomic diversity and conservation of SARS-CoV-2 since the COVID-19 outbreak.

Global sequencing of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has enabled researchers to monitor the genetic information of new pathogenic variants. However, the integrity and accuracy of most genomic information cannot be guaranteed due to limited sample quality and sequencing technology. This restricts researchers from tracking the latest changes in the virus genome. In this study, we aimed to characterize the diversity and conservation balance of complete SARS-CoV-2 genomes worldwide since the outbreak. To achieve this goal, we collected 5,966,490 genome sequences from various parts of the world, excluding those with incomplete or unknown/degenerate sequences. Our methodology included comparing sequences using BLASTN and BLASTP, analyzing RNA secondary structure, and calculating entropy. Our findings provide insights into the characteristics of SARS-CoV-2 and potential mechanisms of pathogenesis. We found that uracil had the highest proportion of all bases among various coronaviruses and cytosine to uracil mutations had the highest proportion among all point mutations. The consistency in the front part (1-26,599 nt) was significantly higher than that in the back part (26,600-29,903 nt) of the genome. For most genes, the similar consistency characteristics were also observed in their protein families. IMPORTANCE Our results indicate that most severe acute respiratory syndrome coronavirus 2 genomes sampled from patients had a mutation rate ≤1.07 ‰ and genome-tail proteins (including S protein) were the main sources of genetic polymorphism. The analysis of the virus-host interaction network of genome-tail proteins showed that they shared some antiviral signaling pathways, especially the intracellular protein transport pathway.

The Boltzmann distributions of molecular structures predict likely changes through random mutations

New folded molecular structures can only evolve after arising through mutations. This aspect is modeled using genotype-phenotype maps, which connect sequence changes through mutations to changes in molecular structures. Previous work has shown that the likelihood of appearing through mutations can differ by orders of magnitude from structure to structure and that this can affect the outcomes of evolutionary processes. Thus, we focus on the phenotypic mutation probabilities φqp, i.e., the likelihood that a random mutation changes structure p into structure q. For both RNA secondary structures and the HP protein model, we show that a simple biophysical principle can explain and predict how this likelihood depends on the new structure q: φqp is high if sequences that fold into p as the minimum-free-energy structure are likely to have q as an alternative structure with high Boltzmann frequency. This generalizes the existing concept of plastogenetic congruence from individual sequences to the entire neutral spaces of structures. Our result helps us understand why some structural changes are more likely than others, may be useful for estimating these likelihoods via sampling and makes a connection to alternative structures with high Boltzmann frequency, which could be relevant in evolutionary processes.

DDX20: A Multifunctional Complex Protein.

DEAD-box decapping enzyme 20 (DDX20) is a putative RNA-decapping enzyme that can be identified by the conserved motif Asp-Glu-Ala-Asp (DEAD). Cellular processes involve numerous RNA secondary structure alterations, including translation initiation, nuclear and mitochondrial splicing, and assembly of ribosomes and spliceosomes. DDX20 reportedly plays an important role in cellular transcription and post-transcriptional modifications. On the one hand, DDX20 can interact with various transcription factors and repress the transcriptional process. On the other hand, DDX20 forms the survival motor neuron complex and participates in the assembly of snRNP, ultimately affecting the RNA splicing process. Finally, DDX20 can potentially rely on its RNA-unwinding enzyme function to participate in microRNA (miRNA) maturation and act as a component of the RNA-induced silencing complex. In addition, although DDX20 is not a key component in the innate immune system signaling pathway, it can affect the nuclear factor kappa B (NF-κB) and p53 signaling pathways. In particular, DDX20 plays different roles in tumorigenesis development through the NF-κB signaling pathway. This process is regulated by various factors such as miRNA. DDX20 can influence processes such as viral replication in cells by interacting with two proteins in Epstein-Barr virus and can regulate the replication process of several viruses through the innate immune system, indicating that DDX20 plays an important role in the innate immune system. Herein, we review the effects of DDX20 on the innate immune system and its role in transcriptional and post-transcriptional modification processes, based on which we provide an outlook on the future of DDX20 research in innate immunity and viral infections.

In vivo-like nearest neighbor parameters improve prediction of fractional RNA base-pairing in cells.

We conducted a thermodynamic analysis of RNA stability in Eco80 artificial cytoplasm, which mimics in vivo conditions, and compared it to transcriptome-wide probing of mRNA. Eco80 contains 80% of Escherichia coli metabolites, with biological concentrations of metal ions, including 2 mM free Mg2+ and 29 mM metabolite-chelated Mg2+. Fluorescence-detected binding isotherms (FDBI) were used to conduct a thermodynamic analysis of 24 RNA helices and found that these helices, which have an average stability of -12.3 kcal/mol, are less stable by ΔΔGo37 ∼1 kcal/mol. The FDBI data was used to determine a set of Watson-Crick free energy nearestneighbor parameters (NNPs), which revealed that Eco80 reduces the stability of three NNPs. This information was used to adjust the NN model using the RNAstructure package. The in vivo-like adjustments have minimal effects on the prediction of RNA secondary structures determined in vitro and in silico, but markedly improve prediction of fractional RNA base pairing in E. coli, as benchmarked with our in vivo DMS and EDC RNA chemical probing data. In summary, our thermodynamic and chemical probing analyses of RNA helices indicate that RNA secondary structures are less stable in cells than in artificially stable in vitro buffer conditions.