On combining Thole's induced point dipole model with fixed charge distributions in molecular mechanics force fields

The Thole induced point dipole model is combined with three different point charge fitting methods, Merz–Kollman (MK), charges from electrostatic potentials using a grid (CHELPG), and restrained electrostatic potential (RESP), and two multipole algorithms, distributed multipole analysis (DMA) and Gaussian multipole model (GMM), which can be used to describe the electrostatic potential (ESP) around molecules in molecular mechanics force fields. This is done to study how the different methods perform when intramolecular polarizability contributions are self‐consistently removed from the fitting done in the force field parametrization. It is demonstrated that the polarizable versions of the partial charge models provide a good compromise between accuracy and computational efficiency in describing the ESP of small organic molecules undergoing conformational changes. For the point charge models, the inclusion of polarizability reduced the the average root mean square error of ESP over the test set by 4–10%. © 2015 Wiley Periodicals, Inc.     

Geometry-dependent atomic charges: Methodology and application to alkanes,aldehydes, ketones,and amides

A general methodology for deriving geometry-dependent atomic charges is presented. The main ingredient of the method is a model that describes the molecular dipole moment in terms of geometry-dependent point charges. The parameters of the model are determined from ab initio calculations of molecular dipole moments and their Cartesian derivatives at various molecular geometries. Transferability of the parameters is built into the model by fitting ab initio calculations for various molecules simultaneously. The results show that charge flux along the bonds is a major contributing factor to the geometry dependence of the atomic charges, with additional contributions from fluxes along valence angles and adjacent bonds. Torsion flux is found to be smaller in magnitude than the bond and valence angle fluxes but is not always unimportant. A set of electrostatic parameters is presented for alkanes, aldehydes, ketones, and amides. Transferability of these parameters for a host of molecules is established to within 3 ?5% error in the predicted dipole moments. A possible extension of the method to include atomic dipoles is outlined. With the inclusion of such atomic dipoles and with the set of transferable point charges and charge flux parameters, it is demonstrated that molecular electrostatic potentials as well as electrostatic forces on nuclei can be reproduced much better than is possible with other models (such as potential derived charges). © 1995 by John Wiley & Sons, Inc.     

Point charge representation of multicenter multipole moments in calculation of electrostatic properties

Summary Distributed Point Charge Models (PCM) for CO, (H₂O)₂, and HS-SH molecules have been computed from analytical expressions using multicenter multipole moments. The point charges (set of charges including both atomic and non-atomic positions) exactly reproduce both molecular and segmental multipole moments, thus constituting an accurate representation of the local anisotropy of electrostatic properties. In contrast to other known point charge models, PCM can be used to calculate not only intermolecular, but also intramolecular interactions. Comparison of these results with more accurate calculations demonstrated that PCM can correctly represent both weak and strong (intramolecular) interactions, thus indicating the merit of extending PCM to obtain improved potentials for molecular mechanics and molecular dynamics computational methods.Dedicated to Prof. Alberte PullmanPacific Northwest Laboratory is operated for the US Department of Energy by Battelle Memorial Institute under contract DE-ACO6-76RLO 1830     

Halogen bonding: a theoretical study based on atomic multipoles derived from quantum theory of atoms in molecules

Atomic multipole moments derived from quantum theory of atoms in molecules are used to study halogen bonds in dihalogens (with general formula YX, in which X refers to the halogen directly interacted with the Lewis base) and some molecules containing C–X group. Multipole expansion is used to calculate the electrostatic potential in a vicinity of halogen atom (which is involved in halogen bonding) in terms of atomic monopole, dipole, and quadrupole moments. In all the cases, the zz component of atomic traceless quadrupole moments (where z axis taken along Y–X or C–X bonds) of the halogens plays a stabilizing role in halogen bond formation. The effects of atomic monopole and dipole moments on the formation of a halogen bond in YX molecules depend on Y and X atoms. In Br₂ and Cl₂, the monopole moment of halogens is zero and has no contribution in electrostatic potential and hence in halogen bonding, while in ClBr, FBr, and FCl it is positive and therefore stabilize the halogen bonds. On the other hand, the negative sign of dipole moment of X in all the YX molecules weakens the corresponding halogen bonds. In the C–X-containing molecules, monopole and dipole moments of X atom are negative and consequently destabilize the halogen bonds. So, in these molecules the quadrupole moment of X atom is the only electrostatic term which strengthens the halogen bonds. In addition, we found good linear correlations between halogen bonds strength and electrostatic potentials calculated from multipole expansion.     

High Resolution Synchrotron Diffraction Study on Charge Density Distribution of Ampicillin Trihydrate

Charge density distribution in ampicillin trihydrate was investigated experimentally. Results were compared with the quantum calculations using density functional theory. The charge derived properties including Mulliken atomic charges, dipole moment, and molecular electrostatic potential were calculated. The multipole analysis was done for the refinement of experimental population parameters. The structure factors obtained from multipole treatment were used for the construction of Fourier maps. Topological properties of the charge distribution were discussed and the characteristics of (3,-1) critical points were analyzed.     

Multipole correction of atomic monopole models of molecular charge distribution. I. Peptides

The defects in atomic monopole models of molecular charge distribution have been analyzed for several model-blocked peptides and compared with accurate quantum chemical values. The results indicate that the angular characteristics of the molecular electrostatic potential around functional groups capable of forming hydrogen bonds can be considerably distorted within various models relying upon isotropic atomic charges only. It is shown that these defects can be corrected by augmenting the atomic point charge models by cumulative atomic multipole moments (CAMMs). Alternatively, sets of off-center atomic point charges could be automatically derived from respective multipoles, providing approximately equivalent corrections. For the first time, correlated atomic multipoles have been calculated for N-acetyl, N'-methylamide-blocked derivatives of glycine, alanine, cysteine, threonine, leucine, lysine, and serine using the MP2 method. The role of the correlation effects in the peptide molecular charge distribution are discussed.     

Transferable atom equivalent multicentered multipole expansion method

The transferability of atomic and functional group properties is an implicit concept in chemistry. The work presented here describes the use of Transferable Atom Equivalents (TAE) to represent molecular electrostatic potential fields through the use of integrated atomic multipole moments that are associated with each TAE atom type used in the reconstruction. TAE molecular surface distributions of electrostatic potentials are compared with analytical ab initio and empirical (Gasteiger) partial charge reference models for several conformations of test peptides. Surface electrostatic potential distributions computed using TAE multipole representations were found to converge at the octopole level, with incremental improvement observed when hexadecapoles were included. Molecular electrostatic potential fields that were produced using the TAE method were observed to be responsive to conformational changes and to compare well with ab initio reference distributions. Generation of TAE atom types and their associated multipoles does not involve fitting to sample electrostatic potential fields, but rather utilizes integrated AIM atomic electron density distributions within representative chemical environments. The RECON program was used for TAE reconstruction. RECON is capable of processing 5,000 drug-sized molecules or 25 proteins per minute per 1.7 GHz P4 Linux processor.     

Atomic forces for geometry‐dependent point multipole and Gaussian multipole models

In standard treatments of atomic multipole models, interaction energies, total molecular forces, and total molecular torques are given for multipolar interactions between rigid molecules. However, if the molecules are assumed to be flexible, two additional multipolar atomic forces arise because of (1) the transfer of torque between neighboring atoms and (2) the dependence of multipole moment on internal geometry (bond lengths, bond angles, etc.) for geometry‐dependent multipole models. In this study, atomic force expressions for geometry‐dependent multipoles are presented for use in simulations of flexible molecules. The atomic forces are derived by first proposing a new general expression for Wigner function derivatives . The force equations can be applied to electrostatic models based on atomic point multipoles or Gaussian multipole charge density. Hydrogen‐bonded dimers are used to test the intermolecular electrostatic energies and atomic forces calculated by geometry‐dependent multipoles fit to the ab initio electrostatic potential. The electrostatic energies and forces are compared with their reference ab initio values. It is shown that both static and geometry‐dependent multipole models are able to reproduce total molecular forces and torques with respect to ab initio, whereas geometry‐dependent multipoles are needed to reproduce ab initio atomic forces. The expressions for atomic force can be used in simulations of flexible molecules with atomic multipoles. In addition, the results presented in this work should lead to further development of next generation force fields composed of geometry‐dependent multipole models. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

SCF Deformation densities and electrostatic potentials of purines and pyrimidines

4-31G wave functions have been computed for five purines and pyrimidines. The calculated deformation densities have been partitioned into atomic fragments, which were integrated to yield atomic multipole moments. The transferability of atomic fragments between related molecules was verified by constructing model maps for uracil and guanine from appropriate fragments of cytosine and adenine. Model electrostatic potentials calculated from the moments of model atoms are similar to the corresponding 4-31G potentials. Comparison of 4-31G and 4-31G** deformation densities of cytosine provides simple rules for estimating the effects of polarization functions on the atomic multipole moments of most atom types occurring in the purines and pyrimidines. These rules were applied to the other molecules and yielded reasonable approximations for their molecular dipole moments. Substituting CH₃ for H has little effect on the deformation density beyond the substitution center.     

Comparison of atomic charges derived via different procedures

Atomic charges were obtained from ab initio molecular orbital calculations using a variety of procedures to compare them and assess their utility. Two procedures based on the molecular orbitals were examined, the Mulliken population analysis and the Weinhold–Reed Natural Population Analysis. Two procedures using the charge density distribution were included; the Hirshfeld procedure and Bader's Atoms in Molecules method. Charges also were derived by fitting the electrostatic potential (CHELPG) and making use of the atomic polar tensors (GAPT). The procedures were first examined for basis set independence, and then applied to a group of hydrocarbons. The dipole moments for these molecules were computed from the various atomic charges and compared to the total SCF dipole moments. This was followed by an examination of a series of substituted methanes, simple hydrides, and a group of typical organic compounds such as carbonyl derivatives, nitriles, and nitro compounds. In some cases, the ability of the charges to reproduce electrostatic potentials was examined. © John Wiley & Sons, Inc.     

Gaussian induced dipole polarization model

A new induced dipole polarization model based on interacting Gaussian charge densities is presented. In contrast to the original induced point dipole model, the Gaussian polarization model is capable of finite interactions at short distances. Aspects of convergence related to the Gaussian model will be explored. The Gaussian polarization model is compared with the damped Thole-induced dipole model and the point dipole model. It will be shown that the Gaussian polarization model performs slightly better than the Thole model in terms of fitting to molecular polarizability tensors. An advantage of the model based on Gaussian charge distribution is that it can be easily generalized to other multipole moments and provide effective damping for both permanent electrostatic and polarization models. Finally, a method of parameterizing polarizabilities is presented. This method is based on probing a molecule with point charges and fitting polarizabilities to electrostatic potential. In contrast to the generic atom type polarizabilities fit to molecular polarizability tensors, probed polarizabilities are significantly more accurate in terms of reproducing molecular polarizability tensors and electrostatic potential, while retaining conformational transferability.     

Atomic multipoles: Electrostatic potential fit,local reference axis systems,and conformational dependence

Currently, all standard force fields for biomolecular simulations use point charges to model intermolecular electrostatic interactions. This is a fast and simple approach but has deficiencies when the electrostatic potential (ESP) is compared to that from ab initio methods. Here, we show how atomic multipoles can be rigorously implemented into common biomolecular force fields. For this, a comprehensive set of local reference axis systems is introduced, which represents a universal solution for treating atom‐centered multipoles for all small organic molecules and proteins. Furthermore, we introduce a new method for fitting atomic multipole moments to the quantum mechanically derived ESP. This methods yields a 50–90% error reduction compared to both point charges fit to the ESP and multipoles directly calculated from the ab initio electron density. It is shown that it is necessary to directly fit the multipole moments of conformational ensembles to the ESP. Ignoring the conformational dependence or averaging over parameters from different conformations dramatically deteriorates the results obtained with atomic multipole moments, rendering multipoles worse than partial charges. © 2012 Wiley Periodicals, Inc.     

Off‐atomic site charges in some anions,metal cations,and complexes: Significance for electrostatic properties and binding

Electronic structures and properties of several anions, metal cations, and their complexes with neutral molecules were investigated at the HF/6‐31G** and B3LYP/6‐31G** levels of theory. Charges shifted from atomic sites due to atomic orbital hybridization called hybridization displacement charges (HDC) were investigated in detail. It has been found that many components of HDC are associated with each atom of ion that are shifted from the atomic sites, those associated with metal cations being shifted by large distances as found previously in electrically neutral systems. It is shown that atomic orbitals are appreciably rehybridized in going from neutral molecules to anions and cations. Molecular dipole moments and surface molecular electrostatic potentials (MEP) are obtained satisfactorily using HDC for the various types of species mentioned above. In the OH^?? H₂O complex, reversal of direction of shift of an HDC component associated with the hydrogen atom of H₂O involved in hydrogen bonding, indicates that the hydrogen bond between OH^? and H₂O would have some covalent character. Other atomic site‐based point charge models cannot provide such information about the nature of bonding. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem 2007     

Base stacking in cytosine dimer. A comparison of correlated ab initio calculations with three empirical potential models and density functional theory calculations

Ab initio MP2/6-31G* interaction energies were calculated for more than 80 geometries of stacked cytosine dimer. Diffuse polarization functions were used to properly cover the dispersion energy. The results of ab initio calculations were compared with those obtained from three electrostatic empirical potential models, constructed as the sum of a Lennard-Jones potential (covering dispersion and repulsion contributions) and the electrostatic term. Point charges and point multipoles of the electrostatic term were also obtained at the MP2/6-31G* level of theory. The point charge MEP model (atomic charges derived from molecular electrostatic potential) satisfactorily reproduced the ab initio data. Addition of π-charges localized below and above the cytosine plane did not affect the calculated energies. The model employing the distributed multipole analysis gave worse agreement with the ab initio data than the MEP approach. The MP2 MEP charges were also derived using larger sets of atomic orbitals: cc-pVDZ, 6-311 + G(2d, p), and aug-cc-pVDZ. Differences between interaction energies calculated using these three sets of point charges and the MP2/6-31G* charges were smaller than 0.8 kcal/mol. The correlated ab initio calculations were also compared with the density functional theory (DFT) method. DFT calculations well reproduced the electrostatic part of interaction energy. They also covered some nonelectrostatic short-range effects which were not reproduced by the empirical potentials. The DFT method does not include the dispersion energy. This energy, approximated by an empirical term, was therefore added to the DFT interaction energy. The resulting interaction energy exhibited an artifact secondary minimum for a 3.9-4.0 vertical separation of bases. This defect is inherent in the DFT functionals, because it is not observed for the Hartree-Fock + dispersion interaction energy.© 1996 John Wiley & Sons, Inc.