Variational calculation of macrostate transition rates

Journal of Chemical Physics

Volumn 109, Issue 5, 1998, Pages 1670-1686

  Ulitsky, Alex ^a   Shalloway, David ^a  

^a Department of Obstetrics Gynecology   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

EID: 0242454276     PISSN: 00219606     EISSN: None     Source Type: Journal    
DOI: 10.1063/1.476882     Document Type: Article

References (101)

1. J. Kaelberer, J. D. Etters, and J. C. Raich, Chem. Phys. Lett. 41, 580 (1976).
2. F. G. Amar and R. S. Berry, J. Chem. Phys. 85, 5943 (1985).
3. A. W. Castleman and K. H. Bowen, Jr., J. Phys. Chem. 100, 12911 (1996).
4. C. J. Jameson, A. K. Jameson, and H.-M. Lim, J. Chem. Phys. 104, 1709 (1996).
5. J. M. Delaye and Y. Limoge, J. Phys. I 3, 2079 (1993).
6. A. Lewis, I. Rousso, E. Khachatryan, I. Brodsky, K. Lieberman, and M. Sheves, Biophys. J. 70, 2380 (1996).
7. S. Crouzy, T. B. Woolf, and B. Roux, Biophys. J. 67, 1370 (1994).
8. F. Bezanilla, E. Perozo, and E. Stefani, Science 254, 679 (1991).
9. B. A. Wallace, Annu. Rev. Biophys. Biophys. Chem. 19, 127 (1990).
10. E. Helfand, Science 226, 647 (1984).
11. J. A. McCammon and S. Harvey, Dynamics of Proteins and Nucleic Acids (Cambridge University Press, New York, 1987).
12. C. L. Brooks III, M. Karplus, and B. M. Pettitt, Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics (Wiley, New York, 1988).
13. Y. Bai, T. R. Sosnick, L. Mayne, and S. W. Englander, Science 269, 192 (1995).
14. F. Parak and H. Frauenfelder, Physica A 201, 332 (1993).
15. R. Kubo, M. Toda, and N. Hashitsume, Statistical Physics II: Nonequilibrium Statistical Mechanics (Springer-Verlag, New York, 1985).
16. 0(here), η(old)→ β/2 η(here), and ξ(old)→ β/2 ξ(here).
17. Throughout the paper, subscripted Greek letters identify macrostates and are unrelated to regular Greek letters such as the inverse temperature β. Summations over these macrostate indices are over all macrostates.
18. N. G. van Kampen, Stochastic Processes in Physics and Chemistry (North Holland, New York, 1992).
19. P. Hänggi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 62, 251 (1990).
20. ' S. Lee, M. Karplus, D. Bashford, and D. Weaver, Biopolymers 26, 481 (1987).
21. A. Heidenreich, J. Jortner, and I. Oref, J. Chem. Phys. 97, 197 (1992).
22. A. Heidenreich, I. Schek, D. Scharf, and J. Jortner, Z. Phys. D 20, 227 (1991).
23. Y. Zhang and R. W. Pastor, Mol. Simul. 13, 25 (1994).
24. P. Ahlstrom, J. Lausmaa, P. Löfgren, and H. Berendsen, in Modeling of Biomolecular Structures and Mechanisms, edited by A. Pullman, J. Jortner, and B. Pullman, Vol. 27 of The Jerusalem Symposia on Quantum Chemistry and Biochemistry (Kluwer Academic, Boston, 1994).
25. R. Elber, Curr. Opin. Struct. Biol. 6, 232 (1996).
26. M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford University Press, New York, 1992).
27. J. A. Board, L. V. Kale, K. Schulten, R. D. Skeel, and T. Schlick, IEEE Comput. Sci. Eng. 1, 19 (1995).
28. J. D. Madura, J. M. Briggs, R. C. Wade, M. E. Davis, B. A. Luty, A. Ilin, J. Antosiewicz, M. K. Gilson, B. Bagheri, L. R. Scott, and J. A. McCamman, Comput. Phys. Commun. 91, 57 (1995).
29. D. G. Truhlar, B. C. Garrett, and S. J. Klippenstein, J. Phys. Chem. 100, 12771 (1996).
30. M. A. Wilson and D. Chandler, Chem. Phys. 149, 11 (1990).
31. D. Brown and J. H. E. Clarke, J. Phys. Chem. 93, 4117 (1990).
32. E. P. Wigner, Trans. Faraday Soc. 34, 29 (1938).
33. D. A. McQuarie, Statistical Mechanics (Harper and Row, New York, 1976).
34. D. Chandler, J. Chem. Phys. 68, 2959 (1978).
35. E. Pollak, in Activated Barrier Crossing: Applications in Physics, Chemistry and Biology, edited by G. R. Fleming and P. Hänggi (World Scientific, Singapore, 1993), p. 5-41.
36. R. Elber and M. Karplus, Chem. Phys. Lett. 139, 375 (1987).
37. R. Czerminski and R. Elber, J. Chem. Phys. 92, 5580 (1990).
38. R. Czerminski and R. Elber, Int. J. Quantum Chem. 24, 167 (1990).
39. A. Ulitsky and R. Elber, J. Chem. Phys. 92, 1510 (1990).
40. H. Risken, The Fokker-Plank Equation (Springer-Verlag, New York, 1984).
41. B. Widom, J. Chem. Phys. 55, 44 (1971).
42. J. T. Bartis and B. Widom, J. Chem. Phys. 60, 3474 (1974).
43. H. Tomita, A. Ito, and H. Kidachi, Prog. Theor. Phys. 56, 786 (1976).
44. H. Tomita, A. Ito, and H. Kidachi, Prog. Theor. Phys. 55, 947 (1976).
45. N. G. van Kampen, J. Stat. Phys. 17, 71 (1977).
46. G. J. Moro, J. Chem. Phys. 103, 7514 (1995).
47. The problem is similar to that of determining the relationship between the molecular wave functions and the localized basis sets in quantum chemistry. However, the quantum chemistry problem is simply solved (e.g., by MO-LCAO approximation) because the localized basis functions are roughly factorizable in terms of atomic orbitals [Hartree-Fock-Ruthaan approximation; R. G. Parr, The Quantum Theory of Molecular Electronic Structure, (W. A. Benjamin, New York, 1964)]. This is not true for complicated multidimensional systems where the localized conformational macrostates are not factorizable in terms of the individual atomic coordinates.
48. α̇ J [P]dR, where P is the projection operator defined in Eq. (24). The current generated by the eigenfunction components that are outside the subspace corresponds to fast relaxation of the probability distribution within, not between, macrostates.
49. M. Orešič and D. Shalloway, J. Chem. Phys. 101, 9844 (1994).
50. β, resulting in an unbounded spectrum of relaxation rates which would lead to a singularity in the rate at t=0.
51. B. W. Church, M. Orešič, and D. Shalloway, in Global Minimization of Nonconvex Energy Functions: Molecular Conformation and Protein Folding: DIMACS Workshop, edited by P. Pardalos, D. Shalloway, and G. Xue, Vol. 23 of DIMACS Series in Discrete Mathematics and Theoretical Computer Science (American Mathematical Society, Providence, RI, 1996), pp. 41-64.
52. obs, so the minimum decay condition can equivalently be imposed in either the occupation or observation representation.
53. P. Lankaster, Theory of Matrices (Academic, New York, 1969).
54. See Ref. 16 for further discussion. If H is invariant under a continuous symmetry group (e.g., rigid body translation and rotation), we will only be interested in group-invariant macrostates; therefore the sums in Eq. (24) can be restricted to group- invariant eigenfunctions, and degeneracies associated with irreducible representations of the group of rank >1 will not enter. If H is invariant under a discrete group of transformations, the natural requirement that all group-transformation-related macrostates be included in the set of macrostates (i.e., that the macrostates form a basis for a representation of the group) implies that the macrostate subspace will include either all or none of the degenerate eigenfunctions belonging to an irreducible representation of the group. Thus, there will be no symmetry-induced ambiguity in the subspace specification.
55. This variational equation for σ is equivalent to the Rayleigh-Ritz variational equation for the first excited eigenfunction; see Ref. 16.
56. The detW term ensures that the node is located between the two macrostates; the denominator will become small if σ is varied so much that the nodal surface does not separate the macrostates. However, this does not occur during normal numerical variation.
57. 0 as R→∞ [A. Ulitsky and D. Shalloway (in preparation)]. Thus, by Eq. (15), the window and transition functions will go to constants.
58. Note that erf(0.83x)-tanh(x)<0.02∀x. Thus the functions provide essentially equivalent parameterizations for the variational condition.
59. For a system with N degrees-of-freedom, both the nodal surface and the transition "plane" are (N-1)-dimensional manifolds.
60. Since η is a vector, Eqs. (40) represents N+1 conditions where N is the number of degrees-of-freedom.
61. H. A. Kramers, Physica (Amsterdam) 7, 284 (1940).
62. X and Z instead of the characteristic harmonic frequency.
63. J. S. Langer, Ann. Phys. (N.Y.) 54, 258 (1969).
64. A. Nitzan and Z. Schuss, in Activated Barrier Crossing: Applications in Physics, Chemistry and Biology, edited by G. R. Fleming and P. Hänggi (World Scientific, Singapore, 1993), pp. 42-81.
65. M. H. Kalos and P. A. Whitlock, Monte Carlo Methods, Volume I: Basics (Wiley, New York, 1986).
66. W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in Fortran (Cambridge University Press, Cambridge, 1992).
67. A. Ulitsky and D. Shalloway, J. Chem. Phys. 106, 10099 (1997).
68. This method (see Ref. 66) was chosen because it uses only the objective function and does not require derivatives which are relatively expensive to accurately compute by Monte Carlo integration.
69. The integrals in the right-hand-sides of Eqs. (47) were computed to 3% accuracy as determined by the variance (Ref. 65).
70. XY<0.03.
71. ξ=0.03.
72. D. Shalloway, J. Chem. Phys. 105, 9986 (1996).
73. c.
74. C. W. Gardiner, Handbook of Stochastic Methods for Physics, Chemistry and Natural Sciences, 2nd ed. (Springer, New York, 1997).
75. 1 to 1% accuracy.
76. S. H. Northrup and J. A. McCammon, J. Chem. Phys. 78, 987 (1983).
77. S. Huo and J. E. Straub, J. Chem. Phys. 107, 5000 (1997).
78. R. D. Etters, R. Danilowicz, and J. B. Kaelberer, J. Chem. Phys. 67, 4145 (1977).
79. S. Sawada and S. Sugano, Z. Phys. D 14, 247 (1989).
80. M. R. Hoare and J. A. McInnes, Adv. Phys. 32, 791 (1983).
81. D. Shalloway (in preparation).
82. W. Kabsch, Acta Crystallogr. 32, 922 (1976).
83. t) from the transition region conformations resulting from the replacement of atoms 3 and 4 with other pairs of adjacent atoms.
84. This procedure described in A. Ulitsky and R. Elber, J. Chem. Phys. 92, 1510 (1990) was used. While multiple reaction paths are anticipated in multidimensional systems only one path connecting the two isoforms of the six-particle cluster was detected in 100 trials starting with different initial conditions.
85. As an additional check, transition rates were also calculated using the mean first-passage time evaluated from Langevin dynamics trajectories. Similar results were obtained.
86. 2 using Eq. (49b) since only the coordinates of one atom enter Eq. (61).
87. E. Pollak, J. Chem. Phys. 85, 865 (1986).
88. A. Berezhkovskii, E. Pollak, and V. Y. Zitserman, J. Chem. Phys. 97, 2422 (1992).
89. G. Gershinsky and E. Pollak, J. Chem. Phys. 101, 7174 (1994).
90. 0exp[-βw(q)]dq/∫exp[-βw(q)]dq, where w(q) is the potential of mean force along an ad hoc preselected direction and q=0 is the maximizer of w(q) in the barrier region.
91. D. Shalloway, J. Global Optim. 2, 281 (1992).
92. T. Coleman, D. Shalloway, and Z. Wu, Comput. Optim. Appl. 2, 145 (1993).
93. T. Coleman, D. Shalloway, and Z. Wu, J. Global Optim. 4, 171 (1994).
94. D. Shalloway, Variable-scale coarse-graining in macromolecular global optimization, in Proceedings of IMA Conference on Large Scale Optimization, edited by L. Biegler, T. Coleman, A. R. Conn, and F. Santosa, (Springer, New York, 1997).
95. J. E. Straub, M. Borkovec, and B. J. Berne, J. Chem. Phys. 89, 4833 (1988).
96. β≈1, but this has not been proven.
97. α to approach zero or one, but this has not been proven.
98. I. S. Gradshteyn and I. M. Ryzik, Table of Integrals, Series, and Products, 5th ed. (Academic, New York, 1994).
99. D. L. Ermak and J. A. McCammon, J. Chem. Phys 69, 1352 (1978).
100. R. W. Pastor, B. R. Brooks, and A. Szabo, Mol. Phys. 65, 1409 (1988).
101. 0).