GARN: Sampling RNA 3D Structure Space with Game Theory and Knowledge-Based Scoring Strategies

Significant improvements have been made in the efficiency and accuracy of RNA 3D structure prediction methods during the succeeding challenges of RNA-Puzzles, a community-wide effort on the assessment of blind prediction of RNA tertiary structures. The RNA-Puzzles contest has shown, among others, that the development and validation of computational methods for RNA fold prediction strongly depend on the benchmark datasets and the structure comparison algorithms. Yet, there has been no systematic benchmark set or decoy structures available for the 3D structure prediction of RNA, hindering the standardization of comparative tests in the modeling of RNA structure. Furthermore, there has not been a unified set of tools that allows deep and complete RNA structure analysis, and at the same time, that is easy to use. Here, we present RNA-Puzzles toolkit, a computational resource including (i) decoy sets generated by different RNA 3D structure prediction methods (raw, for-evaluation and standardized datasets), (ii) 3D structure normalization, analysis, manipulation, visualization tools (RNA_format, RNA_normalizer, rna-tools) and (iii) 3D structure comparison metric tools (RNAQUA, MCQ4Structures). This resource provides a full list of computational tools as well as a standard RNA 3D structure prediction assessment protocol for the community.

RNA molecules with novel functions have revived interest in the accurate prediction of RNA three-dimensional (3D) structure and folding dynamics. However, existing methods are inefficient in automated 3D structure prediction. Here, we report a robust computational approach for rapid folding of RNA molecules. We develop a simplified RNA model for discrete molecular dynamics (DMD) simulations, incorporating base-pairing and base-stacking interactions. We demonstrate correct folding of 150 structurally diverse RNA sequences. The majority of DMD-predicted 3D structures have <4 A deviations from experimental structures. The secondary structures corresponding to the predicted 3D structures consist of 94% native base-pair interactions. Folding thermodynamics and kinetics of tRNA(Phe), pseudoknots, and mRNA fragments in DMD simulations are in agreement with previous experimental findings. Folding of RNA molecules features transient, non-native conformations, suggesting non-hierarchical RNA folding. Our method allows rapid conformational sampling of RNA folding, with computational time increasing linearly with RNA length. We envision this approach as a promising tool for RNA structural and functional analyses.

In order to carry out biological functions, RNA molecules must fold into specific three-dimensional (3D) structures. Current experimental methods to determine RNA 3D structures are expensive and time consuming. With the recent advances in computational biology, RNA structure prediction is becoming increasingly reliable. This chapter describes a recently developed RNA structure prediction software, Vfold, a virtual bond-based RNA folding model. The main features of Vfold are the physics-based loop free energy calculations for various RNA structure motifs and a template-based assembly method for RNA 3D structure prediction. For illustration, we use the yybP-ykoY Orphan riboswitch as an example to show the implementation of the Vfold model in RNA structure prediction from the sequence.

Many recent exciting discoveries have revealed the versatility of RNAs and their importance in a variety of cellular functions which are strongly coupled to RNA structures. To understand the functions of RNAs, some structure prediction models have been developed in recent years. In this review, the progress in computational models for RNA structure prediction is introduced and the distinguishing features of many outstanding algorithms are discussed, emphasizing three-dimensional (3D) structure prediction. A promising coarse-grained model for predicting RNA 3D structure, stability and salt effect is also introduced briefly. Finally, we discuss the major challenges in the RNA 3D structure modeling.

The biosynthesis of several protein cofactors is subject to feedback regulation by riboswitches. Flavin mononucleotide (FMN)-specific riboswitches, also known as RFN elements, direct expression of bacterial genes involved in the biosynthesis and transport of riboflavin (vitamin B(2)) and related compounds. Here we present the crystal structures of the Fusobacterium nucleatum riboswitch bound to FMN, riboflavin and antibiotic roseoflavin. The FMN riboswitch structure, centred on an FMN-bound six-stem junction, does not fold by collinear stacking of adjacent helices, typical for folding of large RNAs. Rather, it adopts a butterfly-like scaffold, stapled together by opposingly directed but nearly identically folded peripheral domains. FMN is positioned asymmetrically within the junctional site and is specifically bound to RNA through interactions with the isoalloxazine ring chromophore and direct and Mg(2+)-mediated contacts with the phosphate moiety. Our structural data, complemented by binding and footprinting experiments, imply a largely pre-folded tertiary RNA architecture and FMN recognition mediated by conformational transitions within the junctional binding pocket. The inherent plasticity of the FMN-binding pocket and the availability of large openings make the riboswitch an attractive target for structure-based design of FMN-like antimicrobial compounds. Our studies also explain the effects of spontaneous and antibiotic-induced deregulatory mutations and provided molecular insights into FMN-based control of gene expression in normal and riboflavin-overproducing bacterial strains.

The close relationship between RNA structure and function underlines the significance of accurately predicting RNA structures from sequence information. Structural topologies such as pseudoknots are of particular interest due to their ubiquity and direct involvement in RNA function, but identifying pseudoknots is a computationally challenging problem and existing heuristic approaches usually perform poorly for RNA sequences of even a few hundred bases. We survey the performance of pseudoknot prediction methods on a data set of full-length RNA sequences representing varied sequence lengths, and biological RNA classes such as RNase P RNA, Group I Intron, tmRNA and tRNA. Pseudoknot prediction methods are compared with minimum free energy and suboptimal secondary structure prediction methods in terms of correct base-pairs, stems and pseudoknots and we find that the ensemble of suboptimal structure predictions succeeds in identifying correct structural elements in RNA that are usually missed in MFE and pseudoknot predictions. We propose a strategy to identify a comprehensive set of non-redundant stems in the suboptimal structure space of a RNA molecule by applying heuristics that reduce the structural redundancy of the predicted suboptimal structures by merging slightly varying stems that are predicted to form in local sequence regions. This reduced-redundancy set of structural elements consistently outperforms more specialized approaches.in data sets. Thus, the suboptimal folding space can be used to represent the structural diversity of an RNA molecule more comprehensively than optimal structure prediction approaches alone.

The prediction of RNA structure with pseudoknots is a nondeterministic polynomial-time hard (NP-hard) problem; according to minimum free energy models and computational methods, we investigate the RNA-pseudoknotted structure. Our paper presents an efficient algorithm for predicting RNA structure with pseudoknots, and the algorithm takes O([Formula: see text]) time and O([Formula: see text]) space, the experimental tests in Rfam10.1 and PseudoBase indicate that the algorithm is more effective and precise. The predicting accuracy, the time complexity and space complexity outperform existing algorithms, such as Maximum Weight Matching (MWM) algorithm, PKNOTS algorithm and Inner Limiting Layer (ILM) algorithm, and the algorithm can predict arbitrary pseudoknots. And there exists a [Formula: see text] ([Formula: see text]) polynomial time approximation scheme in searching maximum number of stackings, and we give the proof of the approximation scheme in RNA-pseudoknotted structure. We have improved several types of pseudoknots considered in RNA folding structure, and analyze their possible transitions between types of pseudoknots.

RNA molecules play essential roles in fundamental biological mechanisms, including the regulation of gene expression and the catalysis of enzymatic reactions. Deciphering RNA structure is critical to understanding its diverse functions, defining its roles within cellular contexts, and enabling the design of RNA-targeted therapeutic strategies. This chapter offers a comprehensive overview of computational techniques developed for predicting RNA secondary and tertiary structures. We begin by introducing the core principles that govern RNA folding, such as base pairing rules, common secondary structure motifs, and tertiary interactions. The discussion then progresses to computational frameworks and comparative sequence analysis methods used to infer RNA secondary structures. The chapter also explores various approaches for modeling RNA three-dimensional architecture, including molecular modeling, coarse-grained simulations, and integrative methodologies that may incorporate experimental data and artificial intelligence-based predictions. Throughout the chapter, we evaluate the advantages, limitations, and practical applications of these tools, with an emphasis on how they contribute to elucidating RNA structure–function relationships.

Bulges in nucleic acids

Highlights of the IMA workshop RNA in Biology, Bioengineering, and Biotechnology are summarized, including recent developments in RNA secondary structure prediction and RNA design, innovative mathematical constructs for RNA structure, bioinformatics advances in RNA structure analysis and prediction, and experimental progress in RNA folding and imaging. In small things considered is the world revealed. Robert Lee Hotz, Feb. 2008.

TriRNASP: A knowledge-based potential with three-body effects for accurate RNA structure evaluation.

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty-two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and x-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as noncanonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

Single-particle cryo-electron microscopy (cryo-EM) can reveal the structures of large and often dynamic molecules, but smaller biomolecules (≤50 kDa) remain challenging targets due to their intrinsic low signal to noise ratio. Methods to help resolve small proteins have been applied but development of similar approaches to aid in structural determination of small, structured RNA elements have lagged. Here, we present a scaffold-based approach that we used to recover maps of sub-25 kDa RNA domains to 4.5-5.0 Å. While lacking the detail of true high-resolution maps, these maps are suitable for model building and preliminary structure determination. We demonstrate this method helped faithfully recover the structure of several RNA elements of known structure, and that it promises to be generalized to other RNAs without disturbing their native fold. This approach may streamline the sample preparation process and reduce the optimization required for data collection. This first-generation scaffold approach provides a robust system to aid in RNA structure determination by cryo-EM and lays the groundwork for further scaffold optimization to achieve higher resolution.

There are some NP-hard problems in the prediction of RNA structures. Prediction of RNA folding structure in RNA nucleotide sequence remains an unsolved challenge. We investigate the computing algorithm in RNA folding structural prediction based on extended structure and basin hopping graph, it is a computing mode of basin hopping graph in RNA folding structural prediction including pseudoknots. This study presents the predicting algorithm based on extended structure, it also proposes an improved computing algorithm based on barrier tree and basin hopping graph, which are the attractive approaches in RNA folding structural prediction. Many experiments have been implemented in Rfam14.1 database and PseudoBase database, the experimental results show that our two algorithms are efficient and accurate than the other existing algorithms.

Untersuchungen zum Mechanismus der katalytischen Aktivierung von Spleißosomen aus Saccharomyces Cerevisiae