Analysis of MHC Class II Antigen Processing by Quantitation of Peptides that Constitute Nested Sets

Transcription

1 This information is current as of August 20, References Subscription Permissions Alerts Analysis of MHC Class II Antigen Processing by Quantitation of Peptides that Constitute Nested Sets John D. Lippolis, Forest M. White, Jarrod A. Marto, Chance J. Luckey, Timothy N. J. Bullock, Jeffrey Shabanowitz, Donald F. Hunt and Victor H. Engelhard J Immunol 2002; 169: ; ; doi: /jimmunol This article cites 32 articles, 12 of which you can access for free at: Why The JI? Submit online. Rapid Reviews! 30 days* from submission to initial decision No Triage! Every submission reviewed by practicing scientists Fast Publication! 4 weeks from acceptance to publication *average Information about subscribing to The Journal of Immunology is online at: Submit copyright permission requests at: Receive free -alerts when new articles cite this article. Sign up at: Downloaded from by guest on August 20, 2018 The Journal of Immunology is published twice each month by The American Association of Immunologists, Inc., 1451 Rockville Pike, Suite 650, Rockville, MD Copyright 2002 by The American Association of Immunologists All rights reserved. Print ISSN: Online ISSN:

2 The Journal of Immunology Analysis of MHC Class II Antigen Processing by Quantitation of Peptides that Constitute Nested Sets John D. Lippolis, 2 * Forest M. White, 2 Jarrod A. Marto, Chance J. Luckey,* Timothy N. J. Bullock,* Jeffrey Shabanowitz, Donald F. Hunt, and Victor H. Engelhard 3 * Peptides associated with class II MHC molecules are of variable length because in contrast to peptides associated with class I MHC molecules, their amino and C termini are not constrained by the structure of the peptide interaction with the binding site. The proteolytic processing events that generate these peptides are still not well understood. To address this question, peptides extracted from HLA-DR*0401 were analyzed using two types of mass spectrometry instrumentation. This enabled identification of >700 candidate peptides in a single analysis and provided relative abundance information on 142 peptides contained in 11 nested sets of 3 36 members each. Peptides of 12 residues or less occurred only at low abundance, despite the fact that they were predicted to fully occupy the HLA-DR*0401 molecule in a single register. Conversely, the relative abundance of longer species suggested that proteolytic events occurring after MHC binding determine the final structure of most class II-associated peptides. Our data suggest that C-terminal residues of these peptides reflect the action of peptidases that cleave at preferred amino acids, while amino termini appear to be determined more by proximity to the class II MHC binding site. Thus, the analysis of abundance information for class II-associated peptides comprising nested sets has offered new insights into proteolytic processing of MHC class IIassociated peptides. The Journal of Immunology, 2002, 169: Class II MHC molecules bind peptides and present them to CD4 T cells. The major sources of such peptides are either external or membrane-associated proteins (1) that are internalized and proteolyzed in endosomes or lysosomes (2 4). The selective presentation by class II MHC molecules of peptides produced in endosomal compartments is due to the involvement of the invariant chain, which associates with the class II MHC - and -chains immediately after synthesis in the endoplasmic reticulum, blocks the peptide binding site, and retargets the complex to the endosomal/lysosomal compartment (5). Endosomal proteases digest most of the invariant chain, leaving a set of peptides called class II-associated invariant chain peptide associated with the class II peptide binding site (2, 4). HLA-DM catalyzes the dissociation of the class II-associated invariant chain peptide peptides, allowing for the subsequent binding of other peptides present in the endosomal compartment (6, 7). Although the assembly, transport, and peptide binding of class II MHC molecules is relatively well understood, little is known about the processes and proteases that form these peptides. Although endosomes contain a variety of proteases (reviewed in Ref. 4), few have been directly implicated in the generation of class II-associated peptides. Cathepsin D has been implicated in the generation of *Department of Microbiology and Carter Immunology Center, University of Virginia, Charlottesville, VA 22908, and Departments of Chemistry and Pathology, University of Virginia, Charlottesville, VA Received for publication June 19, Accepted for publication August 29, The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. 1 This work was supported by U.S. Public Health Service Grants AI20963 (to V.H.E.), AI33993 (to D.F.H.), and AI45199 (to V.H.E. and D.F.H.). C.J.L. was supported by Public Health Service Medical Scientist Training Program Grant (GM07267). 2 J.D.L. and F.M.W. have contributed equally to this work. 3 Address correspondence and reprint requests to Dr. Victor H. Engelhard, Carter Immunology Center, University of Virginia, P.O. Box , Charlottesville, VA address: vhe@virginia.edu T cell epitopes from OVA, hen egg lysozyme, sperm whale myoglobin, and insulin based on in vitro digestion to produce antigenic fragments (8 12). However, splenocytes from cathepsin B or cathepsin D-deficient mice process and present multiple antigenic determinants normally, including those from OVA and hen egg lysozyme (13). Recently, a newly identified asparaginyl endopeptidase was shown to generate a proteolytic intermediate in the processing of a tetanus toxin Ag (14). This study also suggested the existence of at least two processing events: an initial generation of precursor fragments from the intact protein and secondary processing to an optimal size for MHC binding and/or T cell recognition. However, it remains unclear whether processing is completed before binding to class II MHC molecules in the endosome, or whether some processing of MHC-associated peptides occurs subsequent to binding. Binding of peptides to class II MHC molecules involves several hydrogen bonds between residues in the binding site and the peptide backbone, as well as contacts between specific peptide side chains and pockets in the binding site (15 17). In contrast to class I MHC molecules, the ends of the class II molecule binding site are open and do not constrain the length of class II-associated peptides. Indeed, sequence analysis of peptides associated with class II MHC molecules has identified peptides of lengths between 10 and 34 residues (1, 18 20). Furthermore, many of these class II-associated peptides comprise nested sets that contain the same core sequence but vary in length at the amino or C-terminal ends. However, little information has been gathered concerning the relative abundance of related peptides or the factors that determine the preferred length of class II-associated peptides. Previous sequence analyses of class II MHC-associated peptides have relied on either Edman degradation or mass spectrometry (1, 18 22). The utility of Edman degradation is limited by the complexity of the peptide mixture even after HPLC fractionation. Chicz and colleagues (1, 20) used mass spectrometry to correlate masses in individual HPLC fractions with Edman data, and thus inferred the identity of several peptides from different human class Copyright 2002 by The American Association of Immunologists, Inc /02/$02.00

3 5090 ANALYSIS OF MHC CLASS II-ASSOCIATED PEPTIDE NESTED SETS II alleles. Although this improves on the use of Edman degradation alone, the mass accuracy of the instrument used does not allow for truly unambiguous identification of such peptides. Only by using tandem mass spectrometry (MS/MS) 4 is it possible to select an individual peptide ion and subject it to collision-activated dissociation to generate an unambiguous sequence (19, 22). In the present study, we used a high throughput mass spectrometric technique together with automated database searching of collision-activated dissociation spectra to determine candidate sequences for 700 different peptides bound to the human class II MHC molecule HLA-DR*0401. Next, we used Fourier Transform mass spectrometry to analyze 142 peptides that were members of nested sets and determine their relative abundances. Using this information, we provide a comprehensive view of the distribution of peptides within a nested set that are associated with this class II MHC molecule. Our analysis has also allowed us to draw inferences regarding the proteolytic processes that lead to final peptide products displayed by HLA-DR4. Materials and Methods Cell lines The human B lymphoblastoid cell line PRIESS (kindly provided by B. Nepom, Virginia Mason Research Center, Seattle, WA) is homozygous for HLA-DR*0401 (DRA*0101, DRB1*0401). The hybridoma LB3.1 secretes an mab that specifically recognizes human HLA-DR molecules (23). Both cell lines were maintained in RPMI 1640 supplemented with 10% FBS and 2 mm glutamine. Isolation of HLA-DR*0401-associated peptides PREISS cells ( ) were lysed in 5 ml of lysis buffer (20 mm Tris-HCl (ph 8.0), 150 mm NaCl, 1% 3-[(3-cholamidopropyl)dimethylammonio]- 1-propanesulfonate, 5 g/ml aprotinin, 10 g/ml leupeptin, 10 g/ml pepstatin A, 5 mm EDTA, 0.04% sodium azide, and 1 mm PMSF) for 1hat 4 C. The lysate was clarified by centrifugation at 16,000 g for 30 min at 4 C, and the supernatant was precleared by incubation with 100 l recombinant protein A-Sepharose beads (Amersham Pharmacia Biotech, Piscataway, NJ) for 4 h at 4 C. After removal of the beads, the supernatant was incubated overnight at 4 C with 100 l recombinant protein A-Sepharose beads to which 2 mg of LB3.1 Ab had been bound. The beads were subsequently washed twice in lysis buffer, four times in 20 mm Tris-HCl (ph 8.0), 150 mm NaCl, twice in 20 mm Tris (ph 8.0), 1 M NaCl, and three times in 20 mm Tris-HCl (ph 8.0). Peptides were eluted in acid and separated from beads, Ab, and class II MHC molecules by passage through a 10,000-day cutoff ultrafiltration unit (Millipore, Bedford, MA). Mass spectrometry An aliquot corresponding to cell equivalents of the peptide mixture described above was analyzed by nanoflow HPLC-electrospray ionization coupled directly to an ion trap mass spectrometer (24, 25). Following datadependent MS/MS analysis (24, 25), spectra with characteristic features of peptide fragmentation were searched against the nonredundant database maintained at the National Center for Biotechnology Information using the SEQUEST algorithm (26) to identify candidate sequences. Selected sequences were confirmed by manual interpretation of MS/MS spectra and their identification as members of nested sets. A second aliquot corresponding to cell equivalents of the peptide mixture was subsequently analyzed by nanoflow HPLC coupled to electrospray ionization on a homebuilt Fourier transform mass spectrometer (27). Full scan mass spectra (300 m/z 5000) were acquired at approximately one scan per second. Peptide ions were correlated with those identified in ion trap experiments by relative HPLC retention time and accurate mass. Peptide abundance was estimated based on the observed ion current. Because peptides in nested sets are chemically similar to one another, it is reasonable to expect that they will exhibit similar ionization efficiencies. 4 Abbreviations used in this paper: MS/MS, tandem mass spectrometry; Ig, human Ig -chain; BSP, bovine serotransferrin precursor; NSFA, human -soluble NSF attachment protein; LAMP, human lysosomal-associated multitransmembrane protein; HSC71, human heat shock cognate 71-kDa protein. Results Rapid identification and characterization of HLA-DR*0401- associated peptides belonging to nested sets The peptides associated with immunoaffinity-purified HLA- DR*0401 molecules were extracted and analyzed by nanoflow HPLC interfaced to electrospray ionization on an ion-trap mass spectrometer. By using a long, shallow HPLC gradient and by acquiring mass spectra in the data-dependent mode (24), 2000 nonredundant peptide MS/MS spectra were acquired in a single nanoflow HPLC interfaced to electrospray ionization on an iontrap mass spectrometer run. A search of these MS/MS spectra against the NCBI nonredundant database using SEQUEST yielded 701 candidate peptide sequences. Based on estimates of 10,000 peptides presented by class I and II MHC molecules on the surface of a typical lymphoid cell (28), this represents 7% of the total peptide repertoire presented by HLA-DR*0401 on these cells. Among these 701 candidate peptide sequences, 273 were identifiable as members of nested sets (defined as groups of peptides that share core sequences but have distinct amino and C-terminal ends). By manual interpretation of the MS/MS spectra, we confirmed the sequences of 142 peptides that comprised 11 nested sets (Table I). The sample of HLA-DR*0401-associated peptides was also analyzed using a Fourier Transform mass spectrometer, which has enhanced sensitivity and mass accuracy compared with other instruments and also provided quantitative information on the 142 nested set peptides. We used this information to examine the display of related peptides derived from the same protein sequence. We found that the most abundant individual peptide species in the different nested sets varied between 14 and 21 residues long, while the mean length of all peptides in a nested set, after normalizing for the abundance of individual species, varied between 14.5 and 17.8 residues (Table II). Based on visibility in the x-ray crystallographic structure of HLA-DR*0401 associated with the collagen II peptide and involvement in hydrogen bonds, only 12 peptide residues interact with the binding site of this class II MHC molecule (17). However, only 8 of the 142 peptides were 12 residues or smaller, and these were invariably of very low relative abundance (Table I). These data suggested that features of each composite nested set sequence apart from the minimal sequences required for binding to HLA-DR*0401, determined both the optimal length, and distribution of displayed peptides. Alignment of nested set peptides to HLA-DR4 binding motif To facilitate understanding the factors contributing to the presentation of nested set peptides, we used three predictive algorithms to determine the likely binding frames of the peptides in each nested set. The three binding algorithms used different methodology, but predicted the same preferred frame in each case. For human Ig -chain (Ig ) , human -soluble N-ethylaleimidesensitive factor attachment protein , human lysosomalassociated multitransmembrane protein (LAMP) , and HLA- C , these algorithms identified single binding frames that were contained within all peptides comprising these sets and had relative affinities significantly higher than any other predicted frame (Table III). The nested sets from HLA-C , HLA-B , and HLA-E are derived from homologous regions of these three class I MHC molecules, and the sequences contained in any two shares of 19 possible residues. A single homologous high-affinity binding frame accommodated 5/5, 14/16, and 10/11 peptides from these sets, respectively. Similarly, the nested sets HLA-C 52 74, HLA-B 51 72, and HLA-A were also derived from homologous regions of these proteins, and share 20/22 residues. In fact, 12 peptides were derived from a region that is

4 The Journal of Immunology 5091 Table I. Major nested sets extracted from HLA-DR4 a Nested Set Residues Peptide Sequence Abundance ( 10 7 ) b HLA-A DTQFVRFDSDAASQRMEPR VDDTQFVRFDSDAASQR DTQFVRFDSDAASQR VDDTQFVRFDSDAASQRMEP DTQFVRFDSDAASQRMEP DDTQFVRFDSDAASQR DTQFVRFDSDAASQRME TQFVRFDSDAASQR DDTQFVRFDSDAASQRMEPR VDDTQFVRFDSDAASQRME DTQFVRFDSDAASQ VDDTQFVRFDSDAASQ DDTQFVRFDSDAASQRMEP DDTQFVRFDSDAASQRME DTQFVRFDSDAASQRM TQFVRFDSDAASQRME FVRFDSDAASQRM c FVRFDSDAASQR VDDTQFVRFDSDAASQRMEPRAP VDDTQFVRFDSDAASQRM TQFVRFDSDAASQ FVRFDSDAASQRMEP YVDDTQFVRFDSDAASQRMEPRAP DDTQFVRFDSDAASQ VDDTQFVRFDSDAAS DDTQFVRFDSDAASQRM FVRFDSDAASQRME VRFDSDAASQR d TQFVRFDSDAASQRM FVRFDSDAASQ FVRFDSDAASQRMEPR VDDTQFVRFDSDAASQRMEPRA VGYVDDTQFVRFDSDAASQR QFVRFDSDAASQR RFDSDAASQR TQFVRFDSDAASQRMEP 0.1 HLA-B VDDTQFVRFDSDAASPRMAPR VDDTQFVRFDSDAASPRMA VDDTQFVRFDSDAASPRMAP VDDTQFVRFDSDAASPR 14 e DTQFVRFDSDAASPR 12 e VDDTQFVRFDSDAASPRM DTQFVRFDSDAASPRMAPR DTQFVRFDSDAASPRM DTQFVRFDSDAASPRMAP DDTQFVRFDSDAASPR 4.9 e VDDTQFVRFDSDAASP 3.5 e DTQFVRFDSDAASP 2.3 e DDTQFVRFDSDAASPRMAPR DTQFVRFDSDAASPRMA DDTQFVRFDSDAASPRM DDTQFVRFDSDAASPRMAP DDTQFVRFDSDAASPRMA DDTQFVRFDSDAASP 1.3 e TQFVRFDSDAASPR 1.1 e FVRFDSDAASPR 0.4 e YVDDTQFVRFDSDAASPR 0.3 e FVRFDSDAASP 0.2 e VRFDSDAASPR 0.2 e QFVRFDSDAASPR 0.1 e HLA-C VDDTQFVRFDSDAASPRGEPR VDDTQFVRFDSDAASPR 116 e DTQFVRFDSDAASPR 98 e DTQFVRFDSDAASPRGEPR DDTQFVRFDSDAASPRGEPR DDTQFVRFDSDAASPR 41 e VDDTQFVRFDSDAASP 29 e VDDTQFVRFDSDAASPRGEPRAP VDDTQFVRFDSDAASPRGEP DTQFVRFDSDAASP 20 e DTQFVRFDSDAASPRGEP DDTQFVRFDSDAASP 11 e TQFVRFDSDAASPR 8.9 e DDTQFVRFDSDAASPRGE FVRFDSDAASPR 3.4 e DTQFVRFDSDAASPRGE 2.9 (Table continues)

5 5092 ANALYSIS OF MHC CLASS II-ASSOCIATED PEPTIDE NESTED SETS Table I. Continued Nested Set Residues Peptide Sequence Abundance ( 10 7 ) b YVDDTQFVRFDSDAASPR 2.9 e VRFDSDAASPR 2.1 e FVRFDSDAASP 1.4 e QFVRFDSDAASPR 0.6 e HLA-B LSSWTAADTAAQITQ DLSSWTAADTAAQITQ DLSSWTAADTAAQITQR EDLSSWTAADTAAQITQR LSSWTAADTAAQIT LSSWTAADTAAQITQR EDLSSWTAADTAAQITQ DLSSWTAADTAAQITQRK DLSSWTAADTAAQIT EDLSSWTAADTAAQITQRK NEDLSSWTAADTAAQITQ EDLSSWTAADTAAQIT SSWTAADTAAQITQ LNEDLSSWTAADTAAQITQ EDLSSWTAADTAAQ DLSSWTAADTAAQ 0.2 HLA-C DLRSWTAADTAAQITQ DLRSWTAADTAAQITQR NEDLRSWTAADTAAQITQ EDLRSWTAADTAAQITQR LNEDLRSWTAADTAAQITQR 0.6 HLA-E NEDLRSWTAVDTAAQISEQ DLRSWTAVDTAAQISEQ DLRSWTAVDTAAQISEQK DLRSWTAVDTAAQ LRSWTAVDTAAQISE EDLRSWTAVDTAAQISEQK DLRSWTAVDTAAQISE EDLRSWTAVDTAAQISEQ DLRSWTAVDTAAQIS NEDLRSWTAVDTAAQISEQK EDLRSWTAVDTAAQISE 0.1 Ig KVQWKVDNALQSGNS KVQWKVDNALQSGN YPREAKVQWKVDNALQSGN AKVQWKVDNALQSGN YPREAKVQWKVDNALQSGNS KVQWKVDNALQSGNSQ PREAKVQWKVDNALQSGN AKVQWKVDNALQSGNS KVQWKVDNALQSGNSQE EAKVQWKVDNALQSGN YPREAKVQWKVDNALQSGNSQ EAKVQWKVDNALQSGNS EAKVQWKVDNALQSGNSQ AKVQWKVDNALQSGNSQE KVQWKVDNALQSGNSQES AKVQWKVDNALQSGNSQES AKVQWKVDNALQSGNSQ KVQWKVDNALQSG 0.4 BSP DVAFVKDQTVIQNTD GDVAFVKDQTVIQNTD DVAFVKDQTVIQNTDG DVAFVKDQTVIQN GDVAFVKDQTVIQ GDVAFVKDQTVIQN DVAFVKDQTVIQNTDGN DVAFVKDQTVIQ GDVAFVKDQTVIQNTDG GDVAFVKDQTVIQNTDGN GDVAFVKDQTVIQNTDGNN 0.4 NSFA SADYYKGEESNSSANK DYYKGEESNSSANK ADYYKGEESNSSANK 0.6 (Table continues)

6 The Journal of Immunology 5093 Table I. Continued Nested Set Residues Peptide Sequence Abundance ( 10 7 ) b DYYKGEESNSSAN 0.2 LAMP LPSYEEALSLPSKTPE LPSYEEALSLPSKTPEG LPSYEEALSLPSKTP VVLPSYEEALSLPSKTPE VVLPSYEEALSLPSKTPEG VLPSYEEALSLPSKTPEG 0.3 HSC GILNVSAVDKSTGK GILNVSAVDKSTGKE ANGILNVSAVDKSTGKE 0.3 a Peptides were identified in HLA-DRB*0401 extracts from MS/MS spectra that were recorded using an ion trap mass spectrometer and searched against a sequence database using the SEQUEST algorithm. b Peptide abundance is estimated based on the observed ion current. c Sequences in bold do not contain the primary predicted binding frame identified in Table II. d Sequences in italics are shorter than 12 residues. e Because residues of HLA-B and HLA-C are identical, the exact contribution of these two sources to peptides contained within this sequence could not be directly determined. Therefore, the total ion abundances for all HLA-B and all HLA-C peptides that terminated at residue 72 were calculated, and their ratio was then used to calculate the amounts of each common peptide that would arise from each molecule. Peptides terminating at residue R72 were chosen because no peptides for HLA-C were detected with residue G69 as a terminus. Therefore, we could not sum the abundances for residues because the abundance of peptides terminating at residue 69 is not 0 but below detection threshold. In addition, residues 69 and 70 are unique to each molecule. If site-specific endoproteases do generate these peptides, then any difference in amino acid could account for changes in abundance. identical in HLA-B and -C molecules; and therefore, could come from either or both of these molecules. A common homologous high-affinity frame accounted for the binding of 28/32 peptides from the combined HLA-B and -C sets and 27/36 of the peptides from the HLA-A It should be noted that the 16 of 128 peptides from these 9 nested sets that did not contain the complete predicted high affinity binding frame generally fell into two categories (Table I). Five were shorter than 12 residues; and therefore, could not fully contain any frame based on this number of residues. Of the remaining 11, only the 3 peptides of the HLA-B and HLA-E sets that terminated in Q could fully occupy any discernable alternate binding frame, albeit of low affinity (Table III). Regardless, these nonconforming peptides were generally among the lowest abundance species of a particular set. Thus, the predictive algorithm, together with the abundance data provided from the Fourier Transform mass spectrometry experiments, suggests that the single high-affinity binding frames in Table III represent the dominant mode of binding for the peptides in these nested sets. The two nested sets derived from bovine serotransferrin precursor (BSP) and human heat shock cognate 71-kDa protein (HSC71) had two frames of similar predicted affinity that were shifted by one residue with respect to one another. These two frames also accommodated the same sets of peptides. In addition, no predicted alternate binding frame can account for the two peptides of the BSP set that terminate at Q, and thus do not contain the complete primary predicted frame. In these cases, it cannot be determined what proportions of these peptides will be displayed by HLA-DR*0401 in one or both frames. Distribution of peptide length and abundance within a nested set reveals holes in the presented peptide repertoire With the information about likely binding frames in hand, we reexamined in more detail the abundance of individual peptides within two representative nested sets containing a minimum of 11 unique members. The 16-member HLA-B set was dominated by four peptides ( (15 mer), (16 mer), Table II. Length distribution of peptides in nested sets a Set Mean Length Mean Length Normalized for Abundance Most Abundant Species Sequence Length HLA-C ND b ND b VDDTQFVRFDSDAASPRGEPR 21 HLA-B ND b ND b VDDTQFVRFDRDAASPRMAPR 21 HLA-A DTQFVRFDSDAASQRMEPR 19 HLA-E NEDLRSWTAVDTAAQISEQ 19 LAMP LPSYEEALSLPSKTPE 16 HLA-C DLRSWTAADTAAQITQ 16 HLA-B LSSWTAADTAAQITQ 15 Ig KVQWKVDNALQSGNS 15 BSP DVAFVKDQTVIQNTD 15 NSFA DYYKGEESNSSANK 14 HSC GILNVSAVDKSTGK 14 a The mean length was calculated simply by summing the total number of residues in each sequence of a nested set and dividing by the number of sequences. The mean length normalized for abundance was calculated by multiplying the length of each sequence in a nested set by its abundance in Table I, summing these values, and dividing by the total abundance of all members of the set. Determination of most abundant species is based on the data in Table I. b Values for these nested sets were ND, because several peptides could not be unambiguously assigned to one set or another.

7 5094 ANALYSIS OF MHC CLASS II-ASSOCIATED PEPTIDE NESTED SETS Table III. Alignment of nested set peptides to predicted HLA-DR4-binding motifs Set Alignment and Predicted Binding Frame Algorithm P( 2) P( 1) P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 1 a (IC 50 ) 2 b 3 c Peptides Containing Complete Frame d Ig V Q W K V D N A L Q S G /18 NSFA D Y Y K G E E S N S S A /4 A D Y Y K G E E S N S S /4 LAMP P S Y E E A L S L P S K /6 HLA-C R S W T A A D T A A Q I /5 E D L R S W T A A D T A /5 HLA-B S S W T A A D T A A Q I /16 E D L S S W T A A D T A /16 (1) HLA-E R S W T A V D T A A Q I /11 E D L R S W T A V D T A /11 (0) HLA-C T Q F V R F D S D A A S /20 Q F V R F D S D A A S P /20 (1) HLA-B T Q F V R F D S D A A S /24 Q F V R F D S D A A S P /24 (1) HLA-A T Q F V R F D S D A A S /36 Q F V R F D S D A A S Q /36 (1) V G Y V D D T Q F V R F /36 (0) BSP V A F V K D Q T V I Q N /11 A F V K D Q T V I Q N T /11 (0) HSC G I L N V S A V D K S T /3 N G I L N V S A V D K S /3 a IC 50 values were calculated according to Ref. 32. This algorithm does not take into account the effect of an amino terminal P( 2) residue, which has been shown by crystallographic analysis to interact with the class II-binding site by hydrogen bonding interactions with the peptide backbone (17). All other potential binding frames had calculated IC 50 values 40,000 and have been omitted for clarity. b Relative binding affinity values of 9-mer peptides initiating at P1 were calculated according to the algorithm of Southwood et al. (33). Higher values predict higher binding affinity, and values of 0 predict no binding. The average relative binding affinity for all class II-associated peptides is No other binding frames containing the peptides in Table I scored 0. c Binding scores for 15-mer peptides with position 4 corresponding to P1 in the table were obtained using the SYFPEITHI program at Higher values predict higher binding affinity. The maximum possible score is 36 and the minimum is 0. No other binding frames in Table I scored 0. d In cases where the primary binding frame was not contained by all peptides in a nested set, the numbers in parenthesis indicates the additional peptides of that set that contain the complete secondary binding frame (17 mer), and (18 mer)) that collectively accounted for almost 80% of the total peptide abundance (Fig. 1, upper panel). In contrast, eight peptides in this set were each present at 1.5%. The distribution of peptide lengths within this set is relatively regular and continuous, and the dominant peptides were related to one another by the gain or loss of single residues at either the amino or the C termini. A small number of dominant species was also evident in the Ig set, with the most abundant three peptides ( , , and ) accounting for 64% of the total peptide abundance (Fig. 1, lower panel). However, the species was not related to the other dominant peptides by the gain or loss of a single residue at the amino terminal end, but instead involved a difference of five residues. More generally, Ig -derived peptides beginning at residues 140 or 145 were relatively frequent, while those beginning at were much less abundant. At the extreme, no Ig peptides beginning at R 142 were found, leading to a distribution gap and two clusters of peptides based on their amino termini. It is unlikely that an alternate binding frame could account for this discontinuous distribution, since only a single binding frame of any significance was identified (Table III). Thus, these features are most likely due to the specificity of the proteolytic enzymes that produce these peptides during Ag processing. Summed abundance of peptides with common amino and C termini provides evidence for the action of specific proteases The data of Fig. 1, taken together with the differences in dominant and mean peptide lengths among the nested sets and the predicted binding motif information, suggested the possibility that the proteases responsible for the generation of the peptides associated with HLA-DR4 preferentially cleave at certain residues. Alternatively, the HLA-DR4 binding site might sterically inhibit access of a nonspecific exoprotease or endoprotease to the ends of an already bound peptide, even though those ends might not be immobilized in the site. To gain more information about these possibilities, we summed the abundances of all peptides in each of the nested sets that either began or ended at a particular residue. For example, 5 of the 18 peptides in the Ig nested set ended with N 158, and their summed abundance is (Table I). Similarly, six peptides from this set begin with K 145, and their summed relative abundance is By combining this data with the predicted binding frames in Table II, we examined whether the distribution of peptide ends was related simply to binding site proximity, or suggested specificity in proteolysis (Fig. 2). We first examined the relationship of dominant C termini to one another and to the P(-2) through P10 residues predicted to interact with the HLA-DR4 binding site. Only in the case of the Ig set was there significant representation of peptides that terminated at P11: the most abundant C-terminal residues for the remaining 10 sets were usually at the P12 or P13 positions. While the paucity of peptides terminating at P11 could be due to steric inhibition of protease access, we also observed that certain residues were preferred cleavage sites although they were more distant from the apparent end of the binding site. Thus, peptides of the HLA-C 52 74, HLA-B 51 72, and HLA-A sets terminated at R at the P12 or P16 positions with a similar frequency, and peptides terminating at one of these two residues were more abundant than those terminating at intervening residues. R was also found as a strong secondary C terminus at the P13 position for peptides of the HLA- C and HLA-B sets. Similarly, Q was the preferred C-terminal residue for the HLA-C , HLA-B , and HLA-E sets. While Q occupies the P12 position in the first two sets, it occupies the P13 position in the latter, although the sequences of the binding frames for these nested sets are otherwise

8 The Journal of Immunology 5095 FIGURE 1. Abundance of individual peptides in two nested sets associated with HLA-DR*0401. Peptides are identified based on their amino and C-terminal residues as indicated on the X and Y axes, respectively. Abundance data is taken from Table I. Abundance normalized length was calculated by multiplying the length of each sequence in a nested set by its abundance in Table I, summing these values, and dividing by the total abundance of all members of the set. highly homologous. These preferred cleavage sites were not explainable as steric hindrance due to the use of a secondary binding frame, since these frames all would have protected amino terminal, rather than C-terminal residues. The enhanced representation of certain C-terminal residues in multiple sequences at different positions is consistent with the hypothesis that they are due to selective protease activity rather than simply by steric constraint due to proximity to the peptide binding site. In contrast to the C termini, dominant amino termini were often found immediately adjacent to the peptide binding site at the P(-3) position. In two cases (human -soluble NSF attachment protein (NSFA) and HSC ), the dominant terminal residues were at the P(-2) position. In neither case could this observation be accounted for by a secondary binding frame. It is also more difficult to discern evidence for the specificity of proteases that produce the amino termini of the remaining nested sets. These observations are consistent with the idea that steric constraints on protease activity are less severe at this end of the peptide binding site, or that substantially complete processing of peptide amino termini occurs before stable HLA-DR4 binding. Discussion In this study, we used quadrupole ion trap mass spectrometry technology to identify 700 candidate peptides associated with HLA- DR*0401. Of these sequences, 142 were contained within 11 nested sets. The sequences of these peptides were manually confirmed and then quantitated using Fourier transform mass spectrometry. This quantitative information concerning constituent peptide species of nested sets provides several advantages: 1) nested sets derive from a common protein source, thereby eliminating any skew introduced by differences in protein expression; 2) any observed differences in peptide distribution must result from proteolytic processing or binding to the class II MHC molecule; and 3) comparison of multiple nested sets allows us to define general trends in peptide processing. Based on hydrogen bond interactions and visibility in the x-ray crystallographic structure of HLA-DR*0401 (DRA*0101, DRB1*0401) associated with the collagen II peptide (17), as well as other crystal structures of HLA-DR molecules (15), only 12 peptide residues interact with the peptide binding site of this molecule. Furthermore, studies using sequentially smaller peptides have shown that 11 mers bind to HLA-DR4 with affinities comparable to those 13 mers (29). Our work confirms on a large scale that the average naturally occurring class II-associated peptides is four to seven residues longer than necessary for maximal hydrogen bonding between the peptide and the class II molecule. However, perhaps most surprising was that 12 mers corresponding to this ideal binding site were not observed among the peptides of most nested sets, or if observed, were of low relative abundance. Although it is formally possible that these additional residues contribute to the affinity of these particular peptides, these observations suggest that most peptides interact initially with the class II binding site in the context of larger precursors, which are then proteolyzed further at both amino and C-terminal ends. The additional residues beyond the minimum necessary for binding site occupancy are either protected from proteolytic removal after the peptide is bound to the class II molecule or are a reflection of the nonrandom cleavage of peptides by endosomal proteases. However, in a small number of cases, a fraction of the peptides from individual nested sets terminate at residues that lie within the peptide binding site. Such peptide termini likely represent dominant protease cleavage sites that are used early in the degradation process. Several observations are consistent with the idea that the final form of peptides displayed by HLA-DR*0401 reflects the action of proteases with preferences for certain amino acid side chains, rather than simply amino and carboxyl peptidase digestion of exposed ends until steric hindrance prevents further action. First, the average peptide lengths and length of the most dominant species differ among nested sets. Second, when individual species were aligned, some sets showed evidence of gaps or low frequency of occurrence of particular amino or C termini. Third, preferential occurrence of R and Q at the C terminus was observed when different nested sets were compared. This preference was only partly dependent on the absolute position relative to the peptide binding site. The specificities of proteases that have been localized to endosomes are rather broad (4), and so it may be somewhat surprising that any evidence of selectivity in proteolysis was observed in this work. However, with the exception of a recently described asparaginyl endopeptidase (14), none of these proteolytic activities have been directly implicated in the generation of class II-associated peptides other than those derived from the invariant chain. Our work may offer useful insights into the characteristics of proteases that are important in this process. Previous work had established that nested sets were a typical feature of peptides associated with class II MHC molecules (1, 18 22). Indeed, 4 of the 11 sets described in the present work were previously identified with a more limited number of members (1).

9 5096 ANALYSIS OF MHC CLASS II-ASSOCIATED PEPTIDE NESTED SETS FIGURE 2. Summed abundance of peptides with common amino and C termini. For each nested set, the abundances values of all peptides initiating or terminating at a specific position were summed and plotted with respect to the primary predicted binding frame identified in Table III (indicated by shaded area). The secondary predicted binding frame is indicated by a filled horizontal bar only in those instances where that frame can account for the binding of additional peptides in the nested set, as shown in Table III. However, we have shown in this study that while they tend to be dominated by a relatively small number of species, many such sets are quite complex. Previous studies have demonstrated that T cells may differentially recognize different members of such sets. For example, the crystal structure of hen egg lysozyme restricted by I-A k showed that residues 61 and 62 had very weak electron density, consistent with their being considerably solvent exposed (30). However, two-thirds of the T cells that respond to this epitope recognize the peptide 48 63, but not the minimal peptide (31). Therefore, nested sets may add a level of antigenic diversity that increases immunogenicity. It will be interesting to determine how different members of these nested sets are viewed by the immune system. References 1. Chicz, R. M., R. G. Urban, J. C. Gorga, D. A. Vignali, W. S. Lane, and J. L. Strominger Specificity and promiscuity among naturally processed peptides bound to HLA-DR alleles. J. Exp. Med. 178: Cresswell, P Assembly, transport, and function of MHC class II molecules. Annu. Rev. Immunol. 12: Watts, C Capture and processing of exogenous antigens for presentation on MHC molecules. Annu. Rev. Immunol. 15: Chapman, H. A Endosomal proteolysis and MHC class II function. Curr. Opin. Immunol. 10: Roche, P. A., and P. Cresswell Invariant chain association with HLA-DR molecules inhibits immunogenic peptide binding. Nature 345: Sloan, V. S., P. Cameron, G. Porter, M. Gammon, M. Amaya, E. Mellins, and D. M. Zaller Mediation by HLA-DM of dissociation of peptides from HLA-DR. Nature 375: Sherman, M. A., D. A. Weber, and P. E. Jensen DM enhances peptide binding to class II MHC by release of invariant chain-derived peptide. Immunity 3: Rodriguez, G. M., and S. Diment Destructive proteolysis by cysteine proteases in antigen presentation of ovalbumin. Eur. J. Immunol. 25: Rodriguez, G. M., and S. Diment Role of cathepsin D in antigen presentation of ovalbumin. J. Immunol. 149: van Noort, J. M., and M. J. Jacobs Cathepsin D, but not cathepsin B, releases T cell stimulatory fragments from lysozyme that are functional in the context of multiple murine class II MHC molecules. Eur. J. Immunol. 24: van Noort, J. M., J. Boon, A. C. van der Drift, J. P. Wagenaar, A. M. Boots, and C. J. Boog Antigen processing by endosomal proteases determines which sites of sperm-whale myoglobin are eventually recognized by T cells. Eur. J. Immunol. 21: Lang, Y., F. Forquet, E. Speck, J. Blum, and T. L. Delovitch Major histocompatibility complex class II molecules function as a template for the processing of a partially processed insulin peptide into a T-cell epitope. Diabetes 45: Deussing, J., W. Roth, P. Saftig, C. Peters, H. L. Ploegh, and J. A. Villadangos Cathepsins B and D are dispensable for major histocompatibility complex class II-mediated antigen presentation. Proc. Natl. Acad. Sci. USA 95: Manoury, B., E. W. Hewitt, N. Morrice, P. M. Dando, A. J. Barrett, and C. Watts An asparaginyl endopeptidase processes a microbial antigen for class II MHC presentation. Nature 396:695.

10 The Journal of Immunology Stern, L. J., J. H. Brown, T. S. Jardetzky, J. C. Gorga, R. G. Urban, J. L. Strominger, and D. C. Wiley Crystal structure of the human class II MHC protein HLA-DR1 complexed with an influenza virus peptide. Nature 368: Fremont, D. H., W. A. Hendrickson, P. Marrack, and J. Kappler Structures of an MHC class II molecule with covalently bound single peptides. Science 272: Dessen, A., C. M. Lawrence, S. Cupo, D. M. Zaller, and D. C. Wiley X-ray crystal structure of HLA-DR4 (DRA*0101, DRB1*0401) complexed with a peptide from human collagen II. Immunity 7: Rudensky, A. Y., P. Preston-Hurlburt, S. C. Hong, A. Barlow, and C. A. Janeway, Jr Sequence analysis of peptides bound to MHC class II molecules. Nature 353: Hunt, D. F., H. Michel, T. A. Dickinson, J. Shabanowitz, A. L. Cox, K. Sakaguchi, E. Appella, H. M. Grey, and A. Sette Peptides presented to the immune system by the murine class II major histocompatibility complex molecule I-A d. Science 256: Chicz, R. M., R. G. Urban, W. S. Lane, J. C. Gorga, L. J. Stern, D. A. Vignali, and J. L. Strominger Predominant naturally processed peptides bound to HLA-DR1 are derived from MHC-related molecules and are heterogeneous in size. Nature 358: Nelson, C. A., R. W. Roof, D. W. McCourt, and E. R. Unanue Identification of the naturally processed form of hen egg white lysozyme bound to the murine major histocompatibility complex class II molecule I-A k. Proc. Natl. Acad. Sci. USA 89: Sette, A., S. Ceman, R. T. Kubo, K. Sakaguchi, E. Appella, D. F. Hunt, T. A. Davis, H. Michel, J. Shabanowitz, R. Rudersdorf, et al Invariant chain peptides in most HLA-DR molecules of an antigen-processing mutant. Science 258: Gorga, J. C., P. J. Knudsen, J. A. Foran, J. L. Strominger, and S. J. Burakoff Immunochemically purified DR antigens in liposomes stimulate xenogeneic cytolytic T cells in secondary in vitro cultures. Cell. Immunol. 103: Shabanowitz, J., R. E. Settlage, J. A. Marto, R. E. Christian, F. M. White, P. S. Russo, S. E. Martin, and D. F. Hunt Sequencing the primordial soup. In Mass Spectrometry in Biology and Medicine. A. L. Burlingame, S. A. Carr, and M. A. Baldwin, eds. Humana, Towata, NJ, p Nepom, G. T., J. D. Lippolis, F. M. White, S. Masewicz, J. A. Marto, A. Herman, C. J. Luckey, C. Falk, J. Shabanowitz, D. F. Hunt, et al Identification and modulation of a naturally processed T cell epitope from the diabetes-associated autoantigen human glutamic acid decarboxylase 65 (hgad65). Proc. Natl. Acad. Sci. USA 98: Eng, J. K., A. L. McCormack, and J. R. Yates An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein database. J. Am. Soc. Mass Spectrom. 5: Martin, S. E., J. Shabanowitz, D. F. Hunt, and J. A. Marto Sub-femtomole MS and MS/MS peptide sequence analysis using nano-hplc micro-esi Fourier transform ion cyclotron resonance mass spectrometry. Anal. Chem. 72: Engelhard, V. H Structure of peptides associated with class I and class II MHC molecules. Annu. Rev. Immunol. 12: Hill, C. M., A. Liu, K. W. Marshall, J. Mayer, B. Jorgensen, B. Yuan, R. M. Cubbon, E. A. Nichols, L. S. Wicker, and J. B. Rothbard Exploration of requirements for peptide binding to HLA DRB1*0101 and DRB1*0401. J. Immunol. 152: Fremont, D. H., D. Monnaie, C. A. Nelson, W. A. Hendrickson, and E. R. Unanue Crystal structure of I-A k in complex with a dominant epitope of lysozyme. Immunity 8: Carson, R. T., K. M. Vignali, D. L. Woodland, and D. A. Vignali T cell receptor recognition of MHC class II-bound peptide flanking residues enhances immunogenicity and results in altered TCR V region usage. Immunity 7: Marshall, K. W., K. J. Wilson, J. Liang, A. Woods, D. Zaller, and J. B. Rothbard Prediction of peptide affinity to HLA DRB1*0401. J. Immunol. 154: Southwood, S., J. Sidney, A. Kondo, M. F. Del Guercio, E. Appella, S. Hoffman, R. T. Kubo, R. W. Chesnut, H. M. Grey, and A. Sette Several common HLA-DR types share largely overlapping peptide binding repertoires. J. Immunol. 160:3363.