Unravelling Mg${}^{2+}$-RNA binding with atomistic molecular dynamics

Abstract

Interaction with divalent cations is of paramount importance for RNA structural stability and function. We here report a detailed molecular dynamics study of all the possible binding sites for Mg${}^{2+}$ on a RNA duplex, including both direct (inner sphere) and indirect (outer sphere) binding. In order to tackle sampling issues, we develop a modified version of bias-exchange metadynamics which allows us to simultaneously compute affinities with previously unreported statistical accuracy. Results correctly reproduce trends observed in crystallographic databases. Based on this, we simulate a carefully chosen set of models that allows us to quantify the effects of competition with monovalent cations, RNA flexibility, and RNA hybridization. Our simulations reproduce the decrease and increase of Mg${}^{2+}$ affinity due to ion competition and hybridization respectively, and predict that RNA flexibility has a site dependent effect. This suggests a non trivial interplay between RNA conformational entropy and divalent cation binding.

References in corpus (7)

• Canonical sampling through velocity-rescaling
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• Metadynamics with adaptive Gaussians
• Molecular Dynamics Simulations of Solutions at Constant Chemical Potential
• Empirical corrections to the Amber RNA force field with Target Metadynamics
• RNA Folding Pathways in Stop Motion
• Kissing loop interaction in adenine riboswitch: insights from umbrella sampling simulations

Cited by in corpus (8)

• Using the Maximum Entropy Principle to Combine Simulations and Solution Experiments
• Exploring RNA structure and dynamics through enhanced sampling simulations
• Analyzing and biasing simulations with PLUMED
• The mechanism of RNA base fraying: molecular dynamics simulations analyzed with core-set Markov state models
• Refinement of molecular dynamics ensembles using experimental data and flexible forward models
• Effects and limitations of a nucleobase-driven backmapping procedure for nucleic acids using steered Molecular Dynamics
• Exponential Scaling of Water Exchange Rates with Ion Interaction Strength from the Perspective of Dynamic Facilitation Theory
• Molecular dynamics simulations with grand-canonical reweighting suggest cooperativity effects in RNA structure probing experiments