Abstract:RNAs are highly negatively charged chain molecules. Metal ions play a crucial role in RNA folding stability and conformational changes. In this work, we employ the recently developed tightly bound ion (TBI) model, which accounts for the correlation between ions and the fluctuation of ion distributions, to investigate the ion-dependent free energy landscape for the three-way RNA junction in a 16S rRNA domain. The predicted electrostatic free energy landscape suggests that 1), ion-mediated electrostatic interactions cause an ensemble of unfolded conformations narrowly populated around the maximally extended structure; and 2), Mg2+ ion-induced correlation effects help bring the helices to the folded state. Nonelectrostatic interactions, such as noncanonical interactions within the junctions and between junctions and helix stems, might further limit the conformational diversity of the unfolded state, resulting in a more ordered unfolded state than the one predicted from the electrostatic effect. Moreover, the folded state is predominantly stabilized by the coaxial stacking force. The TBI-predicted folding stability agrees well with the experimental measurements for the different Na+ and Mg2+ ion concentrations. For Mg2+ solutions, the TBI model, which accounts for the Mg2+ ion correlation effect, gives more improved predictions than the Poisson-Boltzmann theory, which tends to underestimate the role of Mg2+ in stabilizing the folded structure. Detailed control tests indicate that the dominant ion correlation effect comes from the charge-charge Coulombic correlation rather than the size (excluded volume) correlation between the ions. Furthermore, the model gives quantitative predictions for the ion size effect in the folding energy landscape and folding cooperativity.

Types and Concentrations of Metal Ions Affect Local Structure and Dynamics of RNA.

Dynamics of Metal Ions Around an RNA Molecule

Importance of Diffuse Metal Ion Binding to RNA

RNA Captures More Cations than DNA: Insights from Molecular Dynamics Simulations

Theory and simulations for RNA folding in mixtures of monovalent and divalent cations

Metal Ion Sensing by Tetraloop-Like RNA Fragment: Role of Compact Intermediates with Non-Native Metal Ion-RNA Inner Shell Contacts

Predicting electrostatic forces in RNA folding.

Modeling Loop Composition and Ion Concentration Effects in RNA Hairpin Folding Stability

Non-specific binding of Na$^+$ and Mg$^{2+}$ to RNA determined by force spectroscopy methods

Salt-Dependent Folding Energy Landscape of RNA Three-Way Junction

Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes

Principles of RNA compaction: insights from the equilibrium folding pathway of the P4-P6 RNA domain in monovalent cations.

Predicting Ion Binding Properties for RNA Tertiary Structures

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

Predicting 3D structure, flexibility and stability of RNA hairpins in monovalent and divalent ion solutions

Ion-mediated RNA structural collapse: effect of spatial confinement

Principles of ion binding to RNA inferred from the analysis of a 1.55 Å resolution bacterial ribosome structure – Part I: Mg

Monovalent ions modulate the flux through multiple folding pathways of an RNA pseudoknot

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

Atomistic Simulations Reveal Crucial Role of Metal Ions for Ligand Binding in Gd-I Riboswitch

Ion distributions around left- and right-handed DNA and RNA duplexes: a comparative study