GROMACS Reference Manual

Total Page:16

File Type:pdf, Size:1020Kb

GROMACS Reference Manual GROMACS Groningen Machine for Chemical Simulations Reference Manual Version 2016.3 GROMACS Reference Manual Version 2016.3 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, Erik Marklund, Peiter 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–2017: 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. On-line Resources You can find more documentation and other material at our homepage www.gromacs.org. Among other things there is an on-line reference, several GROMACS mailing lists with archives and contributed topologies/force fields. Citation information When citing this document in any scientific publication please refer to it as: M.J. Abraham, D. van der Spoel, E. Lindahl, B. Hess, and the GROMACS development team, GROMACS User Manual version 2016.3, www.gromacs.org (2017) However, we prefer that you cite (some of) the GROMACS papers [1,2,3,4,5,6,7,8] when you publish your results. Any future development depends on academic research grants, since the package is distributed as free software! GROMACS is Free Software The entire GROMACS package is available under the GNU Lesser General Public License (LGPL), version 2.1. This means it’s free as in free speech, not just that you can use it without pay- ing us money. You can redistribute GROMACS and/or modify it under the terms of the LGPL as published by the Free Software Foundation; either version 2.1 of the License, or (at your option) any later version. For details, check the COPYING file in the source code or consult http://www.gnu.org/licenses/old-licenses/lgpl-2.1.html. The GROMACS source code and and selected set of binary packages are available on our home- page, www.gromacs.org. Have fun. Contents 1 Introduction1 1.1 Computational Chemistry and Molecular Modeling................1 1.2 Molecular Dynamics Simulations..........................2 1.3 Energy Minimization and Search Methods.....................5 2 Definitions and Units7 2.1 Notation.......................................7 2.2 MD units......................................7 2.3 Reduced units....................................9 2.4 Mixed or Double precision............................. 10 3 Algorithms 11 3.1 Introduction..................................... 11 3.2 Periodic boundary conditions............................ 11 3.2.1 Some useful box types........................... 13 3.2.2 Cut-off restrictions............................. 14 3.3 The group concept................................. 14 3.4 Molecular Dynamics................................ 15 3.4.1 Initial conditions.............................. 17 3.4.2 Neighbor searching............................. 18 3.4.3 Compute forces............................... 24 3.4.4 The leap-frog integrator.......................... 26 3.4.5 The velocity Verlet integrator....................... 26 3.4.6 Understanding reversible integrators: The Trotter decomposition..... 27 3.4.7 Multiple time stepping........................... 29 3.4.8 Temperature coupling........................... 29 vi Contents 3.4.9 Pressure coupling.............................. 35 3.4.10 The complete update algorithm...................... 41 3.4.11 Output step................................. 43 3.5 Shell molecular dynamics.............................. 43 3.5.1 Optimization of the shell positions..................... 43 3.6 Constraint algorithms................................ 44 3.6.1 SHAKE................................... 44 3.6.2 LINCS................................... 45 3.7 Simulated Annealing................................ 47 3.8 Stochastic Dynamics................................ 48 3.9 Brownian Dynamics................................ 49 3.10 Energy Minimization................................ 49 3.10.1 Steepest Descent.............................. 49 3.10.2 Conjugate Gradient............................. 50 3.10.3 L-BFGS................................... 50 3.11 Normal-Mode Analysis............................... 51 3.12 Free energy calculations.............................. 51 3.12.1 Slow-growth methods........................... 51 3.12.2 Thermodynamic integration........................ 53 3.13 Replica exchange.................................. 54 3.14 Essential Dynamics sampling............................ 55 3.15 Expanded Ensemble................................. 56 3.16 Parallelization.................................... 56 3.17 Domain decomposition............................... 57 3.17.1 Coordinate and force communication................... 57 3.17.2 Dynamic load balancing.......................... 58 3.17.3 Constraints in parallel........................... 59 3.17.4 Interaction ranges.............................. 60 3.17.5 Multiple-Program, Multiple-Data PME parallelization.......... 61 3.17.6 Domain decomposition flow chart..................... 62 3.18 Implicit solvation.................................. 62 4 Interaction function and force fields 65 4.1 Non-bonded interactions.............................. 65 Contents vii 4.1.1 The Lennard-Jones interaction....................... 66 4.1.2 Buckingham potential........................... 67 4.1.3 Coulomb interaction............................ 67 4.1.4 Coulomb interaction with reaction field.................. 68 4.1.5 Modified non-bonded interactions..................... 69 4.1.6 Modified short-range interactions with Ewald summation......... 71 4.2 Bonded interactions................................. 71 4.2.1 Bond stretching............................... 71 4.2.2 Morse potential bond stretching...................... 72 4.2.3 Cubic bond stretching potential...................... 73 4.2.4 FENE bond stretching potential...................... 74 4.2.5 Harmonic angle potential.......................... 74 4.2.6 Cosine based angle potential........................ 75 4.2.7 Restricted bending potential........................ 75 4.2.8 Urey-Bradley potential........................... 76 4.2.9 Bond-Bond cross term........................... 77 4.2.10 Bond-Angle cross term........................... 77 4.2.11 Quartic angle potential........................... 77 4.2.12 Improper dihedrals............................. 77 4.2.13 Proper dihedrals.............................. 79 4.2.14 Tabulated bonded interaction functions.................. 82 4.3 Restraints...................................... 84 4.3.1 Position restraints.............................. 84 4.3.2 Flat-bottomed position restraints...................... 85 4.3.3 Angle restraints............................... 86 4.3.4 Dihedral restraints............................. 86 4.3.5 Distance restraints............................. 87 4.3.6 Orientation restraints............................ 90 4.4 Polarization..................................... 94 4.4.1 Simple polarization............................. 94 4.4.2 Water polarization............................. 94 4.4.3 Thole polarization............................. 94 4.5 Free energy interactions............................... 94 4.5.1 Soft-core interactions............................ 98 viii Contents 4.6 Methods....................................... 99 4.6.1 Exclusions and 1-4 Interactions....................... 99 4.6.2 Charge Groups............................... 100 4.6.3 Treatment of Cut-offs in the group scheme................ 100 4.7 Virtual interaction sites............................... 101 4.8 Long Range Electrostatics............................. 105 4.8.1 Ewald summation.............................. 105 4.8.2 PME.................................... 106 4.8.3 P3M-AD.................................. 107 4.8.4 Optimizing Fourier transforms and PME calculations........... 107 4.9 Long Range Van der Waals interactions...................... 107 4.9.1 Dispersion correction............................ 107 4.9.2 Lennard-Jones PME............................ 109 4.10 Force field...................................... 112 4.10.1 GROMOS-96................................ 113 4.10.2 OPLS/AA.................................. 114 4.10.3 AMBER.................................. 114 4.10.4 CHARMM................................. 114 4.10.5 Coarse-grained force fields......................... 115 4.10.6 MARTINI.................................. 115 4.10.7 PLUM..................................
Recommended publications
  • MODELLER 10.1 Manual
    MODELLER A Program for Protein Structure Modeling Release 10.1, r12156 Andrej Saliˇ with help from Ben Webb, M.S. Madhusudhan, Min-Yi Shen, Guangqiang Dong, Marc A. Martı-Renom, Narayanan Eswar, Frank Alber, Maya Topf, Baldomero Oliva, Andr´as Fiser, Roberto S´anchez, Bozidar Yerkovich, Azat Badretdinov, Francisco Melo, John P. Overington, and Eric Feyfant email: modeller-care AT salilab.org URL https://salilab.org/modeller/ 2021/03/12 ii Contents Copyright notice xxi Acknowledgments xxv 1 Introduction 1 1.1 What is Modeller?............................................. 1 1.2 Modeller bibliography....................................... .... 2 1.3 Obtainingandinstallingtheprogram. .................... 3 1.4 Bugreports...................................... ............ 4 1.5 Method for comparative protein structure modeling by Modeller ................... 5 1.6 Using Modeller forcomparativemodeling. ... 8 1.6.1 Preparinginputfiles . ............. 8 1.6.2 Running Modeller ......................................... 9 2 Automated comparative modeling with AutoModel 11 2.1 Simpleusage ..................................... ............ 11 2.2 Moreadvancedusage............................... .............. 12 2.2.1 Including water molecules, HETATM residues, and hydrogenatoms .............. 12 2.2.2 Changing the default optimization and refinement protocol ................... 14 2.2.3 Getting a very fast and approximate model . ................. 14 2.2.4 Building a model from multiple templates . .................. 15 2.2.5 Buildinganallhydrogenmodel
    [Show full text]
  • Open Babel Documentation Release 2.3.1
    Open Babel Documentation Release 2.3.1 Geoffrey R Hutchison Chris Morley Craig James Chris Swain Hans De Winter Tim Vandermeersch Noel M O’Boyle (Ed.) December 05, 2011 Contents 1 Introduction 3 1.1 Goals of the Open Babel project ..................................... 3 1.2 Frequently Asked Questions ....................................... 4 1.3 Thanks .................................................. 7 2 Install Open Babel 9 2.1 Install a binary package ......................................... 9 2.2 Compiling Open Babel .......................................... 9 3 obabel and babel - Convert, Filter and Manipulate Chemical Data 17 3.1 Synopsis ................................................. 17 3.2 Options .................................................. 17 3.3 Examples ................................................. 19 3.4 Differences between babel and obabel .................................. 21 3.5 Format Options .............................................. 22 3.6 Append property values to the title .................................... 22 3.7 Filtering molecules from a multimolecule file .............................. 22 3.8 Substructure and similarity searching .................................. 25 3.9 Sorting molecules ............................................ 25 3.10 Remove duplicate molecules ....................................... 25 3.11 Aliases for chemical groups ....................................... 26 4 The Open Babel GUI 29 4.1 Basic operation .............................................. 29 4.2 Options .................................................
    [Show full text]
  • 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]
  • Ee9e60701e814784783a672a9
    International Journal of Technology (2017) 4: 611‐621 ISSN 2086‐9614 © IJTech 2017 A PRELIMINARY STUDY ON SHIFTING FROM VIRTUAL MACHINE TO DOCKER CONTAINER FOR INSILICO DRUG DISCOVERY IN THE CLOUD Agung Putra Pasaribu1, Muhammad Fajar Siddiq1, Muhammad Irfan Fadhila1, Muhammad H. Hilman1, Arry Yanuar2, Heru Suhartanto1* 1Faculty of Computer Science, Universitas Indonesia, Kampus UI Depok, Depok 16424, Indonesia 2Faculty of Pharmacy, Universitas Indonesia, Kampus UI Depok, Depok 16424, Indonesia (Received: January 2017 / Revised: April 2017 / Accepted: June 2017) ABSTRACT The rapid growth of information technology and internet access has moved many offline activities online. Cloud computing is an easy and inexpensive solution, as supported by virtualization servers that allow easier access to personal computing resources. Unfortunately, current virtualization technology has some major disadvantages that can lead to suboptimal server performance. As a result, some companies have begun to move from virtual machines to containers. While containers are not new technology, their use has increased recently due to the Docker container platform product. Docker’s features can provide easier solutions. In this work, insilico drug discovery applications from molecular modelling to virtual screening were tested to run in Docker. The results are very promising, as Docker beat the virtual machine in most tests and reduced the performance gap that exists when using a virtual machine (VirtualBox). The virtual machine placed third in test performance, after the host itself and Docker. Keywords: Cloud computing; Docker container; Molecular modeling; Virtual screening 1. INTRODUCTION In recent years, cloud computing has entered the realm of information technology (IT) and has been widely used by the enterprise in support of business activities (Foundation, 2016).
    [Show full text]
  • Parameterizing a Novel Residue
    University of Illinois at Urbana-Champaign Luthey-Schulten Group, Department of Chemistry Theoretical and Computational Biophysics Group Computational Biophysics Workshop Parameterizing a Novel Residue Rommie Amaro Brijeet Dhaliwal Zaida Luthey-Schulten Current Editors: Christopher Mayne Po-Chao Wen February 2012 CONTENTS 2 Contents 1 Biological Background and Chemical Mechanism 4 2 HisH System Setup 7 3 Testing out your new residue 9 4 The CHARMM Force Field 12 5 Developing Topology and Parameter Files 13 5.1 An Introduction to a CHARMM Topology File . 13 5.2 An Introduction to a CHARMM Parameter File . 16 5.3 Assigning Initial Values for Unknown Parameters . 18 5.4 A Closer Look at Dihedral Parameters . 18 6 Parameter generation using SPARTAN (Optional) 20 7 Minimization with new parameters 32 CONTENTS 3 Introduction Molecular dynamics (MD) simulations are a powerful scientific tool used to study a wide variety of systems in atomic detail. From a standard protein simulation, to the use of steered molecular dynamics (SMD), to modelling DNA-protein interactions, there are many useful applications. With the advent of massively parallel simulation programs such as NAMD2, the limits of computational anal- ysis are being pushed even further. Inevitably there comes a time in any molecular modelling scientist’s career when the need to simulate an entirely new molecule or ligand arises. The tech- nique of determining new force field parameters to describe these novel system components therefore becomes an invaluable skill. Determining the correct sys- tem parameters to use in conjunction with the chosen force field is only one important aspect of the process.
    [Show full text]
  • Computational Chemistry (F14CCH)
    Computational Chemistry (F14CCH) Connecting to Remote Machines Using CCP1GUI CONDOR is a queuing system that can be accessed from the Linux Left Mouse Button: Turn machines of the theoretical research groups and distributes the Hold right + up or down: Zoom submitted jobs to free compute nodes. Since the machines are protected Hold wheel + move: Move by firewalls you have to connect to a specific node, Porthos, on which you were given an account. Click Atom => Select Fragment => Add Assign all elements, also the hydrogens, before creating the GAMESS Use PuTTY to connect to remote machines, porthos in this case. Open input file (click on “All X => H” in the tool panel). If you cannot see one of putty.exe, enter the hostname as porthos.chem.nott.ac.uk. The login the windows of CCP1GUI, it might be hiding outside the desktop area. is your university username and the default password is Right click on the taskbar icon => Move => hold cursor keys to the left ChangeMeSoon. When you first log in on CONDOR you must change until the window appears. your password using kpasswd! It must be strong, including mixed case letters and numbers. Creating a GAMESS Input File / Running a Local Job Compute => GAMESS-UK Copying Files from/to Porthos To exchange files with porthos in order to submit files to CONDOR for Molecule: Options => Title calculation or copy results back, you first have to copy them from the Task => select requested task Windows machine to porthos via the SSHClient. It cannot be installed in Theory: Select SCF Method (and post SCF Method, if required) the CAL but can be accessed at NAL => Accessing the Internet => Optimisation: Runtype => Opt.
    [Show full text]
  • Molecular Simulation and Structure Prediction Using CHARMM and the MMTSB Tool Set Introduction
    Molecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Introduction Charles L. Brooks III MMTSB/CTBP 2006 Summer Workshop MMTSB/CTBP Summer Workshop © Charles L. Brooks III, 2006. CTBP and MMTSB: Who are we? • The CTBP is the Center for Theoretical Biological Physics – Funded by the NSF as a Physics Frontiers Center – Partnership between UCSD, Scripps and Salk, lead by UCSD • The CTBP encompasses a broad spectrum of research and training activities at the forefront of the biology- physics interface. – Principal scientists include • José Onuchic, Herbie Levine, Henry Abarbanel, Charles Brooks, David Case, Mike Holst, Terry Hwa, David Kleinfeld, Andy McCammon, Wouter Rappel, Terry Sejnowski, Wei Wang, Peter Wolynes MMTSB/CTBP Summer Workshop © Charles L. Brooks III, 2006. CTBP and MMTSB: Who are we? MMTSB/CTBP Summer Workshop http://ctbp.ucsd.edu © Charles L. Brooks III, 2006. CTBP and MMTSB: Who are we? • The MMTSB is the Center for Multi-scale Modeling Tools for Structural Biology – Funded by the NIH as a National Research Resource Center – Partnership between Scripps and Georgia Tech, lead by Scripps • The MMTSB aims to develop new tools and theoretical models to aid molecular and structural biologists in interpreting their biological data. – Principal scientists include • Charles Brooks, David Case, Steve Harvey, Jack Johnson, Vijay Reddy, Jeff Skolnick MMTSB/CTBP Summer Workshop © Charles L. Brooks III, 2006. CTBP and MMTSB: Who are we? http://mmtsb.scripps.edu MMTSB/CTBP Summer Workshop © Charles L. Brooks III, 2006. CTBP and MMTSB: Who are we? • Activities – Fundamental research across a broad spectrum – Software and methods development and distribution • MMTSB distributes multiple software packages as well as hosts a variety of web services and databases – Training and research workshops and educational outreach • Both centers have extensive workshop programs – Visitors • Both centers host visitors and collaborators for short and longer term (sabbatical) visits MMTSB/CTBP Summer Workshop © Charles L.
    [Show full text]
  • Trends in Atomistic Simulation Software Usage [1.3]
    A LiveCoMS Perpetual Review Trends in atomistic simulation software usage [1.3] Leopold Talirz1,2,3*, Luca M. Ghiringhelli4, Berend Smit1,3 1Laboratory of Molecular Simulation (LSMO), Institut des Sciences et Ingenierie Chimiques, Valais, École Polytechnique Fédérale de Lausanne, CH-1951 Sion, Switzerland; 2Theory and Simulation of Materials (THEOS), Faculté des Sciences et Techniques de l’Ingénieur, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland; 3National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland; 4The NOMAD Laboratory at the Fritz Haber Institute of the Max Planck Society and Humboldt University, Berlin, Germany This LiveCoMS document is Abstract Driven by the unprecedented computational power available to scientific research, the maintained online on GitHub at https: use of computers in solid-state physics, chemistry and materials science has been on a continuous //github.com/ltalirz/ rise. This review focuses on the software used for the simulation of matter at the atomic scale. We livecoms-atomistic-software; provide a comprehensive overview of major codes in the field, and analyze how citations to these to provide feedback, suggestions, or help codes in the academic literature have evolved since 2010. An interactive version of the underlying improve it, please visit the data set is available at https://atomistic.software. GitHub repository and participate via the issue tracker. This version dated August *For correspondence: 30, 2021 [email protected] (LT) 1 Introduction Gaussian [2], were already released in the 1970s, followed Scientists today have unprecedented access to computa- by force-field codes, such as GROMOS [3], and periodic tional power.
    [Show full text]
  • Application of Various Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens
    International Journal of Molecular Sciences Review Application of Various Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens Anna Helena Mazurek 1 , Łukasz Szeleszczuk 1,* , Thomas Simonson 2 and Dariusz Maciej Pisklak 1 1 Chair and Department of Physical Pharmacy and Bioanalysis, Department of Physical Chemistry, Medical Faculty of Pharmacy, University of Warsaw, Banacha 1 str., 02-093 Warsaw Poland; [email protected] (A.H.M.); [email protected] (D.M.P.) 2 Laboratoire de Biochimie (CNRS UMR7654), Ecole Polytechnique, 91-120 Palaiseau, France; [email protected] * Correspondence: [email protected]; Tel.: +48-501-255-121 Received: 21 July 2020; Accepted: 1 September 2020; Published: 3 September 2020 Abstract: In this review, applications of various molecular modelling methods in the study of estrogens and xenoestrogens are summarized. Selected biomolecules that are the most commonly chosen as molecular modelling objects in this field are presented. In most of the reviewed works, ligand docking using solely force field methods was performed, employing various molecular targets involved in metabolism and action of estrogens. Other molecular modelling methods such as molecular dynamics and combined quantum mechanics with molecular mechanics have also been successfully used to predict the properties of estrogens and xenoestrogens. Among published works, a great number also focused on the application of different types of quantitative structure–activity relationship (QSAR) analyses to examine estrogen’s structures and activities. Although the interactions between estrogens and xenoestrogens with various proteins are the most commonly studied, other aspects such as penetration of estrogens through lipid bilayers or their ability to adsorb on different materials are also explored using theoretical calculations.
    [Show full text]
  • Bringing Open Source to Drug Discovery
    Bringing Open Source to Drug Discovery Chris Swain Cambridge MedChem Consulting Standing on the shoulders of giants • There are a huge number of people involved in writing open source software • It is impossible to acknowledge them all individually • The slide deck will be available for download and includes 25 slides of details and download links – Copy on my website www.cambridgemedchemconsulting.com Why us Open Source software? • Allows access to source code – You can customise the code to suit your needs – If developer ceases trading the code can continue to be developed – Outside scrutiny improves stability and security What Resources are available • Toolkits • Databases • Web Services • Workflows • Applications • Scripts Toolkits • OpenBabel (htttp://openbabel.org) is a chemical toolbox – Ready-to-use programs, and complete programmer's toolkit – Read, write and convert over 110 chemical file formats – Filter and search molecular files using SMARTS and other methods, KNIME add-on – Supports molecular modeling, cheminformatics, bioinformatics – Organic chemistry, inorganic chemistry, solid-state materials, nuclear chemistry – Written in C++ but accessible from Python, Ruby, Perl, Shell scripts… Toolkits • OpenBabel • R • CDK • OpenCL • RDkit • SciPy • Indigo • NumPy • ChemmineR • Pandas • Helium • Flot • FROWNS • GNU Octave • Perlmol • OpenMPI Toolkits • RDKit (http://www.rdkit.org) – A collection of cheminformatics and machine-learning software written in C++ and Python. – Knime nodes – The core algorithms and data structures are written in C ++. Wrappers are provided to use the toolkit from either Python or Java. – Additionally, the RDKit distribution includes a PostgreSQL-based cartridge that allows molecules to be stored in relational database and retrieved via substructure and similarity searches.
    [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]
  • GROMACS Tutorial Lysozyme in Water
    Free Energy Calculations: Methane in Water 1/20 GROMACS Tutorial Lysozyme in water Based on the tutorial created by Justin A. Lemkul, Ph.D. Department of Pharmaceutical Sciences University of Maryland, Baltimore Adapted by Atte Sillanpää, CSC - IT Center for Science Ltd. This example will guide a new user through the process of setting up a simulation system containing a protein (lysozyme) in a box of water, with ions. Each step will contain an explanation of input and output, using typical settings for general use. This tutorial assumes you are using a GROMACS version in the 2018 series. Step One: Prepare the Topology Some GROMACS Basics With the release of version 5.0 of GROMACS, all of the tools are essentially modules of a binary named "gmx" This is a departure from previous versions, wherein each of the tools was invoked as its own command. To get help information about any GROMACS module, you can invoke either of the following commands: Free Energy Calculations: Methane in Water 2/20 gmx help (module) or gmx (module) -h where (module) is replaced by the actual name of the command you're trying to issue. Information will be printed to the terminal, including available algorithms, options, required file formats, known bugs and limitations, etc. For new users of GROMACS, invoking the help information for common commands is a great way to learn about what each command can do. Now, on to the fun stuff! Lysozyme Tutorial We must download the protein structure file with which we will be working. For this tutorial, we will utilize hen egg white lysozyme (PDB code 1AKI).
    [Show full text]