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.     

A finite field method for calculating molecular polarizability tensors for arbitrary multipole rank

A finite field method for calculating spherical tensor molecular polarizability tensors α_lm;l′m′ = ?Δ_lm/??_l′m′* by numerical derivatives of induced molecular multipole Δ_lm with respect to gradients of electrostatic potential ?_l′m′* is described for arbitrary multipole ranks l and l′. Interconversion formulae for transforming multipole moments and polarizability tensors between spherical and traceless Cartesian tensor conventions are derived. As an example, molecular polarizability tensors up to the hexadecapole–hexadecapole level are calculated for water using the following ab initio methods: Hartree–Fock (HF), Becke three‐parameter Lee‐Yang‐Parr exchange‐correlation functional (B3LYP), Møller–Plesset perturbation theory up to second order (MP2), and Coupled Cluster theory with single and double excitations (CCSD). In addition, intermolecular electrostatic and polarization energies calculated by molecular multipoles and polarizability tensors are compared with ab initio reference values calculated by the Reduced Variation Space method for several randomly oriented small molecule dimers separated by a large distance. It is discussed how higher order molecular polarizability tensors can be used as a tool for testing and developing new polarization models for future force fields. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Gaussian multipoles in practice: Electrostatic energies for intermolecular potentials

A method is presented for calculating the total electrostatic interaction energies between molecules from ab initio monomer wave functions. This approach differs from existing methods, such as Stone's distributed multipole analysis (DMA), in including the short-range penetration energy as well as the long-range multipolar energy. The monomer charge densities are expressed as distributed series of atom-centered functions which we call Gaussian multipoles; these are analogous to the distributed point multipoles used in DMA. Our procedure has been encoded in the GMUL program. Calculations have been performed on the formamide/formaldehyde complex, a model system for N? H …? O hydrogen bonding in biological molecules, and also on guanidinium/benzene, modeling amino/aromatic interactions in proteins. We find that the penetration energy can be significant, especially in its contribution to the variation of the electrostatic energy with interaction geometry. A hybrid method, which uses Gaussian multipoles for short-range atom pair interactions and point multipoles for long-range ones, allows the electrostatic energies, including penetration, to be calculated at a much reduced cost. We also note that the penetration energy may provide the best route to an atom–atom anisotropic model for the exchange-repulsion energy in intermolecular potentials. © 1994 by John Wiley & Sons, Inc.     

Rules for minimal atomic multipole expansion of molecular fields

A nonempirical minimal atomic multipole expansion (MAME) defines atomic charges or higher multipoles that reproduce electrostatic potential outside molecules. MAME eliminates problems associated with redundancy and with statistical sampling, and produces atomic multipoles in line with chemical intuition.     

Intramolecular chelation of Zn2+ by α‐ and β‐mercaptocarboxamides. A parallel ab initio and polarizable molecular mechanics investigation. Assessment of the role of multipole transferability

α‐ and β‐mercaptocarboxamides constitute the Zn²⁺‐ligating entity of several highly potent metalloenzyme inhibitors. We have studied their interaction energies with Zn²⁺ using the polarizable molecular mechanics procedure SIBFA, and compared them to the corresponding ab initio supermolecule ones. Such validations are necessary to subsequently undertake simulations on complexes of Zn²⁺–metalloenzymes with inhibitors. If the distributed multipoles and polarizabilities are those derived for each ligand in its appropriate Zn²⁺‐binding conformation, a close reproduction of the ab initio binding energies is afforded. However, this representation is not tractable upon increasing the size of the ligands and/or to explore a continuum of binding conformations. This makes it necessary to construct the ligands by resorting to a library of constitutive fragments, namely in this case methanethiolate, formamide, and methane covalently connected together. A close reproduction of the ab initio interaction energies is enabled, but only if the ligand–ligand interactions are computed simultaneously with those occurring with Zn²⁺. This representation accounts for the nonadditivity occurring in the Zn²⁺–methanethiolate–formamide complex, and justifies the use of the distributed multipoles on the fragments for the construction of larger and flexible molecules. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1038–1047, 2001     

Cumulative coordinates for approximations of high‐order atomic multipole moments in aluminosilicate and aluminophosphate sieves*

A scheme to obtain approximate analytical functions for the atomic distributed multipole moments of the crystallographically different atoms within aluminosilicate and aluminophosphate sieves is discussed. Respective atomic multipole moments are derived within the CRYSTAL95 ab initio periodic Hartree–Fock code with different basis sets, from minimal STO‐3G to 6‐21G*. In order to illustrate the possible applications, distributed analyses are carried out for various structural models from all‐siliceous zeolites and aluminophosphates with ratio Al/P=1 to hydrogen forms of aluminosilicates. Simple approximate forms based on charge and geometry coordinates are proposed for the high‐order moments of each atom, which are further required for the calculation of the electrostatic field within the structures. The possibility to use this analytical approach to evaluate the electrostatic field within embedded cluster models is also shortly discussed. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 83: 70–85, 2001     

Low cost prediction of relative stabilities of hydrogen bonded complexes from atomic multipole moments for overly short intermolecular distances

The relative stability of biologically relevant, hydrogen bonded complexes with shortened distances can be assessed at low cost by the electrostatic multipole term alone more successfully than by ab initio methods. These results imply that atomic multipole moments may help improve ligand–receptor ranking predictions, particularly in cases where accurate structural data are not available. © 2013 Wiley Periodicals, Inc.     

Torque and atomic forces for Cartesian tensor atomic multipoles with an application to crystal unit cell optimization

New equations for torque and atomic force are derived for use in flexible molecule force fields with atomic multipoles. The expressions are based on Cartesian tensors with arbitrary multipole rank. The standard method for rotating Cartesian tensor multipoles and calculating torque is to first represent the tensor with n indexes and 3ⁿ redundant components. In this work, new expressions for directly rotating the unique (n + 1)(n + 2)/2 Cartesian tensor multipole components Θ_pqr are given by introducing Cartesian tensor rotation matrix elements X( R ). A polynomial expression and a recursion relation for X( R ) are derived. For comparison, the analogous rotation matrix for spherical tensor multipoles are the Wigner functions D( R ). The expressions for X( R ) are used to derive simple equations for torque and atomic force. The torque and atomic force equations are applied to the geometry optimization of small molecule crystal unit cells. In addition, a discussion of computational efficiency as a function of increasing multipole rank is given for Cartesian tensors. © 2016 Wiley Periodicals, Inc.     

Accurate modeling of the intramolecular electrostatic energy of proteins

The (?, ψ) energy surface of blocked alanine (N-acetyl–N′-methyl alanineamide) was calculated at the Hartree-Fock (HF)/6-31G* level using ab initio molecular orbital theory. A collection of six electrostatic models was constructed, and the term electrostatic model was used to refer to (1) a set of atomic charge densities, each unable to deform with conformation; and (2) a rule for estimating the electrostatic interaction energy between a pair of atomic charge densities. In addition to two partial charge and three multipole electrostatic models, this collection includes one extremely detailed model, which we refer to as nonspherical CPK. For each of these six electrostatic models, parameters—in the form of partial charges, atomic multipoles, or generalized atomic densities—were calculated from the HF/6-31G* wave functions whose energies define the ab initio energy surface. This calculation of parameters was complicated by a problem that was found to originate from the locking in of a set of atomic charge densities, each of which contains a small polarization-induced deformation from its idealized unpolarized state. It was observed that the collective contribution of these small polarization-induced deformations to electrostatic energy differences between conformations can become large relative to ab initio energy differences between conformations. For each of the six electrostatic models, this contribution was reduced by an averaging of atomic charge densities (or electrostatic energy surfaces) over a large collection of conformations. The ab initio energy surface was used as a target with respect to which relative accuracies were determined for the six electrostatic models. A collection of 42 more complete molecular mechanics models was created by combining each of our six electrostatic models with a collection of seven models of repulsion + dispersion + intrinsic torsional energy, chosen to provide a representative sample of functional forms and parameter sets. A measure of distance was defined between model and ab initio energy surfaces; and distances were calculated for each of our 42 molecular mechanics models. For most of our 12 standard molecular mechanics models, the average error between model and ab initio energy surfaces is greater than 1.5 kcal/mol. This error is decreased by (1) careful treatment of the nonspherical nature of atomic charge densities, and (2) accurate representation of electrostatic interaction energies of types 1—2 and 1—3. This result suggests an electrostatic origin for at least part of the error between standard model and ab initio energy surfaces. Given the range of functional forms that is used by the current generation of protein potential functions, these errors cannot be corrected by compensating for errors in other energy components. © 1995 by John Wiley & Sons, Inc.     

Prediction of conformationally dependent atomic multipole moments in carbohydrates

The conformational flexibility of carbohydrates is challenging within the field of computational chemistry. This flexibility causes the electron density to change, which leads to fluctuating atomic multipole moments. Quantum Chemical Topology (QCT) allows for the partitioning of an “atom in a molecule,” thus localizing electron density to finite atomic domains, which permits the unambiguous evaluation of atomic multipole moments. By selecting an ensemble of physically realistic conformers of a chemical system, one evaluates the various multipole moments at defined points in configuration space. The subsequent implementation of the machine learning method kriging delivers the evaluation of an analytical function, which smoothly interpolates between these points. This allows for the prediction of atomic multipole moments at new points in conformational space, not trained for but within prediction range. In this work, we demonstrate that the carbohydrates erythrose and threose are amenable to the above methodology. We investigate how kriging models respond when the training ensemble incorporating multiple energy minima and their environment in conformational space. Additionally, we evaluate the gains in predictive capacity of our models as the size of the training ensemble increases. We believe this approach to be entirely novel within the field of carbohydrates. For a modest training set size of 600, more than 90% of the external test configurations have an error in the total (predicted) electrostatic energy (relative to ab initio) of maximum 1 kJ mol^?1 for open chains and just over 90% an error of maximum 4 kJ mol^?1 for rings. © 2015 Wiley Periodicals, Inc.     

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.     

Interaction of neutral and zwitterionic glycine with Zn2+ in gas phase: ab initio and SIBFA molecular mechanics calculations

The interaction of Zn²⁺ with glycine (Gly) in the gas phase is studied by a combination of ab initio and molecular mechanics techniques. The structures and energetics of the various isomers of the Gly–Zn²⁺ complex are first established via high‐level ab initio calculations. Two low‐energy isomers are characterized: one in which the metal ion interacts with the carboxylate end of zwitterionic glycine, and another in which it chelates the amino nitrogen and the carbonyl oygen of neutral glycine. These calculations lead to the first accurate value of the gas‐phase affinity of glycine for Zn²⁺. Ab initio calculations were also used to evaluate the performance of various implementations of the SIBFA force field. To assess the extent of transferability of the distributed multipoles and polarizabilities used in the SIBFA computations, two approaches are followed. In the first, approach (a), these quantities are extracted from the ab initio Hartree–Fock wave functions of glycine or its zwitterion in its entirety, and for each individual Zn²⁺‐binding conformation. In the second, approach (b), they are assembled from the appropriate constitutive fragments, namely methylamine and formic acid for neutral glycine, and protonated methylamine and formate for the zwitterion; they undergo the appropriate vector or matrix rotation to be assembled in the conformation studied. The values of the Zn²⁺–glycine interaction energies are compared to those resulting from ab initio SCF and MP2 computations using both the all‐electron 6‐311+G(2d,2p) basis set and an effective core potential together with the valence CEP 4‐31G(2d) basis set. Approach (a) values closely reproduce the ab initio ones, both in terms of the total interaction energies and of the individual components. Approach (b) can provide a similar match to ab initio interaction energies as does approach (a), provided that the two constitutive Gly building blocks are considered as separate entities having mutual interactions that are computed simultaneously with those occurring with Zn²⁺. Thus, the supermolecule is treated as a three‐body rather than a two‐body system. These results indicate that the current implementation of the SIBFA force field should be adequate to undertake accurate studies on zinc metallopeptides. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 963–973, 2000     

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     

Force related atomic multipoles in planar molecules. Derivation of atomic quadrupole and octupole moments

Atomic multipoles as defined by current methods generally do not account for forces in molecules that arise from external electrostatic fields. It is pointed out that such forces and the electrostatic potential that the molecule itself generates are both determined by the molecular multipolar tensors. The latter constitute therefore the fundamental molecular constants that determine the molecular electrostatics apart from polarization. In general the multipolar tensors include contributions from the atomic multipoles and their fluxes. In planar molecules, however, the perpendicular charge flux is zero by symmetry. This gives rise to a (previously introduced) formalism that extracts analytical, force-related, atomic multipoles from the molecular multipolar tensors. This formalism is extended in this work to include force-related (FR) atomic quadrupoles and octupoles in planar molecules. The properties of the FR atomic multipoles, including their perpendicular fluxes, are discussed and some formal theoretical and computational advantages that characterize them are indicated. As an example, the electrostatics of OCS, including the molecular electrostatic potential and the forces on the nuclei due to an external point charge, is discussed.     

Atomic multipoles and perpendicular electrostatic forces in diatomic and planar molecules

Current methods for assigning atomic multipoles focus on reproduction of the molecular electrostatic potential. Another aspect of electrostatic interaction, which is usually overlooked, is the forces that an external electric field exerts on the nuclei of a molecule. In a self-consistent theory, both the electrostatic potential and force should be accounted for. However, in general it is not easy to meet this requirement for the force. For planar molecules, though, a formal solution is available in terms of atomic multipoles that are extracted from the molecular multipolar tensors. These Force-Related (FR) atomic multipoles are discussed in detail for some typical diatomics and planar polyatomics, and are shown to provide a solid uniform framework for treating both aspects of the electrostatics. In contrast, the commonly used potential-derived charges (i.e., the atomic charges obtained by fitting the electrostatic potential) can yield large deviations with respect to electrostatic forces on the nuclei, even when the electrostatic potential is very well reproduced.     

A multipole approximation of the electrostatic potential of molecules

The problem of approximating three-dimensional spatial distributions of quantum-mechanical electrostatic potentials of molecules by analytic potentials on the basis of atomic charges, real dipoles, and atomic multipoles up to quadrupoles inclusive was considered. Real dipole potentials are created by pairs of point charges of opposite signs, and the search for their arrangement in the volume of a molecule is part of the approximation problem. A FitMEP program was developed for the optimization of the parameters of models of the types specified taking into account molecular symmetry. It was shown for the example of several molecules (HF, CO, H₂O, NH₃, CH₄, formaldehyde, methanol, formamide, ethane, cyclopropane, cyclobutane, cyclohexane, tetrahedrane, cubane, adamantane, ethylene, and benzene) that the real dipole and atomic multipole models gave errors in approximated quantum-mechanical electrostatic potential values smaller by one or two orders of magnitude compared with the atomic charge model. The atomic charge model was shown to be virtually inoperative as applied to saturated hydrocarbons. Real dipole models were slightly inferior to atomic multipole models in quality but had all the advantages of the potential of point charges as concerned simplicity and compactness, and their use in potential energy calculations did not require changes in the existing program codes.