Some Accelerated/Enhanced Sampling Techniques
|
|
- Sybil Palmer
- 5 years ago
- Views:
Transcription
1 Some Accelerated/Enhanced Sampling Techniques The Sampling Problem Accelerated sampling can be achieved by addressing one or more of the MD ingredients Back to free energy differences... Examples of accelerated sampling techniques
2 ( G RT ln e Ek E 0 0 S energy kbt k S Sampling: a Problem?! (Free) energy landscape is complicated and huge ) Obtaining a complete picture of the possible conformations and their probability of occurring is hampered by the sheer number of them and the barriers between local minima in the (free) energy landscape co te na di or 2 e1 t a n i coord Pictures adapted from presentation by Sander Pronk and Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006)
3 Enhanced Sampling Opportunities Multiscaling: Coarse-graining: reducing the number of degrees of freedom, preserving the relevant physics reducing detail in the surroundings leading to effective interactions Biasing: Jumping: adapting interactions to exchanging snapshots reduce phase space and/or between conditions to smoothen the free energy overcome barriers landscape Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006)
4 Enhanced Sampling Opportunities Multiscaling: Coarse-graining: THERE IS NO SUCH THING AS A FREE LUNCH! reducing the number of degrees of freedom, preserving the relevant physics reducing detail in the surroundings leading to effective interactions Biasing: Jumping: adapting interactions to exchanging snapshots reduce phase space and/or between conditions to smoothen the free energy overcome barriers landscape Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006)
5 Q: reaction quotient K: equilibrium constant Free Energy Differences from Simulations Direct by counting 1 2 eq p Q12eq = K12 = 2eq = p1 ( ) kbt ( e ) kbt e Ei E 0 i 2 E j E 0 j 1 ( =e G20 G10 ) RT Boltzmann weights! a2 c2 N 2V1 p2 Reliable value for free energy Q12 = =! is obtained only if the a1 c1 N1 V2 p1 difference statistics are good enough: we A thermodynamic state is in fact a collection of configurations! need many transitions between the states and full sampling of each state to capture the entropy
6 ΔG120 = G20 G10 = RT ln K12 Free Energy from Simulations: ONE of the TRICKS Find extra (also called biasing) potential here ΔU to make probabilities equal better statistics 1 p2eq K12 = eq = p1 2 β Ei e i 2 e βej =e 0 βδg12 j 1 p2 =1= p1 e i 2 e j 1 ( β Ei ΔU βej ) = eβδu e β Ei i 2 βej e j 1 A thermodynamic state is in fact a collection of configurations! =e 0 βδu βδg12 e ΔG120 = ΔU 1 E = 0; β = k BT 0
7 ΔG120 = G20 G10 = RT ln K12 Free Energy from Simulations: ONE of the TRICKS More general: make correction for biasing potential ( e β Ei ΔU 1 2 p2 = p1 i 2 e βej ) =e 0 βδu βδg12 e j 1 eq 2 eq 1 p K12 = = p e β Ei e βej i 2 =e 0 βδg12 j 1 p2 ΔG = ΔU RT ln p A thermodynamic state is in fact a collection of configurations! Input Measure with good statistics
8 ΔG 0 12 = RT ln p eq 2 eq p 1 Free energy differences from Simulations Weighted Histogram Analysis Method Apply a restraining potential at different points PMF ( ξ ) = RT ln p ( ξ ) d=1.0 nm d=1.2 nm gmx wham ( ) p ξ 0 Methanol dimer PMF Potential of Mean Force Add harmonic restraining potential to the distance between centers-of-mass ΔU R ( ξ,d) = K ( 2 ξ d ) 2 The PMF is obtained by reweighting the conformations
9 Biasing Metadynamics Adds penalizing potential for conformations already visited Potential energy landscape is changed on the fly Many kinds of variations are possible depending on how the penalizing potentials are implemented or to which degrees of freedom they are applied True Boltzmann probabilities can be determined afterwards by reweighting You pay by losing the true dynamics of the system and possibly sampling utterly uninteresting and/or irrelevant parts of phase space (you need to have an idea of the relevant degree(s) of freedom) From Review by Bernardi et al. BBA 1850, 872 (2015)
10 Local Elevation Example Biasing Gaussian-shaped penalty function as function of torsion angle Glucose ring puckering MD LE 300 K 2 kj mol -1 penalty increments Glucose ring puckering states characterized by dihedral angles. A normal MD simulation samples only one state; local elevation drives the molecule to the other state and later on leads to multiple transitions Christen et al. J. Comput. Chem. 26, 1719 (2005)
11 The Idea of Replica Exchange Run a number of simulations at the same time (replica s of the system) under different conditions Different temperatures is the most popular Different Hamiltonians is a viable option Every so often, attempts are made to swap conformations between (neighboring) conditions The criterion is based on relative Boltzmann weights, this ensures correct overall Boltzmann statistics at one particular condition Jumping p exchange = min 1, e e E 1 k B T 2 + E 2 k B T 1 E 1 k B T 1 + E 2 k B T 2 From Review by Bernardi et al. BBA 1850, 872 (2015)
12 Jumping Replica Exchange Attempts to make trajectories jump between different conditions Temperatures: at higher temperatures, minima are relatively less likely and the free energy landscape is smoother, thus barriers are more easily overcome Hamiltonians: by switching off/on interactions, the landscape may be smoothened or biased towards interesting regions of phase space True Boltzmann probabilities determined by collecting the parts of trajectories at one temperature or hamiltonian You pay by losing the dynamical continuity at one temperature/ hamiltonian and by simulating irrelevant (high temperature) or non-existing (non-physical hamiltonian) systems
13 gmx rms RMSD: structural similarity Root-Mean-Square Deviation average over all particles at one point in time Extensively used in Protein Modeling 1 r k t N p k ( ( ) ref r ) 2 k Here, Np is the number of particles (atoms/beads) in the molecule ref r is the position of particle k k in the reference structure r k ( t) is the position of particle k at time t cf Basic Exercise
14 Replica Exchange Example The Trp-cage peptide is a designed fast folder 64 replica s (13 ns) K T-REMD does not require previous knowledge about interesting degrees of freedom; in fact, it can be used to discover them! Clusters found in REMD simulation (left) characterized by overall rmsd and helix rmsd (rmsdhx) compared to native structures (a: PDB, b: OPLS, c: GROMOS) btw: this is an example of projecting the phase space on a few degrees of freedom (cf PCA) Juraszek and Bolhuis P. N. A. S. 103, (2006)
15 Replica Exchange Example The Trp-cage peptide is a designed fast folder 64 replica s (13 ns) K From folded The free energy landscape is NOT CONVERGED because REMD simulations starting from different conformations do not yield the same probabilities Free energy landscapes characterized by radius of gyration (rg), ratio of native contacts (ρ), solvent accessible surface (sas), helix rmsd (rmsd_hx) as contours From unfolded Juraszek and Bolhuis P. N. A. S. 103, (2006)
16 Coarse-graining Coarse-graining Retaining the relevant degrees of freedom Interesting physics often does not need all details: similar phenomena occur on different length and time scales justifying reduced units and reduced models Coarse-grained models come in many more flavors than atomistic models because they can be specific or generic and tailored to a particular field of research Coarse-graining within a model identifies the slow but overall large motions: these take a system between different states You pay by losing detail and chemical specificity: the model represents one or more (classes of) molecules and you may sometimes miss particularities of the intended (real) system
17 Coarse-graining Example: Martini Model Marrink et al. J. Phys. Chem. B 108, 750 (2004); J. Phys. Chem. B 111, 7812 (2007) Check out
18 Computational gain CG models Effective interactions between beads Reduction in no. degrees of freedom Fewer interactions to be calculated Lower density of particles allows longer lengthscales Fastest motions on ~10 fs time-scale iso ~1 fs Allows larger time-step Speed-up of up to a factor ~1,000 for 4-to-1 heavy-atom mapped particle-based models
19 Comparison Martini to GROMOS Lipid bilayer in liquid-crystalline state Electron density profile across the bilayer shows how the components of the system are distributed (dashed is CG) Snapshots, side view and looking on top of the bilayer Marrink et al. J. Phys. Chem. B 108, 750 (2004)
20 gmx anaeig Characterizing a collection of structures Are conformations realistic? Time scale of sampling Schlitter s formula for configurational entropy Upper bound Approximation for harmonic oscillator S true S = k B 2 ln 1 + k BTe 2! 2 D Mass-weighted covariance matrix Procedure Fit (part of) the structure to remove translation (and rotation) J. Schlitter Chem. Phys. Lett. 215, 617 (1993)
21 Configurational Entropy Build-up S true S = k B 2 ln 1+ k B Te2! 2 D For atomistic hexadecane and lipid tail in a bilayer Baron et al J. Phys. Chem. B 110, (2006)
22 Conf. Entropy Shows Enhanced Sampling For phospholipid DPPC in bilayer Average structures S true S = k B 2 ln 1+ k B Te2 D! 2 Build-up of sampling after mapping sn1 sn2 Baron et al. J. Phys. Chem. B 110, (2006)
23 Multiscaling Multiscaling Combining the best of different level models Concurrent multiscaling: combining two or more levels in a single simulation, e.g. QM-MM, AA-CG Sequential multiscaling: generating snapshots at less detailed or cheaper level and using these at more detailed level to obtain proper statistics or distribution Active developments ongoing in combining different level molecular models You pay by doing a lot of work to make sure that the less detailed model is (thermodynamically) compatible with the more detailed model or accept artifacts around the boundary between the different models
24 Antimicrobial Peptides Induce Pores Atomistic simulations show the pore formation Pore seen in simulation looks unexpected barrel stave torroidal Proposed pore models showed regular arrangements with membrane spanning peptides Leontiadou et al. J. Am. Chem. Soc. 128, (2006)
25 Sequential Multiscaling & Backmapping Use CG model to explore phase space Put in atomistic details in interesting states A membrane is attacked by peptides (UA ~100 ns) and forms a toroidal pore The peptides become membrane-spanning in a 24 µs CG simulation After back-mapping, this pore is seen to be stable also in the atomistic force field 24 µs 50 ns Rzepiela et al. Faraday Discussions 144, 431 (2010)
26 backward.py User Friendly Backmapping Routine Define rules for initial placing of finer detail based on positions of coarser model Cook 3-bead lipids Martini 12-bead lipids GROMOS 50 UAs Being able to easily connect three quite different resolution lipid models, structures that take very long to build or self-assemble at atomistic level can be generated in no time at all... Wassenaar et al. J. Chem. Theory Comput. 10, 676 (2014)
27 Concurrent Multiscaling Example OPLS-AA/L Valine tripeptide in CG polarizable Martini water (Advanced Tutorial) Interaction of peptide with water through virtual sites at center of mass of groups of atoms according to Martini model Forces on virtual sites are redistributed to the atoms Polarizable Martini water (blue 3-particle model) OPLS-AA peptide: licorice Virtual sites: green spheres
28 Adaptive Resolution Simulation (AdResS) Resolution changes depending on the position Simulation done with Espresso++ Available in GROMACS; recent improvements not yet in standard release Protein in water surrounded by Martini water - molecules diffusing into or out of a pre-defined sphere around the protein change between atomistic and Martini models Zavadlav et al. J. Chem. Phys. 140, (2014)
29 Multiscaling is an Exciting Field Cool technology, but many issues still to be resolved or investigated to fully deliver the promise Artifacts around transition regions because of different chemical potentials of the models Putting in atomistic details smoothly and efficiently Combine with multiple time-stepping Combine with QM region
30 Summary Coarse-graining: The Best Way to Enhance Sampling Depends on Your System and Research Question Multiscaling: reducing the number of degrees of freedom, preserving the relevant physics reducing detail in the surroundings leading to effective interactions Biasing: Jumping: adapting interactions to exchanging snapshots reduce phase space and/or between conditions to smoothen the free energy overcome barriers landscape Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006)
31 Thank you for your attention Recent Reviews: Bernardi et al., Enhanced sampling techniques in molecular dynamics simulations of biological systems, BBA 1850, 872 (2015) Cavalli et al., Investigating Drug Target Association and Dissociation Mechanisms Using Metadynamics-Based Algorithms, Acc. Chem. Res. 48, 277 (2015) Ingolfsson et al., The power of coarse graining in biomolecular simulations, WIRE Comput. Mol. Sci., 4, 225 (2014)
32 Hamiltonian Exchange Example A torsional barrier for rotation was varied between 40% (λ=0.6) and 100% (λ=0) of the nominal value 11 replica s 500 molecules all trans start 273 K At 273 K trans-gauche isomerization across the nominal barrier does not occur (the barrier is ~18 kj mol -1 ) Paths of the 11 replica simulations in λ-space are indicated by different colors Christen et al. J. Comput. Chem. 26, 1719 (2005)
33 Hamiltonian Exchange Example A torsional barrier for rotation was varied between 40% (λ=0.6) and 100% (λ=0) of the nominal value 11 replica s 500 molecules all trans start 273 K λ=0 replica RMSD of torsional angle as a function of time shows faster sampling for the λ=0 replica than for reference simulation normal MD Christen et al. J. Comput. Chem. 26, 1719 (2005)
34 D = Mass-weighted covariance matrix 1 K m 1 K ( r k 1 r ) K k k! "! ( r k N r )( N r k 1 r ) 1 # k gmx covar S true S = k B 2 ln 1 + k BTe 2! 2 D ( r k 1 r )( 1 r k N r ) N m N K k ( r k N r ) 2 N Here, K is the number of conformations k r 1 is the position of atom/bead number 1 in frame k m 1 is the mass of atom/bead 1 Note similarity to variance! Baron et al J. Phys. Chem. B 110, (2006)
35 gmx rmsf RMSF: structural mobility/flexibility Root-Mean-Square Fluctuation average over time for each atom (or residue) 1 N f k ( r i r ref ) 2 k i Here, Nf is the number of frames in the trajectory ref r i is the position of particle i in the reference structure i r k is the position of particle i in frame k
36 D = S true S = k B 2 ln 1 + k BTe 2 D! 2 Mass-weighted covariance matrix 1 K m 1 K ( r k 1 r ) K k k! "! ( r k N r )( N r k 1 r ) 1 # k ( r k 1 r )( 1 r k N r ) N m N K k ( r k N r ) 2 N 23 Here, K is the number of conformations k r 1 is the position of atom/bead number 1 in frame k m 1 is the mass of atom/bead 1 Note similarity to variance! Baron et al J. Phys. Chem. B 110, (2006)
Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006)
Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006) Martini Workshop 2015 Coarse Graining Basics Alex de Vries Every word or concept, clear as it may seem to be, has only a limited range of applicability
More informationMARTINI simulation details
S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained
More 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 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 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 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 informationFree energy calculations
Free energy calculations Jochen Hub & David van der Spoel Overview Free energies and Probabilities Thermodynamic cycles (Free energy perturbation (FEP)) Thermodynamic integration (TI) (Jarzynski equality
More 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 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 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 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 informationAnalysis of the simulation
Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities
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 informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More 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 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. 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 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 informationMethods for searching and sampling configurational space
International Spring School Statistical Thermodynamics, Santiago de Chile Tuesday, ovember 28, 2017 Lecture 18 Methods for searching and sampling configurational space Prof. Dr. Wilfred F. van Gunsteren
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 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 informationComputing free energy: Thermodynamic perturbation and beyond
Computing free energy: Thermodynamic perturbation and beyond Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls
More 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 informationEntropyCalculationsonaReversiblyFoldingPeptide: ChangesinSoluteFreeEnergyCannotExplain FoldingBehavior
PROTEINS: Structure, Function, and Genetics 43:45 56 (2001) EntropyCalculationsonaReversiblyFoldingPeptide: ChangesinSoluteFreeEnergyCannotExplain FoldingBehavior HeikoSchä fer, 1 XavierDaura, 1 AlanE.Mark,
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 informationJournal of Physics: Conference Series PAPER. To cite this article: Jiramate Kitjanon et al 2017 J. Phys.: Conf. Ser
Journal of Physics: Conference Series PAPER Transferability of Polymer Chain Properties between Coarse-Grained and Atomistic Models of Natural Rubber Molecule Validated by Molecular Dynamics Simulations
More informationMulti-Scale Hierarchical Structure Prediction of Helical Transmembrane Proteins
Multi-Scale Hierarchical Structure Prediction of Helical Transmembrane Proteins Zhong Chen Dept. of Biochemistry and Molecular Biology University of Georgia, Athens, GA 30602 Email: zc@csbl.bmb.uga.edu
More informationIntroduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012
Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012 K. Kremer Max Planck Institute for Polymer Research, Mainz Overview Simulations, general considerations
More informationMartini Basics. Hands on: how to prepare a Martini
Martini Basics Hands on: how to prepare a Martini gin and vermouth are combined at a ratio of 2:1, stirred in a mixing glass with ice cubes, then strained into a chilled cocktail glass and garnished with
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 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 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 informationComputational Science and Engineering Field of Specialisation: Chemistry and Biology. Contact Person: Prof. Philippe Hünenberger
Computational Science and Engineering Field of Specialisation: Chemistry and Biology Contact Person: Prof. Philippe Hünenberger MOLECULAR MODELING/SIMULATION IN CHEMISTRY & BIOLOGY Understand (and predict)
More informationMixing MARTINI: Electrostatic Coupling in Hybrid Atomistic Coarse- Grained Biomolecular Simulations
pubs.acs.org/jpcb Terms of Use Mixing MARTINI: Electrostatic Coupling in Hybrid Atomistic Coarse- Grained Biomolecular Simulations Tsjerk A. Wassenaar, Helgi I. Ingoĺfsson, Marten Prieß, Siewert J. Marrink,
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 informationMD Simulations of Proteins: Practical Hints and Pitfalls to Avoid
MD Simulations of Proteins: Practical ints and Pitfalls to Avoid Stefan Boresch stefan@mdy.univie.ac.at Department of Computational Biological Chemistry Faculty of Chemistry, University of Vienna Vienna
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 informationBiomolecular modeling. Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010)
Biomolecular modeling Marcus Elstner and Tomáš Kubař Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010) IX. FREE ENERGY SIMULATIONS When searching for a physical quantity that is of most
More 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 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 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 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 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 informationImproved Adaptive Resolution Molecular Dynamics Simulation
Improved Adaptive Resolution Molecular Dynamics Simulation Iuliana Marin, Virgil Tudose, Anton Hadar, Nicolae Goga, Andrei Doncescu To cite this version: Iuliana Marin, Virgil Tudose, Anton Hadar, Nicolae
More informationHamiltonian Replica Exchange Molecular Dynamics Using Soft-Core Interactions to Enhance Conformational Sampling
John von Neumann Institute for Computing Hamiltonian Replica Exchange Molecular Dynamics Using Soft-Core Interactions to Enhance Conformational Sampling J. Hritz, Ch. Oostenbrink published in From Computational
More informationPHASE TRANSITIONS IN SOFT MATTER SYSTEMS
OUTLINE: Topic D. PHASE TRANSITIONS IN SOFT MATTER SYSTEMS Definition of a phase Classification of phase transitions Thermodynamics of mixing (gases, polymers, etc.) Mean-field approaches in the spirit
More informationReaction Dynamics (2) Can we predict the rate of reactions?
Reaction Dynamics (2) Can we predict the rate of reactions? Reactions in Liquid Solutions Solvent is NOT a reactant Reactive encounters in solution A reaction occurs if 1. The reactant molecules (A, B)
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 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 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 information3. Solutions W = N!/(N A!N B!) (3.1) Using Stirling s approximation ln(n!) = NlnN N: ΔS mix = k (N A lnn + N B lnn N A lnn A N B lnn B ) (3.
3. Solutions Many biological processes occur between molecules in aqueous solution. In addition, many protein and nucleic acid molecules adopt three-dimensional structure ( fold ) in aqueous solution.
More informationCharacterization of the free-energy landscapes of proteins by NMR-guided metadynamics
Characterization of the free-energy landscapes of proteins by NMR-guided metadynamics B Results and Discussion BIOPHYSICS ND COMPUTTIONL BIOLOGY SI Text Inset B χ χ Inset SI Text Inset C Folding of GB3
More informationAll-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
More informationSUPPLEMENTAL MATERIAL
SUPPLEMENTAL MATERIAL Systematic Coarse-Grained Modeling of Complexation between Small Interfering RNA and Polycations Zonghui Wei 1 and Erik Luijten 1,2,3,4,a) 1 Graduate Program in Applied Physics, Northwestern
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 informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationMolecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods
Molecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods Charles L. Brooks III MMTSB/CTBP 2006 Summer Workshop CHARMM Simulations The flow of data and information
More informationBiomolecular modeling I
2015, December 15 Biomolecular simulation Elementary body atom Each atom x, y, z coordinates A protein is a set of coordinates. (Gromacs, A. P. Heiner) Usually one molecule/complex of interest (e.g. protein,
More informationPrincipal Component Analysis (PCA)
Principal Component Analysis (PCA) Projection onto 2-dim subspace: v 2 v 1 Trajectory x(t) in 3N-dim configuration space A 100-dimensional data set viewed at random orientation Example: Conformational
More informationProtein dynamics. Folding/unfolding dynamics. Fluctuations near the folded state
Protein dynamics Folding/unfolding dynamics Passage over one or more energy barriers Transitions between infinitely many conformations B. Ozkan, K.A. Dill & I. Bahar, Protein Sci. 11, 1958-1970, 2002 Fluctuations
More informationBasic Ingredients of Free Energy Calculations: A Review
Feature Article Basic Ingredients of Free Energy Calculations: A Review CLARA D. CHRIST, 1 ALAN E. MARK, 2 WILFRED F. van GUNSTEREN 1 1 Laboratory of Physical Chemistry, Swiss Federal Institute of Technology,
More information4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro. Karel Berka, Ph.D.
4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro Karel Berka, Ph.D. UP Olomouc, 21.1.-25.1. 2019 Motto A theory is something nobody believes, except the person who made it An experiment
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 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 informationAtomistic Modeling of Small-Angle Scattering Data Using SASSIE-web
Course Introduction Atomistic Modeling of Small-Angle Scattering Data Using SASSIE-web September 21-23, 2016 Advanced Photon Source Argonne National Laboratory, Argonne, IL ccpsas.org Scatters & Simulators
More informationEuropean Life Science Bootcamp: Case Studies in Enhanced Sampling Methods.
European Life Science Bootcamp: Case Studies in Enhanced Sampling Methods davide.branduardi@schrodinger.com Outline Limitations of molecular dynamics Overcoming them with Enhanced sampling Pragmatic use
More informationUnfolding CspB by means of biased molecular dynamics
Chapter 4 Unfolding CspB by means of biased molecular dynamics 4.1 Introduction Understanding the mechanism of protein folding has been a major challenge for the last twenty years, as pointed out in the
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 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 informationSunyia Hussain 06/15/2012 ChE210D final project. Hydration Dynamics at a Hydrophobic Surface. Abstract:
Hydration Dynamics at a Hydrophobic Surface Sunyia Hussain 6/15/212 ChE21D final project Abstract: Water is the universal solvent of life, crucial to the function of all biomolecules. Proteins, membranes,
More informationWhat makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces: Electronic Supplementary
Electronic Supplementary Material (ESI) for Journal of Materials Chemistry B. This journal is The Royal Society of Chemistry 21 What makes a good graphene-binding peptide? Adsorption of amino acids and
More informationCombining an Elastic Network With a Coarse-Grained Molecular Force Field Periole, Xavier; Cavalli, Marco; Marrink, Siewert; Ceruso, Marco A.
University of Groningen Combining an Elastic Network With a Coarse-Grained Molecular Force Field Periole, Xavier; Cavalli, Marco; Marrink, Siewert; Ceruso, Marco A. Published in: Journal of Chemical Theory
More informationComputing free energies with PLUMED 2.0
Computing free energies with PLUMED 2.0 Davide Branduardi Formerly at MPI Biophysics, Frankfurt a.m. TyOutline of the talk Relevance of free energy computation Free-energy from histograms: issues Tackling
More informationOverview & Applications. T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 04 June, 2015
Overview & Applications T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 4 June, 215 Simulations still take time Bakan et al. Bioinformatics 211. Coarse-grained Elastic
More informationAdaptive resolution molecular-dynamics simulation: Changing the degrees of freedom on the fly
THE JOURNAL OF CHEMICAL PHYSICS 123, 224106 2005 Adaptive resolution molecular-dynamics simulation: Changing the degrees of freedom on the fly Matej Praprotnik, a Luigi Delle Site, and Kurt Kremer b Max-Planck-Institut
More informationSupporting Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Supporting Information Lipid molecules can induce an opening of membrane-facing
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 informationSupplementary Figure S1. Urea-mediated buffering mechanism of H. pylori. Gastric urea is funneled to a cytoplasmic urease that is presumably attached
Supplementary Figure S1. Urea-mediated buffering mechanism of H. pylori. Gastric urea is funneled to a cytoplasmic urease that is presumably attached to HpUreI. Urea hydrolysis products 2NH 3 and 1CO 2
More informationMolecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)
Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments) Fei Gao gaofeium@umich.edu Limitations of MD Time scales Length scales (PBC help a lot) Accuracy of
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 informationto satisfy the large number approximations, W W sys can be small.
Chapter 12. The canonical ensemble To discuss systems at constant T, we need to embed them with a diathermal wall in a heat bath. Note that only the system and bath need to be large for W tot and W bath
More informationInsights on Antibiotic Translocation : Molecular Dynamics Simulations
I Insights on Antibiotic Translocation : Molecular Dynamics Simulations Amit Kumar, Eric Hajjar,Enrico Spiga, Francesa Collu, Paolo Ruggerone & Matteo Ceccarelli Marseille : 11 04 2008 Outline METHODS
More informationCoupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation
Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Bo Li Math Dept & CTBP, UCSD Li-Tien Cheng (Math, UCSD) Zhongming Wang (Math & Biochem, UCSD) Yang Xie (MAE,
More informationGenerating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniques
Generating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniques Peter Bolhuis van t Hoff institute for Molecular Sciences University of Amsterdam, The Netherlands
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 28 Electronic Supplementary Information ph-dependent Cooperativity and Existence of
More informationInfluence of temperature and diffusive entropy on the capture radius of fly-casting binding
. Research Paper. SCIENCE CHINA Physics, Mechanics & Astronomy December 11 Vol. 54 No. 1: 7 4 doi: 1.17/s114-11-4485-8 Influence of temperature and diffusive entropy on the capture radius of fly-casting
More informationFree Radical-Initiated Unfolding of Peptide Secondary Structure Elements
Free Radical-Initiated Unfolding of Peptide Secondary Structure Elements Thesis of the Ph.D. Dissertation by Michael C. Owen, M.Sc. Department of Chemical Informatics Faculty of Education University of
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 informationCoarse-grained Models for Oligomer-grafted Silica Nanoparticles
Modeling So+ Ma-er: Linking Mul3ple Length and Time Scales KITP Conference, Santa Barbara, June 4-8, 2012 Coarse-grained Models for Oligomer-grafted Silica Nanoparticles Bingbing Hong Alexandros Chremos
More informationChapter 6 Cyclic urea - a new central unit in bent-core compounds
82 Chapter 6 Cyclic urea - a new central unit in bent-core compounds A new class of five-ring bent-core molecules with a cyclic urea group as a central unit was synthesized [94]. A significant difference
More informationMolecular Simulation Methods with Gromacs
Molecular Simulation Methods with Gromacs CSC 2017 Alex de Vries with special thanks to Tsjerk Wassenaar Hands-on tutorial Multiscaling Simulation and Back-mapping 1 2 Aim Demonstrate how to combine different
More informationPrecipitation. Size! Shape! Size distribution! Agglomeration!
Precipitation Size! Shape! Size distribution! Agglomeration! Precipitation Four major questions: 1. Why do molecules/ions precipitate? 2. What determines the size? 3. What determines the size distribution?
More informationFile ISM02. Dynamics of Soft Matter
File ISM02 Dynamics of Soft Matter 1 Modes of dynamics Quantum Dynamics t: fs-ps, x: 0.1 nm (seldom of importance for soft matter) Molecular Dynamics t: ps µs, x: 1 10 nm Brownian Dynamics t: ns>ps, x:
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2016 Electronic Supplementary Information Active Site Structure and Absorption Spectrum of Channelrhodopsin-2
More informationElectronic excitations in conjugated. Many-Body Green's Functions. Behnaz Bagheri Varnousfaderani. Behnaz Bagheri Varnousfaderani
Technische Universiteit Eindhoven University of Technology Electronic excitations in conjugated polymers Electronic from excitations QM/MM simulations in conjugated pol based from on Many-Body QM/MM simulations
More informationComputing free energy: Replica exchange
Computing free energy: Replica exchange Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls Macroscopic i a t e regime
More informationSupplementary Figures:
Supplementary Figures: Supplementary Figure 1: The two strings converge to two qualitatively different pathways. A) Models of active (red) and inactive (blue) states used as end points for the string calculations
More informationLecture 20. Chemical Potential
Lecture 20 Chemical Potential Reading: Lecture 20, today: Chapter 10, sections A and B Lecture 21, Wednesday: Chapter 10: 10 17 end 3/21/16 1 Pop Question 7 Boltzmann Distribution Two systems with lowest
More information