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.

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

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

Multipolar Description of Atom–Atom Electrostatic Interaction Energies in Single/Double-Stranded DNAs

Validation of Polarizable Force Field Parameters for Nucleic Acids by Inter-Molecular Interactions

Polarizable force field for RNA based on the classical drude oscillator

A systematic study on RNA NMR chemical shift calculation based on the automated fragmentation QM/MM approach

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

Quantum Chemical Calculations with Machine Learning for Multipolar Electrostatics Prediction in RNA: an Application to Pentose.

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

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

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Extensive Evaluation of Force Fields for G-Quadruplexes

Combined QM/MM, Machine Learning Path Integral Approach to Compute Free Energy Profiles and Kinetic Isotope Effects in RNA Cleavage Reactions.

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

Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

RNA-protein complexes and force field polarizability

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

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

RNA-Specific Force Field Optimization with CMAP and Reweighting

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