Nonstandard peptide binding revealed by crystal structures of HLA-B*5101 complexed with HIV immunodominant epitopes

Class I HLA Alleles Are Associated With an Increased Risk of Osimertinib-Induced Hypersensitivity

BACKGROUND

Osimertinib, a third-generation epidermal growth factor receptor-tyrosine kinase inhibitor (EGFR-TKI), shows superior lung cancer treatment efficacy. However, osimertinib-induced severe hypersensitivity, including Stevens-Johnson syndrome (SJS)/toxic epidermal necrolysis (TEN), is frequently observed in Asian populations and hinders cancer treatment.

OBJECTIVE

We investigated the genetic HLA predisposition and immune pathomechanism of osimertinib-induced hypersensitivity.

METHODS

We enrolled 17 patients with osimertinib-induced delayed hypersensitivity (seven with severe SJS/TEN and 10 with mild maculopapular exanthema), 98 osimertinib-tolerant subjects, and 2,123 general population controls. We performed HLA genotyping, drug-induced lymphocyte activation test, and surface plasmon resonance assay.

RESULTS

HLA-B∗51:02 was present in 83.3% of osimertinib-induced SJS/TEN patients but in only 3.3% of the general population controls (P = 2.8 × 10-7; corrected P = 6.9 × 10-6; odds ratio [OR] = 146), and 0% of osimertinib-tolerant controls (P = 6.5 × 10-8; corrected P = 1.6 × 10-6; OR = 707). The association of HLA-B∗51:01 and HLA-A∗24:02 with osimertinib-induced maculopapular exanthema patients, rather than with osimertinib-tolerant subjects (P = .002, OR = 15.7 for HLA-B∗51:01; and P = .003, OR = 9.5 for HLA-A∗24:02), was identified as a phenotype-specific association. Granulysin, the SJS/TEN-specific cytotoxic protein, was significantly higher in plasma of SJS/TEN patients (39.8 ± 4.5 ng/mL; P < .001) and in in vitro lymphocyte activation test (sensitivity = 83.3%; P < .01) compared with tolerant controls. Patients with osimertinib-induced hypersensitivity appeared to tolerate alternative EGFR-TKIs. Surface plasmon resonance results also confirmed that HLA-B∗51:02 protein has a higher binding affinity for osimertinib and lower or no affinity for other EGFR-TKIs.

CONCLUSIONS

HLA-B∗51:02 frequently occurs in Asian populations and is strongly associated with osimertinib-induced SJS/TEN. Our findings suggest HLA-B∗51:02 screening as a preemptive test to reduce osimertinib-induced severe hypersensitivity.

TripHLApan: predicting HLA molecules binding peptides based on triple coding matrix and transfer learning

Human leukocyte antigen (HLA) recognizes foreign threats and triggers immune responses by presenting peptides to T cells. Computationally modeling the binding patterns between peptide and HLA is very important for the development of tumor vaccines. However, it is still a big challenge to accurately predict HLA molecules binding peptides. In this paper, we develop a new model TripHLApan for predicting HLA molecules binding peptides by integrating triple coding matrix, BiGRU + Attention models, and transfer learning strategy. We have found the main interaction site regions between HLA molecules and peptides, as well as the correlation between HLA encoding and binding motifs. Based on the discovery, we make the preprocessing and coding closer to the natural biological process. Besides, due to the input being based on multiple types of features and the attention module focused on the BiGRU hidden layer, TripHLApan has learned more sequence level binding information. The application of transfer learning strategies ensures the accuracy of prediction results under special lengths (peptides in length 8) and model scalability with the data explosion. Compared with the current optimal models, TripHLApan exhibits strong predictive performance in various prediction environments with different positive and negative sample ratios. In addition, we validate the superiority and scalability of TripHLApan's predictive performance using additional latest data sets, ablation experiments and binding reconstitution ability in the samples of a melanoma patient. The results show that TripHLApan is a powerful tool for predicting the binding of HLA-I and HLA-II molecular peptides for the synthesis of tumor vaccines. TripHLApan is publicly available at https://github.com/CSUBioGroup/TripHLApan.git.

Analysis of the different subpeptidomes presented by the HLA class I molecules of the B7 supertype

MHC-I molecules of the HLA-B7 supertype preferentially bind peptides with proline at position 2. HLA-B*51:01 and B*51:08 present two predominant subpeptidomes, one with Pro2 and hydrophobic residues at P1, and another with Ala2 and Asp enriched at position 1. Here, we present a meta-analysis of the peptidomes presented by molecules of the B7 supertype to investigate the presence of subpeptidomes across different allotypes. Several allotypes presented subpeptidomes differing in the presence of Pro or another residue at P2. The Ala2 subpeptidomes preferred Asp1 except in HLA-B*54:01, where ligands with Ala2 contained Glu1. Sequence alignment and the analysis of crystal structures allowed us to propose positions 45 and 67 of the MHC heavy chain as relevant for the presence of subpeptidomes. Deciphering the principles behind the presence of subpeptidomes could improve our understanding of antigen presentation in other MHC-I molecules. Running title: HLA-B7 supertype subpeptidomes.

Engagement with the TCR induces plasticity in antigenic ligands bound to MHC class I and CD1 molecules

Complementarity-determining regions (CDRs) of αβ T-cell receptors (TCRs) sense peptide-bound MHC (pMHC) complexes via chemical interactions, thereby mediating antigen specificity and MHC restriction. Flexible finger-like movement of CDR loops contributes to the establishment of optimal interactions with pMHCs. In contrast, peptide ligands captured in MHC molecules are considered more static because of the rigid hydrogen-bond network that stabilizes peptide ligands in the antigen-binding groove of MHC molecules. An array of crystal structures delineating pMHC complexes in TCR-docked and TCR-undocked forms is now available, which enables us to assess TCR engagement-induced conformational changes in peptide ligands. In this short review, we overview conformational changes in MHC class I-bound peptide ligands upon TCR docking, followed by those for CD1-bound glycolipid ligands. Finally, we analyze the co-crystal structure of the TCR:lipopeptide-bound MHC class I complex that we recently reported. We argue that TCR engagement-induced conformational changes markedly occur in lipopeptide ligands, which are essential for exposure of a primary T-cell epitope to TCRs. These conformational changes are affected by amino acid residues, such as glycine, that do not interact directly with TCRs. Thus, ligand recognition by specific TCRs involves not only T-cell epitopes but also non-epitopic amino acid residues. In light of their critical function, we propose to refer to these residues as non-epitopic residues affecting ligand plasticity and antigenicity (NR-PA).

Integrative immunoinformatics paradigm for predicting potential B-cell and T-cell epitopes as viable candidates for subunit vaccine design against COVID-19 virulence

Background

The increase in global mortality rates from SARS-COV2 (COVID-19) infection has been alarming thereby necessitating the continual search for viable therapeutic interventions. Due to minimal microbial components, subunit (peptide-based) vaccines have demonstrated improved efficacies in stimulating immunogenic responses by host B- and T-cells.

Methods

Integrative immunoinformatics algorithms were used to determine linear and discontinuous B-cell epitopes from the S-glycoprotein sequence. End-point selection of the most potential B-cell epitope was based on highly essential physicochemical attributes. NetCTL-I and NetMHC-II algorithms were used to predict probable MHC-I and II T-cell epitopes for globally frequent HLA-A∗O2:01, HLA-B∗35:01, HLA-B∗51:01 and HLA-DRB1∗15:02 molecules. Highly probable T-cell epitopes were selected based on their high propensities for C-terminal cleavage, transport protein (TAP) processing and MHC-I/II binding.

Results

Preferential epitope binding sites were further identified on the HLA molecules using a blind peptide-docking method. Phylogenetic analysis revealed close relativity between SARS-CoV-2 and SARS-CoV S-protein. LALHRSYLTPGDSSSGWTAGAA_242→263 was the most probable B-cell epitope with optimal physicochemical attributes. MHC-I antigenic presentation pathway was highly favourable for YLQPRTFLL_269-277 (HLA-A∗02:01), LPPAYTNSF_24-32 (HLA-B∗35:01) and IPTNFTISV7_14-721 (HLA-B∗51:01). Also, LTDEMIAQYTSALLA_865-881 exhibited the highest binding affinity to HLA-DR B1∗15:01 with core interactions mediated by IAQYTSALL_870-878. COVID-19 YLQPRTFLL_269-277 was preferentially bound to a previously undefined site on HLA-A∗02:01 suggestive of a novel site for MHC-I-mediated T-cell stimulation.

Conclusion

This study implemented combinatorial immunoinformatics methods to model B- and T-cell epitopes with high potentials to trigger immunogenic responses to the S protein of SARS-CoV-2.

Mining the Mycobacterium tuberculosis proteome for identification of potential T-cell epitope based vaccine candidates

Identification of protective antigens for designing a high-efficacy tuberculosis vaccine is the need of the hour. Till date only 7% of the Mycobacterium tuberculosis proteome has been explored for discovering antigens capable of activating T-cell responses. Therefore, it becomes crucial to screen the remaining Mycobacterium tuberculosis proteome for more immunodominant T-cell epitopes. An extensive knowledge of the epitopes recognized by our immune system can aid this process of finding potential T cell antigens for development of a better TB vaccine. In the present in-silico study, 237 proteins belonging to the 'virulence, detoxification, and adaptation' category of Mycobacterium tuberculosis proteome were targeted for T-cell epitope screening. 50825 MHC Class I and 49357 MHC Class II epitopes were generated using NetMHC3.4 and IEDB servers respectively and tested for their antigenicity and cytokine stimulation. The highest antigenic epitopes were analyzed for their world population coverage and epitope conservancy. Molecular docking and molecular dynamics simulation studies were performed to corroborate the binding affinities and structural stability of the peptide-MHC complexes. We predicted a total of 3 MHC Class I (ILLKMCWPA, FAVGMNVYV, and SLAGNSAKV) and 7 MHC Class II (DLTIGFFLHIPFPPV, RPDLTIGFFLHIPFP, LTIGFFLHIPFPPVE, VLVFALVVALVYLQF, LVFALVVALVYLQFR, PNLVAARFIQLTPVY, and LVLVFALVVALVYLQ) epitopes that can be promising vaccine candidates. These predicted epitopes belong to 6 distinct proteins: Rv0169 (mce1a), Rv3490 (ostA), Rv3496 (mce4D), Rv1085c, Rv0563 (HtpX), Rv3497c (mce4C). All these proteins are expressed at different stages in the life cycle of Mycobacterium tuberculosis and thus, the predicted epitopes could be employed as candidates for designing a multistage-multiepitopic vaccine.

A comprehensive screening of the whole proteome of hantavirus and designing a multi-epitope subunit vaccine for cross-protection against hantavirus: Structural vaccinology and immunoinformatics study

Hantaviruses are an emerging zoonotic group of rodent-borne viruses that are having serious implications on global public health due to the increase in outbreaks. Since there is no permanent cure, there is increasing interest in developing a vaccine against the hantavirus. This research aimed to design a robust cross-protective subunit vaccine using a novel immunoinformatics approach. After careful evaluation, the best predicted cytotoxic & helper T-cell and B-cell epitopes from nucleocapsid proteins, glycoproteins, RdRp proteins, and non-structural proteins were considered as potential vaccine candidates. Among the four generated vaccine models with different adjuvant, the model with toll-like receptor-4 (TLR-4) agonist adjuvant was selected because of its high antigenicity, non-allergenicity, and structural quality. The selected model was 654 amino acids long and had a molecular weight of 70.5 kDa, which characterizes the construct as a good antigenic vaccine candidate. The prediction of the conformational B-lymphocyte (CBL) epitope secured its ability to induce the humoral response. Thereafter, disulfide engineering improved vaccine stability. Afterwards, the molecular docking confirmed a good binding affinity of -1292 kj/mol with considered immune receptor TLR-4 and the dynamics simulation showed high stability of the vaccine-receptor complex. Later, the in silico cloning confirmed the better expression of the constructed vaccine protein in E. coli K12. Finally, in in silico immune simulation, significantly high levels of immunoglobulin M (IgM), immunoglobulin G1 (IgG1), cytotoxic & helper T lymphocyte (CTL & HTL) populations, and numerous cytokines such as interferon-γ (IFN-γ), interleukin-2 (IL-2) etc. were found as coherence with actual immune response and also showed faster antigen clearance for repeated exposures. Nonetheless, experimental validation can demonstrate the safety and cross-protective ability of the proposed vaccine to fight against hantavirus infection.

Molecular dynamics simulations provide molecular insights into the role of HLA-B51 in Behçet's disease pathogenesis

Behçet's disease is an inflammatory disorder of unknown etiology. Genetic tendency has an important role in its pathogenesis, and HLA-B51, a class I MHC antigen, has been recognized as the strongest susceptibility factor for Behçet's disease. Despite the confirmation of the association of HLA-B51 with Behçet's disease in different populations, its pathogenic mechanisms remain elusive. HLA-B51 differs in only two amino acids from HLA-B52, other split antigen of HLA-B5, which is not associated with Behçet's disease. These two amino acids are located in the B pocket of the antigen-binding groove, which occupies the second amino acids of the bound peptides. To understand the nature of the HLA-peptide interactions, differences in structure and dynamics of two HLA alleles were investigated by molecular dynamics simulations using YAYDGKDYI, LPRSTVINI, and IPYQDLPHL peptides. For HLA-B51, all bound peptides fluctuated to larger extent than HLA-B52. Free energy profiles of unbinding process for YAYDGKDYI by steered molecular dynamics simulations showed that unbinding from HLA-B52 results in greater free energy differences than HLA-B51. These results suggest the possibility of an instability of HLA-B51 associated with the repertoire of peptides, and this finding may provide significant insight to its pathogenic role in Behçet's disease.

HLA-Arena: A Customizable Environment for the Structural Modeling and Analysis of Peptide-HLA Complexes for Cancer Immunotherapy

PURPOSE

HLA protein receptors play a key role in cellular immunity. They bind intracellular peptides and display them for recognition by T-cell lymphocytes. Because T-cell activation is partially driven by structural features of these peptide-HLA complexes, their structural modeling and analysis are becoming central components of cancer immunotherapy projects. Unfortunately, this kind of analysis is limited by the small number of experimentally determined structures of peptide-HLA complexes. Overcoming this limitation requires developing novel computational methods to model and analyze peptide-HLA structures.

METHODS

Here we describe a new platform for the structural modeling and analysis of peptide-HLA complexes, called HLA-Arena, which we have implemented using Jupyter Notebook and Docker. It is a customizable environment that facilitates the use of computational tools, such as APE-Gen and DINC, which we have previously applied to peptide-HLA complexes. By integrating other commonly used tools, such as MODELLER and MHCflurry, this environment includes support for diverse tasks in structural modeling, analysis, and visualization.

RESULTS

To illustrate the capabilities of HLA-Arena, we describe 3 example workflows applied to peptide-HLA complexes. Leveraging the strengths of our tools, DINC and APE-Gen, the first 2 workflows show how to perform geometry prediction for peptide-HLA complexes and structure-based binding prediction, respectively. The third workflow presents an example of large-scale virtual screening of peptides for multiple HLA alleles.

CONCLUSION

These workflows illustrate the potential benefits of HLA-Arena for the structural modeling and analysis of peptide-HLA complexes. Because HLA-Arena can easily be integrated within larger computational pipelines, we expect its potential impact to vastly increase. For instance, it could be used to conduct structural analyses for personalized cancer immunotherapy, neoantigen discovery, or vaccine development.

Epitope-based peptide vaccines predicted against novel coronavirus disease caused by SARS-CoV-2

•

Linear B-cell epitopes in RBD of S protein predicted against SARS-CoV-2.

•

Discontinuous B-cell epitopes from S protein predicted against SARS-CoV-2.

•

T-cell epitopes from S, M and N protein predicted against SARS-CoV-2.

The outbreak of the 2019 novel coronavirus (SARS-CoV-2) has infected millions of people with a large number of deaths across the globe. The existing therapies are limited in dealing with SARS-CoV-2 due to the sudden appearance of the virus. Therefore, vaccines and antiviral medicines are in desperate need. We took immune-informatics approaches to identify B- and T-cell epitopes for surface glycoprotein (S), membrane glycoprotein (M) and nucleocapsid protein (N) of SARS-CoV-2, followed by estimating their antigenicity and interactions with the human leukocyte antigen (HLA) alleles. Allergenicity, toxicity, physiochemical properties analysis and stability were examined to confirm the specificity and selectivity of the epitope candidates. We identified a total of five B cell epitopes in RBD of S protein, seven MHC class-I, and 18 MHC class-II binding T-cell epitopes from S, M and N protein which showed non-allergenic, non-toxic and highly antigenic features and non-mutated in 55,179 SARS-CoV-2 virus strains until June 25, 2020. The epitopes identified here can be a potentially good candidate repertoire for vaccine development.

Informatics investigations into anti-thyroid drug induced agranulocytosis associated with multiple HLA-B alleles

INTRODUCTION

Adverse drug reactions have been linked with HLA alleles in different studies. These HLA proteins play an essential role in the adaptive immune response for the presentation of self and non-self peptides. Anti-thyroid drugs methimazole and propylthiouracil have been associated with drug induced agranulocytosis (severe lower white blood cell count) in patients with B*27:05, B*38:02 and DRB1*08:03 alleles in different populations: Taiwanese, Vietnamese, Han Chinese and Caucasian.

METHODS

In this study, informatics methods were used to investigate if any sequence or structural similarities exist between the two associated HLA-B alleles, compared with a set of "control" alleles assumed not be associated, which could help explain the molecular basis of the adverse drug reaction. We demonstrated using MHC Motif Viewer and MHCcluster that the two alleles do not have a propensity to bind similar peptides, and thus at a gross level the structure of the antigen presentation region of the two alleles are not similar. We also performed multiple sequence alignment to identify polymorphisms shared by the risk but not by the control alleles and molecular docking to compare the predicted binding poses of the drug-allele combinations.

RESULTS

Two residues, Cys67 and Thr80, were identified from the multiple sequence alignments to be unique to these risk alleles alone. The molecular docking showed the poses of the risk alleles to favour the F-pocket of the peptide binding groove, close to the Thr80 residue, with the control alleles generally favouring a different pocket. The data are thus suggestive that Thr80 may be a critical residue in HLA-mediated anti-thyroid drug induced agranulocytosis, and thus can guide future research and risk assessment.

A large peptidome dataset improves HLA class I epitope prediction across most of the human population

Prediction of HLA epitopes is important for the development of cancer immunotherapies and vaccines. However, current prediction algorithms have limited predictive power, in part because they were not trained on high-quality epitope datasets covering a broad range of HLA alleles. To enable prediction of endogenous HLA class I-associated peptides across a large fraction of the human population, we used mass spectrometry to profile >185,000 peptides eluted from 95 HLA-A, -B, -C and -G mono-allelic cell lines. We identified canonical peptide motifs per HLA allele, unique and shared binding submotifs across alleles and distinct motifs associated with different peptide lengths. By integrating these data with transcript abundance and peptide processing, we developed HLAthena, providing allele-and-length-specific and pan-allele-pan-length prediction models for endogenous peptide presentation. These models predicted endogenous HLA class I-associated ligands with 1.5-fold improvement in positive predictive value compared with existing tools and correctly identified >75% of HLA-bound peptides that were observed experimentally in 11 patient-derived tumor cell lines.

Redundancy and Complementarity between ERAP1 and ERAP2 Revealed by their Effects on the Behcet's Disease-associated HLA-B*51 Peptidome

The endoplasmic reticulum aminopeptidases ERAP1 and ERAP2 trim peptides to be loaded onto HLA molecules, including the main risk factor for Behçet's disease HLA-B*51. ERAP1 is also a risk factor among HLA-B*51-positive individuals, whereas no association is known with ERAP2. This study addressed the mutual relationships between both enzymes in the processing of an HLA-bound peptidome, interrogating their differential association with Behçet's disease. CRISPR/Cas9 was used to generate knock outs of ERAP1, ERAP2 or both from transfectant 721.221-HLA-B*51:01 cells. The surface expression of HLA-B*51 was reduced in all cases. The effects of depleting each or both enzymes on the B*51:01 peptidome were analyzed by quantitative label-free mass spectrometry. Substantial quantitative alterations of peptide length, subpeptidome balance, N-terminal residue usage, affinity and presentation of noncanonical ligands were observed. These effects were often different in the presence or absence of the other enzyme, revealing their mutual dependence. In the absence of ERAP1, ERAP2 showed similar and significant processing of B*51:01 ligands, indicating functional redundancy. The high overlap between the peptidomes of wildtype and double KO cells indicates that a large majority of B*51:01 ligands are present in the ER even in the absence of ERAP1/ERAP2. These results indicate that both enzymes have distinct, but complementary and partially redundant effects on the B*51:01 peptidome, leading to its optimization and maximal surface expression. The distinct effects of both enzymes on the HLA-B*51 peptidome provide a basis for their differential association with Behçet's disease and suggest a pathogenetic role of the B*51:01 peptidome.

Identification and Structure of an MHC Class I-Encoded Protein with the Potential to Present N-Myristoylated 4-mer Peptides to T Cells

Similar to host proteins, N-myristoylation occurs for viral proteins to dictate their pathological function. However, this lipid-modifying reaction creates a novel class of "lipopeptide" Ags targeted by host CTLs. The primate MHC class I-encoded protein, Mamu-B*098, was previously shown to bind N-myristoylated 5-mer peptides. Nevertheless, T cells exist that recognize even shorter lipopeptides, and much remains to be elucidated concerning the molecular mechanisms of lipopeptide presentation. We, in this study, demonstrate that the MHC class I allele, Mamu-B*05104, binds the N-myristoylated 4-mer peptide (C14-Gly-Gly-Ala-Ile) derived from the viral Nef protein for its presentation to CTLs. A phylogenetic tree analysis indicates that these classical MHC class I alleles are not closely associated; however, the high-resolution x-ray crystallographic analyses indicate that both molecules share lipid-binding structures defined by the exceptionally large, hydrophobic B pocket to accommodate the acylated glycine (G1) as an anchor. The C-terminal isoleucine (I4) of C14-Gly-Gly-Ala-Ile anchors at the F pocket, which is distinct from that of Mamu-B*098 and is virtually identical to that of the peptide-presenting MHC class I molecule, HLA-B51. The two central amino acid residues (G2 and A3) are only exposed externally for recognition by T cells, and the methyl side chain on A3 constitutes a major T cell epitope, underscoring that the epitopic diversity is highly limited for lipopeptides as compared with that for MHC class I-presented long peptides. These structural features suggest that lipopeptide-presenting MHC class I alleles comprise a distinct MHC class I subset that mediates an alternative pathway for CTL activation.

Variations in HLA-B cell surface expression, half-life and extracellular antigen receptivity

The highly polymorphic human leukocyte antigen (HLA) class I molecules present peptide antigens to CD8+ T cells, inducing immunity against infections and cancers. Quality control mediated by peptide loading complex (PLC) components is expected to ensure the cell surface expression of stable peptide-HLA class I complexes. This is exemplified by HLA-B*08:01 in primary human lymphocytes, with both expression level and half-life at the high end of the measured HLA-B expression and stability hierarchies. Conversely, low expression on lymphocytes is measured for three HLA-B allotypes that bind peptides with proline at position 2, which are disfavored by the transporter associated with antigen processing. Surprisingly, these lymphocyte-specific expression and stability differences become reversed or altered in monocytes, which display larger intracellular pools of HLA class I than lymphocytes. Together, the findings indicate that allele and cell-dependent variations in antigen acquisition pathways influence HLA-B surface expression levels, half-lives and receptivity to exogenous antigens.

General Prediction of Peptide-MHC Binding Modes Using Incremental Docking: A Proof of Concept

The class I major histocompatibility complex (MHC) is capable of binding peptides derived from intracellular proteins and displaying them at the cell surface. The recognition of these peptide-MHC (pMHC) complexes by T-cells is the cornerstone of cellular immunity, enabling the elimination of infected or tumoral cells. T-cell-based immunotherapies against cancer, which leverage this mechanism, can greatly benefit from structural analyses of pMHC complexes. Several attempts have been made to use molecular docking for such analyses, but pMHC structure remains too challenging for even state-of-the-art docking tools. To overcome these limitations, we describe the use of an incremental meta-docking approach for structural prediction of pMHC complexes. Previous methods applied in this context used specific constraints to reduce the complexity of this prediction problem, at the expense of generality. Our strategy makes no assumption and can potentially be used to predict binding modes for any pMHC complex. Our method has been tested in a re-docking experiment, reproducing the binding modes of 25 pMHC complexes whose crystal structures are available. This study is a proof of concept that incremental docking strategies can lead to general geometry prediction of pMHC complexes, with potential applications for immunotherapy against cancer or infectious diseases.

Is Behçet's disease a 'class 1-opathy'? The role of HLA-B*51 in the pathogenesis of Behçet's disease

The association between carriage of the human leucocyte antigen (HLA)-B*51 allele and development of Behçet's disease (BD) has been known since the early 1970s, but the exact mechanisms responsible for its role in pathogenesis remain much-debated. In an effort to explain the disease process, it has been suggested that BD constitutes one of a newly termed group of diseases, the 'MHC-I-opathies'. Other MHC-I-opathies include ankylosing spondylitis and HLA-B*27-associated spondyloarthropathies and HLA-C*0602-associated skin psoriasis. Recent work analysing the peptidome of HLA-B*51 suggests that altered peptide presentation by HLA-B*51 is vital to the disease process. In this review, we argue that immune receptor interactions with HLA-B*51 or the HLA-B*51-peptide complex could lead to development of inflammation in BD. The evidence for CD8⁺ T cell involvement is weak, and based on emerging studies it seems more likely that natural killer (NK) or other cell interactions, perhaps mediated by leucocyte immunoglobulin-like receptor (LILR) or killer immunoglobulin-like receptor (KIR) receptors, are culpable in pathogenesis. HLA misfolding leading directly to inflammation is another hypothesis for BD pathogenesis that deserves greater investigation. Ultimately, greater understanding of HLA-B*51's unique role in BD will probably lead to improved development of therapeutic strategies.

The Behçet's disease-associated variant of the aminopeptidase ERAP1 shapes a low-affinity HLA-B*51 peptidome by differential subpeptidome processing

A low-activity variant of endoplasmic reticulum aminopeptidase 1 (ERAP1), Hap10, is associated with the autoinflammatory disorder Behçet's disease (BD) in epistasis with HLA-B*51, which is the main risk factor for this disorder. The role of Hap10 in BD pathogenesis is unknown. We sought to define the effects of Hap10 on the HLA-B*51 peptidome and to distinguish these effects from those due to HLA-B*51 polymorphisms unrelated to disease. The peptidome of the BD-associated HLA-B*51:08 subtype expressed in a Hap10-positive cell line was isolated, characterized by mass spectrometry, and compared with the HLA-B*51:01 peptidome from cells expressing more active ERAP1 allotypes. We additionally performed synthetic peptide digestions with recombinant ERAP1 variants and estimated peptide-binding affinity with standard algorithms. In the BD-associated ERAP1 context of B*51:08, longer peptides were generated; of the two major HLA-B*51 subpeptidomes with Pro-2 and Ala-2, the former one was significantly reduced, and the latter was increased and showed more ERAP1-susceptible N-terminal residues. These effects were readily explained by the low activity of Hap10 and the differential susceptibility of X-Pro and X-Ala bonds to ERAP1 trimming and together resulted in a significantly altered peptidome with lower affinity. The differences due to ERAP1 were clearly distinguished from those due to HLA-B*51 subtype polymorphism, which affected residue frequencies at internal positions of the peptide ligands. The alterations in the nature and affinity of HLA-B*51·peptide complexes probably affect T-cell and natural killer cell recognition, providing a sound basis for the joint association of ERAP1 and HLA-B*51 with BD.

HLA class I variation in Iranian Lur and Kurd populations: high haplotype and allotype diversity with an abundance of KIR ligands

HLA-A, -B and -C alleles of 285 individuals, representing three Iranian Lur populations and one Iranian Kurd population were sequenced completely, yielding human leukocyte antigen (HLA) class I genotypes at high resolution and filling four fields of the official HLA nomenclature. Each population has 87-99 alleles, evenly distributed between the three HLA class I genes, 145 alleles being identified in total. These alleles were already known, named and deposited in the HLA database. The alleles form 316 different HLA A-B-C haplotypes, with each population having between 80 and 112 haplotypes. The four Iranian populations form a related group that is distinguished from other populations, including other Iranians. All four KIR ligands - the A3/11, Bw4, C1 and C2 epitopes - are well represented, particularly Bw4, which is carried by three high-frequency allotypes: HLA-A*24:02, HLA-A*32:01 and HLA-B*51:01. In the Lur and Kurd populations, between 82% and 94% of individuals have the Bw4 epitope, the ligand for KIR3DL1. HLA-B*51:01 is likely of Neandertal origin and associated with Behcet's disease, also known as the Silk Road disease. The Lordegan Lur have the highest frequency of HLA-B*51:01 in the world. This allele is present on 46 Lur and Kurd haplotypes. Present at lower frequency is HLA-B*51:08, which is also associated with Behcet's disease. In the four Iranian populations, 31 haplotypes encode both Bw4(+) HLA-A and Bw4(+) HLA-B, a dual combination of Bw4 epitopes that is relatively rare in other populations, worldwide. This study both demonstrates and emphasizes the value of studying HLA class I polymorphism at highest resolution in anthropologically well-defined populations.

Genetic Analysis with the Immunochip Platform in Behçet Disease. Identification of Residues Associated in the HLA Class I Region and New Susceptibility Loci

Behcet's disease (BD) is an immuno-mediated vasculitis in which knowledge of its etiology and genetic basis is limited. To improve the current knowledge, a genetic analysis performed with the Immunochip platform was carried out in a population from Spain. A discovery cohort comprising 278 BD cases and 1,517 unaffected controls were genotyped using the Immunochip platform. The validation step was performed on an independent replication cohort composed of 130 BD cases and 600 additional controls. The strongest association signals were observed in the HLA class I region, being HLA-B*51 the highest peak (overall P = 6.82E-32, OR = 3.82). A step-wise conditional logistic regression with classical alleles identified HLA-B*57 and HLA-A*03 as additional independent markers. The amino acid model that best explained the association, includes the position 97 of the HLA-B molecule and the position 66 of the HLA-A. Among the non-HLA loci, the most significant in the discovery analysis were: IL23R (rs10889664: P = 3.81E-12, OR = 2.00), the JRKL/CNTN5 region (rs2848479: P = 5.00E-08, OR = 1.68) and IL12A (rs1874886: P = 6.67E-08, OR = 1.72), which were confirmed in the validation phase (JRKL/CNTN5 rs2848479: P = 3.29E-10, OR = 1.66; IL12A rs1874886: P = 1.62E-08, OR = 1.61). Our results confirm HLA-B*51 as a primary-association marker in predisposition to BD and suggest additional independent signals within the class I region, specifically in the genes HLA-A and HLA-B. Regarding the non-HLA genes, in addition to IL-23R, previously reported in our population; IL12A, described in other populations, was found to be a BD susceptibility factor also in Spaniards; finally, a new associated locus was found in the JRKL/CNTN5 region.

The Peptidome of Behçet's Disease-Associated HLA-B*51:01 Includes Two Subpeptidomes Differentially Shaped by Endoplasmic Reticulum Aminopeptidase 1

OBJECTIVE

To characterize the peptidome of the Behçet's disease-associated HLA-B*51:01 allotype as well as the differential features of major peptide subsets and their distinct endoplasmic reticulum aminopeptidase 1 (ERAP-1)-mediated processing.

METHODS

The endogenous B*51:01-bound peptidome was characterized from 721.221 transfectant cells, after affinity chromatography and acid extraction, by tandem mass spectrometry. Recombinant ERAP-1 variants were used to digest synthetic B*51:01 ligands. HLA and transporter associated with antigen processing (TAP) binding affinities of peptide ligands were calculated with well-established algorithms. ERAP-1 and ERAP-2 from 721.221 cells were characterized by genomic sequencing and Western blotting.

RESULTS

The B*51:01 peptidome consisted of 29.5% octamers, 61.7% nonamers, 4.8% decamers, and 4.0% longer peptides. The major peptide motif consisted of Pro and Ala at position 2, aliphatic/aromatic position 3 residues, and Val and Ile at the C-terminal position. The ligands with Pro or Ala at position 2 constituted 2 distinct subpeptidomes. Peptides with Pro at position 2 showed higher affinity for B*51:01 and lower affinity for TAP than those with Ala at position 2. Most important, both peptide subsets differed drastically in the susceptibility of their position 1 residues to ERAP-1, revealing a distinct influence of this enzyme on both subpeptidomes, which may alter their balance, affecting the global affinity of B*51:01-peptide complexes.

CONCLUSION

ERAP-1 has a significant influence on the B*51:01 peptidome and its affinity. This influence is based on very distinct effects on the 2 subpeptidomes, whereby only peptides in the subpeptidome with Ala at position 2 are extensively destroyed, except when their position 1 residues are ERAP-1 resistant. This pattern provides a mechanism for the epistatic association of ERAP-1 and B*51:01 in Behçet's disease.

Vaccine Design for H5N1 Based on B- and T-cell Epitope Predictions

From 2003 to 2013, Indonesia had the highest number of avian influenza A cases in humans, with 192 cases and 160 fatalities. Avian influenza is caused by influenza virus type A, such as subtype H5N1. This virus has two glycoproteins: hemagglutinin and neuraminidase, which will become the primary target to be neutralized by vaccine. Vaccine is the most effective immunologic intervention. In this study, we use the epitope-based vaccine design from hemagglutinin and neuraminidase of H5N1 Indonesian strain virus by using immunoinformatics approach in order to predict the binding of B-cell and T-cell epitopes (class I and class II human leukocyte antigen [HLA]). BCPREDS was used to predict the B-cell epitope. Propred, Propred I, netMHCpan, and netMHCIIpan were used to predict the T-cell epitope. Two B-cell epitopes of hemagglutinin candidates and one B-cell epitope of neuraminidase candidates were obtained to bind T-cell CD4(+) (class II HLA), and also five T-cell epitope hemagglutinin and four T-cell epitope neuraminidase were obtained to bind T-cell CD8(+) (class I HLA). The visualization of epitopes was done using MOE 2008.10. It shows that the binding affinity of epitope-HLA was based on minimum binding free energy (ΔG binding). Based on this result, visualization, and dynamic simulation, four hemagglutinin epitopes (MEKIVLLLA, CPYLGSPSF, KCQTPMGAI, and IGTSTLNQR) and two neuraminidase epitopes (NPNQKIITI and CYPDAGEIT) were computed as having the best binding affinity from HLA ligand. The results mentioned above are from in silico experiments and need to be validated using wet experiment.

Molecular Dynamics Simulation Reveals the Selective Binding of Human Leukocyte Antigen Alleles Associated with Behçet's Disease

Behçet’s disease (BD), a multi-organ inflammatory disorder, is associated with the presence of the human leukocyte antigen (HLA) HLA-B*51 allele in many ethnic groups. The possible antigen involvement of the major histocompatibility complex class I chain related gene A transmembrane (MICA-TM) nonapeptide (AAAAAIFVI) has been reported in BD symptomatic patients. This peptide has also been detected in HLA-A*26:01 positive patients. To investigate the link of BD with these two specific HLA alleles, molecular dynamics (MD) simulations were applied on the MICA-TM nonapeptide binding to the two BD-associated HLA alleles in comparison with the two non-BD-associated HLA alleles (B*35:01 and A*11:01). The MD simulations were applied on the four HLA/MICA-TM peptide complexes in aqueous solution. As a result, stabilization for the incoming MICA-TM was found to be predominantly contributed from van der Waals interactions. The P2/P3 residue close to the N-terminal and the P9 residue at the C-terminal of the MICA-TM nonapeptide served as the anchor for the peptide accommodated at the binding groove of the BD associated HLAs. The MM/PBSA free energy calculation predicted a stronger binding of the HLA/peptide complexes for the BD-associated HLA alleles than for the non-BD-associated ones, with a ranked binding strength of B*51:01 > B*35:01 and A*26:01 > A*11:01. Thus, the HLAs associated with BD pathogenesis expose the binding efficiency with the MICA-TM nonapeptide tighter than the non-associated HLA alleles. In addition, the residues 70, 73, 99, 146, 147 and 159 of the two BD-associated HLAs provided the conserved interaction for the MICA-TM peptide binding.

HLA-B*59:01: a marker for Stevens-Johnson syndrome/toxic epidermal necrolysis caused by methazolamide in Han Chinese

Methazolamide is an intraocular pressure-lowering drug that is used in the treatment of glaucoma and other ophthalmologic abnormalities. The use of methazolamide has been shown to cause Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) in patients of Asian ancestry. Methazolamide-induced SJS/TEN is associated with the presence of HLA-B59 serotype/HLA-B*59:01 in Korean and Japanese populations. To better understand the genetic risk factors for these adverse reactions in the Han Chinese population, we characterized the HLA class I genotypes of eight Chinese patients with methazolamide-induced SJS/TEN from 2008 to 2014. The frequency of HLA-B*59:01 was 87.5% (7/8) in the case patients, which was significantly different from 0% (0/30) in the methazolamide-tolerant patients (odds ratio (OR)=305.0; P=6.3 × 10(-7)) and 0.35% (1/283) in healthy subjects from the human major histocompatibility complex database (OR=1974.0; P=2.0 × 10(-12)). HLA-C*01:02, which is closely linked to HLA-B*59:01, had a weaker but notable association with methazolamide-induced SJS/TEN compared with the tolerant controls (OR=12.1; P=0.016) and general population (OR=15.5; P=2.0 × 10(-3)). The distribution of the HLA-B*59:01-C*01:02 haplotype was also significantly different in cases and controls. This study demonstrated a strong association between HLA-B*59:01 and methazolamide-induced SJS/TEN in the Han Chinese population for the first time. Pretherapy screening for HLA-B*59:01 would be useful to reduce the risk of methazolamide-induced SJS/TEN.

T cell antigen receptor recognition of antigen-presenting molecules

The Major Histocompatibility Complex (MHC) locus encodes classical MHC class I and MHC class II molecules and nonclassical MHC-I molecules. The architecture of these molecules is ideally suited to capture and present an array of peptide antigens (Ags). In addition, the CD1 family members and MR1 are MHC class I-like molecules that bind lipid-based Ags and vitamin B precursors, respectively. These Ag-bound molecules are subsequently recognized by T cell antigen receptors (TCRs) expressed on the surface of T lymphocytes. Structural and associated functional studies have been highly informative in providing insight into these interactions, which are crucial to immunity, and how they can lead to aberrant T cell reactivity. Investigators have determined over thirty unique TCR-peptide-MHC-I complex structures and twenty unique TCR-peptide-MHC-II complex structures. These investigations have shown a broad consensus in docking geometry and provided insight into MHC restriction. Structural studies on TCR-mediated recognition of lipid and metabolite Ags have been mostly confined to TCRs from innate-like natural killer T cells and mucosal-associated invariant T cells, respectively. These studies revealed clear differences between TCR-lipid-CD1, TCR-metabolite-MR1, and TCR-peptide-MHC recognition. Accordingly, TCRs show remarkable structural and biological versatility in engaging different classes of Ag that are presented by polymorphic and monomorphic Ag-presenting molecules of the immune system.

Behçet disease-associated MHC class I residues implicate antigen binding and regulation of cell-mediated cytotoxicity

The HLA protein, HLA-B*51, encoded by HLA-B in MHC, is the strongest known genetic risk factor for Behçet disease (BD). Associations between BD and other factors within the MHC have been reported also, although strong regional linkage disequilibrium complicates their confident disentanglement from HLA-B*51. In the current study, we examined a combination of directly obtained and imputed MHC-region SNPs, directly obtained HLA-B locus types, and imputed classical HLA types with their corresponding polymorphic amino acid residues for association with BD in 1,190 cases and 1,257 controls. SNP mapping with logistic regression of the MHC identified the HLA-B/MICA region and the region between HLA-F and HLA-A as independently associated with BD (P < 1.7 × 10(-8)). HLA-B*51, -A*03, -B*15, -B*27, -B*49, -B*57, and -A*26 each contributed independently to BD risk. We directly examined rs116799036, a noncoding SNP upstream of HLA-B that was recently suggested to underlie the association of HLA-B*51 with BD, but we were unable to replicate that finding in our collection. Instead, we mapped the BD association to seven MHC class I (MHC-I) amino acid residues, including anchor residues that critically define the selection and binding of peptides to MHC-I molecules, residues known to influence MHC-I-killer immunoglobulin-like receptor interactions, and a residue located in the signal peptide of HLA-B. The locations of these variants collectively implicate MHC-I peptide binding in the pathophysiology of BD. Furthermore, several lines of evidence suggest a role for altered regulation of cellular cytotoxicity in BD pathogenesis.

Distinct HIV-1 escape patterns selected by cytotoxic T cells with identical epitope specificity

Pol283-8-specific, HLA-B*51:01-restricted, cytotoxic T cells (CTLs) play a critical role in the long-term control of HIV-1 infection. However, these CTLs select for the reverse transcriptase (RT) I135X escape mutation, which may be accumulating in circulating HIV-1 sequences. We investigated the selection of the I135X mutation by CTLs specific for the same epitope but restricted by HLA-B*52:01. We found that Pol283-8-specific, HLA-B*52:01-restricted CTLs were elicited predominantly in chronically HIV-1-infected individuals. These CTLs had a strong ability to suppress the replication of wild-type HIV-1, though this ability was weaker than that of HLA-B*51:01-restricted CTLs. The crystal structure of the HLA-B*52:01-Pol283-8 peptide complex provided clear evidence that HLA-B*52:01 presents the peptide similarly to HLA-B*51:01, ensuring the cross-presentation of this epitope by both alleles. Population level analyses revealed a strong association of HLA-B*51:01 with the I135T mutant and a relatively weaker association of HLA-B*52:01 with several I135X mutants in both Japanese and predominantly Caucasian cohorts. An in vitro viral suppression assay revealed that the HLA-B*52:01-restricted CTLs failed to suppress the replication of the I135X mutant viruses, indicating the selection of these mutants by the CTLs. These results suggest that the different pattern of I135X mutant selection may have resulted from the difference between these two CTLs in the ability to suppress HIV-1 replication.

Diverse peptide presentation of rhesus macaque major histocompatibility complex class I Mamu-A 02 revealed by two peptide complex structures and insights into immune escape of simian immunodeficiency virus

Major histocompatibility complex class I (MHC I)-restricted CD8(+) T-cell responses play a pivotal role in anti-human immunodeficiency virus (HIV) immunity and the control of viremia. The rhesus macaque is an important animal model for HIV-related research. Among the MHC I alleles of the rhesus macaque, Mamu-A 02 is prevalent, presenting in ≥20% of macaques. In this study, we determined the crystal structure of Mamu-A 02, the second structure-determined MHC I from the rhesus macaque after Mamu-A 01. The peptide presentation characteristics of Mamu-A 02 are exhibited in complex structures with two typical Mamu-A 02-restricted CD8(+) T-cell epitopes, YY9 (Nef159 to -167; YTSGPGIRY) and GY9 (Gag71 to -79; GSENLKSLY), derived from simian immunodeficiency virus (SIV). These two peptides utilize similar primary anchor residues (Ser or Thr) at position 2 and Tyr at position 9. However, the central region of YY9 is different from that of GY9, a difference that may correlate with the immunogenic variance of these peptides. Further analysis indicated that the distinct conformations of these two peptides are modulated by four flexible residues in the Mamu-A 02 peptide-binding groove. The rare combination of these four residues in Mamu-A 02 leads to a variant presentation for peptides with different residues in their central regions. Additionally, in the two structures of the Mamu-A 02 complex, we compared the binding of rhesus and human β(2) microglobulin (β(2)m) to Mamu-A 02. We found that the peptide presentation of Mamu-A 02 is not affected by the interspecies interaction with human β(2)m. Our work broadens the understanding of CD8(+) T-cell-specific immunity against SIV in the rhesus macaque.

Two distinct conformations of a rinderpest virus epitope presented by bovine major histocompatibility complex class I N*01801: a host strategy to present featured peptides

The presentation of viral peptide epitopes to host cytotoxic T lymphocytes (CTLs) is crucial for adaptive cellular immunity to clear the virus infection, especially for some chronic viral infections. Indeed, hosts have developed effective strategies to achieve this goal. The ideal scenario would be that the peptide epitopes stimulate a broad spectrum of CTL responses with diversified T-cell receptor (TCR) usage (the TCR repertoire). It is believed that a diversified TCR repertoire requires a "featured" peptide to be presented by the host major histocompatibility complex (MHC). A featured peptide can be processed and presented in a number of ways. Here, using the X-ray diffraction method, the crystal structures of an antigenic peptide derived from rinderpest virus presented by bovine MHC class I N*01801 (BoLA-A11) have been solved, and two distinct conformations of the presented peptide are clearly displayed. A detailed analysis of the structure and comparative sequences revealed that the polymorphic amino acid isoleucine 73 (Ile73) is extremely flexible, allowing the MHC groove to adopt different conformations to accommodate the rinderpest virus peptide. This makes the peptide more featured by exposing different amino acids for T-cell recognition. The crystal structures also demonstrated that the N*01801 molecule has an unusually large A pocket, resulting in the special conformation of the P1 residue at the N terminus of the peptide. We propose that this strategy of host peptide presentation might be beneficial for creating a diversified TCR repertoire, which is important for a more-effective CTL response.

Long-term control of HIV-1 in hemophiliacs carrying slow-progressing allele HLA-B*5101

HLA-B*51 alleles are reported to be associated with slow disease progression to AIDS, but the mechanism underlying this association is still unclear. In the present study, we analyzed the effect of HLA-B*5101 on clinical outcome for Japanese hemophiliacs who had been infected with HIV-1 before 1985 and had been recruited in 1998 for this study. HLA-B*5101(+) hemophiliacs exhibited significantly slow progression. The analysis of HLA-B*5101-restricted HIV-1-specific cytotoxic T-lymphocyte (CTL) responses to 4 HLA-B*-restricted epitopes in 10 antiretroviral-therapy (ART)-free HLA-B*5101(+) hemophiliacs showed that the frequency of Pol283-8-specific CD8(+) T cells was inversely correlated with the viral load, whereas the frequencies of CD8(+) T cells specific for 3 other epitopes were positively correlated with the viral load. The HLA-B*5101(+) hemophiliacs whose HIV-1 replication had been controlled for approximately 25 years had HIV-1 possessing the wild-type Pol283-8 sequence or the Pol283-8V mutant, which does not critically affect T-cell recognition, whereas other HLA-B*5101(+) hemophiliacs had HIV-1 with escape mutations in this epitope. The results suggest that the control of HIV-1 over approximately 25 years in HLA-B*5101-positive hemophiliacs is associated with a Pol283-8-specific CD8(+) T-cell response and that lack of control of HIV-1 is associated with the appearance of Pol283-8-specific escape mutants.

Structural basis for T cell alloreactivity among three HLA-B14 and HLA-B27 antigens

The existence of cytotoxic T cells (CTL) cross-reacting with the human major histocompatibility antigens HLA-B14 and HLA-B27 suggests that their alloreactivity could be due to presentation of shared peptides in similar binding modes by these molecules. We therefore determined the crystal structures of the subtypes HLA-B*1402, HLA-B*2705, and HLA-B*2709 in complex with a proven self-ligand, pCatA (peptide with the sequence IRAAPPPLF derived from cathepsin A (residues 2-10)), and of HLA-B*1402 in complex with a viral peptide, pLMP2 (RRRWRRLTV, derived from latent membrane protein 2 (residues 236-244) of Epstein-Barr virus). Despite the exchange of 18 residues within the binding grooves of HLA-B*1402 and HLA-B*2705 or HLA-B*2709, the pCatA peptide is presented in nearly identical conformations. However, pLMP2 is displayed by HLA-B*1402 in a conformation distinct from those previously found in the two HLA-B27 subtypes. In addition, the complexes of HLA-B*1402 with the two peptides reveal a nonstandard, tetragonal mode of the peptide N terminus anchoring in the binding groove because of the exchange of the common Tyr-171 by His-171 of the HLA-B*1402 heavy chain. This exchange appears also responsible for reduced stability of HLA-B14-peptide complexes in vivo and slow assembly in vitro. The studies with the pCatA peptide uncover that CTL cross-reactive between HLA-B14 and HLA-B27 might primarily recognize the common structural features of the bound peptide, thus neglecting amino acid replacements within the rim of the binding grooves. In contrast, structural alterations between the three complexes with the pLMP2 peptide indicate how heavy chain polymorphisms can influence peptide display and prevent CTL cross-reactivity between HLA-B14 and HLA-B27 antigens.

Polymorphic sites away from the Bw4 epitope that affect interaction of Bw4+ HLA-B with KIR3DL1

KIR3DL1 is a polymorphic, inhibitory NK cell receptor specific for the Bw4 epitope carried by subsets of HLA-A and HLA-B allotypes. The Bw4 epitope of HLA-B*5101 and HLA-B*1513 is determined by the NIALR sequence motif at positions 77, 80, 81, 82, and 83 in the alpha(1) helix. Mutation of these positions to the residues present in the alternative and nonfunctional Bw6 motif showed that the functional activity of the Bw4 epitopes of B*5101 and B*1513 is retained after substitution at positions 77, 80, and 81, but lost after substitution of position 83. Mutation of leucine to arginine at position 82 led to loss of function for B*5101 but not for B*1513. Further mutagenesis, in which B*1513 residues were replaced by their B*5101 counterparts, showed that polymorphisms in all three extracellular domains contribute to this functional difference. Prominent were positions 67 in the alpha(1) domain, 116 in the alpha(2) domain, and 194 in the alpha(3) domain. Lesser contributions were made by additional positions in the alpha(2) domain. These positions are not part of the Bw4 epitope and include residues shaping the B and F pockets that determine the sequence and conformation of the peptides bound by HLA class I molecules. This analysis shows how polymorphism at sites throughout the HLA class I molecule can influence the interaction of the Bw4 epitope with KIR3DL1. This influence is likely mediated by changes in the peptides bound, which alter the conformation of the Bw4 epitope.

T cell responses to human endogenous retroviruses in HIV-1 infection

Human endogenous retroviruses (HERVs) are remnants of ancient infectious agents that have integrated into the human genome. Under normal circumstances, HERVs are functionally defective or controlled by host factors. In HIV-1-infected individuals, intracellular defense mechanisms are compromised. We hypothesized that HIV-1 infection would remove or alter controls on HERV activity. Expression of HERV could potentially stimulate a T cell response to HERV antigens, and in regions of HIV-1/HERV similarity, these T cells could be cross-reactive. We determined that the levels of HERV production in HIV-1-positive individuals exceed those of HIV-1-negative controls. To investigate the impact of HERV activity on specific immunity, we examined T cell responses to HERV peptides in 29 HIV-1-positive and 13 HIV-1-negative study participants. We report T cell responses to peptides derived from regions of HERV detected by ELISPOT analysis in the HIV-1-positive study participants. We show an inverse correlation between anti-HERV T cell responses and HIV-1 plasma viral load. In HIV-1-positive individuals, we demonstrate that HERV-specific T cells are capable of killing cells presenting their cognate peptide. These data indicate that HIV-1 infection leads to HERV expression and stimulation of a HERV-specific CD8+ T cell response. HERV-specific CD8+ T cells have characteristics consistent with an important role in the response to HIV-1 infection: a phenotype similar to that of T cells responding to an effectively controlled virus (cytomegalovirus), an inverse correlation with HIV-1 plasma viral load, and the ability to lyse cells presenting their target peptide. These characteristics suggest that elicitation of anti-HERV-specific immune responses is a novel approach to immunotherapeutic vaccination. As endogenous retroviral sequences are fixed in the human genome, they provide a stable target, and HERV-specific T cells could recognize a cell infected by any HIV-1 viral variant. HERV-specific immunity is an important new avenue for investigation in HIV-1 pathogenesis and vaccine design.

The human genome contains a number of remnants or fossils of ancient viral infections referred to as human endogenous retroviruses (HERV). Like fossils, these HERV are considered to be dead or inert in most cases. However, we demonstrate that T cells in the human immune system respond to HERV when a person is infected with the human immunodeficiency virus (HIV). The T cells responding to HERV share characteristics with T cells that effectively control cytomegalovirus, a common chronic viral infection. T cells responding to HERV can also kill target cells carrying HERV protein. For some HIV-positive people, the strength of their response against HERV is related to having a lower HIV viral load. This study has important implications for new directions in HIV vaccine research. One of the key obstacles to creating an effective HIV vaccine is overcoming the ability of some of the viral variants produced when HIV replicates to evade the immune responses that the body mounts to control infections. If T cells that recognize HERV can stably target HIV-infected cells, they could be an important factor in controlling HIV infection.

Lack of KIR3DS1 binding to MHC class I Bw4 tetramers in complex with CD8+ T cell epitopes

In HIV-1 infection, the synergistic association of a subset of Bw4 MHC class I molecules and the activating killer inhibitory receptor (KIR), KIR3DS1, with prolonged AIDS-free survival has been reported. As KIRs represent a diverse group of MHC class I receptors, we questioned whether Bw4 MHC class I molecules expressing isoleucine at position 80 (Bw4Ile80) and in complex with HIV-1-derived T cell epitopes represented KIR3DS1 ligands. MHC class I tetramers are powerful tools for the detection of T cell receptor-MHC class I interactions, and have recently been used to evaluate KIR-MHC class I binding ex vivo. Specifically, this approach has been successfully utilized to assess binding of Bw4 MHC class I tetramers to KIR3DL1, an inhibitory KIR and allele of KIR3DS1. In this study we generated a diverse panel of HIV-1-specific Bw4Ile80 MHC class I tetramers and tested its ability to bind transiently expressed KIR3DS1 on 293-T cells. Using flow cytometry analysis, the expression of KIR3DS1 on 293-T cells was confirmed by anti-FLAG BioM2 staining, prior to incubation with PE-conjugated MHC class I tetramers. Despite choosing a broad array of peptide epitopes and diverse Bw4Ile80 MHC class I molecules, we were unable to detect tetramer binding to KIR3DS1. We speculate that our negative finding may be a consequence of the MHC class I molecules and peptide epitopes chosen, but could also relate to key amino acid differences that distinguish KIR3DS1 from KIR3DL1.

Associations among Epstein-Barr virus subtypes, human leukocyte antigen class I alleles, and the development of posttransplantation lymphoproliferative disorder in bone marrow transplant recipients

The association between Epstein-Barr virus subtype, human leukocyte antigen class I alleles, and the development of posttransplantation lymphoproliferative disorder was examined in a group of 25 bone marrow transplant recipients. A highly statistically significant correlation was observed between the human leukocyte antigen B51 allele and development of posttransplantation lymphoproliferative disorder (P=.0016). This study provides, to our knowledge, the first evidence that the human leukocyte antigen B51 allele might predispose bone marrow transplant recipients to Epstein-Barr virus-associated posttransplantation lymphoproliferative disorder.

Cutting Edge: Allele-specific and peptide-dependent interactions between KIR3DL1 and HLA-A and HLA-B

Although it is clear that KIR3DL1 recognizes Bw4(+) HLA-B, the role of Bw4(+) HLA-A allotypes as KIR3DL1 ligands is controversial. We therefore examined the binding of tetrameric HLA-A and -B complexes, including HLA*2402, a common Bw4(+) HLA-A allotype, to KIR3DL1*001, *005, *007, and *1502 allotypes. Only Bw4(+) tetramers bound KIR3DL1. Three of four HLA-A*2402 tetramers bound one or more KIR3DL1 allotypes and all four KIR3DL1 allotypes bound to one or more HLA-A*2402 tetramers, but with different binding specificities. Only KIR3DL1*005 bound both HLA-A*2402 and HLA-B*5703 tetramers. HLA-A*2402-expressing target cells were resistant to lysis by NK cells expressing KIR3DL1*001 or *005. This study shows that HLA-A*2402 is a ligand for KIR3DL1 and demonstrates how the binding of KIR3DL1 to Bw4(+) ligands depends upon the bound peptide as well as HLA and KIR3DL1 polymorphism.

Sampling of major histocompatibility complex class I-associated peptidome suggests relatively looser global association of HLA-B*5101 with peptides

We have analyzed peptides associated with six human major histocompatibility complex (MHC) class I allomorphs expressed by the U937 cell line. Peptides were isolated by mild acid elution or by MHC class I immunoprecipitation by using W6/32 monoclonal antibody. Eighty-five peptides were sequenced by mass spectrometry, and their putative binding alleles were assigned using bioinformatic tools. Only three peptides isolated by the two approaches were identical, suggesting that the approaches may yield distinct partially overlapping peptide populations. Mild acid treatment-derived peptides manifested overall characteristics suggestive of relatively lower affinity of binding for MHC class I. Interestingly, a large proportion of putative HLA-B*5101-binding peptides was evident among the mild acid treatment-eluted peptides, and to a lesser degree in the affinity-purified peptide pool. These results suggest that HLA-B*5101 may bind a potentially large pool of peptides with relatively lower affinity. We suggest that lower affinity of peptide binding may be the basis for inefficient tolerance to HLA-B*5101-binding self-peptides, a predisposing factor for the development of Behçet disease.

Ab initio prediction of peptide-MHC binding geometry for diverse class I MHC allotypes

Since determining the crystallographic structure of all peptide-MHC complexes is infeasible, an accurate prediction of the conformation is a critical computational problem. These models can be useful for determining binding energetics, predicting the structures of specific ternary complexes with T-cell receptors, and designing new molecules interacting with these complexes. The main difficulties are (1) adequate sampling of the large number of conformational degrees of freedom for the flexible peptide, (2) predicting subtle changes in the MHC interface geometry upon binding, and (3) building models for numerous MHC allotypes without known structures. Whereas previous studies have approached the sampling problem by dividing the conformational variables into different sets and predicting them separately, we have refined the Biased-Probability Monte Carlo docking protocol in internal coordinates to optimize a physical energy function for all peptide variables simultaneously. We also imitated the induced fit by docking into a more permissive smooth grid representation of the MHC followed by refinement and reranking using an all-atom MHC model. Our method was tested by a comparison of the results of cross-docking 14 peptides into HLA-A*0201 and 9 peptides into H-2K(b) as well as docking peptides into homology models for five different HLA allotypes with a comprehensive set of experimental structures. The surprisingly accurate prediction (0.75 A backbone RMSD) for cross-docking of a highly flexible decapeptide, dissimilar to the original bound peptide, as well as docking predictions using homology models for two allotypes with low average backbone RMSDs of less than 1.0 A illustrate the method's effectiveness. Finally, energy terms calculated using the predicted structures were combined with supervised learning on a large data set to classify peptides as either HLA-A*0201 binders or nonbinders. In contrast with sequence-based prediction methods, this model was also able to predict the binding affinity for peptides to a different MHC allotype (H-2K(b)), not used for training, with comparable prediction accuracy.

Role of HLA-B51 and HLA-B52 in susceptibility to pulmonary tuberculosis

MHC class I-restricted CD8+ T cells are important for the generation of protective immune responses in MTB infection. CD8+ CTL (cytotoxic T-lymphocyte)-derived IFN-g may be especially important both for cells lacking MHC class II molecules, e.g. in the lung and for macrophages where mycobacteria can evade recognition during chronic infection by sequestering their antigens away from sensitized CD4+ T cells. This study was designed to detect any association of MHC class I (HLA-B) molecules with pulmonary tuberculosis. A total of 75 individuals, comprising of 33 patients with pulmonary tuberculosis; 12 HIV patients who developed tuberculosis and 30 healthy controls, were included in the study. They were typed for HLA-B by the PCR-SSP method. The results of only HLA-B alleles, which are highly significant, are presented here. The number of healthy individuals with HLA-B52 was significantly high when compared to the patient groups (healthy versus TB: 21.2% versus 0.0%, OR=0.0, P<0.0001, P(c)=0.003; healthy versus HIV-TB: 21.2% versus 16.7%; OR=0.74; P<0.001; P(c)=0.003). In contrast, the number of patients, both TB- and HIV-TB-positive, with HLA-B51 was significantly high when compared to the healthy group of individuals (TB versus healthy: 36.7% versus 3%; OR=18.53; P<0.0001; P(c)=0.001; HIV-TB versus healthy: 41.7% versus 3%; OR=22.86; P<0.0001; P(c)=0.001). Only one healthy control was positive to HLA-B51; however this individual also had HLA-B52. The results of this study suggest that HLA-B52(5) has a negative, i.e. a protective association and HLA-B51(5) has a positive (susceptible) association, for pulmonary tuberculosis. Studies on HLA-B51 and HLA-B52 in a larger population to assess their role in tuberculosis may be useful for TB-vaccination strategies, since HLA profiles are likely to be related to vaccine efficacy.

Structural prediction of peptides bound to MHC class I

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values < or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values < or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.

A biochemical and structural analysis of genetic diversity within the HLA-A*11 subtype

The HLA-A*11 subtype includes 17 naturally occurring variants (-A*1101 to -A*1117) distributed among different ethnic groups worldwide. At present, only HLA-A*1101 has been characterized at the molecular, structural, and immunological level. Developing similar knowledge on other HLA-A*11 alleles is highly important for bone marrow and graft transplantation. This is also important to better understand disease linkages within the HLA-A*11 subtype given that HLA-A*11 molecules are associated with resistance to acquisition of HIV-1 infection and various autoimmune diseases. To broaden our understanding of HLA-A*11 molecules, we have determined the impact of natural polymorphism on the peptide-binding properties of several HLA-A*11 molecules: -A*1103, -A*1106, -A*1108, -A*1110, -A*1111, and -A*1114. We used an approach that combines data from thermal stability studies of recombinant, soluble forms of these molecules in complex with HIV-1 peptides, together with a detailed structural analysis of the resulting HLA-A*11 molecule/peptide complexes based on crystal and molecular model structures. Our analysis shows that natural polymorphism within the HLA-A*11 subtype is distributed along the alpha1 and alpha2 helices of the peptide-binding groove, in marked contrast to the pattern of polymorphism in HLA-A*2 and HLA-B*27 subtypes. Natural polymorphism greatly altered the abilities of individual -A*11 molecules to form stable complexes with HIV-1 peptides. In comparison to -A*1101, natural polymorphism altered the peptide-presenting properties of -A*1103, -A*1108, and -A*1114 and has the potential to affect the peptide-selecting properties of -A*1106, -A*1110, and -A*1111 as well. Overall, our findings suggest that HLA-A*11 molecules may stimulate alloreactive CD8+ cytotoxic T-cell responses.

Crystal structures and KIR3DL1 recognition of three immunodominant viral peptides complexed to HLA-B*2705

We have solved the crystal structures of three HLA-B*2705-peptide complexes with the immunodominant viral peptides: EBV EBNA3C 258-266 (RRIYDLIEL), influenza (flu) nucleoprotein NP383-391 (SRYWAIRTR), and HIV gag 264-273 (KRWIILGLNK). Long-term non-progression during HIV infection has been associated with presentation by HLA-B*2705, and T cell recognition, of the highly immunodominant KRWIILGLNK peptide. The tight hydrogen-bonding network observed between the HLA-B*2705 B-pocket and the peptide P2 arginine guanadinium anchor explains why mutation of this residue during HIV infection results in loss of peptide binding, immune escape and progression to AIDS. Prominent, solvent-exposed structures within these peptides may participate in generating T cell responses to these immunodominant epitopes. In the HLA-B*2705 complex with flu NP383-391, the amino acid side chains of residues 4, 7 and 8 are solvent-exposed whilst in the HIV decamer, the main-chain bulges into the solvent around P7. Thus, HLA-B*2705 presents viral peptides in a range of conformations. Tetrameric complexes of HLA-B*2705 with the HIV and flu but not EBV peptides bound strongly to the killer-Ig-like receptor (KIR)3DL1. Substitution of EBV P8 glutamate to threonine allowed recognition by KIR3DL1. In the HLA-B*2705-EBV structure the P8 glutamate side chain is solvent-exposed and may inhibit KIR3DL1 binding through electrostatic forces.