{'fraction_non_alphanumeric': 0.037669980383681884, 'fraction_numerical': 0.04066230009641919, 'mean_word_length': 4.711735662742119, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Virtual screening and molecular dynamics study of approved drugs as inhibitors of spike protein S1 domain and ACE2 interaction in SARS-CoV-2

Background The receptor binding domain (RBD) of spike protein S1 domain SARS-CoV-2 plays a key role in the interaction with ACE2, which leads to subsequent S2 domain mediated membrane fusion and incorporation of viral RNA into host cells. In this study we tend to repurpose already approved drugs as inhibitors of the interaction between S1-RBD and the ACE2 receptor. Methods 2456 approved drugs were screened against the RBD of S1 protein of SARS-CoV-2 (target PDB ID: 6M17). As the interacting surface between S1-RBD and ACE2 comprises of bigger region, the interacting surface was divided into 3 sites on the basis of interactions (site 1, 2 and 3) and a total of 5 grids were generated (site 1, site 2, site 3, site 1+site 2 and site 2+site 3). A virtual screening was performed using GLIDE implementing HTVS, SP and XP screening. The top hits (on the basis of docking score) were further screened for MM-GBSA. All the top hits were further evaluated in molecular dynamics studies. Performance of the virtual screening protocol was evaluated using enrichment studies. Result and discussion: We performed 5 virtual screening against 5 grids generated. A total of 42 compounds were identified after virtual screening. These drugs were further assessed for their interaction dynamics in molecular dynamics simulation. On the basis of molecular dynamics studies, we come up with 10 molecules with favourable interaction profile, which also interacted with physiologically important residues (residues taking part in the interaction between S1-RBD and ACE2. These are antidiabetic (acarbose), vitamins (riboflavin and levomefolic acid), anti-platelet agents (cangrelor), aminoglycoside antibiotics (Kanamycin, amikacin) bronchodilator (fenoterol), immunomodulator (lamivudine), and anti-neoplastic agents (mitoxantrone and vidarabine). However, while considering the relative side chain fluctuations when compared to the S1-RBD: ACE2 complex riboflavin, fenoterol, cangrelor and vidarabine emerged out as molecules with prolonged relative stability. Conclusion We identified 4 already approved drugs (riboflavin, fenoterol, cangrelor and vidarabine) as possible agents for repurposing as inhibitors of S1:ACE2 interaction. In-vitro validation of these findings are necessary for identification of a safe and effective inhibitor of S1: ACE2 mediated entry of SARS-CoV-2 into the host cell.


Introduction
The spike protein of coronavirus plays a major role in attachment to the host and further series of events e.g., membrane fusion and incorporation of the host RNA genome to the host cell [1,2]. The S1 domain of spike protein comprises of one N terminal domain (NTD) and three c terminal domain (CTD 1 to 3) and the RBD resides within the CTD-1 [3]. The interactions between the S1-RBD with the PD (peptidase domain) of ACE2 leads to attachment of the virus to the host [4,5]. which is followed by S2 domain mediated membrane fusion and internalization of viral ribonucleoprotien complex into the host cell [3].
Compared to SARS-CoV, many changes can be noted in the RBD region of SARS-CoV-2 e.g. at the N-terminus, the following variations, Arg426 to Asn439, Tyr484 to Gln498, and Thr487 to Asn501, at equivalent positions were observed between SARS-CoV and SARS-CoV-2 [4]. Similar findings are also reported by Ortega et al., 2020 [6]. On sequence homology comparison between the SARS-CoV and SARS-CoV-2 spike protein, only 76.04% similarity (SARS-CoV sequence taken from UniProtKB -P59594 (SPIKE_CVHSA) and SARS-CoV-2 sequence I.D GeneID: 43740568) was observed. Hence, there is an urgent need for the development of drugs targeting against SARS-CoV-2.
Although, the combination HCQ þ Azithromycin showed promise [7] in clinical settings for the treatment of 2019-nCoV in the first clinical study, a subsequent study showed negative results [8]. Here comes the need for rapid development of a safe and effective new drug or repurposing the already existing drugs against 2019-nCoV. For the purpose of our study, we have selected S1: ACE2 interaction because, being the first point of contact between the host and the virus, the Spike protein S1 domain represents a potential target for the development of anti-SARS-CoV-2 drugs which might have the potential to inhibit the entry of SARS-CoV-2 into the host cells [4].
A typical new drug development process takes 10e12 years [9]. However, in the settings of this current pandemic of COVID-19, we need a safe and effective drug at the earliest and thus repurposing seems to be the most effective strategy. Use of approved drugs for the repurposing activity carries the advantage that the preclinical, clinical toxicity information and kinetic data of these agents are already available and thus we can expect a shorter regulatory pathway [10]. In this current study, we have repurposed the approved drugs for finding a novel inhibitor of S1: ACE2 interaction, inspired from the theme of relatable studies from Sarma et al. [64]. We did not target the host protein ACE2, as its a major component of the human rennin angiotensin system (RAS) system and is important for maintaining of blood pressure and other important homeostatic mechanisms [11].
As there are only a few thousands of approved drugs are available, a preliminary in-vitro screening of all these agents for the identification of potential inhibitors may be time-consuming and costly. Here comes the need for computational approaches, which can help us in selection potential agents beforehand and help our selection of molecules for further evaluation in in-vitro, [10]. Although a previous paper has already addressed the similar theme [12], they however, took SARS-CoV spike protein (PDB: 2AJF_E) for their virtual screening studies. In this regard, we used the spike protein of SARS-CoV-2. Additionally, we validated the results of our experiments using molecular dynamics studies. So, this is the first study, which purely targeted SARS-CoV-2 spike protein S1 domain and ACE2 interaction for the identification of inhibitors.

Materials and methods
In-silico studies were performed on an ACER Predator Helios 300 laptop running on LINUX Ubuntu OS 18.04.02 LTS using Schrodinger Maestro version 2019e3 for docking and molecular dynamics simulation processes.

The target PDB structure selection
The target protein considered for this study is the Spike protein S1 specifically the RBD (receptor binding domain of S1 protein) as the key target. In the RCSB database [13], we could find 5 PDB IDs (till date 05-04-2020, which are 6M17, 6VSB, 6VW1, 6LZG and 6MOJ) which has a structure of these Spike protein S1 domain of 2019-nCoV.
Among these different structures, 6M17 was screened as it is a cryo-EM structure with a resolution of 2.9 Å, and had the full length of RBD region (229 amino acid residues) and the interaction details between S1 RBD and ACE2 [4]. The rationale for excluding the 6LZG, 6M0J proteins was the non-availability of details of interface interactions, 6VSB PDB was presenting only a small region of RBD which had no ACE-2 binding region and 6VW1 PDB structure was a chimeric structure. Therefore, PDB structure of PDB I.D. 6M17 was considered for this study [14].

Validation of target protein
The validity of the target protein structure was evaluated using PROCHECK tool from SAVES v5.0 server. PROCHECK is a useful tool to examine the stereochemical efficiency of target protein structure. By this server we analyzed the geometry of 'residues by residue' or overall residues and also determined the Ramachandran plot to check the quality of selected target protein model. This tool was helpful to assure the qulity of protein which were used as a model in our study (PDB: 6M17).

Preparation of the protein
The selected protein was prepared for further docking purposes. Protein preparation was performed using the "protein preparation wizard" tool of "Maestro" version 10.2 [62]. The missing side chain and loops were filled in pre-processingat OPLS3 force field minimization algorithm where water molecules beyond 5.0 Å to rule out the hindrence [63].

Retrieval and preparation of ligands
The DrugBank database is a free accessible, comprehensive Details of interaction between spike protein (showed in surface diagram) and ACE2 (showed in ribbon diagram). Three points of interactions can be seen between ACE2 and spike protein. Fig. 1A represents interaction between S1 and ACE2 (first view) and figure and 1.B. represents the same in from 180 from the first view.

Table 1
Roles of important amino acid residues (physiologically important residues) in the S1-RBD: ACE2 interactions.
C terminus of a1 of ACE2 The ACE2: a1 (C terminal) residues Gln 24 and Met 82 interact with Gln 474 and Phe 486 via H-bond and van-der waals forces [22]. Hotspot residues for recognition of host receptor ACE2 In SARS-CoV, two prime residues 479 and 487 of RBD-CoV are involved in the recognition of ACE2 [24,25]. In the case of SARS-CoV-2, the residues corresponding to N479 is Q493 and T487 is N501. Physiologically important changes in SARS-CoV-2 when compared to SARS-CoV The residues val 404 (in SARS-CoV) change into Lys417 in SARS-CoV-2, which may result in stronger association with ACE2. This stronger affinity may be attributed to formation of salt bridge between Lys417 of SARS-CoV-2 and ASP30 of ACE2 [22]. The residues Leu 472 change into Phe486 may shown the stronger van der Waals interaction at met82. How-ever, replacement of Arg426 to Asn439 appears to weaken the interaction by losing one important salt bridge with Asp 329 on ACE2 [22]. Capping loop residues [6] (Two capping loops in the binding domain, which stabilized the interaction with ACE2 by producing an increased amount of electrostatic interactions) database of drugs and drug targets, which also contains the approved small molecules, biologics, nutraceuticals and few thousand of experimental drugs. For the purpose of repurposing, we selected the approved small molecule set (2456 drugs, up to date 10-01-2020) from DrugBank database [18]. SDF files of these 2456 approved drugs were obtained from the drug bank database [19] and prepared using ligprep module. This preparation of ligand is a necessary step for further docking and molecular simulations. OPLS 2003 force field was used for the energy minimization process, which is helpful to produce the low-energy isomers and prepare the ligands for docking [20].

S1-RBD and ACE2 interaction details
The RBD region of SARS-CoV-2 comprises of a core structure and the RBM (receptor binding motif) [21]. The host recognition by SARS-CoV-2 depends upon the interaction between the claw structured peptidase domain of ACE2 and the RBM of RBD region of S1 of spike protein [22]. In SARS-CoV-2, the S1-RBD extends from amino acid residues 331 to 524 [23]. An extended loop region of the RBD spans the arch-shaped a1 helix of the ACE2-PD like a bridge, moreover, these contacts can be arranged into three distinct clusters to apprehend the interacting domain. At the N terminus of a1 helix, Gln498, Thr500, and Asn501 of the S1-RBD forms a network of hydrogen bonds network with Tyr41, Gln42, Lys353, and Arg357 of ACE2. In the mid region, the "bridge", Lys417 and Tyr453 of the S1-RBD can be observed interacting to Asp30 and His34 of ACE2, respectively. The C terminus of a1interacts with ACE2 through Hbond (Gln474 RBD ---Gln 24 ACE2 ), while the van der Waals forces were seen formed between Phe486 of S1-RBD and Met82 of ACE2 [22]. The details of interaction between S1 and ACE-2 are represented in Fig. 1 and Table 1.
The site of interactions can be divided into three sites: the portion of S-RBD interacting with of N terminus of alpha 1 of ACE2 (in our study, we denoted it as site 1), part of S1-RBD interacting with middle of the bridge region of PD of ACE2 (in our study, we denoted it as site 2) and part of S1-RBD interacting with the c terminus of Alpha 1 of PD of ACE2 (in our study, we denoted it as site 3). Details of site 1, site 2 and site 3 are showed in Fig. 1.

Generation of grid
It is the very important and essential step for the binding of ligands to the receptor. The three dimensional boundary were generated around actively participating residues on RBD (which are involved in S1-RBD and ACE2 interactions) for ligand binding at desired site by using of "receptor grid generation" module of GLIDE, Schr€ odinger tool [26].
Site 1: The important interacting residues in the RBD Gln498, Thr500, and Asn501 form a network of hydrogen bonds (Hbonds) with Tyr41, Gln42, Lys353, and Arg357 from ACE2s. Details are showed in Fig. 1. Site 2: The middle regions of RBD have important residues that are Lys417 and Tyr453 interact with Asp30 and His34 of ACE2, respectively. Details are showed in Fig. 1. Site 3: The important residues in third region are Gln474, Phe486 form an interaction with GLN24 and Met82 of ACE2 respectively [22]. Details are showed in Fig. 1.
For the purpose of virtual screening, we made 5 grids (grid ¼ site 1, grid 2 ¼ site 2, grid 3 ¼ site 3, grid 4 ¼ Site 1 þ site 2 and grid , considering their importance in formation of S1-RBD: ACE2 interface and were used as distinct target sites in virtual screenings [ Fig. 2 (A to E)].

Virtual screening
The virtual screening was performed using "Glide" tool of "Maestro". A total of 2456 approved drugs were screened in a tandem HTVS/ SP/XP progression [16]. For virtual screening, a total of 5 parallel virtual screening was performed against the 5 grids generated.

MM-GBSA: Binding free energy calculation
The binding free energy between the ligand and receptor target, which were present in the docked structure, can be measured by using Prime MM-GBSA. Maestro, Schr€ odinger, LLC, New York, NY, USA which computes end-point binding free energy utilizing the receptor -ligand complex pose files [22e24]. The ligand-bound pose-viewer complexes were subjected to MM-GBSA in VSGB solvation within 5.0 Å from flexible residues at OPLS3 force field [29].

ROC curve plot and enrichment study
The performance of the virtual screening protocol was evaluated by performing enrichment studies using parameters area under the receiver operator characteristics curve (ROC) and Boltzmann enhanced discrimination of receiver operating characteristic (BEDROC). This test differentiates whether the used docking protocol can differentiate between actives and decoys. Y axis represents the sensitivity and X axis represents the 1-specificity. An AUC of 1 represents best performance of the virtual screening protocol and a value of 0 signifies worst performance. As five grids were generated for the purpose of our study, we evaluated the performance of each of the of these protocols by performing enrichment studies using actives and decoy molecules (Schrodinger decoy set, n ¼ 1000) [31].

Molecular dynamics
The Desmond module was used to carry out MD simulations of the S1 RBD-drug complexes and apo forms. Firstly, water model were developed and then sodium and chloride ions were added to neutralize the ligand protein complex by using "system builder" tool of Desmond, Schr€ odinger. Liquid simulations optimization was minimized at OPLS3 force field. After creating frames for 100ps interval, the complexes were submitted for 100ns simulations at NPT ensemble at 310 K, where N lied between 34,900e35,000 in the solvation box. Trajectory analysis for the trajectory outputs extracted C a RMSD, side chain RMSF, ligand contact maps and binding profile [15]. We performed 50 ns molecular dynamics studies of all the ligand-target complexes.

Validation of target protein
The stereochemical quality of selected protein structure was checked by PROCHECK method. SAVES v5.0 server was used check the model quality of the 3D structure of target protein PDB: 6M17. This server has PROCHECK tool which helps investigate the obtained crystal structure model of protein. Verifying the 3D structure of 6M17 showed that 87.89% of the residues have averaged 3D-1D score!0.2 (At least 80% of the amino acids have scored!0.2 in the 3D/1D profile), the qulity of model were assessed by Ramachandran plot which revealed that 83.5% of the residues were in the most favored region and 15.3% were in the allowed region. Overall G factor was À0.05. The maximum deviation was 9.2 and planar groups: 100.0% within limits (Fig. 3).

Virtual screening and MM-GBSA results
A total of five virtual screening was performed against the five grids generated. Ligands used were FDA approved drugs. The top performing ligands on the basis of docking score are shown in Table 2.
In our study we have generated 5 grids in target protein and after then virtual screening process were performed. The top XP docking result (on the basis of docking score) were found form each grid. All the ligands which showed good docking score against the respected grids are showed in Table 2.
In case of Grid 1, highest docking score is À9.948 (Acarbose), in grid site 2 have highest docking score is À6.324 of Cangrelor, grid 3 have highest docking score is À9.889 (pentethinin), grid 4 (S1þS2) and grid 5 (S2þS3) were shown highest docking score is À7.488 (diosmin) and À6.375 (norepinephrine) respectively. The results of docking score and glide score of the top performing ligands are showed in Table 2.
In prime MM-GBSA, in grid 1, lowest MM-GBSA score was À58.173 and the highest score was showed À29.972, in the case of grid 2, the lowest MM-GBSA score was À57.793 and highest score was-36.434, In grid 3, lowest MM-GBSA score was À79.537 highest score was À22.296, in the Grid 4 (S1þS2) lowest MM-GBSA score was À59.905 and the highest score was À29.969 and in case of grid 5 (S2þS3), lowest MM-GBSA score was À78.259 and highest score was À18.563. Data showed in Table 2.

Enrichment studies
Overall, the XP docking protocol showed good AUC-ROC with range being 0.71 to 0.97 (for all the grids). For Grid 1, AUC was 0.94 and BEDROC was 0.813 (at alpha ¼ 160.9), Grid 2, AUC was 0.94 and BDROC was 0.946 (at alpha ¼ 160.9). In case of Grid 3, AUC was 0.97 and BEDROC was 0.842 at (at alpha ¼ 160.9). In case of grid 4, AUC was 0.71 and BEDROC was 0.497 (at alpha ¼ 160.9) and in case of Grid 5, AUC was 0.94 and BEDROC was 0.870 (at alpha ¼ 160.9). (Table 3).

Molecular dynamics studies
We performed molecular dynamic simulation experiment at 50ns of all the ligands which showed good docking profile in virtual screening. Data shown in Table 4.
First we performed molecular dynamic simulation study of the apo form of targeted protein to evaluate and analyze the stability of the protein. Here we took two apo forms, first is, S1-RBD: ACE2 complex, and the second is, only RBD region of spike protein (S1-RBD). In our study the RBD region of S1 protein ranged from Cys 336 to Leu 518. We conducted a 100 ns MD simulation study for both the apo forms. The average RMSD value was found to be 2.4 Å for the ACE2: RBD Apo form and 2.7 Å found for the Apo S1e RBD protein. Data shown in Figs. 4 and 5. When comparing the RMSF value of the S1-RBD apo and the S1-RBD: ACE2 complex apo form, we could find that higher fluctuations were noted in the S1-RBD apo, especially in residue index 0e25, 75e100, 125e150 and 150e175.

Grid 1
Among all the ligands evaluated in molecular dynamics studies, favourable interaction profile was seen in cases of acarbose, riboflavin, cangrelor, epinephrine, Kanamycin and fenoterol. However, among these, only acarbose, riboflavin, cangrelor, Kanamycin and fenoterol showed interaction with physiologically important residues. Data of two best performers in terms of RMSD and RMSF value (riboflavin and fenoterol) are showed in Figs. 4 and 5. The RMSF values of riboflavin and fenoterol were found to be 1.2 Å and 0.8 Å. However, we couldn't find much difference in fluctuation between the S1-RBD: ACE2 apo form and the ligand bound complex.

Grid 2
Among all the ligands evaluated in molecular dynamics studies, favourable interaction profile was seen in cases of cangrelor, amikacin, levomefolic acid and fenoterol. However, among these, only cangrelor, amikacin and levomefolic acid showed interaction with physiologically important residues. Among these, in terms of RMSD and RMSF, cangrelor was found to be the best performer in terms of retaining the site. Data showed in Figs. 4 and 5. Compared to the S1-RBD: ACE2 apo form, the cangrelor bound form showed less fluctuation in terms of side chain RMSF values.

Grid 3
Among all the ligands evaluated in molecular dynamics studies, favourable interaction profile was seen in cases of pantethinone, lamivudine, azacitidine and pirbuterol. However RMSF values were higher in case of pantethinone, azacitidine and pirbuterol. Lamivudine showed interaction with physiologically important residue Phe486, however, the RMSD value was 3.6 and a higher fluctuation (RMSF) was noted when compared to the S1-RBD: ACE2 complex apo form. Data showed in Table 4.

Grid 4
Among all the ligands evaluated in molecular dynamics studies, favourable interaction profile was seen in cases of diosmin, mitoxantrone, amikacin, cangrelor and vidarabin. However interaction with physiologically important residue was seen cases of mitoxantrone, amikacin, cangrelor and vidarabin. Among these vidarabine performed better in terms of RMSD (2.4 Å) and RMSF (1 Å) result and fluctuations were less in terms of RMSF vale compared to the S1-RBD: ACE2 apo form. Data showed in Table 4 and Figs. 4 and 5.

Grid 5
Among all the ligands evaluated in molecular dynamics studies, although interaction was seen in cases of nor-epinephrine and sodium ascorbate, however, RMSF value was high for both the cases. Data showed in Table 4.
So, while combining the results of all the experiments, we come up with 10 molecules with favourable interaction profile, which also interacted with physiologically important residues. These are antidiabetic (acarbose), vitamins (riboflavin and levomefolic acid), anti-platelet agents (cangrelor), aminoglycoside antibiotics (Kanamycin, amikacin) bronchodilator (fenoterol), immunomodulator (lamivudine), and anti-neoplastic agents (mitoxantrone and vidarabine). However, while considering the relative fluctuation when compared to the S1-RBD: ACE2 complex riboflavin, fenoterol, cangrelor and vidarabine were the best performers in terms of relative stability. Details of interactions of these ligands with the target protein are showed in Fig. 6.
Acarbose: Acarbose is an antidiabetic. Acarbose is effective in delaying or preventing the progression of prediabetes to type 2 diabetes mellitus (T2DM) [32] In our study, acarbose showed interaction with write all residues Asn501, Gly 496, Gln 493, Ser 494, tyrosine 505 and Gly 502. Hydrogen bonding was seen in cases of Gln 493, Ser 494 and Asn501 Among these, residue Asn501 and Gln493 are physiologically important for the S1: ACE2 interaction. The residue Asn 501 of S1 interacts with the N terminus of alpha 1 of ACE2 and plays an important role in host receptor recognition by the virus [22] and both Asn501 and Gln493 are host-spot residues for host recognition by the virus [24,25]. The 501 residue plays a critical part in human infectivity. One study were predicted that a single N501T mutation (corresponding to the S487T mutation in SARS-CoV) may significantly enhance the binding affinity between 2019-nCoV RBD and human ACE2, highlighting the importance of the residue [21]. Data shown in Table 4.
Riboflavin: Riboflavin is also known as vitamin B2. Riboflavin showed interaction with residues Gln 498, Tyr 505, Arg 403, Tyr 495 and Ser 494. Among these, hydrogen bond formation was seen in cases of Tyr 505, Gln 498, Tyr 495 and Ser 494 (Table 4 and Fig. 6). However, among these only Gln 498 plays a strong physiological role in the interaction S1 and ACE2. Interaction between Gln 498 of S1 and N terminus of a1 of ACE2 plays important role in host receptor recognition and further sequelae [22]. Validity of the findings from our study is also highlighted by the fact that in in-vitro level, in human plasma and whole blood, use of riboflavin was associated with reduction of titre of SARS-CoV-2 to below detection level [33]. This highlights the validity of our methodology. In vero cell plaque assay also, use of riboflavin was associated with reduction of the titre of SARS-CoV-2 in both plasma and platelet products to below the limit of detection in tissue culture [34].
Levomefolic acid: Levomefolic acid is a biologically active form of folic acid [35]. In our study, Levomefolic acid showed interactions with residues of Glu406, Ile 418, Lys417, Tyr453 and H bond interaction was seen with Glu406, Lys417, Tyr453 and Ser494. However, among these only Lys 417 and Tyr 453 plays a strong physiological role in the interaction S1 and ACE. In clinical studies, low serum folate levels were associated with a more severe disease [36]. Homocystine is a predictor of radiological progression of COVID-19 [37] Folic acid supplementation also reduces the level of homocystine [38e40]. In in-silico studies, folic acid was also found to bind to active site of coronavirus main proitease [41]. Another multi target based study also indicated that folic acid may be useful in case of COVID-19 [42]. Folic acid also decreases the level of homocystine in clinical studies [43]. Hyper-homocystinemia is associated with more severe disease in case of COVID-19 [44].
Cangrelor: Cangrelor is an adenosine triphosphate analogue that selectively and specifically blocks P2Y 12 receptor-mediated platelet activation [45]. In this study, we found that cangrelor showed interaction with Arg403, Tyr495, Glu406, Lys417, Arg408, Gln409, Thr415 residues of RBD region of CoV-2 (Table 4 and Fig. 6). However, among these, H bond interaction was seen with Lys417, Gln409, Arg408 residues. However, only Lys 417 residue has a physiological important role in interaction S1 domain and ACE2. Cangrelor is an antiplatelet agent. Coagulopathy is reported in case of COVID-19 and this is a major cause of mortality and morbidity [46]. Another anti-platelet agent "aspirin" is already under evaluation for its protective action against COVID-19 [47]. In this context, being an anti-platelet agent, the effect of cangrelor will be of interest.

Kanamycin
Kanamycin, an aminoglycoside class antibiotic is a widely used antibiotic available as intra-muscular, intra-venous and oral formulation [48]. It is a commonly used antibiotic, which is used in many difficult to treat cases also [49]. Regarding antiviral activity, Kanamycin showed antiviral activity in in-vitro against HSV-2 [50], influenza-A virus [51], respiratory syncytial virus [52] and against foot and mouth disease viruses [53]. Important side effects are nephrotoxicty, ototoxicity and vestibular toxicity and toxicity depends upon cumulative life-time exposure [54]. In our study, Kanamycin showed interaction with physiological important residue Glu484 of S1-RBD. This residue plays important role in the interaction between S1 and ACE2.
Amikacin: Amikacin is an antibiotic of aminoglycoside class [55]. In our study amikacin showed interaction with Tyr495, tyr453, arg403, lys417 and glu406 residues of RBD region. H bond was seen with residues Arg403, Glu406 and Tyr453 of S1 RBD. Among all interaction, physiologically important residues showing interaction with amikacin are Lys417, Tyr453. Amikacin is a semisynthetic aminoglycoside antibiotic that is derived from kanamycin-A [56].
Fenoterol: In our study, fenopterol (inhaled bronchodilator) showed significant interactions with residues Tyr449, tyr495, tyr505 and arg403 of S1 RBD. Among these, H-bond interaction was seen with tyr495 and arg403 residues (Fig. 6). Among all only Tyr449 residue were shown the physiological important residue. Fenoterol is sympatho-mimetic agent. Fenoterol stimulates the beta-2-adrenergic receptors in the lungs, which in turn activates the enzyme adenylyl cyclase which then catalyzes the conversion of ATP to cAMP. The increased cAMP levels relaxes bronchial smooth muscle, relieve bronchial spasms and decline the inflammatory cell mediator release, principally from the mast cells [57].
Lamivudine: Lamivudine have shown intraction with Asn481, val 483, phe486 residues of S1 RBD and H bond intraction was seen with Asn481, gly 482 and phe486 residues. However, among the interacting residues, only phe486 residue was physiologically important in the interaction between S1 RBD and ACE2. Regarding antiviral effect of lamivudine, its a known anti-retroviral agent and also active against hepatitis B virus [58].
Mitoxantrone: It showed interaction with Arg403, Gln498 and Asn501 residue and all residue formed H bond with Mitoxantrone. The physiological important residues among them are gln498, asn501. Mitoxantrone is an anticancer agent and immune suppressor agent used for treating cancer and multiple sclerosis. Interestingly, this drug has also been reported as a potential main protease binder [59]. It also inhibit vaccinia virus replication by blocking virion assembly [60].
Vidarabine: In our study vidarabine have shown H-bond interaction with Gly496, Gln498, and Tyr505 residue of S1-RBD region. Among all interaction, physiological important residues interaction was shown with Gln498 ( Fig. 6). One of recent study were found that vidarabine were useful to overcome the SARS-CoV-2 [65].

Conclusion
Among all the agents, Fenoterol and riboflavin ligand at Grid 1 (site S1), Cangrelor at grid 2 (site 2) and vidarabine at grid 4 (site 1 þ 2) showed good interaction with RBD region of S1. We need invitro proof of concept studies of these identified agents to find out a potential inhibitor of S1: ACE2 interaction.

Declaration of competing interest
The authors declare that they have no known competing financialinterestsor personal relationships that could have appeared to influence the work reported in this paper.