key: cord-0887975-xxkagj6e
authors: Qiao, Zhen; Zhang, Hongtao; Ji, Hai-Feng; Chen, Qian
title: Computational View toward the Inhibition of SARS-CoV-2 Spike Glycoprotein and the 3CL Protease
date: 2020-05-31
journal: Computation (Basel)
DOI: 10.3390/computation8020053
sha: 0a647b38d4d6843432c6bc6afe1b8b1002756c04
doc_id: 887975
cord_uid: xxkagj6e

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.

As of 24 May 2020, over 5 millions people in the world have been confirmed as having the 2019 novel coronavirus disease , 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 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/). 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 Sprotein [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 (3CL pro ) 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 3CL pro .

With these models and crystal data, we performed in silico studies of potential inhibitors of the SARS-CoV-2 S-protein and 3CL pro .

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] .

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.

For the SARS-CoV-2 3 CL 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 (3 CL ) 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 Sprotein 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.

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 Sprotein/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. Table 1 shows the binding affinity of several ligands with SARS-CoV-2 3 CL 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.

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.

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 Jan. 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.

For the inhibition of the SARS-CoV-2 3 CL 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 Sprotein, as shown in Table 2 .

For more active results, a combination of 3 CL protease inhibitors and ergotamine may be considered.

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. An overlay of the modeled structures of amphotericin, ergotamin, indinavir, and ACE-2 with the SARS-CoV-2 spike protein. 

Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding

A pneumonia outbreak associated with a new coronavirus of probable bat origin

Function, and Evolution of Coronavirus Spike Proteins

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

Role of lopinavir/ritonavir in the treatment of SARS: initial virological and clinical findings

Hannongbua S Molecular dynamic simulations analysis of ritonavir and lopinavir as SARS-CoV 3CL(pro) inhibitors

The crytal structure of 2019-nCoV main protease in complex with an inhibitor N3

A programming language for software integration and development

Chimera-A visualization system for exploratory research and analysis

Avogadro: An advanced semantic chemical editor, visualization, and analysis platform

Automated docking with selective receptor flexibility

Improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading

AndrejŠali A; Blundell TL Comparative Protein Modelling by Satisfaction of Spatial Restraints

A program to check the stereochemical quality of protein structures

Antibody structure determination using a combination of homology modeling, energy-based refinement, and loop prediction

Structure-based approach to the prediction of disulfide bonds in proteins

Physics-Based Scoring to Calculate Free Energies of docking for Single Amino Acid Mutations in Protein-Protein Complexes

Friesner RA The VSGB 2.0 model: A next generation energy model for high resolution protein structure modeling

Integrated Modeling Program, Applied Chemical Theory (IMPACT)

Safety and activity of saquinavir in HIV infection

Overview of the Key Routes to the Best Selling 5-membered Ring Heterocyclic Pharmaceuticals

Extended duration rivaroxaban versus short-term enoxaparin for the prevention of venous thromboembolism after total hip arthroplasty: A double-blind, randomised controlled trial

Viagra), a potent and selective inhibitor of type 5 cGMP phosphodiesterase with utility for the treatment of male erectile dysfunction

Dasatinib in imatinib-resistant Philadelphia chromosome-positive leukemias

Kloner RA Pharmacology and drug interaction effects of the phosphodiesterase 5 inhibitors: Focus on alpha-blocker interactions

De Lepeleire I

Decramer M Montelukast causes prolonged, potent leukotriene D4-receptor antagonism in the airways of patients with asthma

Capaldini L Protease inhibitors' metabolic side effects: cholesterol, triglycerides, blood sugar, and "Crix belly

Antifungal and Antiviral Agents In Pharmacology and Therapeutics for Dentistry

Low bioavailability of ergotamine tartrate after oral and rectal administration in migraine sufferers

Control of the leishmaniasis: Report of a meeting of the WHO Expert Committee on the Control of Leishmaniases; World Health Organization

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

Zucker IH Carotid sinus baroreceptor sensitivity in experimental heart failure

Turpie AG New oral anticoagulants in atrial fibrillation

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

SARS-CoV-2 S-protein/Indinavir −61

SARS-CoV-2 S-protein/Zafirlukast 52

SARS-CoV-2 S-protein/Lanicor

SARS-CoV-2 S-protein/Nelfinavir −81

SARS-CoV-2 S-protein/Saquinavir −48

SARS-CoV-2 S-protein/Carfilzomib −88

SARS-CoV-2 S-protein/Lapatinib −83

SARS-CoV-2 S-protein/Atovaquone −68

SARS-CoV-2 S-protein/Celecoxib −74