Limitations of temperature replica exchange (T-REMD) for protein folding simulations
|
|
- Marybeth Booker
- 6 years ago
- Views:
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
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 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 informationCan a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?
Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Ruhong Zhou 1 and Bruce J. Berne 2 1 IBM Thomas J. Watson Research Center; and 2 Department of Chemistry,
More informationPolypeptide 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 informationMany proteins spontaneously refold into native form in vitro with high fidelity and high speed.
Macromolecular Processes 20. Protein Folding Composed of 50 500 amino acids linked in 1D sequence by the polypeptide backbone The amino acid physical and chemical properties of the 20 amino acids dictate
More informationBiology Chemistry & Physics of Biomolecules. Examination #1. Proteins Module. September 29, Answer Key
Biology 5357 Chemistry & Physics of Biomolecules Examination #1 Proteins Module September 29, 2017 Answer Key Question 1 (A) (5 points) Structure (b) is more common, as it contains the shorter connection
More informationApplication of the Markov State Model to Molecular Dynamics of Biological Molecules. Speaker: Xun Sang-Ni Supervisor: Prof. Wu Dr.
Application of the Markov State Model to Molecular Dynamics of Biological Molecules Speaker: Xun Sang-Ni Supervisor: Prof. Wu Dr. Jiang Introduction Conformational changes of proteins : essential part
More informationSupporting 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 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 informationReplica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2)
pubs.acs.org/jpcb Replica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2) Lingle Wang, Richard A. Friesner, and B. J. Berne* Department of Chemistry,
More 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 informationEnergetics 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 informationOutline. The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation. Unfolded Folded. What is protein folding?
The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation By Jun Shimada and Eugine Shaknovich Bill Hawse Dr. Bahar Elisa Sandvik and Mehrdad Safavian Outline Background on protein
More informationRecovering 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 informationFree 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 informationHydrophobic 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 informationarxiv: 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 informationMolecular Interactions F14NMI. Lecture 4: worked answers to practice questions
Molecular Interactions F14NMI Lecture 4: worked answers to practice questions http://comp.chem.nottingham.ac.uk/teaching/f14nmi jonathan.hirst@nottingham.ac.uk (1) (a) Describe the Monte Carlo algorithm
More 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 informationChallenges 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 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 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 informationAlchemical free energy calculations in OpenMM
Alchemical free energy calculations in OpenMM Lee-Ping Wang Stanford Department of Chemistry OpenMM Workshop, Stanford University September 7, 2012 Special thanks to: John Chodera, Morgan Lawrenz Outline
More 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 informationShort 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 informationExploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials
Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials Dr. V. Subramanian Chemical Laboratory, IPC Division CSIR-Central Leather Research
More informationAnalysis 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 informationIntroduction 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 informationChoosing weights for simulated tempering
PHYSICAL REVIEW E 7, 173 7 Choosing weights for simulated tempering Sanghyun Park 1 and Vijay S. Pande 1, 1 Department of Chemistry, Stanford University, Stanford, California 935, USA Department of Structural
More 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 information2008 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 informationThe 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 informationProtein Folding experiments and theory
Protein Folding experiments and theory 1, 2,and 3 Protein Structure Fig. 3-16 from Lehninger Biochemistry, 4 th ed. The 3D structure is not encoded at the single aa level Hydrogen Bonding Shared H atom
More 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 informationPROTEIN-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 informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More informationWhat 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 informationProtein 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 informationTHE 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 informationHOW 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 informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationPHYSICAL 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 informationFibril 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 information3.320 Lecture 23 (5/3/05)
3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials
More informationComputer 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 informationHierarchical 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 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 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 informationTransition 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 informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationReplica 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 informationMethods of Computer Simulation. Molecular Dynamics and Monte Carlo
Molecular Dynamics Time is of the essence in biological processes therefore how do we understand time-dependent processes at the molecular level? How do we do this experimentally? How do we do this computationally?
More informationThermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning
Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning Ben Leimkuhler University of Edinburgh Joint work with C. Matthews (Chicago), G. Stoltz (ENPC-Paris), M. Tretyakov
More informationBiomolecules 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 informationTrpzip-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 informationAb-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 informationEquilibrium 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 informationHigh 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 informationw 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 informationFolding 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 informationALGORITHMS 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 informationCecilia 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 informationPaul 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 informationONETEP 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 informationMolecular 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 informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationDocking. 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 informationExercise 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 informationComputational Studies of the Photoreceptor Rhodopsin. Scott E. Feller Wabash College
Computational Studies of the Photoreceptor Rhodopsin Scott E. Feller Wabash College Rhodopsin Photocycle Dark-adapted Rhodopsin hn Isomerize retinal Photorhodopsin ~200 fs Bathorhodopsin Meta-II ms timescale
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationSupporting 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 informationProtein 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 informationEquilibrium 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 informationRoutine 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 informationThe 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 informationBIOINFORMATICS. 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 informationProtein 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 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 informationPresenter: 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 informationSimulating 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 informationIntroduction 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 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 informationMolecular 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 informationProtein 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 informationProtein 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 informationCS 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 informationProteins 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 informationStructural 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 informationThe 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 informationStatistical Mechanics for Proteins
The Partition Function From Q all relevant thermodynamic properties can be obtained by differentiation of the free energy F: = kt q p E q pd d h T V Q ), ( exp 1! 1 ),, ( 3 3 3 ),, ( ln ),, ( T V Q kt
More informationCalculation 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 informationWhat 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 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 informationSupporting 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 informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationImproving 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 informationDISCRETE 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 information2 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 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 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 information