A Novel Therapeutic Peptide Blocks SARS-CoV-2 Spike Protein Binding with Host Cell ACE2 Receptor.

BACKGROUND AND
OBJECTIVE: Coronavirus disease 2019 is a novel disease caused by the severe acute respiratory syndrome coronavirus (SARS-CoV)-2 virus. It was first detected in December 2019 and has since been declared a pandemic causing millions of deaths worldwide. Therefore, there is an urgent need to develop effective therapeutics against coronavirus disease 2019. A critical step in the crosstalk between the virus and the host cell is the binding of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein to the peptidase domain of the angiotensin-converting enzyme 2 (ACE2) receptor present on the surface of host cells.
METHODS: An in silico approach was employed to design a 13-amino acid peptide inhibitor (13AApi) against the RBD of the SARS-CoV-2 spike protein. Its binding specificity for RBD was confirmed by molecular docking using pyDockWEB, ClusPro 2.0, and HDOCK web servers. The stability of 13AApi and the SARS-CoV-2 spike protein complex was determined by molecular dynamics simulation using the GROMACS program while the physicochemical and ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties of 13AApi were determined using the ExPASy tool and pkCSM server. Finally, in vitro validation of the inhibitory activity of 13AApi against the spike protein was performed by an enzyme-linked immunosorbent assay.
RESULTS: In silico analyses indicated that the 13AApi could bind to the RBD of the SARS-CoV-2 spike protein at the ACE2 binding site with high affinity. In vitro experiments validated the in silico findings, showing that 13AApi could significantly block the RBD of the SARS-CoV-2 spike protein.
CONCLUSIONS: Blockage of binding of the SARS-CoV-2 spike protein with ACE2 in the presence of the 13AApi may prevent virus entry into host cells. Therefore, the 13AApi can be utilized as a promising therapeutic agent to combat coronavirus disease 2019.

Key Points

Introduction

The pandemic of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), with more than 180 million infections and close to four million deaths, can be designated as one of the greatest threats of the century to public health. It has not only posed a significant challenge to human health globally but also caused socio-economic crises for many countries [1, 2]. The clinical spectrum of coronavirus disease 2019 (COVID-19) ranges from mild fever, cough, and shortness of breath, to severe clinical conditions characterized by respiratory failure and other long-term associated conditions [3-7]. Old age, together with pre-existing conditions such as lung or heart disease, diabetes mellitus, or a compromised immune system are known to expedite the infection time and severity [8, 9].

Structurally, the coronavirus is a single-stranded, non-segmented, positive-sense RNA virus that consists of the largest known RNA genome of 26–32 kB compared to other known viruses [10, 11]. The four major structural proteins encoded by its genome include the spike (S) protein, envelope (E) protein, membrane (M) protein, and nucleocapsid (N) protein as well as several non-structural proteins including viral proteases [12, 13]. The surface-bound spike protein, suggested to be highly susceptible to mutations, plays a crucial role in facilitating virus entry by mediating its interaction with transmembrane surface receptors on host cells. The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein interacts with the peptidase domain (PD) of the angiotensin-converting enzyme 2 (ACE2) receptor on host cells, which technically marks virus entry inside the cells [14-18]. The S1 and S2 subunits in the spike protein of SARS-CoV-2 are mainly responsible for the interaction and fusion with the host cells for virus entry. The RBD in the S1 subunit initiates direct binding with the ACE2 PD, whereas the S2 subunit contains essential elements needed for membrane fusion [16, 19–22].

In humans, ACE2 is a type I, transmembrane endothelium-bound metallo-carboxypeptidase, with homology to the angiotensin-converting enzyme, well known for its role in the renin-angiotensin system [23-26]. ACE2  is a target for the treatment of hypertension [27, 28]. It is mainly expressed on endothelial cells of several organs, particularly the cardiovascular system, renal tubular epithelium, Leydig cells in the testes, and alveolar epithelial type II cells in the lungs and brain [23, 29–31]. The full-length structure of ACE2 consists of two main domains: the PD at the N-terminus and the collectrin-like domain at the C-terminus [32-34]. The spike glycoprotein of SARS-CoV-2 is known to bind to a homodimer of ACE2 [16, 17, 32].

In the current SARS-CoV-2 pandemic, a vaccine is rightfully anticipated as the most effective therapy. However, in an ongoing pandemic, vaccine development and production may not be equally efficacious in all countries. Therefore, potent, easy to generate, and cheap therapeutic agents that can effectively curb infection at the early stages also need to be developed. Several approaches, such as decoy-soluble ACE2 proteins, antibodies from the serum of infected patients, epitope-based vaccines, repurposing of drugs, and designing peptide inhibitors have been reported [35-44]. Peptides possess several attractive features compared with small-molecule and protein therapeutics, including high structural compatibility with target proteins, and the ability to specifically disrupt protein–protein interfaces [45, 46].

Effective use of computational tools can ease the way to rapidly reach a therapeutic solution for COVID-19. In the current study, we used computational biology tools to develop a therapeutic strategy utilizing a novel peptide that can inhibit the SARS-CoV-2/ACE2 interaction. We have designed a short and efficient peptide that imitates the binding of ACE2 with RBD of the SARS-CoV-2 spike protein and could significantly block their interaction in an in vitro assay system.

Materials and Methods

Structure Retrieval, Interaction Study, Peptide Design, Structure Preparation, and Molecular Docking Studies

The experimentally solved x-ray crystallographic structure of the SARS-CoV-2 spike protein and human ACE2 protein complex with PDB ID 6M17 was retrieved from the RCSB PDB database (https://www.rcsb.org/). The interaction between the SARS-CoV-2 spike protein and the ACE2 receptor protein was identified and analyzed using the Chimera Version 1.13.1 tool [47]. We focused on the interface residues of both proteins that are involved in the interaction within the 3Å region.

To target the RBD residues of the SARS-CoV-2 spike protein, a decoy peptide from the region of its interaction with ACE2 was designed. A 13-amino acid (13AA) stretch from F28 to L40 (FLDKFNHEAEDLF) from the PD of ACE2 was selected to design a short and effective peptide inhibitor. Modifications of the 13AA peptide sequence were carried out to improve its binding affinity without disturbing its physiochemical properties. Substitutions of residues at the eighth, ninth, and tenth positions gave us a 13AA peptide inhibitor (13AApi) with sequence FLDKFNHDLF, which showed higher binding affinity and specificity. The 13AApi sequence was modeled in a three-dimensional structure and energetically minimized on Chimera using the Amber ff14SB force field [47].

The SARS-CoV-2 spike RBD structure obtained from PDB 6M17 was prepared for docking experiments by removing water, heteroatoms, and ligand groups, followed by adding polar hydrogen atoms to the structure using the Discovery Studio Version 21.1.0 tool [48]. Molecular docking of the 13AApi with SARS-CoV-2 spike RDB was carried out in pyDockWEB [49] and cross-checked with the HDOCK [50] and ClusPro 2.0 [51] web server tools. The docking results were visualized and analyzed in Discovery Studio Visualizer Version 21.1.0 and Chimera Version 1.13.1 tools, respectively [47, 48].

Molecular Dynamics (MD)

The system for MD simulation was prepared by soaking in a dodecahedron box filled with TIP3P water molecules. The simulation system was neutralized by adding Na+ or Cl− ions. The protein, solvent, and ion parameters were assigned using the Amber99SB-ILDN force field [52]. The initial configuration of the system was relaxed through the steepest descent energy minimization using a maximum number of steps value of 50,000. After energy minimization, all simulation systems were equilibrated in two phases: first under NVT ensemble for 1 ns to stabilize the temperature of the system and then under NPT ensemble for another 1 ns to stabilize the pressure of the simulation system. Finally, production MD runs were carried out at a constant temperature of 300 K for 500 ns using GROMACS Version 5.0.4. [53]. A time step of 2 fs was used for the simulations. A cut-off of 10 Å was used for short-range interactions while the particle-mesh Ewald method was used to handle long-range interactions. Coordinates and velocities were saved every 10 ps. The MD trajectory was processed and analyzed using MDAnalysis [54], MDTraj [55], and scikit-learn [56] python libraries. Principal component analysis was performed using MODE-TASK [57]. The MM-PBSA method was employed to evaluate the 13AApi binding energies. MM-PBSA energy calculations were performed using the g_mmpbsa program [58, 59]. Graphics were prepared using Matplotlib python library [60] and R [61].

Physicochemical and ADMET (Adsorption, Distribution, Metabolism, Excretion, and Toxicity) Analysis of 13AApi

The physicochemical and ADMET properties of the 13AApi were analyzed from ExPasy ProtParam (https://web.expasy.org/protparam/) [62] and pkCSM (http://biosig.unimelb.edu.au/pkcsm/prediction) [63] tools available online. For ADMET analysis, the pkCSM tool requires a SMILES file format of the 13AApi, which was generated using the pepSMI (https://www.novoprolabs.com/tools/convert-peptide-to-smiles-string) web tool.

Enzyme-Linked Immunosorbent Assay (ELISA)

An in vitro assay was performed to validate the in silico interaction data with 13AApi. The percentage binding of the SARS-CoV-2 spike protein to ACE2 in the presence of the 13AApi was determined with an ACE2:SARS-CoV-2 spike inhibitor screening assay kit (BPS Bioscience #79936; San Diego, CA, USA) according to the protocol described in the manufacturer’s datasheet. A stock solution of the 13AApi was prepared in distilled water and the indicated working concentrations were used in the assay. A chemiluminescent reading from three biological replicates was used to determine the 13AApi activity graph.

Results and Discussion

Structural Analysis of SARS-CoV-2 Spike Protein Interaction with Host Receptor ACE2 and Design of the Novel Peptide Inhibitor

A crucial step in SARS-CoV-2 infection is the attachment of its surface spike protein with the ACE2 surface receptor on the host cell; this binding has now been widely studied and delineated [16, 17]. In terms of therapeutic strategies whereby antiviral drugs can block this interaction, previous work has mainly demonstrated the host ACE2 domain as a target; very few studies have aimed to target the SARS-CoV-2 spike protein for drug development [64-68]. Targeting the SARS-CoV-2 spike protein instead of the host ACE2 protein ensures direct action on the virus. Of note, no major therapeutic breakthroughs have been reported in targeting the SARS-CoV-2 spike protein to disrupt this primary step of interaction and disable fusion of the virus with the host cell. Several peptide-based therapeutic options to target SARS-CoV-2 and inhibit its interaction with ACE2 receptors have been reported by us and many other laboratories [69-73]. However, these peptide inhibitors are long-stretch peptides that are mainly derived from ACE2 receptors; moreover, very few have been validated in a cell-based assay.

We have taken this opportunity to further explore and develop a short novel peptide inhibitor that can specifically and significantly target the RBD residues of the SARS-CoV-2 spike protein and block its interaction with the host, the ACE2 protein. The N-terminal residues of the ACE2 PD are mainly involved in its interaction with the spike protein. Therefore, a similar stretch of peptide that can mimic this binding was designed with the aim to block the RBD of the SARS-CoV-2 spike protein. The available crystal structure complex of the SARS-CoV-2 spike protein with human ACE2 (PDB ID 6M17) was retrieved from PDB database and explored as the basis of the current study. We analyzed the interaction interface and identified crucial amino acids involved in protein–protein interactions (Fig. 1A, B).

Fig. 1

A Illustration of the interacting interface (encircled) of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) spike receptor-binding domain [RBD] (green) and human angiotensin-converting enzyme 2 [ACE2] (purple) from the crystal structure (PDB ID 6M17). B Identification of the interacting residues within the 3 Å region. The highlighted residues in the SARS-CoV-2 spike RBD (blue) interact with human ACE2 peptidase domain residues (bright green) as demonstrated with the Discovery studio Visualizer tool. The Table shows the interacting residues within a 3 Å region that was analyzed using the Chimera tool

We next prepared a set of peptide sequences similar to the N-terminal region of the ACE2 PD that is known to interact with SARS-CoV-2 within a 3–5 Å region and tested their binding affinity and specificity with the RBD of the spike protein. The residues within the 3 Å region are crucial for a complex formation and blocking such residues could interfere with the interaction of the virus with the host cell. A 13AA stretch from F28 to L40 (FLDKFNHEAEDLF) in the PD of ACE2 was found to effectively bind to the spike protein RBD and was suitable as a template to design a short and effective peptide inhibitor. We performed modifications of the 13AA peptide sequence to improve its binding affinity without disturbing the RBD binding site and physiochemical properties. Substitution of the amino acid residues at the eighth, ninth, and tenth positions (E8N, A9F, and E10K) gave us a 13AA peptide inhibitor with sequence FLDKFNHDLF, which showed higher binding affinity and specificity towards SARS-CoV-2 RBD. This peptide inhibitor was named 13AApi.

A protein-peptide blind docking study indicated that the 13AApi best-fit docking solution could bind to the residues of the spike protein RBD in the same pocket where the ACE2 PD binds (Fig. 2A). The binding energy of the initial peptide inhibitor 13AA (FLDKFNHEAEDLF) was obtained as − 42.177 kcal/mol in pyDockWEB while the docking score for the final peptide inhibitor 13AApi was increased to − 46.967 kcal/mol (Fig. 2A). We next cross-checked the docking on two other platforms: HDOCK and ClusPro 2.0. The dock score for 13AApi was − 187.01 kcal/mol in HDOCK and − 678.2 kcal/mol in ClusPro 2.0 while the scores of 13AA were − 168.91 kcal/mol in HDOCK and − 634.2 kcal/mol in ClusPro 2.0 (data not shown). Therefore, the results from HDOCK and ClusPro 2.0 are in good agreement with the pyDockWEB data, corroborating an enhancement in the docking score of 13AApi in comparison to the initial peptide. The binding position of the 13AApi was also in accord with the pyDockWEB result. The docking complex from pyDockWEB was employed for further analyses on interaction and representation.

Fig. 2

Interaction of the novel 13-amino acid peptide inhibitor (13AApi) with the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) spike receptor-binding domain (RBD) restricts its interaction with host cell receptor angiotensin-converting enzyme 2 (ACE2) peptidase domain. Molecular docking complex of the 13AApi (red) with SARS-CoV-2 spike RBD (green) obtained with pyDockWEB. The interacting residues within the 3 Å region are highlighted in cyan and blue, respectively, while the dock scores (kcal/mol) from all three docking suits (pyDockWEB, HDOCK, and ClusPro 2.0) are listed in the table below

Three residues K417, Y489, and Q493 of the spike protein RBD within the 3 Å region were found to interact with the novel 13AApi; these residues also interact with the ACE2 receptor (Fig. 2A). The docking position of 13AApi in the RBD pocket ensured a high probability to block the interaction with the ACE2 receptor. It is also imperative to understand the stability of a complex at an atomic level. Hence, we proceeded with an MD simulation study to analyze the dynamics and stability components.

MD Simulation Study of 13AApi and the SARS-CoV-2 Spike Protein

The docking complex between the SARS-CoV-2 spike protein and the 13AApi was next subjected to an MD simulation study to determine the stability of the complex. An MD simulation of 500 ns was performed using the GROMACS program. Analysis of the root mean squared deviation and the root mean squared fluctuation of the SARS-CoV-2 spike protein throughout the MD trajectory suggested no significant fluctuations in the ACE2 binding interface of the SARS-CoV-2 spike protein (Fig. 3A, B). Moreover, no significant changes in the radius of gyration and fraction of native contacts were observed, suggesting an overall stability of the SARS-CoV-2 spike protein and the 13AApi protein-peptide complex (Fig. 3C, D).

Fig. 3

MD simulation to study the dynamics of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) spike protein and 13-amino acid peptide inhibitor (13AApi) complex. A Root mean squared deviation (RMSD) of the Cα atoms of the SARS-CoV-2 spike protein (black), 13AApi (blue), and protein-peptide interface (red) plotted against total simulation time. B Root mean squared fluctuation of SARS-CoV-2 spike protein (black) and 13AApi. C The radius of gyration of the whole protein-peptide complex. D Fraction of native contacts of the whole protein-peptide complex. E Distance between Cα atom of K417, Y489, and Q493 with the center of mass of 13AApi at t = 0 ns is shown as dashed lines

The 13AApi initially showed significant fluctuations from the docking predicted binding mode, which was stabilized after 400 ns (Fig. 3B). Conformational and positional rearrangement accounted for most of the deviations and fluctuations as 13AApi remained bound to the SARS-CoV-2 spike protein for the entire 500 ns MD trajectory (Fig. 3E). To further study the dynamics of 13AApi within the SARS-CoV-2 spike protein receptor-binding interface, a principal component analysis using cartesian coordinates and singular value decomposition approach was performed. As shown in Fig. 4A, the first three principal components explained about 85% of the variance. Analysis of the first three principal components (Fig. 4C–E) and snapshots sampled at intervals of 100 ns (Fig. 4E) revealed several conformational states. The 13AApi did not maintain a stable helical conformation as in the case of an ACE2 fragment (FLDKFNHEAEDLF), possibly because of the disruption of intramolecular contacts. Instead, 13AApi exhibited frequent unfolding and folding events (Fig. 4). A figure showing the change in the secondary structure throughout the 500 ns MD simulation is given as supporting information (Fig. S1 of the Electronic Supplementary Material). At the end of the MD simulation, i.e., at around 500 ns, 13AApi appeared to adopt a helical conformation that was necessary for its interaction with the SARS-CoV-2 spike protein. However, it seems that  further optimization of 13AApi is needed to identify a sequence with stable helical conformation. To further estimate the energetics of 13AApi binding to the SARS-CoV-2 spike protein, MM-PBSA analysis was performed that revealed favorable binding energy (Table 1).

Fig. 4

A Principal component analysis (PCA) to study the dynamics of 13-amino acid peptide inhibitor (13AApi) within the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) spike protein receptor-binding interface. B–E Screen plots showing percent variance explained by a PCA

Table 1

13-Amino acid peptide inhibitor binding energies calculated using the MM-PBSA approach. The molecular dynamics trajectory of the last 100 ns was used for the MM-PBSA binding energy calculations

Energy type	Energy (kJ/mol)
Van der Waal (ΔE_vdw)	− 189.89 ± 31.24
Electrostatic (ΔE_elec)	− 397.57 ± 88.48
Polar solvation energy (ΔG_polar)	445.82 ± 90.61
Nonpolar energy (ΔG_nonpolar)	− 25.74 ± 2.70
Binding energy (ΔG_bind)	− 167.38 ± 45.08

In Silico Analysis of Physicochemical and ADMET Properties of 13AApi

To further understand the properties of the 13AApi in the host system, we performed an in silico prediction of its physicochemical and ADMET properties (Table 2). In brief, the 13AApi with an average molecular weight of 1684.91 g/mol and molecular formula C82H113N19O20 contains two each of the positively charged residues (K4 and K10) and negatively charged residues (D3 and D11) with a theoretical isoelectric point (pI) of 7.55. The measure of the instability index suggests that the 13AApi is in the category of a stable peptide with an estimated half-life of 1.1 h. Additionally, the hydropathy data indicated that it has suitable water-solubility properties.

Table 2

In-silico analysis of physicochemical and ADMET (absorption, distribution, metabolism, excretion and toxicity) properties of 13-amino acid peptide inhibitor (13AApi) from ExPasy ProtParam (ExPASy—ProtParam tool) and pkCSM (pkCSM [unimelb.edu.au]) tools

Peptide sequence	FLDKFNHNFKDLF
Physicochemical analysis
Properties
Molecular weight (g/mol)	1684.91
Theoretical pI	7.55
Net charge at pH 7	0.1
Negative + positive residues	2 +2
Molecular formula	C82H113N19O20
Number of atoms	234
Instability Index	− 2.21
Half-life (in h)	1.1
Hydrophobicity (%)	46.15
Acidic + Basic + neutral ratio (%)	15.38 + 23.08 + 15.38
GRAVY	− 0.477
ADMET analysis
Absorption
Water solubility (log mol/L)_	− 2.892
Skin perm (log Kp)	− 2.735
Distribution
Fraction unbound (human) (Fu)	0.369
BBB permeability	− 2.493
CNS permeability (logPS)	− 7.089
Metabolism
Cytochrome P substrate	No
Cyctochrome P inhibitor	No
Excretion
Total clearance (log mL/min/kg)	− 1.052
Toxicity
AMES toxicity	No
Skin sensitisation	No
MRTD human (log mg/kg/day)	0.438
Rat oral LD50 (mol/kg)	2.482

For the ADMET analysis in pkCSM, the required SMILES file format of the 13AApi was created using pepSMI tool (PepSMI: Convert Peptide to SMILES string [novoprolabs.com])

BBB Permeability logBB < − 1 indicates poor distribution to the brain, CNS permeability logPS > − 2 classifies central nervous system penetration and logPS < − 3 classifies no central nervous system penetration, Fraction Unbound defines amount that remains unbound to plasma protein for pharmacological action, GRAVY Grand Average of Hydropathy, Instability Index value below 40 classifies stable protein/peptide, MRTD maximum recommended tolerated dose (should be less than 0.477 mg/kg/day), pI isoelectric point, Skin Perm logkp > − 2.5 classifies low skin permeability, Total Clearance includes both hepatic and renal clearance, Water Solubility logS defines solubility in water at 25 °C

Evaluation of ADMET properties suggested good water solubility with low skin permeability of 13AApi. The ADMET analyses also suggested that 13AApi was unlikely to have blood–brain barrier and central nervous system permeability. In terms of toxicity, the prediction suggested that 13AApi was safe and bore no skin toxicity or carcinogenicity (Ames test). Furthermore, the predicted maximum recommended tolerated dose for humans was below the maximum set value of 0.477 mg/kg/day (Table 2). Other properties of 13AApi are shown in Table 2.

13AApi Binds to SARS-CoV-2 and Blocks the ACE2 interaction

To validate the in silico findings of an efficient inhibition in the interaction between the SARS-CoV-2 spike and ACE2 receptor protein by 13AApi, we used an ELISA-based ACE2: SARS-CoV-2 spike inhibitor screening assay. Our data suggest that 13AApi could significantly inhibit the interaction between the SARS-CoV-2 spike and immobilized ACE2 in a concentration-dependent manner. Greater than 40% inhibition of interaction was observed at a 13AApi concentration of 100µM (Fig. 5A, B). This strongly suggests that the 13AApi can be considered for further investigation as a new modality for the therapeutic intervention of SARS-CoV-2 infection.

Fig. 5

In vitro validation of the inhibitory effect of the 13-amino acid peptide inhibitor (13AApi) on severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) spike receptor-binding domain (RBD) was examined with an enzyme-linked immunosorbent assay (ELISA) test kit. A The image reproduced from the BP Biosciences ELISA kit demonstrates the steps in the assay. B Percent inhibitory activity of the 13AApi at the indicated concentrations on the binding of SARS-CoV-2 spike RBD to angiotensin-converting enzyme 2 (ACE2). The plotted graph represents the mean of three biological replicates

Conclusions

Several advantages of peptides, including ease of synthesis and modification, low toxicity, and high target specificity and selectivity, motivated us to design a potential therapeutic peptide candidate against COVID-19. The 13AApi is a novel potential therapeutic peptide that could inhibit the interaction of the SARS-CoV-2 spike protein with ACE2, thereby blocking cellular entry of the virus (Fig. 6). Our findings suggest that computationally designed inhibitory peptides may be developed as an anti-SARS-CoV-2 agent. Inhibitory activity of the designed 13AApi could be successfully validated with an ELISA. The 13AApi also displayed multiple drug-like properties such as good solubility, stability, and selectivity towards the RBD of the SARS-CoV-2 spike protein and can be further developed as a selective therapeutic candidate for COVID-19. The present study is the first step towards developing candidate drugs that can block SARS-CoV-2 entry into the host cell and further studies are warranted to determine the real-time therapeutic application of these candidate drugs.

Fig. 6

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and host angiotensin-converting enzyme 2 (ACE2) receptor interaction and mechanism of virus entry (left panel). A novel therapeutic peptide 13-amino acid peptide inhibitor (13AApi) blocks the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein from interacting with the host cell surface ACE2 (right panel). The 13AApi is a promising therapeutic candidate to prevent virus entry and infection in the host. The image was created with BioRender

Below is the link to the electronic supplementary material.

Supplementary file1 (DOC 60 kb)

Coronavirus disease 2019 is an ongoing pandemic caused by severe acute respiratory syndrome coronavirus (SARS-CoV)-2. The spike protein of SARS-CoV-2 interacts with the angiotensin-converting enzyme 2 receptor present on the host cell surface, leading to viral entry and infection.

A 13-amino acid novel peptide inhibitor (FLDKFNHNFKDLF) has been designed to block the receptor-binding domain of the SARS-CoV-2 spike protein, which impedes its interaction with the host angiotensin-converting enzyme 2 receptor.

The 13-amino acid peptide inhibitor is expected to restrict virus interaction and fusion with the host cells, thereby making it a promising therapeutic candidate to combat coronavirus disease 2019.