Next Article in Journal
In Silico Prediction and Selection of Target Sequences in the SARS-CoV-2 RNA Genome for an Antiviral Attack
Previous Article in Journal
Effect of Insertion and Deletion in the Meq Protein Encoded by Highly Oncogenic Marek’s Disease Virus on Transactivation Activity and Virulence
Previous Article in Special Issue
Uniqueness of RNA Coliphage Qβ Display System in Directed Evolutionary Biotechnology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phage-Displayed Mimotopes of SARS-CoV-2 Spike Protein Targeted to Authentic and Alternative Cellular Receptors

by
Valery A. Petrenko
1,*,
James W. Gillespie
1,*,
Laura Maria De Plano
2 and
Michael A. Shokhen
3
1
Department of Pathobiology, College of Veterinary Medicine, Auburn University, Auburn, AL 36849, USA
2
Department of Chemical, Biological, Pharmaceutical and Environmental Sciences, University of Messina, 98122 Messina, Italy
3
Department of Chemistry, Bar Ilan University, Ramat Gan 52900, Israel
*
Authors to whom correspondence should be addressed.
Viruses 2022, 14(2), 384; https://doi.org/10.3390/v14020384
Submission received: 11 November 2021 / Revised: 8 February 2022 / Accepted: 10 February 2022 / Published: 14 February 2022

Abstract

:
The evolution of the SARS-CoV-2 virus during the COVID-19 pandemic was accompanied by the emergence of new heavily mutated viral variants with increased infectivity and/or resistance to detection by the human immune system. To respond to the urgent need for advanced methods and materials to empower a better understanding of the mechanisms of virus’s adaptation to human host cells and to the immuno-resistant human population, we suggested using recombinant filamentous bacteriophages, displaying on their surface foreign peptides termed “mimotopes”, which mimic the structure of viral receptor-binding sites on the viral spike protein and can serve as molecular probes in the evaluation of molecular mechanisms of virus infectivity. In opposition to spike-binding antibodies that are commonly used in studying the interaction of the ACE2 receptor with SARS-CoV-2 variants in vitro, phage spike mimotopes targeted to other cellular receptors would allow discovery of their role in viral infection in vivo using cell culture, tissue, organs, or the whole organism. Phage mimotopes of the SARS-CoV-2 Spike S1 protein have been developed using a combination of phage display and molecular mimicry concepts, termed here “phage mimicry”, supported by bioinformatics methods. The key elements of the phage mimicry concept include: (1) preparation of a collection of p8-type (landscape) phages, which interact with authentic active receptors of live human cells, presumably mimicking the binding interactions of human coronaviruses such as SARS-CoV-2 and its variants; (2) discovery of closely related amino acid clusters with similar 3D structural motifs on the surface of natural ligands (FGF1 and NRP1), of the model receptor of interest FGFR and the S1 spike protein; and (3) an ELISA analysis of the interaction between candidate phage mimotopes with FGFR3 (a potential alternative receptor) in comparison with ACE2 (the authentic receptor).

1. Introduction

The evolution of the SARS-CoV-2 virus during the COVID-19 pandemic was accompanied by the emergence of a diverse assortment of new virus variants having multiple mutations in both receptor- and antibody-binding sites located throughout the spike S protein [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25] (Figure 1). Analysis of a variety of coronavirus variants that caused epidemics during the last 60 years showed that their evolution may occur not only through the appearance and selection of immuno-resistant mutants, but also by the selection of variants with altered spike proteins that can use alternative, or “reserve”, human host receptors for entering into target cells [18,26,27,28,29]. Our current understanding of the molecular mechanisms of virus adaptation to human host cells [30] and immuno-resistant human populations is not sufficient to rationally design efficient countermeasures against viral pandemics, such as the SARS-CoV-2 virus. To empower increased insight into the evolutionary mechanisms of the SARS-CoV-2 spike protein, we suggest using recombinant filamentous bacteriophage-based probes [31], displaying on their surface an array of ~4000 foreign peptides, which produces a unique molecular landscape across the viral surface that can mimic the structure of viral receptor-binding sites [32,33,34,35,36,37], as illustrated in Figure 2 and Figure 3.
Phage display, a technique developed over the past 35 years, employs the use of engineered bacteriophages to display a genetically encoded peptide sequence fused onto an exposed structural protein of a bacteriophage particle, such as fd-tet [31]. Landscape phage libraries are composed of a collection of different bacteriophages or “phages” in which a fused protein is displayed on the N-terminus of the p8 major coat protein as opposed to the p3 minor coat protein, as seen with many commonly used phage display libraries [43,44]. Due to the structural modifications introduced by the displayed peptide fusions to the 4000 copies of the p8 protein, each landscape phage particle can be treated as a unique nanomaterial with novel and emergent properties that cannot be observed by use of an individual synthetic peptide alone [45]. In this p8-type phage display system, the array of densely packaged foreign peptides composes a unique landscape in which the structure and function of individual peptides can be influenced by interactions with neighboring peptides and the body of the phage particle [45,46]. For example, it was shown that synthetic peptides corresponding to phage-displayed peptides on landscape phages can demonstrate very different activity in competition ELISA. This property of landscape phages demonstrates the potential of each phage to be a selectable nanomaterial, as opposed to simply a means to discover binding peptides. The discovery of landscape phage mimotopes specific towards viral attachment proteins, such as the SARS-CoV-2 spike protein, provides a unique nanomaterial that can be translated to various applications, such as vaccine products, biosensors, or other analytical devices [44].
Molecular or antigenic mimicry can be used to describe the similarity of different antigens that can cause the cross-recognition of conformational peptides/proteins or amino acid sequence to a receptor, commonly observed between pathogens and immune receptors that triggers an undesired autoimmune response [47,48]. We extend this concept and use in silico computational tools previously developed for studying phage-displayed peptide sequences discovered between antibody-antigen interactions [49].
In this study, we sought to discover phage mimotopes of the SARS-CoV-2 spike S protein using the concepts of phage display and molecular mimicry [31], termed here “phage mimicry.” Three essential steps in the phage mimicry paradigm can be defined: (1) enrichment of a multibillion clone population of landscape phage libraries for thousands of phage variants that interact with cellular receptors exposed on the surface of mammalian normal and cancer cells; (2) identification of computationally predicted amino acid clusters with similar 3D conformations as naturally occurring ligands of the target receptor; and (3) evaluation of candidate phage mimotopes for their ability to interact with a receptor of interest in vitro or in vivo using an ELISA and other diagnostic methods in comparison with phages displaying a non-related peptide.
Here, we used a collection of landscape phages that were previously selected by our research group through enrichment of landscape phage display libraries for cancer cell binding/penetrating phages using a panel of human cancer cell lines, which overexpress a variety of different receptors, including receptors that presumably can be involved in viral pathogenesis. Using the PepSurf algorithm, which maps a set of affinity-selected peptides onto the surface of a solved PDB structure, the linear peptide sequence that mimics the genuine epitope (mimotopes) of the antibody-antigen complex can be identified. The discovered linear peptide sequences can therefore mimic the complex 3D conformational folding of the native ligands. We hypothesized that instead of using peptide sequences discovered against an antibody to discover the molecular mimotopes of the antigen, as was presented in the original manuscript, we could use a library of preselected phages enriched towards cellular receptors to discover viral protein molecular mimotopes involved with viral attachment. Next, we studied the interaction of phage mimotopes mimicking amino acid clusters in the SARS-CoV-2 S1 protein that bind the canonical receptor ACE2 as a proof-of-concept model. We identified amino acid clusters on two different ligands (SARS-CoV-2 S1 protein, FGF1, and NRP1), which presumably interact with the FGFR3 receptor. We predicted conformationally similar domains of phage mimotopes interacting with the same cellular receptor FGFR. Finally, we analyzed candidate phage mimotopes for binding with recombinant FGFR3.

2. Materials and Methods

A panel of phage mimotopes of the SARS-CoV-2 spike S protein that bind certain cell receptors was developed using the phage mimicry technique [31,50,51], supported by a combination of bioinformatic methods [49,50,51,52] and validated by experimental methods.

2.1. Preparation of Libraries of Phages That Display Receptor-Binding Mimotopes Fused to the Major Coat Protein p8

The source of SARS-CoV-2 mimotopes was the p8-type polyvalent peptide phage-displayed library (or, shortly, the landscape phage library) [43,44,53]. Bacteriophage fd (Figure 2), which itself has no natural tropism to mammalian cells, is a suitable vector for generating random peptide phage-displayed libraries—a rich source of ligands for cellular receptors. It was shown that phages selected from these libraries can specifically recognize cellular receptors and penetrate into sub-cellular compartments during their artificial molecular evolution in vitro and in vivo, similarly to evolution observed with naturally evolved viruses. We demonstrated that selection of organ/cell-associated phage variants from their multibillion clone libraries and bioinformatic analysis of their cell-binding peptides in comparison with proteins of natural viruses allow the identification of functional virus-host binding sites that are apparently used during viral pathogenesis. The phages in these libraries have been screened against a variety of human cancer and normal cells harboring cellular proteins in their native structure-functional state [54]. Specifically, in this project we used an assortment of several hundred phage clones isolated from human lung [55], prostate [53], pancreatic [56], and breast cancer [57] cells and characterized for specificity and selectivity towards corresponding cell lines.

2.2. Mapping of Epitopes and Receptor-Binding Sites on the Surface of Spike S Protein

From a variety of computer programs that were prepared for studying antibody-antigen interactions, we chose PepSurf [49], which originally was designed for mapping a set of affinity-selected peptides onto the resolved structure of the antigen, and was adapted in this work for the discovery of receptor-binding mimotopes of SARS-CoV-2 S protein. The important feature of this algorithm is that it assumes that the peptides mimic surface residues (i.e., solvent exposed residues), and buried residues are eliminated from the search. In this application, the S protein played a role of “antigen” and mimotopes interacting with virus receptors were identified by aligning them against all possible paths in the graph and clustering the most significant matches in which a predicted epitope and/or receptor-binding site is implied. To discover phage-displayed mimotopes of epitopes and receptor-binding clusters on the surface of the spike protein, the PepSurf server was used, as recommended for epitope mapping, with the standard input of the Protein Data Bank (PDB) files of the S protein with a solved 3D structure available from the PDB (6VYB and 6M0J) and a set of peptides discovered by screening of landscape phage-displayed libraries against human cancer cells in vitro and in vivo. The first output of the program is the alignment of each peptide to the 3D structure of the spike. When several peptides were aligned, the server implemented a clustering algorithm to detect one or more patches of residues on the surface of the surveyed protein. Thus, the second output is the predicted patches. Such a patch may correspond to a putative epitope site on the spike, or as assumed in our work, a receptor-binding site.

2.3. Molecular Modeling

The 3D structure of a landscape phage displaying the peptide EDYSELVSQ as an N-terminal fusion to all copies of the mature p8 major coat protein was generated by SWISS-MODEL, the homology modeling server [58,59,60,61,62]. The full 55 amino acid residue sequence of p8 (AEDYSELVSQPAKAAFDSLQASATEYIGYAWAMVVVIVGATIGIKLFKKFTSKAS) was used as the input query, where the residues highlighted in bold indicate the displayed N-terminal peptide fusion on every copy of p8. The 3D structure of 1FDM, from the PDB, was identified as a template by a search performed against the SWISS-MODEL template library with BLAST [63] and HHBlits [64]. The resulting 3D structure was a fragment formed by 55 coat proteins of filamentous phage, where the full-length of the latter is ~2700 copies of the coat proteins. It should be noted that the model generated by homology modeling coat proteins had a reduced sequence containing 6–55 residues only from the target model. All 55 reduced-coat proteins were substituted for the full-length intact variant by means of molecular modeling procedures implemented in the YASARA Structure software package [65,66]. Finally, the 3D structure of the phage fragment with corrected protein sequences was optimized under an AMBER ff14SB force field [67] in a periodic solvent box with explicit water molecules and physiological concentrations of Na+ and Cl ions.
The 3D structures of the complexes formed by the fibroblast growth factor receptor 3 (FGFR3) with ligands phage, FGF1, NRP1, and SARS-CoV-2 spike RBD, presented in Figure 8, were predicted by the HADDOCK web server for protein-protein docking [68,69]. The 3D structures of receptor and ligands used as input data for the docking procedure were acquired from the following crystal structures available from the PDB: FGFR3 (1RY7), FGF1 (1RY7), NRP1 (2QQN), and the SARS-CoV-2 spike protein RBD (6M0J).

2.4. Antibodies and Human Recombinant Proteins

The following commercially available antibodies were used in this study: rabbit; the anti-SARS-CoV-2 spike protein mAb (Sino Biological, Wayne, PA, USA #40150-R007, RRID: AB_2827979); and HRP-conjugated goat anti-rabbit IgG (H+L) polyclonal antibody (Jackson ImmunoResearch, West Grove, PA, USA #111-035-045, RRID: AB_2337938).
The following human recombinant proteins were obtained from Sino Biological, Wayne, PA, USA (Table 1): ACE2-His (#10108-H08H), DPP4-His (#10688-H08H), FGFR3-His (#16044-H08H), and the SARS-CoV-2 S1-His spike protein (40591-V08H). All recombinant proteins used in this study were resuspended to a stock concentration of 250 µg/mL in ddH2O, according to the manufacturer’s instructions.

2.5. Phage Indirect ELISA with Human Extracellular Receptors

An indirect ELISA was used to evaluate the interaction between landscape phages displaying recombinant peptide fusions of SARS-CoV-2 spike protein mimotopes and candidate host cellular receptors. Recombinant, extracellular domains of human receptors (hACE2, hNRP1, and hFGFR3) were diluted to a final concentration of 2.0 µg /mL with a coating buffer (1× PBS, pH 7.4). Wells of a 96-well, high-binding microplate (Corning, Corning, NY, USA #9018) were coated with 100 µL of each diluted candidate receptor (200 ng/well) and incubated overnight at 4 °C.
Unbound antigen was removed, and the wells were washed three times with 200 µL of washing buffer (1× PBS, pH 7.4/0.1% Tween 20) at room temperature. Antigen-coated wells were blocked with 200 µL of blocking buffer (1× PBS, pH 7.4/2% BSA/0.1% Tween 20) for 1 h at 37 °C. The wells were then washed three times with 200 µL of washing buffer at room temperature. Candidate phages were diluted to a working concentration of 5.0 × 1011 vir/mL and tenfold serial dilutions were prepared to a final concentration of 5.0 × 105 vir/mL in blocking buffer. The wells were then treated with 100 µL of diluted phages for 1 h at 37 °C. The wells were washed four times with 200 µL of washing buffer at room temperature. Bound phages were detected using 100 µL of rabbit, anti-fd phage IgG [70,71] diluted 1/1000 in blocking buffer for 1 h at 37 °C. The wells were washed four times with 200 µL of washing buffer at room temperature. Bound rabbit IgGs were detected using 100 µL of horseradish peroxidase (HRP)-conjugated goat anti-rabbit secondary antibody diluted 1/40,000 in blocking buffer for 1 h at 37 °C. The wells were washed four times with 200 µL of washing buffer prior to adding 50 µL of 1-Step Ultra TMB-ELISA Substrate (ThermoScientific, Waltham, MA, USA). The reduction of 3,3′,5,5′-tetramethylbenzidine (TMB) substrate by the HRP-conjugated antibody was monitored over a 30 min incubation at room temperature. The reaction was stopped by adding 50 µL of Stop Solution (0.2 M sulfuric acid) to each well. Endpoint absorbances for each well were measured at 450 nm using a Synergy H1 plate reader (BioTek, Winooski, VT, USA).

2.6. Indirect ELISA with Candidate SARS-CoV-2 Spike Protein Receptors

An indirect ELISA was used to evaluate the interaction between a recombinant SARS-CoV-2 spike protein and candidate host cellular receptors. Recombinant, extracellular domains of candidate SARS-CoV-2 S1 receptors (hACE2, hDPP4, hNRP1, and hFGFR3) were diluted to a final concentration of 2.0 µg/mL with a coating buffer (1× PBS, pH 7.4). Bovine serum albumin (BSA), IgG-free and protease-free, (Jackson ImmunoResearch, West Grove, PA, USA #001-000-173, RRID: AB_2336947) was diluted to a final concentration of 2% with the coating buffer and used as a negative control sample. Wells of a 96-well, high-binding microplate (Corning, Corning, NY, USA #9018) were coated with 100 µL of each diluted candidate receptor (200 ng/well) in triplicate and incubated overnight at 4 °C. A set of wells was treated with the coating buffer as a control for plastic binding.
Unbound antigen was removed, and the wells were washed three times with 200 µL of washing buffer (1× PBS, pH 7.4/0.1% Tween 20) at room temperature. Antigen-coated wells were blocked with 200 µL of blocking buffer (1× PBS, pH 7.4/2% BSA/0.1% Tween 20) for 1 h at 37 °C. The wells were then washed three times with 200 µL of washing buffer at room temperature. The recombinant SARS-CoV-2 S1-His spike protein was diluted to a final concentration of 2.0 µg/mL with blocking buffer. The wells were then treated with 100 µL of diluted SARS-CoV-2 S1 spike protein (200 ng/well) for 1 h at 37 °C. The wells were washed four times with 200 µL of washing buffer at room temperature. Bound SARS-CoV-2 spike protein was detected using 100 µL of the rabbit anti-SARS-CoV-2 spike protein monoclonal antibody diluted 1/5000 in blocking buffer for 1 h at 37 °C. The wells were washed four times with 200 µL of washing buffer at room temperature. Bound rabbit IgGs were detected using 100 µL of horseradish peroxidase (HRP)-conjugated goat anti-rabbit secondary antibody diluted 1/80,000 in blocking buffer for 1 h at 37 °C. The wells were washed four times with 200 µL of washing buffer prior to adding 50 µL of 1-Step Ultra TMB-ELISA Substrate (ThermoScientific, Waltham, MA, USA). The reduction of 3,3′,5,5′-tetramethylbenzidine (TMB) substrate by the HRP-conjugated antibody was monitored over a 30 min incubation at room temperature. The reaction was stopped by adding 50 µL of Stop Solution (0.2 M sulfuric acid) to each well. Endpoint absorbances for each well were measured at 450 nm using a Synergy H1 plate reader (BioTek, Winooski, VT, USA).

2.7. Data and Statistical Analysis

Descriptive statistics for endpoint values were calculated for each sample. Due to the small sample size, nonparametric statistical tests were selected to analyze the differences in the data. A Kruskal-Wallis one-way analysis of variance (ANOVA) was performed using R (version 4.1.1) to identify whether there was a statistical difference in mean endpoint absorbances between candidate receptors. A post-hoc Dunnett’s test for multiple comparisons was used to identify statistically significant differences in mean endpoint absorbances when compared to a negative control protein (BSA).

3. Results

3.1. Landscape Phage Libraries

To identify landscape phage mimotopes of the SARS-CoV-2 spike protein that bind a certain cellular receptor (Figure 2) [44], phage collections were prepared by screening the multibillion parental phage libraries against a variety of human cancer cells expressing authentic functional receptors: human lung [55], prostate [53], pancreatic [56], and breast cancer [57] cells, some of which, along with other human cancer cells, have been used for the propagation of SARS-CoV and SARS-CoV-2 in cell culture [30,72,73,74,75,76,77]. We assumed that after treating target cells with the phage library and discarding unbound phages, a number of phage mimics of the SARS-CoV-2 virus remained bound to cellular receptors; for example, previously identified nucleolin [57], integrins [73], and proteinase N [73], as illustrated in Figure 3B. To discover phage-displayed mimotopes that mimic the receptor-binding AA cluster on the surface of the SARS-CoV-2 spike protein (Figure 3C), we used the PepSurf server as recommended for epitope mapping, as discussed above (Section 2.2). We evaluated the performance of phage mimicry using representatives of major classes of cell surface receptors, which mediate entry of coronaviruses into the host cells: (1) metallopeptidases, including angiotensin-converting enzyme 2 (ACE2) [1,26], and (2) growth factor receptors, including fibroblast growth factor receptor 3 (FGFR3) [35,37,78,79,80,81].

3.2. Identification of Receptor-Binding Domain Mimotopes on the SARS-CoV-2 Spike Protein

In the traditional phage-display applications, such as epitope discovery and vaccine development, the final goal is discovery of peptides that can be used as a replacement for natural protein antigens [31]. However, in most proteins, including the spike protein, the epitopes and receptor-binding sites are presented by conformational AA clusters, in which separated AAs are brought together through a specific folding of proteins. That is why the goal of our project was the development of phage mimotopes, which can be used themselves as counterparts or artificial ligands of cellular receptors and antibodies in different applications, including the study of viral evolution and the design of diagnostic systems and elements of molecular vaccines [44,82].
To identify receptor-binding domain mimotopes located on the SARS-CoV-2 spike protein, we collected a panel of over 350 phage displayed peptides that we previously enriched and characterized for binding to various human cancer cell lines [53,55,56,57]. As many receptors and growth factors relevant for increased cell proliferation are overexpressed in cancer cell lines, we hypothesized that the phage-display libraries enriched for binders to these various cancer cell types would provide a representative panel of relevant phages displaying receptor-binding peptides for use in this study. We first analyzed various sets of these cell receptor-binding peptides to generate an alignment cluster of amino acids on the surface of the full-length SARS-CoV-2 spike protein based on the 3D PDB model 6VYB. Over 100 different alignment clusters were identified on the surface of the SARS-CoV-2 spike protein using the PepSurf alignment algorithm, which defines a network of surface-accessible residue paths and aligns peptide inputs to these residue paths. Alignment paths or AA clusters that contain several peptides were scored higher and were hypothesized to be more functional. After an initial round of peptide enrichment to identify SARS-CoV-2 spike protein mimics, we performed another round of mimotope screening, in which individual peptides were used as input into the PepSurf program to identify the best surface alignment of each peptide (AA clusters) on the 6VYB 3D spike model. We then performed a hierarchical clustering of peptides with the 6VYB model of the SARS-CoV-2 spike protein (Figure 4).
Amino acid residues were grouped into 13 potentially functional domains (Figure 5). Cluster 1 appeared to contain amino acids that did not strongly cluster into any functional domain. We speculate that these amino acid residues are the result of non-specific binding of phages in the parental library enrichment procedures. The remaining 12 amino acid clusters appear to contain more residues that cluster into common domains. Some of these clusters overlap (clusters 2, 3, and 5–8) but involve additional residues that may only be functional during specific orientations/conformations of the spike protein. For example, as the spike protein transitions from the closed to open state, additional functional amino acids are exposed, but also previously functional residues may no longer be functional based on their accessibility [83]. However, several of the identified clusters define isolated regions within the spike protein, including clusters 9, 4, and 10. Some of the clusters, i.e., cluster 4, contain residues that are specific for that cluster but cover a wide distribution of positions that do not cluster around a central region. These regions most likely contain functional domains that are highly active or functionally conserved.
Following identification of spike protein functional domains using our bioinformatics approach, we observed that most of the functional amino acid residues were in the spike receptor-binding domain (RBD) as expected and chose to focus our attention on interactions located within this known functional region using the PDB model 6M0J.

3.3. Phage-Displayed Mimotopes of Spike S1 Protein Interacting with ACE2

To validate the phage mimicry approach for the discovery of phage mimotopes that correspond to functionally active sites on the surface of the SARS-CoV-2 spike protein, we focused on the SARS-CoV-2 spike protein RBD. We hypothesized that phage-displayed peptide mimotopes would share common amino acid clusters with the ACE2-interacting S protein’s AA clusters as visualized through the PDB model 6M0J [84]. We narrowed our collection of spike protein mimics to those that contained residues involved in mediating the spike-ACE2 interaction and visually confirming a potential interaction between the two proteins in YASARA. Using these criteria, we narrowed the pool of potential candidates to 55 phages. We hypothesized that peptides containing the most amino acid residues in common with the spike RBD would have highest interaction potential and therefore produce the highest binding in a functional ELISA. The peptides were then ranked based on the number of amino acids shared with the SARS-CoV-2 spike RBD, resulting in a panel of 11 phages that were used for screening in a functional ELISA towards a recombinant ACE2 protein. For example, the phage-displayed peptides DGRADLSYD, VGIDEQRAD, and DGRSIVGDE all contained 9 amino acid residues in common with the spike RBD and were located in the spike-ACE2 interaction site (Figure 6B–D). We then modeled the interaction of a landscape phage displaying the DGRADLSYD peptide on all copies of the p8 major coat protein with ACE2 to demonstrate the various sites of interaction between the phage-ACE2 complex, as illustrated in Figure 6A.
To check the prediction value of the phage mimicry strategy, we characterized the binding of phages displaying the SARS-CoV-2 spike protein RBD-mimicking peptides to a recombinant ACE2 protein in a functional ELISA. Briefly, a recombinant ACE2 protein was adsorbed to the wells of a 96-well plate and blocked with BSA. Serial dilutions of candidate phages were prepared and allowed to interact with bound ACE2 at 37 °C for 1 h. Following extensive washing, bound phages were detected using a rabbit anti-phage antibody and a goat anti-rabbit secondary antibody conjugated with horseradish peroxidase. The signal was generated using a TMB substrate and absorbance data collected over a 30 min interval to calculate the maximum velocity for each reaction well. Spike RBD-mimicking phages produced higher binding in an ELISA when compared to the unrelated parent phage, fd-tet, suggesting that the identified mimotopes serve as functional peptide mimics of the SARS-CoV-2 spike protein (Figure 7). Based on the predicted interaction location provided by the molecular modeling of AA clusters of the SARS-CoV-2 spike RBD and ACE2, the affinity of the phage mimotope to the ACE2 receptor can be modified based on differences in the amino acids involved in the interaction site.

3.4. Identification of Phage-Displayed Mimotopes of the Spike S1 Protein Interacting with FGFR3

Since there is no data in the literature concerning 3D molecular models demonstrating the hypothetical complex between FGFR and the SARS-CoV-2 spike protein that are confirmed by X-ray crystallography or electron microscopy analysis, as in the model described above for the ACE2 and SARS-CoV-2 spike protein complex, we used for our analysis the phage mimicry strategy. The phage mimicry strategy is based on principles of equivalence relations. In mathematics, if object A is equivalent to object B, and object B is equivalent to object C, then object A is also equivalent to object C, according to the transitivity relationship. In virology, molecular mimicry can be defined as a structural similarity between viral proteins and natural ligands of human cellular receptors. In relation to the discovery of phage mimotopes, which are able to compete with RBD for binding to FGFR3, we assumed that they can be found among FGFR-binding phage mimotopes of the natural ligands for FGFR3—neuropilin-1 (NRP1) and fibroblast growth factor 1 (FGF1), which form stable complexes with FGFR3 [81,85,86,87,88]. Since the 3D structure of the complex between NRP1 and FGFR3 is unknown, we used methods of molecular modeling. We found that the CUB2 domain of NRP1 has amino acid clusters with similarity to the SARS-CoV-2 spike RBD and phage mimotopes.
In a similar manner, we determined the structure of the complex between FGFR3 and the SARS-CoV-2 spike protein RBD (Figure 8). The assortment of FGFR3 phage binders was extended using the model of the complex between FGFR3 and FGF1. Mimotopes were identified using the alignment of phage peptides with the FGF1-FGFR3 complex (PDB model 1RY7). Using the same selection criteria as above, we narrowed the pool of candidate phages displaying FGFR3-binding mimotopes to 13 peptide sequences. We prepared PepSurf surface residue alignments of each candidate mimotope to three different ligands, FGF1, NRP1, and the SARS-CoV-2 spike protein RBD (Table 2). From this analysis, we identified amino acid residues that were conserved between the phage-displayed mimotope and the amino acid clusters located on the three ligands. We observed that the identified mimotopes shared several amino acid residues in common with the three ligands. For example, the mimotope DGRMTVYNE contained a shared DGR motif that was common among all three ligands. Other residues in the mimotope were varied in different positions, based on the PepSurf alignment. However, we found that all amino acid clusters identified on the different ligands shared at least five amino acid residues in common, resulting in FGFR3 putative binding mimotopes.
Using the primary protein sequence of the FGFR3 extracellular domain, we identified critical residues that mediated interactions between the receptors and the three ligands (Figure 9). We observed that several critical amino acid clusters around positions 155–175, 244–262, and 310–322 were important for the generation of the FGFR3-binding mimotope. We highlighted the amino acid residues of the phage mimotope EDYSELVSQ, which interact with FGFR3 (Figure 8A), and found that many of these residues, for example D160, E247, R248, and E320, were shared among the three ligands, which can confirm the localization of the ligand-binding site on the surface of FGFR.

3.5. Analysis of Interactions of Candidate Phage Mimotopes with FGFR3

To test the functional activity of the identified FGFR3-binding phage mimotopes, we evaluated them in an indirect ELISA assay to verify binding to the extracellular domain of the FGF3 receptor. As above, the recombinant FGFR3 protein, which contained all of the extracellular domains of the receptor, was adsorbed to the wells of a 96-well plate and blocked with BSA. Serial dilutions of candidate phages were prepared and allowed to interact with bound FGFR3 at 37 °C for 1 h. Phages were detected using a rabbit anti-phage antibody and a goat anti-rabbit secondary antibody conjugated with horseradish peroxidase, as above. FGFR3-mimicking phages produced higher binding in the indirect ELISA, when compared to the unrelated parent phage, fd-tet, suggesting that the identified mimotopes serve as functional peptide mimics of the FGFR3 cellular receptor (Figure 10). Of the identified mimotopes, the phage-displaying peptides EDYSELVSQ and ETRVEPEYD demonstrated the highest binding activity in the FGFR3-binding assay. These peptides had stronger similarities to the NRP1 and spike RBD ligands than the natural FGF1 ligand we identified from the PepSurf alignment to each PDB structure, suggesting that these peptides would mimic the binding for all three proposed ligands.

3.6. Evaluation of Candidate SARS-CoV-2 Spike Protein Receptors

We hypothesized that the SARS-CoV-2 S1 spike protein can interact with several host receptors and serves as a reservoir for additional mutations, which enables increased infection through the use of alternative host cell receptors. Here, we evaluated the binding of a recombinant SARS-CoV-2 S1 spike protein to interact with the extracellular domain of the ACE2 receptor [84], in comparison to the dipeptidyl peptidase 4 (DPP4/CD26) receptor identified as the primary receptor involved in MERS-CoV [89], neuropilin-1 (NRP1), a SARS-CoV-2 co-receptor that may facilitate improved viral entry into host cells [90,91,92], and fibroblast growth factor receptor 3 (FGFR3/CD333) [35]. Recombinant extracellular domains of candidate host receptors were adsorbed to a 96-well microplate before the addition of a recombinant SARS-CoV-2 S1 spike protein, including the receptor-binding domain (319–541), an integrin-binding motif (403–405), and an ACE2 receptor-binding domain (437–508). The bound S1 spike protein was determined using an indirect ELISA with an anti-SARS-CoV-2 spike protein antibody and an appropriate HRP-conjugated secondary. After measurement of endpoint absorbances for each sample, a Kruskal-Wallis one-way ANOVA was performed to identify differences in endpoint means across all samples (p = 8.88 × 10−3). A post-hoc Dunnett’s test was then performed against each candidate receptor to identify statistically significant differences in the average endpoint absorbance among samples, compared to a control (BSA), as summarized in (Figure 11).
As expected, the SARS-CoV-2 S1 spike protein domain demonstrated the highest binding to the canonical receptor ACE2, which contained the single extracellular domain of the host receptor, including a previously identified interaction site with the SARS-CoV-2 spike glycoprotein (30–41, 82–84, 353–357), a substrate-binding domain (345–346), an ADAM17 cleavage site (652–659), and a TMPRSS11D/TMPRSS2 cleavage site (697–716). As opposed to previous reports [89] that predicted that the SARS-CoV-2 S1 protein would have a high affinity to the DPP4 receptor from in silico methods, we demonstrated that binding is higher towards the ACE2 receptor. However, it is speculated that portions of the DPP4 binding site may still be present within the S1 protein of SARS-CoV-2 variants leaving the possibility for a receptor reversion back to using DPP4 as the primary receptor, given the appropriate selective pressure. Finally, we tested the binding of the SARS-CoV-2 S1 protein to both the NRP1 and FGFR3 extracellular domains as an example of an alternative receptor that could be used to broaden the tissue tropism of host-binding or to increase the infectivity/propagation of viral variants. In contrast to the synthetic S1 spike protein, phage mimotopes demonstrated strong specific binding to both ACE2 and FGFR3, while the S1 fragment of the viral spike protein demonstrated strong binding to ACE2 and DPP4 and only background binding to NRP1 and FGFR3. These results can demonstrate the insignificant role of FGFR3 alone as a receptor interacting with the SARS-CoV-2 spike protein, in comparison with ACE2, DPP4, and other alternative receptors of human coronaviruses. However, these findings leave a prospect for proposing the role of FGFR as a co-receptor that enhances the activity of ACE2, in a similar mechanism that was revealed for NRP1, in which overexpression enhances the infectivity of the SARS-CoV-2 virus, although the interaction of the spike RBD-NRP1 was not revealed [92,93,94,95,96]. To study this opportunity, phage mimotopes derived from different spike protein variants can be further studied in live human cell lines. In opposition to the spike-binding antibodies commonly used in studying the interaction of the ACE2 receptor with SARS-CoV-2 variants in vitro, phage spike mimotopes targeted to other cellular receptors would allow the discovery of their role in viral infection in vivo using cell culture, tissue samples, or whole organisms. For example, in novel, emerged variants of SARS-CoV-2, including the recent Omicron strain, different configurations of the spike RBD and the role of alternative receptors can be revealed using more advanced methods for studying viral-host interactions and using phage mimotopes prepared using the phage mimicry strategy presented in this study. The higher performance of phage mimotopes can be explained by the multivalency of landscape phage probes and the involvement of amino acid residues in neighboring p8 subunits of the phage in the target receptor-binding. This remarkable characteristic of the landscape phage can be used in the design of molecular devices to control virus-host interactions.

3.7. Primary Structure Analysis of Receptor-Binding AA Clusters of SARS-CoV-2

To determine the mechanism of virus selectivity towards canonical and alternative receptors, we compared AA clusters and the related phage mimotopes responsible for binding to ACE2 and FGFR3 (Figure 12). It was found that although binding of the spike RBD to FGFR3 is much weaker than binding to ACE2 and CD26 (Dipeptidyl peptidease-4 or DPP4), in a direct ELISA, the presence of amino acids R403, D405, K444, E484, and N501 that determine the interaction of S protein of SARS-CoV-2 with FGFR3, which are absent in preceding coronavirus variant SARS-CoV, may be advantageous for virus infectivity.
These observation are consistent with the hypothesis that the evolution of the virus in the human population can occur not only by escaping immune pressure or increasing the affinity of RBD towards the canonical ACE2 receptor, as commonly assumed [98], but also through adaptation of the virus to alternative receptors and co-receptors of expressed on-host cells. In the retrospective analysis, the discovered cluster of FGFR interacting with SARS-CoV-2 spike amino acids was compared with corresponding amino acids of the SARS-CoV spike. Intriguingly, all of these amino acids mutated to the potential FGFR-binding cluster (R403K, E406D, N439R, K444T, V483P, E484A, and N501T), while members of the AA cluster (D442D, N448N, Y449Y, S494D, Q493N, and L492l) corresponding to phage mimotope EDYSELVSQ were mostly equally presented on both the SARS-CoV-2 and SARS-CoV RBD domains. Considering that both the virus and phage could emerge during their interaction with human growth factor receptors overexpressed on the host Calu-3 cells, one can speculate that the SARS-CoV-2 virus emerged during adaptation of SARS-CoV to Calu-3 or other human cancer cells. These observations are consistent with the hypothesis that the evolution of the virus can occur both by escaping immune pressure in the human population, and through the exchange of viral receptors and co-receptors during virus adaptation to host cells.

4. Discussion

4.1. Key Role of Human Protein Receptors in SARS-CoV-2 Infection and Evolution

The evolution of human coronaviruses during recent decades can be followed by comparing seven variants, starting from HCoV-229E (1962), and continuing to the last outbreak of SARS-CoV-2 (2019) [4] and its mutated variants. In particular, the emergence of the highly pathogenic coronaviruses SARS-CoV and MERS-CoV [29] inspired significant activities in uncovering the molecular mechanisms of viral infectivity and pathogenesis [33,37,99,100,101]. The spike S protein (Figure 1) was identified as a mediator of virus and host cell binding [102], while human angiotensin-converting enzyme 2 (hACE2) was recognized as a receptor of SARS-CoV and SARS-CoV-2. Furthermore, several other receptors and cellular proteases were found to support viral entry in the host cells [1,101]. Intriguingly, the spike RBD and appropriate human protein receptors diverge among human coronaviruses [37]. For example, ACE2 serves as a receptor for NL63, SARS1, and SARS2 strains; aminopeptidase N (CD13) serves as a receptor for the 229E strain; dipeptidyl peptidase-4 (DPP4 or CD26) serves as a receptor for MERS [29]; while peptidases ANPEP, ENPEP, and AGTR2 are related to SARS-CoV-2 infectivity [37].

4.2. Alternative Receptors of Coronaviruses as an “Operational Reserve” for the Evolutionary Escape of Emerging SARS-CoV-2 Variants

Despite minor changes in the spike conformation of human coronaviruses, single mutations in their S proteins can dramatically change the pathogenicity [99]. For example, the enhanced infectivity of SARS-CoV-2 can be attributed to the use of alternative receptors by the virus during the evolution of its predecessors [37]. Thus, the lower infectivity but higher mortality rates of SARS-CoV—the source of the epidemic in 2003 [40]—in comparison with SARS-CoV-2, can be explained in part by the appearance of the mutation K403R in the spike protein of the SARS-CoV-2 variant [84] (Figure 12). This mutation resulted in the creation of an integrin-binding motif, RGD, which is crucial for a successful infection that is known to drive the infectivity of a broad spectrum of viruses [103]. The presence of the RGD site was also thought to be a factor that can boost the affinity of SARS-CoV-2 to ACE2-positive target cell, and ACE2-negative cells. Thus, the RGD motif in the spike protein of SARS-CoV-2 could be critical in infecting cells through the RGD-binding integrins [36]. In their turn, RGD-binding integrins are involved in the activation of transforming growth factor beta (TGF-β), which can be related to complications in COVID-19 patients [36]. This is a good illustration of how a single mutation might be able to redirect a virus to a reserve receptor and assist in the dramatic increase of viral spreading and pathogenicity. A number of other non-canonical receptors were proposed as potential entry points of the virus, such as basigin (CD147) or the tyrosine-protein kinase receptor UFO (AXL) receptor that can mediate the invasion of the SARS-CoV-2 virus into host cells [104,105]. The revealed importance of CD147 in virus infection inspired a search for other potential receptors of the SARS-CoV-2 virus using computational methods [106]. A representative member of the GFR family—fibroblast growth factor receptors (FGFRs)—are TKRs that play an important role in cell proliferation, migration, differentiation, and carcinogenesis [80], and are also found to be relevant in viral infections, for example, during endocytosis of the adeno-associated virus 2, and can be considered as a reserve receptor for SARS-CoV-2.

4.3. Development of Phage-Derived Probes for Monitoring Virus Evolution

Considering the different mechanisms of SARS-CoV-2 infectivity and evolution, including the use of alternative receptors, the knowledge and monitoring of receptor specificity for emerging variants is important in the creation of new antiviral vaccines and medicines. Currently used methods of structural and biochemical analysis, including computer analysis, are not well suited to meet the challenging problem of global screening and monitoring the origins of emerging diseases, such as COVID-19 [106]. Considering the importance of the problem and the lack of simple and reliable approaches to its solution, we thought about phage probes that could compete with pseudoviruses prepared from emerged virus variants in cell culture or animal models and produce a first warning signal regarding a change in virus behavior. In opposition to spike-binding antibodies commonly used to study the interaction of the ACE2 receptor with SARS-CoV-2 variants in vitro, phage spike mimotopes targeted to other cellular receptors would allow the discovery of their role in viral infection in vivo using cell culture, tissue samples, organs, or the whole organism. Phage probes, and specifically landscape phage probes, previously demonstrated several remarkable features, which determined their use in detection of biological threats [44,53,82,107,108,109,110,111,112,113]. The phage mimicry approach, previously suggested [31] for epitope discovery, was essentially modified and adapted for receptor binder discovery.
We evaluated the performance of the phage mimicry approach in the preparation of phage-derived probes targeted to a representative of the GFR family—fibroblast growth factor receptors (FGFR) [80]. Following the phage mimicry strategy, we assumed that if selected cell-associated phage mimotopes are related, for example, to FGF1 or NRP1 (the natural partners of FGFR3 [81,85,86,87,88]), then they can be related to the mimotopes of the spike S protein, which also bind to FGFR3. To realize this idea, we focused on NRP1, which was shown to form a complex with FGFR3. Since the 3D structure of the complex is unknown, we used methods of molecular modelling. In the similar way, we determined a structure of FGFR3 with the SARS-CoV-2 spike RBD and NRP1-FGFR3. We found that the binding AA clusters of NRP1 have a striking similarity to FGF1 and the RBD, as shown in Figure 8 and Table 2. An ELISA was used to confirm our predictions. In this study, a panel of closely related phage mimotopes mimicking the receptor-binding sites on the surface of FGF1, NRP1, and the SARS-CoV-2 spike protein have been identified, and their binding to FGFR3 was shown by an ELISA and confirmed by molecular modelling. Remarkably, it was observed that some of the amino acids of phage mimotopes, and the corresponding spike RBD AA clusters interacting with FGFR3, were absent in the RBD of the SARS-CoV variant [40] (Figure 12). We hypothesize that these residues have been mutated in the emerged variants of the SARS-CoV-2 virus and have granted them a more aggressive phenotype [4]. In relation to these observations, it is interesting to note that phage EDYSELVSQ has been isolated by affinity selection from Calu-3 cells—a representative cell line of non-small cell lung cancer that overexpress human growth factor receptors and other important cellular receptors. Calu-3 cells are sensitive to SARS-CoV-2 infection and, like other human cancer cells, are commonly used for propagation of the SARS-CoV and SARS-CoV-2 viruses in cell culture [30,72,77,114,115,116,117]. These observation are consistent with the hypothesis that the evolution of the virus in the human population can occur not only by escaping immune pressure and increasing the affinity of RBD to ACE2, as commonly assumed [98], but also through exchanging virus receptors and co-receptors during the adaptation to a new host cell [1,10].

5. Conclusions

In this study, a panel of phage mimotopes, closely related to the FGF1, NRP1, and spike receptor-binding AA clusters, were identified. Their binding to ACE2 and FGFR3 was shown by an ELISA and confirmed by molecular modelling. The phage probes that compete with viruses for binding to certain cellular receptors can provide essential information about the mechanisms of virus infectivity. They can also serve as leads in the development of vaccines or drugs, or be used as interfaces in diagnostics. Here, we demonstrated the phage mimicry strategy and developed experimental algorithms that are being used in our ongoing project to evaluate the role of DPP4, NRP1, EGFR, and other human receptors that may be involved in SARS-CoV-2 evolution. We note that this in vitro system is very limited and further studies are required, using more complete pseudoviral and cell receptor expression systems. These additional studies in more advanced systems are beyond the scope of this proof-of-concept project but have led us to identify potential candidate mimotopes for further study.

Author Contributions

Conceptualization, V.A.P.; methodology, J.W.G., M.A.S., L.M.D.P. and V.A.P.; software, J.W.G., L.M.D.P. and M.A.S.; validation, J.W.G., M.A.S., L.M.D.P. and V.A.P.; formal analysis, J.W.G., M.A.S., L.M.D.P. and V.A.P.; investigation, J.W.G. and V.A.P.; resources, V.A.P.; data curation, J.W.G. and V.A.P.; writing—original draft preparation, V.A.P. and J.W.G.; writing—review and editing, V.A.P. and J.W.G.; visualization, J.W.G. and M.A.S.; supervision, V.A.P. and J.W.G.; project administration, V.A.P. and J.W.G.; funding acquisition, V.A.P. All authors have read and agreed to the published version of the manuscript.

Funding

Development of anticancer cell-targeted landscape phages was funded by NIH grants to V.A.P.: 5R01CA125063-01-05, 1R01CA193246-01-05 and NIH (NCI) 1U54CA151881. Project support was provided by the College of Veterinary Medicine (Dean Calvin M. Johnson, DVM, PhD) through the Auburn University Research Initiative in Cancer (AURIC) award to V.A.P. The content is solely the responsibility of the authors and does not necessarily represent the official views of the funding organizations.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or are available from publicly available datasets. Any additional data is available on request from the corresponding author.

Acknowledgments

Anticancer cell-targeted landscape phages used in this project have been developed by students and members of the Petrenko laboratory: Bedi, D., De Plano, L.M., Fagbohun, O.A., Gross, A.L., Gillespie, J.W., Jayanna, P.K., Kuzmicheva, G.A., O’Dell, T., Puzyrev, A.T., Wei, L., Yang, L., and Xu, H.

Conflicts of Interest

The authors declare no competing financial interest.

References

  1. Millet, J.K.; Jaimes, J.A.; Whittaker, G.R. Molecular diversity of coronavirus host cell entry receptors. FEMS Microbiol. Rev. 2021, 45, fuaa057. [Google Scholar] [CrossRef]
  2. Zhou, W.; Wang, W. Fast-spreading SARS-CoV-2 variants: Challenges to and new design strategies of COVID-19 vaccines. Signal Transduct. Target. 2021, 6, 226. [Google Scholar] [CrossRef] [PubMed]
  3. van Oosterhout, C.; Hall, N.; Ly, H.; Tyler, K.M. COVID-19 evolution during the pandemic Implications of new SARS-CoV-2 variants on disease control and public health policies. Virulence 2021, 12, 507–508. [Google Scholar] [CrossRef] [PubMed]
  4. Li, Q.; Wu, J.; Nie, J.; Zhang, L.; Hao, H.; Liu, S.; Zhao, C.; Zhang, Q.; Liu, H.; Nie, L.; et al. The Impact of Mutations in SARS-CoV-2 Spike on Viral Infectivity and Antigenicity. Cell 2020, 182, 1284–1294e1289. [Google Scholar] [CrossRef] [PubMed]
  5. Weisblum, Y.; Schmidt, F.; Zhang, F.; DaSilva, J.; Poston, D.; Lorenzi, J.C.; Muecksch, F.; Rutkowska, M.; Hoffmann, H.H.; Michailidis, E.; et al. Escape from neutralizing antibodies by SARS-CoV-2 spike protein variants. Elife 2020, 9, e61312. [Google Scholar] [CrossRef]
  6. Gan, H.H.; Twaddle, A.; Marchand, B.; Gunsalus, K.C. Structural Modeling of the SARS-CoV-2 Spike/Human ACE2 Complex Interface can Identify High-Affinity Variants Associated with Increased Transmissibility. J. Mol. Biol. 2021, 433, 167051. [Google Scholar] [CrossRef]
  7. Hastie, K.M.; Li, H.; Bedinger, D.; Schendel, S.L.; Dennison, S.M.; Li, K.; Rayaprolu, V.; Yu, X.; Mann, C.; Zandonatti, M.; et al. Defining variant-resistant epitopes targeted by SARS-CoV-2 antibodies: A global consortium study. Science 2021, 374, 472–478. [Google Scholar] [CrossRef]
  8. Ibrahim, I.M.; Elfiky, A.A.; Elgohary, A.M. Recognition through GRP78 is enhanced in the UK, South African, and Brazilian variants of SARS-CoV-2; An in silico perspective. Biochem. Biophys. Res. Commun. 2021, 562, 89–93. [Google Scholar] [CrossRef]
  9. Dos Santos, W.G. Impact of virus genetic variability and host immunity for the success of COVID-19 vaccines. Biomed. Pharm. 2021, 136, 111272. [Google Scholar] [CrossRef]
  10. McCarthy, K.R.; Rennick, L.J.; Nambulli, S.; Robinson-McCarthy, L.R.; Bain, W.G.; Haidar, G.; Duprex, W.P. Recurrent deletions in the SARS-CoV-2 spike glycoprotein drive antibody escape. Science 2021, 371, 1139–1142. [Google Scholar] [CrossRef]
  11. Wang, P.; Casner, R.G.; Nair, M.S.; Wang, M.; Yu, J.; Cerutti, G.; Liu, L.; Kwong, P.D.; Huang, Y.; Shapiro, L.; et al. Increased resistance of SARS-CoV-2 variant P.1 to antibody neutralization. Cell Host. Microbe. 2021, 29, 747–751e744. [Google Scholar] [CrossRef] [PubMed]
  12. Harvey, W.T.; Carabelli, A.M.; Jackson, B.; Gupta, R.K.; Thomson, E.C.; Harrison, E.M.; Ludden, C.; Reeve, R.; Rambaut, A.; Consortium, C.-G.U.; et al. SARS-CoV-2 variants, spike mutations and immune escape. Nat. Rev. Microbiol. 2021, 19, 409–424. [Google Scholar] [CrossRef] [PubMed]
  13. Greaney, A.J.; Starr, T.N.; Gilchuk, P.; Zost, S.J.; Binshtein, E.; Loes, A.N.; Hilton, S.K.; Huddleston, J.; Eguia, R.; Crawford, K.H.D.; et al. Complete Mapping of Mutations to the SARS-CoV-2 Spike Receptor-Binding Domain that Escape Antibody Recognition. Cell Host. Microbe. 2021, 29, 44–57e49. [Google Scholar] [CrossRef] [PubMed]
  14. Stoddard, C.I.; Galloway, J.; Chu, H.Y.; Shipley, M.M.; Sung, K.; Itell, H.L.; Wolf, C.R.; Logue, J.K.; Magedson, A.; Garrett, M.E.; et al. Epitope profiling reveals binding signatures of SARS-CoV-2 immune response in natural infection and cross-reactivity with endemic human CoVs. Cell Rep. 2021, 35, 109164. [Google Scholar] [CrossRef]
  15. Starr, T.N.; Greaney, A.J.; Addetia, A.; Hannon, W.W.; Choudhary, M.C.; Dingens, A.S.; Li, J.Z.; Bloom, J.D. Prospective mapping of viral mutations that escape antibodies used to treat COVID-19. Science 2021, 371, 850–854. [Google Scholar] [CrossRef]
  16. Mascola, J.R.; Graham, B.S.; Fauci, A.S. SARS-CoV-2 Viral Variants-Tackling a Moving Target. JAMA-J. Am. Med. Assoc. 2021, 325, 1261–1262. [Google Scholar] [CrossRef]
  17. Ferraz, M.V.F.; Moreira, E.G.; Coelho, D.F.; Wallau, G.L.; Lins, R.D. Immune evasion of SARS-CoV-2 variants of concern is driven by low affinity to neutralizing antibodies. Chem. Commun. 2021, 57, 6094–6097. [Google Scholar] [CrossRef]
  18. Valdes-Balbin, Y.; Santana-Mederos, D.; Paquet, F.; Fernandez, S.; Climent, Y.; Chiodo, F.; Rodriguez, L.; Ramirez, B.S.; Leon, K.; Hernandez, T.; et al. Molecular Aspects Concerning the Use of the SARS-CoV-2 Receptor Binding Domain as a Target for Preventive Vaccines. ACS Cent. Sci. 2021, 7, 757–767. [Google Scholar] [CrossRef]
  19. Shah, N.N.; Ul Nabi, S.; Rather, M.A.; Kalwar, Q.; Ali, S.I.; Sheikh, W.M.; Ganai, A.; Bashir, S.M. An update on emerging therapeutics to combat COVID-19. Basic Clin. Pharmacol. Toxicol. 2021, 129, 104–129. [Google Scholar] [CrossRef]
  20. Kaur, N.; Singh, R.; Dar, Z.; Bijarnia, R.K.; Dhingra, N.; Kaur, T. Genetic comparison among various coronavirus strains for the identification of potential vaccine targets of SARS-CoV2. Infect. Genet. Evol. 2021, 89, 104490. [Google Scholar] [CrossRef]
  21. Corti, D.; Purcell, L.A.; Snell, G.; Veesler, D. Tackling COVID-19 with neutralizing monoclonal antibodies. Cell 2021, 184, 3086–3108. [Google Scholar] [CrossRef]
  22. Smith, C.C.; Olsen, K.S.; Gentry, K.M.; Sambade, M.; Beck, W.; Garness, J.; Entwistle, S.; Willis, C.; Vensko, S.; Woods, A.; et al. Landscape and selection of vaccine epitopes in SARS-CoV-2. Genome. Med. 2021, 13, 101. [Google Scholar] [CrossRef] [PubMed]
  23. Zhang, Y.N.; Paynter, J.; Sou, C.; Fourfouris, T.; Wang, Y.; Abraham, C.; Ngo, T.; Zhang, Y.; He, L.; Zhu, J. Mechanism of a COVID-19 nanoparticle vaccine candidate that elicits a broadly neutralizing antibody response to SARS-CoV-2 variants. bioRxiv 2021. [Google Scholar] [CrossRef] [PubMed]
  24. Zeyaullah, M.; AlShahrani, A.M.; Muzammil, K.; Ahmad, I.; Alam, S.; Khan, W.H.; Ahmad, R. COVID-19 and SARS-CoV-2 Variants: Current Challenges and Health Concern. Front. Genet. 2021, 12, 693916. [Google Scholar] [CrossRef] [PubMed]
  25. Kreutzberger, A.J.B.; Sanyal, A.; Ojha, R.; Pyle, J.D.; Vapalahti, O.; Balistreri, G.; Kirchhausen, T. Synergistic block of SARS-CoV-2 infection by combined drug inhibition of the host entry factors PIKfyve kinase and TMPRSS2 protease. J. Virol. 2021, 95, e00975-21. [Google Scholar] [CrossRef]
  26. Devarakonda, C.K.V.; Meredith, E.; Ghosh, M.; Shapiro, L.H. Coronavirus Receptors as Immune Modulators. J. Immunol. 2021, 206, 923–929. [Google Scholar] [CrossRef]
  27. Costa, L.B.; Perez, L.G.; Palmeira, V.A.; Cordeiro, T.M.E.; Ribeiro, V.T.; Lanza, K.; Silva, A. Insights on SARS-CoV-2 Molecular Interactions With the Renin-Angiotensin System. Front. Cell. Dev. Biol. 2020, 8, 13. [Google Scholar] [CrossRef]
  28. Li, F. Structure, Function, and Evolution of Coronavirus Spike Proteins. Annu. Rev. Virol. 2016, 3, 237–261. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, N.; Shi, X.; Jiang, L.; Zhang, S.; Wang, D.; Tong, P.; Guo, D.; Fu, L.; Cui, Y.; Liu, X.; et al. Structure of MERS-CoV spike receptor-binding domain complexed with human receptor DPP4. Cell Res. 2013, 23, 986–993. [Google Scholar] [CrossRef] [Green Version]
  30. Wurtz, N.; Penant, G.; Jardot, P.; Duclos, N.; La Scola, B. Culture of SARS-CoV-2 in a panel of laboratory cell lines, permissivity, and differences in growth profile. Eur. J. Clin. Microbiol. Infect. Dis. 2021, 40, 477–484. [Google Scholar] [CrossRef]
  31. Smith, G.P.; Petrenko, V.A. Phage Display. Chem. Rev. 1997, 97, 391–410. [Google Scholar] [CrossRef] [PubMed]
  32. Liu, S.F.; Selvaraj, P.; Lien, C.Z.; Nunez, I.A.; Wu, W.W.; Chou, C.K.; Wang, T.T. The PRRA Insert at the S1/S2 Site Modulates Cellular Tropism of SARS-CoV-2 and ACE2 Usage by the Closely Related Bat RaTG13. J. Virol. 2021, 95, e01751-20. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, J.R.; Yang, Y.Q.; Hou, Z.Y.; Jia, H.; Wang, Y.J. TRPV2-spike protein interaction mediates the entry of SARS-CoV-2 into macrophages in febrile conditions. Theranostics 2021, 11, 7379–7390. [Google Scholar] [CrossRef] [PubMed]
  34. Pieren, M.; Prota, A.E.; Ruch, C.; Kostrewa, D.; Wagner, A.; Biedermann, K.; Winkler, F.K.; Ballmer-Hofer, K. Crystal structure of the Orf virus NZ2 variant of vascular endothelial growth factor-E. Implications for receptor specificity. J. Biol. Chem. 2006, 281, 19578–19587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Hondermarck, H.; Bartlett, N.W.; Nurcombe, V. The role of growth factor receptors in viral infections: An opportunity for drug repurposing against emerging viral diseases such as COVID-19? FASEB Bioadv. 2020, 2, 296–303. [Google Scholar] [CrossRef]
  36. Carvacho, I.; Piesche, M. RGD-binding integrins and TGF-beta in SARS-CoV-2 infections novel targets to treat COVID-19 patients? Clin. Transl. Immunol. 2021, 10, e1240. [Google Scholar] [CrossRef]
  37. Masre, S.F.; Jufri, N.F.; Ibrahim, F.W.; Abdul Raub, S.H. Classical and alternative receptors for SARS-CoV-2 therapeutic strategy. Rev. Med. Virol. 2021, 31, 1–9. [Google Scholar] [CrossRef]
  38. Watanabe, Y.; Allen, J.D.; Wrapp, D.; McLellan, J.S.; Crispin, M. Site-specific glycan analysis of the SARS-CoV-2 spike. Science 2020, 369, 330–333. [Google Scholar] [CrossRef]
  39. Woo, H.; Park, S.J.; Choi, Y.K.; Park, T.; Tanveer, M.; Cao, Y.; Kern, N.R.; Lee, J.; Yeom, M.S.; Croll, T.I.; et al. Developing a Fully Glycosylated Full-Length SARS-CoV-2 Spike Protein Model in a Viral Membrane. J. Phys. Chem. B 2020, 124, 7128–7137. [Google Scholar] [CrossRef]
  40. Zhou, P.; Yang, X.L.; Wang, X.G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.R.; Zhu, Y.; Li, B.; Huang, C.L.; et al. Addendum: A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 588, E6. [Google Scholar] [CrossRef]
  41. Hoffmann, M.; Kleine-Weber, H.; Pohlmann, S. A Multibasic Cleavage Site in the Spike Protein of SARS-CoV-2 Is Essential for Infection of Human Lung Cells. Mol. Cell 2020, 78, 779–784e775. [Google Scholar] [CrossRef] [PubMed]
  42. Jackson, C.B.; Farzan, M.; Chen, B.; Choe, H. Mechanisms of SARS-CoV-2 entry into cells. Nat. Rev. Mol. Cell Biol. 2021, 23, 3–20. [Google Scholar] [CrossRef] [PubMed]
  43. Petrenko, V.A.; Smith, G.P. Chapter 2: Vectors and modes of display. In Chapter 2. Phage Display in Biotechnology and Drug Discovery; Geyer, S.S.S.C.R., Ed.; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2005; pp. 43–74. [Google Scholar]
  44. Petrenko, V.A. Landscape Phage: Evolution from Phage Display to Nanobiotechnology. Viruses 2018, 10, 311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kuzmicheva, G.A.; Jayanna, P.K.; Eroshkin, A.M.; Grishina, M.A.; Pereyaslavskaya, E.S.; Potemkin, V.A.; Petrenko, V.A. Mutations in fd phage major coat protein modulate affinity of the displayed peptide. Protein. Eng. Des. Sel. 2009, 22, 631–639. [Google Scholar] [CrossRef] [Green Version]
  46. Petrenko, V.A.; Smith, G.P.; Gong, X.; Quinn, T. A library of organic landscapes on filamentous phage. Protein. Eng. 1996, 9, 797–801. [Google Scholar] [CrossRef]
  47. Thaper, D.; Prabha, V. Molecular mimicry: An explanation for autoimmune diseases and infertility. Scand. J. Immunol. 2018, 88, e12697. [Google Scholar] [CrossRef]
  48. Wildner, G.; Diedrichs-Mohring, M. Molecular Mimicry and Uveitis. Front. Immunol. 2020, 11, 580636. [Google Scholar] [CrossRef]
  49. Mayrose, I.; Penn, O.; Erez, E.; Rubinstein, N.D.; Shlomi, T.; Freund, N.T.; Bublil, E.M.; Ruppin, E.; Sharan, R.; Gershoni, J.M.; et al. Pepitope: Epitope mapping from affinity-selected peptides. Bioinformatics 2007, 23, 3244–3246. [Google Scholar] [CrossRef] [Green Version]
  50. Huang, J.; Ru, B.; Dai, P. Bioinformatics Resources and Tools for Phage Display. Molecules 2011, 16, 694–709. [Google Scholar] [CrossRef] [Green Version]
  51. Fan, W.; Mencius, J.; Du, W.; Fan, H.; Zhu, H.; Wei, D.; Zhou, M.; Quan, S. Online bioinformatics teaching practice: Comparison of popular docking programs using SARS-CoV-2 spike RBD-ACE2 complex as a benchmark. Biochem. Mol. Biol. Educ. 2021, 49, 833–840. [Google Scholar] [CrossRef]
  52. Huang, J.; Derda, R.; Huang, Y. Phage display informatics. Comput. Math. Methods Med. 2013, 2013, 698395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Jayanna, P.K.; Bedi, D.; Deinnocentes, P.; Bird, R.C.; Petrenko, V.A. Landscape phage ligands for PC3 prostate carcinoma cells. Protein. Eng. Des. Sel. 2010, 23, 423–430. [Google Scholar] [CrossRef] [PubMed]
  54. Andrieu, J.; Re, F.; Russo, L.; Nicotra, F. Phage-displayed peptides targeting specific tissues and organs. J. Drug Target. 2019, 27, 555–565. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Gillespie, J.W.; Wei, L.; Petrenko, V.A. Selection of Lung Cancer-Specific Landscape Phage for Targeted Drug Delivery. Comb. Chem. High Throughput. Screen 2016, 19, 412–422. [Google Scholar] [CrossRef] [PubMed]
  56. Bedi, D.; Gillespie, J.W.; Petrenko, V.A. Selection of pancreatic cancer cell-binding landscape phages and their use in development of anticancer nanomedicines. Protein. Eng. Des. Sel. 2014, 27, 235–243. [Google Scholar] [CrossRef]
  57. Fagbohun, O.A.; Bedi, D.; Grabchenko, N.I.; Deinnocentes, P.A.; Bird, R.C.; Petrenko, V.A. Landscape phages and their fusion proteins targeted to breast cancer cells. Protein. Eng. Des. Sel. 2012, 25, 271–283. [Google Scholar] [CrossRef] [Green Version]
  58. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; de Beer, T.A.P.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic. Acids Res. 2018, 46, W296–W303. [Google Scholar] [CrossRef] [Green Version]
  59. Bienert, S.; Waterhouse, A.; de Beer, T.A.; Tauriello, G.; Studer, G.; Bordoli, L.; Schwede, T. The SWISS-MODEL Repository-new features and functionality. Nucleic. Acids Res. 2017, 45, D313–D319. [Google Scholar] [CrossRef] [Green Version]
  60. Studer, G.; Tauriello, G.; Bienert, S.; Biasini, M.; Johner, N.; Schwede, T. ProMod3-A versatile homology modelling toolbox. PLoS Comput. Biol. 2021, 17, e1008667. [Google Scholar] [CrossRef]
  61. Studer, G.; Rempfer, C.; Waterhouse, A.M.; Gumienny, R.; Haas, J.; Schwede, T. QMEANDisCo-distance constraints applied on model quality estimation. Bioinformatics 2020, 36, 2647. [Google Scholar] [CrossRef] [Green Version]
  62. Bertoni, M.; Kiefer, F.; Biasini, M.; Bordoli, L.; Schwede, T. Modeling protein quaternary structure of homo- and hetero-oligomers beyond binary interactions by homology. Sci. Rep. 2017, 7, 10480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Steinegger, M.; Meier, M.; Mirdita, M.; Vohringer, H.; Haunsberger, S.J.; Soding, J. HH-suite3 for fast remote homology detection and deep protein annotation. BMC Bioinform. 2019, 20, 473. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Krieger, E.; Vriend, G. YASARA View molecular graphics for all devices from smartphones to workstations. Bioinformatics 2014, 30, 2981–2982. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Krieger, E.; Vriend, G. New ways to boost molecular dynamics simulations. J. Comput. Chem. 2015, 36, 996–1007. [Google Scholar] [CrossRef]
  67. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [Green Version]
  68. Honorato, R.V.; Koukos, P.I.; Jimenez-Garcia, B.; Tsaregorodtsev, A.; Verlato, M.; Giachetti, A.; Rosato, A.; Bonvin, A. Structural Biology in the Clouds: The WeNMR-EOSC Ecosystem. Front. Mol. Biosci. 2021, 8, 729513. [Google Scholar] [CrossRef]
  69. van Zundert, G.C.P.; Rodrigues, J.; Trellet, M.; Schmitz, C.; Kastritis, P.L.; Karaca, E.; Melquiond, A.S.J.; van Dijk, M.; de Vries, S.J.; Bonvin, A. The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes. J. Mol. Biol. 2016, 428, 720–725. [Google Scholar] [CrossRef] [Green Version]
  70. Smith, G.P.; Petrenko, V.A.; Matthews, L.J. Cross-linked filamentous phage as an affinity matrix. J. Immunol. Methods 1998, 215, 151–161. [Google Scholar] [CrossRef]
  71. Gillespie, J.W.; Gross, A.L.; Puzyrev, A.T.; Bedi, D.; Petrenko, V.A. Combinatorial synthesis and screening of cancer cell-specific nanomedicines targeted via phage fusion proteins. Front. Microbiol. 2015, 6, 628. [Google Scholar] [CrossRef]
  72. Xia, S.; Lan, Q.; Zhu, Y.; Wang, C.; Xu, W.; Li, Y.; Wang, L.; Jiao, F.; Zhou, J.; Hua, C.; et al. Structural and functional basis for pan-CoV fusion inhibitors against SARS-CoV-2 and its variants with preclinical evaluation. Signal Transduct. Target. 2021, 6, 288. [Google Scholar] [CrossRef] [PubMed]
  73. Gross, A.L.; Gillespie, J.W.; Petrenko, V.A. Promiscuous tumor targeting phage proteins. Protein. Eng. Des. Sel. 2016, 29, 93–103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Drosten, C.; Gunther, S.; Preiser, W.; van der Werf, S.; Brodt, H.R.; Becker, S.; Rabenau, H.; Panning, M.; Kolesnikova, L.; Fouchier, R.A.; et al. Identification of a novel coronavirus in patients with severe acute respiratory syndrome. N. Engl. J. Med. 2003, 348, 1967–1976. [Google Scholar] [CrossRef] [PubMed]
  75. Wu, S.; Du, Z.; Shen, S.; Zhang, B.; Yang, H.; Li, X.; Cui, W.; Cheng, F.; Huang, J. Identification and Validation of a Novel Clinical Signature to Predict the Prognosis in Confirmed Coronavirus Disease 2019 Patients. Clin. Infect. Dis. 2020, 71, 3154–3162. [Google Scholar] [CrossRef]
  76. Zhou, B.P.; Chen, X.C.; Wang, H.S.; Li, M.Z.; Hu, Y.W.; Du, F.; Xu, L.M.; Yang, G.L. Identification and molecular cloning and sequence analysis of a novel coronavirus from patients with SARS by RT-PCR. Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi 2003, 17, 137–139. [Google Scholar]
  77. Kumar, S.; Sarma, P.; Kaur, H.; Prajapat, M.; Bhattacharyya, A.; Avti, P.; Sehkhar, N.; Kaur, H.; Bansal, S.; Mahendiratta, S.; et al. Clinically relevant cell culture models and their significance in isolation, pathogenesis, vaccine development, repurposing and screening of new drugs for SARS-CoV-2: A systematic review. Tissue Cell 2021, 70, 101497. [Google Scholar] [CrossRef]
  78. Khan, I.; Hatiboglu, M.A. Can COVID-19 induce glioma tumorogenesis through binding cell receptors? Med. Hypotheses 2020, 144, 110009. [Google Scholar] [CrossRef]
  79. Shaik, F.; Cuthbert, G.A.; Homer-Vanniasinkam, S.; Muench, S.P.; Ponnambalam, S.; Harrison, M.A. Structural Basis for Vascular Endothelial Growth Factor Receptor Activation and Implications for Disease Therapy. Biomolecules 2020, 10, 25. [Google Scholar] [CrossRef]
  80. Purcaru, O.-S.; Artene, S.-A.; Barcan, E.; Silosi, C.A.; Stanciu, I.; Danoiu, S.; Tudorache, S.; Tataranu, L.G.; Dricu, A. The Interference between SARS-CoV-2 and Tyrosine Kinase Receptor Signaling in Cancer. Int. J. Mol. Sci. 2021, 22, 4830. [Google Scholar] [CrossRef]
  81. Ferguson, H.R.; Smith, M.P.; Francavilla, C. Fibroblast Growth Factor Receptors (FGFRs) and Noncanonical Partners in Cancer Signaling. Cells 2021, 10, 1201. [Google Scholar] [CrossRef]
  82. Petrenko, V.A. Evolution of phage display: From bioactive peptides to bioselective nanomaterials. Expert. Opin. Drug Deliv. 2008, 5, 825–836. [Google Scholar] [CrossRef]
  83. Walls, A.C.; Park, Y.J.; Tortorici, M.A.; Wall, A.; McGuire, A.T.; Veesler, D. Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell 2020, 181, 281–292e286. [Google Scholar] [CrossRef]
  84. Lan, J.; Ge, J.; Yu, J.; Shan, S.; Zhou, H.; Fan, S.; Zhang, Q.; Shi, X.; Wang, Q.; Zhang, L.; et al. Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor. Nature 2020, 581, 215–220. [Google Scholar] [CrossRef] [Green Version]
  85. Abdullah, A.; Akhand, S.S.; Paez, J.S.P.; Brown, W.; Pan, L.; Libring, S.; Badamy, M.; Dykuizen, E.; Solorio, L.; Andy Tao, W.; et al. Epigenetic targeting of neuropilin-1 prevents bypass signaling in drug-resistant breast cancer. Oncogene 2021, 40, 322–333. [Google Scholar] [CrossRef]
  86. Plotnikov, A.N.; Hubbard, S.R.; Schlessinger, J.; Mohammadi, M. Crystal structures of two FGF-FGFR complexes reveal the determinants of ligand-receptor specificity. Cell 2000, 101, 413–424. [Google Scholar] [CrossRef] [Green Version]
  87. Giacomini, A.; Chiodelli, P.; Matarazzo, S.; Rusnati, M.; Presta, M.; Ronca, R. Blocking the FGF/FGFR system as a “two-compartment” antiangiogenic/antitumor approach in cancer therapy. Pharm. Res. 2016, 107, 172–185. [Google Scholar] [CrossRef]
  88. Porta, R.; Borea, R.; Coelho, A.; Khan, S.; Araujo, A.; Reclusa, P.; Franchina, T.; Van Der Steen, N.; Van Dam, P.; Ferri, J.; et al. FGFR a promising druggable target in cancer: Molecular biology and new drugs. Crit. Rev. Oncol. Hematol. 2017, 113, 256–267. [Google Scholar] [CrossRef] [Green Version]
  89. Li, Y.; Zhang, Z.; Yang, L.; Lian, X.; Xie, Y.; Li, S.; Xin, S.; Cao, P.; Lu, J. The MERS-CoV Receptor DPP4 as a Candidate Binding Target of the SARS-CoV-2 Spike. iScience 2020, 23, 101400. [Google Scholar] [CrossRef]
  90. Kyrou, I.; Randeva, H.S.; Spandidos, D.A.; Karteris, E. Not only ACE2-the quest for additional host cell mediators of SARS-CoV-2 infection: Neuropilin-1 (NRP1) as a novel SARS-CoV-2 host cell entry mediator implicated in COVID-19. Signal Transduct. Target. 2021, 6, 21. [Google Scholar] [CrossRef]
  91. Cantuti-Castelvetri, L.; Ojha, R.; Pedro, L.D.; Djannatian, M.; Franz, J.; Kuivanen, S.; van der Meer, F.; Kallio, K.; Kaya, T.; Anastasina, M.; et al. Neuropilin-1 facilitates SARS-CoV-2 cell entry and infectivity. Science 2020, 370, 856–860. [Google Scholar] [CrossRef]
  92. Gudowska-Sawczuk, M.; Mroczko, B. The Role of Neuropilin-1 (NRP-1) in SARS-CoV-2 Infection: Review. J. Clin. Med. 2021, 10, 2772. [Google Scholar] [CrossRef]
  93. Mayi, B.S.; Leibowitz, J.A.; Woods, A.T.; Ammon, K.A.; Liu, A.E.; Raja, A. The role of Neuropilin-1 in COVID-19. PLoS Pathog. 2021, 17, e1009153. [Google Scholar] [CrossRef]
  94. Alnomasy, S.F. Virus-receptor interactions of SARS-CoV-2 spikereceptor-binding domain and human neuropilin-1 b1 domain. Saudi. J. Biol. Sci. 2021, 28, 3926–3928. [Google Scholar] [CrossRef]
  95. Li, Z.L.; Buck, M. Neuropilin-1 Assists SARS-CoV-2 Infection by Stimulating the Separation of Spike Protein Domains S1 and S2. bioRxiv 2021, 120, 2828–2837. [Google Scholar] [CrossRef]
  96. Daly, J.L.; Simonetti, B.; Klein, K.; Chen, K.E.; Williamson, M.K.; Anton-Plagaro, C.; Shoemark, D.K.; Simon-Gracia, L.; Bauer, M.; Hollandi, R.; et al. Neuropilin-1 is a host factor for SARS-CoV-2 infection. Science 2020, 370, 861. [Google Scholar] [CrossRef]
  97. Ghorbani, M.; Brooks, B.R.; Klauda, J.B. Critical Sequence Hotspots for Binding of Novel Coronavirus to Angiotensin Converter Enzyme as Evaluated by Molecular Simulations. J. Phys. Chem. B 2020, 124, 10034–10047. [Google Scholar] [CrossRef]
  98. Chen, C.; Boorla, V.S.; Banerjee, D.; Chowdhury, R.; Cavener, V.S.; Nissly, R.H.; Gontu, A.; Boyle, N.R.; Vandegrift, K.; Nair, M.S.; et al. Computational prediction of the effect of amino acid changes on the binding affinity between SARS-CoV-2 spike RBD and human ACE2. Proc. Natl. Acad. Sci. USA 2021, 118, 42. [Google Scholar] [CrossRef]
  99. Cueno, M.E.; Imai, K. Structural Comparison of the SARS CoV 2 Spike Protein Relative to Other Human-Infecting Coronaviruses. Front. Med. (Lausanne) 2020, 7, 594439. [Google Scholar] [CrossRef]
  100. Xu, C.; Wang, A.; Geng, K.; Honnen, W.; Wang, X.; Bruiners, N.; Singh, S.; Ferrara, F.; D’Angelo, S.; Bradbury, A.R.M.; et al. Human Immunodeficiency Viruses Pseudotyped with SARS-CoV-2 Spike Proteins Infect a Broad Spectrum of Human Cell Lines through Multiple Entry Mechanisms. Viruses 2021, 13, 953. [Google Scholar] [CrossRef]
  101. Li, F. Receptor Recognition Mechanisms of Coronaviruses: A Decade of Structural Studies. J. Virol. 2015, 89, 1954–1964. [Google Scholar] [CrossRef] [Green Version]
  102. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Kruger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 271–280e278. [Google Scholar] [CrossRef]
  103. Hussein, H.A.; Walker, L.R.; Abdel-Raouf, U.M.; Desouky, S.A.; Montasser, A.K.; Akula, S.M. Beyond RGD: Virus interactions with integrins. Arch. Virol. 2015, 160, 2669–2681. [Google Scholar] [CrossRef]
  104. Wang, K.; Chen, W.; Zhang, Z.; Deng, Y.; Lian, J.Q.; Du, P.; Wei, D.; Zhang, Y.; Sun, X.X.; Gong, L.; et al. CD147-spike protein is a novel route for SARS-CoV-2 infection to host cells. Signal Transduct. Target. 2020, 5, 283. [Google Scholar] [CrossRef]
  105. Chen, Z.; Mi, L.; Xu, J.; Yu, J.; Wang, X.; Jiang, J.; Xing, J.; Shang, P.; Qian, A.; Li, Y.; et al. Function of HAb18G/CD147 in invasion of host cells by severe acute respiratory syndrome coronavirus. J. Infect. Dis. 2005, 191, 755–760. [Google Scholar] [CrossRef] [Green Version]
  106. Mobini, S.; Chizari, M.; Mafakher, L.; Rismani, E.; Rismani, E. Structure-based study of immune receptors as eligible binding targets of coronavirus SARS-CoV-2 spike protein. J. Mol. Graph. Model. 2021, 108, 107997. [Google Scholar] [CrossRef]
  107. Petrenko, V.A.; Vodyanoy, V.J. Phage display for detection of biological threat agents. J. Microbiol. Methods 2003, 53, 253–262. [Google Scholar] [CrossRef]
  108. Brigati, J.; Williams, D.D.; Sorokulova, I.B.; Nanduri, V.; Chen, I.H.; Turnbough, C.L., Jr.; Petrenko, V.A. Diagnostic probes for Bacillus anthracis spores selected from a landscape phage library. Clin. Chem. 2004, 50, 1899–1906. [Google Scholar] [CrossRef]
  109. Petrenko, V.A.; Sorokulova, I.B. Detection of biological threats. A challenge for directed molecular evolution. J. Microbiol. Methods 2004, 58, 147–168. [Google Scholar] [CrossRef]
  110. Sorokulova, I.B.; Olsen, E.V.; Chen, I.H.; Fiebor, B.; Barbaree, J.M.; Vodyanoy, V.J.; Chin, B.A.; Petrenko, V.A. Landscape phage probes for Salmonella typhimurium. J. Microbiol. Methods 2005, 63, 55–72. [Google Scholar] [CrossRef]
  111. Brigati, J.R.; Samoylova, T.I.; Jayanna, P.K.; Petrenko, V.A. Phage display for generating peptide reagents. Curr. Protoc. Protein. Sci. 2008, 18, 19. [Google Scholar] [CrossRef]
  112. Han, L.; Xia, H.Q.; Yin, L.; Petrenko, V.A.; Liu, A.H. Selected landscape phage probe as selective recognition interface for sensitive total prostate-specific antigen immunosensor. Biosens. Bioelectron. 2018, 106, 1–6. [Google Scholar] [CrossRef]
  113. Liu, P.; Wang, Y.B.; Han, L.; Cai, Y.Y.; Ren, H.; Ma, T.X.; Li, X.Q.; Petrenko, V.A.; Liu, A.H. Colorimetric Assay of Bacterial Pathogens Based on Co3O4 Magnetic Nanozymes Conjugated with Specific Fusion Phage Proteins and Magnetophoretic Chromatography. ACS Appl. Mater. Interfaces 2020, 12, 9090–9097. [Google Scholar] [CrossRef]
  114. Chu, H.; Chan, J.F.; Yuen, T.T.; Shuai, H.; Yuan, S.; Wang, Y.; Hu, B.; Yip, C.C.; Tsang, J.O.; Huang, X.; et al. Comparative tropism, replication kinetics, and cell damage profiling of SARS-CoV-2 and SARS-CoV with implications for clinical manifestations, transmissibility, and laboratory studies of COVID-19: An observational study. Lancet. Microbe. 2020, 1, e14–e23. [Google Scholar] [CrossRef]
  115. Zhou, Z.; Liu, Z.; Ou, Q.; Wu, X.; Wang, X.; Shao, Y.; Liu, H.; Yang, Y. Targeting FGFR in non-small cell lung cancer: Implications from the landscape of clinically actionable aberrations of FGFR kinases. Cancer Biol. Med. 2021, 18, 490. [Google Scholar] [CrossRef]
  116. Zuo, W.; He, Y.; Li, W.; Wu, H.; Zhao, Z.; Zhang, Y.; Chen, S.; Yin, Y. Landscape of FGF/FGFR Alterations in 12,372 Chinese Cancer Patients. J. Cancer 2020, 11, 6695–6699. [Google Scholar] [CrossRef]
  117. Helsten, T.; Elkin, S.; Arthur, E.; Tomson, B.N.; Carter, J.; Kurzrock, R. The FGFR Landscape in Cancer: Analysis of 4,853 Tumors by Next-Generation Sequencing. Clin. Cancer Res. 2016, 22, 259–267. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Three-dimensional predicted model of the spike (S) glycoprotein of the SARS-CoV-2 virus composed of (A) two well-defined structural domains (S1 and S2), decorated with (B) 22 N-glycan residues as modeled using 6VSB and 6VXX [38,39]. Monomers of the S protein, composed of polypeptide chains of 1273 amino acids, form homotrimer spikes on the virus surface [18]. Spike protein monomers are composed of three major structural domains: head, stalk, and cytoplasmic tail. The head comprises the N-terminal domain (NTD; yellow) and the receptor-binding domain (RBD; orange), which displays the receptor-binding motif (RBM; cyan) that is responsible for interaction with cell receptors [40]. RBDs in non-activated viral S glycoprotein trimers are present in a hidden “down” conformation. The S glycoprotein is cleaved by host proteases (trypsin and furin) at the site between the S1 and the S2 subunits [41,42]. The S2 domain of the S protein consists of fusion peptide (FP; purple), two heptad-repeat domains (HR1 and HR2; red), a transmembrane domain (TM; gray), and a cytoplasm domain (CP; pink). A second proteolytic site (S2′ site), located within the S2 subdomain, is also cut by type II transmembrane serine protease (TMPRSS2) as well as cathepsin B and L (CatB/L) to enable virus-cell fusion by triggering the dissociation of S1 and the irreversible refolding of S2, a conformational change of the S protein and the fusion of the viral envelope and endosomes [42].
Figure 1. Three-dimensional predicted model of the spike (S) glycoprotein of the SARS-CoV-2 virus composed of (A) two well-defined structural domains (S1 and S2), decorated with (B) 22 N-glycan residues as modeled using 6VSB and 6VXX [38,39]. Monomers of the S protein, composed of polypeptide chains of 1273 amino acids, form homotrimer spikes on the virus surface [18]. Spike protein monomers are composed of three major structural domains: head, stalk, and cytoplasmic tail. The head comprises the N-terminal domain (NTD; yellow) and the receptor-binding domain (RBD; orange), which displays the receptor-binding motif (RBM; cyan) that is responsible for interaction with cell receptors [40]. RBDs in non-activated viral S glycoprotein trimers are present in a hidden “down” conformation. The S glycoprotein is cleaved by host proteases (trypsin and furin) at the site between the S1 and the S2 subunits [41,42]. The S2 domain of the S protein consists of fusion peptide (FP; purple), two heptad-repeat domains (HR1 and HR2; red), a transmembrane domain (TM; gray), and a cytoplasm domain (CP; pink). A second proteolytic site (S2′ site), located within the S2 subdomain, is also cut by type II transmembrane serine protease (TMPRSS2) as well as cathepsin B and L (CatB/L) to enable virus-cell fusion by triggering the dissociation of S1 and the irreversible refolding of S2, a conformational change of the S protein and the fusion of the viral envelope and endosomes [42].
Viruses 14 00384 g001
Figure 2. Electron microscopy image of filamentous phage (left) and electron density model (center) of filamentous phage M13 (courtesy of Lee Makowski and Gregory Kishchenko). Blue and red arrows depict the sharp and blunt ends of the phage capsid with attached minor coat proteins pIII/pIV and pVII/pIX, respectively (five copies each). Major coat protein (~2700 copies) forms the tubular capsid around viral single-stranded DNA (scale bar: 100 nm). 3D structure (right) of the complex between phage displaying the peptide EDYSELVSQ (green) with FGFR3 (red). Here, the peptide-displayed phages are designated by the structure of the inserted foreign peptides.
Figure 2. Electron microscopy image of filamentous phage (left) and electron density model (center) of filamentous phage M13 (courtesy of Lee Makowski and Gregory Kishchenko). Blue and red arrows depict the sharp and blunt ends of the phage capsid with attached minor coat proteins pIII/pIV and pVII/pIX, respectively (five copies each). Major coat protein (~2700 copies) forms the tubular capsid around viral single-stranded DNA (scale bar: 100 nm). 3D structure (right) of the complex between phage displaying the peptide EDYSELVSQ (green) with FGFR3 (red). Here, the peptide-displayed phages are designated by the structure of the inserted foreign peptides.
Viruses 14 00384 g002
Figure 3. (A) Schematic of viral evolution (adaptation) leading to increased fitness to cellular receptors through diversification of viral functional domains; (B) schematic of SARS-CoV-2 functional domains and epitopes (or their mimetics) fused to p8 proteins and selected from landscape phage libraries through phage mimicry; (C) the phage-displayed peptide (mimotope) contains the same or similar amino acid (AA) residues as amino acid clusters (AA clusters) on the surface of spike protein, and presumably can interact with viral cellular receptors.
Figure 3. (A) Schematic of viral evolution (adaptation) leading to increased fitness to cellular receptors through diversification of viral functional domains; (B) schematic of SARS-CoV-2 functional domains and epitopes (or their mimetics) fused to p8 proteins and selected from landscape phage libraries through phage mimicry; (C) the phage-displayed peptide (mimotope) contains the same or similar amino acid (AA) residues as amino acid clusters (AA clusters) on the surface of spike protein, and presumably can interact with viral cellular receptors.
Viruses 14 00384 g003
Figure 4. Clusters of amino acids identified by PepSurf surface accessible alignments of cell receptor-binding phage mimotopes on the SARS-CoV-2 spike protein using the PDB model 6VYB. Cluster 1 (dark purple), cluster 2 (green), cluster 3 (bright pink), cluster 4 (turquoise), cluster 5 (pink), cluster 6 (dark blue), cluster 7 (purple), cluster 8 (blue), cluster 9 (dark pink), cluster 10 (light blue), cluster 11 (yellow), cluster 12 (orange), and cluster 13 (red).
Figure 4. Clusters of amino acids identified by PepSurf surface accessible alignments of cell receptor-binding phage mimotopes on the SARS-CoV-2 spike protein using the PDB model 6VYB. Cluster 1 (dark purple), cluster 2 (green), cluster 3 (bright pink), cluster 4 (turquoise), cluster 5 (pink), cluster 6 (dark blue), cluster 7 (purple), cluster 8 (blue), cluster 9 (dark pink), cluster 10 (light blue), cluster 11 (yellow), cluster 12 (orange), and cluster 13 (red).
Viruses 14 00384 g004
Figure 5. Amino acid residues corresponding to each identified functional cluster.
Figure 5. Amino acid residues corresponding to each identified functional cluster.
Viruses 14 00384 g005
Figure 6. Interaction of ACE2 with different ligands. (A) Interaction of a landscape phage displaying the peptide DGRADLSYD on the full-length p8 protein (yellow) with ACE2 (blue) as determined using homology modeling. Here, a segment containing less than 1% of the landscape phage is presented, where the DGRADLSYD peptide is presented as an N-terminal fusion to all copies of the mature p8 major coat protein. Molecular model 6M0J demonstrating the interaction between ACE2 protein (blue) and recombinant SARS-CoV-2 spike RBD (pink) with amino acid clusters corresponding to phage mimotopes. (B) DGRADLSYD; (C) VGIDEQRAD; and (D) DGRSIVGDE, highlighted in red.
Figure 6. Interaction of ACE2 with different ligands. (A) Interaction of a landscape phage displaying the peptide DGRADLSYD on the full-length p8 protein (yellow) with ACE2 (blue) as determined using homology modeling. Here, a segment containing less than 1% of the landscape phage is presented, where the DGRADLSYD peptide is presented as an N-terminal fusion to all copies of the mature p8 major coat protein. Molecular model 6M0J demonstrating the interaction between ACE2 protein (blue) and recombinant SARS-CoV-2 spike RBD (pink) with amino acid clusters corresponding to phage mimotopes. (B) DGRADLSYD; (C) VGIDEQRAD; and (D) DGRSIVGDE, highlighted in red.
Viruses 14 00384 g006
Figure 7. ACE2-binding phages displaying SARS-CoV-2 spike RBD-mimotopes as characterized by an indirect ELISA. Serial dilutions of phages were incubated with a bound ACE2 receptor, followed by incubation with a rabbit anti-phage IgG and an HRP-conjugated goat anti-rabbit IgG. Signal was produced using a TMB substrate and endpoint absorbance at 450 nm was measured after a 30-min incubation. Candidate RBD-binding phages were compared to the wildtype fd-tet phage (dark green).
Figure 7. ACE2-binding phages displaying SARS-CoV-2 spike RBD-mimotopes as characterized by an indirect ELISA. Serial dilutions of phages were incubated with a bound ACE2 receptor, followed by incubation with a rabbit anti-phage IgG and an HRP-conjugated goat anti-rabbit IgG. Signal was produced using a TMB substrate and endpoint absorbance at 450 nm was measured after a 30-min incubation. Candidate RBD-binding phages were compared to the wildtype fd-tet phage (dark green).
Viruses 14 00384 g007
Figure 8. Interaction of fibroblast growth factor receptor 3 (FGFR3) with different ligands: (A) Interaction of a landscape phage displaying the peptide EDYSELVSQ (yellow) on the full-length p8 protein with FGFR3 (cyan), as determined using homology modeling. Here, a segment containing less than 1% of the landscape phage is presented, where the EDYSELVSQ peptide is presented as an N-terminal fusion to all copies of the mature p8 major coat protein. Interaction of FGFR3 (cyan) with (B) FGF1 (gray), (C) NRP1 (magenta), or (D) SARS-CoV-2 spike RBD (blue), with amino acid clusters containing alignments to the EDYSELVSQ mimotope highlighted in red.
Figure 8. Interaction of fibroblast growth factor receptor 3 (FGFR3) with different ligands: (A) Interaction of a landscape phage displaying the peptide EDYSELVSQ (yellow) on the full-length p8 protein with FGFR3 (cyan), as determined using homology modeling. Here, a segment containing less than 1% of the landscape phage is presented, where the EDYSELVSQ peptide is presented as an N-terminal fusion to all copies of the mature p8 major coat protein. Interaction of FGFR3 (cyan) with (B) FGF1 (gray), (C) NRP1 (magenta), or (D) SARS-CoV-2 spike RBD (blue), with amino acid clusters containing alignments to the EDYSELVSQ mimotope highlighted in red.
Viruses 14 00384 g008
Figure 9. Primary structure of the FGFR3 domain. Amino acids involved in the interaction (designated by capital bold letters) between FGF1 (*), NRP1 (^), and RBD (#) were identified using the YASARA Structure and literature data. Promiscuous amino acid residues of FGFR3 involved with the interaction of FGF1, NRP1, and spike RBD are marked with (*^#). Members of AA clusters corresponding to the phage mimotope EDYSELVSQ are marked with (+) and indicated by underlined letters.
Figure 9. Primary structure of the FGFR3 domain. Amino acids involved in the interaction (designated by capital bold letters) between FGF1 (*), NRP1 (^), and RBD (#) were identified using the YASARA Structure and literature data. Promiscuous amino acid residues of FGFR3 involved with the interaction of FGF1, NRP1, and spike RBD are marked with (*^#). Members of AA clusters corresponding to the phage mimotope EDYSELVSQ are marked with (+) and indicated by underlined letters.
Viruses 14 00384 g009
Figure 10. FGFR3-binding phages displaying FGF1 mimotopes as characterized by an indirect ELISA. Serial dilutions of phages were incubated with the bound FGFR3 receptor, followed by incubation with a rabbit anti-phage IgG and an HRP-conjugated goat anti-rabbit IgG. Signal was produced using a TMB substrate and endpoint absorbance at 450 nm was measured after a 30-min incubation. Candidate FGFR3-binding phages were compared to the wildtype fd-tet phage (dark green).
Figure 10. FGFR3-binding phages displaying FGF1 mimotopes as characterized by an indirect ELISA. Serial dilutions of phages were incubated with the bound FGFR3 receptor, followed by incubation with a rabbit anti-phage IgG and an HRP-conjugated goat anti-rabbit IgG. Signal was produced using a TMB substrate and endpoint absorbance at 450 nm was measured after a 30-min incubation. Candidate FGFR3-binding phages were compared to the wildtype fd-tet phage (dark green).
Viruses 14 00384 g010
Figure 11. Evaluation of receptors interacting with the SARS-CoV-2 S1 spike protein by an indirect ELISA. Recombinant SARS-CoV-2 S1-His protein was incubated with bound ACE2, DPP4, FGFR3, NRP1, or BSA protein in triplicate wells for 1 h at 37 °C. The wells were incubated with a rabbit anti-SARS-CoV-2 spike IgG and an HRP-conjugated goat anti-rabbit IgG for 1 h at 37 °C. Signal was produced using a TMB substrate and endpoint absorbance at 450 nm was measured after a 30-min incubation. Mean endpoint absorbances for recombinant protein receptors were compared to the means of an unrelated BSA control using a Dunnett’s test with statistically significant differences (p < 0.05) indicated with a star.
Figure 11. Evaluation of receptors interacting with the SARS-CoV-2 S1 spike protein by an indirect ELISA. Recombinant SARS-CoV-2 S1-His protein was incubated with bound ACE2, DPP4, FGFR3, NRP1, or BSA protein in triplicate wells for 1 h at 37 °C. The wells were incubated with a rabbit anti-SARS-CoV-2 spike IgG and an HRP-conjugated goat anti-rabbit IgG for 1 h at 37 °C. Signal was produced using a TMB substrate and endpoint absorbance at 450 nm was measured after a 30-min incubation. Mean endpoint absorbances for recombinant protein receptors were compared to the means of an unrelated BSA control using a Dunnett’s test with statistically significant differences (p < 0.05) indicated with a star.
Viruses 14 00384 g011
Figure 12. Primary protein structures of the RBD domain of SARS-CoV-2 (upper line) and SARS-CoV (bottom line) aligned using BLASTP. Amino acid residues involved in the interaction (designated by capital bold letters) between FGFR3 (marked with a *) and ACE2 (marked with a ^) were identified using the YASARA Structure and literature data [97]. Promiscuous AA of the SARS-CoV-2 RBD involved in binding both FGFR3 and ACE2 are marked with *^. Members of AA cluster D442, N448, Y449, S494, Q493, L492, corresponding to phage mimotope EDYSELVSQ, are underlined and marked with (+).
Figure 12. Primary protein structures of the RBD domain of SARS-CoV-2 (upper line) and SARS-CoV (bottom line) aligned using BLASTP. Amino acid residues involved in the interaction (designated by capital bold letters) between FGFR3 (marked with a *) and ACE2 (marked with a ^) were identified using the YASARA Structure and literature data [97]. Promiscuous AA of the SARS-CoV-2 RBD involved in binding both FGFR3 and ACE2 are marked with *^. Members of AA cluster D442, N448, Y449, S494, Q493, L492, corresponding to phage mimotope EDYSELVSQ, are underlined and marked with (+).
Viruses 14 00384 g012
Table 1. Characteristics of recombinant human cell receptors and viral spike proteins.
Table 1. Characteristics of recombinant human cell receptors and viral spike proteins.
ProteinVendorCatalog #ResiduesLengthMolecular Weight
(kDa)
StartEnd
hACE2Sino Biological10108-H08HQ18S74073485.1
hDPP4Sino Biological10688-H08HD34P76674486.3
hNRP1Sino Biological10011-H08HF22K64463471.3
hFGFR3Sino Biological16044-H08HE23G37536439.6
SARS-CoV-2 S1 spikeSino Biological40591-V08HV16R68568176.5
Table 2. Amino acid clusters identified on the surface of ligands and their corresponding phage mimotopes. Bold letters highlight the amino acid residues shared with the phage mimotope.
Table 2. Amino acid clusters identified on the surface of ligands and their corresponding phage mimotopes. Bold letters highlight the amino acid residues shared with the phage mimotope.
Phage MimotopeAmino Acid Clusters Identified in:
FGF1NRP1Spike RBD
AGSGQEYLDTG-GQSY--ASSG-E--DAGS-Q-YI-
DGRADLSYDDGR-D-S-DDSRGELNYE-GNSDLSYN
DGRMTVYNEDGRLTV-GDDGR-NM-NEDGR--VYNQ
DGRRDVADDDG-RD-SDTEG-RDFGNDDGQREVGD-
DSFVNAPEDDT-VDLPHDD-F-NG-ED-AFVNGPE-
EDYSELVSQENH-EL---EDFSEFT-QDNYSQLVS-
ETRVEPEYDETQV-PKYN--R-DPEYDEVRAQP-Y-
VNADSASSDVDRDSRSS-VPAKSTSSDLN-DSVSS-
VPTDTDYSD-P-DTDYS-VPSTDYSVPNDSNYS-
VQAFDDSVDTQSYGDTVDTQ-FDQS-DIQAY-ESID
VTALEPGQDI-ALE-GE--SAIEPGQD-TAVEPG--
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Petrenko, V.A.; Gillespie, J.W.; De Plano, L.M.; Shokhen, M.A. Phage-Displayed Mimotopes of SARS-CoV-2 Spike Protein Targeted to Authentic and Alternative Cellular Receptors. Viruses 2022, 14, 384. https://doi.org/10.3390/v14020384

AMA Style

Petrenko VA, Gillespie JW, De Plano LM, Shokhen MA. Phage-Displayed Mimotopes of SARS-CoV-2 Spike Protein Targeted to Authentic and Alternative Cellular Receptors. Viruses. 2022; 14(2):384. https://doi.org/10.3390/v14020384

Chicago/Turabian Style

Petrenko, Valery A., James W. Gillespie, Laura Maria De Plano, and Michael A. Shokhen. 2022. "Phage-Displayed Mimotopes of SARS-CoV-2 Spike Protein Targeted to Authentic and Alternative Cellular Receptors" Viruses 14, no. 2: 384. https://doi.org/10.3390/v14020384

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop