The Force Field Toolkit (fftk)
|
|
- Grace Hood
- 5 years ago
- Views:
Transcription
1 Parameterizing Small Molecules Using: The Force Field Toolkit (fftk) Christopher G. Mayne, Emad Tajkhorshid Beckman Institute for Advanced Science and Technology University of Illinois, Urbana-Champaign Klaus Schulten University of Illinois, Urbana-Champaign James C. Gumbart, Anna Pavlova Georgia Institute of Technology NIH Hands n Workshop Atlanta November 5th, 2014
2 MD Simulations of Biological Systems Molecular Mechanics Force Fields U = Ubonds + Uangles + Udihedrals + UvdW + Ucoulombic bonded non-bonded The CHARMM Force Field U = Σ k bond - ) i (ri r0 2 + Σ ( θ - ) ki angle i 2 + θ0 bonds angles Σ [1 + cos( + )] ki dihedral niɸi δi dihedrals + Σ i Σ j i 4 ij( σ ij rij ) 12 ( σij ) 6 rij - + ΣΣ qiqj i j i rij
3 MD Simulations of Biological Systems Molecular Mechanics Force Fields U = Ubonds + Uangles + Udihedrals + UvdW + Ucoulombic bonded non-bonded The CHARMM Force Field U = Σ k bond - ) i (ri r0 2 + Σ ( θ - ) ki angle i 2 + θ0 bonds angles Σ [1 + cos( + )] ki dihedral niɸi δi dihedrals + Σ i Σ j i 4 ij( σ ij rij ) 12 ( σij ) 6 rij - + ΣΣ qiqj i j i rij
4 Parameter Transferability In Biopolymers Parameter set describes molecular behavior in varied chemical (connectivity) and spatial (conformation) contexts Peptides and Proteins Nucleic Acids H N Me S Me H H N N N N N H H H Me Me NH 2 H H 2 N N N N N P N NH 2 N P HN N Me P N N N NH NH 2 Key Features: R limited set of isolated building blocks repetitive backbone unit
5 Parameterization as an Impasse H N H non-standard or engineered amino acids CH H 2 N CH HN small molecule ligands N HN N N Me Me N Me S H Me N S N Imatinib (Gleevec) Tiotropium (Spiriva) cofactors metal centers HS N H N H Me H Me P P N N NH 2 N N H 2 N S Fe S NH 2 Coenzyme A H P S NH 2 S H 2 N
6 General Parameterization Workflow PSF/PDB System Preparation Charges Bonds & Angles PAR File Dihedrals / Torsions
7 CGenFF Parameterization Workflow PSF/PDB Find Missing Parameters Geometry ptimization (QM) Water Interaction En. (QM) Charge ptimization Hessian Calculation (QM) Bond & Angle ptimization PAR File Torsion Scan (QM) Torsion ptimization K. Vanommeslaeghe et al., J. Comput. Chem. 2010, 31,
8 CGenFF Parameterization Workflow PSF/PDB PAR File Find Missing Parameters Geometry ptimization (QM) Water Interaction En. (QM) Charge ptimization Hessian Calculation (QM) Bond & Angle ptimization Torsion Scan (QM) Torsion ptimization build init PAR update PDB update PSF update PAR update PAR Calculation Action K. Vanommeslaeghe et al., J. Comput. Chem. 2010, 31,
9 CGenFF Parameterization Workflow PSF/PDB PAR File Find Missing Parameters Geometry ptimization (QM) Water Interaction En. (QM) Charge ptimization Hessian Calculation (QM) Bond & Angle ptimization Torsion Scan (QM) Torsion ptimization build init PAR update PDB update PSF update PAR update PAR Calculation Action C.G. Mayne et al., J. Comput. Chem. 2013, 34,
10 fftk Interface
11 fftk Interface
12 fftk Interface tasks organized under tabs
13 fftk Interface standard file dialogs action buttons action menus
14 Functionality Provided by fftk Setup & Perform Multi-dimensional ptimizations Core Functions Abstraction of Gaussian I/ (QM) Assess Performance of Parameters by Visualizing ptimization Data Support Functions Auto-detect Water Interaction Sites Auto-detect Charge Groups Auto-detect Non-redundant Torsions Build & Update Parameter Files Browse Existing Parameter Sets Write Updated Charges to PSF Reset pt. Input from utput Visualize Target Data in VMD Create Graphic bjects in VMD Label Atoms in VMD Read Input Parameters from File Read/Write Data From pt. Logs Export Plot Data to File Monitor ptimization Progress
15 fftk Exemplified by Charge ptimization PSF/PDB PAR File Find Missing Parameters Geometry ptimization (QM) Water Interaction En. (QM) Charge ptimization Hessian Calculation (QM) Bond & Angle ptimization Torsion Scan (QM) Torsion ptimization build init PAR update PDB update PSF update PAR update PAR Calculation Action
16 fftk Exemplified by Charge ptimization PSF/PDB PAR File Find Missing Parameters Geometry ptimization (QM) Water Interaction En. (QM) Charge ptimization Hessian Calculation (QM) Bond & Angle ptimization Torsion Scan (QM) Torsion ptimization build init PAR update PDB update PSF update PAR update PAR Calculation Action
17 Generating Charge ptimization Target Data Load QM optimized geometry Auto-detect interaction sites Genera pyrrolidine fftk GUI VMD main window
18 Generating Charge ptimization Target Data geometry Auto-detect interaction sites Generate Gaussian Input Files Run Q fftk GUI VMD main window
19 Generating Charge ptimization Target Data geometry Auto-detect interaction sites Generate Gaussian Input Files Run Q fftk GUI Donor VMD main window Acceptor
20 Generating Charge ptimization Target Data n sites Generate Gaussian Input Files Run QM Inspect water optimization Compute water position ptimize distance & rotation fftk GUI
21 Generating Charge ptimization Target Data Run QM Inspect water optimization Visually assess QM-optimized water position(s) fftk GUI
22 Generating Charge ptimization Target Data Run QM Inspect water optimization Visually assess QM-optimized water position(s) fftk GUI
23 Charge ptimization Setup ptimization Load QM Target Data Prepare ptimization ptimizer: Assign Charges Compute UMM, dmm, µmm Compute bjective Function Return ptimized Charges Analyze Performance bjective Function Σ wat. int. Σ wat. int. f(u MM -U QM ) + f(d MM -d QM ) + f(μ MM -μ QM ) Write Charges to PSF
24 Assessing MM Water-Interaction Profiles 10 8 Initial Charges Intermediate Charges Final ptimized Charges Literature Charges U MM-QM (kcal/mol) d MM-QM (Å)
25 Sampling MM Water-Interaction Profiles Mode: Simulated Annealing 6 5 H N H N H H N 4000 H 3 H 2 H ΔdMM-QM (Å) Iteration ΔUMM-QM (kcal/mol) 4 H N 4500
26 Plotting Charge ptimization Data
27 Plotting Charge ptimization Data
28 Plotting Charge ptimization Data
29 Plotting Charge ptimization Data
30 Plotting Charge ptimization Data export data ΔE from Target (kcal/mol) ptimization Iteration
31 Parameterization of Erythromycin macrolides - antibiotics that prevent peptide synthesis in bacterial ribosomes small subunit PTC PET large subunit nascent peptide macrolide Divided into four parts
32 Parameterization of Erythromycin macrolides - antibiotics that prevent peptide synthesis in bacterial ribosomes small subunit PTC PET large subunit nascent peptide macrolide Divided into four parts
33 Parameterization of Erythromycin macrolides - antibiotics that prevent peptide synthesis in bacterial ribosomes small subunit PTC PET large subunit nascent peptide macrolide Divided into four parts
34 Parameterization of Erythromycin macrolides - antibiotics that prevent peptide synthesis in bacterial ribosomes small subunit PTC PET large subunit nascent peptide macrolide Divided into four parts
35 Parameterization of Erythromycin macrolides - antibiotics that prevent peptide synthesis in bacterial ribosomes small subunit PTC PET large subunit nascent peptide macrolide Divided into four parts
36 Parameterization of Erythromycin macrolides - antibiotics that prevent peptide synthesis in bacterial ribosomes small subunit PTC PET large subunit nascent peptide macrolide Divided into four parts
37 Parameterization of Erythromycin macrolides - antibiotics that prevent peptide synthesis in bacterial ribosomes small subunit PTC PET large subunit nascent peptide macrolide Divided into four parts
38 Fitting of Partial Charges CHARMM partial charges are derived from interactions with water! Calculated for each atom that is fitted
39 Fitting of Partial Charges CHARMM partial charges are derived from interactions with water! Calculated for each atom that is fitted
40 Fitting of Partial Charges CHARMM partial charges are derived from interactions with water! Calculated for each atom that is fitted Loaded into FFTK
41 Accuracy of Fitting: Charges Partial charges for polar atoms Desosamine charges fitted for red atoms Charge CGenFF FFTK H1 H' FFTK can be used to get better agreement with QM data
42 Accuracy of Fitting: Charges Interaction Partial energies charges with for water polar for atoms polar atoms Desosamine charges fitted for red atoms E (kcal/mol) Charge CGenFF FFTK CGenFF FFTK QME H1 H' H1 H' FFTK can be used to get better agreement with QM data
43 Fitting of Partial Charges! ptimization of water interactions can be difficult for sterically hindered groups! ur self consistent iterative charge optimization scheme! Desosamine) Simpli&ied)Desosamine) H H 3 C N CH 3 CH 3 H CH 3 H H 3 C N H CH 3 H
44 Fitting of Partial Charges! ptimization of water interactions can be difficult for sterically hindered groups! ur self consistent iterative charge optimization scheme! Desosamine) Simpli&ied)Desosamine) H H 3 C N CH 3 CH 3 H CH 3 H H 3 C N H CH 3 H
45 Fitting of Partial Charges! ptimization of water interactions can be difficult for sterically hindered groups! ur self consistent iterative charge optimization scheme! Desosamine) Simpli&ied)Desosamine) H H 3 C N CH 3 CH 3 H CH 3 H H 3 C N H CH 3 H
46 An Alternative to Water Interactions for Charges RESP fitting of charges! Based electrostatic potential!!! Much easier to implement! Tests on compatibility with CHARMM in progress
47 An Alternative to Water Interactions for Charges RESP fitting of charges! Based electrostatic potential!!! Much easier to implement! Tests on compatibility with CHARMM in progress
48 Fitting of Bonds and Angles
49 Fitting of Bonds and Angles input Hessian file
50 Fitting of Bonds and Angles input Hessian file btained! force constants! and values
51 Fitting of Dihedrals K (1 + cos(n ))
52 Fitting of Dihedrals Feed in the dihedral scans K (1 + cos(n ))
53 Fitting of Dihedrals Select multiplicity and phase settings K (1 + cos(n ))
54 Fitting of Dihedrals Select multiplicity and phase settings multiplicity phase multiplicity phase K (1 + cos(n ))
55 Fitting of Dihedrals First fit is not great but better than the starting fit QM PES: black! First fit: blue! Starting fit: red!
56 Fitting of Dihedrals After 7 iterations the fit is much better First fit is not great but better than the starting fit QM PES: black! First fit: blue! Starting fit: red!
57 Two Approaches to Fitting the Dihedrals Several multiplicities! and free phase ne multiplicity and! locked phase Pro: very good fit of QM PES! Cons: possible incorrect behavior! multiple sets of force constants Pros: limits incorrect behavior,! sets of force constants! Con: fit to QM PES not always! possible In practice a tradeoff is needed!
58 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin
59 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin E (kcal/mol) multiplicities, free phase Dihedral Fit QM PES MM Fit Configuration
60 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin
61 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin
62 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin Using too many dihedral multiplicities can leads to distortion of a planar molecule!
63 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin
64 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin planar dihedrals have multiplicity 2 and phase 180
65 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin planar dihedrals have multiplicity 2 and phase 180 E (kcal/mol) Dihedral Fitting QM PES MM Fit Configuration The fit looks even better!
66 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin planar dihedrals have multiplicity 2 and phase 180
67 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin planar dihedrals have multiplicity 2 and phase 180 Planarity is maintained!
68 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin Problems persist eclipsed conformation of the alkane part
69 Example of verfitting imidazole-pyrridine moiety of antibiotic telithromycin Problems persist eclipsed conformation of the alkane part Restraining phase of CH dihedrals to 0 prevents eclipsed confirmations
70 How to know what terms to include?
71 fftk: Find Missing Parameters Geometry ptimization (QM) Water Interaction En. (QM) Charge ptimization Hessian Calculation (QM) Bond & Angle ptimization Torsion Scan (QM) Torsion ptimization Conclusions - Simplifies the parameterization workflow - ffers opportunity for extensive customization - Provides analytical tools to assess parameter performance
72 Future Directions Investigate forcefield performance Support additional force fields (e.g. AMBER) Improve optimization speed Expand analysis tools Support multiple QM packages
73 fftk Mayne et al.; J. Comp. Chem. 2013, 34, pp fftk is available as a VMD Plugin (1.9.1 or newer) Journal of Volume 34 Issues Included in this print edition: Issue 31 (December 5, 2013) Issue 32 (December 15, 2013) CMPUTATINAL rganic CHEMISTRY Inorganic Physical Biological Materials Research in Systems Neuroscience Editors: Charles L. Brooks III Masahiro Ehara Gernot Frenking Peter R. Schreiner Questions?
Force 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 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 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 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 informationPrinciples and Applications of Molecular Dynamics Simulations with NAMD
Principles and Applications of Molecular Dynamics Simulations with NAMD Nov. 14, 2016 Computational Microscope NCSA supercomputer JC Gumbart Assistant Professor of Physics Georgia Institute of Technology
More informationSupplementary Figure 1 Irregular arrangement of E,E-8-mer on TMA. STM height images formed when
Supplementary Figure 1 Irregular arrangement of E,E-8-mer on TMA. STM height images formed when a 5 µl heptanoic acid solution of E,E-8-mer is applied on: (a) a TMA templated HOPG substrate, U s = +1.00
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 informationSUPPLEMENTAL MATERIAL
SUPPLEMENTAL MATERIAL Systematic Coarse-Grained Modeling of Complexation between Small Interfering RNA and Polycations Zonghui Wei 1 and Erik Luijten 1,2,3,4,a) 1 Graduate Program in Applied Physics, Northwestern
More informationAssignment A02: Geometry Definition: File Formats, Redundant Coordinates, PES Scans
Assignment A02: Geometry Definition: File Formats, Redundant Coordinates, PES Scans In Assignments A00 and A01, you familiarized yourself with GaussView and G09W, you learned the basics about input (GJF)
More informationVMD: Visual Molecular Dynamics. Our Microscope is Made of...
VMD: Visual Molecular Dynamics Computational Microscope / Tool to Think amino acid tyrosine enzymatic control BPTI VMD tutorial traficking Ubiquitin case study http://www.ks.uiuc.edu/training/casestudies/
More informationPreparing a PDB File
Figure 1: Schematic view of the ligand-binding domain from the vitamin D receptor (PDB file 1IE9). The crystallographic waters are shown as small spheres and the bound ligand is shown as a CPK model. HO
More informationComputational Structural Biology and Molecular Simulation. Introduction to VMD Molecular Visualization and Analysis
Computational Structural Biology and Molecular Simulation Introduction to VMD Molecular Visualization and Analysis Emad Tajkhorshid Department of Biochemistry, Beckman Institute, Center for Computational
More informationExamples of Protein Modeling. Protein Modeling. Primary Structure. Protein Structure Description. Protein Sequence Sources. Importing Sequences to MOE
Examples of Protein Modeling Protein Modeling Visualization Examination of an experimental structure to gain insight about a research question Dynamics To examine the dynamics of protein structures To
More informationRevised CHARMM Force Field Parameters for Iron- Containing Cofactors of Photosystem II
WWW.C-CHEM.ORG FULL PAPER Revised CHARMM Force Field Parameters for Iron- Containing Cofactors of Photosystem II Suliman Adam, [a] Michaela Knapp-Mohammady, [b] Jun Yi, [c] and Ana-Nicoleta Bondar [a]
More informationAPBS electrostatics in VMD - Software. APBS! >!Examples! >!Visualization! >! Contents
Software Search this site Home Announcements An update on mailing lists APBS 1.2.0 released APBS 1.2.1 released APBS 1.3 released New APBS 1.3 Windows Installer PDB2PQR 1.7.1 released PDB2PQR 1.8 released
More informationComputational 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 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 informationSupporting information for: Redox Potential Tuning Through Differential. Quinone Binding in the Photosynthetic Reaction
Supporting information for: Redox Potential Tuning Through Differential Quinone Binding in the Photosynthetic Reaction Center of Rhodobacter sphaeroides Josh V. Vermaas,,, Alexander T. Taguchi,,,, Sergei
More informationWhy study protein dynamics?
Why study protein dynamics? Protein flexibility is crucial for function. One average structure is not enough. Proteins constantly sample configurational space. Transport - binding and moving molecules
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationSeeSAR 7.1 Beginners Guide. June 2017
SeeSAR 7.1 Beginners Guide June 2017 Part 1: Basics 1 Type a pdb code and press return or Load your own protein or already existing project, or Just load molecules To begin, let s type 2zff and download
More informationRetinal Proteins (Rhodopsins) Vision, Bioenergetics, Phototaxis. Bacteriorhodopsin s Photocycle. Bacteriorhodopsin s Photocycle
Molecular chanisms of Photoactivation and Spectral Tuning in Retinal Proteins Emad Tajkhorshid Theoretical and Computational Biophysics Group Beckman Institute University of Illinois at Urbana-Champaign
More informationThe Molecular Dynamics Simulation Process
The Molecular Dynamics Simulation Process For textbooks see: M.P. Allen and D.J. Tildesley. Computer Simulation of Liquids.Oxford University Press, New York, 1987. D. Frenkel and B. Smit. Understanding
More 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 informationprotein synthesis and the ribosome
protein synthesis and the ribosome Central dogma of biology DNA codes for DNA DNA codes for RNA RNA codes for proteins not surprisingly, many points for regulation of the process RNA codes for proteins
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 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 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 information1 Introduction to Computational Chemistry (Spartan)
1 Introduction to Computational Chemistry (Spartan) Start Spartan by clicking Start / Programs / Spartan Then click File / New Exercise 1 Study of H-X-H Bond Angles (Suitable for general chemistry) Structure
More informationGeneral Chemistry Lab Molecular Modeling
PURPOSE The objectives of this experiment are PROCEDURE General Chemistry Lab Molecular Modeling To learn how to use molecular modeling software, a commonly used tool in chemical research and industry.
More informationEvolution of Translation: Dynamics of Recognition in RNA:Protein Complexes
Evolution of Translation: Dynamics of Recognition in RNA:Protein Complexes Zaida (Zan) Luthey-Schulten Dept. Chemistry, Beckman Institute, Biophysics, Institute of Genomics Biology, & Physics NIH Workshop
More information3. An Introduction to Molecular Mechanics
3. An Introduction to Molecular Mechanics Introduction When you use Chem3D to draw molecules, the program assigns bond lengths and bond angles based on experimental data. The program does not contain real
More informationMolecular 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 informationIdentifying Interaction Hot Spots with SuperStar
Identifying Interaction Hot Spots with SuperStar Version 1.0 November 2017 Table of Contents Identifying Interaction Hot Spots with SuperStar... 2 Case Study... 3 Introduction... 3 Generate SuperStar Maps
More informationTutorial: Structural Analysis of a Protein-Protein Complex
Molecular Modeling Section (MMS) Department of Pharmaceutical and Pharmacological Sciences University of Padova Via Marzolo 5-35131 Padova (IT) @contact: stefano.moro@unipd.it Tutorial: Structural Analysis
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 informationMolecular modeling with InsightII
Molecular modeling with InsightII Yuk Sham Computational Biology/Biochemistry Consultant Phone: (612) 624 7427 (Walter Library) Phone: (612) 624 0783 (VWL) Email: shamy@msi.umn.edu How to run InsightII
More informationEnhanced Sampling and Free Energy Applications in Biomolecular Modeling
Enhanced Sampling and Free Energy Applications in Biomolecular Modeling Emad Tajkhorshid NIH Center for Macromolecular Modeling and Bioinformatics Beckman Institute for Advanced Science and Technology
More informationICM-Chemist-Pro How-To Guide. Version 3.6-1h Last Updated 12/29/2009
ICM-Chemist-Pro How-To Guide Version 3.6-1h Last Updated 12/29/2009 ICM-Chemist-Pro ICM 3D LIGAND EDITOR: SETUP 1. Read in a ligand molecule or PDB file. How to setup the ligand in the ICM 3D Ligand Editor.
More informationMolecular Dynamics Flexible Fitting
Molecular Dynamics Flexible Fitting Ryan McGreevy Research Programmer University of Illinois at Urbana-Champaign NIH Resource for Macromolecular Modeling and Bioinformatics Molecular Dynamics Flexible
More informationHands-on Course in Computational Structural Biology and Molecular Simulation BIOP590C/MCB590C. Course Details
Hands-on Course in Computational Structural Biology and Molecular Simulation BIOP590C/MCB590C Emad Tajkhorshid Center for Computational Biology and Biophysics Email: emad@life.uiuc.edu or tajkhors@uiuc.edu
More informationAMINO ACIDS ARE JOINED TOGETHER BY BONDS. FILE
04 February, 2018 AMINO ACIDS ARE JOINED TOGETHER BY BONDS. FILE Document Filetype: PDF 401.79 KB 0 AMINO ACIDS ARE JOINED TOGETHER BY BONDS. FILE Understanding the chemistry and properties of the peptide
More informationPatrick: An Introduction to Medicinal Chemistry 5e MOLECULAR MODELLING EXERCISES CHAPTER 17
MOLECULAR MODELLING EXERCISES CHAPTER 17 Exercise 17.6 Conformational analysis of n-butane Introduction Figure 1 Butane Me Me In this exercise, we will consider the possible stable conformations of butane
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 informationBiological Macromolecules
Introduction for Chem 493 Chemistry of Biological Macromolecules Dr. L. Luyt January 2008 Dr. L. Luyt Chem 493-2008 1 Biological macromolecules are the molecules of life allow for organization serve a
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 informationCreating a Pharmacophore Query from a Reference Molecule & Scaffold Hopping in CSD-CrossMiner
Table of Contents Creating a Pharmacophore Query from a Reference Molecule & Scaffold Hopping in CSD-CrossMiner Introduction... 2 CSD-CrossMiner Terminology... 2 Overview of CSD-CrossMiner... 3 Features
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 informationIntroduction to Protein Structures - Molecular Graphics Tool
Lecture 1a Introduction to Protein Structures - Molecular Graphics Tool amino acid tyrosine LH2 energetics traficking Ubiquitin enzymatic control BPTI VMD - a Tool to Think Carl Woese Atomic, CG, Particle:
More informationAll-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization
More 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 information3. An Introduction to Molecular Mechanics
3. An Introduction to Molecular Mechanics Introduction When you use Chem3D to draw molecules, the program assigns bond lengths and bond angles based on experimental data. The program does not contain real
More informationFigure 1. Molecules geometries of 5021 and Each neutral group in CHARMM topology was grouped in dash circle.
Project I Chemistry 8021, Spring 2005/2/23 This document was turned in by a student as a homework paper. 1. Methods First, the cartesian coordinates of 5021 and 8021 molecules (Fig. 1) are generated, in
More informationVersion 1.2 October 2017 CSD v5.39
Mogul Geometry Check Table of Contents Introduction... 2 Example 1. Using Mogul to assess intramolecular geometry... 3 Example 2. Using Mogul to explain activity data... 5 Conclusions... 8 Further Exercises...
More informationConformational Analysis of n-butane
Conformational Analysis of n-butane In this exercise you will calculate the Molecular Mechanics (MM) single point energy of butane in various conformations with respect to internal rotation around the
More informationAssignment 2: Conformation Searching (50 points)
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard September 16, 2015 Assignment 2: Conformation Searching (50 points) In this assignment, you will use the Spartan software package to investigate some conformation
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 informationPractical Manual. General outline to use the structural information obtained from molecular alignment
Practical Manual General outline to use the structural information obtained from molecular alignment 1. In order to use the information one needs to know the direction and the size of the tensor (susceptibility,
More informationICM-Chemist How-To Guide. Version 3.6-1g Last Updated 12/01/2009
ICM-Chemist How-To Guide Version 3.6-1g Last Updated 12/01/2009 ICM-Chemist HOW TO IMPORT, SKETCH AND EDIT CHEMICALS How to access the ICM Molecular Editor. 1. Click here 2. Start sketching How to sketch
More informationStructure of the HIV-1 Capsid Assembly by a hybrid approach
Structure of the HIV-1 Capsid Assembly by a hybrid approach Juan R. Perilla juan@ks.uiuc.edu Theoretical and Computational Biophysics Group University of Illinois at Urbana-Champaign Virion HIV infective
More informationDon t forget to bring your MD tutorial. Potential Energy (hyper)surface
Don t forget to bring your MD tutorial Lab session starts at 1pm You will have to finish an MD/SMD exercise on α-conotoxin in oxidized and reduced forms Potential Energy (hyper)surface What is Force? Energy
More informationMolecular Modeling and Conformational Analysis with PC Spartan
Molecular Modeling and Conformational Analysis with PC Spartan Introduction Molecular modeling can be done in a variety of ways, from using simple hand-held models to doing sophisticated calculations on
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 informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationIntroduction to ArcMap
Introduction to ArcMap ArcMap ArcMap is a Map-centric GUI tool used to perform map-based tasks Mapping Create maps by working geographically and interactively Display and present Export or print Publish
More informationImproved Resolution of Tertiary Structure Elasticity in Muscle Protein
Improved Resolution of Tertiary Structure Elasticity in Muscle Protein Jen Hsin and Klaus Schulten* Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois
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 informationFolding WT villin in silico Three folding simulations reach native state within 5-8 µs
Modeling of Cryo-EM Maps Workshop Baylor College of Medicine Klaus Schulten, U. Illinois at Urbana-Champaign Molecular Modeling Flexible Fitting 1: Introduction to Molecular Dynamics MD force field Equation
More informationConformations of Cycloheptadecane. A Comparison of Methods for Conformational Searching.
Conformations of Cycloheptadecane. A Comparison of Methods for Conformational Searching. M. Saunders, K. Houk, Y.-D. Wu, W.C. Still, M. Lipton, G. Chang, W.C. Guida J. Am. Chem. Soc. 1990, 112, 1419-1427
More informationCHEM 463: Advanced Inorganic Chemistry Modeling Metalloproteins for Structural Analysis
CHEM 463: Advanced Inorganic Chemistry Modeling Metalloproteins for Structural Analysis Purpose: The purpose of this laboratory is to introduce some of the basic visualization and modeling tools for viewing
More informationTutorial I: IQ MOL and Basic DFT and MP2 Calculations 1 / 30
Tutorial I: IQ MOL and Basic DFT and MP2 Calculations Q-Chem User Workshop, Denver March 21, 2015 1 / 30 2 / 30 Introduction to IQMOL DFT and MP2 Calculations 3 / 30 IQMOL and Q-CHEM IQMOL is an open-source
More informationSupplementary Figure 1 Preparation of PDA nanoparticles derived from self-assembly of PCDA. (a)
Supplementary Figure 1 Preparation of PDA nanoparticles derived from self-assembly of PCDA. (a) Computer simulation on the self-assembly of PCDAs. Two kinds of totally different initial conformations were
More informationSecondary and sidechain structures
Lecture 2 Secondary and sidechain structures James Chou BCMP201 Spring 2008 Images from Petsko & Ringe, Protein Structure and Function. Branden & Tooze, Introduction to Protein Structure. Richardson, J.
More informationSilica surface - Materials Studio tutorial. CREATING SiO 2 SURFACE
Silica surface - Materials Studio tutorial CREATING SiO 2 SURFACE Our goal surface of SiO2 6.948 Ǻ Import structure The XRD experiment gives us such parameters as: lattice parameters, symmetry group and
More informationIntroduction to Spark
1 As you become familiar or continue to explore the Cresset technology and software applications, we encourage you to look through the user manual. This is accessible from the Help menu. However, don t
More informationTutorial. Getting started. Sample to Insight. March 31, 2016
Getting started March 31, 2016 Sample to Insight CLC bio, a QIAGEN Company Silkeborgvej 2 Prismet 8000 Aarhus C Denmark Telephone: +45 70 22 32 44 www.clcbio.com support-clcbio@qiagen.com Getting started
More informationMolecular Simulation II
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 informationTo learn how to use molecular modeling software, a commonly used tool in the chemical and pharmaceutical industry.
NAME: Lab Day/Time: Molecular Modeling BV 1/2009 Purpose The purposes of this experiment are: To learn how to use molecular modeling software, a commonly used tool in the chemical and pharmaceutical industry.
More informationAssignment 2 Atomic-Level Molecular Modeling
Assignment 2 Atomic-Level Molecular Modeling CS/BIOE/CME/BIOPHYS/BIOMEDIN 279 Due: November 3, 2016 at 3:00 PM The goal of this assignment is to understand the biological and computational aspects of macromolecular
More informationMolecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains
Molecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains Fatemeh Khalili-Araghi, Emad Tajkhorshid, Benoît Roux, and Klaus Schulten Department of Physics, Department of Biochemistry,
More informationLast updated: Copyright
Last updated: 2012-08-20 Copyright 2004-2012 plabel (v2.4) User s Manual by Bioinformatics Group, Institute of Computing Technology, Chinese Academy of Sciences Tel: 86-10-62601016 Email: zhangkun01@ict.ac.cn,
More informationTailoring the Properties of Quadruplex Nucleobases for Biological and Nanomaterial Applications
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Supporting Information for: Tailoring the Properties of Quadruplex Nucleobases
More informationRanjit 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 informationTHE 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 informationThe Computational Microscope
The Computational Microscope Computational microscope views at atomic resolution... Rs SER RER C E M RER N GA L... how living cells maintain health and battle disease John Stone Our Microscope is Made
More informationThis is an author produced version of Privateer: : software for the conformational validation of carbohydrate structures.
This is an author produced version of Privateer: : software for the conformational validation of carbohydrate structures. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/95794/
More informationIon Translocation Across Biological Membranes. Janos K. Lanyi University of California, Irvine
Ion Translocation Across Biological Membranes Janos K. Lanyi University of California, Irvine Examples of transmembrane ion pumps Protein Cofactor, substrate, etc. MW Subunits Mitoch. cytochrome oxidase
More informationMultiscale 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 informationMarvin. Sketching, viewing and predicting properties with Marvin - features, tips and tricks. Gyorgy Pirok. Solutions for Cheminformatics
Marvin Sketching, viewing and predicting properties with Marvin - features, tips and tricks Gyorgy Pirok Solutions for Cheminformatics The Marvin family The Marvin toolkit provides web-enabled components
More informationBuild_model v User Guide
Build_model v.2.0.1 User Guide MolTech Build_model User Guide 2008-2011 Molecular Technologies Ltd. www.moltech.ru Please send your comments and suggestions to contact@moltech.ru. Table of Contents Input
More informationOpenDiscovery: 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 informationGraphical User Interfaces for Emittance and Correlation Plot. Henrik Loos
Graphical User Interfaces for Emittance and Correlation Plot Common GUI Features Overview Files Configs Measure Data Point Export Common GUI Features Save Saves the present data. Auto-generated file name
More information----- Ver October 24, 2014 Bug about reading MOPAC2012 Ver.14 calculations of 1 atom and 2 atoms molecule was fixed.
***** Facio's Release History ***** ----- Ver.18.8.2 ----- October 24, 2014 Bug about reading MOPAC2012 Ver.14 calculations of 1 atom and 2 atoms molecule was fixed. ----- Ver.18.8.1 ----- August 14, 2014
More informationMassHunter TOF/QTOF Users Meeting
MassHunter TOF/QTOF Users Meeting 1 Qualitative Analysis Workflows Workflows in Qualitative Analysis allow the user to only see and work with the areas and dialog boxes they need for their specific tasks
More informationNIH Center for Macromolecular Modeling and Bioinformatics Developer of VMD and NAMD. Beckman Institute
NIH Center for Macromolecular Modeling and Bioinformatics Developer of VMD and NAMD 5 faculty members (2 physics, 1 chemistry, 1 biochemistry, 1 computer science); 8 developers; 1 system admin; 15 post
More informationPerforming a Pharmacophore Search using CSD-CrossMiner
Table of Contents Introduction... 2 CSD-CrossMiner Terminology... 2 Overview of CSD-CrossMiner... 3 Searching with a Pharmacophore... 4 Performing a Pharmacophore Search using CSD-CrossMiner Version 2.0
More informationComputational Molecular Modeling
Computational Molecular Modeling Lecture 1: Structure Models, Properties Chandrajit Bajaj Today s Outline Intro to atoms, bonds, structure, biomolecules, Geometry of Proteins, Nucleic Acids, Ribosomes,
More informationUser Guide for LeDock
User Guide for LeDock Hongtao Zhao, PhD Email: htzhao@lephar.com Website: www.lephar.com Copyright 2017 Hongtao Zhao. All rights reserved. Introduction LeDock is flexible small-molecule docking software,
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 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 information