In silico design of high-affinity antigenic peptides for HLA-B44.

Cancer cell-killing by CD8+ T cells demands effective tumor antigen presentation by human leukocyte antigen class I (HLA-I) molecules. Screening and designing highly immunogenic neoantigens require quantitative computations to reliably predict HLA-peptide binding affinities. Here, with all-atom molecular dynamics (MD) simulations and free energy perturbation (FEP) methods, we design a collection of antigenic peptide candidates through in silico mutagenesis studies on immunogenic neoantigens, yielding enhanced binding affinities to HLA-B*44:02. In-depth structural dissection shows that introducing positively charged residues such as arginine to position 6 or lysine to position 7 of the candidates triggers conformational shifts in both peptides and the antigen-binding groove of the HLA, following the "induced-fit" mechanism. Enhancement in binding affinities compared to the wild-type was found in three out of five mutated candidates. The HLA pocket, capable of accommodating positively charged residues in positions from 5 to 7, is designated as the "dynamic pocket". Taken together, we showcase an effective structure-based binding affinity optimization framework for antigenic peptides of HLA-B*44:02 and underscore the importance of dynamic nature of the antigen-binding groove in concert with the anchoring motifs. This work provides structural insights for rational design of favorable HLA-peptide bindings and future developments in neoantigen-based therapeutics.

Computational investigation of peptide binding stabilities of HLA-B*27 and HLA-B*44 alleles.

Major Histocompatibility Complex (MHC) is a cell surface glycoprotein that binds to foreign antigens and presents them to T lymphocyte cells on the surface of Antigen Presenting Cells (APCs) for appropriate immune recognition. Recently, studies focusing on peptide-based vaccine design have allowed a better understanding of peptide immunogenicity mechanisms, which is defined as the ability of a peptide to stimulate CTL-mediated immune response. Peptide immunogenicity is also known to be related to the stability of peptide-loaded MHC (pMHC) complex. In this study, ENCoM server was used for structure-based estimation of the impact of single point mutations on pMHC complex stabilities. For this purpose, two human MHC molecules from the HLA-B*27 group (HLA-B*27:05 and HLA-B*27:09) in complex with four different peptides (GRFAAAIAK, RRKWRRWHL, RRRWRRLTV and IRAAPPPLF) and three HLA-B*44 molecules (HLA-B*44:02, HLA-B*44:03 and HLA-B*44:05) in complex with two different peptides (EEYLQAFTY and EEYLKAWTF) were analyzed. We found that the stability of pMHC complexes is dependent on both peptide sequence and MHC allele. Furthermore, we demonstrate that allele-specific peptide-binding preferences can be accurately revealed using structure-based computational methods predicting the effect of mutations on protein stability.

Selector function of MHC I molecules is determined by protein plasticity.

The selection of peptides for presentation at the surface of most nucleated cells by major histocompatibility complex class I molecules (MHC I) is crucial to the immune response in vertebrates. However, the mechanisms of the rapid selection of high affinity peptides by MHC I from amongst thousands of mostly low affinity peptides are not well understood. We developed computational systems models encoding distinct mechanistic hypotheses for two molecules, HLA-B*44:02 (B*4402) and HLA-B*44:05 (B*4405), which differ by a single residue yet lie at opposite ends of the spectrum in their intrinsic ability to select high affinity peptides. We used in vivo biochemical data to infer that a conformational intermediate of MHC I is significant for peptide selection. We used molecular dynamics simulations to show that peptide selector function correlates with protein plasticity, and confirmed this experimentally by altering the plasticity of MHC I with a single point mutation, which altered in vivo selector function in a predictable way. Finally, we investigated the mechanisms by which the co-factor tapasin influences MHC I plasticity. We propose that tapasin modulates MHC I plasticity by dynamically coupling the peptide binding region and α3 domain of MHC I allosterically, resulting in enhanced peptide selector function.

Mechanistic Basis for Epitope Proofreading in the Peptide-Loading Complex.

The peptide-loading complex plays a pivotal role in Ag processing and is thus central to the efficient immune recognition of virally and malignantly transformed cells. The underlying mechanism by which MHC class I (MHC I) molecules sample immunodominant peptide epitopes, however, remains poorly understood. In this article, we delineate the interaction between tapasin (Tsn) and MHC I molecules. We followed the process of peptide editing in real time after ultra-fast photoconversion to pseudoempty MHC I molecules. Tsn discriminates between MHC I loaded with optimal and MHC I bound to suboptimal cargo. This differential interaction is key to understanding the kinetics of epitope proofreading. To elucidate the underlying mechanism at the atomic level, we modeled the Tsn/MHC I complex using all-atom molecular dynamics simulations. We present a catalytic working cycle, in which Tsn binds to MHC I with suboptimal cargo and thereby adjusts the energy landscape in favor of MHC I complexes with immunodominant epitopes.

Coupling between side chain interactions and binding pocket flexibility in HLA-B*44:02 molecules investigated by molecular dynamics simulations.

MHC class I molecules present antigenic peptides to cytotoxic T-cells at the cell surface. Peptide loading of class I molecules in the endoplasmatic reticulum can involve interaction with the tapasin chaperone protein. The human class I allotype HLA-B*44:02 with an Asp at position 116 at the floor of the F pocket (which binds the peptide C-terminal residues) depends on tapasin for efficient peptide loading. However, HLA-B*44:05 (identical to B*44:02 except for tyrosine 116) can efficiently load peptides in the absence of tapasin. Both allotypes adopt very similar structures in the presence of the same peptide. Molecular dynamics simulations indicate a significantly higher conformational flexibility of the F pocket in the absence of a peptide for B*44:02 compared to B*44:05. Free energy simulations to open the F pocket indicate a molecular side chain switch mechanism that underlies the global opening motion. This side chain switch involves the rearrangement of salt bridges and hydrogen bonding of the basic arginine 97 with three acidic aspartate residues 114, 116 and 156 near the F pocket. A replica exchange simulation to specifically accelerate side chain motions demonstrates that the same side chain rearrangements induce global opening motions of the F pocket. In case of B*44:05 the free energy barrier for F pocket opening was significantly higher compared to B*44:02 and no associated side chain rearrangement was observed. Such coupling of local side chain rearrangements with global conformational changes might be the basis for allosteric changes in other class I allotypes as well as for allosteric changes in other proteins.

Residue 81 confers a restricted C-terminal peptide binding motif in HLA-B*44:09.

Knowledge about the magnitude of individual polymorphism is a critical part in understanding the complexity of comprehensive mismatching. HLA-B*44:09 differs from the highly frequent HLA-B*44:02 allele by amino acid exchanges at residues 77, 80, 81, 82 and 83. We aimed to identify the magnitude of these mismatches on the features of HLA-B*44:09 bound peptides since residues 77, 80 and 81 comprise part of the F pocket which determines sequence specificity at the pΩ position of the peptide. Using soluble HLA technology we determined >200 individual (nonduplicate) self-peptides from HLA-B*44:09 and compared their features with that of the published peptide features of HLA-B*44:02. Both alleles illustrate an anchor motif of E at p2. In contrast to the C-terminal peptide binding motif of B*44:02 (W, F, Y or L), B*44:09-derived peptides are restricted predominantly to L or F. The source of peptides for both alleles is identical (LCL 721.221 cells) allowing us to identify 23 shared peptides. The majority of these peptides however contained the restricted B*44:09 anchor motif of F or L at the pΩ position. Molecular modelling based on the B*44:02 structure highlights that the differences of the C-terminal peptide anchor between both alleles can be explained primarily by the B*44:02(81Ala) > B*44:09(81Leu) polymorphism which restricts the size of the amino acid that can be accommodated in the F pocket of B*44:09. These results highlight that every amino acid substitution has an impact of certain magnitude on the alleles function and demonstrate how surrounding residues orchestrate peptide specificity.

Specific binding of HLA-B44 to human macrophage MHC receptor 1 on monocytes.

Allograft (H-2D(d)K(d))-induced macrophages (AIM) in C57BL/6 (H-2D(b)K(b)) mice exhibit major histocompatibility complex (MHC) haplotype-specific killing of allografts in a macrophage MHC receptor 1 (MMR1; for H-2D(d))- and MMR2 (for H-2K(d))-dependent manner. Recently, we showed HLA-B62 to be a ligand for the human homologue of mouse MMR2. In the present study, we isolated a cDNA encoding the human homologue of mouse MMR1 and found HLA-B44 to be the sole ligand specific for the human MMR1 by using beads that had been conjugated with 80 kinds of HLA proteins. Flow cytometric analyses revealed that HLA-B44-conjugated beads are specifically bound to HEK293T cells expressing human MMR1, that HLA-B44 tetramers are bound to the human MMR1-transfected HEK293T cells with a dissociation constant of 3.0×10(-9) M, and that the interaction was completely inhibited by the addition of R15 monoclonal antibody specific for mouse MMR1. The MMR1 cDNA (1537-bp) encoded a 473-amino acid polypeptide and was expressed at least in part in the brain and peripheral blood mononuclear cells (PBMCs) or monocytes, but not in granulocytes or lymphocytes. PBMCs from 7 non-H-2D(d) (non-self), but none from 5 H-2D(d) (self), in-bred mice expressed mouse MMR1 specific for H-2D(d). In contrast, PBMCs from none of the 16 human volunteers expressed HLA-B44; whereas those from only 3 of these 16 volunteers expressed human MMR1. These results reveal that human MMR1 on monocytes is a novel receptor specific for HLA-B44.

HLA alleles and amino-acid signatures of the peptide-binding pockets of HLA molecules in vitiligo.

Vitiligo is a depigmenting disorder of the skin that is characterized by the loss of functional melanocytes from the lesional sites. Although the exact etiology is not understood, autoimmunity is thought to be a crucial deterministic factor. A recurring theme of several autoimmune disorders is the aberrant presentation of self-antigens to the immune system, which triggers downstream perturbations. Here we examine the role of alleles of HLA class I and class II loci to delineate vitiligo manifestation in two distinct populations. Our studies have identified three specific alleles, HLA-A*33:01, HLA-B*44:03, and HLA-DRB1*07:01, to be significantly increased in vitiligo patients as compared with controls in both the initial study on North Indians (N=1,404) and the replication study in Gujarat (N=355) cases, establishing their positive association with vitiligo. Both generalized and localized vitiligo have the same predisposing major histocompatibility complex alleles, i.e., B*44:03 and DRB1*07:01, in both the populations studied, beside the differences in the frequencies of other alleles, suggesting that localized vitiligo too may be an autoimmune disorder. Significant differences in the amino-acid signatures of the peptide-binding pockets of HLA-A and HLA-B α-chain and HLA-DR ß-chain were observed between vitiligo patients and unaffected controls.

Position 156 influences the peptide repertoire and tapasin dependency of human leukocyte antigen B*44 allotypes.

BACKGROUND: Polymorphic differences between donor and recipient human leukocyte antigen class I molecules can result in graft-versus-host disease due to distinct peptide presentation. As part of the peptide-loading complex, tapasin plays an important role in selecting peptides from the pool of potential ligands. Class I polymorphisms can significantly alter the tapasin-mediated interaction with the peptide-loading complex and although most class I allotypes are highly dependent upon tapasin, some are able to load peptides independently of tapasin. Several human leukocyte antigen B*44 allotypes differ exclusively at position 156 (B*44:02(156Asp), 44:03(156Leu), 44:28(156Arg), 44:35(156Glu)). From these alleles, only the high tapasin-dependency of human leukocyte antigen B*44:02 has been reported. DESIGN AND METHODS: We investigated the influence of position 156 polymorphisms on both the requirement of tapasin for efficient surface expression of each allotype and their peptide features. Genes encoding human leukocyte antigen B*44 variants bearing all possible substitutions at position 156 were lentivirally transduced into human leukocyte antigen class I-negative LCL 721.221 cells and the tapasin-deficient cell line LCL 721.220. RESULTS: Exclusively human leukocyte antigen B*44:28(156Arg) was expressed on the surface of tapasin-deficient cells, suggesting that the remaining B*44/156 variants are highly tapasin-dependent. Our computational analysis suggests that the tapasin-independence of human leukocyte antigen B*44:28(156Arg) is a result of stabilization of the peptide binding region and generation of a more peptide receptive state. Sequencing of peptides eluted from human leukocyte antigen B*44 molecules by liquid chromatography-electrospray ionization-mass spectrometry (LTQ-Orbitrap) demonstrated that both B*44:02 and B*44:28 share the same overall peptide motif and a certain percentage of their individual peptide repertoires in the presence and/or absence of tapasin. CONCLUSIONS: Here we report for the first time the influence of position 156 on the human leukocyte antigen/tapasin association. Additionally, the results of peptide sequencing suggest that tapasin chaperoning is needed to acquire peptides of unusual length.

Mismatches outside exons 2 and 3 do not alter the peptide motif of the allele group B*44:02P.

Sequence variations outside exons 2 and 3 do not appear to affect the function of human leukocyte antigen (HLA) class I alleles. HLA-B*44:02:01:01 and -B*44:27 are considered functionally identical because they differ by a single amino acid substitution of Val > Ala at position 199, which is located in the α3 domain. To validate that HLA-B*44:02:01:01 and -B*44:27 represent functionally identical alleles that might reflect a permissive mismatch in hematopoetic stem cell transplantation (HSCT), we determined their peptide-binding features. B-lymphoblastic cells were lentivirally transduced with B*44:02 and B*44:27 constructs and soluble recombinant molecules were purified by affinity chromatography. Peptides were isolated and sequencing of single peptides was performed using liquid chromatography-electrospray ionization-tandem mass spectrometry (LTQ-Orbitrap) technology. We demonstrate that the peptide motif of B*44:02(199Val) and B*44:27(199Ala) is identical. Both variants feature E at P2 and Y, F, or W at PΩ in their ligands. Most of the identified peptides are 9 to 11 amino acids in length and approximately 20% of these ligands are shared between the alleles. Our results lead to the conclusion that B*44:02:01:01 and B*44:27 might have the same immune function, validating a theory that is now being used in deciding which donors to select in HSCT when there is no identical donor available.