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.

Force-induced denaturation of RNA

Mechanically probing the folding pathway of single RNA molecules

Force unfolding kinetics of RNA using optical tweezers. I. Effects of experimental variables on measured results

Force unfolding kinetics of RNA using optical tweezers. II. Modeling experiments

Unraveling RNA by Mechanical Unzipping

RNA under Tension: Folding Landscapes, Kinetic Partitioning Mechanism, and Molecular Tensegrity

Force-dependent elasticity of nucleic acids

Atomistic simulations of RNA duplex thermal denaturation: Sequence- and forcefield-dependence

RNA unwinding from reweighted pulling simulations

Non-specific binding of Na$^+$ and Mg$^{2+}$ to RNA determined by force spectroscopy methods

Statistical Physics of RNA-folding

Stem-loop formation drives RNA folding in mechanical unzipping experiments

Force Unfolding Single RNAs: from Equilibrium to Far-from Equilibrium

Can We Ever Develop an Ideal RNA Force Field? Lessons Learned from Simulations of UUCG RNA Tetraloop and Other Systems.

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather Than Electrostatic Repulsion.

Remarkable similarity of force induced dsRNA conformational changes to stretched dsDNA and their detection using electrical measurements

Force Dependence of Proteins' Transition State Position and the Bell-Evans Model

DNA nano-mechanics: how proteins deform the double helix

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

Theory for RNA folding, stretching, and melting including loops and salt

Theoretical models for single-molecule DNA and RNA experiments: from elasticity to unzipping