NEEMP: software for validation, accurate calculation and fast parameterization of EEM charges

BackgroundThe concept of partial atomic charges was first applied in physical and organic chemistry and was later also adopted in computational chemistry, bioinformatics and chemoinformatics. The electronegativity equalization method (EEM) is the most frequently used approach for calculating partial atomic charges. EEM is fast and its accuracy is comparable to the quantum mechanical charge calculation method for which it was parameterized. Several EEM parameter sets for various types of molecules and QM charge calculation approaches have been published and new ones are still needed and produced. Methodologies for EEM parameterization have been described in a few articles, but a software tool for EEM parameterization and EEM parameter sets validation has not been available until now.ResultsWe provide the software tool NEEMP (http://ncbr.muni.cz/NEEMP), which offers three main functionalities: EEM parameterization [via linear regression (LR) and differential evolution with local minimization (DE-MIN)]; EEM parameter set validation (i.e., validation of coverage and quality) and EEM charge calculation. NEEMP functionality is shown using a parameterization and a validation case study. The parameterization case study demonstrated that LR is an appropriate approach for smaller and homogeneous datasets and DE-MIN is a suitable solution for larger and heterogeneous datasets. The validation case study showed that EEM parameter set coverage and quality can still be problematic. Therefore, it makes sense to verify the coverage and quality of EEM parameter sets before their use, and NEEMP is an appropriate tool for such verification. Moreover, it seems from both case studies that new EEM parameterizations need to be performed and new EEM parameter sets obtained with high quality and coverage for key structural databases.ConclusionWe provide the software tool NEEMP, which is to the best of our knowledge the only available software package that enables EEM parameterization and EEM parameter set validation. Additionally, its DE-MIN parameterization method is an innovative approach, developed by ourselves and first published in this work. In addition, we also prepared four high-quality EEM parameter sets tailored to ligand molecules.Graphical abstract.

[1]  R. S. Mulliken Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

[2]  Robert S. Mulliken,et al.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. II. Overlap Populations, Bond Orders, and Covalent Bond Energies , 1955 .

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

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

[5]  Frank Weinhold,et al.  Natural bond orbital analysis of near‐Hartree–Fock water dimer , 1983 .

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

[7]  F. Weinhold,et al.  Natural population analysis , 1985 .

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

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

[10]  R. Bader Atoms in molecules , 1990 .

[11]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[12]  Wilfried J. Mortier,et al.  Probing the reactivity of different sites within a molecule or solid by direct computation of molecular sensitivities via an extension of the electronegativity equalization method , 1991 .

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

[14]  J. Gasteiger,et al.  FROM ATOMS AND BONDS TO THREE-DIMENSIONAL ATOMIC COORDINATES : AUTOMATIC MODEL BUILDERS , 1993 .

[15]  Jorge Gálvez,et al.  Charge Indexes. New Topological Descriptors , 1994, J. Chem. Inf. Comput. Sci..

[16]  Robert P. Sheridan,et al.  Chemical Similarity Using Physiochemical Property Descriptors , 1996, J. Chem. Inf. Comput. Sci..

[17]  C. Lemmen,et al.  FLEXS: a method for fast flexible ligand superposition. , 1998, Journal of medicinal chemistry.

[18]  James Demmel,et al.  LAPACK Users' Guide, Third Edition , 1999, Software, Environments and Tools.

[19]  Jack Dongarra,et al.  LAPACK Users' Guide, 3rd ed. , 1999 .

[20]  D. Rognan,et al.  Protein-based virtual screening of chemical databases. 1. Evaluation of different docking/scoring combinations. , 2000, Journal of medicinal chemistry.

[21]  H. Scheraga,et al.  A Fast Method for Calculating Geometry-Dependent Net Atomic Charges for Polypeptides , 2001 .

[22]  Christopher M. Hadad,et al.  Comparison of different atomic charge schemes for predicting pKa variations in substituted anilines and phenols , 2002 .

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

[24]  Peter Gedeck,et al.  Calculation of Intersubstituent Similarity Using R-Group Descriptors , 2003, J. Chem. Inf. Comput. Sci..

[25]  Haruki Nakamura,et al.  Announcing the worldwide Protein Data Bank , 2003, Nature Structural Biology.

[26]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

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

[28]  Zukang Feng,et al.  Ligand Depot: a data warehouse for ligands bound to macromolecules , 2004, Bioinform..

[29]  Antti Poso,et al.  BRUTUS: optimization of a grid-based similarity function for rigid-body molecular superposition. 1. Alignment and virtual screening applications. , 2005, Journal of medicinal chemistry.

[30]  Milan Hodoscek,et al.  Development and Validation of Empirical Force Field Parameters for Netropsin , 2005, J. Chem. Inf. Model..

[31]  Vladimir A. Palyulin,et al.  Atomic charges via electronegativity equalization: Generalizations and perspectives , 2006 .

[32]  Radka Svobodová Vareková,et al.  Optimized and parallelized implementation of the electronegativity equalization method and the atom‐bond electronegativity equalization method , 2006, J. Comput. Chem..

[33]  Hwangseo Park,et al.  Critical assessment of the automated AutoDock as a new docking tool for virtual screening , 2006, Proteins.

[34]  R. A. Nistor,et al.  A generalization of the charge equilibration method for nonmetallic materials. , 2006, The Journal of chemical physics.

[35]  M. Powell The NEWUOA software for unconstrained optimization without derivatives , 2006 .

[36]  C. Henze,et al.  Fleshing‐out Pharmacophores with Volume Rendering of the Laplacian of the Charge Density and Hyperwall Visualization Technology , 2007 .

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

[38]  Simon Suchomel,et al.  Electronegativity Equalization Method: Parameterization and Validation for Large Sets of Organic, Organohalogene and Organometal Molecule , 2007, International Journal of Molecular Sciences.

[39]  Peter Ertl,et al.  Estimation of pKa for Druglike Compounds Using Semiempirical and Information-Based Descriptors , 2007, J. Chem. Inf. Model..

[40]  D. Mathieu Split charge equilibration method with correct dissociation limits. , 2007, The Journal of chemical physics.

[41]  Yanli Wang,et al.  PubChem: Integrated Platform of Small Molecules and Biological Activities , 2008 .

[42]  A. V. van Duin,et al.  ReaxFF reactive force field for molecular dynamics simulations of hydrocarbon oxidation. , 2008, The journal of physical chemistry. A.

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

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

[45]  Thanyada Rungrotmongkol,et al.  Dynamic Behavior of Avian Influenza A Virus Neuraminidase Subtype H5N1 in Complex with Oseltamivir, Zanamivir, Peramivir, and Their Phosphonate Analogues , 2009, J. Chem. Inf. Model..

[46]  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.

[47]  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..

[48]  R. Bader,et al.  A quantum theory of molecular structure and its applications , 1991 .

[49]  Chris Morley,et al.  Open Babel: An open chemical toolbox , 2011, J. Cheminformatics.

[50]  D. Stalke Meaningful structural descriptors from charge density. , 2011, Chemistry.

[51]  Radka Svobodová Vareková,et al.  Predicting pKa Values of Substituted Phenols from Atomic Charges: Comparison of Different Quantum Mechanical Methods and Charge Distribution Schemes , 2011, J. Chem. Inf. Model..

[52]  Radka Svobodová Vareková,et al.  Predicting pKa values from EEM atomic charges , 2013, Journal of Cheminformatics.

[53]  Radka Svobodová Vareková,et al.  How Does the Methodology of 3D Structure Preparation Influence the Quality of pKa Prediction? , 2015, J. Chem. Inf. Model..

[54]  Radka Svobodová Vareková,et al.  High-quality and universal empirical atomic charges for chemoinformatics applications , 2015, Journal of Cheminformatics.

[55]  Silvio C. E. Tosatto,et al.  Tools and data services registry: a community effort to document bioinformatics resources , 2015, Nucleic Acids Res..