COVID-19 spike-host cell receptor GRP78 binding site prediction

Abstract

Introduction

In late December 2019, it was noticed that several people in Wuhan city of Hubei Province, China, were suffering from SARS-like pneumonia, which the World Health Organization (WHO) would later name COVID-19.^{1 , 2} According to the WHO surveillance draft, in January 2020, any resident or citizen in transit through Wuhan city 14 days before the onset of the symptoms is suspected to be infected by COVID-19.^{2 , 3} Additionally, WHO distributed interim guidance for laboratories that carry out the testing for the newly emerged outbreak and infection prevention and control guidance.^{4 , 5} The COVID-19 virus is suspected of having emerged in an unknown animal (perhaps a bat) and to have subsequently been transmitted to humans in the seafood and wild animal market.¹ All over the world, there are surveillance borders to prevent the spread of the new unknown coronavirus, while some countries stopped flights to and from China.⁶ By the first week of the year 2020, 41 cases were confirmed to be COVID-19 positive, leaving one person dead and seven in critical care.⁷ This number is continuously increasing on a daily basis. The number of confirmed cases at the time of writing this manuscript exceeded 77,000, and there were confirmed more than 2200 deaths, mainly in mainland China. On January 20, 2020, the National Health Commission of China confirmed the human-to-human transmission of the new coronavirus outbreak.⁷ Ten days later, WHO declared COVID-19 as a Public Health Emergency of International Concern (PHEIC). COVID-19 symptoms include fever, malaise, dry cough, shortness of breath, and respiratory distress.¹

COVID-19 is a member of Betacoronaviruses, like the former human coronaviruses SARS and MERS.^{8 , 9} With this novel human coronavirus, there are now seven different strains of Human coronaviruses (HCoVs), namely, 229E and NL63 strains of HCoVs (Alphacoronaviruses), OC43, HKU1, SARS, MERS, and COVID-19 HCoVs (Betacoronaviruses).^{1 , 8 , 10} SARS and MERS HCoV are the widely known strains of coronaviruses, and each has caused about 800 deaths. According to WHO, the mortality rates for SARS and MERS HCoV are 10% and 36%, respectively.^{8 ,} 10, 11, 12, 13 COVID-19 has a 2% mortality rate, but only in a few months will we know how fast the new virus spreads.

HCoVs are positive-sense, long (30,000 bp) single-stranded RNA viruses. Two groups of proteins characterize HCoVs; structural proteins such as Spike (S) that characterize all coronaviruses, Nucleocapsid (N), Matrix (M), and Envelope (E), in addition to the non-structural proteins, such as proteases (nsp3 and nsp5) and RdRp (nsp12).^{8 , 14} The Spike protein is a crucial recognition factor for virus attachment and entry to the host cells. It is present on the virion's outer surface in a homo-trimeric state.^{14 , 15}

The Glucose Regulating Protein 78 (GRP78) or Binding immunoglobulin protein (BiP) is the master chaperone protein of the unfolded protein response (when unfolded or misfolded proteins accumulate).16, 17, 18, 19 Under reasonable conditions, GRP78 is found in the lumen of the Endoplasmic Reticulum (ER) bound to and inactivating three enzymes responsible for cell death or differentiation. These enzymes are Activating Transcription Factor 6 (ATF6), Protein kinase RNA-like Endoplasmic Reticulum Kinase (PERK), and Inositol-requiring Enzyme 1 (IRE1).¹⁴ Above a threshold of accumulated unfolded proteins, GRP78 releases ATF6, PERK, and IRE1, leading to their activation. Inhibition of protein synthesis and enhancement of the refolding is the end result of the enzymes’ activation.^{14 , 20} Overexpression of GRP78 is also initiated upon cell stress, which increases the chance for GRP78 to escape ER retention and translocate to the cell membrane. Once translocated to the cell membrane, GRP78 is susceptible to virus recognition by its substrate-binding domain (SBD), and it can mediate the virus entry in the cell.¹⁴ Pep42 is a cyclic peptide that has been reported to bind the GRP78 overexpressed and expressed at the surface of cancer cells.²¹

In this study, the spike protein of COVID-19 was modeled using solved structures in the protein data bank.²² After model validation, molecular docking was performed to test its binding affinity against GRP78. We hypothesized that GRP78 binds to COVID-19, as it happens in the case of the MERS-CoV coronavirus,²³ and we tried to predict the binding site using the similarity between Pep42 and the COVID-19 Spike protein.¹⁴ Four regions of the spike were predicted to be the binding site to GRP78 based on sequence and structural similarity. The results are promising and suggest the possible recognition of the COVID-19 spike by the cell-surface GRP78 upon cell stress.

Materials and methods

All human coronaviruses (HCoV) spike proteins were downloaded from the National Center for Biotechnology Information (NCBI).²⁴ Multiple Sequence Alignment (MSA) of all of the sequences was performed using the Clustal Omega web server of the European Bioinformatics Institute (EMBL-EBI).²⁵ ESpript 3 software was used to represent the MSA.²⁶ The pairwise sequence identity between COVID-19 spike protein and each of the other HCoV spike proteins was calculated using the same method.

A model was built for the COVID-19 spike using the SWISS-MODEL web server.²⁷ The SARS HCoV spike (PDB ID: 6ACD, chain C), as the closest spike protein to COVID-19 among all of the human coronaviruses, was used as a template for building the COVID-19 spike. The structural superposition of the SARS spike (6ACD, chain C) and COVID-19 spike model was performed using PyMOL software.²⁸

Four regions are selected from the COVID-19 spike protein receptor-binding domain (RBD). These regions are cyclic (starting and ending with Cysteine residues connected by a disulfide bond) and were hypothesized to be a possible binding site to GRP78 based on the alignment with the Pep42 cyclic peptide. Pairwise sequence alignment was performed for the Pep42 peptide against the four predicted regions of the COVID-19 spike utilizing the Clustal Omega web server. ProtScale web server of the ExPASy bioinformatic resource portal was used to compare the Pep42 sequence to the four regions of the COVID-19 spike that we hypothesized to be the binding site to GRP78.²⁹ The Kyte & Doolittle hydrophobicity index was calculated for each residue of the peptides, while the grand average hydrophobicity (GRAVY) was calculated for each region of the spike and Pep42.

Molecular docking was utilized to test the binding affinity of the four regions of the COVID-19 spike against the solved structure of GRP78 Substrate Binding Domain (SBD) β. The only solved structure of wild-type, full-length GRP78 bound to ATP in the open configuration found in the PDB database is 5E84.^{30 , 31} The coordinates of 5E84 were downloaded and prepared for the docking experiment by removing water ions and the ligand. The solvated docking software, HADDOCK, was used in this study to dock the four regions of the spike model for COVID-19 against the solved structure of GRP78.³² The easy interface was utilized since no restraints are defined.³³ GRP78 active residues (I426, T428, V429, V432, T434, F451, S452, V457, and I459) were retrieved from the literature.³⁰ The active residues from the COVID-19 spike protein were chosen for each region of the spike to be the hydrophobic residues. In both proteins, the residues surrounding the active residues were selected as passive in HADDOCK. Active residues are the amino acid residues from the two interacting proteins’ binding sites that take part in direct interaction with the other protein partner while passive residues are the residues that can interact indirectly in HADDOCK.33, 34, 35, 36

Furthermore, Pep42 is tested against the GRP78 solved structure to compare its binding affinity to that of the four regions of the spike protein. The Pep42 3D structure was generated using the I-Tasser web server.³⁷ Pep42 was treated as cyclic during the docking experiment (a distance restraint is added to HADDOCK) since the cyclic form of Pep42 is the selectivity determinant against GRP78 recognition.^{38 , 39}

After docking, the docking complexes were analyzed by the aid of the Protein-Ligand Interaction Profiler (PLIP) web server of Technical University, Dresden.⁴⁰ Two main types of interactions are established upon docking: H-bonding and hydrophobic interactions. PRODIGY software was used to predict the binding affinity for each region of the spike to GRP78. The average docking scores and the residues that take part in the interactions are discussed.

Results and discussion

Sequence and structural alignment

Fig. 1 A shows part of the Multiple Sequence Alignment (MSA) between the seven human coronaviruses (229E, NL63, OC43, HKU1, SARS, MERS, and COVID-19) spike proteins performed by Clustal Omega web server and visualized by ESpript software. The secondary structure for the COVID-19 spike model is displayed at the top of the MSA, and residual surface accessibility is present at the bottom. Alpha helices are shown by helix while arrows show beta-sheets on the top of the MSA. The residues that are surface accessible are in blue, while buried residues are in white at the bottom of the MSA. Identical residues are highlighted in red, while similar residues are highlighted in yellow. The positions of the disulfide bonds are marked by the green numbers below the accessibility rows in the MSA. 13 disulfide bonds are found in the spike protein from which we predict four regions to be the binding site with cell surface GRP78. These four regions of the spike protein, identified with the disulfides numbers 3, 4, 5, and 6 are marked in the MSA with green, blue, magenta, and red dashed lines, respectively. A complete MSA for the spike protein (1273 residues) is found in supplementary figure S1.

Open in a separate window

Fig. 1

(A) Part of the multiple sequence alignment for the spike protein of all of the currently reported human coronaviruses strains (COVID-19, SARS, MERS, NL63, 229E, OC43, and HKU1). The alignment is made using the Clustal Omega web server and is displayed by ESpript 3 software. The red highlighted residues are identical, while yellow highlighted residues are conserved among the seven HCoVs. Secondary structures are represented at the top of the MSA for the COVID-19 spike, while the surface accessibility is shown at the bottom (blue, surface accessible, cyan, partially accessible, and white for buried residues). The four regions of the spike protein are shaded with green, blue, magenta, and red for regions I, II, III, and IV, respectively. (B) Structural superposition of SARS spike structure (green cartoon) and COVID-19 spike model (cyan cartoon). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

The SARS spike protein sequence is the closest to the COVID-19 spike, with 77.38% identity. In contrast, OC43, MERS, HKU1, 229E, and NL63 share only 32.81%, 32.79%, 31.86%, 30.35%, and 28.28%, respectively, with COVID-19 spike. Fig. 1B shows the superposition of the homo-trimeric COVID-19 spike model (cyan cartoon) and SARS spike structure (PDB ID: 6ACD) (green cartoon). Two views are shown with a vertical axis rotation of 180°. The Root Mean Square Deviation (RMSD) between the two structures is only 0.284 Å, while the sequence identity is 77.38%.

Pep42 versus spike regions

Pep42 is reported to specifically target the cell-surface GRP78 in cancer cells.¹⁴ Its selectivity against GRP78 has been reported for its cyclic form but not for the extended form. This may be due to the rigidity of the cyclic structure of the peptide, which causes the stabilization of the hydrophobic patch formed by C1, V3, A4, L5, V10, V12, and C13. These residues become closer to each other by the aid of the disulfide bond, making the cyclic peptide the perfect docking platform for GRP78 SBDβ.

We found 13 disulfide bonds in the COVID-19 spike protein model that form 13 different cyclic regions that may resemble the cyclic Pep42. Four of these disulfides are found in the outer surface of the spike receptor-binding domain that faces the outside part of the virion, a region that has been targeted with neutralizing antibodies against the SARS and MERS spikes. These four regions, namely, the region I C336: C361 (26 residues), region II C379: C432 (54 residues), region III C391: C525 (135 residues), and region IV C480: C488 (9 residues), are marked in Fig. 1A.

Fig. 2 A shows the pairwise sequence alignments between each spike region and Pep42. The percentage of pairwise sequence identity is listed on the right side of each alignment. The percent identity for region III is the most significant (46.15%) compared to other regions (15.38%, 23.08%, and 33.33% for the regions I, II, and IV, respectively). As shown in Fig. 2A, region IV is part of region III. Moreover, regions II and III share some residues. Again, identical residues are highlighted in red, while the conserved residues are highlighted in yellow. The secondary structure is shown at the top of the alignment and the surface accessibility at the bottom. Region IV has all of its residues exposed at the surface (either blue, meaning surface accessible, or cyan, for partially accessible residues). For other regions, some residues are surface exposed (blue or cyan), while others are buried (in white).

Open in a separate window

Fig. 2

(A) Pairwise sequence alignment between Pep42 and four regions of the COVID-19 spike protein model. Secondary structures and surface accessibility are shown at the top and the bottom of the alignment, respectively. The percentage of identity arrows shows beta-sheets on the right side of each alignment. Identical and similar residues are highlighted in red and yellow, respectively. (B) The hydrophobicity index (Kyte & Doolittle) for the four regions of the spike protein model for COVID-19 and the Pep42 peptide. Grand average hydrophobicity (GRAVY) is shown on the right side for each peptide. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Fig. 2B shows the hydrophobicity index (Kyte & Doolittle) for each suggested region and the Pep42. The grand average hydrophobicity index for each region (GRAVY) is listed in front of each peptide. Regions I, II, and III have negative values of GRAVY (−0.24, −0.30, and −0.28, respectively). In contrast, region IV has a positive value (0.08), which means that it has a slightly more hydrophobic character compared to other regions. Pep42 has a highly hydrophobic character (GRAVY value of 1.1) that enables it to be recognized by the cell-surface GRP78.^{14 , 38 , 39}

Fig. 3 A shows the structure model of the COVID-19 spike protein model (homo-trimeric) in a colored cartoon representation. Region IV of the spike (C480: C488) is not only cyclic but also surface-accessible and protrudes to the outer side of the spike, i.e., facing the target cell. It has a slightly hydrophobic character, hence resembling the Pep42 cyclic peptide, and it seems suitable to be the binding site to the cell-surface GRP78. Fig. 3B shows the region III of the spike (black cartoon). As shown in the enlarged panel, region IV is part of region III and it is the most surface-exposed part of the spike receptor-binding domain. The contribution of region IV (C480: C488) in binding region III to GRP78 is high (−9.8 out of −14.0 kcal/mol)

Open in a separate window

Fig. 3

(A) The structure of the spike protein model of COVID-19 in its homotrimer state (colored cartoon). Two chains, A (green) and B (cyan) are in the closed conformation, while chain C (magenta) is the open configuration that makes it able to recognize the host cell receptor. Region IV of the spike (C480-C488), which we suggest is the recognition site for cell-surface GRP78, is shown in the black cartoon in the enlarged panel. (B) The structure of the region III (black carton). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Binding mode of spike-GRP78

GRP78-COVID-19 spike protein docking was performed using the HADDOCK software in four different ways. Each region of the spike (predicted to be the binding site to GRP78) was used as the binding site to GRP78, using its active residues selected to be that have hydrophobic character. The active residues for region I are: C336, F338, V341, F342, A344, R347, F348, S350, A352, I358, and C361. For region II, the active residues are: C379, V382, L387, I390, C391, F392, V395, A397, F400, V401, I402, V407, I410, A411, I418, A419, L425, F429, and C432. For region III, the active residues are: C391, F392, V395, A397, F400, V401, I402, V407, I410, A411, I418, A419, L425, F429, C432, V433, I434, A435, L441, V445, L452, L455, F456, L461, F464, I468, I472, A475, C480, V483, F486, C488, F490, L492, F497, V503, V510, V511, V512, L513, F515, L517, L518, A520, A522, V524, and C52. Finally, for region IV, the active residues are C480, V483, F486, and C488. For GRP78, the active site residues are retrieved from previous work to be I426, T428, V429, V432, T434, F451, S452, V457, and I459.

Figs. 4 A and ​and4B4B show the binding mode for each of the docking trials (the best-formed complexes from each docking experiment) with green cartoon representing GRP78 and yellow cartoon representing the homo-trimeric COVID-19 spike. All the docking trials proved the possibility of fitting the GRP78 SBDβ to the spike with binding affinities (predicted by PRODIGY) ranging from −9.8 up to −14 kcal/mol. In terms of the orientation of the two interacting proteins, regions III and IV are accepted to be the docking platform. In these two trials, the GRP78 and spike can interact in a head-to-head fashion.

Open in a separate window

Fig. 4

The structure of the docking complexes of GRP78 (green cartoon) and COVID-19 spike (yellow cartoon) regions I and II (A) and regions III and IV (B). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Table 1 summarizes the docking trials of the four regions of the COVID-19 spike protein against GRP78. The docking scores are listed, while the interaction pattern is analyzed by PLIP software and listed in the table. As shown from the docking scores, region IV of the spike is the best docking platform to GRP78, with a score of −143.5 ± 4.4. This score is lower (better) than other regions by 18.3%, 15.2%, and 65.5% for the region I, region II, and region III, respectively. The PRODIGY binding affinities are also listed in Table 1. The PLIP analysis partially explains the binding affinity. Region IV of the spike interacts with the substrate-binding domain β of GRP78 with five H-bonds (through P479, N481, E484, and N487) and four hydrophobic interactions (through T478, E484, and F486). The average H-bond length for the docking trial of region IV is 2.26 ± 0.54, while the average hydrophobic contact length is 3.66 ± 0.18. These values are less than other docking trials using other regions.

Table 1

The interactions formed between the spike protein of COVID-19 and cell-surface GRP78 SBDβ upon docking with HADDOCK.

Pep42 conformation at	HADDOCK score	PRODIGY binding affinity (kcal/mol)	H-bonding			Hydrophobic interaction
			number	Residues from the spike	Residues from GRP78	number	Residues from the spike	Residues from GRP78
Region I	−121.3 ±7.2	−12.4	11	R466 N343 N354 N354 D467 E340 G339 I468 I468 K356 T345	T434 G489 S452 V453 T434 V490 G489 T434 K447 G454 G489	7	E465 I468 W353 W353 A344 R346 R466	V432 L436 T428 V429 V490 Q492 T428
Region II	−124.6 ±9.8	−10.1	12	Y396 Y396 R357 R357 R466 R466 R466 F464 E516 E516 K356 N354	S452 S452 Q449 I450 G454 T456 T456 G454 S452 Q492 T428 G430	10	E340 E516 K356 K356 F429 W353 N354 R355 R355 R357	V432 V490 T428 V429 V490 V453 V429 F451 V453 F451
Region III	−86.7 ±4.6	−14.0	9	N460 N487 D420 L455 K417 Y421 Y421 Y473 Y505	Q492 K447 S452 Q449 E427 I450 I450 S448 G430	3	E406 K417 F486	V429 V429 T441
Region IV	−143.5 ±4.4	−9.8	5	E484 E484 N481 P479 N487	T428 V429 S452 S452 T458	4	E484 E484 F486 T478	T428 V429 V457 V490

Open in a separate window

Fig. 5 A shows the predicted binding mode of the GRP78 (cyan surface) to the spike of the newly emerged coronavirus (green surface) using region IV of the spike (red surface) as the docking platform. The interacting residues of GRP78 and the spike proteins are shown in yellow and red, respectively. The enlarged panel shows in more detail the interacting residues of GRP78 (yellow sticks) and the region IV of the spike (red cartoon) labeled with its one-letter codes. This binding mode is acceptable since the two proteins are interacting, as when the virus is approaching the target cell (respiratory system cells) expressing its cell-surface receptor, GRP78. Fig. 5B shows a hypothetical binding model showing the homotrimer spike (red surface) protein of the COVID-19 bound to a respiratory system cell exposing the GRP78 protein (green surface). This scenario could occur in stressed cells when GRP78 is overexpressed and translocated from the Endoplasmic Reticulum (ER) to the cell membrane (M).

Open in a separate window

Fig. 5

(A) The proposed binding mode of the host cell GRP78 (cyan surface) and the COVID-19 spike model (green surface) through region IV (C480-C488) (red surface). The amino acids from the GRP78 SBDβ that interact with the spike protein region IV (red cartoon) are labeled and represented in yellow sticks in the enlarged panel. (B) The proposed recognition mode of the COVID-19 spike (red surface) and cell-surface GRP78 (green surface) through the spike protein region C480-C488. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

The predicted binding site of the Spike protein to GRP7 found in this study is in good agreement with studies that identified the spike receptor-binding domain using antibodies.^{23 , 41} Knowledge of this binding site could open the door for further experimental and simulation studies on the mode of envelope protein recognition by the highly dynamical GRP78 substrate-binding domain.

Conclusion

Spike protein is an essential viral element that helps in the attachment and virus internalization to the host cell. A vast amount of host cell receptors are targets for viruses, including the cell-surface GRP78. Inhibiting the interaction that occurs between the COVID-19 spike protein and the host cell receptor GRP78 would probably decrease the rate of viral infection. Furthermore, a vaccine against the COVID-19 spike protein would likely prevent viral infection. The present in silico perspective suggests the existence of a COVID-19 spike protein-GRP78 binding site, thus paving the route for drug designers to develop suitable inhibitors to prevent the binding and hence the infection. Future work involving the dynamics of GRP78 and the experimental validation is required to suggest potent peptidomimetic inhibitors.

Declaration of Competing Interests

All of the authors declare that there is no competing interest in this work.

Acknowledgments

Footnotes

References