Improved Formulas for the Calculation of the Electrostatic Contribution to the Intermolecular Interaction Energy from Multipolar Expansion of the Electronic Distribution.

We have, within the framework of the molecular mechanics method SIBFA, improved the formulation of the Coulomb (electrostatic) energy contribution to the intermolecular interaction energy. This was done by integrating "overlap-like" terms into two components of the multipolar development used to calculate this contribution in SIBFA. The calibration of the new component is done on five water dimers by fitting this augmented electrostatic contribution to the corresponding Ec term. Several tests are done on (i) representative neutral and ionic hydrogen-bonded complexes; (ii) the complexes of metal cations (Cu(I) and Cu(II)) with a neutral or an anionic ligand; and (iii) a representative stacked complex. The improvement brought by the new formulation reduces the difference between the ab initio (Ec) and molecular mechanics (EMTP*) values by almost an order of magnitude when compared to the values of EMTP calculated using the standard method.

[1]  Kim Palmö,et al.  A polarizable electrostatic model of the N‐methylacetamide dimer , 2001, J. Comput. Chem..

[2]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[3]  E. M. Evleth,et al.  Critical analysis of electric field modeling: Formamide , 1992 .

[4]  Paul S. Bagus,et al.  A new analysis of charge transfer and polarization for ligand–metal bonding: Model studies of Al4CO and Al4NH3 , 1984 .

[5]  Peter Pulay,et al.  Efficient elimination of basis set superposition errors by the local correlation method: Accurate ab initio studies of the water dimer , 1993 .

[6]  A.-L. Derepas,et al.  Can we understand the different coordinations and structures of closed‐shell metal cation‐water clusters? , 2002, J. Comput. Chem..

[7]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[8]  Dennis R. Salahub,et al.  Optimization of Gaussian-type basis sets for local spin density functional calculations. Part I. Boron through neon, optimization technique and validation , 1992 .

[9]  P. Kollman,et al.  Advancing beyond the atom‐centered model in additive and nonadditive molecular mechanics , 1997 .

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

[11]  Claude Millot,et al.  Revised Anisotropic Site Potentials for the Water Dimer and Calculated Properties , 1998 .

[12]  Mark S. Gordon,et al.  The Effective Fragment Potential Method: A QM-Based MM Approach to Modeling Environmental Effects in Chemistry , 2001 .

[13]  Peter A. Kollman,et al.  Application of the multimolecule and multiconformational RESP methodology to biopolymers: Charge derivation for DNA, RNA, and proteins , 1995, J. Comput. Chem..

[14]  Mark S. Gordon,et al.  An effective fragment method for modeling solvent effects in quantum mechanical calculations , 1996 .

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

[16]  P. Kollman,et al.  An SCF partitioning scheme for the hydrogen bond , 1970 .

[17]  Donald E. Williams Representation of the molecular electrostatic potential by atomic multipole and bond dipole models , 1988 .

[18]  Christophe Chipot,et al.  Modeling amino acid side chains. 2. Determination of point charges from electrostatic properties: toward transferable point charge models , 1993 .

[19]  M. Dreyfus,et al.  A non-empirical study of the hydrogen bond between peptide units , 1970 .

[20]  P. Claverie,et al.  The exact multicenter multipolar part of a molecular charge distribution and its simplified representations , 1988 .

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

[22]  Nohad Gresh,et al.  Critical Role of Anisotropy for the Dimerization Energies of Two Protein−Protein Recognition Motifs: cis-N-Methylacetamide versus a β-Sheet Conformer of Alanine Dipeptide. A Joint ab Initio, Density Functional Theory, and Molecular Mechanics Investigation , 1999 .

[23]  Junmei Wang,et al.  Automatic parameterization of force field by systematic search and genetic algorithms , 2001, J. Comput. Chem..

[24]  Visvaldas Kairys,et al.  Evaluation of the charge penetration energy between non-orthogonal molecular orbitals using the Spherical Gaussian Overlap approximation , 1999 .

[25]  Keiji Morokuma,et al.  Molecular Orbital Studies of Hydrogen Bonds. III. C=O···H–O Hydrogen Bond in H2CO···H2O and H2CO···2H2O , 1971 .

[26]  William H. Fink,et al.  Frozen fragment reduced variational space analysis of hydrogen bonding interactions. Application to the water dimer , 1987 .

[27]  Richard A. Friesner,et al.  Pseudospectral localized Mo/ller–Plesset methods: Theory and calculation of conformational energies , 1995 .

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

[29]  György G. Ferenczy,et al.  Transferable net atomic charges from a distributed multipole analysis for the description of electrostatic properties: a case study of saturated hydrocarbons , 1993 .

[30]  Nohad Gresh,et al.  Energetics of Zn2+ binding to a series of biologically relevant ligands: A molecular mechanics investigation grounded on ab initio self‐consistent field supermolecular computations , 1995, J. Comput. Chem..

[31]  Laurent Emmanuel Dardenne,et al.  Reassociation of fragments using multicentered multipolar expansions: peptide junction treatments to investigate electrostatic properties of proteins , 2001, J. Comput. Chem..

[32]  Ulf Berg,et al.  Inter‐ and intramolecular potential for the N‐formylglycinamide‐water system. A comparison between theoretical modeling and empirical force fields , 2003, J. Comput. Chem..

[33]  Richard A. Friesner,et al.  Accurate ab Initio Quantum Chemical Determination of the Relative Energetics of Peptide Conformations and Assessment of Empirical Force Fields , 1997 .