Analysis of a large data base of electrostatic potential derived atomic charges

A large data base of 6‐31G*, MNDO, AM1, and PM3 electrostatic potential (ESP) derived point charges of amino acids and monosaccharides is analyzed. We find that MNDO correlates well with 6‐31G* ESP derived point charges, while AM1 and PM3 do so quite poorly. Furthermore, scaling MNDO ESP derived point charges enhances the ability of MNDO to reproduce 6‐31G* results. We used our data base to attempt to derive a 6‐31G* transferable charge model at an atom‐by‐atom level. We find that it is simple to derive a transferable model for monosaccharides, but for the amino acids statistical difficulties make this a less attractive approach. The transferable charge model for the monosaccharides is slightly better than MNDO, but scaled MNDO charges perform significantly better than the transferable model. We also carried out a QMD simulation on the alanine dipeptide to assess the fluctuations that would be expected in atomic point charges during the course of an MD simulation. Relatively large charge fluctuations are observed and their impact on molecular simulation is addressed. © 1992 by John Wiley & Sons, Inc.

[1]  J. Stewart Optimization of parameters for semiempirical methods II. Applications , 1989 .

[2]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[3]  Peter A. Kollman,et al.  Implementation of nonadditive intermolecular potentials by use of molecular dynamics: development of a water-water potential and water-ion cluster interactions , 1990 .

[4]  S H Kim,et al.  Atomic charges for DNA constituents derived from single-crystal X-ray diffraction data. , 1990, Journal of molecular biology.

[5]  J. Stewart Optimization of parameters for semiempirical methods. III Extension of PM3 to Be, Mg, Zn, Ga, Ge, As, Se, Cd, In, Sn, Sb, Te, Hg, Tl, Pb, and Bi , 1991 .

[6]  Warren J. Hehre,et al.  AB INITIO Molecular Orbital Theory , 1986 .

[7]  Michiel Sprik,et al.  A polarizable model for water using distributed charge sites , 1988 .

[8]  Kenneth B. Wiberg,et al.  Properties of atoms in molecules: Dipole moments and transferability of properties , 1987 .

[9]  M. Eisenstein SCF Deformation densities and electrostatic potentials of purines and pyrimidines , 1988 .

[10]  Peter L. Cummins,et al.  Atomic charges derived from semiempirical electrostatic potentials; an interaction energy method , 1990 .

[11]  R. J. Harrison,et al.  An ab initio distributed multipole study of the electrostatic potential around an undecapeptide cyclosporin derivative and a comparison with point charge electrostatic models , 1989 .

[12]  Walter Thiel,et al.  Ground States of Molecules. 39. MNDO Results for Molecules Containing Hydrogen, Carbon, Nitrogen, and Oxygen , 1977 .

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

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

[15]  Modesto Orozco,et al.  On the use of AM1 and MNDO wave functions to compute accurate electrostatic charges , 1990 .

[16]  J. Stewart Optimization of parameters for semiempirical methods I. Method , 1989 .

[17]  Donald E. Williams,et al.  Alanyl dipeptide potential‐derived net atomic charges and bond dipoles, and their variation with molecular conformation , 1990 .

[18]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

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

[20]  Modesto Orozco,et al.  Effect of electron correlation on the electrostatic potential distribution of molecules , 1991 .

[21]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[22]  György G. Ferenczy,et al.  Semiempirical AM1 electrostatic potentials and AM1 electrostatic potential derived charges: A comparison with ab initio values , 1989 .

[23]  Malcolm E. Davis,et al.  Electrostatics in biomolecular structure and dynamics , 1990 .

[24]  Ji-Min Yan,et al.  Point-Charge Models for Molecules Derived from Least-Squares Fitting of the Electric Potential , 1988 .

[25]  Michael J. S. Dewar,et al.  Analytical first derivatives if the energy in MNDO , 1978, Comput. Chem..

[26]  S. Harvey Treatment of electrostatic effects in macromolecular modeling , 1989, Proteins.

[27]  U. Dinur "Flexible" water molecules in external electrostatic potentials , 1990 .

[28]  S H Kim,et al.  Determinations of atomic partial charges for nucleic acid constituents from x‐ray diffraction data. I. 2′‐Deoxycytidine‐5′‐monophosphate , 1985, Biopolymers.

[29]  H. A. Levy,et al.  α-D-Glucose: Precise Determination of Crystal and Molecular Structure by Neutron-Diffraction Analysis , 1965, Science.

[30]  M. Murcko,et al.  The Response of Electrons to Structural Changes , 1991, Science.

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

[32]  M. L. Connolly Analytical molecular surface calculation , 1983 .

[33]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

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

[35]  A. Kuki,et al.  Partial charges by multipole constraint. Application to the amino acids , 1990 .

[36]  Philip Coppens,et al.  Experimental charge densities in chemistry: what is next? , 1989 .

[37]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[38]  J. Stewart,et al.  Semi‐empirical calculations of molecular trajectories: Method and applications to some simple molecular systems , 1987 .

[39]  F. Allen,et al.  The Cambridge Crystallographic Data Centre: computer-based search, retrieval, analysis and display of information , 1979 .

[40]  Donald E. Williams,et al.  Representation of the molecular electrostatic potential by a net atomic charge model , 1981 .

[41]  M. Dewar,et al.  Ground States of Molecules. 38. The MNDO Method. Approximations and Parameters , 1977 .

[42]  G. Petsko,et al.  Weakly polar interactions in proteins. , 1988, Advances in protein chemistry.

[43]  Martin J. Field,et al.  Free energy perturbation method for chemical reactions in the condensed phase: a dynamic approach based on a combined quantum and molecular mechanics potential , 1987 .

[44]  L. E. Chirlian,et al.  Atomic charges derived from electrostatic potentials: A detailed study , 1987 .

[45]  D. Hohl,et al.  Structure, bonding, and dynamics in heterocyclic sulfur-selenium molecules, SexSy , 1990 .

[46]  J. Almlöf,et al.  Principles for a direct SCF approach to LICAO–MOab‐initio calculations , 1982 .