CROSS-DISCIPLINARY PHYSICS AND RELATED AREAS OF SCIENCE AND TECHNOLOGY |
|
|
|
|
Calculation of Collective Variable-based PMF by Combining WHAM with Umbrella Sampling |
XU Wei-Xin1,2**, LI Yang1, ZHANG John Z. H.1,2,3 |
1State Key Laboratory of Precision Spectroscopy, Department of Physics, East China Normal University, Shanghai 200062 2Institute of Theoretical and Computational Science, Institutes for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062 3Department of Chemistry, New York University, New York, New York 10003, USA |
|
Cite this article: |
XU Wei-Xin, LI Yang, ZHANG John Z. H. 2012 Chin. Phys. Lett. 29 068702 |
|
|
Abstract Potential of mean force (PMF) with respect to localized reaction coordinates (RCs) such as distance is often applied to evaluate the free energy profile along the reaction pathway for complex molecular systems. However, calculation of PMF as a function of global RCs is still a challenging and important problem in computational biology. We examine the combined use of the weighted histogram analysis method and the umbrella sampling method for the calculation of PMF as a function of a global RC from the coarse-grained Langevin dynamics simulations for a model protein. The method yields the folding free energy profile projected onto a global RC, which is in accord with benchmark results. With this method rare global events would be sufficiently sampled because the biased potential can be used for restricting the global conformation to specific regions during free energy calculations. The strategy presented can also be utilized in calculating the global intra- and intermolecular PMF at more detailed levels.
|
|
Received: 05 March 2012
Published: 31 May 2012
|
|
PACS: |
87.15.ap
|
(Molecular dynamics simulation)
|
|
87.15.hm
|
(Folding dynamics)
|
|
87.15.Cc
|
(Folding: thermodynamics, statistical mechanics, models, and pathways)
|
|
87.14.E-
|
(Proteins)
|
|
|
|
|
[1] Kirkwood J G 1935 J. Chem. Phys. 3 300 [2] Wang J and Wang W 1999 Nature Struct. Biol. 6 1033 [3] Kollman P 1993 Chem. Rev. 93 2395 [4] Torrie G M and Valleau J P 1977 J. Comput. Phys. 23 187 [5] Miller T F and Clary D C 2003 J. Chem. Phys. 119 68 [6] Ensing B, De Vivo M, Liu Z W, Moore P and Klein M L 2006 Acc. Chem. Res. 39 73 [7] Chipot C, Kollman P A and Pearlman D A 1996 J. Comput. Chem. 17 1112 [8] Souaille M and Roux B 2001 Comput. Phys. Commum. 135 40 [9] Xu W X, Wang J and Wang W 2005 Proteins 61 777 [10] Zuo G H, Wang J and Wang W 2006 Proteins 63 165 [11] Xu W X, Wang J and Wang W 2005 Chin. Phys. Lett. 22 258 [12] Lu Q and Wang J 2008 J. Am. Chem. Soc. 130 4772 [13] Wang W, Xu W X, Levy Y, Trizacc E and Wolynes P G 2009 Proc. Natl. Acad. Sci. USA 106 5517 [14] Sun L, Wang J and Wang W 2010 Chin. Phys. Lett. 27 038702 [15] Ferrenberg A M and Swendsen R H 1989 Phys. Rev. Lett. 63 1195 [16] Xu W X, Jiang P and Mu Y G 2009 J. Phys. Chem. B 113 7308 [17] Zhang Y J, Zhang J and Wang W 2011 J. Am. Chem. Soc. 133 6882 [18] Oliveira R J, Whitford P C, Chahine J, Leite V B P and Wang J 2010 Methods 52 91 [19] Xu W X, Lai Z Z, Oliveira R J, Leite V B P and Wang J 2012 J. Phys. Chem. B 116 5152 [20] Wang J and Verkhivker G M 2003 Phys. Rev. Lett. 90 188101 [21] Wang J, Zheng X, Yang Y, Drueckhammer D, Yang W, Verkhivker G and Wang E 2007 Phys. Rev. Lett. 99 198101 |
|
Viewed |
|
|
|
Full text
|
|
|
|
|
Abstract
|
|
|
|
|