A coarse-grained model with implicit salt for RNAs: Predicting 3D structure, stability and salt effect

Na+-Mg2+ ion effects on conformation and translocation dynamics of single-stranded RNA: Cooperation and competition

The nanopore-based translocation of a single-stranded RNA (ssRNA) in mixed salt solution has garnered increasing interest for its biological and technological significance. However, it is challenging to comprehensively understand the effects of the mixed ion species on the translocation dynamics due to their cooperation and competition, which can be directly reflected by the ion screening and neutralizing effects, respectively. In this study, Langevin dynamics simulation is employed to investigate the properties of ssRNA conformation and translocation in mixed Na+-Mg2+ ion environments. Simulation results reveal that the ion screening effect dominates the change in the ssRNA conformational size, the ion neutralizing effect controls the capture rate of the ssRNA by the nanopore, and both of them take charge of the different changes in translocation time of the ssRNA under various mixed ion environments. Under high Na+ ion concentration, as Mg2+ concentration increases, the ion neutralizing effect strengthens, weakening the driving force inside the nanopore, leading to longer translocation time. Conversely, at low Na+ concentration, an increase in Mg2+ concentration enhances the ion screening effect, aiding in faster translocation. Furthermore, these simulation results will be explained by quantitative analysis, advancing a deeper understanding of the complicated effects of the mixed Na+-Mg2+ ions.

RNA 3D Structure Prediction: Progress and Perspective

Ribonucleic acid (RNA) molecules play vital roles in numerous important biological functions such as catalysis and gene regulation. The functions of RNAs are strongly coupled to their structures or proper structure changes, and RNA structure prediction has been paid much attention in the last two decades. Some computational models have been developed to predict RNA three-dimensional (3D) structures in silico, and these models are generally composed of predicting RNA 3D structure ensemble, evaluating near-native RNAs from the structure ensemble, and refining the identified RNAs. In this review, we will make a comprehensive overview of the recent advances in RNA 3D structure modeling, including structure ensemble prediction, evaluation, and refinement. Finally, we will emphasize some insights and perspectives in modeling RNA 3D structures.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

RNA Folding Based on 5 Beads Model and Multiscale Simulation

RNA folding prediction is very meaningful and challenging. The molecular dynamics simulation (MDS) of all atoms (AA) is limited to the folding of small RNA molecules. At present, most of the practical models are coarse grained (CG) model, and the coarse-grained force field (CGFF) parameters usually depend on known RNA structures. However, the limitation of the CGFF is obvious that it is difficult to study the modified RNA. Based on the 3 beads model (AIMS_RNA_B3), we proposed the AIMS_RNA_B5 model with three beads representing a base and two beads representing the main chain (sugar group and phosphate group). We first run the all atom molecular dynamic simulation (AAMDS), and fit the CGFF parameter with the AA trajectory. Then perform the coarse-grained molecular dynamic simulation (CGMDS). AAMDS is the foundation of CGMDS. CGMDS is mainly to carry out the conformation sampling based on the current AAMDS state and improve the folding speed. We simulated the folding of three RNAs, which belong to hairpin, pseudoknot and tRNA respectively. Compared to the AIMS_RNA_B3 model, the AIMS_RNA_B5 model is more reasonable and performs better.

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.

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.

Ribonucleic Acid Folding Prediction Based on Iterative Multiscale Simulation

RNA folding prediction is a challenge. Currently, many RNA folding models are coarse-grained (CG) with the potential derived from the known RNA structures. However, this potential is not suitable for modified and entirely new RNA. It is also not suitable for the folding simulation of RNA in the real cellular environment, including many kinds of molecular interactions. In contrast, our proposed model has the potential to address these issues, which is a multiscale simulation scheme based on all-atom (AA) force fields. We fit the CG force field using the trajectories generated by the AA force field and then iteratively perform molecular dynamics (MD) simulations of the two scales. The all-atom molecular dynamics (AAMD) simulation is mainly responsible for the correction of RNA structure, and the CGMD simulation is mainly responsible for efficient conformational sampling. On the basis of this scheme, we can successfully fold three RNAs belonging to a hairpin, a pseudoknot, and a four-way junction.

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.

Computational tools to study RNA-protein complexes

RNA is the key player in many cellular processes such as signal transduction, replication, transport, cell division, transcription, and translation. These diverse functions are accomplished through interactions of RNA with proteins. However, protein–RNA interactions are still poorly derstood in contrast to protein–protein and protein–DNA interactions. This knowledge gap can be attributed to the limited availability of protein-RNA structures along with the experimental difficulties in studying these complexes. Recent progress in computational resources has expanded the number of tools available for studying protein-RNA interactions at various molecular levels. These include tools for predicting interacting residues from primary sequences, modelling of protein-RNA complexes, predicting hotspots in these complexes and insights into derstanding in the dynamics of their interactions. Each of these tools has its strengths and limitations, which makes it significant to select an optimal approach for the question of interest. Here we present a mini review of computational tools to study different aspects of protein-RNA interactions, with focus on overall application, development of the field and the future perspectives.

FebRNA: An automated fragment-ensemble-based model for building RNA 3D structures

Knowledge of RNA three-dimensional (3D) structures is critical to understanding the important biological functions of RNAs. Although various structure prediction models have been developed, the high-accuracy predictions of RNA 3D structures are still limited to the RNAs with short lengths or with simple topology. In this work, we proposed a new model, namely FebRNA, for building RNA 3D structures through fragment assembly based on coarse-grained (CG) fragment ensembles. Specifically, FebRNA is composed of four processes: establishing the library of different types of non-redundant CG fragment ensembles regardless of the sequences, building CG 3D structure ensemble through fragment assembly, identifying top-scored CG structures through a specific CG scoring function, and rebuilding the all-atom structures from the top-scored CG ones. Extensive examination against different types of RNA structures indicates that FebRNA consistently gives the reliable predictions on RNA 3D structures, including pseudoknots, three-way junctions, four-way and five-way junctions, and RNAs in the RNA-Puzzles. FebRNA is available on the Web site: https://github.com/Tan-group/FebRNA.

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.

Diffuse Ions Coordinate Dynamics in a Ribonucleoprotein Assembly

Proper ionic concentrations are required for the functional dynamics of RNA and ribonucleoprotein (RNP) assemblies. While experimental and computational techniques have provided many insights into the properties of chelated ions, less is known about the energetic contributions of diffuse ions to large-scale conformational rearrangements. To address this, we present a model that is designed to quantify the influence of diffuse monovalent and divalent ions on the dynamics of biomolecular assemblies. This model employs all-atom (non-H) resolution and explicit ions, where effective potentials account for hydration effects. We first show that the model accurately predicts the number of excess Mg²⁺ ions for prototypical RNA systems, at a level comparable to modern coarse-grained models. We then apply the model to a complete ribosome and show how the balance between diffuse Mg²⁺ and K⁺ ions can control the dynamics of tRNA molecules during translation. The model predicts differential effects of diffuse ions on the free-energy barrier associated with tRNA entry and the energy of tRNA binding to the ribosome. Together, this analysis reveals the direct impact of diffuse ions on the dynamics of an RNP assembly.

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.

Accurate Estimation of Solvent Accessible Surface Area for Coarse-Grained Biomolecular Structures with Deep Learning

Coarse-grained (CG) models of biomolecules have been widely used in protein/ribonucleic acid (RNA) three-dimensional structure prediction, docking, drug design, and molecular simulations due to their superiority in computational efficiency. Most of these applications strongly depend on the reasonable estimation of solvation free energy, which requires the accurate calculation of solvent accessible surface area (SASA). Although algorithms for SASA calculations with all-atom protein and RNA structures have been well-established, accurately estimating the SASA based on CG structures is extremely challenging. In this work, we developed a deep learning-based SASA estimator (DeepCGSA), which can provide almost perfect SASA estimation based on CG structures of protein and RNA molecules. Extensive testing analysis showed that for three types of widely used CG protein models, including the Cα-based, Cα-Cβ, and Martini models, the correlation coefficients between the predicted values and the reference values can be as high as 0.95-0.99, which perform dramatically better than available methods. In addition, the new method can be used for CG RNA structures and unfolded protein structures with much improved accuracy. We anticipate that DeepCGSA will be highly useful in the protein/RNA structure prediction, drug design, and other applications, in which accurate estimations of SASA for CG biomolecular structures are critically important.

Predicting RNA Scaffolds with a Hybrid Method of Vfold3D and VfoldLA

The ever-increasing discoveries of noncoding RNA functions draw a strong demand for RNA structure determination from the sequence. In recently years, computational studies for RNA structures, at both the two-dimensional and the three-dimensional levels, led to several highly promising new developments. In this chapter, we describe a hybrid method, which combines the motif template-based Vfold3D model and the loop template-based VfoldLA model, to predict RNA 3D structures. The main emphasis is placed on the definition of motifs and loops, the treatment of no-template motifs, and the 3D structure assembly from templates of motifs and loops. For illustration, we use the ZIKV xrRNA1 as an example to show the template-based prediction of RNA 3D structures from the 2D structure. The web server for the hybrid model is freely accessible at http://rna.physics.missouri.edu/vfold3D2 .

In silico prediction of the in vitro behavior of polymeric gene delivery vectors

Non-viral gene delivery vectors have increasingly come under the spotlight, but their performaces are still far from being satisfactory. Therefore, there is an urgent need for forecasting tools and screening methods to enable the development of ever more effective transfectants. Here, coarse-grained (CG) models of gold standard transfectant poly(ethylene imine)s (PEIs) have been profitably used to investigate and highlight the effect of experimentally-relevant parameters, namely molecular weight (2 vs. 10 kDa) and topologies (linear vs. branched), protonation state, and ammine-to-phosphate ratios (N/Ps), on the complexation and the gene silencing efficiency of siRNA molecules. The results from the in vitro screening of cationic polymers and conditions were used to validate the in silico platform that we developed, such that the hits which came out of the CG models were of high practical relevance. We show that our in silico platform enables to foresee the most suitable conditions for the complexation of relevant siRNA-polycation assemblies, thereby providing a reliable predictive tool to test bench transfectants in silico, and foster the design and development of gene delivery vectors.

Salt-Dependent RNA Pseudoknot Stability: Effect of Spatial Confinement

Macromolecules, such as RNAs, reside in crowded cell environments, which could strongly affect the folded structures and stability of RNAs. The emergence of RNA-driven phase separation in biology further stresses the potential functional roles of molecular crowding. In this work, we employed the coarse-grained model that was previously developed by us to predict 3D structures and stability of the mouse mammary tumor virus (MMTV) pseudoknot under different spatial confinements over a wide range of salt concentrations. The results show that spatial confinements can not only enhance the compactness and stability of MMTV pseudoknot structures but also weaken the dependence of the RNA structure compactness and stability on salt concentration. Based on our microscopic analyses, we found that the effect of spatial confinement on the salt-dependent RNA pseudoknot stability mainly comes through the spatial suppression of extended conformations, which are prevalent in the partially/fully unfolded states, especially at low ion concentrations. Furthermore, our comprehensive analyses revealed that the thermally unfolding pathway of the pseudoknot can be significantly modulated by spatial confinements, since the intermediate states with more extended conformations would loss favor when spatial confinements are introduced.

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.

Salt dependent mesoscopic model for RNA at multiple strand concentrations

Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.

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.

Modeling Loop Composition and Ion Concentration Effects in RNA Hairpin Folding Stability

The ability to accurately predict RNA hairpin structure and stability for different loop sequences and salt conditions is important for understanding, modeling, and designing larger RNA folds. However, traditional RNA secondary structure models cannot treat loop-sequence and ionic effects on RNA hairpin folding. Here, we describe a general, three-dimensional (3D) conformation-based computational method for modeling salt concentration-dependent conformational distributions and the detailed 3D structures for a set of three RNA hairpins that contain a variable, 15-nucleotide loop sequence. For a given RNA sequence, the new, to our knowledge, method integrates a Vfold2D two-dimensional structure folding model with IsRNA coarse-grained molecular dynamics 3D folding simulations and Monte Carlo tightly bound ion estimations of ion-mediated electrostatic interactions. The model predicts free-energy landscapes for the different RNA hairpin-forming sequences with variable salt conditions. The theoretically predicted results agree with the experimental fluorescence measurements, validating the strategy. Furthermore, the theoretical model goes beyond the experimental results by enabling in-depth 3D structural analysis, revealing energetic mechanisms for the sequence- and salt-dependent folding stability. Although the computational framework presented here is developed for RNA hairpin systems, the general method may be applied to investigate other RNA systems, such as multiway junctions or pseudoknots in mixed metal ion solutions.

Topological constraints of RNA pseudoknotted and loop-kissing motifs: applications to three-dimensional structure prediction

An RNA global fold can be described at the level of helix orientations and relatively flexible loop conformations that connect the helices. The linkage between the helices plays an essential role in determining the structural topology, which restricts RNA local and global folds, especially for RNA tertiary structures involving cross-linked base pairs. We quantitatively analyze the topological constraints on RNA 3D conformational space, in particular, on the distribution of helix orientations, for pseudoknots and loop-loop kissing structures. The result shows that a viable conformational space is predominantly determined by the motif type, helix size, and loop size, indicating a strong topological coupling between helices and loops in RNA tertiary motifs. Moreover, the analysis indicates that (cross-linked) tertiary contacts can cause much stronger topological constraints on RNA global fold than non-cross-linked base pairs. Furthermore, based on the topological constraints encoded in the 2D structure and the 3D templates, we develop a 3D structure prediction approach. This approach can be further combined with structure probing methods to expand the capability of computational prediction for large RNA folds.

3dRNA: Building RNA 3D structure with improved template library

Most of computational methods of building RNA tertiary structure are template-based. The template-based methods usually can give more accurate 3D structures due to the use of native 3D templates, but they cannot work if the 3D templates are not available. So, a more complete library of the native 3D templates is very important for this type of methods. 3dRNA is a template-based method for building RNA tertiary structure previously proposed by us. In this paper we report improved 3D template libraries of 3dRNA by using two different schemes that give two libraries 3dRNA_Lib1 and 3dRNA_Lib2. These libraries expand the original one by nearly ten times. Benchmark shows that they can significantly increase the accuracy of 3dRNA, especially in building complex and large RNA 3D structures.

Structure folding of RNA kissing complexes in salt solutions: predicting 3D structure, stability, and folding pathway

RNA kissing complexes are essential for genomic RNA dimerization and regulation of gene expression, and their structures and stability are critical to their biological functions. In this work, we used our previously developed coarse-grained model with an implicit structure-based electrostatic potential to predict three-dimensional (3D) structures and stability of RNA kissing complexes in salt solutions. For extensive RNA kissing complexes, our model shows great reliability in predicting 3D structures from their sequences, and our additional predictions indicate that the model can capture the dependence of 3D structures of RNA kissing complexes on monovalent/divalent ion concentrations. Moreover, the comparisons with extensive experimental data show that the model can make reliable predictions on the stability for various RNA kissing complexes over wide ranges of monovalent/divalent ion concentrations. Notably, for RNA kissing complexes, our further analyses show the important contribution of coaxial stacking to the 3D structures and stronger stability than the corresponding kissing-interface duplexes at high salts. Furthermore, our comprehensive analyses for RNA kissing complexes reveal that the thermally folding pathway for a complex sequence is mainly determined by the relative stability of two possible folded states of kissing complex and extended duplex, which can be significantly modulated by its sequence.

3D structure stability of the HIV-1 TAR RNA in ion solutions: A coarse-grained model study

As an extremely common structural motif, RNA hairpins with bulge loops [e.g., the human immunodeficiency virus type 1 (HIV-1) transactivation response (TAR) RNA] can play essential roles in normal cellular processes by binding to proteins and small ligands, which could be very dependent on their three-dimensional (3D) structures and stability. Although the structures and conformational dynamics of the HIV-1 TAR RNA have been extensively studied, there are few investigations on the thermodynamic stability of the TAR RNA, especially in ion solutions, and the existing studies also have some divergence on the unfolding process of the RNA. Here, we employed our previously developed coarse-grained model with implicit salt to predict the 3D structure, stability, and unfolding pathway for the HIV-1 TAR RNA over a wide range of ion concentrations. As compared with the extensive experimental/theoretical results, the present model can give reliable predictions on the 3D structure stability of the TAR RNA from the sequence. Based on the predictions, our further comprehensive analyses on the stability of the TAR RNA as well as its variants revealed that the unfolding pathway of an RNA hairpin with a bulge loop is mainly determined by the relative stability between different states (folded state, intermediate state, and unfolded state) and the strength of the coaxial stacking between two stems in folded structures, both of which can be apparently modulated by the ion concentrations as well as the sequences.

VfoldLA: A web server for loop assembly-based prediction of putative 3D RNA structures

RNA three-dimensional (3D) structures are critical for RNA cellular functions. However, structure prediction for large and complex RNAs remains a challenge, which hampers our understanding of RNA structure-function relationship. We here report a new web server, the VfoldLA server (http://rna.physics.missouri.edu/vfoldLA), for the prediction of RNA 3D structures from nucleotide sequences and base-pair information (2D structure). This server is based on the recently developed VfoldLA, a model that classifies the single-stranded loops (junctions) into four different types and according to the loop-helix connections, assembles RNA 3D structures from the loop/junction templates. The VfoldLA web server provides a user-friendly online interface for a fully automated prediction of putative 3D RNA structures using VfoldLA. With a single-RNA or RNA-RNA complex sequence and 2D structure as input, the server generates structure(s) with the JSmol visualization along with a downloadable PDB file. The output result may serve as useful scaffolds for future structure refinement studies.

Structural Flexibility of DNA-RNA Hybrid Duplex: Stretching and Twist-Stretch Coupling

DNA-RNA hybrid (DRH) duplexes play essential roles during the replication of DNA and the reverse transcription of RNA viruses, and their flexibility is important for their biological functions. Recent experiments indicated that A-form RNA and B-form DNA have a strikingly different flexibility in stretching and twist-stretch coupling, and the structural flexibility of DRH duplex is of great interest, especially in stretching and twist-stretch coupling. In this work, we performed microsecond all-atom molecular dynamics simulations with new AMBER force fields to characterize the structural flexibility of DRH duplex in stretching and twist-stretch coupling. We have calculated all the helical parameters, stretch modulus, and twist-stretch coupling parameters for the DRH duplex. First, our analyses on structure suggest that the DRH duplex exhibits an intermediate conformation between A- and B-forms and closer to A-form, which can be attributed to the stronger rigidity of the RNA strand than the DNA strand. Second, our calculations show that the DRH duplex has the stretch modulus of 834 ± 34 pN and a very weak twist-stretch coupling. Our quantitative analyses indicate that, compared with DNA and RNA duplexes, the different flexibility of the DRH duplex in stretching and twist-stretch coupling is mainly attributed to its apparently different basepair inclination in the helical structure.

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.

Limits in accuracy and a strategy of RNA structure prediction using experimental information

RNA structural complexity and flexibility present a challenge for computational modeling efforts. Experimental information and bioinformatics data can be used as restraints to improve the accuracy of RNA tertiary structure prediction. Regarding utilization of restraints, the fundamental questions are: (i) What is the limit in prediction accuracy that one can achieve with arbitrary number of restraints? (ii) Is there a strategy for selection of the minimal number of restraints that would result in the best structural model? We address the first question by testing the limits in prediction accuracy using native contacts as restraints. To address the second question, we develop an algorithm based on the distance variation allowed by secondary structure (DVASS), which ranks restraints according to their importance to RNA tertiary structure prediction. We find that due to kinetic traps, the greatest improvement in the structure prediction accuracy is achieved when we utilize only 40-60% of the total number of native contacts as restraints. When the restraints are sorted by DVASS algorithm, using only the first 20% ranked restraints can greatly improve the prediction accuracy. Our findings suggest that only a limited number of strategically selected distance restraints can significantly assist in RNA structure modeling.

Competitive Binding of Mg2+ and Na+ Ions to Nucleic Acids: From Helices to Tertiary Structures

Nucleic acids generally reside in cellular aqueous solutions with mixed divalent/monovalent ions, and the competitive binding of divalent and monovalent ions is critical to the structures of nucleic acids because of their polyanionic nature. In this work, we first proposed a general and effective method for simulating a nucleic acid in mixed divalent/monovalent ion solutions with desired bulk ion concentrations via molecular dynamics (MD) simulations and investigated the competitive binding of Mg²⁺/Na⁺ ions to various nucleic acids by all-atom MD simulations. The extensive MD-based examinations show that single MD simulations conducted using the proposed method can yield desired bulk divalent/monovalent ion concentrations for various nucleic acids, including RNA tertiary structures. Our comprehensive analyses show that the global binding of Mg²⁺/Na⁺ to a nucleic acid is mainly dependent on its structure compactness, as well as Mg²⁺/Na⁺ concentrations, rather than the specific structure of the nucleic acid. Specifically, the relative global binding of Mg²⁺ over Na⁺ is stronger for a nucleic acid with higher effective surface charge density and higher relative Mg²⁺/Na⁺ concentrations. Furthermore, the local binding of Mg²⁺/Na⁺ to a phosphate of a nucleic acid mainly depends on the local phosphate density in addition to Mg²⁺/Na⁺ concentrations.

RNA3DCNN: Local and global quality assessments of RNA 3D structures using 3D deep convolutional neural networks

Quality assessment is essential for the computational prediction and design of RNA tertiary structures. To date, several knowledge-based statistical potentials have been proposed and proved to be effective in identifying native and near-native RNA structures. All these potentials are based on the inverse Boltzmann formula, while differing in the choice of the geometrical descriptor, reference state, and training dataset. Via an approach that diverges completely from the conventional statistical potentials, our work explored the power of a 3D convolutional neural network (CNN)-based approach as a quality evaluator for RNA 3D structures, which used a 3D grid representation of the structure as input without extracting features manually. The RNA structures were evaluated by examining each nucleotide, so our method can also provide local quality assessment. Two sets of training samples were built. The first one included 1 million samples generated by high-temperature molecular dynamics (MD) simulations and the second one included 1 million samples generated by Monte Carlo (MC) structure prediction. Both MD and MC procedures were performed for a non-redundant set of 414 RNAs. For two training datasets (one including only MD training samples and the other including both MD and MC training samples), we trained two neural networks, named RNA3DCNN_MD and RNA3DCNN_MDMC, respectively. The former is suitable for assessing near-native structures, while the latter is suitable for assessing structures covering large structural space. We tested the performance of our method and made comparisons with four other traditional scoring functions. On two of three test datasets, our method performed similarly to the state-of-the-art traditional scoring function, and on the third test dataset, our method was far superior to other scoring functions. Our method can be downloaded from https://github.com/lijunRNA/RNA3DCNN.

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.

Radial distribution function of semiflexible oligomers with stretching flexibility

The radial distribution of the end-to-end distance R_ee is crucial for quantifying the global size and flexibility of a linear polymer. For semiflexible polymers, several analytical formulas have been derived for the radial distribution of R_ee ignoring the stretching flexibility. However, for semiflexible oligomers, such as DNA or RNA, the stretching flexibility can be rather pronounced and can significantly affect the radial distribution of R_ee. In this study, we obtained an extended formula that includes the stretch modulus to describe the distribution of R_ee for semiflexible oligomers on the basis of previous formulas for semiflexible polymers without stretching flexibility. The extended formula was validated by extensive Monte Carlo simulations over wide ranges of the stretch modulus and persistence length, as well as all-atom molecular dynamics simulations of short DNAs and RNAs. Additionally, our analyses showed that the effect of stretching flexibility on the distribution of R_ee becomes negligible for DNAs longer than ∼130 base pairs and RNAs longer than ∼240 base pairs.

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.

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.

Hierarchical Assembly of RNA Three-Dimensional Structures Based on Loop Templates

The current RNA structure prediction methods cannot keep up the pace of the rapidly increasing number of sequences and the emerging new functions of RNAs. Template-based RNA three-dimensional structure prediction methods are restricted by the limited number of known RNA structures, and traditional motif-based search for the templates does not always lead to successful results. Here we report a new template search and assembly algorithm, the hierarchical loop template-assembly method (VfoldLA). The method searches for templates for single strand loop/junctions instead of the whole motifs, which often renders no available templates, or short fragments (several nucleotides), which requires a long computational time to assemble and refine. The VfoldLA method has the advantage of accounting for local and nonlocal interloop interactions. Benchmark tests indicate that this new method can provide low-resolution predictions for RNA conformations at different levels of structural complexities. Furthermore, the VfoldLA-predicted conformations may also serve as reliable putative models for further structure prediction and refinements. VfoldLA is accessible at http://rna.physics.missouri.edu/vfoldLA .

Genetic regulation mechanism of the yjdF riboswitch

The yjdF riboswitch resides in potential 5' UTRs of homologues of protein-coding gene yjdF in Firmicutes. Unlike other 30 riboswitch classes previously validated, this riboswitch class, can sense and bind to a broad collection of azaaromatic ligands. Among these compounds, some do activate production of yjdF protein driven by the riboswitch, while others are out of riboswitch-mediated modulation possibly because of the toxicity at high ligand concentrations. By incorporating the structures with pseudoknots and ligand binding kinetics into the co-transcriptional folding theory, we theoretically studied the co-transcriptional folding behaviors of the yjdF riboswitch from Bacillus subtilis at different transcription conditions. Like most riboswitches, the yjdF riboswitch can quickly fold into the aptamer structure without any trapped states during the transcription process. After the aptamer structure is formed, the riboswitch shows two main co-transcriptional folding pathways: aptamer→ON state→OFF state and aptamer → the ligand bound aptamer → the ligand bound ON state. Our results suggested that this translational riboswitch is coupled with the transcription process to exert its biological function and it is kinetically controlled. The threshold concentration for the ligand to activate the riboswitch depends on the transcription rate and the association rate of the ligand binding.

Computational Methods for Modeling Aptamers and Designing Riboswitches

Riboswitches, which are located within certain noncoding RNA region perform functions as genetic "switches", regulating when and where genes are expressed in response to certain ligands. Understanding the numerous functions of riboswitches requires computation models to predict structures and structural changes of the aptamer domains. Although aptamers often form a complex structure, computational approaches, such as RNAComposer and Rosetta, have already been applied to model the tertiary (three-dimensional (3D)) structure for several aptamers. As structural changes in aptamers must be achieved within the certain time window for effective regulation, kinetics is another key point for understanding aptamer function in riboswitch-mediated gene regulation. The coarse-grained self-organized polymer (SOP) model using Langevin dynamics simulation has been successfully developed to investigate folding kinetics of aptamers, while their co-transcriptional folding kinetics can be modeled by the helix-based computational method and BarMap approach. Based on the known aptamers, the web server Riboswitch Calculator and other theoretical methods provide a new tool to design synthetic riboswitches. This review will represent an overview of these computational methods for modeling structure and kinetics of riboswitch aptamers and for designing riboswitches.

Optimization of RNA 3D structure prediction using evolutionary restraints of nucleotide-nucleotide interactions from direct coupling analysis

Direct coupling analysis of nucleotide coevolution provides a novel approach to identify which nucleotides in an RNA molecule are likely in direct contact, and this information obtained from sequence only can be used to predict RNA 3D structures with much improved accuracy. Here we present an efficient method that incorporates this information into current RNA 3D structure prediction methods, specifically 3dRNA. Our method makes much more accurate RNA 3D structure prediction than the original 3dRNA as well as other existing prediction methods that used the direct coupling analysis. In particular our method demonstrates a significant improvement in predicting multi-branch junction conformations, a major bottleneck for RNA 3D structure prediction. We also show that our method can be used to optimize the predictions by other methods. These results indicate that optimization of RNA 3D structure prediction using evolutionary restraints of nucleotide-nucleotide interactions from direct coupling analysis offers an efficient way for accurate RNA tertiary structure predictions.

Understanding the Relative Flexibility of RNA and DNA Duplexes: Stretching and Twist-Stretch Coupling

The flexibility of double-stranded (ds) RNA and dsDNA is crucial for their biological functions. Recent experiments have shown that the flexibility of dsRNA and dsDNA can be distinctively different in the aspects of stretching and twist-stretch coupling. Although various studies have been performed to understand the flexibility of dsRNA and dsDNA, there is still a lack of deep understanding of the distinctive differences in the flexibility of dsRNA and dsDNA helices as pertains to their stretching and twist-stretch coupling. In this work, we have explored the relative flexibility in stretching and twist-stretch coupling between dsRNA and dsDNA by all-atom molecular dynamics simulations. The calculated stretch modulus and twist-stretch coupling are in good accordance with the existing experiments. Our analyses show that the differences in stretching and twist-stretch coupling between dsRNA and dsDNA helices are mainly attributed to their different (A- and B-form) helical structures. Stronger basepair inclination and slide in dsRNA is responsible for the apparently weaker stretching rigidity versus that of dsDNA, and the opposite twist-stretch coupling for dsRNA and dsDNA is also attributed to the stronger basepair inclination in dsRNA than in dsDNA. Our calculated macroscopic elastic parameters and microscopic analyses are tested and validated by different force fields for both dsRNA and dsDNA.

Martini Coarse-Grained Force Field: Extension to RNA

RNA has an important role not only as the messenger of genetic information but also as a regulator of gene expression. Given its central role in cell biology, there is significant interest in studying the structural and dynamic behavior of RNA in relation to other biomolecules. Coarse-grain molecular dynamics simulations are a key tool to that end. Here, we have extended the coarse-grain Martini force field to include RNA after our recent extension to DNA. In the same way DNA was modeled, the tertiary structure of RNA is constrained using an elastic network. This model, therefore, is not designed for applications involving RNA folding but rather offers a stable RNA structure for studying RNA interactions with other (bio)molecules. The RNA model is compatible with all other Martini models and opens the way to large-scale explicit-solvent molecular dynamics simulations of complex systems involving RNA.

Theory and Modeling of RNA Structure and Interactions with Metal Ions and Small Molecules

In addition to continuous rapid progress in RNA structure determination, probing, and biophysical studies, the past decade has seen remarkable advances in the development of a new generation of RNA folding theories and models. In this article, we review RNA structure prediction models and models for ion-RNA and ligand-RNA interactions. These new models are becoming increasingly important for a mechanistic understanding of RNA function and quantitative design of RNA nanotechnology. We focus on new methods for physics-based, knowledge-based, and experimental data-directed modeling for RNA structures and explore the new theories for the predictions of metal ion and ligand binding sites and metal ion-dependent RNA stabilities. The integration of these new methods with theories about the cellular environment effects in RNA folding, such as molecular crowding and cotranscriptional kinetic effects, may ultimately lead to an all-encompassing RNA folding model.

Using 3dRNA for RNA 3-D Structure Prediction and Evaluation

This unit describes how to use 3dRNA to predict RNA 3-D structures from their sequences and secondary (2-D) structures, and how to use 3dRNAscore to evaluate the predicted structures. The predicted RNA 3-D structures can be used to predict or understand their functions and can also be used to find the interactions between the RNA and other molecules. © 2017 by John Wiley & Sons, Inc.

Reversible-Switch Mechanism of the SAM-III Riboswitch

Riboswitches are self-regulatory elements located at the 5' untranslated region of certain mRNAs. The Enterococcus faecalis SAM-III (S_MK) riboswitch regulates downstream gene expression through conformational change by sensing S-adenosylmethionine (SAM) at the translation level. Using the recently developed systematic helix-based computational method, we studied the co-transcriptional folding behavior of the S_MK riboswitch and its shortened construct lacking the first six nucleotides. We find that there are no obvious misfolded structures formed during the transcription and refolding processes for this riboswitch. The full-length riboswitch quickly folds into the ON-state in the absence of SAM, and the coupling between transcription and translation is not required for the riboswitch to function. The potential to form helix P0 is necessary for the riboswitch to function as a switch. For this thermodynamically controlled reversible riboswitch, the fast helix-exchanging transition pathway between the two functional structures guaranteed that this riboswitch can act as a reversible riboswitch.

Predicting 3D Structure, Flexibility, and Stability of RNA Hairpins in Monovalent and Divalent Ion Solutions

A full understanding of RNA-mediated biology would require the knowledge of three-dimensional (3D) structures, structural flexibility, and stability of RNAs. To predict RNA 3D structures and stability, we have previously proposed a three-bead coarse-grained predictive model with implicit salt/solvent potentials. In this study, we further develop the model by improving the implicit-salt electrostatic potential and including a sequence-dependent coaxial stacking potential to enable the model to simulate RNA 3D structure folding in divalent/monovalent ion solutions. The model presented here can predict 3D structures of RNA hairpins with bulges/internal loops (<77 nucleotides) from their sequences at the corresponding experimental ion conditions with an overall improved accuracy compared to the experimental data; the model also makes reliable predictions for the flexibility of RNA hairpins with bulge loops of different lengths at several divalent/monovalent ion conditions. In addition, the model successfully predicts the stability of RNA hairpins with various loops/stems in divalent/monovalent ion solutions.

Structure Prediction of RNA Loops with a Probabilistic Approach

The knowledge of the tertiary structure of RNA loops is important for understanding their functions. In this work we develop an efficient approach named RNApps, specifically designed for predicting the tertiary structure of RNA loops, including hairpin loops, internal loops, and multi-way junction loops. It includes a probabilistic coarse-grained RNA model, an all-atom statistical energy function, a sequential Monte Carlo growth algorithm, and a simulated annealing procedure. The approach is tested with a dataset including nine RNA loops, a 23S ribosomal RNA, and a large dataset containing 876 RNAs. The performance is evaluated and compared with a homology modeling based predictor and an ab initio predictor. It is found that RNApps has comparable performance with the former one and outdoes the latter in terms of structure predictions. The approach holds great promise for accurate and efficient RNA tertiary structure prediction.