Can the Lennard–Jones 6‐12 function replace the 10‐12 form in molecular mechanics calculations?

A protocol to replace “10‐12” hydrogen bonding function with the “6‐12” form to reproduce hydrogen bond distances, energies, and geometries in molecular mechanics calculations is described. The 6‐12 function was least‐squares fit to the normally employed 10‐12 form of the function for the hydrogen bond types of the Weiner et al. force field by iterating over the A and B coefficients. A weighting function was used to fit the curves in the most critical areas. The 6‐12 hydrogen bond model was compared with the Weiner et al. force field, OPLS/AMBER fore field, and quantum mechanical calculations on two simple systems, the water dimer and the chloride‐water interaction. The 6‐12 model produced structures, energies, and geometries that were consistent with the other molecular mechanics calculations and showed reasonable agreement to the quantum mechanical results for the water dimer. The 6‐12 model was also compared with normal calculations using a 10‐12 model on several representative systems. The results indicate that the 6‐12 function, when substituted by the procedure outlined in this work, yields structures and hydrogen bond properties that are similar to the normal 10‐12 model.

[1]  William L. Jorgensen,et al.  Energy component analysis for dilute aqueous solutions of lithium(1+), sodium(1+), fluoride(1-), and chloride(1-) ions , 1984 .

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

[3]  T. R. Dyke,et al.  Partially deuterated water dimers: Microwave spectra and structure , 1980 .

[4]  H. Berendsen,et al.  ALGORITHMS FOR MACROMOLECULAR DYNAMICS AND CONSTRAINT DYNAMICS , 1977 .

[5]  Angelo Vedani,et al.  YETI: An interactive molecular mechanics program for small‐molecule protein complexes , 1988 .

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

[7]  Hideaki Umeyama,et al.  The origin of hydrogen bonding. An energy decomposition study , 1977 .

[8]  S. Lifson,et al.  Energy functions for peptides and proteins. I. Derivation of a consistent force field including the hydrogen bond from amide crystals. , 1974, Journal of the American Chemical Society.

[9]  U. Singh,et al.  Ab initio calculations on the structure and nature of the hydrogen bonded complex H2S ⋅⋅⋅ HF , 1984 .

[10]  E. Clementi,et al.  Study of the structure of molecular complexes. V. Heat of formation for the Li+, Na+, K+, F−, and Cl− ion complexes with a single water molecule , 1973 .

[11]  J. Dannenberg.  An AM1 and ab initio molecular orbital study of water dimer , 1988 .

[12]  H. Schaefer,et al.  Extensive theoretical studies of the hydrogen‐bonded complexes (H2O)2, (H2O)2H+, (HF)2, (HF)2H+, F2H−, and (NH3)2 , 1986 .

[13]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[14]  Peter A. Kollman,et al.  Molecular mechanics studies of enzyme-substrate interactions: the interaction of L- and D-N-acetyltryptophanamide with α-chymotrypsin , 1983 .

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

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

[17]  H. Scheraga,et al.  Structure, energetics, and dynamics of small water clusters , 1975 .

[18]  Peter A. Kollman,et al.  AMBER: Assisted model building with energy refinement. A general program for modeling molecules and their interactions , 1981 .

[19]  M. Whitlow,et al.  An empirical examination of potential energy minimization using the well-determined structure of the protein crambin , 1986 .

[20]  John S. Muenter,et al.  Microwave spectrum and structure of hydrogen bonded water dimer , 1974 .

[21]  M. Frisch,et al.  Molecular orbital study of the dimers (AHn)2 formed from ammonia, water, hydrogen fluoride, phosphine, hydrogen sulfide, and hydrogen chloride , 1985 .

[22]  H. Scheraga,et al.  Energy parameters in polypeptides. 9. Updating of geometrical parameters, nonbonded interactions, and hydrogen bond interactions for the naturally occurring amino acids , 1983 .

[23]  U. Singh,et al.  A water dimer potential based on ab initio calculations using Morokuma component analyses , 1985 .

[24]  Peter A. Kollman,et al.  Molecular mechanics simulation of protein-ligand interactions: binding of thyroid hormone analogs to prealbumin , 1982 .

[25]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[26]  Peter A. Kollman,et al.  A general analysis of noncovalent intermolecular interactions , 1977 .

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

[28]  W. L. Jorgensen.  Quantum and statistical mechanical studies of liquids. 10. Transferable intermolecular potential functions for water, alcohols, and ethers. Application to liquid water , 2002 .