Computationally repurposed drugs and natural products against RNA dependent RNA polymerase as potential COVID-19 therapies
Sakshi Piplani, Puneet Singh, David A. Winkler, Nikolai Petrovsky
CComputationally repurposed drugs and natural products against RNA dependent RNA polymerase as potential COVID-19 therapies
Sakshi Piplani , Puneet Singh , David A. Winkler , Nikolai Petrovsky bstract For fast development of COVID-19, it is only feasible to use drugs (off label use) or approved natural products that are already registered or been assessed for safety in previous human trials. These agents can be quickly assessed in COVID-19 patients, as their safety and pharmacokinetics should already be well understood. Computational methods offer promise for rapidly screening such products for potential SARS-CoV-2 activity by predicting and ranking the affinities of these compounds for specific virus protein targets. The RNA-dependent RNA polymerase (RdRP) is a promising target for SARS-CoV-2 drug development given it has no human homologs making RdRP inhibitors potentially safer, with fewer off-target effects that drugs targeting other viral proteins. We combined robust Vina docking on RdRP with molecular dynamic (MD) simulation of the top 80 identified drug candidates to yield a list of the most promising RdRP inhibitors. Literature reviews revealed that many of the predicted inhibitors had been shown to have activity in in vitro assays or had been predicted by other groups to have activity. The novel hits revealed by our screen can now be conveniently tested for activity in RdRP inhibition assays and if conformed testing for antiviral activity invitro before being tested in human trials. ntroduction
Since December 2019 and the outbreak of a pneumonic disease called COVID-19 caused by a novel coronavirus (SARS-CoV-2) in China, clinicians globally have faced unprecedented challenges in managing the disease caused by this pandemic virus. There has been a massive international research effort to discover effective drugs and vaccines for this and other pathogenic coronaviruses such as SARS and MERS.
Design of potent new coronavirus drugs is very important for future disease preparedness. However, for COVID-19, it is only feasible to use drugs that are already registered (off label use) or approved natural products that been through at least phase 1 clinical trials for safety assessment. These agents can be used in man quickly, as their safety and pharmacokinetics are well understood. Computational methods offer considerable promise for rapidly screening such drugs for SARS-Cov-2 activity by predicting their affinities for relevant virus protein targets. Recent papers in prominent journals have reported, for example, the application of computational de novo drug design based on the structures of the SARS-Cov-2 protease.
The RNA-dependent RNA polymerase (RdRP) is a promising target for SARS-CoV-2 drug development. It has no host cell homologs so selective RdRP inhibitors should be potentially safer, with lower risks of off-target effects. Zhu et al., recently reviewed the biochemical properties of this critical enzyme and cell-based RdRPs assays suitable for high-throughput screening to discover new and repurposed drugs against SARS-CoV-2. Viral RdRP plays a crucial role in the SARS-Cov-2 replicative cycle. Its active site is highly conserved and accessible, making it a promising drug target. Notably, all DNA and RNA viruses employ RdRP proteins for replication and transcription of viral genes and other viral and host factors (Figure 1), suggesting that computational techniques that prove useful for identifying repurposed drugs against COVID-19 RdRP should be broadly applicable to other pathogenic viruses. igure 1 . Interaction of SARS-Cov-2 with cells showing effects of inhibition of RNA replication by an RdRP inhibitor. Created using BioRender template. RdRPs share several sequence motifs and tertiary structures amongst all RNA virus types, including positive-sense RNA, negative-sense RNA and dsRNA viruses. The core structure of RdRP resembles a right-hand, with palm, thumb and finger domains. Five of the seven classical catalytic RdRP motifs (A – E) are in the most preserved palm domain, while the remaining two (F and G) are in the finger domains. Five of the seven classic catalytic RdRP (A-E) motifs are in the most highly conserved palm domain, while the other two (F and G) are in finger domain (Figure 2). The structurally conserved RdRP core and related motifs are important for the catalytic role of viral RdRP and thus represent potential drug targets. While the criteria for the substrates iffer, all known RdRPs share the same catalytic mechanism. On host cell infection, viral RdRP participates in formation of the molecular machinery for genome replication by complexing with other factors. It initiates and regulates the elongation of the RNA strand, which involves the addition of hundreds to thousands of nucleotides. When incorporated into the newly synthesised RNA chain, nucleotide analogues such as remdesivir block the RNA elongation catalysed by RdRP (Figure 2). Figure 2 . Inhibition of RNA-dependent RNA polymerase (RdRp) by repurposed drugs. (A) The SARS-Cov-2 genome composition single strand RNA model. (B) RdRp mediated RNA replication during coronavirus infection. (C) Drugs such as remdesivir inhibit RdRp and block RNA replication. Creative Commons licence figure adapted from Huang et al. In pandemics, computational methods can quickly identify candidate drugs for repurposing, where speed is of utmost importance. Here we show how a well validated molecular docking ollowed by a high throughput molecular dynamics simulation (MDS) drug screening program can screen a large number of drugs and natural products to produce a short list of RdRp-inhibiting drug candidates. MDS calculations were used to predict the optimal binding poses and binding energies for 80 of the top hits from virtual screening on SARS-CoV-2 RdRP. Finally we ranked the top candidates for COVID-19 based on binding affinity and novelty.
Results and discussion
There have been many computational studies aiming to predict which existing drugs and natural products may be inhibit the main protease, M pro , of SARS-Cov-2, but fewer studies have studied drugs for targeting the viral RdRP. The computational workflow used to estimate the relative binding affinities of drugs for the RdRP binding pocket are summarized in Figure 3. The combination of robust Vina docking followed by MDS of the top 80 candidates yields improved performance relative to Vina docking alone. The binding free energies calculated by both computational methods (see Methods) correlated very well (r =0.84) and the free energies calculated by the thermodynamic cycle correlate with the Vina docking scores with r = 0.64. The size of the binding site (area 2920 Å and volume 5335 Å ) tends to select larger ligands, many of which are quite flexible. Binding energy penalties due to ligand entropy are likely to be significant. Hence, substantial correlation between the Vina scores and the binding energies from MMBBSA and thermodynamic cycle are important because these algorithms treat ligand entropy approximately and in different ways. Figure 3 . Computational workflow for repurposing drugs against SARS-Cov-2 RdRP. The binding energies of the 80 top ranked ligands from the docking calculations are listed in Supplementary Table 1. The 20 drugs predicted to have the highest binding to RdRp are summarized in Table 1, together with both their MMPBSA and thermodynamic cycle binding energies. The high correlation between the two methods of calculating binding free energies meant that the ranking of the top 20 compounds was very similar when ordered by either binding free energy estimate. Antiviral drugs, Paritaprevir, Beclabuvir, Remdesivir, Voxilaprevir, Setrobuvir, Galidesvir, Elbasvir, Ciluprevir, Faldaprevir, and Tegobuvir account for half of the list. The remaining 10 hit compounds, which include several natural products, are used to treat a diverse range of afflictions – cancers, infections, cardiac insufficiency, liver damage, circulatory issues, and parasitic infections. Almost all of the drugs in the top 20 have relatively large, complex structures and substantial ligand flexibility. Some have had their in vitro activity against SARS-Cov-2 determined experimentally (Figure 4), further supporting our binding predictions.
Table 1 . Structures and binding energies of 20 top ranked (by MMPBSA score) small molecule ligands to SARS-Cov-2 RDRP. atabase ID C=ChemBL D=Drugbank DrugName Structure D G MMPBSA ( D G thermo) kcal/mol C 3391662 Paritaprevir (antiviral) -54.3 (-67.5) C 1200633 Ivermectin (anti-parasitic) -54.1 (-69.7) C 3126842 Beclabuvir (antiviral) -53.5 (-66.4) C 3809489 Bemcentinib (anticancer) -46.1 (-62.5) D 14761 Remdesivir (antiviral) -44.6 (-56.8) C 1751 Digoxin (cardiac drug) -41.2 (-54.4) D 09298 Silibinin (chemo-protectant -40.2 (-57.2) C 1236524 Galidesvir (antiviral) -40.2 (-49.2) H NHO HNONHON N O NO HN SO O
OOHO OOOH O O O O OHOOO
ONO N NOHNSO ON H N N HN NNNN N
OOHNPOOONHO HONNNH N O
H HH OHOH OOOOHOHOOOHHO O O O
HOO OHHOOHHO O OH OOH
H NH OHOHHO HNNN NH atabase ID C=ChemBL D=Drugbank DrugName Structure D G MMPBSA ( D G thermo) kcal/mol C 1076263 Setrobuvir (antiviral) -40.0 (-51.1) C 3707372 Voxilaprevir (antiviral) -39.5 (-52.8) C 2013174 Vedroprevir (antiviral) -38.4 (-44.6) C 1241348 Faldaprevir (antiviral) -38.1 (-46.4) C 413 Sirolimus (Rapamycin) (immuno-suppressant) -37.4 (-45.2) SO OHNNHN SOOOHNF O
N O NHF FOHN SOOONHO O FFO N ON
NHO N O NH O OOONO Cl NSHN N OHO
OBr O N O NH O OO NH OHONSNHO N
H H O OHO OHHO OOONH OO OHO O atabase ID C=ChemBL D=Drugbank DrugName Structure D G MMPBSA ( D G thermo) kcal/mol C 3301668 Carbetocin (anti-hemorrhagic) -37.4 (-35.1) D 01051 Novobiocin (antibiotic) -37.3 (-45.3) C 1683590 Eribulin (anticancer) -36.5 (-48.3) C 442 Ergotamine (anti-migraine) -36.2 (-45.6) C 1957287 Tegobuvir (antiviral) -34.7 (-46.0) D 12466 Favipiravir (antiviral) -34.2 (-41.6) D 14850 Deleobuvir (antiviral) -34.1 (-39.6) HNHOO NHO S HN O NH OHN O NH OHNOON O HN O NH NH O OO NH OOHOOHNHOHO O O O
OO O OO O OO OO NH HO HNO N O NH O N H NHOOH
FNNN NFF FFFF N NH ONF N OH
NHO O NHNO N N BrN
Figure 4 . Experimental antiviral spectrum for hit compounds. From DrugVirus.info.
Antiviral agents . The calculated binding energy of several of the antiviral drugs, Paritaprevir and Beclabuvir, are very similar, within the calculation uncertainties in energies. Several of the top antiviral agents predicted by this study, have also been identified in other in silico docking studies of small molecule drugs. This provides a degree of validation that our computational methods are appropriate and are yielding similar results to the other published studies for some well-studied antiviral drugs. Beg and Athar used a homology modelled structure for SARS-Cov-2 RdRP and AutoDock Vina to rank the binding of a range of antiviral drugs to the RdRP enzyme. Paritaprevir, Beclabuvir, and Favipiravir had some of the highest docking scores. Cozac et al. reported a combined docking and machine learning study of inhibitors of RdRP from HCV, poliovirus, dengue virus, and influenza virus. They identified Faldaprevir, Vedroprevir, Beclabuvir and Remdesivir as having good docking scores to SARS-Cov-2 RdRP and predicted viral RdRPs IC /EC values below 5µM by multiple machine learning classification models. A homology model for SARS-Cov-2 RdRP and AutoDock was used by Dutta et al. to rank antiviral drugs for potential use in treating COVID-19 patients. One of the most promising rugs for repurposing was Beclabuvir with a predicted docking score of (-10 kcal/mol) and IC of 50nM, although several other antiviral agents had similar docking scores, e.g. tegobuvir (-9.7 kcal/mol), dasabuvir (-9.4 kcal/mol), lomibuvir (-11 kcal/mol), setrobuvir ( -10.5 kcal/mol) and dasabuvir (-10.4 kcal/mol). Remdesivir (-7.4 kcal/mol), radalbuvir (-7.4 kcal/mol) and dasabuvir (-6.7 kcal/mol) had significantly lower docking scores. None of these studies used MDS refinement of docked compounds to improve the estimates of the binding free energy of the drugs. Remdesivir has been predicted to be active against SARS-Cov-2 RdRP in several studies. There is limited in vitro evidence for activity of remdesivir with an IC
31, 32 and recent clinical studies suggest its efficacy in treating COVID-19 infection in man is limted. Elfiki reported two studies aimed at repurposing existing drugs against RdRP.
29, 30
In one study, AutoDock Vina docking, and MDS (50 ns production runs) was used to model viral RdRp and calculate the binding affinity of several drugs and drug candidates. Existing drugs, Sofosbuvir, Ribavirin, Galidesivir, Remdesivir, Favipiravir, Cefuroxime, Tenofovir, Setrobuvir, and Hydroxychloroquine, were all predicted to bind to RdRp. Setrobuvir, YAK, and IDX-184 had even more favourable binding energies, and four novel derivatives of the latter drug showed binding to SARS-CoV-2 RdRp. The second study used only docking to predict the binding of the same drugs to RdRP from SARS-Cov-2, SARS, and HCV. HCV RdRp had greater binding to the antiviral drugs than did SARS-Cov-2 RdRP, with SARS RdRP having the lowest affinity as ranked by Vina docking scores. Ahmed et al. used molecular docking and 100 ns MDS to estimate the relative binding affinities, interactions, and structure-activity-relationships of 76 prescription antiviral drugs for RdRp. MDS on the best docking candidates showed that remdesivir, raltegravir, and simeprevir had the best binding free energies, ranging from –32 to –38 kcal.mole. Aouldate et l. reported virtual screening of 50,000 chemical compounds from the CAS Antiviral COVID19 database against RdRP. They used a combination of 3D-similarity searches, docking, and 20 ns MDS to rank and select lead compounds and reported one compound (833463-19-7) that bound well to RdRP. Banerjee et al. used protein modelling and computational docking techniques to investigate the effects of common mutations in RdRp, 3CLpro and PLpro sequences of Indian patients. Two RdRp mutations occurred in the Indian population with prevalence >5% and these were analysed as possible targets for repurposed drugs. Docking using Autodock Vina predicted Elbasvir as the best inhibitor of RdRp in the Indian population, followed by Remdesivir and Methylprednisolone.
Natural products and analogues . Given the large degree of interest in repurposing antiviral drugs for use in treating COVID-19, natural product hits are of greater interest given their relative novelty. The drug with the strongest binding affinity to RdRP in our study was the natural anti-helminthic product, Ivermectin. Ivermectin and other avermectins and milbemycins are broad spectrum antiparasitic macrocyclic lactones derived from the bacterium
Streptomyces avermitilis . Ivermectin’s mode of action is by enhancing inhibitory neurotransmission by binding to glutamate-gated chloride channels. Ivermectin has been shown to be effective against several positive-sense single-strand RNA viruses, including SARS-CoV-2,
38, 39 and has been suggested by others as a COVID-19 drug repurposing candidate. It has demonstrated broad in vitro antiviral activity, including against positive-sense single-strand RNA viruses such as SARS-CoV-2. It inhibited replication of SARS-CoV-2 in monkey kidney cell culture with an IC of 2.2 - 2.8 µM.
32, 37, 41
It has been predicted to inhibit SARS-Cov-2 RdRP in several computational studies.
42, 43
Parvez et al. used MDS techniques, augmented by very short DS to predict the binding affinities of Ivermectin and several other antiviral agents included in Table 1 to RdRP. Janabi et al. also predicted favourable binding energies for Ivermectin and several milbemycins to RdRp using AutoDock Vina, but without subsequent simulation of the protein-ligand complexes. Figure 5 summarizes the main interactions between Ivermectin and the binding site of RdRP and the docking pose from the MDS.
Figure 5 . LigPlot (left) and hydrophobic protein surface representation (right) of the main interactions between RdRP and ivermectin. Kalhor et al predicted digoxin would inhibit the interaction of the RBD domain of the SARS-CoV-2 with the ACE2 receptor by but no computation studies other than ours have reported the possibility of RdRP inhibition by digoxin. However, very recently digoxin was shown to have potent in vitro antiviral effects against SARS-Cov-2 in Vero cells (IC
37 nM) and by virus growth kinetics with an IC of 190nM in Vero cells. The activity of digoxin in hese assays was substantially better than that of chloroquine or remdesivir. Figure 6 shows the main interactions of digoxin with the RdRP binding site and how the drug binds to the catalytic cavity. The main interactions are also listed in Supplementary Table 2.
Figure 6 . LigPlot (left) and hydrophobic protein surface representation (right) of the main interactions between RdRP and digoxin. Silibinin has been predicted to be a potential inhibitor of SARS-Cov-2 RdRP by two computational docking studies.
46, 47
Bosch-Barrera et al. speculated that Silibinin may have synergistic benefits for treating COVID-19 patients as, apart from its potent antiviral properties against HCV and HIV, it also combated the cytokine storm that causes major problems in COVID-19 patients. Silibinin is the subject of clinical trials planned for the near future. Figure 7 shows the main interactions of Silibinin with the binding site of RdRP and how it binds to the RdRP catalytic cavity. The main interactions are also listed in Supplementary Table 2 for reference.
Figure 7 . LigPlot (left) and hydrophobic protein surface representation (right) of the main interactions between RdRP and silibinin. AutoDock Vina and MDS using the CHARMM forcefield were used by Pokhrel et al. to predict repurposing of existing drugs and natural products. Sirolimus (Rapamycin) was on the top few hits in their screening and simulation study, suggesting a role for this drug in treating COVID-19 and it is now the subject of a trial of COVID-19 treatment (https://clinicaltrials.gov/ct2/show/NCT04461340). Figure 8 shows the main interactions of sirolimus and the RDRP binding site, and the binding pose from the MDS.
Figure 8 . LigPlot (left) and hydrophobic protein surface representation (right) of the main interactions between RdRP and Sirolimus (rapamycin). Carbetocin is a synthetic analogue of the natural hormone, oxytoxin, in which a labile disulfide bond in the macrocycle is replaced by a thioether. Carbetocin was predicted to be one of the top 10 inhibitors of RdRp (single chain) by Ahmad et al. Their studies showed stable binding with a docking score of -9.5 kcal/mol with 8 strong hydrogen bonding interactions with the active site of the enzyme by the Glide docking-scoring function but no MDS was applied to their hits. Eribulin is a fully synthetic macrocyclic ketone analogue of the marine natural product, Halichondrin B, and is a potent anti-mitotic cancer agent. There are no reports of Eribulin binding to SARS-Cov-2 RdRP but Machitani et al. suggested its COVID-19 use by virtue of its known activity against other viral RdRPs. ovobiocin was also selected as one of the top five best RdRP binders by Choudhury et al. using the MoleGro virtual docker software. . Ergotamine and related ergot alkaloids have been predicted to be high binders to SARS-Cov-2 molecular targets, including the main protease, M pro . There are a few literature reports of its binding to the SARS-Cov-2 RdRP enzyme. There is an in silico predicted SARS-Cov-2 IC of 190µM. Figure 9 shows the main interactions between the RdRP active site residues and ergotamine, and also illustrates the binding pose from the MDS. The interactions are also listed in Supplementary Table 2.
Figure 9 . LigPlot (left) and hydrophobic protein surface representation (right) of the main interactions between ergotamine and RdRP. Computational studies on natural products have been reported recently. Singh et al. reported a computational study of the affinity of plant-derived polyphenols as potential inhibitors of ARS-CoV-2 RdRp. They used AutoDock Vina followed by MDS using Amber to predict the binding poses and affinities of 100 polyphenols. The top four polyphenols had Vina binding energies close to –10 kcal/mol, significantly stronger than the predicted binding of remdesivir. Wu et al. generated homology models of 18 SARS-Cov-2 viral proteins and two human targets and used structure-based virtual ligand screening to identify potential inhibitors of these targets. They identified a range of drugs and natural products that their computational methods suggested may inhibit the RdRP, including beclabuvir.
Other drugs
Bemcentinib selectively inhibits AXL kinase activity, which blocks viral entry and enhances the antiviral type I interferon response. It has been identified as a potential inhibitor of the SARS-Cov-2 M pro enzyme in many computational studies. We have not been able to locate any other published computational publications that suggest that bemcentinib may inhibit the SARS-Cov-2RdRP. Bemcentinib exhibits useful in vitro activity against SARS-Cov-2 with Liu et al. reporting 10-40% protection at 50µM in Vero cells. It was also reported to exhibit an IC of 100nM and CC of 4.7µM in human Huh7.5 cells and an IC of 470nM and CC of 1.6µM in Vero cells, Figure 10 . LigPlot (left) and hydrophobic protein surface representation (right) of the main interactions between RdRP and bemcentinib. Computational docking studies identified an unusual arsenic drug, Darinaparsin, as being an active inhibitor of SARS-CoV-2 RdRp. These researchers performed docking calculations on 12 arsenical drugs using iGEMDOCK (http://gemdock.life.nctu.edu.tw/dock/igemdock.php) and AutoDockVina. They screened the arsenical drugs against several other SARS-Cov-2 drugs and additional several promising leads. Neither of these computational studies used MDS to simulate the highest scoring docked structures and calculate more accurate binding energies.
Other novel putative RdRP inhibitors from the short list of 80 drugs part from the drugs discussed above, several other drugs in the Supplementary Table 1 are of interest. The top 80 list is strongly populated by antiviral drugs such as ciluprevir, glecaprevir, indinavir, simeprevir, elbasvir and ruzasvir. There are several kinase inhibitors with good predicted binding affinities to RdRP including imatinib, ponatinib, rebastinib, lonafarnib, tivantinib and entrectinib. Antibiotics also feature in the list of the 80 best binders e.g. quinupristin, dalfopristin, rifapentin and erythromycin. Other drugs with binding energies lower than –25 kcal/mol include hesperidin, eltrombopag (also active against Mpro), dutasteride, etoposide, quarfloxin, epirubicin, telmisartan, brequinar, conivaptan, rifapentine, sertindole, itraconazole, vapreotide, bafilomycin A1, bromocriptine, idarubicin, midostaurin and rutin. As Supplementary Table 1 shows, a substantial percentage of the hits generated by our docking and MDS protocols have already been shown to have in vitro activity against SARS-Cov-2, or have been predicted to bind to the SARS-CoV-2 RdRP. Interesting some of the compounds have also been predicted to bind several other SARS-Cov-2 protein targets suggesting useful multimodal action that may predict a more potent anti-viral effect. Some identified drugs have already been tested in patients or are undergoing clinical trials in COVID-19. For example, ciluprevir inhibits SARS-Cov-2 Mpro with an IC of 21 µM, indinavir inhibits SARS-Cov-2 with an EC >10 µM and CC >50 µM (A549-hACE2 cells) and an EC of 59 μM and CC >81 μM in Vero cells. Simeprevir inhibits SARS-Cov-2 in vitro with an EC of 4 µM and CC
19 µM (Vero cells), and exhibits SARS-Cov-2 Mpro inhibition with an IC = 10 µM. Eltrombopag inhibits SARS-Cov-2 in vitro with an IC of 8 µM and CC >50 µM (Vero cells), and an IC of 8 µM in Vero and Calu-3 cells. Elbasvir also inhibits SARS-Cov-2 in vitro with an EC of 23 µM in Huh7-hACE2 cells. Telmisartan is in clinical trials for COVID-19 treatment (NCT04356495). Similarly, brequinar has an experimental SARS-Cov-2 EC of 0.3 µM and CC > 50 µM in Vero E6 cells, while imatinib produces an 80% reduction in Mpro activity at 10 µM, with an EC = µM (A549 cells), and an experimental SARS-Cov-2 IC of 3-5µM, with a CC >30 µM, Conivaptan also exhibits an IC of 10 µM against SARS-COV-2 in Vero cells, of ~ 10µM against SARS-Cov-2 Mpro in 293T cells. Ponatinib exhibits an experimental SARS-Cov-2 EC of 1 µM and a CC of 9 µM in HEK-293T cells, and grazoprevir also shows experimental SARS-CoV-2 inhibition, with EC of 16 μM and CC of >100 μM in Vero E6 cells. Finally, itraconazole has an experimental SARS-Cov-2 Mpro IC of 110 µM and SARS-Cov-2 EC = 2.3 μM in human Caco-2 cells, idarubicin shows weak in vitro Mpro activity with an IC of 250−600µM, and ciclesonide shows an in vitro EC for SARS-CoV-2 of 5µM in Vero cells, 0.55µM in differentiated human bronchial tracheal epithelial cells, blocks viral RNA replication, and supresses replication of 15 mutants by >90%
62, 73
It is used to treat COVID-19 patients. As was demonstrated when we applied the same protocol to repurposing of drugs and natural products acting against the SARS-Cov-2 main protease, Mpro, our computational methods generate short targeted lists of existing drugs that may be useful for treating COVID-19. s As well as direct antiviral activities, some of the drugs on the hit list may also have synergistic effects on cytokine storm, clotting disorders, secondary bacterial infections, or pulmonary function. These drugs and natural products merit assessment in SAR-Cov-2 assays and RdRP inhibition experiments. Conclusions
Our virtual screening approach which applied Autodock Vina and MDS in tandem to calculate binding energies for repurposed drugs against the RNA-dependent RNA polymerase (RdRP) identified 80 promising compounds for treating SARS-Cov2 infections. The top hits from our study consisted of a mixture of antiviral agents, natural products and drugs with diverse modes f action. The prognostic value of our computational approach was been demonstrated by our earlier studies of drugs and natural products against the viral main protease, mPro. Here we show that the same protocol generates useful predictions of agents against SARS-Cov-2 replication because a substantial number of the diverse range of drugs in the 80 compound short list exhibit useful SARS-Cov-2 antiviral effects in vitro or have been identified in other computational studies on RdRP. The antiviral drugs simeprevir, sofosbuvir, lopinavir, ritonavir and remdesivir exhibit strong antiviral properties and several in in clinical trial or use against SARS-Cov-2. These drugs have also been identified as binding to RdRP by other virtual screening studies and by in vitro assays. Again, this high validation success rate reinforced the view that our virtual screening protocols were able to identify existing drugs and approved natural products, for rapid testing in the clinic against COVID-19. The hits may also have activity against other coronaviruses. The identified drugs may be useful for treating COVID-19 patients and provide a rational computational basis for repurposing drugs for future pandemics and other diseases.
Materials and Methods
Protein structure preparation and grid preparation
Figure 11 AutoDock Tools (ADT) software was used to prepare the required files for Autodock Vina by assigning hydrogen polarities, calculating Gasteiger charges to protein structures and converting protein structures from the .pdb file format to .pdbqt format. . Screening databases
Drugs were downloaded from the DrugBank database (Wishart et al., 2018) and CHEMBL database (FDA approved) (Gaulton et al., 2017). A total of 8773 and 13,308 drugs were retrieved from Drugbank and CHEMBL database, respectively. The drugs were dowloaded in sdf format and converted to .pdbqt format using Raccoon (Forli et al., 2016).
Docking Methodology mall molecule ligand structures were docked against protein structure using the AutoDock Vina (version 1.1.3) package. AutoDock Vina employs gradient-based conformational search approach and an energy-based empirical scoring function. AutoDock Vina is also flexible, easily scripted, extensively validated in many published studies with a variety of proteins and ligands and takes advantage of large multi-CPU or -GPU machines to run many calculations in parallel. The code has also been employed very successfully to dock millions of small molecule drug candidates into a series of protein targets to discover new potent drug leads. The package includes useful scripts for generating modified .pdb files required for grid calculations and for setting up the grid calculations around each protein automatically. The software requires the removal of hydrogens, addition of polar hydrogens, setting of the correct atom types, and calculation of atom charges compatible with the AutoGrid code. The algorithm generates a grid around each protein and calculates the interaction energy of a probe noble gas atom at each grid position outside and within internal cavities of the protein. The grid resolution was set to 1 Å, the maximum number of binding modes to output was fixed at 10, and the exhaustiveness level (controlling the number of independent runs performed) was set at 8. The docking employed a genetic algorithm to optimize the binding conformations of the ligands during docking to the protease site. Drugs were docked individually to the active site of RdRP (refcode 6M71) with the grid coordinates (grid centre) and grid boxes of appropriate sizes generated by the bash script vina_screen.sh (Supplementary Information). The top scored compounds were identified with a python script script1.py (Supplementary Information) and subjected to molecular dynamic simulation. The docked structures were analysed using UCSF Chimera and LigPlot+ software to illustrate hydrogen-bond and hydrophobic interactions. A total of fifty top compounds selected from each of the Drugbank and CHEMBL compounds. Fifteen compounds were common to both database top hits. Molecular dynamics studies were conducted on the unique set of eighty-five compounds from both sets. olecular Dynamics Simulation . Docked complexes of ligands and COVID-19 M pro protein were used as starting geometries for MD simulations. Simulations were carried out using the GPU accelerated version of the program with the CHARMm force field I periodic boundary conditions in ORACLE server. Docked complexes were immersed in a truncated octahedron box of TIP3P water molecules. The solvated box was further neutralized with Na+ or Cl− counter ions using the tleap program. Particle Mesh Ewald (PME) was employed to calculate the long-range electrostatic interactions. The cut-off distance for the long-range van der Waals (VDW) energy term was 12.0 Å. The whole system was minimized without any restraint. The above steps applied 2500 cycles of steepest descent minimization followed by 5000 cycles of conjugate gradient minimization. After system optimization, the MD simulations was initiated by gradually heating each system in the NVT ensemble from 0 to 300 K for 50 ps using a Langevin thermostat with a coupling coefficient of 1.0/ps and with a force constant of 2.0 kcal/mol·Å2 on the complex. Finally, a production run of 20 ns of MD simulation was performed under a constant temperature of 300 K in the NPT ensemble with periodic boundary conditions for each system. During the MD procedure, the SHAKE algorithm was applied for the constraint of all covalent bonds involving hydrogen atoms. The time step was set to 2 fs. The structural stability of the complex was monitored by the RMSD and RMSF values of the backbone atoms of the entire protein. Calculations were also performed for up to 100 ns on few compounds to ensure that 20ns is sufficiently long for convergence. Duplicate production runs starting with different random seeds were also run to allow estimates of binding energy uncertainties to be determined. The binding free energies of the protein‐protein complexes were evaluated in two ways. The traditional method is to calculate the energies of solvated SARS-Cov-2 protease and small molecule ligands and that of the bound complex and derive the binding energy by subtraction. ΔG (binding, aq) = ΔG (complex, aq) – (ΔG (protein, aq) + ΔG (ligand, aq) (1) We also calculated binding energies using the molecular mechanics Poisson Boltzmann surface area (MM/PBSA) tool in GROMACS that is derived from the nonbonded interaction energies of the complex. The method is also widely used method for binding free energy calculations. MMPBSA calculations were conducted by GMXPBSA 2.1 a suite based on Bash/Perl scripts for streamlining MM/PBSA calculations on structural ensembles derived from GROMACS trajectories and to automatically calculate binding free energies for protein–protein or ligand–protein. GMXPBSA 2.1 calculates diverse MM/PBSA energy contributions from molecular mechanics (MM) and electrostatic contribution to solvation (PB) and non-polar contribution to solvation (SA). This tool combines the capability of MD simulations (GROMACS) and the Poisson–Boltzmann equation (APBS) for calculating solvation energy (Baker et., 2001). The g_mmpbsa tool in GROMACS was used after molecular dynamics simulations, the output files obtained were used to post-process binding free energies by the single-trajectory MMPBSA method. In the current study we considered 100 frames at equal distance from 20ns trajectory files. Specifically, for a non-covalent binding interaction in the aqueous phase the binding free energy, ΔG (bind,aq), is: – ΔG (bind,aqu) = ΔG (bind,vac) + ΔG (bind,solv) (2) where ΔG (bind,vac) is the binding free energy in vacuum, and ΔG(bind,solv) is the solvation free energy change upon binding: – G (bind,solv) = ΔG (R:L, solv) - ΔG (R,solv) - ΔG (L,solv) (3) where ΔG (R:L,solv), ΔG (R,solv) and ΔG (L,solv) are solvation free energies of complex, receptor and ligand, respectively. Note added in proof
As we noted in a previous SARS-Cov-2 drug repurposing paper targeting the main protease, Guterres and Im showed how substantial improvement in protein-ligand docking results could be achieved using high-throughput MD simulations. Over 56 protein targets (of 7 different protein classes) and 560 ligands this showed a 22% improvement in the area under receiver operating characteristics curve, from an initial value of 0.68 using AutoDock Vina alone to a final value of 0.83 when the Vina results were refined by MD.
ACKNOWLEDGEMENTS
We would also like to thank Oracle for providing their Cloud computing resources for the modelling studies described herein and in particular, Peter Winn, Dennis Ward, and Alison Derbenwick-Miller in facilitating these studies. NP, SP and PS are supported by National Institutes of Health Contract HHSN272201400053C. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Author contributions
NP conceived project, analysed data, contributed to manuscript, SP and PS performed the computations, analysed data, contributed to the manuscript and DW analysed the data and contributed to manuscript. eferences
1. Zhang, J.; Zeng, H.; Gu, J.; Li, H.; Zheng, L.; Zou, Q. Progress and Prospects on Vaccine Development against SARS-CoV-2.
Vaccines
8, 153. 2. Whitworth, J. COVID-19: a fast evolving pandemic.
Trans. R. Soc. Trop. Med. Hyg.
Nat. Rev. Drug Discov.
19, 305-306. 4. Sohrabi, C.; Alsafi, Z.; O'Neill, N.; Khan, M.; Kerwan, A.; Al-Jabir, A.; Iosifidis, C.; Agha, R. World Health Organization declares global emergency: A review of the 2019 novel coronavirus (COVID-19).
Int. J. Surg.
76, 71-76. 5. Schlagenhauf, P.; Grobusch, M. P.; Maier, J. D.; Gautret, P. Repurposing antimalarials and other drugs for COVID-19.
Travel Med. Infect. Dis. , 101658. 6. Sanders, J. M.; Monogue, M. L.; Jodlowski, T. Z.; Cutrell, J. B. Pharmacologic Treatments for Coronavirus Disease 2019 (COVID-19): A Review.
J. Am. Med. Assoc.
Vaccines
8. 8. Rosa, S. G. V.; Santos, W. C. Clinical trials on drug repositioning for COVID-19 treatment.
Rev. Panam Salud Publica
44, e40. 9. Olsen, M.; Cook, S. E.; Huang, V.; Pedersen, N.; Murphy, B. G. Perspectives: Potential Therapeutic Options for SARS-CoV-2 Patients Based on Feline Infectious Peritonitis Strategies: Central Nervous System Invasion and Drug Coverage.
Int. J. Antimicrob. Agents , 105964. 0. Mendes, A. Research towards treating COVID-19.
Br. J. Commun. Nurs.
25, 204-205. 11. Lu, S. Timely development of vaccines against SARS-CoV-2.
Emerg. Microbes Infect.
9, 542-544. 12. Jiang, S. Don't rush to deploy COVID-19 vaccines and drugs without sufficient safety guarantees.
Nature
Chemother. , 1-9. 14. Berkley, S. COVID-19 needs a big science approach.
Science
Sci China Life Sci
63, 774-776. 17. Yavuz, S.; Unal, S. Antiviral Treatment of Covid-19.
Turk. J. Med. Sci.
50, 611-619. 18. Zhang, L.; Lin, D.; Sun, X.; Curth, U.; Drosten, C.; Sauerhering, L.; Becker, S.; Rox, K.; Hilgenfeld, R. Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors.
Science
Science
Curr. Res. Chem. Biol. , submitted. 21. Zhu, W.; Chen, C. Z.; Gorshkov, K.; Xu, M.; Lo, D. C.; Zheng, W. RNA-Dependent RNA Polymerase as a Target for COVID-19 Drug Discovery.
SLAS Discov. , 2472555220942123. 22. Huang, J.; Song, W.; Huang, H.; Sun, Q. Pharmacological Therapeutics Targeting RNA-Dependent RNA Polymerase, Proteinase and Spike Protein: From Mechanistic Studies to Clinical Trials for COVID-19.
J. Clin. Med.
9, 1131. 23. Guterres, H.; Im, W. Improving Protein-Ligand Docking Results with High-Throughput Molecular Dynamics Simulations.
J. Chem. Inf. Model.
60, 2189-2198. 24. Winkler, D. A. Ligand Entropy Is Hard but Should Not Be Ignored.
J. Chem. Inf. Model.
60, 4421-4423. 25. Beg, M. A.; Athar, F. Anti-HIV and Anti-HCV drugs are the putative inhibitors of RNA-dependent-RNA polymerase activity of NSP12 of the SARS CoV- 2 (COVID-19). .
Pharm. Pharmacol. Int. J.
Antiviral Res.
Life Sci.
J. Biomol. Struct. Dyn. , 1-9. 31. Wang, M.; Cao, R.; Zhang, L.; Yang, X.; Liu, J.; Xu, M.; Shi, Z.; Hu, Z.; Zhong, W.; Xiao, G. Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro.
Cell Res.
30, 269-271. 32. Anastasiou, I. A.; Eleftheriadou, I.; Tentolouris, A.; Tsilingiris, D.; Tentolouris, N. In Vitro Data of Current Therapies for SARS-CoV-2.
Curr. Med. Chem.
27, 4542-4548. 33. Dyer, O. Covid-19: Remdesivir has little or no impact on survival, WHO trial shows.
Br. Med. J.
J. Biomol. Struct. Dyn. , 1-16. 35. Aouidate, A.; Ghaleb, A.; Chtita, S.; Aarjane, M.; Ousaa, A.; Maghat, H.; Sbai, A.; Choukrad, M.; Bouachrine, M.; Lakhlifi, T. Identification of a novel dual-target scaffold for 3CLpro and RdRp proteins of SARS-CoV-2 using 3D-similarity search, molecular docking, molecular dynamics and ADMET evaluation.
J. Biomol. Struct. Dyn. , 1-14. 6. Banerjee, S.; Dey, R.; Seal, S.; Mondal, K. K.; Bhattacharjee, P. Identification of best suitable repurposed drugs considering mutational spectra at RdRp (nsp12), 3CLpro (nsp 5) and PLpro (nsp 3) of SARS-CoV-2 in Indian population. Research Square 2020, rs.3.rs-33879/v1 37. Heidary, F.; Gharebaghi, R. Ivermectin: a systematic review from antiviral effects to COVID-19 complementary regimen.
J. Antibiot.
73, 593-602. 38. Caly, L.; Druce, J. D.; Catton, M. G.; Jans, D. A.; Wagstaff, K. M. The FDA-approved drug ivermectin inhibits the replication of SARS-CoV-2 in vitro.
Antiviral Res.
Ann. Clin. Microbiol. Antimicrob.
19, 23. 40. Dixit, A.; Yadav, R.; Singh, A. V. Ivermectin: Potential Role as Repurposed Drug for COVID-19.
Malays. J. Med. Sci.
27, 154-158. 41. Simsek Yavuz, S.; Unal, S. Antiviral treatment of COVID-19.
Turk. J. Med. Sci.
50, 611-619. 42. Janabi, D.; Hassan, A. Effective Anti-SARS-CoV-2 RNA Dependent RNA Polymerase Drugs Based on Docking Methods: The Case of Milbemycin, Ivermectin, and Baloxavir Marboxil.
Avicenna J. Med. Biotechnol.
12, 246-250. 43. Parvez, M. S. A.; Karim, M. A.; Hasan, M.; Jaman, J.; Karim, Z.; Tahsin, T.; Hasan, M. N.; Hosen, M. J. Prediction of potential inhibitors for RNA-dependent RNA polymerase of SARS-CoV-2 using comprehensive drug repurposing and molecular docking approach.
Int. J. Biol. Macromolecul. npj Sci Rep
10, 16200. 45. Jeon, S.; Ko, M.; Lee, J.; Choi, I.; Byun, S. Y.; Park, S.; Shum, D.; Kim, S. Identification of Antiviral Drug Candidates against SARS-CoV-2 from FDA-Approved Drugs.
Antimicrobial Agents and Chemotherapy
64, e00819-20. 46. Bosch-Barrera, J.; Martin-Castillo, B.; Buxo, M.; Brunet, J.; Encinar, J. A.; Menendez, J. A. Silibinin and SARS-CoV-2: Dual Targeting of Host Cytokine Storm and Virus Replication Machinery for Clinical Management of COVID-19 Patients.
J. Clin. Med.
9, 1770. 47. Wu, C.; Liu, Y.; Yang, Y.; Zhang, P.; Zhong, W.; Wang, Y.; Wang, Q.; Xu, Y.; Li, M.; Li, X.; Zheng, M.; Chen, L.; Li, H. Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods.
Acta Pharm. Sin. B
10, 766-788. 48. Pokhrel, R.; Chapagain, P.; Siltberg-Liberles, J. Potential RNA-dependent RNA polymerase inhibitors as prospective therapeutics against SARS-CoV-2.
J. Med. Microbiol.
69, 864-873. 49. Ahmad, J.; Ikram, S.; Ahmad, F.; Rehman, I. U.; Mushtaq, M. SARS-CoV-2 RNA Dependent RNA polymerase (RdRp) - A drug repurposing study.
Heliyon
6, e04502. 50. Machitani, M.; Yasukawa, M.; Nakashima, J.; Furuichi, Y.; Masutomi, K. RNA-dependent RNA polymerase, RdRP, a promising therapeutic target for cancer and potentially COVID-19.
Cancer Sci.
Med. J. Arm. Forc. India , in press. 2. Nitulescu, G. M.; Paunescu, H.; Moschos, S. A.; Petrakis, D.; Nitulescu, G.; Ion, G. N. D.; Spandidos, D. A.; Nikolouzakis, T. K.; Drakoulis, N.; Tsatsakis, A. Comprehensive analysis of drugs to treat SARS‑CoV‑2 infection: Mechanistic insights into current COVID‑19 therapies (Review).
Int. J. Mol. Med.
46, 467-488. 53. Chandra, A.; Gurjar, V.; Qamar, I.; Singh, N. Identification of Potential Inhibitors of SARS-COV-2 Endoribonuclease (EndoU) from FDA Approved Drugs: A Drug Repurposing Approach to find Therapeutics for COID19.
J. Biomol. Struct. Dynam. , 1-16. 54. Singh, S.; Sk, M. F.; Sonawane, A.; Kar, P.; Sadhukhan, S. Plant-derived natural polyphenols as potential antiviral drugs against SARS-CoV-2 via RNA-dependent RNA polymerase (RdRp) inhibition: an in-silico analysis.
J. Biomol. Struct. Dyn. , 1-16. 55. Liu, S.; Lien, C. Z.; Selvaraj, P.; Wang, T. T. Evaluation of 19 antiviral drugs against SARS-CoV-2 Infection. bioRxiv 2020, 2020.04.29.067983 56. Dittmar, M.; Lee, J. S.; Whig, K.; Segrist, E.; Li, M.; Jurado, K.; Samby, K.; Ramage, H.; Schultz, D.; Cherry, S. Drug repurposing screens reveal FDA approved drugs active against SARS-Cov-2. bioRxiv 2020, 2020.06.19.161042 57. Chowdhury, T.; Roymahapatra, G.; Mandal, S. M. In Silico Identification of a Potent Arsenic Based Approved Drug Darinaparsin against SARS-CoV-2: Inhibitor of RNA Dependent RNA polymerase (RdRp) and Essential Proteases.
Infect. Disord. Drug Targ. , in press. 58. Baker, J. D.; Uhrich, R. L.; Kraemer, G. C.; Love, J. E.; Kraemer, B. C. A drug repurposing screen identifies hepatitis C antivirals as inhibitors of the SARS-CoV-2 main protease. bioRxiv 2020, 2020.07.10.197889 59. Xie, X.; Muruato, A. E.; Zhang, X.; Lokugamage, K. G.; Fontes-Garfias, C. R.; Zou, J.; Liu, J.; Ren, P.; Balakrishnan, M.; Cihlar, T.; Tseng, C.-T. K.; Makino, S.; Menachery, V. .; Bilello, J. P.; Shi, P.-Y. A nanoluciferase SARS-CoV-2 for rapid neutralization testing and screening of anti-infective drugs for COVID-19. bioRxiv 2020, 2020.06.22.165712 60. Yamamoto, N.; Matsuyama, S.; Hoshino, T.; Yamamoto, N. Nelfinavir inhibits replication of severe acute respiratory syndrome coronavirus 2 in vitro. bioRxiv 2020, 2020.04.06.026476 61. Lo, H. S.; Hui, K. P. Y.; Lai, H.-M.; Khan, K. S.; Kaur, S.; Li, Z.; Chan, A. K. N.; Cheung, H. H.-Y.; Ng, K. C.; Ho, J. C. W.; Chen, Y. W.; Ma, B.; Cheung, P. M.-H.; Shin, D.; Wang, K.; Wu, K.-P.; Dikic, I.; Liang, P.-H.; Zuo, Z.; Chan, F. K. L.; Hui, D. S. C.; Mok, V. C. T.; Wong, K.-B.; Ko, H.; Aik, W. S.; Chan, M. C. W.; Ng, W.-L. Simeprevir suppresses SARS-CoV-2 replication and synergizes with remdesivir. bioRxiv 2020, 2020.05.26.116020 62. Jeon, S.; Ko, M.; Lee, J.; Choi, I.; Byun, S. Y.; Park, S.; Shum, D.; Kim, S. Identification of Antiviral Drug Candidates against SARS-CoV-2 from FDA-Approved Drugs.
Antimicrob. Agents Chemother.
64, e00819-20. 63. Ko, M.; Jeon, S.; Ryu, W. S.; Kim, S. Comparative analysis of antiviral efficacy of FDA-approved drugs against SARS-CoV-2 in human lung cells.
J. Med. Virol. , in press. 64. Milani, M.; Donalisio, M.; Bonotto, R. M.; Schneider, E.; Arduino, I.; Boni, F.; Lembo, D.; Marcello, A.; Mastrangelo, E. Combined in silico docking and in vitro antiviral testing for drug repurposing identified lurasidone and elbasvir as SARS-CoV-2 and HCoV-OC43 inhibitors. bioRxiv 2020, 2020.11.12.379958 65. Rothlin, R. P.; Vetulli, H. M.; Duarte, M.; Pelorosso, F. G. Telmisartan as tentative angiotensin receptor blocker therapeutic for COVID-19.
Drug Devel. Res.
81, 768-770. 66. Sales-Medina, D. F.; Ferreira, L. R. P.; Romera, L. M. D.; Gonçalves, K. R.; Guido, R. V. C.; Courtemanche, G.; Buckeridge, M. S.; Durigon, É. L.; Moraes, C. B.; Freitas-Junior, L. H. Discovery of clinically approved drugs capable of inhibiting SARS-CoV-2 in vitro infection sing a phenotypic screening strategy and network-analysis to predict their potential to treat covid-19. bioRxiv 2020, 2020.07.09.196337 67. Drayman, N.; Jones, K. A.; Azizi, S.-A.; Froggatt, H. M.; Tan, K.; Maltseva, N. I.; Chen, S.; Nicolaescu, V.; Dvorkin, S.; Furlong, K.; Kathayat, R. S.; Firpo, M. R.; Mastrodomenico, V.; Bruce, E. A.; Schmidt, M. M.; Jedrzejczak, R.; Muñoz-Alía, M. Á.; Schuster, B.; Nair, V.; Botten, J. W.; Brooke, C. B.; Baker, S. C.; Mounce, B. C.; Heaton, N. S.; Dickinson, B. C.; Jaochimiak, A.; Randall, G.; Tay, S. Drug repurposing screen identifies masitinib as a 3CLpro inhibitor that blocks replication of SARS-CoV-2 in vitro. bioRxiv 2020, 2020.08.31.274639 68. Weston, S.; Haupt, R.; Logue, J.; Matthews, K.; Frieman, M. B. FDA approved drugs with broad anti-coronaviral activity inhibit SARS-CoV-2 in vitro. bioRxiv 2020, 2020.03.25.008482 69. Xiao, X.; Wang, C.; Chang, D.; Wang, Y.; Dong, X.; Jiao, T.; Zhao, Z.; Ren, L.; Dela Cruz, C. S.; Sharma, L.; Lei, X.; Wang, J. Identification of potent and safe antiviral therapeutic candidates against SARS-CoV-2. bioRxiv 2020, 2020.07.06.188953 70. Vatansever, E. C.; Yang, K.; Kratch, K. C.; Drelich, A.; Cho, C.-C.; Mellot, D. M.; Xu, S.; Tseng, C.-T. K.; Liu, W. R. Targeting the SARS-CoV-2 Main Protease to Repurpose Drugs for COVID-19. bioRxiv 2020, 2020.05.23.112235 71. Van Damme, E.; De Meyer, S.; Bojkova, D.; Ciesek, S.; Cinatl, J.; De Jonghe, S.; Jochmans, D.; Leyssen, P.; Buyck, C.; Neyts, J.; Van Loock, M. In Vitro Activity of Itraconazole Against SARS-CoV-2. bioRxiv 2020, 2020.11.13.381194 72. Ghahremanpour, M. M.; Tirado-Rives, J.; Deshmukh, M.; Ippolito, J. A.; Zhang, C.-H.; Cabeza de Vaca, I.; Liosi, M.-E.; Anderson, K. S.; Jorgensen, W. L. Identification of 14 Known Drugs as Inhibitors of the Main Protease of SARS-CoV-2.
ACS Med. Chem. Lett. , in press. 3. Matsuyama, S.; Kawase, M.; Nao, N.; Shirato, K.; Ujike, M.; Kamitani, W.; Shimojima, M.; Fukushi, S. The inhaled steroid ciclesonide blocks SARS-CoV-2 RNA replication by targeting the viral replication-transcription complex in cultured cells.
J. Virol. , JVI.01648-20. 74. Nakajima, K.; Ogawa, F.; Sakai, K.; Uchiyama, M.; Oyama, Y.; Kato, H.; Takeuchi, I. A Case of Coronavirus Disease 2019 Treated With Ciclesonide.
Mayo Clin. Proc.
95, 1296-1297. 75. Pettersen, E. F.; Goddard, T. D.; Huang, C. C.; Couch, G. S.; Greenblatt, D. M.; Meng, E. C.; Ferrin, T. E. UCSF Chimera--a visualization system for exploratory research and analysis.
J. Comput. Chem.
25, 1605-12. 76. Forli, S.; Huey, R.; Pique, M. E.; Sanner, M. F.; Goodsell, D. S.; Olson, A. J. Computational protein-ligand docking and virtual drug screening with the AutoDock suite.
Nat Protoc
11, 905-19. 77. Laskowski, R. A.; Swindells, M. B. LigPlot+: multiple ligand-protein interaction diagrams for drug discovery.
J Chem Inf Model
51, 2778-86. 78. Zoete, V.; Cuendet, M. A.; Grosdidier, A.; Michielin, O. SwissParam: a fast force field generation tool for small organic molecules.
J Comput Chem
32, 2359-68. 79. Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hessa, B.; Lindahlad, E. GROMACS; High performance molecular simulations through multi-level parallelism from laptops to supercomputers.
SoftwareX
Computer Physics Communications omputational screening of repurposed drugs and natural products against SARS-Cov-2 RdRP as potential COVID-19 therapies
Sakshi Piplani , Puneet Singh , Nikolai Petrovsky , David A. Winkler College of Medicine and Public Health, Flinders University, Bedford Park 5046, Australia Vaxine Pty Ltd, 11 Walkley Avenue, Warradale 5046, Australia La Trobe University, Kingsbury Drive, Bundoora 3042, Australia Monash Institute of Pharmaceutical Sciences, Monash University, Parkville 3052, Australia School of Pharmacy, University of Nottingham, Nottingham NG7 2RD. UK CSIRO Data61, Pullenvale 4069, Australia upplementary information
Table S1 . Binding energies and published SARS-Cov-2 data for 80 top ranked small molecule ligands ChemBl (C)) or Drugbank (D) ID Name ΔG
MMPBSA kcal/mol SARS-Cov-2 data 1 C 3391662 Paritaprevir -54.3 Predicted inhibitor of Mpro and RdRP.
1, 2
2 C 1200633 Ivermectin -54.1 IC of 2.2 - 2.8 µM in monkey kidney cells.
3 C 3126842 Beclabuvir -53.6 4 C 3809489 Bemcentinib -46.2 10-40% protection at 50µM in Vero cells. IC of 100nM and CC of 4.7µM in human Huh7.5 cells and an IC of 470nM and CC was 1.6µM in Vero cells,
5 D 14761 Remdesivir -44.7 In clinical trial for COVD-19, results equivocal,
8, 9 many in vitro reports e.g. IC of 1µM and CC of 275 µM in Vero cells, IC50 of 11.4 µM in Vero cells and IC of 1.3 µM in Calu-3 human lung cells, SARS-Cov-2 RdRP EC = 0.007 μM with IC of 1.7 µM (Vero), 0.3 µM (Calu-3) and 0.01 µM (human airway epithelial cells), computational prediction of RdRP inhibition,
13, 14 and Mpro inhibition.
6 C 1751 Digoxin -41.2 Predicted RdRP inhibitor, IC = 0.043 μM and CC >10µM in Vero cells,
7 D 09298 Silibinin -40.3 Predicted RdRP inhibitor,
17, 18
8 C 1236524 Galidesvir -40.3 Clinical trials for COVID-19 and RdRP inhibitor, predicted Mpro inhibitor,
9 C 1076263 Setrobuvir -40.0 Predicted SARS-Cov-2 RdRP
13, 21 and Mpro inhibitor,
10 C 3707372 Voxilaprevir -39.5 Experimental EC >10 µM and CC
16 µM in A549-hACE2 cells, predicted RdRP and Mpro inhibitor.
11 C 2013174 Vedroprevir -38.5 Predicted Mpro and RdRP inhibitor,
12 C 1241348 Faldaprevir -38.1 Predicted Mpro, PLpro, and RdRP inhibitor,
13 C 413 Sirolimus (Rapamycin) -37.5 Clinical trial for COVID-19,
30, 31 predicted Mpro, PLpro and spike inhibition,
32, 33
14 C 3301668 Carbetocin -37.4 Predicted RdRP and Mpro inhibitor, ChemBl (C)) or Drugbank (D) ID Name ΔG
MMPBSA kcal/mol SARS-Cov-2 data 15 D 01051 Novobiocin -37.3 Predicted RdRP inhibition,
16 C 1683590 Eribulin -36.6 Predicted RdRP and 2′-O-ribose methyltransferase inhibitor.
17 C 442 Ergotamine -36.3 Predicted IC of 190µM, predicted Mpro and RdRP inhibitor.
2, 39
18 C 1957287 Tegobuvir -34.7 SARS-Cov-2 EC = >10 and CC = 18 µM, predicted RdRP and Mpro inhibitor,
19 D 12466 Favipiravir -34.3 Human trials for COVID-19,
9, 43, 44 EC = 62 µM and CC =400 µM (Vero cells), predicted RdRP, helicase and Mpro inhibitor.
21, 46, 47
20 D 14850 Deleobuvir -34.2 Predicted Mpro inhibitor. 21 C 297884 Ciluprevir -34.1 Predicted Mpro inhibitor, SARS-Cov-2 Mpro IC = 21 µM,
22 C 1200649 Quinupristin -33.8 Predicted RdRP inhibitor
23 D 01764 Dalfopristin -33.7 Predicted RdRP inhibitor
24 D 04703 Hesperidin -33.4 Predicted Mpro and RdRP inhibitor
25 C 3545363 Glecaprevir -32.9 Predicted Mpro and RdRP inhibitor
26 D 00224 Indinavir -32.6 Predicted Mpro and RdRP inhibitor, SARS-Cov-2 EC >10 µM and CC >50 µM (A549-hACE2 cells) In vitro EC : 59 μM; CC >81 μM (Vero cells),
27 D 06290 Simeprevir -31.7 In vitro SARS-Cov-2 EC = 4 µM and CC
19 µM (Vero cells), in Vitro Mpro inhibition IC = 10 µM, negligible PLpro and RdRP enzyme inhibition, predicted Mpro and RdRP inhibitor.
2, 56
28 C 461101 Eltrombopag -31.6 Predicted Mpro
2, 39 and RdRP inhibitor,
39, 40 in vitro SARS-Cov-2 IC = 8 µM and CC >50 µM (Vero cells), IC = 8 µM (Vero and Calu-3 cells)
29 D 01126 Dutasteride -31.6 Predicted Mpro and RdRP inhibitor and regulates expression of transmembrane serine protease 2 (TMPRSS2),
30 C 3039514 Elbasvir -31.6 SARS-Cov-2 in vitro EC = 23 µM (Huh7-hACE2 cells), predicted Mpro and RdRP inhibition.
31 C 44657 Etoposide -30.4 Predicted Mpro
60, 61
32 D 06638 Quarfloxin -30.2 Predicted Mpro, human ACE2,
2, 62 and PLpro inhibitor.
33 D 00445 Epirubicin -29.9 Predicted Mpro,
34 C 3833385 Ruzasvir -28.9 Predicted Mpro, RdRP inhibitor
1, 2
ChemBl (C)) or Drugbank (D) ID Name ΔG
MMPBSA kcal/mol SARS-Cov-2 data 35 D 11753 Rifamycin -28.5 Predicted Mpro and RdRP inhibitor.
2, 64
36 C 1017 Telmisartan -28.0 Predicted human ACE2 blocker and clinical trial for COVID-19 treatment (NCT04356495). predicted Mpro and RdRP inhibitor.
66, 67
37 D 03523 Brequinar -28.0 Experimental SARS-Cov-2 EC = 0.3 µM, CC > 50 µM (Vero E6 cells),
38 D 00619 Imatinib -27.7 Experimental 80% reduction in Mpro activity at 10 µM, EC = 8 µM (A549 cells), predicted Mpro inhibitor, experimental SARS-Cov-2 IC = 3-5µM, CC >30 µM,
39 D 00872 Conivaptan -27.6 Predicted Mpro and RdRP inhibitor,
39, 72 experimental SARS-COV-2 IC = 10 µM (Vero cells), EC = 4µM (ACE2-A549 cells), SARS-Cov-2 Mpro IC ~ 10µM (293T cells),
40 D 08901 Ponatinib -27.1 Experimental SARS-Cov-2 EC = 1 µM, CC predicted SARS-Cov-2 spike inhibitor.
41 C 1660 Rifapentine -27.1 Predicted RdRP
64, 75 and Mpro inhibitor.
42 D 06144 Sertindole -26.6 Predicted Mpro inhibition, but negligible experimental SARS-Cov-2 Mpro inhibition.
43 C 2063090 Grazoprevir -26.5 Experimental SARS-CoV-2 inhibition with EC = 16 μM (CC value >100 μM, Vero E6 cells), predicted Mpro and RdRP inhibitor.
24, 40
44 D 01167 Itraconazole -26.4 Predicted Mpro and RdRP inhibition, experimental SARS-Cov-2 Mpro EC = 110 µM. and SARS-Cov-2 EC = 2.3 μM (human Caco-2 cells).
45 C 2103975 Vapreotide -26.2 Predicted Mpro inhibition 46 D 06733 Bafilomycin A1 -25.9 Experimental SARS-Cov-2 EC > 50µM and CC > 50 µM,
47 C 1738757 Rebastinib -25.7 Predicted Mpro and 2’-O-ribose methyltransferase
81, 82 inhibitor. 48 D 06448 Lonafarnib -25.3 Predicted RdRP PLpro and Mpro
25, 63 inhibitor. 49 C 493 Bromocriptine -25.3 Predicted Mrpo and NSP14 inhibitor.
50 D 01177 Idarubicin -25.2 Predicted RdRP and Mpro inhibitor.
66, 83
Weak in vitro Mpro activity IC = 250−600µM.
51 C 608533 Midostaurin -25.0 Predicted Mpro, and spike inhibitor.
52 D 01698 Rutin -25.0 Predicted RdRP and Mpro inhibitor.
53 C 2103882 Tivantinib -24.8 Predicted Mpro and 2 2’-O-methyltransferase inhibitor.
42, 82
ChemBl (C)) or Drugbank (D) ID Name ΔG
MMPBSA kcal/mol SARS-Cov-2 data 54 C 490672 Filibuvir -24.6 Predicted RdRP and Mpro inhibitor.
55 C 1983268 Entrectinib -23.8 Predicted SARS-Cov-2 2’-O-methyltransferase inhibitor.
56 C 1429 Desmopressin -23.7 Predicted SARS-Cov-2 spike, helicase, spike, and 2’- O-methyltransferase inhibitor.
57 C 2106409 Elsamitrucin -22.6 Predicted Mpro inhibitor.
58 D 11618 Zorubicin -22.3 Predicted inhibitor of SARS-Cov-2 spike.
59 C 3039525 Golvatinib -22.2 Predicted SARS-Cov-2 RDRP NSP12‐NSP7 interface, NSP14 SAM-dependent N7-methyl transferase, Mpro, and 2'-O-Ribose Methyltransferase Nsp16
60 C 532 Erythromycin -21.6 Predicted SARS-Cov-2 Mpro inhibitor.
61 D 01092 Ouabain -21.3 Predicted SARS-Cov-2 RdRP and Mpro inhibitor.
15, 42
62 C 3318007 Pimodivir -21.2 Predicted SARS-Cov-2 Mpro inhibitor.
42, 93
63 D 04785 Streptolydigin -20.8 Predicted SARS-Cov-2 RdRP inhibition.
64 D 11616 Pirarubicin -20.3 … 65 C 2104415 Moxidectin -20.2 In vitro IC = 3µM in LLC-MK2 cells, predicted inhibitor of Mpro, RdRP and human ACE2 receptor.
2, 94
66 C 3137309 Venetoclax -20.1 Inactive in SARS‐CoV‐2 CPE assay, predicted Mpro inhibitor
67 D 00549 Zafirlukast -20.0 Predicted inhibitor of SARS-Cov-2 Mpro, spike, and 2'-O-methyltransferases.
2, 37, 91
68 C 4093031 BMS-929075 -19.8 … 69 D 01267 Paliperidone -19.4 Predicted inhibitor of SARS-Cov-2 RdRP and 2'-O-methyltransferases.
39, 82
70 C 1951095 Eravacycline -18.8 Predicted inhibitor of SARS-CoV-2 Mpro
2, 96
RdRP (NSP12), and human ACE2 receptor.
71 C 372795 Streptomycin -18.4 … 72 D 09280 Lumacaftor -17.8 Predicted to be a SARS-Cov-2 spike protein, Mpro and helicase inhibitor.
73 C 3989904 Cethromycin -17.0 … 74 D 01254 Dasatinib -16.9 Active against SARS and MERS at low µM in vitro. Predicted to bind to Mpro.
75 D 01410 Ciclesonide -16.7 In vitro EC for SARS-CoV-2 of 5µM in Vero cells and 0.55µM in differentiated human bronchial tracheal epithelial cells, blocks viral RNA replication, supresses replication of 15 mutants by >90%
57, 101
Used to treat COVID-19 patients.
ChemBl (C)) or Drugbank (D) ID Name ΔG
MMPBSA kcal/mol SARS-Cov-2 data 76 C 4297453 Zalypsis -16.0 … 77 C 444172 Zosuquidar -15.7 Predicted to target 2'-O-ribose methyltransferase Nsp16 of SARS-CoV-2 78 C 1471 Aprepitant -15.6 Predicted to bind to SARS-Cov-2 Mpro
79 D 03325 Tyrosyladenylate -13.6 Predicted to bind to SARS-CoV-2 Nsp16 2’-O-MTase
80 C 408 Troglitazone -13.2 Predicted to bind to SARS-Cov-2 spike able S2 . Binding interactions with RdRP binding site for top 10 ranked drugs. No Drug Interacting residues 1 Beclabuvir Arg553, Lys621, Asp618, Tyr619, Pro620, Asp623, Arg624, Ser759, Asp760, Asp761, Lys798, Glu811, Phe812, Ser814, 2 Bemcentinib Tyr455, Arg553, Lys621, Cys622, Asp623, Asp760, Asp761, Lys798, Glu811, Ser814 3 Digoxin Tyr455, Arg553, Trp617, Asp618, Lys621, Cys622, Asp623, Arg624, Asp760, Asp761, Trp800, Glu811, Phe812, Cys813, Ser814. 4 Galidesvir Asp618, Tyr619, Asp760, Asp761, Trp800, Glu811, Phe812, Ser814 5 Ivermectin Arg553, Arg555, Trp617, Asp618, Asp623, Ser682, Thr687, Asp760, Asp761, Glu811, Cys813, Ser814 6 Paritaprivir Tyr455, Arg553, Asp618, Tyr619, Pro620, Lys621, Asp623, Arg624, Asp760, Asp761, Lys798, Ser814 7 Remdesivir Tyr455, Arg553, Tyr619, Trp617, Lys621, Asp623, Cys622, Arg624, Asp760, Asp761, Lys798, Trp800, Glu811 8 Setrobuvir Arg553, Trp617, Asp618, Lys621, Pro620, Asp623, Asp760, Asp761, Glu811 9 Silibinin Arg553, Trp617, Asp618, Tyr619, Asp623, Thr680, Thr687, Asp760, Asp761 10 Voxilaprevir Lys551, Arg553, Arg555, Asp618, Tyr619, Lys621, Asp623, Asp760, Asp761, Lys798, Glu811, Cys813, Ser814 cripts : receptor = 6Y2F.pdbqt center_x= 9.245 center_y= -0.788 center_z = 18.371 size_x = 50 size_y = 50 size_z = 50 num_modes = 10 exhaustiveness = 50 ! /usr/bin/env python import sys import glob def doit(n): file_names = glob.glob('*/*.pdbqt') everything = [] failures = [] print 'Found', len(file_names), 'pdbqt files' for file_name in file_names: file = open(file_name) lines = file.readlines() file.close() try: line = lines[1] result = float(line.split(':')[1].split()[0]) everything.append([result, file_name]) except: failures.append(file_name) everything.sort(lambda x,y: cmp(x[0], y[0])) part = everything[:n] for p in part: print p[1], print if len(failures) > 0: print 'WARNING:', len(failures), 'pdbqt files could not be processed' if __name__ == '__main__': doit(int(sys.argv[1])) eferences
1. Cozac, R.; Medzhidov, N.; Yuki, S. Predicting inhibitors for SARS-CoV-2 RNA-dependent RNA polymerase using machine learning and virtual screening}. arXiv 2020, 2006.06523 2. Piplani, S.; Singh, P.; Petrovsky, N.; Winkler, D. A. Computational screening of repurposed drugs and natural products against SARS-Cov-2 main protease as potential COVID-19 therapies.
Curr. Res. Chem. Biol. , submitted. 3. Heidary, F.; Gharebaghi, R. Ivermectin: a systematic review from antiviral effects to COVID-19 complementary regimen.
J. Antibiot.
73, 593-602. 4. Simsek Yavuz, S.; Unal, S. Antiviral treatment of COVID-19.
Turk. J. Med. Sci.
50, 611-619. 5. Anastasiou, I. A.; Eleftheriadou, I.; Tentolouris, A.; Tsilingiris, D.; Tentolouris, N. In Vitro Data of Current Therapies for SARS-CoV-2.
Curr. Med. Chem.
27, 4542-4548. 6. Liu, S.; Lien, C. Z.; Selvaraj, P.; Wang, T. T. Evaluation of 19 antiviral drugs against SARS-CoV-2 Infection. bioRxiv 2020, 2020.04.29.067983 7. Dittmar, M.; Lee, J. S.; Whig, K.; Segrist, E.; Li, M.; Jurado, K.; Samby, K.; Ramage, H.; Schultz, D.; Cherry, S. Drug repurposing screens reveal FDA approved drugs active against SARS-Cov-2. bioRxiv 2020, 2020.06.19.161042 8. Dyer, O. Covid-19: Remdesivir has little or no impact on survival, WHO trial shows.
Br. Med. J.
ChemBioChem , in press. 0. Pizzorno, A.; Padey, B.; Dubois, J.; Julien, T.; Traversier, A.; Duliere, V.; Brun, P.; Lina, B.; Rosa-Calatrava, M.; Terrier, O. In vitro evaluation of antiviral activity of single and combined repurposable drugs against SARS-CoV-2.
Antiviral Res.
J. Med. Virol. , in press. 12. Pruijssers, A. J.; George, A. S.; Schafer, A.; Leist, S. R.; Gralinksi, L. E.; Dinnon, K. H., 3rd; Yount, B. L.; Agostini, M. L.; Stevens, L. J.; Chappell, J. D.; Lu, X.; Hughes, T. M.; Gully, K.; Martinez, D. R.; Brown, A. J.; Graham, R. L.; Perry, J. K.; Du Pont, V.; Pitts, J.; Ma, B.; Babusis, D.; Murakami, E.; Feng, J. Y.; Bilello, J. P.; Porter, D. P.; Cihlar, T.; Baric, R. S.; Denison, M. R.; Sheahan, T. P. Remdesivir Inhibits SARS-CoV-2 in Human Lung Cells and Chimeric SARS-CoV Expressing the SARS-CoV-2 RNA Polymerase in Mice.
Cell Rep.
32, 107940. 13. Elfiky, A. A. Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir against SARS-CoV-2 RNA dependent RNA polymerase (RdRp): A molecular docking study.
Life Sci.
Med. J. Arm. Forc. India npj Sci Rep
10, 16200. 17. Bosch-Barrera, J.; Martin-Castillo, B.; Buxo, M.; Brunet, J.; Encinar, J. A.; Menendez, J. A. Silibinin and SARS-CoV-2: Dual Targeting of Host Cytokine Storm and Virus Replication Machinery for Clinical Management of COVID-19 Patients.
J. Clin. Med.
9, 1770. 18. Liu, C.; Zhu, X.; Lu, Y.; Zhang, X.; Jia, X.; Yang, T. Potential Treatment of Chinese and Western Medicine Targeting Nsp14 of SARS-CoV-2.
J. Pharm. Anal. , in press. 19. Gil, C.; Ginex, T.; Maestro, I.; Nozal, V.; Barrado-Gil, L.; Cuesta-Geijo, M. A.; Urquiza, J.; Ramirez, D.; Alonso, C.; Campillo, N. E.; Martinez, A. COVID-19: Drug Targets and Potential Treatments.
J. Med. Chem.
63, 12359-12386. 20. Kumar, S.; Sharma, P. P.; Shankar, U.; Kumar, D.; Joshi, S. K.; Pena, L.; Durvasula, R.; Kumar, A.; Kempaiah, P.; Poonam; Rathi, B. Discovery of New Hydroxyethylamine Analogs against 3CLpro Protein Target of SARS-CoV-2: Molecular Docking, Molecular Dynamics Simulation, and Structure-Activity Relationship Studies.
J. Chem. Inf. Mod. , in press. 21. Elfiky, A. A. SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: an in silico perspective.
J. Biomol. Struct. Dyn. , 1-9. 22. Mosquera-Yuqui, F.; Lopez-Guerra, N.; Moncayo-Palacio, E. A. Targeting the 3CLpro and RdRp of SARS-CoV-2 with phytochemicals from medicinal plants of the Andean Region: molecular docking and molecular dynamics simulations.
J. Biomol. Struct. Dyn. , 1-14. 23. Xie, X.; Muruato, A. E.; Zhang, X.; Lokugamage, K. G.; Fontes-Garfias, C. R.; Zou, J.; Liu, J.; Ren, P.; Balakrishnan, M.; Cihlar, T.; Tseng, C.-T. K.; Makino, S.; Menachery, V. D.; ilello, J. P.; Shi, P.-Y. A nanoluciferase SARS-CoV-2 for rapid neutralization testing and screening of anti-infective drugs for COVID-19. bioRxiv 2020, 2020.06.22.165712 24. Indu, P.; Rameshkumar, M. R.; Arunagirinathan, N.; Al-Dhabi, N. A.; Valan Arasu, M.; Ignacimuthu, S. Raltegravir, Indinavir, Tipranavir, Dolutegravir, and Etravirine against main protease and RNA-dependent RNA polymerase of SARS-CoV-2: A molecular docking and drug repurposing approach.
J. Infect. Pub. Health , in press. 25. Ray, A. K.; Gupta, P. S. S.; Panda, S. K.; Biswal, S.; Rana, M. K. Repurposing of FDA Approved Drugs for the Identification of Potential Inhibitors of SARS-CoV-2 Main Protease. ChemRxiv 2020, chemrxiv.12278066.v1 26. Shah, B.; Modi, P.; Sagar, S. R. In silico studies on therapeutic agents for COVID-19: Drug repurposing approach.
Life Sci.
Molecules
Int. J. Mol. Med.
46, 467-488. 1. Wu, R.; Wang, L.; Kuo, H. D.; Shannar, A.; Peter, R.; Chou, P. J.; Li, S.; Hudlikar, R.; Liu, X.; Liu, Z.; Poiani, G. J.; Amorosa, L.; Brunetti, L.; Kong, A. N. An Update on Current Therapeutic Drugs Treating COVID-19.
Curr. Pharmacol. Rep. , 1-15. 32. Mall, R.; Elbasir, A.; Al Meer, H.; Chawla, S.; Ullah, E. Data-driven Drug Repurposing for COVID-19. ChemRxiv chemrxiv.12661103.v1 33. Pokhrel, R.; Chapagain, P.; Siltberg-Liberles, J. Potential RNA-dependent RNA polymerase inhibitors as prospective therapeutics against SARS-CoV-2.
J. Med. Microbiol.
69, 864-873. 34. Ahmad, J.; Ikram, S.; Ahmad, F.; Rehman, I. U.; Mushtaq, M. SARS-CoV-2 RNA Dependent RNA polymerase (RdRp) - A drug repurposing study.
Heliyon
6, e04502. 35. Farag, A.; Wang, P.; Boys, I. N.; Eitson, J. L.; Ohlson, M. B.; Fan, W.; McDougal, M. B.; Ahmed, M.; Schoggins, J. W.; Sadek, H. Identification of Atovaquone, Ouabain and Mebendazole as FDA Approved Drugs Targeting SARS-CoV-2. ChemRxiv 2020, /chemrxiv.12003930.v4 36. Hosseini, M.; Chen, W.; Wang, C. Computational Molecular Docking and Virtual Screening Revealed Promising SARS-CoV-2 Drugs. . ChemRxiv 2020, chemrxiv.12237995.v1 37. Sharma, K.; Morla, S.; Goyal, A.; Kumar, S. Computational guided drug repurposing for targeting 2'-O-ribose methyltransferase of SARS-CoV-2.
Life Sci.
J. Biomol. Struct. Dynam. , 1-16. 39. Gul, S.; Ozcan, O.; Asar, S.; Okyar, A.; Baris, I.; Kavakli, I. H. In silico identification of widely used and well-tolerated drugs as potential SARS-CoV-2 3C-like protease and viral RNA-ependent RNA polymerase inhibitors for direct use in clinical trials.
J. Biomol. Struct. Dynam. , in press. 40. Ruan, Z.; Liu, C.; Guo, Y.; He, Z.; Huang, X.; Jia, X.; Yang, T. SARS-CoV-2 and SARS-CoV: Virtual screening of potential inhibitors targeting RNA-dependent RNA polymerase activity (NSP12).
J. Med. Virol.
Curr. Med. Chem.
27, 4536-4541. 44. Cai, Q.; Yang, M.; Liu, D.; Chen, J.; Shu, D.; Xia, J.; Liao, X.; Gu, Y.; Cai, Q.; Yang, Y.; Shen, C.; Li, X.; Peng, L.; Huang, D.; Zhang, J.; Zhang, S.; Wang, F.; Liu, J.; Chen, L.; Chen, S.; Wang, Z.; Zhang, Z.; Cao, R.; Zhong, W.; Liu, Y.; Liu, L. Experimental Treatment with Favipiravir for COVID-19: An Open-Label Control Study.
Eng. , in press. 45. Wang, M.; Cao, R.; Zhang, L.; Yang, X.; Liu, J.; Xu, M.; Shi, Z.; Hu, Z.; Zhong, W.; Xiao, G. Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro.
Cell Res.
30, 269-271. 6. Borgio, J. F.; Alsuwat, H. S.; Al Otaibi, W. M.; Ibrahim, A. M.; Almandil, N. B.; Al Asoom, L. I.; Salahuddin, M.; Kamaraj, B.; AbdulAzeez, S. State-of-the-art tools unveil potent drug targets amongst clinically approved drugs to inhibit helicase in SARS-CoV-2.
Arch. Med. Sci.
16, 508-518. 47. Sencanski, M.; Perovic, V.; Pajovic, S.; Adzic, M.; Paessler, S.; Glisic, S. Drug Repurposing for Candidate SARS-CoV-2 Main Protease Inhibitors by a Novel in Silico Method.
Molecules
25, 3830. 48. Chakraborti, S.; Bheemireddy, S.; Srinivasan, N. Repurposing drugs against the main protease of SARS-CoV-2: mechanism-based insights supported by available laboratory and clinical data.
Molecular omics . 49. Baker, J. D.; Uhrich, R. L.; Kraemer, G. C.; Love, J. E.; Kraemer, B. C. A drug repurposing screen identifies hepatitis C antivirals as inhibitors of the SARS-CoV-2 main protease. bioRxiv 2020, 2020.07.10.197889 50. Singh, S.; Sk, M. F.; Sonawane, A.; Kar, P.; Sadhukhan, S. Plant-derived natural polyphenols as potential antiviral drugs against SARS-CoV-2 via RNA-dependent RNA polymerase (RdRp) inhibition: an in-silico analysis.
J. Biomol. Struct. Dyn. , 1-16. 51. Das, S.; Sarmah, S.; Lyndem, S.; Singha Roy, A. An investigation into the identification of potential inhibitors of SARS-CoV-2 main protease using molecular docking study.
J. Biomol. Struct. Dynam. . 52. Joshi, R. S.; Jagdale, S. S.; Bansode, S. B.; Shankar, S. S.; Tellis, M. B.; Pandya, V. K.; Chugh, A.; Giri, A. P.; Kulkarni, M. J. Discovery of potential multi-target-directed ligands by targeting host-specific SARS-CoV-2 structurally conserved main protease.
J. Biomol. Struct. Dynam. , 1-16. 3. Fischer, A.; Sellner, M.; Neranjan, S.; Smieško, M.; Lill, M. A. Potential Inhibitors for Novel Coronavirus Protease Identified by Virtual Screening of 606 Million Compounds.
Int. J. Mol. Sci.
21, 3626. 54. Yamamoto, N.; Matsuyama, S.; Hoshino, T.; Yamamoto, N. Nelfinavir inhibits replication of severe acute respiratory syndrome coronavirus 2 in vitro. bioRxiv 2020, 2020.04.06.026476 55. Lo, H. S.; Hui, K. P. Y.; Lai, H.-M.; Khan, K. S.; Kaur, S.; Li, Z.; Chan, A. K. N.; Cheung, H. H.-Y.; Ng, K. C.; Ho, J. C. W.; Chen, Y. W.; Ma, B.; Cheung, P. M.-H.; Shin, D.; Wang, K.; Wu, K.-P.; Dikic, I.; Liang, P.-H.; Zuo, Z.; Chan, F. K. L.; Hui, D. S. C.; Mok, V. C. T.; Wong, K.-B.; Ko, H.; Aik, W. S.; Chan, M. C. W.; Ng, W.-L. Simeprevir suppresses SARS-CoV-2 replication and synergizes with remdesivir. bioRxiv 2020, 2020.05.26.116020 56. Ahmed, S.; Mahtarin, R.; Ahmed, S. S.; Akter, S.; Islam, M. S.; Mamun, A. A.; Islam, R.; Hossain, M. N.; Ali, M. A.; Sultana, M. U. C.; Parves, M. R.; Ullah, M. O.; Halim, M. A. Investigating the binding affinity, interaction, and structure-activity-relationship of 76 prescription antiviral drugs targeting RdRp and Mpro of SARS-CoV-2.
J. Biomol. Struct. Dyn. , 1-16. 57. Jeon, S.; Ko, M.; Lee, J.; Choi, I.; Byun, S. Y.; Park, S.; Shum, D.; Kim, S. Identification of Antiviral Drug Candidates against SARS-CoV-2 from FDA-Approved Drugs.
Antimicrob. Agents Chemother.
64. 58. Milani, M.; Donalisio, M.; Bonotto, R. M.; Schneider, E.; Arduino, I.; Boni, F.; Lembo, D.; Marcello, A.; Mastrangelo, E. Combined in silico docking and in vitro antiviral testing for drug repurposing identified lurasidone and elbasvir as SARS-CoV-2 and HCoV-OC43 inhibitors. bioRxiv 2020, 2020.11.12.379958 9. Beg, M. A.; Athar, F. Anti-HIV and Anti-HCV drugs are the putative inhibitors of RNA-dependent-RNA polymerase activity of NSP12 of the SARS CoV- 2 (COVID-19). .
Pharm. Pharmacol. Int. J.
8, 163‒172. 60. Chen, Y. W.; Yiu, C.-P. B.; Wong, K.-Y. Prediction of the SARS-CoV-2 (2019-nCoV) 3C-like protease (3CL pro) structure: virtual screening reveals velpatasvir, ledipasvir, and other drug repurposing candidates. F1000Research 2020, f1000research.22457.2 61. Anwar, M. U.; Adnan, F.; Abro, A.; Khan, M. R.; Rehman, A. U.; Osama, M.; Javed, S.; Baig, A.; Shabbir, M. R.; Assir, M. Z. Combined Deep Learning and Molecular Docking Simulations Approach Identifies Potentially Effective FDA Approved Drugs for Repurposing Against SARS-CoV-2. ChemRxiv 2020, chemrxiv.12227363.v1 62. Alexpandi, R.; De Mesquita, J. F.; Karutha Pandian, S. K.; Ravi, A. V. Quinolines-Based SARS-CoV-2 3CLpro and RdRp Inhibitors and Spike-RBD-ACE2 Inhibitor for Drug-Repurposing Against COVID-19: An in silico Analysis.
Front. Microbiol.
11, 1796. 63. Murugan, N. A.; Kumar, S.; Jeyakanthan, J.; Srivastava, V. Searching for target-specific and multi-targeting organics for Covid-19 in the Drugbank database with a double scoring approach. npj Sci. Rep.
10, 19125. 64. Elkarhat, Z.; Charoute, H.; Elkhattabi, L.; Barakat, A.; Rouba, H. Potential inhibitors of SARS-cov-2 RNA dependent RNA polymerase protein: molecular docking, molecular dynamics simulations and MM-PBSA analyses.
J. Biomol. Struct. Dynam. , 1-14. 65. Rothlin, R. P.; Vetulli, H. M.; Duarte, M.; Pelorosso, F. G. Telmisartan as tentative angiotensin receptor blocker therapeutic for COVID-19.
Drug Devel. Res.
81, 768-770. 6. Jimenez-Alberto, A.; Ribas-Aparicio, R. M.; Aparicio-Ozores, G.; Castelan-Vega, J. A. Virtual screening of approved drugs as potential SARS-CoV-2 main protease inhibitors.
Comp. Biol. Chem.
88, 107325-107325. 67. Ivanov, J.; Polshakov, D.; Kato-Weinstein, J.; Zhou, Q.; Li, Y.; Granet, R.; Garner, L.; Deng, Y.; Liu, C.; Albaiu, D.; Wilson, J.; Aultman, C. Quantitative Structure–Activity Relationship Machine Learning Models and their Applications for Identifying Viral 3CLpro- and RdRp-Targeting Compounds as Potential Therapeutics for COVID-19 and Related Viral Infections.
ACS Omega
5, 27344-27358. 68. Sales-Medina, D. F.; Ferreira, L. R. P.; Romera, L. M. D.; Gonçalves, K. R.; Guido, R. V. C.; Courtemanche, G.; Buckeridge, M. S.; Durigon, É. L.; Moraes, C. B.; Freitas-Junior, L. H. Discovery of clinically approved drugs capable of inhibiting SARS-CoV-2 in vitro infection using a phenotypic screening strategy and network-analysis to predict their potential to treat covid-19. bioRxiv 2020, 2020.07.09.196337 69. Drayman, N.; Jones, K. A.; Azizi, S.-A.; Froggatt, H. M.; Tan, K.; Maltseva, N. I.; Chen, S.; Nicolaescu, V.; Dvorkin, S.; Furlong, K.; Kathayat, R. S.; Firpo, M. R.; Mastrodomenico, V.; Bruce, E. A.; Schmidt, M. M.; Jedrzejczak, R.; Muñoz-Alía, M. Á.; Schuster, B.; Nair, V.; Botten, J. W.; Brooke, C. B.; Baker, S. C.; Mounce, B. C.; Heaton, N. S.; Dickinson, B. C.; Jaochimiak, A.; Randall, G.; Tay, S. Drug repurposing screen identifies masitinib as a 3CLpro inhibitor that blocks replication of SARS-CoV-2 in vitro. bioRxiv 2020, 2020.08.31.274639 70. Kumar, D.; Chandel, V.; Raj, S.; Rathi, B. In silico identification of potent FDA approved drugs against Coronavirus COVID-19 main protease: A drug repurposing approach.
Int. J. Biol. Macromolecul.
J. Biomol. Struct. Dynam. , 1-19. 77. Vatansever, E. C.; Yang, K.; Kratch, K. C.; Drelich, A.; Cho, C.-C.; Mellot, D. M.; Xu, S.; Tseng, C.-T. K.; Liu, W. R. Targeting the SARS-CoV-2 Main Protease to Repurpose Drugs for COVID-19. bioRxiv 2020, 2020.05.23.112235 8. Hakmi, M.; Bouricha, E.; Kandoussi, I.; El Harti, J.; Ibrahimi, A. Repurposing of known anti-virals as potential inhibitors for SARS-CoV-2 main protease using molecular docking analysis.
Bioinformat.
16, 301-305. 79. Van Damme, E.; De Meyer, S.; Bojkova, D.; Ciesek, S.; Cinatl, J.; De Jonghe, S.; Jochmans, D.; Leyssen, P.; Buyck, C.; Neyts, J.; Van Loock, M. In Vitro Activity of Itraconazole Against SARS-CoV-2. bioRxiv 2020, 2020.11.13.381194 80. Liu, X.; Wang, X. J. Potential inhibitors against 2019-nCoV coronavirus M protease from clinically approved medicines.
J. Genet. Genom.
47, 119-121. 81. Jiang, Y.; Liu, L.; Manning, M.; Bonahoom, M.; Lotvola, A.; Yang, Z.; Yang, Z.-Q. Structural analysis, virtual screening and molecular simulation to identify potential inhibitors targeting 2'-O-ribose methyltransferase of SARS-CoV-2 coronavirus.
J. Biomol. Struct. Dynam. , 1-16. 82. Jiang, Y.; Liu, L.; Manning, M.; Bonahoom, M.; Lotvola, A.; Yang, Z.-Q. Repurposing Therapeutics to Identify Novel Inhibitors Targeting 2'-O-Ribose Methyltransferase Nsp16 of SARS-CoV-2. ChemRxiv 2020, chemrxiv.12252965.v1 83. Wu, C.; Liu, Y.; Yang, Y.; Zhang, P.; Zhong, W.; Wang, Y.; Wang, Q.; Xu, Y.; Li, M.; Li, X.; Zheng, M.; Chen, L.; Li, H. Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods.
Acta Pharm. Sin. B
10, 766-788. 84. Ghahremanpour, M. M.; Tirado-Rives, J.; Deshmukh, M.; Ippolito, J. A.; Zhang, C.-H.; Cabeza de Vaca, I.; Liosi, M.-E.; Anderson, K. S.; Jorgensen, W. L. Identification of 14 Known Drugs as Inhibitors of the Main Protease of SARS-CoV-2.
ACS Med. Chem. Lett. , in press. 5. Trezza, A.; Iovinelli, D.; Santucci, A.; Prischi, F.; Spiga, O. An integrated drug repurposing strategy for the rapid identification of potential SARS-CoV-2 viral inhibitors. npj Sci Rep
10, 13866. 86. da Silva, F. M. A.; da Silva, K. P. A.; de Oliveira, L. P. M.; Costa, E. V.; Koolen, H. H.; Pinheiro, M. L. B.; de Souza, A. Q. L.; de Souza, A. D. L. Flavonoid glycosides and their putative human metabolites as potential inhibitors of the SARS-CoV-2 main protease (Mpro) and RNA-dependent RNA polymerase (RdRp).
Mem. Inst. Oswaldo Cruz
J. Proteome Res.
Nereis , 15-18. 3. Chtita, S.; Belhassan, A.; Aouidate, A.; Belaidi, S.; Bouachrine, M.; Lakhlifi, T. Discovery of Potent SARS-CoV-2 Inhibitors from Approved Antiviral Drugs via Docking Screening.
Comb. Chem. High Throughp. Scr. , in press. 94. Oliveira, M. D.; Oliveira, K. M. Comparative docking of SARS-CoV-2 receptors antagonists from repurposing drugs. ChemRxiv 2020, chemrxiv.12044538.v4 95. Alves, V. M.; Bobrowski, T.; Melo-Filho, C. C.; Korn, D.; Auerbach, S.; Schmitt, C.; Muratov, E. N.; Tropsha, A. QSAR Modeling of SARS-CoV M(pro)Inhibitors Identifies Sufugolix, Cenicriviroc, Proglumetacin, and Other Drugs as Candidates for Repurposing against SARS-CoV-2.
Mol. Info. , in press. 96. Wang, J. Fast Identification of Possible Drug Treatment of Coronavirus Disease-19 (COVID-19) through Computational Drug Repurposing Study.
J. Chem. Inf. Mod.
60, 3277-3286. 97. Beck, B. R.; Shin, B.; Choi, Y.; Park, S.; Kang, K. Predicting commercially available antiviral drugs that may act on the novel coronavirus (SARS-CoV-2) through a drug-target interaction deep learning model.
Comp. Struct. Biotechnol. J.
18, 784-790. 98. White, M. A.; Lin, W.; Cheng, X. Discovery of COVID-19 Inhibitors Targeting the SARS-CoV-2 Nsp13 Helicase.
J. Phys. Chem. Lett. , 9144-9151. 99. Ekins, S.; Mottin, M.; Ramos, P. R.; Sousa, B. K.; Neves, B. J.; Foil, D. H.; Zorn, K. M.; Braga, R. C.; Coffee, M.; Southan, C. Déjà vu: Stimulating open drug discovery for SARS-CoV-2.
Drug Discov. Today
25, 928-941. 100. Qiao, Z.; Zhang, H.; Ji, H. F.; Chen, Q. Computational View toward the Inhibition of SARS-CoV-2 Spike Glycoprotein and the 3CL Protease.
Computat.
8, 53. 01. Matsuyama, S.; Kawase, M.; Nao, N.; Shirato, K.; Ujike, M.; Kamitani, W.; Shimojima, M.; Fukushi, S. The inhaled steroid ciclesonide blocks SARS-CoV-2 RNA replication by targeting the viral replication-transcription complex in cultured cells.
J. Virol. , JVI.01648-20. 102. Nakajima, K.; Ogawa, F.; Sakai, K.; Uchiyama, M.; Oyama, Y.; Kato, H.; Takeuchi, I. A Case of Coronavirus Disease 2019 Treated With Ciclesonide.
Mayo Clin. Proc.2020,