Zgarbová Marie, Jurečka Petr, Šponer Jiří, Otyepka Michal
Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic.
J Chem Theory Comput. 2018 Jan 9;14(1):319-328. doi: 10.1021/acs.jctc.7b00926. Epub 2017 Dec 21.
The A/B transition is a basic element of DNA conformational change. Because of its involvement in the sensing of the ionic conditions by DNA and in specific protein-DNA interactions, this transition is important for biological functions of DNA. Therefore, accurate modeling of the A/B equilibrium by means of empirical force fields is of utmost interest. In this work, we examine the A/B equilibrium in three AMBER force fields, including the recent bsc1 and OL15 modifications, using much longer MD simulations than attempted before. Special attention is paid to the coupling of the A/B equilibrium with the south/north (S/N) transition of the sugar pucker. We found that none of the tested force fields provided a satisfactory description of the A/B equilibrium because the B-form was predicted to be much too stable and the A-form was predicted to be almost absent even in concentrated trifluoroethanol solutions. Based on comparison with NMR data for duplexes and single nucleosides, we hypothesize that this problem arose from the incorrect description of the S/N equilibrium of sugar pucker, where the south conformation is much too stable, thus stabilizing the B-form. Because neither the A/B equilibrium in duplexes nor the S/N equilibrium in nucleosides was described accurately, further refinements of the AMBER DNA force fields are needed.
A/B转变是DNA构象变化的一个基本要素。由于其参与了DNA对离子条件的感知以及特定的蛋白质-DNA相互作用,这种转变对DNA的生物学功能很重要。因此,通过经验力场对A/B平衡进行精确建模极具意义。在这项工作中,我们使用了比以往尝试的长得多的分子动力学(MD)模拟,研究了三种AMBER力场中的A/B平衡,包括最近的bsc1和OL15修正。特别关注了A/B平衡与糖环构象的南/北(S/N)转变之间的耦合。我们发现,所测试的力场均未对A/B平衡给出令人满意的描述,因为预测B型过于稳定,而即使在浓三氟乙醇溶液中,A型也几乎不存在。基于与双链体和单核苷的核磁共振(NMR)数据的比较,我们推测这个问题源于对糖环构象S/N平衡的错误描述,其中南构象过于稳定,从而稳定了B型。由于双链体中的A/B平衡和核苷中的S/N平衡均未得到准确描述,因此需要对AMBER DNA力场进行进一步优化。