Limitations of temperature replica exchange (T-REMD) for protein folding simulations

Size: px
Start display at page:

Download "Limitations of temperature replica exchange (T-REMD) for protein folding simulations"

Transcription

1 Limitations of temperature replica exchange (T-REMD) for protein folding simulations Jed W. Pitera, William C. Swope IBM Research

2 Anomalies in protein folding kinetic thermodynamic 322K 305K Yang & Gruebele, Nature Ma & Gruebele, PNAS 2005 ~320K ~325K Scaled SVD Fl. data, 63 C Fl. fit IR data, 63 C IR fit 0.0 Garcia-Mira et al,science Time (µs) 150

3 Hunting anomalies with a computational model Do we have the correct/converged answer for our model? How do we know? Do we have a model that reflects reality? How do we compare the model against experiment? What is missing from our model?

4 Replica exchange molecular dynamics (REMD) 400K 350K 300K Multiple simulations of the same system are run in parallel at different temperatures (T-REMD), state points or Hamiltonians Monte Carlo moves periodically exchange systems between adjacent temperatures Allows escape from local minima Ideal for cases where we want temperature-dependent properties 350K ensemble { }

5 T-REMD motivations T high realm of perfect sampling Interested in a system at a temperature/state point where sampling is slow (T low ) Long correlation times Broken ergodicity Assume that sampling is fast at some other temperature/state point (T high ) Simulate as many intermediate state points as necessary to bridge T low and T high Trajectories (MD) or Markov Chains (MC) decorrelate at T high, importance sample at T low In many cases T was an interesting variable anyway Experiments often provide A(T) or perturb a system by T T-dependent phenomena of biological interest T low broken ergodicity

6 Model system T-REMD 1-D double well potential, compare MD and T-REMD Energetic barrier; activated process with Arrhenius kinetics (ln(k) linear in 1/T) ln k -6 MD MD transitions MD aggregate transitions in 2.5x10^7 steps vs MD (4.5k/4.5k) 1/T replica Similar rates vs. 1/T All replicas undergo transitions

7 trpzip2 β-hairpin in explicit solvent 12 amino acid peptide trpzip2 in explicit solvent TIP3P, AMBER parm96 or parm99sb 3605 waters, atoms Cubic box (equil.10ns 310K, 1 atm) edge length Å PME electrostatics 9Å Switch for vdw/direct; long range vdw correction Replica exchange molecular dynamics (80 replicas at a range of temperatures) Exchanges every 40 ps; Andersen collisions every 10 ps 2 independent calculations with different initial conditions 80 representative conformations from implicit solvent folded (0.68 µs/replica, aggregate 54 µs) All 80 replicas started in the same fully extended conformation unfolded (1.45 µs/replica, aggregate 116 µs)

8 Effect of exchange period on relaxation from the folded state

9 The hazards of non-thermalized initial conditions Rapid unfolding from folded initial conditions Exchange period shorter than the relaxation time of the potential energy

10 Stability of folded initial conditions

11 Convergence from unfolded initial conditions

12 SASA (Å 2 ) IBM Research trpzip2 thermodynamics SASA PMF vs. T Continuous, weak collapse transition fr from folded temperature (K) from unfolded

13 Cα RMSD (Å) IBM Research trpzip2 thermodynamics CαRMSD PMF vs. T Measure of the backbone deviation from the NMR structure from folded temperature (K) from unfolded

14 heavy atom RMSD (Å) IBM Research trpzip2 thermodynamics heavy atom RMSD PMF vs. T Spurious absence of a barrier from folded temperature (K) from unfolded

15 Comparison of cluster populations Approximate stochastic k-medoid clustering of merged & downsampled data set to produce a set of 40 clusters; metric was distance matrix error of C α, trp C δ and C ζ3

16 The energy landscape Markov model of 425K kinetics (N. Singhal) Compact states are isolated local minima connected by unfolded state Kinetics in/out of these minima are slow

17 Apparent folding rates RMSD criteria folded: < 2.5 Å C α - RMSD from NMR unfolded: > 6 Å C α - RMSD from NMR Track perreplica transitions, record T Exp tl k f 5x10-7, k u 5x10-8 ps 296 K (Snow et al PNAS 2004)

18 Apparent folding rates cluster membership Transitions to/from cluster #1 Order of magnitude difference in rates Different T- dependence

19 Apparent folding rates nonequilibrium data Successive block averages of same data set Started folded, parm99sb Systematic ~2x change in unfolding rate

20 Conclusions T-REMD is useful but not a panacea No increase in aggregate # of transitions Many interesting barriers entropic rather than energetic No T where sampling is infinitely fast Explicit solvent REMD of proteins has limitations Decoupling of D.O.F. of interest (protein) from extended variables (T, U, etc.) Large N small T, limiting replica motion in T Sampling limited by intra-replica correlation time Folding not a simple activated process

21 Acknowledgements William Swope, Hans Horn, Julia Rice (ARC) Robert Germain (YKT), Blue Gene Science & Application Team Martin Gruebele & Wei Yang (UIUC) Vijay Pande, Nina Singhal, Michael Shirts (Stanford) John Chodera & Ken Dill (UCSF)

22 Gordon Research Conference in Computational Chemistry, July 27 Aug Mount Holyoke, MA Chair: Dr. Jed W. Pitera, IBM Research Vice Chair: Prof. Dr. Walter Thiel, MPI-Kohlenforschung Force fields, electronic structure, quantum dynamics, chemical reactions, drug design, docking, coarse-grained simulation accuracy time and space

Convergence of replica exchange molecular dynamics

Convergence of replica exchange molecular dynamics THE JOURNAL OF CHEMICAL PHYSICS 123, 154105 2005 Convergence of replica exchange molecular dynamics Wei Zhang and Chun Wu Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware

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

Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?

Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Ruhong Zhou 1 and Bruce J. Berne 2 1 IBM Thomas J. Watson Research Center; and 2 Department of Chemistry,

More information

Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation

Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation Jakob P. Ulmschneider and William L. Jorgensen J.A.C.S. 2004, 126, 1849-1857 Presented by Laura L. Thomas and

More information

Many proteins spontaneously refold into native form in vitro with high fidelity and high speed.

Many proteins spontaneously refold into native form in vitro with high fidelity and high speed. Macromolecular Processes 20. Protein Folding Composed of 50 500 amino acids linked in 1D sequence by the polypeptide backbone The amino acid physical and chemical properties of the 20 amino acids dictate

More information

Biology Chemistry & Physics of Biomolecules. Examination #1. Proteins Module. September 29, Answer Key

Biology Chemistry & Physics of Biomolecules. Examination #1. Proteins Module. September 29, Answer Key Biology 5357 Chemistry & Physics of Biomolecules Examination #1 Proteins Module September 29, 2017 Answer Key Question 1 (A) (5 points) Structure (b) is more common, as it contains the shorter connection

More information

Application of the Markov State Model to Molecular Dynamics of Biological Molecules. Speaker: Xun Sang-Ni Supervisor: Prof. Wu Dr.

Application of the Markov State Model to Molecular Dynamics of Biological Molecules. Speaker: Xun Sang-Ni Supervisor: Prof. Wu Dr. Application of the Markov State Model to Molecular Dynamics of Biological Molecules Speaker: Xun Sang-Ni Supervisor: Prof. Wu Dr. Jiang Introduction Conformational changes of proteins : essential part

More information

Supporting Information. RNA hairpin

Supporting Information. RNA hairpin Supporting Information RNA hairpin 1. Initial Configurations We started our ST simulations from two different initial configurations as shown in Figure S1: a near-native state and a random coil. The near-native

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

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

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

Energetics and Thermodynamics

Energetics and Thermodynamics DNA/Protein structure function analysis and prediction Protein Folding and energetics: Introduction to folding Folding and flexibility (Ch. 6) Energetics and Thermodynamics 1 Active protein conformation

More information

Outline. The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation. Unfolded Folded. What is protein folding?

Outline. The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation. Unfolded Folded. What is protein folding? The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation By Jun Shimada and Eugine Shaknovich Bill Hawse Dr. Bahar Elisa Sandvik and Mehrdad Safavian Outline Background on protein

More information

Recovering Kinetics from a Simplified Protein Folding Model Using Replica Exchange Simulations: A Kinetic Network and Effective Stochastic Dynamics

Recovering Kinetics from a Simplified Protein Folding Model Using Replica Exchange Simulations: A Kinetic Network and Effective Stochastic Dynamics 11702 J. Phys. Chem. B 2009, 113, 11702 11709 Recovering Kinetics from a Simplified Protein Folding Model Using Replica Exchange Simulations: A Kinetic Network and Effective Stochastic Dynamics Weihua

More information

Free Energy Landscape of Protein Folding in Water: Explicit vs. Implicit Solvent

Free Energy Landscape of Protein Folding in Water: Explicit vs. Implicit Solvent PROTEINS: Structure, Function, and Genetics 53:148 161 (2003) Free Energy Landscape of Protein Folding in Water: Explicit vs. Implicit Solvent Ruhong Zhou* IBM T.J. Watson Research Center, Yorktown Heights,

More information

Hydrophobic Aided Replica Exchange: an Efficient Algorithm for Protein Folding in Explicit Solvent

Hydrophobic Aided Replica Exchange: an Efficient Algorithm for Protein Folding in Explicit Solvent 19018 J. Phys. Chem. B 2006, 110, 19018-19022 Hydrophobic Aided Replica Exchange: an Efficient Algorithm for Protein Folding in Explicit Solvent Pu Liu, Xuhui Huang, Ruhong Zhou,, and B. J. Berne*,, Department

More information

arxiv: v1 [cond-mat.soft] 22 Oct 2007

arxiv: v1 [cond-mat.soft] 22 Oct 2007 Conformational Transitions of Heteropolymers arxiv:0710.4095v1 [cond-mat.soft] 22 Oct 2007 Michael Bachmann and Wolfhard Janke Institut für Theoretische Physik, Universität Leipzig, Augustusplatz 10/11,

More information

Molecular Interactions F14NMI. Lecture 4: worked answers to practice questions

Molecular Interactions F14NMI. Lecture 4: worked answers to practice questions Molecular Interactions F14NMI Lecture 4: worked answers to practice questions http://comp.chem.nottingham.ac.uk/teaching/f14nmi jonathan.hirst@nottingham.ac.uk (1) (a) Describe the Monte Carlo algorithm

More information

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

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

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

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

Alchemical free energy calculations in OpenMM

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

Short Announcements. 1 st Quiz today: 15 minutes. Homework 3: Due next Wednesday.

Short Announcements. 1 st Quiz today: 15 minutes. Homework 3: Due next Wednesday. Short Announcements 1 st Quiz today: 15 minutes Homework 3: Due next Wednesday. Next Lecture, on Visualizing Molecular Dynamics (VMD) by Klaus Schulten Today s Lecture: Protein Folding, Misfolding, Aggregation

More information

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

Analysis of MD trajectories in GROMACS David van der Spoel

Analysis of MD trajectories in GROMACS David van der Spoel Analysis of MD trajectories in GROMACS David van der Spoel What does MD produce? Energy terms E(t) Coordinates x(t) Velocities v(t) Forces f(t) Managing your files trjcat - merging trajectories concatenating

More information

Introduction The gramicidin A (ga) channel forms by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet. Th

Introduction The gramicidin A (ga) channel forms by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet. Th Abstract When conductive, gramicidin monomers are linked by six hydrogen bonds. To understand the details of dissociation and how the channel transits from a state with 6H bonds to ones with 4H bonds or

More information

Choosing weights for simulated tempering

Choosing weights for simulated tempering PHYSICAL REVIEW E 7, 173 7 Choosing weights for simulated tempering Sanghyun Park 1 and Vijay S. Pande 1, 1 Department of Chemistry, Stanford University, Stanford, California 935, USA Department of Structural

More information

Free energy calculations using molecular dynamics simulations. Anna Johansson

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

2008 Biowerkzeug Ltd.

2008 Biowerkzeug Ltd. 2008 Biowerkzeug Ltd. 1 Contents Summary...3 1 Simulation...4 1.1 Setup...4 1.2 Output...4 2 Settings...5 3 Analysis...9 3.1 Setup...9 3.2 Input options...9 3.3 Descriptions...10 Please note that we cannot

More information

The computational Core of Molecular Modeling. (The What? Why? And How? )

The computational Core of Molecular Modeling. (The What? Why? And How? ) The computational Core of Molecular Modeling. (The What? Why? And How? ) CS Seminar. Feb 09. Alexey Onufriev, Dept of Computer Science; Dept. of Physics and the GBCB program. Thanks to: Faculty co-authors:

More information

Protein Folding experiments and theory

Protein Folding experiments and theory Protein Folding experiments and theory 1, 2,and 3 Protein Structure Fig. 3-16 from Lehninger Biochemistry, 4 th ed. The 3D structure is not encoded at the single aa level Hydrogen Bonding Shared H atom

More information

Free energy simulations

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

PROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS

PROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS TASKQUARTERLYvol.20,No4,2016,pp.353 360 PROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS MARTIN ZACHARIAS Physics Department T38, Technical University of Munich James-Franck-Str.

More information

Energy Barriers and Rates - Transition State Theory for Physicists

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

What would happen? + -

What would happen? + - What would happen? + - Lecture Plan Molecules spend more time in low energy conformations Predicting a structure is approximated as the search for the lowest energy conformation on a protein model Modeling

More information

Protein Folding Prof. Eugene Shakhnovich

Protein Folding Prof. Eugene Shakhnovich Protein Folding Eugene Shakhnovich Department of Chemistry and Chemical Biology Harvard University 1 Proteins are folded on various scales As of now we know hundreds of thousands of sequences (Swissprot)

More information

THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION

THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION AND CALIBRATION Calculation of turn and beta intrinsic propensities. A statistical analysis of a protein structure

More information

HOW WELL CAN SIMULATION PREDICT PROTEIN FOLDING KINETICS AND THERMODYNAMICS?

HOW WELL CAN SIMULATION PREDICT PROTEIN FOLDING KINETICS AND THERMODYNAMICS? Annu. Rev. Biophys. Biomol. Struct. 2005. 34:43 69 doi: 10.1146/annurev.biophys.34.040204.144447 Copyright c 2005 by Annual Reviews. All rights reserved First published online as a Review in Advance on

More information

Molecular Dynamics. A very brief introduction

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

PHYSICAL REVIEW LETTERS

PHYSICAL REVIEW LETTERS PHYSICAL REVIEW LETTERS VOLUME 86 28 MAY 21 NUMBER 22 Mathematical Analysis of Coupled Parallel Simulations Michael R. Shirts and Vijay S. Pande Department of Chemistry, Stanford University, Stanford,

More information

Fibril Elongation by Aβ : Kinetic Network Analysis of Hybrid- Resolution Molecular Dynamics Simulations

Fibril Elongation by Aβ : Kinetic Network Analysis of Hybrid- Resolution Molecular Dynamics Simulations pubs.acs.org/jacs Terms of Use Fibril Elongation by Aβ 17 42 : Kinetic Network Analysis of Hybrid- Resolution Molecular Dynamics Simulations Wei Han, and Klaus Schulten*,,, Beckman Institute, Center for

More information

3.320 Lecture 23 (5/3/05)

3.320 Lecture 23 (5/3/05) 3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials

More information

Computer simulation methods (2) Dr. Vania Calandrini

Computer simulation methods (2) Dr. Vania Calandrini Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing

More information

Hierarchical Metastable States and Kinetic Transition Networks: Trajectory Mapping and Clustering

Hierarchical Metastable States and Kinetic Transition Networks: Trajectory Mapping and Clustering Hierarchical Metastable States and Kinetic Transition Networks: Trajectory Mapping and Clustering Xin Zhou xzhou@gucas.ac.cn Graduate University of Chinese Academy of Sciences, Beijing 2012.6.5 KITP, Santa

More information

Orthogonal Space Sampling of Slow Environment Responses

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

Transition Paths of Met-Enkephalin from Markov State Modeling of a Molecular Dynamics Trajectory

Transition Paths of Met-Enkephalin from Markov State Modeling of a Molecular Dynamics Trajectory pubs.acs.org/jpcb Transition Paths of Met-Enkephalin from Markov State Modeling of a Molecular Dynamics Trajectory Rahul Banerjee and Robert I. Cukier* Department of Chemistry, Michigan State University,

More information

Molecular Mechanics. I. Quantum mechanical treatment of molecular systems

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

Methods of Computer Simulation. Molecular Dynamics and Monte Carlo

Methods of Computer Simulation. Molecular Dynamics and Monte Carlo Molecular Dynamics Time is of the essence in biological processes therefore how do we understand time-dependent processes at the molecular level? How do we do this experimentally? How do we do this computationally?

More information

Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning

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

Biomolecules are dynamic no single structure is a perfect model

Biomolecules are dynamic no single structure is a perfect model Molecular Dynamics Simulations of Biomolecules References: A. R. Leach Molecular Modeling Principles and Applications Prentice Hall, 2001. M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids",

More information

Trpzip-based beta hairpin temperature jump IR studies enhanced by sitespecific

Trpzip-based beta hairpin temperature jump IR studies enhanced by sitespecific Trpzip-based beta hairpin temperature jump IR studies enhanced by sitespecific isotope labeling Carsten Kretjschi 1, Karin auser 1, Rong uang 2, Tim Keiderling 2 1 Institute of Biophysics, University of

More information

Ab-initio protein structure prediction

Ab-initio protein structure prediction Ab-initio protein structure prediction Jaroslaw Pillardy Computational Biology Service Unit Cornell Theory Center, Cornell University Ithaca, NY USA Methods for predicting protein structure 1. Homology

More information

Equilibrium sampling of self-associating polymer solutions: A parallel selective tempering approach

Equilibrium sampling of self-associating polymer solutions: A parallel selective tempering approach THE JOURNAL OF CHEMICAL PHYSICS 123, 124912 2005 Equilibrium sampling of self-associating polymer solutions: A parallel selective tempering approach Chakravarthy Ayyagari, a Dmitry Bedrov, and Grant D.

More information

High Temperature Unfolding Simulations of the TRPZ1 Peptide

High Temperature Unfolding Simulations of the TRPZ1 Peptide 4444 Biophysical Journal Volume 94 June 2008 4444 4453 High Temperature Unfolding Simulations of the TRPZ1 Peptide Giovanni Settanni and Alan R. Fersht Centre for Protein Engineering, Cambridge, United

More information

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

Folding of small proteins using a single continuous potential

Folding of small proteins using a single continuous potential JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 17 1 MAY 2004 Folding of small proteins using a single continuous potential Seung-Yeon Kim School of Computational Sciences, Korea Institute for Advanced

More information

ALGORITHMS FOR BUILDING MODELS OF MOLECULAR MOTION FROM SIMULATIONS

ALGORITHMS FOR BUILDING MODELS OF MOLECULAR MOTION FROM SIMULATIONS ALGORITHMS FOR BUILDING MODELS OF MOLECULAR MOTION FROM SIMULATIONS A DISSERTATION SUBMITTED TO THE DEPARTMENT OF COMPUTER SCIENCE AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY IN PARTIAL

More information

Cecilia Clementi s research group.

Cecilia Clementi s research group. Cecilia Clementi s research group http://leonardo.rice.edu/~cecilia/research/ Proteins don t have a folding problem it s we humans that do! Cartoons by Larry Gonick In principle, the laws of physics completely

More information

Paul Sigler et al, 1998.

Paul Sigler et al, 1998. Biological systems are necessarily metastable. They are created, modulated, and destroyed according to a temporal plan that meets the survival needs of the cell, organism, and species...clearly, no biological

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

Molecular Mechanics. Yohann Moreau. November 26, 2015

Molecular Mechanics. Yohann Moreau. November 26, 2015 Molecular Mechanics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Mechanics, Label RFCT 2015 November 26, 2015 1 / 29 Introduction A so-called Force-Field

More information

The Molecular Dynamics Method

The Molecular Dynamics Method The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation

More information

Docking. GBCB 5874: Problem Solving in GBCB

Docking. GBCB 5874: Problem Solving in GBCB Docking Benzamidine Docking to Trypsin Relationship to Drug Design Ligand-based design QSAR Pharmacophore modeling Can be done without 3-D structure of protein Receptor/Structure-based design Molecular

More information

Exercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions

Exercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,

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

Temperature and Pressure Controls

Temperature 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

Supporting Online Material for

Supporting Online Material for www.sciencemag.org/cgi/content/full/330/6002/341/dc1 Supporting Online Material for Atomic-Level Characterization of the Structural Dynamics of Proteins David E. Shaw,* Paul Maragakis, Kresten Lindorff-Larsen,

More information

Protein folding is a fundamental problem in modern structural

Protein folding is a fundamental problem in modern structural Protein folding pathways from replica exchange simulations and a kinetic network model Michael Andrec, Anthony K. Felts, Emilio Gallicchio, and Ronald M. Levy* Department of Chemistry and Chemical Biology

More information

Equilibrium Distribution from Distributed Computing (Simulations of Protein Folding)

Equilibrium Distribution from Distributed Computing (Simulations of Protein Folding) pubs.acs.org/jpcb Equilibrium Distribution from Distributed Computing (Simulations of Protein Folding) Riccardo Scalco and Amedeo Caflisch* Department of Biochemistry, University of Zurich, Winterthurerstrasse

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

The protein folding problem consists of two parts:

The protein folding problem consists of two parts: Energetics and kinetics of protein folding The protein folding problem consists of two parts: 1)Creating a stable, well-defined structure that is significantly more stable than all other possible structures.

More information

BIOINFORMATICS. Markov Dynamic Models for Long-Timescale Protein Motion

BIOINFORMATICS. Markov Dynamic Models for Long-Timescale Protein Motion BIOINFORMATICS Vol. 00 no. 00 2010 Pages 1 9 Markov Dynamic Models for Long-Timescale Protein Motion Tsung-Han Chiang 1, David Hsu 1 and Jean-Claude Latombe 2 1 Department of Computer Science, National

More information

Protein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall How do we go from an unfolded polypeptide chain to a

Protein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall How do we go from an unfolded polypeptide chain to a Lecture 11: Protein Folding & Stability Margaret A. Daugherty Fall 2004 How do we go from an unfolded polypeptide chain to a compact folded protein? (Folding of thioredoxin, F. Richards) Structure - Function

More information

SANDRO BOTTARO, PAVEL BANÁŠ, JIŘÍ ŠPONER, AND GIOVANNI BUSSI

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

Presenter: She Zhang

Presenter: She Zhang Presenter: She Zhang Introduction Dr. David Baker Introduction Why design proteins de novo? It is not clear how non-covalent interactions favor one specific native structure over many other non-native

More information

Simulating Folding of Helical Proteins with Coarse Grained Models

Simulating Folding of Helical Proteins with Coarse Grained Models 366 Progress of Theoretical Physics Supplement No. 138, 2000 Simulating Folding of Helical Proteins with Coarse Grained Models Shoji Takada Department of Chemistry, Kobe University, Kobe 657-8501, Japan

More information

Introduction to Computational Structural Biology

Introduction to Computational Structural Biology Introduction to Computational Structural Biology Part I 1. Introduction The disciplinary character of Computational Structural Biology The mathematical background required and the topics covered Bibliography

More information

Free energy calculations with alchemlyb

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

Molecular dynamics simulations of anti-aggregation effect of ibuprofen. Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov

Molecular dynamics simulations of anti-aggregation effect of ibuprofen. Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov Biophysical Journal, Volume 98 Supporting Material Molecular dynamics simulations of anti-aggregation effect of ibuprofen Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov Supplemental

More information

Protein Dynamics, Allostery and Function

Protein Dynamics, Allostery and Function Protein Dynamics, Allostery and Function Lecture 3. Protein Dynamics Xiaolin Cheng UT/ORNL Center for Molecular Biophysics SJTU Summer School 2017 1 Obtaining Dynamic Information Experimental Approaches

More information

Protein Structure Prediction

Protein Structure Prediction Protein Structure Prediction Michael Feig MMTSB/CTBP 2006 Summer Workshop From Sequence to Structure SEALGDTIVKNA Ab initio Structure Prediction Protocol Amino Acid Sequence Conformational Sampling to

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

Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability

Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability Dr. Andrew Lee UNC School of Pharmacy (Div. Chemical Biology and Medicinal Chemistry) UNC Med

More information

Structural Bioinformatics (C3210) Molecular Mechanics

Structural Bioinformatics (C3210) Molecular Mechanics Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational

More information

The strength of the hydrogen bond in the linking of protein

The strength of the hydrogen bond in the linking of protein Energetics of hydrogen bonds in peptides Sheh-Yi Sheu*, Dah-Yen Yang, H. L. Selzle, and E. W. Schlag *Department of Life Science, National Yang-Ming University, Taipei 112, Taiwan; Institute of Atomic

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

Calculation of the distribution of eigenvalues and eigenvectors in Markovian state models for molecular dynamics

Calculation of the distribution of eigenvalues and eigenvectors in Markovian state models for molecular dynamics THE JOURNAL OF CHEMICAL PHYSICS 126, 244101 2007 Calculation of the distribution of eigenvalues and eigenvectors in Markovian state models for molecular dynamics Nina Singhal Hinrichs Department of Computer

More information

What is Classical Molecular Dynamics?

What is Classical Molecular Dynamics? What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated

More information

Free energy calculations

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

Supporting Online Material for

Supporting Online Material for www.sciencemag.org/cgi/content/full/309/5742/1868/dc1 Supporting Online Material for Toward High-Resolution de Novo Structure Prediction for Small Proteins Philip Bradley, Kira M. S. Misura, David Baker*

More information

Hyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea

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

DISCRETE TUTORIAL. Agustí Emperador. Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING:

DISCRETE TUTORIAL. Agustí Emperador. Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING: DISCRETE TUTORIAL Agustí Emperador Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING: STRUCTURAL REFINEMENT OF DOCKING CONFORMATIONS Emperador

More information

2 Computational methods. 1 Introduction. 2.1 Simulation models

2 Computational methods. 1 Introduction. 2.1 Simulation models Molecular Dynamics Simulations of Folding and Insertion of the Ebola Virus Fusion Peptide into a Membrane Bilayer Mark A. Olson 1, In-Chul Yeh 2, and Michael S. Lee 1,3 1 US Army Medical Research Institute

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

Supporting Information

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