DOCSLIB.ORG

THE MOLECULAR DYNAMICS PACKAGE Q Version 4.20 Revised August 22, 2000

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
THE MOLECULAR DYNAMICS PACKAGE Q Version 4.20 Revised August 22, 2000 MANUAL FOR THE MOLECULAR DYNAMICS PACKAGE Q version 4.20 Revised August 22, 2000 John Marelius Karin Kolmodin Ellen Kindlund Johan Åqvist* <[email protected]> Dept. of Cell and Molecular Biology Uppsala University, Biomedical Centre Box 596 SE-751 24 Uppsala SWEDEN Q web site:site: http://aqvist.bmc.uu.se/Q 1 Contents A. Introduction.............................................................................................................................................. 4 B. User guide ................................................................................................................................................ 5 B.1. Preparing a molecular topology.......................................................................................................... 6 B.1.1. Preparing coordinates................................................................................................................ 6 B.1.2. Selecting a force field................................................................................................................ 7 B.1.3. Adding force field library entries............................................................................................... 7 B.1.4. Running Qprep ......................................................................................................................... 7 B.2. Running dynamics with Qdyn ......................................................................................................... 10 B.2.1. Simulation procedure .............................................................................................................. 10 B.2.2. Output generated by Qdyn....................................................................................................... 11 B.2.3. Preparing Qdyn input files ...................................................................................................... 11 B.2.4. Qdyn input file examples......................................................................................................... 14 B.2.5. FEP file .................................................................................................................................. 16 B.2.6. Monitoring non-bonded interactions ........................................................................................ 21 B.3. Analysis of results........................................................................................................................... 22 B.3.1. Analysing structures from the simulation................................................................................. 22 B.3.2. Free energy calculation using Qfep.......................................................................................... 23 B.3.3. Reaction free energy averaging using avefep............................................................................ 25 B.3.4. Simple energy averaging using Qave........................................................................................ 26 B.3.5. Trajectory averaging using Qavetr............................................................................................ 26 B.4. Useful tips ...................................................................................................................................... 26 C. Tutorials ................................................................................................................................................. 27 C.1. Binding affinity of the FABP-stearate complex from LIE simulations ............................................... 27 C.1.1. Editing the PDB file................................................................................................................ 27 C.1.2. Modelling ionic groups of the protein....................................................................................... 27 C.1.3. Writing the library file ............................................................................................................ 28 C.1.4. Qprep...................................................................................................................................... 28 C.1.5. FEP files................................................................................................................................. 29 C.1.6. Creating input files.................................................................................................................. 29 C.1.7. Using Qdyn3........................................................................................................................... 29 C.1.8. Evaluating the simulation........................................................................................................ 29 D. Reference guide ..................................................................................................................................... 31 D.1. Force field reference information .................................................................................................... 31 D.2. Topology preparation reference ....................................................................................................... 33 D.2.1. Coordinate files for input into Qprep ....................................................................................... 33 D.2.2. Atom masks ............................................................................................................................ 34 D.2.3. Qprep commands .................................................................................................................... 35 D.2.4. Qprep preferences................................................................................................................... 36 D.2.5. Fragment library file format..................................................................................................... 36 D.2.6. Force field parameter file format ............................................................................................. 39 D.2.7. Solvent file format................................................................................................................... 43 D.3. Boundary conditions ....................................................................................................................... 43 D.3.1. Solute boundary restraints....................................................................................................... 43 D.3.2. Water boundary restraints ....................................................................................................... 44 D.4. Units .............................................................................................................................................. 45 D.5. File and format descriptions............................................................................................................ 45 D.5.1. Qdyn input file format............................................................................................................. 45 D.5.2. FEP file format ....................................................................................................................... 49 D.6. Log files.......................................................................................................................................... 55 D.6.1. Qdyn log file ........................................................................................................................... 55 D.6.2. Qfep log file............................................................................................................................ 59 D.7. Differences from earlier versions..................................................................................................... 59 D.7.1. Differences in Qprep............................................................................................................... 59 2 D.7.2. Differences in Qdyn ................................................................................................................ 60 D.8. Utility programs.............................................................................................................................. 60 D.8.1. proq - PDB file analysis and manipulation .............................................................................. 60 D.8.2. mdsh - input preparation shell ................................................................................................ 64 D.8.3. Analysis programs................................................................................................................... 67 D.8.4. Other utility programs ............................................................................................................. 68 E. System guide .......................................................................................................................................... 70 E.1. System requirements ....................................................................................................................... 70 E.1.1. Compilers................................................................................................................................ 70 E.2. Installation...................................................................................................................................... 70 E.2.1. Installing executable images.................................................................................................... 70 E.2.2. Building the programs from source code .................................................................................
Load more

Recommended publications
  • Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development
    processes Review Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development Outi M. H. Salo-Ahen 1,2,* , Ida Alanko 1,2, Rajendra Bhadane 1,2 , Alexandre M. J. J. Bonvin 3,* , Rodrigo Vargas Honorato 3, Shakhawath Hossain 4 , André H. Juffer 5 , Aleksei Kabedev 4, Maija Lahtela-Kakkonen 6, Anders Støttrup Larsen 7, Eveline Lescrinier 8 , Parthiban Marimuthu 1,2 , Muhammad Usman Mirza 8 , Ghulam Mustafa 9, Ariane Nunes-Alves 10,11,* , Tatu Pantsar 6,12, Atefeh Saadabadi 1,2 , Kalaimathy Singaravelu 13 and Michiel Vanmeert 8 1 Pharmaceutical Sciences Laboratory (Pharmacy), Åbo Akademi University, Tykistökatu 6 A, Biocity, FI-20520 Turku, Finland; ida.alanko@abo.fi (I.A.); rajendra.bhadane@abo.fi (R.B.); parthiban.marimuthu@abo.fi (P.M.); atefeh.saadabadi@abo.fi (A.S.) 2 Structural Bioinformatics Laboratory (Biochemistry), Åbo Akademi University, Tykistökatu 6 A, Biocity, FI-20520 Turku, Finland 3 Faculty of Science-Chemistry, Bijvoet Center for Biomolecular Research, Utrecht University, 3584 CH Utrecht, The Netherlands; [email protected] 4 Swedish Drug Delivery Forum (SDDF), Department of Pharmacy, Uppsala Biomedical Center, Uppsala University, 751 23 Uppsala, Sweden; [email protected] (S.H.); [email protected] (A.K.) 5 Biocenter Oulu & Faculty of Biochemistry and Molecular Medicine, University of Oulu, Aapistie 7 A, FI-90014 Oulu, Finland; andre.juffer@oulu.fi 6 School of Pharmacy, University of Eastern Finland, FI-70210 Kuopio, Finland; maija.lahtela-kakkonen@uef.fi (M.L.-K.); tatu.pantsar@uef.fi
    [Show full text]
  • Cygwin User's Guide
    Cygwin User’s Guide Cygwin User’s Guide ii Copyright © Cygwin authors Permission is granted to make and distribute verbatim copies of this documentation provided the copyright notice and this per- mission notice are preserved on all copies. Permission is granted to copy and distribute modified versions of this documentation under the conditions for verbatim copying, provided that the entire resulting derived work is distributed under the terms of a permission notice identical to this one. Permission is granted to copy and distribute translations of this documentation into another language, under the above conditions for modified versions, except that this permission notice may be stated in a translation approved by the Free Software Foundation. Cygwin User’s Guide iii Contents 1 Cygwin Overview 1 1.1 What is it? . .1 1.2 Quick Start Guide for those more experienced with Windows . .1 1.3 Quick Start Guide for those more experienced with UNIX . .1 1.4 Are the Cygwin tools free software? . .2 1.5 A brief history of the Cygwin project . .2 1.6 Highlights of Cygwin Functionality . .3 1.6.1 Introduction . .3 1.6.2 Permissions and Security . .3 1.6.3 File Access . .3 1.6.4 Text Mode vs. Binary Mode . .4 1.6.5 ANSI C Library . .4 1.6.6 Process Creation . .5 1.6.6.1 Problems with process creation . .5 1.6.7 Signals . .6 1.6.8 Sockets . .6 1.6.9 Select . .7 1.7 What’s new and what changed in Cygwin . .7 1.7.1 What’s new and what changed in 3.2 .
    [Show full text]
  • The Focus - Issue 36
    Contents The Focus - Issue 36 A Publication for ANSYS Users Contents Feature Articles ● Linux & ANSYS: Lessons Learned ● Backup Tool ● Design Modeler FAQ On the Web ● APDL Customization course notes now available for purchase ● ANSYS and MathCAD ● ANSYS Acquires Century Dynamics Resources ● PADT Support: How can we help? ● Upcoming Training at PADT ● About The Focus ❍ The Focus Library ❍ Contributor Information ❍ Subscribe / Unsubscribe ❍ Legal Disclaimer http://www.padtinc.com/epubs/focus/common/contents.asp [3/28/2005 9:06:12 AM] Linux & ANSYS: Lessons Learned The Focus - Issue 36 A Publication for ANSYS Users Linux & ANSYS: Lessons Learned by Eric Miller, PADT Every couple of years, the computing picture for analysts gets turned upside down. For a long time now the industry has been moving from Unix workstations to Windows/Intel desktop machines. The wintel price/performance has been fantastic, the IT guys are happier, and all of that productivity software that you spend so much time with runs in the same spot. We have been happy with a stable and known environment. However, accepting the fact that unless you work for a big company that can buy some Unix servers, you just don’t have an easy way to get some extra horsepower other then getting a new box. Then along comes this Finnish guy that may or may not have been named after Lucy’s little brother. With not much of a life and a very large brain, he popped out the majority of a complete and free version of Unix that anyone can use, breaking the stranglehold of (expensive) proprietary Unix OS’s that ran on (expensive) proprietary hardware.
    [Show full text]
  • Development and Application of a Computational Platform for Complex Molecular Design Jaime Rodríguez-Guerra Pedregal
    ADVERTIMENT. Lʼaccés als continguts dʼaquesta tesi queda condicionat a lʼacceptació de les condicions dʼús establertes per la següent llicència Creative Commons: http://cat.creativecommons.org/?page_id=184 ADVERTENCIA. El acceso a los contenidos de esta tesis queda condicionado a la aceptación de las condiciones de uso establecidas por la siguiente licencia Creative Commons: http://es.creativecommons.org/blog/licencias/ WARNING. The access to the contents of this doctoral thesis it is limited to the acceptance of the use conditions set by the following Creative Commons license: https://creativecommons.org/licenses/?lang=en Development and Application of a Computational Platform for Complex Molecular Design a dissertation submitted by Jaime Rodríguez-Guerra Pedregal & directed by Prof. Dr. Jean-Didier Maréchal in fulfillment of the requirements for the degree of Doctor of Biotechnology Tutor: Prof. Dr. Jordi Joan Cairó Badillo Department of Chemical, Biological and Environmental Engineering Universitat Autònoma de Barcelona July 2018 Development and Application of a Computational Platform for Complex Molecular Design a dissertation submitted by & recommended for acceptance by advisor Jaime Rodríguez-Guerra Pedregal Prof. Dr. Jean-Didier Maréchal Tutor: Prof. Dr. Jordi Joan Cairó Badillo Department of Chemical, Biological and Environmental Engineering Universitat Autònoma de Barcelona July 2018 ©2018 – Jaime Rodríguez-Guerra Pedregal Licensed as Creative Commons BY-NC-ND Attribution-NonCommercial-NoDerivs In the beginning, there was nothing. And God said «Let there be light». And there was light. There was still nothing, but you could see it a lot better. —WoodyAllen. Development and Application of a Computational Platform for Complex Molecular Design by Jaime Rodríguez-Guerra Pedregal Abstract In this dissertation, a series of novel computational modeling tools is reported.
    [Show full text]
  • Cygwin User's Guide
    Cygwin User’s Guide i Cygwin User’s Guide Cygwin User’s Guide ii Copyright © 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012 Red Hat, Inc. Permission is granted to make and distribute verbatim copies of this documentation provided the copyright notice and this per- mission notice are preserved on all copies. Permission is granted to copy and distribute modified versions of this documentation under the conditions for verbatim copying, provided that the entire resulting derived work is distributed under the terms of a permission notice identical to this one. Permission is granted to copy and distribute translations of this documentation into another language, under the above conditions for modified versions, except that this permission notice may be stated in a translation approved by the Free Software Foundation. Cygwin User’s Guide iii Contents 1 Cygwin Overview 1 1.1 What is it? . .1 1.2 Quick Start Guide for those more experienced with Windows . .1 1.3 Quick Start Guide for those more experienced with UNIX . .1 1.4 Are the Cygwin tools free software? . .2 1.5 A brief history of the Cygwin project . .2 1.6 Highlights of Cygwin Functionality . .3 1.6.1 Introduction . .3 1.6.2 Permissions and Security . .3 1.6.3 File Access . .3 1.6.4 Text Mode vs. Binary Mode . .4 1.6.5 ANSI C Library . .5 1.6.6 Process Creation . .5 1.6.6.1 Problems with process creation . .5 1.6.7 Signals . .6 1.6.8 Sockets . .6 1.6.9 Select .
    [Show full text]
  • PARCE: Protocol for Amino Acid Refinement Through Computational
    PARCE: Protocol for Amino acid Refinement through Computational Evolution Rodrigo Ochoaa, Miguel A. Solerb, Alessandro Laioc,d, Pilar Cossioa,e,∗ aBiophysics of Tropical Diseases, Max Planck Tandem Group, University of Antioquia, 050010 Medellin, Colombia bItalian Institute of Technology (IIT), Via Melen 83, B Block, 16152, Genova, Italy cInternational School for Advanced Studies (SISSA), Via Bonomea 265, I-34136 Trieste, Italy dThe Abdus Salam International Centre for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy eDepartment of Theoretical Biophysics, Max Planck Institute of Biophysics, 60438 Frankfurt am Main, Germany Abstract The in silico design of peptides and proteins as binders is useful for diagnosis and therapeutics due to their low adverse effects and major specificity. To select the most promising candidates, a key matter is to understand their interactions with protein targets. In this work, we present PARCE, an open source Protocol for Amino acid Refinement through Computational Evolution that implements an advanced and promising method for the design of peptides and proteins. The protocol performs a random mutation in the binder sequence, then samples the bound conformations using molecular dynamics simulations, and evaluates the protein-protein interactions from multiple scoring. Finally, it accepts or rejects the mutation by applying a consensus criterion based on binding scores. The procedure is iterated with the aim to explore efficiently novel sequences with potential better affinities toward their targets.
    [Show full text]
  • Manual-2018.Pdf
    GROMACS Groningen Machine for Chemical Simulations Reference Manual Version 2018 GROMACS Reference Manual Version 2018 Contributions from Emile Apol, Rossen Apostolov, Herman J.C. Berendsen, Aldert van Buuren, Pär Bjelkmar, Rudi van Drunen, Anton Feenstra, Sebastian Fritsch, Gerrit Groenhof, Christoph Junghans, Jochen Hub, Peter Kasson, Carsten Kutzner, Brad Lambeth, Per Larsson, Justin A. Lemkul, Viveca Lindahl, Magnus Lundborg, Erik Marklund, Pieter Meulenhoff, Teemu Murtola, Szilárd Páll, Sander Pronk, Roland Schulz, Michael Shirts, Alfons Sijbers, Peter Tieleman, Christian Wennberg and Maarten Wolf. Mark Abraham, Berk Hess, David van der Spoel, and Erik Lindahl. c 1991–2000: Department of Biophysical Chemistry, University of Groningen. Nijenborgh 4, 9747 AG Groningen, The Netherlands. c 2001–2018: The GROMACS development teams at the Royal Institute of Technology and Uppsala University, Sweden. More information can be found on our website: www.gromacs.org. iv Preface & Disclaimer This manual is not complete and has no pretention to be so due to lack of time of the contributors – our first priority is to improve the software. It is worked on continuously, which in some cases might mean the information is not entirely correct. Comments on form and content are welcome, please send them to one of the mailing lists (see www.gromacs.org), or open an issue at redmine.gromacs.org. Corrections can also be made in the GROMACS git source repository and uploaded to gerrit.gromacs.org. We release an updated version of the manual whenever we release a new version of the software, so in general it is a good idea to use a manual with the same major and minor release number as your GROMACS installation.
    [Show full text]
  • Open Source Molecular Modeling
    Accepted Manuscript Title: Open Source Molecular Modeling Author: Somayeh Pirhadi Jocelyn Sunseri David Ryan Koes PII: S1093-3263(16)30118-8 DOI: http://dx.doi.org/doi:10.1016/j.jmgm.2016.07.008 Reference: JMG 6730 To appear in: Journal of Molecular Graphics and Modelling Received date: 4-5-2016 Accepted date: 25-7-2016 Please cite this article as: Somayeh Pirhadi, Jocelyn Sunseri, David Ryan Koes, Open Source Molecular Modeling, <![CDATA[Journal of Molecular Graphics and Modelling]]> (2016), http://dx.doi.org/10.1016/j.jmgm.2016.07.008 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. Open Source Molecular Modeling Somayeh Pirhadia, Jocelyn Sunseria, David Ryan Koesa,∗ aDepartment of Computational and Systems Biology, University of Pittsburgh Abstract The success of molecular modeling and computational chemistry efforts are, by definition, de- pendent on quality software applications. Open source software development provides many advantages to users of modeling applications, not the least of which is that the software is free and completely extendable. In this review we categorize, enumerate, and describe available open source software packages for molecular modeling and computational chemistry. 1. Introduction What is Open Source? Free and open source software (FOSS) is software that is both considered \free software," as defined by the Free Software Foundation (http://fsf.org) and \open source," as defined by the Open Source Initiative (http://opensource.org).
    [Show full text]
  • Introduction to QM/MM Simulations
    Chapter 3 Introduction to QM/MM Simulations Gerrit Groenhof Abstract Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations have become a popular tool for investigating chemical reactions in condensed phases. In QM/MM methods, the region of the system in which the chemical process takes place is treated at an appropriate level of quantum chemistry theory, while the remainder is described by a molecular mechanics force field. Within this approach, chemical reactivity can be studied in large systems, such as enzymes. In the first part of this contribution, the basic methodol- ogy is briefly reviewed. The two most common approaches for partitioning the two subsystems are presented, followed by a discussion on the different ways of treating interactions between the subsystems. Special attention is given on how to deal with situations in which the boundary between the QM and MM subsystems runs through one or more chemical bonds. The second part of this contribution discusses what properties of larger system can be obtained within the QM/MM framework and how. Finally, as an example of a QM/MM application in practice, the third part presents an overview of recent QM/MM molecular dynamics simulations on photobiological systems. In addition to providing quantities that are experimen- tally accessible, such as structural intermediates, fluorescence lifetimes, quantum yields and spectra, the QM/MM simulations also provide information that is much more difficult to measure experimentally, such as reaction mechanisms and the influence of individual amino acid residues. Key word: Quantum mechanics, Molecular mechanics, QM/MM, Molecular dynamics 1. Introduction In this chapter we present a short introduction into the development and application of computational techniques for modelling chemi- cal reactions in the condensed phase.
    [Show full text]
  • Molecular Dynamics Study of Proton Transfer Reactions in GROMACS
    Molecular Dynamics Study ofProton Transfer Reactions inGROMACS Özge Yoluk Degree project inapplied biotechnology, Master ofScience (2years), 2011 Examensarbete itillämpad bioteknik 45 hp tillmasterexamen, 2011 Biology Education Centre and Department ofCell and Molecular Biology, Uppsala University Supervisors: David van der Spoel and Erik G.Marklund Dedicated to my beloved parents; Emirali and Fidan Yoluk i Figures Figure 1. Grotthuss mechanism. Green dashed lines indicate the transfer. ............................... 1 Figure 2. One dimensional potential energy profile of a proton-transfer reaction. ................... 3 Figure 3. Schematic repsresentation of the regimes in proton-transfer reactions ...................... 4 Figure 4.Schematic representations of the molecules and proton-transfer routes. .................... 8 Figure 5. Flowchart of a step in Q-Hop algorithm..................................................................... 8 Figure 6. A simple flowchart for a molecular dynamics simulation with GROMACS ........... 10 Figure 7. Empirical cumulative distribution function plots with the first (right) and the second (left) set of reactions. ............................................................................................ 12 Figure 8. Radial distribution function of inert hydronium in different water models with or without proton-transfer. .................................................................................................... 13 Figure 9. Rdf of Cl- OW (I) and OW*-OW (II) in simulations with or without
    [Show full text]
  • A Concise Guide to CHARMM and the Analysis of Protein Structure and Function
    A Concise Guide to CHARMM and the Analysis of Protein Structure and Function Robert Schleif Biology Department Johns Hopkins University 3400 N. Charles St. Baltimore, MD 21218 1/8/06 9/17/13 Preface Increasingly, biologists and biochemists are faced with understanding how their favorite proteins work. The structures of many of these proteins have been determined, and the structures of many more will be determined in the next few years. Once a protein's structure has been determined, it becomes possible and also enticing to design experiments probing the protein's mechanism of action. Tools for the graphical display of structure, the manipulation of the structure, and the calculation of various interaction energies all become interesting and important. Additionally, some properties of a protein may best be revealed by modeling the protein in water and simulating its molecular thermal motion at 300 K. A researcher interested in protein structure and function faces the question of whether to use one of the complete, but expensive, computer programs for the manipulation and analysis of protein structure, use a number of the highly specialized but almost completely undocumented programs that are available on the web, or to learn and use a powerful and general program that can perform most of the manipulations and calculations one might need. This book is written for those who decide to follow the latter course and to learn the program CHARMM (Chemistry at Harvard Macromolecular Mechanics) that was initiated in the laboratory of Dr. Martin Karplus. The program has been continuously refined and extended by many workers over the years since the initial publication, "CHARMM: A Program for Macromolecular Energy, Minimization, and Dynamics Calculations", J.
    [Show full text]
  • Tutorial-FEP.Pdf
    Universite´ de Lorraine Centre National de la Recherche Scientifique Laboratoire International Associe´ CNRS-University of Illinois Centre National de la Recherche Scientifique Institut de Biologie Physico-chimique University of Illinois at Urbana-Champaign Beckman Institute for Advanced Science and Technology Theoretical and Computational Biophysics Group In silico alchemy: A tutorial for alchemical free-energy perturbation calculations with NAMD Jer´ omeˆ Henin´ James Gumbart Christophe Chipot September 20, 2017 Current editors: Abhishek Singharoy & Ivan Teo Please visit www.ks.uiuc.edu/Training/Tutorials/ to get the latest version of this tutorial, to obtain more tutorials like this one, or to join the [email protected] mailing list for additional help. Abstract This tutorial explains how NAMD and related tools can be used to setup and perform alchemical free- energy simulations within the free-energy perturbation (FEP) theory. The force-field independent, “zero-sum” transformation of ethane into ethane is used as an introductory, prototypical example. FEP is then used to compute the free energy of charging a naked Lennard-Jones particle into a sodium ion. Next, the variation in solvation free energy upon mutation of a tyrosine residue into alanine is examined in the Ala–Tyr–Ala tripeptide. Last, the concept of standard binding free energy is illustrated in the simple case of a potassium ion binding a ionophore, 18–crown–6. Prior knowledge of NAMD and standard molecular dynamics simulations is assumed. c 2017, Centre National de la Recherche Scientifique, University of Illinois, Urbana–Champaign 2 Alchemical free-energy calculations tutorial 3 Contents 1. Ethane-to-ethane “zero-sum” transformation 5 1.1.
    [Show full text]