The effect of isodensity surface sampling on ESP derived charges and the effect of adding bondcenters on DMA derived charges

The effect of sampling the electrostatic potential around a molecule on the quality of electrostatic potential derived charges is investigated. In addition, the effect of the number of expansion sites in a Distributed Multipole Analysis (DMA) on the quality of charges fitted to the DMA derived electrostatic potential is investigated. Sampling on constant electron density surfaces gives a better fit between the quantum mechanical potential and the potential derived from the fitted charges, compared to sampling on a van der Waals surface composed of intersecting spheres. The fit between the electrostatic potential derived from point charges and the quantum mechanical potential becomes poorer with increasing quality of the employed basis set. The inclusion of bondcenters into the calculations improves the fit between the Quantum Mechanical (QM) electrostatic potential and the DMA derived potential. The number of expansion sites needed for an accurate approximation of the QM electrostatic potential increases with increasing quality of the used basis set.     

Generation of OPLS-like charges from molecular electrostatic potential using restraints

A new method for generating atom-centered charges for use in condensed phase computer simulations is presented, which is based on a restrained electrostatic potential (RESP) procedure. Charges are calculated from a least-squares fit to the quantum mechanical electrostatic potential with a restraint applied to reduce their magnitude. The restraint developed here offers advantages over that used in RESP. The magnitude of the restraint is optimized to yield charges as close to the equivalent OPLS values as possible while still reproducing the molecule's electrostatic potential. A cross-validation analysis is used to show that the restraint is insensitive to the selection of OPLS molecules from which it is derived. Thus, with this method, OPLS-like charges may be produced from the electrostatic potential for atom types not in the OPLS force field. In addition, the restraint is shown to reduce the conformational dependence of the charges. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 483–498, 1999     

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

The atom‐centered partial charges‐approximation is commonly used in current molecular modeling tools as a computationally inexpensive alternative to quantum mechanics for modeling electrostatics. Even today, the use of partial charges remains useful despite significant advances in improving the efficiency of ab initio methods. Here, we report on new parameters for the EEM and SFKEEM electronegativity equalization‐based methods for rapidly determining partial charges that will accurately model the electrostatic potential of flexible molecules. The developed parameters cover most pharmaceutically relevant chemistries, and charges obtained using these parameters reproduce the B3LYP/cc‐pVTZ reference electrostatic potential of a set of FDA‐approved drug molecules at best to an average accuracy of 13 ± 4 kJ mol^?1; thus, equipped with these parameters electronegativity equalization‐based methods rival the current best non‐quantum mechanical methods, such as AM1‐BCC, in accuracy, yet incur a lower computational cost. Software implementations of EEM and SFKEEM, including the developed parameters, are included in the conformer‐generation tool BALLOON , available free of charge at http://web.abo.fi/fak/mnf/bkf/research/johnson/software.php . © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

A practical procedure for the determination of electrostatic charges of large molecules

Summary A practical procedure for the precise determination of electrostatic charges, which are evaluated by fitting the rigorous quantum mechanical molecular electrostatic potential to a monopole-monopole expression, is presented. The proposal of this procedure arises from the study of the minimum requirements necessary to obtain reliable electrostatic charges. Such a study is focused on: (i) the dependence of the electrostatic charges on the set of points where the quantum mechanical and the monopole-monopole molecular electrostatic potentials are fitted; thus, both the influence of the number of points and their distribution in layers located out of the van der Waals radii of the atoms are examined, and (ii) the reliability of the use of fractional models for the evaluation of electrostatic charges of large molecules. Results point out that the optimum number of points is defined by a density of points ranging from 0.45 to 0.60 points/Ų when four layers (separated by 0.2 Å) are considered. Nevertheless, the use of only two layers (separated by 0.4 Å) for large molecules is recommended, thus enabling one to obtain reliable charges at a reduced computational cost. Moreover, results justify the use of fractional models for the determination of electrostatic charges of extremely large molecules, even when aromatic structures exist.     

Conformational dependence of charges in protein simulations

We have studied the conformational dependence of molecular mechanics atomic charges for proteins by calculating the charges fitted to the quantum mechanical (QM) electrostatic potential (ESP) for all atoms in complexes between avidin and seven biotin analogues for 20 snapshots from molecular dynamics simulations. We have studied how various other charge sets reproduce those charges. The QM charges, even if averaged over all snapshots or all residues, in general have a larger magnitude than standard Amber charges, indicating that the restraint toward zero in the restrained ESP method is too strong. This has a significant influence on the electrostatic conformational energies and the interaction energy between the biotin ligand and the protein, giving a difference between the QM and Amber charges of 43 and 8 kJ/mol for the negatively charged and neutral biotin analogues, respectively (3-4%). However, this energy difference is strongly reduced if the solvation energy (calculated by the Poisson-Boltzmann or Generalized Born methods) is added, viz., to 7 kJ/mol for charged and 3 kJ/mol for uncharged ligand. In fact, charges need to be recalculated with a QM method only for residues within 7 or 4 A of the ligand, if the error should be less than 4 kJ/mol. Unfortunately, the QM charges do not give significantly better MM/PBSA estimates of ligand-binding affinities than standard Amber charges.     

Fast tools for calculation of atomic charges well suited for drug design1

Two novel approaches to construct empirical schemes for partial atomic charge calculation were proposed. The charge schemes possess important benefits. First, they produce both topologically symmetrical and environment dependent charges. Second, they can be parameterised to reasonably reproduce ab initio molecular electrostatic potential (MEP), which guarantees their successful use in molecular modelling. To validate the approaches, the parameters of the proposed charge schemes were fitted to best reproduce MEP simultaneously on grids around a set of 227 diverse organic compounds. The residual errors in MEP reproduction due to calculated atomic charges were compared to those due to charges from known charge schemes.     

Consensus bond-charge increments fitted to electrostatic potential or field of many compounds: Application to MMFF94 training set

Bond-charge increments (BCIs) are additive parameters used to assign atomic charges for the MMFF force field. BCI parameters are classified parsimoniously according to two atom types and the bond order. We show how BCIs may be fitted rapidly by linear least squares to the calculated ab initio electrostatic potential (ESP) or to the electrostatic field. When applied simultaneously to a set of compounds or conformations, the method yields consensus values of the BCIs. The method can also derive conventional “ESP-fit” atomic charges with improved numerical stability. The method may be generalized to determine atom multipoles, multicenter charge templates, or electronegativities, but not polarizability or hardness. We determine 65 potential-derived (PD) BCI parameters, which are classified as in MMFF, by fitting the 6-31G* ESP or the electrostatic field of the 45 compounds in the original MMFF94 training set. We compare the consensus BCIs with classified BCIs that were fit to each molecule individually and with “unique-bond” BCIs (ESP-derived atom charges). Consensus BCIs give a satisfactory representation for about half of the structures and are robust to the adjustment of the alkyl CH bond increment to the zero value employed in MMFF94. We highlight problems at three levels: Point approximation: the potential near lone pairs on sulfur and to some extent nitrogen cannot be represented just by atom charges. Bond classification: BCIs classified according to MMFF atom types cannot represent all delocalized electronic effects. The problem is especially severe for bonds between atoms of equivalent MMFF type, whose BCI must be taken as zero. Consensus: discrepancies that occur in forming the consensus across the training set indicate the need for a more detailed classification of BCIs. Contradictions are seen (e.g., between acetic acid and acetone and between guanidine and formaldehydeimine). We then test the three sets of PD-BCIs in energy minimizations of hydrogen-bonded dimers. Unique-bond BCIs used with the MMFF buffered 14–7 potential reproduce unscaled quantum chemical dimer interaction energies within 0.9 kcal/mol root mean square (or 0.5, omitting two N-oxides). These energies are on average 0.7 (or 0.5) kcal/mol too weak to reproduce the scaled quantum mechanical (SQM) results that are a benchmark for MMFF parameterization. Consensus BCIs tend to weaken the dimer energy by a further 0.4–0.6 kcal/mol. Thus, consensus PD-BCIs can serve as a starting point for MMFF parameterization, but they require both systematic and individual adjustments. Used with a “harder” AMBER-like Lennard–Jones potential, unique-bond PD-BCIs without systematic adjustment give dimer energies in fairly good agreement with SQM. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1495–1516, 1999     

Restrained point-charge models for disaccharides

Various methods for deriving atomic partial charges from the quantum chemical electrostatic potential and moments have been tested for the sucrose molecule. We show that if no further information is used, the charges on some carbon atoms become large and charge patterns involving these atoms are badly determined and poorly transferable. Adding lone-pairs on the ether oxygen atoms or dividing the molecule into smaller fragments did not cure the instabilities. We develop a method, CHELP-BOW0, that restrains charges toward zero with different weights for different atoms. These harmonic restraints preserve the linear form of the least-squares equations, which are solved in a single step using singular-value decomposition. CHELP-BOW0 improves the chemical transferability of the charges compared to unrestrained methods, and slightly improves their conformational transferability. It introduces a modest degradation of the fit compared to unrestrained CHELP-BOW (mean average deviation of the potential 0.00016 vs. 0.00010 a.u.). A second new method, CHELP-BOWC, avoids the need for restraints by including several conformations in the fit, weighting each according to its estimated energy in solution. CHELP-BOWC charges are more transferable than CHELP-BOW or CHELP-BOW0 charges to conformations not included in the training set. Restraints to zero charge do not further improve transferability of the CHELP-BOWC charges. We, therefore, recommend CHELP-BOW charges for rigid molecules and CHELP-BOWC charges for flexible molecules.     

Application of error-ranked singular value decomposition for the determination of potential-derived atomic-centered point charges

The present work provides a detailed investigation on the use of singular value decomposition (SVD) to solve the linear least-squares problem (LLS) for the purposes of obtaining potential-derived atom-centered point charges (PD charges) from the ab initio molecular electrostatic potential (V(QM)). Given the SVD of any PD charge calculation LLS problem, it was concluded that (1) all singular vectors are not necessary to obtain the optimal set of PD charges and (2) the most effective set of singular vectors do not necessarily correspond to those with the largest singular values. It is shown that the efficient use of singular vectors can provide statistically well-defined PD charges when compared with conventional PD charge calculation methods without sacrificing the agreement with V(QM). As can be expected, the methodology outlined here is independent of the algorithm for sampling V(QM) as well as the basis set used to calculate V(QM). An algorithm is provided to select the best set of singular vectors used for optimal PD charge calculations. To minimize the subjective comparisons of different PD charge sets, we also provide an objective criterion for determining if two sets of PD charges are significantly different from one another.     

Optimization of a molecular mechanics force field for type-II polyoxometalates focussing on electrostatic interactions: a case study

In this study, we have focussed on type-II polyanions such as [M(7)O(24)](6-), and we have developed and validated optimized force fields that include electrostatic and van der Waals interactions. These contributions to the total steric energy are described by the nonbonded term, which encompasses all interactions between atoms that are not transmitted through the bonds. A first validation of a stochastic technique based on genetic algorithms was previously made for the optimization of force fields dedicated to type-I polyoxometalates. To describe the new nonbonded term added in the functional, a fixed-charged model was chosen. Therefore, one of the main issues was to analyze that which partial atomic charges could be reliably used to describe these interactions in such inorganic compounds. Based on several computational strategies, molecular mechanics (MM) force field parameters were optimized using different types of atomic charges. Moreover, the influence of the electrostatic and van der Waals buffering constants and 1,4-interactions scaling factors used in the force field were also tested, either being optimized as well or fixed with respect to the values of CHARMM force field. Results show that some atomic charges are not well adapted to CHARMM parameters and lead to unrealistic MM-optimized structures or a MM divergence. As a result, a new scaling factor has been optimized for Quantum Theory of Atoms in Molecules charges and charges derived from the electrostatic potential such as ChelpG. The force fields optimized can be mixed with the CHARMM force field, without changing it, to study for the first time hepta-anions interacting with organic molecules.     

Toward a new approach for determination of solute's charge distribution to analyze interatomic electrostatic interactions in quantum mechanical/molecular mechanical simulations

We present an alternative approach to determine "density-dependent property"-derived charges for molecules in the condensed phase. In the case of a solution, it is essential to take into consideration the electron polarization of molecules in the active site of this system. The solute and solvent molecules in this site have to be described by a quantum mechanical technique and the others are allowed to be treated by a molecular mechanical method (QM/MM scheme). For calculations based on this scheme, using the forces and interaction energy as density-dependent property our charges from interaction energy and forces (CHIEF) approach can provide the atom-centered charges on the solute atoms. These charges reproduce well the electrostatic potentials around the solvent molecules and present properly the picture of the electron density of the QM subsystem in the solution system. Thus, the CHIEF charges can be considered as the atomic charges under the conditions of the QM/MM simulation, and then enable one to analyze electrostatic interactions between atoms in the QM and MM regions. This approach would give a view of the QM nuclei and electrons different from the conventional methods.     

Suitability of density functional methods for calculation of electrostatic properties

A systematic analysis was performed on the suitability of the molecular electrostatic potential (MEP) and MEP-derived properties determined by means of density functional (DFT) methods. Attention was paid to the electrostatic potential (ESP) derived charges, the ESP and exact quantum mechanical dipole moments, the depth of MEP minima, and the MEP distribution in layers around the molecule for a large series of molecules. The electrostatic properties were determined at either local or nonlocal DFT levels using different functionals. The results were compared with the values estimated from quantum mechanical calculations performed at Hartree–Fock, Møller–Plesset up to fourth order, and CIPSI levels. The suitability of the MEP-derived properties estimated from DFT methods is discussed for application in different areas of chemical interest. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 980–991, 1997     

Electric potential distributions around discrete charges in a dielectric membrane—electrolyte solution system

Within the framework of the linearized Debye-Hückel theory an exact solution of the problem of calculating the electric potential caused by discrete fixed charges located at arbitrary positions with respect to a dielectric membrane-solution interface is presented. It takes into account the existence of an electrolyte solution on both sides of the membrane. Asymmetric ionic conditions are allowed for. For some interesting typical cases of fixed charge locations and electrolyte ionic strengths electric potential distributions were calculated and discussed. It is shown that, if the fixed charges were at or in front of the membrane surface, the characteristic distance of the electric potential decay was comparable to the Debye-Hückel length. At the opposite membrane surface only very small electric potentials can be observed. If, however, the fixed charge was placed below the membrane surface the electric potential in lateral direction and towards the other membrane surface largely increased. This effect was very sensitive to the position of the fixed charge with respect to the membrane surface.     

Constructing irregular surfaces to enclose macromolecular complexes for mesoscale modeling using the discrete surface charge optimization (DISCO) algorithm

Salt-mediated electrostatics interactions play an essential role in biomolecular structures and dynamics. Because macromolecular systems modeled at atomic resolution contain thousands of solute atoms, the electrostatic computations constitute an expensive part of the force and energy calculations. Implicit solvent models are one way to simplify the model and associated calculations, but they are generally used in combination with standard atomic models for the solute. To approximate electrostatics interactions in models on the polymer level (e.g., supercoiled DNA) that are simulated over long times (e.g., milliseconds) using Brownian dynamics, Beard and Schlick have developed the DiSCO (Discrete Surface Charge Optimization) algorithm. DiSCO represents a macromolecular complex by a few hundred discrete charges on a surface enclosing the system modeled by the Debye-Hückel (screened Coulombic) approximation to the Poisson-Boltzmann equation, and treats the salt solution as continuum solvation. DiSCO can represent the nucleosome core particle (>12,000 atoms), for example, by 353 discrete surface charges distributed on the surfaces of a large disk for the nucleosome core particle and a slender cylinder for the histone tail; the charges are optimized with respect to the Poisson-Boltzmann solution for the electric field, yielding a approximately 5.5% residual. Because regular surfaces enclosing macromolecules are not sufficiently general and may be suboptimal for certain systems, we develop a general method to construct irregular models tailored to the geometry of macromolecules. We also compare charge optimization based on both the electric field and electrostatic potential refinement. Results indicate that irregular surfaces can lead to a more accurate approximation (lower residuals), and the refinement in terms of the electric field is more robust. We also show that surface smoothing for irregular models is important, that the charge optimization (by the TNPACK minimizer) is efficient and does not depend on the initial assigned values, and that the residual is acceptable when the distance to the model surface is close to, or larger than, the Debye length. We illustrate applications of DiSCO's model-building procedure to chromatin folding and supercoiled DNA bound to Hin and Fis proteins. DiSCO is generally applicable to other interesting macromolecular systems for which mesoscale models are appropriate, to yield a resolution between the all-atom representative and the polymer level.     

Fast,efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method

The AM1‐BCC method quickly and efficiently generates high‐quality atomic charges for use in condensed‐phase simulations. The underlying features of the electron distribution including formal charge and delocalization are first captured by AM1 atomic charges for the individual molecule. Bond charge corrections (BCCs), which have been parameterized against the HF/6‐31G* electrostatic potential (ESP) of a training set of compounds containing relevant functional groups, are then added using a formalism identical to the consensus BCI (bond charge increment) approach. As a proof of the concept, we fit BCCs simultaneously to 45 compounds including O‐, N‐, and S‐containing functionalities, aromatics, and heteroaromatics, using only 41 BCC parameters. AM1‐BCC yields charge sets of comparable quality to HF/6‐31G* ESP‐derived charges in a fraction of the time while reducing instabilities in the atomic charges compared to direct ESP‐fit methods. We then apply the BCC parameters to a small “test set” consisting of aspirin, d ‐glucose, and eryodictyol; the AM1‐BCC model again provides atomic charges of quality comparable with HF/6‐31G* RESP charges, as judged by an increase of only 0.01 to 0.02 atomic units in the root‐mean‐square (RMS) error in ESP. Based on these encouraging results, we intend to parameterize the AM1‐BCC model to provide a consistent charge model for any organic or biological molecule. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 132–146, 2000     

Efficient and accurate solvation energy calculation from polarizable continuum models

A new approach is proposed to enhance the efficiency and accuracy for calculation of the long-range electrostatic interaction from implicit solvation models, i.e., the polarizable continuum model (PCM) and its variants, conductorlike PCM/conductorlike screening model and integral equation formalism PCM. In these methods the solvent electrostatics effects are represented by a set of discrete apparent charges distributed on tesserae of the molecular cavity surface embedding the solute. In principle, the accuracy of these methods is improved if the cavity surface is tessellated to finer tesserae; however, the computational time is increased rapidly. We show that such undesired dependency between accuracy and efficiency is a result of the inaccurate treatment of the apparent charge self-contribution to the potential and/or electric field. By taking into account the full effects due to the size and curvature of the segment occupied by each apparent charge, the error in calculated electrostatic solvation free energy is essentially zero for ions (point charge at the center of a sphere) regardless of the degree of tessellation used. For molecules where gradient of apparent charge density is nonzero at the cavity surface, we propose a multiple-sampling technique which significantly lowers the calculated error compared to the original PCM methods, especially when very few numbers of tesserae are used.     

Partition coefficients of methylated DNA bases obtained from free energy calculations with molecular electron density derived atomic charges

Partition coefficients serve in various areas as pharmacology and environmental sciences to predict the hydrophobicity of different substances. Recently, they have also been used to address the accuracy of force fields for various organic compounds and specifically the methylated DNA bases. In this study, atomic charges were derived by different partitioning methods (Hirshfeld and Minimal Basis Iterative Stockholder) directly from the electron density obtained by electronic structure calculations in a vacuum, with an implicit solvation model or with explicit solvation taking the dynamics of the solute and the solvent into account. To test the ability of these charges to describe electrostatic interactions in force fields for condensed phases, the original atomic charges of the AMBER99 force field were replaced with the new atomic charges and combined with different solvent models to obtain the hydration and chloroform solvation free energies by molecular dynamics simulations. Chloroform–water partition coefficients derived from the obtained free energies were compared to experimental and previously reported values obtained with the GAFF or the AMBER‐99 force field. The results show that good agreement with experimental data is obtained when the polarization of the electron density by the solvent has been taken into account, and when the energy needed to polarize the electron density of the solute has been considered in the transfer free energy. These results were further confirmed by hydration free energies of polar and aromatic amino acid side chain analogs. Comparison of the two partitioning methods, Hirshfeld‐I and Minimal Basis Iterative Stockholder (MBIS), revealed some deficiencies in the Hirshfeld‐I method related to the unstable isolated anionic nitrogen pro‐atom used in the method. Hydration free energies and partitioning coefficients obtained with atomic charges from the MBIS partitioning method accounting for polarization by the implicit solvation model are in good agreement with the experimental values. © 2018 Wiley Periodicals, Inc.     

Ab initio calculation of the interaction energy in the P2 binding pocket of HIV‐1 protease

Interactions at the P2 binding pocket of human immunodeficiency virus type 1 (HIV‐1) protease have been studied using calculated interaction energies for model systems that mimic this binding pocket. Models were built for the P2 pocket of HIV‐1 protease in complex with TMC114, nelfinavir, and amprenavir. A two‐step procedure was applied. In the first step, the size of the model system was confined to ～40 atoms, and the interaction energy was calculated at different computational levels. In the second step, the size of the system was increased to 138 atoms, and the calculations were only performed at the HF/6‐31G** level. The interaction energy of the HIV‐1 protease/TMC114 complex was found to be more favorable than the interaction energies of the other complexes because of the additional hydrogen bond interaction this inhibitor is able to make with the HIV‐1 protease backbone. The results of the calculations are supported by stockholder charges and electrostatic potential maps. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005