In silico prediction of cancer immunogens: current state of the art.

Cancer kills 8 million annually worldwide. Although survival rates in prevalent cancers continue to increase, many cancers have no effective treatment, prompting the search for new and improved protocols. Immunotherapy is a new and exciting addition to the anti-cancer arsenal. The successful and accurate identification of aberrant host proteins acting as antigens for vaccination and immunotherapy is a key aspiration for both experimental and computational research. Here we describe key elements of in silico prediction, including databases of cancer antigens and bleeding-edge methodology for their prediction. We also highlight the role dendritic cell vaccines can play and how they can act as delivery mechanisms for epitope ensemble vaccines. Immunoinformatics can help streamline the discovery and utility of Cancer Immunogens.

Integrating in silico and in vitro analysis of peptide binding affinity to HLA-Cw*0102: a bioinformatic approach to the prediction of new epitopes.

BACKGROUND: Predictive models of peptide-Major Histocompatibility Complex (MHC) binding affinity are important components of modern computational immunovaccinology. Here, we describe the development and deployment of a reliable peptide-binding prediction method for a previously poorly-characterized human MHC class I allele, HLA-Cw*0102. METHODOLOGY/FINDINGS: Using an in-house, flow cytometry-based MHC stabilization assay we generated novel peptide binding data, from which we derived a precise two-dimensional quantitative structure-activity relationship (2D-QSAR) binding model. This allowed us to explore the peptide specificity of HLA-Cw*0102 molecule in detail. We used this model to design peptides optimized for HLA-Cw*0102-binding. Experimental analysis showed these peptides to have high binding affinities for the HLA-Cw*0102 molecule. As a functional validation of our approach, we also predicted HLA-Cw*0102-binding peptides within the HIV-1 genome, identifying a set of potent binding peptides. The most affine of these binding peptides was subsequently determined to be an epitope recognized in a subset of HLA-Cw*0102-positive individuals chronically infected with HIV-1. CONCLUSIONS/SIGNIFICANCE: A functionally-validated in silico-in vitro approach to the reliable and efficient prediction of peptide binding to a previously uncharacterized human MHC allele HLA-Cw*0102 was developed. This technique is generally applicable to all T cell epitope identification problems in immunology and vaccinology.

VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines.

BACKGROUND: Vaccine development in the post-genomic era often begins with the in silico screening of genome information, with the most probable protective antigens being predicted rather than requiring causative microorganisms to be grown. Despite the obvious advantages of this approach--such as speed and cost efficiency--its success remains dependent on the accuracy of antigen prediction. Most approaches use sequence alignment to identify antigens. This is problematic for several reasons. Some proteins lack obvious sequence similarity, although they may share similar structures and biological properties. The antigenicity of a sequence may be encoded in a subtle and recondite manner not amendable to direct identification by sequence alignment. The discovery of truly novel antigens will be frustrated by their lack of similarity to antigens of known provenance. To overcome the limitations of alignment-dependent methods, we propose a new alignment-free approach for antigen prediction, which is based on auto cross covariance (ACC) transformation of protein sequences into uniform vectors of principal amino acid properties. RESULTS: Bacterial, viral and tumour protein datasets were used to derive models for prediction of whole protein antigenicity. Every set consisted of 100 known antigens and 100 non-antigens. The derived models were tested by internal leave-one-out cross-validation and external validation using test sets. An additional five training sets for each class of antigens were used to test the stability of the discrimination between antigens and non-antigens. The models performed well in both validations showing prediction accuracy of 70% to 89%. The models were implemented in a server, which we call VaxiJen. CONCLUSION: VaxiJen is the first server for alignment-independent prediction of protective antigens. It was developed to allow antigen classification solely based on the physicochemical properties of proteins without recourse to sequence alignment. The server can be used on its own or in combination with alignment-based prediction methods. It is freely-available online at the URL: http://www.jenner.ac.uk/VaxiJen.

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.

MHCPred 2.0: an updated quantitative T-cell epitope prediction server.

UNLABELLED: The accurate computational prediction of T-cell epitopes can greatly reduce the experimental overhead implicit in candidate epitope identification within genomic sequences. In this article we present MHCPred 2.0, an enhanced version of our online, quantitative T-cell epitope prediction server. The previous version of MHCPred included mostly alleles from the human leukocyte antigen A (HLA-A) locus. In MHCPred 2.0, mouse models are added and computational constraints removed. Currently the server includes 11 human HLA class I, three human HLA class II, and three mouse class I models. Additionally, a binding model for the human transporter associated with antigen processing (TAP) is incorporated into the new MHCPred. A tool for the design of heteroclitic peptides is also included within the server. To refine the veracity of binding affinities prediction, a confidence percentage is also now calculated for each peptide predicted. AVAILABILITY: As previously, MHCPred 2.0 is freely available at the URL http://www.jenner.ac.uk/MHCPred/ CONTACT: Darren R. Flower (darren.flower@jenner.ac.uk).

Synthesis and antiproliferative activity of benzo[d]isothiazole hydrazones.

Several benzo[d]isothiazole hydrazones have been evaluated for their potential antiretroviral activity. Since a number of these compounds were found to be inactive against viruses, but showed cytotoxicity at micromolar concentrations against the human CD4+ lymphocytes (MT-4) that were used to support HIV-1 growth, they were further tested for antiproliferative activity. The compounds resulted as being cytotoxic for MT-4 cells and new derivatives which were rationally designed and synthesized, were tested for antiproliferative activity against several leukaemia and solid tumour cell lines. In addition, these compounds were evaluated against "normal" cell lines. Compound 2h proved to be the most active compound and the fragment -CO-NH-N=CH-2-hydroxyphenyl was identified as being very important for biological activity, suggesting intramolecular hydrogen bond formation or favourable mutual disposition between two important centres in the pharmacophore. 1H-NMR spectra have been explained with the support of a conformational analysis.

Modeling the peptide-T cell receptor interaction by the comparative molecular similarity indices analysis-soft independent modeling of class analogy technique.

A set of 38 epitopes and 183 non-epitopes, which bind to alleles of the HLA-A3 supertype, was subjected to a combination of comparative molecular similarity indices analysis (CoMSIA) and soft independent modeling of class analogy (SIMCA). During the process of T cell recognition, T cell receptors (TCR) interact with the central section of the bound nonamer peptide; thus only positions 4-8 were considered in the study. The derived model distinguished 82% of the epitopes and 73% of the non-epitopes after cross-validation in five groups. The overall preference from the model is for polar amino acids with high electron density and the ability to form hydrogen bonds. These so-called "aggressive" amino acids are flanked by small-sized residues, which enable such residues to protrude from the binding cleft and take an active role in TCR-mediated T cell recognition. Combinations of "aggressive" and "passive" amino acids in the middle part of epitopes constitute a putative TCR binding motif.

In silico prediction of peptide binding affinity to class I mouse major histocompatibility complexes: a comparative molecular similarity index analysis (CoMSIA) study.

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 prediction of peptide-major histocompatibility complex (MHC) affinity. A three-dimensional quantitative structure-activity relationship (3D-QSAR) for the prediction of peptide binding to class I MHC molecules was established using the comparative molecular similarity index analysis (CoMSIA) method. Three MHC alleles were studied: H2-D(b), H2-K(b), and H2-K(k). Models were produced for each allele. Each model consisted of five physicochemical descriptors-steric bulk, electrostatic potentials, hydrophobic interactions, and hydrogen-bond donor and hydrogen-bond acceptor abilities. The models have an acceptable level of predictivity: cross-validation leave-one-out statistical terms q2 and SEP (standard error of prediction) ranged between 0.490 and 0.679 and between 0.525 and 0.889, respectively. The non-cross-validated statistical terms r2 and SEE (standard error of estimate) ranged between 0.913 and 0.979 and between 0.167 and 0.248, respectively. The use of coefficient contour maps, which indicate favored and disfavored areas for each position of the MHC-bound peptides, allowed the binding specificity of each allele to be identified, visualized, and understood. The present study demonstrates the effectiveness of CoMSIA as a method for studying peptide-MHC interactions. The peptides used in this study are available on the Internet (http://www.jenner.ac.uk/AntiJen). The partial least-squares method is available commercially in the SYBYL molecular modeling software package.

New horizons in mouse immunoinformatics: reliable in silico prediction of mouse class I histocompatibility major complex peptide binding affinity.

Quantitative structure-activity relationship (QSAR) analysis is a main cornerstone of modern informatic disciplines. Predictive computational models, based on QSAR technology, of peptide-major histocompatibility complex (MHC) binding affinity have now become a vital component of modern day computational immunovaccinology. Historically, such approaches have been built around semi-qualitative, classification methods, but these are now giving way to quantitative regression methods. The additive method, an established immunoinformatics technique for the quantitative prediction of peptide-protein affinity, was used here to identify the sequence dependence of peptide binding specificity for three mouse class I MHC alleles: H2-D(b), H2-K(b) and H2-K(k). As we show, in terms of reliability the resulting models represent a significant advance on existing methods. They can be used for the accurate prediction of T-cell epitopes and are freely available online ( http://www.jenner.ac.uk/MHCPred).

Coupling in silico and in vitro analysis of peptide-MHC binding: a bioinformatic approach enabling prediction of superbinding peptides and anchorless epitopes.

The ability to define and manipulate the interaction of peptides with MHC molecules has immense immunological utility, with applications in epitope identification, vaccine design, and immunomodulation. However, the methods currently available for prediction of peptide-MHC binding are far from ideal. We recently described the application of a bioinformatic prediction method based on quantitative structure-affinity relationship methods to peptide-MHC binding. In this study we demonstrate the predictivity and utility of this approach. We determined the binding affinities of a set of 90 nonamer peptides for the MHC class I allele HLA-A*0201 using an in-house, FACS-based, MHC stabilization assay, and from these data we derived an additive quantitative structure-affinity relationship model for peptide interaction with the HLA-A*0201 molecule. Using this model we then designed a series of high affinity HLA-A2-binding peptides. Experimental analysis revealed that all these peptides showed high binding affinities to the HLA-A*0201 molecule, significantly higher than the highest previously recorded. In addition, by the use of systematic substitution at principal anchor positions 2 and 9, we showed that high binding peptides are tolerant to a wide range of nonpreferred amino acids. Our results support a model in which the affinity of peptide binding to MHC is determined by the interactions of amino acids at multiple positions with the MHC molecule and may be enhanced by enthalpic cooperativity between these component interactions.

Quantitative online prediction of peptide binding to the major histocompatibility complex.

With its implications for vaccine discovery, the accurate prediction of T cell epitopes is one of the key aspirations of computational vaccinology. We have developed a robust multivariate statistical method, based on partial least squares, for the quantitative prediction of peptide binding to major histocompatibility complexes (MHC), the principal checkpoint on the antigen presentation pathway. As a service to the immunobiology community, we have made a Perl implementation of the method available via a World Wide Web server. We call this server MHCPred. Access to the server is freely available from the URL: http://www.jenner.ac.uk/MHCPred. We have exemplified our method with a model for peptides binding to the common human MHC molecule HLA-B*3501.

JenPep: a novel computational information resource for immunobiology and vaccinology.

JenPep is a relational database containing a compendium of thermodynamic binding data for the interaction of peptides with a range of important immunological molecules: the major histocompatibility complex, TAP transporter, and T cell receptor. The database also includes annotated lists of B cell and T cell epitopes. Version 2.0 of the database is implemented in a bespoke postgreSQL database system and is fully searchable online via a perl/HTML interface (URL: http://www.jenner.ac.uk/JenPep).

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-based 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*0401 and HLA-DRB*0701). MHCPred is available from the URL: http://www.jenner.ac.uk/MHCPred.

A comparative molecular similarity indices (CoMSIA) study of peptide binding to the HLA-A3 superfamily.

Epitope identification is the basis of modern vaccine design. The present paper studied the supermotif of the HLA-A3 superfamily, using comparative molecular similarity indices analysis (CoMSIA). Four alleles with high phenotype frequencies were used: A*1101, A*0301, A*3101 and A*6801. Five physicochemical properties-steric bulk, electrostatic potential, local hydrophobicity, hydrogen-bond donor and acceptor abilities-were considered and 'all fields' models were produced for each of the alleles. The models have a moderate level of predictivity and there is a good correlation between the data. A revised HLA-A3 supermotif was defined based on the comparison of favoured and disfavoured properties for each position of the MHC bound peptide. The present study demonstrated that CoMSIA is an effective tool for studying peptide-MHC interactions.

MHCPred: bringing a quantitative dimension to the online prediction of MHC binding.

The accurate prediction of T cell epitopes is one of the key aspirations of immunoinformatics. We have developed a partial least squares-based, robust multivariate statistical method for the quantitative prediction of peptide binding to major histocompatibility complexes (MHCs), the principal checkpoint on the antigen presentation pathway. As a service to the immunobiology community, we have made a Perl implementation of the method available as a World Wide Web server.

Physicochemical explanation of peptide binding to HLA-A*0201 major histocompatibility complex: a three-dimensional quantitative structure-activity relationship study.

A three-dimensional quantitative structure-activity relationship method for the prediction of peptide binding affinities to the MHC class I molecule HLA-A*0201 was developed by applying the CoMSIA technique on a set of 266 peptides. To increase the self consistency of the initial CoMSIA model, the poorly predicted peptides were excluded from the training set in a stepwise manner and then included in the study as a test set. The final model, based on 236 peptides and considering the steric, electrostatic, hydrophobic, hydrogen bond donor, and hydrogen bond acceptor fields, had q2 = 0.683 and r2 = 0.891. The stability of this model was proven by cross-validations in two and five groups and by a bootstrap analysis of the non-cross-validated model. The residuals between the experimental pIC50 (-logIC50) values and those calculated by "leave-one-out" cross-validation were analyzed. According to the best model, 63.2% of the peptides were predicted with /residuals/ < or = 0.5 log unit; 29.3% with 1.0 < or = /residuals/ < 0.5; and 7.5% with /residuals/ > 1.0 log unit. The mean /residual/ value was 0.489. The coefficient contour maps identify the physicochemical property requirements at each position in the peptide molecule and suggest amino acid sequences for high-affinity binding to the HLA-A*0201 molecule.

Quantitative approaches to computational vaccinology.

This article reviews the newly released JenPep database and two new powerful techniques for T-cell epitope prediction: (i) the additive method; and (ii) a 3D-Quantitative Structure Activity Relationships (3D-QSAR) method, based on Comparative Molecular Similarity Indices Analysis (CoMSIA). The JenPep database is a family of relational databases supporting the growing need of immunoinformaticians for quantitative data on peptide binding to major histocompatibility complexes and to the Transporters associated with Antigen Processing (TAP). It also contains an annotated list of T-cell epitopes. The database is available free via the Internet (http://www.jenner.ac.uk/JenPep). The additive prediction method is based on the assumption that the binding affinity of a peptide depends on the contributions from each amino acid as well as on the interactions between the adjacent and every second side-chain. In the 3D-QSAR approach, the influence of five physicochemical properties (steric bulk, electrostatic potential, local hydrophobicity, hydrogen-bond donor and hydrogen-bond acceptor abilities) on the affinity of peptides binding to MHC molecules were considered. Both methods were exemplified through their application to the well-studied problem of peptides binding to the human class I MHC molecule HLA-A*0201.

JenPep: a database of quantitative functional peptide data for immunology.

MOTIVATION: The compilation of quantitative binding data underlies attempts to derive tools for the accurate prediction of epitopes in cellular immunology and is part of our concerted goal to develop practical computational vaccinology. RESULTS: JenPep is a family of relational databases supporting the growing community of immunoinformaticians. It contains quantitative data on peptide binding to Major Histocompatibility Complexes (MHCs) and to Transmembrane Peptide Transporter (TAP), as well as an annotated list of T-cell epitopes. AVAILABILITY: The database is available via the Internet. An HTML interface allowing searching of the database can be found at the following address: http://www.jenner.ac.uk/JenPep.

Additive method for the prediction of protein-peptide binding affinity. Application to the MHC class I molecule HLA-A*0201.

A method has been developed for prediction of binding affinities between proteins and peptides. We exemplify the method through its application to binding predictions of peptides with affinity to major histocompatibility complex class I molecule HLA-A*0201. The method is named "additive" because it is based on the assumption that the binding affinity of a peptide could be presented as a sum of the contributions of the amino acids at each position and the interactions between them. The amino acid contributions and the contributions of the interactions between adjacent side chains and every second side chain were derived using a partial least squares (PLS) statistical methodology using a training set of 420 experimental IC50 values. The predictive power of the method was assessed using rigorous cross-validation and using an independent test set of 89 peptides. The mean value of the residuals between the experimental and predicted pIC50 values was 0.508 for this test set. The additive method was implemented in a program for rapid T-cell epitope search. It is universal and can be applied to any peptide-protein interaction where binding data is known.