J. W. Caldwell, P. Cieplak, P.A. Kollman
Department of Pharmaceutical Chemistry, University of California San Francisco


P. Cieplak, J. W. Caldwell, P. A. Kollman
J. Comp. Chem. 2001,22,1048
Molecular mechanical models for organic and biological systems going beyond the atom centered two body additive approximation: Aqueous solution free energies of methanol and N-methyl acetamide, nucleic acid base, and amide hydrogen bonding and chloroform/water partition coefficients of the nucleic acid bases
We have developed a methodology to derive RESP charges for molecular mechanical models that include "lone pairs" on lone-pair donor sites and atom-centered polarizabilities. This approach uses a very high level ab initio cc-pVTZ basis set,(1) where the multipole moments of the molecules are as accurate as possible. The partial charges are derived self-consistently so that the model, whose electrostatic potential comes from both partial charges and induced dipoles, reproduces the quantum mechanical electrostatic potential. We then study the ability of such models to reproduce the aqueous solvation free energy of methanol and N-methyl acetamide (NMA), the base pair hydrogen bonding of the 26 base pairs analyzed by Hobza et al, and the chloroform/water partition coefficients of the five N-methyl substituted nucleic acid bases. The base pair H-bond energies are described as well as the atom centered additive model, after modifying the van der Waals parameter on the N-H to give reasonable base pair H-bond distances. The experimental solvation free energies (gas --> water) of methanol and NMA are well described, and the water/CHCl3 partition coefficients are improved over the additive model, without any parameter changes. (C) 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1048-1057, 2001.


all-atom CHCL3 solvent box of 1375 molecules

[ APD Home ] [ Submit Parameters ] [ Amber Home ] [ UK Mirror ] [ Pharmacy Home ] [ Contact ]

Last modified Mon Jan 20 09:00:39 GMT 2003