SARS-CoV-2 Preimmune IgM Epitopes

doi:10.21203/rs.3.rs-21020/v1

While studying the human public IgM igome as represented by a library of 224087 linear mimotopes, three exact matches to peptides in the proteins of SARS-CoV-2 were found: two in the open reading frame 1ab and one in the spike protein. Joining the efforts to fast track SARS-CoV-2 vaccine development, here we describe briefly these potential epitopes in comparison to mimotopes representing peptides of SARS-CoV , HCoV 229E and OC43.

Immunology

Bioinformatics

Virology

SARS-CoV-2

epitope

mimotope

B cell precursors

IgM

The COVID19 pandemic has put to test the capacity of vaccinology to produce as fast as possible relevant vaccines. A number of recent reports predict possible SARS-CoV–2 epitopes for vaccine development but there are no reports on experimentally defined B cell epitopes [1; 2; 3; 4; 5]. The closest to identification of actual epitopes is the finding of pentapeptide sequences from the viral proteome in other known epitopes form IEDB [5]. A library of 224 087 mimotopes corresponding to the human public IgM repertoire as represented in a plasma pool from 10 000 healthy donors was recently designed [6]. The mimotopes were selected from a commercial 7 amino acid random peptide phage display library (Ph.D. 7, New England Biolabs). Conceptually, this mimotope library reflects at a certain level of detail, the repertoire of IgM specificities in the plasma focusing on the recurring ones. The latter can be just natural antibodies or they may represent the product of fast extrafollicularly expanding of IgM clones that may serve as precursors of highly specific, somatically mutated, class-switched B cells. The preimmune repertoire has to be quasi-complete to provide for rapid expansion of clones reactive with any newly encountered antigen. The same may not be true for our library although, due to the polyspecific binding, most of the available public repertoire may be partially represented in it [6]. Here we report that the IgM mimotope library contains heptapeptides identical to peptides in the proteome of SARS-CoV–2. One of them may serve as a potentially neutralizing epitope on the spike protein.

The design and the properties of the mimotope database has been published elsewhere [6]. The available sequences of the genomes of SARS-CoV (NC_004718.3), SARS-CoV–2 (ASM985889v3), HCoV229E (NC_002645.1) and HCoVOC43 (AY391777.1) were split into consecutive overlapping heptamers shifted by one residue and the resultant sequence sets were compared to the sequences in the database of natural mimotopes. Only exact matches were considered.

The homologous sequences in the non-redundant databases of the human proteome and Viridae (taxid:10239) were blast searched using the NCBI blastp suite (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE = Proteins).

As part of an ongoing analysis, the natural mimotope database was represented as a graph by connecting the sequences having at least 5 exact matches (i.e. of maximal Hamming distance 2). The graph had one giant component containing approximately 90% of the sequences which was further considered as the graph of interest. For the present study, the degrees of the vertices representing the natural SARS-CoV–2 epitopes, all of which belonged to the giant component, were used as the number of adjacent mimotopes parameter. For a set of words of length l made by alphabet of L symbols, the theoretical average number of neighbors N at Hamming distance D was calculated using the following formula for the number of neighbors:

For the present study, L = 20, l = 7, and D<3. For the first layer of neighbors N1 = 133 and for the second N2 = 7581. Under the hypothesis that the database is a random sample from the set of heptamer peptides, the probability of the occurrence of each neighbor is:

p = 224087/207≈1.75e–4,

and the expected mean number of distinct neighbors at D<3 was calculated as p.(N1,7+N2,7)≈1.33. The value of p was used subsequently also in a binomial test to calculate the probabilities of finding equal or higher number of adjacent mimotopes (Table 1).

The structure of the spike of SARS-CoV–2 was recently published (6vsb.pdb [2]). It was used to visualize the molecular context of the spike epitope found. The visualization of the structure and the calculation of the relative solvent exposed surface were done using UCSF Chimera, developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, with support from NIH P41-GM103311.

A simple comparison for exact matches to peptides from the SARS-CoV–2 proteome yielded 3 heptapeptides - two in the open reading frame 1ab (³⁵¹⁸AQTGIAV³⁵²⁴ and ⁵¹⁹⁸TKGPHEF⁵²⁰⁴) and one in the spike protein (¹⁰⁸TTLDSKT¹¹⁴). The E value of such a search result is very high and the mere number of sequences is not statistically significant. Yet, this does not refute the fact that 3 heptapeptides which are operationally defined mimotopes of human preimmune antibodies, are part of the viral proteome and, thus, represent (parts of) possible epitopes. On the other hand, the mimotopes in the database sometimes form non-random clusters of homologous sequences much like the mimotopes selected by a single monoclonal antibody. Each one among 224087 randomly selected heptamers should have on the average 1.33 homologous sequences in that same database that differ from it by up to 2 mismatches. As seen from Table 1, all SARS sequences but not those from trivial HCoV were members of clusters significantly greater than random (Binomial test, p<0.05, false discovery rate adjusted). This is an indication that the presence of these sequences is non-random and they represent clusters of mimotopes representing individual (poly)specificities.

An important prerequisite for the functionality of these epitopes is their degree of exposure to the solvent. The recently published structure of the spike (S protein) of SARS-CoV–2 [2] shows that ¹⁰⁸TTLDSKT¹¹⁴ forms a loop exposed to the solvent (Fig 1A.). The relative solvent exposed surface greatly exceeds the threshold of 5% for participating in contacts (Fig.1B). This loop is adjacent to the loop representing the epitope of the neutralizing antibody LCA60 on the SARS-CoV spike [7; 8]. Presumably, it is similarly exposed further in the open conformation of the spike domains. The adjacent N-glycosylation sites are N165 and N234. Dependent on the size of the carbohydrate sidechains, they may partially occlude the epitope.

The closest sequences in the human proteome are ⁵⁴⁰tlTLDSKT⁵⁴⁷ of the prostate-specific transglutaminase (TGM4) and ⁴⁶²TTLDSKi⁴⁶⁸ of mucin–16 (also known as ovarian tumor marker CA125, Q8WXI7.3, [9]). Both are on tumor associated antigens [9; 10]. While TGM4 is an intracellular antigen, mucin–16 is highly accessible on cell surfaces and in a soluble form. The mucin sequence ⁴⁶²TTLDSKI⁴⁶⁸ is T/S biased, represents part of the highly O-glycosylated N-terminal part of mucin–16 and is predicted to be O-glycosylated itself. Normally, such mucin protein core epitopes are occluded by glycosylation and thus, cryptic with respect to immune tolerance. Yet, monoclonals to similar epitopes turned out to bind specifically to tumour expressed mucins [11; 12; 13; 14] which are aberrantly/hypo glycosylated.

The sequences ¹⁰⁸TTLDSKT¹¹⁴has several exact matches in viruses outside the family Coronaviridae in hypothetical proteins of various phages. At least one of them infects L. plantarum which is a common species in the gut microbiome.

It is not surprising that the public IgM repertoire has clones potentially capable of binding to non-conserved regions of novel viruses. Similarly, the IgM igome contained sequences found also in SARS-CoV, although the epidemic was too restricted to be reflected in the antibody repertoires of the donors (Table 1). As expected, no signs of persistent antibody titers after exposure were observed. The representation of clones reactive with the trivial human coronaviruses 229E and OC43 was rather narrower than that of the unknown strains. Some of the epitopes were conserved between SARS-CoV and SARS-CoV–2 (AQTGIAV and TKGPHEF) but they were found in non-structural proteins and are hardly targets for neutralizing antibodies (Table 1). On the other hand, all potential epitopes found could play a role in targeting the viral proteins to specific B1 cells which produce the bulk of natural IgM. The latter are known to be excellent antigen presenting cells able to prime CD4⁺ T cells, and initiate immune responses [15] in antigen specific manner much like activated specific B2 cells [16]. It has been shown that B1 cells secreted IgM is a non-redundant and essential arm of the humoral responses to influenza in mice [17]. This implies that natural antibody epitopes might be essential components of subunit vaccines even though they may not represent typical dominant epitopes. The role of overlapping T and B cell epitopes is not clear except when the B cell receptor has a high enough affinity for the epitope to protect it during processing [18], but it is interesting that one of the SARS-CoV natural epitopes (⁹²²TTSTALG⁹²⁸) is also part of a CD4 T cell epitope in the context of HLA-DR 0401[19].

The current thinking separates the repertoire of natural and induced antibodies [20]. The preimmune IgM epitopes we describe could be epitopes of naïve B clones which are poised to initiate an immune response or epitopes of natural antibodies. As to the capacity of these epitopes to induce fully mature antibody response, it is interesting to note that the two preimmune IgM epitopes found for the spike of SARS-CoV (922TTSTALG928 and 389VKGDDVR395) are proven antibody targets in approximately one fourth of the SARS patients [21]. Furthermore, a recent report indicates the importance of IgM antibodies in the control of the diseases in mild case of COVID19 [22]. Thus, it is quite possible that the SARS-CoV–2 spike epitope TTLDSKT is bound by B cells that will contribute to the induced immune response. The correlation between the actual reactivities in SARS-CoV patients’ sera and the Bepipred score (Fig. 1C) confirms the low power of linear B cell epitope predicting algorithms, and underlies the necessity to base the proposals of new epitopes as much as possible on actual binding data.

These considerations make the novel SARS-CoV–2 epitopes valid targets in the search for a vaccine for COVID–19. The whole paradigm followed here focuses exclusively on the relatively rare linear epitopes. A lot more information about conformational epitopes may be hidden in the natural mimotope database but the approaches for sorting out clusters of mimotopes defining a conformational epitope are still being developed. The proposed actual preimmune IgM epitopes of SARS-CoV–2 can be instrumental both as parts of subunit vaccines or in the design of nanoparticle-based vaccines but also in the development of therapeutic monoclonal antibodies.

Table 1. Human public IgM repertoire (igome) selected mimotopes and their exact matches in the proteomes of SARS-CoV-2, SRAS-CoV, HCoV 229E and HCoV OC43.

	Mimotopes	Protein_ID	Protein	Strain	Starting_pos	Number of Adjacent³ Mimotopes
1	AQTGIAV	YP_009724389.1	orf1ab¹ nsp5²	SARS-CoV-2	3518	6^*
2	TKGPHEF	YP_009724389.1	orf1ab nsp12	SARS-CoV-2	5198	22^*
3	TTLDSKT	YP_009724390.1	spike	SARS-CoV-2	108	7^*
4	HSYGIDL	NP_828849	orf1ab nsp2	SARS-CoV	134	5^*
5	TTYNGYL	NP_828849	orf1ab nsp3	SARS-CoV	1460	6^*
6	AQTGIAV	NP_828849	orf1ab nsp5	SARS-CoV	3495	6^*
7	TKGPHEF	NP_828849	orf1ab nsp12	SARS-CoV	5175	22^*
8	VKGDDVR	NP_828851	spike	SARS-CoV	389	8^*
9	TTSTALG	NP_828851	spike	SARS-CoV	922	6^*
10	QLSLSMA	NP_828859	hypothetical	SARS-CoV	52	11^*
11	GAGDAGH	NP_073549	orf1ab nsp3	229E	2483	5^*
12	QTSQALQ	NP_073551	spike	229E	818	2
13	ANSFRLF	NP_073555	M protein	229E	96	2
14	NGSWVLN	AAR01012	orf1ab nsp3	OC43	2987	2

¹ orf – open reading frame,

²nsp – non-structural protein,

³Adjacent mimotopes are considered those that have no more than 2 mismatches (Hamming distance<3), ^*p<0.05, Binomial test, fdr adjusted.

Author Contributions

VS – conceptualizing, manuscript preparation, data analysis

PP – data analysis

AP - data analysis, conceptualizing, manuscript preparation,

Funding

This work was funded partially by the Bulgarian Fund for Scientific Research Grant D01-11/2016.

Acknowledgments

The authors wish to thank Prof. Angel Galabov for critical reading of the manuscript

Data Availability Statement

The datasets analyzed and the scripts for this study can be found in the GitHub Repository https://github.com/ansts/SARS-CoV-2.

[1] M. Zheng, and L. Song, Novel antibody epitopes dominate the antigenicity of spike glycoprotein in SARS-CoV-2 compared to SARS-CoV. Cellular & molecular immunology (2020).

[2] A.C. Walls, Y.-J. Park, M.A. Tortorici, A. Wall, A.T. McGuire, and D. Veesler, Structure, function and antigenicity of the SARS-CoV-2 spike glycoprotein. bioRxiv (2020) 2020.02.19.956581.

[3] A. Grifoni, J. Sidney, Y. Zhang, R.H. Scheuermann, B. Peters, and A. Sette, Candidate targets for immune responses to 2019-Novel Coronavirus (nCoV): sequence homology- and bioinformatic-based predictions. bioRxiv (2020) 2020.02.12.946087.

[4] S.F. Ahmed, A.A. Quadeer, and M.R. McKay, Preliminary Identification of Potential Vaccine Targets for the COVID-19 Coronavirus (SARS-CoV-2) Based on SARS-CoV Immunological Studies. Viruses 12 (2020) 254.

[5] G. Lucchese, Epitopes for a 2019-nCoV vaccine. Cellular & molecular immunology (2020).

[6] A. Pashov, V. Shivarov, M. Hadzhieva, V. Kostov, D. Ferdinandov, K.-M. Heintz, S. Pashova, M. Todorova, T. Vassilev, T. Kieber-Emmons, L.A. Meza-Zepeda, and E. Hovig, Diagnostic Profiling of the Human Public IgM Repertoire With Scalable Mimotope Libraries. Frontiers in Immunology 10 (2019).

[7] D. Wrapp, N. Wang, K.S. Corbett, J.A. Goldsmith, C.-L. Hsieh, O. Abiona, B.S. Graham, and J.S. McLellan, Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science (2020) eabb2507.

[8] A.C. Walls, X. Xiong, Y.J. Park, M.A. Tortorici, J. Snijder, J. Quispe, E. Cameroni, R. Gopal, M. Dai, A. Lanzavecchia, M. Zambon, F.A. Rey, D. Corti, and D. Veesler, Unexpected Receptor Functional Mimicry Elucidates Activation of Coronavirus Fusion. Cell 176 (2019) 1026-1039.e15.

[9] M. Felder, A. Kapur, J. Gonzalez-Bosquet, S. Horibata, J. Heintz, R. Albrecht, L. Fass, J. Kaur, K. Hu, H. Shojaei, R.J. Whelan, and M.S. Patankar, MUC16 (CA125): tumor biomarker to cancer therapy, a work in progress. Molecular cancer 13 (2014) 129.

[10] Z. Cao, Y. Wang, Z.-Y. Liu, Z.-S. Zhang, S.-C. Ren, Y.-W. Yu, M. Qiao, B.-B. Zhai, and Y.-H. Sun, Overexpression of transglutaminase 4 and prostate cancer progression: a potential predictor of less favourable outcomes. Asian J Androl 15 (2013) 742-746.

[11] J. Burchell, and J. Taylor-Papadimitriou, Effect of modification of carbohydrate side chains on the reactivity of antibodies with core-protein epitopes of the MUC1 gene product. Epithelial Cell Biol 2 (1993) 155-162.

[12] J. Burchell, S. Gendler, J. Taylor-Papadimitriou, A. Girling, A. Lewis, R. Millis, and D. Lamport, Development and characterization of breast cancer reactive monoclonal antibodies directed to the core protein of the human milk mucin. Cancer research 47 (1987) 5476-5482.

[13] E. Petrakou, A. Murray, and M.R. Price, Epitope mapping of anti-MUC1 mucin protein core monoclonal antibodies. Tumour biology : the journal of the International Society for Oncodevelopmental Biology and Medicine 19 Suppl 1 (1998) 21-29.

[14] D. Zhou, L. Xu, W. Huang, and T. Tonn, Epitopes of MUC1 Tandem Repeats in Cancer as Revealed by Antibody Crystallography: Toward Glycopeptide Signature-Guided Therapy. Molecules (Basel, Switzerland) 23 (2018) 1326.

[15] N. Baumgarth, A two-phase model of B-cell activation. Immunol Rev 176 (2000) 171-80.

[16] A.F. Popi, I.M. Longo-Maugéri, and M. Mariano, An Overview of B-1 Cells as Antigen-Presenting Cells. Frontiers in immunology 7 (2016) 138-138.

[17] N. Baumgarth, O.C. Herman, G.C. Jager, L.E. Brown, L.A. Herzenberg, and J. Chen, B-1 and B-2 cell-derived immunoglobulin M antibodies are nonredundant components of the protective response to influenza virus infection. J Exp Med 192 (2000) 271-280.

[18] P.D. Simitsek, D.G. Campbell, A. Lanzavecchia, N. Fairweather, and C. Watts, Modulation of antigen processing by bound antibodies can boost or suppress class II major histocompatibility complex presentation of different T cell determinants [see comments]. J Exp Med 181 (1995) 1957-63.

[19] J. Yang, E. James, M. Roti, L. Huston, J.A. Gebe, and W.W. Kwok, Searching immunodominant epitopes prior to epidemic: HLA class II-restricted SARS-CoV spike protein epitopes in unexposed individuals. International immunology 21 (2009) 63-71.

[20] N. Baumgarth, J.W. Tung, and L.A. Herzenberg, Inherent specificities in natural antibodies: a key to immune defense against pathogen invasion. Springer Semin Immunopathol 26 (2005) 347-62.

[21] Y. He, Y. Zhou, H. Wu, B. Luo, J. Chen, W. Li, and S. Jiang, Identification of immunodominant sites on the spike protein of severe acute respiratory syndrome (SARS) coronavirus: implication for developing SARS diagnostics and vaccines. Journal of immunology (Baltimore, Md. : 1950) 173 (2004) 4050-4057.

[22] I. Thevarajan, T.H.O. Nguyen, M. Koutsakos, J. Druce, L. Caly, C.E. van de Sandt, X. Jia, S. Nicholson, M. Catton, B. Cowie, S.Y.C. Tong, S.R. Lewin, and K. Kedzierska, Breadth of concomitant immune responses prior to patient recovery: a case report of non-severe COVID-19. Nature medicine (2020).

SARS-CoV-2 Preimmune IgM Epitopes

Status:

Version 1

Abstract

Figures

1. Introduction

2. Methods

3. Results and Discussion

Table

Declarations

References

Status:

Version 1