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.

CEITEC Central European Institute of Technology Masaryk University Campus Bohunice Kamenice 5 625 00 Brno Czech Republic

Department of Biomedical Engineering The University of Texas at Austin Austin Texas 78712 United States

Institut Universitaire de France Paris Cedex 05 75231 France

Institute of Biophysics Academy of Sciences of the Czech Republic Kralovpolska 135 612 65 Brno Czech Republic

Laboratoire de Chimie Théorique Sorbonne Universités UPMC UMR7616 CNRS 75006Paris France

Citace poskytuje Crossref.org

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview