BEROLECMI: a novel prediction method to infer circRNA-miRNA interaction from the role definition of molecular attributes and biological networks.

Circular RNA (CircRNA)-microRNA (miRNA) interaction (CMI) is an important model for the regulation of biological processes by non-coding RNA (ncRNA), which provides a new perspective for the study of human complex diseases. However, the existing CMI prediction models mainly rely on the nearest neighbor structure in the biological network, ignoring the molecular network topology, so it is difficult to improve the prediction performance. In this paper, we proposed a new CMI prediction method, BEROLECMI, which uses molecular sequence attributes, molecular self-similarity, and biological network topology to define the specific role feature representation for molecules to infer the new CMI. BEROLECMI effectively makes up for the lack of network topology in the CMI prediction model and achieves the highest prediction performance in three commonly used data sets. In the case study, 14 of the 15 pairs of unknown CMIs were correctly predicted.

Tradeoffs in the design of RNA thermometers.

The synthesis of RNA thermometers is aimed at achieving temperature responses with desired thresholds and sensitivities. Although previous works have generated thermometers with a variety of thresholds and sensitivities as well as guidelines for design, possible constraints in the achievable thresholds and sensitivities remain unclear. We addressed this issue using a two-state model and its variants, as well as melt profiles generated from thermodynamic computations. In the two-state model, we found that the threshold was inversely proportional to the sensitivity, in the case of a fixed energy difference between the two states. Notably, this constraint could persist in variations of the two-state model with sequentially unfolding states and branched parallel pathways. Furthermore, the melt profiles generated from a library of thermometers exhibited a similar constraint. These results should inform the design of RNA thermometers as well as other responses that are mediated in a similar fashion.

Atomistic insights into the reentrant phase-transitions in polyuracil and polylysine mixtures.

The phase separation of protein and RNA mixtures underpins the assembly and regulation of numerous membraneless organelles in cells. The ubiquity of protein-RNA condensates in cellular regulatory processes is in part due to their sensitivity to RNA concentration, which affects their physical properties and stability. Recent experiments with poly-cationic peptide-RNA mixtures have revealed closed-loop phase diagrams featuring lower and upper critical solution temperatures. These diagrams indicate reentrant phase transitions shaped by biomolecular interactions and entropic forces such as solvent and ion reorganization. We employed atomistic simulations to study mixtures with various RNA-polylysine stoichiometries and temperatures to elucidate the microscopic driving forces behind reentrant phase transitions in protein-RNA mixtures. Our findings reveal an intricate interplay between hydration, ion condensation, and specific RNA-polylysine hydrogen bonding, resulting in distinct stoichiometry-dependent phase equilibria governing stabilities and structures of the condensate phase. Our simulations show that reentrant transitions are accompanied by desolvation around the phosphate groups of RNA, with increased contacts between phosphate and lysine side chains. In RNA-rich systems at lower temperatures, RNA molecules can form an extensive pi-stacking and hydrogen bond network, leading to percolation. In protein-rich systems, no such percolation-induced transitions are observed. Furthermore, we assessed the performance of three prominent water force fields-Optimal Point Charge (OPC), TIP4P-2005, and TIP4P-D-in capturing reentrant phase transitions. OPC provided a superior balance of interactions, enabling effective capture of reentrant transitions and accurate characterization of changes in solvent reorganization. This study offers atomistic insights into the nature of reentrant phase transitions using simple model peptide and nucleotide mixtures. We believe that our results are broadly applicable to larger classes of peptide-RNA mixtures exhibiting reentrant phase transitions.

Identifying small-molecules binding sites in RNA conformational ensembles with SHAMAN.

The rational targeting of RNA with small molecules is hampered by our still limited understanding of RNA structural and dynamic properties. Most in silico tools for binding site identification rely on static structures and therefore cannot face the challenges posed by the dynamic nature of RNA molecules. Here, we present SHAMAN, a computational technique to identify potential small-molecule binding sites in RNA structural ensembles. SHAMAN enables exploring the conformational landscape of RNA with atomistic molecular dynamics simulations and at the same time identifying RNA pockets in an efficient way with the aid of probes and enhanced-sampling techniques. In our benchmark composed of large, structured riboswitches as well as small, flexible viral RNAs, SHAMAN successfully identifies all the experimentally resolved pockets and ranks them among the most favorite probe hotspots. Overall, SHAMAN sets a solid foundation for future drug design efforts targeting RNA with small molecules, effectively addressing the long-standing challenges in the field.

Using Modified Synthetic Oligonucleotides to Assay Nucleic Acid-Metabolizing Enzymes.

The availability of a range of modified synthetic oligonucleotides from commercial vendors has allowed the development of sophisticated assays to characterize diverse properties of nucleic acid metabolizing enzymes that can be run in any standard molecular biology lab. The use of fluorescent labels has made these methods accessible to researchers with standard PAGE electrophoresis equipment and a fluorescent-enabled imager, without using radioactive materials or requiring a lab designed for the storage and preparation of radioactive materials, i.e., a Hot Lab. The optional addition of standard modifications such as phosphorylation can simplify assay setup, while the specific incorporation of modified nucleotides that mimic DNA damages or intermediates can be used to probe specific aspects of enzyme behavior. Here, the design and execution of assays to interrogate several aspects of DNA processing by enzymes using commercially available synthetic oligonucleotides are demonstrated. These include the ability of ligases to join or nucleases to degrade different DNA and RNA hybrid structures, differential cofactor usage by the DNA ligase, and evaluation of the DNA-binding capacity of enzymes. Factors to consider when designing synthetic nucleotide substrates are discussed, and a basic set of oligonucleotides that can be used for a range of nucleic acid ligase, polymerase, and nuclease enzyme assays are provided.

Advances in the investigation of N6-isopentenyl adenosine i6A RNA modification.

Prenylation (isopentenylation), a key post-transcriptional modification with a hydrophobic prenyl group onto the biomacromolecules such as RNA and proteins, influences their localization and function. Prenyltransferases mediate this process, while cytokinin oxidases degrade the prenylated adenosine in plants. This review summarizes current progress in detecting prenylation modifications in RNA across species and their effects on protein synthesis. Advanced methods have been developed to label and study these modifications in vitro and in vivo, despite challenges posed by the inert chemical properties of prenyl groups. Continued advancements in bioorthogonal chemistry promise new tools for understanding the precise biological functions of prenylated RNA modifications and other related proteins.

Information Gradient among Nucleotide Sequences of Essential RNAs from an Evolutionary Perspective.

We hypothesize that the first ancestral "protocell" molecular structures, i.e., the first RNAs and peptides that gradually transformed into real cells once the Earth had cooled sufficiently for organic molecules to appear, have left traces in the RNAs and the genes in present cells. We propose a circular RNA that could have been one of these ancestral structures whose vestigial pentameric subsequences would mark the evolution from this key moment when the protocells began to join with living organisms. In particular, we propose that, in present RNAs (ribosomal or messenger), which play an important role in the metabolism of current cells, we look for traces of the proposed primitive structure in the form of pentamers (or longer fragments) that belong to their nucleotide sequence. The result obtained can be summarized in the existence of a gradient of occurrence of such pentamers, with a high frequency for the most vital functions (protein synthesis, nucleic synthesis, cell respiration, etc.). This gradient is also visible between organisms, from the oldest (Archaea) to the most recent (Eukaryotes) in the evolution of species.

Phase Separation Modulates the Thermodynamics and Kinetics of RNA Hybridization.

Biomolecular condensates can influence cellular function in a number of ways, including by changing the structural dynamics and conformational equilibria of the molecules partitioned within them. Here we use methyl transverse relaxation optimized spectroscopy (methyl-TROSY) NMR in conjunction with 2'-O-methyl labeling of RNA to characterize the thermodynamics and kinetics of RNA-RNA base pairing in condensates formed by the C-terminal intrinsically disordered region of CAPRIN1, an RNA-binding protein involved in RNA transport, translation, and stability. CAPRIN1 condensates destabilize RNA-RNA base pairing, resulting from a ∼270-fold decrease and a concomitant ∼15-fold increase in the on- and off-rates for duplex formation, respectively. The ∼30-fold slower diffusion of RNA single strands within the condensed phase partially accounts for the reduced on-rate, but the further ∼9-fold reduction likely reflects shedding of CAPRIN1 chains that are interacting with the RNA prior to hybridization. Our study emphasizes the important role of protein solvation in modulating nucleic acid recognition processes inside condensates.

If the 5' cap fits (wear it) - Non-canonical RNA capping.

RNA capping is a prominent RNA modification that influences RNA stability, metabolism, and function. While it was long limited to the study of the most abundant eukaryotic canonical m7G cap, the field recently went through a large paradigm shift with the discovery of non-canonical RNA capping in bacteria and ultimately all domains of life. The repertoire of non-canonical caps has expanded to encompass metabolite caps, including NAD, FAD, CoA, UDP-Glucose, and ADP-ribose, alongside alarmone dinucleoside polyphosphate caps, and methylated phosphate cap-like structures. This review offers an introduction into the field, presenting a summary of the current knowledge about non-canonical RNA caps. We highlight the often still enigmatic biological roles of the caps together with their processing enzymes, focusing on the most recent discoveries. Furthermore, we present the methods used for the detection and analysis of these non-canonical RNA caps and thus provide an introduction into this dynamic new field.

Ratiometric sandwich-type assays for RNAs with a point mutation using benzo[a]pyrene-modified probes.

Benzo[a]pyrene-modified oligonucleotides were developed for the detection of RNAs with a point mutation. The probes produced two distinct fluorescence signals in response to single nucleotide differences in the RNA sequences, allowing for discrimination between the matched and single base mismatched RNA sequences in colorimetric and ratiometric manners.

RNA Adduction Resulting from the Metabolic Activation of Myristicin by P450 Enzymes and Sulfotransferases.

Myristicin (MYR) mainly occurs in nutmeg and belongs to alkoxy-substituted allylbenzenes, a class of potentially toxic natural chemicals. RNA interaction with MYR metabolites in vitro and in vivo has been investigated in order to gain a better understanding of MYR toxicities. We detected two guanosine adducts (GA1 and GA2), two adenosine adducts (AA1 and AA2), and two cytosine adducts (CA1 and CA2) by LC-MS/MS analysis of total RNA extracts from cultured primary mouse hepatocytes and liver tissues of mice after exposure to MYR. An order of nucleoside adductions was found to be GAs > AAs > CAs, and the result of density functional theory calculations was in agreement with that detected by the LC-MS/MS-based approach. In vitro and in vivo studies have shown that MYR was oxidized by cytochrome P450 enzymes to 1'-hydroxyl and 3'-hydroxyl metabolites, which were then sulfated by sulfotransferases (SULTs) to form sulfate esters. The resulting sulfates would react with the nucleosides by SN1 and/or SN2 reactions, resulting in RNA adduction. The modification may alter the biochemical properties of RNA and disrupt RNA functions, perhaps partially contributing to the toxicities of MYR.

Polyamines (PAs) but not small peptides with closely spaced positively charged groups interact with DNA and RNA, but they do not represent a relevant buffer system at physiological pH values.

Polyamines (PAs) including putrescine (PUT), spermidine (SPD) and spermine (SPM) are small, versatile molecules with two or more positively charged amino groups. Despite their importance for almost all forms of life, their specific roles in molecular and cellular biology remain partly unknown. The molecular structures of PAs suggest two presumable biological functions: (i) as potential buffer systems and (ii) as interactants with poly-negatively charged molecules like nucleic acids. The present report focuses on the question, whether the molecular structures of PAs are essential for such functions, or whether other simple molecules like small peptides with closely spaced positively charged side chains might be suitable as well. Consequently, we created titration curves for PUT, SPD, and SPM, as well as for oligolysines like tri-, tetra-, and penta-lysine. None of the molecules provided substantial buffering capacity at physiological intracellular pH values. Apparently, the most important mechanism for intracellular pH homeostasis in neurons is not a buffer system but is provided by the actions of the sodium-hydrogen and the bicarbonate-chloride antiporters. In a similar approach we investigated the interaction with DNA by following the extinction at 260 nm when titrating DNA with the above molecules. Again, PUT and tri-lysine were not able to interact with herring sperm DNA, while SPD and SPM were. Obviously, the presence of several positively charged groups on its own is not sufficient for the interaction with nucleic acids. Instead, the precise spacing of these groups is necessary for biological activity.

pyRBDome: a comprehensive computational platform for enhancing RNA-binding proteome data.

High-throughput proteomics approaches have revolutionised the identification of RNA-binding proteins (RBPome) and RNA-binding sequences (RBDome) across organisms. Yet, the extent of noise, including false positives, associated with these methodologies, is difficult to quantify as experimental approaches for validating the results are generally low throughput. To address this, we introduce pyRBDome, a pipeline for enhancing RNA-binding proteome data in silico. It aligns the experimental results with RNA-binding site (RBS) predictions from distinct machine-learning tools and integrates high-resolution structural data when available. Its statistical evaluation of RBDome data enables quick identification of likely genuine RNA-binders in experimental datasets. Furthermore, by leveraging the pyRBDome results, we have enhanced the sensitivity and specificity of RBS detection through training new ensemble machine-learning models. pyRBDome analysis of a human RBDome dataset, compared with known structural data, revealed that although UV-cross-linked amino acids were more likely to contain predicted RBSs, they infrequently bind RNA in high-resolution structures. This discrepancy underscores the limitations of structural data as benchmarks, positioning pyRBDome as a valuable alternative for increasing confidence in RBDome datasets.

Excited-Ground-State Transition of the RNA Strand Slippage Mechanism Captured by the Base-Specific Force Field.

Excited-ground-state transition and strand slippage of RNA play key roles in transcription and translation of central dogma. Due to limitation of current experimental techniques, the dynamic structure ensembles of RNA remain inadequately understood. Molecular dynamics simulations offer a promising complementary approach, whose accuracy depends on the force field. Here, we develop the new version of RNA base-specific force field (BSFF2) to address underestimation of base pairing stability and artificial backbone conformations. Extensive evaluations on typical RNA systems have comprehensively confirmed the accuracy of BSFF2. Furthermore, BSFF2 demonstrates exceptional efficiency in de novo folding of tetraloops and reproducing base pair reshuffling transition between RNA excited and ground states. Then, we explored the RNA strand slippage mechanism with BSFF2. We conducted a comprehensive three-dimensional structural investigation into the strand slippage of the most complex r(G4C2)9 repeat element and presented the molecular details in the dynamic transition along with the underlying mechanism. Our results of capturing the strand slippage, excited-ground transition, de novo folding, and simulations for various typical RNA motifs indicate that BSFF2 should be one of valuable tools for dynamic conformation research and structure prediction of RNA, and a future contribution to RNA-targeted drug design as well as RNA therapy development.

A Cyanine Dye for Highly Specific Recognition of Parallel G-Quadruplex Topology and Its Application in Clinical RNA Detection for Cancer Diagnosis.

G-quadruplex (G4), an unconventional nucleic acid structure, shows polymorphism in its topological morphology. The parallel G4 topology is the most prevalent form in organisms and plays a regulatory role in many biological processes. Designing fluorescent probes with high specificity for parallel G4s is important but challenging. Herein, a supramolecular assembly of the anionic cyanine dye SCY-5 is reported, which selectively identifies parallel G4 topology. SCY-5 can clearly distinguish parallel G4s from other G4s and non-G4s, even including hybrid-type G4s with parallel characteristics. The high specificity mechanism of SCY-5 involves a delicate balance between electrostatic repulsion and π-π interaction between SCY-5 and G4s. Using SCY-5, cellular RNA extracted from peripheral venous blood was quantitatively detected, and a remarkable increase in RNA G4 content in cancer patients compared to healthy volunteers was confirmed for the first time. This study provides new insights for designing specific probes for parallel G4 topology and opens a new path for clinical cancer diagnosis using RNA G4 as a biomarker.

Template-based copying in chemically fuelled dynamic combinatorial libraries.

One of science's greatest challenges is determining how life can spontaneously emerge from a mixture of molecules. A complicating factor is that life and its molecules are inherently unstable-RNA and proteins are prone to hydrolysis and denaturation. For the de novo synthesis of life or to better understand its emergence at its origin, selection mechanisms are needed for unstable molecules. Here we present a chemically fuelled dynamic combinatorial library to model RNA oligomerization and deoligomerization and shine new light on selection and purification mechanisms under kinetic control. In the experiments, oligomers can only be sustained by continuous production. Hybridization is a powerful tool for selecting unstable molecules, offering feedback on oligomerization and deoligomerization rates. Moreover, we find that templation can be used to purify libraries of oligomers. In addition, template-assisted formation of oligomers within coacervate-based protocells changes its compartment's physical properties, such as their ability to fuse. Such reciprocal coupling between oligomer production and physical properties is a key step towards synthetic life.

Protocell Effects on RNA Folding, Function, and Evolution.

ConspectusCreating a living system from nonliving matter is a great challenge in chemistry and biophysics. The early history of life can provide inspiration from the idea of the prebiotic "RNA World" established by ribozymes, in which all genetic and catalytic activities were executed by RNA. Such a system could be much simpler than the interdependent central dogma characterizing life today. At the same time, cooperative systems require a mechanism such as cellular compartmentalization in order to survive and evolve. Minimal cells might therefore consist of simple vesicles enclosing a prebiotic RNA metabolism.The internal volume of a vesicle is a distinctive environment due to its closed boundary, which alters diffusion and available volume for macromolecules and changes effective molecular concentrations, among other considerations. These physical effects are mechanistically distinct from chemical interactions, such as electrostatic repulsion, that might also occur between the membrane boundary and encapsulated contents. Both indirect and direct interactions between the membrane and RNA can give rise to nonintuitive, "emergent" behaviors in the model protocell system. We have been examining how encapsulation inside membrane vesicles would affect the folding and activity of entrapped RNA.Using biophysical techniques such as FRET, we characterized ribozyme folding and activity inside vesicles. Encapsulation inside model protocells generally promoted RNA folding, consistent with an excluded volume effect, independently of chemical interactions. This energetic stabilization translated into increased ribozyme activity in two different systems that were studied (hairpin ribozyme and self-aminoacylating RNAs). A particularly intriguing finding was that encapsulation could rescue the activity of mutant ribozymes, suggesting that encapsulation could affect not only folding and activity but also evolution. To study this further, we developed a high-throughput sequencing assay to measure the aminoacylation kinetics of many thousands of ribozyme variants in parallel. The results revealed an unexpected tendency for encapsulation to improve the better ribozyme variants more than worse variants. During evolution, this effect would create a tilted playing field, so to speak, that would give additional fitness gains to already-high-activity variants. According to Fisher's Fundamental Theorem of Natural Selection, the increased variance in fitness should manifest as faster evolutionary adaptation. This prediction was borne out experimentally during in vitro evolution, where we observed that the initially diverse ribozyme population converged more quickly to the most active sequences when they were encapsulated inside vesicles.The studies in this Account have expanded our understanding of emergent protocell behavior, by showing how simply entrapping an RNA inside a vesicle, which could occur spontaneously during vesicle formation, might profoundly affect the evolutionary landscape of the RNA. Because of the exponential dynamics of replication and selection, even small changes to activity and function could lead to major evolutionary consequences. By closely studying the details of minimal yet surprisingly complex protocells, we might one day trace a pathway from encapsulated RNA to a living system.

Biophysical modeling with variational autoencoders for bimodal, single-cell RNA sequencing data.

Here we present biVI, which combines the variational autoencoder framework of scVI with biophysical models describing the transcription and splicing kinetics of RNA molecules. We demonstrate on simulated and experimental single-cell RNA sequencing data that biVI retains the variational autoencoder's ability to capture cell type structure in a low-dimensional space while further enabling genome-wide exploration of the biophysical mechanisms, such as system burst sizes and degradation rates, that underlie observations.

FURNA: A database for functional annotations of RNA structures.

Despite the increasing number of 3D RNA structures in the Protein Data Bank, the majority of experimental RNA structures lack thorough functional annotations. As the significance of the functional roles played by noncoding RNAs becomes increasingly apparent, comprehensive annotation of RNA function is becoming a pressing concern. In response to this need, we have developed FURNA (Functions of RNAs), the first database for experimental RNA structures that aims to provide a comprehensive repository of high-quality functional annotations. These include Gene Ontology terms, Enzyme Commission numbers, ligand-binding sites, RNA families, protein-binding motifs, and cross-references to related databases. FURNA is available at https://seq2fun.dcmb.med.umich.edu/furna/ to enable quick discovery of RNA functions from their structures and sequences.