ModeRNA: a tool for comparative modeling of RNA 3D structure

State-of-the-RNArt: benchmarking current methods for RNA 3D structure prediction

RNAs are essential molecules involved in numerous biological functions. Understanding RNA functions requires the knowledge of their 3D structures. Computational methods have been developed for over two decades to predict the 3D conformations from RNA sequences. These computational methods have been widely used and are usually categorised as either ab initio or template-based. The performances remain to be improved. Recently, the rise of deep learning has changed the sight of novel approaches. Deep learning methods are promising, but their adaptation to RNA 3D structure prediction remains difficult. In this paper, we give a brief review of the ab initio, template-based and novel deep learning approaches. We highlight the different available tools and provide a benchmark on nine methods using the RNA-Puzzles dataset. We provide an online dashboard that shows the predictions made by benchmarked methods, freely available on the EvryRNA platform: https://evryrna.ibisc.univ-evry.fr/evryrna/state_of_the_rnart/.

Prediction of 3D RNA Structures from Sequence Using Energy Landscapes of RNA Dimers: Application to RNA Tetraloops

Access to the three-dimensional structure of RNA enables an ability to gain a more profound understanding of its biological mechanisms, as well as the ability to design RNA-targeting drugs, which can take advantage of the unique chemical environment imposed by a folded RNA structure. Due to the dynamic and structurally complex properties of RNA, both experimental and traditional computational methods have difficulty in determining RNA's 3D structure. Herein, we introduce TAPERSS (Theoretical Analyses, Prediction, and Evaluation of RNA Structures from Sequence), a physics-based fragment assembly method for predicting 3D RNA structures from sequence. Using a fragment library created using discrete path sampling calculations of RNA dinucleoside monophosphates, TAPERSS can sample the physics-based energy landscapes of any RNA sequence with relatively low computational complexity. We have benchmarked TAPERSS on 21 RNA tetraloops, using a combinatorial algorithm as a proof-of-concept. We show that TAPERSS was successfully able to predict the apo-state structures of all 21 RNA hairpins, with 16 of those structures also having low predicted energies as well. We demonstrate that TAPERSS performs most accurately on GNRA-like tetraloops with mostly stacked loop-nucleotides, while having limited success with more dynamic UNCG and CUYG tetraloops, most likely due to the influence of the RNA force field used to create the fragment library. Moreover, we show that TAPERSS can successfully predict the majority of the experimental non-apo states, highlighting its potential in anticipating biologically significant yet unobserved states. This holds great promise for future applications in drug design and related studies. With discussed improvements and implementation of more efficient sampling algorithms, we believe TAPERSS may serve as a useful tool for a physics-based conformational sampling of large RNA structures.

Dynamics of Protein-RNA Interfaces Using All-Atom Molecular Dynamics Simulations

Facing the current challenges posed by human health diseases requires the understanding of cell machinery at a molecular level. The interplay between proteins and RNA is key for any physiological phenomenon, as well protein-RNA interactions. To understand these interactions, many experimental techniques have been developed, spanning a very wide range of spatial and temporal resolutions. In particular, the knowledge of tridimensional structures of protein-RNA complexes provides structural, mechanical, and dynamical pieces of information essential to understand their functions. To get insights into the dynamics of protein-RNA complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on nine different protein-RNA complexes with different functions and interface size by taking into account the bound and unbound forms. First, we characterized structural changes upon binding and, for the RNA part, the change in the puckering. Second, we extensively analyzed the interfaces, their dynamics and structural properties, and the structural waters involved in the binding, as well as the contacts mediated by them. Based on our analysis, the interfaces rearranged during the simulation time showing alternative and stable residue-residue contacts with respect to the experimental structure.

Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation

Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.

RNA threading with secondary structure and sequence profile

MOTIVATION

RNA threading aims to identify remote homologies for template-based modeling of RNA 3D structure. Existing RNA alignment methods primarily rely on secondary structure alignment. They are often time- and memory-consuming, limiting large-scale applications. In addition, the accuracy is far from satisfactory.

RESULTS

Using RNA secondary structure and sequence profile, we developed a novel RNA threading algorithm, named RNAthreader. To enhance the alignment process and minimize memory usage, a novel approach has been introduced to simplify RNA secondary structures into compact diagrams. RNAthreader employs a two-step methodology. Initially, integer programming and dynamic programming are combined to create an initial alignment for the simplified diagram. Subsequently, the final alignment is obtained using dynamic programming, taking into account the initial alignment derived from the previous step. The benchmark test on 80 RNAs illustrates that RNAthreader generates more accurate alignments than other methods, especially for RNAs with pseudoknots. Another benchmark, involving 30 RNAs from the RNA-Puzzles experiments, exhibits that the models constructed using RNAthreader templates have a lower average RMSD than those created by alternative methods. Remarkably, RNAthreader takes less than two hours to complete alignments with ∼5000 RNAs, which is 3-40 times faster than other methods. These compelling results suggest that RNAthreader is a promising algorithm for RNA template detection.

AVAILABILITY AND IMPLEMENTATION

https://yanglab.qd.sdu.edu.cn/RNAthreader.

Inconsistencies in the Classification of the Family Cydnidae (Hemiptera: Heteroptera: Pentatomoidea) Revealed by Molecular Apomorphies in the Secondary and Tertiary Structures of 18S rRNA Length-Variable Region L (LVR L)

The SSU nuclear rDNA (encoding 18S ribosomal RNA) is one of the most frequently sequenced genes in the molecular analysis of insects. Molecular apomorphies in the secondary and tertiary structures of several 18S rRNA length-variable regions (LVRs) located within the V2, V4, and V7 hypervariable regions can be good indicators for recovering monophyletic groups within some heteropteran families. Among the LVRs that have been analysed, the LVR L in the V4 hypervariable region is the longest and most crucial for such assessments. We analysed the 18S rRNA V4 hypervariable region sequences of 45 species from the family Cydnidae, including all 6 subfamilies (Amaurocorinae, Amnestinae, Cephalocteinae, Cydninae, Garsauriinae, and Sehirinae) and three pentatomoid families (Parastrachiidae, Thaumastellidae, and Thyreocoridae), which have often been included in the broadly defined Cydnidae family. This is the first time that representatives of all Cydnidae subfamilies have been included in a molecular analysis. Only taxa from two subfamilies, Sehirinae and Cydninae, have been used in previous molecular studies. The secondary and tertiary structures of the LVR L were predicted for each species using the two-step procedure already accepted for such analyses to recover any molecular apomorphy essential for determining monophyly. The results of our comparative studies contradict the current understanding of the relationships among burrowing bugs and the current family classification.

RNA G-quadruplex folding is a multi-pathway process driven by conformational entropy

The kinetics of folding is crucial for the function of many regulatory RNAs including RNA G-quadruplexes (rG4s). Here, we characterize the folding pathways of a G-quadruplex from the telomeric repeat-containing RNA by combining all-atom molecular dynamics and coarse-grained simulations with circular dichroism experiments. The quadruplex fold is stabilized by cations and thus, the ion atmosphere forming a double layer surrounding the highly charged quadruplex guides the folding process. To capture the ionic double layer in implicit solvent coarse-grained simulations correctly, we develop a matching procedure based on all-atom simulations in explicit water. The procedure yields quantitative agreement between simulations and experiments as judged by the populations of folded and unfolded states at different salt concentrations and temperatures. Subsequently, we show that coarse-grained simulations with a resolution of three interaction sites per nucleotide are well suited to resolve the folding pathways and their intermediate states. The results reveal that the folding progresses from unpaired chain via hairpin, triplex and double-hairpin constellations to the final folded structure. The two- and three-strand intermediates are stabilized by transient Hoogsteen interactions. Each pathway passes through two on-pathway intermediates. We hypothesize that conformational entropy is a hallmark of rG4 folding. Conformational entropy leads to the observed branched multi-pathway folding process for TERRA25. We corroborate this hypothesis by presenting the free energy landscapes and folding pathways of four rG4 systems with varying loop length.

RNA tertiary structure modeling with BRiQ potential in CASP15

We describe the modeling method for RNA tertiary structures employed by team AIchemy_RNA2 in the 15th Critical Assessment of Structure Prediction (CASP15). The method consists of the following steps. Firstly, secondary structure information was derived from various manually-verified sources. With this information, the full length RNA was fragmented into structural modules. The structures of each module were predicted and then assembled into the full structure. To reduce the searching conformational space, an RNA structure was organized into an optimal base folding tree. And to further improve the sampling efficiency, the energy surface was smoothed at high temperatures during the Monte Carlo sampling to make it easier to move across the energy barrier. The statistical potential energy function BRiQ was employed during Monte Carlo energy optimization.

RNA tertiary structure prediction in CASP15 by the GeneSilico group: Folding simulations based on statistical potentials and spatial restraints

Ribonucleic acid (RNA) molecules serve as master regulators of cells by encoding their biological function in the ribonucleotide sequence, particularly their ability to interact with other molecules. To understand how RNA molecules perform their biological tasks and to design new sequences with specific functions, it is of great benefit to be able to computationally predict how RNA folds and interacts in the cellular environment. Our workflow for computational modeling of the 3D structures of RNA and its interactions with other molecules uses a set of methods developed in our laboratory, including MeSSPredRNA for predicting canonical and non-canonical base pairs, PARNASSUS for detecting remote homology based on comparisons of sequences and secondary structures, ModeRNA for comparative modeling, the SimRNA family of programs for modeling RNA 3D structure and its complexes with other molecules, and QRNAS for model refinement. In this study, we present the results of testing this workflow in predicting RNA 3D structures in the CASP15 experiment. The overall high score of the computational models predicted by our group demonstrates the robustness of our workflow and its individual components in terms of predicting RNA 3D structures of acceptable quality that are close to the target structures. However, the variance in prediction quality is still quite high, and the results are still too far from the level of protein 3D structure predictions. This exercise led us to consider several improvements, especially to better predict and enforce stacking interactions and non-canonical base pairs.

Structural and dynamic mechanisms for coupled folding and tRNA recognition of a translational T-box riboswitch

T-box riboswitches are unique riboregulators where gene regulation is mediated through interactions between two highly structured RNAs. Despite extensive structural insights, how RNA-RNA interactions drive the folding and structural transitions of T-box to achieve functional conformations remains unclear. Here, by combining SAXS, single-molecule FRET and computational modeling, we elaborate the folding energy landscape of a translational T-box aptamer consisting of stems I, II and IIA/B, which Mg2+-induced global folding and tRNA binding are cooperatively coupled. smFRET measurements reveal that high Mg2+ stabilizes IIA/B and its stacking on II, which drives the pre-docking of I and II into a competent conformation, subsequent tRNA binding promotes docking of I and II to form a high-affinity tRNA binding groove, of which the essentiality of IIA/B and S-turn in II is substantiated with mutational analysis. We highlight a delicate balance among Mg2+, the intra- and intermolecular RNA-RNA interactions in modulating RNA folding and function.

trRosettaRNA: automated prediction of RNA 3D structure with transformer network

RNA 3D structure prediction is a long-standing challenge. Inspired by the recent breakthrough in protein structure prediction, we developed trRosettaRNA, an automated deep learning-based approach to RNA 3D structure prediction. The trRosettaRNA pipeline comprises two major steps: 1D and 2D geometries prediction by a transformer network; and 3D structure folding by energy minimization. Benchmark tests suggest that trRosettaRNA outperforms traditional automated methods. In the blind tests of the 15th Critical Assessment of Structure Prediction (CASP15) and the RNA-Puzzles experiments, the automated trRosettaRNA predictions for the natural RNAs are competitive with the top human predictions. trRosettaRNA also outperforms other deep learning-based methods in CASP15 when measured by the Z-score of the Root-Mean-Square Deviation. Nevertheless, it remains challenging to predict accurate structures for synthetic RNAs with an automated approach. We hope this work could be a good start toward solving the hard problem of RNA structure prediction with deep learning.

The structure of a hibernating ribosome in a Lyme disease pathogen

The spirochete bacterial pathogen Borrelia (Borreliella) burgdorferi (Bbu) affects more than 10% of the world population and causes Lyme disease in about half a million people in the US annually. Therapy for Lyme disease includes antibiotics that target the Bbu ribosome. Here we present the structure of the Bbu 70S ribosome obtained by single particle cryo-electron microscopy at 2.9 Å resolution, revealing a bound hibernation promotion factor protein and two genetically non-annotated ribosomal proteins bS22 and bL38. The ribosomal protein uL30 in Bbu has an N-terminal α-helical extension, partly resembling the mycobacterial bL37 protein, suggesting evolution of bL37 and a shorter uL30 from a longer uL30 protein. Its analogy to proteins uL30m and mL63 in mammalian mitochondrial ribosomes also suggests a plausible evolutionary pathway for expansion of protein content in mammalian mitochondrial ribosomes. Computational binding free energy predictions for antibiotics reflect subtle distinctions in antibiotic-binding sites in the Bbu ribosome. Discovery of these features in the Bbu ribosome may enable better ribosome-targeted antibiotic design for Lyme disease treatment.

When will RNA get its AlphaFold moment?

The protein structure prediction problem has been solved for many types of proteins by AlphaFold. Recently, there has been considerable excitement to build off the success of AlphaFold and predict the 3D structures of RNAs. RNA prediction methods use a variety of techniques, from physics-based to machine learning approaches. We believe that there are challenges preventing the successful development of deep learning-based methods like AlphaFold for RNA in the short term. Broadly speaking, the challenges are the limited number of structures and alignments making data-hungry deep learning methods unlikely to succeed. Additionally, there are several issues with the existing structure and sequence data, as they are often of insufficient quality, highly biased and missing key information. Here, we discuss these challenges in detail and suggest some steps to remedy the situation. We believe that it is possible to create an accurate RNA structure prediction method, but it will require solving several data quality and volume issues, usage of data beyond simple sequence alignments, or the development of new less data-hungry machine learning methods.

NuFold: A Novel Tertiary RNA Structure Prediction Method Using Deep Learning with Flexible Nucleobase Center Representation

RNA is not only playing a core role in the central dogma as mRNA between DNA and protein, but also many non-coding RNAs have been discovered to have unique and diverse biological functions. As genome sequences become increasingly available and our knowledge of RNA sequences grows, the study of RNA's structure and function has become more demanding. However, experimental determination of three-dimensional RNA structures is both costly and time-consuming, resulting in a substantial disparity between RNA sequence data and structural insights. In response to this challenge, we propose a novel computational approach that harnesses state-of-the-art deep learning architecture NuFold to accurately predict RNA tertiary structures. This approach aims to offer a cost-effective and efficient means of bridging the gap between RNA sequence information and structural comprehension. NuFold implements a nucleobase center representation, which allows it to reproduce all possible nucleotide conformations accurately.

Integrating end-to-end learning with deep geometrical potentials for ab initio RNA structure prediction

RNAs are fundamental in living cells and perform critical functions determined by their tertiary architectures. However, accurate modeling of 3D RNA structure remains a challenging problem. We present a novel method, DRfold, to predict RNA tertiary structures by simultaneous learning of local frame rotations and geometric restraints from experimentally solved RNA structures, where the learned knowledge is converted into a hybrid energy potential to guide RNA structure assembly. The method significantly outperforms previous approaches by >73.3% in TM-score on a sequence-nonredundant dataset containing recently released structures. Detailed analyses showed that the major contribution to the improvements arise from the deep end-to-end learning supervised with the atom coordinates and the composite energy function integrating complementary information from geometry restraints and end-to-end learning models. The open-source DRfold program with fast training protocol allows large-scale application of high-resolution RNA structure modeling and can be further improved with future expansion of RNA structure databases.

RNA 3D structure modeling by fragment assembly with small-angle X-ray scattering restraints

SUMMARY

Structure determination is a key step in the functional characterization of many non-coding RNA molecules. High-resolution RNA 3D structure determination efforts, however, are not keeping up with the pace of discovery of new non-coding RNA sequences. This increases the importance of computational approaches and low-resolution experimental data, such as from the small-angle X-ray scattering experiments. We present RNA Masonry, a computer program and a web service for a fully automated modeling of RNA 3D structures. It assemblies RNA fragments into geometrically plausible models that meet user-provided secondary structure constraints, restraints on tertiary contacts, and small-angle X-ray scattering data. We illustrate the method description with detailed benchmarks and its application to structural studies of viral RNAs with SAXS restraints.

AVAILABILITY AND IMPLEMENTATION

The program web server is available at http://iimcb.genesilico.pl/rnamasonry. The source code is available at https://gitlab.com/gchojnowski/rnamasonry.

Modeling the 3D structure and conformational dynamics of very large RNAs using coarse-grained molecular simulations

We describe a computational approach to building and simulating realistic 3D models of very large RNA molecules (>1000 nucleotides) at a resolution of one "bead" per nucleotide. The method starts with a predicted secondary structure and uses several stages of energy minimization and Brownian dynamics (BD) simulation to build 3D models. A key step in the protocol is the temporary addition of a 4 ^th spatial dimension that allows all predicted helical elements to become disentangled from each other in an effectively automated way. We then use the resulting 3D models as input to Brownian dynamics simulations that include hydrodynamic interactions (HIs) that allow the diffusive properties of the RNA to be modelled as well as enabling its conformational dynamics to be simulated. To validate the dynamics part of the method, we first show that when applied to small RNAs with known 3D structures the BD-HI simulation models accurately reproduce their experimental hydrodynamic radii (Rh). We then apply the modelling and simulation protocol to a variety of RNAs for which experimental Rh values have been reported ranging in size from 85 to 3569 nucleotides. We show that the 3D models, when used in BD-HI simulations, produce hydrodynamic radii that are usually in good agreement with experimental estimates for RNAs that do not contain tertiary contacts that persist even under very low salt conditions. Finally, we show that sampling of the conformational dynamics of large RNAs on timescales of 100 µs is computationally feasible with BD-HI simulations.

In silico methods for predicting functional synonymous variants

Single nucleotide variants (SNVs) contribute to human genomic diversity. Synonymous SNVs are previously considered to be "silent," but mounting evidence has revealed that these variants can cause RNA and protein changes and are implicated in over 85 human diseases and cancers. Recent improvements in computational platforms have led to the development of numerous machine-learning tools, which can be used to advance synonymous SNV research. In this review, we discuss tools that should be used to investigate synonymous variants. We provide supportive examples from seminal studies that demonstrate how these tools have driven new discoveries of functional synonymous SNVs.

P-FARFAR2: A multithreaded greedy approach to sampling low-energy RNA structures in Rosetta FARFAR2

RNA (ribonucleic acid) structure prediction finds many applications in health science and drug discovery due to its importance in several life regulatory processes. But despite significant advances in the close field of protein prediction, RNA 3D structure still poses a tremendous challenge to predict, especially for large sequences. In this regard, the approach unfolded by Rosetta FARFAR2 (Fragment Assembly of RNA with Full-Atom Refinement, version 2) has shown promising results, but the algorithm is non-deterministic by nature. In this paper, we develop P-FARFAR2: a parallel enhancement of FARFAR2 that increases its ability to assemble low-energy structures via multithreaded exploration of random configurations in a greedy manner. This strategy, appearing in the literature under the term "parallel mechanism", is made viable through two measures: first, the synchronization window is coarsened to several Monte Carlo cycles; second, all but one of the threads are differentiated as auxiliary and set to perform a weakened version of the problem. Following empirical analysis on a diverse range of RNA structures, we report achieving statistical significance in lowering the energy levels of ensuing samples. And consequently, despite the moderate-to-weak correlation between energy levels and prediction accuracy, this achievement happens to propagate to accuracy measurements.

RNAJP: enhanced RNA 3D structure predictions with non-canonical interactions and global topology sampling

RNA 3D structures are critical for understanding their functions. However, only a limited number of RNA structures have been experimentally solved, so computational prediction methods are highly desirable. Nevertheless, accurate prediction of RNA 3D structures, especially those containing multiway junctions, remains a significant challenge, mainly due to the complicated non-canonical base pairing and stacking interactions in the junction loops and the possible long-range interactions between loop structures. Here we present RNAJP ('RNA Junction Prediction'), a nucleotide- and helix-level coarse-grained model for the prediction of RNA 3D structures, particularly junction structures, from a given 2D structure. Through global sampling of the 3D arrangements of the helices in junctions using molecular dynamics simulations and in explicit consideration of non-canonical base pairing and base stacking interactions as well as long-range loop-loop interactions, the model can provide significantly improved predictions for multibranched junction structures than existing methods. Moreover, integrated with additional restraints from experiments, such as junction topology and long-range interactions, the model may serve as a useful structure generator for various applications.

The structure of a hibernating ribosome in a Lyme disease pathogen

The spirochete bacterial pathogen Borrelia ( Borreliella) burgdorferi ( Bbu ) affects more than 10% of the world population and causes Lyme disease in about half a million people in the US annually. Therapy for Lyme disease includes antibiotics that target the Bbu ribosome. We determined the structure of the Bbu 70S ribosome by single particle cryo-electron microscopy (cryo-EM) at a resolution of 2.9 Å, revealing its distinctive features. In contrast to a previous study suggesting that the single hibernation promoting factor protein present in Bbu (bbHPF) may not bind to its ribosome, our structure reveals a clear density for bbHPF bound to the decoding center of the small ribosomal 30S subunit. The 30S subunit has a non-annotated ribosomal protein, bS22, that has been found only in mycobacteria and Bacteroidetes so far. The protein bL38, recently discovered in Bacteroidetes, is also present in the Bbu large 50S ribosomal subunit. The protein bL37, previously seen only in mycobacterial ribosomes, is replaced by an N-terminal α-helical extension of uL30, suggesting that the two bacterial ribosomal proteins uL30 and bL37 may have evolved from one longer uL30 protein. The longer uL30 protein interacts with both the 23S rRNA and the 5S rRNA, is near the peptidyl transferase center (PTC), and could impart greater stability to this region. Its analogy to proteins uL30m and mL63 in mammalian mitochondrial ribosomes also suggests a plausible evolutionary pathway for expansion of protein content in mammalian mitochondrial ribosomes. Computational binding free energies are predicted for antibiotics, bound to the decoding center or PTC and are in clinical use for Lyme disease, that account for subtle distinctions in antibiotic-binding regions in the Bbu ribosome structure. Besides revealing unanticipated structural and compositional features for the Bbu ribosome, our study thus provides groundwork to enable ribosome-targeted antibiotic design for more effective treatment of Lyme disease.

Recognizing the power of machine learning and other computational methods to accelerate progress in small molecule targeting of RNA

RNA structures regulate a wide range of processes in biology and disease, yet small molecule chemical probes or drugs that can modulate these functions are rare. Machine learning and other computational methods are well poised to fill gaps in knowledge and overcome the inherent challenges in RNA targeting, such as the dynamic nature of RNA and the difficulty of obtaining RNA high-resolution structures. Successful tools to date include principal component analysis, linear discriminate analysis, k-nearest neighbor, artificial neural networks, multiple linear regression, and many others. Employment of these tools has revealed critical factors for selective recognition in RNA:small molecule complexes, predictable differences in RNA- and protein-binding ligands, and quantitative structure activity relationships that allow the rational design of small molecules for a given RNA target. Herein we present our perspective on the value of using machine learning and other computation methods to advance RNA:small molecule targeting, including select examples and their validation as well as necessary and promising future directions that will be key to accelerate discoveries in this important field.

Targeting Tumor Necrosis Factor Receptor 1 with Selected Aptamers for Anti-Inflammatory Activity

Tumor necrosis factor-α (TNFα) inhibitors are widely used in treating autoimmune diseases like rheumatoid arthritis (RA). These inhibitors can presumably alleviate RA symptoms by blocking TNFα-TNF receptor 1 (TNFR1)-mediated pro-inflammatory signaling pathways. However, the strategy also interrupts the survival and reproduction functions conducted by TNFα-TNFR2 interaction and causes side effects. Thus, it is urgently needed to develop inhibitors that can selectively block TNFα-TNFR1 but not TNFα-TNFR2. Here, nucleic acid-based aptamers against TNFR1 are explored as potential anti-RA candidates. Through the systematic evolution of ligands by exponential enrichment (SELEX), two types of TNFR1-targeting aptamers were obtained, and their K_D values are approximately 100-300 nM. In silico analysis shows that the binding interface of aptamer-TNFR1 highly overlapped with natural TNFα-TNFR1 binding. On the cellular level, the aptamers can exert TNFα inhibitory activity by binding to TNFR1. The anti-inflammatory efficiencies of aptamers were assessed and further enhanced using divalent aptamer constructs. These findings provide a new strategy to block TNFR1 for potential anti-RA treatment precisely.

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.

3D Structural Analysis of Long Noncoding RNA by Small Angle X-ray Scattering and Computational Modeling

Small angle X-ray scattering (SAXS) has been widely applied as an enabling integrative technique for comprehensive analysis of the structure of biomacromolecules by multiple, complementary techniques in solution. SAXS in combination with computational modeling can be a powerful strategy bridging the secondary and 3D structural analysis of large RNAs, including the long noncoding RNAs (lncRNA). Here, we outline the major procedures and techniques in the combined use of SAXS and computational modeling for 3D structural characterization of a lncRNA, the subgenomic flaviviral RNA from Zika virus. lncRNA production and purification, RNA buffer and sample preparation for SAXS experiments, SAXS data collection and analysis, SAXS-aided RNA 3D structure prediction, and computational modeling are described.

Combining Experimental Restraints and RNA 3D Structure Prediction in RNA Nanotechnology

Precise RNA tertiary structure prediction can aid in the design of RNA nanoparticles. However, most existing RNA tertiary structure prediction methods are limited to small RNAs with relatively simple secondary structures. Large RNA molecules usually have complex secondary structures, including multibranched loops and pseudoknots, allowing for highly flexible RNA geometries and multiple stable states. Various experiments and bioinformatics analyses can often provide information about the distance between atoms (or residues) in RNA, which can be used to guide the prediction of RNA tertiary structure. In this chapter, we will introduce a platform, iFoldNMR, that can incorporate non-exchangeable imino protons resonance data from NMR as restraints for RNA 3D structure prediction. We also introduce an algorithm, DVASS, which optimizes distance restraints for better RNA 3D structure prediction.

Sequence-sensitive elastic network captures dynamical features necessary for miR-125a maturation

The Elastic Network Contact Model (ENCoM) is a coarse-grained normal mode analysis (NMA) model unique in its all-atom sensitivity to the sequence of the studied macromolecule and thus to the effect of mutations. We adapted ENCoM to simulate the dynamics of ribonucleic acid (RNA) molecules, benchmarked its performance against other popular NMA models and used it to study the 3D structural dynamics of human microRNA miR-125a, leveraging high-throughput experimental maturation efficiency data of over 26 000 sequence variants. We also introduce a novel way of using dynamical information from NMA to train multivariate linear regression models, with the purpose of highlighting the most salient contributions of dynamics to function. ENCoM has a similar performance profile on RNA than on proteins when compared to the Anisotropic Network Model (ANM), the most widely used coarse-grained NMA model; it has the advantage on predicting large-scale motions while ANM performs better on B-factors prediction. A stringent benchmark from the miR-125a maturation dataset, in which the training set contains no sequence information in common with the testing set, reveals that ENCoM is the only tested model able to capture signal beyond the sequence. This ability translates to better predictive power on a second benchmark in which sequence features are shared between the train and test sets. When training the linear regression model using all available data, the dynamical features identified as necessary for miR-125a maturation point to known patterns but also offer new insights into the biogenesis of microRNAs. Our novel approach combining NMA with multivariate linear regression is generalizable to any macromolecule for which relatively high-throughput mutational data is available.

Advances in RNA 3D Structure Prediction

RNA molecules carry out various cellular functions, and understanding the mechanisms behind their functions requires the knowledge of their 3D structures. Different types of computational methods have been developed to model RNA 3D structures over the past decade. These methods were widely used by researchers although their performance needs to be further improved. Recently, along with these traditional methods, machine-learning techniques have been increasingly applied to RNA 3D structure prediction and show significant improvement in performance. Here we shall give a brief review of the traditional methods and recent related advances in machine-learning approaches for RNA 3D structure prediction.

Structures of the eukaryotic ribosome and its translational states in situ

Ribosomes translate genetic information into primary structure. During translation, various cofactors transiently bind to the ribosome that undergoes prominent conformational and structural changes. Different translational states of ribosomes have been well characterized in vitro. However, to which extent the known translational states are representative of the native situation inside cells has thus far only been addressed in prokaryotes. Here, we apply cryo-electron tomography to cryo-FIB milled Dictyostelium discoideum cells combined with subtomogram averaging and classification. We obtain an in situ structure that is locally resolved up to 3 Angstrom, the distribution of eukaryotic ribosome translational states, and unique arrangement of rRNA expansion segments. Our work demonstrates the use of in situ structural biology techniques for identifying distinct ribosome states within the cellular environment.

Computer Aided Development of Nucleic Acid Applications in Nanotechnologies

Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.

HDOCK update for modeling protein-RNA/DNA complex structures

Protein-nucleic acid interactions are involved in various cellular processes. Therefore, determining the structures of protein-nucleic acid complexes can provide insights into the mechanisms of the interactions and thus guide the rational drug design to modulate these interactions. Due to the high cost and technical difficulties of solving complex structures experimentally, computational modeling such as molecular docking has been playing an important role in the study of molecular interactions. In order to make it easier for researchers to obtain biomolecular complex structures through molecular docking, we developed the HDOCK server for protein-protein and protein-RNA/DNA docking (accessed at http://hdock.phys.hust.edu.cn/). Since its first release in 2017, HDOCK has been widely used in the scientific community. As nucleic acids may include single-stranded (ss) RNA/DNA and double-stranded (ds) RNA/DNA, we now present an updated version of HDOCK, which offers new options for structural modeling of ssRNA, ssDNA, dsRNA, and dsDNA. We hope this update will better help the scientific community solve important biological problems, thereby advancing the field. In this article, we describe the general protocol of HDOCK with emphasis on the new functions on RNA/DNA modeling. Several application examples are also given to illustrate the usage of the new functions.

NACDDB: Nucleic Acid Circular Dichroism Database

The Nucleic Acid Circular Dichroism Database (NACDDB) is a public repository that archives and freely distributes circular dichroism (CD) and synchrotron radiation CD (SRCD) spectral data about nucleic acids, and the associated experimental metadata, structural models, and links to literature. NACDDB covers CD data for various nucleic acid molecules, including DNA, RNA, DNA/RNA hybrids, and various nucleic acid derivatives. The entries are linked to primary sequence and experimental structural data, as well as to the literature. Additionally, for all entries, 3D structure models are provided. All entries undergo expert validation and curation procedures to ensure completeness, consistency, and quality of the data included. The NACDDB is open for submission of the CD data for nucleic acids. NACDDB is available at: https://genesilico.pl/nacddb/.

rRNA expansion segment 7 in eukaryotes: from Signature Fold to tentacles

The ribosomal core is universally conserved across the tree of life. However, eukaryotic ribosomes contain diverse rRNA expansion segments (ESs) on their surfaces. Sites of ES insertions are predicted from sites of insertion of micro-ESs in archaea. Expansion segment 7 (ES7) is one of the most diverse regions of the ribosome, emanating from a short stem loop and ranging to over 750 nucleotides in mammals. We present secondary and full-atom 3D structures of ES7 from species spanning eukaryotic diversity. Our results are based on experimental 3D structures, the accretion model of ribosomal evolution, phylogenetic relationships, multiple sequence alignments, RNA folding algorithms and 3D modeling by RNAComposer. ES7 contains a distinct motif, the 'ES7 Signature Fold', which is generally invariant in 2D topology and 3D structure in all eukaryotic ribosomes. We establish a model in which ES7 developed over evolution through a series of elementary and recursive growth events. The data are sufficient to support an atomic-level accretion path for rRNA growth. The non-monophyletic distribution of some ES7 features across the phylogeny suggests acquisition via convergent processes. And finally, illustrating the power of our approach, we constructed the 2D and 3D structure of the entire LSU rRNA of Mus musculus.

Visualizing translation dynamics at atomic detail inside a bacterial cell

Translation is the fundamental process of protein synthesis and is catalysed by the ribosome in all living cells1. Here we use advances in cryo-electron tomography and sub-tomogram analysis2,3 to visualize the structural dynamics of translation inside the bacterium Mycoplasma pneumoniae. To interpret the functional states in detail, we first obtain a high-resolution in-cell average map of all translating ribosomes and build an atomic model for the M. pneumoniae ribosome that reveals distinct extensions of ribosomal proteins. Classification then resolves 13 ribosome states that differ in their conformation and composition. These recapitulate major states that were previously resolved in vitro, and reflect intermediates during active translation. On the basis of these states, we animate translation elongation inside native cells and show how antibiotics reshape the cellular translation landscapes. During translation elongation, ribosomes often assemble in defined three-dimensional arrangements to form polysomes4. By mapping the intracellular organization of translating ribosomes, we show that their association into polysomes involves a local coordination mechanism that is mediated by the ribosomal protein L9. We propose that an extended conformation of L9 within polysomes mitigates collisions to facilitate translation fidelity. Our work thus demonstrates the feasibility of visualizing molecular processes at atomic detail inside cells.

3dDNA: A Computational Method of Building DNA 3D Structures

Considerable progress has been made in the prediction methods of 3D structures of RNAs. In contrast, no such methods are available for DNAs. The determination of 3D structures of the latter is also increasingly needed for understanding their functions and designing new DNA molecules. Since the number of experimental structures of DNA is limited at present, here, we propose a computational and template-based method, 3dDNA, which combines DNA and RNA template libraries to predict DNA 3D structures. It was benchmarked on three test sets with different numbers of chains, and the results show that 3dDNA can predict DNA 3D structures with a mean RMSD of about 2.36 Å for those with one or two chains and fewer than 4 Å with three or more chains.

Protein DEK and DTA Aptamers: Insight Into the Interaction Mechanisms and the Computational Aptamer Design

By blocking the DEK protein, DEK-targeted aptamers (DTAs) can reduce the formation of neutrophil extracellular traps (NETs) to reveal a strong anti-inflammatory efficacy in rheumatoid arthritis. However, the poor stability of DTA has greatly limited its clinical application. Thus, in order to design an aptamer with better stability, DTA was modified by methoxy groups (DTA_OMe) and then the exact DEK–DTA interaction mechanisms were explored through theoretical calculations. The corresponding 2′-OCH3-modified nucleotide force field was established and the molecular dynamics (MD) simulations were performed. It was proved that the 2′-OCH3-modification could definitely enhance the stability of DTA on the premise of comparative affinity. Furthermore, the electrostatic interaction contributed the most to the binding of DEK–DTA, which was the primary interaction to maintain stability, in addition to the non-specific interactions between positively-charged residues (e.g., Lys and Arg) of DEK and the negatively-charged phosphate backbone of aptamers. The H-bond network analysis reminded that eight bases could be mutated to probably enhance the affinity of DTA_OMe. Therein, replacing the 29th base from cytosine to thymine of DTA_OMe was theoretically confirmed to be with the best affinity and even better stability. These research studies imply to be a promising new aptamer design strategy for the treatment of inflammatory arthritis.

Recent advances in RNA structurome

RNA structures are essential to support RNA functions and regulation in various biological processes. Recently, a range of novel technologies have been developed to decode genome-wide RNA structures and novel modes of functionality across a wide range of species. In this review, we summarize key strategies for probing the RNA structurome and discuss the pros and cons of representative technologies. In particular, these new technologies have been applied to dissect the structural landscape of the SARS-CoV-2 RNA genome. We also summarize the functionalities of RNA structures discovered in different regulatory layers-including RNA processing, transport, localization, and mRNA translation-across viruses, bacteria, animals, and plants. We review many versatile RNA structural elements in the context of different physiological and pathological processes (e.g., cell differentiation, stress response, and viral replication). Finally, we discuss future prospects for RNA structural studies to map the RNA structurome at higher resolution and at the single-molecule and single-cell level, and to decipher novel modes of RNA structures and functions for innovative applications.

rna-tools.online: a Swiss army knife for RNA 3D structure modeling workflow

Significant improvements have been made in the efficiency and accuracy of RNA 3D structure prediction methods in recent years; however, many tools developed in the field stay exclusive to only a few bioinformatic groups. To perform a complete RNA 3D structure modeling analysis as proposed by the RNA-Puzzles community, researchers must familiarize themselves with a quite complex set of tools. In order to facilitate the processing of RNA sequences and structures, we previously developed the rna-tools package. However, using rna-tools requires the installation of a mixture of libraries and tools, basic knowledge of the command line and the Python programming language. To provide an opportunity for the broader community of biologists to take advantage of the new developments in RNA structural biology, we developed rna-tools.online. The web server provides a user-friendly platform to perform many standard analyses required for the typical modeling workflow: 3D structure manipulation and editing, structure minimization, structure analysis, quality assessment, and comparison. rna-tools.online supports biologists to start benefiting from the maturing field of RNA 3D structural bioinformatics and can be used for educational purposes. The web server is available at https://rna-tools.online.

Combining Coarse-Grained Simulations and Single Molecule Analysis Reveals a Three-State Folding Model of the Guanidine-II Riboswitch

Riboswitch RNAs regulate gene expression by conformational changes induced by environmental conditions and specific ligand binding. The guanidine-II riboswitch is proposed to bind the small molecule guanidinium and to subsequently form a kissing loop interaction between the P1 and P2 hairpins. While an interaction was shown for isolated hairpins in crystallization and electron paramagnetic resonance experiments, an intrastrand kissing loop formation has not been demonstrated. Here, we report the first evidence of this interaction in cis in a ligand and Mg²⁺ dependent manner. Using single-molecule FRET spectroscopy and detailed structural information from coarse-grained simulations, we observe and characterize three interconvertible states representing an open and kissing loop conformation as well as a novel Mg²⁺ dependent state for the guanidine-II riboswitch from E. coli. The results further substantiate the proposed switching mechanism and provide detailed insight into the regulation mechanism for the guanidine-II riboswitch class. Combining single molecule experiments and coarse-grained simulations therefore provides a promising perspective in resolving the conformational changes induced by environmental conditions and to yield molecular insights into RNA regulation.

Evaluation of the stereochemical quality of predicted RNA 3D models in the RNA-Puzzles submissions

In silico prediction is a well-established approach to derive a general shape of an RNA molecule based on its sequence or secondary structure. This paper reports an analysis of the stereochemical quality of the RNA three-dimensional models predicted using dedicated computer programs. The stereochemistry of 1052 RNA 3D structures, including 1030 models predicted by fully automated and human-guided approaches within 22 RNA-Puzzles challenges and reference structures, is analyzed. The evaluation is based on standards of RNA stereochemistry that the Protein Data Bank requires from deposited experimental structures. Deviations from standard bond lengths and angles, planarity, or chirality are quantified. A reduction in the number of such deviations should help in the improvement of RNA 3D structure modeling approaches.

Prediction of the RNA Tertiary Structure Based on a Random Sampling Strategy and Parallel Mechanism

Background: Macromolecule structure prediction remains a fundamental challenge of bioinformatics. Over the past several decades, the Rosetta framework has provided solutions to diverse challenges in computational biology. However, it is challenging to model RNA tertiary structures effectively when the de novo modeling of RNA involves solving a well-defined small puzzle. Methods: In this study, we introduce a stepwise Monte Carlo parallelization (SMCP) algorithm for RNA tertiary structure prediction. Millions of conformations were randomly searched using the Monte Carlo algorithm and stepwise ansatz hypothesis, and SMCP uses a parallel mechanism for efficient sampling. Moreover, to achieve better prediction accuracy and completeness, we judged and processed the modeling results. Results: A benchmark of nine single-stranded RNA loops drawn from riboswitches establishes the general ability of the algorithm to model RNA with high accuracy and integrity, including six motifs that cannot be solved by knowledge mining-based modeling algorithms. Experimental results show that the modeling accuracy of the SMCP algorithm is up to 0.14 Å, and the modeling integrity on this benchmark is extremely high. Conclusion: SMCP is an ab initio modeling algorithm that substantially outperforms previous algorithms in the Rosetta framework, especially in improving the accuracy and completeness of the model. It is expected that the work will provide new research ideas for macromolecular structure prediction in the future. In addition, this work will provide theoretical basis for the development of the biomedical field.

RNAStat: An Integrated Tool for Statistical Analysis of RNA 3D Structures

The 3D architectures of RNAs are essential for understanding their cellular functions. While an accurate scoring function based on the statistics of known RNA structures is a key component for successful RNA structure prediction or evaluation, there are few tools or web servers that can be directly used to make comprehensive statistical analysis for RNA 3D structures. In this work, we developed RNAStat, an integrated tool for making statistics on RNA 3D structures. For given RNA structures, RNAStat automatically calculates RNA structural properties such as size and shape, and shows their distributions. Based on the RNA structure annotation from DSSR, RNAStat provides statistical information of RNA secondary structure motifs including canonical/non-canonical base pairs, stems, and various loops. In particular, the geometry of base-pairing/stacking can be calculated in RNAStat by constructing a local coordinate system for each base. In addition, RNAStat also supplies the distribution of distance between any atoms to the users to help build distance-based RNA statistical potentials. To test the usability of the tool, we established a non-redundant RNA 3D structure dataset, and based on the dataset, we made a comprehensive statistical analysis on RNA structures, which could have the guiding significance for RNA structure modeling. The python code of RNAStat, the dataset used in this work, and corresponding statistical data files are freely available at GitHub (https://github.com/RNA-folding-lab/RNAStat).

rsRNASP: A residue-separation-based statistical potential for RNA 3D structure evaluation

Knowledge-based statistical potentials have been shown to be rather effective in protein 3-dimensional (3D) structure evaluation and prediction. Recently, several statistical potentials have been developed for RNA 3D structure evaluation, while their performances are either still at a low level for the test datasets from structure prediction models or dependent on the "black-box" process through neural networks. In this work, we have developed an all-atom distance-dependent statistical potential based on residue separation for RNA 3D structure evaluation, namely rsRNASP, which is composed of short- and long-ranged potentials distinguished by residue separation. The extensive examinations against available RNA test datasets show that rsRNASP has apparently higher performance than the existing statistical potentials for the realistic test datasets with large RNAs from structure prediction models, including the newly released RNA-Puzzles dataset, and is comparable to the existing top statistical potentials for the test datasets with small RNAs or near-native decoys. In addition, rsRNASP is superior to RNA3DCNN, a recently developed scoring function through 3D convolutional neural networks. rsRNASP and the relevant databases are available to the public.

Constrained peptides mimic a viral suppressor of RNA silencing

The design of high-affinity, RNA-binding ligands has proven very challenging. This is due to the unique structural properties of RNA, often characterized by polar surfaces and high flexibility. In addition, the frequent lack of well-defined binding pockets complicates the development of small molecule binders. This has triggered the search for alternative scaffolds of intermediate size. Among these, peptide-derived molecules represent appealing entities as they can mimic structural features also present in RNA-binding proteins. However, the application of peptidic RNA-targeting ligands is hampered by a lack of design principles and their inherently low bio-stability. Here, the structure-based design of constrained α-helical peptides derived from the viral suppressor of RNA silencing, TAV2b, is described. We observe that the introduction of two inter-side chain crosslinks provides peptides with increased α-helicity and protease stability. One of these modified peptides (B3) shows high affinity for double-stranded RNA structures including a palindromic siRNA as well as microRNA-21 and its precursor pre-miR-21. Notably, B3 binding to pre-miR-21 inhibits Dicer processing in a biochemical assay. As a further characteristic this peptide also exhibits cellular entry. Our findings show that constrained peptides can efficiently mimic RNA-binding proteins rendering them potentially useful for the design of bioactive RNA-targeting ligands.

CoCoNet-boosting RNA contact prediction by convolutional neural networks

Co-evolutionary models such as direct coupling analysis (DCA) in combination with machine learning (ML) techniques based on deep neural networks are able to predict accurate protein contact or distance maps. Such information can be used as constraints in structure prediction and massively increase prediction accuracy. Unfortunately, the same ML methods cannot readily be applied to RNA as they rely on large structural datasets only available for proteins. Here, we demonstrate how the available smaller data for RNA can be used to improve prediction of RNA contact maps. We introduce an algorithm called CoCoNet that is based on a combination of a Coevolutionary model and a shallow Convolutional Neural Network. Despite its simplicity and the small number of trained parameters, the method boosts the positive predictive value (PPV) of predicted contacts by about 70% with respect to DCA as tested by cross-validation of about eighty RNA structures. However, the direct inclusion of the CoCoNet contacts in 3D modeling tools does not result in a proportional increase of the 3D RNA structure prediction accuracy. Therefore, we suggest that the field develops, in addition to contact PPV, metrics which estimate the expected impact for 3D structure modeling tools better. CoCoNet is freely available and can be found at https://github.com/KIT-MBS/coconet.

Structure of the translating Neurospora ribosome arrested by cycloheximide

Ribosomes translate RNA into proteins. The protein synthesis inhibitor cycloheximide (CHX) is widely used to inhibit eukaryotic ribosomes engaged in translation elongation. However, the lack of structural data for actively translating polyribosomes stalled by CHX leaves unanswered the question of which elongation step is inhibited. We elucidated CHX's mechanism of action based on the cryo-electron microscopy structure of actively translating Neurospora crassa ribosomes bound with CHX at 2.7-Å resolution. The ribosome structure from this filamentous fungus contains clearly resolved ribosomal protein eL28, like higher eukaryotes but unlike budding yeast, which lacks eL28. Despite some differences in overall structures, the ribosomes from Neurospora, yeast, and humans all contain a highly conserved CHX binding site. We also sequenced classic Neurospora CHX-resistant alleles. These mutations, including one at a residue not previously observed to affect CHX resistance in eukaryotes, were in the large subunit proteins uL15 and eL42 that are part of the CHX-binding pocket. In addition to A-site transfer RNA (tRNA), P-site tRNA, messenger RNA, and CHX that are associated with the translating N. crassa ribosome, spermidine is present near the CHX binding site close to the E site on the large subunit. The tRNAs in the peptidyl transferase center are in the A/A site and the P/P site. The nascent peptide is attached to the A-site tRNA and not to the P-site tRNA. The structural and functional data obtained show that CHX arrests the ribosome in the classical PRE translocation state and does not interfere with A-site reactivity.

Modeling Noncanonical RNA Base Pairs by a Coarse-Grained IsRNA2 Model

Noncanonical base pairs contribute crucially to the three-dimensional architecture of large RNA molecules; however, how to accurately model them remains an open challenge in RNA 3D structure prediction. Here, we report a promising coarse-grained (CG) IsRNA2 model to predict noncanonical base pairs in large RNAs through molecular dynamics simulations. By introducing a five-bead per nucleotide CG representation to reserve the three interacting edges of nucleobases, IsRNA2 accurately models various base-pairing interactions, including both canonical and noncanonical base pairs. A benchmark test indicated that IsRNA2 achieves a comparable performance to the atomic model in de novo modeling of noncanonical RNA structures. In addition, IsRNA2 was able to refine the 3D structure predictions for large RNAs in RNA-puzzle challenges. Finally, the graphics processing unit acceleration was introduced to speed up the sampling efficiency in IsRNA2 for very large RNA molecules. Therefore, the CG IsRNA2 model reported here offers a reliable approach to predict the structures and dynamics of large RNAs.