期刊
JOURNAL OF CHEMICAL THEORY AND COMPUTATION
卷 10, 期 3, 页码 1326-1340出版社
AMER CHEMICAL SOC
DOI: 10.1021/ct4009969
关键词
-
资金
- Czech Science Foundation [P208/11/1822]
- European Regional Development Fund [CZ.1.05/1.1.00/02.0068]
- EU Seventh Framework Programme under the Capacities specific programme [286154]
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.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据