key: cord-0697492-mwf0fjua
authors: Ekins, Sean; Mottin, Melina; Ramos, Paulo R.P.S.; Sousa, Bruna K.P.; Neves, Bruno Junior; Foil, Daniel H.; Zorn, Kimberley M.; Braga, Rodolpho C.; Coffee, Megan; Southan, Christopher; Puhl, Ana C.; Andrade, Carolina Horta
title: Déjà vu: Stimulating open drug discovery for SARS-CoV-2
date: 2020-04-19
journal: Drug Discov Today
DOI: 10.1016/j.drudis.2020.03.019
sha: 643cc95b17a80db19513d6943567196d082b3e0d
doc_id: 697492
cord_uid: mwf0fjua

In the past decade we have seen two major Ebola virus outbreaks in Africa, the Zika virus in Brazil and the Americas and the current pandemic of coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). There is a strong sense of déjà vu because there are still no effective treatments. In the COVID-19 pandemic, despite being a new virus, there are already drugs suggested as active in in vitro assays that are being repurposed in clinical trials. Promising SARS-CoV-2 viral targets and computational approaches are described and discussed. Here, we propose, based on open antiviral drug discovery approaches for previous outbreaks, that there could still be gaps in our approach to drug discovery.

Q2 The two major Ebola virus outbreaks in Africa [1] , the Zika virus in Brazil [2] and the current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [3] , which was declared by the World Health Organization (WHO) as a pandemic [4, 5] , have all occurred in the past decade. In each case there is a strong sense of déjà vu for those scientists involved in drug discovery. A lack of available antivirals to treat the infected patients leads to a clamoring to test anything available, and some pharmaceutical companies charge in to offer their drugs. We also seem to see similar patterns in response across outbreaks. There is a rush to be 'first' and this sense of priority might not lead to the best or even any outcome for patients. Always, there is the immediate proposal to create a vaccine and pronouncements that this will be available in a short time or by the end of the year of the actual outbreak in question; and it never happens within these optimistic artificial deadlines. Again, we have experienced this with the current outbreak. For Ebola the vaccine was ready for the second outbreak and has now been approved [6] .

Governments are out to calm their populations while, at the same time, needing to be seen to do something that will vanquish the virus. In the case of SARS-CoV-2 it results in pneumonia [7] and shares aspects of pathology and pathogenesis with severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) [8] . SARS-CoV-2, SARS-CoV and MERS-CoV belong to the same family Coronaviridae and genus Betacoronavirus. Several drug companies immediately made their antivirals available. For example, Abbvie made their HIV approved drug lopinavir/ritonavir combination available to the Chinese government [9] . This drug had previously been suggested as promising based on combination study results from a marmoset study [10] . Publications have already appeared describing the use of this drug in a patient in Korea, significantly decreasing the viral load [11] . Simultaneously, Gilead made available remdesivir, which had previously failed in clinical trials for Ebola [12] but had recently shown activity against MERS in rhesus macaques [13] . Remdesivir, chloroquine and several other drugs were tested in vitro against SARS-CoV-2 (Table 1 ) and shown to be active [14] . Some researchers had also suggested as early as January 2020 what treatment options might be most likely and these included lopinavir/ritonavir, remdesivir, favilavir, arbidol, as well as a broad array of nucleoside analogs, neuraminidase inhibitors, peptides, RNA synthesis inhibitors, anti-inflammatory drugs and traditional Chinese medicines [15] [16] [17] .

The need for other treatment options will be important as we learn more about its transmission [18] and as we see increases in the number of travel-related infections [19] . Already in the past weeks we have seen the rapid spread to Italy and from there to many other countries including Brazil and beyond. Such virus outbreaks are greatly straining the healthcare systems of affected countries and will test every aspect of our infrastructure, supply chains and preparedness [20] . This raises the question of what we can learn from our past drug discovery efforts in an attempt to develop additional drugs for SARS-CoV-2 more efficiently. The acceleration of global SARS-Cov-2 data generation, effectively just in a couple of months, has been remarkable. However, scientific information gaps that impede progress toward small-molecule treatments can still be perceived. Many of these persist from when they were first outlined for Ebola back in 2015 [21] . We can thus make some 'then-and-now' comparisons (see Supplementary material Figs. S1-S4 and Table S1 online) as well as share our own experiences of research on these virus outbreaks.

Our efforts in Ebola drug discovery initially started in 2014 with a computational pharmacophore analysis of a small number of in vitro Ebola active compounds [22] . This was followed by identifying and summarizing all the FDA-approved drugs that could be used against the virus [23, 24] and highlighting various strategies for the next virus outbreak [21] . These steps preceded a pivotal point for us in developing machine learning models for the Ebola virus derived from in vitro data [25, 26] .

An early drug [22] identified by these screens was the antimalarial amodiaquine, which was subsequently shown to be associated with decreased mortality [27] , as the drug (artesunate/ amodiaquine) was used for malaria treatment in some Ebola patients, whereas others took a different malaria medicine (artemether/lumefantrine) [27] . Our Ebola machine learning models were used to select three molecules for in vitro testing [28] . We identified pyronaridine, tilorone and quinacrine as having good in vitro activity (nM) against the Ebola virus [28] . These preliminary data enabled us to obtain funding from the NIH to take one of the compounds through in vivo testing. We also leveraged NIH support to a collaborator to test the other two compounds as well. These molecules were then each tested in the mouse model for Ebola infection and demonstrated significant efficacy [29] [30] [31] . Pyronaridine is currently being pursued in larger animal models of Ebola virus infection. Working on computational models through to in vitro testing happened in the space of a few months, whereas it took several years to obtain funding for our first mouse studies. This is by no means a streamlined approach to drug discovery but it was cost effective for the amount of data ultimately generated and led to revitalized interest in these molecules.

In 2016 we pulled together a team of researchers in Brazil and the USA to provide some suggestions for an open drug discovery effort for the Zika virus. These included various computational strategies to repurpose molecules and docking into the Zika proteins [32] . We also described resources and molecules that could be prioritized for testing. This was followed by our homology modeling of every Zika virus protein [33] weeks before the first crystal structures for Zika proteins were released. All of this work was performed in the open and it attracted the attention of the IBM World Community Grid team who requested we submit a proposal to them for a Zika project using their distributed computing facilities. This became the OpenZika project led by our team in Brazil and the USA and it involved the docking of millions of molecules from the ZINC database in all Zika virus and related flavivirus protein crystal structures [34] . We have since also performed further modeling of Zika proteins [35] and an analysis of the first 2 years following the outbreak [36] . Our most recent efforts have led to the discovery of an inhibitor for NS3 helicase [37] . What we can learn from the OpenZika project to streamline future open science projects would be as follows. We would pick fewer targets to dock molecules into as the OpenZika project created billions of docking results for tens of targets that will take many years to process. We would also suggest that it is best to work with high-resolution crystal structures versus homology models as soon as possible. We would also start with a library of commercially available drugs and antivirals, this would massively limit the number of molecules docked to a few thousand because we found that docking of the ZINC library available at that time led to compounds that were not readily available. Focused libraries are interesting options to perform target-directed virtual screening. Examples are ChemDiv [38] and Asinex [39] databases, which have specific subsets for protease inhibitors, for example facilitating and accelerating the experimental screening. During OpenZika we learnt that having multiple biology labs committed to doing testing in vitro is also very important to have access to experimental results more quickly to validate the computational approach. Working with the IBM World Community Grid did not include any funding. We found it very difficult to obtain NIH funding to support this project from the very beginning (because it was computational) and also attempted small business grants with academic collaborators without success, even though we were leveraging massive computational resources through the World Community Grid. We would recommend that open efforts like this need some funding for in vitro validation. Millions of dollars have been invested by the NIH to fund other Zika research projects since the 2016 outbreak and there is as yet no small molecule of note in the drug discovery pipeline [36] . There is still therefore an unmet medical need for treatments for Zika virus and computational efforts could aid this and be cost-effective.

The whole genome sequences for SARS-CoV-2 were isolated from patients from several countries, such as Brazil, Canada, China, Germany, Korea and USA, and were made rapidly available at the Global Initiative on Sharing All Influenza Data (GSAID) platform [40] . There were 135 genomes sampled between December 2019 and February 2020. The samples were found to be closely related with few mutations relative to a common ancestor. Preliminary genetic analysis indicates that the Brazilian genome differs by three mutations to the Wuhan reference strain [41] . Two of these mutations are shared with its closest sequence-the Germany strain [41] .

In the space of a few weeks after the acknowledgement of the latest coronavirus outbreak, several laboratories started to deposit homology models for the main viral protease [42, 43] , and then for all the SARS-CoV-2 proteins using the server I-TASSER [44] . Others performed docking in these models [45] or developed their own models to understand how the virus enters cells by modeling the Spike (S) protein and human angiotensin-converting enzyme 2 (ACE2) protein interaction [46] . At the time of writing there were structures available for the SARS-CoV-2 main protease in complex with an inhibitor N3 (PDB ID 6LU7) [47] and S1 and S2 subunits of Spike glycoprotein receptor-binding domain up (6VSB), the post fusion core subunit (6LXT) and the HR2 domain (6LVN) [48, 49] . A useful overview is provided by searching the entire virus sequence against the Protein Data Bank (PDB) (see Supplementary material Fig. S2 online). In comparison, for Zika it took several months from the beginning of the outbreak until there was a single protein available in the PDB [50] .

The Spike glycoprotein is responsible for the fusion of the viral membrane with the host-cell membrane. During virus maturation, it is cleaved to its subunits: S1, which attaches the virion to the cell membrane by interacting with the host receptor; and S2, which mediates the fusion of the virion and cellular membranes. A recent study reports that SARS-CoV-2 Spike protein binds ACE2 with higher affinity than SARS-CoV [48] .

Because there was only one Cryo-EM structure of SARS-CoV-2 Spike protein with low resolution (3.46 Å ) without the human ACE2 protein [48] at the time of writing this paper, we preferred to build a model by homology modeling (see Supplementary material methods online). Using our SARS-CoV-2 Spike protein model and the available SARS-CoV Spike protein crystal structure (PDB: 2AJF), we aligned both structures to compare the residues in the Spike-ACE2 interface residues for both viruses. As can be seen in Fig. 1 , the residues for SARS-CoV are Asn473, Tyr475, Tyr440, Tyr436, Tyr484, Leu472, Tyr442 and Asn479 and for SARS-CoV-2 are Asn167, Tyr169, Tyr133, Tyr129, Gln178, Phe166, Leu135 and Gln173. Moreover, some SARS-CoV-2 residues, such as Gln178, Phe166, Leu135 and Gln173, are not conserved but have the same physical-chemical features, except the residues SARS-CoV-2 Gln178 (polar) and SARS Tyr484 (hydrophobic) (Fig. 1) .

The main protease (M pro ) of coronaviruses is a cysteine protease that participates in the viral replication, cleaving the polyprotein. M pro is organized into three domains and is active in a homodimer form [51] . We aligned the crystal structures of M pro of SARS-CoV, MERS-CoV and SARS-CoV-2 to compare the structures and the residues of the catalytic site. The catalytic site for SARS-CoV and SARS-CoV-2 are the same: His41-Cys145 [52] , whereas for MERS-Cov it is His41-Cys148 [52] (Fig. 2) . This catalytic site is localized between the interface of domain I (residues 8-101) and II (102-184). The third domain (residues 201-303) is C-terminal, and consists of a loop region, required for the dimerization [53] [54] [55] [56] .

Recently, it has been reported that the structures of M pro of SARS-CoV and SARS-CoV-2 have high identity of 96% [43] . The function and the structural equivalence of M pro of SARS-CoV, MERS-CoV and SARS-CoV-2 demonstrate that the coronavirus proteases are promising targets for development of broad-spectrum antivirals in these viruses.

Drug repurposing offers a potentially faster approach to identify drugs that are already approved for other uses and might work for the disease of interest [57] [58] [59] . This type of strategy was attempted for MERS-CoV and SARS-CoV [60, 61] and the repurposing of hostbased therapies has also been proposed for coronaviruses [62] . Several studies have screened a few hundred molecules and tested the activity against MERS-CoV and SARS-CoV (Table 2) [63, 64] .

The first screen of 348 FDA-approved drugs identified four FDAapproved drugs including chloroquine, chlorpromazine, loperamide and lopinavir [63] . At the same time, a second study screened 290 compounds and identified 27 FDA-and non-FDA-approved drugs including chloroquine among the many hits [64] , and the Abl2 kinase inhibitors such as imatinib which impacts the early stages of infection after internalization and endosomal trafficking. This latter molecule appears to work by inhibiting fusion at the viral membrane [65] . Chloroquine is widely used for malaria and loperamide is used for diarrhea. These represent drugs that are readily available in many locations around the world and for which physicians are also familiar with the side-effect profile. The ligand datasets from these and other screens can be used to apply machine learning methods to produce predictive models that can then be used to screen compound libraries and suggest Superposition of SARS-CoV and SARS-CoV-2 spike proteins highlighting the spike-ACE2 interface binding pocket. Zoom at the residues of SARS-CoV-2 spike S1 Asn167, Tyr169, Tyr133, Tyr129, Gln178, Phe166, Leu135 and Gln173 and the correspondent SARS-CoV spike S1 Asn473, Tyr475, Tyr440, Tyr436, Tyr484, Leu472, Tyr442 and Asn479.

Superposition of SARS-CoV, MERS-CoV and SARS-CoV-2 main proteases, and respective catalytic biding site. SARS-CoV (PDB ID 1Q2W) and SARS-CoV-2 (PDB ID 6LU7) share the same residues at the catalytic site: His41-Cys145, and MERS-CoV (PDB ID 5WKL) His41 and Cys148, situated between domains I and II. (Fig. 3 , the experimental details are described in the Supplementary material online; the Bayesian models are also available at http://models.collabchem.com/8/). The advantage of this approach is that the method is fast, does not require crystal structures and it enables small-molecule structures to be scored against many models simultaneously.

As a well-established principle any SARS-CoV-2 relevant enzyme assay, protein-protein interaction antagonism or viral infectivity determination can be run against a collection of approved drugs. However, it is neither trivial to get a definitive list of these structures in silico prepared for docking or to acquire these as a physical collection ready for screening (even though there are various repurposing subsets available from vendors). Notwithstanding, an approximation to a definitive approved set can be 

An alternative to screening every approved drug is to take one class of compounds. For example we have now taken the opportunity to dock FDA-approved HIV-1 protease inhibitors (amprenavir, atazanavir, darunavir, fosamprenavir, indinavir, lopinavir, nelfinavir, ritonavir, saquinavir and tipranavir) and hepatitis C NS3/4A protease inhibitors (asunaprevir, boceprevir, grazoprevir, paritaprevir, simeprevir and telaprevir) into the M pro structure of SARS-CoV-2.

The docking was performed with two different commercial docking software products: Biovia Discovery Studio and Schrödinger Glide [70] [71] [72] [73] , and the best docking results were compared (see Supplementary material Table S2 online). The results of both docking approaches indicated that azatanavir and lopinavir could be of most interest for in vitro or in vivo testing (Fig. 4) . Atazanavir can interact with M pro by H-bonding at Cys145 and the amide group. The aromatic ring of atazanavir (LigandFit score of 191.196) also formed hydrophobic interactions with residues Cys44, Met49, Leu50, Tyr54 and Cys145. In relation to lopinavir (LigandFit score of 175.497), H-bonds were established between residues Cys145 and Glu166 and nitrogen atoms of carbamate groups. Similar hydrophobic interactions were observed with lopinavir (Fig. 4) . In China and South Korea lopinavir 200 mg/ ritonavir 50 mg has already been used for the treatment of COVID-19. From the day after lopinavir/ritonavir administration, the viral load started to decrease and no detectable or little coronavirus titers have been observed since then [11] . Lopinavir, co-formulated with ritonavir is commercially available and known as Kaletra1 around the world-it is on formularies for HIV. By contrast, atazanavir is less commonly used. Kaletra1 can be widely produced and distributed quickly through existing channels but it will be important to ensure HIV patients do not lose access to the medications they need. It is also a medication that clinicians are familiar with in terms of side effects, unlike other less commonly used medications. A broader analysis of all FDA-approved drugs docked in this protease structure could point to further drugs worth testing. We have also performed docking calculations of SARS-CoV-2 spike glycoprotein complexed with the host ACE2 protein against all FDA-approved drugs ($2400 molecules) using Schrödinger Glide software [70] [71] [72] [73] . Because there was only one Cryo-EM structure of SARS-CoV-2 Spike protein with low resolution (3.46 Å ) without the human ACE2 protein [48] at the time of doing this analysis, we preferred to build a model by homology modeling, because lower-resolution structures (3 Å or higher) show the basic contours of the protein chain and the atomic structure must be inferred. We built the SARS-CoV-2 Spike subunit S1 model using the SwissModel server [74] . The template used was the SARS-CoV Spike S1 protein (PDB ID 2GHV) [75] , which has resolution of 2.2 Å and 73.22% of sequence identity with SARS-CoV-2 Spike protein. After the protein refinement, we added the human ACE2 protein to the model, based on the SARS-CoV Spike-ACE2 structure (PDB ID 2AJF) [76] , and built the full SARS-CoV-2 Spike S1-ACE2 model (see Fig. S5 in Supplementary material online). The quality statistics of the Intermolecular 3D and 2D interactions of atazanavir (a) and lopinavir (b) with the main protease structure of SARS-CoV-2 (PDB ID 6LU7) obtained by docking. In the 3D representation (left), H-bonds are presented as magenta dashed lines and the color code of carbon, oxygen, nitrogen and hydrogen atoms are grey, red, blue and white, respectively. In the 2D interaction diagrams (right) H-bonds are presented as magenta arrows. The Pymol software [79] was used for the visual inspection of 3D docking poses and to render the 3D molecular images. 3D and 2D intermolecular interactions of rizatriptan (a), dasabuvir (b), pravastatin (c) and empagliflozin (d) with the spike-ACE2 interface (light blue and green, respectively). In 3D representation (left), H-bonds are presented as magenta dashed lines and the color code of carbon, oxygen, nitrogen and hydrogen atoms is grey, red, blue and white, respectively. The SARS-CoV-2 Spike (S1 subunit) is colored in blue and human ACE2 in light blue. In the 2D ligand-interaction diagram (right), H-bonds, p-p stacking and p-cation interactions are presented as magenta, green and red lines, respectively. Glide software [71, 80] was used to dock the ligands in a grid centered at Asn167, Tyr169, Tyr133, Tyr129, Gln178, Phe166, Leu135 and Gln173 residues.

www.drugdiscoverytoday.com 9 modeled SARS-CoV-2 Spike-ACE2 were good: 96.53% of the residues lie in the most favorable regions of the Ramachandran plot and the Molprobity score was 1.66. Our SARS-CoV-2 Spike-ACE2 model is available (see Supplemental material PDB file online). Subsequently, we performed docking-based virtual screening with a grid centered at the Spike-ACE2 interface. As can be seen in Table 3 , docking scores indicated that rizatriptan (an antimigraine agent, Glide score = À7.49 kcal mol À1 ), dasabuvir (an antiviral for hepatitis C, Glide score = À7.92 kcal mol À1 ), pravastatin (a lipid-lowering agent, Glide score = À7.81 kcal mol À1 ) and empagliflozin (a hypoglycemic agent, Glide score = À7.09 kcal mol À1 ) could be of interest for in vitro testing. Details of intermolecular interactions of these drugs with the interface of Spike-ACE2 are shown in Fig. 5 . In summary, these drugs form p-cation interactions with Lys31a, pp interactions with Phe170b and H-bonds with His34a, Asp38a, Gly176b and Ser174b. In good agreement with our results, a previous study suggests that an early and high dose of a statin (such as pravastatin) could be crucial to surviving MERS-CoV infection [77] . These suggested molecules could be validated in future in vitro experiments.

Modern drug discovery efforts try to eliminate inefficiency by increasing knowledge, through extensive utilization of generated data to guide new experiments. Therefore, we have also taken a drug repurposing approach based on polypharmacology using a method called the 'molecular and biological signature' [78] . This is a profiling approach, guided by artificial intelligence, capable of presenting the complex relationships between the chemical substance with molecular targets, cells and/or organs. This could create a more complete understanding of the 'multi-biological profile' caused by the drug analyzed in response to a disease. As an example, we have applied the drugs ritonavir [9] , favipiravir, ribavirin and chloroquine as scaffolds to find new potential molecular and biological signature profiles, to search for non-obvious, mechanism-of-action-based associations between drugs, molecular targets and diseases ( Fig. 6; and Supplementary material Fig. S6 and Table S3 online). Comparing the molecular and biological signature of ritonavir with other drugs, we found that the drugs cobicistat, carfilzomib and ombitasvir presented the highest matching scores (81-71%), suggesting that those drugs could be of further interest for SARS-CoV-2. 

Molecular and biological signature profiles generated for ritonavir compared with other drugs, using artificial intelligence to create thousands of decisionmakers that mimic a biological response through the interaction of a small molecule with multiple enzyme targets, off-targets and in vitro and in vivo endpoints. The drugs cobicistat, carfilzomib and ombitasvir presented the higher matching scores (81-71%) with ritonavir.

Although there are already several drugs being assessed clinically for SARS-CoV-2, there is still a need to identify additional treatments as alternatives (to offset drug shortages) and possibly discover molecules with increased efficacy, safety and costeffectiveness. We have previously proposed and applied computational drug discovery approaches such as machine learning and docking to aid in the discovery of new candidate molecules for Ebola and Zika, respectively. What we have learnt in these efforts can be applied to SARS-CoV-2 as we have demonstrated here. However, the research in this case has moved at a much faster pace initially, delivering a first crystal structure in just a few weeks after first public notification of the virus. Companies have mobilized to bring their drugs to patients and clinicians have used prior in vitro research that proposed compounds active against MERS-CoV and SARS-CoV to suggest treatments for SARS-CoV-2. We currently do not clinically use any broad-spectrum antivirals in the USA but they would be invaluable in situations like this for SARS-CoV-2 and other viruses. Identifying medications with side-effect profiles known by clinicians that are already produced at scale and available at local pharmacies around the world will allow more patients to be treated, especially as travel restrictions and quarantines could reduce production and distribution of certain drugs. It will remain important that those who need these medications at baseline, such as HIV patients who take ritonavir/lopinavir, are able to continue to access the medications they need or else there is a risk of poor clinical outcomes and antiviral resistance in these patients. Although other viruses such as the flu might have killed more patients in the USA this season so far, the crippling financial impact of SARS-CoV-2 remains to be seen globally. If past Ebola, ZIKV, MERS-CoV and SARS-CoV outbreaks are anything to go by this SARS-CoV-2 pandemic will be far more significant. Having multiple effective treatments available for SARS-CoV-2 that can be rapidly delivered to patients could save lives. We can also supplement these treatments by a combination of computational screening with machine learning and docking to select additional drugs to test as proposed herein and this can be repeated for future virus outbreaks. The sense of déjà vu will probably be repeated again and again with the next virus outbreaks until we develop broaderspectrum antivirals. Each time we hope to be better prepared. Hope however does not save lives. Effective antiviral treatments do, and we need to have them ready in the future. 

Ebola virus: pathogenesis and countermeasure development

Zika virus infection -after the pandemic

Severe acute respiratory syndrome-related coronavirusthe species and its viruses, a statement of the Coronavirus Study Group

First case of 2019 novel coronavirus in the United States

WHO Director-General's Opening Remarks at the Media Briefing on COVID-19

First FDA-approved vaccine for the prevention of Ebola virus disease, marking a critical milestone in public health preparedness and response

Initial CT findings and temporal changes in patients with the novel coronavirus pneumonia (2019-nCoV): a study of 63 patients in Wuhan, China

Overlapping and discrete aspects of the pathology and pathogenesis of the emerging human pathogenic coronaviruses SARS-CoV, MERS-CoV, and 2019-nCoV

AbbVie Statement on Q5 Coronavirus and Lopinavir/Ritonavir | AbbVie News Center

A systematic review of therapeutic agents for the treatment of the Middle East respiratory syndrome coronavirus (MERS-CoV)

Case of the index patient who caused tertiary transmission of COVID-19 infection in Korea: the application of lopinavir/ritonavir for the treatment of COVID-19 infected pneumonia monitored by quantitative RT-PCR

A randomized, controlled trial of Ebola virus disease therapeutics

Prophylactic and therapeutic remdesivir (GS-5734) treatment in the rhesus macaque model of MERS-CoV infection

Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro

Drug treatment options for the 2019-new coronavirus (2019-nCoV)

Potential interventions for novel coronavirus in China: a systemic review

Potential antiviral therapeutics for 2019 novel coronavirus

Initial cluster of novel coronavirus (2019-nCoV) infections in Wuhan, China is consistent with substantial human-to-human transmission

Persons evaluated for 2019 novel coronavirus -United States

Straining the system: novel coronavirus (COVID-19) and preparedness for concomitant disasters

Finding small molecules for the

A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus

Small molecules with antiviral activity against the Ebola virus. F1000Research

FDA approved drugs as potential Ebola treatments

Evaluation of Ebola virus inhibitors for drug repurposing

A systematic screen of FDA-approved drugs for inhibitors of biological threat agents

Effect of artesunate-amodiaquine on mortality related to Ebola virus disease

Machine learning models identify molecules active against the Ebola virus in vitro

Efficacy of tilorone dihydrochloride against Ebola virus infection

Repurposing quinacrine against Ebola virus infection in vivo

Repurposing the antimalarial pyronaridine tetraphosphate to protect against Ebola virus infection

Open drug discovery for the Zika virus

Illustrating and homology modeling the proteins of the Zika virus

OpenZika: an IBM World Community Grid project to accelerate Zika virus drug discovery

Molecular dynamics simulations of Zika virus NS3 helicase: insights into RNA binding site activity

The A-Z of Zika drug discovery

A diarylamine derived from anthranilic acid inhibits ZIKV replication

GISAID Initiative the Genome Sequences of the Newly Emerging Coronavirus

First Report of COVID-19 in South America

Comparative Model of Novel Coronavirus 2019-nCoV Protease Mpro

Homology models of Wuhan coronavirus 3CLpro protease

Structure Models of All Mature Peptides in COVID-19 Genome by C-I-TASSER

Deep learning based drug screening for novel coronavirus 2019-nCov

Evolution of the novel coronavirus from the ongoing Wuhan outbreak and modeling of its spike protein for risk of human transmission

The Crystal Structure of COVID-19 Main Protease in Complex with an Inhibitor N3

Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation

Crystal Structure of HR2 Domain of 2019-nCoV S2 Subunit

Zika virus NS1 structure reveals diversity of electrostatic surfaces among flaviviruses

Quaternary structure of the severe acute respiratory syndrome (SARS) coronavirus main protease

Structure-guided design of potent and permeable inhibitors of MERS coronavirus 3CL protease that utilize a piperidine moiety as a novel design element

The 3C-like proteinase of an invertebrate nidovirus links coronavirus and potyvirus homologs

The crystal structures of severe acute respiratory syndrome virus main protease and its complex with an inhibitor

Dissection study on the severe acute respiratory syndrome 3C-like protease reveals the critical role of the extra domain in dimerization of the enzyme

Structures of two coronavirus main proteases: implications for substrate binding and antiviral drug design

In silico chemogenomics drug repositioning strategies for neglected tropical diseases

New uses for old drugs

In silico repositioning of approved drugs for rare and neglected diseases

Screening of FDA-approved drugs for treatment of emerging pathogens

Middle East respiratory syndrome and severe acute respiratory syndrome: current therapeutic options and potential targets for novel therapies

Repurposing host-based therapeutics to control coronavirus and influenza virus

Screening of an FDA-approved compound library identifies four small-molecule inhibitors of Middle East respiratory syndrome coronavirus replication in cell culture

Repurposing of clinically developed drugs for treatment of Middle East respiratory syndrome coronavirus infection

Abelson kinase inhibitors are potent inhibitors of severe acute respiratory syndrome coronavirus and Middle East respiratory syndrome coronavirus fusion

Exploiting machine learning for end-to-end drug discovery and development

Comparing multiple machine learning algorithms and metrics for estrogen receptor binding prediction

Multiple machine learning comparisons of HIV cell-based and reverse transcriptase data sets

Ebola virus Bayesian machine learning models enable new in vitro leads

Extra Precision Glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes

Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening

Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy

SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information

Structural basis of neutralization by a human anti-severe acute respiratory syndrome Spike protein antibody, 80R

Structure of SARS coronavirus Spike receptor-binding domain complexed with receptor

Statins may decrease the fatality rate of Middle East respiratory syndrome infection

Drug repurposing: progress, challenges and recommendations

The PyMOL Molecular Graphics System, Version 1

Schrödinger Release 2019-4: Glide

Supplementary data associated with this article can be found, in the online version, at https://doi.org/10.1016/j.drudis.2020.03.019.