Selection of vaccine-candidate peptides from Mycobacterium avium subsp. paratuberculosis by in silico prediction, in vitro T-cell line proliferation, and in vivo immunogenicity.

Mycobacterium avium subspecies paratuberculosis (MAP) is a global concern in modern livestock production worldwide. The available vaccines against paratuberculosis do not offer optimal protection and interfere with the diagnosis of bovine tuberculosis. The aim of this study was to identify immunogenic MAP-specific peptides that do not interfere with the diagnosis of bovine tuberculosis. Initially, 119 peptides were selected by either (1) identifying unique MAP peptides that were predicted to bind to bovine major histocompatibility complex class II (MHC-predicted peptides) or (2) selecting hydrophobic peptides unique to MAP within proteins previously shown to be immunogenic (hydrophobic peptides). Subsequent testing of peptide-specific CD4+ T-cell lines from MAP-infected, adult goats vaccinated with peptides in cationic liposome adjuvant pointed to 23 peptides as being most immunogenic. These peptides were included in a second vaccine trial where three groups of eight healthy goat kids were vaccinated with 14 MHC-predicted peptides, nine hydrophobic peptides, or no peptides in o/w emulsion adjuvant. The majority of the MHC-predicted (93%) and hydrophobic peptides (67%) induced interferon-gamma (IFN-γ) responses in at least one animal. Similarly, 86% of the MHC-predicted and 89% of the hydrophobic peptides induced antibody responses in at least one goat. The immunization of eight healthy heifers with all 119 peptides formulated in emulsion adjuvant identified more peptides as immunogenic, as peptide specific IFN-γ and antibody responses in at least one heifer was found toward 84% and 24% of the peptides, respectively. No peptide-induced reactivity was found with commercial ELISAs for detecting antibodies against Mycobacterium bovis or MAP or when performing tuberculin skin testing for bovine tuberculosis. The vaccinated animals experienced adverse reactions at the injection site; thus, it is recommend that future studies make improvements to the vaccine formulation. In conclusion, immunogenic MAP-specific peptides that appeared promising for use in a vaccine against paratuberculosis without interfering with surveillance and trade tests for bovine tuberculosis were identified by in silico analysis and ex vivo generation of CD4+ T-cell lines and validated by the immunization of goats and cattle. Future studies should test different peptide combinations in challenge trials to determine their protective effect and identify the most MHC-promiscuous vaccine candidates.

Allergen-specific IgG+ memory B cells are temporally linked to IgE memory responses.

BACKGROUND: IgE is the least abundant immunoglobulin and tightly regulated, and IgE-producing B cells are rare. The cellular origin and evolution of IgE responses are poorly understood. OBJECTIVE: The cellular and clonal origin of IgE memory responses following mucosal allergen exposure by sublingual immunotherapy (SLIT) were investigated. METHODS: In a randomized double-blind, placebo-controlled, time course SLIT study, PBMCs and nasal biopsy samples were collected from 40 adults with seasonal allergic rhinitis at baseline and at 4, 8, 16, 28, and 52 weeks. RNA was extracted from PBMCs, sorted B cells, and nasal biopsy samples for heavy chain variable gene repertoire sequencing. Moreover, mAbs were derived from single B-cell transcriptomes. RESULTS: Combining heavy chain variable gene repertoire sequencing and single-cell transcriptomics yielded direct evidence of a parallel boost of 2 clonally and functionally related B-cell subsets of short-lived IgE+ plasmablasts and IgG+ memory B cells. Mucosal grass pollen allergen exposure by SLIT resulted in highly diverse IgE and IgGE repertoires. These were extensively mutated and appeared relatively stable as per heavy chain isotype, somatic hypermutations, and clonal composition. Single IgGE+ memory B-cell and IgE+ preplasmablast transcriptomes encoded antibodies that were specific for major grass pollen allergens and able to elicit basophil activation at very low allergen concentrations. CONCLUSION: For the first time, we have shown that on mucosal allergen exposure, human IgE memory resides in allergen-specific IgG+ memory B cells. These cells rapidly switch isotype, expand into short-lived IgE+ plasmablasts, and serve as a potential target for therapeutic intervention.

Diverse and highly cross-reactive T-cell responses in ragweed allergic patients independent of geographical region.

BACKGROUND: Ragweed frequently causes seasonal allergies in North America and Europe. In the United States, several related ragweed species with diverse geographical distribution cause allergic symptoms. Cross-reactivity towards related ragweed species of IgE and treatment-induced IgG4 has been demonstrated previously. However, less is known about the underlying T-cell cross-reactivity. METHODS: The allergen content of ragweed extracts was determined by mass spectrometry and related to T-cell epitopes of Amb a allergens (group 1, 3, 4, 5, 8 and 11) in 20 American ragweed allergic patients determined by FluoroSpot and proliferation assays. T-cell responses to 50 frequently recognized Amb a-derived T-cell epitopes and homologous peptides from western ragweed (Amb p), giant ragweed (Amb t) and mugwort (Art v) were investigated in an additional 11 American and 14 Slovakian ragweed allergic donors. RESULTS: Ragweed extracts contained all known allergens and isoallergens thereof. Donor T-cell responses were diverse and directed against all Amb a 1 isoallergens and to most minor allergens investigated. Similar response patterns were seen in American and Slovakia donors. Several epitopes were cross-reactive between isoallergens and ragweed species, some even including mugwort. T-cell cross-reactivity generally correlated with allergen sequence homology. CONCLUSION: T-cell epitopes of multiple allergens/isoallergens are involved in the diverse T-cell responses in ragweed allergic individuals. T-cell lines were highly cross-reactive to epitopes of related ragweed species without any apparent geographical response bias. These data support that different ragweed species can be considered an allergen homology group with Amb a as the representative species regarding diagnosis as well as allergy immunotherapy.

MHCcluster, a method for functional clustering of MHC molecules.

The identification of peptides binding to major histocompatibility complexes (MHC) is a critical step in the understanding of T cell immune responses. The human MHC genomic region (HLA) is extremely polymorphic comprising several thousand alleles, many encoding a distinct molecule. The potentially unique specificities remain experimentally uncharacterized for the vast majority of HLA molecules. Likewise, for nonhuman species, only a minor fraction of the known MHC molecules have been characterized. Here, we describe a tool, MHCcluster, to functionally cluster MHC molecules based on their predicted binding specificity. The method has a flexible web interface that allows the user to include any MHC of interest in the analysis. The output consists of a static heat map and graphical tree-based visualizations of the functional relationship between MHC variants and a dynamic TreeViewer interface where both the functional relationship and the individual binding specificities of MHC molecules are visualized. We demonstrate that conventional sequence-based clustering will fail to identify the functional relationship between molecules, when applied to MHC system, and only through the use of the predicted binding specificity can a correct clustering be found. Clustering of prevalent HLA-A and HLA-B alleles using MHCcluster confirms the presence of 12 major specificity groups (supertypes) some however with highly divergent specificities. Importantly, some HLA molecules are shown not to fit any supertype classification. Also, we use MHCcluster to show that chimpanzee MHC class I molecules have a reduced functional diversity compared to that of HLA class I molecules. MHCcluster is available at www.cbs.dtu.dk/services/MHCcluster-2.0.

In silico peptide-binding predictions of passerine MHC class I reveal similarities across distantly related species, suggesting convergence on the level of protein function.

The major histocompatibility complex (MHC) genes are the most polymorphic genes found in the vertebrate genome, and they encode proteins that play an essential role in the adaptive immune response. Many songbirds (passerines) have been shown to have a large number of transcribed MHC class I genes compared to most mammals. To elucidate the reason for this large number of genes, we compared 14 MHC class I alleles (α1-α3 domains), from great reed warbler, house sparrow and tree sparrow, via phylogenetic analysis, homology modelling and in silico peptide-binding predictions to investigate their functional and genetic relationships. We found more pronounced clustering of the MHC class I allomorphs (allele specific proteins) in regards to their function (peptide-binding specificities) compared to their genetic relationships (amino acid sequences), indicating that the high number of alleles is of functional significance. The MHC class I allomorphs from house sparrow and tree sparrow, species that diverged 10 million years ago (MYA), had overlapping peptide-binding specificities, and these similarities across species were also confirmed in phylogenetic analyses based on amino acid sequences. Notably, there were also overlapping peptide-binding specificities in the allomorphs from house sparrow and great reed warbler, although these species diverged 30 MYA. This overlap was not found in a tree based on amino acid sequences. Our interpretation is that convergent evolution on the level of the protein function, possibly driven by selection from shared pathogens, has resulted in allomorphs with similar peptide-binding repertoires, although trans-species evolution in combination with gene conversion cannot be ruled out.

Bioinformatics identification of antigenic peptide: predicting the specificity of major MHC class I and II pathway players.

Bioinformatics methods for immunology have become increasingly used over the last decade and now form an integrated part of most epitope discovery projects. This wide usage has led to the confusion of defining which of the many methods to use for what problems. In this chapter, an overview is given focusing on the suite of tools developed at the Technical University of Denmark.

Characterization of HIV-specific CD4+ T cell responses against peptides selected with broad population and pathogen coverage.

CD4+ T cells orchestrate immunity against viral infections, but their importance in HIV infection remains controversial. Nevertheless, comprehensive studies have associated increase in breadth and functional characteristics of HIV-specific CD4+ T cells with decreased viral load. A major challenge for the identification of HIV-specific CD4+ T cells targeting broadly reactive epitopes in populations with diverse ethnic background stems from the vast genomic variation of HIV and the diversity of the host cellular immune system. Here, we describe a novel epitope selection strategy, PopCover, that aims to resolve this challenge, and identify a set of potential HLA class II-restricted HIV epitopes that in concert will provide optimal viral and host coverage. Using this selection strategy, we identified 64 putative epitopes (peptides) located in the Gag, Nef, Env, Pol and Tat protein regions of HIV. In total, 73% of the predicted peptides were found to induce HIV-specific CD4+ T cell responses. The Gag and Nef peptides induced most responses. The vast majority of the peptides (93%) had predicted restriction to the patient's HLA alleles. Interestingly, the viral load in viremic patients was inversely correlated to the number of targeted Gag peptides. In addition, the predicted Gag peptides were found to induce broader polyfunctional CD4+ T cell responses compared to the commonly used Gag-p55 peptide pool. These results demonstrate the power of the PopCover method for the identification of broadly recognized HLA class II-restricted epitopes. All together, selection strategies, such as PopCover, might with success be used for the evaluation of antigen-specific CD4+ T cell responses and design of future vaccines.

Immune epitope database analysis resource.

The immune epitope database analysis resource (IEDB-AR: http://tools.iedb.org) is a collection of tools for prediction and analysis of molecular targets of T- and B-cell immune responses (i.e. epitopes). Since its last publication in the NAR webserver issue in 2008, a new generation of peptide:MHC binding and T-cell epitope predictive tools have been added. As validated by different labs and in the first international competition for predicting peptide:MHC-I binding, their predictive performances have improved considerably. In addition, a new B-cell epitope prediction tool was added, and the homology mapping tool was updated to enable mapping of discontinuous epitopes onto 3D structures. Furthermore, to serve a wider range of users, the number of ways in which IEDB-AR can be accessed has been expanded. Specifically, the predictive tools can be programmatically accessed using a web interface and can also be downloaded as software packages.

NetMHCcons: a consensus method for the major histocompatibility complex class I predictions.

A key role in cell-mediated immunity is dedicated to the major histocompatibility complex (MHC) molecules that bind peptides for presentation on the cell surface. Several in silico methods capable of predicting peptide binding to MHC class I have been developed. The accuracy of these methods depends on the data available characterizing the binding specificity of the MHC molecules. It has, moreover, been demonstrated that consensus methods defined as combinations of two or more different methods led to improved prediction accuracy. This plethora of methods makes it very difficult for the non-expert user to choose the most suitable method for predicting binding to a given MHC molecule. In this study, we have therefore made an in-depth analysis of combinations of three state-of-the-art MHC-peptide binding prediction methods (NetMHC, NetMHCpan and PickPocket). We demonstrate that a simple combination of NetMHC and NetMHCpan gives the highest performance when the allele in question is included in the training and is characterized by at least 50 data points with at least ten binders. Otherwise, NetMHCpan is the best predictor. When an allele has not been characterized, the performance depends on the distance to the training data. NetMHCpan has the highest performance when close neighbours are present in the training set, while the combination of NetMHCpan and PickPocket outperforms either of the two methods for alleles with more remote neighbours. The final method, NetMHCcons, is publicly available at www.cbs.dtu.dk/services/NetMHCcons , and allows the user in an automatic manner to obtain the most accurate predictions for any given MHC molecule.

Reliable B cell epitope predictions: impacts of method development and improved benchmarking.

The interaction between antibodies and antigens is one of the most important immune system mechanisms for clearing infectious organisms from the host. Antibodies bind to antigens at sites referred to as B-cell epitopes. Identification of the exact location of B-cell epitopes is essential in several biomedical applications such as; rational vaccine design, development of disease diagnostics and immunotherapeutics. However, experimental mapping of epitopes is resource intensive making in silico methods an appealing complementary approach. To date, the reported performance of methods for in silico mapping of B-cell epitopes has been moderate. Several issues regarding the evaluation data sets may however have led to the performance values being underestimated: Rarely, all potential epitopes have been mapped on an antigen, and antibodies are generally raised against the antigen in a given biological context not against the antigen monomer. Improper dealing with these aspects leads to many artificial false positive predictions and hence to incorrect low performance values. To demonstrate the impact of proper benchmark definitions, we here present an updated version of the DiscoTope method incorporating a novel spatial neighborhood definition and half-sphere exposure as surface measure. Compared to other state-of-the-art prediction methods, Discotope-2.0 displayed improved performance both in cross-validation and in independent evaluations. Using DiscoTope-2.0, we assessed the impact on performance when using proper benchmark definitions. For 13 proteins in the training data set where sufficient biological information was available to make a proper benchmark redefinition, the average AUC performance was improved from 0.791 to 0.824. Similarly, the average AUC performance on an independent evaluation data set improved from 0.712 to 0.727. Our results thus demonstrate that given proper benchmark definitions, B-cell epitope prediction methods achieve highly significant predictive performances suggesting these tools to be a powerful asset in rational epitope discovery. The updated version of DiscoTope is available at www.cbs.dtu.dk/services/DiscoTope-2.0.

Predictions versus high-throughput experiments in T-cell epitope discovery: competition or synergy?

Prediction methods as well as experimental methods for T-cell epitope discovery have developed significantly in recent years. High-throughput experimental methods have made it possible to perform full-length protein scans for epitopes restricted to a limited number of MHC alleles. The high costs and limitations regarding the number of proteins and MHC alleles that are feasibly handled by such experimental methods have made in silico prediction models of high interest. MHC binding prediction methods are today of a very high quality and can predict MHC binding peptides with high accuracy. This is possible for a large range of MHC alleles and relevant length of binding peptides. The predictions can easily be performed for complete proteomes of any size. Prediction methods are still, however, dependent on good experimental methods for validation, and should merely be used as a guide for rational epitope discovery. We expect prediction methods as well as experimental validation methods to continue to develop and that we will soon see clinical trials of products whose development has been guided by prediction methods.

Human leukocyte antigen (HLA) class I restricted epitope discovery in yellow fewer and dengue viruses: importance of HLA binding strength.

Epitopes from all available full-length sequences of yellow fever virus (YFV) and dengue fever virus (DENV) restricted by Human Leukocyte Antigen class I (HLA-I) alleles covering 12 HLA-I supertypes were predicted using the NetCTL algorithm. A subset of 179 predicted YFV and 158 predicted DENV epitopes were selected using the EpiSelect algorithm to allow for optimal coverage of viral strains. The selected predicted epitopes were synthesized and approximately 75% were found to bind the predicted restricting HLA molecule with an affinity, K(D), stronger than 500 nM. The immunogenicity of 25 HLA-A*02:01, 28 HLA-A*24:02 and 28 HLA-B*07:02 binding peptides was tested in three HLA-transgenic mice models and led to the identification of 17 HLA-A*02:01, 4 HLA-A*2402 and 4 HLA-B*07:02 immunogenic peptides. The immunogenic peptides bound HLA significantly stronger than the non-immunogenic peptides. All except one of the immunogenic peptides had K(D) below 100 nM and the peptides with K(D) below 5 nM were more likely to be immunogenic. In addition, all the immunogenic peptides that were identified as having a high functional avidity had K(D) below 20 nM. A*02:01 transgenic mice were also inoculated twice with the 17DD YFV vaccine strain. Three of the YFV A*02:01 restricted peptides activated T-cells from the infected mice in vitro. All three peptides that elicited responses had an HLA binding affinity of 2 nM or less. The results indicate the importance of the strength of HLA binding in shaping the immune response.

Prediction of epitopes using neural network based methods.

In this paper, we describe the methodologies behind three different aspects of the NetMHC family for prediction of MHC class I binding, mainly to HLAs. We have updated the prediction servers, NetMHC-3.2, NetMHCpan-2.2, and a new consensus method, NetMHCcons, which, in their previous versions, have been evaluated to be among the very best performing MHC:peptide binding predictors available. Here we describe the background for these methods, and the rationale behind the different optimization steps implemented in the methods. We go through the practical use of the methods, which are publicly available in the form of relatively fast and simple web interfaces. Furthermore, we will review results obtained in actual epitope discovery projects where previous implementations of the described methods have been used in the initial selection of potential epitopes. Selected potential epitopes were all evaluated experimentally using ex vivo assays.

NetTurnP--neural network prediction of beta-turns by use of evolutionary information and predicted protein sequence features.

UNLABELLED: ß-turns are the most common type of non-repetitive structures, and constitute on average 25% of the amino acids in proteins. The formation of ß-turns plays an important role in protein folding, protein stability and molecular recognition processes. In this work we present the neural network method NetTurnP, for prediction of two-class ß-turns and prediction of the individual ß-turn types, by use of evolutionary information and predicted protein sequence features. It has been evaluated against a commonly used dataset BT426, and achieves a Matthews correlation coefficient of 0.50, which is the highest reported performance on a two-class prediction of ß-turn and not-ß-turn. Furthermore NetTurnP shows improved performance on some of the specific ß-turn types. In the present work, neural network methods have been trained to predict ß-turn or not and individual ß-turn types from the primary amino acid sequence. The individual ß-turn types I, I', II, II', VIII, VIa1, VIa2, VIba and IV have been predicted based on classifications by PROMOTIF, and the two-class prediction of ß-turn or not is a superset comprised of all ß-turn types. The performance is evaluated using a golden set of non-homologous sequences known as BT426. Our two-class prediction method achieves a performance of: MCC=0.50, Qtotal=82.1%, sensitivity=75.6%, PPV=68.8% and AUC=0.864. We have compared our performance to eleven other prediction methods that obtain Matthews correlation coefficients in the range of 0.17-0.47. For the type specific ß-turn predictions, only type I and II can be predicted with reasonable Matthews correlation coefficients, where we obtain performance values of 0.36 and 0.31, respectively. CONCLUSION: The NetTurnP method has been implemented as a webserver, which is freely available at http://www.cbs.dtu.dk/services/NetTurnP/. NetTurnP is the only available webserver that allows submission of multiple sequences.

NetMHCIIpan-2.0 - Improved pan-specific HLA-DR predictions using a novel concurrent alignment and weight optimization training procedure.

BACKGROUND: Binding of peptides to Major Histocompatibility class II (MHC-II) molecules play a central role in governing responses of the adaptive immune system. MHC-II molecules sample peptides from the extracellular space allowing the immune system to detect the presence of foreign microbes from this compartment. Predicting which peptides bind to an MHC-II molecule is therefore of pivotal importance for understanding the immune response and its effect on host-pathogen interactions. The experimental cost associated with characterizing the binding motif of an MHC-II molecule is significant and large efforts have therefore been placed in developing accurate computer methods capable of predicting this binding event. Prediction of peptide binding to MHC-II is complicated by the open binding cleft of the MHC-II molecule, allowing binding of peptides extending out of the binding groove. Moreover, the genes encoding the MHC molecules are immensely diverse leading to a large set of different MHC molecules each potentially binding a unique set of peptides. Characterizing each MHC-II molecule using peptide-screening binding assays is hence not a viable option. RESULTS: Here, we present an MHC-II binding prediction algorithm aiming at dealing with these challenges. The method is a pan-specific version of the earlier published allele-specific NN-align algorithm and does not require any pre-alignment of the input data. This allows the method to benefit also from information from alleles covered by limited binding data. The method is evaluated on a large and diverse set of benchmark data, and is shown to significantly out-perform state-of-the-art MHC-II prediction methods. In particular, the method is found to boost the performance for alleles characterized by limited binding data where conventional allele-specific methods tend to achieve poor prediction accuracy. CONCLUSIONS: The method thus shows great potential for efficient boosting the accuracy of MHC-II binding prediction, as accurate predictions can be obtained for novel alleles at highly reduced experimental costs. Pan-specific binding predictions can be obtained for all alleles with know protein sequence and the method can benefit by including data in the training from alleles even where only few binders are known. The method and benchmark data are available at http://www.cbs.dtu.dk/services/NetMHCIIpan-2.0.

State of the art and challenges in sequence based T-cell epitope prediction.

Sequence based T-cell epitope predictions have improved immensely in the last decade. From predictions of peptide binding to major histocompatibility complex molecules with moderate accuracy, limited allele coverage, and no good estimates of the other events in the antigen-processing pathway, the field has evolved significantly. Methods have now been developed that produce highly accurate binding predictions for many alleles and integrate both proteasomal cleavage and transport events. Moreover have so-called pan-specific methods been developed, which allow for prediction of peptide binding to MHC alleles characterized by limited or no peptide binding data. Most of the developed methods are publicly available, and have proven to be very useful as a shortcut in epitope discovery. Here, we will go through some of the history of sequence-based predictions of helper as well as cytotoxic T cell epitopes. We will focus on some of the most accurate methods and their basic background.

Mice, men and MHC supertypes.

Mice are still the most used model organism in initial phases of vaccine design. Bioinformatics is becoming increasingly more important in vaccine development, both for the design of novel simplified epitope-based vaccines and also in order to understand the specific immune response of selected vaccine formulations. Toxoplasma gondii, an intracellular parasite, causes severe neurologic and ocular disease in congenitally infected and immunocompromised individuals. No protective vaccine exists against human toxoplasmosis. However, studies with mice have revealed immunodominant cytotoxic T lymphocyte epitopes originating from type II strains. These verified epitopes might be useful in human vaccines as the peptide binding repertoire of H-2L(d) MHC and MHCs belonging to the HLA-B07 supertype are very similar. Here, the results obtained using these epitopes in BALB/c as well as transgenic HLA-B*0702 mice are discussed. The stunning results obtained from the use of refined computational methods for the prediction of cytotoxic T lymphocyte epitopes are also discussed. The results highlight some important issues regarding both the use of mice but also the choice of bioinformatics methods in vaccine development.

CPHmodels-3.0--remote homology modeling using structure-guided sequence profiles.

CPHmodels-3.0 is a web server predicting protein 3D structure by use of single template homology modeling. The server employs a hybrid of the scoring functions of CPHmodels-2.0 and a novel remote homology-modeling algorithm. A query sequence is first attempted modeled using the fast CPHmodels-2.0 profile-profile scoring function suitable for close homology modeling. The new computational costly remote homology-modeling algorithm is only engaged provided that no suitable PDB template is identified in the initial search. CPHmodels-3.0 was benchmarked in the CASP8 competition and produced models for 94% of the targets (117 out of 128), 74% were predicted as high reliability models (87 out of 117). These achieved an average RMSD of 4.6 A when superimposed to the 3D structure. The remaining 26% low reliably models (30 out of 117) could superimpose to the true 3D structure with an average RMSD of 9.3 A. These performance values place the CPHmodels-3.0 method in the group of high performing 3D prediction tools. Beside its accuracy, one of the important features of the method is its speed. For most queries, the response time of the server is <20 min. The web server is available at http://www.cbs.dtu.dk/services/CPHmodels/.

Major histocompatibility complex class I binding predictions as a tool in epitope discovery.

SUMMARY: Over the last decade, in silico models of the major histocompatibility complex (MHC) class I pathway have developed significantly. Before, peptide binding could only be reliably modelled for a few major human or mouse histocompatibility molecules; now, high-accuracy predictions are available for any human leucocyte antigen (HLA) -A or -B molecule with known protein sequence. Furthermore, peptide binding to MHC molecules from several non-human primates, mouse strains and other mammals can now be predicted. In this review, a number of different prediction methods are briefly explained, highlighting the most useful and historically important. Selected case stories, where these 'reverse immunology' systems have been used in actual epitope discovery, are briefly reviewed. We conclude that this new generation of epitope discovery systems has become a highly efficient tool for epitope discovery, and recommend that the less accurate prediction systems of the past be abandoned, as these are obsolete.