Improving replica exchange using driven scaling
J. Chem. Phys. 131, 174113 (2009); doi:10.1063/1.3259049
Published 6 November 2009
You are not logged in to this journal. Log in
Replica exchange is a powerful simulation method in which simulations are run at a series of temperatures, with the highest temperature chosen so phase space can be sampled efficiently. In order for swaps to be accepted, the energy distributions of adjacent replicas must have some overlap. This can create the need for many replicas for large systems. In this paper, we present a new method in which the potential energy is scaled by a parameter, which has an explicit time dependence. Scaling the potential energy broadens the distribution of energy and reduces the number of replicas necessary to span a given temperature range. We demonstrate that if the system is driven by the time-dependent potential sufficiently slowly, then equilibrium is maintained and energetic and structural properties are identical to those of conventional replica exchange. The method is tested using two systems, the alanine dipeptide and the trpzip2 polypeptide, both in water.
©2009 American Institute of Physics
| History: | Received 20 August 2009; accepted 15 October 2009; published 6 November 2009 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/131/174113/1 |
KEYWORDS and PACS
PUBLICATION DATA
0021-9606 (print)
1089-7690 (online)
REFERENCES (51)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- D. J. Earl and M. W. Deem,
Phys. Chem. Chem. Phys. 7, 3910 (2005) . - H. Fukunishi, O. Watanabe, and S. Takada, J. Chem. Phys. 116, 9058 (2002).
- J. W. Pitera, I. Haque, and W. C. Swope, J. Chem. Phys. 124, 141102 (2006).
- R. Zhou, B. Berne, and R. Germain,
Proc. Natl. Acad. Sci. U.S.A. 98, 14931 (2001) . - J. Zhang, M. Qin, and W. Wang,
Proteins 62, 672 (2006) . - S. Yang, J. N. Onuchic, A. E. Garcia, and H. Levine,
J. Mol. Biol. 372, 756 (2007) . - R. Zhou and B. J. Berne, J. Chem. Phys. 107, 9185 (1997).
- I. Andricioaei, J. E. Straub, and A. F. Voter, J. Chem. Phys. 114, 6994 (2001).
- S. Brown and T. Head-Gordon,
J. Comput. Chem. 24, 68 (2003) . - X. Li, C. P. O'Brien, G. Collier, N. A. Vellore, F. Wang, R. A. Latour, D. A. Bruce, and S. J. Stuart, J. Chem. Phys. 127, 164116 (2007).
- E. Lyman and D. M. Zuckerman, J. Chem. Phys. 127, 065101 (2007).
- Y. Sugita and Y. Okamoto,
Chem. Phys. Lett. 329, 261 (2000) . - F. Calvo and J. P. K. Doye, Phys. Rev. E 63, 010902 (2000).
- R. Faller, Q. Yan, and J. J. de Pablo, J. Chem. Phys. 116, 5419 (2002).
- M. K. Fenwick and F. A. Escobedo, J. Chem. Phys. 119, 11998 (2003).
- A. Mitsutake and Y. Okamoto, J. Chem. Phys. 121, 2491 (2004).
- T. W. Whitfield, L. Bu, and J. E. Straub,
Physica A 305, 157 (2002) . - S. Jang, S. Shin, and Y. Pak, Phys. Rev. Lett. 91, 058305 (2003).
- P. Liu, B. Kim, R. A. Friesner, and B. J. Berne,
Proc. Natl. Acad. Sci. U.S.A. 102, 13749 (2005) . - P. Liu, X. Huang, R. Zhou, and B. J. Berne,
J. Phys. Chem. B 110, 19018 (2006) . - S. Kannan and M. Zacharias,
Proteins 66, 697 (2007) . - S. W. Rick, J. Chem. Phys. 126, 054102 (2007).
- X. Kong and C. L. Brooks III, J. Chem. Phys. 105, 2414 (1996).
- R. M. Lynden-Bell and J. C. Rasaiah, J. Chem. Phys. 107, 1981 (1997).
- Z. Guo, C. L. Brooks III, and X. Kong,
J. Phys. Chem. B 102, 2032 (1998) . - N. F. A. van der Vegt and W. J. Briels, J. Chem. Phys. 109, 7578 (1998).
- S. W. Rick and R. E. Cachau, J. Chem. Phys. 112, 5230 (2000).
- U. Börjesson and P. H. Hünenberger, J. Chem. Phys. 114, 9706 (2001).
- M. S. Lee, F. R. Salsbury, and C. L. Brooks III,
Proteins 56, 738 (2004) . - J. Khandogin and C. L. Brooks III,
Biophys. J. 89, 141 (2005) . - J. B. Abrams, L. Rosso, and M. E. Tuckerman, J. Chem. Phys. 125, 074115 (2006).
- C. Predescu, M. Predescu, and C. V. Ciobanu,
J. Phys. Chem. B 109, 4189 (2005) . - S. Trebst, M. Troyer, and U. H. E. Hansmann, J. Chem. Phys. 124, 174903 (2006).
- N. Rathore, M. Chopra, and J. J. de Pablo, J. Chem. Phys. 122, 024111 (2005).
- A. Kone and D. A. Kofke, J. Chem. Phys. 122, 206101 (2005).
- A. G. Cochran, N. J. Skelton, and M. A. Starovasnik,
Proc. Natl. Acad. Sci. U.S.A. 98, 5578 (2001) . - N. Metropolis, A. Rosenbluth, M. Rosenbluth, A. Teller, and E. Teller, J. Chem. Phys. 21, 1087 (1953).
- W. L. Jorgensen, D. S. Maxwell, and J. Tirado-Rives,
J. Am. Chem. Soc. 118, 11225 (1996) . - G. A. Kaminski, R. A. Friesner, J. Tirado-Rives, and W. L. Jorgensen,
J. Phys. Chem. B 105, 6474 (2001) . - W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 79, 926 (1983).
- Y. Sugita and Y. Okamoto,
Chem. Phys. Lett. 314, 141 (1999) . - J. Wang, P. Cieplak, and P. Kollman,
J. Comput. Chem. 21, 1049 (2000) . - V. Hornak, R. Abel, A. Okur, B. Strockbine, A. Roitberg, and C. Simmerling,
Proteins 65, 712 (2006) . - D. A. Case, T. A. Darden, T. E. Cheatham III, C. L. Simmerling, J. Wang, R. E. Duke, R. Luo, K. M. Merz, D. A. Pearlman, M. Crowley et al., AMBER9, University of California, San Francisco, 2006.
- In Fig. 1 of Ref. 22, the x-axis is incorrectly labeled. It should read 0.02, 0.04, 0.06, 0.08, and 0.10.
- S. B. Opps and J. Schofield, Phys. Rev. E 63, 056701 (2001).
- W. Zhang, C. Wu, and Y. Duan, J. Chem. Phys. 123, 154105 (2005).
- X. Periole and A. E. Mark, J. Chem. Phys. 126, 014903 (2007).
- D. Sindhikara, Y. Meng, and A. E. Roitberg, J. Chem. Phys. 128, 024103 (2008).
- M. J. Abraham and J. E. Gready,
J. Chem. Theory Comput. 4, 1119 (2008) . - The C7eq/C5 region is defined as
<0 or
>130 and
<−110 or
>80;
R/
2 is
<0 or
>130 and −110<
<80;
L is 0<
<130 and −60<
<90; and C7ax is 0<
<130 and
<−60 or
>90.








