Optimal scaling factors for CM1 and CM3 atomic charges in RM1-based aqueous simulations

Jonah Z. Vilseck, Somisetti V. Sambasivarao, Orlando Acevedo

Research output: Contribution to journalArticlepeer-review

12 Scopus citations

Abstract

Scaling factors for atomic charges derived from the RM1 semiempirical quantum mechanical wavefunction in conjunction with CM1 and CM3 charge models have been optimized by minimizing errors in absolute free energies of hydration, ΔGhyd, for a set of 40 molecules. Monte Carlo statistical mechanics simulations and free energy perturbation theory were used to annihilate the solutes in gas and in a box of TIP4P water molecules. Lennard-Jones parameters from the optimized potentials for liquid simulations-all atom (OPLS-AA) force field were utilized for the organic compounds. Optimal charge scaling factors have been determined as 1.11 and 1.14 for the CM1R and CM3R methods, respectively, and the corresponding unsigned average errors in ΔGhyd relative to experiment were 2.05 and 1.89 kcal/mol. Computed errors in aniline and two derivatives were particularly large for RM1 and their removal from the data set lowered the overall errors to 1.61 and 1.75 kcal/mol for CM1R and CM3R. Comparisons are made to the AM1 method which yielded total errors in ΔGhyd of 1.50 and 1.64 kcal/mol for CM1A1.14 and CM3A1.15, respectively. This work is motivated by the need for a highly efficient yet accurate quantum mechanical (QM) method to study condensed-phase and enzymatic chemical reactions via mixed QM and molecular mechanical (QM/MM) simulations. As an initial test, the Menshutkin reaction between NH3 and CH3Cl in water was computed using a RM1/TIP4P-Ew/CM3R procedure and the resultant ΔG, ΔGrxn, and geometries were in reasonable accord with other computational methods; however, some potentially serious shortcomings in RM1 are discussed.

Original languageEnglish (US)
Pages (from-to)2836-2842
Number of pages7
JournalJournal of Computational Chemistry
Volume32
Issue number13
DOIs
StatePublished - Oct 1 2011
Externally publishedYes

Keywords

  • Menshutkin reaction
  • QM/MM
  • RM1
  • atomic charges
  • free energy of hydration

ASJC Scopus subject areas

  • Chemistry(all)
  • Computational Mathematics

Fingerprint Dive into the research topics of 'Optimal scaling factors for CM1 and CM3 atomic charges in RM1-based aqueous simulations'. Together they form a unique fingerprint.

Cite this