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

[1]  J. Gasteiger,et al.  ITERATIVE PARTIAL EQUALIZATION OF ORBITAL ELECTRONEGATIVITY – A RAPID ACCESS TO ATOMIC CHARGES , 1980 .

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

[3]  P. Geerlings,et al.  Conceptual density functional theory. , 2003, Chemical reviews.

[4]  Ulf Ryde,et al.  Comparison of methods for deriving atomic charges from the electrostatic potential and moments , 1998 .

[5]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[6]  David M. Gange,et al.  Charges fit to electrostatic potentials. II. Can atomic charges be unambiguously fit to electrostatic potentials? , 1996 .

[7]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[8]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[9]  Chang-Sheng Wang,et al.  Atom−Bond Electronegativity Equalization Method. 1. Calculation of the Charge Distribution in Large Molecules , 1997 .

[10]  Patrick Bultinck,et al.  Algebraic relationships between conceptual DFT quantities and the electronegativity equalization hardness matrix , 2002 .

[11]  Michel Waroquier,et al.  The Electronegativity Equalization Method I: Parametrization and Validation for Atomic Charge Calculations , 2002 .

[12]  Kalyanmoy Deb,et al.  A fast and elitist multiobjective genetic algorithm: NSGA-II , 2002, IEEE Trans. Evol. Comput..

[13]  Yizeng Liang,et al.  A modified electronegativity equalization method for fast and accurate calculation of atomic charges in large biological molecules. , 2009, Physical chemistry chemical physics : PCCP.

[14]  Ralph G. Pearson,et al.  HARD AND SOFT ACIDS AND BASES , 1963 .

[15]  Y. Mo,et al.  Energy decomposition analysis of intermolecular interactions using a block-localized wave function approach , 2000 .

[16]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[17]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[18]  Patrick Bultinck,et al.  High-speed calculation of AIM charges through the electronegativity equalization method , 2004 .

[19]  David S. Wishart,et al.  DrugBank: a knowledgebase for drugs, drug actions and drug targets , 2007, Nucleic Acids Res..

[20]  Mark S. Johnson,et al.  Generating Conformer Ensembles Using a Multiobjective Genetic Algorithm , 2007, J. Chem. Inf. Model..

[21]  F. L. Hirshfeld.  Bonded-atom fragments for describing molecular charge densities , 1977 .

[22]  Michael K. Gilson,et al.  Fast Assignment of Accurate Partial Atomic Charges: An Electronegativity Equalization Method that Accounts for Alternate Resonance Forms , 2003, J. Chem. Inf. Comput. Sci..

[23]  Patrick Bultinck,et al.  Toward an Alternative Hardness Kernel Matrix Structure in the Electronegativity Equalization Method (EEM) , 2006, J. Chem. Inf. Model..

[24]  Paola Gramatica,et al.  The Importance of Being Earnest: Validation is the Absolute Essential for Successful Application and Interpretation of QSPR Models , 2003 .

[25]  Robert A. Schoonheydt,et al.  INVESTIGATION OF SUPRAMOLECULAR SYSTEMS BY A COMBINATION OF THE ELECTRONEGATIVITY EQUALIZATION METHOD AND A MONTE CARLO SIMULATION TECHNIQUE , 1995 .

[26]  Patrick Bultinck,et al.  Critical analysis and extension of the Hirshfeld atoms in molecules. , 2007, The Journal of chemical physics.

[27]  Wilfried J. Mortier,et al.  Electronegativity-equalization method for the calculation of atomic charges in molecules , 1986 .

[28]  Victor Guallar,et al.  Importance of accurate charges in molecular docking: Quantum mechanical/molecular mechanical (QM/MM) approach , 2005, J. Comput. Chem..

[29]  Radka Svobodová Vareková,et al.  Software news and updates electronegativity equalization method: Parameterization and validation for organic molecules using the Merz‐Kollman‐Singh charge distribution scheme , 2009, J. Comput. Chem..

[30]  T. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[31]  Patrick Bultinck,et al.  The Electronegativity Equalization Method II: Applicability of Different Atomic Charge Schemes , 2002 .

[32]  R. Parr,et al.  Absolute hardness: companion parameter to absolute electronegativity , 1983 .

[33]  S. A. Pisarev,et al.  Fast tools for calculation of atomic charges well suited for drug design , 2008, SAR and QSAR in environmental research.

[34]  Alexander D. MacKerell.  Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..

[35]  Alfredo Mayall Simas,et al.  RM1: A reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I , 2006, J. Comput. Chem..

[36]  Johann Gasteiger,et al.  A new model for calculating atomic charges in molecules , 1978 .

[37]  David Weininger,et al.  SMILES. 2. Algorithm for generation of unique SMILES notation , 1989, J. Chem. Inf. Comput. Sci..

[38]  Eszter Czinki,et al.  A combined electronegativity equalization and electrostatic potential fit method for the determination of atomic point charges , 2007, J. Comput. Chem..

[39]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[40]  Binne Zwanenburg,et al.  A two-step chirality transfer from (−)-endo- to (−)-exo-tricyclo[5.2.1.026]deca-4,8-dien-3-one , 1993 .

[41]  G. Scuseria,et al.  Gaussian 03, Revision E.01. , 2007 .

[42]  Jianfen Fan,et al.  Extending and simplifying the electronegativity equalization method , 1998 .

[43]  David Weininger,et al.  SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules , 1988, J. Chem. Inf. Comput. Sci..

[44]  P. Kollman,et al.  An approach to computing electrostatic charges for molecules , 1984 .

[45]  Patrick Bultinck,et al.  Electrostatic Potentials from Self-Consistent Hirshfeld Atomic Charges. , 2009, Journal of chemical theory and computation.

[46]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[47]  Kenneth Levenberg.  A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[48]  Peter Politzer,et al.  The fundamental nature and role of the electrostatic potential in atoms and molecules , 2002 .

[49]  Richard F. W. Bader.  A quantum theory of molecular structure and its applications , 1991 .

[50]  Konstantin S. Smirnov,et al.  Consistent implementation of the electronegativity equalization method in molecular mechanics and molecular dynamics , 1996 .