Differences in Ion-RNA binding modes due to charge density variations explain the stability of RNA in monovalent salts

The stability of RNA increases as the charge density of the alkali metal cations increases. The molecular mechanism for the stability does not exist. To fill this gap, we performed all-atom molecular dynamics (MD) pulling simulations to dissect the microscopic origin of this phenomenon. We first established that the free energy landscape obtained in the simulations is in excellent agreement with the single-molecule optical tweezer experiments. The origin of the stronger stability in Na+ compared to K+ is found to be due to the differences in the charge-density related binding modes. The smaller hydrated Na+ ion preferentially binds to the highly charged phosphates. In contrast, the larger K+ ions interact with the major grooves. As a result, the electrostatic repulsion between the phosphate groups is reduced more effectively by Na+ ions. Because the proposed mechanism is generic, we predict that the same conclusions are valid for divalent alkaline earth metal cations.