Computational View toward the Inhibition of SARS-CoV-2 Spike Glycoprotein and the 3CL Protease
1
Department of Orthopaedics, The Warren Alpert Medical School of Brown University, Rhode Island Hospital, Providence, RI 02903, USA
2
Departments of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19301, USA
3
Department of Chemistry, Drexel University, Philadelphia, PA 19104, USA
*
Author to whom correspondence should be addressed.
Computation 2020, 8(2), 53; https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020053
Submission received: 16 April 2020
/
Revised: 25 May 2020
/
Accepted: 27 May 2020
/
Published: 31 May 2020
(This article belongs to the Special Issue Computation to Fight SARS-CoV-2 (CoVid-19))
Abstract
:Since the outbreak of the 2019 novel coronavirus disease (COVID-19), the medical research community is vigorously seeking a treatment to control the infection and save the lives of severely infected patients. The main potential candidates for the control of viruses are virally targeted agents. In this short letter, we report our calculations on the inhibitors for the SARS-CoV-2 3CL protease and the spike protein for the potential treatment of COVID-19. The results show that the most potent inhibitors of the SARS-CoV-2 3CL protease include saquinavir, tadalafil, rivaroxaban, sildenafil, dasatinib, etc. Ergotamine, amphotericin b, and vancomycin are most promising to block the interaction of the SARS-CoV-2 S-protein with human ACE-2.
1. Introduction
As of 24 May 2020, over 5 millions people in the world have been confirmed as having the 2019 novel coronavirus disease (COVID-19), an infection with Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) (initially called 2019-nCoV before 11 February 2020) which is part of the Coronaviridae family of positive-sense single-stranded RNA viruses that includes SARS-CoV and MERS-CoV (Middle East Respiratory Syndrome coronavirus), both of which also cause severe respiratory infections. The death count in China so far has been over 1700, but the number is expected to go higher with the increasing number of confirmed and non-confirmed cases. The medical research community is vigorously seeking a treatment to control the infection and save the lives of severely infected patients.
Just a few weeks after the COVID-19 outbreak, the complete genome of SARS-CoV-2 was determined and reported to GenBank (accession MN908947). Viruses were also isolated from patients to understand the genomic characteristics and mechanism of the viral infection. As revealed by the analysis, the SARS-CoV-2 shared 79% sequence identity to SARS-CoV. In one study, SARS-CoV-2 was found to be closely related to two bat-derived Severe Acute Respiratory Syndrome (SARS)-like coronaviruses, with 87.5% and 87.6% shared identity [1]. In another study, SARS-CoV-2 was 96% identical at the whole-genome level to a bat coronavirus [2].
Despite the high sequence identity between the SARS-CoV-2 and the SARS-CoV in the open reading frame regions, the envelop spike protein (S-protein) [3], which mediates the infection of SARS-CoV via the human host protein ACE-2, has only about 80% shared sequence identity between the SARS-CoV and SARS-CoV-2 [1]. Within the S-protein, the receptor docking domain has a higher divergence, with four out of five critical ACE-2 interacting amino acid residues replaced in the SARS-CoV-2. However, structural modeling indicated that the four residues in the SARS-CoV-2 retain a structural conformation similar to that of SARS-CoV, and the SARS-CoV-2 S-protein should be able to bind ACE-2 with reasonable affinity [4]. Indeed, studies by Zhou et al. using cells expressing human ACE-2 confirmed that the SARS-CoV-2 could infect cells via the same protein on ACE-2 as SARS-CoV did [2]. Thus, one option to treat the infection is to search for an inhibitor that can prevent the interaction of the SARS-CoV-2 S-protein with human ACE-2. The availability of the genome sequence of SARS-CoV-2 allows us to establish structural models for the S-protein [4].
The RNA of coronaviruses encodes polyproteins that can be processed by viral proteases to yield mature proteins. The same mechanism is shared by picornaviruses and retroviruses. Patients treated with protease inhibitors appeared to have much better clinical outcomes than without using the inhibitors (SARS death: 28.8% vs. 2.4%) [5]. Molecular dynamics simulations have revealed that, by molecular docking to the active site of the main protease 3CL of SARS-CoV, both lopinavir and ritonavir could induce conformation changes and potentially interfere with infection by SARS virus [6]. We expect the same will apply for SARS-CoV-2. The crystal structure of the SARS-CoV-2 protease (3CLpro) was just recently reported by Liu et al. [7]. Thus, another option to treat the SARS-CoV-2 infection is to search for inhibitors of the SARS-CoV-2 3CLpro.
With these models and crystal data, we performed in silico studies of potential inhibitors of the SARS-CoV-2 S-protein and 3CLpro.
2. Computational Methods
All calculations were operated on Dell PowerEdge C6220 servers. The chemical structures were prepared by AutoDockTools-1.5.6 [8], Chimera 1.14 [9], and Avogadro [10]. The docking studies were performed with Autodock 4.2.6, Autodock4, AutoDockTools4 [11], and Autodock Vina 1.1.2 [12].
2.1. Preparation of Receptor and Ligands
The 3CL protease’s three-dimensional crystal structure was retrieved from the Protein Data Bank (PDB ID: 6LU7), and it was applied as the receptor for molecular docking after a cleaning with Chimera. The ligands observed, i.e., FDA-approved drugs (2454 structures in total), were retrieved from the BindingDB (https://www.bindingdb.org), and the structures of the ligands were further optimized with Avogadro. The force field applied for geometry optimization was MMFF94.
The SARS-CoV-2/ACE-2 structure was retrieved using the function of the comparative modeling of the Chimera interface with the modeler (version 9.23) [13]. For the preparation of the SARS-CoV-2/ACE-2 structure, the target template sequence was retrieved from Zhang et al.’s work and the SARS-CoV/ACE-2 (PDB ID: 6ACD) served as a template, as it was also the top candidate from Basic Local Alignment Search Tool (BLAST) results. Because SARS-CoV and SARS-CoV-2 have an 88% similarity, the 3D structure can be predicted with a high accuracy. Next, the sequence alignments were performed using SARS-CoV as a template. Then, the model was built followed by refining the loops, side chain optimization, and model optimization. When the homology model was generated, it was further validated using the WHATCHECK/PROCHECK program [14] for basic parameters like torsion angle, rotational angle, bond length, etc. Finally, this model was used as receptor for docking purposes. The loop refinement and side chain optimization were performed using Chimera 1.14 by selecting the active region; all the parameters were the default of the version.
It is noteworthy that this calculated work was performed before the crystal structure of the COVID-19 S-protein was released (6LZG, 6VW1, etc.). After the crystal structures were released, their structures were compared with ours and the structures overlaid well (Figure 1), with 93.22% of its residues in the allowed region and a minor difference on the top right loop, which was not a site that interacts with the ACE-2, so a re-calculation was not conducted using the new crystal structures.
2.2. Molecular Docking with Autodock Vina
For the SARS-CoV-2 3CL inhibition calculation, the input files for Autodock Vina were prepared in the receptor’s original file (PDB format) and ligands files (SDF format) using AutoDockTools-1.5.6. After minimizing, the grid box was set at 22.00 Å × 22.00 Å ×22.00 Å along the x, y, and z axis, respectively. The docking site was defined at 1.00 Å when using the Autodock Vina. The grid box was set into the docking site at the H41, C145, and E166 regions according to the docking site of the coronavirus main proteinase (3CL) of Severe Acute Respiratory Syndrome (SARS). Then, the receptor file (PDBQT format, for docking purposes) was prepared by the addition of polar only hydrogen atoms, the removal of all water molecules, and the calculation of the Gasteiger charge. The instructed command prompts were used for the docking process. The docking output file includes the docking energy (in kcal/mol, which is an indication of the binding affinity/efficiency of one specific ligand to the receptor molecule) and the interaction of the ligands with the receptor (hydrogen bond, pi-pi stacking, etc.).
For the SARS-CoV-2 S-protein inhibition calculation, the PDB files of the SARS-CoV-2 S-protein were generated using the homology modeling method in Chimera; the template used for this was the SARS-CoV S-protein. After minimization, the input file was prepared using AutoDockTools-1.5.6. The grid box, which was a rectangular shaped area that covered all the possible docking sites of the SARS-CoV-2 S-protein with its receptor ACE-2, was chosen as 22.00 Å × 42.00 Å ×22.00 Å along the x, y, and z axis, respectively. The docking site was defined at 1.00 Å when using the Autodock Vina. Then, the receptor file (PDBQT format, for docking purposes) was prepared by the addition of polar only hydrogen atoms, the removal of all water molecules, and the calculation of the Gasteiger charge.
2.3. Analyzing the Docking Results with Chimera and BioLuminate
The docking results were ranked in the order from high to low in different modes according to the docking scores (docking energy, kcal/mol). The ligands with the most negative docking scores—i.e., the highest affinities—were selected for the visualization of the docked complexes using Chimera [9].
The docking energies of the SARS-CoV-2 S-protein and human ACE-2 were calculated using BioLuminate [15,16,17], and then compared to the docking energies of the SARS-CoV S-protein and human ACE-2. To verify whether those ligands can be used for blocking the interaction of the S-protein with human ACE-2, the docking energies of the SARS-CoV-2 S-protein/ligands and human ACE-2 were also calculated. The solvation model used was VSGB [18], and the force field chosen was OPLS_2005 [19] for all the docking energy predictions.
3. Results
3.1. Results of the SARS-CoV-2 3CL Protease
Table 1 shows the binding affinity of several ligands with SARS-CoV-2 3CL protease sorted according to the docking scores (binding affinities) calculated from the Autodock Vina; Figure 2 shows the docking of those with high docking scores—Tadalafil, Dasatinib, and Saquinavir—with the protease in the docking sites of the protease.
3.2. Results of SARS-CoV-2 S-Protein
We modeled ligands that may bind at a large docking area on the top of the S-protein that interacts with ACE-2 (red cycle in Figure 3a). Table 2 shows the binding affinities of several ligands with the highest docking scores toward the top docking side of the S-protein. Figure 3b–i shows the dockings of several ligands with the SARS-CoV-2 S-protein.
To understand whether these ligands are reasonably good inhibitors that block the interaction of the SARS-CoV-2 S-protein with ACE-2, the docking energy of the S-protein/ligand complex with ACE-2 was calculated and the results are listed in Table 3. For comparison, the docking energy between the SARS-CoV S-protein and ACE-2 was also calculated and the score was −92.7 kcal/mol, which was close to the −78.6 kcal/mol reported by Xu et al.’s work [4]. The docking energy between the SARS-CoV-2 S-protein and ACE-2 was calculated to be −82.2 kcal/mol, suggesting a slightly weaker interaction than that of the SARS-CoV S-protein with ACE-2. The observation is similar to that reported by Xu et al.’s work [4]. Table 3 shows that more than half of those ligands docking onto the SARS-CoV-2 S-protein do not significantly change the interaction of the SARS-CoV-2 S-protein with ACE-2—i.e., they are not inhibitors to block the interaction of the SARS-CoV-2 S-protein with ACE-2. However, ergotamine, amphotericin B, vancomycin, zafirlukast, and lanicor showed that once they were bound to the S-protein, the interactions of these complex with ACE-2 were no longer energetically favored interactions—i.e., these ligands acted as desired inhibitors that can efficiently block the interaction of the SARS-CoV S-protein with ACE-2. Among these, ergotamine and amphotericin b are most promising, since they demonstrate the highest docking energy to the SARS-CoV-2 S-protein (Table 2). Thus, they are strongly suggested as the core drugs for clinical trials to treat COVID-19 patients. Considering the severe and potentially lethal side effects of amphotericin b [28], ergotamine and vancomycin seem be the top choices.
Table 2.
Different docking scores of ligands for the SARS-CoV-2 S-protein.
| Drug Name | Docking Score (kcal/mol) | Usage |
|---|---|---|
| Ergotamine | −8.8 | For treatment of acute migraine attacks [29] |
| Amphotericin b | −8.3 | An antifungal medication used for serious fungal infections and leishmaniasis [30] |
| Indinavir | −8.1 | A component antiretroviral therapy to treat HIV/AIDS |
| Vancomycin | −7.7 | For treatment bacterial infections [31] |
| Lonpinavir | −7.7 | An antiretroviral, often used against HIV infections |
| Zafirlukast | −7.6 | For the chronic treatment of asthma |
| Lanicor | −7.5 | Used to treat heart conditions [32] |
| PubChem ID: 54098557 | −7.5 | -- |
| Digitaline Nativelle | −7.5 | For treatment of congestive heart failure, also used as angiotensin-converting enzyme (ACE) inhibitor |
| Rivaroxaban | −7.5 | To treat and prevent blood clots [33] |
| Tadalafil | −7.5 | To treat erectile dysfunction |
| Nelfinavir | −7.3 | The treatment of HIV |
| Montelukast | −7.2 | Treatment of asthma |
| Saquinavir | −7.1 | The treatment of HIV |
| Carfilzomib | −7.1 | Anti-cancer drug as proteasome inhibitor |
| Lapatinib | −7.0 | Anti-cancer drug |
| Atovaquone | −7.0 | To treat pneumocystis pneumonia, toxoplasmosis, malaria and babesia |
| Celecoxib | −7.0 | An anti-inflammation drug |
| Vardenafil | −6.9 | For treatment of erectile dysfunction |
| Dasatinib | −6.8 | To treat certain cases of chronic myelogenous leukemia |
| Cortisone | −6.6 | Can be used for a variety of conditions |
4. Discussion
Disulfiram, lopinavir, and ritonavir are the three approved and active protease inhibitors against SARS and MERS. Indeed, lopinavir and ritonavir were successfully used to treat a patient in Thailand in January 2020. Our results show that among these ligands, saquinavir, tadalafil, rivaroxaban, sildenafil, dasatinib, vardenafil, montelukast are most promising due to their higher docking scores (<−8.5 kcal/mol, which coresponds to < 1 μM IC50) than others. All of these scores appear better than that of the antiviral drug Lopinavir (−8.2 kcal/mol). As a comparison, the docking scores reported for lopinavir with the viral RNA polymerase is −8.3 kcal/mol [18]. It is a remarkable observation that some SARS-CoV-2 inhibitors such as indinavir could not block the interaction of the SARS-CoV-2 S-protein with ACE-2, while other inhibitors, such as ergotamine and amphotericin B, can effectively inhibit such interaction. This is somewhat confusing, since all of these three compounds dock on the same docking site that is marked by the red circles in Figure 3—the grove between an extended insertion that contains the β5/β6 strands and the receptor-binding motif (RBM) loop [28]. To comprehend what caused the significant difference, we overlaid the structures of the three docked compounds and ACE-2 on the SARS-CoV-2 spike protein in Figure 4. The comparison clearly shows that ergotamine (red) and amphotericin B (blue) extend further out toward the ACE-2 and thus effectively block the interaction of the SARS-CoV-2 spike protein with ACE-2 while indinavir (green) clings to the SARS-CoV-2 spike protein, leaving room for ACE-2 to interact with the spike protein.
5. Conclusions
For the inhibition of the SARS-CoV-2 3CL protease, saquinavir, tadalafil, rivaroxaban, sildenafil, dasatinib, vardenafil, and montelukast are most promising due to their high docking scores (<−8.5 kcal/mol), which were more negative than those of other ligands.
Among these that showed an excellent inhibiting property to block the interaction of SARS-CoV-2 S-protein with ACE-2 in Table 3, ergotamine, amphotericin b, and vancomycin are the most promising since they are also among the highest to bind to the SARS-CoV-2 S-protein, as shown in Table 2.
For more active results, a combination of 3CL protease inhibitors and ergotamine may be considered.
Author Contributions
Conceptualization, Z.Q., and H.-F.J.; methodology, Z.Q.; software, Z.Q.; validation, Z.Q., H.Z. and H.-F.J.; formal analysis, Z.Q.; writing—original draft preparation, Z.Q. and H.Z.; writing—review and editing, Z.Q., H.Z., Q.C., H.-F.J. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by NIH grant number P30 GM122732, H. Zhang was funded by NIH R01HL128895.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Lu, R.; Zhao, X.; Li, J.; Niu, P.; Yang, B.; Wu, H.; Wang, W.; Song, H.; Huang, B.; Zhu, N.; et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding. Lancet 2020, 20, 30251–30258. [Google Scholar] [CrossRef] [Green Version]
- Zhou, P.; Yang, X.; Wang, X.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.; Zhu, Y.; Li, B.; Huang, C.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Li, F. Structure, Function, and Evolution of Coronavirus Spike Proteins. Annu. Rev. Virol. 2016, 3, 237–261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Xu, X.; Chen, P.; Wang, J.; Feng, J.; Zhou, H.; Li, X.; Zhong, W.; Hao, P. Evolution of the novel coronavirus from the ongoing Wuhan outbreak and modeling of its spike protein for risk of human transmission. Sci. China Life Sci. 2020, 63, 457–460. [Google Scholar] [CrossRef] [Green Version]
- Chu, C.M.; Cheng, V.C.; Hung, I.F.; Wong, M.M.; Chan, K.H.; Chan, K.S.; Kao, R.Y.; Poon, L.L.; Wong, C.L.; Guan, Y.; et al. Role of lopinavir/ritonavir in the treatment of SARS: Initial virological and clinical findings. Thorax 2004, 59, 252–256. [Google Scholar] [CrossRef] [Green Version]
- Nukoolkarn, V.; Lee, V.S.; Malaisree, M.; Aruksakulwong, O.; Hannongbua, S. Molecular dynamic simulations analysis of ritonavir and lopinavir as SARS-CoV 3CL(pro) inhibitors. J. Theor. Biol. 2008, 254, 861–867. [Google Scholar] [CrossRef]
- Liu, X.; Zhang, B.; Jin, Z.; Yang, H.; Rao, Z. The crytal structure of 2019-nCoV main protease in complex with an inhibitor N3. Nature 2020. [Google Scholar] [CrossRef] [Green Version]
- Sanner, M.F. Python: A programming language for software integration and development. J. Mol. Graph. Model. 1999, 17, 57–61. [Google Scholar]
- 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. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [Green Version]
- Hanwell, M.D.; Curtis, D.E.; Lonie, D.C.; Vandermeersch, T.; Zurek, E.; Hutchison, G.R. Avogadro: An advanced semantic chemical editor, visualization, and analysis platform. J. Cheminform. 2012, 4. [Google Scholar] [CrossRef] [Green Version]
- Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 16, 2785–2791. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [Green Version]
- AndrejŠali, A.; Blundell, T.L. Comparative Protein Modelling by Satisfaction of Spatial Restraints. J. Mol. Biol. 1993, 234, 779–815. [Google Scholar]
- Laskowski, R.A.; MacArthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Cryst. 1993, 26, 283–291. [Google Scholar] [CrossRef]
- Zhu, K.; Day, T.; Warshaviak, D.; Murrett, C.; Friesner, R.; Pearlman, D. Antibody structure determination using a combination of homology modeling, energy-based refinement, and loop prediction. Proteins 2014, 82, 1646–1655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Salam, N.K.; Adzhigirey, M.; Sherman, W.; Pearlman, D.A. Structure-based approach to the prediction of disulfide bonds in proteins. Protein Eng. Des. Sel. 2014, 27, 365–374. [Google Scholar] [CrossRef] [Green Version]
- Beard, H.; Cholleti, A.; Pearlman, D.; Sherman, W.; Loving, K.A. Applying Physics-Based Scoring to Calculate Free Energies of docking for Single Amino Acid Mutations in Protein-Protein Complexes. PLoS ONE 2013, 8, e82849. [Google Scholar] [CrossRef] [Green Version]
- Li, J.; Abel, R.; Zhu, K.; Cao, Y.; Zhao, S.; Friesner, R.A. The VSGB 2.0 model: A next generation energy model for high resolution protein structure modeling. Proteins 2011, 79, 2794–2812. [Google Scholar] [CrossRef] [Green Version]
- Banks, J.L.; Beard, H.S.; Cao, Y.; Cho, A.E.; Damm, W.; Farid, R.; Felts, A.K.; Halgren, T.A.; Mainz, D.T.; Maple, J.R.; et al. Integrated Modeling Program, Applied Chemical Theory (IMPACT). J. Comput. Chem. 2005, 26, 1752–1780. [Google Scholar] [CrossRef] [Green Version]
- Kitchen, V.S.; Skinner, C.; Ariyoshi, K.; Lane, E.A.; Duncan, I.B.; Burckhardt, J.; Burger, H.U.; Bragman, K.; Pinching, A.J.; Weber, J.N. Safety and activity of saquinavir in HIV infection. Lancet 1995, 345, 952–955. [Google Scholar] [CrossRef]
- Baumann, M. An Overview of the Key Routes to the Best Selling 5-membered Ring Heterocyclic Pharmaceuticals. Beilstein J. Org. Chem. 2011, 7, 442–495. [Google Scholar] [CrossRef] [PubMed]
- Kakkar, A.K.; Brenner, B.; Dahl, O.E.; Eriksson, B.I.; Mouret, P.; Muntz, J.; Soglian, A.G.; Pap, A.F.; Misselwitz, F.; Haas, S.; et al. Extended duration rivaroxaban versus short-term enoxaparin for the prevention of venous thromboembolism after total hip arthroplasty: A double-blind, randomised controlled trial. Lancet 2008, 372, 31–39. [Google Scholar] [CrossRef]
- Terrett, N.K.; Bell, A.S.; Brown, D.; Ellis, P. Sildenafil (Viagra), a potent and selective inhibitor of type 5 cGMP phosphodiesterase with utility for the treatment of male erectile dysfunction. Bioorg. Med. Chem. Lett. 1996, 6, 1819–1824. [Google Scholar] [CrossRef]
- Talpaz, M.; Shah, N.P.; Kantarjian, H.; Donato, N.; Nicoll, J.; Paquette, R.; Cortes, J.; O’Brien, S.; Nicaise, C.; Bleickardt, E.; et al. Dasatinib in imatinib-resistant Philadelphia chromosome-positive leukemias. N. Engl. J. Med. 2006, 354, 2531–2541. [Google Scholar] [CrossRef] [Green Version]
- Kloner, R.A. Pharmacology and drug interaction effects of the phosphodiesterase 5 inhibitors: Focus on alpha-blocker interactions. Am. J. Cardiol. 2005, 96, 42M–46M. [Google Scholar] [CrossRef]
- De Lepeleire, I.; Reiss, T.F.; Rochette, F.; Botto, A.; Zhang, J.; Kundu, S.; Decramer, M. Montelukast causes prolonged, potent leukotriene D4-receptor antagonism in the airways of patients with asthma. Clin. Pharmacol. Ther. 1997, 61, 83–92. [Google Scholar] [CrossRef]
- Capaldini, L. Protease inhibitors’ metabolic side effects: Cholesterol, triglycerides, blood sugar, and “Crix belly. Interview with Lisa Capaldini, M.D. Interview by John S. James”. AIDS Treatment News 1997, 277, 1–4. [Google Scholar]
- Park, N.H.; Shin, K.H.; Kang, M.K. Chapter 34. Antifungal and Antiviral Agents. In Pharmacology and Therapeutics for Dentistry, 7th ed.; Elsevier: Amsterdam, The Netherlands, 2017; pp. 488–503. [Google Scholar]
- Ibraheem, J.J.; Paalzow, L.; Tfelt-Hansen, P. Low bioavailability of ergotamine tartrate after oral and rectal administration in migraine sufferers. Br. J. Clin. Pharmacol. 1983, 16, 695–699. [Google Scholar] [CrossRef] [Green Version]
- World Health Organization. Control of the Leishmaniasis: Report of a Meeting of the WHO Expert Committee on the Control of Leishmaniases; World Health Organization: Geneva, Switzerland, 22–26 March 2010; pp. 55–186. [Google Scholar]
- Liu, C.; Bayer, A.; Cosgrove, S.E.; Daum, R.S.; Fridkin, S.K.; Gorwitz, R.J.; Kaplan, S.L.; Karchmer, A.W.; Levine, D.P.; Murray, B.E.; et al. Clinical practice guidelines by the infectious diseases society of America for the treatment of methicillin-resistant Staphylococcus aureus infections in adults and children: Executive summary. Clin. Infect. Dis. 2011, 52, 285–292. [Google Scholar] [CrossRef]
- Wang, W.; Chen, J.S.; Zucker, I.H. Carotid sinus baroreceptor sensitivity in experimental heart failure. Circulation 1990, 81, 1959–1966. [Google Scholar] [CrossRef] [Green Version]
- Turpie, A.G. New oral anticoagulants in atrial fibrillation. Eur. Heart J. 2008, 29, 155–165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1.
The comparison of the crystal structures of the SARS-2 spike protein (6VW1, pink color) with our homology modeling (light brown color) using the SARS-2 template.
Figure 1.
The comparison of the crystal structures of the SARS-2 spike protein (6VW1, pink color) with our homology modeling (light brown color) using the SARS-2 template.
Figure 2.
Different ligands in the docking site of SARS-CoV-2 protease. (a,b) Tadalafil; (c,d) Saquinavir; (e,f) Dasatinib.
Figure 2.
Different ligands in the docking site of SARS-CoV-2 protease. (a,b) Tadalafil; (c,d) Saquinavir; (e,f) Dasatinib.
Figure 3.
Potential active site selection and ligand-receptor interaction. (a) The docking site (inside the red frame) was chosen between the ACE-2 (light brown color) and SARS-CoV-2 S-protein (Cyan color); (b,c) Amphotericin b docks onto the SARS-CoV-2 S-protein; (d,e) Saquinavir docks onto the SARS-CoV-2 S-protein; (f,g) Indinavir docks onto the SARS-CoV-2 S-protein; (h,i) Ergotamine binds onto the SARS-CoV-2 S-protein.
Figure 3.
Potential active site selection and ligand-receptor interaction. (a) The docking site (inside the red frame) was chosen between the ACE-2 (light brown color) and SARS-CoV-2 S-protein (Cyan color); (b,c) Amphotericin b docks onto the SARS-CoV-2 S-protein; (d,e) Saquinavir docks onto the SARS-CoV-2 S-protein; (f,g) Indinavir docks onto the SARS-CoV-2 S-protein; (h,i) Ergotamine binds onto the SARS-CoV-2 S-protein.
Figure 4.
An overlay of the modeled structures of amphotericin, ergotamin, indinavir, and ACE-2 with the SARS-CoV-2 spike protein.
Figure 4.
An overlay of the modeled structures of amphotericin, ergotamin, indinavir, and ACE-2 with the SARS-CoV-2 spike protein.
Table 1.
Different docking scores (binding affinities) of the tested drugs for SARS-CoV-2 proteinase.
Table 1.
Different docking scores (binding affinities) of the tested drugs for SARS-CoV-2 proteinase.
| Drug Name | Docking Score (kcal/mol) | Usage |
|---|---|---|
| Saquinavir | −9.5 | Antiretroviral drug to treat or prevent HIV/AIDS [20] |
| Tadalafil | −9.3 | A medication used to treat erectile dysfunction (ED), benign prostatic hyperplasia (BPH), and pulmonary arterial hypertension [21] |
| Rivaroxaban | −9.2 | An anticoagulant medication used to treat and prevent blood clots [22] |
| Sildenafil | −8.9 | A medication used to treat erectile dysfunction and pulmonary arterial hypertension [23] |
| Dasatinib | −8.8 | A targeted therapy used to treat certain cases of chronic myelogenous leukemia (CML) and acute lymphoblastic leukemia (ALL) [24] |
| Vardenafil | −8.7 | A PDE5 inhibitor used to treat erectile dysfunction [25] |
| Montelukast | −8.5 | To treat seasonal and year-round allergies [26] |
| Indinavir | −8.3 | A component antiretroviral therapy to treat HIV/AIDS [27] |
| Lopinavir | −8.2 | Protease inhibitor |
| Cortisone | −8.2 | Can be used for a variety of conditions |
| celecoxib | −8.1 | An anti-inflammation drug |
| Atazanavir | −8.1 | An antiretroviral drug used for HIV treatment |
| Iressa | −7.9 | A drug for cancer treatment |
| Darunavir | −7.7 | An antiretroviral drug used for HIV treatment |
| Sorafenib | −7.5 | A drug for cancer treatment |
Table 3.
Docking energy of the SARS-CoV S-protein with and without ligands to human ACE-2.
| Interaction of S-Protein and S-Protein/Drug Complex with ACE-2 | Docking Energy (kcal/mol) |
|---|---|
| SARS-CoV S-protein (for comparison) | −92.7 |
| SARS-CoV-2 S-protein | −82.2 |
| SARS-CoV-2 S-protein/Ergotamine | 56.4 |
| SARS-CoV-2 S-protein/Amphotericin b | 78.6 |
| SARS-CoV-2 S-protein/Indinavir | −61.9 |
| SARS-CoV-2 S-protein/Vancomycin | 81.7 |
| SARS-CoV-2 S-protein/Zafirlukast | 52.6 |
| SARS-CoV-2 S-protein/Lanicor | 4.2 |
| SARS-CoV-2 S-protein/Nelfinavir | −81.5 |
| SARS-CoV-2 S-protein/Montelukast | −71.3 |
| SARS-CoV-2 S-protein/Saquinavir | −48.2 |
| SARS-CoV-2 S-protein/Carfilzomib | −88.1 |
| SARS-CoV-2 S-protein/Lapatinib | −83.1 |
| SARS-CoV-2 S-protein/Atovaquone | −68.2 |
| SARS-CoV-2 S-protein/Celecoxib | −74.2 |
| SARS-CoV-2 S-protein/Dasatinib | −42.3 |
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Qiao, Z.; Zhang, H.; Ji, H.-F.; Chen, Q. Computational View toward the Inhibition of SARS-CoV-2 Spike Glycoprotein and the 3CL Protease. Computation 2020, 8, 53. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020053
AMA Style
Qiao Z, Zhang H, Ji H-F, Chen Q. Computational View toward the Inhibition of SARS-CoV-2 Spike Glycoprotein and the 3CL Protease. Computation. 2020; 8(2):53. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020053
Chicago/Turabian StyleQiao, Zhen, Hongtao Zhang, Hai-Feng Ji, and Qian Chen. 2020. "Computational View toward the Inhibition of SARS-CoV-2 Spike Glycoprotein and the 3CL Protease" Computation 8, no. 2: 53. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020053
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
