Bioinformatic Techniques for Vaccine Development: Epitope Prediction and Structural Vaccinology

Structural vaccinology involves characterizing the interactions between an antigen and antibodies or host immune receptors. Central to this is the task of epitope prediction, which involves describing the binding affinity and interactions of a given peptide typically to the major histocompatibility complex in the case of T-cells or to the antibodies in the case of B-cells. Several computational models exist for this purpose which we will review here. Generally, epitope predictions for MHC-I and MHC-II are substantially different tasks as well as epitope prediction for continuous versus discontinuous B-cell epitopes. Overall, these models suffer from overprediction of epitopes although general themes support both the use of neural networks as well as the incorporation of more abundant and more varied experimental annotation into model training as valuable in improving predictive performance.

In Silico Prediction of Human Leukocytes Antigen (HLA) Class II Binding Hepatitis B Virus (HBV) Peptides in Botswana

Hepatitis B virus (HBV) is the primary cause of liver-related malignancies worldwide, and there is no effective cure for chronic HBV infection (CHB) currently. Strong immunological responses induced by T cells are associated with HBV clearance during acute infection; however, the repertoire of epitopes (epi) presented by major histocompatibility complexes (MHCs) to elicit these responses in various African populations is not well understood. In silico approaches were used to map and investigate 15-mers HBV peptides restricted to 9 HLA class II alleles with high population coverage in Botswana. Sequences from 44 HBV genotype A and 48 genotype D surface genes (PreS/S) from Botswana were used. Of the 1819 epi bindings predicted, 20.2% were strong binders (SB), and none of the putative epi bind to all the 9 alleles suggesting that multi-epitope, genotype-based, population-based vaccines will be more effective against HBV infections as opposed to previously proposed broad potency epitope-vaccines which were assumed to work for all alleles. In total, there were 297 unique epi predicted from the 3 proteins and amongst, S regions had the highest number of epi (n = 186). Epitope-densities (Depi) between genotypes A and D were similar. A number of mutations that hindered HLA-peptide binding were observed. We also identified antigenic and genotype-specific peptides with characteristics that are well suited for the development of sensitive diagnostic kits. This study identified candidate peptides that can be used for developing multi-epitope vaccines and highly sensitive diagnostic kits against HBV infection in an African population. Our results suggest that viral variability may hinder HBV peptide-MHC binding, required to initiate a cascade of immunological responses against infection.

Software tools for simultaneous data visualization and T cell epitopes and disorder prediction in proteins

We have developed EpDis and MassPred, extendable open source software tools that support bioinformatic research and enable parallel use of different methods for the prediction of T cell epitopes, disorder and disordered binding regions and hydropathy calculation. These tools offer a semi-automated installation of chosen sets of external predictors and an interface allowing for easy application of the prediction methods, which can be applied either to individual proteins or to datasets of a large number of proteins. In addition to access to prediction methods, the tools also provide visualization of the obtained results, calculation of consensus from results of different methods, as well as import of experimental data and their comparison with results obtained with different predictors. The tools also offer a graphical user interface and the possibility to store data and the results obtained using all of the integrated methods in the relational database or flat file for further analysis. The MassPred part enables a massive parallel application of all integrated predictors to the set of proteins. Both tools can be downloaded from http://bioinfo.matf.bg.ac.rs/home/downloads.wafl?cat=Software. Appendix A includes the technical description of the created tools and a list of supported predictors.

Epitope distribution in ordered and disordered protein regions - part A. T-cell epitope frequency, affinity and hydropathy

Highly disordered protein regions are prevalently hydrophilic, extremely sensitive to proteolysis in vitro, and are expected to be under-represented as T-cell epitopes. The aim of this research was to find out whether disorder and hydropathy prediction methods could help in understanding epitope processing and presentation. According to the pan-specific T-cell epitope predictors NetMHCpan and NetMHCIIpan and 9 publicly available disorder predictors, frequency of epitopes presented by human leukocyte antigens (HLA) class-I or -II was found to be more than 2.5 times higher in ordered than in disordered protein regions (depending on the disorder predictor). Both HLA class-I and HLA class-II binding epitopes are prevalently hydrophilic in disordered and prevalently hydrophobic in ordered protein regions, whereas epitopes recognized by HLA class-II alleles are more hydrophobic than those recognized by HLA class-I. As regards both classes of HLA molecules, high-affinity binding epitopes display more hydrophobicity than low affinity-binding epitopes (in both ordered and disordered regions). Epitopes belonging to disordered protein regions were not predicted to have poor affinity to HLA class-II molecules, as expected from disorder intrinsic proteolytic instability. The relation of epitope hydrophobicity and order/disorder location was also valid if alleles were grouped according to the HLA class-I and HLA class-II supertypes, except for the class-I supertype A3 in which the main part of recognized epitopes was prevalently hydrophilic. Regarding specific supertypes, the affinity of epitopes belonging to ordered regions varies only slightly (depending on the disorder predictor) compared to the affinity of epitopes in corresponding disordered regions. The distribution of epitopes in ordered and disordered protein regions has revealed that the curves of order-epitope distribution were convex-like while the curves of disorder-epitope distribution were concave-like. The percentage of prevalently hydrophobic epitopes increases with the enhancement of epitope promiscuity level and moving from disordered to ordered regions. These data suggests that reverse vaccinology, oriented towards promiscuous and high-affinity epitopes, is also oriented towards prevalently hydrophobic, ordered regions. The analysis of predicted and experimentally evaluated epitopes of cancer-testis antigen MAGE-A3 has confirmed that the majority of T-cell epitopes, particularly those that are promiscuous or naturally processed, was located in ordered and disorder/order boundary protein regions overlapping hydrophobic regions.

Sequence-based in silico analysis of well studied hepatitis C virus epitopes and their variants in other genotypes (particularly genotype 5a) against South African human leukocyte antigen backgrounds

Background

Host genetics influence the outcome of HCV disease. HCV is also highly mutable and escapes host immunity. HCV genotypes are geographically distributed and HCV subtypes have been shown to have distinct repertoires of HLA-restricted viral epitopes which explains the lack of cross protection across genotypes observed in some studies. Despite this, immune databases and putative epitope vaccines concentrate almost exclusively on HCV genotype 1 class I-epitopes restricted by the HLA-A*02 allele. While both genotype and allele predominate in developed countries, we hypothesise that HCV variation and population genetics will affect the efficacy of proposed epitope vaccines in South Africa. This in silico study investigates HCV viral variability within well-studied epitopes identified in genotype 1 and uses algorithms to predict the immunogenicity of their variants from other less studied genotypes and thus rate the most promising vaccine candidates for the South African population. Six class I- and seven class II- restricted epitope sequences within the core, NS3, NS4B and NS5B regions were compared across the six HCV genotypes using local genotype 5a sequence data together with global data. Common HLA alleles in the South African population are A30:01, A02:01, B58:02, B07:02; DRB1*13:01 and DRB1*03:01. Epitope binding to 13 class I- and 8 class –II alleles were described using web-based prediction servers, Immune Epitope Database, (IEDB) and Propred. Online population coverage tools were used to assess vaccine efficacy.

Results

Despite the homogeneity of genotype 1 and genotype 5 over the epitopes, there was limited promiscuity to local HLA-alleles.Host differences will make a putative vaccine less effective in South Africa. Of the 6 well-characterized class I- epitopes, only 2 class I- epitopes were promiscuous and 3 of the 7 class-II epitopes were better conserved and promiscuous. By fine tuning the putative vaccine using an optimal cocktail of genotype 1 and 5a epitopes and local HLA data, the coverage was raised from 65.85% to 91.87% in South African Blacks.

Conclusion

While in vivo and in vitro studies are needed to confirm immunogenic epitopes, in silico HCV epitope vaccine design which takes into account HCV variation and host allele frequency will maximize population coverage in different ethnic groups.