Evaluation of different generic in silico methods for predicting HLA class I binding peptide vaccine candidates using a reverse approach

Rare variants and HLA haplotypes associated in patients with neuromyelitis optica spectrum disorders

Neuromyelitis optica spectrum disorders (NMOSD) are rare, debilitating autoimmune diseases of the central nervous system. Many NMOSD patients have antibodies to Aquaporin-4 (AQP4). Prior studies show associations of NMOSD with individual Human Leukocyte Antigen (HLA) alleles and with mutations in the complement pathway and potassium channels. HLA allele associations with NMOSD are inconsistent between populations, suggesting complex relationships between the identified alleles and risk of disease. We used a retrospective case-control approach to identify contributing genetic variants in patients who met the diagnostic criteria for NMOSD and their unaffected family members. Potentially deleterious variants identified in NMOSD patients were compared to members of their families who do not have the disease and to existing databases of human genetic variation. HLA sequences from patients from Belgrade, Serbia, were compared to the frequency of HLA haplotypes in the general population in Belgrade. We analyzed exome sequencing on 40 NMOSD patients and identified rare inherited variants in the complement pathway and potassium channel genes. Haplotype analysis further detected two haplotypes, HLA-A*01, B*08, DRB1*03 and HLA-A*01, B*08, C*07, DRB1*03, DQB1*02, which were more prevalent in NMOSD patients than in unaffected individuals. In silico modeling indicates that HLA molecules within these haplotypes are predicted to bind AQP4 at several sites, potentially contributing to the development of autoimmunity. Our results point to possible autoimmune and neurodegenerative mechanisms that cause NMOSD, and can be used to investigate potential NMOSD drug targets.

Cancer systems immunology

Tumor immunology is undergoing a renaissance due to the recent profound clinical successes of tumor immunotherapy. These advances have coincided with an exponential growth in the development of -omics technologies. Armed with these technologies and their associated computational and modeling toolsets, systems biologists have turned their attention to tumor immunology in an effort to understand the precise nature and consequences of interactions between tumors and the immune system. Such interactions are inherently multivariate, spanning multiple time and size scales, cell types, and organ systems, rendering systems biology approaches particularly amenable to their interrogation. While in its infancy, the field of 'Cancer Systems Immunology' has already influenced our understanding of tumor immunology and immunotherapy. As the field matures, studies will move beyond descriptive characterizations toward functional investigations of the emergent behavior that govern tumor-immune responses. Thus, Cancer Systems Immunology holds incredible promise to advance our ability to fight this disease.

GradDock: rapid simulation and tailored ranking functions for peptide-MHC Class I docking

Motivation

The identification of T-cell epitopes has many profound translational applications in the areas of transplantation, disease diagnosis, vaccine/therapeutic protein development and personalized immunotherapy. While data-driven methods have been widely used for the prediction of peptide binders with notable successes, the structural modeling of peptide binding to MHC molecules is crucial for understanding the underlying molecular mechanism of the immunological processes.

Results

We developed GradDock, a structure-based method for the rapid and accurate modeling of peptide binding to MHC Class I (pMHC-I). GradDock explicitly models diverse unbound peptides in vacuo and inserts them into the MHC-I groove through a steered gradient descent with a topological correction process. The simulation process yields diverse structural conformations including native-like peptides. We completely revised the Rosetta score terms and developed a new ranking function specifically for pMHC-I. Using the diverse peptides, a linear programming approach is applied to find the optimal weights for the individual Rosetta score terms. Our examination revealed that a refinement of the dihedral angles and a modification of the repulsion can dramatically improve the modeling quality. GradDock is five-times faster than a Rosetta-based docking approach for pMHC-I. We also demonstrate that the predictive capability of GradDock with the re-weighted Rosetta ranking function is consistently more accurate than the Rosetta-based method with the standard Rosetta score (approximately three-times better for a cross-docking set).

Availability and implementation

GradDock is freely available for academic purposes. The program and the ranking score weights for Rosetta are available at http://bel.kaist.ac.kr/research/GradDock.

Contact

hskim76@kaist.ac.kr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Dynamics of T Cells Repertoire During Trypanosoma cruzi Infection and its Post-Treatment Modulation

Chagas disease courses with different clinical phases and has a variable clinical presentation and progression. The acute infection phase mostly exhibits a non-specific symptomatology. In the absence of treatment, the acute phase is followed by a chronic phase, which is initially asymptomatic. This chronic asymptomatic phase of the disease is characterized by a fragile balance between the host's immune response and the parasite replication. The loss of this balance is crucial for the progression of the sickness. The virulence and tropism of the T. cruzi infecting strain together to the inflammation processes in the cardiac tissue are the main factors for the establishment and severity of the cardiomyopathy. The efficacy of treatment in chronic Chagas disease patients is controversial. However, several studies carried out in chronic patients demonstrated that antiparasitic treatment reduces parasite load in the bloodstream and leads to an improvement in the immune response against the Trypanosoma cruzi parasite. The present review is mainly focused on the cellular patterns associated to the clinical status and the evolution of the disease in chronic patients, as well as the effectiveness of the treatment related to T. cruzi infection control. Therefore, an emphasis is placed on the dynamics of specific-antigens T cell subpopulations, their memory and activation phenotypes, their functionality and their contribution to pathogenesis or disease control, as well as their association with risk of congenital transmission of the parasite.

Differential phenotypic and functional profiles of TcCA-2 -specific cytotoxic CD8+ T cells in the asymptomatic versus cardiac phase in Chagasic patients

It has been reported that the immune response mediated by T CD8⁺ lymphocytes plays a critical role in the control of Trypanosoma cruzi infection and that the clinical symptoms of Chagas disease appear to be related to the competence of the CD8⁺ T immune response against the parasite. Herewith, in silico prediction and binding assays on TAP-deficient T2 cells were used to identify potential HLA-A*02:01 ligands in the T. cruzi TcCA-2 protein. The TcCA-2-specific CD8⁺ T cells were functionality evaluated by Granzyme B and cytokine production in peripheral blood mononuclear cells (PBMC) from Chagas disease patients stimulated with the identified HLA-A*02:01 peptides. The specific cells were phenotypically characterized by flow cytometry using several surface markers and HLA-A*02:01 APC-labeled dextramer loaded with the peptides. In the T. cruzi TcCA-2 protein four T CD8⁺ epitopes were identified which are processed and presented during Chagas disease. Interestingly, a differential cellular phenotypic profile could be correlated with the severity of the disease. The TcCA-2-specific T CD8⁺ cells from patients with cardiac symptoms are mainly effector memory cells (T_EM and T_EMRA) while, those present in the asymptomatic phase are predominantly naive cells (T_NAIVE). Moreover, in patients with cardiac symptoms the percentage of cells with senescence features is significantly higher than in patients at the asymptomatic phase of the disease. We consider that the identification of these new class I-restricted epitopes are helpful for designing biomarkers of sickness pathology as well as the development of immunotherapies against T. cruzi infection.

Bioinformatic and empirical analysis of novel hypoxia-inducible targets of the human antituberculosis T cell response

We analyzed whole genome-based transcriptional profiles of Mycobacterium tuberculosis subjected to prolonged hypoxia to guide the discovery of novel potential Ags, by a combined bioinformatic and empirical approach. We analyzed the fold induction of the 100 most highly induced genes at 7 d of hypoxia, as well as transcript abundance, peptide-binding prediction (ProPred) adjusted for population-specific MHC class II allele frequency, and by literature search. Twenty-six candidate genes were selected by this bioinformatic approach and evaluated empirically using IFN-γ and IL-2 ELISPOT using immunodominant Ags (Acr-1, CFP-10, ESAT-6) as references. Twenty-three of twenty-six proteins induced an IFN-γ response in PBMCs of persons with active or latent tuberculosis. Five novel immunodominant proteins-Rv1957, Rv1954c, Rv1955, Rv2022c, and Rv1471-were identified that induced responses similar to CFP-10 and ESAT-6 in both magnitude and frequency. IL-2 responses were of lower magnitude than were those of IFN-γ. Only moderate evidence of infection stage-specific recognition of Ags was observed. Reconciliation of bioinformatic and empirical hierarchies of immunodominance revealed that Ags could be predicted, providing transcriptomic data were combined with peptide-binding prediction adjusted by population-specific MHC class II allele frequency.

Trypanosoma cruzi paraflagellar rod proteins 2 and 3 contain immunodominant CD8(+) T-cell epitopes that are recognized by cytotoxic T cells from Chagas disease patients

The protozoan parasite Trypanosoma cruzi is the etiological agent of Chagas disease. To date, no vaccine is available for protection against T. cruzi infection. The CD8(+) T cells immune response against specific antigens has shown to efficiently control the spread of the parasite in murine experimental infection. However, data concerning CD8(+) response in Chagas patients are still restricted to a few epitopes. We have studied the existence of immunodominant CD8(+) T cell epitopes in the paraflagellar rod proteins 2 and 3 (PFR2 and PFR3) from T. cruzi in a mouse model, and analyzed their recognition by cytotoxic T lymphocytes from Chagas disease patients. Immunization of C57BL/6-A2/K(b) transgenic mice with plasmids coding for the fusion proteins PFR2-HSP70 and PFR3-HSP70 induced a specific CTL response against two PFRs epitopes (PFR2(449-457) and PFR3(481-489)), and showed specific lysis percentages of 24 and 12, respectively. Moreover, the PFR2(19-28), PFR2(156-163), PFR2(449-457), PFR3(428-436), PFR3(475-482) and PFR3(481-489) peptides were observed to have a high binding affinity to the HLA-A*02:01 molecule. Remarkably, these HLA-A*02:01-binding peptides are successfully processed and presented during natural infection by T. cruzi in the context of MHC class I as evidenced by using peptide-pulsed K562-A2 cells as antigen presenting cells. The T cells from Chagas disease chronic patients specific for PFR2/PFR3 selected CD8(+) epitopes showed a pro-inflammatory cytokine secretion profile (IFN-γ, TNF-α and IL-6). A positive Granzime B secretion was observed in three out of 16 patients in response to PFR2(156-163) and PFR2(449-457) peptides, two out of 11 patients in response to PFR2(19-28) peptide and one out of 14 and 11 patients in response to PFR3(428-436) and PFR3(481-489) peptides, respectively. The PFRs-specific cytotoxic activity in purified PBMC was only detected in patients in the indeterminate phase of the disease.

Potential T cell epitopes of Mycobacterium tuberculosis that can instigate molecular mimicry against host: implications in autoimmune pathogenesis

Background

Molecular mimicry between microbial antigens and host-proteins is one of the etiological enigmas for the occurrence of autoimmune diseases. T cells that recognize cross-reactive epitopes may trigger autoimmune reactions. Intriguingly, autoimmune diseases have been reported to be prevalent in tuberculosis endemic populations. Further, association of Mycobacterium tuberculosis (M. tuberculosis) has been implicated in different autoimmune diseases, including rheumatoid arthritis and multiple sclerosis. Although, in silico analyses have identified a number of M. tuberculosis specific vaccine candidates, the analysis on prospective cross-reactive epitopes, that may elicit autoimmune response, has not been yet attempted. Here, we have employed bioinformatics tools to determine T cell epitopes of homologous antigenic regions between M. tuberculosis and human proteomes.

Results

Employing bioinformatics tools, we have identified potentially cross-reactive T cell epitopes restricted to predominant class I and II alleles of human leukocyte antigens (HLA). These are similar to peptides of mycobacterial proteins and considerable numbers of them are promiscuous. Some of the identified antigens corroborated with established autoimmune diseases linked with mycobacterial infection.

Conclusions

The present study reveals many target proteins and their putative T cell epitopes that might have significant application in understanding the molecular basis of possible T cell autoimmune reactions during M. tuberculosis infections.

Identification of HLA-A∗02:01-restricted CTL epitopes in Trypanosoma cruzi heat shock protein-70 recognized by Chagas disease patients

CD8(+) cytotoxic T lymphocyte (CTL) response is critical for controlling the infection of the protozoan parasite Trypanosoma cruzi, the causative agent of Chagas disease. Since only a few CD8 antigens have been described in Chagas disease patients, the identification of new class I-restricted epitopes is urgently needed for the development of immunotherapies against T. cruzi infection. In this study, bioinformatic methods were used to predict HLA-A∗02:01-binders, and 30 peptides were selected, synthesized and tested for HLA-A∗02:01 binding. Among them, sixteen peptides with medium-to-high affinity were assayed for their recognition by CTL from HSP70-immunized or T. cruzi-infected transgenic B6-A2/K(b) mice. Our results show that four immunodominant epitopes (HSP70(210-8), HSP70(255-63), HSP70(316-24) and HSP70(345-53)) are contained in the T. cruzi HSP70 antigen. Indeed two of them (HSP70(210-8) and HSP70(316-24)) were also recognized by CTL of HLA-A∗02:01(+) Chagas disease patients, indicating that these peptides are processed and displayed as MHC class I epitopes during the natural history of T. cruzi infection. The HLA-A∗02:01 restriction was evidenced using peptide-pulsed K562-A2 cells as antigen-presenting cells. Both cytotoxic and cytokine-secreting activities were detected in response to the former two peptides and, moreover, 10/12 patients (83%) recognized at least one of these two HSP70-derived CD8(+) epitopes.

In silico identification of novel protective VSG antigens expressed by Trypanosoma brucei and an effort for designing a highly immunogenic DNA vaccine using IL-12 as adjuvant

African trypanosomiasis continues to be a major health problem, with more adults dying from this disease world-wide. As the sequence diversity of Trypanosoma brucei is extreme, with VSGs having 15-25% identity with most other VSGs, hence it displays a huge diversity of adaptations and host specificities. Therefore the need for an improved vaccine has become an international priority. The highly conserved and specific epitopes acting as both CD8+ and CD4+ T-cell epitopes (FLINKKPAL and FTALCTLAA) were predicted from large bunch of VSGs of T. brucei. Besides, some other potential epitopes with very high affinity for MHC I and II molecules were also determined while taking consideration on the most common HLA in the general population which accounts for major ethnicities. The vaccine candidates were found to be effective even for non-african populations as predicted by population coverage analysis. Hence the migrating travelers acting as a spread means of the infection can probably also be treated successfully after injection of such a multiepitopic vaccine. Exploiting the immunoinformatics approaches, we designed a potential vaccine by using the consensus epitopic sequence of 388 VSG proteins of T. brucei and performed in silico cloning of multiepitopic antigenic DNA sequence in pBI-CMV1 vector. Moreover, various techniques like codon adaptation, CpG optimization, removal of self recognized epitopes, use of adjuvant and co-injection with plasmids expressing immune-stimulatory molecules were implemented to enhance the immunogenicity of the proposed in silico vaccine.

MULTIPRED2: a computational system for large-scale identification of peptides predicted to bind to HLA supertypes and alleles

MULTIPRED2 is a computational system for facile prediction of peptide binding to multiple alleles belonging to human leukocyte antigen (HLA) class I and class II DR molecules. It enables prediction of peptide binding to products of individual HLA alleles, combination of alleles, or HLA supertypes. NetMHCpan and NetMHCIIpan are used as prediction engines. The 13 HLA Class I supertypes are A1, A2, A3, A24, B7, B8, B27, B44, B58, B62, C1, and C4. The 13 HLA Class II DR supertypes are DR1, DR3, DR4, DR6, DR7, DR8, DR9, DR11, DR12, DR13, DR14, DR15, and DR16. In total, MULTIPRED2 enables prediction of peptide binding to 1077 variants representing 26 HLA supertypes. MULTIPRED2 has visualization modules for mapping promiscuous T-cell epitopes as well as those regions of high target concentration – referred to as T-cell epitope hotspots. Novel graphic representations are employed to display the predicted binding peptides and immunological hotspots in an intuitive manner and also to provide a global view of results as heat maps. Another function of MULTIPRED2, which has direct relevance to vaccine design, is the calculation of population coverage. Currently it calculates population coverage in five major groups in North America. MULTIPRED2 is an important tool to complement wet-lab experimental methods for identification of T-cell epitopes. It is available at http://cvc.dfci.harvard.edu/multipred2/.