Abstract:The electrostatic potential (ESP) is an important property of interactions within and between macromolecules, including those of importance in the life sciences. Semiempirical quantum chemical methods and classical Coulomb calculations fail to provide even qualitative ESP for many of these biomolecules. A new empirical ESP calculation method, namely, EM-ESP, is developed in this study, in which the traditional approach of point atomic charges and the classical Coulomb equation is discarded. In its place, the EM-ESP generates a three-dimensional electrostatic potential V(EM)(r) in molecular space that is the sum of contributions from all component atoms. The contribution of an atom k is formulated as a Gaussian function g(r(k);alpha(k),beta(k)) = alpha(k)/r(k)(betak) with two parameters (alpha(k) and beta(k)). The benchmark for the parameter optimization is the ESP obtained by using higher-level quantum chemical approaches (e.g., CCSD/TZVP). A set of atom-based parameters is optimized in a training set of common organic molecules. Calculated examples demonstrate that the EM-ESP approach is a vast improvement over the Coulombic approach in producing the molecular ESP contours that are comparable to the results obtained with higher-level quantum chemical methods. The atom-based parameters are shown to be transferrable between one part of closely related aromatic molecules. The atom-based ESP formulization and parametrization strategy can be extended to biological macromolecules, such as proteins, DNA, and RNA molecules. Since ESP is frequently used to rationalize and predict intermolecular interactions, we expect that the EM-ESP method will have important applications for studies of protein-ligand and protein-protein interactions in numerous areas of chemistry, molecular biology, and other life sciences.

Energetic Molecule Specific Polarizable Force Field

Development of a Polarizable Force Field Using Multiple Fluctuating Charges per Atom

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

PhyNEO: A Neural-Network-Enhanced Physics-Driven Force Field Development Workflow for Bulk Organic Molecule and Polymer Simulations.

First-Principles-Based Force Field for 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

Physically motivated, robust, ab initio force fields for CO2 and N2.

A fast, convenient, polarizable electrostatic model for molecular dynamics

An Electrostatic Energy-Based Charge Model For Molecular Dynamics Simulation

A Polarizable Valence Electron Density Based Force Field for High-Energy Interactions between Atoms and Molecules

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Molecular Forcefield Methods for Describing Energetic Molecular Crystals: A Review

Development of an anisotropic polarizable model for the all-atom AMOEBA force field

Accurate conformation-dependent molecular electrostatic potentials for high-throughput in silico drug discovery

A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids

A Simple Electron-Density Based Force Field Model for High-Energy Interactions between Atoms and Molecules

AUTOMATED FORCE FIELD PARAMETERIZATION FOR NON-POLARIZABLE AND POLARIZABLE ATOMIC MODELS BASED ON AB INITIO TARGET DATA.

Empirical and Accurate Method for the Three-Dimensional Electrostatic Potential (EM-ESP) of Biomolecules

Polarizable Multipole-Based Force Field for Dimethyl and Trimethyl Phosphate

Benchmarking polarizable and non-polarizable force fields for Ca2+-peptides against a comprehensive QM dataset

All-atom, Non-Empirical, and Tailor-Made Force Field for Α-Rdx from First Principles