title: Paradynamics: an effective and reliable model for ab initio QM/MM free-energy calculations and related tasks.
abstract: Recent years have seen tremendous effort in the development of approaches with which to obtain quantum mechanics/molecular mechanics (QM/MM) free energies for reactions in the condensed phase. Nevertheless, there remain significant challenges to address, particularly, the high computational cost involved in performing proper configurational sampling and, in particular, in obtaining ab initio QM/MM (QM(ai)/MM) free-energy surfaces. One increasingly popular approach that seems to offer an ideal way to progress in this direction is the elegant metadynamics (MTD) approach. However, in the current work, we point out the subtle efficiency problems associated with this approach and illustrate that we have at hand what is arguably a more powerful approach. More specifically, we demonstrate the effectiveness of an updated version of our original idea of using a classical reference potential for QM(ai)/MM calculations [J. Phys. Chem. 1995, 99, 17516)], which we refer to as paradynamics (PD). This approach is based on the use of an empirical valence bond (EVB) reference potential, which is already similar to the real ab initio potential. The reference potential is fitted to the ab initio potential by an iterative and, to a great degree, automated refinement procedure. The corresponding free-energy profile is then constructed using the refined EVB potential, and the linear response approximation (LRA) is used to evaluate the QM(ai)/MM activation free-energy barrier. The automated refinement of the EVB surface (and thus the reduction of the difference between the reference and ab initio potentials) is a key factor in accelerating the convergence of the LRA approach. We apply our PD approach to a test reaction, namely, the S(N)2 reaction between a chloride ion and methyl chloride, and demonstrate that, at present, this approach is far more powerful and cost-effective than the metadynamics approach (at least in its current implementation). We also discuss the general features of the PD approach in terms of its ability to explore complex systems and clarify that it is not a specialized approach limited to only accelerating QM(ai)/MM calculations with proper sampling, but rather can be used in a wide variety of applications. In fact, we point out that the use of a reference (CG) potential coupled with its PD refinement, as well as our renormalization approach, provides very general and powerful strategies that can be used very effectively to explore any property that has been studied by the MTD approach.
authors: Department of Chemistry (SGM418), University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089, United States.
authors: Nikolay V
authors: NV
authors: Plotnikov
authors: Shina C L
authors: SC
authors: Kamerlin
authors: Arieh
authors: A
authors: Warshel
references: J Am Chem Soc. 2010 Feb 10;132(5):1526-8
references: 20070076
references: Adv Protein Chem. 2003;66:249-313
references: 14631821
references: Proc Natl Acad Sci U S A. 2009 Oct 13;106(41):17359-64
references: 19805169
references: Annu Rev Phys Chem. 2002;53:291-318
references: 11972010
references: J Am Chem Soc. 2005 May 18;127(19):6940-1
references: 15884923
references: J Phys Chem A. 2005 Sep 1;109(34):7682-7
references: 16834142
references: Proc Natl Acad Sci U S A. 2008 Sep 16;105(37):13877-82
references: 18779576
references: J Phys Chem B. 2009 Nov 26;113(47):15692-8
references: 19888735
references: Annu Rev Phys Chem. 2011;62:41-64
references: 21034218
references: J Am Chem Soc. 2003 Apr 30;125(17):5025-39
references: 12708852
references: J Phys Chem B. 2006 Feb 16;110(6):2934-41
references: 16471904
references: Chemistry. 2007;13(24):6828-40
references: 17566132
references: Biochemistry. 2007 Dec 25;46(51):15076-88
references: 18052079
references: Nature. 1975 Feb 27;253(5494):694-8
references: 1167625
references: J Phys Chem B. 2008 May 8;112(18):5680-92
references: 18412414
references: J Phys Chem B. 2009 Jun 4;113(22):7867-73
references: 19432459
references: Phys Rev Lett. 2005 Apr 8;94(13):138302
references: 15904045
references: Proc Natl Acad Sci U S A. 2010 Nov 23;107(47):20340-5
references: 21059937
references: J Chem Phys. 2008 Jan 7;128(1):014109
references: 18190187
references: Annu Rev Phys Chem. 2008;59:573-601
references: 18393679
references: J Phys Chem B. 2009 Feb 5;113(5):1253-72
references: 19055405
references: Faraday Discuss. 2010;145:71-106
references: 25285029
references: J Mol Biol. 2002 Jan 25;315(4):713-36
references: 11812142
references: Proteins. 2010 Apr;78(5):1212-27
references: 20052756
references: Biochemistry. 2008 Mar 25;47(12):3725-35
references: 18307312
references: Proc Natl Acad Sci U S A. 2002 Oct 1;99(20):12562-6
references: 12271136
references: J Biomol NMR. 2007 Dec;39(4):265-73
references: 17929172
references: J Phys Chem B. 2005 Oct 20;109(41):19516-22
references: 16853521
references: J Comput Aided Mol Des. 1994 Dec;8(6):695-708
references: 7738605
references: Proteins. 2003 Aug 15;52(3):412-26
references: 12866052
references: Proc Natl Acad Sci U S A. 2009 May 5;106(18):7449-54
references: 19383795
references: Acc Chem Res. 2006 Feb;39(2):73-81
references: 16489726
references: J Chem Phys. 2005 Mar 1;122(9):091102
references: 15836102
references: J Phys Chem B. 2010 Oct 7;114(39):12720-8
references: 20836533
references: Chemphyschem. 2008 Aug 25;9(12):1767-73
references: 18666265
references: J Phys Chem B. 2006 Oct 5;110(39):19570-4
references: 17004821