enhanced sampling methods recovering <a> Accelerated MD Introduction to Accelerated Molecular Dynamics an enhanced sampling method

Size: px
Start display at page:

Download "enhanced sampling methods recovering <a> Accelerated MD Introduction to Accelerated Molecular Dynamics an enhanced sampling method"

Transcription

1 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 the potential of mean force using computer simulations", Comp. Phys. Comm. 91, 1-8, Conformational flooding O. F. Lange, L. V. Schäfer, and H. Grubmüller. Flooding in GROMACS: Accelerated Barrier Crossings in Molecular Dynamics. J. Comput. Chem. 27: (2006) Adaptive biasing force Implementation in NAMD: Research/namd/2.6/ug/node35.html Metadynamics Laio, A.; Parrinello, M. (2002). "Escaping free-energy minima". Proceedings of the National Academy of Sciences of the United States of America 99 (20): and more.. non-boltzmann sampling reweighting Boltzmann distribution Accelerated MD recovering <a> non-boltzmann sampling Boltzmann distribution Donald Hamelberg, John Mongan and J. Andrew McCammon, J. Chem. Phys., 120: (2004).

2 amd modes Boosting the total potential algorithm The algorithm of amd Boosting the dihedral potential only Dual boost: protein->dihedral, rest->total Donald Hamelberg, John Mongan and J. Andrew McCammon, J. Chem. Phys., 120: (2004). Implementing amd in namd Parallelization Strategy Space decomposition Force decomposition Basic Components Patches Compute Objects A typical MD step under the hood The velocity verlet algorithm James C. Phillips, et. al., JCC, 26: (2005). time

3 implementation of amd A typical MD step --- on each patch implementation of amd Total-potential amd A new reduction class Run Compute Objects Report Energy for amd Report/Receive Energy time Y Wang, C.B. Harrison, K. Schulten, J.A. McCammon, 2011 Comput. Sci. Disc Dihedral-only amd Store dihedral force in a new amd force array The amd force array is only initiated and used when amddihe is on. Y Wang, C.B. Harrison, K. Schulten, J.A. McCammon, 2011 Comput. Sci. Disc running amd in namd Benchmark & Output Input parameters ( accelmd Is accelerated molecular dynamics active? Acceptable Values: on or off Default Value: off Description: Specifies if accelerated MD is active. accelmddihe Apply boost to dihedrals? Acceptable Values: on or off Default Value: on Description: Only applies boost to the dihedral potential. By default, accelmddihe is turned on and the boost energy is applied to the dihedral potential of the simulated system. WhenaccelMDdihe is turned off, amd switches to the accelmdt mode, and the boost is applied to the total potential. accelmde Threshold energy Acceptable Values: Real number Description: Specifies the threshold energy in the amd equations. accelmdalpha Acceleration factor Acceptable Values: Positive real number Description: Specifies the acceleration factor in the amd equations. accelmddual Use dual boost mode? Acceptable Values: on or off Default Value: off Description: When accelmddual is on, amd switches to the dual boost mode. Two independent boost potentials will be applied: one to the dihedral potential that is controlled by the parameters accelmde and accelmdalpha, and a second to the (Total - Dihedral) potential that is controlled by the accelmdte and accelmdtalpha parameters described below. amdd: 7% amdt: 12% TCL: Running for steps ACCELERATED MD: STEP 0 dv dvavg BOND ANGL DIHE IMPR ELEC VDW POT PRESSURE: GPRESSURE: ETITLE: TS BOND ANGLE DIHED IMPRP ELECT VDW BOUNDARY MISC KINETIC TOTAL TEMP POTENTIAL TOTAL3 TEMPAVG

4 Choice of parameters an example Two free parameters: E and! Larger E and smaller! lead to higher acceleration levels. As! increases, V*(r) asymptotically approaches V(r).! cannot be zero. E = <V> + N * 0.2! = 1/5* (E - <V>) Alanine dipeptide alanine dipeptide choice of parameters Empirical equation: E = <V> + N * (0.2 to 0.5)! = 1/5* (E - <V>) Classic MD (cmd) (-159, 165) Accelerated MD (amd) (-159, 165) (78, -56) The optimal choice of E/! is systemspecific and may require trial-and-error runs for each system.

5 effect of amd parameters pros & cons Advantages Only two free parameters. amd requires no prior knowledge of the system, i.e., no reaction coordinate needs to be defined a priori. Limitations The exponential form used in reweighting can produce large statistical noise. smaller! Parameter scan can be expensive for large systems. bigger E development of amd More to read Donald Hamelberg, John Mongan and J. Andrew McCammon, JCP, 120: , Replica-exchange amd Mikolai Fajer, Donald Hamelberg and J. Andrew McCammon. JCTC, 4: , Selective boost Donald Hamelberg, César Augusto F. de Oliveira and J. Andrew McCammon, JCP, 127:155102, Mikolai Fajer, Donald Hamelberg and J. Andrew McCammon, JCTC, 4: , Markwick, P., C. Cervantes, B. Abel, E. Komives, M. Blackledge, J.A. McCammon. JACS., 132 (4), , Y Wang, C.B. Harrison, K. Schulten, J.A. McCammon, 2011 Comput. Sci. Disc J. Wereszczynski and J. Andrew McCammon. JCTC, 6: , Adaptive amd Phineus R. L. Markwick, Levi C. T. Pierce, David B. Goodin, J. Andrew McCammon, JPCL, 2: , 2011.

6 acknowledgement San Diego: the McCammon group prof. Andy McCammon Phineus Markwick Cesar de Oliveira Urbana-Champaign: the TCBG group Chris Harrison Jim Phillips Peter Freddolino

w REXAMD: A Hamiltonian Replica Exchange Approach to Improve Free Energy Calculations for Systems with Kinetically Trapped Conformations

w REXAMD: A Hamiltonian Replica Exchange Approach to Improve Free Energy Calculations for Systems with Kinetically Trapped Conformations pubs.acs.org/jctc Downloaded via 148.251.232.83 on March 8, 2019 at 14:33:02 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles. w REXAMD: A Hamiltonian

More information

Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation

Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. pubs.acs.org/jctc Improved

More information

Routine access to millisecond timescale events with accelerated molecular dynamics

Routine access to millisecond timescale events with accelerated molecular dynamics Routine access to millisecond timescale events with accelerated molecular dynamics Levi C.T. Pierce, Romelia Salomon-Ferrer, Cesar Augusto F. de Oliveira #, J. Andrew McCammon #, Ross C. Walker * SUPPORTING

More information

Computing free energies with PLUMED 2.0

Computing 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 information

Protecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations

Protecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations pubs.acs.org/jctc Protecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations William Sinko,*,, Cesar Augusto F. de Oliveira,*,,, Levi C. T. Pierce,

More information

accelerated Molecular Dynamics (amd) Tutorial Levi Pierce 2012 NBCR Summer InsAtute

accelerated Molecular Dynamics (amd) Tutorial Levi Pierce 2012 NBCR Summer InsAtute accelerated Molecular Dynamics (amd) Tutorial Levi Pierce 2012 NBCR Summer InsAtute Time Scales Accessible with Molecular Dynamics Pierce, L.C.T.; Salomon- Ferrer, R.; de Oliveira C.A.; McCammon, J.A.;

More information

Enhancing Amber for Use on the Blue Waters High- Performance Computing Resource

Enhancing Amber for Use on the Blue Waters High- Performance Computing Resource Enhancing Amber for Use on the Blue Waters High- Performance Computing Resource Daniel R. Roe a, Jason M. Swails b, Romelia Salomon-Ferrer c, Adrian E. Roitberg b, and Thomas E. Cheatham III a * a. Department

More information

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 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 information

Performing Metadynamics Simulations Using NAMD

Performing 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 information

Supporting Information

Supporting Information Supporting Information for Activation and Dynamic Network of the M2 Muscarinic Receptor by Yinglong Miao, Sara E. Nichols, Paul M. Gasper, Vincent T. Metzger, and J. Andrew McCammon. Methods Accelerated

More information

Free energy calculations and the potential of mean force

Free 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 information

Using Spectral Clustering to Sample Molecular States and Pathways

Using Spectral Clustering to Sample Molecular States and Pathways Using Spectral Clustering to Sample Molecular States and Pathways Surl-Hee Ahn 1, a) 2, b) and Johannes Birgmeier 1) Chemistry Department, Stanford University, Stanford, California 94305, USA 2) Computer

More information

Computing free energy: Thermodynamic perturbation and beyond

Computing 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 information

Sampling the free energy surfaces of collective variables

Sampling 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 information

Multi-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018

Multi-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 information

Essential dynamics sampling of proteins. Tuorial 6 Neva Bešker

Essential 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 information

Characterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi

Characterizing 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 information

GROMACS Implementation of Metadynamics in Essential Coordinates

GROMACS 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 information

Variational Implicit Solvation of Biomolecules: From Theory to Numerical Computations

Variational 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 information

Determining 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 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 information

Biomolecular modeling III

Biomolecular 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 information

Level-Set Variational Implicit-Solvent Modeling of Biomolecular Solvation

Level-Set Variational Implicit-Solvent Modeling of Biomolecular Solvation Level-Set Variational Implicit-Solvent Modeling of Biomolecular Solvation Bo Li Department of Mathematics and Quantitative Biology Graduate Program UC San Diego The 7 th International Congress of Chinese

More information

Controlling fluctuations

Controlling 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 information

The Molecular Dynamics Method

The Molecular Dynamics Method H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular

More information

ONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole

ONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole ONETEP PB/SA: Application to G-Quadruplex DNA Stability Danny Cole Introduction Historical overview of structure and free energy calculation of complex molecules using molecular mechanics and continuum

More information

Challenges in generation of conformational ensembles for peptides and small proteins

Challenges in generation of conformational ensembles for peptides and small proteins Challenges in generation of conformational ensembles for peptides and small proteins Carlos Simmerling Stony Brook University What could (and does) go wrong? 1. Sampling: difficult to obtain converged

More information

Molecular 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 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 information

Principles and Applications of Molecular Dynamics Simulations with NAMD

Principles 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 information

Thomas E. Cheatham III

Thomas E. Cheatham III NSF OCI-1036208: PRAC Hierarchical molecular dynamics sampling for assessing pathways and free energies of RNA catalysis, ligand binding, and conformational change. NEIS-P2 update, May 2013 Thomas E. Cheatham

More information

Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers

Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers Jaewoon Jung (RIKEN, RIKEN AICS) Yuji Sugita (RIKEN, RIKEN AICS, RIKEN QBiC, RIKEN ithes) Molecular Dynamics

More information

Advanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)

Advanced 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 information

Improved Resolution of Tertiary Structure Elasticity in Muscle Protein

Improved Resolution of Tertiary Structure Elasticity in Muscle Protein Improved Resolution of Tertiary Structure Elasticity in Muscle Protein Jen Hsin and Klaus Schulten* Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois

More information

What makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces: Electronic Supplementary

What 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 information

Local Hyperdynamics. Arthur F. Voter Theoretical Division Los Alamos National Laboratory Los Alamos, NM, USA

Local 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 information

arxiv:physics/ v1 [physics.atm-clus] 21 Jun 2004

arxiv:physics/ v1 [physics.atm-clus] 21 Jun 2004 Conformational properties of neutral and charged alanine and glycine chains arxiv:physics/0406093v1 [physics.atm-clus] 21 Jun 2004 Alexander V. Yakubovitch, Ilia A. Solov yov, Andrey V. Solov yov, and

More information

Computational Chemistry - MD Simulations

Computational 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 information

Molecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains

Molecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains Molecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains Fatemeh Khalili-Araghi, Emad Tajkhorshid, Benoît Roux, and Klaus Schulten Department of Physics, Department of Biochemistry,

More information

Replica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2)

Replica 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 information

The Molecular Dynamics Simulation Process

The Molecular Dynamics Simulation Process The Molecular Dynamics Simulation Process For textbooks see: M.P. Allen and D.J. Tildesley. Computer Simulation of Liquids.Oxford University Press, New York, 1987. D. Frenkel and B. Smit. Understanding

More information

Exploring the Free Energy Surface of Short Peptides by Using Metadynamics

Exploring 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 information

Molecular dynamics (MD) and the Monte Carlo simulation

Molecular 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 information

Level-Set Variational Solvation Coupling Solute Molecular Mechanics with Continuum Solvent

Level-Set Variational Solvation Coupling Solute Molecular Mechanics with Continuum Solvent Level-Set Variational Solvation Coupling Solute Molecular Mechanics with Continuum Solvent Bo Li Department of Mathematics and Center for Theoretical Biological Physics (CTBP) University of California,

More information

Advanced Molecular Dynamics

Advanced 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 information

Analysis of MD Results Using Statistical Mechanics Methods. Molecular Modeling

Analysis of MD Results Using Statistical Mechanics Methods. Molecular Modeling Analysis of MD Results Using Statistical Mechanics Methods Ioan Kosztin eckman Institute University of Illinois at Urbana-Champaign Molecular Modeling. Model building. Molecular Dynamics Simulation 3.

More information

The Mechanism of the Sarco/Endoplasmic Re7culum ATP- Driven Calcium Pump

The Mechanism of the Sarco/Endoplasmic Re7culum ATP- Driven Calcium Pump The Mechanism of the Sarco/Endoplasmic Re7culum ATP- Driven Calcium Pump Blue Waters Symposium Champaign, May 13, 2014 Avisek Das Benoît Roux Department of Biochemistry and Molecular Biology 1 Sarco/endoplasmic

More information

Polymer in Water. Davide Bochicchio and Giovanni M. Pavan* Switzerland, Galleria 2, CH-6928 Manno, Switzerland.

Polymer in Water. Davide Bochicchio and Giovanni M. Pavan* Switzerland, Galleria 2, CH-6928 Manno, Switzerland. Accelerated Atomistic Simulations of a Supramolecular Polymer in Water Davide Bochicchio and Giovanni M. Pavan* Department of Innovative Technologies, University of Applied Sciences and Arts of Southern

More information

Metadynamics with adaptive Gaussians

Metadynamics 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 information

Dielectric Boundary in Biomolecular Solvation Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego

Dielectric Boundary in Biomolecular Solvation Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego Dielectric Boundary in Biomolecular Solvation Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego International Conference on Free Boundary Problems Newton Institute,

More information

Accepted Manuscript. Correcting Mesh-Based Force Calculations to Conserve Both Energy and Momentum in Molecular Dynamics Simulations

Accepted Manuscript. Correcting Mesh-Based Force Calculations to Conserve Both Energy and Momentum in Molecular Dynamics Simulations Accepted Manuscript Correcting Mesh-Based Force Calculations to Conserve Both Energy and Momentum in Molecular Dynamics Simulations Robert D. Seel, David J. Hardy, James C. Phillips PII: S001-1(0)001-1

More information

Ab initio molecular dynamics

Ab 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 information

From Dynamics to Thermodynamics using Molecular Simulation

From Dynamics to Thermodynamics using Molecular Simulation From Dynamics to Thermodynamics using Molecular Simulation David van der Spoel Computational Chemistry Physical models to describe molecules Software to evaluate models and do predictions - GROMACS Model

More information

The Computational Microscope

The Computational Microscope The Computational Microscope Computational microscope views at atomic resolution... Rs SER RER C E M RER N GA L... how living cells maintain health and battle disease John Stone Our Microscope is Made

More information

Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation

Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Bo Li Math Dept & CTBP, UCSD Li-Tien Cheng (Math, UCSD) Zhongming Wang (Math & Biochem, UCSD) Yang Xie (MAE,

More information

Molecular 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 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 information

Replica exchange methodology. John Karanicolas June 2003

Replica exchange methodology. John Karanicolas June 2003 Replica exchange methodology John Karanicolas MMTSB@PSC, June 2003 Outline o Motivation o Theory o Practical considerations o MMTSB Tool Set Why? o Sampling on rugged potential energy surfaces is difficult

More information

Ifs and Buts of DEER. Likai Song, Ilker Sen, Marco Bonora, Florida State University, NHMFL

Ifs and Buts of DEER. Likai Song, Ilker Sen, Marco Bonora, Florida State University, NHMFL Ifs and Buts of DEER Likai Song, Ilker Sen, Marco Bonora, P. Fajer Florida State University, NHMFL 1. X-band v. W-band: is bigger better? DEER STEPR 2. DEER Analysis: how not to over interpret the data?

More information

Dynamics and Calcium Association to the N-Terminal Regulatory Domain of Human Cardiac Troponin C: A Multiscale Computational Study

Dynamics and Calcium Association to the N-Terminal Regulatory Domain of Human Cardiac Troponin C: A Multiscale Computational Study pubs.acs.org/jpcb Dynamics and Calcium Association to the N-Terminal Regulatory Domain of Human Cardiac Troponin C: A Multiscale Computational Study Steffen Lindert,*, Peter M. Kekenes-Huskey, Gary Huber,

More information

A 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 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 information

Monte Carlo (MC) Simulation Methods. Elisa Fadda

Monte 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

Free energy calculations

Free 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 information

Lecture 14: Advanced Conformational Sampling

Lecture 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 information

Modeling 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 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 information

Characterization of the free-energy landscapes of proteins by NMR-guided metadynamics

Characterization 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 information

The role of water and steric constraints in the kinetics of cavity-ligand unbinding

The role of water and steric constraints in the kinetics of cavity-ligand unbinding The role of water and steric constraints in the kinetics of cavity-ligand unbinding However, apart from one recent work [12], to the best of our knowledge there is no reported study in which the timescales

More information

Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006)

Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006) Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006) Martini Workshop 2015 Coarse Graining Basics Alex de Vries Every word or concept, clear as it may seem to be, has only a limited range of applicability

More information

Two-Metal Ion Catalysis by Ribonuclease H

Two-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 information

Electronic structure simulations of water solid interfaces

Electronic 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 information

Semi Empirical Force Fields and Their Limitations. Potential Energy Surface (PES)

Semi Empirical Force Fields and Their Limitations. Potential Energy Surface (PES) Semi Empirical Force Fields and Their Limitations Ioan Kosztin Beckman Institute University of Illinois at Urbana-Champaign Potential Energy Surface (PES) Schrödinger equation: H T Ψ( r, = E Ψ( r, H =

More information

Folding WT villin in silico Three folding simulations reach native state within 5-8 µs

Folding WT villin in silico Three folding simulations reach native state within 5-8 µs Modeling of Cryo-EM Maps Workshop Baylor College of Medicine Klaus Schulten, U. Illinois at Urbana-Champaign Molecular Modeling Flexible Fitting 1: Introduction to Molecular Dynamics MD force field Equation

More information

Carlo Camilloni, Andrea Cavalli,, and Michele Vendruscolo INTRODUCTION

Carlo Camilloni, Andrea Cavalli,, and Michele Vendruscolo INTRODUCTION pubs.acs.org/jpcb Assessment of the Use of NMR Chemical Shifts as Replica-Averaged Structural Restraints in Molecular Dynamics Simulations to Characterize the Dynamics of Proteins Carlo Camilloni, Andrea

More information

Competing sources of variance reduction in parallel replica Monte Carlo, and optimization in the low temperature limit

Competing sources of variance reduction in parallel replica Monte Carlo, and optimization in the low temperature limit Competing sources of variance reduction in parallel replica Monte Carlo, and optimization in the low temperature limit Paul Dupuis Division of Applied Mathematics Brown University IPAM (J. Doll, M. Snarski,

More information

Testing Convergence of Different Free-Energy Methods in a Simple Analytical System with Hidden Barriers

Testing 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 information

Molecular sampling is difficult!

Molecular sampling is difficult! Physics of trajectories and the weighted ensemble method Daniel M. Zuckerman Department of Computational & Systems Biology University of Pittsburgh School of Medicine Molecular sampling is difficult! BPTI

More information

Density Functional Theory: from theory to Applications

Density 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 information

Computational Studies of the Photoreceptor Rhodopsin. Scott E. Feller Wabash College

Computational 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 information

Cy3-DNA Stacking Interactions Strongly Depend on the Identity of the Terminal Basepair

Cy3-DNA Stacking Interactions Strongly Depend on the Identity of the Terminal Basepair Biophysical Journal Volume 100 February 2011 1049 1057 1049 Cy3-DNA Stacking Interactions Strongly Depend on the Identity of the Terminal Basepair Justin Spiriti, Jennifer K. Binder, Marcia Levitus, *

More information

Markov State Models. Gregory R. Bowman Miller Fellow University of California, Berkeley

Markov State Models. Gregory R. Bowman Miller Fellow University of California, Berkeley Markov State Models Gregory R. Bowman Miller Fellow University of California, Berkeley Tuesday, August 16, 2011 Defining The Protein Folding Problem Defining The Protein Folding Problem Why Mechanism?

More information

CO 2 molecule. Morse Potential One of the potentials used to simulate chemical bond is a Morse potential of the following form: O C O

CO 2 molecule. Morse Potential One of the potentials used to simulate chemical bond is a Morse potential of the following form: O C O CO 2 molecule The aim of this project is a numerical analysis of adsorption spectra of CO2 molecule simulated by a double Morse potential function. In the project you should achieve following tasks: 1.

More information

Simulating biomolecular function from motions across multiple scales (I) Peter J. Bond (BII)

Simulating biomolecular function from motions across multiple scales (I) Peter J. Bond (BII) Simulating biomolecular function from motions across multiple scales (I) Peter J. Bond (BII) peterjb@bii.a-star.edu.sg Structural Biology: Why the Need for Simulation? 2017 year 1972 0 RCSB PDB: RCSB Protein

More information

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid Force Fields for Classical Molecular Dynamics simulations of Biomolecules Emad Tajkhorshid Beckman Institute Departments of Biochemistry Center for Biophysics and Computational Biology University of Illinois

More information

CS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C.

CS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C. CS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring 2006 Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C. Latombe Scribe: Neda Nategh How do you update the energy function during the

More information

Probing the Flexibility of Tropomyosin and Its Binding to Filamentous Actin Using Molecular Dynamics Simulations

Probing the Flexibility of Tropomyosin and Its Binding to Filamentous Actin Using Molecular Dynamics Simulations 1882 Biophysical Journal Volume 105 October 2013 1882 1892 Probing the Flexibility of Tropomyosin and Its Binding to Filamentous Actin Using Molecular Dynamics Simulations Wenjun Zheng, * Bipasha Barua,

More information

Statistical Mechanics for Proteins

Statistical 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 information

Dihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769

Dihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769 Dihedral Angles Homayoun Valafar Department of Computer Science and Engineering, USC The precise definition of a dihedral or torsion angle can be found in spatial geometry Angle between to planes Dihedral

More information

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields Department of Chemistry and Biochemistry, Concordia University! page 1 of 6 CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields INTRODUCTION The goal of this tutorial

More information

Introduction to molecular dynamics

Introduction to molecular dynamics 1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.

More information

Equilibrium Molecular Thermodynamics from Kirkwood Sampling

Equilibrium 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 information

Improving Protein Function Prediction with Molecular Dynamics Simulations. Dariya Glazer Russ Altman

Improving Protein Function Prediction with Molecular Dynamics Simulations. Dariya Glazer Russ Altman Improving Protein Function Prediction with Molecular Dynamics Simulations Dariya Glazer Russ Altman Motivation Sometimes the 3D structure doesn t score well for a known function. The experimental structure

More information

Signaling Proteins: Mechanical Force Generation by G-proteins G

Signaling Proteins: Mechanical Force Generation by G-proteins G Signaling Proteins: Mechanical Force Generation by G-proteins G Ioan Kosztin Beckman Institute University of Illinois at Urbana-Champaign Collaborators: Robijn Bruinsma (Leiden & UCLA) Paul O Lague (UCLA)

More information

Equilibrium Molecular Thermodynamics from Kirkwood Sampling

Equilibrium 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 information

Supplemental Material for Global Langevin model of multidimensional biomolecular dynamics

Supplemental 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 information

MOLECULAR DYNAMICS PREDICTS THE SOLUTION CONFORMATIONS OF POLY-L-LYSINE IN SALT SOLUTIONS. by Liqi Feng B.S. in Chemistry, Nankai University, 2011

MOLECULAR DYNAMICS PREDICTS THE SOLUTION CONFORMATIONS OF POLY-L-LYSINE IN SALT SOLUTIONS. by Liqi Feng B.S. in Chemistry, Nankai University, 2011 MOLECULAR DYNAMICS PREDICTS THE SOLUTION CONFORMATIONS OF POLY-L-LYSINE IN SALT SOLUTIONS by Liqi Feng B.S. in Chemistry, Nankai University, 2011 Submitted to the Graduate Faculty of the Kenneth P. Dietrich

More information

XXL-BIOMD. Large Scale Biomolecular Dynamics Simulations. onsdag, 2009 maj 13

XXL-BIOMD. Large Scale Biomolecular Dynamics Simulations. onsdag, 2009 maj 13 XXL-BIOMD Large Scale Biomolecular Dynamics Simulations David van der Spoel, PI Aatto Laaksonen Peter Coveney Siewert-Jan Marrink Mikael Peräkylä Uppsala, Sweden Stockholm, Sweden London, UK Groningen,

More information

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields Department of Chemistry and Biochemistry, Concordia University page 1 of 6 CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields INTRODUCTION The goal of this tutorial

More information

The Force Field Toolkit (fftk)

The Force Field Toolkit (fftk) Parameterizing Small Molecules Using: The Force Field Toolkit (fftk) Christopher G. Mayne, Emad Tajkhorshid Beckman Institute for Advanced Science and Technology University of Illinois, Urbana-Champaign

More information

Quantum Chemistry for Supramolecular Complexes

Quantum Chemistry for Supramolecular Complexes Quantum Chemistry for Supramolecular Complexes Stefan Grimme Mulliken Center for Theoretical Chemistry, University of Bonn COSMO Symposium, March 2015 Supramolecular chemistry: chemistry beyond the covalent

More information

Accelerated Quantum Molecular Dynamics

Accelerated 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 information

Multiscale simulations of complex fluid rheology

Multiscale simulations of complex fluid rheology Multiscale simulations of complex fluid rheology Michael P. Howard, Athanassios Z. Panagiotopoulos Department of Chemical and Biological Engineering, Princeton University Arash Nikoubashman Institute of

More information

Energy 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 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 information

Parallelization of Molecular Dynamics (with focus on Gromacs) SeSE 2014 p.1/29

Parallelization of Molecular Dynamics (with focus on Gromacs) SeSE 2014 p.1/29 Parallelization of Molecular Dynamics (with focus on Gromacs) SeSE 2014 p.1/29 Outline A few words on MD applications and the GROMACS package The main work in an MD simulation Parallelization Stream computing

More information