Vulnerabilities in coronavirus glycan shields despite extensive glycosylation

Severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) coronaviruses (CoVs) are zoonotic pathogens with high fatality rates and pandemic potential. Vaccine development focuses on the principal target of the neutralizing humoral immune response, the spike (S) glycoprotein. Coronavirus S proteins are extensively glycosylated, encoding around 66–87 N-linked glycosylation sites per trimeric spike. Here, we reveal a specific area of high glycan density on MERS S that results in the formation of oligomannose-type glycan clusters, which were absent on SARS and HKU1 CoVs. We provide a comparison of the global glycan density of coronavirus spikes with other viral proteins including HIV-1 envelope, Lassa virus glycoprotein complex, and influenza hemagglutinin, where glycosylation plays a known role in shielding immunogenic epitopes. Overall, our data reveal how organisation of glycosylation across class I viral fusion proteins influence not only individual glycan compositions but also the immunological pressure across the protein surface.

C oronaviruses (CoVs) are enveloped pathogens responsible for multiple respiratory disorders of varying severity in humans 1 . Certain CoVs represent a significant threat to global human health, as illustrated by outbreaks of severe acute respiratory syndrome coronavirus (SARS-CoV) in 2003 2 , Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012 3 , and most recently of SARS-CoV-2 4 . Given their mortality rates, the current lack of targeted treatments and licensed vaccines, and their capacity to transmit between humans and across species barriers 5,6 , there is an urgent need for effective countermeasures to combat these pathogens. Ongoing vaccine development efforts focus on the spike (S) proteins that protrude from the viral envelope and constitute the main target of neutralizing antibodies 7,8 .
Previous analyses of viral glycan shields have revealed the presence of underprocessed oligomannose-type glycans that seemingly arise due to steric constraints that prevent access of glycan processing enzymes to substrate glycans 24,27,28 , especially when the viral glycoprotein has evolved to mask immunogenic epitopes with a particularly dense array of host-derived glycans 26,[29][30][31][32][33][34] . Restricted access to these glycan sites or interference with surrounding protein surface or neighbouring glycan residues can render glycan processing enzymes ineffective in specific regions 27,28,35 . Glycan processing on soluble glycoproteins has also been shown to be a strong reporter of native-like protein architecture and thus immunogen integrity [36][37][38] ; and glycan processing on a successful immunogen candidate should therefore mimic, as closely as possible, the structural features observed on the native virus 39,40 .
Here, we provide global and site-specific analyses of N-linked glycosylation on soluble SARS, MERS and HKU1 CoV S glycoproteins and reveal extensive heterogeneity, ranging from oligomannose-type glycans to highly-processed complex-type glycosylation. The structural mapping of glycans of trimeric S proteins revealed that some of these glycans contribute to the formation of a cluster of oligomannose-type glycans at specific regions of high glycan density on MERS-CoV S. Molecular evolution analysis of SARS and MERS S genes also reveals a higher incidence of amino-acid diversity on the exposed surfaces of the S proteins that are not occluded by N-linked glycans. In addition, we compare the structures of the respective glycan coats of SARS and HIV-1 envelope proteins using cryo-electron microscopy (cryo-EM) and computational modelling, which delineate a sparse glycan shield exhibited on SARS S compared with other viral glycoproteins. We therefore undertook a comparative analysis of viral glycan shields from characterized class I fusion proteins to highlight how glycosylation density influences oligomannose-type glycan abundance, and the relationship between effective glycan shields and viral evasion ability.
Together, these data underscore the importance of glycosylation in viral immune evasion.

Results and discussion
Glycan processing of trimeric SARS and MERS spike proteins. To generate a soluble mimic of the viral S proteins, we used the 2P-stabilised native-like SARS and MERS S protein antigens, the design and structures of which have been described previously by Pallesen et al. 41 . SARS, MERS and HKU1 S genes encode many N-linked glycan sequons; 23, 23 and 29, respectively (Fig. 1a). We initially sought to quantitatively assess the composition of the carbohydrate structures displayed on the S glycoproteins. Nlinked glycans were enzymatically released, fluorescently labelled, and subjected to hydrophilic interaction chromatography-ultraperformance liquid chromatography (HILIC-UPLC). Treatment with endoglycosidase H (Endo H) revealed a population (SARS 32.2%; MERS 33.8%, HKU1 25.0%) of underprocessed oligomannose-type glycans (Fig. 1b). This observation of both complex and oligomannose-type glycans reveals that the majority of N-linked glycans can be processed, although there is limited processing at specific sites across the S proteins. It is also interesting to note that the distribution of oligomannose-type glycans was broad, with Man 5 GlcNAc 2 to Man 9 GlcNAc 2 glycans all present, without one particular dominant peak, as is the case for some viral glycoproteins, such as HIV-1 Env 36 . The proportion of oligomannose-type glycans on recombinant coronavirus S proteins is consistent with previous studies performed on virally derived MERS and SARS coronavirus S proteins 17,42 . Coronaviruses have been previously been reported to form virions by budding into the lumen of endoplasmic reticulum-Golgi intermediate compartments (ERGIC) 43,44 . Observations of hybridand complex-type glycans on virally derived material 17,42 would, however, suggest that it is likely that coronavirus virions travel through the Golgi apparatus after virion formation in the ERGIC en route to the cell surface, thus supporting recombinant immunogens as models of viral glycoproteins.
To ascertain the precise structures of N-linked glycans, glycan pools of each coronavirus S protein were analysed by negative-ion ion-mobility-electrospray ionisation mass spectrometry (IM-ESI MS) ( Supplementary Fig. 1). Consistent with the UPLC data, IM-ESI MS confirmed an array of complex-type glycans ranging from mono-to tetra-antennary, but also oligomannose-and hybridtype glycans. The glycan compositions characterised in the spectra were largely invariant among the coronaviruses with no major structural differences observed.
Clustering of underprocessed glycans on MERS S. We subsequently performed glycopeptide analysis to ascertain the compositions of glycans at all of the potential N-linked glycosylation sites (PNGs). MERS, SARS and HKU1 recombinant S proteins were reduced, alkylated and digested with an assortment of proteases to yield glycopeptides, which were subjected to in-line liquid chromatography-mass spectrometry (LC-MS). This revealed differential levels of oligomannose, hybrid, and complextype glycan populations (Fig. 2a, b). Using structures of the trimeric MERS and SARS S proteins (PDB ID: 5X59 and 5X58, respectively), we generated models of fully glycosylated coronavirus spikes using experimentally determined glycan compositions (Fig. 3a, b). This revealed that oligomannose-type glycans on MERS S co-localize to specific clusters on the head of the S protein, consisting of glycans at Asn155, Asn166, and Asn236 (Fig. 3a). We hypothesized that the fully oligomannose-type glycan population in this cluster arises due to the hindered accessibility of glycan processing enzymes to access the substrate glycan 28 . As such, we performed mutagenesis to knock out glycosylation sites with N155A, N166A, and N236A mutations. Site-specific analysis of these glycan-KO mutants revealed enhanced trimming of mannose residues, i.e. increased processing, when glycan clustering was reduced (SI Fig. 4). The presence of clustered oligomannose-type glycans is reminiscent of that found on other viral glycoproteins, including HIV-1 Env and LASV GPC 24,31,34,36,45,46 .
Interestingly, SARS and HKU1 (SI Fig. 2) S proteins did not exhibit specific mannose clusters that contribute to the overall mannose abundance, but only isolated glycans were underprocessed. We speculate that the oligomannose-type glycans here arise from protein-directed inhibition of glycan processing, as opposed to the glycan-influenced processing observed on MERS. Importantly, oligomannose-type glycans has also been implicated in innate immune recognition of coronaviruses by lectins 47,48 that recognise these underprocessed glycans as pathogen-associated molecular patterns.
Given that the receptor-binding domain is the main target of neutralising antibodies 8 , it is surprising that the DPP4 receptorbinding site of MERS S was not occluded by glycans (Fig. 3a), as observed for other receptor-binding sites of class I viral fusion proteins, including SARS S (Fig. 3b), HIV-1 Env 49 , LASV GPC 24 and influenza HA 50 . We suggest that this is likely due to the intrinsic functionality of the receptor-binding domain of MERS S, that would be sterically hindered by the presence of N-linked glycans, whereas other viruses are able to accommodate the posttranslational modifications, without greatly perturbing functionality.
Sequence diversification of CoV spikes. We hypothesized that solvent-accessible, amino-acid residues on S proteins would be undergoing higher rates of mutations compared with buried residues and regions that are occluded by glycans, which are unable to be targeted by host immune responses. To that end, we performed an evaluation of amino-acid diversification on a residue-specific level, using publicly available gene sequences of SARS and MERS S, which was calculated as the number of observed pairwise differences divided by the total number of pairwise comparisons. Firstly, we found that amino-acid diversity was elevated at known epitopes targeted by neutralizing antibodies, such as the N-terminal domain and the receptor-binding domains, and reduced in the regions in the S2 domain, such as the fusion peptide, heptad repeat one, and the central helix domains, which are likely subject to greater functional constraints (Fig. 4a).
Analysis of the relative ratio of non-synonymous to synonymous nucleotide substitutions (i.e. dN/dS ratios) revealed that exposed residues exhibited significantly higher dN/dS values (Fig. 4b). Buried residues on SARS had mean dN/dS ratios of 0.31 compared with 2.82 for exposed resides. Likewise, the buried residues on MERS had a calculated dN/dS ratio of 0.10 compared with exposed residues with a value of 0.45. Furthermore, when per-site amino-acid diversities were mapped onto the fully glycosylated structural model of the respective CoV S proteins ( Fig. 4c), hotspots of mutations were highlighted on the protein surface throughout the trimer revealing extensive vulnerabilities permeating through the glycan shield of SARS and MERS CoVs. It is interesting to note the lack of amino-acid diversity on the receptor-binding domains of MERS S proteins that protrude away from the glycans. We would suggest that this may result from the intrinsic receptor-binding functionality of these domains.
Although dN/dS estimates are comparable within each viral outbreak, they are not directly comparable between viral families as they can only be considered in the environment in which they are measured (i.e. multiple differences in transmission ecology and host-virus interactions disallow meaningful comparisons). For example, differences in the epidemic behaviour and host immune environment of MERS and SARS outbreaks likely contribute to the observed genetic diversity and thus dN/dS. MERS was characterized by repeated spillover events from camels into humans, where it circulated transiently. In contrast, the SARS outbreak corresponded to a single zoonotic event followed by extensive human-to-human transmission. Consequently, inferring the degree of selection acting upon MERS and SARS from dN/dS analysis is extremely difficult. Importantly, while similar analyses of SARS-CoV-2 is desirable, due to the low genetic variation among the current SARS-CoV-2 sequences (as of 17 March 2020), which likely include deleterious mutations a b S1 S1 S2 S1/S2 S2' S1/S2 S2' S1 S2 S1/S2 S2' Oligomannose-type glycans (M5 to M9; Man 5 GlcNAc 2 -Man 9 GlcNAc 2 ) (green) and complex-type glycans (magenta) were identified by Endo H digestion, with quantification of major glycan types summarised as a pie chart. Oligomannose-type glycans are schematically annotated with mannose residues as green circles and GlcNAc residues as blue squares.
Visualising the HIV-1 and SARS glycan shields by cryo-EM. HIV-1 Env is a prototypic viral class I fusion protein that exhibits extensive surface glycosylation, resulting in an effective glycan shield to aid evasion from the host adaptive immune response 21,31 . In order to visualize the structure of the respective glycan "shields" of HIV-1 and SARS coronavirus we used singleparticle cryo-electron microscopy (cryo-EM). The results for HIV-1 Env were reproduced directly from Berndsen et al. 51 while the previously published SARS 2P dataset 52 was reprocessed for this study. Although cryo-EM datasets of fully glycosylated MERS S 41 and chimpanzee simian immunodeficiency virus (SIVcpz) 53 are also available, only the HIV and SARS data were of sufficient quality (Fig. 5). We recently showed 51 that dynamics in surface exposed glycans HIV-1 Env leads to an extensive network of interactions that drive higher-order structuring in the glycan shield. This structure defines diffuse boundaries between buried and exposed surface protein surface, which can serve to define potential sites of vulnerability. Cryo-EM captures the ensembleaverage structure of biomolecules and therefore glycan dynamics results in blurred density at the resolutions necessary for building atomic structure. However, we showed how a simple combination of low-pass filtering and auto-thresholding, as well as 3D variability analysis, can reveal the previously hidden structure of the SARS glycan shield and compare it with the HIV-1 Env glycan shield 51 (Fig. 5). We observe the nearly all-encompassing glycan density on HIV-1 Env and evidence for extensive glycan-glycan interactions, especially in the oligomannose patch regions, whereas the glycans on SARS S appear more isolated and lack the wide-ranging glycan networks that are the hallmark of an effective glycan shield 54,55 . The 3D variability maps are more sensitive to low intensity signal and reveal additional glycan-glycan interactions in both maps, however the S1 receptor-binding domains in the SARS dataset were shown to exist in both up and down conformations 52 , leading to poor resolution and significant 2D-variability which is convolved with the variability coming from glycans and limits the interpretability of glycan shielding effects in this region of the map.
Disparate shielding efficacies of viral glycosylation. Viral envelope proteins are glycosylated to varying degrees, but depending on their overall mass, surface area, and volume, the overall density of glycan shielding may differ significantly. For example, both LASV GPC and coronavirus S proteins consist of 25% glycan by molecular weight. However, given the significantly larger protein surface area and volume of coronavirus S proteins, coverage of the glycan "shield" over the proteinaceous surface is considerably sparser in comparison to the smaller LASV GPC, which occludes a far greater proportion of the protein surface with fewer glycans. To demonstrate that the presence of  glycosylation plays a major role in the immune response to these different glycoproteins, we studied the glycome of several biomedically important coronaviruses and compared their glycan compositions in a structural context. We then investigated the glycan shield densities of seven viral class I fusion proteins using a global structural approach which was calculated by dividing the number of amino-acids that interact with glycans by the number of solvent-accessible aminoacid residues of each respective glycoprotein and plotted this against oligomannose abundance. A strong correlation was observed (Fig. 6) and viruses historically classified as "evasion strong" 56 had significantly elevated glycan shield densities and oligomannose abundance, which underscores the importance of glycan shielding in immune evasion.
Whether the restricted glycan shielding observed on coronaviruses is linked to the zoonosis of the pathogens is unknown. However, it is tempting to speculate, for example, that MERS has not evolved a dense shield since it would not offer as much of a protective advantage against camel nanobodies (also known as single-domain antibodies) which could more easily penetrate it. Investigation of the host immune response to viruses in their natural reservoirs may offer a route to understanding why coronavirus glycosylation does not reach the density of other viruses such as HIV-1. In addition, it may be that functional constraints, such as maintaining flexibility of the receptorbinding domains, limit the accretion of glycans on coronavirus spikes, which would render it incapable of performing its primary functions, including receptor-binding and membrane fusion. This phenomenon has been observed on other viral glycoproteins, including influenza HAs, where there is a limit to the accumulation of glycosylation sites that can be incorporated in vivo 57,58 , compared with in vitro 59 , with H3N2 and H1N1 HAs replacing existing PNGs rather than continually adding them upon the glycoprotein 22,58 . The importance of glycosylation in modulation of viral infectivity and immune responses have also been investigated during influenza vaccine research 22,60 and should be considered in coronavirus vaccine research.
More topically, it is interesting to note the conservation of Nlinked glycosylation sites on S proteins from the SARS-CoV-2 and SARS (SI Fig. 6). SARS-CoV-2 possesses a total of 22 Nlinked glycan sites compared with 23 on SARS, with 18 of these sites being in common. As such, it is likely that these glycans on this novel coronavirus would shield similar immunogenic epitopes that are observed on SARS S. As expected, most of the differences between the two viruses are observed on the S1 subunit, due to its amenability to substitutions while still remaining functionally competent. Furthermore, likely targets for the majority of antibodies targeting the spike are located on S1, resulting in greater levels of immune pressure upon this subunit. This notion is further reflected in terms of glycosylation, with all of the glycan sites conserved on the S2 subunit between SARS and SARS-CoV-2, whereas the S1 subunit exhibits glycan site additions and deletions (SI Fig. 7). Bioinformatic analysis of current SARS-CoV-2 spike genes (n = 566 as of 17 March 2020)  11 and SARS S (PDB ID 5X58) 11 glycoproteins. The glycans are colored according to oligomannose content, as defined by the key. DPP4 receptor-binding sites and ACE2 receptor-binding sites for MERS and SARS, respectively, are indicated in light blue. The S1 and S2 subunits are colored light grey and dark grey, respectively. NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-020-16567-0 ARTICLE from nextstrain 61 (https://nextstrain.org/ncov) revealed low sequence diversity and no changes in glycosylation sites (SI Fig. 8).
Although it is difficult to directly compare viruses in terms of immunogenic responses, on the one hand, SARS and MERS coronaviruses readily elicit neutralizing antibodies following infection or immunization [62][63][64][65] . Indeed, many potential MERS-CoV vaccine candidates are able to elicit high titres of serum IgG upon immunization but fail to produce sufficient mucosal immunity 65 . In contrast, the high mutation rate 66 and the evolving glycan shield of HIV-1 20,39 , which firmly exemplifies it as "evasion strong" virus, hinders the development of broadly neutralizing antibodies 67 .Viruses classified as "evasion strong" 26,56 may then differ due to varied efficacies of protein surface shielding by glycans.
Overall, this study adds further evidence suggesting that extensive N-linked glycan modifications of SARS and MERS CoV S proteins do not constitute an effective shield, in comparison to glycan shields of certain other viruses, which is reflected by the overall structure, density and oligomannose abundance across the corresponding trimeric glycoproteins. We also demonstrate that amino-acid diversification indeed occurs at antibody accessible c Mapping of the per residue amino-acid diversity shown in A onto the structures of SARS and MERS S (PDB ID 5X58 and 5X59, respectively) 11 . S proteins are presented as backbone traces with residues colored according to amino-acid diversity. Residues with elevated diversity are colored in red, and N-linked glycans are presented as white surfaces.
regions on the trimer, which confirms that glycans play a role in occluding specific regions if vulnerability on the glycoprotein. Furthermore, comparisons between glycan shields from a number of viruses highlight the importance of a glycan shield in immune evasion and reveal structural principles that govern glycosylation status.

Methods
Expression and purification of coronavirus spike glycoproteins. Human embryonic kidney 293 Freestyle (HEK293F) cells were transfected with mammalian-codon-optimised genes encoding 2P-stabilised SARS MERS and HKU1 S proteins containing a C-terminal T4 fibritin trimerization domain, an HRV3C cleavage site, an 8xHis-tag and a Twin-Strep-tag 41 . H3N2 Victoria 2011 hemagglutinin was also expressed in the HEK293F cells. The 200 ml cultures were harvested 6 days after transfection, filtered and purified by nickel-affinity chromatography and size exclusion chromatography using a Superdex TM 16/600 75 pg column (GE Healthcare).
Release and labelling of N-linked glycans. Excised coronavirus S gel bands were washed alternately with acetonitrile and water before drying in a vacuum centrifuge. The bands were rehydrated with 100 μL of water and incubated with PNGase F at 37°C overnight. Aliquots of released N-linked glycans were also fluorescently labelled with procainamide, by adding 100 μL of labelling mixture (110 mg/mL procainamide and 60 mg/mL sodium cyanoborohydrate in 70% DMSO and 30% glacial acetic acid) and incubating for 4 h at 65°C. Procainamide labelled glycans were purified using Spe-ed Amide 2 columns (Applied Separations).  4) and B was acetonitrile. Quantification of oligomannose-type glycans was achieved by digestion of fluorescently labelled glycans with Endo H, and clean-up using a PVDF protein-binding membrane (Millipore). Empower 3 software (Waters) was used for data processing.
Mass spectrometry of glycans. Prior to ion-mobility electrospray ionisation MS and tandem MS analysis, PNGase F released N-linked glycans were purified on a Nafion ® 117 membrane (Sigma-Aldrich) and a trace amount of ammonium phosphate was added to promote phosphate adduct formation. Glycans were analyzed by direct infusion using a Synapt G2Si instrument (Waters) with the following settings: capillary voltage, 0. Mass spectrometry of glycopeptides. Aliquots of 30-50 μg of coronavirus spikes were denatured, reduced and alkylated as described previously 36 . Proteins were proteolytically digested with trypsin (Promega), chymotrypsin (Promega), alphalytic protease (Sigma-Aldrich) and Glu-C (Promega). Reaction mixtures were dried and peptides/glycopeptides were extracted using C18 Zip-tip (MerckMilipore) following the manufacturer's protocol. Samples were resuspended in 0.1% formic acid prior to analysis by liquid chromatography-mass spectrometry using an Easy-nLC 1200 system coupled to an Orbitrap Fusion mass spectrometer (Thermo Fisher Scientific

SARS-CoV S 2P
HIV Env BG505 SOSIP.664 Low-pass filtered cryo-EM map Sparx 3D variability map Precursor and fragmentation detection were performed using an Orbitrap at the following resolution: MS1 = 100,000 and MS2 = 30,000. The automatic gain control (AGC) targets were MS1 = 4e 5 and MS2 = 5e 4 , and injection times were MS1 = 50 and MS2 = 54. The following cleavage sites were used for the respective proteases; trypsin=R/K, chymotrypsin=F/Y/W, alpha lytic protease=T/A/S/V, Glu C=E/D. Number of missed cleavages were set at 3. The following modifications were also included: Carbamidomethyl (+57.021464, target=C, fine control=fixed), Oxidation (+15.994915, target=M, fine control=variable rare 1), Glu to pyro-Glu (−18.010565, target=peptide N-term E, fine control=variable rare 1), and Gln to pyro-Glu (−17.026549, target peptide N-term Q, fine control=variable rare 1). Glycopeptide fragmentation data were extracted form raw files using Byonic TM (Version 3.5.0) and Byologic TM (Version 3.5-15; Protein Metrics Inc.). Glycopeptide fragmentation data were manually evaluated with true-positive assignments given when correct b-and y-fragments and oxonium ions corresponding to the peptide and glycan, respectively, were observed. The precursor mass tolerance was set at 4 ppm for precursor ions and 10 ppm for fragment ions. MS data were searched using a glycan library (SI Fig. 9) with the identical peptide sequence. A 1% false discovery rate (FDR) was applied. The extracted ion chromatographic areas for each true-positive glycopeptide, with the same amino-acid sequence, were compared to determine the relative quantitation of glycoforms at each specific Nlinked glycan site.
Model construction. Structural models of N-linked glycan presentation on SARS, MERS and HKU1 S were created using electron microscopy structures (PDB ID 5X58, 5X59, and 5I08, respectively) 9,11 , along with complex-, hybrid-, and oligomannose-type N-linked glycans (PDB ID 4BYH, 4B7I, and 2WAH  72 and a codon-structured substitution model 73 . Multiple independent MCMC runs of 10-20 million steps were executed to ensure that stationarity and convergence had been achieved. Empirical distributions of time-scaled phylogenies were obtained by combining (after the removal of burnin) the posterior tree distributions from the separate runs, which were subsequently used to estimate dN/dS ratios using the renaissance counting approach 74,75 implemented in BEAST v 1.8.4. We also estimated per-site amino-acid diversity, which was calculated as the average number of amino-acid difference between two sequences at an amino-acid position in all possible pairs in the sequence alignment.
Cryo-EM data analysis and visualization. Single-particle cryo-EM data analysis of BG505 SOSIP.664 in complex with RM20A3 Fab was reproduced directly from Berndsen et al. 51 . Data for the SARS-CoV S 2P ectodomain was previously published 52 and the final particle stack and alignment parameters from the published reconstruction were used for 3D variability analysis in the SPARX software package 76,77 . All metadata for these reconstructions along with raw data images and FSC resolution curves can be found in the original publications. In summary, both datasets were acquired on a FEI Titan Krios (Thermo Fisher) operating at 300 KeV equipped with a K2 Summit Direct Electron Detector (Gatan). Movie micrographs were aligned and dose weighted with MotionCor2 78 and CTF estimation was performed with Gctf 79 . Single-particle data processing was performed using CryoSparc v.2 80 and Relion v.3 81 . Maps were low-pass filtered using a Gaussian kernel and visualized in UCSF chimera 82 . A detailed description of the auto-thresholding method used to set the isosurface value for visualisation of lowpass filtered maps can be found in Berndsen et al. 51 .
Clustering analysis of viral glycan shields. Solvent-accessible residues and interactions between N-linked glycans and amino-acid residues were calculated using Proteins, Interfaces, Structures and Assemblies (PISA) European Bioinformatics Institute (EBI) 83 . Glycan shield density was calculated by the number of amino-acid residues interacting with glycans divided by the total number of solvent-accessible amino-acid residues.
Reporting summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
In addition to the data reported in this paper and accompanying supplementary materials, the raw mass spectrometric data that support the findings of this study  5). The number of amino-acid residues interacting with N-linked glycans was divided by the number of solvent-accessible amino-acid residues of the glycoprotein as a measure for global glycan shield density. All viral glycoproteins analysed were expressed as trimers in HEK293F cells apart from LASV GPC, which was derived from virus-like particles from Madin-Darby canine kidney II cells.