Energy Landscapes and Accelerated Molecular- Dynamical Techniques for the Study of Protein Folding
|
|
- Jade Short
- 6 years ago
- Views:
Transcription
1 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 November 3,
2 What is a Protein? See 2
3 Peptide Units Combine to Form Proteins Space Filling Model Stick Model Ribbon Model See 3
4 The Native Structure of Proteins is Unique and Complex Bovine F1 ATP Synthase See 4
5 So How Do Proteins Fold? Denatured state Folded (native) state? Thanks to John Chodera, Dill Biophysics Group, UCSF 5
6 The Protein Folding Problem The holy grail : Compute structure from sequence ab initio This involves exploring the conformational space of the protein using first-principle physics methods in order to identify the native structure. 6
7 The Energy Funnel José Onuchic, Workshop IIII Structural Proteomics, IPAM-UCLA, May
8 A Realistic Energy Funnel Landscape Denatured state Folding Funnel Native state 8
9 How to Theoretically Explore Folding? The Molecular Dynamics Approach The AMBER forcefield: Thanks to John Chodera, Dill Biophysics Group, UCSF 9
10 Molecular Dynamics Numerically solve Newton s 2 nd law over small time increments to compute the position and momentum of the atoms in the protein over time. m i N dvi = dt = j 1 j i F i, j q = p = p H; H q 10
11 All-Atom MD Calculations of Folding These calculations include all atoms in the protein, implicit or explicit solvent, and account for all forces between the various atoms. Each MD timestep is no more than 2 fs (10-15 s). A 4 µs calculation would require 2 x 10 9 timesteps to simulate observed real life folding times. At 10 ms per timestep, such a simulation would require 231 single CPU days to compute a single trajectory on a modern high-end workstation. Brute force, single trajectory full-folding folding calculations are thus in general not practical at this time. 11
12 Statistical Molecular Dynamics Consider a protein consisting of N atoms whose instantaneous positions and momenta at time t are given by the 3N dimensional vectors r t and p ( t). The time average of a property A is defined as: 1 A dt A r t p t t τ τ = lim (), () Numerically, we calculate this as: τ 0 ( ) M 1 A A r t p t t ( ), ( ) m m M = m 1 12
13 Statistical Molecular Dynamics (2) Note that the ergodic hypothesis says that A t = where the ensemble average of property A is given by A A dp dr A r, p ρ r, p e = e ( ) ( ) and the probability density in the canonical (PVT) ensemble is given by E ( r, p) ρ ( r, p) = exp Q kt B 13
14 Folding Times are Often Dominated by the Time to Escape Kinetic Traps Low energy trajectories tend to get trapped in metastable local minima. Higher energy (temperature) trajectories avoid kinetic traps. By varying temperatures to kick a trajectory out of a kinetic trap, we can accelerate MD folding simulations. 14
15 The Replica Exchange Method Simulations are conducted in parallel at temperatures T, T,, T, T 1 2 Λ 1 Λ T 1 T 2 P = { ( ) } min 1,exp T Λ 1 T Λ Because each replica can undergo cycles of heating and cooling, conformation space is broadly sampled. Thanks to John Chodera, Dill Biophysics Group, UCSF 15
16 Details of the REM (1) ( ) ( () ) [ i Let (1) ] [ i ( Λ ) ] [] 1 [ Λ X = x ] 1,, xλ = x,, x λ 1 λ( Λ) stand for a state in this generalized ensemble. The superscript and subscript label the replica and temperature respectively. Each state X is specified by Λ sets of coordinates and momenta of N atoms in replica i(λ) at temperature T λ. We assume each replica is in the canonical (PVT) ensemble. i( λ) i( λ) i( λ) ( ), λ x r p λ 16
17 Details of the REM (2) The replicas are non-interacting. Therefore, the weight factor for the state X in this generalized ensemble is given by the product of the Boltzmann factors for each replica (each at a single temperature): where Λ WREM ( X) = exp β () H r, p λ i i= 1 ( [] i [] i ) Λ = exp βλ H r, λ= 1 1 β = kt B ( i( λ) i( λ) p ) 17
18 Details of the REM (3) Now consider exchanging a pair of replicas in the generalized ensemble. Specifically, exchange the temperatures between replicas i and j. Thus: ( ) ( [] i [ j] ) [ j] [ i] λ γ λ γ X =, x,, x, X =, x,, x, where ( ) [] [] [] x r p x r p [ ] [ ] [ ] x r p x r p ( ) [] [] [], λ, i i i i i i λ λ ( ) ( ) [ ] [ ] [ ], γ, j j j j j j γ γ γ λ 18
19 Details of the REM (4) The post-exchange momenta are defined by: p p [] i γ [] i [ j] [ j] This specific scaling of momenta correctly rescales the kinetic energy to be consistent with the equipartition of energy theorem K p ( ) where K is the kinetic energy. T T T λ T T λ γ p p N 2 pk 3 = = Nk T B 2m 2 k = 1 k T 19
20 Details of the REM (5) In order for this exchange process to converge toward an equilibrium distribution in the canonical ensemble (for each temperature), it is sufficient to impose the detailed balance condition on the transition probability w( X X ): W X w X X = W X w X X REM ( ) ( ) ( ) ( ) REM This can be satisfied by the usual Metropolis criterion: [] [ ] ( ) ( ) i j 1 for 0 w X X w xλ xγ = exp( ) otherwise where U is the potential energy and = βγ β λ ( ( []) ( [ ])) i j U r U r 20
21 Details of the REM (6) The canonical expectation value of a quantity A at temperature T λ λ = 1,, Λ can then be calculated by ( ) 1 M f m m= 1 ( ( λ, ) ( )) λ m A = A x t Tλ M Expectations values for intermediate temperatures can be obtained using the multiple-histogram reweighting technique (WHAM). 21
22 The Replica Exchange Method Simulations are conducted in parallel at temperatures T, T,, T, T 1 2 Λ 1 Λ T 1 T 2 P = { ( ) } min 1,exp T Λ 1 T Λ Because each replica can undergo cycles of heating and cooling, conformation space is broadly sampled. Thanks to John Chodera, Dill Biophysics Group, UCSF 22
23 Advantages of REM Exchanging temperatures kicks replicas out of kinetic traps and accelerates the simulated folding process. Each replica calculation can be done on a different processor on a parallel computer. There is little interprocessor communication required, so the problem is virtually embarrassingly parallel. Thus, on a parallel computer with sufficient processors, there is no wall clock penalty for the additional work required in REM while the simulation time to reach the native state is reduced because the sampling of phase space is improved. 23
24 Results for Met-enkephallin Time series of temperature exchange and total potential energy for one replica in an 8 replica REM calculation of Met-enkephalin (Tyr-Gly-Gly-Phe-Met). Yuji Sugita and Yuko Okamoto, Chemical Physics Letters 314 (1999)
25 Results for Met-enkephallin Yuji Sugita and Yuko Okamoto, Chemical Physics Letters 314 (1999)
26 MREM Method Multiplexed REM (MREM) is closely related to REM, but uses multiple replicas at the same temperature in addition to replicas at different temperatures. This improves sampling and further reduces folding times. Rhee and Pande performed MREM calculations for the 23 amino acid model protein BBA5. The calculations used 20 temperatures and 200 replicas per temperature. The MREM calculations reached the native state approximately 10 times faster than a constant temperature molecular dynamics calculation. Rhee and Pande, Biophysical Journal 84 (2003)
27 Results for BBA5 Using MREM Rhee and Pande, Biophysical Journal 84 (2003)
28 Results for BBA5 Using MREM Rhee and Pande, Biophysical Journal 84 (2003)
29 Results for BBA5 Using MREM Rhee and Pande, Biophysical Journal 84 (2003)
30 Another Way to Accelerate Folding t 1/ k Note that transition rates merely reflect the average time required for a transition. Some trajectories will transition faster than others and some slower. The parallel replica method exploits this as another method for accelerating simulations. 30
31 The Parallel Replica Method Pande, et. al., Biopolymers, vol 68, Consider N identical conformation replicas at the same temperature, but with randomized velocities. The probability that a particular simulation has crossed the barrier is: P t k kt () 1 = exp( ) For M simulations, the probability of the first stimulation crossing the barrier is t PM t = MP t 1 dτp = Mk exp Mkt () () ( ) ( ) τ M 1 31
32 Parallel replica requires infrequent exchange of information between replicas. It is therefore nearly embarrassingly parallel. See distributes the computation of each replica to individual PCs across the World Wide Web. To date, over 1 million computers have participated in folding@home. This has made possible the calculation of realistic complete folds for moderate sized proteins out to the equivalent of milliseconds. 32
33 Conclusions Folding rates are dominated by time spent in kinetic traps. In silico brute-force folding of proteins of biological interest via a single molecular dynamics trajectory is not in general feasible with current computers. The replica exchange and parallel replica techniques avoid this problem by artificially minimizing the time spent in kinetic traps (at the expense of destroying information about the kinetics). Speedups of 10x in calculation of native structures have been realized for MREM compared to constant temperature MD. Speedups of several orders of magnitude for calculation of native structures have been realized using parallel replica. 33
34 Acknowledgements I would like to acknowledge the invaluable assistance in preparing this talk of John Chodera and Ken Dill of the Dill Biophysics Group, UCSF. 34
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 informationLimitations of temperature replica exchange (T-REMD) for protein folding simulations
Limitations of temperature replica exchange (T-REMD) for protein folding simulations Jed W. Pitera, William C. Swope IBM Research pitera@us.ibm.com Anomalies in protein folding kinetic thermodynamic 322K
More informationFree-Energy Calculations in Protein Folding by Generalized-Ensemble Algorithms
Free-Energy Calculations in Protein Folding by Generalized-Ensemble Algorithms Yuji Sugita 1 and Yuko Okamoto 2 Department of Theoretical Studies Institute for Molecular Science Okazaki, Aichi 444-8585,
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 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 informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationSmart walking: A new method for Boltzmann sampling of protein conformations
Smart walking: A new method for Boltzmann sampling of protein conformations Ruhong Zhou and B. J. Berne Department of Chemistry and the Center for Biomolecular Simulation, Columbia University, New York,
More informationJavier Junquera. Statistical mechanics
Javier Junquera Statistical mechanics From the microscopic to the macroscopic level: the realm of statistical mechanics Computer simulations Thermodynamic state Generates information at the microscopic
More informationGeneralized-Ensemble Algorithms for Molecular Simulations of Biopolymers
Generalized-Ensemble Algorithms for Molecular Simulations of Biopolymers Ayori Mitsutake, a,b Yuji Sugita, a,b and Yuko Okamoto a,b,1 a Department of Theoretical Studies Institute for Molecular Science
More informationConvergence 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 informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More 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 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 informationParallel Tempering Algorithm for. Conformational Studies of Biological Molecules
Parallel Tempering Algorithm for Conformational Studies of Biological Molecules Ulrich H.E. Hansmann 1 Department of Theoretical Studies, Institute for Molecular Science Okazaki, Aichi 444, Japan ABSTRACT
More informationAnalysis 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 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 informationStatistical thermodynamics (mechanics)
Statistical thermodynamics mechanics) 1/15 Macroskopic quantities are a consequence of averaged behavior of many particles [tchem/simplyn.sh] 2/15 Pressure of ideal gas from kinetic theory I Molecule =
More informationMolecular Dynamics and Accelerated Molecular Dynamics
Molecular Dynamics and Accelerated Molecular Dynamics Arthur F. Voter Theoretical Division National Laboratory Lecture 3 Tutorial Lecture Series Institute for Pure and Applied Mathematics (IPAM) UCLA September
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 informationWhy Proteins Fold. How Proteins Fold? e - ΔG/kT. Protein Folding, Nonbonding Forces, and Free Energy
Why Proteins Fold Proteins are the action superheroes of the body. As enzymes, they make reactions go a million times faster. As versatile transport vehicles, they carry oxygen and antibodies to fight
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 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 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 informationEffective stochastic dynamics on a protein folding energy landscape
THE JOURNAL OF CHEMICAL PHYSICS 125, 054910 2006 Effective stochastic dynamics on a protein folding energy landscape Sichun Yang, a José N. Onuchic, b and Herbert Levine c Center for Theoretical Biological
More informationMolecular Structure Prediction by Global Optimization
Molecular Structure Prediction by Global Optimization K.A. DILL Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, CA 94118 A.T. PHILLIPS Computer Science
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 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 informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
More informationDominant Paths in Protein Folding
Dominant Paths in Protein Folding Henri Orland SPhT, CEA-Saclay France work in collaboration with P. Faccioli, F. Pederiva, M. Sega University of Trento Henri Orland Annecy meeting 2006 Outline Basic notions
More informationCalculation of point-to-point short time and rare. trajectories with boundary value formulation
Calculation of point-to-point short time and rare trajectories with boundary value formulation Dov Bai and Ron Elber Department of Computer Science, Upson Hall 413 Cornell University, Ithaca, NY 14853
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 informationEquilibrium Molecular Thermodynamics from Kirkwood Sampling
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source
More 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 informationEquilibrium Molecular Thermodynamics from Kirkwood Sampling
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source
More informationEvolutionary design of energy functions for protein structure prediction
Evolutionary design of energy functions for protein structure prediction Natalio Krasnogor nxk@ cs. nott. ac. uk Paweł Widera, Jonathan Garibaldi 7th Annual HUMIES Awards 2010-07-09 Protein structure prediction
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 informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationThe 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 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 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 informationA new combination of replica exchange Monte Carlo and histogram analysis for protein folding and thermodynamics
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 3 15 JULY 2001 A new combination of replica exchange Monte Carlo and histogram analysis for protein folding and thermodynamics Dominik Gront Department of
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
More informationEnergy Minimization of Protein Tertiary Structure by Parallel Simulated Annealing using Genetic Crossover
Minimization of Protein Tertiary Structure by Parallel Simulated Annealing using Genetic Crossover Tomoyuki Hiroyasu, Mitsunori Miki, Shinya Ogura, Keiko Aoi, Takeshi Yoshida, Yuko Okamoto Jack Dongarra
More informationEfficient 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 informationProtein Structure Analysis
BINF 731 Protein Modeling Methods Protein Structure Analysis Iosif Vaisman Ab initio methods: solution of a protein folding problem search in conformational space Energy-based methods: energy minimization
More informationAccelerated Quantum Molecular Dynamics
Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended
More informationSimulations with MM Force Fields. Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi
Simulations with MM Force Fields Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi Some slides taken with permission from Howard R. Mayne Department of Chemistry University of New Hampshire Walking
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 6 Aug 2003
Metropolis Algorithms in Generalized Ensemble Yuko Okamoto arxiv:cond-mat/0308119v1 [cond-mat.stat-mech] 6 Aug 2003 Department of Theoretical Studies Institute for Molecular Science Okazaki, Aichi 444-8585,
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
More informationSupplementary Information Intrinsic Localized Modes in Proteins
Supplementary Information Intrinsic Localized Modes in Proteins Adrien Nicolaï 1,, Patrice Delarue and Patrick Senet, 1 Department of Physics, Applied Physics and Astronomy, Rensselaer Polytechnic Institute,
More information9.3: Separable Equations
9.3: Separable Equations An equation is separable if one can move terms so that each side of the equation only contains 1 variable. Consider the 1st order equation = F (x, y). dx When F (x, y) = f (x)g(y),
More informationMilestoning: Extending time scales of molecular simulations
Milestoning: Extending time scales of molecular simulations Anthony West, Peter Majek, and RE Department of Chemistry and Biochemistry and Institute of Computational Sciences and Engineering (ICES) University
More informationOptimization Methods via Simulation
Optimization Methods via Simulation Optimization problems are very important in science, engineering, industry,. Examples: Traveling salesman problem Circuit-board design Car-Parrinello ab initio MD Protein
More informationMulticanonical parallel tempering
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 13 1 APRIL 2002 Multicanonical parallel tempering Roland Faller, Qiliang Yan, and Juan J. de Pablo Department of Chemical Engineering, University of Wisconsin,
More informationLecture 8: Computer Simulations of Generalized Ensembles
Lecture 8: Computer Simulations of Generalized Ensembles Bernd A. Berg Florida State University November 6, 2008 Bernd A. Berg (FSU) Generalized Ensembles November 6, 2008 1 / 33 Overview 1. Reweighting
More informationProtein Folding. I. Characteristics of proteins. C α
I. Characteristics of proteins Protein Folding 1. Proteins are one of the most important molecules of life. They perform numerous functions, from storing oxygen in tissues or transporting it in a blood
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 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 informationSerial Replica Exchange
1416 J. Phys. Chem. B 2007, 111, 1416-1423 Serial Replica Exchange Morten Hagen, Byungchan Kim, Pu Liu, Richard A. Friesner, and B. J. Berne* Department of Chemistry and Center for Biomolecular Simulation,
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 informationBiological Thermodynamics
Biological Thermodynamics Classical thermodynamics is the only physical theory of universal content concerning which I am convinced that, within the framework of applicability of its basic contents, will
More informationApplication of Computational Modelling to Protein Folding and Aggregation Studies
Application of Computational Modelling to Protein Folding and Aggregation Studies NEVENA TODOROVA, ANDREW HUNG AND IRENE YAROVSKY School of Applied Sciences RMIT University GPO Box 2476V Melbourne Victoria
More informationParallel Tempering Algorithm in Monte Carlo Simulation
Parallel Tempering Algorithm in Monte Carlo Simulation Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL) Metropolis Algorithm on Ising Model
More informationMonte Carlo simulation of proteins through a random walk in energy space
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 16 22 APRIL 2002 Monte Carlo simulation of proteins through a random walk in energy space Nitin Rathore and Juan J. de Pablo a) Department of Chemical Engineering,
More information4/18/2011. Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania
1. Introduction Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania Bibliography Computer experiments Ensemble averages and time averages Molecular
More informationUsing 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 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 informationStatistical Mechanics in a Nutshell
Chapter 2 Statistical Mechanics in a Nutshell Adapted from: Understanding Molecular Simulation Daan Frenkel and Berend Smit Academic Press (2001) pp. 9-22 11 2.1 Introduction In this course, we will treat
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 Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
More informationMolecular Dynamics in Principal Component Space
pubs.acs.org/jpcb Molecular Dynamics in Principal Component Space Servaas Michielssens,*, Titus S. van Erp, Carsten Kutzner, Arnout Ceulemans, and Bert L. de Groot Department of Chemistry and Centre for
More informationAb initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
More 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 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 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 informationClassical Monte-Carlo simulations
Classical Monte-Carlo simulations Graduate Summer Institute Complex Plasmas at the Stevens Insitute of Technology Henning Baumgartner, A. Filinov, H. Kählert, P. Ludwig and M. Bonitz Christian-Albrechts-University
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 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 informationLecture 11: Protein Folding & Stability
Structure - Function Protein Folding: What we know Lecture 11: Protein Folding & Stability 1). Amino acid sequence dictates structure. 2). The native structure represents the lowest energy state for a
More informationProtein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall Protein Folding: What we know. Protein Folding
Lecture 11: Protein Folding & Stability Margaret A. Daugherty Fall 2003 Structure - Function Protein Folding: What we know 1). Amino acid sequence dictates structure. 2). The native structure represents
More informationOptimized statistical ensembles for slowly equilibrating classical and quantum systems
Optimized statistical ensembles for slowly equilibrating classical and quantum systems IPAM, January 2009 Simon Trebst Microsoft Station Q University of California, Santa Barbara Collaborators: David Huse,
More informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More informationIntroduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou
Introduction to Simulation - Lectures 17, 18 Molecular Dynamics Nicolas Hadjiconstantinou Molecular Dynamics Molecular dynamics is a technique for computing the equilibrium and non-equilibrium properties
More information8.333: Statistical Mechanics I Problem Set # 3 Due: 10/18/13. Kinetic Theory
8.333: Statistical Mechanics I Problem Set # 3 Due: 10/18/13 1. Poisson Brackets: Kinetic Theory (a) Show that for observable O(p(µ),q(µ)), do/dt = {O,H}, along the time trajectory of any micro state µ,
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 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 informationCHE210D Principles of Modern Molecular Simulation Methods
CHE210D Principles of Modern Molecular Simulation Methods Instructor: M. Scott Shell Spring Quarter 2012 TuTh 12:30-1:45pm in Eng II 1519 www.engr.ucsb.edu/~shell/che210d The goals of this course This
More informationOrganization of NAMD Tutorial Files
Organization of NAMD Tutorial Files .1.1. RMSD for individual residues Objective: Find the average RMSD over time of each residue in the protein using VMD. Display the protein with the residues colored
More informationSampling the conformation space of complex systems, such as
Replica exchange with solute tempering: A method for sampling biological systems in explicit water Pu Liu*, Byungchan Kim*, Richard A. Friesner, and B. J. Berne Department of Chemistry and Center for Biomolecular
More informationThermodynamics. Entropy and its Applications. Lecture 11. NC State University
Thermodynamics Entropy and its Applications Lecture 11 NC State University System and surroundings Up to this point we have considered the system, but we have not concerned ourselves with the relationship
More informationProtein Folding In Vitro*
Protein Folding In Vitro* Biochemistry 412 February 29, 2008 [*Note: includes computational (in silico) studies] Fersht & Daggett (2002) Cell 108, 573. Some folding-related facts about proteins: Many small,
More informationRare Events. Transition state theory Bennett-Chandler Approach 16.2 Transition path sampling16.4
Rare Events Transition state theory 16.1-16.2 Bennett-Chandler Approach 16.2 Transition path sampling16.4 Outline Part 1 Rare event and reaction kinetics Linear Response theory Transition state theory
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
More informationLecture 21 (11/3/17) Protein Stability, Folding, and Dynamics Hydrophobic effect drives protein folding
Reading: Ch4; 142-151 Problems: Ch4 (text); 14, 16 Ch6 (text); 1, 4 NEXT (after exam) Reading: Ch8; 310-312, 279-285, 285-289 Ch24; 957-961 Problems: Ch8 (text); 1,2,22 Ch8 (study-guide:facts); 1,2,3,4,5,9,10
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 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 informationIsothermal experiments characterize time-dependent aggregation and unfolding
1 Energy Isothermal experiments characterize time-dependent aggregation and unfolding Technical ote Introduction Kinetic measurements have, for decades, given protein scientists insight into the mechanisms
More informationComputer simulations of protein folding with a small number of distance restraints
Vol. 49 No. 3/2002 683 692 QUARTERLY Computer simulations of protein folding with a small number of distance restraints Andrzej Sikorski 1, Andrzej Kolinski 1,2 and Jeffrey Skolnick 2 1 Department of Chemistry,
More informationExchange frequency in replica exchange molecular dynamics
THE JOURNAL OF CHEMICAL PHYSICS 128, 024103 2008 Exchange frequency in replica exchange molecular dynamics Daniel Sindhikara Department of Physics, University of Florida, Gainesville, Florida 32611-8435,
More informationAppendix C extra - Concepts of statistical physics Cextra1-1
Appendix C extra - Concepts of statistical physics Cextra1-1 Appendix C extra - Concepts of statistical physics A conventional introductory course in physics concentrates on one-body motion: the change
More information