Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation

Journal of Computational Chemistry - Tập 23 Số 16 - Trang 1623-1641 - 2002
Araz Jakalian1, D. B. Jack, Christopher I. Bayly
1Boehringer Ingelheim (Canada) Ltd Research and Development, 2100 Rue Cunard, Laval, Quebec, Canada, H7S 2G5.

Tóm tắt

AbstractWe present the first global parameterization and validation of a novel charge model, called AM1‐BCC, which quickly and efficiently generates high‐quality atomic charges for computer simulations of organic molecules in polar media. The goal of the charge model is to produce atomic charges that emulate the HF/6‐31G* electrostatic potential (ESP) of a molecule. Underlying electronic structure features, including formal charge and electron delocalization, are first captured by AM1 population charges; simple additive bond charge corrections (BCCs) are then applied to these AM1 atomic charges to produce the AM1‐BCC charges. The parameterization of BCCs was carried out by fitting to the HF/6‐31G* ESP of a training set of >2700 molecules. Most organic functional groups and their combinations were sampled, as well as an extensive variety of cyclic and fused bicyclic heteroaryl systems. The resulting BCC parameters allow the AM1‐BCC charging scheme to handle virtually all types of organic compounds listed in The Merck Index and the NCI Database. Validation of the model was done through comparisons of hydrogen‐bonded dimer energies and relative free energies of solvation using AM1‐BCC charges in conjunction with the 1994 Cornell et al. forcefield for AMBER.13 Homo‐ and hetero‐dimer hydrogen‐bond energies of a diverse set of organic molecules were reproduced to within 0.95 kcal/mol RMS deviation from the ab initio values, and for DNA dimers the energies were within 0.9 kcal/mol RMS deviation from ab initio values. The calculated relative free energies of solvation for a diverse set of monofunctional isosteres were reproduced to within 0.69 kcal/mol of experiment. In all these validation tests, AMBER with the AM1‐BCC charge model maintained a correlation coefficient above 0.96. Thus, the parameters presented here for use with the AM1‐BCC method present a fast, accurate, and robust alternative to HF/6‐31G* ESP‐fit charges for general use with the AMBER force field in computer simulations involving organic small molecules. © 2002 Wiley Periodicals, Inc. J Comput Chem 23: 1623–1641, 2002

Từ khóa


Tài liệu tham khảo

Williams D. E., 1991

Francl M. M., 2000

10.1021/ja00074a030

10.1002/(SICI)1096-987X(20000130)21:2<132::AID-JCC5>3.0.CO;2-P

10.1021/ja00299a024

(b)Stewart J. J. P.MOPAC 6.0 (SGI version) Quantum Program Chemistry Exchange (QCPE) #455. AM1 atomic charges are the standard population charges produced by default by MOPAC 6.0 and are derived from the Coulson density matrix (see pp.4–8 notes 18 and 20 of the MOPAC 6.0 manual).

10.1002/jcc.540110404

10.1002/jcc.540140706

10.1007/BF00117280

10.1021/jp972682r

10.1021/j100142a004

Frisch M. J., 1992, Gaussian 92

10.1063/1.453069

10.1107/S0021889883010985

Frisch M. J., 1995, Gaussian 94

McCammon J. A., 1989, Dynamics of Proteins and Nucleic Acids

10.1021/ja00124a002

Case D. A., 1997

10.1021/jp952760f

10.1002/jcc.540141112

http://www.ccl.net/cca/data/ff_evaluation_suite(accessed 18 November1998).

The CD‐ROM of the Merck Index 12thEdition is commercially available.

The NCI database of ca. 250 000 molecules was downloaded fromhttp://cactus.cit.nih.gov/ncidb/download.html

10.1002/(SICI)1096-987X(199710)18:13<1632::AID-JCC5>3.0.CO;2-S

10.1021/ci00015a015

10.1002/(SICI)1096-987X(199604)17:5/6<490::AID-JCC1>3.0.CO;2-P

10.1002/(SICI)1096-987X(199604)17:5/6<520::AID-JCC2>3.0.CO;2-W

10.1021/ja00085a033

10.1021/ja00105a060

10.1021/ja950302e

10.1023/A:1007902804814

10.1007/BF00126219

10.1021/ja984106u

10.1021/ja993663t

Thông tin
Thông tin xuất bản

Journal of Computational Chemistry

Tập 23 Số 16

1623-1641

Thông tin tác giả