Abstract:Molecular dynamics (MD) simulations that rely on force field methods has been widely used to explore the structure and function of RNAs. However, the current commonly used force fields are limited by the electrostatic description offered by atomic charge, dipole and at most quadrupole moments, failing to capture the anisotropic picture of electronic features. Actually, the distribution of electrons around atomic nuclei is not spherically symmetric but is geometry dependent. A multipolar electrostatic model based on high rank multipole moments is described in this work, which allows us to combine polarizability and anisotropy of electron density. RNA secondary structure was taken as a research system, and its substructures including stem, loops (hairpin loop, bulge loop, internal loop, and multi-branch loop), and pseudoknots (H-type and K-type) were investigated, respectively, as well as the hairpin. First, the atom-atom electrostatic properties derived from one chain of a duplex RNA 2MVY in our previous work (Ref. 58) were measured by the pilot RNA systems of hairpin, hairpin loop, stem, and H-type pseudoknot, respectively. The prediction results were not satisfactory. Consequently, to obtain a general set of electrostatic parameters for RNA force fields, the convergence behavior of the atom-atom electrostatic interactions in the pilot RNA systems was explored using high rank atomic multipole moments. The pilot RNA systems were cut into four types of different-sized molecular fragments, and the single nucleotide fragment and nucleotide-paired fragment proved to be the most reasonable systems for base-unpairing regions and base-pairing regions to investigate the convergence behavior of all types of atom-atom electrostatic interactions, respectively. Transferability of the electrostatic properties drawn from the pilot RNA systems to the corresponding test systems was also investigated. Furthermore, the convergence behavior of atomic electrostatic interactions in other substructures including bulge loop, internal loop, multi-branch loop, and K-type pseudoknot was expected to be modeled via the hairpin.

RNA force field with accuracy comparable to state-of-the-art protein force fields

Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields

Base-specific RNA Force Field Improving the Dynamics Conformation of Nucleotide.

Molecular Dynamics Simulations of Protein RNA Complexes by Using an Advanced Electrostatic Model

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA–Protein Complexes

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

RNA-Specific Force Field Optimization with CMAP and Reweighting

Benchmarking RNA all-atom force fields using hairpin motifs

Polarizable force field for RNA based on the classical drude oscillator

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

Assessing RNA atomistic force fields via energy landscape explorations in implicit solvent

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Simulations of the molecular dynamics of nucleic acids

RNA-protein complexes and force field polarizability

Multipolar Electrostatics for Hairpin and Pseudoknots in RNA: Improving the Accuracy of Force Field Potential Energy Function

Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics

Assessing Force Field Potential Energy Function Accuracy Via a Multipolar Description of Atomic Electrostatic Interactions in RNA.

Adjusting the energy profile for CH–O interactions leads to improved stability of RNA stem-loop structures in MD simulations

The development of nucleic acids force fields: From an unchallenged past to a competitive future

All-Atom Knowledge-Based Potential for RNA Structure Discrimination Based on the Distance-Scaled Finite Ideal-Gas Reference State.

RNA stability under different combinations of amber force fields and solvation models.