Lymphocyte proliferation induced by high-affinity peptides for HLA-B*51:01 in Behçet’s uveitis
Toshikatsu Kaburaki aff001; Hisae Nakahara aff001; Rie Tanaka aff001; Kimiko Okinaga aff001; Hidetoshi Kawashima aff002; Youichiro Hamasaki aff003; Thanyada Rungrotmongkol aff004; Supot Hannongbua aff005; Hiroshi Noguchi aff006; Makoto Aihara aff001; Fujio Takeuchi aff007
Authors place of work:
Department of Ophthalmology, The University of Tokyo Graduate School of Medicine, Tokyo, Japan
aff001; Department of Ophthalmology, Jichi Medical University, Tochigi, Japan
aff002; Department of Dermatology, Dokkyo Medical University, Tochigi, Japan
aff003; Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
aff004; Computational Chemistry Unit Cell, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
aff005; Department of Pharmacognosy, Nihon Pharmaceutical University, Saitama, Japan
aff006; Department of Pharmacology, University of Shizuoka, Shizuoka, Japan
aff007; Department of Health and Nutrition, Tokyo Seiei University, Tokyo, Japan
Published in the journal:
PLoS ONE 14(9)
Several proteins have been proposed as candidate auto-antigens in the pathogenesis of Behçet’s disease (BD). In this study, we aimed to confirm the cellular responses to candidate peptide autoantigens with high affinity for the HLA-B*51:01 molecule using computerized binding predictions and molecular dynamics simulations. We identified two new candidate peptides (HSP65PD, derived from heat shock protein-65, and B51PD, derived from HLA-B*51:01) with high-affinity to the HLA-B*51:01 binding pocket using the Immune Epitope Database for Major Histocompatibility Complex-I Binding Prediction and molecular dynamics simulations. The peptide-induced proliferation of lymphocytes from patients with BD, sarcoidosis, Vogt–Koyanagi–Harada disease (VKH) with panuveitis, systemic scleroderma (SSc) without uveitis, and healthy controls (HC) was investigated using the bromodeoxyuridine assay. The proliferative response of leukocytes to HSP65PD was significantly higher in BD (SI 1.92 ± 0.65) than that in sarcoidosis (SI 1.38 ± 0.46), VKH (SI 1.40 ± 0.33), SSc (SI 1.32 ± 0.31), and HC (SI 1.27 ± 0.28) (P = 0.0004, P = 0.0007, P < 0.0001, P < 0.0001, respectively, Mann-Whitney’s U-test). The proliferative response of leukocytes to B51PD was also higher in BD than that in sarcoidosis, VKH, SSc without uveitis, and HC, whereas no significant differences were observed among the five groups in response to a control peptide derived from topoisomerase 1. A significantly higher response to HPS65PD and B51PD was observed in the HLA-B*51:01-positive patients with BD than in the HLA-B*51:01-negative patients. In conclusion, two peptides that had high affinity to HLA-B*51:01 in computerized binding prediction showed significantly higher response in HLA-B*51:01-positive patients with BD, indicating the usefulness of computerized simulations for identifying autoreactive peptides to HLAs.
Medicine and health sciences – Inflammatory diseases – Sarcoidosis – Rheumatology – Ophthalmology – Uveitis – Immunology – Lymphocyte proliferation – Pathology and laboratory medicine – Pathogenesis – Biology and life sciences – Cell biology – Cellular types – Animal cells – Blood cells – White blood cells – T cells – Lymphocytes – Immune cells – Anatomy – Body fluids – Blood – Physiology
Behçet’s disease (BD), a multi-organ inflammatory disease characterized by recurrent uveitis, oral aphthous stomatitis, genital ulcers, and skin lesions, is more prevalent in countries along the ancient “Silk Road” from Japan to the Mediterranean region than elsewhere in the world. Although the pathogenesis of BD is not well understood, genetic factors, especially human leukocyte antigens (HLA), are believed to play an important role in the pathogenesis of this disease [1, 2]. T-cell mediated immune reactions may also be associated with this disease as patients with active BD had significantly more IFNγ-producing CD4+ cells  and activated CD8+ and γδ T cells  than those in inactive cases and control subjects. Moreover, the frequencies of Th1/Th17 cells in peripheral blood were increased in BD with active uveitis  and the alterations of T cell balance, especially increase of T helper 1 (Th1) and Th17 cells and decrease of regulatory T (Treg) cells, are considered to have a significant role in the pathogenesis of BD .
T cells recognize a complex composed of a major histocompatibility complex (MHC) molecule and an antigenic peptide or epitope. The identification of pathogenic antigens or antigenic epitopes is crucial for facilitating the study of immune-mediated inflammatory diseases. Several studies have demonstrated that peripheral blood lymphocytes (PBL) from patients with BD and panuveitis are sensitized to retinal S-antigen (S-Ag), and autoreactive lymphocytes to S-Ag or its peptide motif were present in the peripheral blood of patients with BD [7,8]. Furthermore, several candidate antigens or their peptide motifs , such as major histocompatibility complex class I chain-related gene A (MICA) , heat shock protein-65 (HSP65) , retinal S-antigen , and alpha-tropomyosin , have been reported to elicit response from PBL of patients with BD. In addition, peptides derived from the HLA-B*51:01 molecule itself have a predicted binding capacity to HLA-B*51:01 [9, 13]. Although the high prevalence of HLA-B51 in BD suggested the importance of this HLA type in the pathogenesis of BD, the putative role of these candidate antigens and HLA-B51 in the pathogenesis of BD is still obscure.
The accuracy of predicting peptide-MHC binding is useful for candidate T cell epitope selection as it minimizes the number of experiments required for their identification. The binding groove of class I MHC can accommodate 8–11 amino acid-long short peptides . Sequence motifs are still the best-known tools for predicting the peptide specificity of allele-dependent MHC peptide binding. Currently, the MHC-I binding for any peptide sequences is predicted using the Immune Epitope Database (IEDB) analysis resource Consensus tool , which combines predictions from ANN aka NetMHC (4.0) [16–18], SMM , and Comblib , all of which are freely available on the IEDB website (http://www.immuneepitope.org).
Recently, using molecular dynamics simulations, we demonstrated that a 9-mer peptide, MICA-TM, which had been reported to induce proliferation of PBL from patients with BD, had a significantly stronger total binding free energy with the BD-associated HLA alleles (HLA-B*51:01 and A*26:01) than with the non-associated alleles (HLA-B*35:01 and A*11:01) . This result suggests that structure, dynamics, and energetics may explain the recognition and selective binding of immunogenic peptides to specific HLAs. Thus, computerized binding predictions of peptide motifs with high affinity to HLA  or computational dynamics simulation of the bonding energy between HLA and candidate peptides  may select and provide T cell epitopes to induce T cell proliferation in BD effectively.
In the present study, we investigated the proliferation of PBL from patients with BD patients and panuveitis along with that in non-BD patients induced by newly determined candidate peptides that have high affinity to HLA-B*51:01 using the IEDB analysis resource Consensus tool .
Materials and methods
A total of 47 patients with BD and panuveitis who met the diagnostic criteria established by the Behçet’s Disease Research Committee of Japan  and were followed at the outpatient uveitis clinic of the University of Tokyo Hospital were studied. Patients with other forms of non-infectious panuveitis (n = 37; 19 cases of sarcoidosis and 17 cases of Vogt–Koyanagi–Harada disease (VKH)) and 15 patients with systemic inflammatory disease without uveitis and systemic scleroderma (SSc) acted as controls. Furthermore, 17 age- and sex-matched healthy individuals (healthy control, HC) were enrolled in this study as negative controls. The diagnosis of sarcoidosis , VKH , and SSc  was based on the respective specific diagnostic criteria. The inclusion criteria for BD and other non-infectious panuveitis in this study were as follows; (1) patients with panuveitis who had recurrence of uveitis within 6 months, or (2) patients with active systemic inflammatory signs owing to the disease (i.e., oral aphthae in BD and granulomatous skin lesions in sarcoidosis). The patients were enrolled in this study irrespective of their current medical therapy and disease severity. This study was approved by the ethical committees of the University of Tokyo Hospital (G10136) and written informed consent for participating in this study was obtained from all participants prior to blood collection.
The severity of ocular inflammation of uveitis in the inflamed eyes and use of systemic corticosteroid and immunosuppressants at blood collection were retrospectively examined from the clinical records. The severity of ocular inflammation was graded according to the criteria established by the Standardization of Uveitis Nomenclature (SUN) Working Group . Briefly, cells in the anterior chamber (maximum 4 points) were graded using a semi-quantitative scoring system of 6 grades (0, 0.5+, 1+, 2+, 3+, and 4+). Vitreous haze (maximum 4 points) was evaluated using a semi-quantitative scoring system of 6 grades (0, trace, 1+, 2+, 3+, and 4+) based on the clarity of the optic disc, retinal vessels, and nerve fiber layers in fundus examination reported by Nussenblatt et al. . For convenience, vitreous haze graded as trace in the Nussenblatt grading  was noted as 0.5+ in this study. The severity of uveitis was scored in the severe inflamed eye if both eyes were inflamed at blood collection.
Selection and synthesis of antigenic peptides
Several candidate peptide sequences with potentially high binding capacity to HLA-B*51:01 determined using bioinformatics programs  have been previously reported . In the present study, we focused on two peptides, HSP65PD and B51PD, derived from HSP65 and HLA-B*51:01, respectively, as they were speculated to have high binding capacity to HLA-B*5101  and the proliferation of lymphocytes from patients with BD against these peptide sequences had not been previously investigated. A well-known peptide for MICA (MICA-TM), which induced significant proliferation of T lymphocytes from patients with BD carrying HLA-B*51:01 was used as the positive control [4, 21]. A peptide derived from topoisomerase 1 (Topo1PD), a specific autoreactive antigen in patients with SSc, was used as the negative control [28, 29]. Recently, molecular dynamics simulation revealed that a 20-mer peptide sequence of topoisomerase 1 showed significantly high binding to HLA-DR subtypes associated with anti-topoisomerase 1 (or anti scl-70) antibody-positive SSc (HLA-DRB1*08:02, HLA-DRB1*11:01, and HLA-DRB1*11:04) . The peptide sequences in the present study and their affinity to HLA-B*51:01 calculated using the IEDB analysis resource Consensus tool  are shown in Table 1.
Separation of PBL and proliferation assay
PBL from patients and controls were isolated via centrifugation on a Ficoll-Hypaque gradient for 20 min, washed twice, and resuspended in Roswell Park Memorial Institute 1640 medium with HEPES, supplemented with 2 mM glutamine (Sigma-Aldrich Japan, Tokyo, Japan), 100 units/ml penicillin G (Sigma-Aldrich Japan, Tokyo, Japan), 100 μg/ml streptomycin sulfate, and 10% heat-inactivated human AB serum (Biowest, Nuaillé, France). Cells (1 × 105/well) were incubated in a total volume of 100 ml with the peptides or concanavalin A (4 μg/ml, Fujifilm Wako Pure Chemical Co., Osaka, Japan) in triplicates in 96-well round-bottom culture plates (Nunc, Roskilde, Denmark) at 37°C in a humidified atmosphere with 5% CO2. The cultures were incubated for 4 days and then pulsed with bromodeoxyuridine (BrdU) during the last 16 h of incubation. Cellular proliferation was assessed using a BrdU enzyme-linked immunosorbent assay kit (Cell Proliferation ELISA BrdU, Sigma), according to the manufacturer’s instructions. After co-culture with Brdu for 16 h, the BrdU labeling solution was removed and anti-BrdU antibody was added. The OD at 450 nm for BrdU was measured using a luminescence plate reader (ARVO X3, Perkin Elmer Japan, Yokohama, Japan). The results were expressed as stimulation index (SI: mean OD450nm in stimulated cultures / mean OD450 nm in unstimulated control cultures).
Lymphocyte proliferation was analyzed using the Mann–Whitney U-test. P value < 0.05 was considered statistically significant. Mean and standard deviation (SD) values are shown in the text. Statistical analysis was performed using SPSS for Windows version 23.0 (IBM Corp., Armonk, NY, USA).
The demographics of the participants of this study are shown in Table 2. The ages of the patients ranged between 21 and 81 years, and the mean age for each disease varied as the peak age of onset for each disease was different. Sex distribution was also considerably different for the same reason. Male sex was more frequent in patients with BD than in those with control uveitis (sarcoidosis and VKH) and SSc (P < 0.0001, Chi-square test). However, there were no significant differences in the age and sex distribution between BD and HC (P = 0.226, Chi-square test). The duration of the disease was longer in patients with BD, probably because patients with relatively severe disease showed recurrence of uveoretinitis despite administration of systemic immunosuppressants.
The grades of ocular inflammation in anterior chamber cells in the inflamed eyes at blood collection were 0.28 ± 0.62 in BD, 0.37 ± 0.57 in sarcoidosis, and 0.47 ± 0.86 in VKH. Furthermore, the severity scores of vitreous haze were 0.28 ± 0.62 in BD, 0.39 ± 0.58 in sarcoidosis, and 0.21 ± 0.40 in VKH (Table 2). There were no significant differences in the scores in anterior chamber cells and vitreous haze at blood collection among the three uveitis groups (P = 0.345 and P = 0.245, Kruskal-Wallis test). Systemic corticosteroid was administrated at blood sampling in 8 cases of BD (17.0%, 0.89 ± 2.19 mg), 7 cases of sarcoidosis (36.8%, 2.53 ± 4.25 mg), 10 cases of VKH (58.8%, 4.38 ± 4.25 mg), and 4 cases of SSc (26.7%, 2.93 ± 5.83 mg). Moreover, usage of TNF inhibitor (infliximab or adalimumab) was more frequent in BD (48.9%) than in sarcoidosis (10.5%), VKH (29.4%), and SSc (0%).
Responses to candidate peptides
The proliferation reactions of PBL from patients with BD after stimulation by candidate peptides (HSP65PD and B51PD) and control peptides (MICA-TM and Topo1PD) were compared to those of PBL from disease controls (sarcoidosis, VKH, SSc) and HC (Fig 1). Patients with BD had significantly higher lymphocyte proliferation responses to HSP65PD (SI 1.92 ± 0.65) than did HC (SI 1.27 ± 0.28) (P = 0.0004, Mann–Whitney’s U-test), whereas the lymphocyte responses of those with sarcoidosis (SI 1.38 ± 0.46), VKH (SI 1.40 ± 0.33), and SSc (SI 1.32 ± 0.31) to HSP65PD were not significantly higher than the response of HC (SI 1.27 ± 0.28). The responses to B51PD were also significantly stronger in patients with BD (SI 1.85 ± 0.55) than those in HC (SI 1.27 ± 0.28) (P = 0.001, Mann–Whitney’s U-test), whereas responses in patients with sarcoidosis (SI 1.40 ± 0.45), VKH (SI 1.40 ± 0.33), and SSc (SI 1.32 ± 0.31) were similar to those in HC (SI 1.27 ± 0.28) (n.s. in HC vs. sarcoidosis, VKH, and SSc, respectively, Mann–Whitney’s U-test).
As a positive control, the proliferative response to MICA-TM was assessed. The SIs for BD, sarcoidosis, VKH, SSc, and HC were 1.85 ± 0.55, 1.40 ± 0.45, 1.40 ± 0.33, 1.32 ± 0.31, and 1.27 ± 0.28, respectively (P < 0.0001, n.s., n.s., n.s. for HC vs. BD, sarcoidosis, VKH, and SSc, respectively, Mann–Whitney’s U-test). The corresponding SIs for the negative control peptide, Topo1PD, were 1.37 ± 0.31, 1.27 ± 0.34, 1.35 ± 0.34, 1.41 ± 0.26, and 1.24 ± 0.20, respectively. The peptide Topo1PD did not significantly increase the response of PBL from patients with BD, sarcoidosis, VKH, and SSc, compared to those from HC (n.s., in HC vs. BD, sarcoidosis, VKH, and SSc, respectively, Mann–Whitney’s U-test). We also examined the effect of sex on proliferation of PBL in response to the peptides in patients with BD. However, there were no significant differences according to sex on lymphocyte proliferation (S1 Fig).
Association of HLA-B*51:01 and T cell proliferation
Next, we investigated the association of HLA-B*51:01 in patients with BD and T cell proliferation induced by candidate peptides (Fig 2). The SI was significantly higher in the HLA-B*51:01-positive BD group than that in the HLA-B*51:01-negative BD group when the PBL were stimulated by HSP65, B51PD, and MICA-TM (P < 0.05, P < 0.05, P < 0.05, respectively, Mann–Whitney’s U-test). However, there was no significant difference in PBL proliferation between HLA-B*51:01-positive and -negative BD groups when the PBLs were stimulated by Topo1PD, the control peptide (n.s., Mann–Whitney’s U-test).
Association of HLA-A*26 and T cell proliferation
We assessed the association of HLA-A*26 and T cell proliferation induced by candidate peptides (Fig 3). There were no significant differences in PBL proliferation between HLA-A*26-positive and -negative BD groups when the PBL were stimulated by all four peptides (n.s., Mann–Whitney’s U-test).
Although the etiology of BD is not clear, the association of autoreactive T lymphocytes and the increase of Th1/Th17 cells and decrease of regulatory Treg cells was speculated as one of the causes of recurrent uveitis . Recently, several antigens and their peptides have been assessed as candidate autoantigens associated with BD with uveitis [4,7–12].
Recently, using molecular dynamics simulation, we demonstrated that MICA-TM, one of the candidate peptides that induced proliferation of PBL from patients with BD, had a significantly stronger total binding free energy with the BD-associated HLA alleles (HLA-B*51:01 and A*26:01) than with the non-associated alleles (HLA-B*35:01 and A*11:01) . Because molecular dynamics simulation can suggest actual molecular binding energies between peptides and HLA receptors, we speculate that this method will be a useful tool for selecting candidate peptides that can induce higher proliferative response in PBL from patients with T cell-mediated inflammatory diseases associated with specific HLA alleles. However, further analysis is required to confirm this prediction.
In the present study, we examined the proliferative response of PBL from patients with BD and other diseases to HSP65PD and B51PD. These peptide sequences had high-affinity to the HLA-B*51:01 binding pocket as observed using computerized binding predictions (IEDB analysis resource Consensus tool) . The results of this study clearly demonstrated that these two peptides can induce significantly higher proliferation of PBL from patients with BD than from HC, whereas proliferation of PBL from patients with other diseases (sarcoidosis, VKH, and SSc) was similar to that of HC (Fig 1). Notably, the proliferative reaction of PBL was HLA-B*51:01 dependent; it was significantly stronger in HLA-B*51:01-positive BD patients than in HLA-B*51:01-negative BD patients (Fig 2). In contrast, there were no significant differences in lymphocyte proliferation between HLA-A*26-positive BD patients and HLA-A*26-negative BD patients (Fig 3). This suggested that the higher affinity of the peptide to the HLA binding pocket might cause stronger lymphocyte proliferation. As this study investigated only one HLA epitope (HLA-B*51:01), further studies on other HLA epitopes and diseases are necessary to confirm this hypothesis. Furthermore, sequencing the peptides eluted from HLA-B51 positive antigen-presenting cells in patients with BD and control groups may help confirm the current results, although this may be challenging.
The control peptide Topo1PD was a 9-mer peptide from topoisomerase 1, an antigen highly associated with SSc, as anti-topoisomerase 1 antibody (anti-Scl-70) was detected in the sera of 20−70% of patients with SSc . In a previous report, a 15-mer peptide for topoisomerase 1 induced significantly strong proliferation of PBL from patients with anti-topoisomerase 1 antibody-positive SSc . However, the Topo1PD peptide, a partial 9-mer peptide of the 15-mer peptide used in the present study, did not induce a significant increase in PBL proliferation in either BD or SSc. Although the underlying reason for this observation was unclear, one of the main reasons may be that only 26.7% (4/15) cases of SSc were anti-topoisomerase 1-positive in the present study. Further studies with higher numbers of anti-topoisomerase 1-positive SSc patients are required to clarify this issue.
The reason why both HLA- B51 negative and B51 positive subjects responded in this study is unclear and may be difficult to explain. Potential reasons include: (1) there exist HLA class I types other than HLA-B51 with high binding ability to these peptides, (2) HLA class 2 (HLA-DR) might play a role in lymphocyte proliferation by these peptides, or (3) receptors other than HLA, such as Toll-like receptors, might be involved in lymphocyte proliferation.
The major strength of this study was that all experiments were performed at the University of Tokyo, which allowed maintenance of a stable experimental system and yielded reproducible results. However, this study has several limitations. First, the number of patients with sarcoidosis, VKH, SSc, and HC was low. Second, the age and sex distribution of HC and patients with sarcoidosis, VKH, and SSc were slightly different; the age and sex distribution of HC matched with those of patients with BD in this study, whereas the patients with BD and HC were relatively younger than those with sarcoidosis, VKH, and SSc. Furthermore, the female to male ratio was high among patients with SSc in Japan. Third, no inclusion criteria were set regarding current medical therapy or severity of disease, which may affect the individual SI values. Fourth, lymphocyte proliferation in response to an antigen does not prove the contribution of the antigen to pathogenesis. Despite these limitations, we believe that the present findings provide new and important insights regarding the benefits of the recent computerized binding predictions  for the selection of peptide epitopes specific for HLA molecule, which might induce higher proliferation of PBL.
In conclusion, two peptides with high affinity to HLA-B*51:01 were identified using computerized binding predictions . These peptides showed significantly high response in HLA-B*51:01-positive patients with BD. The recent computerized simulations for predicting affinity between peptides and HLA pockets might be useful for determining autoreactive peptides for HLAs.
S1 Fig [tif]
Difference in lymphocyte proliferation by sex in patients with Behçet’s uveitis.
1. Meguro A, Inoko H, Ota M, Katsuyama Y, Oka A, Okada E, et al. Genetics of Behçet disease inside and outside the MHC. Ann Rheum Dis. 2010;69: 747–754. doi: 10.1136/ard.2009.108571 19684014
2. de Menthon M, Lavalley MP, Maldini C, Guillevin L, Mahr A. HLA-B51/B5 and the risk of Behçet's disease: a systematic review and meta-analysis of case-control genetic association studies. Arthritis Rheumatol. 2009;61: 1287–1296. doi: 10.1002/art.24642 19790126
4. Yasuoka H, Yamaguchi Y, Mizuki N, Nishida T, Kawakami Y, Kuwana M. Preferential activation of circulating CD8+ and gamma delta T cells in patients with active Behçet's disease and HLA-B51. Clin Exp Rheumatol. 2008;26: S59–S63. 19026117
5. Chi W, Zhu X, Yang P, Liu X, Lin X, Zhou H, et al. Upregulated IL-23 and IL-17 in Behçet patients with active uveitis. Invest Ophthalmol Vis Sci. 2008;49: 3058–3064. doi: 10.1167/iovs.07-1390 18579762
6. Greco A, De Virgilio A, Ralli M, Ciofalo A, Mancini P, Attanasio G, et al. Behçet’s disease: New insights into pathophysiology, clinical features and treatment options. Autoimmun Rev. 2018;17: 567–575. doi: 10.1016/j.autrev.2017.12.006 29631062
7. Yamamoto JH, Minami M, Inaba G, Masuda K, Mochizuki M. Cellular autoimmunity to retinal specific antigens in patients with Behçet's disease. Br J Ophthalmol. 1993;77: 584–589. doi: 10.1136/bjo.77.9.584 8218058
8. de Smet MD, Bitar G, Mainigi S, Nussenblatt RB. Human S-antigen determinant recognition in uveitis. Invest Ophthalmol Vis Sci. 2001;42: 3233–3238. 11726628
9. Baharav E, Weinberger A. The HLA-B*5101 molecule-binding capacity to antigens used in animal models of Behçet's disease: a bioinformatics study. Isr Med Assoc J. 2012;14: 424–428. 22953618
10. Direskeneli H, Ekşioğlu-Demiralp E, Yavuz Ş, Ergun T, Shinnick T, Lehner T, et al. T cell responses to 60/65 kDa heat shock protein derived peptides in Turkish patients with Behçet’s disease. J Rheumatol. 2000;27: 708–713. 10743813
11. Kurhan-Yavuz S, Direskeneli H, Bozkurt N, Ozyazgan Y, Bavbek T, Kazokoglu H, et al. Anti-MHC autoimmunity in Behçet's disease: T cell responses to an HLA-B-derived peptide cross-reactive with retinal-S antigen in patients with uveitis. Clin Exp Immunol. 2000;120: 162–166. doi: 10.1046/j.1365-2249.2000.01176.x 10759778
12. Mahesh SP, Li Z, Buggage R, Mor F, Cohen IR, Chew EY, et al. Alpha tropomyosin as a self-antigen in patients with Behçet’s disease. Clin Exp Immunol. 2005;140: 368–375. doi: 10.1111/j.1365-2249.2005.02760.x 15807864
13. Wildner G, Thurau SR. Cross-reactivity between an HLA-B27-derived peptide and a retinal autoantigen peptide: a clue to major histocompatibility complex association with autoimmune disease. Eur J Immunol. 1994;24: 2579–2585. doi: 10.1002/eji.1830241103 7957552
14. Lundegaard C, Lund O, Nielsen M. Accurate approximation method for prediction of class I MHC affinities for peptides of length 8, 10 and 11 using prediction tools trained on 9mers. Bioinformatics. 2008;24: 1397–1398. doi: 10.1093/bioinformatics/btn128 18413329
15. Kim Y, Ponomarenko J, Zhu Z, Tamang D, Wang P, Greenbaum J, et al. Immune epitope database analysis resource. Nucleic Acids Res. 2012;4: W525–W530. doi: 10.1093/nar/gks438 22610854
16. Nielsen M, Lundegaard C, Worning P, Lauemøller SL, Lamberth K, Buus S, et al. Reliable prediction of T-cell epitopes using neural networks with novel sequence representations. Protein Sci. 2003;12: 1007–1017. doi: 10.1110/ps.0239403 12717023
17. Lundegaard C, Lamberth K, Harndahl M, Buus S, Lund O, Nielsen M. NetMHC-3.0: Accurate web accessible predictions of human, mouse, and monkey MHC class I affinities for peptides of length 8–11. Nucleic Acids Res. 2008;36: W509–W512. doi: 10.1093/nar/gkn202 18463140
18. Andreatta M, Nielsen M. Gapped sequence alignment using artificial neural networks: application to the MHC class I system. Bioinformatics. 2016; 32: 511–517. doi: 10.1093/bioinformatics/btv639 26515819
19. Peters B, Sette A. Generating quantitative models describing the sequence specificity of biological processes with the stabilized matrix method. BMC Bioinformatics. 2005;6: 132. doi: 10.1186/1471-2105-6-132 15927070
20. Sidney J, Assarsson E, Moore C, Ngo S, Pinilla C, Sette A, et al. Quantitative peptide binding motifs for 19 human and mouse MHC class I molecules derived using positional scanning combinatorial peptide libraries. Immunome Res. 2008;4: 2. doi: 10.1186/1745-7580-4-2 18221540
21. Kongkaew S, Yotmanee P, Rungrotmongkol T, Kaiyawet N, Meeprasert A, Kaburaki T, et al. Molecular dynamics simulation reveals the selective binding of human leukocyte antigen alleles associated with Behçet's disease. PLoS One. 2015;10: e0135575. doi: 10.1371/journal.pone.0135575 26331842
22. Suzuki-Kurokawa M, Suzuki N. Behçet’s disease. Clin Exp Med. 2004;4: 10–20. 15598081
23. Shijubo N, Yamaguchi T. Diagnosis criteria and classification of disease severity for sarcoidosis in Japan. Jpn J Sarcoid Granulomatous Disord. 2015;35: 3–8. (in Japanese, with English abstract). doi: 10.7878/jjsogd.35.3
24. Read RW, Holland GN, Rao NA, Tabbara KF, Ohno S, Arellanes-Garcia L, et al. Revised diagnostic criteria for Vogt-Koyanagi Harada disease: report of an international committee on nomenclature. Am J Ophthalmol. 2001;131: 647–652. doi: 10.1016/s0002-9394(01)00925-4 11336942
25. Asano Y, Jinnin M, Kawaguchi Y, Kuwana M, Goto D, Sato S, et al. Diagnostic criteria, severity classification and guidelines of systemic sclerosis. J Dermatol. 2018;45: 633–691. doi: 10.1111/1346-8138.14162 29687465
26. Jabs DA, Nussenblatt RB, Rosenbaum JT, Standardization of Uveitis Nomenclature (SUN) Working Group. Standardization of uveitis nomenclature for reporting clinical data. Results of the First International Workshop. Am J Ophthalmol. 2005;140: 509–516. doi: 10.1016/j.ajo.2005.03.057 16196117
27. Nussenblatt RB, Palestine AG, Chan CC, Roberge F. Standardization of vitreal inflammatory activity in intermediate and posterior uveitis. Ophthalmology. 1985;92: 467–471. doi: 10.1016/s0161-6420(85)34001-0 4000641
28. Rizou C, Ioannidis JP, Panou-Pomonis E, Sakarellos-Daitsiotis M, Sakarellos C, Moutsopoulos HM, et al. B-cell epitope mapping of DNA topoisomerase I defines epitopes strongly associated with pulmonary fibrosis in systemic sclerosis. Am J Respir Cell Mol Biol. 2000;22: 344–351. doi: 10.1165/ajrcmb.22.3.3850 10696071
29. Veeraraghavan S, Renzoni EA, Jeal H, Jones M, Hammer J, Wells AU, et al. Mapping of the immunodominant T cell epitopes of the protein topoisomerase I. Ann Rheum Dis. 2004;63: 982–987. doi: 10.1136/ard.2003.008037 15249326
30. Kongkaew S, Rungrotmongkol T, Punwong C, Noguchi H, Takeuchi F, Kungwan N, et al. Interactions of HLA-DR and topoisomerase I epitope modulated genetic risk for systemic sclerosis. Sci Rep. 2019;9: 745. doi: 10.1038/s41598-018-37038-z 30679605