Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

Open Access Research article

PeptX: Using Genetic Algorithms to optimize peptides for MHC binding

Bernhard Knapp1*, Verena Giczi12, Reiner Ribarics1 and Wolfgang Schreiner1

Author Affiliations

1 Center for Medical Statistics, Informatics and Intelligent Systems, Department for Biosimulation and Bioinformatics, Medical University of Vienna, Vienna, Austria

2 University of Applied Sciences, FH Campus Wien, Department of Bioengineering, Vienna, Austria

For all author emails, please log on.

BMC Bioinformatics 2011, 12:241  doi:10.1186/1471-2105-12-241


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/12/241


Received:13 January 2011
Accepted:17 June 2011
Published:17 June 2011

© 2011 Knapp et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The binding between the major histocompatibility complex and the presented peptide is an indispensable prerequisite for the adaptive immune response. There is a plethora of different in silico techniques for the prediction of the peptide binding affinity to major histocompatibility complexes. Most studies screen a set of peptides for promising candidates to predict possible T cell epitopes. In this study we ask the question vice versa: Which peptides do have highest binding affinities to a given major histocompatibility complex according to certain in silico scoring functions?

Results

Since a full screening of all possible peptides is not feasible in reasonable runtime, we introduce a heuristic approach. We developed a framework for Genetic Algorithms to optimize peptides for the binding to major histocompatibility complexes. In an extensive benchmark we tested various operator combinations. We found that (1) selection operators have a strong influence on the convergence of the population while recombination operators have minor influence and (2) that five different binding prediction methods lead to five different sets of "optimal" peptides for the same major histocompatibility complex. The consensus peptides were experimentally verified as high affinity binders.

Conclusion

We provide a generalized framework to calculate sets of high affinity binders based on different previously published scoring functions in reasonable runtime. Furthermore we give insight into the different behaviours of operators and scoring functions of the Genetic Algorithm.

Background

Antigen presenting cells (APCs) present peptides via their major histocompatibility complex (MHC) to T cell receptors (TCRs) of T cells which play an essential role in the adaptive immune system [1]. Before any recognition between T cell and APC can take place the peptides need to be processed within the APC and afterwards presented in a stable way on the cell surface of the APC. For most of these steps ample prediction methods exist [2-5]. In this context the binding prediction between MHC and the presented peptides is usually referred as T cell epitope prediction. In a usual workflow one wants to test different peptides or even possible fragments of a whole protein for its binding affinity to a given MHC. After this rough pre-selection step the most promising candidates are then tested in wet-lab experiments for their definitive binding affinity and applicability. The success rate of these approaches is discussed abundantly in the literature: While it is known that there is still much space for improvement of B cell epitope prediction methods [6] and that MHC class I T cell epitope prediction works quite well [2,7-9], the opinions on the success of MHC class II T cell epitope prediction differ. Various reviews conclude that it lacks reliability [10-13], though others are more enthusiastic about the results obtainable with MHC class II prediction tools [3].

While most methods are sequence-based and it is generally believed that sequence-based approaches are more successful [8], there were also some structure-based approaches reported [14]. They range from molecular dynamics (MD) methods [15-17] over quantitative structure-activity relationship (QSAR) [18] methods to empirical scoring methods [19-21].

By contrast, when asking the question vice versa the scientific challenge changes: One does not want to test a specific set of peptides for their binding affinities, but knows a specific MHC allele and wants to determine a set of peptides which have a very high binding affinity to this MHC allele. Applications could be found for example in peptide immunotherapy (PIT), in allergy [22-24], cancer therapy [25], or peptide vaccines [26] by indentifying altered peptide ligands with favourable properties for the respective purpose.

Given the fact that for MHC class I 209 and for MHC class II 2018 or even more different peptides exist, it is obvious that not each individual peptide can be tested. This limitation also includes in silico techniques since even if we assume that modern computers can test 10 peptides per second it would still take more than 1 600 years (209/(10×60×60×24×365)) to predict the binding affinity for all theoretically possible peptides binding to one single MHC class I allele. The runtime would further increase if the effect on the interaction between peptide/MHC (pMHC) and TCR would additionally be predicted.

Several approaches to address the challenge of optimization for pMHC/TCR interaction were reported: In early methods Alexander et al. employed the main HLA-DR anchors to increase the peptide/MHC binding affinity as well as the use of bulky and charged residues to increase T cell recognition [27,28]. Shang et al. used computational alanine scanning to indentify hotspots which were then systematically substituted and scored to optimize a tumor immunodominant epitope [29]. Also for the peptide and MHC interaction several approaches were reported: Reche et al. published a webserver for the formulation of multi-epitope vaccines [30]. Toussaint et al. developed a mathematical framework on the basis of integer linear programming to obtain good candidates for epitope-based vaccines. In this context they define good as a combination of mutation tolerance, allele coverage, antigen coverage, and antigen processing [31]. Contrary, Parker et al. implemented an optimization algorithm to find point mutations to reduce the immunogenicity of a protein while maintaining stability and function [32]. Lazar et al. presented an optimization approach for antibody humanization [33]. Bhasin et al. reported quantitative matrices on the basis of binding and non binding peptides sets [34]. Similarly, in the studies of Doytchinova et al. and Walshe et al. the preferred amino acids for peptide positions were selected based on experimental data sets. This leads to a set of a few hundred peptides which were then evaluated in silico and experimentally [35,36]. Guan et al. reported amino acid descriptors to characterize the interaction between peptide and MHC. On this basis of the defined binding model high-affinity binders were designed [37].

However, in these approaches the search space for optimal binders is frequently reduced by an initial selection of preferred amino acids for several positions. Given this limitation and the fact that the total search space would be 209 combinations we propose the use of a Genetic Algorithm (GA) to be able to investigate the whole search space in justifiable runtime. Hereby the task reformulates to finding the most efficient and reliable modes of GAs. Hence, we implemented a framework "PeptX" for optimizing GAs for the prediction of sets of peptides with high binding affinities within reasonable runtime. By means of the framework "PeptX" we evaluated which combination of parameters and operators yields the most rapid convergence of the GA towards a set of high binding affinity peptides.

It is not the major aim of this study to compare different scoring functions since this issue was already addressed several times in the recent literature e.g. MHC class I [9,38-41] and MHC class II [3,10-13]. Even a "machine learning in immunology competition" was organized for this purpose (http://www.kios.org.cy/ICANN09/MLI.html webcite). Since there is already a plethora of benchmarks for prediction methods our study is focused on a generalized approach on how to predict huge sets of high affinity binders on the basis of these previously published methods.

Methods

We implemented the C++ framework "PeptX" for GAs to optimize a set of peptides for MHC binding and tested it for HLA-A*02:01. We implemented the whole framework in an object-oriented way which makes it easily maintainable and extendable. The workflow of the GA is illustrated in Figure 1 and described in detail in the subsequent sections.

thumbnailFigure 1. Workflow of the GA.

Initialization of random start generation and removal of identical individuals

We generated a random set of peptides as initial generation. Each individual consisted of a string of 9 characters, each of them representing one position in the peptide. Therefore, we randomly selected for each position in the peptide one of the 20 essential amino acids. To avoid too dominant individuals (genetic drift) we purged the whole population from identical peptides at each generation. The removed individuals were replaced by new random ones. This purging step ensures that the algorithm will not get stuck in a local optimum and the genetic diversity will remain over the whole runtime of the GA.

Fitness evaluation of the individuals

In the next step we assigned a fitness value to each of the peptides by employing a scoring function for peptide/MHC binding. We utilized the following sequence-based scoring functions: Immune Epitope Database (IEDB) [42], SYFPEITHI [43], SVMHC [44], SVRMHC [45]. Additionally, we employed the structure-based scoring function XScore [46] in combination with the threading technique of SCATD [47] for the construction of the 3D peptide/MHC model. This combination was found most appropriate for structural peptide/MHC binding predictions [40,48].

Selection of parent individuals

Individuals were selected as parents for the next generation based on the previously described fitness evaluation. We implemented 7 different selection methods according to the literature which are briefly summarized in the following. For details of the single algorithms we refer to [49-51]:

1. Proportional (roulette wheel) selection: To each individual an area on a roulette wheel is assigned depending on its fitness value. Individuals with higher value have a higher probability to be selected than individuals with lower fitness value.

2. Linear rank selection: The individuals obtain a rank in correspondence to their fitness value. The selection is performed on the basis of this rank.

3. Binary tournament selection: Two randomly selected individuals compete for their selection where the one with the higher fitness value wins without any other stochastic influences.

4. Random selection: The individuals are selected in a completely random way.

5. Best percent selection: The n best percent of the population are chosen straight forwardly.

6. Q-tournament selection: All individuals participate in q tournaments, where the individuals with the most victories are selected.

7. Stochastic universal sampling: Similar to the proportional selection, every individual obtains a segment on a roulette wheel according to its fitness value. However, it is turned only one time with n-balls where n is the number of individuals in the population.

Recombination (crossover) of parent individuals

To recombine the previously selected parent individuals we implemented 6 different crossover operators according to [51,52]:

1. Single point crossover: The amino acid sequence of parent 1 (P1) and parent 2 (P2) are cleaved at one randomly chosen position. The four parts of P1 and P2 are recombined at the cleavage point in a way that each child contains one part of P1 and the other part of P2.

2. Double point crossover: P1 and P2 are recombined at two randomly chosen positions.

3. Distance bisector crossover: P1 and P2 are recombined in the middle (e.g. for MHC class I: position 5 of the peptide).

4. Multi point crossover: P1 and P2 are recombined at r randomly chosen positions where r itself is also a random number (r ≤ 9).

5. Uniform crossover: Related to the multi-point crossover, however, each position within the sequence obtains a randomly assigned probability for recombination. If this probability exceeds a certain threshold, a recombination at this position occurs.

6. Shuffle crossover: Similar to the double point crossover P1 and P2 are recombined at two randomly selected positions. However, before recombination the amino acids are shuffled in both parents. After recombination the amino acids are unshuffled.

Mutations in the offspring

To maintain genetic diversity we performed mutations on each offspring generation. For this purpose we implemented 2 different operators [51]:

1. "Single point amino acid mutation": One amino acid is substituted randomly by another amino acid which belongs to another amino acid class (unpolar, polar, basic, or acid).

2. "Single point nucleobase mutation": A mutation is introduced in the triplet of the genetic code. Stop codons are avoided.

We applied the mutation with a probability of 3% for each amino acid position.

Stop criterion

Although it turned out that the most reasonable combinations of parameters let the GA converge much earlier (see result section) we always simulated 50 generations, so as to allow for comparability between the different parameterizations.

Results and Discussion

Operators and their combinations

In total we evaluated 7 selection operators, 6 recombination operators, 5 fitness functions, and 2 mutation operators in all possible combinations. This yields 420 (7 × 6 × 5 × 2) independent runs of the GA. In each run we calculated 50 generations, each generation consisting of 100 individuals (peptides). Therefore, in total we performed 2 100 000 (420 × 50 × 100) different peptide evaluations. The results are illustrated in Figure 2 and will be summarized in the subsequent sections grouped by selection operator, recombination operator and scoring function. Since Figure 2 contains a plethora of data we additionally depict several numerical details of these data in additional file 1.

thumbnailFigure 2. Median fitness score of the individuals over generations. For each generation we calculate the median of the fitness scores of the 100 individuals per generation. At generation 1 the individuals are random while at generation 50 the individuals should be optimized according to the corresponding scoring function: (a) IEDB, (b) SYFPEITHI, (c) SVRMHC, (d) SVMHC and, (e) XSCORE. For all scoring functions except IEDB the optimum is a maximum (higher values indicate a higher binding affinity, e.g. pIC50-value). In IEDB the optimum is a minimum (lower values indicate a higher binding affinity, e.g. IC50-value). The selection operators are grouped by colour while the recombination operators are grouped by marker symbols. Only the data for the "single point amino acid mutation" is shown.

Additional file 1. Numerical details of individual operator combinations. This file contains the numerical details of Figure 2. The median value over the first 10 generations, the median value over the last 10 generations, the maximum median, and the minimum median is shown.

Format: DOC Size: 276KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Selection operators have a strong influence on the convergence

It can be seen that most operator combinations work well and converge within the first 10 to 20 generations. However, there are a few exceptions: Obviously, a random selection does not improve fitness in an evolutionary process. We included this operator only as negative control. Additionally, the roulette wheel selection operator is, in its basic implementation, not able to generate enough evolutionary pressure to lead the population to convergence. This can be explained by the fact that, in the present setting, the whole roulette wheel was divided into small sections for each individual where the area of the section is proportional to the fitness value. The theoretical maximum for IC50 values is at around 50 000 while binding peptides usually have IC50 values below 500. Therefore, the area for a binding individual is not much larger than for non binders closely above 500. These proportions of the roulette wheel seem to hamper the evolutionally pressure on the population. A possibility to improve the performance of the roulette wheel selection operator would be a transformation of the distribution of the IC50 values. In contrast, the best percent selection usually converges fastest, however, it seems that the genetic pressure is too strong and the population loses the genetic diversity after reaching a certain optimum and cannot deal with the random mutations and removal of identical individuals anymore. Therefore, after a very fast convergence of the best percent selection between generation 1 and 10 the results are becoming worse between generations 10 and 20. Quite good results are usually achieved by the tournament selection operators. Although they are converging slightly slower than the best percent selection, they reach a better and more stable optimum.

Recombination operators have minor influence on the convergence

The influence of recombination operators on the convergence of the population is limited. The results strongly depend on the interaction with the selection operator. There are examples for a good convergence as well for a bad one for nearly each of the recombination operators. The only exception is the distance bisector cross over whose convergence is never among the top performers. This result was to be expected since if the peptide is always cleaved exactly in the middle it takes much longer to find the global optimum.

Different mutation operators did not have a significant influence on the population

The usage of "single point amino acid mutation" versus "single point nucleobase mutation" method did not influence the convergence of the populations. Operators converging with the "single point amino acid mutation" did also converge with the "single point nucleobase mutation". The same also applies for non converging operators. Hence, we show only the data for the "single point amino acid mutation" in Figure 2.

Different scoring functions lead to diverse sets of optimized peptides

The most severe influence is created by the scoring function used for the calculation of the fitness values. Although it is not purposeful to compare the convergence between different scoring functions, it is interesting to investigate individuals which are predicted as optimal for HLA-A*02:01 by different scoring functions (see Figure 3). In the subsequent paragraphs they are discussed and compared to experimentally known anchor residues determined by Falk et al. [53]. These residues are described in the SYFPEITHI [43] database as L or M for position 2, L or V for position 9, and V as auxiliary anchor for position 6.

thumbnailFigure 3. Sequence logos of the best individuals. Logos created by WebLogo [56] on the basis of the best individuals of generation number 50 of each operator combination grouped by scoring function: (a) IEDB, (b) SYFPEITHI, (c) SVRMHC, (d) SVMHC and, (e) XSCORE.

IEDB

At positions 1, 7, 8 and 9 the scoring function exhibits a clear preference for one specific amino acid. The preferred residues for the anchor positions 2 and 9 are in agreement with the experimental data. In contrast, the algorithm also proposed Q, a rather polar amino acid, for position 2. Furthermore, it is noteworthy that there is a tendency towards amino acids with a higher Van-der-Waals volume. The consensus sequence would exhibit a volume of 1 194 Å3.

SYFPEITHI

The SYFPEITHI scoring function tends to prefer small, polar, and apolar amino acids while retaining a rather low background noise. The anchors at positions 2, 6 and 9 are clearly the same as described in the literature, however, this is not surprising since the SYFPEITHI database is partly based on the data of Falk et al. [53]. The Van-der-Waals volume of the consensus sequence is 877 Å3.

SVRMHC

This scoring function exhibits a rather promiscuous impression. There is no strong preference for a single residue. However, also here L and M are strong at position 2 while L and V are present in position 9. Besides, most of the remaining top performing residues for position 9 fall into the same class as V or L. The Van-der-Waals volume of the consensus sequence is 1 066 Å3.

SVMHC

While most of the time the SVMHC scoring function prefers one single very dominant amino acid per position, the algorithm shows almost never a clear "second place". One residue dominates and the remaining ones are more or less random. Interesting is that at position 3, the algorithm mixes an apolar (W) residue with a charged one (H). Despite this fact, there arises a pattern, because both amino acids carry aromatic residues that seem to be important for position 3 according to SVMHC. Also the other algorithms, with the exception of SYFPEITHI, support this preference. The consensus sequence of SVMHC matches the experimental anchor residues. However, V at position 6 is not present. The Van-der-Waals volume of the consensus sequence is 1 017 Å3.

XSCORE

XSCORE shows a strong tendency to assert high scores to big, hydrophobic amino acids. The dominating residue at position 1, 2, 3, 5, and 9 is W, whereas at position 5 and 7 F and Y are preferred. The corresponding consensus sequence would be repetitive and the probability of occurrence in the genetic code would be very low. The consensus sequence yields a Van-der-Waals volume of 1 331 Å3. The theoretical maximum (poly W-peptide) for a nonamer would only be slightly above at 1 467 Å3. One explanation for this exceptional behaviour could be the fact that XSCORE is the only structural scoring function used. Results suggest that XSCORE tries to fill the MHC binding groove as tightly as possible.

Comparison of the optimized peptides with public experimental data

One question which arises when having a look at the optimized peptides is: Are the predicted peptides related to training sets used by the respective scoring functions? While this is an interesting question, it is hard to give a definitive answer since even if the training data set is given in the publication of the scoring function, it is likely that the scoring function was re-calibrated with new data in the meantime since publication. For this reason we compared the resulting top 210 peptides (7×6×5) with the 14828 experimental peptide binding affinities available from the IEDB [54]. Interestingly, no peptide was identical. Only 2 peptides had a Levenshtein-distance of 2. Those 2 peptides are experimentally determined high affinity binders. All other peptides are more different from the public data. This would suggest that the scoring functions predict novel peptides by the use of the GA. However, as already mentioned above we cannot guarantee that the scoring functions were not trained on further "private" datasets.

Validation of the 5 consensus peptides via experimental binding assays

Since a validation on the basis of existing public experimental data was not possible we were able to find a co-operation partner who tested our consensus peptides for experimental peptide/MHC binding affinity. This was done at the revision level of the manuscript hence it is evident that the whole prediction was validated.

Although at first sight it seemed that the consensus peptides will be difficult to synthesise because of the multiple W and M residues, all peptides could be synthesised experimentally. All 5 consensus peptides are high affinity binders. The detailed results are shown in Table 1.

Table 1. Experimental validation of the consensus peptides.

Conclusion

We developed the framework "PeptX" for the evaluation of different operators for GAs in the context of peptide/MHC binding optimization. At the current time the focus of this study is mainly on the information theoretical part and may not be directly applicable in the wet-lab e.g. for questions like "what is the set of peptides (from a relevant pathogen), that bind an MHC with high affinity". However, in a next step further applications can be found in peptide immunotherapy: One wants to find high affinity binding peptides, but with certain constraints in the sequence to avoid allergic reactions (see Introduction).To be further able to address this issue we evaluated different parameter-sets of Genetic Algorithms in relation to MHC. On this basis further studies with direct relation to the wet-lab can be carried out.

The work most similar to our study was published by Wisniewska and co-workers. They combined an ant colony optimization algorithm with an artificial neural network classifier to iteratively adapt octapeptides for MHC class I stabilization [55]. However, to our knowledge our study is the first study which investigates the operators of a GA in relation to maximizing peptide/MHC binding affinity.

On the basis of our study we found two remarkable characteristics in the evolutionary process of the individuals (peptides). Firstly, it is intriguing that although selection operators have a strong influence on the convergence of the population while recombination operators have only minor influence, most reasonable operator combinations lead to convergence of the population. Long before generation 50 an optimum for most of the populations is found (see Figure 2). The fastest convergence is usually achieved by the best percent selection; however, the tournament selections often provide a slightly better optimized median score. Secondly, the most crucial step is the choice of a scoring function appropriate for the particular investigation. All other operators have minor influence on what is finally found optimal for binding to a specific MHC. Although there are similarities between the consensus sequences of the scoring functions (see Figure 3), one obtains different sets of optimal peptides by each scoring function. Hence, choosing a scoring function which is suboptimal for the purpose of the respective study renders almost all other parameters of the GA irrelevant. One should choose the employed scoring function with caution on the basis of previously published benchmarks (see Introduction). However, the GA itself can be utilized with an arbitrary scoring function and the convergence will mainly depend on the parameters evaluated in this study.

The PeptX framework is available for download including C++ source code for Linux at http://www.meduniwien.ac.at/msi/md/sourceCodes/peptX/peptX.htm webcite The download and usage is for free for academic researchers.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

BK designed the study, carried out programming work and drafted the manuscript. VG carried out programming work and helped to draft the manuscript. RR helped in the analysis of data and revised the manuscript critically. WS supervised the research and revised the manuscript critically. All authors read and approved the final manuscript.

Acknowledgements and Funding

The authors thank J. Sidney for testing the IC50 binding affinities of the consensus peptides.

This study was supported by the Austrian Science Fund (FWF P22258-B12) and the Austrian Grid Project of the Austrian Ministry of Education, Science and Culture (GZ 4003/2-VI/4c/2004).

References

  1. Rudolph MG, Stanfield RL, Wilson IA: How TCRs bind MHCs, peptides, and coreceptors.

    Annu Rev Immunol 2006, 24:419-466. OpenURL

  2. Tong JC, Tan TW, Ranganathan S: Methods and protocols for prediction of immunogenic epitopes.

    Brief Bioinform 2007, 8:96-108. OpenURL

  3. Dimitrov I, Garnev P, Flower DR, Doytchinova IA: MHC Class II Binding Prediction-A Little Help from a Friend.

    J Biomed Biotechnol 2010, 2010:705821. OpenURL

  4. Sieker F, May A, Zacharias M: Predicting affinity and specificity of antigenic peptide binding to major histocompatibility class I molecules.

    Curr Protein Pept Sci 2009, 10:286-296. OpenURL

  5. Saxova P, Buus S, Brunak S, Kesmir C: Predicting proteasomal cleavage sites: a comparison of available methods.

    Int Immunol 2003, 15:781-787. OpenURL

  6. Tong JC, Ren EC: Immunoinformatics: Current trends and future directions.

    Drug Discov Today 2009, 14:684-689. OpenURL

  7. Lafuente EM, Reche PA: Prediction of MHC-peptide binding: a systematic and comprehensive overview.

    Curr Pharm Des 2009, 15:3209-3220. OpenURL

  8. Flower DR, Macdonald IK, Ramakrishnan K, Davies MN, Doytchinova IA: Computer aided selection of candidate vaccine antigens.

    Immunome Res 2010, 6(Suppl 2):S1. OpenURL

  9. Lundegaard C, Lund O, Buus S, Nielsen M: Major histocompatibility complex class I binding predictions as a tool in epitope discovery.

    Immunology 2010, 130:309-318. OpenURL

  10. Gowthaman U, Agrewala JN: In silico tools for predicting peptides binding to HLA-class II molecules: more confusion than conclusion.

    J Proteome Res 2008, 7:154-163. OpenURL

  11. Wang P, Sidney J, Dow C, Mothe BR, Sette A, Peters B: A systematic assessment of MHC class II peptide binding predictions and evaluation of a consensus approach.

    PLoS Comput Biol 2008, 4:e1000048. OpenURL

  12. El-Manzalawy Y, Dobbs D, Honavar V: On evaluating MHC-II binding peptide prediction methods.

    PLoS ONE 2008, 3:e3268. OpenURL

  13. Lin HH, Zhang GL, Tongchusak S, Reinherz EL, Brusic V: Evaluation of MHC-II peptide binding prediction servers: applications for vaccine research.

    BMC Bioinformatics 2008, 9(Suppl 12):S22. OpenURL

  14. Tsurui H, Takahashi T: Prediction of T-cell epitope.

    J Pharmacol Sci 2007, 105:299-316. OpenURL

  15. Davies MN, Sansom CE, Beazley C, Moss DS: A novel predictive technique for the MHC class II peptide-binding interaction.

    Mol Med 2003, 9:220-225. OpenURL

  16. Knapp B, Omasits U, Bohle B, Maillere B, Ebner C, Schreiner W, Jahn-Schmid B: 3-Layer-based analysis of peptide-MHC-interaction: in silico prediction, peptide binding affinity and T cell activation in a relevant allergen-specific model.

    Molecular Immunology 2009, 46:1839-1844. OpenURL

  17. Knapp B, Omasits U, Schreiner W, Epstein MM: A comparative approach linking molecular dynamics of altered peptide ligands and MHC with in vivo immune responses.

    PLoS ONE 2010, 5:e11653. OpenURL

  18. Davies MN, Hattotuwagama CK, Moss DS, Drew MG, Flower DR: Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity.

    BMC Struct Biol 2006, 6:5. OpenURL

  19. Rognan D, Lauemoller SL, Holm A, Buus S, Tschinke V: Predicting binding affinities of protein ligands from three-dimensional models: application to peptide binding to class I major histocompatibility proteins.

    J Med Chem 1999, 42:4650-4658. OpenURL

  20. Logean A, Rognan D: Recovery of known T-cell epitopes by computational scanning of a viral genome.

    J Comput Aided Mol Des 2002, 16:229-243. OpenURL

  21. Schueler-Furman O, Altuvia Y, Sette A, Margalit H: Structure-based prediction of binding peptides to MHC class I molecules: application to a broad range of MHC alleles.

    Protein Sci 2000, 9:1838-1846. OpenURL

  22. Larche M: Update on the current status of peptide immunotherapy.

    J Allergy Clin Immunol 2007, 119:906-909. OpenURL

  23. Larche M, Wraith DC: Peptide-based therapeutic vaccines for allergic and autoimmune diseases.

    Nat Med 2005, 11:S69-S76. OpenURL

  24. Larche M: Inhibition of human T-cell responses by allergen peptides.

    Immunology 2001, 104:377-382. OpenURL

  25. Mishra S, Sinha S: Immunoinformatics and modeling perspective of T cell epitope-based cancer immunotherapy: a holistic picture.

    J Biomol Struct Dyn 2009, 27:293-306. OpenURL

  26. Purcell AW, McCluskey J, Rossjohn J: More than one reason to rethink the use of peptides in vaccine design.

    Nat Rev Drug Discov 2007, 6:404-414. OpenURL

  27. Alexander J, Fikes J, Hoffman S, Franke E, Sacci J, Appella E, Chisari FV, Guidotti LG, Chesnut RW, Livingston B, Sette A: The optimization of helper T lymphocyte (HTL) function in vaccine development.

    Immunol Res 1998, 18:79-92. OpenURL

  28. Alexander J, Sidney J, Southwood S, Ruppert J, Oseroff C, Maewal A, Snoke K, Serra HM, Kubo RT, Sette A: Development of high potency universal DR-restricted helper epitopes by modification of high affinity DR-blocking peptides.

    Immunity 1994, 1:751-761. OpenURL

  29. Shang X, Wang L, Niu W, Meng G, Fu X, Ni B, Lin Z, Yang Z, Chen X, Wu Y: Rational optimization of tumor epitopes using in silico analysis-assisted substitution of TCR contact residues.

    Eur J Immunol 2009, 39:2248-2258. OpenURL

  30. Reche PA, Reinherz EL: PEPVAC: a web server for multi-epitope vaccine development based on the prediction of supertypic MHC ligands.

    Nucleic Acids Res 2005, 33:W138-W142. OpenURL

  31. Toussaint NC, Donnes P, Kohlbacher O: A mathematical framework for the selection of an optimal set of peptides for epitope-based vaccines.

    PLoS Comput Biol 2008, 4:e1000246. OpenURL

  32. Parker AS, Zheng W, Griswold KE, Bailey-Kellogg C: Optimization algorithms for functional deimmunization of therapeutic proteins.

    BMC Bioinformatics 2010, 11:180. OpenURL

  33. Lazar GA, Desjarlais JR, Jacinto J, Karki S, Hammond PW: A molecular immunology approach to antibody humanization and functional optimization.

    Mol Immunol 2007, 44:1986-1998. OpenURL

  34. Bhasin M, Raghava GP: Prediction of promiscuous and high-affinity mutated MHC binders.

    Hybrid Hybridomics 2003, 22:229-234. OpenURL

  35. Doytchinova IA, Walshe VA, Jones NA, Gloster SE, Borrow P, Flower DR: Coupling in silico and in vitro analysis of peptide-MHC binding: a bioinformatic approach enabling prediction of superbinding peptides and anchorless epitopes.

    J Immunol 2004, 172:7495-7502. OpenURL

  36. Walshe VA, Hattotuwagama CK, Doytchinova IA, Wong M, Macdonald IK, Mulder A, Claas FH, Pellegrino P, Turner J, Williams I, Turnbull EL, Borrow P, Flower DR: Integrating in silico and in vitro analysis of peptide binding affinity to HLA-Cw*0102: a bioinformatic approach to the prediction of new epitopes.

    PLoS ONE 2009, 4:e8095. OpenURL

  37. Guan P, Doytchinova IA, Walshe VA, Borrow P, Flower DR: Analysis of peptide-protein binding using amino acid descriptors: prediction and experimental verification for human histocompatibility complex HLA-A0201.

    J Med Chem 2005, 48:7418-7425. OpenURL

  38. Roomp K, Antes I, Lengauer T: Predicting MHC class I epitopes in large datasets.

    BMC Bioinformatics 2010, 11:90. OpenURL

  39. Peters B, Bui HH, Frankild S, Nielson M, Lundegaard C, Kostem E, Basch D, Lamberth K, Harndahl M, Fleri W, Wilson SS, Sidney J, Lund O, Buus S, Sette A: A community resource benchmarking predictions of peptide binding to MHC-I molecules.

    PLoS Comput Biol 2006, 2:e65. OpenURL

  40. Knapp B, Omasits U, Frantal S, Schreiner W: A critical cross-validation of high throughput structural binding prediction methods for pMHC.

    J Comput Aided Mol Des 2009, 23:301-307. OpenURL

  41. Lin HH, Ray S, Tongchusak S, Reinherz EL, Brusic V: Evaluation of MHC class I peptide binding prediction servers: applications for vaccine research.

    BMC Immunol 2008, 9:8. OpenURL

  42. Zhang Q, Wang P, Kim Y, Haste-Andersen P, Beaver J, Bourne PE, Bui HH, Buus S, Frankild S, Greenbaum J, Lund O, Lundegaard C, Nielsen M, Ponomarenko J, Sette A, Zhu Z, Peters B: Immune epitope database analysis resource (IEDB-AR).

    Nucleic Acids Res 2008, 36:W513-W518. OpenURL

  43. Rammensee HG, Bachmann J, Emmerich NP, Bachor OA, Stevanovic S: SYFPEITHI: database for MHC ligands and peptide motifs.

    Immunogenetics 1999, 50:213-219. OpenURL

  44. Donnes P, Elofsson A: Prediction of MHC class I binding peptides, using SVMHC.

    BMC Bioinformatics 2002, 3:25. OpenURL

  45. Wan J, Liu W, Xu Q, Ren Y, Flower DR, Li T: SVRMHC prediction server for MHC-binding peptides.

    BMC Bioinformatics 2006, 7:463. OpenURL

  46. Wang R, Lai L, Wang S: Further development and validation of empirical scoring functions for structure-based binding affinity prediction.

    Journal of Computer-Aided Molecular Design 2002, 16:11-26. OpenURL

  47. Xu J: Rapid side-chain prediction via tree decomposition.

    RECOMB 2005, 3500:423-439. OpenURL

  48. Knapp B, Omasits U, Schreiner W: Side chain substitution benchmark for peptide/MHC interaction.

    Protein Sci 2008, 17:977-982. OpenURL

  49. Baeck T: Selective Pressure in Evolutionary Algorithms: A Characterization of Selection Mechanisms. In Proceedings of the First IEEE Conference on Evolutionary Computation. IEEE Press; 1994:57-62. OpenURL

  50. Baker JE: Reducing Bias and Inefficiency in the Selection Algorithm.

    New Jersey: Proceedings of the Second International Conference on Genetic Algorithms and their Application 1987, 14-21. OpenURL

  51. Holland JH: Adaptation in Natural and Artificial Systems. Ann Arbor: University of Michigan Press; 1975.

  52. Ackley DH: A connectionist machine for genetic hillclimbing. Springer Verlag Gmbh; 1987.

  53. Falk K, Rotzschke O, Stevanovic S, Jung G, Rammensee HG: Allele-specific motifs revealed by sequencing of self-peptides eluted from MHC molecules.

    Nature 1991, 351:290-296. OpenURL

  54. Vita R, Zarebski L, Greenbaum JA, Emami H, Hoof I, Salimi N, Damle R, Sette A, Peters B: The immune epitope database 2.0.

    Nucleic Acids Res 2010, 38:D854-D862. OpenURL

  55. Wisniewska JM, Jager N, Freier A, Losch FO, Wiesmuller KH, Walden P, Wrede P, Schneider G, Hiss JA: MHC I stabilizing potential of computer-designed octapeptides.

    J Biomed Biotechnol 2010, 2010:396847. OpenURL

  56. Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator.

    Genome Res 2004, 14:1188-1190. OpenURL