Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.

Size: px
Start display at page:

Download "Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar."

Transcription

1 Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002

2 Molecular Dynamics: Introduction At physiological conditions, the biomolecules undergo several movements and changes The time-scales of the motions are diverse, ranging from few femtoseconds to few seconds These motions are crucial for the function of the biomolecules

3 Molecular Dynamics: Introduction Newton s second law of motion

4 Molecular Dynamics: Introduction We need to know The motion of the atoms in a molecule, x(t) and therefore, the potential energy, V(x)

5 Molecular Dynamics: Introduction How do we describe the potential energy V(x) for a molecule? Potential Energy includes terms for Bond stretching Angle Bending Torsional rotation Improper dihedrals

6 Molecular Dynamics: Introduction Potential energy includes terms for (contd.) Electrostatic Interactions van der Waals Interactions

7 Molecular Dynamics: Introduction Equation for covalent terms in P.E. ( ) )] cos( [1 ) ( ) ( ) ( φ φ ω ω θ θ ω θ = n A k k l l k R V torsions n impropers angles bonds l bonded

8 Molecular Dynamics: Introduction Equation for non-bonded terms in P.E. V nonbonded ( R) min min rij r 12 ij 6 = ( ε [( ) 2( ) ] + ij r r i< j ij ij q q i 4πε r j ε r 0 ij

9 Molecular Dynamics: Introduction Each of these interactions exerts a force onto a given atom of the molecule The total resulting force on each atom is calculated using the PE function Knowing the force on an atom, its movement due to the force is then calculated:

10 Molecular Dynamics: Introduction To do this, we should know at given time t, initial position of the atom x 1 its velocity v 1 = dx 1 /dt and the acceleration a 1 = d 2 x 1 /dt 2 = m -1 F(x 1 )

11 Molecular Dynamics: Introduction The position x 2, of the atom after time interval t would be, x 2 = x 1 + v 1 t and the velocity v 2 would be, v = v + a t = v + 1 m F x t = v ( m ) 1 1 dv dx x 1 t

12 How a molecule changes during MD

13 Molecular Dynamics: Introduction In general, given the values x 1, v 1 and the potential energy V(x), the molecular trajectory x(t) can be calculated, using, t dx x dv m v v t v x x x i i i i i i = + = 1 ) (

14 Generalizing these ideas, the trajectories for all the atoms of a molecule can be calculated.

15 The Necessary Ingredients Description of the structure: atoms and connectivity Initial structure: geometry of the system Potential Energy Function: force field AMBER CVFF CFF95 Universal

16 Protein-specific Applications of MD Calculation of thermodynamic properties such as internal energy, free energy Studying the protein folding / unfolding process Studying conformational properties and transitions due to environmental conditions Studying conformational distributions in molecular system.

17 An overview of various motions in proteins (1) Motion Relative vibration of bonded atoms Elastic vibration of globular region Rotation of side chains at surface Spatial extent (nm) 0.2 to to to 1 Log 10 of characteristic time (s) -14 to to to 10

18 An overview of various motions in proteins (2) Motion Relative motion of different globular regions (hinge bending) Rotation of medium-sized side chains in interior Protein Folding Spatial Extent (nm) 1 to 2 0.5??? Log 10 of characteristic time (s) -11 to 7-4 to 0-5 to 2

19 A typical MD simulation protocol Initial random structure generation Initial energy minimization Equilibration Dynamics run with capture of conformations at regular intervals Energy minimization of each captured conformation

20 Essential Parameters for MD (to be set by user) Temperature Pressure Time step Dielectric constant Force field Durations of equilibration and MD run ph effect (addition of ions)

21 WHAT IS AMBER? AMBER (Assisted Model Building with Energy Refinement). Allows users to carry out molecular dynamics simulations Updated forcefield for proteins and nucleic acids Parallelized dynamics codes Ewald sum periodicity New graphical and text-based tools for building molecules Powerful tools for NMR spectral simulations New dynamics and free energy program

22 CASE STUDY Type II restriction endonucleases recognize DNA sequences of 4 to 8 base pairs in length and require Mg 2+ to hydrolyse DNA. The recognition of DNA sequences by endonucleases is still an open question. PvuII endonuclease, recognizes the sequence 5 - CAGCTG-3 and cleaves between the central G and C bases in both strands. Though crystal structure of the PvuII-DNA complex have been reported, very little is known about the steps involved in the recognition of the cleavage site by the PvuII enzyme. Molecular dynamics (MD) simulation is a powerful computational approach to study the macromolecular structure and motions.

23 CASE STUDY: METHODS (MD Simulations) Simulations were carried out on the sequence 5 -TGACCAGCTGGTC-3 Rectangular box (60 X 48 X 54 Å 3 ) containing 24 Na +, using PBC SHAKE algorithm Integration time step of 1 fs 283 K with Berendsen coupling Particle Mesh Ewald (PME) method 9.0 Å cutoff was applied to the Lennard-Jones interaction term. Equilibration was performed by slowly raising the temperature from 100 to 283 K. Production run was initiated for ns and the structures were saved at intervals of one picosecond. The trajectory files were imaged using the RDPARM program and viewed and analysed using the MOIL-VIEW and CURVES packages respectively.

24 STARTING DNA MODEL

25 DNA MODEL WITH IONS

26 DNA in a box of water

27 SNAPSHOTS

28 SNAPSHOTS

29 DOCKING

30 DOCKING

31 Computing Potentials Boundary Conditions Truncating Potentials Approximating Potentials

32 Computing Potentials Fast decaying potentials can be truncated (for example, Lennart-Jones Potentials) Other potentials such as Coloumbic potentials must be approximated.

33 Fast Multipole Method What is the gravitational potential where we are sitting because of the sun? Can we treat the sun as a point mass or should we integrate over the volume of the sun. Fast Multipole Method works on the principle of hierarchical aggregation.

34 Fast Multipole Method Collect spatially proximate objects and represent it with a series -- a multipole series. For an evaluation point far from this group, compute the series as opposed to point-wise evaluation. Do this hierarchically. Reduces complexity from O(n 2 ) to O(n).

35 Ewald Summations Decompose the potential into two terms: an infinite term that whose FFT can be computed. a fast converging local term that can be truncated. Reduces complexity to O(n log n)

36 Software CHARMM NAMD many other in-grown codes.

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

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

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

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 Molecular Dynamics Method

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

Lecture 11: Potential Energy Functions

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

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

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

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

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

Introduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA

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

Non-bonded interactions

Non-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 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

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

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 Allowed phi-psi angles Red areas are preferred, yellow areas are allowed, and white is avoided 2.3a Hamiltonian

More information

Protein Dynamics. The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron.

Protein Dynamics. The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron. Protein Dynamics The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron. Below is myoglobin hydrated with 350 water molecules. Only a small

More information

Non-bonded interactions

Non-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 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

Molecular dynamics simulation of Aquaporin-1. 4 nm

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

Molecular Dynamics. A very brief introduction

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

Advanced Molecular Molecular Dynamics

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

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

All-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010

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

Computational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data

Computational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data Computational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data ABSTRACT Keyword Lei Shi 1 Advisor: Gianluigi Veglia 1,2 Department of Chemistry 1, & Biochemistry, Molecular

More information

Modeling Biological Systems Opportunities for Computer Scientists

Modeling Biological Systems Opportunities for Computer Scientists Modeling Biological Systems Opportunities for Computer Scientists Filip Jagodzinski RBO Tutorial Series 25 June 2007 Computer Science Robotics & Biology Laboratory Protein: πρώτα, "prota, of Primary Importance

More information

3rd 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. 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 information

Applications of Molecular Dynamics

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

Computation of non-bonded interactions: Part 1

Computation of non-bonded interactions: Part 1 Computation of non-bonded interactions: Part 1 Computation of non-local (non-bonded) interactions in the potential function E p non-local scales with the number of atoms as N. For the large molecular systems

More information

Gromacs Workshop Spring CSC

Gromacs 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 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

An FPGA Implementation of Reciprocal Sums for SPME

An FPGA Implementation of Reciprocal Sums for SPME An FPGA Implementation of Reciprocal Sums for SPME Sam Lee and Paul Chow Edward S. Rogers Sr. Department of Electrical and Computer Engineering University of Toronto Objectives Accelerate part of Molecular

More information

Computational Molecular Biophysics. Computational Biophysics, GRS Jülich SS 2013

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

Introduction to molecular dynamics

Introduction 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 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

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

Molecular Dynamics Simulation of a Nanoconfined Water Film

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

Using Molecular Dynamics to Compute Properties CHEM 430

Using Molecular Dynamics to Compute Properties CHEM 430 Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system

More information

Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers

Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers Jaewoon Jung (RIKEN, RIKEN AICS) Yuji Sugita (RIKEN, RIKEN AICS, RIKEN QBiC, RIKEN ithes) Molecular Dynamics

More 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

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

Subject of the Lecture:

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

Contents. xiii. Preface v

Contents. xiii. Preface v Contents Preface Chapter 1 Biological Macromolecules 1.1 General PrincipIes 1.1.1 Macrornolecules 1.2 1.1.2 Configuration and Conformation Molecular lnteractions in Macromolecular Structures 1.2.1 Weak

More information

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

Principles and Applications of Molecular Dynamics Simulations with NAMD

Principles 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 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

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

Scuola di Chimica Computazionale

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

Can 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? 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 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

A Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane

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

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid

Force 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 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

Molecular Dynamics Studies of Human β-glucuronidase

Molecular 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 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 Outline Energy functions for proteins (or biomolecular systems more generally)

More information

Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland

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

Protein-mineral interactions: molecular dynamics simulations capture importance of variations in mineral surface composition and structure

Protein-mineral interactions: molecular dynamics simulations capture importance of variations in mineral surface composition and structure Supplemental Information for: Protein-mineral interactions: molecular dynamics simulations capture importance of variations in mineral surface composition and structure Amity Andersen, *,1 Patrick N. Reardon,

More information

Computational Modeling of Protein Dynamics. Yinghao Wu Department of Systems and Computational Biology Albert Einstein College of Medicine Fall 2014

Computational Modeling of Protein Dynamics. Yinghao Wu Department of Systems and Computational Biology Albert Einstein College of Medicine Fall 2014 Computational Modeling of Protein Dynamics Yinghao Wu Department of Systems and Computational Biology Albert Einstein College of Medicine Fall 2014 Outline Background of protein dynamics Basic computational

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

Lecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects

Lecture 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 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

Supporting Information: Improved Parametrization. of Lithium, Sodium, Potassium, and Magnesium ions. for All-Atom Molecular Dynamics Simulations of

Supporting Information: Improved Parametrization. of Lithium, Sodium, Potassium, and Magnesium ions. for All-Atom Molecular Dynamics Simulations of Supporting Information: Improved Parametrization of Lithium, Sodium, Potassium, and Magnesium ions for All-Atom Molecular Dynamics Simulations of Nucleic Acid Systems. Jejoong Yoo and Aleksei Aksimentiev,,

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

Atomistic Modeling of Small-Angle Scattering Data Using SASSIE-web

Atomistic Modeling of Small-Angle Scattering Data Using SASSIE-web Course Introduction Atomistic Modeling of Small-Angle Scattering Data Using SASSIE-web September 21-23, 2016 Advanced Photon Source Argonne National Laboratory, Argonne, IL ccpsas.org Scatters & Simulators

More information

Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials

Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials Dr. V. Subramanian Chemical Laboratory, IPC Division CSIR-Central Leather Research

More 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

Biomolecular modeling I

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

Biomolecular modeling I

Biomolecular modeling I 2016, December 6 Biomolecular structure Structural elements of life Biomolecules proteins, nucleic acids, lipids, carbohydrates... Biomolecular structure Biomolecules biomolecular complexes aggregates...

More information

Introduction to MD simulations of DNA systems. Aleksei Aksimentiev

Introduction to MD simulations of DNA systems. Aleksei Aksimentiev Introduction to MD simulations of DNA systems Aleksei Aksimentiev Biological Modeling at Different Scales spanning orders of magnitude in space and time The computational microscope Massive parallel computer

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

Computer simulation methods (2) Dr. Vania Calandrini

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

An Informal AMBER Small Molecule Force Field :

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

Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia

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

MARTINI simulation details

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

Exploring the Sequence Dependent Structure and Dynamics of DNA with Molecular Dynamics Simulation

Exploring the Sequence Dependent Structure and Dynamics of DNA with Molecular Dynamics Simulation Exploring the Sequence Dependent Structure and Dynamics of DNA with Molecular Dynamics Simulation Sarah Harris School of Physics and Astronomy University of Leeds Introduction Calculations of the charge

More information

PH 548 Atomistic Simulation Techniques

PH 548 Atomistic Simulation Techniques PH 548 Atomistic Simulation Techniques Lectures: Lab: 2-0-2-6 Tuesday 12-1 PM Thursday 12-1 PM Monday 2-5 PM P. K. Padmanabhan PH548: Two Excellent Books How to do it? 1 2 F 12 F 23 3 F i = F ij j i F

More information

Molecular Mechanics / ReaxFF

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

The Molecular Dynamics Method

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

Analysis of the simulation

Analysis of the simulation Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities

More 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

Towards accurate calculations of Zn 2+ binding free energies in zinc finger proteins

Towards accurate calculations of Zn 2+ binding free energies in zinc finger proteins Towards accurate calculations of Zn 2+ binding free energies in zinc finger proteins UNDERGRADUATE HONORS RESEARCH THESIS Presented in Partial Fulfillment of the Requirements for the Bachelor of Science

More information

Course Introduction. SASSIE CCP-SAS Workshop. January 23-25, ISIS Neutron and Muon Source Rutherford Appleton Laboratory UK

Course Introduction. SASSIE CCP-SAS Workshop. January 23-25, ISIS Neutron and Muon Source Rutherford Appleton Laboratory UK Course Introduction SASSIE CCP-SAS Workshop January 23-25, 2017 ISIS Neutron and Muon Source Rutherford Appleton Laboratory UK ccpsas.org Scatters & Simulators Develop a community of users and developers

More information

Lecture 3: methods and terminology, part I

Lecture 3: methods and terminology, part I Eritis sicut Deus, scientes bonum et malum. Lecture 3: methods and terminology, part I Molecular dynamics, internal coordinates, Z-matrices, force fields, periodic boundary conditions, Hartree-Fock theory,

More information

The Molecular Dynamics Simulation Process

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

Force fields, thermo- and barostats. Berk Hess

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

Structural and mechanistic insight into the substrate. binding from the conformational dynamics in apo. and substrate-bound DapE enzyme

Structural and mechanistic insight into the substrate. binding from the conformational dynamics in apo. and substrate-bound DapE enzyme Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 215 Structural and mechanistic insight into the substrate binding from the conformational

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

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

Molecular dynamics simulations of a single stranded (ss) DNA

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

Supporting Material for. Microscopic origin of gating current fluctuations in a potassium channel voltage sensor

Supporting Material for. Microscopic origin of gating current fluctuations in a potassium channel voltage sensor Supporting Material for Microscopic origin of gating current fluctuations in a potassium channel voltage sensor J. Alfredo Freites, * Eric V. Schow, * Stephen H. White, and Douglas J. Tobias * * Department

More information

Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization

Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization Richard H. Gee * Naida Lacevic and Laurence E. Fried University of California Lawrence

More information

Supporting Information

Supporting Information Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 207 Supporting Information Carbon Nanoscroll- Silk Crystallite Hybrid Structure with Controllable Hydration

More information

Routine access to millisecond timescale events with accelerated molecular dynamics

Routine access to millisecond timescale events with accelerated molecular dynamics Routine access to millisecond timescale events with accelerated molecular dynamics Levi C.T. Pierce, Romelia Salomon-Ferrer, Cesar Augusto F. de Oliveira #, J. Andrew McCammon #, Ross C. Walker * SUPPORTING

More information

MECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS. Vinu Unnithan and J. N. Reddy

MECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS. Vinu Unnithan and J. N. Reddy MECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS Vinu Unnithan and J. N. Reddy US-South American Workshop: Mechanics and Advanced Materials Research and Education

More information

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

Atomistic Modeling of DNA and Protein Structures

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

Supporting Information. for An atomistic modeling of F-actin mechanical responses and determination of. mechanical properties

Supporting Information. for An atomistic modeling of F-actin mechanical responses and determination of. mechanical properties Supporting Information for An atomistic modeling of F-actin mechanical responses and determination of mechanical properties Authors: Si Li 1, Jin Zhang 2, Chengyuan Wang 1, Perumal Nithiarasu 1 1. Zienkiewicz

More information