Home Browse In silico screening of known small molecules to bind ACE2 specific...
ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
▬
✚
Research Article

In silico screening of known small molecules to bind ACE2 specific RBD on Spike glycoprotein of SARS-CoV-2 for repurposing against COVID-19

[version 1; peer review: 2 approved]
Bharath BR
https://orcid.org/0000-0002-7405-0268
1, Hrishikesh Damle2, Shiban Ganju2, Latha Damle1
Bharath BR
https://orcid.org/0000-0002-7405-0268
1, Hrishikesh Damle2, Shiban Ganju2, Latha Damle1
PUBLISHED 01 Jul 2020
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

This article is included in the Cheminformatics gateway.

This article is included in the Coronavirus collection.

Abstract

Background: Human coronavirus (SARS-CoV-2) is causing a pandemic with significant morbidity and mortality. As no effective novel drugs are available currently, drug repurposing is an alternative intervention strategy. Here we present an in silico drug repurposing study that implements successful concepts of computer-aided drug design (CADD) technology for repurposing known drugs to interfere with viral cellular entry via the spike glycoprotein (SARS-CoV-2-S), which mediates host cell entry via the hACE2 receptor.
Methods: A total of 4015 known and approved small molecules were screened for interaction with SARS-CoV-2-S through docking studies and 15 lead molecules were shortlisted. Additionally, streptomycin, ciprofloxacin, and glycyrrhizic acid (GA) were selected based on their reported anti-viral activity, safety, availability and affordability. The 18 molecules were subjected to molecular dynamics (MD) simulation.
Results: The MD simulation results indicate that GA of plant origin may be repurposed for SARS-CoV-2 intervention, pending further studies.
Conclusions: Repurposing is a beneficial strategy for treating COVID-19 with existing drugs. It is aimed at using docking studies to screen molecules for clinical application and investigating their efficacy in inhibiting SARS-CoV-2-S. SARS-CoV-2-S is a key pathogenic protein that mediates pathogen-host interaction. Hence, the molecules screened for inhibitory properties against SARS-CoV-2-S can be clinically used to treat COVID-19 since the safety profile is already known.

Keywords

SARS-CoV-2, Pathogen host interaction, SARS-CoV-2-Spike Glycoprotein, ACE2, Small Molecules

Corresponding author: Bharath BR
Competing interests: Latha Damle is the founder of Atrimed Biotech LLP and holds equity in Atrimed Pharmaceuticals. Shiban Ganju and Hrishikesh Damle hold equity shares in Atrimed Pharmaceuticals.
Grant information: The study was sponsored by Atrimed Pharmaceuticals Pvt. Ltd and no external funding was received for the study.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Copyright:  © 2020 BR B et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
How to cite: BR B, Damle H, Ganju S and Damle L. In silico screening of known small molecules to bind ACE2 specific RBD on Spike glycoprotein of SARS-CoV-2 for repurposing against COVID-19 [version 1; peer review: 2 approved]. F1000Research 2020, 9:663 (https://doi.org/10.12688/f1000research.24143.1) First published: 01 Jul 2020, 9:663 (https://doi.org/10.12688/f1000research.24143.1) Latest published: 01 Jul 2020, 9:663 (https://doi.org/10.12688/f1000research.24143.1)

Introduction

The complete genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is 82% identical to SARS-CoV, both viruses share a common clade encompassing the genus Betacoronavirus as the root node1,2. Currently no novel antivirals exist that are effective against either of the viruses3–6.

Drug repurposing is a commercially viable strategy, as it exploits existing drugs, thus significantly reducing the cost and time involved in developing effective therapeutics7–9. Experimental approaches, however, at pre-clinical and clinical stages for drug repurposing involve high cost and time10. Computational approaches can offer quick, considerable, and novel testable hypotheses for systematic drug repositioning8.

Current drugs in different phases of clinical trials are being investigated for inhibitory activity against viral targets that play a significant role in the coronavirus infection lifecycle. The drug targets might be involved in entry into the host (e.g. umifenovir and chloroquine), replication (e.g. lopinavir/ritonavir), or RNA synthesis (e.g. remdesivir/favipiravir). Among these, targeting SARS-CoV-2 cellular entry via the spike glycoprotein (SARS-CoV-2-S) has emerged as the leading option for repurposing9. As SARS-CoV-2-S is a surface protein involved in adhesion/fusion and entry into host cells, it has been identified as a potential drug target for both biologics and small molecules10.

The entry of COVID-19 pathogen is mediated by the homotrimeric transmembrane protein SARS-CoV-2-S. It is comprised of two functional subunits, S1 and S2, which are non-covalently bound in the pre-fusion conformation. The S1 subunit interacts with the human ACE2 receptor through the receptor binding domain (RBD), while the S2 subunit is one of the components of viral envelope11–19.

Apart from interacting with the ACE2 receptor, the RBD site also contributes to the stabilisation of the prefusion state of the S2 subunit equipped with fusion machinery18–24. In CoVs, the S-protein is cleaved by host proteases at the S20 site located above the fusion peptide16,25. This activates the protein via extensive irreversible conformational changes11,16,17,23,26. It is well understood that the entry of CoV into the susceptible host is a complex process that requires the vigorous actions of receptor binding and proteolytic processing of the S-protein to promote fusion with the pathogen27.

Hence, the current study aims to predict and validate the structure of SARS-CoV-2-S protein using computer-aided homology modelling tools and screen a library of small molecules for their interaction with the SARS-CoV-2-S protein.

Methods

Sequence analysis

The whole genome of SARS-CoV-2 (GenBank accession number: MT159721.1, length: 29882 bp) was retrieved from NCBI and used as a query to perform a sequence similarity search using NCBI-BLAST28. The BLAST search revealed 96.04%, 91.64% and 82.30% identity with Bat coronavirus RaTG13 (GenBank accession number: MN996532.1, length: 29855 bp), Pangolin coronavirus isolate MP789 (GenBank accession number: MT084071.1, length: 27213 bp) and SARS-CoV (GenBank accession number: JX163927, length: 29646bp), respectively.

Topological analysis of pathogen-host interactome for target validation

Drug target identification and validation were carried out using a network-based topological analysis method using the web-based application Pathogen-Host Interaction Search Tool (PHISTO) by setting pathogen type to virus, family to coronaviridae, species to SARS-Related coronavirus and strain to SARS-Cov. The node properties like the degree of connectivity (k) and betweenness centrality (BC) were assessed29,30. The statistical significance of k and BC values were assessed by the Fligner-Killeen (median) test.

Prediction of ligand binding site

The similarity between RBD domains of S-protein from SARS-CoV-2 (accession number QII57328, length:1273aa), SARS-CoV (accession number: AFR58728, length:1255aa) and RatG13 (accession number: QHR63300, length:1269aa) was evaluated by using the multiple sequence alignment (MSA) tool Clustal Omega from EMBL-EBI. Conservation in ACE2 receptor interaction was seen among all the three sequences aligned. This conservation aided in the active binding site prediction.

The protein-protein interaction between SARS-CoV-2-S and host ACE2 receptor complex was studied using the crystal structure from Protein Data Bank (PDB ID: 6CS2). The amino acids involved in the interaction were identified as ligand binding sites for inhibitor molecules.

Homology modelling of SARS-CoV-2-S protein

Homology modelling was performed with SWISS-MODEL for the protein sequence of SARS-CoV-2-S using the crystal structure of SARS-CoV-S and ACE2 complex (PDB ID: 6ACD) as a template. The modelled protein was validated for quality using a Ramachandran plot and prepared for molecular docking studies using the Protein Preparation Wizard feature of the Schrodinger Small Molecule Suite31. This analysis could also have been performed using open source software such as AutoDock32 or SwissDock33.

Protein preparation

The modelled receptor was processed for docking studies by deleting crystallographic water molecules with less than three H-bonds. This could also be done manually by editing the .PDB file in a text editor. Next, hydrogen atoms corresponding to neutral pH were added in consideration of ionisation states of amino acids. Following this, coordinates for any missing side-chain atoms were added using Prime v4.0, Schrödinger 2019-234. Finally, the energy of the modelled structure was minimised using the OPLS_2005 force field35. This analysis could also have been performed using open source software such as AutoDock32.

Ligand preparation

The three-dimensional conformations of the 4015 small molecule drugs already in use to treat various diseases and as nutritional supplements were downloaded from the DrugCentral database and subjected to ligand minimisation using Ligprip (LigPrep, version 2.3, Schrödinger, LLC, New York, NY, 2009). This analysis could also have been performed using open source software such as AutoDock32. The compounds were minimised by assigning force field OPLS_2005 and stereoisomers were calculated after retaining specific chiralities. The absorption, distribution, metabolism and excretion (ADME) predictions were performed for all ligands using the QikProp package36. This analysis could also have been performed using open source software such as SWISS-ADME37.

Molecular docking

The active site on the prepared receptor was defined around the selected residues (Arg426, Tyr436, Pro462, Thr486, Gly488, and Tyr491) with a 10Å radius. This generated a grid box measuring 20X20X20Å. The docking of small molecules over SARS-CoV-2-S was performed using Glide v7.8, This analysis could also have been performed using open source software such as AutoDock32. Schrödinger 2019-238 in different modes sequentially with defined and incremental precision, and computational time differences. The best-docked conformer with minimum Glide energy and E model energy was selected and lowest-energy docked complex of three known molecules streptomycin, ciprofloxacin, and glycyrrhizic acid (GA) in complex with SARS-CoV-2-S were selected for molecular dynamic simulations.

Molecular dynamics (MD) simulations

The MD39 of shortlisted complexes were studied using the OPLS_2005 force field40 in a plane TIP3P water model41. MD simulations were performed using Desmond version 4.242. This analysis could also have been performed using open source software such as GROMACS43. The system was built by dissolving the streptomycin/SARS-CoV-2-S, ciprofloxacin/SARS-CoV-2-S, and GA/SARS-CoV-2-S complexes in an orthorhombic box containing water molecules, allowing a buffer region of 10Å between atoms and box peripherals. The system was further minimised using the L-BFGS algorithm for a minimum of 10 steepest descent steps and a maximum of 2000 iterations until a gradient threshold of 25 kcal/mol/Å and convergence threshold of 1.0 kcal/mol/Å was reached. For short-range electrostatic interactions, the solid-phase microextraction44 method was employed at 1e-09 tolerance and 9Å cut-off radius. The built systems were gradually warmed up to 300K in the NPT ensemble with a time step of 2fs. A 100ns MD simulation in the NPT ensemble was performed using a Nose–Hoover thermostat40. Resulting root mean square deviation (RMSD) and root mean square fluctuation (RMSF) values were analysed.

Results and discussion

Multiple sequence alignment of complete genomes

The complete genome of coronavirus SARS-CoV-2, Bat coronavirus RaTG13, Pangolin coronavirus isolate MP789, and SARS-CoV obtained as BLAST hits were aligned and a phylogenetic tree was constructed (Figure 1).

The MSA demonstrated the molecular similarities between the organisms. RaTG13 has been identified as a neighbour genome for SARS-CoV-2 and this justifies the hypothesis that the infection may be transmitted from bats. Meanwhile, the subsequent neighbours were Pangolin MP789 and SARS-CoV. This preliminary sequence alignment enabled the understanding of sequence similarities and evolutionary information, which is deeply fundamental to the process of drug discovery.

159e66cc-ae9a-4e51-a245-a815f777c650_figure1.gif

Figure 1. Cladogram constructed using neighbour joining tree construction method, depicting the distance between the complete genome sequences subjected for multiple sequence alignment.

This figure was generated using EMBL-EBI Clustal Omega.

Topological analysis and target validation

A detailed investigation of the pathogen-host interactome can shed clear insights on the mechanism of viral infection and the pathology involved. Due to a lack of interaction data on SARS-CoV-2, the SARS-CoV proteome was considered and the SARS-CoV/human interactome was built by screening domain interactions between SARS-CoV/human protein-protein interactions, and then the network distribution, topological and functional analyses were performed (Figure 2). The circular shapes correspond to proteins (nodes) which are labelled by Uniprot_IDs and details about the nodes are listed in Table 1.

Among 14 proteins of SARS-CoV, the majority of SARS-CoV/human interaction involves five non-structural proteins (NS3B, NS6, NS7A, NS7B and NS8A with four, three, eight, three and one human proteins, respectively), three open reading frame (ORF) polyproteins (ORF9B, A7J8L3 and A7J8L2 with four, one and one human proteins, respectively), two replicase proteins (R1A and R1AB with two and one human proteins, respectively), the Membrane protein (VME1 with human IKKB), Envelope membrane protein (VEMP with human B2CL1), Nucleoprotein (NCAP with four human proteins) and Spike glycoprotein (SPIKE with human ACE2). With these observations, we determine the high specificity of Membrane, Envelope and Spike glycoprotein interactions with the host through specific entry points. Hence, these three SARS-CoV proteins can be a potential target to inhibit the pathogen-host interaction specifically, while other interactions are more versatile. To ensure the impact of inhibition of IKKB, B2CL1, and ACE2 mediated interaction, the landscape of the SARS-CoV/human interaction was further analysed for degree and betweenness centrality distributions of the host, as shown in Table 2.

159e66cc-ae9a-4e51-a245-a815f777c650_figure2.gif

Figure 2. The SARS-CoV/human pathogen-host interactome illustrating the protein-protein interactions between proteins.

The 28 blue nodes correspond to human proteins and 14 red nodes are SARS-CoV-2 proteins. The nodes are connected by 35 edges. This figure was generated using Pathogen-Host Interaction Search Tool (PHISTO).

Table 1. The list of proteins involved in SARS-CoV/human interactions with experimental methods and Pubmed ID as reference.

Uniprot IDPathogen proteinUniprot IDHuman proteinExperimental methodPubmed ID
P59596VME1O14920IKKBCoimmunoprecipitation17705188
P59637VEMPQ07817B2CL1Coimmunoprecipitation16048439
P59637VEMPQ07817B2CL1Other Methods16048439
P59594SPIKEQ9BYF1ACE2Experimental interaction detection18343844
P59594SPIKEQ9BYF1ACE2Pull down23486063
P0C6U8R1AP62753RS6Anti tag coimmunoprecipitation19838190
P0C6U8R1AA2A3R6A2A3R6Anti tag coimmunoprecipitation19838190
P0C6X7R1ABO75348VATG1Anti bait coimmunoprecipitation16226257
P0C6X7R1ABO75348VATG1Anti tag coimmunoprecipitation16226257
P0C6X7R1ABO75348VATG1Experimental interaction detection16226257
P0C6X7R1ABO75348VATG1Phage display16226257
P0C6X7R1ABO75348VATG1ELISA16226257
P59636ORF9BP08708RS17Two hybrid22046132
P59636ORF9BP60866RS20Two hybrid22046132
P59636ORF9BP27448MARK3Two hybrid22046132
P59636ORF9BP25787PSA2Two hybrid22046132
P59636ORF9BP25787PSA2LUMIER22046132
P59636ORF9BP27448MARK3LUMIER22046132
P59636ORF9BP60866RS20LUMIER22046132
P59636ORF9BP08708RS17LUMIER22046132
Q7TFA0NS8AP46379BAG6Two hybrid22046132
Q7TFA0NS8AP46379BAG6LUMIER22046132
Q7TFA1NS7BP49069CAMLGTwo hybrid22046132
Q7TFA1NS7BP05155IC1Two hybrid22046132
Q7TFA1NS7BE9KL26E9KL26Two hybrid22046132
Q7TFA1NS7BP05155IC1LUMIER22046132
Q7TFA1NS7BP49069CAMLGLUMIER22046132
Q7TFA1NS7BE9KL26E9KL26LUMIER22046132
P59635NS7AQ9BQB6VKOR1Two hybrid22046132
P59635NS7AQ9H4F8SMOC1Two hybrid22046132
P59635NS7AQ92843B2CL2Coimmunoprecipitation17428862
P59635NS7AQ07820MCL1Coimmunoprecipitation17428862
P59635NS7AP10415BCL2Coimmunoprecipitation17428862
P59635NS7AQ16548B2LA1Coimmunoprecipitation17428862
P59635NS7AO43765SGTACoimmunoprecipitation16580632
P59635NS7AQ07817B2CL1Coimmunoprecipitation17428862
P59635NS7AQ9BQB6VKOR1LUMIER22046132
P59635NS7AQ9H4F8SMOC1LUMIER22046132
P59635NS7AO43765SGTAOther Methods16580632
P59634NS6Q13561DCTN2Two hybrid22046132
P59634NS6Q92994TF3BTwo hybrid22046132
P59634NS6P52292IMA1Two hybrid17596301
P59634NS6Q13561DCTN2LUMIER22046132
P59634NS6Q92994TF3BLUMIER22046132
P59633NS3BQ13561DCTN2Two hybrid22046132
P59633NS3BP02768ALBUTwo hybrid22046132
P59633NS3BP05155IC1Two hybrid22046132
P59633NS3BE9KL26E9KL26Two hybrid22046132
P59633NS3BP05155IC1LUMIER22046132
P59633NS3BQ13561DCTN2LUMIER22046132
P59633NS3BP02768ALBULUMIER22046132
P59633NS3BE9KL26E9KL26LUMIER22046132
P59595NCAPP63165SUMO1Anti tag coimmunoprecipitation15848177
P59595NCAPP09651ROA1Two hybrid15862300
P59595NCAPP63279UBC9Two hybrid16998888
P59595NCAPP63279UBC9Coimmunoprecipitation16998888
P59595NCAPP63165SUMO1Coimmunoprecipitation15848177
P59595NCAPP62937PPIACoimmunoprecipitation15688292
P59595NCAPP09651ROA1Pull down15862300
P59595NCAPP09651ROA1Surface plasmon resonance15862300
P59595NCAPP62937PPIASurface plasmon resonance15688292
P59595NCAPP63279UBC9Other Methods16998888
P59595NCAPP63279UBC9Other Methods17037517
P59595NCAPP09651ROA1Other Methods15862300
A7J8L3A7J8L3A2A3R6A2A3R6Anti tag coimmunoprecipitation19838190
A7J8L2A7J8L2O75348VATG1Anti bait coimmunoprecipitation16226257
A7J8L2A7J8L2O75348VATG1Anti tag coimmunoprecipitation16226257
A7J8L2A7J8L2O75348VATG1Experimental interaction detection16226257
A7J8L2A7J8L2O75348VATG1Phage display16226257
A7J8L2A7J8L2O75348VATG1ELISA16226257

LUMIER, luminescence-based mammalian interactome mapping; ELISA, enzyme-linked immunosorbent assay.

Table 2. Topological analysis of the human proteins and of the human proteins targeted by SARS-CoV in the human interactome.

Cumulative degree and betweenness centrality distributions.

Sl.NoHuman
Uniprot ID		Human
protein		DegreeBetweenness
centrality
1A2A3R6A2A3R610.0
2Q9BYF1ACE242271.0
3Q9H4F8SMOC14115.53
4Q9BQB6VKOR184431.9
5P49069CAMLG113788.8
6P05155IC1138518.7
7O75348VATG11517898.0
8Q16548B2LA11510821.0
9Q92994TF3B197621.0
10Q92843B2CL22012060.0
11P08708RS172415475.0
12Q07820MCL13523003.0
13Q13561DCTN24432832.0
14P27448MARK34932544.0
15P62937PPIA5790417.0
16P25787PSA26826732.0
17Q07817B2CL184113360.0
18P10415BCL297171060.0
19O43765SGTA133545880.0
20O14920IKKB13499858.0
21P52292IMA1143310540.0
22P63165SUMO1164334120.0
23P60866RS2017037806.0
24P09651ROA1174176900.0
25P02768ALBU178735080.0
26P46379BAG6205684320.0
27P62753RS623867566.0
28P63279UBC9252817850.0

The degree of connectivity estimates the number of directly connecting neighbours to a particular node, while betweenness centrality estimates the frequency of nodes occurring on the shortest paths in the context of other nodes. In the protein interactomes, a node with a high degree of connectivity is identified as hub protein and a node with maximum betweenness centrality is identified as bottleneck protein. In the current topological analysis, the node with the lowest degree of distribution 1 and betweenness centrality 0.0 was A2A3R6. However, the molecular function of A2A3R6 (Uniprot ID: A2A3R6) is not well understood in both human physiology or pathology. Hence, ACE2, with the degree of distribution 4 and betweenness centrality 2271, was the next most significant node, as shown in Figure 3, and it was identified as a key node or key player in the SARS-CoV/host interaction. Hence, the SARS-CoV-S interaction with host ACE2 was identified as a potential drug target. As information about the SARS-CoV-2/human interaction is not available, the SARS-CoV/human interaction data was used. We studied the similarity between SARS-CoV and SARS-CoV-2 by sequence analysis and RBD prediction.

159e66cc-ae9a-4e51-a245-a815f777c650_figure3.gif

Figure 3. Cumulative degree and betweenness centrality distributions among host proteins involved in the SARS-CoV interaction.

Host proteins that are targeted by SARS-CoV have an approximate degree and betweenness centrality. A) Degree distribution highlighting ACE2 position and B) betweenness centrality distribution highlighting ACE2 position. These findings are statistically significant by the Fligner-Killeen (median) test. This figure was generated using Pathogen-Host Interaction Search Tool (PHISTO).

Sequence analysis and RBD predictions

As depicted in Figure 4, the alignment between the S-protein of SARS-CoV-2 and that of Bat coronavirus RaTG13 was closer than with the S-protein of SARS-CoV. The alignment at RBD site residues 317 to 569 was found to be more than 80% similar to SARS-CoV and RaTG13, particularly at major residues including Tyr436, Thr486, Gly488 and Tyr491 but excluding Arg426 and Pro462, as shown in Figure 5. Considering the evolution, the available elucidated structure of the SARS-CoV/ACE2 complex (PDB ID: 6CS2) was used as a template for homology modelling.

159e66cc-ae9a-4e51-a245-a815f777c650_figure4.gif

Figure 4. Multiple sequence alignment of S-protein from SARS_CoV-2, SARS-CoV and RatG13 to understand the conservation of amino acids as generated using EMBL-EBI Clustal Omega.

159e66cc-ae9a-4e51-a245-a815f777c650_figure5.gif

Figure 5. Aligned region of S-protein from SARS_CoV-2, SARS-CoV and RatG13 to understand the conservation of amino acids at the receptor-binding domain region as generated using EMBL-EBI Clustal Omega.

The residues involved in the interaction of SARS-CoV with ACE2 were predicted using the Prime module available in Schrodinger Small Molecule Suite and the major interactions are tabulated in Table 3 and shown in Figure 6. A very strong interaction was seen between the smallest amino acid, Gly488, with Lys353, Gly354, and Asp355. This interaction is facilitated by two features; one hydrogen bond and 94.9% buried solvent accessible surface area. Remaining residues also showed significantly strong interactions with ACE2. Hence, the same residues were made centric to generate the grid.

Table 3. List of residues involved in SARS-CoV-ACE2 interaction.

The residues with hydrogen bonding, better surface complementarity and buried SASA are tabulated and a few other weak interactions are ignored.

ResidueClosestDistanceInteractions# HB# SB# vWCSCB_SASA
B:426:ArgD:329:Glu2.7 A2x hb, 1x salt bridge, 1x clash to D:329:Glu2110.4372.0%
B:436:TyrD:38:Asp
D:42:Gln		3.2 A
3.2 A		1x hb to D:38:Asp1000.7536.9%
B:462:ProD:19:Ser
D:24:Gln		2.8 A
3.7 A		1x hb to D:19:Ser1000.7387.3%
B:486:ThrD:41:Tyr
D:330:Asn
D:355:Asp
D:357:Arg		2.8 A
2.9 A
3.1 A
3.8 A		1x hb to D:41:Tyr1000.8779.8%
B:488:GlyD:353:Lys
D:354:Gly
D:355:Asp		3.0 A
3.9 A
4.0 A		1x hb to D:353:Lys1000.6594.9%
B:491:TyrD:353:Lys
D:354:Gly
D:393:Arg		3.2 A
3.6 A
3.7 A		0000.8570.0%

HB, hydrogen bond; DS, disulfides, vWC, van der Wsals clash; SC, surface complementarity; B_SASA, buried solvent accessible surface area (SASA).

159e66cc-ae9a-4e51-a245-a815f777c650_figure6.gif

Figure 6. Amino acids involved in the SARS-CoV B-chain and ACE2 D-chain interaction.

The residues highlighted in blue correspond to the B-Chain of SARS-CoV and residues in green correspond to the ACE2 D-Chain. This figure was generated using Schrodinger’s Maestro visualizer. As an alternative, the open source visualizer Python Molecule Viewer (PMV)48 could also be used.

Homology modelling and validation of SARS-CoV-2-S protein

The modelling of SARS-CoV-2 was performed using the crystal structure of SARS-CoV-S as a template, which was 97% identical to the query. The modelled protein shown in Figure 7A was validated for quality and preparedness. The Ramachandran plot generated using the protein preparation wizard confirmed the quality of modelled structure by plotting >95% residues in the allowed region, as shown in Figure 7B.

159e66cc-ae9a-4e51-a245-a815f777c650_figure7.gif

Figure 7.

A) The modelled structure of SARS-CoV-2 with Arg426, Tyr436, Pro462, Thr486, Gly488 and Tyr491 centric grid box generated for molecular docking, B) Ramachandran plot has >95% of amino acids plotted in the allowed region (red and yellow). This figure was generated using Schrodinger’s Maestro visualizer. As an alternative, the open source visualizer Python Molecule Viewer (PMV)48 could also be used.

Molecular docking of small molecules with SARS-CoV-2-S

The repurposing of small molecules as therapeutics to treat COVID-19 requires knowledge of the interaction of the therapeutic molecule with SARS-CoV-2-S. Initial high-throughput virtual screening suggested 142 molecules that exhibit reasonable interaction with SARS-CoV-2-S. Following this, the shortlisted molecules were docked in SP mode where the accuracy of prediction was improved. The docking in SP mode suggested 15 top molecules, listed in Table 4 as lead molecules. As hydroxychloroquine has been identified as a possible treatment for COVID-19, it was also subjected to subsequent docking in XP mode. All 15 molecules showed better interaction than hydroxychloroquine with SARS-CoV-2-S. The three molecules streptomycin, ciprofloxacin, and GA had low interaction penalties and displayed better interactions with the ACE2 binding site on the RBD of SARS-CoV-2-S, as shown in Figure 8A–C, respectively. The three molecules were selected based on their reported anti-viral activity, safety, availability, and affordability45–47.

Table 4. List of top 15 lead molecules and hydroxychloroquine subjected for molecular docking in XP mode.

Sl. NoMolecule nameXP_ScoreMol_StructurePharmacological activity
1Fytic acid-10.296159e66cc-ae9a-4e51-a245-a815f777c650_t1.gifHypocalcemic agent
2Diquafosol-9.109159e66cc-ae9a-4e51-a245-a815f777c650_t2.gifPurinoceptor P2Y(2) receptor agonist
3Oxiglutatione-7.869159e66cc-ae9a-4e51-a245-a815f777c650_t3.gifIntraocular Irrigating Solution
4Ceftolozane-8.07159e66cc-ae9a-4e51-a245-a815f777c650_t4.gifAntibiotic
5Glycyrrhizic acid-7.474159e66cc-ae9a-4e51-a245-a815f777c650_t5.gifAnti-inflammatory
6Leucovorin-7.39159e66cc-ae9a-4e51-a245-a815f777c650_t6.gifFolic acid antagonist
7Cefpimizole-6.51159e66cc-ae9a-4e51-a245-a815f777c650_t7.gifAntibiotic
8Streptomycin-6.509159e66cc-ae9a-4e51-a245-a815f777c650_t8.gifAntibiotic
9Ertapenem-5.758159e66cc-ae9a-4e51-a245-a815f777c650_t9.gifAntibiotic
10Enoxacin-5.539159e66cc-ae9a-4e51-a245-a815f777c650_t10.gifAntibacterial
11Acetyldigitoxin-5.407159e66cc-ae9a-4e51-a245-a815f777c650_t11.gifCardioactive
12Pefloxacin-5.399159e66cc-ae9a-4e51-a245-a815f777c650_t12.gifAntibacterial
13Faropenem-5.354159e66cc-ae9a-4e51-a245-a815f777c650_t13.gifAntibiotic
14Ciprofloxacin-5.318159e66cc-ae9a-4e51-a245-a815f777c650_t14.gifAntibacterial
15Ceforanide-5.302159e66cc-ae9a-4e51-a245-a815f777c650_t15.gifAntibiotic
16Hydroxychloroquine-5.302159e66cc-ae9a-4e51-a245-a815f777c650_t16.gifAntimalerial

For SARS-CoV-2-S, the Glide generated docking model showed that streptomycin could bind to SARS-CoV-2-S in a manner highly similar to the SARS-CoV-2-S and ACE2 interaction. The binding pocket of streptomycin was in the RBD site, which has been observed to be an acceptor for ACE2. Streptomycin was well-fitted with the shape of the pocket, as shown in Figure 8A and Figure 9A, with an XP score of -6.5, where it formed a total five hydrogen bonds, among which two hydrogen bonds were formed by donating electrons from N31 and N32 atoms to the Glu493 side-chain atoms. Simultaneously, two other hydrogen bonds were observed between the backbone atoms of Leu501 by receiving electrons from hydroxyl groups at the 5th and 6th carbon atoms of the S1 six-carbon ring of streptomycin. The remaining H-bonds were formed between the backbone atom of Ser503 and the hydroxyl group at the 6th carbon atom at the G3 group of streptomycin. However, the stability of the interaction cannot be pronounced without molecular dynamic simulations.

159e66cc-ae9a-4e51-a245-a815f777c650_figure8.gif

Figure 8. Two-dimensional (2D) illustration of the small molecule interaction in the ACE2 binding site of the receptor-binding domain of SARS-CoV-2-S.

A) Interaction of streptomycin facilitated by five H-bonds, shown as pink arrows. B) Interaction of ciprofloxicine facilitated by two H-bonds, shown as pink arrows. C) Interaction of glycyrrhizic acid, facilitated by give H-bonds, shown as pink arrows. This figure was generated using Schrodinger’s Maestro visualizer. As an alternative, the open source visualizer Python Molecule Viewer (PMV)48 could also be used.

159e66cc-ae9a-4e51-a245-a815f777c650_figure9.gif

Figure 9. Three-dimensional (3D) illustration of the small molecule interaction in the ACE2 binding site of the receptor-binding domain of SARS-CoV-2-S.

A) Interaction of streptomycin facilitated by five H-bonds. B) Interaction of ciprofloxicine facilitated by two H-bonds. C) Interaction of glycyrrhizic acid facilitated by five H-bonds. This figure was generated using Schrodinger’s Maestro visualizer. As an alternative, the open source visualizer Python Molecule Viewer (PMV)48 could also be used.

The docking model of ciprofloxacin illustrated its binding mode on the RBD site, which has been observed to be a key interference site for virus-host interaction. The ciprofloxacin fit with reasonable steric complementarity into the RBD pocket, as shown in Figure 8B and Figure 9B, with an XP score of -5.31. The interaction of ciprofloxacin with SARS-CoV-2-S was facilitated by two hydrogen bonds between Val492 and Phe499; each bond being formed by receiving and donating electrons from hydroxyl and ketone groups, respectively.

The docking model of GA illustrated its binding mode on the RBD site, which has been observed to be a key site for interference of the virus-host interaction. The GA fit with steric complementarity in the RBD pocket, as shown in Figure 8C and Figure 9C, with an XP score of -7.474. The docking of GA with SARS-CoV-2-S was facilitated by three hydrogen bonds with Leu464, Val492, and Glu493 by receiving electrons from the hydroxyl groups of GA. Additionally, the ketone group of GA formed a hydrogen bond, with the backbone atoms of Phe499 receiving the electrons.

Molecular dynamics simulation of protein-ligand complex

As the SARS-CoV-2-S receptor has 1273aa, it requires enormous computational time to perform MD simulation for the whole range of protein, hence, we confined this study only to the RBD portion, ranging from 317th residue to 569th residue, for MD Simulation.

RMSD analysis of protein-ligand complex. The RMSD can illustrate the average difference in the displacement of selected atoms in a particular frame compared to its reference frame. The plots in Figure 10 illustrate the evolution of a protein (left Y-axis) and ligand (right Y-axis) RMSD. Post simulation, the protein and ligand frames are initially aligned over the backbone atom coordinates of the reference frame, and then the RMSD is extrapolated. The information on protein-ligand RMSD can dissect and demonstrate the conformational differences that occurred throughout the simulation. The RMSD of between 1–3 Å is fairly acceptable for small, globular proteins. An RMSD exceeding this indicates a major conformational change during the simulation and pronounces the instability of the complex.

159e66cc-ae9a-4e51-a245-a815f777c650_figure10.gif

Figure 10. The root mean square deviation (RMSD) plot for protein ligand complexes generated after 100 ns molecular dynamics simulation in the NPT ensemble.

A) streptomycin/SARS-CoV-2-S complex; B) ciprofloxacin SARS-CoV-2-S complex; C) glycyrrhizic acid/SARS-CoV-2-S complex. This figure was generated using Desmond version 4.242. This analysis could also have been performed using open source software such as GROMACS43.

The RMSD plot for the streptomycin/SARS-CoV-2-S complex, shown in Figure 10A, attained equilibrium at 5ns and thereafter showed stability with a maximum RMSD of 1 Å (peaks of 2.5 Å - 3.0 Å) up to 55ns. After 55ns a change in the equilibrium state was observed. However, the RMSD was within 1.5 Å, which is acceptable. Similarly, the streptomycin RMSD (right Y-axis) was observed to be significantly higher than the RMSD of the receptor at the RBD site. Thus, it is likely that streptomycin diffuses from its initial binding site after 48ns.

The RMSD plot for the ciprofloxacin/SARS-CoV-2-S complex, shown in Figure 10B, attained equilibrium at 2ns and thereafter showed stability with a maximum RMSD of 1.8 Å (peaks of 2.4 Å - 4.2 Å) up to 58ns. After 58ns a sudden change in equilibrium state was observed. However, the RMSD was within 2 Å, which is acceptable. On the other hand, the RMSD values for ciprofloxacin were observed to be significantly in alignment with the RMSD of SARS-CoV-2-S at the RBD site. Thus, it is likely that ciprofloxacin can retain its initial binding site up to 100ns.

The RMSD plot for the GA/SARS-CoV-2-S complex, shown in Figure 10C, attained equilibrium until 100ns. Compared to SARS-CoV-2-S complexes with streptomycin and ciprofloxacin, it was found that the SARS-CoV-2-S complex with GA was stable until the end of the simulation without any drift in equilibrium. On the other hand, the RMSD values for GA were observed to be significantly in alignment with the RMSD of the SARS-CoV-2-S RBD domain in almost all the frames. Hence, it is likely that it remains in its initial binding site up to 100ns. It is predicted to inhibit SARS-CoV-2-S at the RBD domain comparatively better than streptomycin and ciprofloxacin and for a longer duration, but its contact with key ligands has to be confirmed through RMSF and protein-ligand contact analysis.

RMSF analysis. The RMSF helps characterise minute differences in the protein chain during the simulation. In RMSF plots, peaks correspond to the residues on the protein that fluctuate more during the course of a simulation. Usually, terminals and loop regions fluctuate more than other secondary structures like alpha-helices and beta-strands. The secondary structure of the RBD of SARS-CoV-2-S has the same secondary structural elements as the RBD from SARS-CoV, with 74% homologous residues. These residues are majorly formed of loops and are highly flexible. A unique Phe486 residue in the loop plays a key role in ACE2 interaction by occupying a deep hydrophobic pocket in ACE2. In the trimmed RBD structure this loop starts from 148th residue and ends at 172nd residue. Since the ligand-binding site is located in this loop region, a higher RMSD was noticed. In the RMSF plot for the RBD domain of the streptomycin/SARS-CoV-2-S complex, shown in Figure 11A, the RMSF at the loop region was 5.6Å with many ligand contacts (green-coloured vertical bars). This was on par with molecular docking interactions. In the RMSF plot for the RBD domain of the ciprofloxacin/SARS-CoV-2-S complex, shown in Figure 11B, the RMSF at the loop region was 5.6Å with a few ligand contacts (green-coloured vertical bars). This justifies the interactions seen in molecular docking. Further, in the RMSF plot for the RBD domain of the GA/SARS-CoV-2-S complex, shown in Figure 11C, the RMSF at loop region was 6.3Å with a high number of ligand contacts (green-coloured vertical bars), justifying the interactions seen in molecular docking. Though the ligand contacts are seen in interactions, the simulation time coverage determines their stability.

Protein-ligand contacts. Protein-ligand interactions can be traced throughout the simulation and can be categorised into four types: hydrogen bonds, hydrophobic, ionic, and water bridges, as summarised in Figure 12A–12C. The stacked bars in the plots are normalised over the course of the trajectory and help us to understand the retention of contact throughout the simulation time. The contacts with a value of more than 0.7 are expected to be retained for over 70% of the total simulation time.

159e66cc-ae9a-4e51-a245-a815f777c650_figure11.gif

Figure 11. The root mean square fluctuation (RMSF) plot for protein ligand complexes generated after 100 ns molecular dynamics simulation in the NPT ensemble.

A) streptomycin/SARS-CoV-2-S complex; B) ciprofloxacin SARS-CoV-2-S complex; C) glycyrrhizic acid/SARS-CoV-2-S complex. Desmond version 4.242. This analysis could also have been performed using open source software such as GROMACS43.

In the protein-ligand contact plot for the streptomycin/SARS-CoV-2-S complex, shown in Figure 12A, residues Glu493 and Lys544 showed maximum interactions fractions, i.e. 0.20 facilitated by hydrogen bonds and water bridges. This suggests that the specific interaction is maintained for 20% of the simulation time, and such short interactions are not promising. Hence, streptomycin cannot be a potential inhibitor of SARS-CoV-2-S to offer anti-COVID-19 activity.

159e66cc-ae9a-4e51-a245-a815f777c650_figure12.gif

Figure 12. Histogram of protein ligand contacts generated for residues that occur for more than 64.0% of the simulation time in the selected trajectory.

A) streptomycin/SARS-CoV-2-S complex; B) ciprofloxacin SARS-CoV-2-S complex; C) glycyrrhizic acid/SARS-CoV-2-S complex. Desmond version 4.242. This analysis could also have been performed using open source software such as GROMACS43.

In the protein-ligand contact plot for the ciprofloxacin/SARS-CoV-2-S complex shown in Figure 12B, residues Phe465, Tyr482, Tyr498, and Phe499 were seen to have the interactions fractions 0.75, 0.6, 0.35 and 0.39 respectively facilitated by hydrophobic, hydrogen bonds and water bridges. This suggests that for 70%, 60%, 35% and 39% of the simulation time, the specific interaction is maintained by respective residues and such interactions are considered good. Hence, ciprofloxacin may be a potential inhibitor of SARS-CoV-2-S and may offer anti-COVID19 activity.

In the protein-ligand contact plot for the GA/SARS-CoV-2-S complex shown in Figure 12C, residues Val492, Glu493, Asn496, Cys497, and Phe499 had the interactions fractions 0.78, 1.12, 0.80, 0.60 and 0.80, respectively, facilitated by hydrophobic, hydrogen bonds and water bridges. This suggests that for 78%,100%, 80%, 60% and 80% of the simulation time, the specific interaction is maintained by respective residues and such interactions are excellent and promising. Hence, GA can be a potential inhibitor of SARS-CoV-2-S and can offer anti-COVID19 activity.

Conclusion

Through our topological analysis, we have determined the degree of distribution for viral proteins, and we show that, due to its low degree of distribution, ACE2 is likely to be targeted by viruses like SARS-CoV. Hence, the interaction between the viral protein SARS-CoV-2-S and the host ACE2 receptor is a potential drug target for the repurposing of known drugs. Further, sequence alignment and domain analysis suggest that the RBD is the ligand-binding site. Molecular docking studies have suggested streptomycin, ciprofloxacin, and GA as possible leads to inhibit SARS-CoV-2-S. Molecular dynamic simulation analysis has indicated that GA is a promising small molecule that could be repurposed as a potential inhibitor of SARS-CoV-2-S to offer anti-COVID19 activity.

Data availability

All data underlying the results are available as part of the article and no additional source data are required.

References

  • 1.  Zhou P, Yang XL, Wang XG, et al.: A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020; 579(7798): 270–273. PubMed Abstract | Publisher Full Text | Free Full Text
  • 2.  Gorbalenya AE, Baker SC, Baric RS, et al.: Severe acute respiratory syndrome-related coronavirus: The species and its viruses - a statement of the Coronavirus Study Group. Nat Microbiol. 2020. Publisher Full Text
  • 3.  World Health Organization: Clinical management of severe acute respiratory infection when novel coronavirus (CoV) infection is suspected. Reference Source
  • 4.  CDC. Reference Source
  • 5.  FDA. Reference Source
  • 6.  Korea Biomedical Review. Reference Source
  • 7.  Cheng F: In silico oncology drug repositioning and polypharmacology. Methods Mol Biol. 2019; 1878: 243–261. PubMed Abstract | Publisher Full Text
  • 8.  Cheng F, Hong H, Yang S, et al.: Individualized network-based drug repositioning infrastructure for precision oncology in the panomics era. Brief Bioinform. 2017; 18(4): 682–697. PubMed Abstract | Publisher Full Text
  • 9.  Nature Biotechnology. visited on 20th March. Reference Source
  • 10.  Walls AC, Park YJ, Tortorici MA, et al.: Structure, Function, and Antigenicity of the SARS-CoV-2-Spike Glycoprotein. Cell. 2020; 281–292.e6. PubMed Abstract | Publisher Full Text | Free Full Text
  • 11.  Belouzard S, Chu VC, Whittaker GR: Activation of the SARS coronavirus spike protein via sequential proteolytic cleavage at two distinct sites. Proc Natl Acad Sci U S A. 2009; 106(14): 5871–5876. PubMed Abstract | Publisher Full Text | Free Full Text
  • 12.  Bosch BJ, van der Zee R, de Haan CA, et al.: The coronavirus spike protein is a class I virus fusion protein: structural and functional characterization of the fusion core complex. J Virol. 2003; 77(16): 8801–8811. PubMed Abstract | Publisher Full Text | Free Full Text
  • 13.  Burkard C, Verheije MH, Wicht O, et al.: Coronavirus cell entry occurs through the endo-/lysosomal pathway in a proteolysis-dependent manner. PLoS Pathog. 2014; 10(11): e1004502. PubMed Abstract | Publisher Full Text | Free Full Text
  • 14.  Kirchdoerfer RN, Cottrell CA, Wang N, et al.: Pre-fusion structure of a human coronavirus spike protein. Nature. 2016; 531(7592): 118–121. PubMed Abstract | Publisher Full Text | Free Full Text
  • 15.  Millet JK, Whittaker GR: Host cell entry of Middle East Respiratory Syndrome coronavirus after two-step, furin-mediated activation of the spike protein. Proc Natl Acad Sci U S A. 2014; 111(42): 15214–15219. PubMed Abstract | Publisher Full Text | Free Full Text
  • 16.  Millet JK, Whittaker GR: Host cell proteases: Critical determinants of coronavirus tropism and pathogenesis. Virus Res. 2015; 202: 120–134. PubMed Abstract | Publisher Full Text | Free Full Text
  • 17.  Park JE, Li K, Barlan A, et al.: Proteolytic processing of Middle East respiratory syndrome coronavirus spikes expands virus tropism. Proc Natl Acad Sci U S A. 2016; 113(43): 12262–12267. PubMed Abstract | Publisher Full Text | Free Full Text
  • 18.  Walls AC, Tortorici MA, Bosch BJ, et al.: Cryo-electron microscopy structure of a coronavirus spike glycoprotein trimer. Nature. 2016; 531(7592): 114–117. PubMed Abstract | Publisher Full Text | Free Full Text
  • 19.  Walls AC, Tortorici MA, Frenz B, et al.: Glycan shield and epitope masking of a coronavirus spike protein observed by cryo-electron microscopy. Nat Struct Mol Biol. 2016; 23(10): 899–905. PubMed Abstract | Publisher Full Text | Free Full Text
  • 20.  Gui M, Song W, Zhou H, et al.: Cryo-electron microscopy structures of the SARS-CoV spike glycoprotein reveal a prerequisite conformational state for receptor binding. Cell Res. 2017; 27(1): 119–129. PubMed Abstract | Publisher Full Text | Free Full Text
  • 21.  Kirchdoerfer RN, Cottrell CA, Wang N, et al.: Pre-fusion structure of a human coronavirus spike protein. Nature. 2016; 531(7592): 118–121. PubMed Abstract | Publisher Full Text | Free Full Text
  • 22.  Pallesen J, Wang N, Corbett KS, et al.: Immunogenicity and structures of a rationally designed prefusion MERS-CoV spike antigen. Proc Natl Acad Sci U S A. 2017; 114(35): E7348–E7357. PubMed Abstract | Publisher Full Text | Free Full Text
  • 23.  Walls AC, Tortorici MA, Snijder J, et al.: Tectonic conformational changes of a coronavirus spike glycoprotein promote membrane fusion. Proc Natl Acad Sci U S A. 2017; 114(42): 11157–11162. PubMed Abstract | Publisher Full Text | Free Full Text
  • 24.  Yuan Y, Cao D, Zhang Y, et al.: Cryo-EM structures of MERS-CoV and SARS-CoV spike glycoproteins reveal the dynamic receptor binding domains. Nat Commun. 2017; 8: 15092. PubMed Abstract | Publisher Full Text | Free Full Text
  • 25.  Madu IG, Roth SL, Belouzard S, et al.: Characterization of a highly conserved domain within the severe acute respiratory syndrome coronavirus spike protein S2 domain with characteristics of a viral fusion peptide. J Virol. 2009; 83(15): 7411–7421. PubMed Abstract | Publisher Full Text | Free Full Text
  • 26.  Heald-Sargent T, Gallagher T: Ready, Set, Fuse! The Coronavirus Spike Protein and Acquisition of Fusion Competence. Viruses. 2012; 4(4): 557–580. PubMed Abstract | Publisher Full Text | Free Full Text
  • 27.  Morris DE, Cleary DW, Clarke SC: Secondary Bacterial Infections Associated with Influenza Pandemics. Front Microbiol. 2017; 8: 1041. PubMed Abstract | Publisher Full Text | Free Full Text
  • 28.  Zhang Z, Schwartz S, Wagner L, et al.: A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000; 7(1–2): 203–14. PubMed Abstract | Publisher Full Text
  • 29.  Raman K: Construction and analysis of protein-protein interaction networks. Autom Exp. 2010; 2(1): 2–10. PubMed Abstract | Publisher Full Text | Free Full Text
  • 30.  Hwang S, Son SW, Kim SC, et al.: A protein interaction network associated with asthma. J Theor Biol. 2008; 252(4): 722–731. PubMed Abstract | Publisher Full Text
  • 31.  Sastry, GM, Adzhigirey M, Day T, et al.: Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments. J Comput Aid Mol Des. 2013; 27(3): 221–234. PubMed Abstract | Publisher Full Text
  • 32.  Morris GM, Huey R, Lindstrom W, et al.: Autodock4 and AutoDockTools4: automated docking with selective receptor flexiblity. J Comput Chem. 2009; 30(16): 2785–91. PubMed Abstract | Publisher Full Text | Free Full Text
  • 33.  Grosdidier A, Zoete V, Michielin O: SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 2011; 39(Web Server issue): W270–W277. PubMed Abstract | Publisher Full Text | Free Full Text
  • 34.  Jacobson MP, Pincus DL, Rapp CS, et al.: A Hierarchical Approach to All-Atom Protein Loop Prediction. Proteins. 2004; 55(2): 351–367. PubMed Abstract | Publisher Full Text
  • 35.  Shivakumar D, Williams J, Wu Y, et al.: Prediction of absolute solvation free energies using molecular dynamics free energy perturbation and the OPLS force field. J Chem Theory Comput. 2010; 6(5): 1509–1519. PubMed Abstract | Publisher Full Text
  • 36.  Dagan-Wiener A, et al.: Bitter or not? BitterPredict, a tool for predicting taste from chemical structure. Sci Rep. 2017; 7(1): 12074. PubMed Abstract | Publisher Full Text | Free Full Text
  • 37.  Daina A, Michielin O, Zoete V: SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci Rep. 2017; 7: 42717. PubMed Abstract | Publisher Full Text | Free Full Text
  • 38.  Friesner RA, Murphy RB, Repasky MP, et al.: Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. J Med Chem. 2006; 49(21): 6177–6196. PubMed Abstract | Publisher Full Text
  • 39.  Guo Z, Mohanty U, Noehre J, et al.: Probing the α-helical structural stability of stapled p53 peptides: molecular dynamics simulations and analysis. Chem Biol Drug Des. 2010; 75(4): 348–359. PubMed Abstract | Publisher Full Text
  • 40.  Martyna GJ, Klein ML, Tuckerman M: Nose-Hoover chains-the canonical ensemble via continuous dynamics. J Chem Phys. 1992; 97(4): 2635–2643. Publisher Full Text
  • 41.  Jorgensen WL, Chandrasekhar J: Comparison of simple potential functions for simulating liquid water. J Chem Phys. 1983; 79(2): 926–935. Publisher Full Text
  • 42.  Bowers KJ, Chow E, Xu H, et al.: Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters. Proceedings of the ACM/IEEE Conference on Supercomputing (SC06), Tampa, Florida, 2006. Reference Source
  • 43.  Abraham MJ, Murtola T, Schulz R, et al.: GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX. 2015; 1-2: 19–25. Publisher Full Text
  • 44.  Essmann U, Perera L, Berkowit ML, et al.: A smooth particle mesh Ewald method. J Chem Phys. 1995; 103: 8577–8593. Publisher Full Text
  • 45.  Thomas DB, Susan OW: Streptomycin as an Antiviral Agent: Mode of Action. Science. 1963; 141(3585): 1065–1067. PubMed Abstract | Publisher Full Text
  • 46.  Ilaria U, Elena A, Francesca T, et al.: Possible Antiviral Effect of Ciprofloxacin Treatment on Polyomavirus BK Replication and Analysis of Non-Coding Control Region Sequences. Virol J. 2013; 10: 274. PubMed Abstract | Publisher Full Text | Free Full Text
  • 47.  Pompei P, Pani A, Flore O, et al.: Antiviral activity of glycyrrhizic acid. Experientia. 1980; 36(3): 304. PubMed Abstract | Publisher Full Text
  • 48.  Sanner MF: Python: A Programming Language for Software Integration and Development. J Mol Graph Model. 1999; 17(1): 57–61. PubMed Abstract

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 01 Jul 2020
Comment
Author details Author details
Competing interests
Grant information
Article Versions (1)
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
 - -
Citations
SEE MORE DETAILS
CITE
how to cite this article
BR B, Damle H, Ganju S and Damle L. In silico screening of known small molecules to bind ACE2 specific RBD on Spike glycoprotein of SARS-CoV-2 for repurposing against COVID-19 [version 1; peer review: 2 approved] F1000Research 2020, 9:663 (https://doi.org/10.12688/f1000research.24143.1)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 1
VERSION 1
PUBLISHED 01 Jul 2020
Views
12
Cite
Reviewer Report 27 Jul 2020
Raghavendra HL, College of Medical and Health Sciences, Wollega University, Nekemte, Ethiopia 
Approved
VIEWS 12
https://doi.org/10.5256/f1000research.26631.r67868
The manuscript entitled “In silico screening of known small molecules to bind ACE2 specific RBD on Spike glycoprotein of SARS-CoV-2 for repurposing against COVID-19” addresses the repurposability of known drugs against COVID-19. Spike glycoprotein is the most popular drug target for ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
HL R. Reviewer Report For: In silico screening of known small molecules to bind ACE2 specific RBD on Spike glycoprotein of SARS-CoV-2 for repurposing against COVID-19 [version 1; peer review: 2 approved]. F1000Research 2020, 9:663 (https://doi.org/10.5256/f1000research.26631.r67868)
The direct URL for this report is:
https://f1000research.com/articles/9-663/v1#referee-response-67868
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Report a concern
Views
15
Cite
Reviewer Report 14 Jul 2020
S.M. Vidya, Department of Biotechnology, NMAM Institute of Technology, Udupi, Karnataka, India 
Approved
VIEWS 15
https://doi.org/10.5256/f1000research.26631.r66947
In the present study authors have made an effort to repurpose the known drugs against SARS-CoV-2. Although the study is well organized and presented, a few minor changes are suggested for betterment. 
  1. The name of the molecules
... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Vidya SM. Reviewer Report For: In silico screening of known small molecules to bind ACE2 specific RBD on Spike glycoprotein of SARS-CoV-2 for repurposing against COVID-19 [version 1; peer review: 2 approved]. F1000Research 2020, 9:663 (https://doi.org/10.5256/f1000research.26631.r66947)
The direct URL for this report is:
https://f1000research.com/articles/9-663/v1#referee-response-66947
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Report a concern

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 01 Jul 2020
Comment

Open Peer Review

Reviewer Status

Alongside their report, reviewers assign a status to the article:

Approved
The paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations
A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved
Fundamental flaws in the paper seriously undermine the findings and conclusions

Reviewer Reports

Invited Reviewers
1 2
Version 1
01 Jul 20 		read read
  1. S.M. Vidya, NMAM Institute of Technology, Udupi, India
  2. Raghavendra HL, Wollega University, Nekemte, Ethiopia

Comments on this article

All Comments(0)

Add a comment
Sign up for content alerts

Browse by related subjects

    keyboard_arrow_left Back to all reports
    keyboard_arrow_left Back to all reports
    Alongside their report, reviewers assign a status to the article:
    Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
    Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
    Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
    Sign In
    If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

    The email address should be the one you originally registered with F1000.
    Email address not valid, please try again

    You registered with F1000 via Google, so we cannot reset your password.

    To sign in, please click here.

    If you still need help with your Google account password, please click here.

    You registered with F1000 via Facebook, so we cannot reset your password.

    To sign in, please click here.

    If you still need help with your Facebook account password, please click here.

    Code not correct, please try again
    Email us for further assistance.
    Server error, please try again.