Accuracy and tractability of a kriging model of intramolecular polarizable multipolar electrostatics and its application to histidine

We propose a generic method to model polarization in the context of high‐rank multipolar electrostatics. This method involves the machine learning technique kriging, here used to capture the response of an atomic multipole moment of a given atom to a change in the positions of the atoms surrounding this atom. The atoms are malleable boxes with sharp boundaries, they do not overlap and exhaust space. The method is applied to histidine where it is able to predict atomic multipole moments (up to hexadecapole) for unseen configurations, after training on 600 geometries distorted using normal modes of each of its 24 local energy minima at B3LYP/apc‐1 level. The quality of the predictions is assessed by calculating the Coulomb energy between an atom for which the moments have been predicted and the surrounding atoms (having exact moments). Only interactions between atoms separated by three or more bonds (“1, 4 and higher” interactions) are included in this energy error. This energy is compared with that of a central atom with exact multipole moments interacting with the same environment. The resulting energy discrepancies are summed for 328 atom–atom interactions, for each of the 29 atoms of histidine being a central atom in turn. For 80% of the 539 test configurations (outside the training set), this summed energy deviates by less than 1 kcal mol−1. © 2013 Wiley Periodicals, Inc.

[1]  Laurent Joubert,et al.  Convergence of the electrostatic interaction based on topological atoms , 2001 .

[2]  Sarah L. Price,et al.  A TRANSFERABLE DISTRIBUTED MULTIPOLE MODEL FOR THE ELECTROSTATIC INTERACTIONS OF PEPTIDES AND AMIDES , 1990 .

[3]  Bernd Engels,et al.  Accurate Intermolecular Potentials with Physically Grounded Electrostatics. , 2011, Journal of chemical theory and computation.

[4]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole-based Molecular Mechanics for Organic Molecules. , 2011, Journal of chemical theory and computation.

[5]  Paul L A Popelier,et al.  Aqueous imidazole solutions: a structural perspective from simulations with high-rank electrostatic multipole moments. , 2011, The journal of physical chemistry. B.

[6]  Donald R. Jones,et al.  Efficient Global Optimization of Expensive Black-Box Functions , 1998, J. Glob. Optim..

[7]  Noel A. C. Cressie,et al.  Statistics for Spatial Data: Cressie/Statistics , 1993 .

[8]  P. Popelier,et al.  A convergent multipole expansion for 1,3 and 1,4 Coulomb interactions. , 2006, The Journal of chemical physics.

[9]  Benoît Roux,et al.  Modeling induced polarization with classical Drude oscillators: Theory and molecular dynamics simulation algorithm , 2003 .

[10]  P. Popelier,et al.  Intramolecular polarisable multipolar electrostatics from the machine learning method Kriging , 2011 .

[11]  Hans Ågren,et al.  Recipe of Polarized One-Electron Potential Optimization for Development of Polarizable Force Fields. , 2007, Journal of chemical theory and computation.

[12]  A. Forrester,et al.  Design and analysis of 'noisy' computer experiments , 2006 .

[13]  Paul Bratley,et al.  Algorithm 659: Implementing Sobol's quasirandom sequence generator , 1988, TOMS.

[14]  Paul L. A. Popelier,et al.  Simulation of liquid water using a high-rank quantum topological electrostatic potential , 2004 .

[15]  Teodoro Laino,et al.  Semiempirical self-consistent polarization description of bulk water, the liquid-vapor interface, and cubic ice. , 2011, The journal of physical chemistry. A.

[16]  Frances Y. Kuo,et al.  Constructing Sobol Sequences with Better Two-Dimensional Projections , 2008, SIAM J. Sci. Comput..

[17]  Mihaly Mezei,et al.  Results from an Early Polarization Model Based on Maxwell's Invariant Multipole Form. , 2007, Journal of chemical theory and computation.

[18]  Frederick E. Petry,et al.  Principles and Applications , 1997 .

[19]  A. Warshel,et al.  What are the dielectric “constants” of proteins and how to validate electrostatic models? , 2001, Proteins.

[20]  Peter J. Winn,et al.  Towards improved force fields: III. Polarization through modified atomic charges , 1999, Journal of computational chemistry.

[21]  Steven J. Stuart,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994 .

[22]  Nohad Gresh,et al.  Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy. , 2007, Journal of chemical theory and computation.

[23]  Klaus Schulten,et al.  Polarizable intermolecular potentials for water and benzene interacting with halide and metal ions. , 2009, Journal of chemical theory and computation.

[24]  Paul L. A. Popelier,et al.  High-rank quantum topological electrostatic potential: Molecular dynamics simulation of liquid hydrogen fluoride , 2003 .

[25]  J. Behler,et al.  Construction of high-dimensional neural network potentials using environment-dependent atom pairs. , 2012, The Journal of chemical physics.

[26]  G. Matheron Principles of geostatistics , 1963 .

[27]  Kai Zhu,et al.  Improved Methods for Side Chain and Loop Predictions via the Protein Local Optimization Program:  Variable Dielectric Model for Implicitly Improving the Treatment of Polarization Effects. , 2007, Journal of chemical theory and computation.

[28]  P. Popelier,et al.  Atomic properties of selected biomolecules: quantum topological atom types of carbon occurring in natural amino acids and derived molecules. , 2003, Journal of the American Chemical Society.

[29]  Mark S. Gordon,et al.  Evaluation of Charge Penetration Between Distributed Multipolar Expansions , 2000 .

[30]  Alexander D. MacKerell,et al.  Evaluating CHARMM Additive and Drude Polarizable Force Fields Through QM/MM Computations of the Vibrational Stark Effect , 2011 .

[31]  Samuel Krimm,et al.  Group Polarizability Model for Molecular Mechanics Energy Functions. , 2007, Journal of chemical theory and computation.

[32]  William L Jorgensen,et al.  Special Issue on Polarization. , 2007, Journal of chemical theory and computation.

[33]  F. Escudero,et al.  Atoms in molecules , 1982 .

[34]  David J. Wales,et al.  Topography of potential-energy surfaces for Van der Waals complexes , 1994 .

[35]  Jianshu Cao,et al.  Theory and simulation of polar and nonpolar polarizable fluids , 1993 .

[36]  C. Cramer,et al.  Self-Consistent Reaction Field Model for Aqueous and Nonaqueous Solutions Based on Accurate Polarized Partial Charges. , 2007, Journal of chemical theory and computation.

[37]  Paul L. A. Popelier,et al.  Polarisable multipolar electrostatics from the machine learning method Kriging: an application to alanine , 2012, Theoretical Chemistry Accounts.

[38]  Frank Jensen,et al.  Polarization consistent basis sets. III. The importance of diffuse functions , 2002 .

[39]  Paul L. A. Popelier,et al.  Geometrically faithful homeomorphisms between the electron density and the bare nuclear potential , 2009 .

[40]  Paul L. A. Popelier,et al.  Atomic properties of selected biomolecules. Part 1. The interpretation of atomic integration errors , 2000 .

[41]  W. Jonathan Lederer,et al.  STIM1 in Rat Ventricular Myocytes , 2011 .

[42]  Uwe Koch,et al.  Conformational dependence of atomic multipole moments , 1995 .

[43]  Paul L. A. Popelier,et al.  Atom–atom partitioning of intramolecular and intermolecular Coulomb energy , 2001 .

[44]  Frances Y. Kuo,et al.  Remark on algorithm 659: Implementing Sobol's quasirandom sequence generator , 2003, TOMS.

[45]  P. Popelier,et al.  Atomic properties of amino acids: computed atom types as a guide for future force-field design. , 2003, Chemphyschem : a European journal of chemical physics and physical chemistry.

[46]  Sarah L Price,et al.  Modelling organic crystal structures using distributed multipole and polarizability-based model intermolecular potentials. , 2010, Physical chemistry chemical physics : PCCP.

[47]  Anna-Pitschna E. Kunz,et al.  A simple, efficient polarizable molecular model for liquid carbon tetrachloride , 2011 .

[48]  Paul L A Popelier,et al.  Properties and 3D Structure of Liquid Water:  A Perspective from a High-Rank Multipolar Electrostatic Potential. , 2008, Journal of chemical theory and computation.

[49]  Paul L. A. Popelier,et al.  FORMULAS FOR THE 1ST AND 2ND DERIVATIVES OF ANISOTROPIC POTENTIALS WITH RESPECT TO GEOMETRICAL PARAMETERS , 1994 .

[50]  Glenn I. Hawe,et al.  A water potential based on multipole moments trained by machine learning—Reducing maximum energy errors , 2010 .

[51]  Jonas Boström,et al.  A NEMO potential that includes the dipole–quadrupole and quadrupole–quadrupole polarizability , 2010, J. Comput. Chem..

[52]  P. Popelier,et al.  Dynamically Polarizable Water Potential Based on Multipole Moments Trained by Machine Learning. , 2009, Journal of chemical theory and computation.

[53]  Douglas B Kell,et al.  Optimal construction of a fast and accurate polarisable water potential based on multipole moments trained by machine learning. , 2009, Physical chemistry chemical physics : PCCP.

[54]  S. Price,et al.  Computed crystal energy landscapes for understanding and predicting organic crystal structures and polymorphism. , 2009, Accounts of chemical research.

[55]  Nohad Gresh,et al.  Many‐body exchange‐repulsion in polarizable molecular mechanics. I. orbital‐based approximations and applications to hydrated metal cation complexes , 2011, J. Comput. Chem..

[56]  Laurent Joubert,et al.  Improved convergence of the ‘atoms in molecules’ multipole expansion of electrostatic interaction , 2002 .

[57]  Donald R. Jones,et al.  A Taxonomy of Global Optimization Methods Based on Response Surfaces , 2001, J. Glob. Optim..

[58]  J. R. Carl,et al.  Atom dipole interaction model for molecular polarizability. Application to polyatomic molecules and determination of atom polarizabilities , 1972 .

[59]  R. Eberhart,et al.  Comparing inertia weights and constriction factors in particle swarm optimization , 2000, Proceedings of the 2000 Congress on Evolutionary Computation. CEC00 (Cat. No.00TH8512).

[60]  Anthony J. Stone,et al.  Distributed multipole analysis, or how to describe a molecular charge distribution , 1981 .

[61]  P. Popelier,et al.  Convergence of the multipole expansion for 1,2 Coulomb interactions: The modified multipole shifting algorithm. , 2010, The Journal of chemical physics.

[62]  PENG ZHANG,et al.  Dipole preserving and polarization consistent charges , 2011, J. Comput. Chem..

[63]  Paul L A Popelier,et al.  Beyond Point Charges: Dynamic Polarization from Neural Net Predicted Multipole Moments. , 2008, Journal of chemical theory and computation.

[64]  James Kennedy,et al.  Particle swarm optimization , 1995, Proceedings of ICNN'95 - International Conference on Neural Networks.

[65]  Christophe Chipot,et al.  Derivation of Distributed Models of Atomic Polarizability for Molecular Simulations. , 2007, Journal of chemical theory and computation.

[66]  R. T. Sanderson,et al.  An Interpretation of Bond Lengths and a Classification of Bonds. , 1951, Science.