PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of mcpMolecular & Cellular Proteomics : MCP
 
Mol Cell Proteomics. 2014 December; 13(12): 3367–3380.
Published online 2014 September 3. doi:  10.1074/mcp.M114.039214
PMCID: PMC4256490

Peptide Handling by HLA-B27 Subtypes Influences Their Biological Behavior, Association with Ankylosing Spondylitis and Susceptibility to Endoplasmic Reticulum Aminopeptidase 1 (ERAP1)*An external file that holds a picture, illustration, etc.
Object name is sbox.jpg

Abstract

HLA-B27 is strongly associated with ankylosing spondylitis (AS). We analyzed the relationship between structure, peptide specificity, folding, and stability of the seven major HLA-B27 subtypes to determine the role of their constitutive peptidomes in the pathogenicity of this molecule. Identification of large numbers of ligands allowed us to define the differences among subtype-bound peptidomes and to elucidate the peptide features associated with AS and molecular stability. The peptides identified only in AS-associated or high thermostability subtypes with identical A and B pockets were longer and had bulkier and more diverse C-terminal residues than those found only among non-AS-associated/lower-thermostability subtypes. Peptides sequenced from all AS-associated subtypes and not from non-AS-associated ones, thus strictly correlating with disease, were very rare. Residue 116 was critical in determining peptide binding, thermodynamic properties, and folding, thus emerging as a key feature that unified HLA-B27 biology. HLA-B27 ligands were better suited to TAP transport than their N-terminal precursors, and AS-associated subtype ligands were better than those from non-AS-associated subtypes, suggesting a particular capacity of AS-associated subtypes to bind epitopes directly produced in the cytosol. Peptides identified only from AS-associated/high-thermostability subtypes showed a higher frequency of ERAP1-resistant N-terminal residues than ligands found only in non-AS-associated/low-thermostability subtypes, reflecting a more pronounced effect of ERAP1 on the former group. Our results reveal the basis for the relationship between peptide specificity and other features of HLA-B27, provide a unified view of HLA-B27 biology and pathogenicity, and suggest a larger influence of ERAP1 polymorphism on AS-associated than non-AS-associated subtypes.

The current ideas concerning the pathogenetic role of HLA-B27 in ankylosing spondylitis (AS) emphasize specific antigen presentation (1), misfolding (2), or immunomodulation mediated by heavy chain homodimers (3) expressed at the cell surface upon endosomal recycling (4). Recent research provided evidence that both misfolded HLA-B27 heavy chains and surface expressed B27 homodimers may activate the IL-23/IL-17 axis, a key inflammatory pathway in spondyloarthropathies, through distinct mechanisms, namely the unfolded protein response (5) and the stimulation of IL-17-producing T cells (6). In contrast, the fact that CD8+ T cells are not required for the HLA-B27-associated disease in transgenic rats (7, 8), and the failure to identify specific arthritogenic peptides, point out to a pathogenetic role of HLA-B27 based on its folding and/or non-canonical forms, rather than to an autoimmune mechanism based on molecular mimicry between foreign and self-derived peptides. Yet, on the basis of genetic and immunological studies (9, 10), an involvement of CD8+T cells in the human disease cannot be ruled out.

Beyond the pathogenetic relevance of specific peptides, the constitutive HLA-B27-bound peptidome is related to the folding and stability of HLA-B27, because both features are peptide-dependent (11). This is strongly supported by the association of ERAP1, an aminopeptidase that trims peptides to their optimal size for MHC-I binding (12, 13), with ankylosing spondylitis (AS)1 among HLA-B27-positive individuals (14), and by the demonstration that AS-associated ERAP1 polymorphism has a substantial effect on the HLA-B27 peptidome in live cells (15).

Any pathogenetic mechanism must account for the differential association of HLA-B27 subtypes with AS. Whereas B*27:02, B*27:04 and B*27:05 are clearly associated with this disease, B*27:06 and B*27:09 are not (16, 17). B*27:07, a subtype present in multiple populations, is generally associated with AS, with one reported exception (18, 19). All these subtypes have the same structure in the A and B pockets of their peptide binding site, which accommodate the two N-terminal residues of their peptide ligands, but they differ in one or more positions in the F pocket, which binds the C-terminal peptide residue, as well as in other positions of the peptide binding site. In contrast, B*27:03, a subtype prevalent only in populations of Sub-Saharan African ancestry, differs from the B*27:05 prototype by a single Y59H change in the A pocket (20, 21), a difference that also sets it apart from all other subtypes (supplemental Table S1) and affects the binding preferences for N-terminal peptide residues (2224). The nature of B*27:03 as a putative susceptibility factor for AS is unclear (19). In African populations in which this subtype is prevalent, neither this subtype nor B*27:05 are associated with this disease (25), presumably because of concurrent protective factor(s).

In this study we carried out an extensive sequence analysis of HLA-B27 subtype-bound peptidomes to define their differential features as well as the extent and nature of peptide sharing among subtypes. The results revealed the basis for the intimate relationship between peptide specificity, folding, and stability of HLA-B27, provided a unified explanation on how subtype polymorphism alters the molecular biology of HLA-B27 and its association with AS, and demonstrated a differential influence of TAP and ERAP1 on AS-associated and non-AS-associated subtypes.

EXPERIMENTAL PROCEDURES

Cell Lines and Antibodies

HMy2.C1R (C1R) is a human lymphoid cell line expressing low levels of its endogenous HLA-I antigens (26). C1R transfectants expressing HLA-B*27:02, 03, 04, 05, 06, 07, or 09 (24, 2730) were used. Cells were cultured in RPMI 1640 medium supplemented with 10% fetal bovine serum (FBS) (both from Invitrogen, Paisley, UK). The W6/32 monoclonal antibody (IgG2a; specific for a monomorphic HLA class I determinant) (31) was used for immunopurification of HLA-B27.

Isolation of HLA-B27-bound Peptidomes

This was carried out as previously described (32) with some modifications. About 5 × 109 cells were lysed in 1% Igepal CA-630 (Sigma-Aldrich, St Louis, MO) in the presence of a mixture of protease inhibitors (Roche, Mannheim, Germany). HLA-B27/peptide complexes were purified from the soluble fraction by affinity chromatography using the W6/32 antibody and eluted with 0.1% trifluoroacetic acid (TFA). HLA-B27 ligands were isolated by filtering with Centricon 3 (Amicon, Beberly, MA), concentrated in SpeedVac and subjected to reverse phase HPLC chromatography in a Waters Alliance 2690 instrument (Waters, Milford, MA) using a Vydac 218TP52-C18 column (Vydac, Hesperia, CA) as previously described (33). The eluted peptides were collected and stored at −20 °C.

Electrospray-Orbitrap Mass Spectrometry Analysis

Peptide sequencing by Electrospray-Orbitrap was performed exactly as described elsewhere (34) by μLC-MS/MS using an Orbitrap XL mass spectrometer (Thermo Fisher, San Jose, CA) fitted with a capillary HPLC (Eksigent, Dublin, CA).

Electrospray-Q-TOF Mass Spectrometry Analysis

The peptide pools eluted from B*27:03, B*27:06, or B*27:09 were additionally analyzed in a nano-LC Ultra HPLC (Eksigent, AB Sciex, Redwood City, CA) coupled online with a 5600 triple TOF mass spectrometer (AB Sciex) through a nanospray III source (AB Sciex) exactly as described elsewhere (35).

Database Searches

Peptide sequences obtained by Electrospray-Orbitrap were identified as follows. Pep-Miner (36) was used for generation of the peak-lists based on the μLC-MS/MS data. The peptides were identified using multiple search engines: Pep-Miner, Sequest (Thermo-Fisher) (37) and Mascot (server 2.2, Matrix Science Inc. Boston, USA) (38), and searched against the human part of the Uniprot database (http://www.uniprot.org, May 2011) including 20,381 proteins. The Sequest and Mascot search results were combined into one report by Proteome Discoverer 1.3 (Thermo-Fisher). The search was not limited by enzymatic specificity, the peptide tolerance was set to 0.01 Da and the fragment ion tolerance was set to 0.5 Da. Oxidized M was searched as a variable modification. Peptide identifications were accepted if their masses were below 1500 Da and mass accuracy was better than 0.005 Da. Other required criteria were: Pep-Miner Score ≥85, Sequest Xcore 2.0, Sequest Probability ≤0.05. Some peptides were also observed with different charge-states and as overlapping, longer or shorter peptides by one amino acid, listed in the table as variants. The false discovery rate (FDR) of the peptides adhering to four or more of these criteria was <5%. Similar data analysis was also carried out with the MaxQuant software (39). MaxQuant (version 1.4.1.2) was used with the Andromeda search engine and the human section of the UniProt-proteome/Swiss-Prot database (release 2013_12) including 69,001 proteins. Peptide precursor and fragment mass tolerance were set at 6 ppm and 0.5 Da, respectively. The FDR was set to <5%.

In the experiments carried out using the 5600 triple TOF mass spectrometer, raw data were processed with the Protein Pilot software (version 4.5) to generate a recalibrated mgf file that was used in a MS/MS search with MASCOT (version 2.4) as search engine. A concatenated target-decoy protein database was generated by combining the Uniprot complete human proteome set (downloaded on May 23, 2011) with its corresponding reversed database generated with the DBToolKit software (version 4.1.4). Search parameters were set as follows: No enzyme, peptide tolerance: 10 ppm, MS/MS tolerance: 25 mDa, instrument: ESI-QUAD-TOF, and variable modifications: dioxidation of C, oxidation of M and pyroglutamic acid formation from N-terminal Q. Only the peptide sequences that fitted the canonical B27 binding motif (R or Q at P2) with a FDR < 5% at the peptide level were considered. FDR was estimated by decoy hit counting as previously described (40, 41).

Statistical Analyses

Except when otherwise stated the χ2 test with Bonferroni's correction was used. Statistical significance was assessed at p < 0.05.

Analysis of Protein Features

Molecular weights and isoelectric point (pI) of the precursor proteins and the human proteome (UniProtKB database release 2011.1: 2012/11/28) were obtained with the Expasy compute pI/Mw tool (http://web.expasy.org/compute_pi/). Subcellular locations were obtained by using the Protein Information and Knowledge Extractor tool (42).

TAP-binding Affinity

This was estimated for HLA-B27 ligands and their N-terminally extended precursors using the TAPREG tool (43), which makes use of support vector machines trained with large peptide sets to predict the TAP binding affinity of peptides of eight to 16 residues. Statistical significance between peptide sets was assessed with the Mann-Whitney test.

Estimation of ERAP1 Sensitivity

ERAP1 susceptibility scores for N-terminal flanking and P1 residues were calculated as previously described (15). The method assigned a score, ranging from 0 to 100, to each amino acid, based on its susceptibility to ERAP1, as established in experimental trimming assays (44).

RESULTS

HLA-B27-bound Peptidomes and their Parental Proteins

The constitutive peptidomes of B*27:02, B*27:03, B*27:04, B*27:05, B*27:06, B*27:07, and B*27:09 were subjected to MS/MS sequencing using Electrospray-Orbitrap MS and both Proteome Discoverer and MaxQuant-based peptide identifications. The peptidomes of B*27:03, B*2706, and B*27:09 were additionally subjected to Triple TOF MS analyses. Only peptides up to 18 residues long and showing either R or Q at position 2 were included in this study. With these data 2229 peptide sequences (supplemental Table S2) were compiled from the following four sources: 1) electrospray-orbitrap MS analyzed using Pep-Miner, Sequest, Mascot, and the Protein Discoverer software (supplemental Table S3), 2) electrospray-orbitrap MS data re-analyzed using MaxQuant, which significantly increased the number of assignments both of new ligands and of known ligands to additional subtypes (supplemental Table S4), 3) B*27:03, B*27:06, and B*27:09 ligands identified by electrospray-Q-TOF MS (supplemental Table S5), and 4) previously published ligands from multiple subtypes. This approach precludes quantitative comparisons and leaves open the possibility that peptides not found in a given subtype may still bind to it in some amount. Yet, we were able for the first time to identify sufficiently high numbers of peptides from multiple HLA-B27 subtypes to perform a reliable comparison of their peptidomes and of their main distinctive features.

This total set, which was used for all subsequent analyses, encompassed 1523 novel HLA-B27 ligands and many others that were previously reported from some of the subtypes and were found now in additional ones. These ligands arose from at least 2092 proteins, whose MW, pI, and subcellular distribution were similar to the human proteome, except for a sub-representation of ligands arising from membrane proteins (supplemental Fig. S1). The MW of this whole peptide set showed a Gaussian distribution centered on 1000–1200 Da and somewhat deformed at the high MW end, reflecting a higher tolerance for peptides of high MW, relative to small ones. A majority of the peptides were nonamers, followed by decamers and longer peptides (Fig. 1A). This pattern was similar to that found among HLA-B*07:02 ligands and only slightly different from HLA-A*02:01 (Fig. 1B1C) in a previous study (45).

Fig. 1.
Mass range (upper panel) and length distribution (lower panel) of 2229 HLA-B27-bound and other peptides. A, The HLA-B27 ligands are listed in supplemental Table S2. Five peptides with MW>1700 Da are not included in the upper panel. B, Mass range ...

Subtype-bound Peptidomes: General Features

The total number of peptide sequences from each subtype is shown in Table I. Differences were found in the mean MW (ranging from 1071.3 to 1149.7 Da) and length (9.3 to 9.8 residues) of subtype-bound peptidomes. The length of peptides from individual subtypes showed some significant differences. For instance, 9-mers ranged from 53.6/54.6% to 79.8% (p: 1.2 × 10−24/4.8 × 10−23), 10-mers from 27.0/25.0 to 11.5% (p: 4.0 × 10−12/8.2 × 10−10), and peptides >10-mers from 19.0/19.2% to 6.7% (p: 1.6 × 10−10/9.6 × 10−11) for B*27:02/B*27:05 and B*27:06, respectively. Other subtypes showed intermediate values. The percent of peptides found only in one subtype ranged from 9.1 to 24.9% (B*27:07/B*27:04), indicating a large overlap among subtype-bound peptidomes.

Table I
HLA-B27 subtype-bound ligands

Subtype-bound Peptidomes: Residue Usage

The comparative analysis of residue frequencies (RF) among subtype-bound peptide repertoires was carried out separately for main anchor (P2 and PC), secondary anchor (P1, P3, and PC-2) and “nonanchor” positions (P4-P6, and P8). N-fold and statistical significance of relevant differences are shown in supplemental Table S6.

Main Anchor Positions (Fig. 2)
Fig. 2.
Residue frequencies (% RF) among subtype-bound HLA-B27 ligands at their main anchor positions P2 and PC. AS-associated (B*27:02, 04, 05, and 07) and non-AS-associated (B*27:06 and 09) subtypes are shown, in that order, with white and black bars, respectively. ...

P2

As expected from the identity of the B pocket, the overwhelming majority of peptides from any given subtype had R2. In addition, Q2, previously reported only in B*27:05 (4648), was found with low frequency (1.8 to 3.3%) in all subtypes.

PC

F and L were among the most frequent C-terminal residues in all subtypes (RF: 15.4–29.8% and 23.3–57.1%, respectively). Other aliphatic residues were less frequent. Y was frequent in most AS-associated subtypes (B*27:02, 04, 05) and B*27:03 (17.8–30.1%), but rare in B*27:06, 07, and 09 (1.9–2.4%). Basic residues (K, R) were frequent only in B*27:03 and 05 (7.3–10.7%). Other residues were rare or absent altogether. Significant differences in both the chemical nature and size of the C-terminal peptide residues were observed among subtypes, as a function of residue 116 (Table II and supplemental Table S6). Subtypes lacking D116 (B*27:06, 07, and 09) showed an increased frequency of aliphatic residues relative to subtypes with D116 (67.9–78.7% versus 35.2–50.2%), an almost absolute restriction (>90%) to nonpolar aliphatic residues and F, and a preference for small residues (104–115 Da: 52.1–63.7%). Subtypes with D116 were permissive to C-terminal Y (17.7–30.1%) and, some of them (B*27:03, 05), to basic residues (21.9 and 17.4%, respectively), and showed an increased preference (45.6–62.6%) for bulky (>138 Da) C-terminal residues.

Table II
Polarity and size of C-terminal residues among HLA-B27 subtype-bound peptides. The nature of residue 116 in each subtype is indicated. Figures are residue frequencies (%) among the peptides sequenced from the corresponding subtype. Values selectively ...

There was a close correlation between the structure, folding, and thermostability of HLA-B27 subtypes with identical A and B pockets, therefore excluding B*2703, and their C-terminal motifs (Table III). These subtypes could be classified in two groups: (1) D116/slow folding/high thermostability/allowance for Y/preference for bulky C-terminal residues (B*27:02, 04, and 05), and (2) non-D116/fast folding/low thermostability/high restriction to nonpolar C-terminal residues, almost excluding Y/preference for C-terminal residues of moderate size (B*27:06, 07, and 09).

Table III
Summary of HLA-B27 subtype features

Secondary Anchor Positions (Fig. 3)
Fig. 3.
Residue frequencies (% RF) among subtype-bound HLA-B27 ligands at their secondary anchor positions P1, P3, and PC-2. B*27:02, 04, 05 and 07 are shown, in that order, with white bars, B*27:06 and 09 are shown with black bars and B*27:03 in shown with a ...

P1

The most frequent residues at this position were A, G, H, K, and S. They showed comparable frequencies in all subtypes except B*27:03, which had statistically increased H and K, and decreased G and A frequencies, relative to other subtypes. R was also frequent and most abundant in B*27:03 and 04 (supplemental Table S6). Other residues were used with smaller and similar frequencies in all subtypes.

P3

The most frequent residues at this position were L and F, although many other residues were allowed. B*27:04 and B*27:06 showed slightly increased frequencies of basic P3 residues (R, K), relative to other subtypes. In contrast, B*27:03 showed a particularly low allowance for K and H (supplemental Table S6).

PC-2

Restrictions at this position were low, allowing for both polar and nonpolar side chains. Residue usage was generally similar among subtypes, with some exceptions. For instance, L and K were statistically increased in B*27:04 and B*27:06, respectively, R was increased and E was decreased in both of these subtypes, V was decreased in B*27:06, relative to most other subtypes (supplemental Table S6).

“Nonanchor” Positions (supplemantal Fig. S2)

Because a comparison of the central peptide positions requires a precise alignment of peptide sequences, this analysis can only be properly done with peptides of the same length and was carried out only for positions P4–P6 and P8 of nonamers. Restrictions were low at all four positions, but each one showed some statistically significant bias toward particular residues as well as individual subtype preferences (supplemental Table S6). The most frequent residues at P4 were P (15.9–20.8%) and G (12.4–14.0%) for all subtypes. At P5 the most frequent residues were G (7.2–11.3%), L (8.2–11.4%), and V (5.8–12.6%); the two latter residues were less frequent in B*27:06 and B*27:04, respectively. In contrast, R5 was more frequent in these (2.7–4.6%) than in other subtypes (0.3–1.8%). At P6 the most frequent residues were I, L, V, and polar residues (i.e. Q, S), with RF≈10% each in all subtypes at P8 the most frequent residues were A, E, G, K, L, N, Q, S, and T. Some differences among subtypes were observed, most notably concerning the lower frequency of N (0.5–0.6-fold) and higher frequency of V (about twofold) in B*27:04 and B*27:06. These results indicate that subtype polymorphism has an influence on residue usage also at nonanchor positions.

Peptide Sharing Among HLA-B27 Subtypes

In previous studies we carried out a series of pairwise comparisons between subtype-bound peptide repertoires (29, 30, 49, 50). These studies were based on identity of MW, chromatographic retention time, and limited sequencing. Here, we re-assessed the overlap among subtype-bound peptidomes based on the many novel ligands sequenced from distinct subtypes (supplemental Table S7). All subtype pairs showed similar levels of peptide sharing (about 60–74%). Of 2229 HLA-B27 ligands, 1008 (45.2%) were found only in one subtype and the remaining peptides showed highly diverse patterns of subtype distribution (supplemental Table S2).

Peptides Selectively Presented by AS-associated or non-AS-associated Subtypes

A total of 1167 ligands were identified from one or more AS-associated subtypes, but were not found in the non-AS-associated ones. Likewise, 227 peptides were found only in the non-AS-associated B*27:06 and/or B*27:09 subtypes (Table IV). The mean length and MW of the former peptide set were bigger by 0.5 residues and more than 100 Da, respectively, relative to the latter one. A clear shift was observed in the modal distribution of MW and length between both peptide sets: Heavier and longer ligands were more abundant in the AS-associated subtypes (Fig. 4A-B). Of note, 9-mers accounted for 80.2% of the peptides selectively found among non-AS-associated subtypes compared with only 49.5% in AS-associated ones (p = 1.6 × 10−16).

Table IV
Length and MW of peptide ligands from subtypes with differential AS association or thermostablity
Fig. 4.
MW and length of peptide subsets correlating with AS or molecular stability. A, Mass range distribution of peptides found exclusively in AS-associated (dashed line) or non-AS-associated (solid line) HLA-B27 subtypes. B, Length distribution of peptides ...

Significant differences in residue usage at the C-terminal position, concerning both polarity and size, were also found between both subsets (Fig. 5A5C). Bulky aromatic and basic residues were much more frequent in the AS-associated set. In contrast, small and aliphatic residues, especially L, prevailed among the peptides found only in non-AS-associated subtypes.

Fig. 5.
C-terminal residue usage among peptide subsets correlating with AS or molecular stability. A, Residue frequencies (% RF) among peptides found exclusively in the AS-associated subtypes B*27:02, 04, 05, and 07 (white bars) or the non-AS-associated ones, ...

These results indicate that AS-associated and non-AS associated subtypes selectively bind peptides of distinct length and MW, which are distinguished by differential use of their prevalent C-terminal residues.

Peptides Selectively Presented by HLA-B27 Subtypes with High or Low Thermostability

A total of 988 and 313 ligands were found only among high (B*27:02, 04, and 05) or low (B*27:06, 07, and 09) thermostability subtypes, respectively (Table IV). B*27:03 was not included in this analysis because its low thermostability is due to a change in the A pocket and, therefore, has a different basis relative to the other subtypes. The peptides found only in the high thermostability subtypes had higher mean length (by 0.4 residues) and MW (by 105 Da) than those found only in the low thermostability subtypes (Table IV) and the MW and length distribution of both subsets (Fig. 4C4D) indicated that peptides heavier and longer were more abundant in the former one. Bulky aromatic and basic C-terminal residues were preferentially found in high thermostability subtypes, whereas L and other small aliphatic residues were predominant in the low thermostability subset (Fig. 5D5F).

These results indicate that the differential features between peptides associated with high or low thermostability are similar to those distinguishing the peptides found only among AS-associated or non-AS-associated subtypes. High/low thermostability and AS-/non-AS-associated subtype groups differ only by B*27:07, which is AS-associated but shows low thermostability. Because of this, the peptides found only among high thermostability subtypes were less than those found only among AS-associated subtypes (988 compared with 1167), and the peptides found only among the low thermostability subtypes were increased relative to non-AS-associated subtypes (313 compared with 227). Yet, this was not sufficient to alter the similarity of peptide features between AS-associated/high-thermostability or non-AS-associated/low-thermostability subtypes.

Peptides Strictly Correlating with AS Association or Thermostability are Rare

A total of 26 peptides, mostly with C-terminal F, were found in all four AS-associated subtypes and not in non-AS-associated ones. In addition, 22 peptides, all with C-terminal aliphatic residues, were found in both B*27:06 and 09 and not among AS-associated subtypes. As many as 125 ligands were exclusively found in all three high thermostability subtypes. Their predominant C-terminal motif (88%) was Y. Only two ligands, both with C-terminal L, were exclusively shared by the three subtypes with low thermostability (supplemental Table S8).

HLA-B27 Ligands are Better Suited to TAP than their Precursors

To determine how the last stages of antigen processing contribute to shape HLA-B27-bound peptidomes, the N-terminal flanking residues of these ligands in their parental proteins were analyzed for their contribution to TAP-mediated transport and to ERAP1 trimming.

The affinity of peptides for TAP depends on their whole sequence (43), but mainly on the three N-terminal and the C-terminal residues (51). The TAP-binding affinity of HLA-B27 ligands and their N-terminally extended precursors, by 1 to 4 residues, was compared using a predictive algorithm that takes into account the contribution of all peptide residues (43). The final HLA-B27 ligands were significantly more suited to TAP than any of their N-terminally extended precursors. Among these, the best ones were those extended by only one residue (Fig. 6A). This pattern is different to that found among HLA-B*07:02 or HLA-A*02:01 ligands (Fig. 6C6D) reported in a previous study (45). Our results strongly suggest that natural HLA-B27 ligands that are generated in their final form in the cytosol will reach the ER more efficiently than their precursors.

Fig. 6.
Predicted affinity of HLA-B27 and other ligands and their precursors for TAP. A, TAP binding affinity of the natural HLA-B27 ligands in supplemental Table S2 and their N-terminally extended precursors (Nt) by 1 to 4 residues. Bars indicate the geometrical ...

Peptides Found Only Among AS-associated Subtypes are the Best TAP Ligands

The natural ligands selectively found among AS-associated subtypes were significantly better suited to TAP than those found only in non-AS-associated ones (Fig. 6B). Because the former set has more diverse C-terminal motifs (Fig. 5A), and these are major contributors to TAP affinity, we considered the possibility that this diversity could account for the better suitability to TAP. Thus, the peptides found only among AS-associated subtypes were classified on the basis of their C-terminal residues and their TAP affinity was estimated separately for each subset (Fig. 6B). The results showed that the higher affinity of the AS-associated peptide set to TAP is mainly caused by the peptides with aromatic C-terminal residues, although peptides with aliphatic or basic residues were also better suited to TAP than peptides found only among non-AS-associated subtypes. A large majority of these peptides had aliphatic C-terminal residues (Fig. 5A). Thus, the fact that AS-associated peptides with aliphatic C-terminal residues were still better TAP binders indicates that residues at other peptide positions contribute to the higher affinity of this peptide subset.

Distinct Influence of ERAP1 on the Peptidomes of Subtypes with Differential AS Association and Thermostability

ERAP1 polymorphism is associated with AS (14, 52) and lack of this enzyme affects the surface expression of AS-associated, but not of non-AS-associated HLA-B27 subtypes (53). Thus, we compared the global susceptibility to ERAP1 of the two N-terminal flanking (P-2 and P-1) and the P1 residues of peptide ligands found only among AS-associated subtypes with those found only in non-AS-associated ones. This comparison focused on the frequency of the residues highly susceptible to ERAP1 (score >50: Y, M, L, A, and C) and those with the lowest susceptibility (score <10: E, D, V, W, K, R, and P). The joint frequency of susceptible flanking and P1 residues was similar in both peptide subsets (Fig. 7A). However, peptides found only in AS-associated subtypes showed a statistically significant increase of ERAP1-resistant P1 residues relative to non-AS-associated subtypes (Fig. 7B). When the same analysis was carried out for peptides selectively found among either high or low thermostability subtypes a similar tendency was found (Fig. 7C7D). These differences were due mainly to the 9-mers and shorter peptides (Fig. 7E7F).

Fig. 7.
ERAP1 susceptibility of the flanking and P1 residues of peptides selectively found among HLA-B27 subtypes with differential AS-association or thermostability. A, Joint frequency of highly susceptible residues (score>50) at the P-2, P-1, and P1 ...

In addition, there was an increased frequency of resistant versus susceptible P1 residues among peptides selectively found in AS-associated and high thermostability subtypes, (ratio 1.5). The differences were increased (ratio 1.9) when only peptides ≤9-mers were considered. In contrast, peptides from non-AS-associated or low thermostability subtypes, or the corresponding subsets of peptides ≤9-mers, showed no increased frequencies (ratio 0.7–0.8) of resistant versus susceptible P1 residues (Table V).

Table V
ERAP1 susceptibility of N-terminal residues selectively found among HLA-B27 subtypes with differential AS association or thermostability. The data represent the mean frequency (%) of residues with trimming susceptibility scores <10 (resistant), ...

Taken together, these results reveal a differential influence of ERAP1 on AS-associated/high-thermostability and non-AS-associated/low-thermostability subtypes.

DISCUSSION

The present study is, by far, the most extensive analysis of HLA-B27 subtype-bound peptidomes reported to date and allowed us, for the first time, to define the relationship between peptide specificity and other pathogenetic features of this molecule. Because subtype-bound peptides were identified with various methods and compiled together in a single database there is a heterogeneity that is likely to affect the completeness of subtype assignments. It is obvious that peptides not found in a given subtype cannot be excluded as ligands of that subtype. A more uniform and sensitive methodology would certainly show higher peptide sharing among subtypes, which should then require a quantitative assessment, because shared peptides can be presented by different subtypes at widely different levels. A quantitative comparison of subtype-bound peptidomes was not the purpose and is out of the scope of this study.

An important issue is the use of a 5% FDR threshold for peptide assignments. It can be argued that it is a rather relaxed one and may introduce an excessive number of false positives. It has to be noted, however, that in MHC peptidomics MS/MS searches are typically conducted without enzyme restriction. This leads to a huge increase of the search space compared with a standard tryptic workflow. Under these circumstances, the chances of random matching increase and high-scoring decoys lead to a huge loss in sensitivity if peptide matches are filtered at 1% FDR, an issue discussed at length elsewhere (45). Additionally, in this scenario, a FDR of 1% may be reached with a very low number of decoy hits, compromising the accuracy of FDR estimations. For these reasons FDR filters are typically less stringent, and a 5% FDR is currently standard, in MHC peptidomics. We are aware that this threshold is rather lean and will introduce a certain fraction of false positives. Nevertheless, this should be largely compensated by avoiding the great loss of true positives occurring at more stringent FDR. Improvements in the accuracy of MS instruments, and in search and scoring algorithms, should facilitate in a near future a more robust FDR assessment, in line with that of trypsin-based proteomics.

With the above limitations in mind, HLA-B27 subtypes showed three peptide specificity patterns according to their preferences for C-terminal residues: (1) B*27:03 and B*27:05 preferentially bound peptides with basic, aliphatic and aromatic residues, (2) B*27:02 and B*27:04 were restricted to C-terminal aliphatic/aromatic residues with a low allowance for basic ones, and (3) B*27:06, 07, and 09 were restricted to nonpolar aliphatic and F residues, almost excluding Y. Subtypes from groups 1 and 2 showed a preference for large C-terminal residues, whereas those from group 3 preferred smaller ones. The exclusion of polar and large C-terminal residues in B*27:06, 07, and 09 is because of the substitution of D116 in these subtypes for non-acidic and bulkier residues (Y, H) that decrease the polarity and size of the F pocket. The low preference for C-terminal basic residues, but not Y, in B*27:02 and 04 is presumably determined by the change of D77 to small polar residues (N, S) in these subtypes (supplemental Table S1).

The comparable peptide sharing among subtypes binding either very distinct (i.e. B*27:05/07) or similar (i.e. B*27:06/09) C-terminal residues indicates that differences at other peptide positions, which are determined by polymorphism outside the F pocket, introduce significant diversity among subtype-bound peptidomes. For instance, the increased frequencies of basic P3 and PC-2 residues, and Y3, in B*27:04 and B*27:06, are presumably because of the V152E change in both subtypes and to the H114D change in B*27:06. The high frequency of N-terminal basic residues in B*27:03, which extends previous observations (22, 23), is thought to be a compensatory effect of interactions between basic P1 side chains and E163, for the altered hydrogen bonding of the peptidic N terminus in the A pocket of this subtype (54).

The relationship between the peptidomes and other features of HLA-B27 subtypes revealed the key role of residue 116 in unifying HLA-B27 biology. Among subtypes with identical A and B pockets, absence of D116 correlated strictly with decreased molecular stability, fast folding and almost exclusive binding of peptides with C-terminal nonpolar residues. This is likely the result of an increased hydrophobicity of the F pocket because of abrogation of the acidic charge. The thermodynamic requirements to hide a hydrophobic F pocket from water during the biosynthesis of HLA-B27 would promote fast peptide binding, decrease tapasin dependence, as shown for B*27:06 (55), and lead to fast folding with a suboptimal peptide repertoire. A similar situation was observed among HLA-B44 subtypes (56).

The correspondence of D116/slow folding/high thermostability with AS is incomplete because B*27:07 lacks D116 and shows fast folding and decreased thermostability (57, 58). Yet, the features of the peptides segregating with AS susceptibility or differential thermostability were very similar. This might be related the fact that the B*27:07 has a comparable peptide sharing with B*27:05 and B*27:06/B*27:09. Like B*27:07, B*14:03, which is associated with AS in African populations (19, 59, 60), has lower thermostability than AS-associated HLA-B27 subtypes and equally low as B*14:02, which is not associated with AS (61). The F pocket of B*14:03 is also hydrophobic and almost exclusively binds peptides with nonpolar C-terminal residues (62). Thus, the association of B*27:07 and B*14:03 with AS cannot be explained just as a consequence of their folding and stability features, and its basis remains unknown. This conclusion raises the possibility that, in a more general way, the association of HLA-B27 with AS may be dependent on additional peptide-dependent mechanisms besides those based on folding and stability.

The small number of peptides found in all four AS-associated subtypes and not otherwise might be compatible with a contribution of specific peptides to AS pathogenesis. Yet, the many more peptides found only among AS-associated subtypes, but with diverse distribution patterns, are presumably more relevant because they could globally affect other features of the molecule, such as folding and stability.

It was previously noted that HLA-B27 ligands are particularly well suited to TAP preferences (63). Now we showed that these ligands, irrespective of the subtype(s) in which they were found, are better suited to TAP than their respective N-terminally extended precursors. In addition, peptides with N-terminal dibasic sequences, which are frequent in HLA-B27, are relatively resistant to cytosolic aminopeptidases (64). Thus, HLA-B27 is particularly well suited to bind peptides directly produced in the cytosol. The observation that HLA-B27 ligands selectively found among AS-associated subtypes were better suited to TAP than those from non-AS-associated ones raises the possibility of a relationship between the direct import of fully processed HLA-B27 ligands from the cytosol and AS association, and is consistent with the reported involvement of the cytosolic puromycin-sensitive aminopeptidase in AS (65). Yet, the significant influence of ERAP1 on the HLA-B27 peptidome (15, 66) implies that many natural ligands are generated and/or destroyed in the ER.

Two features distinguishing the peptides selectively found among AS-associated/high thermostability or non-AS-associated/low thermostability subtypes revealed a differential influence of ERAP1 on both subtype groups: (1) the distinct susceptibility of their P1 residues to ERAP1 trimming, and (2) their distinct length distribution. As discussed above, the fast folding and export rates of non-AS-associated/low-thermostability subtypes presumably result from their fast and suboptimal loading of peptides in the ER. Thus, these subtypes could load peptides with both short and long half-life, depending on their affinity and abundance. Because the half-life of peptides in the ER is strongly influenced by the susceptibility of their P1 residues to ERAP1, this explains that resistant and susceptible P1 residues are similarly frequent among the peptides selectively found in these subtypes. Moreover, in the absence of tapasin-mediated editing, peptides with optimal length might be preferred, which explains the prevalence of 9-mers in these subtypes.

In contrast, AS-associated/high-thermostability subtypes fold more slowly and bind more optimized peptidomes (57, 58) after tapasin editing (55, 56). Thus, they may preferentially bind peptides with longer half-life in the ER, which are those with P1 residues more resistant to ERAP1, in agreement with our results. Peptide editing would also allow selecting longer peptides if these bind with higher stability than 9-mers, resulting in a more balanced size distribution within the range compatible with MHC-I binding, as also observed here.

The previous observation that ERAP1 knockdown affected the expression levels of AS-associated, but not of non-AS-associated subtypes at the cell surface (53) is most likely explained by the differences in stability and ERAP1 dependence of these subtypes.

Our observations suggest a differential effect of ERAP1 polymorphism on AS- and non-AS-associated subtypes. We previously reported that B*27:04 ligands prevalent in the context of a highly active ERAP1 variant had an increased frequency of ERAP1-resistant P1 residues, relative to the predominant peptides expressed in a less active ERAP1 context (15). This was interpreted as the result of a more extensive destruction of ligands with susceptible P1 residues, relative to those with resistant ones. Thus, the lower discrimination of non-AS-associated/low-thermostability subtypes for peptides differing in the susceptibility of their P1 residues to ERAP1 suggests that these subtypes may be less influenced by ERAP1 polymorphism than AS-associated/high-thermostability ones. If so, AS-associated ERAP1 alleles might not favor disease development in B*27:06 and B*27:09 positive individuals.

In conclusion, the extensive characterization of HLA-B27 subtype-bound peptidomes reveals the structural basis of their peptide specificity and its close relationship to folding, stability and association with AS. Residue 116 is the key that interconnects all the molecular features of HLA-B27, but other polymorphisms significantly influence the diversity of subtype-bound peptide repertoires. D116 confers an increased polarity to the F pocket and a higher plasticity in its interaction with chemically diverse and bulkier C-terminal peptide residues. These features allow for a proper optimization of the peptidome through tapasin-mediated editing (56), leading to highly stable HLA-B27 molecules and a significant dependence on ERAP1 and its functional polymorphism. Optimization of the peptide cargo is achieved through slow folding, a byproduct of which is increased heavy chain misfolding (67, 68). This behavior would be critical for the disease-predisposing capacity of HLA-B27 and explains that subtypes lacking D116 are biochemically distinct and, except B*27:07, not associated with AS. The pathogenetic relevance of heavy chain homodimers (69), whose expression at the cell surface presumably follows endosomal dissociation of the heterodimer/peptide complexes (4), might seem difficult to reconcile with the lower thermostability of non-AS-associated subtypes, unless one assumes that the longer half-life and/or increased endosomal recycling of high-thermostability subtypes might favor homodimer formation. That B*27:09 shows less surface homodimers than B*27:05 (70) is consistent with this possibility. Thus, the close relationship between peptide specificity, folding, and stability unifies the biology of HLA-B27 and makes the peptidome, its alterations by subtype polymorphism and the differential ERAP1 dependence of subtype-bound peptidomes, critical features of the pathogenetic role of HLA-B27 in AS.

Supplementary Material

Supplemental Data:

Acknowledgments

We thank Sergio Ciordia, Silvia Juarez and Juan P. Albar (ProteoRed-ISCIII, Centro Nacional de Biotecnología, Madrid) for their help.

Footnotes

Contributed by

Author contributions: J.A.L. designed research; N.G., A.S., C.A., P.G., E.B., M.M., and A.A. performed research; P.G. and A.A. contributed new reagents or analytic tools; N.G., A.S., C.A., P.G., E.B., M.M., A.A., and J.A.L. analyzed data; N.G. and J.A.L. wrote the paper; A.S. and C.A. contributed fundamental ideas; A.A. provided critical immunoproteomics input; J.A.L. supervision and direction of the research work.

* This work was supported by grants SAF2011/25681 (Plan Nacional de I+D+i), and RD08/0075 (RIER-ISCIII) to JALC, a Binational Science Foundation Grant 2009393 to AA, and an institutional grant of the Fundación Ramón Areces to the Centro de Biología Molecular Severo Ochoa. The work of MM was supported by grant PT13/0001 (Instituto de Salud Carlos III) to Juan P. Albar.

An external file that holds a picture, illustration, etc.
Object name is sbox.jpg This article contains supplemental Figs. S1 and S2 and Tables S1 to S8.

1 The abbreviations used are:

AS
ankylosing spondylitis
C1R
HMy2.C1R
FDR
false discovery rate
RF
residue frequency
TFA
trifluoroacetic acid.

REFERENCES

1. Benjamin R., Parham P. (1990) Guilt by association: HLA-B27 and ankylosing spondylitis. Immunol. Today 11, 137–142 [PubMed]
2. Colbert R. A., Tran T. M., Layh-Schmitt G. (2014) HLA-B27 misfolding and ankylosing spondylitis. Mol. Immunol. 57, 44–51 [PMC free article] [PubMed]
3. Shaw J., Hatano H., Kollnberger S. (2014) The biochemistry and immunology of non-canonical forms of HLA-B27. Mol. Immunol. 57, 52–58 [PubMed]
4. Bird L. A., Peh C. A., Kollnberger S., Elliott T., McMichael A. J., Bowness P. (2003) Lymphoblastoid cells express HLA-B27 homodimers both intracellularly and at the cell surface following endosomal recycling. Eur. J. Immunol. 33, 748–759 [PubMed]
5. DeLay M. L., Turner M. J., Klenk E. I., Smith J. A., Sowders D. P., Colbert R. A. (2009) HLA-B27 misfolding and the unfolded protein response augment interleukin-23 production and are associated with Th17 activation in transgenic rats. Arthritis Rheum. 60, 2633–2643 [PMC free article] [PubMed]
6. Bowness P., Ridley A., Shaw J., Chan A. T., Wong-Baeza I., Fleming M., Cummings F., McMichael A., Kollnberger S. (2011) Th17 cells expressing KIR3DL2+ and responsive to HLA-B27 homodimers are increased in ankylosing spondylitis. J. Immunol. 186, 2672–2680 [PMC free article] [PubMed]
7. Breban M., Fernandez-Sueiro J. L., Richardson J. A., Hadavand R. R., Maika S. D., Hammer R. E., Taurog J. D. (1996) T cells, but not thymic exposure to HLA-B27, are required for the inflammatory disease of HLA-B27 transgenic rats. J. Immunol. 156, 794–803 [PubMed]
8. May E., Dorris M. L., Satumtira N., Iqbal I., Rehman M. I., Lightfoot E., Taurog J. D. (2003) CD8αβ T cells are not essential to the pathogenesis of arthritis or colitis in HLA-B27 transgenic rats. J. Immunol. 170, 1099–1105 [PubMed]
9. Robinson P. C., Brown M. A. (2014) Genetics of ankylosing spondylitis. Mol. Immunol. 57, 2–11 [PubMed]
10. Sorrentino R., Bockmann R. A., Fiorillo M. T. (2014) HLA-B27 and antigen presentation: at the crossroads between immune defense and autoimmunity. Mol. Immunol. 57, 22–27 [PubMed]
11. Marcilla M., Lopez de Castro J. A. (2008). Peptides: the cornerstone of HLA-B27 biology and pathogenetic role in spondyloarthritis. Tissue Antigens 71, 495–506 [PubMed]
12. Serwold T., Gonzalez F., Kim J., Jacob R., Shastri N. (2002) ERAAP customizes peptides for MHC class I molecules in the endoplasmic reticulum. Nature 419, 480–483 [PubMed]
13. Saric T., Chang S. C., Hattori A., York I. A., Markant S., Rock K. L., Tsujimoto M., Goldberg A. L. (2002) An IFN-gamma-induced aminopeptidase in the ER, ERAP1, trims precursors to MHC class I-presented peptides. Nat. Immunol. 3, 1169–1176 [PubMed]
14. The TASK and WTCCC2 Consortia (2011). Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility. Nat. Genet. 43, 761–767 [PMC free article] [PubMed]
15. Garcia-Medel N., Sanz-Bravo A., Van Nguyen D., Galocha B., Gomez-Molina P., Martin-Esteban A., Alvarez-Navarro C., López de Castro J. A. (2012). Functional interaction of the ankylosing spondylitis-associated endoplasmic reticulum aminopeptidase 1 polymorphism and HLA-B27 in vivo. Mol. Cell. Proteomics 11, 1416–1429 [PMC free article] [PubMed]
16. Lopez-Larrea C., Sujirachato K., Mehra N. K., Chiewsilp P., Isarangkura D., Kanga U., Dominguez O., Coto E., Peña M., Setien F., Gonzalez-Roces S. (1995) HLA-B27 subtypes in Asian patients with ankylosing spondylitis. Evidence for new associations. Tissue Antigens 45, 169–176 [PubMed]
17. D'Amato M., Fiorillo M. T., Carcassi C., Mathieu A., Zuccarelli A., Bitti P. P., Tosi R., Sorrentino R. (1995) Relevance of residue 116 of HLA-B27 in determining susceptibility to ankylosing spondylitis. Eur. J. Immunol. 25, 3199–3201 [PubMed]
18. Varnavidou-Nicolaidou A., Karpasitou K., Georgiou D., Stylianou G., Kokkofitou A., Michalis C., Constantina C., Gregoriadou C., Kyriakides G. (2004) HLA-B27 in the Greek Cypriot population: distribution of subtypes in patients with ankylosing spondylitis and other HLA-B27-related diseases. The possible protective role of B*2707. Hum. Immunol. 65, 1451–1454 [PubMed]
19. Diaz-Pena R., Lopez-Vazquez A., Lopez-Larrea C. (2012) Old and new HLA associations with ankylosing spondylitis. Tissue Antigens 80, 205–213 [PubMed]
20. Rojo S., Aparicio P., Hansen J. A., Choo S. Y., Lopez de Castro J. A. (1987). Structural analysis of an HLA-B27 functional variant, B27d, detected in American blacks. J. Immunol. 139, 3396–3401 [PubMed]
21. Choo S. Y., St. John T., Orr H. T., Hansen J. A. (1988) Molecular analysis of the variant alloantigen HLA-B27d (HLA- B*2703) identifies a unique single amino acid substitution. Hum. Immunol. 21, 209–219 [PubMed]
22. Colbert R. A., Rowland Jones S. L., McMichael A. J., Frelinger J. A. (1994) Differences in peptide presentation between B27 subtypes: the importance of the P1 side chain in maintaining high affinity peptide binding to B*2703. Immunity 1, 121–130 [PubMed]
23. Boisgérault F., Tieng V., Stolzenberg M. C., Dulphy N., Khalil I., Tamouza R., Charron D., Toubert A. (1996) Differences in endogenous peptides presented by HLA-B*2705 and B*2703 allelic variants: implications for susceptibility to spondylarthropathies. J. Clin. Invest. 98, 2764–2770 [PMC free article] [PubMed]
24. Villadangos J. A., Galocha B., Lopez D., Calvo V., Lopez de Castro J. A. (1992). Role of binding pockets for amino-terminal peptide residues in HLA-B27 allorecognition. J. Immunol. 149, 505–510 [PubMed]
25. Brown M., Jepson A., Young A., Whittle H., Greenwood B., Wordsworth P. (1997) Spondyloarthritis in West Africa: evidence for a non-B27 protective effect. Ann. Rheum. Dis. 56, 68–70 [PMC free article] [PubMed]
26. Zemmour J., Little A. M., Schendel D. J., Parham P. (1992) The HLA-A,B “negative” mutant cell line C1R expresses a novel HLA-B35 allele, which also has a point mutation in the translation initiation codon. J. Immunol. 148, 1941–1948 [PubMed]
27. Calvo V., Rojo S., Lopez D., Galocha B., Lopez de Castro J. A. (1990). Structure and diversity of HLA-B27-specific T cell epitopes. Analysis with site-directed mutants mimicking HLA-B27 subtype polymorphism. J. Immunol. 144, 4038–4045 [PubMed]
28. Villadangos J. A., Galocha B., Lopez de Castro J. A. (1994). Unusual topology of an HLA-B27 allospecific T cell epitope lacking peptide specificity. J. Immunol. 152, 2317–2323 [PubMed]
29. Ramos M., Paradela A., Vazquez M., Marina A., Vazquez J., Lopez de Castro J. A. (2002). Differential association of HLA-B*2705 and B*2709 to ankylosing spondylitis correlates with limited peptide subsets but not with altered cell surface stability. J. Biol. Chem. 277, 28749–28756 [PubMed]
30. Gomez P., Montserrat V., Marcilla M., Paradela A., López de Castro J. A. (2006). B*2707 differs in peptide specificity from B*2705 and B*2704 as much as from HLA-B27 subtypes not associated to spondyloarthritis. Eur. J. Immunol. 36, 1867–1881 [PubMed]
31. Barnstable C. J., Bodmer W. F., Brown G., Galfre G., Milstein C., Williams A. F., Ziegler A. (1978) Production of monoclonal antibodies to group A erythrocytes, HLA and other human cell surface antigens. New tools for genetic analysis. Cell 14, 9–20 [PubMed]
32. Paradela A., Garcia-Peydro M., Vazquez J., Rognan D., Lopez de Castro J. A. (1998). The same natural ligand is involved in allorecognition of multiple HLA-B27 subtypes by a single T cell clone: role of peptide and the MHC molecule in alloreactivity. J. Immunol. 161, 5481–5490 [PubMed]
33. Paradela A., Alvarez I., Garcia-Peydro M., Sesma L., Ramos M., Vazquez J., Lopez de Castro J. A. 2000. Limited diversity of peptides related to an alloreactive T cell epitope in the HLA-B27-bound peptide repertoire results from restrictions at multiple steps along the processing-loading pathway. J. Immunol. 164, 329–337 [PubMed]
34. Garcia-Medel N., Sanz-Bravo A., Barnea E., Admon A., Lopez de Castro J. A. (2012). The origin of proteasome-inhibitor resistant HLA class I peptidomes: a study with HLA-A*68:01. Mol. Cell. Proteomics 11, 1–15 [PMC free article] [PubMed]
35. Marcilla M., Alpizar A., Lombardia M., Ramos-Fernandez A., Ramos M., Albar J. P. (2014) Increased diversity of the HLA-B40 ligandome by the presentation of peptides phosphorylated at their main anchor residue. Mol. Cell. Proteomics 13, 462–474 [PMC free article] [PubMed]
36. Beer I., Barnea E., Ziv T., Admon A. (2004) Improving large-scale proteomics by clustering of mass spectrometry data. Proteomics 4, 950–960 [PubMed]
37. Eng J. K., McCormack A. L., Yates J. R. (1994) An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein database. J. Amer. Soc. Mass. Spect. 5, 976–989 [PubMed]
38. Perkins D. N., Pappin D. J., Creasy D. M., Cottrell J. S. (1999) Probability-based protein identification by searching sequence databases using mass spectrometry data. Electrophoresis 20, 3551–3567 [PubMed]
39. Cox J., Mann M. (2008) MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 26, 1367–1372 [PubMed]
40. Elias J. E., Gygi S. P. (2007) Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry. Nat. Methods 4, 207–214 [PubMed]
41. Elias J. E., Gygi S. P. (2010) Target-decoy search strategy for mass spectrometry-based proteomics. Methods Mol. Biol. 604, 55–71 [PMC free article] [PubMed]
42. Medina-Aunon J. A., Paradela A., Macht M., Thiele H., Corthals G., Albar J. P. (2010) Protein information and knowledge extractor: discovering biological information from proteomics data. Proteomics 10, 3262–3271 [PubMed]
43. Diez-Rivero C. M., Chenlo B., Zuluaga P., Reche P. A. (2010) Quantitative modeling of peptide binding to TAP using support vector machine. Proteins 78, 63–72 [PubMed]
44. Hearn A., York I. A., Rock K. L. (2009) The specificity of trimming of MHC class I-presented peptides in the endoplasmic reticulum. J. Immunol. 183, 5526–5536 [PMC free article] [PubMed]
45. Hassan C., Kester M. G., de Ru A. H., Hombrink P., Drijfhout J. W., Nijveen H., Leunissen J. A., Heemskerk M. H., Falkenburg J. H., van Veelen P. A. (2013) The human leukocyte antigen-presented ligandome of B lymphocytes. Mol. Cell Proteomics 12, 1829–1843 [PMC free article] [PubMed]
46. Alvarez I., Martí M., Vazquez J., Camafeita E., Ogueta S., Lopez de Castro J. A. (2001). The Cys-67 residue of HLA-B27 influences cell surface stability, peptide specificity, and T-cell antigen presentation. J. Biol. Chem. 276, 48740–48747 [PubMed]
47. Ben Dror L., Barnea E., Beer I., Mann M., Admon A. (2010) The HLA-B*2705 peptidome. Arthritis Rheum. 62, 420–429 [PubMed]
48. Infantes S., Lorente E., Barnea E., Beer I., Barriga A., Lasala F., Jimenez M., Admon A., Lopez D. (2013) Natural HLA-B*2705 protein ligands with glutamine as anchor motif: implications for HLA-B27 association with spondyloarthropathy. J. Biol. Chem. 288, 10882–10889 [PMC free article] [PubMed]
49. Sesma L., Montserrat V., Lamas J. R., Marina A., Vazquez J., Lopez de Castro J. A. (2002). The peptide repertoires of HLA-B27 subtypes differentially associated to spondyloarthropathy (B*2704 and B*2706) differ by specific changes at three anchor positions. J. Biol. Chem. 277, 16744–16749 [PubMed]
50. Marti M., Alvarez I., Montserrat V., Lopez de Castro J. A. (2001). Large sharing of T-cell epitopes and natural ligands between HLA-B27 subtypes (B*2702 and B*2705) associated with spondyloarthritis. Tissue Antigens 58, 351–362 [PubMed]
51. Van Endert P. M., Riganelli D., Greco G., Fleischhauer K., Sidney J., Sette A., Bach J. F. (1995) The peptide-binding motif for the human transporter associated with antigen processing. J. Exp. Med. 182, 1883–1895 [PMC free article] [PubMed]
52. WTCC Consortium (2007) Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants. Nat. Genet. 39, 1329–1337 [PMC free article] [PubMed]
53. Haroon N., Tsui F. W., Uchanska-Ziegler B., Ziegler A., Inman R. D. (2012) Endoplasmic reticulum aminopeptidase 1 (ERAP1) exhibits functionally significant interaction with HLA-B27 and relates to subtype specificity in ankylosing spondylitis. Ann. Rheum. Dis. 71, 589–595 [PubMed]
54. Villadangos J. A., Galocha B., Garcia F., Albar J. P., Lopez de Castro J. A. (1995). Modulation of peptide binding by HLA-B27 polymorphism in pockets A and B, and peptide specificity of B*2703. Eur. J. Immunol. 25, 2370–2377 [PubMed]
55. Montserrat V., Galocha B., Marcilla M., Vazquez M., Lopez de Castro J. A. (2006). HLA-B*2704, an allotype associated with ankylosing spondylitis, is critically dependent on transporter associated with antigen processing and relatively independent of tapasin and immunoproteasome for maturation, surface expression, and T cell recognition: relationship to B*2705 and B*2706. J. Immunol. 177, 7015–7023 [PubMed]
56. Williams A. P., Peh C. A., Purcell A. W., McCluskey J., Elliott T. S. (2002) Optimization of the MHC class I peptide cargo is dependent on tapasin. Immunity. 16, 509–520 [PubMed]
57. Galocha B., López de Castro J. A. (2008). Folding of HLA-B27 subtypes is determined by the global effect of polymorphic residues and shows incomplete correspondence to ankylosing spondylitis. Arthritis Rheum. 58, 401–412 [PubMed]
58. Galocha B., López de Castro J. A. (2010). Mutational analysis reveals a complex interplay of peptide binding and multiple biological features of HLA-B27. J. Biol. Chem. 285, 39180–39190 [PMC free article] [PubMed]
59. Lopez-Larrea C., Mijiyawa M., Gonzalez S., Fernandez-Morera J. L., Blanco-Gelaz M. A., Martinez-Borra J., Lopez-Vazquez A. (2002) Ankylosing spondylitis is associated with HLA-B*1403 in a West African population. Arthritis Rheum. 46, 2968–2971 [PubMed]
60. Diaz-Pena R., Blanco-Gelaz M. A., Njobvu P., Lopez-Vazquez A., Suarez-Alvarez B., Lopez-Larrea C. (2008) Influence of HLA-B*5703 and HLA-B*1403 on susceptibility to spondyloarthropathies in the Zambian population. J. Rheumatol. 35, 2236–2240 [PubMed]
61. Merino E., Galocha B., Vazquez M. N., Lopez de Castro J. A. (2008). Disparate folding and stability of the ankylosing spondylitis-associated HLA-B*1403 and B*2705 proteins. Arthritis Rheum. 58, 3693–3704 [PubMed]
62. Merino E., Montserrat V., Paradela A., Lopez de Castro J. A. (2005). Two HLA-B14 subtypes (B*1402 and B*1403) differentially associated with ankylosing spondylitis differ substantially in peptide specificity, but have limited peptide and T-cell epitope sharing with HLA-B27. J. Biol. Chem. 280, 35868–35880 [PubMed]
63. Daniel S., Brusic V., Caillat-Zucman S., Petrovsky N., Harrison L., Riganelli D., Sinigaglia F., Gallazzi F., Hammer J., Van Endert P. M. (1998) Relationship between peptide selectivities of human transporters associated with antigen processing and HLA class I molecules. J. Immunol. 161, 617–624 [PubMed]
64. Herberts C. A., Neijssen J. J., de Haan J., Janssen L., Drijfhout J. W., Reits E. A., Neefjes J. J. (2006) Cutting edge: HLA-B27 acquires many N-terminal dibasic peptides: coupling cytosolic peptide stability to antigen presentation. J. Immunol. 176, 2697–2701 [PubMed]
65. Cortes A., Hadler J., Pointon J. P., Robinson P. C., Karaderi T., Leo P., Cremin K., Pryce K., Harris J., Lee S., Joo K. B., Shim S. C., Weisman M., Ward M., Zhou X., Garchon H. J., Chiocchia G., Nossent J., Lie B. A., Forre O., Tuomilehto J., Laiho K., Jiang L., Liu Y., Wu X., Bradbury L. A., Elewaut D., Burgos-Vargas R., Stebbings S., Appleton L., Farrah C., Lau J., Kenna T. J., Haroon N., Ferreira M. A., Yang J., Mulero J., Fernandez-Sueiro J. L., Gonzalez-Gay M. A., Lopez-Larrea C., Deloukas P., Donnelly P., Bowness P., Gafney K., Gaston H., Gladman D. D., Rahman P., Maksymowych W. P., Xu H., Crusius J. B., van der Horst-Bruinsma I. E., Chou C. T., Valle-Onate R., Romero-Sanchez C., Hansen I. M., Pimentel-Santos F.M., Inman R. D., Videm V., Martin J., Breban M., Reveille J. D., Evans D. M., Kim T. H., Wordsworth B. P., Brown M. A. (2013) Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci. Nat. Genet. 45, 730–738 [PMC free article] [PubMed]
66. Chen L., Fischer R., Peng Y., Reeves E., McHugh K., Ternette N., Hanke T., Dong T., Elliott T., Shastri N., Kollnberger S., James E., Kessler B., Bowness P. (2014) Critical role of endoplasmic reticulum aminopeptidase 1 in determining the length and sequence of peptides bound and presented by HLA-B27. Arthritis Rheumatol. 66, 284–294 [PubMed]
67. Mear J. P., Schreiber K. L., Münz C., Zhu X., Stevanovic S., Rammensee H. G., Rowland-Jones S. L., Colbert R. A. (1999) Misfolding of HLA-B27 as a result of its B pocket suggests a novel mechanism for its role in susceptibility to spondyloarthropathies. J. Immunol. 163, 6665–6670 [PubMed]
68. Dangoria N. S., DeLay M. L., Kingsbury D. J., Mear J. P., Uchanska-Ziegler B., Ziegler A., Colbert R. A. (2002) HLA-B27 misfolding is associated with aberrant intermolecular disulfide bond formation (dimerization) in the endoplasmic reticulum. J. Biol. Chem. 277, 23459–23468 [PubMed]
69. Bowness P. (2002) HLA B27 in health and disease: a double-edged sword? Rheumatology 41, 857–868 [PubMed]
70. Cauli A., Shaw J., Giles J., Hatano J., Rysnik O., Payeli S., McHugh K., Dessole G., Porru G., Desogus E., Fiedler S., Holper S., Carette A., Blanco-Gelaz M. A., Vacca A., Piga M., Ibba V., Garau P., La Nasa G., Lopez-Larrea C., Mathieu A., Renner C., Bowness P., Kollnberger S. (2013) The arthritis-associated HLA-B*27:05 allele forms more cell surface B27 dimer and free heavy chain ligands for KIR3DL2 than HLA-B*27:09. Rheumatology 52, 1952–1962 [PMC free article] [PubMed]

Articles from Molecular & Cellular Proteomics : MCP are provided here courtesy of American Society for Biochemistry and Molecular Biology