Computing Best Transition Pathways in High-Dimensional Dynamical Systems: Application to the AlphaL \leftrightharpoons Beta \leftrightharpoons AlphaR Transitions in Octaalanine

The direct computation of rare transitions in high-dimensional dynamical systems such as biomolecules via numerical integration or Monte Carlo is limited by the sampling problem. Alternatively, the dynamics of these systems can be modeled by transition networks (TNs) which are weighted graphs whose edges represent transitions between stable states of the system. The computation of the globally best transition paths connecting two selected stable states is straightforward with available graph-theoretical methods. However, these methods require that the energy barriers of all TN edges be determined, which is often computationally infeasible for large systems. Here, we introduce energy-bounded TNs, in which the transition barriers are specified in terms of lower and upper bounds. We present algorithms permitting the determination of the globally best paths on these TNs while requiring the computation of only a small subset of the true transition barriers. Several variants of the algorithm are given which ach...

[1]  M. Blue,et al.  Unified approach to fuzzy graph problems , 2002, Fuzzy Sets Syst..

[2]  David J Wales,et al.  Finding pathways between distant local minima. , 2005, The Journal of chemical physics.

[3]  J. Straub,et al.  The MaxFlux algorithm for calculating variationally optimized reaction paths for conformational transitions in many body systems at finite temperature , 1997 .

[4]  David J. Wales,et al.  Energy Landscape of a Model Protein , 1999, cond-mat/9904304.

[5]  Frank Noé,et al.  Automated computation of low‐energy pathways for complex rearrangements in proteins: Application to the conformational switch of Ras p21 , 2005, Proteins.

[6]  D. Wales.  Discrete path sampling , 2002 .

[7]  T Schlick,et al.  Biomolecular dynamics at long timesteps: bridging the timescale gap between simulation and experimentation. , 1997, Annual review of biophysics and biomolecular structure.

[8]  Ron Elber,et al.  A method for determining reaction paths in large molecules: application to myoglobin , 1987 .

[9]  R. Elber,et al.  Reaction path study of conformational transitions in flexible systems: Applications to peptides , 1990 .

[10]  O. Becker,et al.  Solvent effects on the energy landscapes and folding kinetics of polyalanine , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[11]  Laxmikant V. Kale,et al.  Algorithmic Challenges in Computational Molecular Biophysics , 1999 .

[12]  D. Frenkel.  Introduction to Monte Carlo Methods , 2004 .

[13]  Daniele Frigioni,et al.  Fully Dynamic Algorithms for Maintaining Shortest Paths Trees , 2000, J. Algorithms.

[14]  O. Becker,et al.  Energy landscapes of conformationally constrained peptides , 2001 .

[15]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[16]  Joshua Jortner,et al.  Dynamics of hierarchical folding on energy landscapes of hexapeptides , 2001 .

[17]  D. Wales,et al.  From Topographies to Dynamics on Multidimensional Potential Energy Surfaces of Atomic Clusters , 1996, Science.

[18]  Thomas A. Weber,et al.  Hidden structure in liquids , 1982 .

[19]  Wilhelm Huisinga,et al.  Biomolecular Conformations can be Identified as Metastable Sets of Molecular Dynamics , 2003 .

[20]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[21]  R Elber,et al.  Reaction path study of conformational transitions and helix formation in a tetrapeptide. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[22]  M. Karplus,et al.  The topology of multidimensional potential energy surfaces: Theory and application to peptide structure and kinetics , 1997 .

[23]  David J. Wales,et al.  The free energy landscape and dynamics of met-enkephalin , 2003 .

[24]  P. Hänggi,et al.  Reaction-rate theory: fifty years after Kramers , 1990 .

[25]  Jonathan P. K. Doye,et al.  Stationary points and dynamics in high-dimensional systems , 2003 .

[26]  H. Jónsson,et al.  Nudged elastic band method for finding minimum energy paths of transitions , 1998 .

[27]  David J Wales,et al.  Archetypal energy landscapes: dynamical diagnosis. , 2005, The Journal of chemical physics.

[28]  Yaakov Levy,et al.  Effect of Conformational Constraints on the Topography of Complex Potential Energy Surfaces , 1998 .

[29]  David J Wales,et al.  Folding of the GB1 hairpin peptide from discrete path sampling. , 2004, The Journal of chemical physics.

[30]  C L Brooks,et al.  Taking a Walk on a Landscape , 2001, Science.

[31]  Edsger W. Dijkstra,et al.  A note on two problems in connexion with graphs , 1959, Numerische Mathematik.

[32]  Jing-Shing Yao,et al.  Fuzzy Shortest-Path Network Problems With Uncertain Edge Weights , 2003, J. Inf. Sci. Eng..

[33]  M. Karplus,et al.  A Comprehensive Analytical Treatment of Continuum Electrostatics , 1996 .

[34]  David J. Wales,et al.  Free energy landscapes of model peptides and proteins , 2003 .

[35]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[36]  Berry,et al.  Topography and Dynamics of Multidimensional Interatomic Potential Surfaces. , 1995, Physical review letters.

[37]  John D. Morgan,et al.  Diffusion‐controlled reactions: A variational formula for the optimum reaction coordinate , 1983 .

[38]  David J. Wales,et al.  Energy landscapes of model polyalanines , 2002 .

[39]  M. Karplus,et al.  Conjugate peak refinement: an algorithm for finding reaction paths and accurate transition states in systems with many degrees of freedom , 1992 .

[40]  David Eppstein,et al.  Finding the k Shortest Paths , 1999, SIAM J. Comput..

[41]  F. Stillinger,et al.  A Topographic View of Supercooled Liquids and Glass Formation , 1995, Science.

[42]  F. Stillinger,et al.  Packing Structures and Transitions in Liquids and Solids , 1984, Science.