Abstract:The computer-aided folding of biomolecules, particularly RNAs, is one of the most difficult challenges in computational structural biology. RNA tetraloops are fundamental RNA motifs playing key roles in RNA folding and RNA-RNA and RNA-protein interactions. Although state-of-the-art Molecular Dynamics (MD) force fields correctly describe the native state of these tetraloops as a stable free-energy basin on the microsecond time scale, enhanced sampling techniques reveal that the native state is not the global free energy minimum, suggesting yet unidentified significant imbalances in the force fields. Here, we tested our ability to fold the RNA tetraloops in various force fields and simulation settings. We employed three different enhanced sampling techniques, namely, temperature replica exchange MD (T-REMD), replica exchange with solute tempering (REST2), and well-tempered metadynamics (WT-MetaD). We aimed to separate problems caused by limited sampling from those due to force-field inaccuracies. We found that none of the contemporary force fields is able to correctly describe folding of the 5'-GAGA-3' tetraloop over a range of simulation conditions. We thus aimed to identify which terms of the force field are responsible for this poor description of TL folding. We showed that at least two different imbalances contribute to this behavior, namely, overstabilization of base-phosphate and/or sugar-phosphate interactions and underestimated stability of the hydrogen bonding interaction in base pairing. The first artifact stabilizes the unfolded ensemble, while the second one destabilizes the folded state. The former problem might be partially alleviated by reparametrization of the van der Waals parameters of the phosphate oxygens suggested by Case et al., while in order to overcome the latter effect we suggest local potentials to better capture hydrogen bonding interactions.

Parametric analysis of RNA folding

Geometric combinatorics and computational molecular biology: branching polytopes for RNA sequences

RNA fold prediction by Monte Carlo in graph space and the statistical mechanics of tertiary interactions

Polynomial-time parametric optimisation

RNA Folding on the 3D Triangular Lattice

RNA folding algorithm based on fuzzy dynamic programming

RNA secondary structure formation: a solvable model of heteropolymer folding

LinearFold: linear-time approximate RNA folding by 5'-to-3' dynamic programming and beam search

The four ingredients of single-sequence RNA secondary structure prediction. A unifying perspective

Helix-Based RNA Landscape Partition and Alternative Secondary Structure Determination

Similarity analysis of RNA secondary structures based on 4D representation

Ab initio RNA folding

Statistical Physics of RNA-folding

RNA folding using quantum computers

Detecting Ribosnitches with RNA Folding Algorithms: a Genome-Wide Benchmark

Novel Representation of RNA Secondary Structure Used to Improve Prediction Algorithms.

RNA Secondary Structure Thermodynamics

Bayesian Analysis of Folding and Unfolding Time Series of Single-Forced RNAs

Monte Carlo Simulation for Single RNA Unfolding by Force.

Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies

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