In Silico Strategy for Targeting the Mtor Kinase at Rapamycin Binding Site by Small Molecules

Total Page:16

File Type:pdf, Size:1020Kb

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. Furthermore, the stability of the resulting Adornato, I.; Russo, E.; De Luca, L. In complexes was studied by means of MD simulation which revealed that the selected compounds Silico Strategy for Targeting the were able to form a stable ternary complex with FKBP12 and FRB domain, thus underlining their mTOR Kinase at Rapamycin Binding potential ability to inhibit mTOR with a rapamycin-like mechanism. Site by Small Molecules. Molecules 2021, 26, 1103. https://doi.org/ Keywords: mTOR; molecular dynamics simulation; ternary complex FKBP12-rapamycin-FRB; molec- 10.3390/molecules26041103 ular docking and structure-based virtual screening; pharmacophore modeling Academic Editor: Rita Guedes Received: 23 December 2020 1. Introduction Accepted: 17 February 2021 Published: 19 February 2021 The mammalian target of rapamycin (mTOR) protein is a highly conserved eukary- otic Serine/Threonine (Ser/Thr) protein kinase, which belongs to phosphatidylinositol Publisher’s Note: MDPI stays neutral 3-kinase-related kinase (PI3K) family. mTOR plays a central role in various cellular pro- with regard to jurisdictional claims in cesses including protein biosynthesis by controlling transcription, translation and ribosomal published maps and institutional affil- biosynthesis thus regulating cell growth and size. Since the mTOR inhibition causes cell iations. cycle arrest at G1 phase, the intervention on mTOR pathway might provide a significant direction to new drug development aimed to the identification of anticancer, antiprolifera- tive, anti-restenosis and immunosuppressor agents. In addition, recent evidence suggests that an altered mTOR pathway is implicated in several neurological disorders [1–3]. Copyright: © 2021 by the authors. The mTOR is composed of 2549 amino acidic residues (289 kDa) and organized into Licensee MDPI, Basel, Switzerland. five distinct functional domains, (i) 32 replicates called “Huntington, elongation factor This article is an open access article 3,PR65/A, TOR” (HEAT) (from amino acid residue 16 to 1345); (ii) “FRAP, ATM, TRRAP” distributed under the terms and (FAT) (from amino acid residue 1382 to1982); (iii) “FKBP12 rapamycin-binding” region conditions of the Creative Commons (FRB) (from aminoacid residue 2012 to 2144); (iv) catalytic “kinase domain”(KD) and Attribution (CC BY) license (https:// “regulatory domain” (RD) (from aminoacid residue 2182 to 2516); and (v) “FRAP ATM creativecommons.org/licenses/by/ TRRAP carboxyterminus” (FATC) (from aminoacid residue 2517 to 2549). 4.0/). Molecules 2021, 26, 1103. https://doi.org/10.3390/molecules26041103 https://www.mdpi.com/journal/molecules Molecules 2021, 26, x FOR PEER REVIEW 2 of 14 Molecules 2021, 26, 1103 (FRB) (from aminoacid residue 2012 to 2144); (iv) catalytic “kinase domain”(KD) and “reg-2 of 13 ulatory domain” (RD) (from aminoacid residue 2182 to 2516); and (v) “FRAP ATM TRRAP carboxyterminus” (FATC) (from aminoacid residue 2517 to 2549) FromFrom a functional a functional point point of of view, view, mTOR existsexists inin two two distinct distinct multi-protein multi-protein complexes com- plexescommonly commonly known known as mTOR as complexmTOR 1complex (mTORC1) 1 and(mTORC1) mTOR complex and mTOR 2 (mTORC2) complex located 2 (mTORC2)in different located cellular in compartments.different cellular It has comp beenartments. well established It has been that well the rapamycin-sensitive established that themTORC1 rapamycin-sensitive complex transmits mTORC1 signals complex of nutrient transmits availability signals of to nutrient control availability various cellular to controlfunctions, various while cellular the rapamycin-insensitivefunctions, while the rapamycin-insensitive mTORC2 complex regulatesmTORC2 thecomplex assembly reg- of ulatesthe cytoskeletonthe assembly byof actin.the cytoskeleton by actin. It isIt well-known is well-known that that for for the the microbial-derived microbial-derived macrolide macrolide rapamycin rapamycin (sirolimus) (sirolimus) and and its itsanalogues analogues (Scheme (Scheme 1) 1the) the inhibition inhibition mechanism mechanism of ofmTOR mTOR pathway pathway occurs occurs through through the the bindingbinding with with FKBP12 FKBP12 domain domain to to form form a a binary binary complex; inin turn,turn, this this complex complex binds binds to to FRB FRBdomain domain of of mTOR mTOR only only when when is associatedis associated to otherto other mTORC1 mTORC1 components; components; this this association asso- ciationinduces induces changes changes in the in mTOR the mTOR conformation, conformation, thus thus inhibiting inhibiting the crucial the crucial kinase kinase activity ac- of tivitythe of mTORC1 the mTORC1 [4,5]. Experimental [4,5]. Experimental studies studies highlighted highlighted that in absencethat in absence of FKBP12, of FKBP12, rapamycin rapamycinpossesses possesses a modest a affinitymodestfor affinity the FRB for domain,the FRB domain, suggesting suggesting that the FKBP12-rapamycinthat the FKBP12- rapamycincomplex complex plays a pivotal plays a role pivotal for the role binding for the tobinding mTOR to [6 mTOR]. [6]. Scheme 1. Chemical structures of Rapamycin and rapalogs as well-established mTOR/FKBP12 ligands. Scheme 1. Chemical structures of Rapamycin and rapalogs as well-established mTOR/FKBP12 ligands. Rapamycin is 31-membered macrocyclic lactone produced by Streptomyces hygro- scopicusRapamycinwhich is was 31-membered developed asmacr anocyclic immunosuppressant lactone produced agent by as Streptomyces allosteric inhibitor hygro- of scopicusmTORC1. which The was crystal developed structure as an of theimmunosuppressant FKBP12–rapamycin–mTOR agent as allosteric ternary complexinhibitor(PDB of mTORC1.code 1FAP) The crystal unveiled structure the protein of the interactions. FKBP12–rapamycin–mTOR It has been found ternary that complex the pipecolyl (PDB α- codeketoamide 1FAP) unveiled of rapamycin the protein anchored interactions it into the. It proline-bindinghas been found pocket,that thewhereas pipecolyl the α-ke- triene toamidesystem of was rapamycin exposed anchored for interactions it into withthe pr mTOR.oline-binding Rapamycin pocket, displays whereas low the water-solubility triene sys- temand was poor exposed stability, for sointeractions that rapamycin with analoguesmTOR. Rapamycin (also named displays rapalogs) low withwater-solubility improved bio- andpharmaceutical poor stability, properties so that rapamycin have been analogues developed [(also7,8] andnamed approved rapalogs) by FDA with (see improved Scheme 1) biopharmaceuticalas the first-generation properties of mTOR have inhibitors been developed to fight cancer [7,8] malignanciesand approved and by other FDA diseases. (see SchemeApart 1) from as the the first-generation weakness in poor of mTOR druglike inhibitors properties, to fight the rapalogscancer malignancies possess a complex and otherchemical diseases. structure Apart [from5]; therefore, the weakness the structural in poor druglike modifications properties of macrolide, the rapalogs ring were possess gener- a complexally limited. chemical structure [5]; therefore, the structural modifications of macrolide ring were generallyFurther limited. allosteric mTOR inhibitors belonging to rapalog series are modified at C- 7,Further C-22, C-27 allosteric and C-42 mTOR positions inhibitors as well belonging as the to C-1/C4 rapalog fragment. series are A modified carefully at analysis C-7, C-22,of structure-activityC-27 and C-42 positions relationships as well as of the rapalogs C-1/C4 hasfragment. been recentlyA carefully reported analysis [5 of]; struc- the best ture-activityresults were relationships obtained for of structural rapalogs has optimization been recently carried reported out addressing
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]
  • 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]
  • 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]
  • 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]
  • 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]
  • Pelvic Fin Flexibility in Tree Climbing Fish
    G Model ZOOL-25524; No. of Pages 7 ARTICLE IN PRESS Zoology xxx (2016) xxx–xxx Contents lists available at ScienceDirect Zoology journal homepage: www.elsevier.com/locate/zool The significance of pelvic fin flexibility for tree climbing fish a b a,b c Adhityo Wicaksono , Saifullah Hidayat , Yudithia Damayanti , Desmond Soo Mun Jin , a b,∗ a,∗ Erly Sintya , Bambang Retnoaji , Parvez Alam a Laboratory of Paper Coating and Converting, Centre for Functional Materials, Abo Akademi University, Porthaninkatu 3, 20500 Turku, Finland b Laboratory of Animal Embryology, Faculty of Biology, Universitas Gadjah Mada, Yogyakarta, Indonesia c Rapid Gain Global Corporation, Singapore a r t i c l e i n f o a b s t r a c t Article history: In this article, we compare the characteristics of biomechanical attachment exhibited by two morpholog- Received 1 March 2016 ically different mudskipper species, Boleophthalmus boddarti (with fused pelvic fins) and Periophthalmus Received in revised form 14 April 2016 variabilis (with unfused pelvic fins). P. variabilis is a tree and rock climber while B. boddarti dwells in the Accepted 17 June 2016 muddy shallows and is unable to climb. Our aim in this article is to determine whether it is predominantly Available online xxx chemical or morphological properties of the pelvic fins from each species that may allow P. variabilis to climb trees whilst preventing B. boddarti from doing the same. To fulfil our objective we perform friction Keywords: and suction resistance tests, Fourier transform infrared spectroscopy of the mucosal secretions under Mudskipper the fins, direct geometrical measurements and finite element modelling.
    [Show full text]
  • Chemistry Packages at CHPC
    Chemistry Packages at CHPC Anita M. Orendt Center for High Performance Computing [email protected] Fall 2014 Purpose of Presentation • Identify the computational chemistry software and related tools currently available at CHPC • Present brief overview of these packages • Present how to access packages on CHPC • Information on usage of Gaussian09 http://www.chpc.utah.edu Ember Lonepeak /scratch/lonepeak/serial 141 nodes/1692 cores 16 nodes/256 cores Infiniband and GigE Large memory, GigE General 70 nodes/852 cores 12 GPU nodes (6 General) Kingspeak 199 nodes/4008 cores Ash cluster (417 nodes) Infiniband and GigE Switch General 44 nodes/832 cores Apex Arch Great Wall PFS Telluride cluster (93 nodes) /scratch/ibrix meteoXX/atmosXX nodes NFS Home NFS Administrative Directories & /scratch/kingspeak/serial Group Nodes Directories 10/1/14 http://www.chpc.utah.edu Slide 3 Brief Overview CHPC Resources • Computational Clusters – kingspeak, ember – general allocation and owner-guest on owner nodes – lonepeak – no allocation – ash and telluride – can run as smithp-guest and cheatham-guest, respectively • Home directory – NFS mounted on all clusters – /uufs/chpc.utah.edu/common/home/<uNID> – generally not backed up (there are exceptions) • Scratch systems – /scratch/ibrix/chpc_gen – kingspeak, ember, ash, telluride – 56TB – /scratch/kingspeak/serial – kingspeak, ember, ash – 175TB – /scratch/lonepeak/serial – lonepeak – 33 TB – /scratch/local on compute nodes • Applications – /uufs/chpc.utah.edu/sys/pkg (use for lonepeak, linux desktops if chpc
    [Show full text]
  • Download Article (PDF)
    International Journal of Networked and Distributed Computing, Vol. 5, No. 1 (January 2017) 52–61 ___________________________________________________________________________________________________________ An FPGA Accelerator for Molecular Dynamics Simulation Using OpenCL Hasitha Muthumala Waidyasooriya 1, Masanori Hariyama 1, Kota Kasahara 2 1 Graduate School of Information Sciences, Tohoku University, Aoba 6-6-05, Aramaki, Aoba, Sendai, Miyagi 980-8579, Japan E-mail: {hasitha, hariyama}@tohoku.ac.jp 2 College of Life Sciences, Ritsumeikan University, 1-1-1, Noji-higashi, Kusatsu, Shiga 525-8577, Japan E-mail: [email protected] Abstract Molecular dynamics (MD) simulations are very important to study physical properties of the atoms and molecules. However, a huge amount of processing time is required to simulate a few nano-seconds of an actual experiment. Although the hardware acceleration using FPGAs provides promising results, huge design time and hardware design skills are required to implement an accelerator successfully. In this paper, we use a heterogeneous computing system for MD simulations, that can be used in C-based programming environment. We propose an FPGA accelerator designed using C-based OpenCL for the heterogeneous environment. We achieved over 4.6 times of speed-up compared to CPU-based processing, by using only 36% of the Stratix V FPGA resources. We also evaluate the processing times of different tasks in the heterogeneous environment. Keywords: OpenCL for FPGA, molecular dynamics simulation, hardware acceleration, scientific com- puting. 1. Introduction tation results of one iteration are used as the inputs in the next iteration. Each iteration consists of two major phases: force computation and motion update 1 Molecular dynamics (MD) simulations are very as shown in Fig.1.
    [Show full text]
  • Download PDF 137.14 KB
    Why should biochemistry students be introduced to molecular dynamics simulations—and how can we introduce them? Donald E. Elmore Department of Chemistry and Biochemistry Program; Wellesley College; Wellesley, MA 02481 Abstract Molecular dynamics (MD) simulations play an increasingly important role in many aspects of biochemical research but are often not part of the biochemistry curricula at the undergraduate level. This article discusses the pedagogical value of exposing students to MD simulations and provides information to help instructors consider what software and hardware resources are necessary to successfully introduce these simulations into their courses. In addition, a brief review of the MD based activities in this issue and other sources are provided. Keywords molecular dynamics simulations, molecular modeling, biochemistry education, undergraduate education, computational chemistry Published in Biochemistry and Molecular Biology Education, 2016, 44:118-123. Several years ago, I developed a course on the modeling of biochemical systems at Wellesley College, which was described in more detail in a previous issue of this journal [1]. As I set out to design that course, I decided to place a particular emphasis on molecular dynamics (MD) simulations. In full candor, this decision was likely related to my experience integrating these simulations into my own research, but I also believe it was due to some inherent advantages of exposing students to this particular computational approach. Notably, even as this initial biochemical molecular modeling course evolved into a more general (i.e. not solely biochemical) course in computational chemistry, a significant emphasis on MD has remained, emphasizing the central role I and my other colleague who teaches the current course felt it played in students’ education.
    [Show full text]
  • Lessons Learned from Comparing Molecular Dynamics Engines on the SAMPL5 Dataset
    JCAMD manuscript No. (will be inserted by the editor) Lessons learned from comparing molecular dynamics engines on the SAMPL5 dataset Michael R. Shirts · Christoph Klein · Jason M. Swails · Jian Yin · Michael K. Gilson · David L. Mobley · David A. Case · Ellen D. Zhong Received: date / Accepted: date Abstract We describe our efforts to prepare common starting structures and models for the SAMPL5 blind prediction challenge. We generated the starting input files and single configuration potential energies for the host-guest in the SAMPL5 blind prediction challenge for the GROMACS, AMBER, LAMMPS, DESMOND and CHARMM molecular simulation programs. All conversions were fully automated from the originally prepared AMBER input files using a combination of the ParmEd and InterMol conversion programs. We find that the energy calculations for all molecular dynamics engines for this molecular set agree to a better than 0.1% relative absolute energy for all energy components, and in most cases an order of magnitude better, when Michael R. Shirts Department of Chemical and Biological Engineering, University of Colorado Boulder, Boul- der, Colorado, USA E-mail: [email protected] Christoph Klein Department of Chemical Engineering, Vanderbilt University, Nashville, TN, USA Jason M. Swails Department of Chemistry and Chemical Biology, Rutgers University, Rutgers, NJ, USA Jian Yin Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, CA, USA Michael K. Gilson Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, CA, USA David L. Mobley Departments of Pharmaceutical Sciences and Chemistry, University of California, Irvine, Irvine, CA, USA David A. Case Department of Chemistry and Chemical Biology, Rutgers University, Rutgers, NJ, USA Ellen D.
    [Show full text]
  • NVIDIA Cryo-EM May 2019
    Molecular Dynamics (MD) on GPUs March 2019 Accelerating Discoveries Using a supercomputer powered by the Tesla Platform with over 3,000 Tesla accelerators, University of Illinois scientists performed the first all-atom simulation of the HIV virus and discovered the chemical structure of its capsid — “the perfect target for fighting the infection.” Without gpu, the supercomputer would need to be 5x larger for similar performance. 2 Overview of Life & Material Accelerated Apps MD QC All key codes are GPU-accelerated All key codes are ported or optimizing Great multi-GPU, multi-node (dense) performance GPU-accelerated math libraries, OpenACC directives GPU-accelerated apps GPU-accelerated apps ACEMD*, AMBER*, BAND, CHARMM, DESMOND, ESPResso, Folding@Home, ABINIT, ACES III, ADF, BigDFT, CP2K, GAMESS, GAMESS-UK, GPAW, LATTE, LSDalton, LSMS, GPUgrid.net, GROMACS, HALMD, HOOMD-Blue*, LAMMPS, Lattice Microbes*, mdcore, MOLCAS, MOPAC2012, NWChem, OCTOPUS*, PEtot, QUICK, Q-Chem, QMCPack, Quantum MELD, miniMD, NAMD, OpenMM, PolyFTS, SOP-GPU* & more Espresso/PWscf, QUICK, TeraChem* Active acceleration projects CASTEP, GAMESS, Gaussian, ONETEP, Quantum Supercharger Library*, VASP & more green* >90% of the workload is on GPU 3 MD vs. QC on GPUs Molecular Dynamics Quantum Chemistry Properties - electronic properties, Simulates atomic positions over time Calculations ground state, excitation, spectra Chemical-biological or chemical-material Examples: MO, PW, DFT, semi-emp Simple empirical formulas Electron wave function Forces No bond rearrangements
    [Show full text]