Toward the prediction of class I and II mouse major histocompatibility complex-peptide-binding affinity:in silico bioinformatic step-by-step guide using quantitative structure-activity relationships

Quantitative structure-activity relationship (QSAR) analysis is a cornerstone of modern informatics. Predictive computational models of peptide-major histocompatibility complex (MHC)-binding affinity based on QSAR technology have now become important components of modern computational immunovaccinology. Historically, such approaches have been built around semiqualitative, classification methods, but these are now giving way to quantitative regression methods. We review three methods--a 2D-QSAR additive-partial least squares (PLS) and a 3D-QSAR comparative molecular similarity index analysis (CoMSIA) method--which can identify the sequence dependence of peptide-binding specificity for various class I MHC alleles from the reported binding affinities (IC50) of peptide sets. The third method is an iterative self-consistent (ISC) PLS-based additive method, which is a recently developed extension to the additive method for the affinity prediction of class II peptides. The QSAR methods presented here have established themselves as immunoinformatic techniques complementary to existing methodology, useful in the quantitative prediction of binding affinity: current methods for the in silico identification of T-cell epitopes (which form the basis of many vaccines, diagnostics, and reagents) rely on the accurate computational prediction of peptide-MHC affinity. We have reviewed various human and mouse class I and class II allele models. Studied alleles comprise HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203, HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3101, HLA-A*6801, HLA-A*6802, HLA-B*3501, H2-K(k), H2-K(b), H2-D(b) HLA-DRB1*0101, HLA-DRB1*0401, HLA-DRB1*0701, I-A(b), I-A(d), I-A(k), I-A(S), I-E(d), and I-E(k). In this chapter we show a step-by-step guide into predicting the reliability and the resulting models to represent an advance on existing methods. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made are freely available online at the URL http://www.jenner.ac.uk/MHCPred.

Predicting class I major histocompatibility complex (MHC) binders using multivariate statistics:comparison of discriminant analysis and multiple linear regression

The accurate in silico identification of T-cell epitopes is a critical step in the development of peptide-based vaccines, reagents, and diagnostics. It has a direct impact on the success of subsequent experimental work. Epitopes arise as a consequence of complex proteolytic processing within the cell. Prior to being recognized by T cells, an epitope is presented on the cell surface as a complex with a major histocompatibility complex (MHC) protein. A prerequisite therefore for T-cell recognition is that an epitope is also a good MHC binder. Thus, T-cell epitope prediction overlaps strongly with the prediction of MHC binding. In the present study, we compare discriminant analysis and multiple linear regression as algorithmic engines for the definition of quantitative matrices for binding affinity prediction. We apply these methods to peptides which bind the well-studied human MHC allele HLA-A*0201. A matrix which results from combining results of the two methods proved powerfully predictive under cross-validation. The new matrix was also tested on an external set of 160 binders to HLA-A*0201; it was able to recognize 135 (84%) of them.

Toward prediction of class II mouse major histocompatibility complex peptide binding affinity:in silico bioinformatic evaluation using partial least squares, a robust multivariate statistical technique

The accurate identification of T-cell epitopes remains a principal goal of bioinformatics within immunology. As the immunogenicity of peptide epitopes is dependent on their binding to major histocompatibility complex (MHC) molecules, the prediction of binding affinity is a prerequisite to the reliable prediction of epitopes. The iterative self-consistent (ISC) partial-least-squares (PLS)-based additive method is a recently developed bioinformatic approach for predicting class II peptide−MHC binding affinity. The ISC−PLS method overcomes many of the conceptual difficulties inherent in the prediction of class II peptide−MHC affinity, such as the binding of a mixed population of peptide lengths due to the open-ended class II binding site. The method has applications in both the accurate prediction of class II epitopes and the manipulation of affinity for heteroclitic and competitor peptides. The method is applied here to six class II mouse alleles (I-Ab, I-Ad, I-Ak, I-As, I-Ed, and I-Ek) and included peptides up to 25 amino acids in length. A series of regression equations highlighting the quantitative contributions of individual amino acids at each peptide position was established. The initial model for each allele exhibited only moderate predictivity. Once the set of selected peptide subsequences had converged, the final models exhibited a satisfactory predictive power. Convergence was reached between the 4th and 17th iterations, and the leave-one-out cross-validation statistical terms - q2, SEP, and NC - ranged between 0.732 and 0.925, 0.418 and 0.816, and 1 and 6, respectively. The non-cross-validated statistical terms r2 and SEE ranged between 0.98 and 0.995 and 0.089 and 0.180, respectively. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made freely available online (http://www.jenner.ac.uk/MHCPred).

MHCPred : a server for quantitative prediction of peptide-MHC binding

Accurate T-cell epitope prediction is a principal objective of computational vaccinology. As a service to the immunology and vaccinology communities at large, we have implemented, as a server on the World Wide Web, a partial least squares-base multivariate statistical approach to the quantitative prediction of peptide binding to major histocom-patibility complexes (MHC), the key checkpoint on the antigen presentation pathway within adaptive,cellular immunity. MHCPred implements robust statistical models for both Class I alleles (HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203,HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3301, HLA-A*6801, HLA-A*6802 and HLA-B*3501) and Class II alleles (HLA-DRB*0401, HLA-DRB*0401and HLA-DRB* 0701).

EpiDOCK:a molecular docking-based tool for MHC class II binding prediction

Cellular peptide vaccines contain T-cell epitopes. The main prerequisite for a peptide to act as a T-cell epitope is that it binds to a major histocompatibility complex (MHC) protein. Peptide MHC binder identification is an extremely costly experimental challenge since human MHCs, named human leukocyte antigen, are highly polymorphic and polygenic. Here we present EpiDOCK, the first structure-based server for MHC class II binding prediction. EpiDOCK predicts binding to the 23 most frequent human, MHC class II proteins. It identifies 90% of true binders and 76% of true non-binders, with an overall accuracy of 83%. EpiDOCK is freely accessible at http://epidock.ddg-pharmfac. net. © The Author 2013. Published by Oxford University Press. All rights reserved.

Protein expression of Myt272-3 recombinant clone and in silico prediction of a possible vaccine candidate against Mycobacterium tuberculosis

Purpose: To investigate the expression of Myt272-3 recombinant protein and also to predict a possible protein vaccine candidate against Mycobacterium tuberculosis . Methods: Myt272-3 protein was expressed in pET30a+-Myt272-3 clone. The purity of the protein was determined using Dynabeads® His-Tag Isolation & Pulldown. Protein sequence was analysed in silico using bioinformatics software for the prediction of allergenicity, antigenicity, MHC-I and MHC-II binding, and B-cell epitope binding. Results: The candidate protein was a non-allergen with 15.19 % positive predictive value. It was also predicted to be antigenic, with binding affinity to MHC-I and MHC-II, as well as B-cell epitope binding. Conclusion: The predicted results obtained in this study provide a guide for practical design of a new tuberculosis vaccine.

The effect of the HLA B27 allele on the immune response to acute HCV in HIV infected patients

In mono-infected individuals, the HLA-B27 allele is strongly associated with spontaneous clearance of HCV in association with a strong CD8+ response targeted against a single epitope within the HCV RNA-dependent RNA polymerase (NS5B). We studied variation across the whole HCV genome and T cell responses over time in a rare cohort of HLA-B27+ patients with acute HCV and HIV co-infection, the majority of whom progressed to chronicity. We used next generation sequencing to detect changes within and outwith the immuno-dominant HLA-B27 restricted HCV-specific CD8+ T cell epitope NS5B2841-2849 (ARMILMTHF) during evolving progression of early HCV infection. Within the Acute HCV UK cohort, 10 patients carried the HLA B27 allele. Of these, 3/8 patients (37.5%) with HIV infection and 2/2 (100%) without HIV spontaneously cleared HCV (p=0.44). Sequential samples from nine HLA-B27+ patients (2 with monoinfection and 7 with HIV co-infection) were available for analysis (four spontaneous clearers and five evolving progressors). Mutations identified using NGS were assessed using a replicon genotype 1a system to evaluate viral fitness. Multiple mutations within the HLA-B27 restricted NS5B2841-2849 epitope were associated with progression to chroncity whereas patients who cleared the HCV infection spontaneously had no or only one mutation at this site (p=0.03). A triple NS5B2841-2849 mutant observed during progression to chronicity was associated with restored replication when compared to wild-type virus while single or double mutants were significantly associated with impaired replication (p=0.0495). T cell responses measured in these patients using ELISpot and flow cytometry. HLA-B27+ patients had significantly higher IFN-γ responses than patients who were HLA-B27- (p=0.0014). Those who progressed to chronicity had lower IFN-γ responses than those who cleared HCV (p=0.0011). Mono-infected patients had higher IFN-γ responses compared to co-infected patients (p=0.0015). HIV co-infection is associated with a lower likelihood of spontaneous clearance of HCV in HLA B27+ patients and this is associated with impaired T cell function in this group.

Consequences of immunodominant epitope deletion for minor influenza virus-specific CD8+-T-cell responses

The extent to which CD8+ T cells specific for other antigens expand to compensate for the mutational loss of the prominent DbNP366 and DbPA224 epitopes has been investigated using H1N1 and H3N2 influenza A viruses modified by reverse genetics. Significantly increased numbers of CD8+ KbPB1703+ , CD8+ KbNS2114+, and CD8+ DbPB1-F262+ T cells were found in the spleen and in the inflammatory population recovered by bronchoalveolar lavage from mice that were first given the -NP-PA H1N1 virus intraperitoneally and then challenged intranasally with the homologous H3N2 virus. The effect was less consistent when this prime-boost protocol was reversed. Also, though the quality of the response measured by cytokine staining showed some evidence of modification when these minor CD8+-T-cell populations were forced to play a more prominent part, the effects were relatively small and no consistent pattern emerged. The magnitude of the enhanced clonal expansion following secondary challenge suggested that the prime-boost with the -NP-PA viruses gave a response overall that was little different in magnitude from that following comparable exposure to the unmanipulated viruses. This was indeed shown to be the case when the total response was measured by ELISPOT analysis with virus-infected cells as stimulators. More surprisingly, the same effect was seen following primary challenge, though individual analysis of the CD8+ KbPB1703+ , CD8+ KbNS2114+, and CD8+ DbPB1-F262+ sets gave no indication of compensatory expansion. A possible explanation is that novel, as yet undetected epitopes emerge following primary exposure to the -NP-PA deletion viruses. These findings have implications for both natural infections and vaccines.

Potent T cell response to a class I-binding 13-mer viral epitope and the influence of HLA micropolymorphism in controlling epitope length

The BZLF1 antigen of Epstein-Barr virus includes three overlapping sequences of different lengths that conform to the binding motif of human leukocyte antigen (HLA) B*3501. These 9-mer ((56)LPOGQLTAy(64)), 11-mer ((54)EPLPQGQLTAy(64)), and 13-mer ((52)LPEPLPQGQLTAY(64)) peptides all bound well to B*3501; however, the CTL response in individuals expressing this HILA allele was directed strongly and exclusively towards the 11-mer peptide. In contrast, EBV-exposed donors expressing HLA B*3503 showed no significant CTL response to these peptides because the single amino acid difference between B*3501 and B*3503 within the F pocket inhibited HLA binding by these peptides. The extraordinarily long 13-mer peptide was the target for the CTL response in individuals expressing B*3508, which differs from B*3501 at a single position within the D pocket (B*3501, 156 Leucine; B*3508, 156 Arginine). This minor difference was shown to enhance binding of the 13-mer peptide, presumably through a stabilizing interaction between the negatively charged glutamate at position 3 of the peptide and the positively charged arginine at HLA position 156. The 13-mer epitope defined in this study represents the longest class I-binding viral epitope identified to date as a minimal determinant. Furthermore, the potency of the response indicates that peptides of this length do not present a major structural barrier to CTL recognition.

TCR{alpha} Genes Direct MHC Restriction in the Potent Human T Cell Response to a Class I-Bound Viral Epitope

The underlying generic properties of {alpha}β TCRs that control MHC restriction remain largely unresolved. To investigate MHC restriction, we have examined the CTL response to a viral epitope that binds promiscuously to two human leukocyte Ags (HLAs) that differ by a single amino acid at position 156. Individuals expressing either HLA-B*3501 (156Leucine) or HLA-B*3508 (156Arginine) showed a potent CTL response to the 407HPVGEADYFEY417 epitope from EBV. Interestingly, the response was characterized by highly restricted TCR β-chain usage in both HLA-B*3501+ and HLA-B*3508+ individuals; however, this conserved TRBV9+ β-chain was associated with distinct TCR {alpha}-chains depending upon the HLA-B*35 allele expressed by the virus-exposed host. Functional assays confirmed that TCR {alpha}-chain usage determined the HLA restriction of the CTLs. Structural studies revealed significant differences in the mobility of the peptide when bound to HLA-B*3501 or HLA-B*3508. In HLA-B*3501, the bulged section of the peptide was disordered, whereas in HLA-B*3508 the bulged epitope adopted an ordered conformation. Collectively, these data demonstrate not only that mobile MHC-bound peptides can be highly immunogenic but can also stimulate an extremely biased TCR repertoire. In addition, TCR {alpha}-chain usage is shown to play a critical role in controlling MHC restriction between closely related allomorphs.

MHC class I bound to an immunodominant Theileria parva epitope demonstrates unconventional presentation to T cell receptors

T cell receptor (TCR) recognition of peptide-MHC class I (pMHC) complexes is a crucial event in the adaptive immune response to pathogens. Peptide epitopes often display a strong dominance hierarchy, resulting in focusing of the response on a limited number of the most dominant epitopes. Such T cell responses may be additionally restricted by particular MHC alleles in preference to others. We have studied this poorly understood phenomenon using Theileria parva, a protozoan parasite that causes an often fatal lymphoproliferative disease in cattle. Despite its antigenic complexity, CD8+ T cell responses induced by infection with the parasite show profound immunodominance, as exemplified by the Tp1(214-224) epitope presented by the common and functionally important MHC class I allele N*01301. We present a high-resolution crystal structure of this pMHC complex, demonstrating that the peptide is presented in a distinctive raised conformation. Functional studies using CD8+ T cell clones show that this impacts significantly on TCR recognition. The unconventional structure is generated by a hydrophobic ridge within the MHC peptide binding groove, found in a set of cattle MHC alleles. Extremely rare in all other species, this feature is seen in a small group of mouse MHC class I molecules. The data generated in this analysis contribute to our understanding of the structural basis for T cell-dependent immune responses, providing insight into what determines a highly immunogenic p-MHC complex, and hence can be of value in prediction of antigenic epitopes and vaccine design.

Human Kallikrein 4 signal peptide induces cytotoxic T cell responses in healthy donors and prostate cancer patients.

Immunotherapy is a promising new treatment for patients with advanced prostate and ovarian cancer, but its application is limited by the lack of suitable target antigens that are recognized by CD8+ cytotoxic T lymphocytes (CTL). Human kallikrein 4 (KLK4) is a member of the kallikrein family of serine proteases that is significantly overexpressed in malignant versus healthy prostate and ovarian tissue, making it an attractive target for immunotherapy. We identified a naturally processed, HLA-A*0201-restricted peptide epitope within the signal sequence region of KLK4 that induced CTL responses in vitro in most healthy donors and prostate cancer patients tested. These CTL lysed HLA-A*0201+ KLK4 + cell lines and KLK4 mRNA-transfected monocyte-derived dendritic cells. CTL specific for the HLA-A*0201-restricted KLK4 peptide were more readily expanded to a higher frequency in vitro compared to the known HLA-A*0201-restricted epitopes from prostate cancer antigens; prostate-specific antigen (PSA), prostate-specific membrane antigen (PSMA) and prostatic acid phosphatase (PAP). These data demonstrate that KLK4 is an immunogenic molecule capable of inducing CTL responses and identify it as an attractive target for prostate and ovarian cancer immunotherapy.

The enigma of IgE+ B-cell memory in human subjects

Our understanding of the origin and fate of the IgE-switched B cell has been markedly improved by studies in mouse models. The immediate precursor of the IgE-switched B cell is either a relatively naive nonswitched B cell or a mature IgG-switched B cell. These 2 routes are referred to as the direct and indirect pathways, respectively. IgE responses derived from each pathway differ significantly, largely reflecting the difference in time spent in a germinal center and thus time for clonal expansion, somatic hypermutation, affinity maturation, and acquisition of a memory phenotype. The clinical and therapeutic implications for IgE responses in human subjects are still a matter of debate, largely because the immunization procedures used in the animal models are significantly different from classical atopic sensitization to allergens from pollen and mites. On the basis of the limited information available, it seems likely that these atopic IgE responses are characterized by a relatively low IgG/IgE ratio, low B-cell memory, and modest affinity maturation, which fits well with the direct switching pathway. It is still unresolved how the IgE response evolves to cover a wide epitope repertoire involving many epitopes per allergen, as well as many different allergens from a single allergen source. © 2013 American Academy of Allergy, Asthma & Immunology.

Introduction: Epitope mimicry as a component cause of autoimmune disease

The causes of autoimmune diseases have yet to be fully elucidated. Autoantibodies, autoreactive T cell responses, the presence of a predisposing major histocompatibility complex (MHC) haplotype and responsiveness to corticosteroids are features, and some are possibly contributory causes of autoimmune disease. The most challenging question is how autoimmune diseases are triggered. Molecular mimicry of host cell determinants by epitopes of infectious agents with ensuing cross-reactivity is one of the most popular yet still controversial theories for the initiation of autoimmune diseases [1]. Throughout the 1990s, hundreds of research articles focusing to various extents on epitope mimicry, as it is more accurately described in an immunological context, were published annually. Many of these articles presented data that were consistent with the hypothesis of mimicry but that did not actually prove the theory. Other equally convincing reports indicated that epitope mimicry was not the cause of the autoimmune disease despite sequence similarity between molecules of infectious agents and the host. Some 20 years ago, Rothman [2] proposed a model for disease causation and I have used this as a framework to examine the role of epitope mimicry in the development of autoimmune disease. The thesis of Rothman’s model is that an effect, in this instance autoimmune disease, arises as a result of a cause. In most cases, multiple-component causes contribute synergistically to yield the effect, and each of these components alone is insufficient as a cause. Logically, some component causes, such as the presence of a particular autoimmune response, are also necessary causes.