Position-specific N- and O-glycosylation of the reactive center loop impacts neutrophil elastase-mediated proteolysis of corticosteroid-binding globulin.

Corticosteroid-binding globulin (CBG) delivers anti-inflammatory cortisol to inflamed tissues through proteolysis of an exposed reactive center loop (RCL) by neutrophil elastase (NE). We previously demonstrated that RCL-localized Asn347-linked N-glycans impact NE proteolysis, but a comprehensive structure-function characterization of the RCL glycosylation is still required to better understand CBG glycobiology. Herein, we first performed RCL-centric glycoprofiling of serum-derived CBG to elucidate the Asn347-glycans and then used molecular dynamics simulations to study their impact on NE proteolysis. Importantly, we also identified O-glycosylation (di/sialyl T) across four RCL sites (Thr338/Thr342/Thr345/Ser350) of serum CBG close to the NE-targeted Val344-Thr345 cleavage site. A restricted N- and O-glycan co-occurrence pattern on the RCL involving exclusively Asn347 and Thr338 glycosylation was experimentally observed and supported in silico by modeling of a CBG-GalNAc-transferase (GalNAc-T) complex with various RCL glycans. GalNAc-T2 and GalNAc-T3 abundantly expressed by liver and gall bladder, respectively, showed in vitro a capacity to transfer GalNAc (Tn) to multiple RCL sites suggesting their involvement in RCL O-glycosylation. Recombinant CBG was then used to determine roles of RCL O-glycosylation through longitudinal NE-centric proteolysis experiments, which demonstrated that both sialoglycans (disialyl T) and asialoglycans (T) decorating Thr345 inhibit NE proteolysis. Synthetic RCL O-glycopeptides expanded on these findings by showing that Thr345-Tn and Thr342-Tn confer strong and moderate protection against NE cleavage, respectively. Molecular dynamics substantiated that short Thr345-linked O-glycans abrogate NE interactions. In conclusion, we report on biologically relevant CBG RCL glycosylation events, which improve our understanding of mechanisms governing cortisol delivery to inflamed tissues.

Sensitive and Specific Global Cell Surface N-Glycoproteomics Shows Profound Differences Between Glycosylation Sites and Subcellular Components.

Cell surface glycans are essential for establishing cell communication, adhesion, and migration. However, it remains challenging to obtain cell surface-specific information about glycoconjugate structures. Acquiring this information is essential for unraveling the functional role of glycans and for exploiting them as clinical targets. To specifically analyze the N-glycoprotein forms expressed at the cell surface, we developed a C18 liquid chromatography (LC)-mass spectrometry (MS)-based glycoproteomics method in combination with highly specific cell surface protein labeling and enrichment using a biotin label. The surface-specificity of the method was validated by MS-based proteomics of subcellular component marker proteins. Using the human keratinocytes N/TERT-1 as a model system, we identified and quantified the glycosylation of hundreds of cell surface N-glycosylation sites. This approach allowed us to study the glycoforms present at the functional relevant cell surface, omitting immaturely glycosylated proteins present in the secretory pathway. Interestingly, the different stages of N-glycan processing at individual sites displayed at the cell surface were found to correlate with their accessibility for ER-residing processing enzymes, as investigated through molecular dynamics simulations. Using the new approach, we compared N-glycosylation sites of proteins expressed on the cell surface to their counterparts in a total cell lysate, showing profound differences in glycosylation between the subcellular components and indicating the relevance of the method for future studies in understanding contextual glycan functions.

Hyper-truncated Asn355- and Asn391-glycans modulate the activity of neutrophil granule myeloperoxidase.

Myeloperoxidase (MPO) plays essential roles in neutrophil-mediated immunity via the generation of reactive oxidation products. Complex carbohydrates decorate MPO at discrete sites, but their functional relevance remains elusive. To this end, we have characterised the structure-biosynthesis-activity relationship of neutrophil MPO (nMPO). Mass spectrometry demonstrated that nMPO carries both characteristic under-processed and hyper-truncated glycans. Occlusion of the Asn355/Asn391-glycosylation sites and the Asn323-/Asn483-glycans, located in the MPO dimerisation zone, was found to affect the local glycan processing, thereby providing a molecular basis of the site-specific nMPO glycosylation. Native mass spectrometry, mass photometry and glycopeptide profiling revealed significant molecular complexity of diprotomeric nMPO arising from heterogeneous glycosylation, oxidation, chlorination and polypeptide truncation variants and a previously unreported low-abundance monoprotomer. Longitudinal profiling of maturing, mature, granule-separated and pathogen-stimulated neutrophils demonstrated that nMPO is dynamically expressed during granulopoiesis, unevenly distributed across granules and degranulated upon activation. We also show that proMPO-to-MPO maturation occurs during early/mid-stage granulopoiesis. While similar global MPO glycosylation was observed across conditions, the conserved Asn355-/Asn391-sites displayed elevated glycan hyper-truncation, which correlated with higher enzyme activities of MPO in distinct granule populations. Enzymatic trimming of the Asn355-/Asn391-glycans recapitulated the activity gain and showed that nMPO carrying hyper-truncated glycans at these positions exhibits increased thermal stability, polypeptide accessibility and ceruloplasmin-mediated inhibition potential relative to native nMPO. Finally, molecular modelling revealed that hyper-truncated Asn355-glycans positioned in the MPO-ceruloplasmin interface are critical for uninterrupted inhibition. Here, through an innovative and comprehensive approach, we report novel functional roles of MPO glycans, providing new insight into neutrophil-mediated immunity.

PD-L1 Glycosylation and Its Impact on Binding to Clinical Antibodies.

Immune checkpoint inhibitors, including PD-L1/PD-1, are key regulators of the immune response and promising targets in cancer immunotherapy. N-glycosylation of PD-L1 affects its interaction with PD-1, but little is known about the distribution of glycoforms at its four NXS/T sequons. We optimized LC-MS/MS methods using collision energy modulation for the site-specific resolution of specific glycan motifs. We demonstrate that PD-L1 on the surface of breast cancer cell line carries mostly complex glycans with a high proportion of polyLacNAc structures at the N219 sequon. Contrary to the full-length protein, the secreted form of PD-L1 expressed in breast MDA-MB-231 or HEK293 cells demonstrated minimum N219 occupancy and low contribution of the polyLacNAc structures. Molecular modeling of PD-L1/PD-1 interaction with N-glycans suggests that glycans at the N219 site of PD-L1 and N74 and N116 of PD-1 may be involved in glycan-glycan interactions, but the impact of this potential interaction on the protein function remains at this point unknown. The interaction of PD-L1 with clinical antibodies is also affected by glycosylation. In conclusion, PD-L1 expressed in the MDA-MB-231 breast cancer cell line carries polyLacNAc glycans mostly at the N219 sequon, which displays the highest variability in occupancy and is most likely to influence the interaction with PD-1.

O-glycosylation on cerebrospinal fluid and plasma apolipoprotein E differs in the lipid-binding domain.

The O-glycoprotein apolipoprotein E (APOE), the strongest genetic risk factor for Alzheimer's disease, associates with lipoproteins. Cerebrospinal fluid (CSF) APOE binds only high-density lipoproteins (HDLs), while plasma APOE attaches to lipoproteins of diverse sizes with binding fine-tuned by the C-terminal loop. To better understand the O-glycosylation on this critical molecule and differences across tissues, we analyzed the O-glycosylation on APOE isolated from the plasma and CSF of aged individuals. Detailed LC-MS/MS analyses allowed the identification of the glycosite and the attached glycan and site occupancy for all detectable glycosites on APOE and further three-dimensional modeling of physiological glycoforms of APOE. APOE is O-glycosylated at several sites: Thr8, Thr18, Thr194, Ser197, Thr289, Ser290 and Ser296. Plasma APOE held more abundant (20.5%) N-terminal (Thr8) sialylated core 1 (Neu5Acα2-3Galß1-3GalNAcα1-) glycosylation compared to CSF APOE (0.1%). APOE was hinge domain glycosylated (Thr194 and Ser197) in both CSF (27.3%) and plasma (10.3%). CSF APOE held almost 10-fold more abundant C-terminal (Thr289, Ser290 and Ser296) glycosylation (36.8% of CSF peptide283-299 was glycosylated, 3.8% of plasma peptide283-299), with sialylated and disialylated (Neu5Acα2-3Galß1-3(Neu5Acα2-6) GalNAcα1-) core 1 structures. Modeling suggested that C-terminal glycosylation, particularly the branched disialylated structure, could interact across domains including the receptor-binding domain. These data, although limited by sample size, suggest that there are tissue-specific APOE glycoforms. Sialylated glycans, previously shown to improve HDL binding, are more abundant on the lipid-binding domain of CSF APOE and reduced in plasma APOE. This indicates that APOE glycosylation may be implicated in lipoprotein-binding flexibility.

Rapidly Display Glycan Symbols in 3D Structures: 3D-SNFG in LiteMol.

The representation of carbohydrates in 3D space using symbols is a powerful visualization method, but such representations are lacking in currently available visualization software. The work presented here allows researchers to display carbohydrate 3D structures as 3D-SNFG symbols using LiteMol from a web browser (e.g., v.litemol.org/?loadFromCS=5T3X ). Any PDB ID can be substituted at the end of the URL. Alternatively, the user may enter a PDB ID or upload a structure. LiteMol is available at https://v.litemol.org and automatically depicts any carbohydrate residues as 3D-SNFG symbols. To embed LiteMol in a webpage, visit https://github.com/dsehnal/LiteMol .

Galectin binding to cells and glycoproteins with genetically modified glycosylation reveals galectin-glycan specificities in a natural context.

Galectins compose a protein family defined by a conserved sequence motif conferring affinity for ß-galactose-containing glycans. Moreover, galectins gain higher affinity and fine-tune specificity by glycan interactions at sites adjacent to their ß-galactoside-binding site, as revealed by extensive testing against panels of purified glycans. However, in cells, galectins bind glycans on glycoproteins and glycolipids in the context of other cellular components, such as at the cell surface. Because of difficulties in characterizing natural cellular environments, we currently lack a detailed understanding of galectin-binding specificities in the cellular context. To address this challenge, we used a panel of genetically stable glycosylation mutated CHO cells that express defined glycans to evaluate the binding affinities of 10 different carbohydrate-recognition domains in galectins to N-glycans and mucin-type O-glycans. Using flow cytometry, we measured the cell-surface binding of the galectins. Moreover, we used fluorescence anisotropy to determine the galectin affinities to recombinant erythropoietin used as a reporter glycoprotein produced by the glycoengineered cells and to synthetic N-glycans with defined branch structures. We found that all galectins, apart from galectin-8N, require complex N-glycans for high-affinity binding. Galectin-8N targeted both N- and O-linked glycans with high affinity, preferring 2,3-sialylated N-acetyllactosamine (LacNAc) structures. Furthermore, we found that 2,3-sialylation suppresses high-affinity binding of select galectins, including galectin-2, -3, -4N, and -7. Structural modeling provided a basis for interpreting the observed binding preferences. These results underscore the power of a glycoengineered platform to dissect the glycan-binding specificities of carbohydrate-binding proteins.

Enhanced Human-Type Receptor Binding by Ferret-Transmissible H5N1 with a K193T Mutation.

All human influenza pandemics have originated from avian influenza viruses. Although multiple changes are needed for an avian virus to be able to transmit between humans, binding to human-type receptors is essential. Several research groups have reported mutations in H5N1 viruses that exhibit specificity for human-type receptors and promote respiratory droplet transmission between ferrets. Upon detailed analysis, we have found that these mutants exhibit significant differences in fine receptor specificity compared to human H1N1 and H3N2 and retain avian-type receptor binding. We have recently shown that human influenza viruses preferentially bind to α2-6-sialylated branched N-linked glycans, where the sialic acids on each branch can bind to receptor sites on two protomers of the same hemagglutinin (HA) trimer. In this binding mode, the glycan projects over the 190 helix at the top of the receptor-binding pocket, which in H5N1 would create a stearic clash with lysine at position 193. Thus, we hypothesized that a K193T mutation would improve binding to branched N-linked receptors. Indeed, the addition of the K193T mutation to the H5 HA of a respiratory-droplet-transmissible virus dramatically improves both binding to human trachea epithelial cells and specificity for extended α2-6-sialylated N-linked glycans recognized by human influenza viruses.IMPORTANCE Infections by avian H5N1 viruses are associated with a high mortality rate in several species, including humans. Fortunately, H5N1 viruses do not transmit between humans because they do not bind to human-type receptors. In 2012, three seminal papers have shown how these viruses can be engineered to transmit between ferrets, the human model for influenza virus infection. Receptor binding, among others, was changed, and the viruses now bind to human-type receptors. Receptor specificity was still markedly different compared to that of human influenza viruses. Here we report an additional mutation in ferret-transmissible H5N1 that increases human-type receptor binding. K193T seems to be a common receptor specificity determinant, as it increases human-type receptor binding in multiple subtypes. The K193T mutation can now be used as a marker during surveillance of emerging viruses to assess potential pandemic risk.

Three mutations switch H7N9 influenza to human-type receptor specificity.

The avian H7N9 influenza outbreak in 2013 resulted from an unprecedented incidence of influenza transmission to humans from infected poultry. The majority of human H7N9 isolates contained a hemagglutinin (HA) mutation (Q226L) that has previously been associated with a switch in receptor specificity from avian-type (NeuAcα2-3Gal) to human-type (NeuAcα2-6Gal), as documented for the avian progenitors of the 1957 (H2N2) and 1968 (H3N2) human influenza pandemic viruses. While this raised concern that the H7N9 virus was adapting to humans, the mutation was not sufficient to switch the receptor specificity of H7N9, and has not resulted in sustained transmission in humans. To determine if the H7 HA was capable of acquiring human-type receptor specificity, we conducted mutation analyses. Remarkably, three amino acid mutations conferred a switch in specificity for human-type receptors that resembled the specificity of the 2009 human H1 pandemic virus, and promoted binding to human trachea epithelial cells.

Integrated Omics and Computational Glycobiology Reveal Structural Basis for Influenza A Virus Glycan Microheterogeneity and Host Interactions.

Despite sustained biomedical research effort, influenza A virus remains an imminent threat to the world population and a major healthcare burden. The challenge in developing vaccines against influenza is the ability of the virus to mutate rapidly in response to selective immune pressure. Hemagglutinin is the predominant surface glycoprotein and the primary determinant of antigenicity, virulence and zoonotic potential. Mutations leading to changes in the number of HA glycosylation sites are often reported. Such genetic sequencing studies predict at best the disruption or creation of sequons for N-linked glycosylation; they do not reflect actual phenotypic changes in HA structure. Therefore, combined analysis of glycan micro and macro-heterogeneity and bioassays will better define the relationships among glycosylation, viral bioactivity and evolution. We present a study that integrates proteomics, glycomics and glycoproteomics of HA before and after adaptation to innate immune system pressure. We combined this information with glycan array and immune lectin binding data to correlate the phenotypic changes with biological activity. Underprocessed glycoforms predominated at the glycosylation sites found to be involved in viral evolution in response to selection pressures and interactions with innate immune-lectins. To understand the structural basis for site-specific glycan microheterogeneity at these sites, we performed structural modeling and molecular dynamics simulations. We observed that the presence of immature, high-mannose type glycans at a particular site correlated with reduced accessibility to glycan remodeling enzymes. Further, the high mannose glycans at sites implicated in immune lectin recognition were predicted to be capable of forming trimeric interactions with the immune-lectin surfactant protein-D.

Structural Analysis of the Glycosylated Intact HIV-1 gp120-b12 Antibody Complex Using Hydroxyl Radical Protein Footprinting.

Glycoprotein gp120 is a surface antigen and virulence factor of human immunodeficiency virus 1. Broadly neutralizing antibodies (bNAbs) that react to gp120 from a variety of HIV isolates offer hope for the development of broadly effective immunogens for vaccination purposes, if the interactions between gp120 and bNAbs can be understood. From a structural perspective, gp120 is a particularly difficult system because of its size, the presence of multiple flexible regions, and the large amount of glycosylation, all of which are important in gp120-bNAb interactions. Here, the interaction of full-length, glycosylated gp120 with bNAb b12 is probed using high-resolution hydroxyl radical protein footprinting (HR-HRPF) by fast photochemical oxidation of proteins. HR-HRPF allows for the measurement of changes in the average solvent accessible surface area of multiple amino acids without the need for measures that might alter the protein conformation, such as mutagenesis. HR-HRPF of the gp120-b12 complex coupled with computational modeling shows a novel extensive interaction of the V1/V2 domain, probably with the light chain of b12. Our data also reveal HR-HRPF protection in the C3 domain caused by interaction of the N330 glycan with the b12 light chain. In addition to providing information about the interactions of full-length, glycosylated gp120 with b12, this work serves as a template for the structural interrogation of full-length glycosylated gp120 with other bNAbs to better characterize the interactions that drive the broad specificity of the bNAb.

Asn347 Glycosylation of Corticosteroid-binding Globulin Fine-tunes the Host Immune Response by Modulating Proteolysis by Pseudomonas aeruginosa and Neutrophil Elastase.

Corticosteroid-binding globulin (CBG) delivers anti-inflammatory cortisol to inflamed tissues upon elastase-based proteolysis of the exposed reactive center loop (RCL). However, the molecular mechanisms that regulate the RCL proteolysis by co-existing host and bacterial elastases in inflamed/infected tissues remain unknown. We document that RCL-localized Asn(347) glycosylation fine-tunes the RCL cleavage rate by human neutrophil elastase (NE) and Pseudomonas aeruginosa elastase (PAE) by different mechanisms. NE- and PAE-generated fragments of native and exoglycosidase-treated blood-derived CBG of healthy individuals were monitored by gel electrophoresis and LC-MS/MS to determine the cleavage site(s) and Asn(347) glycosylation as a function of digestion time. The site-specific (Val(344)-Thr(345)) and rapid (seconds to minutes) NE-based RCL proteolysis was significantly antagonized by several volume-enhancing Asn(347) glycan features (i.e. occupancy, triantennary GlcNAc branching, and α1,6-fucosylation) and augmented by Asn(347) NeuAc-type sialylation (all p < 0.05). In contrast, the inefficient (minutes to hours) PAE-based RCL cleavage, which occurred equally well at Thr(345)-Leu(346) and Asn(347)-Leu(348), was abolished by the presence of Asn(347) glycosylation but was enhanced by sialoglycans on neighboring CBG N-sites. Molecular dynamics simulations of various Asn(347) glycoforms of uncleaved CBG indicated that multiple Asn(347) glycan features are modulating the RCL digestion efficiencies by NE/PAE. Finally, high concentrations of cortisol showed weak bacteriostatic effects toward virulent P. aeruginosa, which may explain the low RCL potency of the abundantly secreted PAE during host infection. In conclusion, site-specific CBG N-glycosylation regulates the bioavailability of cortisol in inflamed environments by fine-tuning the RCL proteolysis by endogenous and exogenous elastases. This study offers new molecular insight into host- and pathogen-based manipulation of the human immune system.

Induction of Antibodies Directed Against Branched Core O-Mannosyl Glycopeptides-Selectivity Complimentary to the ConA Lectin.

Mammalian protein O-mannosylation, initiated by attachment of α-mannopyranose to Ser or Thr residues, comprise a group of post-translational modifications (PTMs) involved in muscle and brain development. Recent advances in glycoproteomics methodology and the "SimpleCell" strategy have enabled rapid identification of glycoproteins and specific glycosylation sites. Despite the enormous progress made, the biological impact of the mammalian O-mannosyl glycoproteome remains largely unknown to date. Tools are still needed to investigate the structure, role, and abundance of O-mannosyl glycans. Although O-mannosyl branching has been shown to be of relevance in integrin-dependent cell migration, and also plays a role in demyelinating diseases, such as multiple sclerosis, a broader understanding of the biological roles of branched O-mannosyl glycans is lacking in part due to the paucity of detection tools. In this work, a glycopeptide vaccine construct was synthesized and used to generate antibodies against branched O-mannosyl glycans. Glycopeptide microarray screening revealed high selectivity of the induced antibodies for branched glycan core structures presented on different peptide backbones, with no cross-reactivity observed with related linear glycans. For comparison, microarray screening of the mannose-binding lectin concanavalin A (ConA), which is commonly used in glycoproteomics workflows to enrich tryptic O-mannosyl peptides, showed that the ConA lectin did not recognize branched O-mannosyl glycans. The binding preference of ConA for short linear O-mannosyl glycans was rationalized in terms of molecular structure using crystallographic data augmented by molecular modeling. The contrast between the ConA binding specificity and that of the new antibodies indicates a novel role for the antibodies in studies of protein O-mannosylation.

Combining 3D structure with glycan array data provides insight into the origin of glycan specificity.

Defining how a glycan-binding protein (GBP) specifically selects its cognate glycan from among the ensemble of glycans within the cellular glycome is an area of intense study. Powerful insight into recognition mechanisms can be gained from 3D structures of GBPs complexed to glycans; however, such structures remain difficult to obtain experimentally. Here an automated 3D structure generation technique, called computational carbohydrate grafting, is combined with the wealth of specificity information available from glycan array screening. Integration of the array data with modeling and crystallography allows generation of putative co-complex structures that can be objectively assessed and iteratively altered until a high level of agreement with experiment is achieved. Given an accurate model of the co-complexes, grafting is also able to discern which binding determinants are active when multiple potential determinants are present within a glycan. In some cases, induced fit in the protein or glycan was necessary to explain the observed specificity, while in other examples a revised definition of the minimal binding determinants was required. When applied to a collection of 10 GBP-glycan complexes, for which crystallographic and array data have been reported, grafting provided a structural rationalization for the binding specificity of >90% of 1223 arrayed glycans. A webtool that enables researchers to perform computational carbohydrate grafting is available at www.glycam.org/gr (accessed 03 March 2016).

Analysis of site-specific N-glycan remodeling in the endoplasmic reticulum and the Golgi.

The hallmark of N-linked protein glycosylation is the generation of diverse glycan structures in the secretory pathway. Dynamic, non-template-driven processes of N-glycan remodeling in the endoplasmic reticulum and the Golgi provide the cellular setting for structural diversity. We applied newly developed mass spectrometry-based analytics to quantify site-specific N-glycan remodeling of the model protein Pdi1p expressed in insect cells. Molecular dynamics simulation, mutational analysis, kinetic studies of in vitro processing events and glycan flux analysis supported the defining role of the protein in N-glycan processing.

Presentation, presentation, presentation! Molecular-level insight into linker effects on glycan array screening data.

Changes in cell-surface glycan patterns are markers of the presence of many different disease and cancer types, offering a relatively untapped niche for glycan-targeting reagents and therapeutics in diagnosis and treatment. Of paramount importance for the success of any glycan-targeting reagent is the ability to specifically recognize the target among the plethora of different glycans that exist in the human body. The preeminent technique for defining specificity is glycan array screening, in which a glycan-binding protein (GBP) can be simultaneously screened against multiple glycans. Glycan array screening has provided unparalleled insight into GBP specificity, but data interpretation suffers from difficulties in identifying false-negative binding arising from altered glycan presentation, associated with the linker used to conjugate the glycan to the surface. In this work, we model the structure and dynamics of the linkers employed in the glycan arrays developed by the Consortium for Functional Glycomics. The modeling takes into account the physical presence and surface polarity of the array, and provides a structure-based rationalization of false-negative results arising from the so-called "linker effect." The results also serve as a guide for interpreting glycan array screening data in a biological context; in particular, we show that attempts to employ natural amino acids as linkers may be prone to unexpected artifacts compromising glycan recognition.

Structure and binding analysis of Polyporus squamosus lectin in complex with the Neu5Ac{alpha}2-6Gal{beta}1-4GlcNAc human-type influenza receptor.

Glycan chains that terminate in sialic acid (Neu5Ac) are frequently the receptors targeted by pathogens for initial adhesion. Carbohydrate-binding proteins (lectins) with specificity for Neu5Ac are particularly useful in the detection and isolation of sialylated glycoconjugates, such as those associated with pathogen adhesion as well as those characteristic of several diseases including cancer. Structural studies of lectins are essential in order to understand the origin of their specificity, which is particularly important when employing such reagents as diagnostic tools. Here, we report a crystallographic and molecular dynamics (MD) analysis of a lectin from Polyporus squamosus (PSL) that is specific for glycans terminating with the sequence Neu5Acα2-6Galß. Because of its importance as a histological reagent, the PSL structure was solved (to 1.7 Å) in complex with a trisaccharide, whose sequence (Neu5Acα2-6Galß1-4GlcNAc) is exploited by influenza A hemagglutinin for viral adhesion to human tissue. The structural data illuminate the origin of the high specificity of PSL for the Neu5Acα2-6Gal sequence. Theoretical binding free energies derived from the MD data confirm the key interactions identified crystallographically and provide additional insight into the relative contributions from each amino acid, as well as estimates of the importance of entropic and enthalpic contributions to binding.

Landscape and selection of vaccine epitopes in SARS-CoV-2.

BACKGROUND: Early in the pandemic, we designed a SARS-CoV-2 peptide vaccine containing epitope regions optimized for concurrent B cell, CD4+ T cell, and CD8+ T cell stimulation. The rationale for this design was to drive both humoral and cellular immunity with high specificity while avoiding undesired effects such as antibody-dependent enhancement (ADE). METHODS: We explored the set of computationally predicted SARS-CoV-2 HLA-I and HLA-II ligands, examining protein source, concurrent human/murine coverage, and population coverage. Beyond MHC affinity, T cell vaccine candidates were further refined by predicted immunogenicity, sequence conservation, source protein abundance, and coverage of high frequency HLA alleles. B cell epitope regions were chosen from linear epitope mapping studies of convalescent patient serum, followed by filtering for surface accessibility, sequence conservation, spatial localization near functional domains of the spike glycoprotein, and avoidance of glycosylation sites. RESULTS: From 58 initial candidates, three B cell epitope regions were identified. From 3730 (MHC-I) and 5045 (MHC-II) candidate ligands, 292 CD8+ and 284 CD4+ T cell epitopes were identified. By combining these B cell and T cell analyses, as well as a manufacturability heuristic, we proposed a set of 22 SARS-CoV-2 vaccine peptides for use in subsequent murine studies. We curated a dataset of ~ 1000 observed T cell epitopes from convalescent COVID-19 patients across eight studies, showing 8/15 recurrent epitope regions to overlap with at least one of our candidate peptides. Of the 22 candidate vaccine peptides, 16 (n = 10 T cell epitope optimized; n = 6 B cell epitope optimized) were manually selected to decrease their degree of sequence overlap and then synthesized. The immunogenicity of the synthesized vaccine peptides was validated using ELISpot and ELISA following murine vaccination. Strong T cell responses were observed in 7/10 T cell epitope optimized peptides following vaccination. Humoral responses were deficient, likely due to the unrestricted conformational space inhabited by linear vaccine peptides. CONCLUSIONS: Overall, we find our selection process and vaccine formulation to be appropriate for identifying T cell epitopes and eliciting T cell responses against those epitopes. Further studies are needed to optimize prediction and induction of B cell responses, as well as study the protective capacity of predicted T and B cell epitopes.