Molecular Dynamics, Monte Carlo and Docking. Lecture 21. Introduction to Bioinformatics MNW2

Size: px
Start display at page:

Download "Molecular Dynamics, Monte Carlo and Docking. Lecture 21. Introduction to Bioinformatics MNW2"

Transcription

1 Molecular Dynamics, Monte Carlo and Docking Lecture 21 Introduction to Bioinformatics MNW2

2

3

4 Allowed phi-psi angles Red areas are preferred, yellow areas are allowed, and white is avoided

5

6

7

8

9 2.3a Hamiltonian equations of motion Hamiltonian equations (one degree of freedom): H Hamiltonian function, Hamiltonian, q, p Canonical variables: generalized coordinate (q) and momentum [impulses] (p).

10 q = coordinates p = momentum

11 v(t) = (r(t + t) - r(t - t))/2 t

12

13 Molecular Dynamics Knowledge of the atomic forces and masses can be used to solve the position of each atom along a series of extremely small time steps (on the order of femtoseconds = seconds). The resulting series of snapshots of structural changes over time is called a trajectory. The use of this method to compute trajectories can be more easily seen when Newton's equation is expressed in the following form: The "leapfrog" method is a common numerical approach to calculating trajectories based on Newton's equation. The steps can be summarized as follows:

14

15 Force field The potential energy of a system can be expressed as a sum of valence (or bond), crossterm, and nonbond interactions: The energy of valence interactions comprises bond stretching (E bond ), valence angle bending (E angle ), dihedral angle torsion (E torsion ), and inversion (also called out-ofplane interactions) (E inversion or E oop ) terms, which are part of nearly all forcefields for covalent systems. A Urey-Bradley term (E UB ) may be used to account for interactions between atom pairs involved in 1-3 configurations (i.e., atoms bound to a common atom): E valence = E bond + E angle + E torsion + E oop + E UB Modern (second-generation) forcefields include cross terms to account for such factors as bond or angle distortions caused by nearby atoms. Crossterms can include the following terms: stretch-stretch, stretch-bend-stretch, bend-bend, torsion-stretch, torsion-bend-bend, bend-torsion-bend, stretch-torsion-stretch. The energy of interactions between nonbonded atoms is accounted for by van der Waals (E vdw ), electrostatic (E Coulomb ), and (in some older forcefields) hydrogen bond (E hbond ) terms: E nonbond = E vdw + E Coulomb + E hbond

16 Force field

17 energy distance f = a/r 6 -b/r 12 Van der Waals forces

18 F = kq1q1/r2

19

20 Figure: Snapshots of ubiquitin pulling with constant velocity at three different time steps.

21

22

23

24

25

26

27

28

29

30 antibody HyHEL-63 (cyan) complexed with Hen Egg White Lysozyme The X-ray structure of the antibody HyHEL-63 (cyan) uncomplexed and complexed with Hen Egg White Lysozyme (yellow) has shown that there are small but significant, local conformational changes in the antibody paratope on binding. The structure also reveals that most of the charged epitope residues face the antibody. Details are in Li YL, Li HM, Smith-Gill SJ and Mariuzza RA (2000) The conformations of the X-ray structure Three-dimensional structures of the free and antigen-bound Fab from monoclonal antilysozyme antibody HyHEL-63. Biochemistry 39: Salt links and electrostatic interactions provide much of the free energy of binding. Most of the charged residues face in interface in the X-ray structure. The importance of the salt link between Lys97 of HEL and Asp27 of the antibody heavy chain is revealed by molecular dynamics simulations. After 1NSec of MD simulation at 100 C the overall conformation of the complex has changed, but the salt link persists. Details are described in Sinha N and Smith-Gill SJ (2002) Electrostatics in protein binding and function. Current Protein & Peptide Science 3:

Molecular Dynamics, Monte Carlo and Docking. Lecture 21. Introduction to Bioinformatics MNW2

Molecular 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 information

Today s lecture. Molecular Mechanics and docking. Lecture 22. Introduction to Bioinformatics Docking - ZDOCK. Protein-protein docking

Today s lecture. Molecular Mechanics and docking. Lecture 22. Introduction to Bioinformatics Docking - ZDOCK. Protein-protein docking C N F O N G A V B O N F O M A C S V U Molecular Mechanics and docking Lecture 22 ntroduction to Bioinformatics 2007 oday s lecture 1. Protein interaction and docking a) Zdock method 2. Molecular motion

More information

Experimentally solving protein structures, protein-protein interactions and simulating protein dynamics. Lecture 15

Experimentally solving protein structures, protein-protein interactions and simulating protein dynamics. Lecture 15 Experimentally solving protein structures, protein-protein interactions and simulating protein dynamics C E N T R F O R I N T B I O I N F E Lecture 15 Introduction to Bioinformatics 2007 E G R A T I V

More information

Structural Bioinformatics (C3210) Molecular Mechanics

Structural 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 information

Why 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. 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 information

Why study protein dynamics?

Why 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 information

Dihedral 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 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 information

THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION

THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION AND CALIBRATION Calculation of turn and beta intrinsic propensities. A statistical analysis of a protein structure

More information

Differences in Electrostatic Properties at Antibody Antigen Binding Sites: Implications for Specificity and Cross-Reactivity

Differences in Electrostatic Properties at Antibody Antigen Binding Sites: Implications for Specificity and Cross-Reactivity 2946 Biophysical Journal Volume 83 December 2002 2946 2968 Differences in Electrostatic Properties at Antibody Antigen Binding Sites: Implications for Specificity and Cross-Reactivity Neeti Sinha, Srinivasan

More information

Gerd Krause, Structural bioinformatics and protein design Leibniz-Institute of molecular Pharmacology

Gerd Krause, Structural bioinformatics and protein design Leibniz-Institute of molecular Pharmacology Basics of Molecular modelling _2 History, Interdisciplinary classification, Goals, Model definition, Molecular mechanics forcefields, Energy minimization algorithms Molecular dynamics simulation Comparative

More information

Molecular Dynamics Simulation of HIV-1 Reverse. Transcriptase

Molecular Dynamics Simulation of HIV-1 Reverse. Transcriptase Molecular Dynamics Simulation of HIV-1 Reverse Transcriptase Abderrahmane Benghanem 1 Maria Kurnikova 2 1 Rensselaer Polytechnic Institute, Troy, NY 2 Carnegie Mellon University, Pittsburgh, PA June 16,

More information

DISCRETE 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: 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 information

Molecular 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 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 information

Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,

Biochemistry,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 information

Computational Biology & Computational Medicine

Computational Biology & Computational Medicine Computational Biology & Computational Medicine Homayoun Valafar Outline Why proteins? What are proteins? How do we compute them? How do we use computational approaches? Why Proteins? Molecular basis of

More information

CE 530 Molecular Simulation

CE 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 information

Semi Empirical Force Fields and Their Limitations. Potential Energy Surface (PES)

Semi 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 information

Biochemistry Prof. S. DasGupta Department of Chemistry Indian Institute of Technology Kharagpur. Lecture - 06 Protein Structure IV

Biochemistry Prof. S. DasGupta Department of Chemistry Indian Institute of Technology Kharagpur. Lecture - 06 Protein Structure IV Biochemistry Prof. S. DasGupta Department of Chemistry Indian Institute of Technology Kharagpur Lecture - 06 Protein Structure IV We complete our discussion on Protein Structures today. And just to recap

More information

Homology modeling. Dinesh Gupta ICGEB, New Delhi 1/27/2010 5:59 PM

Homology 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 information

Molecular Modelling. part of Bioinformatik von RNA- und Proteinstrukturen. Sonja Prohaska. Leipzig, SS Computational EvoDevo University Leipzig

Molecular Modelling. part of Bioinformatik von RNA- und Proteinstrukturen. Sonja Prohaska. Leipzig, SS Computational EvoDevo University Leipzig part of Bioinformatik von RNA- und Proteinstrukturen Computational EvoDevo University Leipzig Leipzig, SS 2011 Protein Structure levels or organization Primary structure: sequence of amino acids (from

More information

Molecular 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 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 information

Molecular Mechanics, Dynamics & Docking

Molecular 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 information

An introduction to Molecular Dynamics. EMBO, June 2016

An 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 information

Potential Energy (hyper)surface

Potential 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 information

Other Cells. Hormones. Viruses. Toxins. Cell. Bacteria

Other Cells. Hormones. Viruses. Toxins. Cell. Bacteria Other Cells Hormones Viruses Toxins Cell Bacteria ΔH < 0 reaction is exothermic, tells us nothing about the spontaneity of the reaction Δ H > 0 reaction is endothermic, tells us nothing about the spontaneity

More information

HOMOLOGY MODELING. The sequence alignment and template structure are then used to produce a structural model of the target.

HOMOLOGY MODELING. The sequence alignment and template structure are then used to produce a structural model of the target. HOMOLOGY MODELING Homology modeling, also known as comparative modeling of protein refers to constructing an atomic-resolution model of the "target" protein from its amino acid sequence and an experimental

More information

Bioengineering 215. An Introduction to Molecular Dynamics for Biomolecules

Bioengineering 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 information

Conformational Geometry of Peptides and Proteins:

Conformational Geometry of Peptides and Proteins: Conformational Geometry of Peptides and Proteins: Before discussing secondary structure, it is important to appreciate the conformational plasticity of proteins. Each residue in a polypeptide has three

More information

Multiscale Materials Modeling

Multiscale Materials Modeling Multiscale Materials Modeling Lecture 09 Quantum Mechanics/Molecular Mechanics (QM/MM) Techniques Fundamentals of Sustainable Technology These notes created by David Keffer, University of Tennessee, Knoxville,

More information

The protein folding problem consists of two parts:

The 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 information

Protein Folding experiments and theory

Protein 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 information

Biomolecules: lecture 9

Biomolecules: lecture 9 Biomolecules: lecture 9 - understanding further why amino acids are the building block for proteins - understanding the chemical properties amino acids bring to proteins - realizing that many proteins

More information

Biomolecules are dynamic no single structure is a perfect model

Biomolecules 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 information

Protein Structure Basics

Protein Structure Basics Protein Structure Basics Presented by Alison Fraser, Christine Lee, Pradhuman Jhala, Corban Rivera Importance of Proteins Muscle structure depends on protein-protein interactions Transport across membranes

More information

Molecular Mechanics. I. Quantum mechanical treatment of molecular systems

Molecular 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 information

Supplement information

Supplement information Electronic Supplementary Material (ESI) for Physil Chemistry Chemil Physics. This journal is the Owner Societies 216 Supplement information Fullerenol C 6 (OH) 16 prevents amyloid fibrillization of Aβ

More information

Assessing the potential of atomistic molecular dynamics simulations to probe reversible. protein-protein recognition and binding

Assessing the potential of atomistic molecular dynamics simulations to probe reversible. protein-protein recognition and binding Supporting information for Assessing the potential of atomistic molecular dynamics simulations to probe reversible protein-protein recognition and binding Luciano A. Abriata and Matteo Dal Peraro Laboratory

More information

Bulk behaviour. Alanine. FIG. 1. Chemical structure of the RKLPDA peptide. Numbers on the left mark alpha carbons.

Bulk behaviour. Alanine. FIG. 1. Chemical structure of the RKLPDA peptide. Numbers on the left mark alpha carbons. Bulk behaviour To characterise the conformational behaviour of the peptide, first we looked at the statistics of alpha carbons and the torsion angles. Then they were correlated with positions of positively

More information

The Structure and Functions of Proteins

The Structure and Functions of Proteins Wright State University CORE Scholar Computer Science and Engineering Faculty Publications Computer Science and Engineering 2003 The Structure and Functions of Proteins Dan E. Krane Wright State University

More information

Useful background reading

Useful background reading Overview of lecture * General comment on peptide bond * Discussion of backbone dihedral angles * Discussion of Ramachandran plots * Description of helix types. * Description of structures * NMR patterns

More information

1 Contents. 1. Contents i. 1. Introduction Forcefields 9

1 Contents. 1. Contents i. 1. Introduction Forcefields 9 1 Contents 1. Contents i 1. Introduction 1 Who should use this documentation.................. 2 What can simulation engines do?.................... 2 Energy minimization........................... 2 Molecular

More information

Lecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability

Lecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Lecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions Van der Waals Interactions

More information

NWChem: Molecular Dynamics and QM/MM

NWChem: 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 information

Introduction The gramicidin A (ga) channel forms by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet. Th

Introduction The gramicidin A (ga) channel forms by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet. Th Abstract When conductive, gramicidin monomers are linked by six hydrogen bonds. To understand the details of dissociation and how the channel transits from a state with 6H bonds to ones with 4H bonds or

More information

Lecture 15: Enzymes & Kinetics. Mechanisms ROLE OF THE TRANSITION STATE. H-O-H + Cl - H-O δ- H Cl δ- HO - + H-Cl. Margaret A. Daugherty.

Lecture 15: Enzymes & Kinetics. Mechanisms ROLE OF THE TRANSITION STATE. H-O-H + Cl - H-O δ- H Cl δ- HO - + H-Cl. Margaret A. Daugherty. Lecture 15: Enzymes & Kinetics Mechanisms Margaret A. Daugherty Fall 2004 ROLE OF THE TRANSITION STATE Consider the reaction: H-O-H + Cl - H-O δ- H Cl δ- HO - + H-Cl Reactants Transition state Products

More information

Molecular Mechanics. Yohann Moreau. November 26, 2015

Molecular 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 information

Elbow Flexibility and Ligand-Induced Domain Rearrangements in Antibody Fab NC6.8: Large Effects of a Small Hapten

Elbow Flexibility and Ligand-Induced Domain Rearrangements in Antibody Fab NC6.8: Large Effects of a Small Hapten 614 Biophysical Journal Volume 79 August 2000 614 628 Elbow Flexibility and Ligand-Induced Domain Rearrangements in Antibody Fab NC6.8: Large Effects of a Small Hapten Christoph A. Sotriffer, Bernd M.

More information

Supplementary Information

Supplementary Information Supplementary Information Resveratrol Serves as a Protein-Substrate Interaction Stabilizer in Human SIRT1 Activation Xuben Hou,, David Rooklin, Hao Fang *,,, Yingkai Zhang Department of Medicinal Chemistry

More information

Molecular Interactions F14NMI. Lecture 4: worked answers to practice questions

Molecular 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 information

Multi-scale approaches in description and design of enzymes

Multi-scale approaches in description and design of enzymes Multi-scale approaches in description and design of enzymes Anastassia Alexandrova and Manuel Sparta UCLA & CNSI Catalysis: it is all about the barrier The inside-out protocol: Big Aim: development of

More information

Using Higher Calculus to Study Biologically Important Molecules Julie C. Mitchell

Using Higher Calculus to Study Biologically Important Molecules Julie C. Mitchell Using Higher Calculus to Study Biologically Important Molecules Julie C. Mitchell Mathematics and Biochemistry University of Wisconsin - Madison 0 There Are Many Kinds Of Proteins The word protein comes

More information

Dominant Paths in Protein Folding

Dominant Paths in Protein Folding Dominant Paths in Protein Folding Henri Orland SPhT, CEA-Saclay France work in collaboration with P. Faccioli, F. Pederiva, M. Sega University of Trento Henri Orland Annecy meeting 2006 Outline Basic notions

More information

Destruction of Amyloid Fibrils by Graphene through Penetration and Extraction of Peptides

Destruction of Amyloid Fibrils by Graphene through Penetration and Extraction of Peptides Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2015 Destruction of Amyloid Fibrils by Graphene through Penetration and Extraction of Peptides Zaixing

More information

Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation

Polypeptide 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 information

Proteins polymer molecules, folded in complex structures. Konstantin Popov Department of Biochemistry and Biophysics

Proteins polymer molecules, folded in complex structures. Konstantin Popov Department of Biochemistry and Biophysics Proteins polymer molecules, folded in complex structures Konstantin Popov Department of Biochemistry and Biophysics Outline General aspects of polymer theory Size and persistent length of ideal linear

More information

April, The energy functions include:

April, The energy functions include: REDUX A collection of Python scripts for torsion angle Monte Carlo protein molecular simulations and analysis The program is based on unified residue peptide model and is designed for more efficient exploration

More information

Molecular Simulation III

Molecular 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 information

Analysis of MD trajectories in GROMACS David van der Spoel

Analysis 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 information

Methods of Computer Simulation. Molecular Dynamics and Monte Carlo

Methods 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 information

ALL LECTURES IN SB Introduction

ALL LECTURES IN SB Introduction 1. Introduction 2. Molecular Architecture I 3. Molecular Architecture II 4. Molecular Simulation I 5. Molecular Simulation II 6. Bioinformatics I 7. Bioinformatics II 8. Prediction I 9. Prediction II ALL

More information

Packing of Secondary Structures

Packing of Secondary Structures 7.88 Lecture Notes - 4 7.24/7.88J/5.48J The Protein Folding and Human Disease Professor Gossard Retrieving, Viewing Protein Structures from the Protein Data Base Helix helix packing Packing of Secondary

More information

Chemistry Problem Set #9 Due on Thursday 11/15/18 in class.

Chemistry Problem Set #9 Due on Thursday 11/15/18 in class. Chemistry 391 - Problem Set #9 Due on Thursday 11/15/18 in class. Name 1. There is a real enzyme called cocaine esterase that is produced in bacteria that live at the base of the coca plant. The enzyme

More information

MOLECULAR MODELING IN BIOLOGY (BIO 3356) SYLLABUS

MOLECULAR MODELING IN BIOLOGY (BIO 3356) SYLLABUS New York City College of Technology School of Arts and Sciences Department of Biological Sciences MOLECULAR MODELING IN BIOLOGY (BIO 3356) SYLLABUS Course Information Course Title: Molecular Modeling in

More information

Protein Structure Analysis

Protein 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 information

Modelação e Simulação de Sistemas para Micro/Nano Tecnologias

Modelação e Simulação de Sistemas para Micro/Nano Tecnologias Modelação e Simulação de Sistemas para Micro/Nano Tecnologias http://gec.di.uminho.pt/mmnt/modsim/ Alberto José Proença, António Joaquim Esteves 2011/12 Mestrado em Micro/Nano Tecnologias ESCOLA DE ENGENHARIA

More information

Molecular Simulation II. Classical Mechanical Treatment

Molecular 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 information

Why Proteins Fold. How Proteins Fold? e - ΔG/kT. Protein Folding, Nonbonding Forces, and Free Energy

Why Proteins Fold. How Proteins Fold? e - ΔG/kT. Protein Folding, Nonbonding Forces, and Free Energy Why Proteins Fold Proteins are the action superheroes of the body. As enzymes, they make reactions go a million times faster. As versatile transport vehicles, they carry oxygen and antibodies to fight

More information

Supporting Online Material for

Supporting Online Material for www.sciencemag.org/cgi/content/full/309/5742/1868/dc1 Supporting Online Material for Toward High-Resolution de Novo Structure Prediction for Small Proteins Philip Bradley, Kira M. S. Misura, David Baker*

More information

The Molecular Dynamics Method

The 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 information

Modeling for 3D structure prediction

Modeling for 3D structure prediction Modeling for 3D structure prediction What is a predicted structure? A structure that is constructed using as the sole source of information data obtained from computer based data-mining. However, mixing

More information

Free-Energy Relationships for the Interactions of Tryptophan with Phosphocholines. Georg Blaser, a John M. Sanderson* a and Mark R.

Free-Energy Relationships for the Interactions of Tryptophan with Phosphocholines. Georg Blaser, a John M. Sanderson* a and Mark R. Free-Energy Relationships for the Interactions of Tryptophan with Phosphocholines Georg Blaser, a John M. Sanderson* a and Mark R. Wilson a S1 Supplementary Information Contents (click to jump to the corresponding

More information

Ranjit P. Bahadur Assistant Professor Department of Biotechnology Indian Institute of Technology Kharagpur, India. 1 st November, 2013

Ranjit P. Bahadur Assistant Professor Department of Biotechnology Indian Institute of Technology Kharagpur, India. 1 st November, 2013 Hydration of protein-rna recognition sites Ranjit P. Bahadur Assistant Professor Department of Biotechnology Indian Institute of Technology Kharagpur, India 1 st November, 2013 Central Dogma of life DNA

More information

k θ (θ θ 0 ) 2 angles r i j r i j

k θ (θ θ 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 information

Kd = koff/kon = [R][L]/[RL]

Kd = koff/kon = [R][L]/[RL] Taller de docking y cribado virtual: Uso de herramientas computacionales en el diseño de fármacos Docking program GLIDE El programa de docking GLIDE Sonsoles Martín-Santamaría Shrödinger is a scientific

More information

Protein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall How do we go from an unfolded polypeptide chain to a

Protein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall How do we go from an unfolded polypeptide chain to a Lecture 11: Protein Folding & Stability Margaret A. Daugherty Fall 2004 How do we go from an unfolded polypeptide chain to a compact folded protein? (Folding of thioredoxin, F. Richards) Structure - Function

More information

Principles of Physical Biochemistry

Principles of Physical Biochemistry Principles of Physical Biochemistry Kensal E. van Hold e W. Curtis Johnso n P. Shing Ho Preface x i PART 1 MACROMOLECULAR STRUCTURE AND DYNAMICS 1 1 Biological Macromolecules 2 1.1 General Principles

More information

Energy functions and their relationship to molecular conformation. CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror

Energy functions and their relationship to molecular conformation. CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror Energy functions and their relationship to molecular conformation CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror Yesterday s Nobel Prize: single-particle cryoelectron microscopy 2 Outline Energy

More information

Lecture 10 (10/4/17) Lecture 10 (10/4/17)

Lecture 10 (10/4/17) Lecture 10 (10/4/17) Lecture 10 (10/4/17) Reading: Ch4; 125, 138-141, 141-142 Problems: Ch4 (text); 7, 9, 11 Ch4 (study guide); 1, 2 NEXT Reading: Ch4; 125, 132-136 (structure determination) Ch4; 12-130 (Collagen) Problems:

More information

BSc and MSc Degree Examinations

BSc and MSc Degree Examinations Examination Candidate Number: Desk Number: BSc and MSc Degree Examinations 2018-9 Department : BIOLOGY Title of Exam: Molecular Biology and Biochemistry Part I Time Allowed: 1 hour and 30 minutes Marking

More information

L718Q mutant EGFR escapes covalent inhibition by stabilizing. a non-reactive conformation of the lung cancer drug. osimertinib

L718Q mutant EGFR escapes covalent inhibition by stabilizing. a non-reactive conformation of the lung cancer drug. osimertinib Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2018 Electronic Supplementary Information (ESI) for L718Q mutant EGFR escapes covalent inhibition

More information

Protein Structure & Motifs

Protein Structure & Motifs & Motifs Biochemistry 201 Molecular Biology January 12, 2000 Doug Brutlag Introduction Proteins are more flexible than nucleic acids in structure because of both the larger number of types of residues

More information

Biophysics II. Hydrophobic Bio-molecules. Key points to be covered. Molecular Interactions in Bio-molecular Structures - van der Waals Interaction

Biophysics II. Hydrophobic Bio-molecules. Key points to be covered. Molecular Interactions in Bio-molecular Structures - van der Waals Interaction Biophysics II Key points to be covered By A/Prof. Xiang Yang Liu Biophysics & Micro/nanostructures Lab Department of Physics, NUS 1. van der Waals Interaction 2. Hydrogen bond 3. Hydrophilic vs hydrophobic

More information

Lecture 5: Electrostatic Interactions & Screening

Lecture 5: Electrostatic Interactions & Screening Lecture 5: Electrostatic Interactions & Screening Lecturer: Prof. Brigita Urbanc (brigita@drexel.edu) PHYS 461 & 561, Fall 2009-2010 1 A charged particle (q=+1) in water, at the interface between water

More information

Advanced in silico drug design

Advanced 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 information

Protein Structure Analysis

Protein 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 information

= (-22) = +2kJ /mol

= (-22) = +2kJ /mol Lecture 8: Thermodynamics & Protein Stability Assigned reading in Campbell: Chapter 4.4-4.6 Key Terms: DG = -RT lnk eq = DH - TDS Transition Curve, Melting Curve, Tm DH calculation DS calculation van der

More information

ENERGY MINIMIZATION AND CONFORMATION SEARCH ANALYSIS OF TYPE-2 ANTI-DIABETES DRUGS

ENERGY MINIMIZATION AND CONFORMATION SEARCH ANALYSIS OF TYPE-2 ANTI-DIABETES DRUGS Int. J. Chem. Sci.: 6(2), 2008, 982-992 EERGY MIIMIZATI AD CFRMATI SEARC AALYSIS F TYPE-2 ATI-DIABETES DRUGS R. PRASAA LAKSMI a, C. ARASIMA KUMAR a, B. VASATA LAKSMI, K. AGA SUDA, K. MAJA, V. JAYA LAKSMI

More information

PROTEIN STRUCTURE AMINO ACIDS H R. Zwitterion (dipolar ion) CO 2 H. PEPTIDES Formal reactions showing formation of peptide bond by dehydration:

PROTEIN STRUCTURE AMINO ACIDS H R. Zwitterion (dipolar ion) CO 2 H. PEPTIDES Formal reactions showing formation of peptide bond by dehydration: PTEI STUTUE ydrolysis of proteins with aqueous acid or base yields a mixture of free amino acids. Each type of protein yields a characteristic mixture of the ~ 20 amino acids. AMI AIDS Zwitterion (dipolar

More information

Conformational Analysis of the Crystal Structure for MDI/ BDO Hard Segments of Polyurethane Elastomers

Conformational Analysis of the Crystal Structure for MDI/ BDO Hard Segments of Polyurethane Elastomers Conformational Analysis of the Crystal Structure for MDI/ BDO Hard Segments of Polyurethane Elastomers CHRIS W. PATTERSON, DAVID HANSON, ANTONIO REDONDO, STEPHEN L. SCOTT, NEIL HENSON Theoretical Division,

More information

Molecular Mechanics. C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology. January 2001

Molecular Mechanics. C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology. January 2001 Molecular Mechanics C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology January 2001 Introduction Molecular Mechanics uses classical type models to predict the energy

More information

OpenDiscovery: Automated Docking of Ligands to Proteins and Molecular Simulation

OpenDiscovery: Automated Docking of Ligands to Proteins and Molecular Simulation OpenDiscovery: Automated Docking of Ligands to Proteins and Molecular Simulation Gareth Price Computational MiniProject OpenDiscovery Aims + Achievements Produce a high-throughput protocol to screen a

More information

Protein Structures. 11/19/2002 Lecture 24 1

Protein Structures. 11/19/2002 Lecture 24 1 Protein Structures 11/19/2002 Lecture 24 1 All 3 figures are cartoons of an amino acid residue. 11/19/2002 Lecture 24 2 Peptide bonds in chains of residues 11/19/2002 Lecture 24 3 Angles φ and ψ in the

More information

Protein Bioinformatics Computer lab #1 Friday, April 11, 2008 Sean Prigge and Ingo Ruczinski

Protein Bioinformatics Computer lab #1 Friday, April 11, 2008 Sean Prigge and Ingo Ruczinski Protein Bioinformatics 260.655 Computer lab #1 Friday, April 11, 2008 Sean Prigge and Ingo Ruczinski Goals: Approx. Time [1] Use the Protein Data Bank PDB website. 10 minutes [2] Use the WebMol Viewer.

More information

Parallel Utility for Modeling of Molecular Aggregation

Parallel 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 information

Statistical Mechanics for Proteins

Statistical Mechanics for Proteins The Partition Function From Q all relevant thermodynamic properties can be obtained by differentiation of the free energy F: = kt q p E q pd d h T V Q ), ( exp 1! 1 ),, ( 3 3 3 ),, ( ln ),, ( T V Q kt

More information

Exercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions

Exercise 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 information

Multiscale Materials Modeling

Multiscale Materials Modeling Multiscale Materials Modeling Lecture 02 Capabilities of Classical Molecular Simulation These notes created by David Keffer, University of Tennessee, Knoxville, 2009. Outline Capabilities of Classical

More information

Electronic Supplementary Information for A molecular dynamics investigation on the cross-linking and physical properties of epoxy-based materials

Electronic Supplementary Information for A molecular dynamics investigation on the cross-linking and physical properties of epoxy-based materials Electronic Supplementary Material (ESI) for RSC Advances. This journal is The Royal Society of Chemistry 2014 Electronic Supplementary Information for A molecular dynamics investigation on the cross-linking

More information

Physiochemical Properties of Residues

Physiochemical Properties of Residues Physiochemical Properties of Residues Various Sources C N Cα R Slide 1 Conformational Propensities Conformational Propensity is the frequency in which a residue adopts a given conformation (in a polypeptide)

More information

1 of 43 Boardworks Ltd Chemistry 11. Chemical Bonding

1 of 43 Boardworks Ltd Chemistry 11. Chemical Bonding 1 of 43 Boardworks Ltd 2009 Chemistry 11 Chemical Bonding 2 of 43 Boardworks Ltd 2009 Electrostatic Forces An electrostatic force is a forces existing as a result of the attraction or repulsion between

More information