DOCSLIB.ORG

Desmond User's Guide

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Desmond User's Guide Desmond User’s Guide Desmond Version 2.0 / Document Version 0.1 June 2008 D. E. Shaw Research Notice The Desmond User’s Guide and the information it contains is offered solely for educa‐ tional purposes, as a service to users. It is subject to change without notice, as is the software it describes. D. E. Shaw Research assumes no responsibility or liability regarding the correctness or completeness of the information provided herein, nor for damages or loss suffered as a result of actions taken in accordance with said informa‐ tion. No part of this guide may be reproduced, displayed, transmitted, or otherwise copied in any form without written authorization from D. E. Shaw Research. The software described in this guide is copyrighted and licensed by D. E. Shaw Research under separate agreement. This software may be used only according to the terms and conditions of such agreement. Copyright © 2008 by D. E. Shaw Research. All rights reserved. Trademarks Ethernet is a trademark of Xerox Corporation. InfiniBand is a registered trademark of systemI/O Inc. Intel and Pentium are trademarks of Intel Corporation in the U.S. and other countries. Linux is the registered trademark of Linus Torvalds in the U.S. and other countries. All other trademarks are the property of their respective owners. June 2008 Preface Intended audience This guide is intended for computational scientists using Desmond to prepare configura‐ tion and structure files for molecular dynamics simulations. It assumes a broad familiarity with the concepts and techniques of molecular dynamics simulation. Prerequisites Desmond runs on Intel™‐based Linux® systems with Pentium™ 4 or more recent proces‐ sors; running CentOS 3 (RHEL3) or later. Linux clusters can be networked with either Ethernet™ or Infiniband®. To build the source code, Desmond requires gcc Version 4.1.2 and glibc Version 2.3. Viparr requires a recent version of Python; we recommend Version 2.5.1. This guide assumes that someone has prepared Desmond executables for you, either by installing a binary release or by building an executable. About this guide This manual contains the following sections: Key Concepts explains the basic ideas underlying Desmond and describes how Desmond fits into a workflow. June 2008 D. E. Shaw Research i Desmond User’s Guide Running explains the basics of working with configuration files; describes how to invoke the vari‐ Desmond ous Desmond executables, including in parallel; and describes how to configure the Des‐ mond executables and built‐in plugins, as well as the optional profiling mechanism. The Global Cell discusses Desmond’s parallelization strategies and describes how to configure the global cell. Preparing a describes how Desmond interprets structure file format, and explains how to use the structure file scripts provided to add force field and constraint information to the structure file. Calculating provides a high‐level overview of configuring force fields; then discusses the computa‐ force and tions involved in, and how to configure, the various interactions. It also describes addi‐ energy tional off‐atom interaction sites. Constraints describes the constraints available to eliminate the fastest vibrational motions and how to configure them. Dynamics summarizes the basic concepts of particle dynamics and describes how to configure the migration interval, timestep scheduling, pressure, and temperature. It also describes each of the available dynamical systems, and how to configure them. Free Energy explains the concepts necessary to configure ligand‐binding and alchemical free‐energy Simulations simulations, as well as those using the Bennett acceptance ratio method, as well as describing how to prepare a structure file for free energy simulations. Extending provides practical details for implementing a plugin, describes the C++ base classes Desmond available for you to subclass, and explains how to load your plugin, what happens when it begins execution, and what happens when it stops. Units explains how numbers provided as configuration parameters are interpreted. Configuration describes configuration file syntax. syntax Clone Radius provides the full set of restrictions on the size of the clone radius, for those who need Restrictions more than the practical guidelines given in Chapter 3. References provides bibliographic information for references given in the text. License explains the licensing restrictions for using Desmond and its companion software. Information Format conventions Command lines appear in a monospaced font; in some cases, bolding draws your atten‐ tion to a particular part of the command: mdsim --include equil.cfg Placeholders intended to be replaced by actual values are obliqued: mdsim --tpp 4 --restore checkpoint_file Configuration file examples also appear in a monospaced font: mdsim = { title = w last_time = t1 checkpt = { ... } plugins = [ ... ] ii D. E. Shaw Research June 2008 Format conventions ... } Configuration files are divided into sections, which can in turn contain other sections; parameters occur at all levels. When discussed in the context of their particular section, configuration parameters appear by name in a monospace font, thus: plugins. When discussed outside of the context of their sections, however, configuration parame‐ ters appear as a keypath, in which the name of each enclosing section appears in order from outermost to innermost, separated by periods. For example, force.non- bonded.far.sigma refers to the sigma configuration parameter in the far subsection of the nonbonded subsection of the force section of the configuration file. About the equations The equations in this document are concerned with scalars, vectors, and matrices of vari‐ ous sorts. To help clarify the type of a quantity, equations in this manual use the following conventions: •An upper‐ or lowercase letter without bolding or arrows, such as A or a , is a scalar. •An arrow over a variable, such as A or a , indicates three variables as a three‐dimen‐ sional vector. •A boldfaced lowercase letter, such as a ,is a vector of unspecified dimension, th withai indicating the i element of the vector. •A boldfaced uppercase letter is a matrix of unspecified dimensions, though usually 33× , with Aij being the element of row i and column j in matrix A . Certain quantities that are ‐3n dimensional vectors, such as r , the positions of n particles, are indexed differently. The manual does not use ri to refer to one of its 3n components, th but instead ri denotes the i three‐dimensional vector in r , which is the position of the th i particle in this case. June 2008 D. E. Shaw Research iii Desmond User’s Guide iv D. E. Shaw Research June 2008 Contents Preface i Intended audience . i Prerequisites . i About this guide . i Format conventions. ii About the equations. iii 1. Key Concepts. 1 What is Desmond? . 1 Forces . 2 Particles. 3 Force fields . 4 Space . 4 Time. 5 Dynamics. 5 Using Desmond . 6 Input . 7 Executables and scripts . 7 Output . 7 Workflow. 8 Customizing Desmond . 9 June 2008 D. E. Shaw Research v Desmond User Guide 2. Running Desmond . 11 About configuration . 11 Invoking Desmond . 12 Restoring from a checkpoint . 14 Using plugins . 14 Running Desmond in parallel . 15 Configuring Desmond applications . 16 mdsim. 16 minimize . 18 vrun . 19 Naming output files . 19 Configuring the built‐in plugins. 20 compute_forces . 20 energy_groups . 21 eneseq . 21 maeff_output . 22 posre_schedule. 23 randomize_velocities . 23 status . 24 trajectory . 24 Configuring optional sections . 25 profile . 25 3. The Global Cell. 27 Parallelization . 27 Configuration . 30 4. Preparing a structure file. 33 Format . 33 Global cell . 35 Particles and pseudoparticles . 35 Force field sections . 36 Adding force field information . 37 Available force fields . 37 Input and output files . 38 Running Viparr . 39 Options to Viparr . 39 Specifying multiple force fields . 40 User‐defined force fields . 41 Adding constraints . 41 Running the build_constraints script . 42 5. Calculating force and energy . 43 Configuring force fields . 43 Bonded, pair, and excluded interactions . 45 vi D. E. Shaw Research June 2008 Stretch terms . 45 Angle terms . 46 Proper and improper dihedral terms . 46 CMAP torsion terms . ..
Load more

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]
  • 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]
  • 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]
  • Maestro 10.2 User Manual
    Maestro User Manual Maestro 10.2 User Manual Schrödinger Press Maestro User Manual Copyright © 2015 Schrödinger, LLC. All rights reserved. While care has been taken in the preparation of this publication, Schrödinger assumes no responsibility for errors or omissions, or for damages resulting from the use of the information contained herein. Canvas, CombiGlide, ConfGen, Epik, Glide, Impact, Jaguar, Liaison, LigPrep, Maestro, Phase, Prime, PrimeX, QikProp, QikFit, QikSim, QSite, SiteMap, Strike, and WaterMap are trademarks of Schrödinger, LLC. Schrödinger, BioLuminate, and MacroModel are registered trademarks of Schrödinger, LLC. MCPRO is a trademark of William L. Jorgensen. DESMOND is a trademark of D. E. Shaw Research, LLC. Desmond is used with the permission of D. E. Shaw Research. All rights reserved. This publication may contain the trademarks of other companies. Schrödinger software includes software and libraries provided by third parties. For details of the copyrights, and terms and conditions associated with such included third party software, use your browser to open third_party_legal.html, which is in the docs folder of your Schrödinger software installation. This publication may refer to other third party software not included in or with Schrödinger software ("such other third party software"), and provide links to third party Web sites ("linked sites"). References to such other third party software or linked sites do not constitute an endorsement by Schrödinger, LLC or its affiliates. Use of such other third party software and linked sites may be subject to third party license agreements and fees. Schrödinger, LLC and its affiliates have no responsibility or liability, directly or indirectly, for such other third party software and linked sites, or for damage resulting from the use thereof.
    [Show full text]
  • Kepler Gpus and NVIDIA's Life and Material Science
    LIFE AND MATERIAL SCIENCES Mark Berger; [email protected] Founded 1993 Invented GPU 1999 – Computer Graphics Visual Computing, Supercomputing, Cloud & Mobile Computing NVIDIA - Core Technologies and Brands GPU Mobile Cloud ® ® GeForce Tegra GRID Quadro® , Tesla® Accelerated Computing Multi-core plus Many-cores GPU Accelerator CPU Optimized for Many Optimized for Parallel Tasks Serial Tasks 3-10X+ Comp Thruput 7X Memory Bandwidth 5x Energy Efficiency How GPU Acceleration Works Application Code Compute-Intensive Functions Rest of Sequential 5% of Code CPU Code GPU CPU + GPUs : Two Year Heart Beat 32 Volta Stacked DRAM 16 Maxwell Unified Virtual Memory 8 Kepler Dynamic Parallelism 4 Fermi 2 FP64 DP GFLOPS GFLOPS per DP Watt 1 Tesla 0.5 CUDA 2008 2010 2012 2014 Kepler Features Make GPU Coding Easier Hyper-Q Dynamic Parallelism Speedup Legacy MPI Apps Less Back-Forth, Simpler Code FERMI 1 Work Queue CPU Fermi GPU CPU Kepler GPU KEPLER 32 Concurrent Work Queues Developer Momentum Continues to Grow 100M 430M CUDA –Capable GPUs CUDA-Capable GPUs 150K 1.6M CUDA Downloads CUDA Downloads 1 50 Supercomputer Supercomputers 60 640 University Courses University Courses 4,000 37,000 Academic Papers Academic Papers 2008 2013 Explosive Growth of GPU Accelerated Apps # of Apps Top Scientific Apps 200 61% Increase Molecular AMBER LAMMPS CHARMM NAMD Dynamics GROMACS DL_POLY 150 Quantum QMCPACK Gaussian 40% Increase Quantum Espresso NWChem Chemistry GAMESS-US VASP CAM-SE 100 Climate & COSMO NIM GEOS-5 Weather WRF Chroma GTS 50 Physics Denovo ENZO GTC MILC ANSYS Mechanical ANSYS Fluent 0 CAE MSC Nastran OpenFOAM 2010 2011 2012 SIMULIA Abaqus LS-DYNA Accelerated, In Development NVIDIA GPU Life Science Focus Molecular Dynamics: All codes are available AMBER, CHARMM, DESMOND, DL_POLY, GROMACS, LAMMPS, NAMD Great multi-GPU performance GPU codes: ACEMD, HOOMD-Blue Focus: scaling to large numbers of GPUs Quantum Chemistry: key codes ported or optimizing Active GPU acceleration projects: VASP, NWChem, Gaussian, GAMESS, ABINIT, Quantum Espresso, BigDFT, CP2K, GPAW, etc.
    [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]
  • In Silico Strategy for Targeting the Mtor Kinase at Rapamycin Binding Site by Small Molecules
    molecules Article In Silico Strategy for Targeting the mTOR Kinase at Rapamycin Binding Site by Small Molecules Serena Vittorio 1 , Rosaria Gitto 1 , Ilenia Adornato 1, Emilio Russo 2 and Laura De Luca 1,* 1 Department of Chemical, Biological, Pharmaceutical and Environmental Sciences, University of Messina, Viale Palatucci 13, I-98168 Messina, Italy; [email protected] (S.V.); [email protected] (R.G.); [email protected] (I.A.) 2 Science of Health Department, School of Medicine, University “Magna Graecia” of Catanzaro, Viale Europa e Germaneto, I-88100 Catanzaro, Italy; [email protected] * Correspondence: [email protected] or [email protected]; Tel.: +39-090-676-6410 Abstract: Computer aided drug-design methods proved to be powerful tools for the identification of new therapeutic agents. We employed a structure-based workflow to identify new inhibitors targeting mTOR kinase at rapamycin binding site. By combining molecular dynamics (MD) simulation and pharmacophore modelling, a simplified structure-based pharmacophore hypothesis was built starting from the FKBP12-rapamycin-FRB ternary complex retrieved from RCSB Protein Data Bank (PDB code 1FAP). Then, the obtained model was used as filter to screen the ZINC biogenic compounds library, containing molecules derived from natural sources or natural-inspired compounds. The resulting hits were clustered according to their similarity; moreover, compounds showing the highest pharmacophore fit-score were chosen from each cluster. The selected molecules were subjected to docking studies to clarify their putative binding mode. The binding free energy of the obtained complexes was calculated by MM/GBSA method and the hits characterized by the lowest DGbind Citation: Vittorio, S.; Gitto, R.; values were identified as potential mTOR inhibitors.
    [Show full text]
  • Lawrence Berkeley National Laboratory Recent Work
    Lawrence Berkeley National Laboratory Recent Work Title From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape. Permalink https://escholarship.org/uc/item/4sm897jh Journal Chemical reviews, 121(8) ISSN 0009-2665 Authors Kowalski, Karol Bair, Raymond Bauman, Nicholas P et al. Publication Date 2021-04-01 DOI 10.1021/acs.chemrev.0c00998 Peer reviewed eScholarship.org Powered by the California Digital Library University of California From NWChem to NWChemEx: Evolving with the computational chemistry landscape Karol Kowalski,y Raymond Bair,z Nicholas P. Bauman,y Jeffery S. Boschen,{ Eric J. Bylaska,y Jeff Daily,y Wibe A. de Jong,x Thom Dunning, Jr,y Niranjan Govind,y Robert J. Harrison,k Murat Keçeli,z Kristopher Keipert,? Sriram Krishnamoorthy,y Suraj Kumar,y Erdal Mutlu,y Bruce Palmer,y Ajay Panyala,y Bo Peng,y Ryan M. Richard,{ T. P. Straatsma,# Peter Sushko,y Edward F. Valeev,@ Marat Valiev,y Hubertus J. J. van Dam,4 Jonathan M. Waldrop,{ David B. Williams-Young,x Chao Yang,x Marcin Zalewski,y and Theresa L. Windus*,r yPacific Northwest National Laboratory, Richland, WA 99352 zArgonne National Laboratory, Lemont, IL 60439 {Ames Laboratory, Ames, IA 50011 xLawrence Berkeley National Laboratory, Berkeley, 94720 kInstitute for Advanced Computational Science, Stony Brook University, Stony Brook, NY 11794 ?NVIDIA Inc, previously Argonne National Laboratory, Lemont, IL 60439 #National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37831-6373 @Department of Chemistry, Virginia Tech, Blacksburg, VA 24061 4Brookhaven National Laboratory, Upton, NY 11973 rDepartment of Chemistry, Iowa State University and Ames Laboratory, Ames, IA 50011 E-mail: [email protected] 1 Abstract Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems.
    [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]
  • Elucidating Binding Sites and Affinities of Erα Agonist and Antagonist to Human Alpha-Fetoprotein by in Silico Modeling and Point Mutagenesis Nurbubu T
    Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 29 December 2019 doi:10.20944/preprints201912.0374.v1 Peer-reviewed version available at Int. J. Mol. Sci. 2020, 21, 893; doi:10.3390/ijms21030893 Elucidating Binding Sites and Affinities of ERα Agonist and Antagonist to Human Alpha-Fetoprotein by in silico Modeling and Point Mutagenesis Nurbubu T. Moldogazieva1*, Daria S. Ostroverkhova1,2, Nikolai N. Kuzmich1,3, Vladimir V. Kadochnikov4, Alexander A. Terentiev5, Yuri B. Porozov1,4 1 Laboratory of Bioinformatics, I.M. Sechenov First Moscow State Medical University (Sechenov University); Moscow, Russia; 2Department of Bioengineering, M.V. Lomonosov Moscow State University; Moscow, Russia; 3Department of Drug Safety, I.M. Smorodintsev Research Institute of Influenza, WHO National Influenza Centre of Russia; Saint Petersburg, Russia 4Department of Information and Communication Technologies, Saint Petersburg National Research University of Information Technologies, Mechanics and Optics; Saint Petersburg, Russia 5Deparment of Biochemistry and Molecular Biology, N.I. Pirogov Russian National Research Medical University; Moscow Russia *Corresponding author: Nurbubu T. Moldogazieva, e-mail: [email protected] (NM) ¶ These authors contributed equally to this work Short title: Binding of ERα ligands to human alpha-fetoprotein Key words: alpha-fetoprotein, estrogens, selective estrogen receptor modulators, homology-based modeling, molecular docking, protein-ligand interaction, amino acid substitutions. Abbreviations: HAFP, human alpha-fetoprotein; HSA, human serum albumin; E2, 17β-estradiol; DES, diethylstilbestrol. Author contribution: Conceptualization, Nurbubu T. Moldogazieva, Alexander A. Terentiev and Yuri B. Porozov; Data curation, Nurbubu T. Moldogazieva; Formal analysis, Nurbubu T. Moldogazieva and Nikolai N. Kuzmich; Investigation, Daria S. Ostroverkhova and Vladimir V. Kadochnikov; Methodology, Nikolai N. Kuzmich; Supervision, Alexander A.
    [Show full text]