Rare Events. Transition state theory Bennett-Chandler Approach 16.2 Transition path sampling16.4
|
|
- Arline Short
- 5 years ago
- Views:
Transcription
1 Rare Events Transition state theory Bennett-Chandler Approach 16.2 Transition path sampling16.4
2 Outline Part 1 Rare event and reaction kinetics Linear Response theory Transition state theory Free energy methods Bennet Chandler approach Example zeolites Part 2 Two ended methods Transition path sampling Rate constants Reaction coordinate analysis Application to biomolecules
3 crystallisation catalysis solution reactions complex fluids Rare Events enzyme reactions solvent effects folding & binding isomerization
4 Rare events Interesting transitions in complex systems solution chemistry protein folding F A B enzymatic reactions complex surface reactions diffusion in porous media nucleation q q B These reactions happen on a long time scale compared to the molecular timescale A dominated by collective, rare events Straightforward MD very inefficient mol t stable
5 Example: Diffusion in porous material
6 Phenomenological reaction kinetics A rare event can be seen as a chemical reaction between reactant A and product B A B A B The change in population c(t) is (0<c<1) dc A ( t) dt = k A B c ( A t) + k B A c ( B t) dc B dt ( t)!# A B " $ = 0 dt Equilibrium: c A t Total number change in population d c ( t) + c ( t) ( ) ( ) =+ k c t k c t ( ) = c B ( t) = 0 A B A B A B This gives a relation between equilibrium population and reaction rates c A c B = k B A k A B
7 Let us make a perturbation of the equilibrium populations. The dynamics of such an equilibrium fluctuation gives the desired information on the response of the system to an external field ( ) = +Δ ( ) For state A c t c c t For state B: A A A We can rewrite the kinetics in terms of the perturbation: 1 dδc A dt ( t) ( ) ( ) = k Δc t k Δc t A B A B A A ( ) ( 0exp ) ( ) Δ ca t = ΔcA "& ka B+ kb A t# ' With τ A ( 0exp ) [ t τ ] = Δc ( k k ) 1 1 = + ( ) 1 A B B A Relaxation time population = k 1+ c c = A B A B c A ( t) + c B ( t) =1 ( ) = Δ ( ) c t c c t B B A C A (t) C B (t) k c B A B time
8 Microscopic theory Microscopic description of the progress of a reaction q Reaction coordinate: in this case the z-coordinate of the particle We need to write the kinetics of the reaction in terms of this microscopic reaction coordinate q
9 A B Reaction coordinate Reactant A: q< q* Product B: q > q* Transition state: q = q* Let us introduce the function g A : Heaviside θ-function # 0 q q* < 0 θ ( q q* ) = $ % 1 q q* > 0 g q q* = 1 θ q q* = θ q* q A ( ) ( ) ( ) With this function we write for the probability c A (t) the system is in state A: c A ( t) = g ( A t)
10 Microscopic theory Is going to give us the macroscopic relaxation in terms of a microscopic time correlation function exp [ t τ ] = A ( 0) Δ ( ) Δg g t c A c B A
11 Perturbed Hamiltonian Let us consider the effect of a static perturbation: H = H g q q 0 ε A * ( ) This external potential increases the concentration of A For the equilibrium concentration as a function of ε: Δ c = c c = g A ε g A 0 A A ε A 0 We need to compute the ensemble average in the form of : A 0 = dγaexp dγexp [ β H0 ] [ β H ] 0
12 Linear Response theory (static) The original Hamiltonian (H 0 ) is perturbed by εd: H = H 0 ε D This gives as change in the expectation value of A: Δ A = A A 0 with A = ( ) dγ Aexp β H 0 ε D dγ exp β ( H 0 ε D) A 0 = dγaexp dγexp [ β H0 ] [ β H ] 0 If the perturbation is small we can write A = A + A 0 ε 0 ε
13 For such a small perturbation ΔA = A 0 ε = ε A ε 0 ε with A ε ( ) ( ) = A ε = dγβ ADexp! " β H 0 εd # $ dγexp! " β H 0 εd { dγexp! " β ( H 0 εd # ) } 2 $ ( ) dγ Aexp! " β H 0 εd ( ) # $ dγβdexp! " β H 0 εd ( ) { # } 2 dγexp! " β H 0 εd $ # $ # $ Evaluated for ε= 0 A = dγ β ADexp β H 0 ε 0 { dγexp β H 0 } dγ Aexp β H 0 dγβ Dexp β H 0 dγexp β H 0 dγexp β H 0 Giving: A = β { AD ε 0 A 0 D } 0 0
14 If we apply this result for c A : A = β { AD ε 0 A 0 D } 0 0 with H = H g q q 0 ε A * ( ) Δ c = g g A A ε A ( ) ( ) 2 2 ga ga 0 Δ ca = β 0 ε 0 Since g A =0 or 1: g A (x) g A (x) =g A (x) ( g ( 1 )) A ga = β 0 0 ( ( 1 )) = β c c = β c c A 0 A 0 A 0 B 0 Giving: Δc A = β c A 0 c B 0 ε
15 Linear Response theory (dynamic) Let us now switch off the perturbation at t=0 H = H 0 ε D H = H 0 at t>0: Let us see how the system relaxes to equilibrium (dynamical perturbation) We take <A> 0 =0 ΔA( t) = A( t) A 0 = A( t) Similar as for the static case for small values of ε, we have A( t) ε 0 = dγ β A t ( ) Dexp β H 0 { } dγ exp β H 0 = β D( 0) A( t) Giving: ΔA( t) = βε D( 0) A( t)
16 ΔA( t) = βε D( 0) A( t) If we apply this result to D = Δg A and A = Δg A We obtain: Δc A From static perturbation: Δc A ( t) = βε Δg ( A 0)Δg ( A t) ( t) = Δc ( A 0) Δg 0 A ( )Δg A t c A c B ( ) ( ) βε = Δc A 0 c A c B Compare linear response expression with the macroscopic expression Δc A ( t) = Δc A 0 ( )exp t τ
17 exp Derivative Microscopic rate expression [ t τ ] = A ( 0) Δ ( ) Δg g t 1 τ exp t τ = g 0 A c A c B A ( ) g ( A t) c A c B Stationary (t is arbitrary, only depends on τ) Δ has disappeared because of the derivative = g ( 0 A )g A t c A c B ( ) d dt A( t)b( t + τ ) = 0 A( t) B ( t + τ ) + A ( t) B( t + τ ) = 0 A( 0) B ( τ ) = A ( 0) B( τ )
18 We have 1 τ exp " # t τ $ % = g 0 A ( ) g ( A t) c A c B Using 1 τ = k A B 1+ c A c B ( ) 1 = c B k A B For sufficiently short t, we obtain Using the definition of g A we can write k A B ( t) = g ( 0 A )g A t c A ( ) ( ) = q g ( q q* ) A g A q q* q = q g ( q q* ) B q k A B ( t) = ( ) g q ( 0 ) q* B q 0 ( ) q c A g B ( t) We now have an expression that relates the macroscopic reaction rate to microscopic properties
19 k A B g B ( t) = ( ) g q ( 0 ) q* B q 0 ( t) = θ q t) ( ( ) q * ) g ( B q( 0) q *) q ( ) q c A = Θ q 0 g B ( t) ( ( ) q *) q ( ) = δ q( 0) q * Let us look at the different terms in this equation Only when the system is in the product state we get a contribution to the ensemble average Only when the system starts at the transition state, we get a contribution to the ensemble average ( ) Velocity at t=0 q 0 c A = Θ( q * q) Concentration of A k A B ( t) = q 0 ( )δ ( q( 0) q *)θ q t ( ( ) q *) ( ) θ q * q
20 k A B Transition state theory ( t) = q 0 ( )δ ( q( 0) q *)θ q t ( ( ) q *) ( ) θ q * q Let us consider the limit: t 0 + lim t 0 + ( ) = θ ( q ( t) ) = θ q( t) q* This gives: k TST A B = q ( 0)δ ( q( 0) q* )θ q θ ( q* q) ( ) Eyring s transition state theory
21 C(t) k A B k TST Decay of rate expression ( t) = q 0 ( )δ ( q( 0) q *)θ q t ( ( ) q *) ( ) θ q * q lower value because of recrossings t k(t) k AB ~exp(-t/! rxn )! mol t
22 k A B ( t) = Transition state theory q ( 0)δ ( q( 0) q* )θ q t ( ( ) q* ) ( ) θ q* q We can rewrite this expression as a product k A B ( t) = q ( 0)δ ( q( 0) q* )θ q t ( ) δ q( 0) q* ( ( ) q* ) ( ( ) q* ) ( ) δ q 0 θ q* q k A B Conditional probability to find a particle on the top of the barrier with a positive velocity ( t) = q 0 ( )θ q t ( ( ) q* ) q=q* Ratio of probabilities to find particle on top of the barrier and in the state A δ ( q( 0) q* ) θ ( q* q)
23 Free energy barrier ( ) δ q( 0) q * ( ) θ q * q Ratio of the probabilities to find a particle on top of the barrier and in the state A Probability to be on top of the barrier: ( ) = C dqδ q q* δ q* q Probability to be in state A: Θ q * q ( ) ( ( )) ( )exp βf ( q) = C exp βf q* ( ) = C dqθ q q * This gives: ( ) ( ( )) ( )exp β F ( q) = C dqexp β F q ( ) δ q( 0) q* ( ) θ q* q We need to determine the free energy as a function of the order parameter = q<q* q<q* ( ( )) ( ) exp βf q* ( ) dqexp βf q
24 q ( 0)θ q t ( ( ) q* ) q=q* Conditional probability to find a particle on the top of the barrier with a positive velocity q ( 0) Assume that on top of the barrier the particle is in equilibrium: use Maxwell-Boltzmann distribution to generate this velocity q ( 0)θ q t ( ( ) q* ) Only particles with a positive velocity end up in the lim t 0 + q ( 0)θ q t product state. We assume that once in the product state they stay there. ( ( ) q* ) = q 0 ( ( )) = 0.5 q ( 0) ( )θ q 0 k TST = lim q ( 0 A B t 0 )θ q t + k TST A B = 0.5 q ( 0) q<q* ( ( ) q* ) q=q* ( ) ( ( )) exp βf ( q* ) dqexp βf q δ ( q( 0) q* ) θ ( q* q) Eyring s TST
25 1-D ideal gas particle on a hill k TST A B = 0.5 q ( 0) q<q* ( ) ( ( )) exp βf ( q* ) dqexp βf q Maxwell-Boltzmann: q ( 0) = 2k T B π m This gives for the hopping rate k TST A B = k B T 2πm q<q* ( ) ( ( )) exp βu ( q* ) dqexp βu q
26 Ideal gas particle on a not-so-ideal hill q * is the true transition state q 1 is the estimated transition state
27 For this case transition state theory will overestimate the hopping rate
28 Transition state theory One has to know the free energy accurately (MC/MD) C(t) Gives only an upper bound to the reaction rate Assumptions underlying transition theory should hold: no recrossings k A B k TST k(t) ( t) = q 0 ( )δ ( q( 0) q *)θ q t ( ( ) q *) ( ) θ q * q k AB t lower value because of recrossings ~exp(-t/! rxn )! mol t
29 Outline Part 1 Rare event and reaction kinetics Linear Response theory Transition state theory Free energy methods Bennet Chandler approach Example zeolites Part 2 Two ended methods Transition path sampling Rate constants Reaction coordinate analysis Application to biomolecules
30 Free energy barriers in complex systems Straightforward MD or MC and then use is highly inefficient for high barriers βf(q) = ln δ ( q(r) q) Many tricks have been proposed to overcome and sample barriers Temperature enhanced sampling: simulated tempering, parallel tempering, Temperature accelerated molecular dynamics ) Constraint dynamics: thermodynamic integration, blue moon sampling... Flat histogram sampling: umbrella sampling, hyperdynamics,... history dependent search: Wang-Landau, adaptive biasing force, metadynamics, non-equilibrium methods: steered MD, targeted MD,... trajectory-based methods: nudged elastic band, action minimization, string method, transition path sampling, forward flux sampling,...
31 Replica exchange/parallel tempering F! room temperature! Q! F! high temperature! Q!
32 Two replicas T = 360K 2 e 2U 2 (r N ) Total Boltzmann weight! e 1U 1 (r N ) e 2U 2 (r N ) T = 290K 1 e 1U 1 (r N )
33 Switching temperatures T = 290K 2 e 1U 2 (r N ) Total Boltzmann weight! e 1U 2 (r N ) e 2U 1 (r N ) T = 360K 1 e 2U 1 (r N )
34 The ratio of the new Boltzmann factor over the old one is: N (n) N (o) = e( 2 1)[U 2 (r N ) U 2 (r N )] the rule for switching temperatures should obey detailed balance Metropolis Monte Carlo scheme acc(1 2) = min 1,e ( 2 1)[U 2 (r N ) U 1 (r N )]
35 Set of replicas T = 360K R T = 293K 2 e 2U 2 (r N ) T = 290K 1 e 1U 1 (r N )
36 Overlap in potential energy
37 N Replica Exchange MD (REMD) 2 e 2U 2 (r N ) 1 e 1U 1 (r N ) Hansmann Chem Phys Lett 1997! Sugita & Okamoto Chem Phys Lett 1999!
38 Replica Exchange Exchange as a function of time. Free energy follows from Advantage: no order parameters needed Disadvantage: convergence of free energy landscape can be still slow, especially around phase transition: many replicas needed.
39 Metadynamics method to obtain free energy in a single simulation similar idea as Wang Landau sampling: add history dependent biasing potential to forcefield Laio and Parrinello, PNAS (2002)! s = predefined order parameters w = height of hills σ = width of gaussians F(s) w is reduced every cycle F (s) = lim V (s; t) t more controlled version: well tempered MetaD!! Barducci, Bussi, Parrinello, PRL, (2008).! s
40 Link to bernds animation
41 SN2 reaction between Cl - and CH3Cl Reactant Transition Product Complex State Complex S(R) = r C-Cl r C-Cl` S 1 (R) = r C-Cl S 2 (R) = r C-Cl` Meta-dynamics can relax the requirement of choosing a good reaction coordinate Bernd Ensing, Alessandro Laio, Michele Parrinello and Michael L. Klein, J. Phys. Chem. B 109 (2005),
42 Outline Part 1 Rare event and reaction kinetics Linear Response theory Transition state theory Free energy methods Bennet Chandler approach Example zeolites Part 2 Two ended methods Transition path sampling Rate constants Reaction coordinate analysis Application to biomolecules
43 Problem with TST There are recrossings that cause overestimation of the rate constant trajectories that seem to overcome the barrier but in fact bounce back
44 Bennett-Chandler approach k A B k A B ( t) = ( t) = q ( 0)δ ( q( 0) q* )θ q t ( ( ) q* ) ( ) θ q* q q ( 0)δ ( q( 0) q* )θ q t ( ) δ q( 0) q* ( ( ) q* ) ( ( ) q* ) ( ) δ q 0 θ q* q Computational scheme: 1. Determine the probability from the free energy using MC or MD, e.g. by umbrella sampling, thermodynamic integration or other free energy methods 2. Compute the conditional average from a MD simulation
45 TST k A B k A B ( t) = ( t) = q ( 0)δ ( q( 0) q 1 )θ q δ q 0 Transmission coefficient ( ) k ( t ) A B κ t Bennett-Chandler approach = ( ) q 1 ( ) q 1 ( ) q ( 0)δ ( q 0 )θ q t TST k A B ( ) δ q( 0) q 1 q ( 0)δ ( q( 0) q 1 )θ q t 0.5 q 0 ( ) δ ( q ( 0 ) q ) 1 θ ( q 1 q) ( ( ) q ) 1 δ ( q ( 0 ) q ) 1 θ ( q 1 q) ( ( ) q ) 1 ( ) MD simulation to correct the transition state result! MD simulation: 1. At t=0 q=q 1 2. Determine fraction at product state weighted with initial velocity
46 Example diffusion in zeolite Zeolites important class of materials Diffusion of alkanes in matrix is poorly described Approach molecular simulation of alkanes in fixed zeolite frame Unified atom FF by Dubbeldam et al. D. Dubbeldam, et al., J. Phys. Chem. B, 108, 12301, 2004
47 t : θ=1 For both q 0 ( ) and q 0 ( ) q ( 0)δ ( q( 0) q* )θ q t ( ( ) q* ) Low value of κ
48 Reaction coordinate cage window cage q * cage window cage q * βf(q) βf(q) q q
49
50
51 Outline Part 1 Rare event and reaction kinetics Linear Response theory Transition state theory Free energy methods Bennet Chandler approach Example zeolites Part 2 Two ended methods Transition path sampling Rate constants Reaction coordinate analysis Application to biomolecules
52 Barriers on smooth and rough energy landscapes Clearly, barrier is most important for rare event But how to obtain this barrier? In multidimensional energy landscapes barrier is saddle point # saddle points limited determined by potential energy use eigenvectors or Hessian to find them # saddle points uncountable entropy important, many pathways determined by free energy exploring requires sampling schemes Dellago logo TM
53 Breakdown of BC approach kappa can become immeasurable low if the reaction coordinate is at the wrong value the reaction coordinate is wrongly chosen W ( q) = kt ln dq" exp{ βe( q, q" )} If the reaction coordinate is not known, the wrong order parameter can lead to wrong transition states, mechanism and rates
54 Two ended methods Methods that take the entire path and fix the begin and end point Many methods proposed: Action minimization Nudged elastic band String method Path metadynamics Milestoning Transition path sampling...
55 Transition path sampling Samples the path ensemble: all trajectories that lead over barrier Sampling by Monte Carlo Requires definition of stable states A,B only Results in ensemble of pathways Reaction coordinate is a result of simulation not an input Allows for calculation of rate constants C. Dellago, P.G. Bolhuis, P.L. Geissler Adv. Chem. Phys. 123, Apply when process of interest is a rare event is complex and reaction coordinate is not known Examples: nucleation, reactions in solution, protein folding
56 Path probability density Path = Sequence of states x iδt
57 Transition path ensemble h A =1 h B =1
58 Metropolis MC of pathways 1. Generate new path from old one 2. Accept new path according to detailed balance: 3. Satisfy detailed balance with the Metropolis rule:
59 Shooting moves ) ( ) ( )] ( ) ( [ ) ( ) ( 0 ) ( ) ( n T B n A n o acc x h x h T x T x P = accept reject! " # = A x if A x if t h t t A 0 1 ) (
60 Standard TPS algorithm take existing path choose random time slice t change momenta slightly at t integrate forward and backward in time to create new path of length L accept if A and B are connected, otherwise reject and retain old path calculate averages repeat
61 Definition of the stable states
62 Classical nucleation (1926) ΔG 0 R Liquid R * R Crystal nucleus surface bulk ΔG = 4πR 2 γ 4 3 πr3 ρδµ ls How does the crystal form? What is the structure of the critical nucleus Is classical nucleation theory correct? What is the barrier? Rate constant γ : surface tension Δµ : chem. pot difference ρ: density
63 Path sampling of nucleation TIS in NPH ensemble, as density and temperature change N=10000, P=5.68 H=1.41 (25 % undercooling) D. Moroni, P. R. ten Wolde, and P. G. Bolhuis, Phys. Rev. Lett. 94, (2005)
64 Sampling paths is only the beginning Eugene Wigner: "It is nice to know that the computer understands the problem. But I would like to understand it too. Path ensemble needs to be further explored to obtain: Rate constants Free energy Transition state ensembles Mechanistic picture Reaction coordinate Illustrative example: crystal nucleation
65 Transition interface sampling B A Overall states in phase space: directly coming from A directly coming from B T. S. van Erp, D. Moroni and P. G. Bolhuis, J. Chem. Phys. 118, 7762 (2003) T. S. van Erp and P. G. Bolhuis, J. Comp. Phys. 205, 157 (2005)
66 A B λ i-1 λ i λ i+1 P A (λ i+1 λ i ) = probability that path crossing i for first time after leaving A reaches i+1 before A k TIS AB = φ AB h A n 1 = φ AB P A (λ i+1 λ i ) = Φ A P A (λ i+1 λ i ) h A i=1 n 1 i=1
67 TIS results for nucleation Free energy follows directly Moroni, van Erp, Bolhuis, PRE, 2005 Structural analysis?
68 Committor (aka p-fold, splitting probability) Probability that a trajectory initiated at r relaxes into B r B A L. Onsager, Phys. Rev. 54, 554 (1938). M. M. Klosek, B. J. Matkowsky, Z. Schuss, Ber. Bunsenges. Phys. Chem. 95, 331 (1991) V. Pande, A. Y. Grosberg, T. Tanaka, E. I. Shaknovich, J. Chem. Phys. 108, 334 (1998)
69 Transition state ensemble r is a transition state (TS) if p B (r) = p A (r) =0.5 A B A B TSE: Intersections of transition pathways with the p B =1/2 surface
70 Committor distributions
71 Committor distribution N=243 Clearly, n is not entire story
72 Small and structured Structure Big and unstructured Committor analysis gives valuable insight
73 Outline Part 1 Rare event and reaction kinetics Linear Response theory Transition state theory Free energy methods Bennet Chandler approach Example zeolites Part 2 Two ended methods Transition path sampling Rate constants Reaction coordinate analysis Application to biomolecules
74 All-atom force fields for biomolecules Potential energy for protein V (r) = k r (r r eq ) 2 + k θ (θ θ eq ) i< j ' ) ( bonds a ij r ij 12 b ij r ij 6 + q i q j εr ij *, + angles 1 2 v (1+ cos(nφ φ ) n 0 dihedrals θ( 3 r 2,3 φ( 4 4 vdw interactions only between non-bonded i-j >4
75 Currently available empirical force fields CHARMm (MacKerrel et 96) AMBER (Cornell et al. 95) GROMOS (Berendsen et al 87) OPLS-AA (Jorgensen et al 95) ENCAD (Levitt et al 83) Subtle differences in improper torsions, scale factors 1-4 bonds, united atom rep. Partial charges based on empirical fits to small molecular systems Amber & Charmm also include ab-initio calculations Not clear which FF is best : top 4 mostly used Water models also included in description TIP3P, TIP4P SPC/E Current limit: 10 6 atoms, microseconds ( with Anton ms)
76 Folding of Trp-cage 20-residue protein NLYIQ WLKDG GPSSG RPPPS 2-state folder, experimental rate 4 µs (Andersen et al, Nature 2002, Zhou et al. PNAS 2004, others) System: 1L2Y in 2800 waters OPLSAA, PME, Nose-Hoover, GROMACS one time per 4 s! What is folding mechanism and kinetics in explicit water at 300K? Strategy: - Stable states by replica exchange MD - Mechanism by path sampling - Rate calculation Free! Energy! about equal population of unfolded and folded!
77 TPS of Trp-cage folding loop-helix 80% Rates by TIS! folding! helix-loop 20% unfolding! Jarek Juraszek! J. Juraszek, PGB PNAS 2006! Biophys. J (2008)! Rates! TIS! exp! k fol! 0.2 s -1! 0.24 s -1!! k unf! 0.8 s -1! 0.08 s -1!
78 Advanced path sampling results Very recent single replica multiple state transition interface sampling revealed much more complex network with many short-lived intermediates 1.0" U" Lo".2" Commi5or" other" SN" LSN" I".1" Pd" LN".001" meta" 0" N" W.Du & PGB, J. Chem. Phys 2014
79 Photoactive Yellow Protein bacterial blue-light sensor Tyr42 pca Glu46 Cys69 N-terminal domain PAS domain Thr50 Arg52 Absorption of a blue-light photon triggers the photo cycle J. Vreede et al. Biophys. J. 2005
80 Partial unfolding Loss of α-helical structure Exposure of hydrophobic groups Increased flexibility in parts of the protein backbone cis-chroh + Glu46- Force field MD (gromacs) Gromos96 - SPC water - PME Replica Exchange H 2 O H 2 O - H 2 O H 2 O H 2 O H 2 O
81 Partial unfolding: REMD free energy Unfolding happens on millisecond timescale (1000 times longer than possible by direct MD) Helix unfolded! Chromophore and Glu46 are solvent exp Agreement with experiment! Vreede et al. Biophys. J. 2005! Vreede et al. Proteins 2008! Bernard et al. Structure 2005! Jocelyne Vreede!
82 Partial unfolding: REMD free energy Unfolding happens on sub-millisecond timescale (1000 times longer than possible by direct MD) Helix unfolded! Chromophore and Glu46 are solvent exp Agreement with experiment! Vreede et al. Biophys. J. 2005! Vreede et al. Proteins 2008! Bernard et al. Structure 2005! Jocelyne Vreede!
83 Transitions in the partial unfolding
84 Solvent exposure transitions Vreede, Juraszek, PGB, PNAS 2010!
85 Juraszek, Vreede, PGB, Chem Phys 2011!
86 The end
Generating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniques
Generating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniques Peter Bolhuis van t Hoff institute for Molecular Sciences University of Amsterdam, The Netherlands
More informationOptimal reaction coordinates, kinetics, and free energy of rare events by reweighted path sampling
Optimal reaction coordinates, kinetics, and free energy of rare events by reweighted path sampling Peter Bolhuis van t Hoff institute for Molecular Sciences University of Amsterdam, The Netherlands Outline
More informationRare Event Simulations
Rare Event Simulations Homogeneous nucleation is a rare event (e.g. Liquid Solid) Crystallization requires supercooling (µ solid < µ liquid ) Crystal nucleus 2r Free - energy gain r 3 4! 3! GBulk = ñr!
More informationOn the Reaction Rate Calculation Methods using Path Sampling
On the Reaction Rate Calculation Methods using Path Sampling Paris, France France 28 November, 2006 Titus S. van Erp Centre for Surface Chemistry and Catalysis, KULeuven Leuven, Belgium Introduction Straightforward
More informationLecture 14: Advanced Conformational Sampling
Lecture 14: Advanced Conformational Sampling Dr. Ronald M. Levy ronlevy@temple.edu Multidimensional Rough Energy Landscapes MD ~ ns, conformational motion in macromolecules ~µs to sec Interconversions
More informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationComputing free energies with PLUMED 2.0
Computing free energies with PLUMED 2.0 Davide Branduardi Formerly at MPI Biophysics, Frankfurt a.m. TyOutline of the talk Relevance of free energy computation Free-energy from histograms: issues Tackling
More informationFree energy simulations
Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational
More informationMajor Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1
Major Concepts Onsager s Regression Hypothesis Relaxation of a perturbation Regression of fluctuations Fluctuation-Dissipation Theorem Proof of FDT & relation to Onsager s Regression Hypothesis Response
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate heory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationComputing free energy: Thermodynamic perturbation and beyond
Computing free energy: Thermodynamic perturbation and beyond Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls
More informationDynamics of Nucleation in the Ising Model
J. Phys. Chem. B 2004, 108, 19681-19686 19681 Dynamics of Nucleation in the Ising Model Albert C. Pan and David Chandler* Department of Chemistry, UniVersity of California, Berkeley, California 94720-1460
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationMolecular Dynamics and Accelerated Molecular Dynamics
Molecular Dynamics and Accelerated Molecular Dynamics Arthur F. Voter Theoretical Division National Laboratory Lecture 3 Tutorial Lecture Series Institute for Pure and Applied Mathematics (IPAM) UCLA September
More informationFree energy calculations and the potential of mean force
Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationComputing free energy: Replica exchange
Computing free energy: Replica exchange Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls Macroscopic i a t e regime
More informationDynamical Rare event simulation techniques for equilibrium and non-equilibrium systems arxiv: v2 [cond-mat.stat-mech] 31 Mar 2011
Dynamical Rare event simulation techniques for equilibrium and non-equilibrium systems arxiv:1101.0927v2 [cond-mat.stat-mech] 31 Mar 2011 Contents Titus S. van Erp Centrum voor Oppervlaktechemie en Katalyse,
More informationExploring the Free Energy Surface of Short Peptides by Using Metadynamics
John von Neumann Institute for Computing Exploring the Free Energy Surface of Short Peptides by Using Metadynamics C. Camilloni, A. De Simone published in From Computational Biophysics to Systems Biology
More informationMulti-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018
Multi-Ensemble Markov Models and TRAM Fabian Paul 21-Feb-2018 Outline Free energies Simulation types Boltzmann reweighting Umbrella sampling multi-temperature simulation accelerated MD Analysis methods
More informationSome Accelerated/Enhanced Sampling Techniques
Some Accelerated/Enhanced Sampling Techniques The Sampling Problem Accelerated sampling can be achieved by addressing one or more of the MD ingredients Back to free energy differences... Examples of accelerated
More informationBiomolecular modeling. Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010)
Biomolecular modeling Marcus Elstner and Tomáš Kubař Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010) IX. FREE ENERGY SIMULATIONS When searching for a physical quantity that is of most
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationKinetic Pathways of Ion Pair Dissociation in Water
3706 J. Phys. Chem. B 1999, 103, 3706-3710 Kinetic Pathways of Ion Pair Dissociation in Water Phillip L. Geissler, Christoph Dellago, and David Chandler* Department of Chemistry, UniVersity of California,
More informationChoosing weights for simulated tempering
PHYSICAL REVIEW E 7, 173 7 Choosing weights for simulated tempering Sanghyun Park 1 and Vijay S. Pande 1, 1 Department of Chemistry, Stanford University, Stanford, California 935, USA Department of Structural
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
More informationLimitations of temperature replica exchange (T-REMD) for protein folding simulations
Limitations of temperature replica exchange (T-REMD) for protein folding simulations Jed W. Pitera, William C. Swope IBM Research pitera@us.ibm.com Anomalies in protein folding kinetic thermodynamic 322K
More informationMilestoning: Extending time scales of molecular simulations
Milestoning: Extending time scales of molecular simulations Anthony West, Peter Majek, and RE Department of Chemistry and Biochemistry and Institute of Computational Sciences and Engineering (ICES) University
More information3.320 Lecture 23 (5/3/05)
3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials
More informationMany proteins spontaneously refold into native form in vitro with high fidelity and high speed.
Macromolecular Processes 20. Protein Folding Composed of 50 500 amino acids linked in 1D sequence by the polypeptide backbone The amino acid physical and chemical properties of the 20 amino acids dictate
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationPerforming Metadynamics Simulations Using NAMD
Performing Metadynamics Simulations Using NAMD Author: Zhaleh Ghaemi Contents 1 Introduction 2 1.1 Goals of this tutorial.......................... 2 2 Standard metadynamics simulations 3 2.1 Free energy
More informationThe PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics
The PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics Davide Branduardi, Theoretical Molecular Biophysics Group Max Planck for Biophysics, Frankfurt am Main, Germany
More informationControlling fluctuations
Controlling fluctuations Michele Parrinello Department of Chemistry and Applied Biosciences ETH Zurich and ICS, Università della Svizzera Italiana, Lugano, Switzerland Today s menu Introduction Fluctuations
More informationOutline. The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation. Unfolded Folded. What is protein folding?
The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation By Jun Shimada and Eugine Shaknovich Bill Hawse Dr. Bahar Elisa Sandvik and Mehrdad Safavian Outline Background on protein
More informationEnergy Landscapes and Accelerated Molecular- Dynamical Techniques for the Study of Protein Folding
Energy Landscapes and Accelerated Molecular- Dynamical Techniques for the Study of Protein Folding John K. Prentice Boulder, CO BioMed Seminar University of New Mexico Physics and Astronomy Department
More informationReaction Dynamics (2) Can we predict the rate of reactions?
Reaction Dynamics (2) Can we predict the rate of reactions? Reactions in Liquid Solutions Solvent is NOT a reactant Reactive encounters in solution A reaction occurs if 1. The reactant molecules (A, B)
More informationarxiv: v1 [q-bio.bm] 9 Jul 2007
Exploring the Protein G Helix Free Energy Surface by Solute Tempering Metadynamics Carlo Camilloni, 1, 2 Davide Provasi, 1, 2 Guido Tiana, 1, 2, and Ricardo A. Broglia 1,2,3 1 Department of Physics, University
More informationMolecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)
Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments) Fei Gao gaofeium@umich.edu Limitations of MD Time scales Length scales (PBC help a lot) Accuracy of
More informationMethods of Computer Simulation. Molecular Dynamics and Monte Carlo
Molecular Dynamics Time is of the essence in biological processes therefore how do we understand time-dependent processes at the molecular level? How do we do this experimentally? How do we do this computationally?
More informationLocal Hyperdynamics. Arthur F. Voter Theoretical Division Los Alamos National Laboratory Los Alamos, NM, USA
Local Hyperdynamics Arthur F. Voter Theoretical Division National Laboratory, NM, USA Summer School in Monte Carlo Methods for Rare Events Division of Applied Mathematics Brown University Providence, RI
More informationFree energy calculations
Free energy calculations Berk Hess May 5, 2017 Why do free energy calculations? The free energy G gives the population of states: ( ) P 1 G = exp, G = G 2 G 1 P 2 k B T Since we mostly simulate in the
More informationMethods for predicting phase transition: VCNEB and TPS approaches
Methods for predicting phase transition: VCNEB and TPS approaches Guangrui Qian Aug 7, 2013 Department of Geosciences, Stony Brook University Materials at High Pressure New physics and new chemistry New
More informationAtomistic Modeling of Small-Angle Scattering Data Using SASSIE-web
Course Introduction Atomistic Modeling of Small-Angle Scattering Data Using SASSIE-web September 21-23, 2016 Advanced Photon Source Argonne National Laboratory, Argonne, IL ccpsas.org Scatters & Simulators
More informationProtein Folding experiments and theory
Protein Folding experiments and theory 1, 2,and 3 Protein Structure Fig. 3-16 from Lehninger Biochemistry, 4 th ed. The 3D structure is not encoded at the single aa level Hydrogen Bonding Shared H atom
More informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationCurrent address: Department of Chemistry, Hong Kong Baptist University, Kowloon Tong, Hong Kong,
Hydrolysis of Cisplatin - A Metadynamics Study Supporting Information Justin Kai-Chi Lau a and Bernd Ensing* b Department of Chemistry and Applied Bioscience, ETH Zurich, USI Campus, Computational Science,
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationIdentifying the molecular mechanism of a reaction in solution,
Reaction coordinates and rates from transition paths Robert B. Best and Gerhard Hummer* Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes
More informationOn the efficient sampling of pathways in the transition path ensemble
On the efficient sampling of pathways in the transition path ensemble Thijs J. H. Vlugt{ and Berend Smit* Department of Chemical Engineering, University of Amsterdam, Nieuwe Achtergracht 166, 1018 WV Amsterdam,
More informationComputational Chemistry - MD Simulations
Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated
More informationCan a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?
Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Ruhong Zhou 1 and Bruce J. Berne 2 1 IBM Thomas J. Watson Research Center; and 2 Department of Chemistry,
More informationComputational Studies of the Photoreceptor Rhodopsin. Scott E. Feller Wabash College
Computational Studies of the Photoreceptor Rhodopsin Scott E. Feller Wabash College Rhodopsin Photocycle Dark-adapted Rhodopsin hn Isomerize retinal Photorhodopsin ~200 fs Bathorhodopsin Meta-II ms timescale
More informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationObtaining reaction coordinates by likelihood maximization
THE JOURNAL OF CHEMICAL PHYSICS 125, 054108 2006 Obtaining reaction coordinates by likelihood maximization Baron Peters and Bernhardt L. Trout a Department of Chemical Engineering, Massachusetts Institute
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
More informationUnderstanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry
JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore
More informationWhat makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces: Electronic Supplementary
Electronic Supplementary Material (ESI) for Journal of Materials Chemistry B. This journal is The Royal Society of Chemistry 21 What makes a good graphene-binding peptide? Adsorption of amino acids and
More informationMonte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects
Monte Carlo Methods Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Lecture 1 2 Lecture 1&2 3 Lecture 1&3 4 Different Ensembles Ensemble ame Constant (Imposed) VT Canonical,V,T P PT
More informationEnhanced sampling via molecular dynamics II: Unbiased approaches
Enhanced sampling via molecular dynamics II: Unbiased approaches Mark E. Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Sciences New York University, 100 Washington Square East, NY
More informationGeneralized Ensembles: Multicanonical Simulations
Generalized Ensembles: Multicanonical Simulations 1. Multicanonical Ensemble 2. How to get the Weights? 3. Example Runs and Re-Weighting to the Canonical Ensemble 4. Energy and Specific Heat Calculation
More informationAb initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
More informationTRANSITION PATH SAMPLING: Throwing Ropes Over Rough Mountain Passes, in the Dark
Annu. Rev. Phys. Chem. 2002. 53:291 318 DOI: 10.1146/annurev.physchem.53.082301.113146 Copyright c 2002 by Annual Reviews. All rights reserved TRANSITION PATH SAMPLING: Throwing Ropes Over Rough Mountain
More informationFree energy calculations using molecular dynamics simulations. Anna Johansson
Free energy calculations using molecular dynamics simulations Anna Johansson 2007-03-13 Outline Introduction to concepts Why is free energy important? Calculating free energy using MD Thermodynamical Integration
More informationMolecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods
Molecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods Charles L. Brooks III MMTSB/CTBP 2006 Summer Workshop CHARMM Simulations The flow of data and information
More informationWater models in classical simulations
Water models in classical simulations Maria Fyta Institut für Computerphysik, Universität Stuttgart Stuttgart, Germany Water transparent, odorless, tasteless and ubiquitous really simple: two H atoms attached
More informationComputer simulations as concrete models for student reasoning
Computer simulations as concrete models for student reasoning Jan Tobochnik Department of Physics Kalamazoo College Kalamazoo MI 49006 In many thermal physics courses, students become preoccupied with
More informationNucleation rate (m -3 s -1 ) Radius of water nano droplet (Å) 1e+00 1e-64 1e-128 1e-192 1e-256
Supplementary Figures Nucleation rate (m -3 s -1 ) 1e+00 1e-64 1e-128 1e-192 1e-256 Calculated R in bulk water Calculated R in droplet Modified CNT 20 30 40 50 60 70 Radius of water nano droplet (Å) Supplementary
More informationApplication of the Markov State Model to Molecular Dynamics of Biological Molecules. Speaker: Xun Sang-Ni Supervisor: Prof. Wu Dr.
Application of the Markov State Model to Molecular Dynamics of Biological Molecules Speaker: Xun Sang-Ni Supervisor: Prof. Wu Dr. Jiang Introduction Conformational changes of proteins : essential part
More informationMolecular Dynamics Simulation of Methanol-Water Mixture
Molecular Dynamics Simulation of Methanol-Water Mixture Palazzo Mancini, Mara Cantoni University of Urbino Carlo Bo Abstract In this study some properties of the methanol-water mixture such as diffusivity,
More informationBiology Chemistry & Physics of Biomolecules. Examination #1. Proteins Module. September 29, Answer Key
Biology 5357 Chemistry & Physics of Biomolecules Examination #1 Proteins Module September 29, 2017 Answer Key Question 1 (A) (5 points) Structure (b) is more common, as it contains the shorter connection
More informationMolecular Interactions F14NMI. Lecture 4: worked answers to practice questions
Molecular Interactions F14NMI Lecture 4: worked answers to practice questions http://comp.chem.nottingham.ac.uk/teaching/f14nmi jonathan.hirst@nottingham.ac.uk (1) (a) Describe the Monte Carlo algorithm
More information(Crystal) Nucleation: The language
Why crystallization requires supercooling (Crystal) Nucleation: The language 2r 1. Transferring N particles from liquid to crystal yields energy. Crystal nucleus Δµ: thermodynamic driving force N is proportional
More informationOptimized statistical ensembles for slowly equilibrating classical and quantum systems
Optimized statistical ensembles for slowly equilibrating classical and quantum systems IPAM, January 2009 Simon Trebst Microsoft Station Q University of California, Santa Barbara Collaborators: David Huse,
More informationMulti-scale approaches in description and design of enzymes
Multi-scale approaches in description and design of enzymes Anastassia Alexandrova and Manuel Sparta UCLA & CNSI Catalysis: it is all about the barrier The inside-out protocol: Big Aim: development of
More informationWang-Landau Monte Carlo simulation. Aleš Vítek IT4I, VP3
Wang-Landau Monte Carlo simulation Aleš Vítek IT4I, VP3 PART 1 Classical Monte Carlo Methods Outline 1. Statistical thermodynamics, ensembles 2. Numerical evaluation of integrals, crude Monte Carlo (MC)
More informationEssential dynamics sampling of proteins. Tuorial 6 Neva Bešker
Essential dynamics sampling of proteins Tuorial 6 Neva Bešker Relevant time scale Why we need enhanced sampling? Interconvertion between basins is infrequent at the roomtemperature: kinetics and thermodynamics
More informationMarkov Processes. Stochastic process. Markov process
Markov Processes Stochastic process movement through a series of well-defined states in a way that involves some element of randomness for our purposes, states are microstates in the governing ensemble
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationMolecular Dynamics in Principal Component Space
pubs.acs.org/jpcb Molecular Dynamics in Principal Component Space Servaas Michielssens,*, Titus S. van Erp, Carsten Kutzner, Arnout Ceulemans, and Bert L. de Groot Department of Chemistry and Centre for
More informationCourse Introduction. SASSIE CCP-SAS Workshop. January 23-25, ISIS Neutron and Muon Source Rutherford Appleton Laboratory UK
Course Introduction SASSIE CCP-SAS Workshop January 23-25, 2017 ISIS Neutron and Muon Source Rutherford Appleton Laboratory UK ccpsas.org Scatters & Simulators Develop a community of users and developers
More informationAnalysis of the simulation
Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities
More informationKinetic Monte Carlo: from transition probabilities to transition rates
Kinetic Monte Carlo: from transition probabilities to transition rates With MD we can only reproduce the dynamics of the system for 100 ns. Slow thermallyactivated processes, such as diffusion, cannot
More informationReplica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2)
pubs.acs.org/jpcb Replica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2) Lingle Wang, Richard A. Friesner, and B. J. Berne* Department of Chemistry,
More informationMolecular dynamics (MD) and the Monte Carlo simulation
Escaping free-energy minima Alessandro Laio and Michele Parrinello* Centro Svizzero di Calcolo Scientifico, Via Cantonale, CH-6928 Manno, Switzerland; and Department of Chemistry, Eidgenössische Technische
More informationEquilibrium Molecular Thermodynamics from Kirkwood Sampling
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source
More informationExploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials
Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials Dr. V. Subramanian Chemical Laboratory, IPC Division CSIR-Central Leather Research
More informationDominant Paths in Protein Folding
Dominant Paths in Protein Folding Henri Orland SPhT, CEA-Saclay France work in collaboration with P. Faccioli, F. Pederiva, M. Sega University of Trento Henri Orland Annecy meeting 2006 Outline Basic notions
More informationEquilibrium Molecular Thermodynamics from Kirkwood Sampling
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source
More informationFree energy, electrostatics, and the hydrophobic effect
Protein Physics 2016 Lecture 3, January 26 Free energy, electrostatics, and the hydrophobic effect Magnus Andersson magnus.andersson@scilifelab.se Theoretical & Computational Biophysics Recap Protein structure
More informationFree energy calculations with alchemlyb
Free energy calculations with alchemlyb Oliver Beckstein Arizona State University SPIDAL Teleconference 2019-02-01 Binding free energy measure of how strong a protein P and a ligand X stick together key
More informationWhy study protein dynamics?
Why study protein dynamics? Protein flexibility is crucial for function. One average structure is not enough. Proteins constantly sample configurational space. Transport - binding and moving molecules
More informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More informationComputing stationary distributions in equilibrium and nonequilibrium systems with forward flux sampling
THE JOURNAL OF CHEMICAL PHYSICS 127, 114109 2007 Computing stationary distributions in equilibrium and nonequilibrium systems with forward flux sampling Chantal Valeriani a FOM Institute for Atomic and
More informationMonte Carlo (MC) Simulation Methods. Elisa Fadda
Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is
More information