Introduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA
|
|
- Phillip Hopkins
- 5 years ago
- Views:
Transcription
1 Introduction to Classical Molecular Dynamics Giovanni Chillemi HPC department, CINECA
2 MD ingredients Coordinates Velocities Force field Topology MD Trajectories Input parameters Analysis
3 Equation of motion The equations that describe the temporal evolution of a physical system is called equation of motion. There are different equations of motions, which characterize the motion with different levels of approximation: Time-dependent Schrödinger's Equation for quantum-mechanical system Newton s Equation for classical-mechanical system Langevin's Equation for stochastic system
4 Newton s Equation of motion Molecules are quantum-mechanical systems whose motion should be described by Schrödinger's Equation. However, technical difficulties make solving Schrödinger's Equation for large systems impractical. Therefore the motion of a molecule is usually approximated by the laws of classical mechanics and by Newton's equation of Motion. In its most simplistic form Newton's second law of motion states: fi = m i ai where mi is the mass of particle i, ai is its acceleration. The force fi is given as the derivative of the potential energy function V: where ri is the position of particle i V fi = ri
5 Potential energy function 2 2 V r, r2,..., rn kbn bn b0 n k n n 0 n bond 2 angle 2 2 k k n cos mn n n n n 0n improper 2 dihedral dihedral Cij 2 Cij rij rij nonbonded pairs ij qi q j + 4 r 0 r ij bonded interactions non bonded interactions
6 Model Degree of freedom Example of predicted properties Considered Removed Quantum mechanic Nucleus, electrons Nucleons Chemistry reaction Polarizable atoms Atoms, dipoles Electrons Binding of charged substrates Non polarizable atoms Solute atoms, solvent atoms Dipoles Conformational transitions Hydration Implicit solvent Solute atoms Solvent atoms Folding topology of macromolecules Classical Molecular Dynamics
7 A Brief History Year System 964 Liquid Argon (Rahman Phys Rev) 974 Water (Rahman J. Chem Phys. ) 977 Small protein in vacuo (Mc Cammon Karplus Nature) 988 First Protein in explicit water (Levitt PNAS ) From 995 Protein-DNA Complexes Membrane Proteins- Complex Systems
8 Bond Stretching Energy 2 kbn bn b0 n.. bond 2 Unique kb and b0 assigned for each bond pair, i.e. C-C, O-H kb is the spring constant of the bond b0 is the bond length at equilibrium
9 Bond Stretching Force rij V bond V bond rij fi kb rij b0 ri rij ri rij f j f i If atom i and j are closer than b0, the bond force separates them If atom i and j are farther than b0, the bond force draws them nearer
10 Bending Energy 2 k n n 0n angle 2 k is the spring constant of the bending 0 is the angle bending at equilibrium Unique parameters for angle bending are assigned to each bonded triplet of atoms based on their types (e.g. C-C-C, C-O-C, C-CH, etc.)
11 Torsional o Dihedral Energy k cos m dihedral n n n n = angle = phase m = number of peaks in a full rotation
12 Improper Dihedral Energy The energy required to deform a group of atoms from its equilibrium angle, xo. Used for tetrahedral or planar groups 2 k n n 0 n improper 2 dihedral Again this system can be modeled by a spring, and the energy is given by the Hookean potential with respect to the planar angle
13 The Hookean potential kb and k broaden or steepen the slope of the parabola The larger the value of k, the more energy is required to deform an angle (or bond) from its equilibrium value
14 Lennard Jones (Van der Waals) interactions Sir John Lennard Jones Johannes Diderik Van der Waals LJ interactions Cij 2 Cij rij rij nonbonded pairs ij non bonded interactions qi q j + 4 r 0 r ij
15 Electrostatic interactions qi and qj are the partial atomic charges for atoms i and j, separated by a distance rij r is the relative dielectric constant Electrostatic Cij 2 Cij rij rij nonbonded pairs ij non bonded interactions qi q j + 4 r 0 r ij Charles Augustin de Coulomb
16 Periodic boundary conditions and cut-off radius To simulate our finite system in liquid conditions, we apply the pbc: i.e. the system box is virtually surrounded in all directions by copy of itself. An atom close to a box border interacts with the atoms in another pbc image. The non-bonded interactions are only calculated between atom pairs closer than a spherical cut-off
17 BOX dimension The edge of cubic box must be large enough to avoid interactions of the solute with itself. Its minimal dimension therefore depends by the chosen cut-off for the non bonded interactions cut-off Edge of the box = 3.2 nm (2.4 nm of water plus the solute size)
18 Electrostatic interactions: Particle Mesh Ewald (PME) The cut-off radius method for electrostatic interactions is particularly inaccurate for charged molecules such as DNA of for dipolar groups such as alpha helices - Short range in the real space - Long range in the Fourier space PME corrects these errors and it helps maintaining short the cut-off in the real space: i.e. the number of atom pairs is reduced Darden et al., Particle mesh Ewald-an NlogN method for Ewald sums in large systems. J. Chem. Phys, 98 (993), 0089
19 Charge groups and atom types A charge group is a neutral charge group composed by several partially charged atoms of a chemical group. Electrostatics can be calculated between charge groups instead that atom pairs Charge groups were first introduced to reduce artifacts in the electrostatics calculation but they can also speed up the calculations; given a pair of water molecules for instance, we only need to determine one atom distance instead of nine (or sixteen for a four-site water model) Note that an atom type is not a physical feature. O8 is defined with a different atom type than O9. In fact, their bond constants with C7 and atomic charges are different
20 Force field V r, r2,..., rn 2 2 k b b k bn n 0n n n 0n bond 2 angle 2 2 k k n cos mn n n n n 0n improper 2 dihedral dihedral Cij 2 Cij rij rij nonbonded pairs ij qi q j + 4 r 0 r ij The potential energy function, together with the parameters required to describe the behavior of different kinds of atoms and bonds (kb, k, k, Cij, ), is called a force field. Several force fields are currently used and the choice depends from the studied system. Some force field are better suited for nucleic acids, for example, while others for membrane proteins
21 Available forcefield in Gromacs (4.6.5). AMBER03 protein, nucleic AMBER94 (Duan et al., J. Comp. Chem. 24, , 2003) 2. AMBER94 force field (Cornell et al., JACS 7, , 995) 3. AMBER96 protein, nucleic AMBER94 (Kollman et al., Acc. Chem. Res. 29, , 996) 4. AMBER99 protein, nucleic AMBER94 (Wang et al., J. Comp. Chem. 2, , 2000) 5. AMBER99SB protein, nucleic AMBER94 (Hornak et al., Proteins 65, , 2006) 6. AMBER99SB-ILDN protein, nucleic AMBER94 (Lindorff-Larsen et al., Proteins 78, , 200) 7. AMBERGS force field (Garcia & Sanbonmatsu, PNAS 99, , 2002) 8. CHARMM27 all-atom force field (with CMAP) - version GROMOS96 43a force field 0. GROMOS96 43a2 force field (improved alkane dihedrals). GROMOS96 45a3 force field (Schuler JCC ) 2. GROMOS96 53a5 force field (JCC 2004 vol 25 pag 656) 3. GROMOS96 53a6 force field (JCC 2004 vol 25 pag 656) 4. GROMOS96 54a7 force field (Eur. Biophys. J. (20), 40,, ) 5. OPLS-AA/L all-atom force field (200 aminoacid dihedrals) 6. [DEPRECATED] Encad all-atom force field, using full solvent charges 7. [DEPRECATED] Encad all-atom force field, using scaled-down vacuum charges 8. [DEPRECATED] Gromacs force field (see manual) 9. [DEPRECATED] Gromacs force field with hydrogens for NMR 2
22 Integration of the equation of motion Numeric integration of Newton's equation of motion is typically done step by step using methods that are called Finite Difference methods. These methods use the information available at time t to predict the system's coordinates and velocities at a time t + t, where t is a short time interval and are based on a Taylor expansion of the position at time t + t r (t t ) r (t ) v(t ) t a (t ) t
23 Integration of the equation of motion r (t t ) 2r (t ) r (t t ) a (t ) t 2 Verlet integrator r (t t ) r (t ) v(t t ) t 2 Leap-frog integrator v(t t ) v(t t ) a(t ) t 2 2 r (t t ) r (t ) v(t ) t a (t ) t 2 Velocity 2 t Verlet v(t t ) v(t ) a (t ) a (t t ) 2
24 Choice of the timestep The length of the timestep must be small compared to the period of the highest frequency motions being simulated Force characteristics Relaxation time (fs) Time step (fs) High frequency motion: bond stretching vibrations Medium frequency motion: angle bending, proper and improper dihedral angle deformation, LJ and short range Coulombian interactions Low frequency motion: long range coulombian interactions The bond stretching vibrations are generally of minimal interest in the study of biomolecular structure and function. Therefore this degree of freedom is usually kept frozen with constraint algorithms (2 x 0-5 s)
25 Constraints The application of geometrical constraints to maintain fixed all the bond lengths, during a simulation, allows the use of a time step up to 2 fs In SHAKE (the first constraints algorithm to be implemented in a MD code) changes a set of unconstrained coordinates to a set of coordinates that fulfill a list of distance constraints, solving a set of Lagrange multipliers in the constrained equations of motion (iterative)
26 Constraints The SETTLE algorithm (developed in 992) is an analytical solution of SHAKE, specifically for water (non-iterative) The LINCS algorithm (developed in 997, twenty years after SHAKE) solves the bond length constraints, always in two steps (non-iterative)
27 Constraints
28 Timescale Protein Folding milliseconds/seconds (0-3-s) Ligand Binding micro/milliseconds ( s) Enzyme catalysis micro/milliseconds ( s) Conformational transitions pico/nanoseconds ( s) Collective vibrations picosecond (0-2 s) Bond vibrations femtosecond (0-5 s)
29 Topology The topology file describes the atoms composing a molecule and their bond connections Es: flexspc.itp in gromacs [ moleculetype ] ; molname nrexcl SOL 2 [ atoms ] ; id at type OW_spc 2 HW_spc 3 HW_spc res nr res name at name cg nr charge mass SOL OW SOL HW SOL HW [ bonds ] ;i j funct length force.c [ angles ] ;i j k funct angle force.c
30 Constraints in Topology Only in case of water, the constraint algorithm can be selected in the topology file flexspc.itp spc.itp
31 Topology
32 Initial velocities The initial velocity of each atom is random assigned through a Maxwell-Boltzmann distribution that is function of the temperature
33 To recapitulate.. Coordinates Velocities Force field Topolog y MD Trajectories Input parameters Analysis
34
35 Molecular Dynamics ensembles The method discussed above is appropriate for the microcanonical ensemble: constant N (number of particles), V (volume) and ET (total energy = E + Ekin) Note that if time step is short enough, the system loses/gains no net energy (potential + kinetic) when running MD in the NVE ensemble When simulating biological macromolecules, it might be more appropriate to simulate under constant Temperature (T) or constant Pressure (P): Canonical ensemble: NVT Isothermal-isobaric: NPT
36 Simulating at constant T: the Berendsen scheme system Heat bath Bath supplies or removes heat from the system as appropriate dt t T0 T t dt T where determines how strong the bath influences the system Exponentially scale the velocities at each time step by the factor : t T0 T T t 2 T: kinetic temperature Berendsen et al. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 8:3684 (984)
37 Simulating at constant T: the Berendsen scheme system Heat bath A small, close to the timestep (strong thermostat), is useful in the equilibration phase, when the quick decreasing of the potential energy could increase too much the kinetic energy of the protein A bigger, e.g. equal to ten times the timestep (weak thermostat), is useful in the production phase, when we want to keep at minimum the perturbation to the conformational sampling Berendsen et al. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 8:3684 (984)
38 Simulating at constant P: the Berendsen scheme system Pressure bath dp t P0 P t dt P A change in pressure P is related to a change in volume V Couple the system to a pressure bath To regulate pressure: exponentially scale the volume of the simulation box at each time step by a factor t Po P t where kt : isothermal compressibility t kt P P : coupling constant Berendsen et al. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 8:3684 (984) 3
39 Sample input file of gromacs title = Yo cpp = /lib/cpp include = -I../top define = integrator = md dt = nsteps = nstxout = 5000 nstvout = 5000 nstlog = 5000 nstenergy = 250 nstxout-compressed = 250 compressed-x-grps = Protein energygrps = Protein SOL nstlist = 0 ns-type = grid rlist = 0.8 coulombtype = cut-off rcoulomb =.4 rvdw = 0.8 tcoupl = Berendsen tc-grps = Protein SOL tau-t = ref-t = Pcoupl = Berendsen tau-p =.0 compressibility = 4.5e-5 ref-p =.0 gen-vel = yes gen-temp = 300 gen-seed = constraints = all-bonds
40 Conformational sampling
41 Energy minimization The potential energy surface of a molecule is defined by only a global minumum and a great number of local minima: i.e. conformations where all the first derivative of the potential energy function with respect to the coordinates are zero and all second derivatives are non-negative
42 The energy minimization algorithms finds the nearest local minimum, i.e. the minimum that can be reached by systematically moving down the steepest local gradient Usually they cannot found the global minumum Two algorithms very used in MD codes are steepest descent conjugate gradient Both of the first order: they use the first derivative of the energy potential function with respect to the coordinates
43 The steepest descent uses only the gradient of the potential energy function in the local position to calculate the coordinate movement It is quicker in the single iteration but less precise in finding the local minumum. Therefore is useful in the first steps of minimizations The conjugate gradient uses also the gradient of the potential energy function in the previous step It is more accurate in finding the local minumum but it is slower than the steepest descent. Therefore is usually applied after some steps of steepest descent
44 Starting from A, the steepest descents goes through A-B-C (or A-D-F) While the conjugate gradient, weighting the A-B and B-C gradients, goes through A-B-O
45 GROMACS USER MANUAL
46 -4 interactions Atoms covalently bound are defined as first neighbours second neighbours and so on. LJ and electrostatic interactions are not calculated among first and second neighbours since they are considered in the stretching (first) or in the bending potential (second) The standard non-bonding interactions are too strong for the third neighbours and are reduced (interactions -4; list -4)
Molecular Dynamics in practice with GROMACS
Molecular Dynamics in practice with GROMACS GROMACS is one of the wold s fastest software package for molecular dynamics simulations. One can find many helpful materials, manual as well as to download
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 informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
More informationGeneration of topology files of a protein chain and simulations of a dipeptide
Chapter-36 Generation of topology files of a protein chain and simulations of a dipeptide After studying this chapter, you will be able to 1) Generate a topology file of a dipeptide molecule. 2) Generate
More informationfreely available at
freely available at www.gromacs.org Generally 3 to 10 times faster than other Molecular Dynamics programs Very user-friendly: issues clear error messages, no scripting language is required to run the programs,
More informationWhy study protein dynamics?
Why study protein dynamics? Protein flexibility is crucial for function. One average structure is not enough. Proteins constantly sample configurational space. Transport - binding and moving molecules
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationChapter 32. Computations of radial distributions functions, PMFs and diffusion constants
Chapter 32 Computations of radial distributions functions, PMFs and diffusion constants After studying this chapter, you will be able to perform the 1) Calculations of radial distribution functions between
More informationForce fields, thermo- and barostats. Berk Hess
Force fields, thermo- and barostats Berk Hess What is a force field? A force field usually consists of three parts: a set of functional forms parameters for the functional forms that, usually, depend on
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 informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
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 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 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 informationGromacs Workshop Spring CSC
Gromacs Workshop Spring 2007 @ CSC Erik Lindahl Center for Biomembrane Research Stockholm University, Sweden David van der Spoel Dept. Cell & Molecular Biology Uppsala University, Sweden Berk Hess Max-Planck-Institut
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 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 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 informationPotential Energy (hyper)surface
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 informationRunning MD on HPC architectures I. Hybrid Clusters
Running MD on HPC architectures I. Hybrid Clusters Alessandro Grottesi Cineca Bologna PATC 2015 Today's lecture You will learn: Gromacs @ CINECA: set up and launch of simulations Launch MD code (GROMACS,
More informationk θ (θ θ 0 ) 2 angles r i j r i j
1 Force fields 1.1 Introduction The term force field is slightly misleading, since it refers to the parameters of the potential used to calculate the forces (via gradient) in molecular dynamics simulations.
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 informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More 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 informationMolecular Modeling Lab #2: Using Molecular Dynamics Simulations To Compare the Stability of Polylysine in Different Protonation States
Molecular Modeling Lab #2: Using Molecular Dynamics Simulations To Compare the Stability of Polylysine in Different Protonation States For this lab report, answer the following questions. Be sure to support
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 informationPrinciples and Applications of Molecular Dynamics Simulations with NAMD
Principles and Applications of Molecular Dynamics Simulations with NAMD Nov. 14, 2016 Computational Microscope NCSA supercomputer JC Gumbart Assistant Professor of Physics Georgia Institute of Technology
More informationComputational Chemistry - MD Simulations
Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated
More 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 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 informationThe GROMOS Software for (Bio)Molecular Simulation
The GROMOS Software for (Bio)Molecular Simulation Volume 2: Algorithms and Formulae for Modelling of Molecular Systems September 28, 2017 Contents Chapter 1. Introduction 2-1 Chapter 2. Basic choices
More informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
More informationMolecular Mechanics, Dynamics & Docking
Molecular Mechanics, Dynamics & Docking Lawrence Hunter, Ph.D. Director, Computational Bioscience Program University of Colorado School of Medicine Larry.Hunter@uchsc.edu http://compbio.uchsc.edu/hunter
More informationUnit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement
Unit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement Xin Yin, Yeqi Shi, Yanbing Wei, Yongho Joo, Padma Gopalan, Izabela Szlufarska,
More informationBiomolecular modeling I
2016, December 6 Biomolecular structure Structural elements of life Biomolecules proteins, nucleic acids, lipids, carbohydrates... Biomolecular structure Biomolecules biomolecular complexes aggregates...
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 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 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 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 informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
More informationMolecular Simulation II. Classical Mechanical Treatment
Molecular Simulation II Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationAdvanced Molecular Molecular Dynamics
Advanced Molecular Molecular Dynamics Technical details May 12, 2014 Integration of harmonic oscillator r m period = 2 k k and the temperature T determine the sampling of x (here T is related with v 0
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 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 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 informationThe Molecular Dynamics Simulation Process
The Molecular Dynamics Simulation Process For textbooks see: M.P. Allen and D.J. Tildesley. Computer Simulation of Liquids.Oxford University Press, New York, 1987. D. Frenkel and B. Smit. Understanding
More informationHomology modeling. Dinesh Gupta ICGEB, New Delhi 1/27/2010 5:59 PM
Homology modeling Dinesh Gupta ICGEB, New Delhi Protein structure prediction Methods: Homology (comparative) modelling Threading Ab-initio Protein Homology modeling Homology modeling is an extrapolation
More informationIntroduction to Molecular Dynamics Simulations for Biological Systems
Introduction to Molecular Dynamics Simulations for Biological Systems Ángel Piñeiro 27th November 2009 Contents 1 Introduction 2 2 Historical Evolution 2 3 Theory 3 3.1 Initial Conditions...............................
More informationDihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769
Dihedral Angles Homayoun Valafar Department of Computer Science and Engineering, USC The precise definition of a dihedral or torsion angle can be found in spatial geometry Angle between to planes Dihedral
More informationComputational Molecular Biophysics. Computational Biophysics, GRS Jülich SS 2013
Computational Molecular Biophysics Computational Biophysics, GRS Jülich SS 2013 Computational Considerations Number of terms for different energy contributions / E el ~N ~N ~N ~N 2 Computational cost for
More informationMolecular Mechanics / ReaxFF
Molecular Dynamics simulations Lecture 09: Molecular Mechanics / ReaxFF Dr. Olli Pakarinen University of Helsinki Fall 2012 Lecture notes based on notes by Dr. Jani Kotakoski, 2010 CONTENTS Molecular mechanics
More informationUniversity of Groningen. Characterization of oil/water interfaces van Buuren, Aldert Roelf
University of Groningen Characterization of oil/water interfaces van Buuren, Aldert Roelf IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it.
More informationForce Field for Water Based on Neural Network
Force Field for Water Based on Neural Network Hao Wang Department of Chemistry, Duke University, Durham, NC 27708, USA Weitao Yang* Department of Chemistry, Duke University, Durham, NC 27708, USA Department
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the
More informationApplications of Molecular Dynamics
June 4, 0 Molecular Modeling and Simulation Applications of Molecular Dynamics Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo Tohru
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 informationThe Computational Microscope
The Computational Microscope Computational microscope views at atomic resolution... Rs SER RER C E M RER N GA L... how living cells maintain health and battle disease John Stone Our Microscope is Made
More 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 informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationThe micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction
ISBN 978-1-84626-081-0 Proceedings of the 2010 International Conference on Application of Mathematics and Physics Volume 1: Advances on Space Weather, Meteorology and Applied Physics Nanjing, P. R. China,
More informationLecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects
Statistical Thermodynamics Lecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects Dr. Ronald M. Levy ronlevy@temple.edu Bare Molecular Mechanics Atomistic Force Fields: torsion
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 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 informationMolecular Dynamics, Monte Carlo and Docking. Lecture 21. Introduction to Bioinformatics MNW2
Molecular Dynamics, Monte Carlo and Docking Lecture 21 Introduction to Bioinformatics MNW2 If you throw up a stone, it is Physics. If you throw up a stone, it is Physics. If it lands on your head, it is
More informationSupporting Information
Supporting Information Structure and Dynamics of Uranyl(VI) and Plutonyl(VI) Cations in Ionic Liquid/Water Mixtures via Molecular Dynamics Simulations Katie A. Maerzke, George S. Goff, Wolfgang H. Runde,
More informationMolecular Dynamics. Coherency of Molecular Motions. Local Minima & Global Minimum. Theoretical Basis of Molecular Dynamic Calculations
Molecular Dynamics Molecular Dynamics provide a good approach to determine the preferred conformers and global minimum of a molecule Theoretical Basis of Molecular Dynamic Calculations Molecular dynamics
More informationHands-on : Model Potential Molecular Dynamics
Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration
More informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
More informationThe potential energy function is typically called the force field (although it is really an energy function, force is a negative gradient of energy).
7. Molecular Mechanics In molecular mechanics (MM) a molecule is described as a collection of atoms that obey the laws of classical mechanics. The potential energy of a molecule is a simple analytical
More informationExploring the Free Energy Surface of Short Peptides by Using Metadynamics
John von Neumann Institute for Computing Exploring the Free Energy Surface of Short Peptides by Using Metadynamics C. Camilloni, A. De Simone published in From Computational Biophysics to Systems Biology
More informationAdvanced in silico drug design
Advanced in silico drug design RNDr. Martin Lepšík, Ph.D. Lecture: Advanced scoring Palacky University, Olomouc 2016 1 Outline 1. Scoring Definition, Types 2. Physics-based Scoring: Master Equation Terms
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 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 informationThe Potential Energy Surface (PES) And the Basic Force Field Chem 4021/8021 Video II.iii
The Potential Energy Surface (PES) And the Basic Force Field Chem 4021/8021 Video II.iii Fundamental Points About Which to Be Thinking It s clear the PES is useful, so how can I construct it for an arbitrary
More informationMolecular Dynamics. Molecules in motion
Molecular Dynamics Molecules in motion 1 Molecules in mo1on Molecules are not sta1c, but move all the 1me Source: h9p://en.wikipedia.org/wiki/kine1c_theory 2 Gasses, liquids and solids Gasses, liquids
More information1 Bulk Simulations in a HCP lattice
1 Bulk Simulations in a HCP lattice 1.1 Introduction This project is a continuation of two previous projects that studied the mechanism of melting in a FCC and BCC lattice. The current project studies
More informationSupporting Information for: Physics Behind the Water Transport through. Nanoporous Graphene and Boron Nitride
Supporting Information for: Physics Behind the Water Transport through Nanoporous Graphene and Boron Nitride Ludovic Garnier, Anthony Szymczyk, Patrice Malfreyt, and Aziz Ghoufi, Institut de Physique de
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 informationMedical Research, Medicinal Chemistry, University of Leuven, Leuven, Belgium.
Supporting Information Towards peptide vaccines against Zika virus: Immunoinformatics combined with molecular dynamics simulations to predict antigenic epitopes of Zika viral proteins Muhammad Usman Mirza
More informationA Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane
A Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane Supplementary Information: Simulation Procedure and Physical Property Analysis Simulation Procedure The molecular
More informationNWChem: Molecular Dynamics and QM/MM
NWChem: Molecular Dynamics and QM/MM Molecular Dynamics Functionality Target systems: General features: biomolecules (proteins, DNA/RNA, biomembranes) energy evaluation (SP) energy minimization (EM) molecular
More informationGenerate topology (Gromacs)
Running MD code What for...? 1. Minimizzation 2. Molecular Dynamics (classic, brownian, Langevin) 3. Normal Mode Analysis 4. Essential Dynamics and Sampling 5. Free Energy calculations (FEP, Umbrella sampling,
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 information3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D.
3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D. Thierry Langer, Ph.D. Jana Vrbková, Ph.D. UP Olomouc, 23.1.-26.1. 2018
More informationMolecular dynamics simulations of a single stranded (ss) DNA
Molecular dynamics simulations of a single stranded (ss) DNA Subhasish Chatterjee 1, Bonnie Gersten 1, Siddarth Thakur 2, Alexander Burin 2 1 Department of Chemistry, Queens College and the Graduate Center
More informationWater models in classical simulations
Water models in classical simulations Maria Fyta Institut für Computerphysik, Universität Stuttgart Stuttgart, Germany Water transparent, odorless, tasteless and ubiquitous really simple: two H atoms attached
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
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 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 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 informationSimulation of molecular systems by molecular dynamics
Simulation of molecular systems by molecular dynamics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Dynamics, Label RFCT 2015 November 26, 2015 1 / 35 Introduction
More informationSupplemental Material for Global Langevin model of multidimensional biomolecular dynamics
Supplemental Material for Global Langevin model of multidimensional biomolecular dynamics Norbert Schaudinnus, Benjamin Lickert, Mithun Biswas and Gerhard Stock Biomolecular Dynamics, Institute of Physics,
More 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 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 informationMolecular Dynamics Studies of Human β-glucuronidase
American Journal of Applied Sciences 7 (6): 83-88, 010 ISSN 1546-939 010Science Publications Molecular Dynamics Studies of Human β-lucuronidase Ibrahim Ali Noorbatcha, Ayesha Masrur Khan and Hamzah Mohd
More informationParallel Utility for Modeling of Molecular Aggregation
Parallel Utility for Modeling of Molecular Aggregation Peter Spijker, Bart Markvoort and Peter Hilbers BioModeling and bioinformatics Eindhoven University of Technology Eindhoven,The Netherlands p.spijker@pumma.nl
More informationParticle-Based Simulation of Bio-Electronic Systems
Particle-Based Simulation of Bio-Electronic Systems Alex Smolyanitsky, and Marco Saraniti Arizona State University Outline Particle-based Brownian dynamics simulations for bioelectronic systems Complex-field
More informationCitation for published version (APA): Hess, B. (2002). Stochastic concepts in molecular simulation Groningen: s.n.
University of Groningen Stochastic concepts in molecular simulation Hess, Berk IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check
More informationMolecular Dynamics Simulation of Methanol-Water Mixture
Molecular Dynamics Simulation of Methanol-Water Mixture Palazzo Mancini, Mara Cantoni University of Urbino Carlo Bo Abstract In this study some properties of the methanol-water mixture such as diffusivity,
More information