Importance of explicit smeared lone-pairs in anisotropic polarizable molecular mechanics. Torture track angular tests for exchange-repulsion and charge transfer contributions

5 years ago

Importance of explicit smeared lone-pairs in anisotropic polarizable molecular mechanics. Torture track angular tests for exchange-repulsion and charge transfer contributions

Jean-Philip Piquemal, Léa El Khoury, Richard G. Maroun, Sehr Naseem-Khan, Zeina Hobaika, Karolina Kwapien, Nohad Gresh

A correct representation of the short-range contributions such as exchange-repulsion (Erep) and charge-transfer (Ect) 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, Erep and Ect 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 sp2 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. We report a further refinement in the representation of the conjugated sp and sp2 lone-pairs (lp). These are, in the context of the polarizable SIBFA procedure, involved in the expression of the short-range contributions Erep and Ect. This is done by partially ‘smearing’ such lp's above and below the atom bearer. A significant improvement of the out-of-plane angular dependencies of Erep and Ect results. Further validations are reported on a diverse series of polyligated complexes of Zn(II), Mg(II), and Ca(II).

Publisher URL: http://onlinelibrary.wiley.com/resolve/doi

DOI: 10.1002/jcc.24830