Effect of Water Models on The Stability of RNA: Role of Counter-Ions

Abstract

Various force fields and water model potentials influence significantly RNA conformations. The polyanionic nature of RNA attracts the water molecules and the counter ions which in turn affects their stability. The interfacial water's structural and dynamic aspects affect the RNA's base-pair opening and denaturation by breaking or making inter/intra-hydrogen bonds. Herein, we employed an MD simulations study using SPC/E and modified TIP3P water models in combination with different force fields CHARMM and AMBER to find their influence on the hydration shell of the SARS-CoV-2 RNA genome at different temperatures. AMBER-mTIP3P model was found to give more dynamic and transient conformations for RNA. The lower dielectric constant of the SPC/E model helps in the formation of the ion-contact pair near the negatively charged phosphate group (Na+-PO4−) leading to strong RNA-ion interaction and strong hydration shells having higher hydrogen bond lifetime. The Na+ ion survival probability at the interface was found to be more in the SPC/E model. At lower temperatures, the water molecules inside these hydration shells were found to be inhomogeneous, with lower void space, higher-coordinated, and non-tetrahedral. The higher dielectric constant of the mTIP3P model screened out the attraction between the ion pairs leading to a more homogenous solvation shell having a lesser hydrogen bond lifetime and more diffusive water. The distribution of the ions near the RNA structure is confirmed by metadynamics simulations. Both water models were found to disrupt the base pair orientation due to the formation of water bridges between the O2ʹ group of RNA and the water molecules.