Inter‐ and intramolecular potential for the N‐formylglycinamide‐water system. A comparison between theoretical modeling and empirical force fields

An intramolecular NEMO potential is presented for the N‐formylglycinamide molecule together with an intermolecular potential for the N‐formylglycinamide‐water system. The intramolecular N‐formylglycinamide potential can be used as a building block for the backbone of polypeptides and proteins. Two intramolecular minima have been obtained. One, denoted as C5, is stabilized by a hydrogen bonded five member ring, and the other, denoted as C7, corresponds to a seven membered ring. The interaction between one water molecule and the N‐formylglycinamide system is also studied and compared with Hartree‐Fock SCF calculations and with the results obtained for some of the more commonly used force fields. The agreement between the NEMO and SCF energies for the complexes is in general superior to that of the other force fields. In the C7 region the surfaces obtained from the intramolecular part of the commonly used force fields are too flat compared to the NEMO potential and the ab initio calculations. We further analyze the possibility of using a charge distribution obtained from one conformation to describe the charge distribution of other conformations. We have found that the use of polarizabilities and generic dipoles can model most of the changes in charge density due to the different geometry of the new conformations, but that one can expect additional errors in the interaction energies that are of the order of 1 kcal/mol. © 2002 Wiley Periodicals, Inc. J Comput Chem 24: 161–176, 2003

[1]  Gunnar Karlström,et al.  The inclusion of electron correlation in intermolecular potentials: applications to the formamide dimer and liquid formamide , 2000 .

[2]  B. Thole.  Molecular polarizabilities calculated with a modified dipole interaction , 1981 .

[3]  P. Åstrand,et al.  NONEMPIRICAL INTERMOLECULAR POTENTIALS FOR UREA-WATER SYSTEMS , 1994 .

[4]  Peter A. Kollman,et al.  Calculation of the Φ-Ψ maps for alanyl and glycyl dipeptides with different additive and non-additive molecular mechanical models , 1997 .

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

[6]  G. Karlström.  Local polarizabilities in molecules, based on ab initio Hartree-Fock calculations , 1982 .

[7]  E. Cabaleiro-Lago,et al.  Development of an intermolecular potential function for interactions in formamide clusters based on ab initio calculations , 1999 .

[8]  I. H. Hillier,et al.  A molecular dynamics study of the conformation of the alanine dipeptide in aqueous solution using a quantum mechanical potential , 1997 .

[9]  Keith E. Gubbins,et al.  Theory of molecular fluids , 1984 .

[10]  A. Stone,et al.  Towards an accurate intermolecular potential for water , 1992 .

[11]  V. Barone,et al.  SOLVENT EFFECTS ON THE CONFORMATIONAL BEHAVIOR OF MODEL PEPTIDES. A COMPARISON BETWEEN DIFFERENT CONTINUUM MODELS , 1996 .

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

[13]  Ian R. Gould,et al.  A quantum Mechanical Investigation of the Conformational Energetics of the Alanine and Glycine Dipeptides in the Gas Phase and in Aqueous Solution , 1994 .

[14]  Sarah L. Price,et al.  Electrostatic models for polypeptides: can we assume transferability? , 1992 .

[15]  R. Bader,et al.  Effect of twisting a polypeptide on its geometry and electron distribution , 1994 .

[16]  Roland Lindh,et al.  The water dimer interaction energy: Convergence to the basis set limit at the correlated level , 1997 .

[17]  K. Tang,et al.  An improved simple model for the van der Waals potential based on universal damping functions for the dispersion coefficients , 1984 .

[18]  G. Chang,et al.  Macromodel—an integrated software system for modeling organic and bioorganic molecules using molecular mechanics , 1990 .

[19]  Per-Olof Widmark,et al.  Density matrix averaged atomic natural orbital (ANO) basis sets for correlated molecular wave functions , 1995 .

[20]  Ola Engkvist,et al.  Intermolecular Potential for the 1,2-Dimethoxyethane−Water Complex , 1996 .

[21]  G. Karlström,et al.  MODELING OF THE EXCHANGE REPULSION ENERGY , 1998 .

[22]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[23]  Peter A. Kollman,et al.  Ion solvation in polarizable water: molecular dynamics simulations , 1991 .

[24]  Teresa Head-Gordon,et al.  Stabilization of helices in glycine and alanine dipeptides in a reaction field model of solvent , 1994 .

[25]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[26]  F. B. van Duijneveldt,et al.  Convergence to the basis‐set limit in ab initio calculations at the correlated level on the water dimer , 1992 .

[27]  G. Karlström.  A new method for approximate estimates of the dispersion interaction between two molecules , 1980 .

[28]  Charles L. Brooks,et al.  Theoretical study of blocked glycine and alanine peptide analogs , 1991 .