| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Endocrinology Division (H.I., S.Q., L.J.D.G.), Department of Medicine, and TB/HIV Research Lab (A.S.D.G.), Brown University, Providence, Rhode Island 02903; and EpiVax, Inc. (W.M., A.S.D.G.), Providence, Rhode Island 02903
Address all correspondence and requests for reprints to: Leslie J. De Groot, M.D., Brown University/Medicine/Endocrinology, Box G, Room E-308, 70 Ship Street, Providence, Rhode Island 02903. E-mail: leslie_degroot{at}brown.edu.
| Abstract |
|---|
|
|
|---|
Objective: We sought to identify T cell epitopes in TSHR that initiate this immune response and their interaction with human histocompatibility leukocyte antigen (HLA) molecules predisposing to GD.
Design: We examined the affinity of 31 overlapping peptides spanning the TSHR extracellular domain for binding in vitro to five purified HLA-DR molecules; DRB1*0101 (DR1), DRB1*1501 (DR2), DRB1*0301 (DR3), DRB1*1101 (DR5), and DRB1*0701 (DR7). We scanned the TSHR extracellular domain using a T cell epitope-mapping algorithm, EpiMatrix. We compared these results with clinical studies of GD patients measuring in vitro T cell responses to the peptides.
Setting: The study was conducted at a university laboratory.
Patients: Patients included 200 serial adult clinic patients with GD.
Intervention: There were no interventions.
Main Outcome Measurements: Binding affinity of epitopes, predicted affinity, and reported T cell stimulation data were measured.
Results: Most peptides bound with intermediate or high affinity to one or more HLA-DR molecule. Peptides binding to HLA-DR3 and HLA-DR5, which predispose to GD, exhibited moderate binding affinities overall, whereas most peptides binding to GD-protective HLA-DR7 bound with high affinity. These differences may relate to T cell selection in the thymus. Binding affinity of peptides correlated strongly with EpiMatrix-predicted affinity for HLA-DRB1*0101, DRB1*1501, DR3, and DRB1*0701 but not HLA-DR5. Average IC50 values correlated significantly with clinical T cell stimulation data.
Conclusions: Three different methods for identifying immunogenic peptides did not provide a uniform picture of important TSHR epitopes. However, peptide 132150 (GIFNTGLKMFPDLTKVYST) was identified by three methods as an important epitope in GD; the possible importance of peptides 145163, 158176, 207222, 248263, 272291, and 343362 was also identified.
| Introduction |
|---|
|
|
|---|
Patients with GD present with hyperthyroidism due to uncontrolled stimulation of the thyroid by antibodies that bind to the TSH receptor (TSHR) (1). CD4+ T cells recognizing TSHR peptides presented by human histocompatibility leukocyte antigen (HLA) class II molecules play a central role in the development of this autoimmunity. Attention has focused on the extracellular domain (ECD) of the TSHR because it is exposed on the surface of the cell and because a large portion of the ECD appears to be shed into surrounding tissue fluid during receptor processing (2). The ECD of the TSHR is composed of 418 amino acids linked to the transmembrane domain of 346 residues. The transmembrane domain contains seven transmembrane segments. The location of CD4+ T cell epitopes within the TSHR-ECD has been the subject of a number of studies.
Genetic factors clearly play an important role in GD. An association between certain HLA class II alleles and GD has been documented. Recently the CTLA-4 gene region has also been identified as a second locus conferring susceptibility to GD (3). This gene encodes an important negative regulator of the immune system.
Certain HLA class II genes such as HLA-DR3 and HLA-DQA1*0501 predispose to GD among Caucasian patients (4, 5). In contrast, HLA-DR7 or HLA-DQA1*0201 may provide protection against GD (6, 7). The role of HLA class II molecules is to selectively bind and present peptides, derived from processed protein captured by endocytosis, in the HLA class II epitope presenting cleft (8). Immunogenic proteins such as TSHR-ECD are endocytosed by antigen-presenting cell (APC), processed in an endosomal pathway in which they are trimmed to 1020 amino acid peptides, displace the common CLIP peptide from its association with HLA class II molecules, and are then transported to the cell surface in which they are displayed to T cell receptors.
As has been documented for other autoimmune diseases, it is not unusual to find multiple reactive epitopes within the same autoantigen. This has been observed in the case of glutamic acid decarboxylase peptides in type 1 diabetes mellitus with HLA-DR4 (9) and for type IV collagen peptides in Goodpastures disease with HLA-DR15 (10) and acetylcholine receptor peptides in myasthenia gravis with HLA-DR3 (11). Many studies have partially defined the epitopes important in development of GD. We have previously reported sequences containing amino acid residues 145163, 158176, 207222, 248263, 272291, and 343362 considered to be important epitopes recognized by T cells from patients with GD (12, 13, 14, 15, 16, 17). Martin et al. (18) found TSHR peptides 5271, 142161, 202221, and 247266 to be frequently recognized by CD4+ T cells from patients with GD. Tandon et al. (19) found that TSHR 146165, 160179, and 202221 were also possibly relevant.
One explanation of the association between HLA and GD is that high-risk HLA alleles present autoantigenic epitopes more efficiently than low-risk HLA alleles (20). Alternatively, both protective and predisposing alleles may be capable of presenting TSHR-derived epitopes but differ with respect to the functional characteristics of the CD4+ T cells that respond to the epitopes. Competition between high- and low-risk alleles for binding of autoantigenic fragments derived from TSHR could also affect the development of GD. Due to a higher affinity for specific fragments, protective alleles might prevent binding and presentation of crucial epitopes by high-risk alleles. High-affinity binding during thymic T cell selection might lead to deletion of cells reactive to specific epitopes. On the other hand, CD4+CD8CD25+ regulatory T cells generated in thymus may play a protective role in GD. And finally, some HLA alleles may be less likely to present key epitopes that promote TSHR antibodies.
To elucidate the mechanism underlying the role of HLA class II molecules, we purified HLA-DR molecules from Epstein-Barr virus-transformed human B-lymphoblastoid cell lines (BLCL) and evaluated peptide binding to the HLA-DR molecules in vitro. We also scanned the TSHR-ECD using the EpiMatix T cell epitope-mapping algorithm. Then we compared these two methods with clinical T cell stimulation data.
| Materials and Methods |
|---|
|
|
|---|
L243 (anti-HLA-DR monoclonal antibody) hybridoma cells, HB-55, were purchased from American Type Culture Collection (Manassas, VA) and expanded in serum-free medium (Life Technologies, Inc., Grand Island, NY). This antibody is reported to bind any DR and produce minimal contamination with DQ protein (21). Antibody secreted by these cells was purified by using protein A-Sepharose chromatography. An L243-Sepharose column was prepared by coupling 10 mg of purified L243 with 2 ml of cyanogen bromide-activated Sepharose 4B (Pharmacia, Uppsala, Sweden).
Cell lines
The following homozygous Epstein-Barr virus-transformed BLCLs were kindly provided by Gerard Nepom and Susan Masewicz (Virginia Mason Research Center, Seattle, WA): LG2 (HLA-DRB1*0101, DR1), MGAR (HLA-DRB1*1501, DR2), QBL (HLA-DRB1*0301, DR3), Sweig (HLA-DRB1*1101, DR5), and PLH (HLA-DRB1*0701, DR7). These transformed cells were expanded in RPMI 1640 with 2 mM L-glutamine, 10 mM HEPES buffer, 50 µg/ml streptomycin, 50 U/ml penicillin, and 10% heat-inactivated fetal bovine serum (Life Technologies).
Purification of HLA-DR molecules
HLA-DR molecules were purified from 108 BLCLs expressing a known DR protein in the homozygous state, as previously described (12). All of the cell manipulations were performed at 4 C. Cells were lysed in 20 ml of lysis buffer [10 mM Tris, 150 mM NaCl, 0.5% Triton X-100, 1 mM EDTA, 1 mM Pefabloc (pH 7.4)]. The lysates were cleared by centrifugation at 100,000 x g for 60 min at 4 C. Lysates were passed through a protein A column and then through an L243 (anti-HLA-DR)-Sepharose affinity column. The column was washed extensively with PBS and 0.1% Triton X-100 (pH 7.4) and eluted with 10 ml of 50 mM glycine-NaOH, 0.1% Triton X-100 (pH 11.5). Eluates were immediately neutralized with 2 M glycine-HCl (pH 3.0), dialyzed against 50-fold volume of 10 mM Tris, 137 mM NaCl, 0.1% Triton X-100 (pH8), and concentrated with centricon (10 kDa; Amicon, Beverly, MA). The HLA-DR molecules were then evaluated on 10% SDS-PAGE for purity. The concentration of HLA-DR molecules was determined using BCA protein concentration assay kit (Pierce, Rockford, IL).
Peptides
Thirty-one individual TSHR-ECD peptides 1621 amino acids in length were prepared by solid-phase peptide synthesis, each one overlapping the previous sequence by five to six amino acids (22). Peptides were purified by reverse-phase HPLC.
Mycobacterium tuberculosis 65-kDa heat shock protein peptide 313 (hsp313) (restricted by HLA-DR3) and influenza hemagglutinin peptide 307319 (HA307319) (restricted by HLA-DR1, HLA-DR2, HLA-DR5, and HLA-DR7) were used as controls, respectively (23, 24).
Both the hsp313 and HA307319 peptides were biotinylated and used as competitor peptides.
All synthetic peptides were checked for purity by amino acid analysis and mass spectrometry. Purity was greater than 90% by these methods.
Binding and competition assays of HLA-DR and TSHR-ECD peptides
Serial dilutions (0.01100 µM) of nonbiotinylated TSHR-ECD peptides and hsp313 (or HA307319) were incubated in a 96-well plate with a given purified HLA-DR molecule (100 ng) for 45 min, followed by the addition of a single concentration (0.2 µM) of biotinylated hsp313 or HA307317 peptide in 150 mM citrate phosphate buffer containing 0.75% n-octyl-ß-D-glucopyranoside and 1 mM Pefabloc (pH 5.4). The plate was incubated at 37 C overnight. Each well of a 96-well plate (Coaster, Cambridge, MA) was coated with 100 µl of anti-HLA-DR capture antibody (L243, 10 µg/ml) in 0.2M borate buffer (pH 8.1) and incubated overnight at 4 C. The antibody plate was washed five times with PBS, containing 0.05% Tween 20, blocked with PBS blocking solution [including 5% fetal bovine serum/fetal calf serum (1:1)] at room temperature for 3 h, and washed. Tris buffer (50 µl) [50 mM (pH 8.0)] containing 0.75% n-octyl-ß-D-glucopyranoside was added to each well. Bound major histocompatibility complex (MHC)-peptide complexes were transferred from the first incubation plate to this plate, which was incubated at 4 C, overnight. After washing, europium-labeled streptavidin (Wallac, Gaithersburg, MD) was added to each well and incubated at room temperature for 60 min, followed by enhancement buffer (Wallac) for 20 min at room temperature. Fluorescence was measured with a Delfia 1232 fluorometer (Wallac). Each assay of each peptide was done in triplicate, and at least three complete studies were done of binding of all peptides to each HLA-DR molecule. The results of the three separate assays were in close agreement and are presented as an average of the results. IC50 of greater than 100 were not determined and are reported as greater than 100.
Sum of T cells score
To determine T cell proliferation in response to TSHR peptide, peripheral blood mononuclear cells (PBMCs) (3 x 105/well) were cultured with or without antigen and pulsed with [3H]thymidine (1 µCi/well; ICN Radiochemicals, Irvine, CA) for 16 h as previously described (13, 15). [3H]thymidine incorporation was measured by liquid scintillation counting. Results were expressed as a stimulation index (SI), the ratio of [3H]thymidine uptake in the presence of antigen to [3H]thymidine uptake in culture medium alone.
A summary of clinically derived T cell stimulation data were defined as follows and is reported below. For each of our previous studies of T cell stimulation by TSHR epitopes, we scored statistically significant positive responses for each epitope as 1 and amalgamated the results. Epitopes positive in stimulation of T cell lines or T cell clones were also given a score of 1. We included in the same manner the data by Martin et al. (18) and Tandon et al. (19) for epitopes recognized in their studies. These scores thus represented an arbitrary but useful compilation of available reported studies evaluating T call reactivity to TSHR peptides. HLA typing data were available on only one small subset of 13 patients in our studies
EpiMatrix Z-score
EpiMatrix is a T cell epitope-mapping algorithm that is used to identify putative HLA ligands/T cell epitopes contained within protein sequences. Computation is performed by comparing peptides sequences with a set of HLA allele-specific coefficient matrices. Each matrix contains a set of 180 coefficients, one for each of nine positions, or pockets, contained in the floor of the MHC binding groove and for each of the 20 amino acids. To complete an analysis, target protein sequences are parsed into overlapping 9-mer frames in which each frame overlaps the last by eight amino acids. For any given frame, each amino acid is assigned a coefficient based on its type and position within the HLA binding groove. Coefficients are summed to produce a raw score. Raw scores are then normalized with respect to a distribution derived from a large set of randomly generated peptide sequences. The resulting Z-scores from this distribution are directly comparable across prediction for different alleles. For this study the sequences of each of the synthesized peptides were screened against the same alleles used in the binding assays. For each peptide we selected the highest scoring 9-mer peptide, within the synthetic peptide sequence, to represent the score for the sequence. EpiMatrix has been used in a large number of studies (25, 26) to predict T cell epitopes. The algorithm has been shown to be an accurate predictor of class I and class II epitopes (25, 26)
Statistical analysis
Pearson product moment correlation coefficients (r values) were calculated for IC50 values vs. Z-scores, IC50 values vs. sum of T cells, and Z-scores vs. sum of T cells. Pearson correlation coefficients were also calculated for IC50 values vs. SI of DR3+ patients and Z-scores vs. SI of DR3+ patients. Probabilities for the correlations of less than 0.05 were considered statistically significant.
For comparison between in vitro peptide binding assays and EpiMatrix predictions, IC50 cut-off values for positive and negative binders were set to 50 µM, and Z-score cut-off values were set to 1.44 to compare sensitivity and specificity.
| Results |
|---|
|
|
|---|
|
An EpiMatrix analysis of predicted epitopes within the TSHR-ECD amino acid sequence was performed for each HLA-DR allele using progressive 9-mer sequences (overlapping by eight amino acids). The 9-mer core epitope sequence, within the specific 1420 amino acid synthetic peptide, predicted by EpiMatrix to have the highest Z-score for a DR molecule, is reported in the columns labeled "P" in Table 2
. Each DR molecule was separately analyzed. The predicted binding affinity for the 9-mer sequence is given as a Z-score in the columns in Table 2
headed by the specific DR allele. Peptides that receive Z-scores greater than 1.64 are generally considered to be in the top 5% of all binding scores, and Z-scores greater than 1.28 are in the top 10% of all binding scores. Peptides receiving Z-scores above 1.64 are considered to be moderately likely to bind, whereas peptides scoring 1.281.64 are considered only weakly likely to bind. Lower scoring peptides generally do not bind. The 9-mers with the highest Z-score within each of the sequences labeled 131 in Table 1
were chosen to represent the predicted binding affinity of each HLA-DR allele for each of the 31 peptides in Table 1
.
|
|
|
|
| Discussion |
|---|
|
|
|---|
Different interpretations of the relation of binding affinity to importance in disease process have been offered. A common theme is that high-affinity binding to HLA-DR correlates with both efficiency in antigen presentation and T cell deletion in the thymic T cell selection process (9, 10, 11). Alternatively, some reports indicate that peptides with high-affinity binding to HLA-DR molecules lead to thymic deletion of the cognate T cells and that peptides that exhibit only moderate binding affinities are more apt to enter in the circulation and participate in autoimmune disease (28, 31, 32).
It must also be noted that the interaction of epitope, APC, and T cell involves numerous other factors, including variations in the specific type of APC involved and the nature of the responding T cell. As a result, development of tolerance occurs at both central and peripheral locations (33, 34). Because of immunological tolerance, severe autoimmune reactions are generally suppressed. But surprisingly, 1030% of females and 515% of males develop autoantibodies against thyroid-related antigens, such as TSHR, thyroglobulin, or thyroid peroxidase (35, 36, 37, 38).
The current studies do not evaluate the affinity of binding of the T cell receptor complex to the epitope-HLA complex, which must in fact be a crucial interaction. We do have data on the T cell responses to the epitopes, presented herein as sum of T cell scores. It is known that amino acids in the third, fifth, and eighth positions of the core nonamer tend to point up, out of the binding cleft, and make contact with the T cell receptor (39). This aspect of epitope function can be more thoroughly analyzed when epitope-specific T cell lines or clones are available. T cell receptor contact modeling algorithms are not currently available. B cell epitope mapping algorithms have been described, but none have been shown to be accurate predictors of B cell immunogenicity. Finally, whereas some T cell epitopes overlap B cell epitopes (40), in most cases there is no good correlation between B cell epitopes and T cell epitopes. Also, and perhaps more importantly in this study, the relatively arbitrary start and stop points chosen for synthetic peptides may, or may not, interrupt the linear sequence of a naturally processed epitope.
In this study, we determined binding affinities of a set of overlapping 1621 mer synthetic peptides covering the TSHR-ECD to five kinds of purified HLA-DR molecules. This approach is valuable because it defines affinity of peptide binding to a specific HLA allele. However, it does not exactly mimic the complexities of antigen presentation in patients, who would generally be heterozygous for DR genes and would also carry functional DQ alleles. We also determined the probable epitopes by using a T cell epitope-mapping algorithm. Nearly all peptides bound to at least one of the HLA-DR molecules, which is expected and in fact represents the presumed function of having different HLA-DR molecules in the population by providing protection against pathogens to at least some subset of individuals. One interesting observation that we made was that many of the peptides (11 of 31) that we selected for this study bound with high affinity to HLA-DR7, which has been associated with protection from GD. Only two of 31 bound with high affinity to DR3. One explanation for this observation may be that high-affinity binding contributed to the removal of disease-causing T cells in the process of thymic deletion in subjects bearing DR. There is good evidence that TSHR is present in the thymus (41). In addition, Schmidt et al. (42) have shown that part of the mechanism for MHC-linked resistance to autoimmunity is negative selection of pathogenic autoreactive T cells, implying high affinity of some immunodominant peptides for disease-protective MHC molecules.
In contrast, only two peptides (109-124, 132-150) bound with high affinity to HLA-DR3, an allele that has been associated with thyroid disease (Table 1
). One of the two (peptide 132150) also binds strongly to GD-related HLA-DR5. These two peptides are not among those showing high affinity for HLA-DR7. Four peptides also displayed high affinity binding to HLA-DR5, another GD-associated allele. The implication of this observation is not entirely clear because high-affinity binding has been associated with thymic deletion. Eighteen of 31 peptides bound with moderate affinity to DR3. The lack of many high-affinity binders and the predominance of moderate affinity binders for this allele might lead to retention of autoreactive TSHR-ECD-specific T cells in HLA-DR3+ patients
In a previous study, we examined binding of 13 TSHR-ECD peptides to HLA-DR3 using a different assay (12). The peptides in that study were classified as immunodominant or not. Peptides classified as nonimmunodominant had similar IC50 values in the two studies. Almost all nonimmunodominant had IC50 values greater than 40 in both studies. The peptides designated as immunodominant in the prior study had IC50 values of approximately 220, whereas in the present study, all but one had IC50 values of 2040. The reasons for this small shift are uncertain but may relate to the different techniques used. Peptide 109124 demonstrated the same very high affinity in both studies.
We previously performed several T cell immunogenicity studies using blood samples from patients with GD and reported several sequences that appeared to be immunodominant epitopes (peptides 145163, 158176, 207222, 248263, 272291, and 343362). In this study each of these peptides bound with moderate affinity (IC50 1450) to DR3.
We evaluated each of the peptides by comparing the results of our binding affinity studies with those of a T cell epitope-mapping algorithm and found a significant correlation between the predictions and HLA-DR binding for four of the five HLA-DR molecules evaluated (Table 3A
). There was no correlation between HLA-DR5 predictions and in vitro binding studies. These observations support the use of T cell epitope mapping using the EpiMatrix algorithm for the prediction of HLA ligands, except for HLA-DR5. The reason for the discrepancy between HLA-DR5 prediction and binding is under evaluation. The peptides that were synthesized often truncated the predicted epitope, which could have led to reductions in the correlation between predictions and immunogenicity. Furthermore, flanking amino acids act to stabilize the MHC molecule outside the binding groove. Indeed, amino terminal truncation of the flanking residues (Table 2
) seemed to reduce the likelihood that peptides predicted to bind would bind in vitro, whereas COOH terminus truncation did not affect the accuracy of the binding predictions. This is due to the relative importance of the amino terminal amino acids in MHC/peptide interactions. Shorter amino terminal flanks could reduce the ability of the important pocket-binding side chain from the key amino acid to interact with the P1 binding pocket in the HLA binding groove.
We also compared the IC50 values and T cell epitope predictions with results from T cell stimulation assays using GD patients blood samples and found a poor correlation with the T cell stimulation assays. Although this is disappointing, perhaps it should be anticipated. First, the sequences we chose may not correspond to those produced in vivo in APCs. Also, eptitope binding assays and predictions relate to specific DR structures, but there is effectively no information about the HLA alleles in the patient populations studied. Furthermore, whereas binding of epitope to HLA is required for T cell response, it is not sufficient. The T cell receptor must recognize the epitope in the context of the HLA molecule, and amino acids in positions 3, 5, and 8 in the core sequence are important in this regard. Epitopes can bind effectively but not stimulate a T cell, depending on single amino acid differences in the structure. However, the correlation between IC50 values and T cell immunogenicity was significant when values for all HLA-DRs were averaged, in a manner that might be more representative of a population with diverse DRs.
Evaluating these TSHR epitopes by three distinct methods does not provide a uniform picture of their importance in development of GD. The results are summarized in Table 6
by criteria of high in vitro binding affinity to DR molecules, high EpiMatrix Z-score, and clinical T cell stimulation data. Only peptide 132150 satisfies all criteria, which confirms its possible immunogenic importance. Peptide 158176, which is most active in clinical assays, fails by other criteria. Peptides 11 (145163), 16 (207222), 20 (248263), 22 (272291), and 27 (343362) are identified again as possibly important epitopes, based on observed activity in clinical studies and moderate binding affinity.
|
| Acknowledgments |
|---|
| Footnotes |
|---|
First Published Online April 4, 2006
Abbreviations: APC, Antigen-presenting cell; BLCL, B-lymphoblastoid cell line; ECD, extracellular domain; GD, Graves disease; HLA, histocompatibility leukocyte antigen; hsp313, heat shock protein peptide 313; MHC, major histocompatibility complex; PBMC, peripheral blood mononuclear cell; SI, stimulation index; TSHR, TSH receptor.
Received November 22, 2005.
Accepted March 29, 2006.
| References |
|---|
|
|
|---|
Read all eLetters
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |