Free energy calculations and the potential of mean force
|
|
- Ursula Adams
- 6 years ago
- Views:
Transcription
1 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 100 Washington Square East New York University, New York, NY 10003
2 Free Energy Canonical Ensemble (Helmholtz free energy): 1 QNVT (,, ) = d p d r e = e N! h N N βh( pr, ) βa( N, V, T) 3 N DV ( ) A( NVT,, ) = ktln QNVT (,, ) Isothermal-Isobaric Ensemble (Gibbs free energy): 1 ( N, P, T) dv d d e = e VN! h = 0 N N β( H( pr, ) + PV) βg( N, P, T) p r 3N 0 DV ( ) GNPT (,, ) = ktln ( NPT,, )
3 Free Energy (cont d) State function: A21 = A2 A1 P 2 1 V T
4 Free energy and work If an amount of work W is required to change the thermodynamic state of the system from 1 to 2, then W A Equality holds when the work is performed infinitely slowly or reversibly. Jarzynski s equality [PRL, , (1997)] shows how to relate irreversible work to the free energy difference. Let W 21 (x) be a microscopic function whose ensemble average is the thermodynamic work W 21. e βw β A = e
5 d 1 d 2 Free energy profiles A q= d1 d2 k = κe β A
6 Protein Folding Energetics From G. Bussi, et al. JACS 128, (2006)
7 Binding Free Energies Inhibition constant: U ( r, r ) = U ( r ) + U ( r ) 1 U ( r, r ) = U ( r ) + U ( r ) + U ( r, r ) 2 K i [ E][ I] [ EI ] = = E I E E I I G bind E I E E I I EI E I e β Thermodynamic state potentials: Meta-potential: U( r, r, λ) = f( λ) U ( r, r ) + g( λ) U ( r, r ) E I 1 E I 2 E I f(0) = 1 f(1) = 0 g(0) = 0 g(1) = 1 G = bind 1 0 dλ U λ Thermodynamic integration (Kirkwood, 1935) λ
8 Binding free energies: Thermodynamic perturbation 1 N N β H(, ) Z( NVT,, ) QNVT (,, ) = d pr p d r e = 3N ( ) 3N N! h DV N! λ Z N V T = d r e λ = βh πm ( ) 2 (,, ) N βu r / 2 DV ( ) Free energy difference related to partition function ratio: Perturbation formula: Q A21 = ktln = ktln Q β ( U U ) 2 1 Z Z Z2 1 1 = dr e = dr e e Z Z Z = e βu ( r) βu ( r) β( U ( r) U ( r)) Need sufficient overlap between two ensembles
9 λ dynamics methods Use molecular dynamics to sample λ via a Hamiltonian: H λ 2 2 pλ pi = + + U( r1,..., rn, λ) 2m 2m λ i i 2 2 pλ pi = + + f( λ) U1( r1,..., rn ) + g( λ) U2( r1,..., rn ) 2m 2m λ i Free energy from probability distribution of λ: i P( λ) = dr e β λ U ( r, ) A( λ) = ktln P( λ) A = A(1) A(0) 21 Need to have best sampling at the endpoints of the λ-path, which are normally the most difficult to sample.
10 Aim for a profile with a barrier: A( λ) λ dynamics methods λ = 0 λ = 1 In order to generate such a profile, we need: 1. A high temperature T λ >> T to ensure barrier crossing 2. An adiabatic decoupling between λ and other degrees of freedom 3. Choose m λ >> m i.
11 λ dynamics methods Under adiabatic conditions, we generate a free energy profile at T λ β β λ λ βλ A( λ; β) βa( λ; β) β βu( r, λ) β e = e = dr e Free energy profile at temperature T from probability distribution generated under adiabatic conditions: A( λ; β) = kt ln P ( λ; β, β ) λ adb λ
12 Chemical Potential of Lennard-Jones Argon ( ) 4 2 λ = [ λ 1] 4 2 f g( λ) = [( λ 1) 1] m = 2000 m T = 200T λ λ TI
13 Nbins 1 ς () t = P( xi;) t Pexact ( xi) N bins i= 1 [ ]
14 H HO H C H Backbone (Serine) 0.09 H 0.09 H C Backbone (Alanine) Solvation free energies of amino acid side-chain analogs H HO H C H H 0.06 (Methanol) 0.09 H 0.09 H C H 0.09 (Methane) 1 Solute (CHARMm22 Parameters) 256 TIP3P Water molecules Cubic Simulation Box (L = A) Periodic Boundary Conditions Ewald Summation Technique for charges System Temperature: 298 K NVT via GGMT Thermostats (Liu,MET 2000) λ-afed Parameters: kt λ = 12,000 K = 50 kt m λ = 16,000 g.mol-1 a. Wolfenden, et al. Biochem. 20, 849 (1981) b. Shirts, et al. JCP 119, 5740 (2003); JCP 122, (2005) c. Yin and Mackerell J. Comp. Chem. 19, 334 (1998)
15 The Free Energy Profile Probability distribution function: 1 P q = d d e q q Q N N β H( pr, ) ( ) p r δ ( ( r1,..., rn ) ) Free energy profile in q: F( q) = ktln P( q)
16 Bluemoon Ensemble Approach E. A. Carter, et al.chem. Phys. Lett (1989); M. Sprik and G. Ciccotti, J. Chem. Phys. 109, 7737 (1998). Impose a constraint of the form: q( r,..., r ) 1 N = q However, constraints also require: q q q ( r,..., r, p,..., p ) = ir = i = 0 p i 1 N 1 N i i ri i ri mi But in constrained MD, what we are actually computing is: 1 Pq = d p d r e δ qr r q δ q r r p p Q N N β H( pr, ) ( ) ( ( 1,..., N) ) ( ( 1,..., N, 1,..., N))
17 Bluemoon Ensemble df dq = Z 1/2 11 [ λ + kti] Z 1/2 11 Using the Lagrange multiplier to compute free energy: df dq = z 1/2 [ λ + kti] z 1/2 constr constr z = i I= i i m r r z mm r r r r i 2 q q q q q 2 i i i i, j i j i i j j M. Sprik and G. Ciccotti, J. Chem. Phys. 109, 7737 (1998). **When using SHAKE/RATTLE, λ must be the SHAKE multiplier!!
18 d 1 d 2 Free energy profiles A δ = d d 1 2
19 From D. Marx and MET, PRL 86, 4946 (2001).
20 Variable transformations and statistical mechanics
21 Adiabatic dynamics and free energy profiles L. Rosso, P. Minary, Z. Zhu and MET, J. Chem. Phys. 116, 4389 (2000); Maragliano, Vanden Eijnden CPL 426, 168 (2006) Hamiltonian from transformation: n 2 3N 2 pi pi H = + + Uq ( 1,..., q3n ) 2m 2m i= 1 i i= n+ 1 i Adiabatic conditions: m T q k T m k Free energy surface: Aq (,..., q) = kt ln Pq (,..., q) 1 n q 1 n
22
23
24
25
26 Conformational sampling of the solvated alanine dipeptide [L Rosso, J. B. Abrams and MET J. Phys. Chem. B 109, 2099 (2005)] Time FF β α R C7 ax α L φ ψ Meta 1 5ns CHm US 2 400ns CHm AFED 5ns CHm AFED T φ,ψ = 5T, M φ,ψ = 50M C Umbrella Sampling 50 ns 4.7 ns 1. Ensing, et al. ACR 39, 73 (2005) 2. Smith, JCP 111, 5568 (1999) α R C7 ax α L CHARMm22 β
27 NATMA (gas phase) (N-acetyl-Tryptophan-methyl-amide) F(φ,ψ) kt(φ,ψ)=20kt, m(φ,ψ)=600m C t=2.5 ns N H H H CH 3 N ψ N φ O CH 3 O φ C7eq ψ C5(AΦ) Why NATMA? AFED Small Minima compared Predictions: to actual proteins and can be easily studied C5(AP) Ab initio energies 1 Conf. Tryptophan Energy side-chain gives Location rise to a (φ,ψ) free energy Conf. landscape endemic of Rel. actual Energy proteins C5(AP) 0.0 kcal/mol (-160, 160) C5 (AP) 0.0 kcal/mol C5(AΦ) Experimental kcal/mol and DFT 1 Minimization (-140, 140) data available C5(AΦ) kcal/mol C7eq kcal/mol (-100, 100) C7eq kcal/mol 1. Dian, B.C., et al. Science, 296, 2369 (2002); J. Chem. Phys. 117, (2002)
28 Add bias potential to Hamiltonian: Metadynamics A. Laio and M. Parrinello, PNAS 99, (2002); A. Laio, et al. JPCB 109, 6714 (2005) 2 i H = p + U( r1,..., rn) + UG( r1,..., rn, t) 2m i i UG( r,...,, t) W exp G G ( q () q ( () t ) 2 n k k G 1 rn = 2 t= τ,2 τ,... k= 1 2 r r Free energy is negative of bias potential: Aq (,..., q) W exp G G ( q q ( r ( t) ) 2 n 1 n = k k G 2 t= τ,2 τ,... k= 1 2
29
30
31 REPSWA (Reference Potential Spatial Warping Algorithm)
32 V No Transformation Transformation 5kT 10kT 10kT
33 How it works Forces: V =kt V =5kT ( V Vref ) x Fu ( ) = x u = F( x) F ( x) ( ) ref x u x / u becomes large in the barrier region! V =10kT
34 Barrier Crossing Transformations (cont d)
35 V ( φ) ref
36 P. Minary, G. J. Martyna and MET SIAM J. Sci. Comp. (accepted) V ( φ,{}) r = V ( φ) S ( φ) + αv ( r ( φ,{}) r r ) S ( r ( φ,{}) r r ) ref tors 1 inter
37
38
39 Comparison for 50-mer using TraPPE with all interactions PT replicas = 10; PT exchange prob. = 5%, REPSWA α = 0.8; Every 10 th dihedral not transformed
40 Comparison to parallel tempering and CBMC Siepmann and Frenkel, Mol. Phys. 75, 59 (1992) End-to-end distance fluctuations
41 Comparison for 50-mer using CHARMM22 all interactions
42 Comparison of 50-mer using CHARMM22 all interactions
43 Honeycutt and Thirumalai
44 No Transformation Parallel Tempering SDC-REPSWA PT replicas = 16; PT exchange prob. = 5% Model sheet protein β
Enhanced 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 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 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 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 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 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 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 informationA Brief Introduction to Statistical Mechanics
A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade
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 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 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 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 informationEntropy and Free Energy in Biology
Entropy and Free Energy in Biology Energy vs. length from Phillips, Quake. Physics Today. 59:38-43, 2006. kt = 0.6 kcal/mol = 2.5 kj/mol = 25 mev typical protein typical cell Thermal effects = deterministic
More informationAspects of nonautonomous molecular dynamics
Aspects of nonautonomous molecular dynamics IMA, University of Minnesota, Minneapolis January 28, 2007 Michel Cuendet Swiss Institute of Bioinformatics, Lausanne, Switzerland Introduction to the Jarzynski
More informationBiomolecular modeling III
2016, January 5 Déjà vu Enhanced sampling Biomolecular simulation Each atom x, y, z coordinates Déjà vu Enhanced sampling Expression for energy the force field = 1 2 + N i i k i (r i r 0 i ) 2 + 1 2 N
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
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 informationEntropy and Free Energy in Biology
Entropy and Free Energy in Biology Energy vs. length from Phillips, Quake. Physics Today. 59:38-43, 2006. kt = 0.6 kcal/mol = 2.5 kj/mol = 25 mev typical protein typical cell Thermal effects = deterministic
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 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 informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 9 I. OVERVIEW Our treatment of the classical ensembles makes clear that the free energy is a quantity of particular importance in statistical mechanics.
More informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 11 I. PHYSICAL AND CHEMICAL RELEVANCE OF FREE ENERGY In this section, we will consider some examples showing the significance of free energies.
More informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationFree Energy Simulation Methods
Free Energy Simulation Methods Free energy simulation methods Many methods have been developed to compute (relative) free energies on the basis of statistical mechanics Free energy perturbation Thermodynamic
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 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 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 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 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 informationLecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods
Statistical Thermodynamics Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods Dr. Ronald M. Levy ronlevy@temple.edu Free energy calculations Free energy
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 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 informationMolecular Dynamics Lecture 3
Molecular Dynamics Lecture 3 Ben Leimkuhler the problem of the timestep in MD constraints - SHAKE/RATTLE multiple timestepping stochastic methods for holonomic constraints stochastic multiple timestepping
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 informationThermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning
Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning Ben Leimkuhler University of Edinburgh Joint work with C. Matthews (Chicago), G. Stoltz (ENPC-Paris), M. Tretyakov
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 informationApproach to Thermal Equilibrium in Biomolecular
Approach to Thermal Equilibrium in Biomolecular Simulation Eric Barth 1, Ben Leimkuhler 2, and Chris Sweet 2 1 Department of Mathematics Kalamazoo College Kalamazoo, Michigan, USA 49006 2 Centre for Mathematical
More informationStatistical Mechanics. Atomistic view of Materials
Statistical Mechanics Atomistic view of Materials What is statistical mechanics? Microscopic (atoms, electrons, etc.) Statistical mechanics Macroscopic (Thermodynamics) Sample with constrains Fixed thermodynamics
More informationCrossing the barriers - simulations of activated processes
Crossing the barriers - simulations of activated processes Mgr. Ján Hreha for 6 th Student Colloquium and School on Mathematical Physics Faculty of Mathematics, Physics and Informatics Comenius University
More information5th CCPN Matt Crump. Thermodynamic quantities derived from protein dynamics
5th CCPN 2005 -Matt Crump Thermodynamic quantities derived from protein dynamics Relaxation in Liquids (briefly!) The fluctuations of each bond vector can be described in terms of an angular correlation
More informationComputational Predictions of 1-Octanol/Water Partition Coefficient for Imidazolium based Ionic Liquids.
Computational Predictions of 1-Octanol/Water Partition Coefficient for Imidazolium based Ionic Liquids. Ganesh Kamath,* a Navendu Bhatnagar b, Gary A. Baker a, Sheila N. Baker c and Jeffrey J. Potoff b
More informationCharacterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi
Characterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi NCSA Blue Waters Symposium for Petascale Science and Beyond Sunriver, Oregon May 11, 2015 Outline Introduction
More informationalthough Boltzmann used W instead of Ω for the number of available states.
Lecture #13 1 Lecture 13 Obectives: 1. Ensembles: Be able to list the characteristics of the following: (a) icrocanonical (b) Canonical (c) Grand Canonical 2. Be able to use Lagrange s method of undetermined
More informationOrthogonal Space Sampling of Slow Environment Responses
IMA University of Minnesota 2015 Orthogonal Space Sampling of Slow Environment Responses Lianqing Zheng, Chao Lv, Dongsheng Wu, William Harris, Xubin Li, Erick Aitchison, and Wei Yang Institute of Molecular
More informationStatistical Mechanics for Proteins
The Partition Function From Q all relevant thermodynamic properties can be obtained by differentiation of the free energy F: = kt q p E q pd d h T V Q ), ( exp 1! 1 ),, ( 3 3 3 ),, ( ln ),, ( T V Q kt
More informationEnergy and Forces in DFT
Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density
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 informationEXTENDING THE BOUNDARIES OF STRUCTURAL MODELING EPIGENETIC EFFECTS ON NUCLEOSOME POSITIONING
EXTENDING THE BOUNDARIES OF STRUCTURAL MODELING EPIGENETIC EFFECTS ON NUCLEOSOME POSITIONING Peter Minary Computational Structural Biology Group Stanford University Stanford, CA 94305 ~ 0.3 nm TO MODEL
More informationNonequilibrium thermodynamics at the microscale
Nonequilibrium thermodynamics at the microscale Christopher Jarzynski Department of Chemistry and Biochemistry and Institute for Physical Science and Technology ~1 m ~20 nm Work and free energy: a macroscopic
More informationStatistical thermodynamics for MD and MC simulations
Statistical thermodynamics for MD and MC simulations knowing 2 atoms and wishing to know 10 23 of them Marcus Elstner and Tomáš Kubař 22 June 2016 Introduction Thermodynamic properties of molecular systems
More information2. Thermodynamics. Introduction. Understanding Molecular Simulation
2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular
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 informationSupporting Information
Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos
More informationSet the initial conditions r i. Update neighborlist. r i. Get new forces F i
Set the initial conditions r i t 0, v i t 0 Update neighborlist Get new forces F i r i Solve the equations of motion numerically over time step t : r i t n r i t n + v i t n v i t n + Perform T, P scaling
More information(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble
Recall from before: Internal energy (or Entropy): &, *, - (# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble & = /01Ω maximized Ω: fundamental statistical quantity
More informationSampling the free energy surfaces of collective variables
Sampling the free energy surfaces of collective variables Jérôme Hénin Enhanced Sampling and Free-Energy Calculations Urbana, 12 September 2018 Please interrupt! struct bioinform phys chem theoretical
More informationElectronic structure simulations of water solid interfaces
Electronic structure simulations of water solid interfaces Angelos Michaelides London Centre for Nanotechnology & Department of Chemistry, University College London www.chem.ucl.ac.uk/ice Main co-workers:
More informationAdvanced in silico drug design
Advanced in silico drug design RNDr. Martin Lepšík, Ph.D. Lecture: Advanced scoring Palacky University, Olomouc 2016 1 Outline 1. Scoring Definition, Types 2. Physics-based Scoring: Master Equation Terms
More informationThe Dominant Interaction Between Peptide and Urea is Electrostatic in Nature: A Molecular Dynamics Simulation Study
Dror Tobi 1 Ron Elber 1,2 Devarajan Thirumalai 3 1 Department of Biological Chemistry, The Hebrew University, Jerusalem 91904, Israel 2 Department of Computer Science, Cornell University, Ithaca, NY 14853
More informationChE 210B: Advanced Topics in Equilibrium Statistical Mechanics
ChE 210B: Advanced Topics in Equilibrium Statistical Mechanics Glenn Fredrickson Lecture 1 Reading: 3.1-3.5 Chandler, Chapters 1 and 2 McQuarrie This course builds on the elementary concepts of statistical
More informationFree energy calculations
Free energy calculations Jochen Hub & David van der Spoel Overview Free energies and Probabilities Thermodynamic cycles (Free energy perturbation (FEP)) Thermodynamic integration (TI) (Jarzynski equality
More informationSupporting Information. for. Influence of Cononsolvency on the. Aggregation of Tertiary Butyl Alcohol in. Methanol-Water Mixtures
Supporting Information for Influence of Cononsolvency on the Aggregation of Tertiary Butyl Alcohol in Methanol-Water Mixtures Kenji Mochizuki,, Shannon R. Pattenaude, and Dor Ben-Amotz Research Institute
More informationRare Events. Transition state theory Bennett-Chandler Approach 16.2 Transition path sampling16.4
Rare Events Transition state theory 16.1-16.2 Bennett-Chandler Approach 16.2 Transition path sampling16.4 Outline Part 1 Rare event and reaction kinetics Linear Response theory Transition state theory
More informationHow Accurately do Current Force Fields Predict. Experimental Peptide Conformations? An Adiabatic Free Energy Dynamics Study:
How Accurately do Current Force Fields Predict Experimental Peptide Conformations? An Adiabatic Dynamics Study: Supplementary Information (SI) Alexandar T. Tzanov, Michel A. Cuendet,, and Mark E. Tuckerman,,
More informationA self-learning algorithm for biased molecular dynamics
A self-learning algorithm for biased molecular dynamics Tribello, G. A., Ceriotti, M., & Parrinello, M. (2010). A self-learning algorithm for biased molecular dynamics. Proceedings of the National Academy
More informationAlchemical free energy calculations in OpenMM
Alchemical free energy calculations in OpenMM Lee-Ping Wang Stanford Department of Chemistry OpenMM Workshop, Stanford University September 7, 2012 Special thanks to: John Chodera, Morgan Lawrenz Outline
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 informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
More informationIonic Liquids simulations : obtention of structural and transport properties from molecular dynamics. C. J. F. Solano, D. Beljonne, R.
Ionic Liquids simulations : obtention of structural and transport properties from molecular dynamics C. J. F. Solano, D. Beljonne, R. Lazzaroni Ionic Liquids simulations : obtention of structural and transport
More informationAccelerated Quantum Molecular Dynamics
Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended
More informationTheory of electron transfer. Winterschool for Theoretical Chemistry and Spectroscopy Han-sur-Lesse, Belgium, December 2011
Theory of electron transfer Winterschool for Theoretical Chemistry and Spectroscopy Han-sur-Lesse, Belgium, 12-16 December 2011 Electron transfer Electrolyse Battery Anode (oxidation): 2 H2O(l) O2(g) +
More informationLow-temperature isomers of the water hexamer
Low-temperature isomers of the water hexamer Volodymyr Babin and Francesco Paesani Department of Chemistry & Biochemistry, University of California, San Diego, La Jolla, CA 92093-0314 Hundreds of years
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 informationSupplemental Material for Global Langevin model of multidimensional biomolecular dynamics
Supplemental Material for Global Langevin model of multidimensional biomolecular dynamics Norbert Schaudinnus, Benjamin Lickert, Mithun Biswas and Gerhard Stock Biomolecular Dynamics, Institute of Physics,
More informationVariational Implicit Solvation of Biomolecules: From Theory to Numerical Computations
Variational Implicit Solvation of Biomolecules: From Theory to Numerical Computations Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego CECAM Workshop: New Perspectives
More informationIntroduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA
Introduction to Classical Molecular Dynamics Giovanni Chillemi g.chillemi@cineca.it HPC department, CINECA MD ingredients Coordinates Velocities Force field Topology MD Trajectories Input parameters Analysis
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 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 informationLecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability
Lecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions: Van der Waals Interactions
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More informationSupporting Information
Supporting Information ph-responsive self-assembly of polysaccharide through a rugged energy landscape Brian H. Morrow, Gregory F. Payne, and Jana Shen Department of Pharmaceutical Sciences, School of
More informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
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 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 informationCalculating free energies using average force
Center for Turbulence Research Annual Research Briefs 200 27 Calculating free energies using average force By Eric Darve AND Andrew Pohorille. Introduction Many molecular-dynamics computer simulations
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 informationElectronic excitations in conjugated. Many-Body Green's Functions. Behnaz Bagheri Varnousfaderani. Behnaz Bagheri Varnousfaderani
Technische Universiteit Eindhoven University of Technology Electronic excitations in conjugated polymers Electronic from excitations QM/MM simulations in conjugated pol based from on Many-Body QM/MM simulations
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 informationPHYS 352 Homework 2 Solutions
PHYS 352 Homework 2 Solutions Aaron Mowitz (, 2, and 3) and Nachi Stern (4 and 5) Problem The purpose of doing a Legendre transform is to change a function of one or more variables into a function of variables
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 informationSelf-referential Monte Carlo method for calculating the free energy of crystalline solids
Self-referential Monte Carlo method for calculating the free energy of crystalline solids M. B. Sweatman* Department of Chemical and Process Engineering, University of Strathclyde, Glasgow, G1 1XJ, United
More informationClose agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations
Close agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations Alexandre V. Morozov, Tanja Kortemme, Kiril Tsemekhman, David Baker
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationNWChem: Molecular Dynamics and QM/MM
NWChem: Molecular Dynamics and QM/MM Molecular Dynamics Functionality Target systems: General features: biomolecules (proteins, DNA/RNA, biomembranes) energy evaluation (SP) energy minimization (EM) molecular
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 informationEnhanced Free Energy Based Structure Prediction in Materials Science
Enhanced Free Energy Based Structure Prediction in Materials Science Mark E. Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Sciences New York University, 100 Washington Square East,
More informationPrinciples and Applications of Molecular Dynamics Simulations with NAMD
Principles and Applications of Molecular Dynamics Simulations with NAMD Nov. 14, 2016 Computational Microscope NCSA supercomputer JC Gumbart Assistant Professor of Physics Georgia Institute of Technology
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 informationAb initio statistical mechanics
Ab initio statistical mechanics Luca M. Ghiringhelli Fritz-Haber-Institut der MPG, Berlin Hands-on Workshop Density-Functional Theory and Beyond: Accuracy, Efficiency and Reproducibility in Computational
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 information