title: Ion Binding to Quadruplex DNA Stems. Comparison of MM and QM Descriptions Reveals Sizable Polarization Effects Not Included in Contemporary Simulations.
abstract: Molecular mechanical (MM) force fields are commonly employed for biomolecular simulations. Despite their success, the nonpolarizable nature of contemporary additive force fields limits their performance, especially in long simulations and when strong polarization effects are present. Guanine quadruplex D(R)NA molecules have been successfully studied by MM simulations in the past. However, the G-stems are stabilized by a chain of monovalent cations that create sizable polarization effects. Indeed, simulation studies revealed several problems that have been tentatively attributed to the lack of polarization. Here, we provide a detailed comparison between quantum chemical (QM) DFT-D3 and MM potential energy surfaces of ion binding to G-stems and assess differences that may affect MM simulations. We suggest that MM describes binding of a single ion to the G-stem rather well. However, polarization effects become very significant when a second ion is present. We suggest that the MM approximation substantially limits accuracy of description of energy and dynamics of multiple ions inside the G-stems and binding of ions at the stem-loop junctions. The difference between QM and MM descriptions is also explored using symmetry-adapted perturbation theory and quantum theory of atoms in molecules analyses, which reveal a delicate balance of electrostatic and induction effects.
authors.affiliation: CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic.
authors.forename: Konstantinos
authors.initials: K
authors.lastname: Gkionis
authors.affiliation: CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic.
authors.forename: Holger
authors.initials: H
authors.lastname: Kruse
authors.affiliation: School of Chemistry, Cardiff University , Park Place, Cardiff CF10 3AT, United Kingdom.
authors.forename: James A
authors.initials: JA
authors.lastname: Platts
authors.affiliation: Institute of Biophysics , Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.
authors.forename: Arnošt
authors.initials: A
authors.lastname: Mládek
authors.affiliation: CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic.
authors.forename: Jaroslav
authors.initials: J
authors.lastname: Koča
authors.affiliation: CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic.|Institute of Biophysics , Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.
authors.forename: Jiří
authors.initials: J
authors.lastname: Šponer