200 likes | 357 Vues
This presentation discusses innovative methods for simulating protein folding, utilizing Molecular Dynamics and Wang-Landau techniques to analyze complex trajectories. We explore strategies for constructing transition matrices that reveal transition states and patterns within protein folding processes. By studying a tripeptide system through GROMACS, we investigate the probabilities of state transitions and analyze the resulting Ramachandran Plot for insights into protein conformations. This work aims to bridge molecular simulations with statistical analysis to enhance our understanding of protein behavior.
E N D
Simulating Protein Folding - Some Ideas Christian Hedegaard Jensen Dmitry Nerukh
Outline • Introduction • Transition Matrix • Wang-Landau
Outline • Introduction • Transition Matrix • Wang-Landau
Introduction • We are simulation protein folding using Molecular Dynamics, and understanding using complexity analysis. • Try to do statistics on trajectories, to see if there are any patterns in where they go. (Transition Matrix) • Try and look at the paths that a protein can take when it folds. (Wang-Landau)
Outline • Introduction • Transition Matrix • Wang-Landau
Transition Matrix • Find states in trajectory • Find transitions between states in trajectory • Write up transition matrix for system • Possibly include some memory into the system.
Transition Matrix • We have a molecular dynamics simulation of a tripeptide. • Trianaline • 300 K, 1 atm • 12 ps (2fs) • OPLSAA-2001 force field • GROMACS
Transition Matrix Ramachandran Plot
Transition Matrix • Gives useful information about what is going on in the system. • Find out where we end and with what probabilities. • Look at transitions between clusters of states.
Outline • Introduction • Transition Matrix • Wang-Landau
Wang-Landau • Wang-Landau simulations, using a combination of Molecular Dynamics and Monte Carlo.
Wang-Landau • Wang-Landau gives the (relative) density of states. • The dos of a simple protein was obtained, d(E). • Try to calculate d(,E) or d(, ,E) .
dos q 0 2p Wang-Landau
Wang-Landau • Define reaction coordinates • Sample as function of these coordinates • Method does not say what the system does, but is says what is possible.
Wang-Landau • Matlab