20th-banner-rev.png

Albany 2019: 20th Conversation - Abstracts

category image Albany 2019
Conversation 20
June 11-15 2019
Adenine Press (2019)

A Revisiting of the RESP Charge Derivation Model

The representation of the electrostatic interactions by Coulombic interactions between the atom-centered partial charges is a fundamental part of the molecular mechanics and empirical force field methods. The broad success of the AMBER force field family originates mainly in the RESP charge model (Bayly et al., 1993) that derives the partial charges to reproduce the electrostatic field around the molecules.

In this study we revisited the RESP charge derivation model in order to improve its description of the electrostatic potential around the molecules and thus the description of the electrostatic interactions in the force field. In particular, we re-optimized the atomic radii used for definition of the grid points for evaluation of the electrostatic field around the molecule. Some grid points introduced by the standard RESP procedure using Singh-Merz-Kollman radii (Singh et al., 1984) are deeply buried into the electron density, especially in case of aromatic molecules such as nucleobases, so that they significantly bias the fitted charges in an artificial way. In addition, we redefined the restraining scheme, in which we replaced the hyperbolic restraints toward zero charge values introduced by Bayly et al. by weighted parabolic restraints toward CM5 charges.(Marenich et al., 2012) These restraints are able to selectively eliminate the poor statistical of the ESP charges of the atoms buried in the molecule. On the other hand, they are negligibly small whenever the ESP charges are statistically well-defined. The redefined charges are able to reproduce the electrostatic potential around the molecules more accurately than standard RESP charges, especially for aromatic systems, and are able, e.g., to entirely eliminate underestimation of the base pairing interactions between nucleobases (Banas et al., 2012), which precludes accurate description of the nucleic acids by empirical force fields.

The authors gratefully acknowledge the support by the Ministry of Education, Youth and Sports of the Czech Republic and Czech Science Foundation no. 18-25349S and by the Operational Programme Research, Development and Education – European Regional Development Fund project no. CZ.02.1.01/0.0/0.0/16_019/0000754.

References

Bayly, C. I.; Cieplak, P.; Cornell, W. D.; Kollman, P. A. J. Phys. Chem. 1993, 97, 10269-10280.

Singh, U. C.; Kollman, P. A. J. Comput. Chem. 1984, 5, 129-145.

Pavel Banáš1,2
Michal Janeček1
Petra Kührová1
Michal Otyepka1,2
Jiří Šponer 1,2

banas_photo.gif

Pavel, University of Palacky will provide a short oral from the platform.

1 Regional Centre of Advanced Technologies and Materials
Department of Physical Chemistry
Faculty of Science
Palacký University, tř. 17 listopadu 12, 771 46
Olomouc, Czech Republic

2 Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic

Ph: (+420) 585 634 756
Fx: (+420) 585 634 761
Email: pavel.banas@upol.cz