Computing free energies with PLUMED 2.0
|
|
- Anabel Woods
- 5 years ago
- Views:
Transcription
1 Computing free energies with PLUMED 2.0 Davide Branduardi Formerly at MPI Biophysics, Frankfurt a.m.
2 TyOutline of the talk Relevance of free energy computation Free-energy from histograms: issues Tackling the sampling problem Tackling the order-parameter problem 2
3 TyRelevance of free energy calculations Free energy is a measure of probability Chemical reactions Phase transitions Free Energy find elusive structures and mechanisms! React Drug binding A A Prod s (RC) A(s 0 )= k B T ln P (s 0 ) Biophysics Molecular recognition 3
4 TyRare event: troubles with ergodic hypothesis Free Energy React A stop? Prod 1 2 k BT RC k(t )= k BT hc 0 e A /RT Now I sampled all the two states. Is it sufficient? property Free energy requires statistics! A = 1 ln N 1 N 2 property [3] Eyring, Chem. Rev vol. 17 pp 65 time time 4
5 Ty A test on ATP-Mg 2+ in water ATP 4- -Mg 2+ is of fundamental importance in bioenergetics and regulation Free energy of α β γ to β γ α β γ β γ OW3 O2α OW2 + OW4 O1α OW3 Mg OW4 OW1 + OW2 O2α O1α [1] Liao, Sun, Chandler, Oster Eur. Biophys. J. 2004, vol. 33, p. 29 5
6 TyErgodic hypothesis: am I visiting all the states? Classical MD, CHARMM27 force field, 915 TIP3P water molecules, GROMACS code 6
7 TyPossible solutions Replica exchange methods: multiple simulations at different temperatures or Hamiltonians Biased dynamics: add a potential function along the degree(s) of freedom you want to accelerate Biased dynamics + Replica exchange 7
8 Ty Thermodynamic-Integration-like approaches The ancestor is mean force calculation through harmonic potential F (s) = k B T ln P (s) F (s) s s 0 = k B T 1 R V (x) e k B T (s 0 (x) s)dx s Z e V (x) k B T (s 0 (x) s)dx s 0 Dirac delta function Gaussian function (I like this trick!) F (s) s s 0 ' k B T F (s) s s 0 ' R e V (x) k B T e 1 R e V (x)+k/2(s 0 (x) s 0 ) 2 k B T dx 1 k 2k B T (s0 (x) s) 2 dx s It is an average of a biased simulation! Z k(s 0 (x) Z e V (x) k B T e k 2k B T (s0 (x) s) 2 dx s 0 )e V (x)+k/2(s 0 (x) s0 ) 2 k B T dx s 0 = khs 0 (x) s 0 i bias,s0 (opposite of the mean force ) 8
9 TyThermodynamic Integration 101 Make a constrained simulation (go over barrier!) Acquire mean force Integrate probability descriptor descriptor or +WHAM etc... descriptor 9
10 TyAn example with ATP-Mg Free energy (kcal/mol) O2α -6 Mg Distance (Å) O2α Mg 10
11 TyThermodynamic-Integration-like approaches Thermodynamic integration WHAM Roux, Comp Phys Comm 1995, vol. 91, pp Free energy perturbation Beveridge, DiCapua, Annu Rev Biophys Biophys Chem 1989, vol 18, pp Beveridge, DiCapua, Annu Rev Biophys Biophys Chem 1989, vol 18, pp Jarzynski-equation based approaches (steered-md) Jarzynski Phys Rev Lett 1997, vol. 78, pp Crooks-equation based approaches (two directions steered-md) Crooks J Stat Phys 1998, vol. 90, pp
12 Ty Avoid recoding with MD positions time=t forces time=t evolution time=t+ t Enhanced sampling calculate additional (history dependent) forces Each single engine had only few enhanced sampling techniques (according developers interests) in NAMD, GROMACS, AMBER, LAMMPS, QUANTUM-ESPRESSO Bonomi, Branduardi, Bussi, Camilloni, Provasi, Raiteri, Donadio, Marinelli, Pietrucci, Broglia, Parrinello, M. Comp. Phys. Commun. 2009, 180, Tribello, Bonomi, Branduardi, Camilloni, Bussi Comp. Phys. Commun. (2013) in press. 12
13 TyNow in ESPResSo: 13
14 TyHow to set this up in PLUMED? Mg O2α m(x,t)= k 2 (d(x) x 0(t)) 2 c(x) = SOD X i PHO X j 1 1 NN xi x j d 0 r 0 MM xi x j d 0 r 0 c(x) > 0.2 : U(X) =k (c(x) 0.2) 4 c(x) < 0.2 : U(X) =0. 14
15 TyAn example with ATP-Mg 2+ Free energy (kcal/mol) Hysteresis! Forward Backward O2α -6 Mg Distance (Å) O2α Mg 15
16 TyAdaptive methods: Torrie & Valleau [3] Free energy: Free energy for a biased system: use biasing potential to overcome barriers and retrieve free energies [3] Torrie, Valleau J. Comp. Phys, 1977 vol. 23, pp
17 TyLessons from Torrie & Valleau Free Energy Biasing potentials: React biasing pot 1 2 k BT resulting potential Prod 1 2 k BT RC allow to measure free energies: useful allow to overcome barriers: cheap React Prod need only an additional energy term to standard force field: easy Free Energy Metastabilities are removed RC [3] Torrie, Valleau J. Comp. Phys, 1977 vol. 23, pp
18 TyAdaptive sampling methods Adaptive Biasing Force Darve & Pohorille, J. Chem. Phys. 2001, vol. 115, pp Adaptive Umbrella Sampling Mezei, M. J Comput Phys 1987, vol. 68, pp Self Healing Umbrella Sampling Marsili et. al J. Chem. Phys. B vol. 110, pp Metadynamics Laio & Parrinello, PNAS 2002, vol. 20, pp Conformational Flooding Grubmüller, Phys Rev E 1995, vol. 52, pp
19 Ty Enhanced sampling: metadynamics Apply a repulsive gaussian potential up to compensate the underlying free energy: V g (t, s(x)) = X N hills i=1 We (s(x) s i ) 2 2 eh(t, x) =H(x)+V g (t, s(x)) Produces a non markovian dynamics. The negative of the potential deposited is proven to converge to the correct free energy surface [6] Transition between metastable states are enhanced by artificial removal of metastable states [4] Laio & Parrinello, PNAS 2002, vol. 20, p [5] Bussi & Laio, Parrinello Phys. Rev. Lett. 2006, vol 96, p [6] Branduardi, Bussi, Parrinello, J. Chem. Theory Comput. 2012, vol 8, p
20 TyHow to set this up in PLUMED? V g (t, s(x)) = NX hills i=1 We (s(x) s i ) 2 2 eh(t, x) =H(x)+V g (t, s(x)) d(x) 0.65 c(x) > 0.65 : U(X) = c(x) < 0.65 : U(X) =
21 TyA test on ATP-Mg 2+ in water: meta with distance Choose the CV(s): maximum 3! Free energy estimate React. Prod. AReact.-AProd. 21
22 TyAdaptive sampling vs TI-like methods Adaptive sampling: The problem is multidimensional but probably the accessible phase space is limited: they explore only what you need The free energy landscape has competitive, parallel reactive paths (solvent degrees of freedom, rotations) that can be overcome at some point TI-like: In 1-d (max 2-d, via WHAM or, better with rbf fitting schemes) very effective Sometimes trivially parallel (many umbrellas at same time) 22
23 TyMisconceptions in free energy calculations Technique (TI, WHAM, SMD, MetaD, ABF?) Choice of order parameter 23
24 TyMisconceptions in free energy calculations Technique (TI, WHAM, SMD, MetaD, ABF?) Choice of order parameter 24
25 Ty The importance of the order parameter Different order parameters give different landscapes HighE LowE hidden var random Transition state is mixed pre post Poor convergence of free energy free energy enhanced A CV CV Unreliable transition rates k(t )= k BT hc 0 e A /RT Hard to compare with Expt! [7] Eyring, Chem. Rev vol. 17 p 65 25
26 Ty How to choose the order parameter Different order parameters give different landscapes HighE LowE hidden var Combination of CVs hidden var free energy A CV CV free energy enhanced PRO: the best perspective has good convergence and rates! CONS: it is a trial and error procedure! enhanced A CV CV 26
27 TyPath free energy calculations Monodimensional nonlinear mapping of the transition in terms of many coordinates Progress Distance along from the path hidden var selected path path snapshots in a in given a given representation z(x) s(x) = actual MD snapshot actual MD snapshot P P i ie x x i P 1 PX ln P i e e x x i metrics (cartesian? MSD?) CV i x x i CV CV Adaptive (it uses information from a previous transition) Preserves low dimensional description for efficiency Allows extensions (distance from path) Beyond metadynamics: Steered molecular dynamics and others [8] Branduardi, Gervasio, Parrinello,J. Chem. Phys. 2007, vol.126, p
28 Ty Path collective variables: ATP-Mg 2+ Free energy estimate AReact.-AProd. React. Prod. the shape is conserved, the free energy difference has smaller fluctuations! 28
29 TyPLUMED provides many CVs: are they enough? 29
30 Ty Plumed provides many CVs : combine with functions 30
31 Ty Plumed provides many CVs : matheval c(t) =cos(0.5t 1.028) s(t) =sin(0.5t 1.028) y x U(X,t) = 100 [(s(t) sin(x)) 2 +(c(t) cos(x)) 2 ] 31
32 Ty Plumed provides many CVs : matheval 32
33 TyPlumed is also a set of tools 33
34 TySummary Enhanced sampling method makes you cross barriers Free energies can be calculated from biased dynamics Free energies depends on the technique and from the order parameter Order parameters are as important as the accelerated sampling method PLUMED offers a flexible tool to access a variety of methods and collective variables Beyond what I shown here: mixed replica exchange/biased dynamics (PT-WTE, PT-MetaD, Bias Exchange) Even more: String methods 34
35 TyAcknowledgements Giovanni Bussi (SISSA, Trieste) Max Bonomi (UCSF, San Francisco) Gareth Tribello (Queen s University, Belfast) Carlo Camilloni (Cambridge University) Michele Parrinello (ETH, Zurich) 35
The 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 informationa micro PLUMED2.0 Tutorial for ESPResSo users PLUMED2.0: portable plugin for free-energy calculations with molecular dynamics
a micro PLUMED2.0 Tutorial for ESPResSo users PLUMED2.0: portable plugin for free-energy calculations with molecular dynamics ESPResSo summer school 2013 Organized by CECAM, SimTech, ICP and University
More informationSupporting Information
Supporting Information The effect of strong acid functional groups on electrode rise potential in capacitive mixing by double layer expansion Marta C. Hatzell, a Muralikrishna Raju, b Valerie J. Watson,
More informationa micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics
a micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics FHI-aims Developers and Users Meeting Organized by Fritz Haber Institute, CECAM-mm1p,
More informationa micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics
a micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics Density functional theory and beyond: Computational materials science for real materials
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 informationSupporting Information
1 Supporting Information 2 3 4 Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds 5 6 7 Fabio Pietrucci 1, José C. Aponte 2,3,*, Richard Starr 2,4, Andrea Pérez-Villa
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 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 informationMETAGUI A VMD EXTENSION TO ANALYZE AND VISUALIZE METADYNAMICS SIMULATIONS
METAGUI A VMD EXTENSION TO ANALYZE AND VISUALIZE METADYNAMICS SIMULATIONS Alessandro Laio SISSA & DEMOCRITOS, Trieste Coworkers: Xevi Biarnes Fabio Pietrucci Fabrizio Marinelli Metadynamics (Laio A. and
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 informationA Tunable, Strain-Controlled Nanoporous MoS 2 Filter for Water Desalination
Supporting Information A Tunable, Strain-Controlled Nanoporous MoS 2 Filter for Water Desalination Weifeng Li 1, Yanmei Yang 1, Jeffrey K. Weber 2, Gang Zhang* 3, Ruhong Zhou* 1,2,4 1. School for Radiological
More informationMetadynamics with adaptive Gaussians
Metadynamics with adaptive Gaussians Davide Branduardi Theoretical Molecular Biophysics Group, Max Planck Institute for Biophysics, Max-von-Laue strasse 5, 60438, Frankfurt am Main, Germany Giovanni Bussi
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 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 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 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 informationMaking the best of a bad situation: a multiscale approach to free energy calculation
Making the best of a bad situation: arxiv:1901.04455v2 [physics.comp-ph] 15 Jan 2019 a multiscale approach to free energy calculation Michele Invernizzi,, and Michele Parrinello, Department of Physics,
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 informationQuick reference guide on PLUMED with Quantum ESPRESSO
Quick reference guide on PLUMED with Quantum ESPRESSO Changru Ma SISSA, Trieste March 30, 2011 Contents 1 Introduction 2 1.1 Overview................................... 2 1.2 Collective variables.............................
More informationarxiv: v3 [physics.comp-ph] 30 Apr 2009
PLUMED: a portable plugin for free-energy calculations with molecular dynamics arxiv:0902.0874v3 [physics.comp-ph] 30 Apr 2009 Massimiliano Bonomi, 1 Davide Branduardi, 1 Giovanni Bussi, 2 Carlo Camilloni,
More informationFlexMD Manual. ADF Program System Release 2013
FlexMD Manual ADF Program System Release 2013 Scientific Computing & Modelling NV Vrije Universiteit, Theoretical Chemistry De Boelelaan 1083; 1081 HV Amsterdam; The Netherlands WWW: www.scm.com E-mail:
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 informationenhanced sampling methods recovering <a> Accelerated MD Introduction to Accelerated Molecular Dynamics an enhanced sampling method
enhanced sampling methods Introduction to Accelerated Molecular Dynamics an enhanced sampling method Yi Wang McCammon group Aug 2011, NBCR Summer Institute Umbrella sampling B. ROUX, "The calculation of
More informationPLUMED: a portable plugin for free-energy calculations with molecular dynamics
PLUMED: a portable plugin for free-energy calculations with molecular dynamics Massimiliano Bonomi,a, Davide Branduardi a, Giovanni Bussi b,carlo Camilloni c,davideprovasi d, Paolo Raiteri e, Davide Donadio
More informationCharacterization of the free-energy landscapes of proteins by NMR-guided metadynamics
Characterization of the free-energy landscapes of proteins by NMR-guided metadynamics B Results and Discussion BIOPHYSICS ND COMPUTTIONL BIOLOGY SI Text Inset B χ χ Inset SI Text Inset C Folding of GB3
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 Jochen Hub & David van der Spoel Overview Free energies and Probabilities Thermodynamic cycles (Free energy perturbation (FEP)) Thermodynamic integration (TI) (Jarzynski equality
More informationA self-learning algorithm for biased molecular dynamics
A self-learning algorithm for biased molecular dynamics Gareth A. Tribello, Michele Ceriotti and Michele Parrinello Computational Science, Department of Chemistry and Applied Biosciences, ETHZ Zurich USI-Campus,
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 informationA general overview of. Free energy methods. Comer et al. (2014) JPCB 119:1129. James C. (JC) Gumbart Georgia Institute of Technology, Atlanta
A general overview of Free energy methods Comer et al. (2014) JPCB 119:1129. James C. (JC) Gumbart Georgia Institute of Technology, Atlanta Computational Biophysics Workshop Georgia Tech Nov. 16 2016 Outline
More informationTesting Convergence of Different Free-Energy Methods in a Simple Analytical System with Hidden Barriers
computation Article Testing Convergence of Different Free-Energy Methods in a Simple Analytical System with Hidden Barriers S. Alexis Paz 1, ID and Cameron F. Abrams 3, * ID 1 Departamento de Química Teórica
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 informationarxiv: v2 [physics.comp-ph] 4 Feb 2016
Sampling Free Energy Surfaces as Slices by Combining Umbrella Sampling and Metadynamics Shalini Awasthi, Venkat Kapil and Nisanth N. Nair 1 arxiv:1507.06764v2 [physics.comp-ph] 4 Feb 2016 Department of
More informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
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 informationEncarna Pucheta-Martinez+, Nicola D Amelio+, Moreno Lelli, Jorge L. Martinez-Torrecuadrada, Marius Sudol, Giorgio Saladino* and Francesco L.
Supplementary Information for: Changes in the folding landscape of the WW domain provide a molecular mechanism for an inherited genetic syndrome Encarna Pucheta-Martinez+, Nicola D Amelio+, Moreno Lelli,
More informationDetermining Energy Barriers and Selectivities of a Multi-Pathway System With Infrequent Metadynamics
Determining Energy Barriers and Selectivities of a Multi-Pathway System With Infrequent Metadynamics Christopher D. Fu 1, Luiz F.L. Oliveira 1 and Jim Pfaendtner 1,a) 1 Department of Chemical Engineering,
More informationarxiv: v1 [physics.comp-ph] 12 Dec 2014
Accurate multiple time step in biased molecular simulations Marco Jacopo Ferrarotti, Sandro Bottaro, Andrea Pérez-Villa, and Giovanni Bussi Scuola Internazionale Superiore di Studi Avanzati (SISSA), via
More informationBacterial Outer Membrane Porins as Electrostatic Nanosieves: Exploring Transport Rules of Small Polar Molecules
Bacterial Outer Membrane Porins as Electrostatic Nanosieves: Exploring Transport Rules of Small Polar Molecules Harsha Bajaj, Silvia Acosta Gutiérrez, Igor Bodrenko, Giuliano Malloci, Mariano Andrea Scorciapino,
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 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 informationFree-energy calculations
Measuring free-energy differences using computer simulations Theoretical and Computational Biophysics Group University of Illinois, Urbana Champaign and Équipe de dynamique des assemblages membranaires,
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 informationHow wet should be the reaction coordinate for ligand unbinding?
THE JOURNAL OF CHEMICAL PHYSICS 145, 054113 (2016) How wet should be the reaction coordinate for ligand unbinding? Pratyush Tiwary a) and B. J. Berne b) Department of Chemistry, Columbia University, New
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 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 informationGROMETA 2.0. Carlo Camilloni and Davide Provasi. Doing metadynamics with GROMACS
Carlo Camilloni and Davide Provasi GROMETA 2.0 Doing metadynamics with GROMACS May 8, 2008 Protein Physics Research Group Physics Department, University of Milano via Celoria, 16-20133 Milano - Italy Table
More informationMARTINI simulation details
S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained
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 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 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 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 informationSupporting Information for. Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite. Edge Water Interface
Supporting Information for Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite Edge Water Interface Zhen Jiang, Konstantin Klyukin, and Vitaly Alexandrov,, Department
More informationComparison of Sampling Methods via Robust. Free Energy Inference: Application to. Calmodulin
Comparison of Sampling Methods via Robust Free Energy Inference: Application to arxiv:1704.00343v1 [q-bio.bm] 2 Apr 2017 Calmodulin Annie M. Westerlund, Tyler J. Harpole, Christian Blau,, and Lucie Delemotte,
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 informationMathematical analysis of an Adaptive Biasing Potential method for diffusions
Mathematical analysis of an Adaptive Biasing Potential method for diffusions Charles-Edouard Bréhier Joint work with Michel Benaïm (Neuchâtel, Switzerland) CNRS & Université Lyon 1, Institut Camille Jordan
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 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 informationEnhanced sampling of transition states
Enhanced sampling of transition states arxiv:1812.09032v1 [physics.chem-ph] 21 Dec 2018 Jayashrita Debnath,, Michele Invernizzi,, and Michele Parrinello,,, Department of Chemistry and Applied Biosciences,
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 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 informationGROMACS Implementation of Metadynamics in Essential Coordinates
GROMACS implementation of metadynamics in essential coordinates 1 GROMACS Implementation of Metadynamics in Essential Coordinates Contact: Vojtěch Spiwok Centre for Glycomics Chemical Department, Slovak
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 informationProbing α/β Balances in Modified Amber Force Fields from a Molecular Dynamics Study on a ββα Model Protein (1FSD)
Benchmark Study of Modified Amber Force Fields on ββα Protein Bull. Korean Chem. Soc. 214, Vol. 35, No. 6 1713 http://dx.doi.org/1.512/bkcs.214.35.6.1713 Probing α/β Balances in Modified Amber Force Fields
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 2 Feb 2005
Topological defects and bulk melting of hexagonal ice Davide Donadio, Paolo Raiteri, and Michele Parrinello Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus,
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 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 informationMetadynamics. Day 2, Lecture 3 James Dama
Metadynamics Day 2, Lecture 3 James Dama Metadynamics The bare bones of metadynamics Bias away from previously visited configuraaons In a reduced space of collecave variables At a sequenaally decreasing
More informationBayesian reconstruction of free energy profiles from umbrella samples
Bayesian reconstruction of free energy profiles from umbrella samples Thomas Stecher with Gábor Csányi and Noam Bernstein 21st November 212 Motivation Motivation and Goals Many conventional methods implicitly
More informationarxiv: v1 [cond-mat.stat-mech] 27 Nov 2015
The solid-liquid interfacial free energy out of equilibrium arxiv:1511.8668v1 [cond-mat.stat-mech] 27 Nov 215 Bingqing Cheng, 1 Gareth A. Tribello, 2 and Michele Ceriotti 1, 1 Laboratory of Computational
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 informationMetashooting: A Novel Tool for Free Energy Reconstruction from Polymorphic Phase Transition Mechanisms. Abstract
Metashooting: A Novel Tool for Free Energy Reconstruction from Polymorphic Phase Transition Mechanisms Samuel Alexander Jobbins, a Salah Eddine Boulfelfel, b and Stefano Leoni a a School of Chemistry,
More informationInsights into Ligand Protein Binding from Local Mechanical Response
pubs.acs.org/jctc Insights into Ligand Protein Binding from Local Mechanical Response Jagdish Suresh Patel, Davide Branduardi,*, Matteo Masetti,, Walter Rocchia, and Andrea Cavalli*,, Department of Drug
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 informationA Parallel-pulling Protocol for Free Energy Evaluation INTRODUCTION B A) , p
A Parallel-pulling Protocol for Free Energy Evaluation Van Ngo Collaboratory for Advance Computing and Simulations (CACS, Department of Physics and Astronomy, University of Southern California, 3651 Watt
More informationarxiv: v1 [cond-mat.stat-mech] 7 Dec 2017
Learning Free Energy Landscapes Using Artificial Neural Networks arxiv:72.02840v [cond-mat.stat-mech] 7 Dec 207 Hythem Sidky and Jonathan K. Whitmer, Department of Chemical and Biomolecular Engineering,
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 informationAdvanced Molecular Dynamics
Advanced Molecular Dynamics Introduction May 2, 2017 Who am I? I am an associate professor at Theoretical Physics Topics I work on: Algorithms for (parallel) molecular simulations including GPU acceleration
More informationProbing the Mechanism of ph-induced Large-Scale Conformational Changes in Dengue Virus Envelope Protein Using Atomistic Simulations
588 Biophysical Journal Volume 99 July 21 588 594 Probing the Mechanism of ph-induced Large-Scale Conformational Changes in Dengue Virus Envelope Protein Using Atomistic Simulations Meher K. Prakash,*
More informationPerforming 2D Hamiltonian-Exchange method for Biological Free Energy Calculations on Petascale Supercomputer
Performing D Hamiltonian-Echange method for Biological Free Energy Calculations on Petascale Supercomputer Yun Lyna Luo, Wei Jiang, Benoît Rou* Argonne Leadership Computing Facility, Argonne National Laboratory
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More informationFree energy calculation from steered molecular dynamics simulations using Jarzynski s equality
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 6 8 AUGUST 2003 Free energy calculation from steered molecular dynamics simulations using Jarzynski s equality Sanghyun Park and Fatemeh Khalili-Araghi Beckman
More informationSANDRO BOTTARO, PAVEL BANÁŠ, JIŘÍ ŠPONER, AND GIOVANNI BUSSI
FREE ENERGY LANDSCAPE OF GAGA AND UUCG RNA TETRALOOPS. SANDRO BOTTARO, PAVEL BANÁŠ, JIŘÍ ŠPONER, AND GIOVANNI BUSSI Contents 1. Supplementary Text 1 2 Sample PLUMED input file. 2. Supplementary Figure
More informationMetadynamics-biased ab initio Molecular Dynamics Study of
Metadynamics-biased ab initio Molecular Dynamics Study of Heterogeneous CO 2 Reduction via Surface Frustrated Lewis Pairs Mireille Ghoussoub 1, Shwetank Yadav 1, Kulbir Kaur Ghuman 1, Geoffrey A. Ozin
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 informationIntroduction to Numerical Simulations and High Performance Computing: From Materials Science to Biochemistry
Introduction to Numerical Simulations and High Performance Computing: From Materials Science to Biochemistry Mauro Boero Institut de Physique et Chimie des Matériaux de Strasbourg University of Strasbourg
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 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 informationTwo-Metal Ion Catalysis by Ribonuclease H
Two-Metal Ion Catalysis by Ribonuclease H Edina Rosta Department of Chemistry King s College London Phosphate Groups as Building Blocks Biological importance: Reproduction: DNA and RNA hydrolysis, synthesis
More informationMulticanonical parallel tempering
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 13 1 APRIL 2002 Multicanonical parallel tempering Roland Faller, Qiliang Yan, and Juan J. de Pablo Department of Chemical Engineering, University of Wisconsin,
More informationElectronic Supporting Information Topological design of porous organic molecules
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2017 Electronic Supporting Information Topological design of porous organic molecules Valentina Santolini,
More informationRepresenting High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks
Energy Landscapes, Chemnitz 200 Representing High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks Jörg Behler Lehrstuhl für Theoretische Chemie Ruhr-Universität Bochum D-44780 Bochum,
More informationAu-C Au-Au. g(r) r/a. Supplementary Figures
g(r) Supplementary Figures 60 50 40 30 20 10 0 Au-C Au-Au 2 4 r/a 6 8 Supplementary Figure 1 Radial bond distributions for Au-C and Au-Au bond. The zero density regime between the first two peaks in g
More informationSupporting Information
Prebiotic NH 3 formation: Insights from simulations András Stirling 1, Tamás Rozgonyi 2, Matthias Krack 3, Marco Bernasconi 4 1 Institute of Organic Chemistry, Research Centre for Natural Sciences of the
More informationLarge Scale Electronic Structure Calculations
Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)
More informationOptimal Thermodynamic Control and the Riemannian Geometry of Ising magnets
Optimal Thermodynamic Control and the Riemannian Geometry of Ising magnets Gavin Crooks Lawrence Berkeley National Lab Funding: Citizens Like You! MURI threeplusone.com PRE 92, 060102(R) (2015) NSF, DOE
More informationGenerating, 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 informationIs there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich
Is there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich Chemistry Chemistry is the science of atomic matter, especially its chemical reactions,
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 information