Probabilistic Roadmaps: A Tool for Computing Ensemble Properties of Molecular Motions Serkan Apaydin, Doug Brutlag 1 Carlos Guestrin, David Hsu 2 Jean-Claude Latombe, Chris Varma Computer Science Department Stanford University 1 Department of Biochemistry, Stanford University 2 Computer Science Department, University of North Carolina
42
Embed
Probabilistic Roadmaps: A Tool for Computing Ensemble Properties of Molecular Motions Serkan Apaydin, Doug Brutlag 1 Carlos Guestrin, David Hsu 2 Jean-Claude.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Probabilistic Roadmaps: A Tool for Computing
Ensemble Properties of Molecular Motions
Serkan Apaydin, Doug Brutlag1
Carlos Guestrin, David Hsu2
Jean-Claude Latombe, Chris Varma
Computer Science DepartmentStanford University
1 Department of Biochemistry, Stanford University2 Computer Science Department, University of North Carolina
Goal of our ResearchGoal of our Research
Develop efficient computational representations and algorithms to study molecular pathways for protein folding and ligand-protein binding
Protein folding RECOMB ’02Ligand-protein binding ECCB ‘02
AcknowledgementsAcknowledgements
People: Leo Guibas Michael Levitt, Structural BiologyItay Lotan Vijay Pande, ChemistryFabian Schwarzer Amit SinghRohit Singh
Funding: NSF-ITR ACI-0086013Stanford’s Bio-X and Graduate Fellowship
The probability that a roadmap fails to correctly capture the connectivity of the
free space goes to 0 exponentially in the number of milestones (~ running time).
Random sampling is convenient incremental
scheme for approximating the free space
Computed ExamplesComputed Examples
Biology Biology Robotics Robotics
Energy field, instead of joint controlContinuous energy field, instead of binary free and in-collision spacesMultiple pathways, instead of single collision-free pathPotentially many more degrees of freedomRelation to real world is more complex
Study of ligand-protein bindingProbabilistic roadmaps with edges weighted by energetic plausibilitySearch of most plausible pathsStudy of energy profiles along such paths
Study of ligand-protein bindingProbabilistic roadmaps with edges weighted by energetic plausibilitySearch of most plausible pathsStudy of energy profiles along such pathsExtensions to protein folding[Song and Amato, 01] [Apaydin et al., 01]
New Idea: New Idea: Capture the stochastic nature of molecular Capture the stochastic nature of molecular motion by assigning probabilities to edgesmotion by assigning probabilities to edges
vi
vj
Pij
Why is this a good idea?Why is this a good idea?
1) We can approximate Monte Carlo simulation as closely as we wish
2) Unlike with MC simulation, we avoid the local-minima problem
3) We can consider all pathways in the roadmap at once to compute ensemble properties
Edge probabilitiesEdge probabilities
Follow Metropolis criteria:
otherwise. ,
1
;0 if ,)/exp(
i
iji
Bij
ij
N
EN
TkE
P
Self-transition probability:
ijijii PP 1
vj
vi
Pij
Pii
Stochastic simulation on roadmap and Monte Carlo simulation converge to same Boltzmann distribution
Chemical environment of LDH-NADH-substrate complex (pyruvate) (catalyzes conversion of pyruvate to lactate in the presence of NADH
CH3
Some amino acids aredeleted entirely, replaced by other amino acids, or sidechains altered
Binding of Pyruvate to LDHBinding of Pyruvate to LDH
ASP-195HIS-193
ASP-166
ARG-169
+
+
+
THR-245
C
C
OO
O
CH3
NADH
GLN-101
ARG-106Loop
ResultsResults
C
C
OO
O
GLN-101
ARG-106
ASP-195HIS-193
ASP-166
ARG-169
NADH
+
+
+
Loop
CH3
THR-245
Mutant Escape Time
Change
Wildtype 3.216E6 N/A
ResultsResults
C
C
OO
O
GLN-101
ALA-106
ASP-195ALA-193
ASP-166
ARG-169
NADH
+
Loop
CH3
Mutant Escape Time
Change
Wildtype 3.216E6 N/A
His193 AlaArg106 Ala
4.126E2
ResultsResults
Mutant Escape Time Change
Wildtype 3.216E6 N/A
His193 AlaArg106 Ala
4.126E2
His193 Ala 3.381E3
Arg106 Ala 2.550E2
Asp195 Asn 5.221E7
Gln101 Arg 1.669E6 No change
Thr245 Gly 4.607E5
C
C
OO
O
GLN-101
ARG-106
ASP-195HIS-193
ASP-166
ARG-169
NADH
+
+
+
Loop
CH3
GLY-245
ConclusionConclusion
Probabilistic roadmaps are a promising computational tool for studying ensemble properties of molecular pathwaysCurrent and future work: Better kinetic/energetic models Experimentally verifiable tests Non-uniform sampling strategies Encoding MD simulation
Stochastic simulation on a roadmap and MC simulation converge to the same distribution (Boltzman):For any set S, >0, >0,>0, there exists N such that a roadmap with N milestones has error bounded by:
• DOF = 10 – 3 coordinates to position root atom;– 2 angles to specify first bond;– Angles for all remaining non-terminal atoms;– Bond angles are assumed constant;
• Protein assumed rigid[Singh, Latombe and Brutlag `99]
x,y,z
Energy of InteractionEnergy of Interaction
Ev
Rij
Ec
Rij
Ev = 0.2[(R0/Rij)12 - 2(R0/Rij)
6 ]Ec = 332 QiQj/(Rij)
Energy = van der Waals interaction (Ev)
+ electrostatic interaction (Ec)
Solvent Effects
• Is only valid for an infinite medium of uniform dielectric;• Dielectric discontinuities result in induced surface
charges;
Solution: Poisson-Boltzman equation
Ec = 332 QiQj/(Rij)
Use Delphi [Rocchia et al `01] Finite Difference solution is based on discretizing