Classic techniques for simulating molecular motion, such as Monte Carlo method and molecular dynamics, generate individual motion pathways one at a time and are inefficient if applied in a naive fashion to deal with many pathways. In this paper, we introduce stochastic roadmap simulation (SRS), a new approach for exploring the kinetics of molecular motion by examining multiple pathways simultaneously. In SRS, we compactly encode many pathways in a graph, called a roadmap. Every path in the roadmap represents a potential motion pathway and is associated with a probability indicating the likelihood that a molecule may follow the path. By viewing the roadmap as a Markov chain, we can efficiently compute kinetic properties of molecular motion over the entire energy landscape. Furthermore we prove that in the limit, SRS converges to the same distribution as Monte Carlo simulation. To test the effectiveness of our method, we applied it in the computation of the transmission coefficient for protein folding, which is an important order parameter that measures the “kinetic distance ” of a conformation to the folded state of a protein. Our computational studies demonstrate that compared with Monte Carlo method, SRS obtains more accurate results and achieves several orders-of-magnitude reduction in running time. 1
|
1120
|
Equation of state calculations by fast computing machines
– Metropolis, Rosenbluth, et al.
- 1953
|
|
721
|
Iterative Methods for Sparse Linear Systems
– Saad
- 2003
|
|
716
|
Probability inequalities for sums of bounded random variables
– Hoeffding
- 1963
|
|
493
|
Probabilistic Roadmaps for Path Planning in High-Dimensional Configuration Spaces
– Kavraki, Svetska, et al.
- 1996
|
|
152
|
Cooling schedules for optimal annealing
– Hajek
- 1988
|
|
113
|
der Stappen. The gaussian sampling strategy for probabilistic roadmap planners
– Boor, Overmars, et al.
- 1999
|
|
108
|
The evolution of the minimum degree ordering algorithm
– George, Liu
- 1989
|
|
98
|
Principles that govern folding of protein chains
– Anfinsen
- 1973
|
|
89
|
An introduction to stochastic modeling
– Taylor, Karlin
- 1998
|
|
87
|
Sparse matrices in Matlab: Design and implementation
– GILBERT, MOLER, et al.
- 1991
|
|
69
|
Molecular Modelling – Principle and Applications
– Leach
- 1991
|
|
44
|
Hierarchical Protein Structure Superposition using both Secondary Structure and Atomic Representations
– Singh, Brutlag
- 1997
|
|
42
|
Using motion planning to study protein folding pathways
– Amato, Song
|
|
41
|
A Motion Planning Approach to Flexible Ligand Binding
– Singh, Latombe, et al.
- 1999
|
|
37
|
A simple protein folding algorithm using a binary code and secondary structure constraints
– Sun, Thomas, et al.
- 1995
|
|
36
|
Stochastic roadmap simulation: An efficient representation and algorithm for analyzing molecular motion
– Apaydin, Brutlag, et al.
- 2002
|
|
36
|
Using motion planning to map protein folding landscapes and analyze folding kinetics of known native structures
– Amato, Song
- 2001
|
|
36
|
Molecular Dynamics Simulation: Elementary Methods
– Haile
- 1997
|
|
36
|
Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function
– Morris, Goodsell, et al.
- 1998
|
|
35
|
Funnels, pathways and the energy landscape of protein folding: A synthesis,” Proteins
– Bryngelson, Onuchic, et al.
- 1995
|
|
32
|
Electrostatic interactions in macromolecules: Theory and applications
– Sharp, Honig
- 1990
|
|
28
|
Structure and Mechanism in Protein Science: A Guide to Enzyme Catalysis and Protein Folding
– Fersht
- 1999
|
|
26
|
B.: Numerical Recipes in C. Cambridge University Press
– Press, Teukolsky, et al.
- 1993
|
|
20
|
Capturing molecular energy landscapes with probabilistic conformational roadmaps
– Apaydin, Singh, et al.
- 2001
|
|
19
|
Protein docking along smooth association pathways
– Camacho, Vajda
|
|
16
|
From Levinthal to Pathways to Funnels
– Dill, Chan
- 1997
|
|
13
|
Ligand binding with OBPRM and haptic user input: Enhancing automatic motion planning with virtual touch
– Bayazit, Song, et al.
- 2001
|
|
12
|
On the transition coordinate for protein folding
– Du, Pande, et al.
- 1998
|
|
11
|
Monte Carlo Methods, volume 1
– Kalos, Whitlock
- 1986
|
|
10
|
The folding thermodynamics and kinetics of crambin using an all-atom Monte Carlo simulation
– Shimada, Kussell, et al.
- 2001
|
|
10
|
Screen savers of the world, unite
– Shirts, Pande
|
|
8
|
Stochastic roadmap simulation for the study of ligand-protein interactions
– Apaydin, Guestrin, et al.
- 2002
|
|
8
|
et al., “The protein databank: a computer-based archival file for macromolecular structures
– Bernstein, Koetzle, et al.
- 1977
|
|
7
|
The rates of defined changes in protein structure during the catalytic cycle of lactate dehydrogenase
– Clarke, Waldman, et al.
- 1985
|
|
7
|
Lactate dehydrogenase
– Holbrook, Liljas, et al.
- 1975
|
|
7
|
The fundamentals of protein folding: bringing together theory and experiment
– Dobson, Karplus
- 1999
|
|
7
|
Pathways for protein folding: Is a new view needed
– Pande, Grosberg, et al.
- 1998
|
|
6
|
Lattice Models of Protein Folding, Dynamics and Thermodynamics. Chapmann
– Kolinski, Skolnick
- 1996
|
|
5
|
An investigation of the contribution made by the carboxylate group of an active site histidine-aspartate couple to binding and catalysis in lactate dehydrogenase
– Clarke, Barstow, et al.
- 1988
|
|
5
|
Site-directed mutagenesis reveals the role of a mobile arginine residue in lactate dehydrogenase catalysis
– Clarke, Wigley, et al.
- 1986
|
|
5
|
Design and synthesis of new enzymes based on the lactate dehydrogenase framework
– Dunn, Wilks, et al.
- 1991
|
|
5
|
A strong carboxylate-arginine interaction is important in substrate orientation and recognition in lactate dehydrogenase
– Hart, Clarke, et al.
- 1987
|
|
4
|
A dimensionality reduction approach to modeling protein flexibility
– Phillips, Kavraki
- 2002
|
|
4
|
A specific, highly active malate dehydrogenase by redesign of a lactate dehydrogenase framework
– Wilks, Hart, et al.
- 1988
|
|
3
|
Blue Gene Team. Blue gene: A vision for protein science using a petaflop supercomputer
– IBM
|
|
1
|
Funnels, pathways, and the energy landscape of protein folding: A synthesis
– Socci, Wolynes
- 1995
|
|
1
|
The evolution of the minimum degree ordering algorithm
– Company
|
|
1
|
The folding thermodynamics and kinetics of crambin using an allatom monte carlo simulation
– Shakhnovich
|
|
1
|
on Intelligent Systems for Molecular Biology
– Sun, Thomas, et al.
- 1999
|