Sialic Acid and Fucose Residues on the SARS-CoV-2 Receptor-Binding Domain Modulate IgG Antibody Reactivity.

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein is a conserved domain and a target for neutralizing antibodies. We defined the carbohydrate content of the recombinant RBD produced in different mammalian cells. We found a higher degree of complex-type N-linked glycans, with less sialylation and more fucosylation, when the RBD was produced in human embryonic kidney cells compared to the same protein produced in Chinese hamster ovary cells. The carbohydrates on the RBD proteins were enzymatically modulated, and the effect on antibody reactivity was evaluated with serum samples from SARS-CoV-2 positive patients. Removal of all carbohydrates diminished antibody reactivity, while removal of only sialic acids or terminal fucoses improved the reactivity. The RBD produced in Lec3.2.8.1-cells, which generate carbohydrate structures devoid of sialic acids and with reduced fucose content, exhibited enhanced antibody reactivity, verifying the importance of these specific monosaccharides. The results can be of importance for the design of future vaccine candidates, indicating that it is possible to enhance the immunogenicity of recombinant viral proteins.

Introduction

The adaptive immune response to SARS-COV-2 depends on T-cells that direct the immune responses and contribute to killing of infected cells and on antibody-producing B-cells.[1] Seroconversion has been detected in 93–99% of patients with diagnosed SARS-CoV-2 infection, with disease severity correlating with antibody titres.[2−4] Neutralizing antibodies (NAb) are a key component in the response toward viruses, and an important aspect after immunization is whether the generated antibodies possess neutralizing capabilities.[5] In the case of SARS-CoV-2, many neutralizing antibodies recognize the receptor-binding motif (RBM) within the receptor-binding domain (RBD) of the spike (S) protein[6] and execute their neutralizing capacity by sterically hindering viral binding to the angiotensin converting enzyme 2 (ACE2) receptor[7] or by keeping the RBD in its “down-conformation”.[6] However, there are reports about neutralizing antibodies targeting epitopes also outside the RBM, such as 47D11, which binds to the conserved core of the RBD,[8] or S309, which recognizes a conserved epitope involving interactions with the fucose and other glycan moieties of the N343-glycan within the RBD.[9] Neutralizing antibodies have also been found to target the N-terminal domain (NTD) of the S protein. For example, the NAb 4A8 targets residues within the NTD, including the N147-glycosite, and neutralizes possibly by inhibiting conformational changes of the S protein.[10] Thus, the glycosylation profile of the S protein appears to be important for antibody recognition and neutralization.

The S protein is inserted in the viral envelope as a trimer, forming the characteristic “spikes” protruding from the viral surface. The S protein is cleaved by host proteases to form S1 and S2. The S1 domain contains the RBM and mediates binding to the ACE2 receptor, while fusion with the host cell membrane is mediated by S2.[11−13] The S protein ectodomain contains 22 consensus sites for N-linked glycosylation (Asn-X-Ser/Thr where X is any amino acid except Pro). Most of the sites have been reported as glycosylated, carrying complex and high-mannose glycans for recombinant S proteins, expressed in cell culture.[14−17] Although most consensus sites for N-linked glycosylation in the S protein appear to be occupied, the composition and structure of glycans at respective sites appear highly variable.[15,16] The O-linked glycosylation pattern of the S protein is not entirely established, although the presence of several O-linked glycans has been identified within the RBD.[14,17−19]

The vaccines against SARS-CoV-2 induce antibodies that after immunization target specific domains of the S protein. The vector-based DNA vaccines and mRNA vaccines utilize the human glycosylation profile on the produced protein, while the glycosylation profile of protein-based subunit vaccines is dependent on the cell type used for production.[20] It has been shown that it is possible to alter the glycan content of the S protein without affecting the serological properties,[21] and by truncating the glycans, it was possible to enhance protection in an animal model.[22] However, the glycan content is also important for correct folding of the S protein and may influence the antigen stability and could therefore indirectly affect the distribution and presentation of the antigen.[23−25] Thus, the glycosylation profile of the S protein may directly affect the antibody recognition by shielding specific epitopes, or indirectly by altering the architecture of the protein, leading to variability in the effectivity of a potential vaccine candidate.[26,27]

In this work, we have characterized the glycan content of a recombinant RBD protein expressed in three different mammalian cell lines and showed a diverse glycan composition at each site. The N- and O-linked glycans were stepwise modulated using enzymatic degradation. Serum samples from patients previously infected with SARS-CoV-2 were used to assess the impact of glycan composition on antibody reactivity. A glycan hot spot within the RBD was found to be essential for antibody reactivity. In addition, modulation of the glycan content revealed specific monosaccharides that were able to enhance the antibody reactivity.

Results

Glycosylation Pattern of the Recombinant RBD Produced in CHO-S and HEK293F Cells

The recombinant RBD, produced in HEK293F- and CHO-S-cells, respectively, was subjected to nano-LC–MS/MS analysis. We defined the level of occupancies and composition of the N-linked and O-linked glycans present in the RBD. In addition, we suggest N- and O-linked glycan structures based on the observed glycan compositions and the knowledge of the mammalian glycan biosynthesis pathways. The HEK293F-produced RBD showed a nearly complete occupancy for both N-linked sites (99.1 and 100%), while the CHO-S-produced construct presented a partial occupancy of 93.3% for site N331 and a full occupancy for N343 (Figure A,B). Complex-type N-linked glycans were the most abundant structure in both cell lines; still, a higher degree of glycans processed to complex type was associated with the HEK293F-cell line, while oligomannose structures were relatively more abundant for the CHO-S-produced protein. The observed CHO-S-produced oligomannose glycans were different at the two sites with N331 displaying higher levels of oligomannose-6-phosphate glycans (Figure A,B and Table ).

Figure 1

Schematic presentation of glycan distribution at the respective sites of the CHO-S- and HEK293F-produced RBD. The glycan structures are inferred from the obtained mass spectrometry data as well as from previously reported glycan structures and known glycan biosynthesis pathways in mammalian cells. (A) Degree of glycosylation, distribution between the glycan types, degree of sialylation, and degree of fucosylation at site N331 as detected on the CHO-S- and HEK293F-produced RBD. (B) Degree of glycosylation, distribution between the glycan types, degree of sialylation, and degree of fucosylation at site N343 as detected on the CHO-S- and HEK293F-produced RBD. (C) Degree of glycosylation, distribution between glycan compositions, and degree of sialylation at site T323/S325 as detected on the CHO-S- and HEK293F-produced RBD. (D) Degree of glycosylation, distribution between glycan compositions, and degree of sialylation at site T523 as detected on the CHO-S- and HEK293F-produced RBD. (E) Glycosylation of the recombinant RBD produced in CHO-S cells with the most prevalent glycans drawn at the respective site. The yellow circle highlights the glycan hotspot. (F) Glycosylation of the recombinant RBD produced in HEK293F cells with the most prevalent glycans drawn at the respective site. The yellow circle highlights the glycan hotspot.

Table 1

Percentage Distribution of Glycan Types at Sites N331 and N343 When Produced in CHO-S, HEK293F, and Lec3.2.8.1 Cells

	N331			N343
	CHO-S	HEK293F	Lec3.2.8.1	CHO-S	HEK293F	Lec3.2.8.1
unoccupied	7.8	0.9	2.5	0.0	0.0	0.0
oligomannose	1.3	0.4	76.1	12.1	4.9	79.8
oligomannose-6-P	19.1	2.7	10.5	1.4	0.2	2.3
paucimannose	2.1	2.0	7.4	0.6	0.6	12.4
hybrid	2.5	2.1	2.2	1.6	0.5	5.5
complex	67.3	91.8	1.3	84.3	93.8	0.0

Among the complex-type N-linked glycans, biantennary structures were most frequently found at both positions (Table S1). Despite similar glycan compositions in both cell lines, the fragment spectral evaluation identified prominent differences for glycans produced in CHO-S and HEK293F cells. The major difference was the prominent LacDiNAc-containing structures in the HEK293F-produced RBD, while those were absent in the CHO-S-produced protein (Table S1).

Within a given cell line, the frequency of fucose residues was similar for both sites, while the overall fucosylation was higher for HEK293F compared to CHO-S (Figure A,B and Table ). For HEK293F and CHO-S, fucosylation was observed for paucimannose, hybrid, and complex structures. The Lec3.2.8.1 cell mainly produced the oligomannose structures with fucosylation observed on HexNAc(2)Hex(5). As all glycan groups were observed being fucosylated, all of them were included in calculation of the total fucosylation level. Not only the total fucosylation level but also the degree of fucosylation (number of fucose residues per glycan) differs between the cell lines (Table S2). CHO-S cells predominantly produced monofucosylated structures with the fucose placed at the core, as based on the fragment ion analysis. Multiple fucosylation, with up to four fucose residues per glycan, was observed for the RBD produced in the HEK293F cell line. The attachment of fucose to LacNAc and LacDiNAc was observed based on the fragment spectral evaluation.

Table 2

Percentage of Detected Glycans Carrying at Least One Fucose or Sialic Acid at Sites N331 and N343 When Produced in CHO-S, HEK293F, and Lec3.2.8.1 Cells

	N331			N343
	CHO-S	HEK293F	Lec3.2.8.1	CHO-S	HEK293F	Lec3.2.8.1
fucosea	67.5	95.7	31.4	78.4	93.5	6.3
fucoseb	73.2	96.6	31.4	78.4	93.5	6.3
sialic acida	38.9	35.2	0.0	40.8	13.3	0.0
sialic acidc	65.0	38.4	0.0	47.9	14.1	0.0

The percentage is calculated in relation to all glycosylated forms sharing the same peptide sequence, including nonglycosylated peptides.

The percentage is calculated in relation of a total amount of all observed glycosylated forms, and unoccupied peptides were excluded when performing the calculations.

The percentage is calculated in relation of a total amount of all observed hybrid- and complex-type glycoforms, and the unoccupied peptides and all oligomannose glycoforms were excluded from the calculations.

In contrast to fucosylation, the sialylation level, calculated including all glycan groups that can be sialylated, i.e., hybrid- and complex-type glycans, was lower for HEK293F compared to CHO-S (Figure A,B and Table ). Also, the degree of sialylation (number of sialic acid residues per glycan) differed between the cell lines (Table S3). CHO-S cells produced multiple sialylated forms in contrast to HEK293F, where mainly monosialylated structures were observed. The lower sialylation level in glycans produced by HEK293F is likely a result of the extensive fucosylation in this cell type. In both cell lines, the major N-glycan type carrying sialic acid was the complex type and a difference between the sites was noted, with N-linked glycans at position N331 displaying higher sialylation levels.

The O-linked glycans were similar for the two cell types (Figure C,D and Table S4). The O-linked glycan close to the NTD of the RBD could not be defined to a single amino acid, due to the absence of fragment ions between the two adjacent potential sites, and thus could be placed either at amino acid position T323 or S325. The site T323/S325 was glycosylated to a high degree (97 and 91%, for CHO-S and HEK293F, respectively), while T523 was scarcely decorated and mainly remained nonglycosylated in both the CHO-S- and HEK293F-produced RBD (5 and 1%, respectively). Comparison of O-linked glycans at the individual sites revealed more extensive processing in the HEK293F-produced protein, while the CHO-S-produced O-linked glycans almost exclusively consisted of core 1 structures (Figure C,D). The degree of sialylated structures at site T323/S325 was similar between the cell lines (80 and 82% for CHO-S and HEK293F, respectively), while the degree of monosialylated structures (66 and 35%, respectively) and disialylated structures (15 and 46%, respectively) differed. Similarly, the frequency of sialylated structures at site T523 was similar between CHO-S (80%) and HEK293F (79%) cells. The degree of monosialylation was higher in the CHO-S-produced RBD, as compared to the HEK293F-produced protein (38 and 14%, respectively), while a higher degree of disialylation was seen on the HEK293F-produced RBD (42 and 62%, respectively) (Figure C,D).

In summary, the RBD produced in CHO-S cells carried O-linked glycans at two positions although only position T323/S325 appeared to be glycosylated with a high frequency. The main type of O-linked glycan found at this position was a core 1 structure with a single sialic acid at the distal galactose (Figure E). This RBD protein also carried two N-linked glycans at positions N331 and N343. The predominant type of N-linked glycan was the biantennary complex type, although many variants of complex-type glycans were found. The RBD produced in HEK293F cells predominantly carried core 1 O-linked glycans with two sialic acids, one attached to the distal galactose and one to the innermost GalNAc residue. The N-linked glycans on the HEK293F RBD were almost exclusively of complex type with a high degree of fucosylation (Figure F). To note, all glycan structures presented in Figure are inferred from the obtained mass spectrometry data as well as from previously reported glycan structures and known glycan biosynthesis pathways in mammalian cells. Also, the measurements and calculations of abundances used for values presented in Figure , Tables and 2 are described in detail in the Materials and Methods Section, The Analysis of RBD Glycosylation Section, and Tables S7 and S8.

Evaluation of Convalescent Sera from COVID-19 Patients

Serum samples were collected from 24 individuals previously infected with SARS-CoV-2 as determined by a PCR-positive nasopharyngeal sample. Blood samples were collected 25–100 days following positive diagnosis. All sera were characterized with respect to anti-RBD IgG levels and the capability to neutralize a DE-Gbg20 strain of SARS-CoV-2 grown in VERO-cells. Based on the neutralization capability, the sera were divided to three groups: non-neutralizing (NT negative, n = 7), weakly neutralizing (NT titre 3–6, n = 7), and highly neutralizing (NT titre 48–96, n = 10) (Figure and Table S5). High neutralization capability correlated well with high levels of IgG targeting the RBD, as all highly neutralizing sera also were anti-RBD IgG positive, while six of the seven serum samples in the weakly neutralizing group were anti-RBD IgG negative. Interestingly, four out of seven serum samples in the NT negative group were anti-RBD IgG positive.

Figure 2

NT-titre (left y axis, gray bars) and anti-RBD IgG-levels (right y axis, transparent circles) for the 24 characterized serum samples. The sera were divided into three groups based on the neutralizing capability: non-neutralizing (NT-negative, n = 7), weakly neutralizing (NT titre 3–6, n = 7), and highly neutralizing (NT titre 48–96, n = 10). Anti RBD-IgG value ≥50 AU/mL is considered positive.

Impact of Glycan Structures on Antibody Reactivity against the RBD

To assess the impact of the different types of glycan structures found within the RBD, we removed the N-linked, the O-linked, or a combination of both glycans using enzymatic treatment. Removal of glycans was verified by a size shift on an SDS-page gel, visualized by silver staining (Figure S1A). Based on the glycan structures identified in the present study, the used enzyme combination would in theory result in complete removal of O-linked glycans from the CHO-S- and Lec3.2.8.1-produced RBD and removal of 71.4% of the O-linked glycans from the HEK293F-produced RBD.

The effect of glycan removal on the antibody reactivity against the recombinant RBD was tested using the defined serum samples described above. The RBD produced in CHO-S cells elicited a strong reactivity to the highly neutralizing sera, with reduced reactivity following removal of N-linked glycans, O-linked glycans, or a combination of both (Figure A). The highly neutralizing sera showed reduced reactivity against the RBD with removed N-linked glycans and the RBD lacking both N- and O-linked glycans produced in HEK293F-cells, while no effect on reactivity against the RBD lacking only O-linked glycans was observed (Figure B). The enzymatic treatment utilized to remove O-linked glycans will also remove sialic acid residues from the remaining N-linked glycans. Thus, the observed effects on antibody binding to the RBD lacking O-linked glycans may reflect changes in the N-linked glycans. Weakly neutralizing sera did not show any reactivity to the CHO-S- or the HEK293F-produced RBD regardless of the glycosylation profile (Figure S2A,C). The non-neutralizing serum samples displayed low reactivity against all recombinant RBD with only a minor difference depending on the glycosylation status (Figure S2B,D). The intensity of the reactivity of individual serum samples against the recombinant RBD correlated well with the anti-RBD IgG levels detected in each serum (Figure S3). Only heat treatment for 24 h without addition of enzymes (mock) resulted in significantly higher sera reactivity for the CHO-S-produced RBD compared to the untreated variant, while no difference was noted for the HEK293F- and Lec3.2.8.1-produced RBD (Figure S4).

Figure 3

Reactivity of highly neutralizing sera (NT titre 48–96, n = 10) against the fully glycosylated RBD (mock treated) and against the deglycosylated RBD produced in CHO-S and HEK293F cells. (A) Removal of both N-linked and O-linked glycans, removal of N-linked glycans alone, or removal of O-linked glycans alone from the RBD produced in CHO-S cells. (B) Removal of both N-linked and O-linked glycans, removal of N-linked glycans alone, or removal of O-linked glycans alone from the RBD produced in HEK293F cells. (C) Removal of sialic acids alone from the RBD produced in CHO-S and HEK293F cells. (D) Removal of fucose alone from the RBD produced in CHO-S and HEK293F cells. Data information: dark red color symbolizes a serum with high levels of anti-RBD IgG, and white color indicates anti-RBD IgG-negative serum (<50 AU/mL). Statistical analysis was performed with the Wilcoxon matched-pair signed rank test, ** = p < 0.001.

To further assess the impact of specific glycan residues on antibody reactivity, sialic acids and fucose groups were enzymatically removed from the CHO-S- and HEK293F-produced RBD. SDS-page gel electrophoresis with silver stain and lectin blots were used to confirm the removal of sialic acids or fucose groups. Small, but distinct, size shifts were evident after the enzymatic treatments (Figure S1B,C), and lectin blots using MAL II further verified removal of sialic acids from both the CHO-S- and HEK293F-produced RBD (Figure S5A). Removal of fucose was verified by the use of UEA I-lectin blots (Figure S5B). Removal of sialic acids from the CHO-S-produced RBD resulted in a significant increase in the serum reactivity, as compared to the fully glycosylated RBD. A similar effect was seen following removal of sialic acids from the HEK293F-produced RBD, but the difference was less prominent (Figure C). Removal of fucose groups also resulted in a significant increase in serum reactivity for both the CHO-S- and HEK293F-produced RBD, with the HEK293F-produced construct showing a more prominent increase (Figure D).

To confirm the impact of sialic acids and fucose groups on the antibody reactivity, the RBD was produced in Lec3.2.8.1 cells deficient in synthesis of complex-type glycans. Highly neutralizing serum samples showed a significantly higher reactivity against the RBD produced in Lec3.2.8.1 cells, as compared to the CHO-S- or HEK293F-produced RBD constructs (Figure A). As expected, enzymatic removal of sialic acids and fucose from the Lec3.2.8.1-produced RBD did not confer any detectable size shift on an SDS-page gel (Figure S1D) or a change in antibody reactivity by highly neutralizing sera (Figure B).

Figure 4

Antibody reactivity of highly neutralizing sera (NT titre 48–96, n = 10). (A) Reactivity against the fully glycosylated RBD (untreated) expressed in CHO-S, HEK293F, and Lec3.2.8.1 cells. (B) Reactivity against the fully glycosylated RBD (mock treated) produced in Lec-3.2.8.1 cells, or against the Lec3.2.8.1-produced RBD following enzymatic removal of sialic acids and fucose. Data information: dark red color symbolizes a serum with high levels of anti-RBD IgG, and white color indicates anti-RBD IgG-negative serum (<50 AU/mL). Statisitical analysis was performed with the Wilcoxon matched-pair signed rank test, ** = p < 0.001.

Glycosylation of the Recombinant RBD Produced in Lec3.2.8.1 Cells

In order to verify that the recombinant RBD produced in Lec3.2.8.1 cells lacked complex-type glycans, it was subjected to nano-LC–MS/MS analysis. Both N-linked sites were found to be glycosylated to a high degree (98 and 88%, respectively). The structural distribution was similar between the sites, with high mannose as the dominating glycan type (Table ). No sialic acid or end-fucose was found; however, 6% of the structures at site N331 and 31% of the structures at site N343 carried core fucose (Table ). Position T323/S325 was frequently (98%) decorated with an O-linked glycan, while site T523 was more sparsely decorated (15%). A single HexNAc was the most frequent structure at both O-linked sites (Table S4).

Discussion

Viruses that infect humans do not carry their own glycosyltransferases but instead rely on the enzymes of the host cell in the processing of their glycans. Hence, viral envelope proteins are therefore glycosylated by the host glycosylation machinery that initially folds and then add glycans to them. This often results in a viral glycosylation profile where glycans protecting antibody epitopes are selected, which does not stimulate a strong immune response. Thus, in order to circumvent the antibody responses, many viruses utilize the host cell glycosylation machinery to cover B-cell epitopes with a dense network of N-linked glycans, which due to physical hindrance shield the epitopes and prevent binding by neutralizing antibodies.[28] With a few exceptions, N-linked glycans alone rarely act as antibody epitopes,[29] but the opposite effect has been observed with small O-linked glycans. Olofsson et al. showed that 70% of tested sera against herpes simplex virus (HSV) type 2 contains antibodies targeting a peptide decorated with a single O-linked GalNAc residue. Removal of the glycan moiety diminished this response.[30] Similarly, our group has previously shown that a single GalNAc-residue added to a naked peptide can alter the antibody binding toward specific domains of the glycoprotein E from the varicella zoster virus (VZV).[31]

It is established that viral glycoproteins are heterogeneously glycosylated when they are expressed either as recombinant proteins in cell culture or after natural infection in cells.[31−33] This could imply that the antibody response to a viral glycoprotein is more diverse than previously thought. Hence, serum from infected individuals contains a polyclonal antibody pool, which could recognize multiple epitopes, and various glycoforms can constitute parts of these epitopes. The S protein of SARS-CoV-2 is highly glycosylated, with 17 to 22 previously identified sites carrying N-linked glycosylation that can shield B cell epitopes.[15,17,19] Of these, 2 N-linked glycans are present in the RBD, and Yang et al. identified as many as 10 O-linked glycans in this region, although most of them appeared to be of low abundance and their biological significance is therefore uncertain.[34] We demonstrate that enzymatic removal of all N-linked and/or O-linked glycans resulted in decreased antibody reactivity of the recombinant SARS-CoV-2 RBD produced in both CHO-S and HEK293F cells. This could indicate (i) that the glycans of the RBD constitute parts of antibody epitopes, (ii) that glycosylation is required for proper protein folding and in maintaining the protein conformation,[35] or (iii) a combination of both.

The strong immunoreactivity of RBD, as verified by the identification of a multitude of NAbs targeting this domain,[6,8,9] could be explained by the presence of multiple structural epitopes. If large glycan moieties are lacking, as after enzymatic removal, the folding of the protein could be compromised and thereby also transform structural epitopes, which would render them inaccessible to antibodies present in the serum samples. No linear B-cell epitopes have been identified within the RBD.[36] However, screening for B cell epitopes is performed with synthetic peptides lacking glycans and would then not be able to identify epitopes that are dependent on the presence of N- or O-linked glycans.

We performed a structural screening of the glycan profile of the recombinant RBD. Overall, our screening is in line with previous studies,[14,16] but we found significant differences in the amount of sialic acid and fucose content when comparing the RBD produced in CHO-S cells and HEK293F cells, respectively. Interestingly, the CHO-S-produced RBD is also presented with a high degree of mannose-6-phosphate (M-6-P). This structure has previously been observed in the SARS-CoV-2 spike protein when expressed in cell lines but also when isolated from intact viral particles.[33,37] Mannose-6-phosphate is recognized by the M-6-P receptor present in the trans-Golgi compartment, and it directs tagged proteins to late endosomes/lysosomes. Lysosomal egress dependent on M-6-P has been described for both HSV (26) and VZV (27). Also, SARS-CoV-2 egress mediated by lysosomes has been proposed.[38] However, we observed only a minor fraction of the peptides carrying M-6-P and to what extent they potentially could contribute to viral particle egress remains to be clarified.

While each glycosite on the recombinant RBD was glycosylated at a similar frequency independent of the production cell line, the RBD produced in HEK293F cells had a higher degree of fucosylation compared to CHO-S cells. Selective removal of the fucose groups resulted in a significantly increased antibody reactivity. While abundant fucosylation was a trait of the HEK293F-produced construct, the RBD produced in CHO-S cells had a larger content of sialic acid moieties. Selective removal of sialic acids enhanced antibody reactivity for both constructs, but the effect was more prominent for the CHO-S produced construct. The use of the Lec3.2.8.1-cell line resulted in RBD with a glycosylation profile completely deficient in sialic acids and end-fucose. The observation of oligomannose structures with core fucosylation in Lec3.2.8.1-cells, although unclear how this type of structure could occur given the known pathways for N-linked glycosylation, has previously been described.[39] The antibody reactivity toward the Lec3.2.8.1-produced RBD was enhanced, and additional treatment to remove the core fucose did not result in any change in the antibody reactivity. Altogether, these results point to an important function of specific terminal-sugar residues in the antibody reactivity against glycosylated viral antigens and suggest that core fucosylation is of minor importance, despite the report of an NAb that specifically interacts with the core fucose of the N-linked glycan situated on position N343.[9]

In line with our findings that removal of sialic acids leads to increased antibody reactivity, the nonsialylated glycan structures of yeast-cell-produced proteins could possibly be part of the explanation of the highly efficient yeast-produced vaccines against HBV.[40,41] This suggests that it is possible to optimize recombinantly expressed RBD or S proteins in order to generate effective vaccine candidates. However, important to note is the possibility that immunization with a recombinant expressed subunit vaccine directs the humoral immune response toward B-cell epitopes with species-specific glycosylation profiles. This can possibly result in skewed immunodominance, directing the antibody response toward epitopes that are not exposed after a natural infection with the virus, resulting in disturbed efficiency of the vaccine.[42] The data presented in this work confirm the necessity of correct glycosylation and show that also small differences in the glycosylation profile of a viral antigen can have a large impact on the reactivity by antibodies generated after a natural infection with SARS-CoV-2. A conscious decision regarding the glycosylation traits of the production cell line could hence affect the antibody response triggered by a recombinant protein. We suggest that the glycosylation characteristics should be considered during the production of recombinant vaccines toward SARS-CoV-2 but also other enveloped viruses, which carry glycoproteins.

Materials and Methods

Expression of Recombinant S Protein Constructs

The RBD of the SARS-CoV-2 spike protein (amino acids 319–541) was produced in three cell lines using an expression vector obtained through BEI Resources, NIAID, NIH, which is vector pCAGGS containing SARS-CoV-2, and Wuhan-Hu-1 spike glycoprotein gene RBD with the C-terminal Hexa-Histidine tag (NR-52309) (Table S6).

CHO-S cells (Cat nr R80007, Thermo Fisher Scientific, Waltham, MA) were adapted to grow in suspension in the FectoCHO medium (Polyplus transfection, Illkirch-Graffenstaden, France) at 37 °C in 5% CO2 in Optimum Growth flasks (Thomson instrument company, Oceanside, CA) at 130 rpm in a Multitron 4 incubator (Infors, Bottmingen, Schweiz). Lec3.2.8.1 cells (a mutated CHO cell line kindly received from Prof. P Stanley[43]) were cultured under the same conditions. The HEK293 derivate HEK293F cell line (Cat nr R79007, Thermo Fisher Scientific) was cultured in the Freestyle 293 medium. Cells were transfected at 2 × 106 cells/mL using a FectoPro transfection reagent (Polyplus transfection). The temperature was reduced to 32 °C (Lec3.2.8.1) or 31 °C (CHO-S) 4 h post transfection, while transfected HEK293F cells were kept at 37 °C. Protein-containing culture supernatants (800 mL–1 L) were harvested when cell viability was below 80%, which was after 168 h (CHO-S), 74 h (Lec3.2.8.1), or 90 h (HEK293F), filtered using Polydisc AS 0.45 μm (Whatman, Maidstone, UK) and loaded onto a 5 mL HisExcel column (Cytiva, Marlborough, MA). After sample loading, the column was washed with 20 mM sodium phosphate, 0.5 M NaCl, and 30 mM imidazole before elution of the protein using the same buffer with 500 mM imidazole (Lec3.2.8.1-produced RBD) or 300 mM imidazole (CHO-S- and HEK293F-produced RBD). Pooled fractions were concentrated using 10 kDa Vivaspin concentrators (MWCO 10 kDa, Sartorius, Göttingen, Germany), passed over a HiPrep 26/10 desalting column (Cytiva) in phosphate-buffered saline, and finally concentrated again. The Lec3.2.8.1-produced RBD was further purified by gel filtration using a Superose 200 Increase 16/300 GL column (Cytiva) in phosphate-buffered saline. Integrity and purity of the different RBD preparations were checked by SDS-PAGE and Western blot.

Sample Preparation Prior to Assessment of the Position and Structure of Glycans

The purified RBD preparations from CHO-S, HEK293F, and Lec3.2.8.1 (20 μg each) were diluted with digestion buffer (DB), 1% sodium deoxycholate (SDC) in 50 mM triethylammonium bicarbonate (TEAB) pH 8.0 (Sigma Aldrich, St. Louis, MO), to give protein concentrations of 0.5 μg/μL. The RBD preparations were reduced with 4.5 mM dithiothreitol (DTT) at 56 °C for 30 min and alkylated with 9 mM 2-iodoacetamide in the dark for 30 min at room temperature (RT). The alkylation reactions were then quenched by incubation with DTT (9 mM final concentration) for 15 min at RT. Additional 20 μL of DB was added prior to the proteolytic digest with Pierce MS grade trypsin and Glu-C (overnight at 37 °C, 0.2 and 0.3 μg, respectively). The digested samples were purified using a High Protein and Peptide Recovery Detergent Removal Spin Column (Thermo Fisher Scientific) according to the manufacturer’s instructions. SDC was removed by acidification with 10% trifluoroacetic acid and subsequent centrifugation.

The supernatants were further purified using Pierce peptide desalting spin columns (Thermo Fisher Scientific) according to the manufacturer’s instructions. Each of the purified RBD preparations was divided into three parts: (1) 7.5 μg for nano-LC–MS/MS analysis, (2) 7.5 μg for neuraminidase treatment, and (3) 5 μg for PNGaseF treatment.

For sialic acid removal, RBD preparations were incubated with 1 μL of Sialidase A (GK80040, Agilent, Santa Clara, CA) in 50 μL of provided buffer, overnight at 37 °C. For N-glycan removal, samples were dissolved in 50 μL of 50 mM TEAB and treated with 1 μL of recombinant PNGaseF (Promega, Madison, WI) overnight at 37 °C. All preparations were desalted using Pierce peptide desalting spin columns (Thermo Fisher Scientific) prior to nano-LC–MS/MS analysis.

Nano-LC–MS/MS Analysis of the Recombinant RBD

The RBD proteolytic preparations were analyzed using a QExactive HF mass spectrometer interfaced with an Easy-nLC1200 liquid chromatography system (Thermo Fisher Scientific). Peptides were trapped using an Acclaim Pepmap 100 C18 trap column (100 μm × 2 cm, particle size 5 μm, Thermo Fischer Scientific) and separated with an in-house packed analytical column (75 μm × 300 mm, particle size 3 μm, Reprosil-Pur C18, Dr. Maisch) using a gradient from 7 to 50% of solvent B over 75 min, followed by an increase to 100% of solvent B for 5 min at a flow of 300 nL/min, where solvent A was 0.2% formic acid (FA) and solvent B was 80% acetonitrile in 0.2% FA. The precursor ion mass spectra were acquired in either 600–2000 m/z or 375–1500 m/z ranges at a resolution of 120,000. For nano-LC–MS/MS analysis, the instrument operated in data-dependent mode with the 10 most intense ions with charge states 2 to 5 being selected for fragmentation using higher-energy collision dissociation (HCD). The isolation window was set to 3 m/z and dynamic exclusion to 20 s. MS/MS spectra were recorded at a resolution of 30,000 with the maximum injection time set to 110 ms. To facilitate glycosylated peptide characterization, multiple injections were acquired with precursor detection in the 600–2000 m/z range and different settings for the normalized HCD energies of 22, 28, and 34.

Glycan Database Search and Data Processing

The acquired data were analyzed using Proteome Discoverer version 2.4 (Thermo Fisher Scientific). Database searches were performed with either Byonic (Protein Metrics, Cupertino, CA) or Sequest as search engines. To evaluate the protein preparation purity, the data were initially searched against the custom database consisting of the Uniprot_Chinese hamster_CHO-K1 cell line database (24,147 proteins), SwissProt_human database (20,342 proteins), and the sequence of expressed RBD protein. For later searches aimed to identify the available glycoforms, the raw data acquired with different HCD energies were searched with Proteome Discoverer/Byonic, with Minora Feature Detector node, against the single RBD protein sequence. Precursor mass tolerance was set to 10 ppm and fragment mass tolerance to 30 ppm. Proteolytic peptides with up to two missed cleavages (combined Trypsin Glu-C cleavage sites) were accepted together with variable modification of methionine oxidation and fixed cysteine alkylation. Several different N-glycan databases were used during the data processing.

The initial N-glycan database contained 227 glycan compositions, where 224 were COVID-19-associated N-glycan compositions reported in GlyConnect Compozitor version: 1.0.0 at SIB Swiss Institute of Bioinformatics | Expasy site plus 3 additional mannose-phosphate containing compositions. This database was used for the analysis of neuraminidase-treated samples. The 45 curated nonsialylated compositions, retrieved from the analysis of neuraminidase-treated preparation, were used to create a new glycan database consisting of 153 glycan compositions for the follow-up analysis of native preparations. An O-glycan database consisted of 6 reported, in GlyConnect Compozitor COVID-19 O-glycan, compositions and was used for the analysis of PNGaseF-treated samples.

All glycopeptide identifications were manually evaluated prior to the final assignment of the observed glycosylation forms. The data, acquired with the normalized HCD energy of 22, were used for oxonium ion evaluation to suggest glycan structures for the observed compositions. The extracted ion chromatogram (EIC) peak intensities of the observed glycoforms were used to calculate their relative abundances. For the N-glycopeptide analysis, the relative abundances were calculated using average EIC values from three injections and are expressed as percent of the total signal for all modified and nonmodified forms. For the O-glycopeptide analysis, two injections were acquired for each PNGAseF-treated sample with precursor ions measured at 375–1500 and 600–2000 m/z. The acquisition in the 375–1500 m/z range provided the best detection of O-glycosylated peptides and was used for further data evaluation The O-glycopeptide intensities were calculated based on single injection, and no CV (%) calculations are available for their measurements.

Analysis of RBD Glycosylation

The RBD domain contains two canonical NXS/TN-glycosylation sites, N331 and N343. Proteolytic digest using a combination of Trypsin and Glu-C was selected to access these sites. No glycopeptide enrichment was performed prior to LC–MS/MS, to avoid selective enrichment of specific glycoforms. Instead, to facilitate detection of N-glycopeptides, precursor ions were acquired in the 600–2000 m/z range. The acquired LC–MS/MS data were first evaluated for the presence of oxonium ions to confirm that glycopeptides were observed for each site and to evaluate their potential glycan compositions. Both sites were found to be glycosylated in all three expression systems. The expected HexNAc and HexNAcHex oxonium ions, m/z 204 and 365, were found for the RBD produced in all three cell lines. The NeuAc oxonium ions (m/z 274 and 292) were observed for the RBD expressed in CHO-S and HEK293F but not in Lec3.2.8.1 cell lines. The presence of oligomanose-6-P structures was suggested by the presence of oxonium ions at m/z 243, observed for all three RBD preparations.

The initial data evaluation revealed high heterogeneity at both N-glycosylation sites with possibilities of both multifucosylated and multisialylated structures. Therefore, native RBD preparations from CHO-S and HEK293F preparations were further treated with neuraminidase, prior to the first round of the site-specific glycosylation analysis. The removal of sialic acid served several purposes. (i) It decreased site heterogeneity and thus improved the detection of the existing glycoforms. (ii) It simplified identification of multifucosylated structures. (iii) In addition, it facilitated relative quantification of the site microheterogeneity. In the positive ion mode, the ionization of glycopeptide occurs at the polypeptide chain, thus allowing one to use the observed signal intensities for relative glycoform quantification at the same polypeptide base.[44] The use of glycopeptide signal intensities for both neutral and sialylated glycoforms detected by LC-MS/MS has been shown to have good correlation with quantification carried on AB-labeled glycans.[45] We commonly observed different charge state distributions as well as later retention times for sialylated glycoforms compared to neutral glycoforms, and therefore, we always perform glycan profile evaluation on both native and neuraminidase-treated preparations.

To minimize the effect of ionization efficiencies, we first used LC-MS/MS data acquired for the neuraminidase-treated samples to compare site-specific glycan profiles of the different RBD preparations. This included the evaluation of glycan-type distribution at each site, the antenna distribution within complex glycans, and the degree of fucosylation. Each RBD preparation was analyzed at least three times with identical MS1 settings but different fragmentation energies in MS2 to facilitate glycoform identifications. The EIC peak intensities were used to determine the glycoform abundances. The average values from three injections were used to calculate glycoform abundances expressed as percent of total signal for all modified and nonmodified peptides sharing the same amino acid sequence (Table S7). These data were used to calculate the values for glycan-type distribution at each site, the antenna distribution within complex glycans, and the degree of fucosylation presented in Figure , Tables , 2, and S3. The MS2 data were acquired at different collision energies and were used to evaluate the fucose position and antenna compositions including the suggested presence of LacDiNac. The native RBD preparations were then used to evaluate the degree of sialylation for all hybrid and complex structures observed in neuraminidase-treated preparations. The glycoform abundances were calculated as average of three injections (Table S8) and were used to calculate sialylation levels for each site presented in Tables and S4.

To evaluate the presence of O-glycans, the RBD preparations were treated with PNGaseF. For each sample, two injections were acquired with precursor ions measured at 375–1500 and 600–2000 m/z. The acquisition in the 375–1500 m/z range provided the best detection of O-glycosylated peptides and was used for further data evaluation and calculations presented in Figure and Table S4. The presented values were calculated based on a single injection, and no CV (%) values are available for their measurements.

Deglycosylation of the Recombinant RBD Using Glycosidase Treatment

Removal of N-linked glycans was performed using PNGaseF (New England Biolabs, Ipswich, USA) at a concentration of 125 U/μg protein. Removal of O-linked glycans was performed using O-glycosidase (New England Biolabs, 20,000 U/μg protein), α2-3,6,8 neuraminidase (New England Biolabs, 25 U/μg protein), and α-N-acetyl-galactosaminidase (New England Biolabs, 10 U/μg protein). Removal of both N-linked and O-linked glycans was performed using PNGaseF (New England Biolabs, 125 U/μg protein), O-glycosidase (New England Biolabs, 20,000 U/μg protein), α2-3,6,8 neuraminidase (New England Biolabs, 25 U/μg protein), and α-N-acetyl-galactosaminidase (New England Biolabs, 10 U/μg protein). Removal of sialic acids was performed using the α2-3,6,8 neuraminidase (New England Biolabs, 50 U/μg protein). Removal of fucose was performed using α1-2,4,5,6 fucosidase O (New England Biolabs, 2 U/μg protein) and α1-3,4 fucosidase (New England Biolabs, 4 U/μg protein). All enzymatic reactions were performed as a 1-step reaction with 1× Glycobuffer 2 (New England Biolabs) and 10 μg of RBD produced in CHO-S-, HEK293F-, or Lec3.2.8.1 cells and incubation at 37 °C for 24 h. As heat-treated controls, peptides were incubated at 37 °C for 24 h but without additional enzymes.

Gel Electrophoresis and In-Gel Staining of the Glycosidase-Treated Recombinant RBD

To control efficiency of the enzymatic treatment, 5 μg of the enzyme-treated products or controls were run on a NuPage 4–12% Bis-Tris gel (Invitrogen, Carlsbad, USA) at 100 V for 60 min using an EI9001-XCELL II Mini Cell (Novex, San Diego, CA) together with a Powerease 500 (Novex) and subsequently stained with a SilverQuest Stain kit (Invitrogen) according to the instructions from the manufacturer.

Lectin Blot of the Glycosidase-Treated Recombinant RBD

First, 2 μg of RBD produced in CHO-S or HEK293F cells was treated with neuraminidase, fucosidase, or mock-treated and separated in NuPAGE Bis-Tris 4–12% gels (Invitrogen) using MOPS buffer (Invitrogen) at 100 V. Separated proteins were blotted to the polyvinylidene difluoride membrane (Immobilon-FL 0.45 μm, Merck Millipore, Burlington, MA) using SemiDry Transblot SD (Bio-Rad Laboratories, Hercules, CA). Membranes were blocked with 2% BSA Factor V (Sigma-Aldrich) and 0.1% Tween-20 (VWR chemicals, Radnor, PA) in phosphate-buffered saline (PBS, Medicago, Uppsala, Sweden) (PBS-BSA-T) and then incubated with biotinylated lectin, either 5 μg/mL MALII or 3 μg/mL UEA (both from Vector laboratories, Burlingame, CA), in PBS-BSA-T at 4 °C for 16–20 h. Membranes were then washed three times with PBS + 0.1% Tween-20, incubated with Streptavidin-alkaline phosphatase diluted 1/2000 (Southern Biotech) for 1 h at 20–22 °C, and washed again. Membranes were developed using BCIP/NBT (Sigma-Aldrich) for 1 min.

Levels of Human Anti-SARS-CoV-2 IgG Antibodies in Convalescent Serum Samples

Serum samples from SARS-CoV-2 convalescent individuals (n = 24) were obtained from the department of Clinical Microbiology, Sahlgrenska University Hospital, Gothenburg, Sweden. Samples were collected between 06-03-2020 and 08-27-2020, 25–100 days following a positive PCR-test. Serum was stored at −80 °C until use. The SARS-CoV-2 IgG II Quant assay is a chemiluminescent microparticle immunoassay used for quantitative determination of IgG antibodies to SARS-CoV-2 in human serum and plasma on an ARCHITECT System (Abbott Laboratories, Chicago, IL). The assay measures IgG binding to the RBD of the S-protein. IgG concentrations ≥50 antibody units (AU)/mL were defined as positive.

Viral CPE Neutralization Assay

The titre of neutralizing antibodies against SARS-CoV-2 in the patient sera was determined against the DE-Gbg20 viral strain (NCBI GenBank ID: MW092768) at a titre of 10–6. A 50% tissue culture infectious dose (TCID50) assay was performed as defined by Reed and Muench.[46] All sera were heat inactivated at 56 °C for 30 min before twofold serial dilution in serum-free Dulbecco’s Modified Eagle Medium (DMEM) with 100TCID50 DEGbg20 followed by incubation for 2 h at 37 °C. The virus–antibody mixture was added to a monolayer of VERO CCL-81 cells grown in 96 well-plates in DMEM supplemented with 2% penicillin–streptomycin and 2% fetal calf serum. The plates were incubated for 72 h at 37 °C with 5% CO2. The neutralizing titre for each serum was defined at the highest serum dilution at which 50% of the added virus was neutralized.

Anti-SARS-CoV-2 Antibody Reactivity Assay

The antibody reactivities toward glycosidase-treated proteins were assessed using an enzyme-linked immunosorbent assay (ELISA). Briefly, Nunc Maxisorp 96-well plates (Thermo Fischer Scientific) were coated with 0.1 μg of glycosidase-treated peptides or heat-treated controls diluted in carbonate buffer (pH 9.6). Coating was performed overnight at 4 °C followed by washing three times with 0.05% tween20 in PBS. The plates were blocked in 2% milk for 30 min at RT prior to addition of sera (diluted 1:100 in 1% milk in PBS with 0.05% tween20) and 1.5 h incubation at 37 °C. The plates were washed three times before addition of alkaline phosphatase-conjugated goat antihuman IgG (Jackson ImmunoResearch, Cambridgeshire, UK) diluted 1:1000 in 1% milk in PBS with 0.05% tween. After 1.5 h incubation at 37 °C, the plate was washed six times and 1 mg/mL p-nitrophenylphosphate (Medicago, Danmarks-Berga, Sweden) in diethanolamine substrate buffer was added. The plates were incubated in the dark for 30 min before spectrophotometric measurement at 405 nm.

Statistics

For the comparison of antibody reactivity as determined by ELISA, the Wilcoxon matched-pair signed rank test was used. The comparison between anti-RBD IgG levels and antibody reactivity toward recombinant RBD was done with Pearson correlation coefficients, assuming normal distribution. All statistical analyses were performed using the Graphpad Prism software version 9.3.1 (GraphPad Software Inc., San Diego, CA, USA).

Ethical Statement

The study was approved by the ethical review board in Gothenburg (Dnr: 2021-02252).