Channeling through Two Stacked Guanine Quartets of One and Two Alkali Cations in the Li+, Na+, K+, and Rb+ Series. Assessment of the Accuracy of the SIBFA Anisotropic Polarizable Molecular Mechanics Potential
Language English Country United States Media print-electronic
Document type Journal Article, Research Support, Non-U.S. Gov't
- MeSH
- Alkalies chemistry MeSH
- Anisotropy MeSH
- Potassium chemistry MeSH
- G-Quadruplexes * MeSH
- Guanine chemistry MeSH
- Cations, Monovalent chemistry MeSH
- Lithium chemistry MeSH
- Models, Molecular MeSH
- Rubidium chemistry MeSH
- Sodium chemistry MeSH
- Thermodynamics MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
- Names of Substances
- Alkalies MeSH
- Potassium MeSH
- Guanine MeSH
- Cations, Monovalent MeSH
- Lithium MeSH
- Rubidium MeSH
- Sodium MeSH
Stacking of guanine quartets (GQs) can trigger the formation of DNA or RNA quadruple helices, which play numerous biochemical roles. The GQs are stabilized by alkali cations, mainly K+ and Na+, which can reside in, or channel through, the central axis of the GQ stems. Further, ion conduction through GQ wires can be leveraged for nanochemistry applications. G-quadruplex systems have been extensively studied by classical molecular dynamics (MD) simulations using pair-additive force fields or by quantum-chemical (QC) calculations. However, the non-polarizable force fields are very approximate, while QC calculations lack the necessary sampling. Thus, ultimate description of GQ systems would require long-enough simulations using advanced polarizable molecular mechanics (MM). However, to perform such calculations, it is first mandatory to evaluate the method's accuracy using benchmark QC. We report such an evaluation for SIBFA polarizable MM, bearing on the channeling (movement) of an alkali cation (Li+, Na+, K+, or Rb+) along the axis of two stacked G quartets interacting with either one or two ions. The QC energy profiles display markedly different features depending upon the cation but can be retrieved in the majority of cases by the SIBFA profiles. An appropriate balance of first-order (electrostatic and short-range repulsion) and second-order (polarization, charge-transfer, and dispersion) contributions within ΔE is mandatory. With two cations in the channel, the relative weights of the second-order contributions increase steadily upon increasing the ion size. In the G8 complexes with two K+ or two Rb+ cations, the sum of polarization and charge-transfer exceeds the first order terms for all ion positions.
Institut Universitaire de France Paris Cedex 05 75231 France
Laboratoire de Chimie Théorique Sorbonne Universités UPMC UMR7616 CNRS 75006Paris France
References provided by Crossref.org
RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview
