Abstract
We urgently need to identify drugs to treat patients suffering from COVID-19 infection. Drugs rarely act at single molecular targets. Off-target effects are responsible for undesirable side effects and beneficial synergy between targets for specific illnesses. They have provided blockbuster drugs, e.g., Viagra for erectile dysfunction and Minoxidil for male pattern baldness. Existing drugs, those in clinical trials, and approved natural products constitute a rich resource of therapeutic agents that can be quickly repurposed, as they have already been assessed for safety in man. A key question is how to screen such compounds rapidly and efficiently for activity against new pandemic pathogens such as SARS-CoV-2. Here, we show how a fast and robust computational process can be used to screen large libraries of drugs and natural compounds to identify those that may inhibit the main protease of SARS-CoV-2. We show that the shortlist of 84 candidates with the strongest predicted binding affinities is highly enriched (≥25%) in compounds experimentally validated in vivo or in vitro to have activity in SARS-CoV-2. The top candidates also include drugs and natural products not previously identified as having COVID-19 activity, thereby providing leads for experimental validation. This predictive in silico screening pipeline will be valuable for repurposing existing drugs and discovering new drug candidates against other medically important pathogens relevant to future pandemics.
Introduction
The devastating impact of the COVID-19 pandemic caused by SARS coronavirus-2 (SARS-CoV-2) has stimulated unprecedented international activity to discover effective drugs for this and other pathogenic coronaviruses such as SARS and MERS CoV (; Zhang J. et al., 2020; ; Zhang J.-J. et al., 2020; Zhang T. et al., 2020; ; ; ; ; ; ; ; ; ; ; ; ). Computational methods are useful, fast approaches to determine the affinities of small drug-like molecules for SARS-CoV-2 protein targets. Recent papers in Science have reported effective computational de novo drug design based on the structures of the SARS-CoV-2 protease (Zhang L. et al., 2020; ). Clearly, the design of potent new drugs for coronaviruses is very important for future pandemic preparedness, given that the last three serious epidemics have been caused by coronaviruses. However, to make an impact on the current COVID-19 pandemic, given the 10- to 15-year time frame required to take drug leads from lab to clinic, it is only feasible to repurpose drugs that are already registered (off label use), have been through at least phase 1 clinical trials to establish initial human safety, or are approved natural products. Any COVID-19 drug candidates identified in this way can then be used very quickly, as their safety and pharmacokinetics should be already well understood. Drugs that reduce viral replication primarily by targeting viral proteases and polymerases are classified as direct-acting antivirals and are the focus of the current work. Other studies have explored host-targeted drugs that inhibit cellular functions required for viral replication and thereby inhibit SARS-CoV-2 infection, albeit with more potential for host side effects().
The SARS-CoV-2 genome encodes >20 proteins, many of which are potential antiviral drug targets (Figure 1). Two proteases (PLpro and 3CLpro) are essential for virus replication. These enzymes cleave the PP1A and PP1AB polyproteins into functional components. 3-Chymotrypsin-like protease (3CLpro) catalytically self-cleaves a peptide bond between a glutamine at position P1 and a small amino acid (serine, alanine, or glycine) at position P1’. This protease corresponds to non-structural protein 5 (nsp5), the main protease (aka main protease, Mpro) in coronaviruses. 3CL protease is crucial to the processing of the coronavirus replicase polyprotein (P0C6U8), cleaving it at 11 conserved sites. It employs a Cys-His catalytic dyad in its active site, where the cysteine sulfur is the nucleophile, and the histidine imidazole ring acts as a general base. Mpro is a conserved drug target present in all Coronavirinae. It does not have a human homolog, reducing the risk that drugs inhibiting it will exhibit side effects (). Very recent research has shown that strong Mpro inhibitors can substantially reduce SARS-CoV-2 virus titers, reduce weight loss, and improve survival in mice (), making Mpro a promising drug target for structure-based drug discovery.
FIGURE 1
Computational methods can rapidly and efficiently identify candidate drugs for repurposing in pandemic situations where speed is of utmost importance. A very recent paper by Llanos et al. analyzed the strengths and weaknesses of docking simulations for SARS-CoV-2 drug repurposing for Mpro (). This study disclosed that most published studies do not check the ability of the docking method to accurately redock ligands from protein structures and do not account for protein and ligand flexibility using MD calculations, and only a tiny percentage validate predictions using experimental measurements of virus activity. To address these shortcomings, here we used validated molecular docking followed by high-throughput molecular dynamics simulations to prioritize, from an initial large number of licensed or clinical trial drugs and natural products, a short list of the most promising candidates.
Results
Molecular dynamics calculations were used to predict the optimal binding poses and binding energies for 84 of the top hits from docking-based virtual screening of ∼12,000 drug candidates against the SARS-CoV-2 Mpro. The docking protocols were validated by redocking ligands from 10 x-ray structures. The top candidates were ranked for COVID-19 repurposing based on binding affinity and novelty. Conspicuously, we found that ∼30% of the computationally repurposed drug candidates have experimentally validated activity against the Mpro target protein, the SARS-CoV-2 virus, or both. Several of the drugs we identified are currently in clinical trials for COVID-19.
The binding energies of the 84 top ranked ligands from the docking calculations are listed in Supplementary Table S1. Note that calculating accurate absolute binding energies is difficult, and the approach we have taken provides good estimates of the relative binding energies of repurposing candidates. The ten drugs with the tightest binding to Mpro are summarized in Table 1, together with their GMXPBSA binding energies. The binding energies of several of the antiviral drugs, namely, simeprevir, sofosbuvir, lopinavir, and ritonavir, are very similar, within the uncertainties in calculated energies. Some of the antivirals were also identified in other in silico docking studies or wet-lab SARS-CoV-2 activity studies, as we discuss below. This, together with a subsequent extensive search of the literature for experimental data, provides strong validation of the utility of our computational methods to find leads consistent with other studies that also inhibit SARS-CoV-2 or the relevant protein target. It strongly suggests that the computational protocols we have adopted are very capable of generating a list of repurposing candidates, many of which are likely to exhibit useful experimental in vitro activity at least.
TABLE 1
| ID | Structure | Description | ΔGMMPBSA (ΔGbind) (kcal/mol) |
|---|---|---|---|
| C3809489 bemcentinib |  | Inhibitor of the kinase domain of AXL receptor | −34.7 ± 2.6 (−30.7) |
| C4291143 PC786 |  | Respiratory syncytial virus (RSV) L protein polymerase inhibitor | −33.1 ± 0.3 (−29.2) |
| C787 Montelukast |  | Leukotriene receptor antagonist used with cortico-steroids for asthma therapy | −32.7 ± 0.2 (−20.6) |
| C442 Ergotamine |  | Alpha-1 selective adrenergic agonist used in migraine treatment | −31.5 ± 0.3 (−28.7) |
| D06290 simeprevir |  | Hepatitis C virus (HCV) NS3/4A protease inhibitor | −31.4 ± 0.2 (−29.2) |
| D08934 sofosbuvir |  | Nucleotide prodrug and HCV NS5B polymerase inhibitor | −31.0 ± 0.5 (−22.8) |
| D01601 lopinavir |  | Antiretroviral protease inhibitor for treatment of HIV-1 | −30.7 ± 0.3 (−20.4) |
| D00503 ritonavir |  | Peptidomimetic inhibitor of HIV-1 and HIV-2 proteases | −30.5 ± 0.5 (−21.3) |
| C2105887 Mergocriptine |  | Synthetic ergot derivative, dopamine receptor agonist | −30.0 ± 0.3 (−17.9) |
| D14761 remdesivir |  | Viral RNA-dependent RNA polymerase inhibitor | −30.0 ± 0.2 (−27.1) |
Binding energies of 10 top ranked small-molecule ligands for SARS-CoV-2 Mpro.
In this paper, we have focused particularly on tightly binding drugs with novel structures, such as ergot compounds, bemcentinib, PC786, and montelukast.
Although the main focus of the paper is to show that appropriate computational methods can make useful predictions of the repurposing potential of drugs, we also provide a preliminary analysis of the binding of candidate drugs to the Mpro active site. Mpro achieves protein cleavage via the catalytic dyad His41and Cys145. The main active site residues that have previously been implicated in drug binding are His41, Gly143, Cys145, His163, Glu166, and Glu166. All of the drugs whose interactions with the Mpro binding site are summarized below interact with these six residues (see Supplementary Table S2). Most form strong hydrogen binds to one or more of Gly143, Cys145, and His163. All docked and MD simulated structures of Mpro with the repurposed drug candidates were also deposited in open access data archives. Supplementary Figure S2 shows a superimposition of the top 10 drugs bound to the Mpro site.
Ergotamine and mergocriptine, a synthetic long-acting ergot derivative, are α1 selective adrenergic agonist vasoconstrictors and an agonist of dopamine receptors, respectively. Figure 2 shows a LigPlot representation of the interactions of key functional groups in ergotamine and mergocriptine with protease active site residues. These, together with the accompanying Mpro binding site molecular surface plots encoded for lipophilicity, illustrate how these drugs bind in the protease binding site. The specific interactions between these drugs and the residues in the binding site are summarized in Supplementary Table S2. Both drugs make strong and multiple interactions with 20 active site residues, notably hydrogen bonds with Gly143, His164, Met165, Cys145, and Thr190.
FIGURE 2
Montelukast is a cysteinyl leukotriene receptor antagonist used to treat asthma and allergic rhinitis. It reduces pulmonary responses to antigen, tissue eosinophilia and IL-5 expression in inflammatory cells and decreases elevated levels of IL-1β and IL8 in viral upper respiratory tract infections (). Figure 3 shows a LigPlot representation of the interactions of key functional groups in montelukast with protease active site residues and a representation of how this drug binds in the active site of Mpro. The specific interactions between montelukast and the active site residues are also summarized in Supplementary Table S2. The drug interacts extensively with the active site, binding to 23 residues, forming strong hydrogen bonds with Ser144 and Cys145. Montelukast spans the relatively broad binding pocket of the enzyme.
FIGURE 3
Bemcentinib selectively inhibits AXL kinase activity, which blocks viral entry and enhances the antiviral type I interferon response. Figure 4 provides a LigPlot representation of the interactions of key functional groups in bemcentinib with protease active site residues, which are also summarized in detail in Supplementary Table S2. It forms strong hydrogen-bonding interactions with Val186, Arg188, andGln192. Figure 4 also shows the docking pose of bemcentinib in the protease active site after simulation by MD. The hydrophobic benzocycloheptapyridazine moiety occupies a relative hydrophobic pocket, while the hydrophilic triazolyldiamine moiety binds strongly to the polar pocket formed partially by Asp187 and Arg188.
FIGURE 4
PC786 targets the respiratory syncytial virus (RSV) L protein and is designed to be a topical inhalation treatment, a likely route of infection for SARS-CoV-2. Figure 5 shows a LigPlot representation of the interactions of key functional groups in PC786 with protease active site residues, with the specific interactions listed in Supplementary Table S2. It forms a hydrogen bond network with Gly143, Ser144, and Cys145. Figure 5 illustrates the binding pose of PC786 in the Mpro binding site after MD simulations based on the structure obtained from Vina docking calculations. The hydrophobic phenyl ring of the benzazepine moiety projects into a hydrophobic pocket formed partly by Thr25 and Thr26.
FIGURE 5
The predicted binding energies of the 84 drugs in the short list are summarized in Supplementary Table S1, along with details of any experiments to determine their activities against Mpro or SARS-CoV-2 in vitro or in vivo. This suggests that our screening and MD simulation methods are sufficiently robust and accurate to identify drugs for repurposing against SARS-CoV-2 and, more broadly, other coronaviruses. The 33% of drugs in the hit list that have not been reported before are clearly of potential interest as novel drugs for treating COVID-19. We discuss below some of the more interesting and novel hit compounds with stronger binding affinities.
Eltrombopag is a thrombopoietin (TPO) receptor agonist that acts at the transmembrane domain of its cognate receptor C-Mpl via a histidine residue that occurs only in humans and apes. It scored highly in the docking studies, suggesting that it could inhibit the Mpro and exhibit antiviral activity. Figure 6 shows a LigPlot representation of the interactions of key functional groups in eltrombopag with protease active site residues. The binding pose of eltrombopag in the active site of Mpro from the MD simulations is also shown in Figure 6. Close analysis of the binding mode shows that eltrombopag occupies the main part of the Mpro binding pocket, with the hydrophilic biphenyl moiety binding to the hydrophobic pocket formed partly by Cys145. The hydrophilic pyrazolone lies in a polar cleft bounded by Glu189 and Glu166, with the terminal dimethyl phenyl ring undergoing a hydrophobic interaction with Thr190.
FIGURE 6
Eltrombopag is of particular interest as an Mpro inhibitor lead because it is novel and is also a member of a large class of small molecular TPO receptor agonists that may also exhibit activity against the viral protease, and potentially the spike protein and human ACE2.(). However, given the clotting disorders that SAR-CoV-2 generates, the TPOR agonist activities would need to be minimized to prevent platelet enhancement, while retaining or enhancing the antiviral activities.
Apart from the drugs discussed above, several other drugs in Supplementary Table S1 are of interest. There are several other ergot derivatives with good predicted binding affinities to Mpro. Metergotamine and dihydroergocristine were predicted to have ΔGbind of –29 and –24 kcal/mol, respectively.
Discussion
Our virtual screening approach, using Autodock Vina and MD simulation in tandem to calculate binding poses and energies for repurposed drugs, identified 84 compounds with potential for treating COVID-19. The top hits from our study consisted of a mixture of antiviral agents, natural products and drugs developed for other applications and that have additional models of action. We now discuss the results of our computational screening in the context of other computational studies of Mpro in the literature.
Relevant Computational Drug Repurposing Modeling Studies
We reviewed the literature for other in silico studies that also identified some of these hit compounds as potential Mpro inhibitors and SARS-CoV-2 antiviral agents. Many drugs on the list in Supplementary Table S1 are predicted by published computational studies to be potential inhibitors of SARS-CoV-2 target proteins, largely Mpro but also RNA-dependent RNA polymerase (RdRp), spike, helicase, 2′-O-methyltransferase, nsp16/nsp10 complex, nsp1, PLpro, nsp3, and nsp12, and human angiotensin converting enzyme 2 (ACE2). Satisfyingly, those with the best predicted binding affinity from our study have also been of greatest interest clinically, with a larger number of in vitro assay results and clinical trials for drugs with the highest binding affinities (see below).
Modeling Studies Related to the Top 10 Predicted Drugs for Repurposing
Simeprevir was reported to be an inhibitor of the Mpro by ) They used a pharmacophore search followed by grid-based ligand docking (GLIDE, Schrodinger) and binding energy estimates from the MMGBSA method of −81.7 kcal/mol. However, they did not use MD to simulate the interaction of simeprevir in the Mpro binding site. Similarly, sofosbuvir was reported to be a strong inhibitor of the protease by .
The potential protease inhibition properties of lopinavir and ritonavir were reported by Bolcato et al., who used supervised MD to calculate the trajectories of the ligands in the protease binding site (). Muralidharan et al. also used AutoDock (another docking program similar to Vina produced by the Scripps group) followed by MD simulations using the Generalized Amber Force Field (GAFF) in Amber16 to screen for repurposed drugs (). They reported AutoDock binding energies for lopinavir, oseltamivir, and ritonavir of −4.1 kcal/mol, −4.65 kcal/mol, and −5.11 kcal/mol, respectively, but did not provide the binding energies from the MD calculations. The best-known antiviral drug, which has been the subject of several clinical trials for COVID-19, is remdesivir (). The potential inhibition of Mpro by this drug has been reported in several computational screening studies. For example, Al-Khafaji and colleagues reported a combined computational docking and MD study of a range of antiviral drugs to the viral protease (). They calculated a binding energy for remdesivir of −65.19 kcal/mol from a gromacs simulation and a MMGBSA binding energy calculation. Beck et al. reported a Kd for binding of remdesivir to 3CLPro of 113 nM using a deep learning model.
Novel potential Mpro inhibitors that emerged from our study included the ergot alkaloids ergotamine, mergocriptine, the thrombopoietin receptor agonist eltrombopag (ranked 13 with ΔGMMPBSA = –28.2 kcal/mol, see Supplementary Table S1), bemcentinib, PC786, and montelukast. These drugs were predicted to have better binding energies than the antiviral drugs discussed above and were not previously known to be antiviral.
Gurung et al. reported potential binding of ergotamine to the SARS-CoV-2 main protease in a preprint (). They employed AutoDock Vina but without subsequent MD simulation of the complex. They reported the binding energy as −9.4 kcal/mol for dihydroergotamine and −9.3 kcal/mol for ergotamine. Mevada et al. also reported in silico estimates of the binding of ergotamine to the protease using AutoDock Vina for the virtual screening (). They found the drug bound with an energy of −10.2 kcal/mol, calculated using Vina (no subsequent MD simulation). Gul et al. used a similar docking approach, this time with MD simulation, and identified ergotamine and its derivatives, dihydroergotamine and bromocriptine, as having high binding affinity to SARS-CoV-2 Mpro. Ergotamine is an alpha-1 selective adrenergic agonist and vasoconstrictor and exhibited a favorable docking binding energy against SARS-CoV-2 Mpro of −8.6 kcal/mol. Dihydroergotamine, the 9,10-alpha-dihydro derivative of ergotamine, showed a similar high affinity of −8.6 kcal/mol, and bromocriptine had a high affinity of −9.2 kcal/mol. Ergotamine has also been predicted to bind tightly to the SARS-CoV-2 spike (S) protein ().
Montelukast has been shown to inhibit at least one other protease, eosinophil protease (). Mansoor and colleagues deduced that it may bind to Mpro on the basis of a simple molecular docking study (). Wu et al. also reported putative binding of montelukast to Mpro in a computational study using the same Internal Coordinate Mechanics modeling methods (). No accurate binding affinities were reported in either study.
There is very little published work on the PC786 SARS-CoV-2 efficacy or predicted binding affinity to Mpro. Panda and coworkers reported a binding energy ΔGbind for PC786 of −179.79, tighter binding than calculated for lopinavir (−131.49 kJ/mol), using a combined docking and MD approach (). Like our study, they employed Autodock Vina to dock a molecular library into the active site of Mpro, followed by MD simulation using GROMACS.
Relevant Modeling Studies of the Drugs From the List of 84 Drugs
Several in silico screening studies have identified eltrombopag as a potential SARS-CoV-2 drug. Feng et al.‘s study suggested that eltrombopag bound not only to the Mpro active site but also to the viral spike protein and to human ACE2 (). This potential synergistic polypharmacy could be particularly beneficial for treating COVID-19. Eltrombopag has also been proposed as a useful drug against SARS-CoV-2 spike protein on the basis of its predicted strong binding to a pocket in the fusion cores of S2 domain (). Eltrombopag was identified as having a high binding affinity to human ACE2, the primary binding site for the SARS-CoV-2 spike protein. This virtual screening study also used Autodock Vina, but no subsequent MD simulation was used for the top hit compounds from the screen. Surface Plasmon Resonance (SPR) was used to assess the binding of the drug to Mpro.
Other drugs with binding energies stronger than –25 kcal/mol include galicaftor (in clinical trial for cystic fibrosis), rolitetracycline (a broad spectrum antibiotic), disogluside (a natural product from Dioscorea nipponica Makino that reduces liver chronic inflammation and fibrosis), zafirlukast (a leukotriene receptor antagonist for asthma), diosmin (a natural flavone for treating venous disease), AZD-5991 (in clinical trial for relapsed or refractory hematologic malignancies), and ruzasvir (in clinical trials for treatment of hepatitis C). Li et al. also reported predicted Mpro binding for galicaftor (). These drugs and natural products merit assessment in SAR-CoV-2 assays and Mpro inhibition experiments.
As we stated in the introduction, few studies have used computational docking followed by MD simulation of the best repurposing candidates, and fewer still have reported experimental validation of the computational predictions. Here, we report a comprehensive review of experimentally determined protein target, in vitro, or in vivo activities of our 84 top binding drug candidates.
Experimental Validation of Biological Activity of Computational Repurposing Candidates
Clearly, blind computational predictions of likely Mpro activity and, potentially, SARS-CoV-2 activities are of limited use if the predictions are not validated experimentally. Supplementary Table S1 shows that 70% of the top 10 hit compounds have confirmed experimental activity against SARS-CoV-2, and 25% of the 84 entries in this table also have confirmed experimental activity either against SARS-CoV-2 or Mpro. Of the remaining 75% of putative repurposed drugs, most have not been studied experimentally so they may have relevant antiviral activity, at least in vitro. The relatively high experimental validation rate of compounds predicted to be strong binders to Mpro suggests that our computational paradigm is useful for selecting drugs for repurposing against SARS-CoV-2. Clearly, computational studies are investigating molecular interactions at the target, which will be most useful for identifying candidates with strong target activity, Mpro in this instance. The high experimental validation rate also strongly suggests that the drugs not yet experimentally tested should at least be screened in an in vitro antiviral assay.
Validation of SARS-CoV-2 Activity of Top 10 Predicted Drugs for Repurposing
The website DrugVirus.info provides a concise picture of the broad-spectrum antiviral activity of a range of drugs. A summary for four of the top 10 antiviral hits (Table 1) from our in silico screens is provided in Figure 7. Here, we discuss the experimental SARS-CoV-2 or molecular target activity of seven of the top ten repurposing drug candidates identified by our computational studies—bemcentinib, montelukast, simeprevir, sofosbuvir, lopinavir, ritonavir, and remdesivir.
FIGURE 7
Bemcentinib selectively inhibits AXL kinase activity, which blocks viral entry and enhances the antiviral type I interferon response. Its in vitro activity against SARS-CoV-2 has been assessed by several groups. In a Vero cell assay, Liu et al. reported 10–40% protection at 50 µM (). However, in an alternative assay using human Huh7.5 cells (), bemcentinib exhibited an IC50 of 100 nM and CC50 of 4.7 µM. These authors also developed an assay in Vero cells and reported an IC50 of 470 nM and CC50 of 1.6 µM, considerably higher activity than that reported by Liu et al. As a result, it is an investigational treatment for COVID-19 (www.clinicaltrialsregister.eu), with a phase 2 trial underway (). Dittmar and co-workers also reported an ED50 for bemcentinib of 0.1 µM (Huh7.5 cells), 0.47 µM (Vero cells), and 2.1 µM (Calu3 cells) (). Six of the top 10 drugs (Table 1) are antiviral agents. Using a Vero E6 cellular infection model, they also reported that simeprevir was the only drug among their prioritized candidates that suppressed SARS-CoV-2 replication at below 10 μM. Dose–response studies showed that simeprevir had an EC50 of 4 μM and a CC50 of 20 μM, similar to remdesivir in their experiments. Simeprevir had an experimental in vitro EC50 activity of 4.08 μM. Ma et al. developed a fluorescence resonance energy transfer (FRET)-based enzymatic assay for the SARS-CoV-2 Mpro and used it to screen a library of Mpro inhibitors (). In their assay, simeprevir exhibited an IC50 of 14 ± 3 µM.
Sofosbuvir was reported to be a strong inhibitor of the main protease by . It has in vitro antiviral EC50 values of 6.2 and 9.5 μM (). Lopinavir exhibits an antiviral in vitro EC50 of 5.7 µM in one study (), and an EC50 of 26.6 μM in another study (). It is also the subject of multiple single-agent and combination human trials (e.g., ; ). Ritonavir has an experimental in vitro EC50 of 8.6 µM (), and it too is being assessed in multiple single-agent and combination human trials (e.g., ; ). Costanzo and colleagues likewise reported high protease binding for these two antiviral drugs (). They also reported updates on experimental drugs successfully employed in the treatment of the disease caused by SARS-CoV-2 coronavirus. Patient recovery has been reported after treatment with lopinavir/ritonavir (used to treat HIV infection) in combination with the anti-flu drug oseltamivir.
Remdesivir has also been assessed in multiple human trials (e.g., ; ), and it has reported antiviral in vitro EC50 values of 23.2 μM () and 0.77 μM and a CC50 > 100 μM in another study (). It had a SARS-CoV-2 EC50 in Vero cells of 6.6 μM and CC50 > 100 µM (). Beck et al. reported a Kd for binding of remdesivir to Mpro of 113 nM using a deep learning model. Liu et al. reported an in vitro assay that exploited the pronounced cytopathic effects of SAR-CoV-2 on Vero cells and the ability of a range of antiviral drugs to protect cells against the virus (). In their assay, remdesivir exhibited an IC50 of 2.5 µM and a CC50 of 175 µM, while sofosbuvir, lopinavir, and ritonavir were inactive.
Montelukast has been shown to produce a significant reduction in SARS-CoV-2 infection in the treated elderly asthmatic patients. Kumar et al. also reported in vitro SARS-CoV-2 inhibition, with an IC50 of 18.8 µM and CC50 > 20 µM (; ).
Other Novel Putative SARS-CoV-2 Drugs From the List of 84 Drugs
We also highlight some novel and interesting drugs for repurposing in the list of 74 (Supplementary Table S1) with weaker predicted binding affinity than the top 10 listed in Table 1. As stated above, 25% of the drugs in Supplementary Table S1 have reported experimental data that support the validity of predictions from our computational experiments.
Eltrombopag, a thrombopoietin receptor agonist used to treat thrombocytopenia, has a reported IC50 8.3 µM for SARS-CoV-2 infection in Vero and Calu-3 cells (). Recently, Vogel et al. reported direct inhibition of cytomegalovirus (CMV) by therapeutic doses of eltrombopag used to treat thrombocytopenia (). They showed that eltrombopag inhibits the late stages of the HCMV replication cycle and reduces virus titers by 1.8 × 104-fold at 10 µM and by 15-fold at 500 nM.
Saquinavir, an HIV protease inhibitor used in combination with other antiretroviral agents for the treatment of HIV-1, displays an in vitro EC50 of 8.8 µM (); zafirlukast, a leukotriene receptor antagonist used for prophylaxis and chronic treatment of asthma, exhibits an in vitro SARS-CoV-2 IC50 value of 3.6 µM (). Zhu and coworkers also measured the SARS-CoV-2 and Mpro inhibition of zafirlukast (Zhu et al., 2020). The IC50 for Mpro was 24 µM and the EC50 for the virus >20 µM.
Eravacycline, a tetracycline antibiotic used to treat complicated intra-abdominal infections, has a reported in vitro activity against recombinant SARS-CoV-2, SARS-CoV, and MERS-CoV main proteases, with IC50 values of 1.7, 10.0, and 16.4 µM, respectively. It also inhibits SARS-CoV-2 infection in VeroE6 cells with an IC50 = 30.6 µM (). Umifenovir (Arbidol), exerts antiviral effects through multiple pathways that see its use against a variety of enveloped and non-enveloped RNA and DNA viruses. It inhibits coronavirus OC43 with an IC50 of 4.4 µM and SARS-CoV-2 with an IC50 of 10 µM (). It has also been reported to inhibit SARS-CoV-2 infection at 10–30 μM in vitro (). Multiple clinical trials show a larger negative rate of PCR on day 14 in adult COVID-19 patients (), and it shortens the viral shedding interval (). Atazanavir, an antiretroviral drug of the protease inhibitor (PI) class, displays SARS-CoV-2 inhibition (in Vero cells) with an EC50 of 2 μM and a CC50 of 312 µM. It is also active against SARS-CoV-2 infection in a human epithelial pulmonary cell line (A549) with an EC50 of 0.22 µM (). Zhu and corkers measured the SARs-CoV-2 and Mpro inhibition of zafirlukast (Zhu et al., 2020). The IC50 for Mpro was 24 µM and the EC50 for the virus is >20 µM.
The protease binding of rolitetracycline has been reported by ) and ) The potential of the natural product diosmin as an antiviral agent targeting Mpro has also been reported in several recent computational studies (; ; ; ). Chakraborti et al. reported the potential of ruzasvir as a drug against SARS-CoV-2, although no data were provided ().
The prognostic value of our computational approach has been demonstrated by the fact that it identified a diverse range of drugs that have been reported in other computational studies or that exhibit useful SARS-CoV-2 antiviral effects in vitro. The antiviral drugs simeprevir, sofosbuvir, lopinavir, ritonavir, and remdesivir exhibit strong antiviral properties, and several are in clinical trial or used against SARS-CoV-2. These drugs have also been reported as binding to Mpro by numerous virtual screening studies, and by in vitro assays. The more interesting and least studied hit drugs among our candidate list, bemcentinib, PC786, montelukast, ergotamine, and mergocriptine, were predicted to have binding affinities equal to or greater than the antiviral drugs, and have also been shown to have in vitro antiviral activity against SARS-CoV-2. A few computational studies mostly using less rigorous methods than those employed here have also suggested that these drugs may bind to Mpro.
This high validation success rate strongly suggests that this type of virtual screening approach is capable of identifying compounds with potentially useful activity against SARS-CoV-2 and, by analogy, other coronaviruses. In particular, the 28 drugs for which no SARS-CoV-2 activity has been yet reported may be of particular interest for in vitro screening. The results of the current drug repurposing study provide information that could be useful to identify additional candidate drugs for testing for use in the current pandemic, as well as a rational computational paradigm for identifying therapeutic agents for future viral pandemics.
Materials and Methods
Protein Structure Preparation and Grid Preparation
The crystal structure of the COVID-19 Mpro was downloaded from the RCSB PDB (http://www.rcsb.org; refcode 6Y2F) (Zhang L. et al., 2020).
Protein preparation and removal of non-essential and non-bridging water molecules for docking studies were performed using the UCSF Chimera package (https://www.cgl.ucsf.edu/chimera/) () 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 (). The surface area of the 3CLPro binding pocket is 335 Å2, and the volume is 364.101 Å3 ().
As recommended by Llanos et al., the ability of Vina to redock known ligands from x-ray structures was assessed to determine the reliability of the algorithm for this target. Table 2 shows the RMSD values for redocking the ligands for 10 experimental structures of Mpro with bound ligands. The relatively low RMSD values show that Vina can recapitulate the experimental binding poses well.
TABLE 2
| Crystal ligand | Docking score (kcal/mol) | RMSD (Å) |
|---|---|---|
| 5R7Y | −5.0 | 0.44 |
| 5R81 | −5.0 | 0.51 |
| 5RE4 | −4.2 | 0.62 |
| 5REJ | −5.3 | 0.52 |
| 5RG0 | −4.5 | 0.55 |
| 6LU7 | −6.6 | 0.44 |
| 6W63 | −7.2 | 0.39 |
| 5R7Z | −5.3 | 0.57 |
| 5REL | −5.6 | 0.51 |
| 5R83 | −5.4 | 0.57 |
RMSD errors for redocking small molecule ligands into the binding site of Mpro.
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 8,773 and 13,308 drugs were retrieved from DrugBank and CHEMBL database, respectively. The drugs were downloaded in sdf format and converted to pdbqt format using Raccoon ().
Docking Methodology Small-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 that includes an approximate correction for ligand conformational entropy. AutoDock Vina is also flexible, easily scripted, and 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 Mpro (3CLPro, refcode 6Y2F) with the grid coordinates (grid center) 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 1 py (Supplementary Information) and subjected to molecular dynamic simulation. The docked structures were analyzed using UCSF Chimera () and LigPlot + software () to illustrate hydrogen-bond and hydrophobic interactions. A total of fifty top compounds were selected from each of the DrugBank and CHEMBL compounds. Sixteen compounds were common to both database top hits. Molecular dynamics studies were conducted on the unique set of eighty-four compounds from both sets.
Molecular Dynamics Simulation The top screened compound complexes with protease were minimized with CHARMm force field. The topology files of the ligands were prepared from Swissparam (http://www.swissparam.ch/) (Zoete et al., 2011) and minimized in Gromacs 2020 (http://www.gromacs.org/) (). Docked complexes of ligands and COVID-19 Mpro 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 cutoff 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 2,500 cycles of steepest descent minimization followed by 5,000 cycles of conjugate gradient minimization. After equilibration at 300 K using Langevin thermostat NVT ensemble for 50 ps, the system was then equilibrated at 1 atm pressure using Berendsen thermostat NPT ensemble for 50 ps. After the system was fully equilibrated at the desired temperature and pressure (NVT/NPT ensembles), we used Parrinelo–Rahman pressure coupling to run MD for data collection. Duplicate production runs starting with different random seeds were also run to allow estimates of binding energy uncertainties to be determined. Finally, a production run of 20 ns of MD simulation was performed.
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 20 ns is sufficiently long for convergence. We checked the RMSD of MPro and drug during this time and it was within the range of 1.5 Å. The RMSF graph revealed minimal fluctuations and relatively stable conformations of SARS-CoV-2 Mpro bound to screened drugs.
The protein–ligand binding affinities were evaluated in two ways. One calculates the energies of solvated SARS-CoV-2 protease and small-molecule ligands and the other calculates that of the bound complex and derive the binding energy by subtraction.
We also calculated binding free energies using the molecular mechanics Poisson Boltzmann surface area (MM/PBSA) tool in GROMACS that uses the nonbonded interaction energies of the complex. The method is also a widely used method for binding free energy calculations (). However, accurate calculation of absolute binding energies requires very extensive sampling, so the methods we employed provide accurate relative binding energies of ligands that are useful for ranking them, as we have done in this work.
We used GMXPBSA2.1 program to perform MM/PBSA calculations on selected docked complexes derived from GROMACS trajectories (). It is a suite of 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, which provides the freedom to calculate free binding energy of complexes with any force field, 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, and 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 20-ns trajectory files.
Specifically, for a non-covalent binding interaction in the aqueous phase, the binding free energy, ΔG (bind,aq), is:where ΔG (bind, vac) is the binding free energy in vacuum, and ΔG (bind, solv) is the solvation free energy change upon binding:where ΔG (R:L,solv), ΔG (R,solv), and ΔG (L,solv) are solvation free energies of complex, receptor, and ligand, respectively.
Method Note Added in Proof
Guterres and Im recently showed how substantial improvements in protein–ligand docking results could be achieved using high-throughput MD simulations (). As with our study, they also employed AutoDock Vina for docking, followed by MD simulation using CHARMM. The MD parameters they advocated were very similar to those used in our study. Proteins were solvated in a box of TIP3P water molecules extending 10 Å beyond the proteins and the particle-mesh Ewald method was used for electrostatic interactions. Nonbonded interactions over 10 and 12 Å were truncated. Their systems were minimized for 5,000 steps using the steepest descent method followed by 1 ns of equilibration with an NVT setting. For each protein–ligand complex, they ran 3 × 100-ns production runs from the same initial structure using different initial velocity random seeds and an integration step size of 2 fs. Over 56 protein targets (of seven different protein classes) and 560 ligands, this shows 22% improvement in the area under the receiver operating characteristic 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.
Statements
Data availability statement
The data have been deposited in the OPAL repository at La Trobe University and are available at DOI 10.26181/19235004.
Author contributions
NP conceived the project, analyzed data, and contributed to the manuscript; SP and PS performed the computations, analyzed data, and contributed to the manuscript; DW analyzed data and contributed to the manuscript.
Funding
We would like to thank Oracle for providing their Cloud computing resources for the modeling studies described herein.
Acknowledgments
We would like to thank Harinda Rajapaksha for assistance in optimizing gromacs for this project. We would also like to thank Oracle for providing their Cloud computing resources for the modeling studies described herein. In particular, we wish to thank Peter Winn, Dennis Ward, and Alison Derbenwick Miller from Oracle in facilitating these studies. The opinions expressed herein are solely those of the individual authors and should not be inferred to reflect the views of their affiliated institutions, funding bodies, or Oracle corporation.
Conflict of interest
SP, PS, and NP were employed by Vaxine Pty Ltd.
The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.781039/full#supplementary-material
References
1
AbhithajJ.DileepF.SharanyaC. S.ArunK. G.SadasivanC.JayadeviV. (2022). Repurposing Simeprevir, Calpain Inhibitor IV and a Cathepsin F Inhibitor Against SARS-CoV-2 and Insights into Their Interactions With Mpro. J. Biomol. Struct. Dyn.40 (1), 325–336. 10.1080/07391102.2020.1813200
2
AbrahamM. J.MurtolaT.SchulzR.PállS.SmithJ. C.HessB.et al (2015). GROMACS: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX1-2, 19–25. 10.1016/j.softx.2015.06.001
3
Al-KhafajiK.Al-DuhaidahawiD.Taskin TokT. (2021). Using Integrated Computational Approaches to Identify Safe and Rapid Treatment for SARS-CoV-2. J. Biomol. Struct. Dyn.39, 3387–3395. 10.1080/07391102.2020.1764392
4
AlmerieM. Q.KerriganD. D. (2020). The Association between Obesity and Poor Outcome after COVID-19 Indicates a Potential Therapeutic Role for Montelukast. Med. Hypotheses143, 109883. 10.1016/j.mehy.2020.109883
5
ArunK. G.SharanyaC. S.AbhithajJ.FrancisD.SadasivanC. (2020). Drug Repurposing against SARS-CoV-2 Using E-Pharmacophore Based Virtual Screening, Molecular Docking and Molecular Dynamics with Main Protease as the Target. J. Biomol. Struct. Dyn., 1–12. 10.1080/07391102.2020.1779819
6
BerkleyS. (2020). COVID-19 Needs a Big Science Approach. Science367, 1407. 10.1126/science.abb8654
7
BolcatoG.BissaroM.PavanM.SturleseM.MoroS. (2020). Targeting the Coronavirus SARS-CoV-2: Computational Insights into the Mechanism of Action of the Protease Inhibitors Lopinavir, Ritonavir and Nelfinavir. Sci. Rep.10, 20927. 10.1038/s41598-020-77700-z
8
BozekA.WintersteinJ.GaluszkaB.JarzabJ. (2020). Different Development Forms of Local Allergic Rhinitis towards Birch. Biomed. Res. Int.2020, 1–9. 10.1155/2020/3408561
9
CaoB.WangY.WenD.LiuW.WangJ.FanG.et al (2020). A Trial of Lopinavir-Ritonavir in Adults Hospitalized with Severe Covid-19. N. Engl. J. Med.382, 1787–1799. 10.1056/nejmoa2001282
10
ChakrabortiS.BheemireddyS.SrinivasanN. (2020). Repurposing Drugs against the Main Protease of SARS-CoV-2: Mechanism-Based Insights Supported by Available Laboratory and Clinical Data. Mol. Omics16, 474–491. 10.1039/d0mo00057d
11
ChoyK.-T.WongA. Y.-L.KaewpreedeeP.SiaS. F.ChenD.HuiK. P. Y.et al (2020). Remdesivir, Lopinavir, Emetine, and Homoharringtonine Inhibit SARS-CoV-2 Replication In Vitro. Antiviral Res.178, 104786. 10.1016/j.antiviral.2020.104786
12
CiottiM.AngelettiS.MinieriM.GiovannettiM.BenvenutoD.PascarellaS.et al (2019). COVID-19 Outbreak: An Overview. Chemotherapy64, 215–223. 10.1159/000507423
13
CostanzoM.De GiglioM. A. R.RovielloG. N. (2020). SARS-CoV-2: Recent Reports on Antiviral Therapies Based on Lopinavir/Ritonavir, Darunavir/Umifenovir, Hydroxychloroquine, Remdesivir, Favipiravir and Other Drugs for the Treatment of the New Coronavirus. Curr. Med. Chem.27, 4536–4541. 10.2174/0929867327666200416131117
14
DaiW.ZhangB.JiangX.-M.SuH.LiJ.ZhaoY.et al (2020). Structure-based Design of Antiviral Drug Candidates Targeting the SARS-CoV-2 Main Protease. Science368, 1331–1335. 10.1126/science.abb4489
15
DittmarM.LeeJ. S.WhigK.SegristE.LiM.KamaliaB.et al (2021). Drug Repurposing Screens Reveal Cell-Type-Specific Entry Pathways and FDA-Approved Drugs Active Against SARS-Cov-2. Cell Rep.35 (1), 108959. 10.1016/j.celrep.2021.108959
16
DurdagiS. (2020). Virtual Drug Repurposing Study against SARS-CoV-2 TMPRSS2 Target. Turk J. Biol.44, 185–191. 10.3906/biy-2005-112
17
FengS.LuanX.WangY.WangH.ZhangZ.WangY.et al (2020). Eltrombopag Is a Potential Target for Drug Intervention in SARS-CoV-2 Spike Protein. Infect. Genet. Evol.85, 104419. 10.1016/j.meegid.2020.104419
18
Fintelman-RodriguesN.SacramentoC. Q.Ribeiro LimaC.Souza Da SilvaF.FerreiraA. C.MattosM.et al (2020). Atazanavir, Alone or in Combination with Ritonavir, Inhibits SARS-CoV-2 Replication and Proinflammatory Cytokine Production. Antimicrob. Agents Chemother.64, e00825–00820. 10.1128/AAC.00825-20
19
ForliS.HueyR.PiqueM. E.SannerM. F.GoodsellD. S.OlsonA. J. (2016). Computational Protein-Ligand Docking and Virtual Drug Screening with the AutoDock Suite. Nat. Protoc.11, 905–919. 10.1038/nprot.2016.051
20
GaultonA.HerseyA.NowotkaM.BentoA. P.ChambersJ.MendezD.et al (2017). The ChEMBL Database in 2017. Nucleic Acids Res.45 (D1), D945–D954. 10.1093/nar/gkw1074
21
GulS.OzcanO.AsarS.OkyarA.BarısI.KavakliI. H. (2020). In Silico identification of Widely Used and Well-Tolerated Drugs as Potential SARS-CoV-2 3C-like Protease and Viral RNA-dependent RNA Polymerase Inhibitors for Direct Use in Clinical Trials. J. Biomol. Struct. Dyn., 39(17):6772–6791. 10.1080/07391102.2020.1802346
22
GurungA. B.AliM. A.LeeJ.Abul FarahM.Al-AnaziK. M. (2020). In Silico screening of FDA Approved Drugs Reveals Ergotamine and Dihydroergotamine as Potential Coronavirus Main Protease Enzyme Inhibitors. Saudi J. Biol. Sci.27, 2674–2682. 10.1016/j.sjbs.2020.06.005
23
GuterresH.ImW. (2020). Improving Protein-Ligand Docking Results with High-Throughput Molecular Dynamics Simulations. J. Chem. Inf. Model.60, 2189–2198. 10.1021/acs.jcim.0c00057
24
HendausM. A. (2021). Remdesivir in the Treatment of Coronavirus Disease 2019 (COVID-19): a Simplified Summary. J. Biomol. Struct. Dyn.39, 3787–3792. 10.1080/07391102.2020.1767691
25
HuangD.YuH.WangT.YangH.YaoR.LiangZ. (2021). Efficacy and Safety of Umifenovir for Coronavirus Disease 2019 (COVID‐19): A Systematic Review and Meta‐analysis. J. Med. Virol.93, 481–490. 10.1002/jmv.26256
26
HuangH.GuanL.YangY.Le GrangeJ. M.TangG.XuY.et al (2020). Chloroquine, Arbidol (Umifenovir) or Lopinavir/ritonavir as the Antiviral Monotherapy for COVID-19 Patients: a Retrospective Cohort Study. Research Square. 10.21203/rs.3.rs-24667/v1
27
JiangS. (2020). Don't rush to Deploy COVID-19 Vaccines and Drugs without Sufficient Safety Guarantees. Nature579, 321. 10.1038/d41586-020-00751-9
28
KoM.JeonS.RyuW. S.KimS. (2021). Comparative Analysis of Antiviral Efficacy of FDA‐approved Drugs against SARS‐CoV‐2 in Human Lung Cells. J. Med. Virol.93, 1403–1408. 10.1002/jmv.26397
29
KumarS.SinghB.KumariP.KumarP. V.AgnihotriG.KhanS.et al (2021). Identification of Multipotent Drugs for COVID-19 Therapeutics with the Evaluation of Their SARS-CoV2 Inhibitory Activity. Comput. Struct. Biotechnol. J.19, 1998–2017. 10.1016/j.csbj.2021.04.014
30
LangloisA.FerlandC.TremblayG. M.LavioletteM. (2006). Montelukast Regulates Eosinophil Protease Activity through a Leukotriene-independent Mechanism. J. Allergy Clin. Immunol.118, 113–119. 10.1016/j.jaci.2006.03.010
31
LaskowskiR. A.SwindellsM. B. (2011). LigPlot+: Multiple Ligand-Protein Interaction Diagrams for Drug Discovery. J. Chem. Inf. Model.51, 2778–2786. 10.1021/ci200227u
32
LiY.ZhangY.HanY.ZhangT.DuR. (2020). Prioritization of Potential Drugs Targeting the SARS-CoV-2 Main Protease. ChemRxiv. 10.26434/chemrxiv.12629858.v1
33
LiuS.LienC. Z.SelvarajP.WangT. T. (2020). Evaluation of 19 Antiviral Drugs against SARS-CoV-2 Infection. bioRxiv. 10.1101/2020.04.29.067983
34
LlanosM. A.GantnerM. E.RodriguezS.AlbercaL. N.BelleraC. L.TaleviA.et al (2021). Strengths and Weaknesses of Docking Simulations in the SARS-CoV-2 Era: the Main Protease (Mpro) Case Study. J. Chem. Inf. Model.61, 3758–3770. 10.1021/acs.jcim.1c00404
35
LoH. S.HuiK. P.LaiH.-M.KhanK. S.KaurS.HuangJ.et al (2021). Simeprevir Potently Suppresses SARS-CoV-2 Replication and Synergizes with Remdesivir. ACS Central Science7 (5), 792–802. 10.1021/acscentsci.0c01186
36
LuS. (2020). Timely Development of Vaccines against SARS-CoV-2. Emerging Microb. Infect.9, 542–544. 10.1080/22221751.2020.1737580
37
MaC.SaccoM. D.HurstB.TownsendJ. A.HuY.SzetoT.et al (2020). Boceprevir, GC-376, and Calpain Inhibitors II, XII Inhibit SARS-CoV-2 Viral Replication by Targeting the Viral Main Protease. Cell Res30, 678–692. 10.1038/s41422-020-0356-z
38
MansoorS.SaadatS.AminA.AliI.GhaffarM. T.AminU.et al (2020). A Case for Montelukast in COVID-19: The Use of Computational Docking to Estimate the Effects of Montelukast on Potential Viral Main Protease Catalytic Site. Research Square. [Preprint]. 10.21203/rs.3.rs-27079/v1
39
MendesA. (2020). Research towards Treating COVID-19. Br. J. Community Nurs.25, 204–205. 10.12968/bjcn.2020.25.4.204
40
MengistH. M.FanX.JinT. (2020). Designing of Improved Drugs for COVID-19: Crystal Structure of SARS-CoV-2 Main Protease Mpro. Sig Transduct Target. Ther.5, 67. 10.1038/s41392-020-0178-y
41
MevadaV.DudhagaraP.GandhiH.VaghamshiN.BeladiyaU.PatelR. (2020). Drug Repurposing of Approved Drugs Elbasvir, Ledipasvir, Paritaprevir, Velpatasvir, Antrafenine and Ergotamine for Combating COVID19. ChemRxiv. 10.26434/chemrxiv.12115251.v2
42
MuralidharanN.SakthivelR.VelmuruganD.GromihaM. M. (2021). Computational Studies of Drug Repurposing and Synergism of Lopinavir, Oseltamivir and Ritonavir Binding with SARS-CoV-2 Protease against COVID-19. J. Biomol. Struct. Dyn.39, 2673–2678. 10.1080/07391102.2020.1752802
43
NgoS. T.Quynh Anh PhamN.Thi LeL.PhamD.-H.VuV. V. (2020). Computational Determination of Potential Inhibitors of SARS-CoV-2 Main Protease. J. Chem. Inf. Model.60, 5771–5780. 10.1021/acs.jcim.0c00491
44
OlenderS. A.PerezK. K.GoA. S.BalaniB.Price-HaywoodE. G.ShahN. S.et al (2021). Remdesivir for Severe COVID-19 versus a Cohort Receiving Standard of Care. Clin. Infect. Dis.73 (11), e4166–e4174. ciaa1041. 10.1093/cid/ciaa1041
45
OlsenM.CookS. E.HuangV.PedersenN.MurphyB. G. (2020). 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. Agents55, 105964. 10.1016/j.ijantimicag.2020.105964
46
PaissoniC.SpiliotopoulosD.MuscoG.SpitaleriA. (2015). GMXPBSA 2.1: A GROMACS Tool to Perform MM/PBSA and Computational Alanine Scanning. Computer Phys. Commun.186, 105–107. 10.1016/j.cpc.2014.09.010
47
PandaP. K.ArulM. N.PatelP.VermaS. K.LuoW.RubahnH. G.et al (2020). Structure-based Drug Designing and Immunoinformatics Approach for SARS-CoV-2. Sci. Adv.6, eabb8097. 10.1126/sciadv.abb8097
48
PetersonL. (2020a). COVID-19 and Flavonoids: In Silico Molecular Dynamics Docking to the Active Catalytic Site of SARS-CoV and SARS-CoV-2 Main Protease. SSRN Electronic Journal. 10.2139/ssrn.3599426
49
PetersonL. (2020b). In Silico Molecular Dynamics Docking of Drugs to the Inhibitory Active Site of SARS-CoV-2 Protease and Their Predicted Toxicology and ADME. SSRN Electronic Journal. 10.2139/ssrn.3580951
50
PettersenE. F.GoddardT. D.HuangC. C.CouchG. S.GreenblattD. M.MengE. C.et al (2004). UCSF Chimera?A Visualization System for Exploratory Research and Analysis. J. Comput. Chem.25, 1605–1612. 10.1002/jcc.20084
51
PirzadaR. H.HaseebM.BatoolM.KimM.ChoiS. (2021). Remdesivir and Ledipasvir Among the FDA-Approved Antiviral Drugs Have Potential to Inhibit SARS-CoV-2 Replication. Cells10, 1052. 10.3390/cells10051052
52
QiaoZ.ZhangH.JiH.-F.ChenQ. (2020). Computational View toward the Inhibition of SARS-CoV-2 Spike Glycoprotein and the 3CL Protease. Computat.8, 53. 10.3390/computation8020053
53
RathnayakeA. D.ZhengJ.KimY.PereraK. D.MackinS.MeyerholzD. K.et al (2020). 3C-like Protease Inhibitors Block Coronavirus Replication In Vitro and Improve Survival in MERS-CoV-Infected Mice. Sci. Transl Med.12, eabc5332. 10.1126/scitranslmed.abc5332
54
ReigN.ShinD.-H. (2020). 560. Repurposing Eravacycline for the Treatment of SARS-CoV-2 Infections. Open Forum Infect. Dis.7, S345. 10.1093/ofid/ofaa439.754
55
RosaS. G. V.SantosW. C. (2020). Clinical Trials on Drug Repositioning for COVID-19 Treatment. Rev. Panam Salud Publica44, e40. 10.26633/RPSP.2020.40
56
Rosales-MendozaS.Márquez-EscobarV. A.González-OrtegaO.Nieto-GómezR.Arévalo-VillalobosJ. I. (2020). What Does Plant-Based Vaccine Technology Offer to the Fight against COVID-19?Vaccines8, 183. 10.3390/vaccines8020183
57
SacramentoC. Q.Fintelman-RodriguesN.TemerozoJ. R.de Paula Dias Da SilvaA.Da Silva Gomes DiasS.dos Santos da SilvaC.et al (2021). In vitro Antiviral Activity of the Anti-HCV Drugs Daclatasvir and Sofosbuvir Against SARS-CoV-2, the Aetiological Agent of COVID-19. J Antimicrob Chemother.76 (7), 1874–1885. 10.1093/jac/dkab072
58
SandersJ. M.MonogueM. L.JodlowskiT. Z.CutrellJ. B. (2020). Pharmacologic Treatments for Coronavirus Disease 2019 (COVID-19): A Review. JAMA323, 1824–1836. 10.1001/jama.2020.6019
59
SaulS.EinavS. (2020). Old Drugs for a New Virus: Repurposed Approaches for Combating COVID-19. ACS Infect. Dis.6, 2304–2318. 10.1021/acsinfecdis.0c00343
60
SchlagenhaufP.GrobuschM. P.MaierJ. D.GautretP. (2020). Repurposing Antimalarials and Other Drugs for COVID-19. Trav. Med Infect Dis34, 101658. 10.1016/j.tmaid.2020.101658
61
Sheik AmamuddyO.VerkhivkerG. M.Tastan BishopÖ. (2020). Impact of Early Pandemic Stage Mutations on Molecular Dynamics of SARS-CoV-2 Mpro. J. Chem. Inf. Model.60, 5080–5102. 10.1021/acs.jcim.0c00634
62
SohrabiC.AlsafiZ.O'neillN.KhanM.KerwanA.Al-JabirA.et al (2020). World Health Organization Declares Global Emergency: A Review of the 2019 Novel Coronavirus (COVID-19). Int. J. Surg.76, 71–76. 10.1016/j.ijsu.2020.02.034
63
SpiliotopoulosD.SpitaleriA.MuscoG. (2012). Exploring PHD Fingers and H3K4me0 Interactions with Molecular Dynamics Simulations and Binding Free Energy Calculations: AIRE-PHD1, a Comparative Study. PLoS One7, e46902. 10.1371/journal.pone.0046902
64
TarasovaA.WinklerD. A. (2009). Modelling Atypical Small-Molecule Mimics of an Important Stem Cell Cytokine, Thrombopoietin. ChemMedChem4, 2002–2011. 10.1002/cmdc.200900340
65
Thanh LeT.AndreadakisZ.KumarA.Gómez RománR.TollefsenS.SavilleM.et al (2020). The COVID-19 Vaccine Development Landscape. Nat. Rev. Drug Discov.19, 305–306. 10.1038/d41573-020-00073-5
66
TianW.ChenC.LeiX.ZhaoJ.LiangJ. (2018). CASTp 3.0: Computed Atlas of Surface Topography of Proteins. Nucleic Acids Res.46, W363–W367. 10.1093/nar/gky473
67
VafaeiS.RazmiM.MansooriM.Asadi-LariM.MadjdZ. (2020). Spotlight of Remdesivir in Comparison with Ribavirin, Favipiravir, Oseltamivir and Umifenovir in Coronavirus Disease 2019 (COVID-19) Pandemic. SSRN Electronic Journal. 10.2139/ssrn.3569866
68
Verdugo-PaivaF.IzcovichA.RagusaM.RadaG. (2020). Lopinavir-ritonavir for COVID-19: A Living Systematic Review. Medwave20, e7967. 10.5867/medwave.2020.06.7966
69
VogelJ.-U.SchmidtS.SchmidtD.RothweilerF.KochB.BaerP.et al (2019). The Thrombopoietin Receptor Agonist Eltrombopag Inhibits Human Cytomegalovirus Replication via Iron Chelation. Cells9, 31. 10.3390/cells9010031
70
WangM.CaoR.ZhangL.YangX.LiuJ.XuM.et al (2020a). Remdesivir and Chloroquine Effectively Inhibit the Recently Emerged Novel Coronavirus (2019-nCoV) In Vitro. Cell Res30, 269–271. 10.1038/s41422-020-0282-0
71
WangY.ZhangD.DuG.DuR.ZhaoJ.JinY.et al (2020b). Remdesivir in Adults with Severe COVID-19: a Randomised, Double-Blind, Placebo-Controlled, Multicentre Trial. Lancet395, 1569–1578. 10.1016/s0140-6736(20)31022-9
72
WhitworthJ. (2020). COVID-19: a Fast Evolving Pandemic. Trans. R. Soc. Trop. Med. Hyg.114, 241–248. 10.1093/trstmh/traa025
73
WilkinsonT.DixonR.PageC.CarrollM.GriffithsG.HoL.-P.et al (2020). ACCORD: A Multicentre, Seamless, Phase 2 Adaptive Randomisation Platform Study to Assess the Efficacy and Safety of Multiple Candidate Agents for the Treatment of COVID-19 in Hospitalised Patients: A Structured Summary of a Study Protocol for a Randomised Controlled Trial. Trials21, 691. 10.1186/s13063-020-04584-9
74
WishartD. S.FeunangY. D.GuoA. C.LoE. J.MarcuA.GrantJ. R.et al (2018). DrugBank 5.0: A Major Update to the DrugBank Database for 2018. Nucleic Acids Res.46 (D1), D1074–D1082. 10.1093/nar/gkx1037
75
WuC.LiuY.YangY.ZhangP.ZhongW.WangY.et al (2020). Analysis of Therapeutic Targets for SARS-CoV-2 and Discovery of Potential Drugs by Computational Methods. Acta Pharmaceut. Sin.10, 766–788. 10.1016/j.apsb.2020.02.008
76
XiaoX.WangC.ChangWangD.WangY.DongX.JiaoT.et al (2020). Identification of Potent and Safe Antiviral Therapeutic Candidates against SARS-CoV-2. Front. Immunol.11, 586572. 10.3389/fimmu.2020.586572
77
YamamotoN.MatsuyamaS.HoshinoT.YamamotoN. (2020). Nelfinavir Inhibits Replication of Severe Acute Respiratory Syndrome Coronavirus 2 in Vitro. bioRxiv. 10.1101/2020.04.06.026476
78
YavuzS. S.ÜnalS. (2020). Antiviral Treatment of Covid-19. Turk. J. Med. Sci.50, 611–619. 10.3906/sag-2004-145
79
ZengJ.WeissmannF.BertolinA. P.PosseV.CanalB.UlfertsR.et al (2021). Identifying SARS-CoV-2 Antiviral Compounds by Screening for Small Molecule Inhibitors of Nsp13 Helicase. Biochem. J.478, 2405–2423. 10.1042/bcj20210201
80
ZhangJ.-J.ShenX.YanY.-M.WangY.ChengY.-X. (2020b). Discovery of Anti-SARS-CoV-2 Agents from Commercially Available Flavor via Docking Screening. 10.31219/osf.io/vjch2
81
ZhangJ.ZengH.GuJ.LiH.ZhengL.ZouQ. (2020a). Progress and Prospects on Vaccine Development against SARS-CoV-2. Vaccines8, 153. 10.3390/vaccines8020153
82
ZhangL.LinD.SunX.CurthU.DrostenC.SauerheringL.et al (2020c). Crystal Structure of SARS-CoV-2 Main Protease Provides a Basis for Design of Improved α-ketoamide Inhibitors. Science368, 409–412. 10.1126/science.abb3405
83
ZhangT.HeY.XuW.MaA.YangY.XuK.-F. (2020d). Clinical Trials for the Treatment of Coronavirus Disease 2019 (COVID-19): A Rapid Response to Urgent Need. Sci. China Life Sci.63, 774–776. 10.1007/s11427-020-1660-2
84
ZhuW.XuM.ChenC. Z.GuoH.ShenM.HuX.et al (2020). Identification of SARS-CoV-2 3CL Protease Inhibitors by a Quantitative High-Throughput Screening. ACS Pharmacol. Transl. Sci.3 (5), 1008–1016. 10.1021/acsptsci.0c00108
85
ZoeteV.CuendetM. A.GrosdidierA.MichielinO. (2011). SwissParam: a Fast Force Field Generation Tool for Small Organic Molecules. J. Comput. Chem.32, 2359–2368. 10.1002/jcc.21816
Summary
Keywords
SARS-CoV-2, binding affinity, main protease, 3CL, computational chemistry, docking, molecular dynamics
Citation
Piplani S, Singh P, Petrovsky N and Winkler DA (2022) Computational Repurposing of Drugs and Natural Products Against SARS-CoV-2 Main Protease (Mpro) as Potential COVID-19 Therapies. Front. Mol. Biosci. 9:781039. doi: 10.3389/fmolb.2022.781039
Received
22 September 2021
Accepted
28 January 2022
Published
14 March 2022
Volume
9 - 2022
Edited by
James Leland Olds, George Mason University, United States
Reviewed by
Yu Kang, Zhejiang University, China
Debsindhu Bhowmik, Oak Ridge National Laboratory (DOE), United States
Lance Hellman, Nevada State College, United States
Richard J. Bingham, University of Huddersfield, United Kingdom
Updates
Copyright
© 2022 Piplani, Singh, Petrovsky and Winkler.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nikolai Petrovsky, nikolai.petrovsky@flinders.edu.au; David A. Winkler, d.winkler@latrobe.edu.au
† These authors share first authorship
‡ These authors share senior authorship
This article was submitted to Biological Modeling and Simulation, a section of the journal Frontiers in Molecular Biosciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.