El Khoury Léa, Naseem-Khan Sehr, Kwapien Karolina, Hobaika Zeina, Maroun Richard G, Piquemal Jean-Philip, Gresh Nohad
Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC, UMR7616 CNRS, Paris, France.
Centre d'Analyses et de Recherche, UR EGFEM, LSIM, Faculté des Sciences, Saint Joseph University of Beirut, BP 11-514, Riad El Solh, Beirut, 1116-2050, Lebanon.
J Comput Chem. 2017 Aug 15;38(22):1897-1920. doi: 10.1002/jcc.24830. Epub 2017 May 30.
A correct representation of the short-range contributions such as exchange-repulsion (E ) and charge-transfer (E ) is essential for the soundness of separable, anisotropic polarizable molecular mechanics potentials. Within the context of the SIBFA procedure, this is aimed at by explicit representations of lone pairs in their expressions. It is necessary to account for their anisotropic behaviors upon performing not only in-plane, but also out-of-plane, variations of a probe molecule or cation interacting with a target molecule or molecular fragment. Thus, E and E have to reproduce satisfactorily the corresponding anisotropies of their quantum chemical (QC) counterparts. A significant improvement of the out-of-plane dependencies was enabled when the sp and sp localized lone-pairs are, even though to a limited extent, delocalized on both sides of the plane, above and below the atom bearer but at the closely similar angles as the in-plane lone pair. We report calibration and validation tests on a series of monoligated complexes of a probe Zn(II) cation with several biochemically relevant ligands. Validations are then performed on several polyligated Zn(II) complexes found in the recognition sites of Zn-metalloproteins. Such calibrations and validations are extended to representative monoligated and polyligated complexes of Mg(II) and Ca(II). It is emphasized that the calibration of all three cations was for each ΔE contribution done on a small training set bearing on a limited number of representative N, O, and S monoligated complexes. Owing to the separable nature of ΔE, a secure transferability is enabled to a diversity of polyligated complexes. For these the relative errors with respect to the target ΔE(QC) values are generally < 3%. Overall, the article proposes a full set of benchmarks that could be useful for force field developers. © 2017 Wiley Periodicals, Inc.
对于可分离的各向异性极化分子力学势的合理性而言,正确表示诸如交换排斥(E )和电荷转移(E )等短程贡献至关重要。在SIBFA程序的背景下,这通过在其表达式中明确表示孤对电子来实现。不仅要考虑与目标分子或分子片段相互作用的探针分子或阳离子在平面内变化时孤对电子的各向异性行为,还需考虑其在平面外变化时的情况。因此,E 和E 必须令人满意地再现其量子化学(QC)对应物的相应各向异性。当sp和sp局域孤对电子即使在有限程度上在平面两侧、原子载体上方和下方离域,但与平面内孤对电子的角度非常相似时,平面外依赖性有了显著改善。我们报告了对一系列探针Zn(II)阳离子与几种生物化学相关配体形成的单配位络合物的校准和验证测试。然后对在锌金属蛋白识别位点发现的几种多配位Zn(II)络合物进行验证。此类校准和验证扩展到Mg(II)和Ca(II)的代表性单配位和多配位络合物。需要强调的是,对所有三种阳离子的校准都是针对每个ΔE贡献在一个包含有限数量代表性N、O和S单配位络合物的小训练集上进行的。由于ΔE的可分离性质,能够安全地转移到各种多配位络合物。对于这些络合物,相对于目标ΔE(QC)值的相对误差通常<3%。总体而言,本文提出了一整套对力场开发者可能有用的基准。© 2017威利期刊公司