Capturing RNA Folding Free Energy with Coarse-Grained Molecular Dynamics Simulations

rCGMM: A Coarse-Grained Force Field Embedding Elastic Network for Studying Small Noncoding RNA Dynamics

Short noncoding RNA molecules play significant roles in catalysis, biological regulation, and disease pathways. Their assessment through sequence-based approaches has been a challenge, compounded by the significant structural flexibility accrued from six free backbone torsions per nucleotide. To efficiently study the structure and dynamics of an extensive repertoire of these molecules in a high throughput mode, we have built a coarse-grained force field using one, two, three, and four pseudoatoms to represent the phosphate, sugar, pyrimidines, and purines, respectively. The Boltzmann inversion method was applied to structures of 5 piRNA, 8 miRNA, and 13 siRNA from the Nucleic Acid Database (NDB) to estimate the initial force field parameters and iteratively optimized through 1 μs molecular dynamics run by comparing against an equivalent all-atom simulation using the CHARMM36 force field. We applied an elastic net to model the hydrogen bond network stabilizing the local structure for double-stranded cases. A spine using pseudoatoms was calculated for the same from the coarse-grain beads, and all beads within a threshold radial distance were constrained using soft distance potentials. Lennard-Jones and Coulomb's potential function modeled the nonbonded interaction. Benchmarks on 26 molecules compared through root-mean-square deviation graphs against all-atom simulation show close concurrence for single- and double-stranded small noncoding RNA molecules. The rCGMM force field is available for download at https://github.com/majumderS/rCGMM.

The Unpaved Road of Non-Coding RNA Structure-Function Relationships: Current Knowledge, Available Methodologies, and Future Trends

The genomes from complex eukaryotes are enriched in non-coding genes whose transcription products (non-coding RNAs) are involved in the regulation of genomic output at different levels. Non-coding RNA action is predominantly driven by sequence and structural motifs that interact with specific functional partners. Despite the exponential growth in primary RNA sequence data facilitated by next-generation sequencing studies, the availability of tridimensional RNA data is comparatively more limited. The subjacent reasons for this relative lack of information regarding RNA structure are related to the specific chemical nature of RNA molecules and the limitations of the currently available methods for structural characterization of biomolecules. In this review, we describe and analyze the different structural motifs involved in non-coding RNA function and the wet-lab and computational methods used to characterize their structure-function relationships, highlighting the current need for detailed structural studies to explore the molecular determinants of non-coding RNA function.

Computational biology and artificial intelligence in mRNA vaccine design for cancer immunotherapy

Messenger RNA (mRNA) vaccines offer an adaptable and scalable platform for cancer immunotherapy, requiring optimal design to elicit a robust and targeted immune response. Recent advancements in bioinformatics and artificial intelligence (AI) have significantly enhanced the design, prediction, and optimization of mRNA vaccines. This paper reviews technologies that streamline mRNA vaccine development, from genomic sequencing to lipid nanoparticle (LNP) formulation. We discuss how accurate predictions of neoantigen structures guide the design of mRNA sequences that effectively target immune and cancer cells. Furthermore, we examine AI-driven approaches that optimize mRNA-LNP formulations, enhancing delivery and stability. These technological innovations not only improve vaccine design but also enhance pharmacokinetics and pharmacodynamics, offering promising avenues for personalized cancer immunotherapy.

Conformational Properties of Poly(A)-Binding Protein Complexed with Poly(A) RNA

Microscopic understanding of protein-RNA interactions is important for different biological activities, such as RNA transport, translation, splicing, silencing, etc. Polyadenine (Poly(A)) binding proteins (PABPs) make up a class of regulatory proteins that play critical roles in protecting the poly(A) tails of cellular mRNAs from nuclease degradation. In this work, we performed molecular dynamics simulations to investigate the conformational modifications of human PABP protein and poly(A) RNA that occur during complexation. It is demonstrated that the intermediate linker domain of the protein transforms from a disordered coil-like structure to a helical form during the recognition process, leading to the formation of the complex. On the other hand, disordered collapsed coil-like RNA on complexation has been found to transform into a rigid extended conformation. Importantly, the binding free energy calculation showed that the thermodynamic stability of the complex is primarily guided by favorable hydrophobic interactions between the protein and the RNA.

Quantifying RNA structures and interactions with a unified reduced chain representation model

RNA is a pivotal molecule that plays critical roles in various cellular processes. Quantifying RNA structures and interactions is essential to understanding RNA function and developing RNA-based therapeutics. Using a unified five-bead model and a non-redundant database, this paper investigates the structural features and interactions of five commonly occurring RNA motifs, i.e., double-stranded helices, hairpin loops, internal/bulge loops, multi-branched junctions, and single-stranded terminal tails. Analyzing detailed distributions of RNA local structural features and base-base interactions reveals a preference for helical structures in both local backbone structures and base orientations. The interactions between adjacent bases exhibit motif-specific and sequence-dependent characteristics, reflecting the distinct topological constraints imposed by different loop-helix connection modes and the varying pairing and stacking interactions among different sequences. These findings shed light on the stability of RNA helices, emphasizing their significance in providing dominant base pairing and stacking interactions for RNA structures and stability. The four non-helix motifs encompass unpaired nucleotide loops and exhibit diverse base-base interactions, contributing to the structural diversity observed in RNA. Overall, the complexity of RNA structure arises from the intricate interplay of base-base interactions.

Predicting the Effect of Single Mutations on Protein Stability and Binding with Respect to Types of Mutations

The development of methods and algorithms to predict the effect of mutations on protein stability, protein-protein interaction, and protein-DNA/RNA binding is necessitated by the needs of protein engineering and for understanding the molecular mechanism of disease-causing variants. The vast majority of the leading methods require a database of experimentally measured folding and binding free energy changes for training. These databases are collections of experimental data taken from scientific investigations typically aimed at probing the role of particular residues on the above-mentioned thermodynamic characteristics, i.e., the mutations are not introduced at random and do not necessarily represent mutations originating from single nucleotide variants (SNV). Thus, the reported performance of the leading algorithms assessed on these databases or other limited cases may not be applicable for predicting the effect of SNVs seen in the human population. Indeed, we demonstrate that the SNVs and non-SNVs are not equally presented in the corresponding databases, and the distribution of the free energy changes is not the same. It is shown that the Pearson correlation coefficients (PCCs) of folding and binding free energy changes obtained in cases involving SNVs are smaller than for non-SNVs, indicating that caution should be used in applying them to reveal the effect of human SNVs. Furthermore, it is demonstrated that some methods are sensitive to the chemical nature of the mutations, resulting in PCCs that differ by a factor of four across chemically different mutations. All methods are found to underestimate the energy changes by roughly a factor of 2.

Predicting 3D structures and stabilities for complex RNA pseudoknots in ion solutions

RNA pseudoknots are a kind of important tertiary motif, and the structures and stabilities of pseudoknots are generally critical to the biological functions of RNAs with the motifs. In this work, we have carefully refined our previously developed coarse-grained model with salt effect through involving a new coarse-grained force field and a replica-exchange Monte Carlo algorithm, and employed the model to predict structures and stabilities of complex RNA pseudoknots in ion solutions beyond minimal H-type pseudoknots. Compared with available experimental data, the newly refined model can successfully predict 3D structures from sequences for the complex RNA pseudoknots including SARS-CoV-2 programming-1 ribosomal frameshifting element and Zika virus xrRNA, and can reliably predict the thermal stabilities of RNA pseudoknots with various sequences and lengths over broad ranges of monovalent/divalent salts. In addition, for complex pseudoknots including SARS-CoV-2 frameshifting element, our analyses show that their thermally unfolding pathways are mainly dependent on the relative stabilities of unfolded intermediate states, in analogy to those of minimal H-type pseudoknots.

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.

cgRNASP: coarse-grained statistical potentials with residue separation for RNA structure evaluation

Knowledge-based statistical potentials are very important for RNA 3-dimensional (3D) structure prediction and evaluation. In recent years, various coarse-grained (CG) and all-atom models have been developed for predicting RNA 3D structures, while there is still lack of reliable CG statistical potentials not only for CG structure evaluation but also for all-atom structure evaluation at high efficiency. In this work, we have developed a series of residue-separation-based CG statistical potentials at different CG levels for RNA 3D structure evaluation, namely cgRNASP, which is composed of long-ranged and short-ranged interactions by residue separation. Compared with the newly developed all-atom rsRNASP, the short-ranged interaction in cgRNASP was involved more subtly and completely. Our examinations show that, the performance of cgRNASP varies with CG levels and compared with rsRNASP, cgRNASP has similarly good performance for extensive types of test datasets and can have slightly better performance for the realistic dataset-RNA-Puzzles dataset. Furthermore, cgRNASP is strikingly more efficient than all-atom statistical potentials/scoring functions, and can be apparently superior to other all-atom statistical potentials and scoring functions trained from neural networks for the RNA-Puzzles dataset. cgRNASP is available at https://github.com/Tan-group/cgRNASP.

TPP Riboswitch Populates Holo-Form-like Structure Even in the Absence of Cognate Ligand at High Mg2+ Concentration

Riboswitches are noncoding RNA that regulate gene expression by folding into specific three-dimensional structures (holo-form) upon binding by their cognate ligand in the presence of Mg²⁺. Riboswitch functioning is also hypothesized to be under kinetic control requiring large cognate ligand concentrations. We ask the question under thermodynamic conditions, can the riboswitches populate structures similar to the holo-form only in the presence of Mg²⁺ and absence of cognate ligand binding. We addressed this question using thiamine pyrophosphate (TPP) riboswitch as a model system and computer simulations using a coarse-grained model for RNA. The folding free energy surface (FES) shows that with the initial increase in Mg²⁺ concentration ([Mg²⁺]), the aptamer domain (AD) of TPP riboswitch undergoes a barrierless collapse in its dimensions. On further increase in [Mg²⁺], intermediates separated by barriers appear on the FES, and one of the intermediates has a TPP ligand-binding competent structure. We show that site-specific binding of the Mg²⁺ aids in the formation of tertiary contacts. For [Mg²⁺] greater than physiological concentration, AD folds into a structure similar to the crystal structure of the TPP holo-form even in the absence of the TPP ligand. The folding kinetics shows that TPP AD populates an intermediate due to the misalignment of two arms present in the structure, which acts as a kinetic trap, leading to larger folding timescales. The predictions of the intermediate structures from the simulations are amenable for experimental verification.

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.

RNA Secondary Structures Regulate Adsorption of Fragments onto Flat Substrates

RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.

Biasing RNA Coarse-Grained Folding Simulations with Small-Angle X-ray Scattering Data

RNA molecules can easily adopt alternative structures in response to different environmental conditions. As a result, a molecule's energy landscape is rough and can exhibit a multitude of deep basins. In the absence of a high-resolution structure, small-angle X-ray scattering data (SAXS) can narrow down the conformational space available to the molecule and be used in conjunction with physical modeling to obtain high-resolution putative structures to be further tested by experiments. Because of the low resolution of these data, it is natural to implement the integration of SAXS data into simulations using a coarse-grained representation of the molecule, allowing for much wider searches and faster evaluation of SAXS theoretical intensity curves than with atomistic models. We present here the theoretical framework and the implementation of a simulation approach based on our coarse-grained model HiRE-RNA combined with SAXS evaluations "on-the-fly" leading the simulation toward conformations agreeing with the scattering data, starting from partially folded structures as the ones that can easily be obtained from secondary structure prediction-based tools. We show on three benchmark systems how our approach can successfully achieve high-resolution structures with remarkable similarity with the native structure recovering not only the overall shape, as imposed by SAXS data, but also the details of initially missing base pairs.

RNA 3D Structure Prediction Using Coarse-Grained Models

The three-dimensional (3D) structures of Ribonucleic acid (RNA) molecules are essential to understanding their various and important biological functions. However, experimental determination of the atomic structures is laborious and technically difficult. The large gap between the number of sequences and the experimentally determined structures enables the thriving development of computational approaches to modeling RNAs. However, computational methods based on all-atom simulations are intractable for large RNA systems, which demand long time simulations. Facing such a challenge, many coarse-grained (CG) models have been developed. Here, we provide a review of CG models for modeling RNA 3D structures, compare the performance of the different models, and offer insights into potential future developments.

IsRNA1: De Novo Prediction and Blind Screening of RNA 3D Structures

Modeling structures and functions of large ribonucleic acid (RNAs) especially with complicated topologies is highly challenging due to the inefficiency of large conformational sampling and the presence of complicated tertiary interactions. To address this problem, one highly promising approach is coarse-grained modeling. Here, following an iterative simulated reference state approach to decipher the correlations between different structural parameters, we developed a potent coarse-grained RNA model named as IsRNA1 for RNA studies. Molecular dynamics simulations in the IsRNA1 can predict the native structures of small RNAs from a sequence and fold medium-sized RNAs into near-native tertiary structures with the assistance of secondary structure constraints. A large-scale benchmark test on RNA 3D structure prediction shows that IsRNA1 exhibits improved performance for relatively large RNAs of complicated topologies, such as large stem-loop structures and structures containing long-range tertiary interactions. The advantages of IsRNA1 include the consideration of the correlations between the different structural variables, the appropriate characterization of canonical base-pairing and base-stacking interactions, and the better sampling for the backbone conformations. Moreover, a blind screening protocol was developed based on IsRNA1 to identify good structural models from a pool of candidates without prior knowledge of the native structures.

A Bayes-inspired theory for optimally building an efficient coarse-grained folding force field

Because of their potential utility in predicting conformational changes and assessing folding dynamics, coarse-grained (CG) RNA folding models are appealing for rapid characterization of RNA molecules. Previously, we reported the iterative simulated RNA reference state (IsRNA) method for parameterizing a CG force field for RNA folding, which consecutively updates the simulation force field to reflect marginal distributions of folding coordinates in the structure database and extract various energy terms. While the IsRNA model was validated by showing close agreement between the IsRNA-simulated and experimentally observed distributions, here, we expand our theoretical understanding of the model and, in doing so, improve the parameterization process to optimize the subset of included folding coordinates, which leads to accelerated simulations. Using statistical mechanical theory, we analyze the underlying, Bayesian concept that drives parameterization of the energy function, providing a general method for developing predictive, knowledge-based, polymer force fields on the basis of limited data. Furthermore, we propose an optimal parameterization procedure, based on the principal of maximum entropy.

Mg2+ Sensing by an RNA Fragment: Role of Mg2+-Coordinated Water Molecules

RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg²⁺, Ca²⁺, and K⁺) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg²⁺ solution. We show that the rigid network of Mg²⁺-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K⁺ and Ca²⁺ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.

In silico design and validation of high-affinity RNA aptamers targeting epithelial cellular adhesion molecule dimers

Nucleic acid aptamers hold great promise for therapeutic applications due to their favorable intrinsic properties, as well as high-throughput experimental selection techniques. Despite the utility of the systematic evolution of ligands by the exponential enrichment (SELEX) method for aptamer determination, complementary in silico aptamer design is highly sought after to facilitate virtual screening and increased understanding of important nucleic acid-protein interactions. Here, with a combined experimental and theoretical approach, we have developed two optimal epithelial cellular adhesion molecule (EpCAM) aptamers. Our structure-based in silico method first predicts their binding modes and then optimizes them for EpCAM with molecular dynamics simulations, docking, and free energy calculations. Our isothermal titration calorimetry experiments further confirm that the EpCAM aptamers indeed exhibit enhanced affinity over a previously patented nanomolar aptamer, EP23. Moreover, our study suggests that EP23 and the de novo designed aptamers primarily bind to EpCAM dimers (and not monomers, as hypothesized in previous published works), suggesting a paradigm for developing EpCAM-targeted therapies.

Small RNA Biosensor Design Strategy To Mitigate Off-Analyte Response

Several bottlenecks in the design of current sensor technologies for small noncoding RNA must be addressed. The small size of the sensors and the large number of other nucleotides that may have sequence similarity makes selectivity a real concern. Many of the current sensors have one strand with an exposed region called a toehold. The toehold serves as a place for the analyte nucleic acid strand to bind and initiate competitive displacement of sensors' secondary strands. Since the toehold region is not protected, any endogenous oligonucleotide sequences that are similar or only different by a few nucleic acids will interact with the toehold and cause false signals. To address sensor selectivity, we investigated how the toehold location in the sensor impacts the sensitivity and selectivity for the analyte of interest. We will discuss the differences in sensitivity and selectivity for a miR-146a-5p biosensor in the presence of different naturally occurring mismatch sequences. We found that altering the toehold location lowered the rate of the false signal from off-analyte microRNA by upward of 20 percentage points. Detection limits as low as 56 pM were observed when the sensor concentration was 5 nM. The findings herein are broadly applicable to other small and large RNAs as well as other types of sensing platforms.

Brief Comparison between Experimental and Computationally Generated Ensembles of RNA Dinucleotides

The computational modeling of RNA and its interactions is of crucial importance for the understanding of the wide variety of biological functions it performs. Among these approaches, several coarse-grained models employ structural databases to derive their energy functions or to define scoring functions for structure prediction purposes. In many cases, the parametrization is done by using as a reference a set of experimentally determined structures or data obtained from Molecular Dynamics simulations. Since the two choices are clearly different, we present here a brief comparison of the essential spaces of a set of conformations of two topologically connected nucleotides generated by these means. We find that when the nucleotides are embedded into a duplex, the essential spaces of both ensembles are quite restricted and do not differ substantially. Nevertheless, when the conformational space of a free dinucleoside monophosphate simulation is compared against a similar set obtained from the structural database, the differences of the essential spaces are considerable. In addition, we show a brief comparison of a specific distance between the nucleotides which correlates with the sugar pucker of the first nucleotide and analyze its distribution under similar conditions.

Theory and simulations for RNA folding in mixtures of monovalent and divalent cations

RNA molecules cannot fold in the absence of counterions. Experiments are typically performed in the presence of monovalent and divalent cations. How to treat the impact of a solution containing a mixture of both ion types on RNA folding has remained a challenging problem for decades. By exploiting the large concentration difference between divalent and monovalent ions used in experiments, we develop a theory based on the reference interaction site model (RISM), which allows us to treat divalent cations explicitly while keeping the implicit screening effect due to monovalent ions. Our theory captures both the inner shell and outer shell coordination of divalent cations to phosphate groups, which we demonstrate is crucial for an accurate calculation of RNA folding thermodynamics. The RISM theory for ion-phosphate interactions when combined with simulations based on a transferable coarse-grained model allows us to predict accurately the folding of several RNA molecules in a mixture containing monovalent and divalent ions. The calculated folding free energies and ion-preferential coefficients for RNA molecules (pseudoknots, a fragment of the rRNA, and the aptamer domain of the adenine riboswitch) are in excellent agreement with experiments over a wide range of monovalent and divalent ion concentrations. Because the theory is general, it can be readily used to investigate ion and sequence effects on DNA properties.

What is the best reference state for building statistical potentials in RNA 3D structure evaluation?

Knowledge-based statistical potentials have been shown to be efficient in protein structure evaluation/prediction, and the core difference between various statistical potentials is attributed to the choice of reference states. However, for RNA 3D structure evaluation, a comprehensive examination on reference states is still lacking. In this work, we built six statistical potentials based on six reference states widely used in protein structure evaluation, including averaging, quasi-chemical approximation, atom-shuffled, finite-ideal-gas, spherical-noninteracting, and random-walk-chain reference states, and we examined the six reference states against three RNA test sets including six subsets. Our extensive examinations show that, overall, for identifying native structures and ranking decoy structures, the finite-ideal-gas and random-walk-chain reference states are slightly superior to others, while for identifying near-native structures, there is only a slight difference between these reference states. Our further analyses show that the performance of a statistical potential is apparently dependent on the quality of the training set. Furthermore, we found that the performance of a statistical potential is closely related to the origin of test sets, and for the three realistic test subsets, the six statistical potentials have overall unsatisfactory performance. This work presents a comprehensive examination on the existing reference states and statistical potentials for RNA 3D structure evaluation.

Coarse-grained dynamic RNA titration simulations

Electrostatic interactions play a pivotal role in many biomolecular processes. The molecular organization and function in biological systems are largely determined by these interactions. Owing to the highly negative charge of RNA, the effect is expected to be more pronounced in this system. Moreover, RNA base pairing is dependent on the charge of the base, giving rise to alternative secondary and tertiary structures. The equilibrium between uncharged and charged bases is regulated by the solution pH, which is therefore a key environmental condition influencing the molecule's structure and behaviour. By means of constant-pH Monte Carlo simulations based on a fast proton titration scheme, coupled with the coarse-grained model HiRE-RNA, molecular dynamic simulations of RNA molecules at constant pH enable us to explore the RNA conformational plasticity at different pH values as well as to compute electrostatic properties as local pK _a values for each nucleotide.

Computational modeling of RNA 3D structure based on experimental data

RNA molecules are master regulators of cells. They are involved in a variety of molecular processes: they transmit genetic information, sense cellular signals and communicate responses, and even catalyze chemical reactions. As in the case of proteins, RNA function is dictated by its structure and by its ability to adopt different conformations, which in turn is encoded in the sequence. Experimental determination of high-resolution RNA structures is both laborious and difficult, and therefore the majority of known RNAs remain structurally uncharacterized. To address this problem, predictive computational methods were developed based on the accumulated knowledge of RNA structures determined so far, the physical basis of the RNA folding, and taking into account evolutionary considerations, such as conservation of functionally important motifs. However, all theoretical methods suffer from various limitations, and they are generally unable to accurately predict structures for RNA sequences longer than 100-nt residues unless aided by additional experimental data. In this article, we review experimental methods that can generate data usable by computational methods, as well as computational approaches for RNA structure prediction that can utilize data from experimental analyses. We outline methods and data types that can be potentially useful for RNA 3D structure modeling but are not commonly used by the existing software, suggesting directions for future development.

Bayesian selection for coarse-grained models of liquid water

The necessity for accurate and computationally efficient representations of water in atomistic simulations that can span biologically relevant timescales has born the necessity of coarse-grained (CG) modeling. Despite numerous advances, CG water models rely mostly on a-priori specified assumptions. How these assumptions affect the model accuracy, efficiency, and in particular transferability, has not been systematically investigated. Here we propose a data driven comparison and selection for CG water models through a Hierarchical Bayesian framework. We examine CG water models that differ in their level of coarse-graining, structure, and number of interaction sites. We find that the importance of electrostatic interactions for the physical system under consideration is a dominant criterion for the model selection. Multi-site models are favored, unless the effects of water in electrostatic screening are not relevant, in which case the single site model is preferred due to its computational savings. The charge distribution is found to play an important role in the multi-site model’s accuracy while the flexibility of the bonds/angles may only slightly improve the models. Furthermore, we find significant variations in the computational cost of these models. We present a data informed rationale for the selection of CG water models and provide guidance for future water model designs.

Chemically Accurate Relative Folding Stability of RNA Hairpins from Molecular Simulations

To benchmark RNA force fields, we compared the folding stabilities of three 12-nucleotide hairpin stem loops estimated by simulation to stabilities determined by experiment. We used umbrella sampling and a reaction coordinate of end-to-end (5' to 3' hydroxyl oxygen) distance to estimate the free energy change of the transition from the native conformation to a fully extended conformation with no hydrogen bonds between non-neighboring bases. Each simulation was performed four times using the AMBER FF99+bsc0+χ_OL3 force field, and each window, spaced at 1 Å intervals, was sampled for 1 μs, for a total of 552 μs of simulation. We compared differences in the simulated free energy changes to analogous differences in free energies from optical melting experiments using thermodynamic cycles where the free energy change between stretched and random coil sequences is assumed to be sequence-independent. The differences between experimental and simulated ΔΔ G° are, on average, 0.98 ± 0.66 kcal/mol, which is chemically accurate and suggests that analogous simulations could be used predictively. We also report a novel method to identify where replica free energies diverge along a reaction coordinate, thus indicating where additional sampling would most improve convergence. We conclude by discussing methods to more economically perform these simulations.

Modeling Structure, Stability, and Flexibility of Double-Stranded RNAs in Salt Solutions

Double-stranded (ds) RNAs play essential roles in many processes of cell metabolism. The knowledge of three-dimensional (3D) structure, stability, and flexibility of dsRNAs in salt solutions is important for understanding their biological functions. In this work, we further developed our previously proposed coarse-grained model to predict 3D structure, stability, and flexibility for dsRNAs in monovalent and divalent ion solutions through involving an implicit structure-based electrostatic potential. The model can make reliable predictions for 3D structures of extensive dsRNAs with/without bulge/internal loops from their sequences, and the involvement of the structure-based electrostatic potential and corresponding ion condition can improve the predictions for 3D structures of dsRNAs in ion solutions. Furthermore, the model can make good predictions for thermal stability for extensive dsRNAs over the wide range of monovalent/divalent ion concentrations, and our analyses show that the thermally unfolding pathway of dsRNA is generally dependent on its length as well as its sequence. In addition, the model was employed to examine the salt-dependent flexibility of a dsRNA helix, and the calculated salt-dependent persistence lengths are in good accordance with experiments.

Predicting 3D structure and stability of RNA pseudoknots in monovalent and divalent ion solutions

RNA pseudoknots are a kind of minimal RNA tertiary structural motifs, and their three-dimensional (3D) structures and stability play essential roles in a variety of biological functions. Therefore, to predict 3D structures and stability of RNA pseudoknots is essential for understanding their functions. In the work, we employed our previously developed coarse-grained model with implicit salt to make extensive predictions and comprehensive analyses on the 3D structures and stability for RNA pseudoknots in monovalent/divalent ion solutions. The comparisons with available experimental data show that our model can successfully predict the 3D structures of RNA pseudoknots from their sequences, and can also make reliable predictions for the stability of RNA pseudoknots with different lengths and sequences over a wide range of monovalent/divalent ion concentrations. Furthermore, we made comprehensive analyses on the unfolding pathway for various RNA pseudoknots in ion solutions. Our analyses for extensive pseudokonts and the wide range of monovalent/divalent ion concentrations verify that the unfolding pathway of RNA pseudoknots is mainly dependent on the relative stability of unfolded intermediate states, and show that the unfolding pathway of RNA pseudoknots can be significantly modulated by their sequences and solution ion conditions.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

IsRNA: An Iterative Simulated Reference State Approach to Modeling Correlated Interactions in RNA Folding

Coarse-grained RNA folding models promise great potential for RNA structure prediction. A key component in a coarse-grained folding model is the force field. One of the challenges in the coarse-grained force field calculation is how to treat the correlation between the different degrees of freedoms. Here, we describe a new approach (IsRNA) to extract the correlated energy functions from the known structures. Through iterative molecular dynamics simulations, we build the correlation effects into the reference states, from which we extract the energy functions. The validity of IsRNA is supported by the close agreement between the simulated Boltzmann-like probability distributions for all the structure parameters and those observed from the experimentally determined structures. The correlated energy functions derived here may provide a new tool for RNA 3D structure prediction.