1
|
Acharya A, Agarwal R, Baker M, Baudry J, Bhowmik D, Boehm S, Byler KG, Chen S, Coates L, Cooper C, Demerdash O, Daidone I, Eblen J, Ellingson S, Forli S, Glaser J, Gumbart JC, Gunnels J, Hernandez O, Irle S, Kneller D, Kovalevsky A, Larkin J, Lawrence T, LeGrand S, Liu SH, Mitchell J, Park G, Parks J, Pavlova A, Petridis L, Poole D, Pouchard L, Ramanathan A, Rogers D, Santos-Martins D, Scheinberg A, Sedova A, Shen Y, Smith J, Smith M, Soto C, Tsaris A, Thavappiragasam M, Tillack A, Vermaas J, Vuong V, Yin J, Yoo S, Zahran M, Zanetti-Polzi L. Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19. J Chem Inf Model 2020; 60:5832-5852. [PMID: 33326239 PMCID: PMC7754786 DOI: 10.1021/acs.jcim.0c01010] [Citation(s) in RCA: 109] [Impact Index Per Article: 27.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2020] [Indexed: 01/18/2023]
Abstract
We present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.
Collapse
Affiliation(s)
- A. Acharya
- School of Physics, Georgia Institute of Technology, Atlanta, GA 30332, USA
| | - R. Agarwal
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830, USA
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996, USA
- Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN, 37996, USA
| | - M. Baker
- Computer Science and Mathematics Division, Oak Ridge National Lab, Oak Ridge, TN 37830, USA
| | - J. Baudry
- The University of Alabama in Huntsville, Department of Biological Sciences. 301 Sparkman Drive, Huntsville, AL 35899, USA
| | - D. Bhowmik
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
| | - S. Boehm
- Computer Science and Mathematics Division, Oak Ridge National Lab, Oak Ridge, TN 37830, USA
| | - K. G. Byler
- The University of Alabama in Huntsville, Department of Biological Sciences. 301 Sparkman Drive, Huntsville, AL 35899, USA
| | - S.Y. Chen
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - L. Coates
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
| | - C.J. Cooper
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830, USA
- Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN, 37996, USA
| | - O. Demerdash
- Biosciences Division, Oak Ridge National Lab, Oak Ridge, TN 37830, USA
| | - I. Daidone
- Department of Physical and Chemical Sciences, University of L’Aquila, I-67010 L’Aquila, Italy
| | - J.D. Eblen
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830, USA
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996, USA
| | - S. Ellingson
- University of Kentucky, Division of Biomedical Informatics, College of Medicine, UK Medical Center MN 150, Lexington KY, 40536, USA
| | - S. Forli
- Scripps Research, La Jolla, CA, 92037, USA
| | - J. Glaser
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| | - J. C. Gumbart
- School of Physics, Georgia Institute of Technology, Atlanta, GA 30332, USA
| | - J. Gunnels
- HPC Engineering, Amazon Web Services, Seattle, WA 98121, USA
| | - O. Hernandez
- Computer Science and Mathematics Division, Oak Ridge National Lab, Oak Ridge, TN 37830, USA
| | - S. Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
- Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN 37996, USA
| | - D.W. Kneller
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
| | - A. Kovalevsky
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
| | - J. Larkin
- NVIDIA Corporation, Santa Clara, CA 95051, USA
| | - T.J. Lawrence
- Biosciences Division, Oak Ridge National Lab, Oak Ridge, TN 37830, USA
| | - S. LeGrand
- NVIDIA Corporation, Santa Clara, CA 95051, USA
| | - S.-H. Liu
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830, USA
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996, USA
| | - J.C. Mitchell
- Biosciences Division, Oak Ridge National Lab, Oak Ridge, TN 37830, USA
| | - G. Park
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - J.M. Parks
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830, USA
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996, USA
- Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN, 37996, USA
| | - A. Pavlova
- School of Physics, Georgia Institute of Technology, Atlanta, GA 30332, USA
| | - L. Petridis
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830, USA
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996, USA
| | - D. Poole
- NVIDIA Corporation, Santa Clara, CA 95051, USA
| | - L. Pouchard
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - A. Ramanathan
- Data Science and Learning Division, Argonne National Lab, Lemont, IL 60439, USA
| | - D. Rogers
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| | | | | | - A. Sedova
- Biosciences Division, Oak Ridge National Lab, Oak Ridge, TN 37830, USA
| | - Y. Shen
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830, USA
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996, USA
- Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN, 37996, USA
| | - J.C. Smith
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830, USA
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996, USA
| | - M.D. Smith
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830, USA
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996, USA
| | - C. Soto
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - A. Tsaris
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| | | | | | - J.V. Vermaas
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| | - V.Q. Vuong
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
- Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN 37996, USA
| | - J. Yin
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| | - S. Yoo
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - M. Zahran
- Department of Biological Sciences, New York City College of Technology, The City University of New York (CUNY), Brooklyn, NY 11201, USA
| | | |
Collapse
|
2
|
Acharya A, Agarwal R, Baker M, Baudry J, Bhowmik D, Boehm S, Byler KG, Coates L, Chen SY, Cooper CJ, Demerdash O, Daidone I, Eblen JD, Ellingson S, Forli S, Glaser J, Gumbart JC, Gunnels J, Hernandez O, Irle S, Larkin J, Lawrence TJ, LeGrand S, Liu SH, Mitchell JC, Park G, Parks JM, Pavlova A, Petridis L, Poole D, Pouchard L, Ramanathan A, Rogers D, Santos-Martins D, Scheinberg A, Sedova A, Shen S, Smith JC, Smith MD, Soto C, Tsaris A, Thavappiragasam M, Tillack AF, Vermaas JV, Vuong VQ, Yin J, Yoo S, Zahran M, Zanetti-Polzi L. Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19. ChemRxiv 2020:12725465. [PMID: 33200117 PMCID: PMC7668744 DOI: 10.26434/chemrxiv.12725465] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Revised: 07/29/2020] [Indexed: 01/18/2023]
Abstract
We present a supercomputer-driven pipeline for in-silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. We also describe preliminary results obtained for 23 systems involving eight protein targets of the proteome of SARS CoV-2. THe MD performed is temperature replica-exchange enhanced sampling, making use of the massively parallel supercomputing on the SUMMIT supercomputer at Oak Ridge National Laboratory, with which more than 1ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to ten configurations of each of the 23 SARS CoV-2 systems using AutoDock Vina. We also demonstrate that using Autodock-GPU on SUMMIT, it is possible to perform exhaustive docking of one billion compounds in under 24 hours. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and AI methods to cluster MD trajectories and rescore docking poses.
Collapse
Affiliation(s)
- A Acharya
- School of Physics, Georgia Institute of Technology, Atlanta, GA 30332
| | - R Agarwal
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996
- Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN, 37996
| | - M Baker
- Computer Science and Mathematics Division, Oak Ridge National Lab, Oak Ridge, TN 37830
| | - J Baudry
- The University of Alabama in Huntsville, Department of Biological Sciences. 301 Sparkman Drive, Huntsville, AL 35899
| | - D Bhowmik
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831
| | - S Boehm
- Computer Science and Mathematics Division, Oak Ridge National Lab, Oak Ridge, TN 37830
| | - K G Byler
- The University of Alabama in Huntsville, Department of Biological Sciences. 301 Sparkman Drive, Huntsville, AL 35899
| | - L Coates
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831
| | - S Y Chen
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973
| | - C J Cooper
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830
- Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN, 37996
| | - O Demerdash
- Biosciences Division, Oak Ridge National Lab, Oak Ridge, TN 37830
| | - I Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, I-67010 L'Aquila, Italy
| | - J D Eblen
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996
| | - S Ellingson
- University of Kentucky, Division of Biomedical Informatics, College of Medicine, UK Medical Center MN 150, Lexington KY, 40536
| | - S Forli
- Scripps Research, La Jolla, CA, 92037
| | - J Glaser
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830
| | - J C Gumbart
- School of Physics, Georgia Institute of Technology, Atlanta, GA 30332
| | - J Gunnels
- HPC Engineering, Amazon Web Services, Seattle, WA 98121
| | - O Hernandez
- Computer Science and Mathematics Division, Oak Ridge National Lab, Oak Ridge, TN 37830
| | - S Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831
- Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN 37996
| | - J Larkin
- NVIDIA Corporation, Santa Clara, CA 95051
| | - T J Lawrence
- Biosciences Division, Oak Ridge National Lab, Oak Ridge, TN 37830
| | - S LeGrand
- NVIDIA Corporation, Santa Clara, CA 95051
| | - S-H Liu
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996
| | - J C Mitchell
- Biosciences Division, Oak Ridge National Lab, Oak Ridge, TN 37830
| | - G Park
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973
| | - J M Parks
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996
- Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN, 37996
| | - A Pavlova
- School of Physics, Georgia Institute of Technology, Atlanta, GA 30332
| | - L Petridis
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996
| | - D Poole
- NVIDIA Corporation, Santa Clara, CA 95051
| | - L Pouchard
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973
| | - A Ramanathan
- Data Science and Learning Division, Argonne National Lab, Lemont, IL 60439
| | - D Rogers
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830
| | | | | | - A Sedova
- Biosciences Division, Oak Ridge National Lab, Oak Ridge, TN 37830
| | - S Shen
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996
- Graduate School of Genome Science and Technology, University of Tennessee, Knoxville, TN, 37996
| | - J C Smith
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996
| | - M D Smith
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, TN, 37830
- The University of Tennessee, Knoxville. Department of Biochemistry & Cellular and Molecular Biology, 309 Ken and Blaire Mossman Bldg. 1311 Cumberland Avenue Knoxville, TN, 37996
| | - C Soto
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973
| | - A Tsaris
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830
| | | | | | - J V Vermaas
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830
| | - V Q Vuong
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831
- Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN 37996
| | - J Yin
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37830
| | - S Yoo
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973
| | - M Zahran
- Department of Biological Sciences, New York City College of Technology, The City University of New York (CUNY), Brooklyn, NY 11201
| | | |
Collapse
|
3
|
Guzzo T, Mandaliti W, Nepravishta R, Aramini A, Bodo E, Daidone I, Allegretti M, Topai A, Paci M. Conformational Change in the Mechanism of Inclusion of Ketoprofen in β-Cyclodextrin: NMR Spectroscopy, Ab Initio Calculations, Molecular Dynamics Simulations, and Photoreactivity. J Phys Chem B 2016; 120:10668-10678. [PMID: 27689813 DOI: 10.1021/acs.jpcb.6b07913] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Inclusion of drugs in cyclodextrins (CDs) is a recognized tool for modifying several properties such as solubility, stability, bioavailability, and so on. The photoreactive behavior of the β-CD/ketoprofen (KP) complex upon UV exposure showed a significant increase in photodecarboxylation, whereas the secondary degradation products by hydroxylation of the benzophenone moiety were inhibited. The results may account for an improvement of KP photophysical properties upon inclusion, thus better fostering its topical use. To correlate the structural details of the inclusion with these results, an NMR spectroscopic study of KP upon inclusion in β-CD was performed. Effects of the magnetically anisotropic centers of KP, changing their orientations upon inclusion and giving chemical shift variations, were specifically correlated with the results of the molecular dynamic simulations and ab initio calculations. In the large variety of papers focusing on the structural analysis of β-CD complexes, this work represents one of the few examples in which a detailed analysis of these simultaneous upfield-downfield NMR shifts of the same aromatic molecule upon inclusion is reported. Interestingly, the results demonstrate that the observed upfield and downfield shifts upon inclusion are not related to any direct magnetic role of β-CD. The conformational change of KP upon the inclusion process consists of a slight reduction in the angle between the two phenyl rings and in a remarkable reduction in the mobility of the carboxyl group, the latter being one of the main contributions to the NMR resonance shifts. These structural details help in understanding the features of the inclusion complex and, eventually, the driving force for its formation.
Collapse
Affiliation(s)
- T Guzzo
- C4T (Colosseum Combinatorial Chemistry Centre for Technology) S.C.a r.l. , Via della Ricerca Scientifica s.n.c., 00133 Rome, Italy
| | - W Mandaliti
- Department of Chemical Science and Technology, University of Rome "Tor Vergata" , 00133 Rome, Italy
| | - R Nepravishta
- Department of Chemical Science and Technology, University of Rome "Tor Vergata" , 00133 Rome, Italy.,Department of Chemical Pharmaceutical and Biomolecular Technologies, Faculty of Pharmacy Catholic University "Our Lady of Good Counsel" , Rr. D. Hoxha, 1000 Tirane, Albania
| | - A Aramini
- Department of Discovery, Dompé Farmaceutici SpA Research Center , Via Campo di Pile, 67100 L'Aquila, Italy
| | - E Bodo
- Department of Chemistry, University of Rome "Sapienza" , 00133 Rome, Italy
| | - I Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila , Via Vetoio (Coppito 1), 67010 L'Aquila, Italy
| | - M Allegretti
- Department of Discovery, Dompé Farmaceutici SpA Research Center , Via Campo di Pile, 67100 L'Aquila, Italy
| | - A Topai
- C4T (Colosseum Combinatorial Chemistry Centre for Technology) S.C.a r.l. , Via della Ricerca Scientifica s.n.c., 00133 Rome, Italy
| | - M Paci
- Department of Chemical Science and Technology, University of Rome "Tor Vergata" , 00133 Rome, Italy
| |
Collapse
|