An optimized charge penetration model for use with the AMOEBA force field.

Journal article

Joshua A. Rackers, Qiantaoa Wang, Chengwen Liu, Jean‐Philip Piquemal, Pengyu Y. Ren, J. Ponder
Physical Chemistry, Chemical Physics - PCCP, 2016

Semantic Scholar DOI PubMed

Cite

APA Click to copy
Rackers, J. A., Wang, Q., Liu, C., Piquemal, J. P., Ren, P. Y., & Ponder, J. (2016). An optimized charge penetration model for use with the AMOEBA force field. Physical Chemistry, Chemical Physics - PCCP.

Chicago/Turabian Click to copy
Rackers, Joshua A., Qiantaoa Wang, Chengwen Liu, Jean‐Philip Piquemal, Pengyu Y. Ren, and J. Ponder. “An Optimized Charge Penetration Model for Use with the AMOEBA Force Field.” Physical Chemistry, Chemical Physics - PCCP (2016).

MLA Click to copy
Rackers, Joshua A., et al. “An Optimized Charge Penetration Model for Use with the AMOEBA Force Field.” Physical Chemistry, Chemical Physics - PCCP, 2016.

BibTeX Click to copy

@article{joshua2016a,
  title = {An optimized charge penetration model for use with the AMOEBA force field.},
  year = {2016},
  journal = {Physical Chemistry, Chemical Physics - PCCP},
  author = {Rackers, Joshua A. and Wang, Qiantaoa and Liu, Chengwen and Piquemal, Jean‐Philip and Ren, Pengyu Y. and Ponder, J.}
}

Abstract

The principal challenge of using classical physics to model biomolecular interactions is capturing the nature of short-range interactions that drive biological processes from nucleic acid base stacking to protein-ligand binding. In particular most classical force fields suffer from an error in their electrostatic models that arises from an ability to account for the overlap between charge distributions occurring when molecules get close to each other, known as charge penetration. In this work we present a simple, physically motivated model for including charge penetration in the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) force field. With a function derived from the charge distribution of a hydrogen-like atom and a limited number of parameters, our charge penetration model dramatically improves the description of electrostatics at short range. On a database of 101 biomolecular dimers, the charge penetration model brings the error in the electrostatic interaction energy relative to the ab initio SAPT electrostatic interaction energy from 13.4 kcal mol-1 to 1.3 kcal mol-1. The model is shown not only to be robust and transferable for the AMOEBA model, but also physically meaningful as it universally improves the description of the electrostatic potential around a given molecule.