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.

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

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

Immuno-potentiating activity of surfactant vesicles in DNA and subunit vaccination

Enhanced immune responses for DNA and subunit vaccines potentiated by surfactant vesicle based delivery systems outlined in the present study, provides proof of principle for the beneficial aspects of vesicle mediated vaccination. The dehydration-rehydration technique was used to entrap plasmid DNA or subunit antigens into lipid-based (liposomes) or non-ionic surfactant-based (niosomes) dehydration-rehydration vesicles (DRV). Using this procedure, it was shown that both these types of antigens can be effectively entrapped in DRV liposomes and DRV niosomes. The vesicle size of DRV niosomes was shown to be twice the diameter (~2µm) of that of their liposome counterparts. Incorporation of cryoprotectants such as sucrose in the DRV procedure resulted in reduced vesicle sizes while retaining high DNA incorporation efficiency (~95%). Transfection studies in COS 7 cells demonstrated that the choice of cationic lipid, the helper lipid, and the method of preparation, all influenced transfection efficiency indicating a strong interdependency of these factors. This phenomenon has been further reinforced when 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE): cholesteryl 3b- [N-(N’ ,N’ -dimethylaminoethane)-carbamoyl] cholesterol (DC-Chol)/DNA complexes were supplemented with non-ionic surfactants. Morphological analysis of these complexes using transmission electron microscopy and environmental scanning electron microscopy (ESEM) revealed the presence of heterogeneous structures which may be essential for an efficient transfection in addition to the fusogenic properties of DOPE. In vivo evaluation of these DNA incorporated vesicle systems in BALB/c mice showed weak antibody and cell-mediated immune (CMI) responses. Subsequent mock challenge with hepatitis B antigen demonstrated that, 1-monopalmitoyl glycerol (MP) based DRV, is a more promising DNA vaccine adjuvant. Studying these DRV systems as adjuvants for the Hepatitis B subunit antigen (HBsAg) revealed a balanced antibody/CMI response profile on the basis of the HBsAg specific antibody and cytokine responses which were higher than unadjuvated antigen. The effect of addition of MP, cholesterol and trehalose 6,6’-dibehenate (TDB) on the stability and immuno-efficacy of dimethyldioctadecylammonium bromide (DDA) vesicles was investigated. Differential scanning calorimetry showed a reduction in transition temperature of DDA vesicles by ~12°C when incorporated with surfactants. ESEM of MP based DRV system indicated an increased vesicle stability upon incorporation of antigen. Adjuvant activity of these systems tested in C57BL/6j mice against three subunit antigens i.e., mycobacterial fusion protein- Ag85B-ESAT-6, and two malarial antigens - merozoite surface protein-1, (MSP1), and glutamate rich protein, (GLURP) revealed that while MP and DDA based systems induced comparable antibody responses, DDA based systems induced powerful CMI responses.

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.

Class I T-cell epitope prediction:improvements using a combination of proteasome cleavage, TAP affinity, and MHC binding

Cleavage by the proteasome is responsible for generating the C terminus of T-cell epitopes. Modeling the process of proteasome cleavage as part of a multi-step algorithm for T-cell epitope prediction will reduce the number of non-binders and increase the overall accuracy of the predictive algorithm. Quantitative matrix-based models for prediction of the proteasome cleavage sites in a protein were developed using a training set of 489 naturally processed T-cell epitopes (nonamer peptides) associated with HLA-A and HLA-B molecules. The models were validated using an external test set of 227 T-cell epitopes. The performance of the models was good, identifying 76% of the C-termini correctly. The best model of proteasome cleavage was incorporated as the first step in a three-step algorithm for T-cell epitope prediction, where subsequent steps predicted TAP affinity and MHC binding using previously derived models.

Comparison of vesicle based antigen delivery systems for delivery of hepatitis B surface antigen

There is a clinical need for a more effective vaccine against hepatitis B, and in particular vaccines that may be suitable for therapeutic administration. This study assesses the potential of cationic surfactant vesicle based formulations using two agents; the cationic amine containing [N-(N′,N′-dimethylaminoethane)-carbamoyl] cholesterol (DC-Chol) or dimethyl dioctadecylammonium bromide (DDA) with hepatitis B surface antigen (HBsAg). Synthetic mycobacterial cord factor, trehalose 6,6′-dibehenate (TDB) has been used as an adjuvant and the addition of 1-monopalmitoyl glycerol (C16:0) (MP) and cholesterol (Chol) to DDA-TDB is assessed for its potential to facilitate formation of dehydration-rehydration vesicles (DRV) at room temperature, and the effect of this on immune responses. A DRV formulation is directly compared to an adsorbed formulation of the same composition and preparation protocol (MP:dioleoyl phosphoethanolamine (DOPE):Chol:DC-Chol) and the direct substitution of MP with phosphatidylcholine (PC) is also compared in DRV antigen-entrapped formulations. MP and Chol were shown to facilitate the use of DDA-TDB in DRV formulations prepared at room temperature, whilst there was marginal alteration of immunogenicity (a reduction in HBsAg-specific IL-2). The HBsAg adsorbed DRV formulation was not significantly different from the HBsAg entrapped DRV formulation. Overall, DDA formulations incorporating TDB showed markedly increased antigen specific splenocyte proliferation and elicited cytokine production concomitant with a strong T cell driven response, delineating formulations that may be useful for further evaluation of their clinical potential. © 2007 Elsevier B.V. All rights reserved.

Novel formulations for antigen delivery using biodegradable polymers: new approaches for the use of new and established adjuvants

The use of immunological adjuvants has been established since 1924 and ever since many candidates have been extensively researched in vaccine development. The controlled release of vaccine is another area of biotechnology research, which is advancing rapidly with great potential and success. Encapsulation of peptide and protein drugs within biodegradable microspheres has been amongst the most successful of approaches within the past decade. The present studies have focused on combining the advantages of microsphere delivery systems composed of biodegradable polylactide (PLLA) and polylactide-co-glycolide (PLGA) polymers with that of safe and effective adjuvants. The research efforts were directed to the development of single-dose delivery vehicles which, can be manufactured easily, safely, under mild and favourable conditions to the encapsulated antigens. In pursuing this objective non ionic block copolymers (NIBCs) (Pluronics@ LI01 and L121) were incorporated within poly-dl-lactide (PDLA) micorospheres prepared with emulsification-diffusion method. LI0I and L121 served both as adjuvants and stabilising agents within these vaccine delivery vehicles. These formulations encapsulating the model antigens lysozyme, ovalbumin (OVA) and diphtheria toxoid (DT) resulted in high entrapment efficiency (99%), yield (96.7%) and elicited high and sustained immune response (IgG titres up to 9427) after one single administration over nine months. The structural integrity of the antigens was preserved within these formulations. In evaluating new approaches for the use of well-established adjuvants such as alum, these particles were incorporated within PLLA and PLGA microspheres at much lesser quantities (5-10 times lower) than those contained within conventional alum-adsorbed vaccines. These studies focused on the incorporation of the clinically relevant tetanus toxoid (TT) antigen within biodegradable microspheres. The encapsulation of both alum particles and TT antigen within these micropheres resulted in preparations with high encapsulation efficiency (95%) and yield (91.2%). The immune response to these particles was also investigated to evaluate the secretion of serum IgG, IgG1, IgG2a and IgG2b after a single administration of these vaccines. The Splenic cells proliferation was also investigated as an indication for the induction of cell mediated immunity. These particles resulted in high and sustained immune response over a period of 14 months. The stability of TT within particles was also investigated under dry storage over a period of several months. NIBC microspheres were also investigated as potential DNA vaccine delivery systems using hepatitis B plasmid. These particles resulted in micro spheres of 3-5 μm diameter and were shown to preserve the integrity of the encapsulated (27.7% entrapment efficiency) hepatitis B plasmid.

Identification of the HLA-DM/HLA-DR interface

Human leukocyte antigen (HLA)-DM is a critical participant in antigen presentation that catalyzes the dissociation of the Class II-associated Invariant chain-derived Peptide (CLIP) from the major histocompatibility complex (MHC) Class II molecules. There is competition amongst peptides for access to an MHC Class II groove and it has been hypothesised that DM functions as a 'peptide editor' that catalyzes the replacement of one peptide for another within the groove. It is established that the DM catalyst interacts directly with the MHC Class II but the precise location of the interface is unknown. Here, we combine previously described mutational data with molecular docking and energy minimisation simulations to identify a putative interaction site of >4000A2 which agrees with known point mutational data for both the DR and DM molecule. The docked structure is validated by comparison with experimental data and previously determined properties of protein-protein interfaces. A possible dissociation mechanism is suggested by the presence of an acidic cluster near the N terminus of the bound peptide.

The classification of HLA supertypes by GRID/CPCA and hierarchical clustering methods

Biological experiments often produce enormous amount of data, which are usually analyzed by data clustering. Cluster analysis refers to statistical methods that are used to assign data with similar properties into several smaller, more meaningful groups. Two commonly used clustering techniques are introduced in the following section: principal component analysis (PCA) and hierarchical clustering. PCA calculates the variance between variables and groups them into a few uncorrelated groups or principal components (PCs) that are orthogonal to each other. Hierarchical clustering is carried out by separating data into many clusters and merging similar clusters together. Here, we use an example of human leukocyte antigen (HLA) supertype classification to demonstrate the usage of the two methods. Two programs, Generating Optimal Linear Partial Least Square Estimations (GOLPE) and Sybyl, are used for PCA and hierarchical clustering, respectively. However, the reader should bear in mind that the methods have been incorporated into other software as well, such as SIMCA, statistiXL, and R.

Transporter associated with antigen processing preselection of peptides binding to the MHC:a bioinformatic evaluation

TAP is responsible for the transit of peptides from the cytosol to the lumen of the endoplasmic reticulum. In an immunological context, this event is followed by the binding of peptides to MHC molecules before export to the cell surface and recognition by T cells. Because TAP transport precedes MHC binding, TAP preferences may make a significant contribution to epitope selection. To assess the impact of this preselection, we have developed a scoring function for TAP affinity prediction using the additive method, have used it to analyze and extend the TAP binding motif, and have evaluated how well this model acts as a preselection step in predicting MHC binding peptides. To distinguish between MHC alleles that are exclusively dependent on TAP and those exhibiting only a partial dependence on TAP, two sets of MHC binding peptides were examined: HLA-A*0201 was selected as a representative of partially TAP-dependent HLA alleles, and HLA-A*0301 represented fully TAP-dependent HLA alleles. TAP preselection has a greater impact on TAP-dependent alleles than on TAP-independent alleles. The reduction in the number of nonbinders varied from 10% (TAP-independent) to 33% (TAP-dependent), suggesting that TAP preselection is an important component in the successful in silico prediction of T cell epitopes.

The administration route is decisive for the ability of the vaccine adjuvant CAF09 to induce antigen-specific CD8+ T-cell responses:the immunological consequences of the biodistribution profile

A prerequisite for vaccine-mediated induction of CD8+ T-cell responses is the targeting of dendritic cell (DC) subsets specifically capable of cross-presenting antigen epitopes to CD8+ T cells. Administration of a number of cationic adjuvants via the intraperitoneal (i.p.) route has been shown to result in strong CD8+ T-cell responses, whereas immunization via e.g. the intramuscular (i.m.) or subcutaneous (s.c.) routes often stimulate weak CD8+ T-cell responses. The hypothesis for this is that self-drainage of the adjuvant/antigen to the lymphoid organs, which takes place upon i.p. immunization, is required for the subsequent activation of cross-presenting lymphoid organ-resident CD8α+ DCs. In contrast, s.c. or i.m. immunization usually results in the formation of a depot at the site of injection (SOI), which hinders the self-drainage and targeting of the vaccine to cross-presenting CD8α+ DCs. We investigated this hypothesis by correlating the biodistribution pattern and the adjuvanticity of the strong CD8+ T-cell inducing liposomal cationic adjuvant formulation 09 (CAF09), which is composed of dimethyldioctadecylammonium bromide/monomycoloyl glycerol liposomes with polyinosinic:polycytidylic acid electrostatically adsorbed to the surface. Biodistribution studies with radiolabeled CAF09 and a surface-adsorbed model antigen [ovalbumin (OVA)] showed that a significantly larger fraction of the vaccine dose localized in the draining lymph nodes (dLNs) and the spleen 6 h after i.p. immunization, as compared to after i.m. immunization. Studies with fluorescently labelled OVA + CAF09 demonstrated a preferential association of OVA + CAF09 to DCs/monocytes, as compared to macrophages and B cells, following i.p. immunization. Administration of OVA + CAF09 via the i.p. route did also result in DC activation, whereas no DC activation could be measured within the same period with unadjuvanted OVA and OVA + CAF09 administered via the s.c. or i.m. routes. In the dLNs, the highest level of activated, cross-presenting CD8α+ DCs was detected at 24 h post immunization, whereas an influx of activated, migrating and cross-presenting CD103+ DCs to the dLNs could be measured after 48 h. This suggests that the CD8α+ DCs are activated by self-draining OVA + CAF09 in the lymphoid organs, whereas the CD103+ DCs are stimulated by the OVA + CAF09 at the SOI. These results support the hypothesis that the self-drainage of OVA + CAF09 to the draining LNs is required for the activation of CD8α+ DCs, while the migratory CD103+ DCs may play a role in sustaining the subsequent induction of strong CD8+ T-cell responses.