ProPred analysis and experimental evaluation of promiscuous T-cell epitopes of three major secreted antigens of Mycobacterium tuberculosis

Optimization algorithms for functional deimmunization of therapeutic proteins

Background

To develop protein therapeutics from exogenous sources, it is necessary to mitigate the risks of eliciting an anti-biotherapeutic immune response. A key aspect of the response is the recognition and surface display by antigen-presenting cells of epitopes, short peptide fragments derived from the foreign protein. Thus, developing minimal-epitope variants represents a powerful approach to deimmunizing protein therapeutics. Critically, mutations selected to reduce immunogenicity must not interfere with the protein's therapeutic activity.

Results

This paper develops methods to improve the likelihood of simultaneously reducing the anti-biotherapeutic immune response while maintaining therapeutic activity. A dynamic programming approach identifies optimal and near-optimal sets of conservative point mutations to minimize the occurrence of predicted T-cell epitopes in a target protein. In contrast with existing methods, those described here integrate analysis of immunogenicity and stability/activity, are broadly applicable to any protein class, guarantee global optimality, and provide sufficient flexibility for users to limit the total number of mutations and target MHC alleles of interest. The input is simply the primary amino acid sequence of the therapeutic candidate, although crystal structures and protein family sequence alignments may also be input when available. The output is a scored list of sets of point mutations predicted to reduce the protein's immunogenicity while maintaining structure and function. We demonstrate the effectiveness of our approach in a number of case study applications, showing that, in general, our best variants are predicted to be better than those produced by previous deimmunization efforts in terms of either immunogenicity or stability, or both factors.

Conclusions

By developing global optimization algorithms leveraging well-established immunogenicity and stability prediction techniques, we provide the protein engineer with a mechanism for exploring the favorable sequence space near a targeted protein therapeutic. Our mechanism not only helps identify designs more likely to be effective, but also provides insights into the interrelated implications of design choices.

Mycobacterial antigen-induced T helper type 1 (Th1) and Th2 reactivity of peripheral blood mononuclear cells from diabetic and non-diabetic tuberculosis patients and Mycobacterium bovis bacilli Calmette-Guérin (BCG)-vaccinated healthy subjects

Patients with diabetes mellitus are more susceptible to tuberculosis (TB), and the clinical conditions of diabetic TB patients deteriorate faster than non-diabetic TB patients, but the immunological basis for this phenomenon is not understood clearly. Given the role of cell-mediated immunity (CMI) in providing protection against TB, we investigated whether CMI responses in diabetic TB patients are compromised. Peripheral blood mononuclear cells (PBMC) obtained from diabetic TB patients, non-diabetic TB patients and Mycobacterium bovis bacilli Calmette-Guérin (BCG)-vaccinated healthy subjects were cultured in the presence of complex mycobacterial antigens and pools of M. tuberculosis regions of difference (RD)1, RD4, RD6 and RD10 peptides. The PBMC were assessed for antigen-induced cell proliferation and secretion of T helper 1 (Th1) [interferon (IFN)-gamma, interleukin (IL)-2, tumour necrosis factor (TNF)-beta], and Th2 (IL-4, IL-5, IL-10) cytokines as CMI parameters. All the complex mycobacterial antigens and RD1(pool) stimulated strong proliferation of PBMC of all groups, except moderate responses to RD1(pool) in healthy subjects. In response to complex mycobacterial antigens, both IFN-gamma and TNF-beta were secreted by PBMC of all groups whereas diabetic TB patients secreted IL-10 with concentrations higher than the other two groups. Furthermore, in response to RD peptides, IFN-gamma and IL-10 were secreted by PBMC of diabetic TB patients only. The analyses of data in relation to relative cytokine concentrations showed that diabetic TB patients had lower Th1 : Th2 cytokines ratios, and a higher Th2 bias. The results demonstrate a shift towards Th2 bias in diabetic TB patients which may explain, at least in part, a faster deterioration in their clinical conditions.

Short communication: association of HLA-A*1101 with resistance and B*4006 with susceptibility to HIV and HIV-TB: an in silico analysis of promiscuous T cell epitopes

We have shown the association of HLA-A*11 with resistance and HLA-B*40 and -DR2 with susceptibility to HIV and HIV-TB. In the present study, we performed high-resolution subtyping of HLA-A*11 and -B*40 to identify the subtype level association, using the polymerase chain reaction-based sequence-specific oligonucleotide probe method. Underrepresentation of HLA-A*1101 was observed in overall HIV [p(c) = 0.012, OR 0.42 (95% confidence interval (CI) 0.24-0.72)] and HIV(+)TB(+) [p(c) = 0.001, OR 0.18 (95% CI 0.06-0.46)] compared to healthy controls. Significantly higher frequencies of HLA-B*4006 were observed in overall HIV [p = 0.0001, p(c) = 0.004, OR 2.71 (95% CI 1.58-4.75)], HIV(+)TB(-) [p = 0.0003, p(c) = 0.008, OR 2.82 (95% CI 1.56-5.17)], and HIV(+)TB(+) [p = 0.003, p(c) = 0.086, OR 2.56 (95% CI 1.33-4.95)] compared to healthy controls. An in silico analysis of potential T cell epitopes of consensus Gag and Pol sequences of HIV-1 subtype C Indian strains revealed relatively higher number of promiscuous HLA-B40, HLA-DRB1*1501, and -DRB1*1502 (HLA-DR2)-restricted epitopes in contrast to limited numbers of promiscuous binders restricted by HLA-A*1101. The results suggest that HLA-A*1101 may be associated with protection against HIV and the development of TB in HIV patients while HLA-B*4006 may be associated with susceptibility to HIV and TB development in HIV patients. The present study also suggests that the extent of promiscuity of T cell epitopes of HIV-1 subtype C restricted by HLA alleles exerting opposing effects might differ.

Characterization of human cellular immune responses to novel Mycobacterium tuberculosis antigens encoded by genomic regions absent in Mycobacterium bovis BCG

Comparative genomics has identified several regions of differences (RDs) between the infectious Mycobacterium tuberculosis and the vaccine strains of Mycobacterium bovis BCG. We aimed to evaluate the cellular immune responses induced by antigens encoded by genes predicted in 11 RDs. Synthetic peptides covering the sequences of RD1, RD4 to RD7, RD9 to RD13, and RD15 were tested for antigen-induced proliferation and secretion of Th1 cytokine, gamma interferon (IFN-gamma), by peripheral blood mononuclear cells (PBMC) obtained from culture-proven pulmonary tuberculosis (TB) patients and M. bovis BCG-vaccinated healthy subjects. Among the peptide pools, RD1 induced the best responses in both donor groups and in both assays. In addition, testing of TB patients' PBMC for secretion of proinflammatory cytokines (tumor necrosis factor alpha [TNF-alpha], interleukin 6 [IL-6], IL-8, and IL-1beta), Th1 cytokines (IFN-gamma, IL-2, and TNF-beta), and Th2 cytokines (IL-4, IL-5, and IL-10) showed differential effects of RD peptides in the secretion of IFN-gamma and IL-10, with high IFN-gamma/IL-10 ratios (32 to 5.0) in response to RD1, RD5, RD7, RD9, and RD10 and low IFN-gamma/IL-10 ratios (<1.0) in response to RD12, RD13, and RD15. Peptide-mixing experiments with PBMC from healthy subjects showed that secretion of large quantities of IL-10 in response to RD12 and RD13 correlated with inhibition of Th1 responses induced by RD1 peptides. In conclusion, our results suggest that M. tuberculosis RDs can be divided into two major groups--one group that activates PBMC to preferentially secrete IFN-gamma and another group that activates preferential secretion of IL-10--and that these two groups of RDs may have roles in protection against and pathogenesis of TB, respectively.

Overexpression of calcineurin B subunit (CnB) enhances the oncogenic potential of HEK293 cells

Calcineurin (CaN) is a Ca(2+)/calmodulin (CaM)-stimulated protein phosphatase. It is a heterodimeric enzyme consisting of a catalytic A subunit (CnA) and a Ca(2+)-binding regulatory B subunit (CnB). CaN's key role in vivo is well known, while the function of CnB keeps unclear except that it acts as a regulator of CaN. The present paper shows that CnB overexpression promotes proliferation of human embryonic kidney HEK293 cells by comparing with vector control cells in the complete or serum reduced medium. Furthermore, stable CnB transfectants showed dramatically improved growth in soft agar. And the migration ability of CnB overexpressors also was enhanced remarkably. But in the progress of transformation, the activity of CaN remained unchanged between CnB overexpressors and controls. Therefore, CnB, rather than CaN, is involved in the proliferation promotion of HEK293 cells. Subsequently, 11 proteins with different expression levels between CnB transfectants and controls were identified using two-dimensional gel electrophoresis and electrospray ionization time-of-fight mass spectrometry. Therein, the expression of heat shock protein 27 (HSP27) and protein DJ-1 increased along with CnB overexpression. The important role of CnB in cell neoplastic transformation was found and the possible mechanism was analyzed.

Whole genome identification of Mycobacterium tuberculosis vaccine candidates by comprehensive data mining and bioinformatic analyses

Background

Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), infects ~8 million annually culminating in ~2 million deaths. Moreover, about one third of the population is latently infected, 10% of which develop disease during lifetime. Current approved prophylactic TB vaccines (BCG and derivatives thereof) are of variable efficiency in adult protection against pulmonary TB (0%–80%), and directed essentially against early phase infection.

Methods

A genome-scale dataset was constructed by analyzing published data of: (1) global gene expression studies under conditions which simulate intra-macrophage stress, dormancy, persistence and/or reactivation; (2) cellular and humoral immunity, and vaccine potential. This information was compiled along with revised annotation/bioinformatic characterization of selected gene products and in silico mapping of T-cell epitopes. Protocols for scoring, ranking and prioritization of the antigens were developed and applied.

Results

Cross-matching of literature and in silico-derived data, in conjunction with the prioritization scheme and biological rationale, allowed for selection of 189 putative vaccine candidates from the entire genome. Within the 189 set, the relative distribution of antigens in 3 functional categories differs significantly from their distribution in the whole genome, with reduction in the Conserved hypothetical category (due to improved annotation) and enrichment in Lipid and in Virulence categories. Other prominent representatives in the 189 set are the PE/PPE proteins; iron sequestration, nitroreductases and proteases, all within the Intermediary metabolism and respiration category; ESX secretion systems, resuscitation promoting factors and lipoproteins, all within the Cell wall category. Application of a ranking scheme based on qualitative and quantitative scores, resulted in a list of 45 best-scoring antigens, of which: 74% belong to the dormancy/reactivation/resuscitation classes; 30% belong to the Cell wall category; 13% are classical vaccine candidates; 9% are categorized Conserved hypotheticals, all potentially very potent T-cell antigens.

Conclusion

The comprehensive literature and in silico-based analyses allowed for the selection of a repertoire of 189 vaccine candidates, out of the whole-genome 3989 ORF products. This repertoire, which was ranked to generate a list of 45 top-hits antigens, is a platform for selection of genes covering all stages of M. tuberculosis infection, to be incorporated in rBCG or subunit-based vaccines.

Efficient testing of large pools of Mycobacterium tuberculosis RD1 peptides and identification of major antigens and immunodominant peptides recognized by human Th1 cells

Comparative genomics has identified several regions of difference (RDs) of Mycobacterium tuberculosis that are deleted or absent in Mycobacterium bovis BCG vaccines. To determine their relevance for diagnostic and vaccine applications, it is imperative that efficient methods are developed to test the encoded proteins for immunological reactivity. In this study, we have used 220 synthetic peptides covering sequences of 12 open reading frames (ORFs) of RD1 and tested them as a single pool (RD1(pool)) with peripheral blood mononuclear cells obtained from pulmonary tuberculosis (TB) patients and M. bovis BCG-vaccinated healthy subjects in Th1 cell assays that measure antigen-induced proliferation and IFN-gamma secretion. The results showed that RD1(pool) induced strong responses in both TB patients and BCG-vaccinated healthy subjects. The subsequent testing of peptide pools of individual ORFs revealed that all ORFs induced positive responses in a portion of donors, but PPE68, CFP10, and ESAT6 induced strong responses in TB patients and PPE68 induced strong responses in BCG-vaccinated healthy subjects. In addition, HLA-DR and -DQ typing of donors and HLA-DR binding prediction analysis of proteins suggested HLA-promiscuous presentation of PPE68, CFP10, and ESAT6. Further testing of individual peptides showed that a single peptide of PPE68 (121-VLTATNFFGINTIPIALTEMDYFIR-145) was immunodominant. The search for sequence homology revealed that a part of this peptide, 124-ATNFFGINTIPIAL-137, was present in several PPE family proteins of M. tuberculosis and M. bovis BCG vaccines. Further experiments limited the promiscuous and immunodominant epitope region to the 10-amino-acid cross-reactive sequence 127-FFGINTIPIA-136.

Whole blood assays to identify Th1 cell antigens and peptides encoded by Mycobacterium tuberculosis-specific RD1 genes

OBJECTIVE

To identify Th1 cell-stimulating antigens/peptides encoded by the genes predicted in the Mycobacterium tuberculosis-specific genomic region of difference (RD)1, deleted in Mycobacterium bovis Bacille Calmette-Guérin(BCG), by using synthetic peptides and whole blood from tuberculosis (TB) patients.

MATERIALS AND METHODS

Heparinized peripheral blood was obtained from culture-proven pulmonary TB patients (n = 16) attending the Chest Disease Hospital, Kuwait. Whole blood was diluted with tissue culture medium RPMI-1640 and tested for Th1 cell stimulation using antigen-induced proliferation and interferon-gamma (IFN-gamma) secretion assays. The antigens included a peptide pool of 220 peptides covering the sequence of 12 open reading frames (ORFs) of RD1 (RD1(mix)), peptide pools of RD1 ORF5 (ORF5(mix)), ORF6 (ORF6(mix)) and ORF7 (ORF7(mix)), and individual peptides of ORF6 (P6.1-P6.6) and ORF7 (P7.1-P7.6). M. tuberculosis culture filtrate, cell walls and whole-cell M. bovis BCG were used as complex mycobacterial antigens. The results obtained with different antigens and peptides were statistically analyzed for significant differences using Z test.

RESULTS

The complex mycobacterial antigens (culture filtrate, cell walls and M.bovis BCG) and RD1(mix) induced comparable (p > 0.05) positive antigen-induced proliferation and IFN-gamma responses with whole blood from TB patients. However, the positive IFN-gamma responses induced by ORF6(mix) and ORF7(mix) were higher than ORF5(mix). Among the individual peptides, P6.4 and P7.1 of ORF6 and ORF7, respectively, induced the highest IFN-gamma responses, suggesting that these peptides represented the immunodominant Th1 cell epitopes of RD1 ORF6 and ORF7 in the patients tested.

CONCLUSION

The whole blood assays with synthetic peptides are useful to identify Th1 cell antigens/peptides encoded by genes located in M. tuberculosis-specific genomic regions.

Immunoinformatics: an overview of computational tools and techniques for understanding immune function

In recent years, there has been a rapid expansion in the application of information technology to biological data. Although the use of information science techniques is less common for the discipline of immunology, this field has seen great strides in recent years. This review addresses why in silico modeling is needed in immunology research, highlights some of the major areas of research and suggests what may be important for the future of immunoinformatics.

Demonstration of In vivo Expression of a Hypothetical Open Reading Frame (ORF-14) Encoded by the RD1 Region of Mycobacterium tuberculosis

Previously we identified a novel antigenic open reading frame (ORF), designated as ORF-14, on the RD1 region of Mycobacterium tuberculosis that was not originally predicted by Mahairas or by annotation of the M. tuberculosis H37 Rv genome. Here we show that anti-ORF-14 antibodies either from mice immunized with recombinant ORF-14 protein or isolated from serum samples from tuberculosis patients, react with a protein in culture filtrate but not in cytoplasmic or cell wall fractions from M. tuberculosis. Our data indicate that the ORF-14 protein is expressed as a secreted protein, representing one more secreted protein antigen not previously identified by genomics.

Direct ex vivo analyses of HLA-DR1 transgenic mice reveal an exceptionally broad pattern of immunodominance in the primary HLA-DR1-restricted CD4 T-cell response to influenza virus hemagglutinin

The recent threat of an avian influenza pandemic has generated significant interest in enhancing our understanding of the events that dictate protective immunity to influenza and in generating vaccines that can induce heterosubtypic immunity. Although antigen-specific CD4 T cells are known to play a key role in protective immunity to influenza through the provision of help to B cells and CD8 T cells, little is known about the specificity and diversity of CD4 T cells elicited after infection, particularly those elicited in humans. In this study, we used HLA-DR transgenic mice to directly and comprehensively identify the specificities of hemagglutinin (HA)-specific CD4 T cells restricted to a human class II molecule that were elicited following intranasal infection with a strain of influenza virus that has been endemic in U.S. human populations for the last decade. Our results reveal a surprising degree of diversity among influenza virus-specific CD4 T cells. As many as 30 different peptides, spanning the entire HA protein, were recognized by CD4 T cells, including epitopes genetically conserved among H1, H2, and H5 influenza A viruses. We also compared three widely used major histocompatibility class II algorithms to predict HLA-DR binding peptides and found these as yet inadequate for identifying influenza virus-derived epitopes. The results of these studies offer key insights into the spectrum of peptides recognized by HLA-DR-restricted CD4 T cells that may be the focus of immune responses to infection or to experimental or clinical vaccines in humans.

Epitope mapping of HSV-1 glycoprotein K (gK) reveals a T cell epitope located within the signal domain of gK

Glycoprotein K (gK) is a virion envelope component of herpes simplex virus types 1 (HSV-1) and 2 (HSV-2), which plays an important role in virion morphogenesis and egress. We previously demonstrated that immunization of mice with gK, but not with any of the 10 other HSV-1 glycoproteins, resulted in exacerbation of corneal scarring and herpetic dermatitis following ocular HSV-1 infection. However, little is known about the gK epitope(s) that is (are) involved in T cell activities in vitro or in vivo. Thus, epitope mapping of gK was performed using a panel of 15-mer peptides with five-amino acid overlaps spanning the full-length gK, and four expressed gK recombinant proteins representing different regions of gK. Epitope mapping within the gK polypeptide defined the amino acid sequence STVVLITAYGLVLVW as the predominant CD4(+) and CD8(+) T cell stimulatory region both in vitro and in vivo. IFN-gamma expression by CD4(+) T cells was CD8(+) T cells-dependent. This immunodominant epitope is located within the signal sequence of the gK polypeptide and is highly conserved in HSV-1 and HSV-2 strains. Using prediction algorithms, the peptide is predicted to bind to numerous MHC class I and class II molecules.

Field evaluation of a novel differential diagnostic reagent for detection of Mycobacterium bovis in cattle

In the search for improved tools with which to control bovine tuberculosis, the development of enhanced immunodiagnostic reagents is a high priority. Such reagents are required to improve the performance of tuberculin-based reagents and allow the discrimination of vaccinated cattle from those infected with Mycobacterium bovis. In this study, we identified the immunodominant, frequently recognized peptides from Rv3873, Rv3879c, Rv0288, and Rv3019c, which, together with peptides comprising the current lead diagnostic antigens, ESAT-6 and CFP-10, were formulated into a peptide cocktail. In a test of naturally infected cattle, this cocktail was significantly better than tuberculin was for identifying skin test-negative animals with confirmed bovine tuberculosis. In addition, the specificity of this cocktail was not compromised by Mycobacterium bovis BCG vaccination. In summary, our results prioritize this peptide-based, fully synthetic reagent for assessment in larger trials.