All-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
|
|
- Ada Washington
- 5 years ago
- Views:
Transcription
1 All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE /27/2010
2 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization Motivations and applications Biomolecular Force Fields
3 Molecular Models Quantum Mechanics physical, but expensive Schrödinger equation: H Ψ = E Ψ wave functions defines electron density Molecular Mechanics less physical --> empirical parameterization cheap and accurate
4
5 Molecular Mechanics Force Field Connectivity (topology) Coordinates Potential energy Parameters Potential Energy Function Functional Form V top, parm r Function of coordinates ( r ) Force Field Dependent on topology and parameters
6 Molecular Mechanics Every atom is a sphere with a radius (Lennard Jones) Point charge is located at each atomic center (Coulomb s law) Bonds and angles are held by springs to ideal lengths eg. V ( ) 2 bond = kb r r 0 Hooke's Law, K b : spring constant, r 0 : ideal length Dihedrals are represented by sigmoidal function which has energy wells at favorable angles. Improper torsions force atoms to be a defined angle to plane. Class I Potential Energy function
7 The "Tinker-toy Model" Bonded terms bonds angles dihedrals (torsions) V bond = k ( r r ) 2 b 0 V = k ( θ ) 2 angle θ = K ( 1+ cos( nχ δ )) θ 0 V dihidral χ improper V coul = q i q ε r i Through space interactions V improper ( ϕ ) 2 = k ϕ ϕ 0 R min i, = V R min LJ i = ε i, + R min 2 R min r i, i, i, = 12 ε ε * ε R min 2 ri, i i, 6
8 Dihedral Term V dihidral χ ( 1+ cos( χ δ )) = K n δ 2/n*π 2*K π (3.14) radians = 180 degrees Molecular Modelling Principles and applications, Leach Pearson Prentice hall second edition (chapter 4) Class I Potential Energy function
9 Lennard-Jones Equation V LJ = ε i, R min r i, i, 12 2 R min r i, i, 6 Energy (kcal/mol (σ, 0) (R min, ε) R min = 2 1/ 6 *σ distance (Å) Molecular Modelling Principles and applications, Leach Pearson Prentice hall second edition (chapter 4) Class I Potential Energy function
10 Potential Energy function V = bonds k b 2 2 ( b b ) + k ( θ θ ) + k ( ϕ ϕ ) 0 angles θ 0 ϕ impropers K dihedrals χ ( 1+ cos( nχ δ )) + N N ε i= 1 = i+ 1 i, R min r i, i, 12 2 R min r i, i, 6 + q i q εr i Different Force-Field CHARMM AMBER GROMOS OPLS Parameterization Experimental observables Quantum Mechanical calculations Interdependences among parameters Molecular Modelling Principles and applications, Leach Pearson Prentice hall second edition (chapter 4) Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry Class I Potential Energy function
11 Applications Dynamics of Molecules how do proteins work Energy: affinity, specificity how do proteins interact with one another and other molecules (drugs, substrates, DNA or RNA) Protein design can proteins be designed to have a specific function
12 Applications Angew. Chem. Int. Ed. 2006, 45,
13 Things to consider size of the system Energy Function sampling method Angew. Chem. Int. Ed. 2006, 45,
14 Limitations to Molecular Mechanics?
15 Limitations to MM MM cannot model the following easily: Chemical reactions Formation and breaking of bonds Polarizability Protonation states
16 Potential Energy Function Electronic Polarizability Combining Rules 1,4 Interactions Lone Pairs All-Atom vs. United Atom Force Fields Treatment of Solvation Treatment of Long-Range Interactions Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
17 Potential Energy Function Electronic Polarizability Combining Rules 1,4 Interactions Lone Pairs All-Atom vs. United Atom Force Fields Treatment of Solvation Treatment of Long-Range Interactions Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
18 Electronic Polarizability Polarization + Drude particle induced dipole models + U pol 1 = µ i E 2 µ i is the dipole moment of atom i E i is the electrostatic field at atom i i i J. Phys. Chem. B 2007, 111, Class II Potential Energy function Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
19 Electronic Polarizability Common force fields use a fixed charge model for every atom Polarizable or non-additive force fields allow for atoms to vary their charges in the context of the electric field Drude particle (Drude oscillator method) every atom has two point charges (atom center and Drude particle) connected by a spring (hooke's law)
20 Potential Energy Function Electronic Polarizability Combining Rules 1,4 Interactions Lone Pairs All-Atom vs. United Atom Force Fields Treatment of Solvation Treatment of Long-Range Interactions Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
21 Combining Rules force fields combine parameters differently arithmetic mean R min i, = R min i + R min 2 0.1, 0.9 geometric mean 0.5, 0.5 ε ε * ε i, = i A1 = A2 V LJ = ε i, R min r i, i, 12 R min r i, i, 6 ε i A1 ε ε i,, A2 ε i,,
22 Combining Rules Arithmetic mean Well depth i, Well depth Radius Radius i, Geometric mean Radius i Well depth ( ) / 2 = 0.5 ( ) / 2 = 0.5 (0.1 * 0.9)^(1 / 2) = 0.3 (0.5 * 0.5)^(1 / 2) = 0.5
23 Combining Rules Harmonic mean Euclidean relationship 2 / (1 / / 0.9) = / (1 / / 0.5) = 0.5 (0.1^ ^2)^(1 / 2) = 0.9 (0.5^ ^2)^(1 / 2) = 0.7
24 Potential Energy Function Electronic Polarizability Combining Rules 1,4 Interactions Lone Pairs All-Atom vs. United Atom Force Fields Treatment of Solvation Treatment of Long-Range Interactions Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
25 1,4 (non-bonded) Interactions 1,2 and 1,3 are normally ignored by force-fields 1,4 are scaled CHARMM no scaling AMBER and OPLS scale 0.83 Bonded terms dominate 2 3 4
26 1,5 interactions Not scaled (s = 1)
27 1,4 interactions Maybe scaled (amber: s = 0.83) Note that 1 / 1.2 = 0.83
28 Potential Energy Function Electronic Polarizability Combining Rules 1,4 Interactions Lone Pairs All-Atom vs. United Atom Force Fields Treatment of Solvation Treatment of Long-Range Interactions Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
29 Lone Pairs Force field without lone pairs do well Adding more parameters Fitting charges from quantum mechanical electrostatic potentials TIP5P and ST2 Water models
30 Potential Energy Function Electronic Polarizability Combining Rules 1,4 Interactions Lone Pairs All-Atom vs. United Atom Force Fields Treatment of Solvation Treatment of Long-Range Interactions Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
31 All-Atom vs. United Atom Force Fields Hydrogen atoms are not represented explicitly. incorporated in to the bonded heavy atom Methyl group will be represented by only one sphere Corse grain models
32 Potential Energy function Electronic Polarizability Combining Rules 1,4 Interactions Lone Pairs All-Atom vs. United Atom Force Fields Treatment of Solvation Treatment of Long-Range Interactions Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
33 Treatment of Solvation Explicit water models TIP3P, TIP4P, TIP5P, SPC,ST2 Implicit water models (Continuum with fixed dielectric) Poisson-Boltzmann (PB) model Generalized Born (GB) model This will be discussed in detail in a following class
34 Treatment of Solvation Explicit water boundaries Angew. Chem. Int. Ed. 2006, 45,
35 Potential Energy Function Electronic Polarizability Combining Rules 1,4 Interactions Lone Pairs All-Atom vs. United Atom Force Fields Treatment of Solvation Treatment of Long-Range Interactions Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
36 Treatment of Long-Range Cut-offs Interactions Partial Mesh Ewald (PME) Smoothing functions Re-parameterization with different treatments of long-range interactions Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
37 Solvation and long range electrostatics Angew. Chem. Int. Ed. 2006, 45,
38 Biomolecular Force Fields Protein Nucleic Acid Lipid Carbohydrates Drug-Like Molecules Heterogeneous Biomolecular Systems (combining different types of biomolecules) Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
39 Conclusions Introduced Molecular Mechanics Described functional form Discussed parameterization Considered dependences and other issues Discussed applications
40 Thank you for your attention. Questions?
41 Protein Force Field Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
42 Force Field Optimization Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
43 Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
44 Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
45 Carbohydrate Force Field Mackerell, Vol. 25, No. 13, Journal of Computational Chemistry
Lecture 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 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 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 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 informationForce Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid
Force Fields for Classical Molecular Dynamics simulations of Biomolecules Emad Tajkhorshid Theoretical and Computational Biophysics Group, Beckman Institute Departments of Biochemistry and Pharmacology,
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 informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
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 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 informationComputer simulation of Biological Systems
UNIVERSITY OF TRENTO Faculty of Mathematical, Physical and Natural Sciences Physics Department Ph.D. Thesis in Physics Computer simulation of Biological Systems Supervisors: Dr. Maurizio Dapor Dr. Giovanni
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 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 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 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 informationForce Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid
Force Fields for Classical Molecular Dynamics simulations of Biomolecules Emad Tajkhorshid Beckman Institute Departments of Biochemistry Center for Biophysics and Computational Biology University of Illinois
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
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 informationBiochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,
Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Course,Informa5on, BIOC%530% GraduateAlevel,discussion,of,the,structure,,func5on,,and,chemistry,of,proteins,and, nucleic,acids,,control,of,enzyma5c,reac5ons.,please,see,the,course,syllabus,and,
More informationSemi Empirical Force Fields and Their Limitations. Potential Energy Surface (PES)
Semi Empirical Force Fields and Their Limitations Ioan Kosztin Beckman Institute University of Illinois at Urbana-Champaign Potential Energy Surface (PES) Schrödinger equation: H T Ψ( r, = E Ψ( r, H =
More informationAtomic and molecular interaction forces in biology
Atomic and molecular interaction forces in biology 1 Outline Types of interactions relevant to biology Van der Waals interactions H-bond interactions Some properties of water Hydrophobic effect 2 Types
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 informationIntroduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA
Introduction to Classical Molecular Dynamics Giovanni Chillemi g.chillemi@cineca.it HPC department, CINECA MD ingredients Coordinates Velocities Force field Topology MD Trajectories Input parameters Analysis
More informationClose agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations
Close agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations Alexandre V. Morozov, Tanja Kortemme, Kiril Tsemekhman, David Baker
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 informationComputational protein design
Computational protein design There are astronomically large number of amino acid sequences that needs to be considered for a protein of moderate size e.g. if mutating 10 residues, 20^10 = 10 trillion sequences
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 informationDeveloping Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model. John Dood Hope College
Developing Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model John Dood Hope College What are MD simulations? Model and predict the structure and dynamics of large macromolecules.
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 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 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 informationMolecular Simulation III
Molecular Simulation III Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Molecular Dynamics Jeffry D. Madura Department of Chemistry & Biochemistry Center
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 informationSubject of the Lecture:
Subject of the Lecture: Conceptual basis for the development of force fields. Implementation/validation Water - a worked example Extensions - combining molecular mechanics and quantum mechanics (QM/MM)
More informationThis semester. Books
Models mostly proteins from detailed to more abstract models Some simulation methods This semester Books None necessary for my group and Prof Rarey Molecular Modelling: Principles and Applications Leach,
More informationChapter 11 Molecular Mechanics
Chapter 11 Molecular Mechanics Molecular Mechanics uses an analytical, differentiable, and relatively simple potential energy function, (R), for describing the interactions between a set of atoms specified
More informationNon-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationNon-bonded interactions
speeding up the number-crunching Marcus Elstner and Tomáš Kubař May 8, 2015 why care? key to understand biomolecular structure and function binding of a ligand efficiency of a reaction color of a chromophore
More informationForce fields in computer simulation of soft nanomaterials
Lecture 2 Part B Force fields in computer simulation of soft nanomaterials Recommended reading: Leach Chapter 4 1 Force Field Methods Force field methods are simulation methods that use classical force
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 informationAtomistic Modeling of DNA and Protein Structures
Atomistic Modeling of DNA and Protein Structures George C. Schatz Introduction In this lecture we introduce one of the most commonly used computational tools for studying protein and DNA structure, the
More informationScuola di Chimica Computazionale
Societa Chimica Italiana Gruppo Interdivisionale di Chimica Computazionale Scuola di Chimica Computazionale Introduzione, per Esercizi, all Uso del Calcolatore in Chimica Organica e Biologica Modellistica
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 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 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 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 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 informationForce Fields for MD simulations
Force Fields for MD simulations Topology/parameter files Where do the numbers an MD code uses come from? ow to make topology files for ligands, cofactors, special amino acids, ow to obtain/develop missing
More informationW tot (z) = W (z) k B T ln e i q iφ mp(z i )/k B T, [1]
Supporting Text Simulation Details. Simulations were carried out with the program charmm (1) using the PARAM27 force field with standard lipid (2), protein (3), and TIP3P water model (4) parameters. The
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 informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz January 17, 2011 CP-PAW/COSMO Interface Mixed quantum/ classical molecular dynamics (QM/MM) Quantum mechanics is computationally expensive
More informationConformational sampling of macrocycles in solution and in the solid state
Conformational sampling of macrocycles in solution and in the solid state Paul Hawkins, Ph.D. Head of Scientific Solutions Stanislaw Wlodek, Ph.D. Senior Scientific Developer 6/6/2018 https://berkonomics.com/?p=2437
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationMichael W. Mahoney Department of Physics, Yale University, New Haven, Connecticut 06520
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 23 15 DECEMBER 2001 Quantum, intramolecular flexibility, and polarizability effects on the reproduction of the density anomaly of liquid water by simple potential
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 informationSolvation and reorganization energies in polarizable molecular and continuum solvents
Solvation and reorganization energies in polarizable molecular and continuum solvents Joel S. Bader CuraGen Corporation, 322 East Main Street, Branford, Connecticut 06405 Christian M. Cortis Department
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 informationGENERALIZED BORN MODELS OF MACROMOLECULARSOLVATION EFFECTS
Annu. Rev. Phys. Chem. 2000. 51:129 52 Copyright c 2000 by Annual Reviews. All rights reserved GENERALIZED BORN MODELS OF MACROMOLECULARSOLVATION EFFECTS Donald Bashford and David A. Case Department of
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 informationTOPOLOGIES AND FORCE FIELD PARAMETERS FOR NITROXIDE SPIN LABELS
TOPOLOGIES AND FORCE FIELD PARAMETERS FOR NITROXIDE SPIN LABELS 1. Introduction The force field used in these studies is the CHARMM19 (Chemistry at HARvard Molecular Mechanics) extended atom force field
More informationAn Informal AMBER Small Molecule Force Field :
An Informal AMBER Small Molecule Force Field : parm@frosst Credit Christopher Bayly (1992-2010) initiated, contributed and lead the efforts Daniel McKay (1997-2010) and Jean-François Truchon (2002-2010)
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 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 informationCHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields
Department of Chemistry and Biochemistry, Concordia University! page 1 of 6 CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields INTRODUCTION The goal of this tutorial
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 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 informationMolecular mechanics. classical description of molecules. Marcus Elstner and Tomáš Kubař. April 29, 2016
classical description of molecules April 29, 2016 Chemical bond Conceptual and chemical basis quantum effect solution of the SR numerically expensive (only small molecules can be treated) approximations
More informationThe split-charge method:
Split charge equilibration: A charge transfer method for electrolytes and other non-metallic materials The split-charge method: Martin H. Müser Materialwissenschaften und Werkstoffwissenschaften Motivation
More informationT6.2 Molecular Mechanics
T6.2 Molecular Mechanics We have seen that Benson group additivities are capable of giving heats of formation of molecules with accuracies comparable to those of the best ab initio procedures. However,
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 informationBIOC 530 Fall, 2011 BIOC 530
Fall, 2011 Course Information Course Description Graduate-level discussion of the structure, function, and chemistry of proteins and nucleic acids, control of enzymatic reactions. Please see the course
More informationStructural and dynamical properties of Polyethylenimine in explicit water at different protonation states: A Molecular Dynamics Study
This journal is The Royal Society of Chemistry Structural and dynamical properties of Polyethylenimine in explicit water at different protonation states: A Molecular Dynamics Study Chandan Kumar Choudhury
More informationVariational Implicit Solvation of Biomolecules: From Theory to Numerical Computations
Variational Implicit Solvation of Biomolecules: From Theory to Numerical Computations Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego CECAM Workshop: New Perspectives
More informationVan 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 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 informationSpecific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate
Specific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate A. Arslanargin, A. Powers, S. Rick, T. Pollard, T. Beck Univ Cincinnati Chemistry Support: NSF, OSC TSRC 2016 November 2, 2016 A. Arslanargin,
More informationpk a Calculations in Solution and Proteins with QM/MM Free Energy Perturbation Simulations: A Quantitative Test of QM/MM Protocols
J. Phys. Chem. B 2005, 109, 17715-17733 17715 pk a Calculations in Solution and Proteins with QM/MM Free Energy Perturbation Simulations: A Quantitative Test of QM/MM Protocols Demian Riccardi, Patricia
More informationPerformance Comparison of Generalized Born and Poisson Methods in the Calculation of Electrostatic Solvation Energies for Protein Structures
Performance Comparison of Generalized Born and Poisson Methods in the Calculation of Electrostatic Solvation Energies for Protein Structures MICHAEL FEIG, 1 ALEXEY ONUFRIEV, 1 MICHAEL S. LEE, 2 WONPIL
More informationPROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS
TASKQUARTERLYvol.20,No4,2016,pp.353 360 PROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS MARTIN ZACHARIAS Physics Department T38, Technical University of Munich James-Franck-Str.
More informationBchem 675 Lecture 9 Electrostatics-Lecture 2 Debye-Hückel: Continued Counter ion condensation
Bchem 675 Lecture 9 Electrostatics-Lecture 2 Debye-Hückel: Continued Counter ion condensation Ion:ion interactions What is the free energy of ion:ion interactions ΔG i-i? Consider an ion in a solution
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 informationResearch Statement. Shenggao Zhou. November 3, 2014
Shenggao Zhou November 3, My research focuses on: () Scientific computing and numerical analysis (numerical PDEs, numerical optimization, computational fluid dynamics, and level-set method for interface
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 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 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 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 informationLevel-Set Variational Solvation Coupling Solute Molecular Mechanics with Continuum Solvent
Level-Set Variational Solvation Coupling Solute Molecular Mechanics with Continuum Solvent Bo Li Department of Mathematics and Center for Theoretical Biological Physics (CTBP) University of California,
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationSolvent Scales. ε α β α: solvent's ability to act as a hydrogen bond-donor to a solute
Solvent Scales ε α β α: solvent's ability to act as a hydrogen bond-donor to a solute Water 78 1.17 0.47 DMS 47 0.00 0.76 DM 37 0.00 0.76 Methanol 33 0.93 0.66 MPA 29 0.00 1.05 Acetone 21 0.08 0.43 Methylene
More informationAvailable online at
Available online at www.sciencedirect.com Chemical Physics Letters 452 (2008) 78 83 www.elsevier.com/locate/cplett Geometry and excitation energy fluctuations of NMA in aqueous solution with CHARMM, AMBER,
More informationBiochemistry 530: Introduction to Structural Biology. Autumn Quarter 2014 BIOC 530
Biochemistry 530: Introduction to Structural Biology Autumn Quarter 2014 Course Information Course Description Graduate-level discussion of the structure, function, and chemistry of proteins and nucleic
More informationThe generalized Born model: its foundation, applications, and limitations
The generalized Born model: its foundation, applications, and limitations Alexey Onufriev September 8, 2010 Abstract The generalized Born (GB) approximation is introduced in the context of the implicit
More informationfile:///biology Exploring Life/BiologyExploringLife04/
Objectives Identify carbon skeletons and functional groups in organic molecules. Relate monomers and polymers. Describe the processes of building and breaking polymers. Key Terms organic molecule inorganic
More informationPolypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation
Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation Jakob P. Ulmschneider and William L. Jorgensen J.A.C.S. 2004, 126, 1849-1857 Presented by Laura L. Thomas and
More informationAn Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. pubs.acs.org/cr An Empirical
More informationG pol = 1 2 where ρ( x) is the charge density at position x and the reaction electrostatic potential
EFFICIENT AND ACCURATE HIGHER-ORDER FAST MULTIPOLE BOUNDARY ELEMENT METHOD FOR POISSON BOLTZMANN ELECTROSTATICS CHANDRAJIT BAJAJ SHUN-CHUAN CHEN Abstract. The Poisson-Boltzmann equation is a partial differential
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 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 informationBack to molecular interac/ons Quantum theory and molecular structure
Back to molecular interac/ons Quantum theory and molecular structure Atoms are arranged in 3D to cons1tute a molecule. Atoms in one molecule are connected by strong covalent bonds that are not easily broken
More informationForce fields and molecular modelling
Force fields and molecular modelling Carine Clavaguéra Laboratoire de Chimie Physique, CNRS & Université Paris Sud, Université Paris Saclay carine.clavaguera@u-psud.fr Carine Clavaguéra Force fields and
More information