z-logo
open-access-imgOpen Access
A Comparative Study of Transferable Aspherical Pseudoatom Databank and Classical Force Fields for Predicting Electrostatic Interactions in Molecular Dimers
Author(s) -
Prashant Kumar,
Sławomir A. Bojarowski,
Katarzy. Jarzembska,
Sławomir Domagała,
Kenno Vanommeslaeghe,
Alexander D. MacKerell,
Paulina M. Dominiak
Publication year - 2014
Publication title -
journal of chemical theory and computation
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.001
H-Index - 185
eISSN - 1549-9626
pISSN - 1549-9618
DOI - 10.1021/ct4011129
Subject(s) - van der waals force , electrostatics , force field (fiction) , intermolecular force , point particle , perturbation theory (quantum mechanics) , macromolecule , range (aeronautics) , electric potential energy , field (mathematics) , wave function , computational chemistry , chemistry , quantum , statistical physics , chemical physics , physics , quantum mechanics , molecule , materials science , energy (signal processing) , mathematics , biochemistry , pure mathematics , composite material
Accurate and fast evaluation of electrostatic interactions in molecular systems is one of the most challenging tasks in the rapidly advancing field of macromolecular chemistry and drug design. Electrostatic interactions are of crucial importance in biological systems. They are well represented by quantum mechanical methods; however, such calculations are computationally expensive. In this study, we have evaluated the University of Buffalo Pseudoatom Databank (UBDB)1,2 approach for approximation of electrostatic properties of macromolecules and their complexes. We selected the S663 and JSCH-20054 data sets (208 molecular complexes in total) for this study. These complexes represent a wide range of chemical and biological systems for which hydrogen bonding, electrostatic, and van der Waals interactions play important roles. Reference electrostatic energies were obtained directly from wave functions at the B3LYP/aug-cc-pVTZ level of theory using the SAPT (Symmetry-Adapted Perturbation Theory) scheme for calculation of electrostatic contributions to total intermolecular interaction energies. Electrostatic energies calculated on the basis of the UBDB were compared with corresponding reference results. Results were also compared with energies computed using a point charge model from popular force fields (AM1-BCC and RESP used in AMBER and CGenFF from CHARMM family). The energy trends are quite consistent ( R 2 ≈ 0.98) for the UBDB method as compared to the AMBER5 and CHARMM force field methods6( R 2 ≈ 0.93 on average). The RSMEs do not exceed 3.2 kcal mol -1 for the UBDB and are in the range of 3.7-7.6 kcal mol -1 for the point charge models. We also investigated the discrepancies in electrostatic potentials and magnitudes of dipole moments among the tested methods. This study shows that estimation of electrostatic interaction energies using the UBDB databank is accurate and reasonably fast when compared to other known methods, which opens potential new applications to macromolecules.

The content you want is available to Zendy users.

Already have an account? Click here to sign in.
Having issues? You can contact us here