Inhibitors for Novel Coronavirus Protease Identified by Virtual Screening of 687 Million Compounds

The rapid outbreak of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in China followed by its spread around the world poses a serious global concern for public health with almost 90000 people infected and thousands of fatalities To conquer viral infections, the inhibition of proteases essential for proteolytic processing of viral polyproteins is a conventional therapeutic strategy To this date, no specific drugs or vaccines are available to treat SARS-CoV-2 despite its close relation to the SARS-CoV-1 virus that caused a similar epidemic in 2003 Thus, there remains an urgent need for the development of specific antiviral therapeutics to conquer SARS-CoV-2 In order to find novel inhibitors, we computationally screened a compound library of over 687 million compounds for binding at the recently solved crystal structure of the main protease of SARS-CoV-2 A screening of such a vast chemical space for SARS-CoV-2 protease inhibitors has not been reported before After shape screening, two docking protocols were applied followed by the determination of pharmacokinetically relevant molecular descriptors to narrow down the initial hits Next, molecular dynamics simulations were conducted to validate the stability of docked binding modes and comprehensively quantify ligand binding energies After evaluation of off-target binding, we report a list of 11 drug-like compounds with improved binding free energy to the target protease in contrast to the cocrystallized peptidomimetic lead compound that suffers from poor pharmacokinetic properties Furthermore, we identified one potent binder with comparable properties from the natural compound library

[1]  Arun K. Ghosh,et al.  Recent Progress in the Development of HIV-1 Protease Inhibitors for the Treatment of HIV/AIDS. , 2016, Journal of medicinal chemistry.

[2]  Arun K. Ghosh,et al.  Second‐Generation Approved HIV Protease Inhibitors for the Treatment of HIV/AIDS , 2011 .

[3]  R. M. Owen,et al.  An analysis of the attrition of drug candidates from four major pharmaceutical companies , 2015, Nature Reviews Drug Discovery.

[4]  Yu-An Huang,et al.  Synthesis, crystal structure, structure-activity relationships, and antiviral activity of a potent SARS coronavirus 3CL protease inhibitor. , 2006, Journal of medicinal chemistry.

[5]  Ting Yu,et al.  Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study , 2020, The Lancet.

[6]  Kwok-Yin Wong,et al.  Prediction of the SARS-CoV-2 (2019-nCoV) 3C-like protease (3CL pro) structure: virtual screening reveals velpatasvir, ledipasvir, and other drug repurposing candidates , 2020, F1000Research.

[7]  Y. Hayashi,et al.  An Overview of Severe Acute Respiratory Syndrome–Coronavirus (SARS-CoV) 3CL Protease Inhibitors: Peptidomimetics and Small Molecule Chemotherapy , 2016, Journal of medicinal chemistry.

[8]  Yan Li,et al.  Therapeutic Drugs Targeting 2019-nCoV Main Protease by High-Throughput Screening , 2020, bioRxiv.

[9]  David Ryan Koes,et al.  Lessons Learned in Empirical Scoring with smina from the CSAR 2011 Benchmarking Exercise , 2013, J. Chem. Inf. Model..

[10]  Guo-Wei Wei,et al.  Machine intelligence design of 2019-nCoV drugs , 2020, bioRxiv.

[11]  A. Vedani,et al.  VirtualToxLab - a platform for estimating the toxic potential of drugs, chemicals and natural products. , 2012, Toxicology and applied pharmacology.

[12]  V. Zarubaev,et al.  Virus-inhibiting activity of dihydroquercetin, a flavonoid from Larix sibirica, against coxsackievirus B4 in a model of viral pancreatitis , 2016, Archives of Virology.

[13]  Weiliang Zhu,et al.  Nelfinavir was predicted to be a potential inhibitor of 2019-nCov main protease by an integrative approach combining homology modelling, molecular docking and binding free energy calculation , 2020, bioRxiv.

[14]  David S. Goodsell,et al.  AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility , 2009, J. Comput. Chem..

[15]  Haitao Yang,et al.  Structure of Main Protease from Human Coronavirus NL63: Insights for Wide Spectrum Anti-Coronavirus Drug Design , 2016, Scientific Reports.

[16]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[17]  Evan Bolton,et al.  PubChem 2019 update: improved access to chemical data , 2018, Nucleic Acids Res..

[18]  Richard D. Smith,et al.  The economic impact of SARS: How does the reality match the predictions? , 2008, Health Policy.

[19]  G. Gottlieb,et al.  Four Amino Acid Changes in HIV-2 Protease Confer Class-Wide Sensitivity to Protease Inhibitors , 2015, Journal of Virology.

[20]  Young-Min Kim,et al.  Flavonoid-mediated inhibition of SARS coronavirus 3C-like protease expressed in Pichia pastoris , 2012, Biotechnology Letters.

[21]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[22]  Olga Golosova,et al.  Unipro UGENE: a unified bioinformatics toolkit , 2012, Bioinform..

[23]  Z. Memish,et al.  Reducing mortality from 2019-nCoV: host-directed therapies should be an option , 2020, The Lancet.

[24]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[25]  A. Vedani,et al.  OpenVirtualToxLab--a platform for generating and exchanging in silico toxicity data. , 2015, Toxicology letters.

[26]  J. P. Grossman,et al.  Anton 2: Raising the Bar for Performance and Programmability in a Special-Purpose Molecular Dynamics Supercomputer , 2014, SC14: International Conference for High Performance Computing, Networking, Storage and Analysis.

[27]  John J. Irwin,et al.  ZINC 15 – Ligand Discovery for Everyone , 2015, J. Chem. Inf. Model..

[28]  Yu Wai Chen,et al.  Prediction of the 2019-nCoV 3C-like Protease (3CLpro) Structure: Virtual Screening Reveals Velpatasvir, Ledipasvir, and Other Drug Repurposing Candidates , 2020 .

[29]  W Patrick Walters,et al.  Going further than Lipinski's rule in drug design , 2012, Expert opinion on drug discovery.

[30]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[31]  R. Wade,et al.  The role of hydrogen-bonds in drug binding. , 1989, Progress in clinical and biological research.

[32]  Z. Memish,et al.  The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health — The latest 2019 novel coronavirus outbreak in Wuhan, China , 2020, International Journal of Infectious Diseases.

[33]  Yanjie Wei,et al.  Deep Learning Based Drug Screening for Novel Coronavirus 2019-nCov , 2020, Interdisciplinary Sciences: Computational Life Sciences.

[34]  Kai Zhao,et al.  A pneumonia outbreak associated with a new coronavirus of probable bat origin , 2020, Nature.

[35]  E. Lionta,et al.  Structure-Based Virtual Screening for Drug Discovery: Principles, Applications and Recent Advances , 2014, Current topics in medicinal chemistry.

[36]  Xiu-Jie Wang,et al.  Potential inhibitors against 2019-nCoV coronavirus M protease from clinically approved medicines , 2020, Journal of Genetics and Genomics.

[37]  F. Maione,et al.  The Potential Therapeutic Application of Peptides and Peptidomimetics in Cardiovascular Disease , 2017, Front. Pharmacol..

[38]  Sanjivanjit K. Bhal,et al.  The Rule of Five revisited: applying log D in place of log P in drug-likeness filters. , 2007, Molecular pharmaceutics.

[39]  Zhènglì Shí,et al.  The crystal structure of COVID-19 main protease in complex with an inhibitor N3 , 2020 .