Magnitude of structural changes of the T-cell receptor binding regions determine the strength of T-cell antagonism: molecular dynamics simulations of HLA-DR4 (DRB1*0405) complexed with analogue peptide

In our model system, we generated T cell clones specific for the HLA-DR4 (DRB1*0405)-index peptide (YWALEAAAD) complex. Based on response patterns of the T cell clones, analogue peptides containing single amino acid substitutions of the index peptide were classified into three types, agonists, antagonists or null peptides (non-agonistic and non-antagonistic peptides). Subtle structural changes induced by the antagonists in the T-cell receptor (TCR) binding regions have already been explained using the root mean square (r.m.s.) deviations from the DR4-index peptide complex in the molecular dynamics (MD) trajectory. In this work, we performed additional MD simulations at 300 K with explicit solvent molecules to reveal the structural character of the HLA-DR4 complexed with the analogue peptides. We examined the r.m.s. deviations of the TCR-binding sites and the exposed areas of the bound peptides. Remarkable differences of the r.m.s. deviations among the DR4-antagonist complexes, together with our previous data, suggest that the magnitude of structural changes of TCR-binding regions would determine the strength of TCR antagonism. The simulations also indicate that TCR could discriminate null peptides from other ligands mainly through the changes of exposed side chains of the bound peptide, rather than the conformational changes of TCR-binding surfaces on HLA molecule.     

Modeling alternative binding registers of a minimal immunogenic peptide on two class II major histocompatibility complex (MHC II) molecules predicts polarized T-cell receptor (TCR) contact positions.

Several major histocompatibility complex class II (MHC II) complexes with known minimal immunogenic peptides have now been solved by X-ray crystallography. Specificity pockets within the MHC II binding groove provide distinct peptide contacts that influence peptide conformation and define the binding register within different allelic MHC II molecules. Altering peptide ligands with respect to the residues that contact the T-cell receptor (TCR) can drastically change the nature of the ensuing immune response. Here, we provide an example of how MHC II (I-A) molecules may indirectly effect TCR contacts with a peptide and drive functionally distinct immune responses. We modeled the same immunogenic 12-amino acid peptide into the binding grooves of two allelic MHC II molecules linked to distinct cytokine responses against the peptide. Surprisingly, the favored conformation of the peptide in each molecule was distinct with respect to the exposure of the N- or C-terminus of the peptide above the MHC II binding groove. T-cell clones derived from each allelic MHC II genotype were found to be allele-restricted with respect to the recognition of these N- vs. C-terminal residues on the bound peptide. Taken together, these data suggest that MHC II alleles may influence T-cell functions by restricting TCR access to specific residues of the I-A-bound peptide. Thus, these data are of significance to diseases that display genetic linkage to specific MHC II alleles, e.g. type 1 diabetes and rheumatoid arthritis.     

Structures of MART-126/27-35 Peptide/HLA-A2 complexes reveal a remarkable disconnect between antigen structural homology and T cell recognition

Small structural changes in peptides presented by major histocompatibility complex (MHC) molecules often result in large changes in immunogenicity, supporting the notion that T cell receptors are exquisitely sensitive to antigen structure. Yet there are striking examples of TCR recognition of structurally dissimilar ligands. The resulting unpredictability of how T cells will respond to different or modified antigens impacts both our understanding of the physical bases for TCR specificity as well as efforts to engineer peptides for immunomodulation. In cancer immunotherapy, epitopes and variants derived from the MART-1/Melan-A protein are widely used as clinical vaccines. Two overlapping epitopes spanning amino acid residues 26 through 35 are of particular interest: numerous clinical studies have been performed using variants of the MART-1 26-35 decamer, although only the 27-35 nonamer has been found on the surface of targeted melanoma cells. Here, we show that the 26-35 and 27-35 peptides adopt strikingly different conformations when bound to HLA-A2. Nevertheless, clonally distinct MART-1(26/27-35)-reactive T cells show broad cross-reactivity towards these ligands. Simultaneously, however, many of the cross-reactive T cells remain unable to recognize anchor-modified variants with very subtle structural differences. These dichotomous observations challenge our thinking about how structural information on unligated peptide/MHC complexes should be best used when addressing questions of TCR specificity. Our findings also indicate that caution is warranted in the design of immunotherapeutics based on the MART-1 26/27-35 epitopes, as neither cross-reactivity nor selectivity is predictable based on the analysis of the structures alone.     

CTL recognition of a bulged viral peptide involves biased TCR selection

MHC class I molecules generally present peptides of 8-10 aa long, forming an extended coil in the HLA cleft. Although longer peptides can also bind to class I molecules, they tend to bulge from the cleft and it is not known whether the TCR repertoire has sufficient plasticity to recognize these determinants during the antiviral CTL response. In this study, we show that unrelated individuals infected with EBV generate a significant CTL response directed toward an HLA-B*3501-restricted, 11-mer epitope from the BZLF1 Ag. The 11-mer determinant adopts a highly bulged conformation with seven of the peptide side chains being solvent-exposed and available for TCR interaction. Such a complex potentially creates a structural challenge for TCR corecognition of both HLA-B*3501 and the peptide Ag. Surprisingly, unrelated B*3501 donors recognizing the 11-mer use identical or closely related alphabeta TCR sequences that share particular CDR3 motifs. Within the small number of dominant CTL clonotypes observed, each has discrete fine specificity for the exposed side chain residues of the peptide. The data show that bulged viral peptides are indeed immunogenic but suggest that the highly constrained TCR repertoire reflects a limit to TCR diversity when responding to some unusual MHC peptide ligands.     

I-Ad-binding peptides derived from unrelated protein antigens share a common structural motif

Purified Ia molecules can specifically bind many unrelated peptide Ag, and such binding appears to be a necessary, albeit not sufficient, prerequisite for the immunogenicity of the proteins from which such peptides are derived. We have recently analyzed the affect of single amino acid substitutions on the I-Ad binding of the immunogenic peptide OVA 323-339. The results obtained demonstrated the very permissive nature of Ag-Ia interaction. We also showed that unrelated peptides that are good I-Ad binders share a common structural motif and speculated that recognition of such motifs could represent a mechanism to achieve a very permissive type of interaction that yet retained some degree of specificity. In the present set of experiments we analyzed the I-Ad binding pattern of a series of overlapping peptides derived from sperm whale myoglobin (residues 102-125) and influenza hemagglutinin (residues 121-146) to determine whether the peptide regions predicted on the basis of structural similarity to be involved in I-Ad binding were in fact involved. In both cases, the I-Ad-interacting determinants were found to contain the sequence motif postulated to be important for I-Ad binding. These data support the hypothesis that I-Ad molecules recognize a large library of Ag by virtue of common structural motifs present in peptides derived from phylogenetically unrelated proteins.     

MHC allele-specific molecular features determine peptide/HLA-A2 conformations that are recognized by HLA-A2-restricted T cell receptors

The structures of alphabeta TCRs bound to complexes of class I MHC molecules and peptide show that the TCRs make multiple contacts with the alpha1 and alpha2 helixes of the MHC. Previously we have shown that the A6 TCR in complex with the HLA-A2/Tax peptide has 15 contact sites on HLA-A2. Single amino acid mutagenesis of these contact sites demonstrated that mutation of only three amino acids clustered on the alpha1 helix (R65, K66, A69) disrupted recognition by the A6 TCR. In the present study we have asked whether TCRs that recognize four other peptides presented by HLA-A2 interact with the MHC in identical, similar, or different patterns as the A6 TCR. Mutants K66A and Q155A had the highest frequency of negative effects on lysis. A subset of peptide-specific CTL also selectively recognized mutants K66A or Q155A in the absence of exogenous cognate peptides, indicating that these mutations affected the presentation of endogenous peptide/HLA-A2 complexes. These findings suggest that most HLA-A2-restricted TCRs recognize surfaces on the HLA-A2/peptide complex that are dependent upon the side chains of K66 and Q155 in the central portion of the peptide binding groove. Crystallographic structures of several peptide/HLA-A2 structures have shown that the side chains of these critical amino acids that make contact with the A6 TCR also contact the bound peptide. Collectively, our results indicate that the generalized effects of changes at these critical amino acids are probably due to the fact that they can be directly contacted by TCRs as well as influence the binding and presentation of the bound peptides.     

Solid-phase epitope recovery: a high throughput method for antigen identification and epitope optimization

Self tolerance to MHC class I-restricted nonmutated self Ags is a significant hurdle to effective cancer immunotherapy. Compelling evidence is emerging that altered peptide ligands can be far more immunogenic than their corresponding native epitopes; however, there is no way to reliably predict which modifications will lead to enhanced native epitope-specific immune responses. We reasoned that this limitation could be overcome by devising an empirical screen in which the nearly complete combinatorial spectrum of peptides of optimal length can be rapidly assayed for reactivity with a MHC class I-restricted cytotoxic T cell clone. This method, solid-phase epitope recovery, quantitatively ranks all reactive peptides in the library and allows selection of altered peptide ligands having desirable immunogenic properties of interest. In contrast to rationally designed MHC anchor-modified peptides, peptides identified by the present method are highly substituted in predicted TCR contact residues and can reliably activate and expand effector cell populations in vitro which lyse target cells presenting the wild-type epitope. We demonstrate that solid-phase epitope recovery peptides corresponding to a poorly immunogenic epitope of the melanoma Ag, gp100, can reliably induce wild-type peptide-specific CTL using normal donor T cells in vitro. Furthermore, these peptides can complement one another to induce these responses in an overwhelming majority of normal individuals in vitro. These data provide a rationale for the design of superior vaccines comprising a mixture of structurally diverse yet functionally convergent peptides.     

Stability of peptide-HLA-I complexes and tapasin folding facilitation--tools to define immunogenic peptides

Only a small fraction of the peptides generated inside the cell end up being presented by HLA-I on the cell surface. High stability of peptide–HLA-I complexes and a low HLA-I tapasin-facilitation have been proposed to predict immunogenicity. We here set out to investigate if these parameters correlated and defined immunogenic peptides. Both peptide–HLA–B¹08:01 and peptide–HLA–A¹02:01 complexes showed small differences in tapasin-facilitation and larger differences in stability. This suggests that the stability of immunogenic peptide–HLA-I complexes vary above an HLA-I allomorph dependent lower limit (e.g. >2 h for HLA–A¹02:01), immunogenicity predicted by tapasin-facilitation may be defined by an equally allomorph unique upper value (e.g. tapasin-facilitation <1.5 for HLA–A¹02:01), and variation above the stability-threshold does not directly reflect a variation in tapasin-facilitation.     

POPI: predicting immunogenicity of MHC class I binding peptides by mining informative physicochemical properties

MOTIVATION: Both modeling of antigen-processing pathway including major histocompatibility complex (MHC) binding and immunogenicity prediction of those MHC-binding peptides are essential to develop a computer-aided system of peptide-based vaccine design that is one goal of immunoinformatics. Numerous studies have dealt with modeling the immunogenic pathway but not the intractable problem of immunogenicity prediction due to complex effects of many intrinsic and extrinsic factors. Moderate affinity of the MHC-peptide complex is essential to induce immune responses, but the relationship between the affinity and peptide immunogenicity is too weak to use for predicting immunogenicity. This study focuses on mining informative physicochemical properties from known experimental immunogenicity data to understand immune responses and predict immunogenicity of MHC-binding peptides accurately. RESULTS: This study proposes a computational method to mine a feature set of informative physicochemical properties from MHC class I binding peptides to design a support vector machine (SVM) based system (named POPI) for the prediction of peptide immunogenicity. High performance of POPI arises mainly from an inheritable bi-objective genetic algorithm, which aims to automatically determine the best number m out of 531 physicochemical properties, identify these m properties and tune SVM parameters simultaneously. The dataset consisting of 428 human MHC class I binding peptides belonging to four classes of immunogenicity was established from MHCPEP, a database of MHC-binding peptides (Brusic et al., 1998). POPI, utilizing the m = 23 selected properties, performs well with the accuracy of 64.72% using leave-one-out cross-validation, compared with two sequence alignment-based prediction methods ALIGN (54.91%) and PSI-BLAST (53.23%). POPI is the first computational system for prediction of peptide immunogenicity based on physicochemical properties. AVAILABILITY: A web server for prediction of peptide immunogenicity (POPI) and the used dataset of MHC class I binding peptides (PEPMHCI) are available at http://iclab.life.nctu.edu.tw/POPI     

Peptide recognition by two HLA-A2/Tax11-19-specific T cell clones in relationship to their MHC/peptide/TCR crystal structures

The crystal structures of two human TCRs specific for a HTLV-I Tax peptide bound to HLA-A2 were recently determined, for the first time allowing a functional comparison of TCRs for which the MHC/peptide/TCR structures are known. Extensive amino acid substitutions show that the native Tax residues are optimal at each peptide position. A prominent feature of the TCR contact surface is a deep pocket that accommodates a tyrosine at position 5 of the peptide. For one of these TCRs, this pocket is highly specific for aromatic residues. In the other TCR structure, this pocket is larger, allowing many different residues to be accommodated. The CTL clones also show major differences in the specificity for several other peptide residues, including side chains that are not directly contacted by the TCR. Despite the specificity of these clones, peptides that are distinct at five or six positions from Tax11-19 induce CTL activity, indicating that substantial changes of the peptide surface are tolerated. Human peptides with limited sequence homology to Tax11-19 represent partial TCR agonists for these CTL clones. The distinct functional properties of these CTL clones highlight structural features that determine TCR specificity and cross-reactivity for MHC-bound peptides.     

Molecular dynamics simulations of helical antimicrobial peptides in SDS micelles: what do point mutations achieve?

We report long time scale simulations of the 18-residue helical antimicrobial peptide ovispirin-1 and its analogs novispirin-G10 and novispirin-T7 in SDS micelles. The SDS micelle serves as an economical and effective model for a cellular membrane. Ovispirin, which is initially placed along a micelle diameter, diffuses out to the water-SDS interface and stabilizes to an interface-bound steady state in 16.35 ns of simulation. The final conformation, orientation, and the structure of ovispirin are in good agreement with the experimentally observed properties of the peptide in presence of lipid bilayers. The simulation succeeds in capturing subtle differences of the membrane-bound peptide structure as predicted by solid state NMR. The novispirins also undergo identical diffusion patterns and similar final conformations. Although the final interface-bound states are similar, the simulations illuminate the structural and binding properties of the mutant peptides which make them less toxic compared to ovispirin. Based on previous data and the current simulations, we propose that introduction of a bend/hinge at the center of helical antimicrobial peptides (containing a specific C-terminal motif), without disrupting the helicity of the peptides might attenuate host-cell toxicity as well as improve membrane binding properties to bacterial cellular envelopes.     

BiPPred: Combined sequence‐ and structure‐based prediction of peptide binding to the Hsp70 chaperone BiP

Substrate binding to Hsp70 chaperones is involved in many biological processes, and the identification of potential substrates is important for a comprehensive understanding of these events. We present a multi‐scale pipeline for an accurate, yet efficient prediction of peptides binding to the Hsp70 chaperone BiP by combining sequence‐based prediction with molecular docking and MMPBSA calculations. First, we measured the binding of 15mer peptides from known substrate proteins of BiP by peptide array (PA) experiments and performed an accuracy assessment of the PA data by fluorescence anisotropy studies. Several sequence‐based prediction models were fitted using this and other peptide binding data. A structure‐based position‐specific scoring matrix (SB‐PSSM) derived solely from structural modeling data forms the core of all models. The matrix elements are based on a combination of binding energy estimations, molecular dynamics simulations, and analysis of the BiP binding site, which led to new insights into the peptide binding specificities of the chaperone. Using this SB‐PSSM, peptide binders could be predicted with high selectivity even without training of the model on experimental data. Additional training further increased the prediction accuracies. Subsequent molecular docking (DynaDock) and MMGBSA/MMPBSA‐based binding affinity estimations for predicted binders allowed the identification of the correct binding mode of the peptides as well as the calculation of nearly quantitative binding affinities. The general concept behind the developed multi‐scale pipeline can readily be applied to other protein‐peptide complexes with linearly bound peptides, for which sufficient experimental binding data for the training of classical sequence‐based prediction models is not available. Proteins 2016; 84:1390–1407. © 2016 Wiley Periodicals, Inc.     

NMR structure of Plasmodium falciparum malaria peptide correlates with protective immunity

Apical membrane antigen-1 is an integral Plasmodium falciparum malaria parasite membrane protein. High activity binding peptides (HABPs) to human red blood cells (RBCs) have been identified in this protein. One of them (peptide 4313), for which critical binding residues have already been defined, is conserved and nonimmunogenic. Its critical binding residues were changed for amino acids having similar mass but different charge to change such immunological properties; these changes generated peptide analogues. Some of these peptide analogues became immunogenic and protective in Aotus monkeys.Three-dimensional models of peptide 4313 and three analogues having different immune characteristics, were calculated from nuclear magnetic resonance (NMR) experiments with distance geometry and restrained molecular dynamic methods. All peptides contained a beta-turn structure spanning amino acids 7 to 10, except randomly structured 4313. When analysing dihedral angle phi and psi values, distorted type III or III' turns were identified in the protective and/or immunogenic peptides, whilst classical type III turns were found for the nonimmunogenic nonprotective peptides. This data shows that some structural modifications may lead to induction of immunogenicity and/or protection, suggesting a new way to develop multicomponent, subunit-based malarial vaccines.     

Structural features of a chimeric peptide inducing cytotoxic T lymphocyte responses in saline.

Little information is available correlating the structural properties of peptides with their immunogenicity in terms of responses via cytotoxic T lymphocytes (CTLs). The TT-NP6 chimeric peptide, consisting of two copies of a promiscuous T-helper epitope (T: residues 288-302 from the fusion protein of the measles virus) linked to the NP6 T-cytotoxic epitope (NP6: residues 52-60 from the nucleoprotein of measles virus) was able to induce virus-specific CTL responses in the absence of any adjuvant and hydrophobic component. The present work was undertaken to gain insight into structural features of the TT-NP6 peptide that may be important in optimizing the CTL immunogenicity of the peptide. Circular dichroism data, obtained in a buffer of physiological ionic strength and pH, strongly suggest a self-associated state for the peptide, which was confirmed by a sedimentation velocity experiment. However, helix association is accompanied by loss of overall helical content. Thermal-dependence studies show that the unfolding of self-associated alpha-helices is significantly more pronounced than the unfolding of isolated alpha-helices. Circular dichroism data, together with tryptic limited proteolysis, suggest the presence of a charged amino acid within the hydrophobic core. This study should provide a basis for engineering more effective immunogenic peptides against the measles virus by increasing the stability of the TT-NP6 peptide.     

Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.     

Probing the conformational and dynamical effects of O-glycosylation within the immunodominant region of a MUC1 peptide tumor antigen.

MUC1 mucin is a large transmembrane glycoprotein, the extracellular domain of which is formed by a repeating 20 amino acid sequence, GVTSAPDTRPAPGSTAPPAH. In normal breast epithelial cells, the extracellular domain is densely covered with highly branched complex carbohydrate structures. However, in neoplastic breast tissue, the extracellular domain is under-glycosylated, resulting in the exposure of a highly immunogenic core peptide epitope (PDTRP in bold above), as well as in the exposure of normally cryptic core Tn (GalNAc), STn (sialyl alpha2-6 GalNAc) and TF (Gal beta1-3 GalNAc) carbohydrates. Here, we report the results of 1H NMR structural studies, natural abundance 13C NMR relaxation measurements and distance-restrained MD simulations designed to probe the structural and dynamical effects of Tn-glycosylation within the PDTRP core peptide epitope. Two synthetic peptides were studied: a nine-residue MUC1 peptide of the sequence, Thr1-Ser2-Ala3-Pro4-Asp5-Thr6-Arg7-Pro8-Ala9, and a Tn-glycosylated version of this peptide, Thr1-Ser2-Ala3-Pro4-Asp5-Thr6(alphaGalNAc)-Arg7-Pro8-Ala9. The results of these studies show that a type I beta-turn conformation is adopted by residues PDTR within the PDTRP region of the unglycosylated MUC1 sequence. The existence of a similar beta-turn within the PDTRP core peptide epitope of the under-glycosylated cancer-associated MUC1 mucin protein might explain the immunodominance of this region in vivo, as the presence of defined secondary structure within peptide epitope regions has been correlated with increased immunogenicity in other systems. Our results have also shown that Tn glycosylation at the central threonine within the PDTRP core epitope region shifts the conformational equilibrium away from the type I beta-turn conformation and toward a more rigid and extended state. The significance of these results are discussed in relation to the possible roles that peptide epitope secondary structure and glycosylation state may play in MUC1 tumor immunogenicity.     

The 3A6-TCR/superagonist/HLA-DR2a complex shows similar interface and reduced flexibility compared to the complex with self-peptide

T-cell receptor (TCR) recognition of the myelin basic protein (MBP) peptide presented by major histocompatibility complex (MHC) protein HLA-DR2a, one of the MHC class II alleles associated with multiple sclerosis, is highly variable. Interactions in the trimolecular complex between the TCR of the MBP83-99-specific T cell clone 3A6 with the MBP-peptide/HLA-DR2a (abbreviated TCR/pMHC) lead to substantially different proliferative responses when comparing the wild-type decapeptide MBP90-99 and a superagonist peptide, which differs mainly in the residues that point toward the TCR. Here, we investigate the influence of the peptide sequence on the interface and intrinsic plasticity of the TCR/pMHC trimolecular and pMHC bimolecular complexes by molecular dynamics simulations. The intermolecular contacts at the TCR/pMHC interface are similar for the complexes with the superagonist and the MBP self-peptide. The orientation angle between TCR and pMHC fluctuates less in the complex with the superagonist peptide. Thus, the higher structural stability of the TCR/pMHC tripartite complex with the superagonist peptide, rather than a major difference in binding mode with respect to the self-peptide, seems to be responsible for the stronger proliferative response.     

Hierarchy of simulation models in predicting molecular recognition mechanisms from the binding energy landscapes: structural analysis of the peptide complexes with SH2 domains.

Computer simulations using the simplified energy function and simulated tempering dynamics have accurately determined the native structure of the pYVPML, SVLpYTAVQPNE, and SPGEpYVNIEF peptides in the complexes with SH2 domains. Structural and equilibrium aspects of the peptide binding with SH2 domains have been studied by generating temperature-dependent binding free energy landscapes. Once some native peptide-SH2 domain contacts are constrained, the underlying binding free energy profile has the funnel-like shape that leads to a rapid and consistent acquisition of the native structure. The dominant native topology of the peptide-SH2 domain complexes represents an extended peptide conformation with strong specific interactions in the phosphotyrosine pocket and hydrophobic interactions of the peptide residues C-terminal to the pTyr group. The topological features of the peptide-protein interface are primarily determined by the thermodynamically stable phosphotyrosyl group. A diversity of structurally different binding orientations has been observed for the amino-terminal residues to the phosphotyrosine. The dominant native topology for the peptide residues carboxy-terminal to the phosphotyrosine is tolerant to flexibility in this region of the peptide-SH2 domain interface observed in equilibrium simulations. The energy landscape analysis has revealed a broad, entropically favorable topology of the native binding mode for the bound peptides, which is robust to structural perturbations. This could provide an additional positive mechanism underlying tolerance of the SH2 domains to hydrophobic conservative substitutions in the peptide specificity region.