In silico identification of common epitopes from pathogenic mycobacteria

  • Bárbara de la Caridad Addine Ramírez1,

    Affiliated with

    • Reynel Marrón1,

      Affiliated with

      • Rommel Calero1,

        Affiliated with

        • Mayelin Mirabal1,

          Affiliated with

          • Juan Carlos Ramírez1,

            Affiliated with

            • María E Sarmiento1,

              Affiliated with

              • Mohd Nor Norazmi2, 3 and

                Affiliated with

                • Armando Acosta1Email author

                  Affiliated with

                  BMC Immunology201314(Suppl 1):S6

                  DOI: 10.1186/1471-2172-14-S1-S6

                  Published: 25 February 2013

                  Abstract

                  An in silico study was carried out to identify antigens for their possible collective use as vaccine candidates against diseases caused by different classes of pathogenic mycobacteria with significant clinical relevance. The genome sequences of the relevant causative agents were used in order to search for orthologous genes among them. Bioinformatics tools permitted us to identify several conserved sequences with 100% identity with no possibility of cross-reactivity to the normal flora and human proteins. Nine different proteins were characterized using the strain H37Rv as reference and taking into account their functional category, their in vivo expression and subcellular location. T and B cell epitopes were identified in the selected sequences. Theoretical prediction of population coverage was calculated for individual epitopes as well as their combinations. Several identical sequences, belonging to six proteins containing T and B cell epitopes which are not present in selected microorganisms of the normal microbial flora or in human proteins were obtained.

                  Introduction

                  There are different species of mycobacteria that cause many human infections with high mortality and morbility [1, 2]. Three of the most common diseases caused by them are tuberculosis (TB), leprosy and Buruli ulcer.

                  Atypical mycobacteria are the causative agents of diseases caused by mycobacteria other than M. tuberculosis (MTB) and M. leprae. Infections with these species have become more frequent nowadays. HIV infection, socioeconomic conditions and antibiotic resistance are some of the reasons for the increased prevalence of these diseases [3]. This fact may also be associated with the waning protective effect of BCG [4].

                  The absence of an effective vaccine against mycobacterial diseases [57] and its high prevalence in poor countries with low accessibility to health services justified the search for vaccine candidates that can be used to simultaneously protect humans against diseases caused by the different pathogenic mycobacteria [8].

                  Vaccine design can be expedited via the application of in silico techniques combined with immunological methods. Bioinformatics tools enable researchers to move rapidly from genomic sequence to vaccine design. These new tools allow the selection of regions of microbial genomes that are predicted to trigger protective immune responses to be used as components for vaccine constucts [9, 10].

                  The present research aims to identify antigens using in silico studies for the development of vaccine candidates that can simultaneously protect humans against diseases caused by pathogenic mycobacteria.

                  Materials and methods

                  The mycobacteria genomes were retrieved from Gene Bank database (http://​www.​ncbi.​nlm.​nih.​gov/​).

                  A whole–genome alignment was performed using Mauve tools in order to search for orthologous genes with identities ranging from 98 to 100% homology.

                  A local alignment among all protein sequences for each orthologous genes was perfomed with Clustalw program, and all sequences having identity of more than 20 continuous amino acids were chosen.

                  Selection of similar genetic regions among mycobacteria and genetic regions from selected bacterial flora of the normal microbiome was also performed.

                  The amino acid sequences obtained were compared with the available sequences of selected microbiome using the local alignment tool FASTA 36.3.4. The parameters for the program were ktup=2, DNA STRAND=N/A and SEQUENCE= PROTEIN.

                  Nine agents were selected for the comparison:

                  Bacteroides fragilis 638R

                  Bifidobacterium bifidum PRL2010

                  Clostridium dificile

                  Escherichia coli 'BL21-Gold (DE3)pLysS AG'

                  Lactobacillus acidophilus NCFM

                  Mycobacterium smegmatis str. MC2 155

                  Staphylococcus epidermidis ATCC 12228,

                  Streptococcus agalactiae 2603V/R

                  Streptococcus sanguinis SK36.

                  Finally, sequences below 70% homology and those without theoretical possibilities of forming linear B cell epitopes were selected.

                  Prediction of subcellular localization of proteins

                  The subcelular localization of the selected proteins was defined using the report of the identification and localization of 1044 MTB proteins [11].

                  Subcelular localization of proteins that did not appear in the report was predicted using three servers: PSORTb [12] , TBpred, [13] and SignalP [14].

                  Identification and selection of in-vivo expressed genes

                  A bibliographic search was performed to look for reports of in-vivo studies of MTB gene expression in humans and animals using google scholar, PubMed and PubMed Central.

                  HLAPred server was used for prediction of T-cell epitopes in this study

                  Thirty six HLA class I and 51 HLA class II were selected for the prediction. The results were displayed in HTML Mapping form and the threshold of 3% were use as default for the prediction parameters [11]. T cell epitopes similar to human epitopes were eliminated.

                  B-cell epitope prediction

                  Two servers, Bcepred and ABCpred, were combined for the prediction of linear B cell epitomes. For Bcepred server: Seven physico-chemical properties of amino acids (hydrophilicity, flexibility, accessibility, polarity, exposed surface and turns and antigenic propensity) were combined with a threshold at 2.38. For ABCpred server: A threshold of 0.5 and the predicted B cell epitope length of 16 amino acids were used [12]. The regions containing both T and B cell epitopes were selected by BioHelper tools [15].

                  The epitopes analyzed were compared with the report of Iñakis et al. [16] for eliminating hyperconserved epitopes.

                  Population coverage calculation of individual epitopes and their combinations

                  The theoretical prediction of the presentation of the individual epitopes and their combinations for MHC alleles of Cuban, Malaysian, Brazilian, Mexican, Australian, African, and North American populations were calculated using the Population Coverage Calculation program [18].

                  The percentage of coverage greater than 70% was accepted as good consistent with reports in the literature [18, 10, 19, 11, 14] (figure 1).
                  http://static-content.springer.com/image/art%3A10.1186%2F1471-2172-14-S1-S6/MediaObjects/12865_2013_Article_503_Fig1_HTML.jpg
                  Figure 1

                  Flowchart of the procedure –bioinformatics analysis

                  Results

                  Thirteen-whole mycobacteria genome alignments were performed and 26 orthologous genes were determined. A local alignment among each protein sequence for each orthologous gene was performed and 71 sequences were obtained with more than 20 identical amino acids between them, corresponding to 14 proteins.

                  In order to avoid the selection of epitopes share with the host microbiome only, only 20 epitopes were retained in this study. It was found that 13 of 14 proteins were expressed in the membrane; with six only expressed in the membrane, two in the membrane and cell wall; three in the membrane and cytosol and two in all of these compartments. One protein was only expressed in the cytosol [20, 21].

                  Table 1 shows a characterization of predicted proteins containing sequences showing 100% identity among pathogenic mycobacteria, and not shared with human proteins or the normal human microbiome. These proteins have also been predicted to contain T and/ or B cell epitopes, and are expressed during infection.
                  Table 1

                  Characterization of predicted proteins

                  Identification (ORF)

                  T epitopes

                  B epitopes

                  Up-regulation during infection

                  Rv0667

                  YES

                  YES

                  YES

                  Rv1299

                  YES

                  NO

                  YES

                  Rv1315

                  YES

                  YES

                  NO

                  Rv1420

                  YES

                  NO

                  NO

                  Rv1421

                  YES

                  YES

                  NO

                  Rv1547

                  YES

                  YES

                  YES

                  Rv0701

                  NO

                  YES

                  NO

                  Rv1308

                  NO

                  YES

                  NO

                  Rv1384

                  NO

                  YES

                  YES

                  Four of the identified genes (Rv1269, Rv0667, Rv1547 and Rv1384) have been reported to be up-regulated during infection. Rv1299 is over-expressed during infection in human macrophages, mice and guinea pig lungs [22]. Rv0667, Rv1547and Rv1384 are over-expressed only during infection in guinea pig lungs [23].

                  The selected sequences from Rv0701, Rv1384 and Rv1308 are not predicted to contain T cell epitopes but were chosen based on reports of their potential role in the humoral immune response against MTB [2325].

                  Hyperconserved epitopes are not present in the selected sequences. These sequences were excluded to avoid regions, which may be selected during the co-evolution of MTB and humans and hence may provide an advantage to the pathogen [17].

                  Five sequences that contain epitopes similar to human epitopes were eliminated [13] to avoid potential autoimmune response.

                  In this research, the identified sequences belonging to 15 MTB proteins and some of its combinations were analyzed. The combination of all sequences were predicted to give more than 70% population coverage when presented via the MHC class I plus MHC class II molecules in the selected geographical regions. In contrast, the predicted population coverage of individual sequences was highly variable in the different geographical regions and was always bellow 70 % (data not shown). However, these sequences were not excluded because of their potential use in the specific geographical regions.

                  The combination of all sequences shows population coverage over 70 % in all geographical regions studied (Figure 2).
                  http://static-content.springer.com/image/art%3A10.1186%2F1471-2172-14-S1-S6/MediaObjects/12865_2013_Article_503_Fig2_HTML.jpg
                  Figure 2

                  Population coverage for combination of all sequences

                  Evaluation of the identified sequences in in vitro and in animal models should be carried out to verify their utility to be included in vaccine formulations.

                  Conclusions

                  We used bioinformatics tools to predict sequences, which are common among different mycobacteria belonging to the MTB complex and are absent from selected microorganisms of the normal flora. The predicted sequences comprise those that are expressed primarily in the membrane and contain B and T cell epitopes that do not match with the human epitopic database. The combination of the identified sequences gives appropriated theoretical population coverage in the different geographic regions studied. These predicted sequences are potential vaccine candidates but functional/biological assays should be performed to verify whether they are indeed appropriate to be included in a vaccine formulation.

                  Declarations

                  Acknowledgements

                  This work was supported by the Ministry of Science, Technology and Innovation, Malaysia (10-01-05-MEB002) and the Ministry of Higher Education, Malaysia LRGS Grant (203.PSK.6722001) and Ministry of Science and Technology, Cuba

                  Declarations

                  This article has been published as part of BMC Immunology Volume 14 Supplement 1, 2013: Proceedings from Delivery Systems and Current strategies to drug design. The full contents of the supplement are available online at http://​www.​biomedcentral.​com/​bmcimmunol/​supplements/​14/​S1

                  Authors’ Affiliations

                  (1)
                  Finlay Institute
                  (2)
                  School of Health Sciences Universiti Sains Malaysia
                  (3)
                  Institute for Research in Molecular Medicine, Universiti Sains Malaysia

                  References

                  1. Smith I: Mycobacterium tuberculosis Pathogenesis and Molecular Determinants of Virulence. Clin. Microbiol. Rev. 2003, 16 (3): 463-496. 10.1128/CMR.16.3.463-496.2003.PubMed CentralView ArticlePubMed
                  2. Martín C, Bigi F, Gicquel B: New Vaccines against Tuberculosis. Tuberculosis 2007: From basic science to patient care. 2007, 341-360. First
                  3. Vargas García R: Epidemiología de las micobacteriosis atípicas en la salud humana. BNF Emerge. 2008, 2 (4): 220-230.
                  4. Marcel AB: Mycobacterium du jour: what's on tomorrow's menu?. Microbes and Infection. 2008, 10 (9): 968-972. 10.1016/j.micinf.2008.07.001.View Article
                  5. Young DB, Perkins MD, Duncan K, Barry CE: Confronting the scientific obstacles to global control of tuberculosis. J Clin Invest. 2008, 118 (4): 1255-65. 10.1172/JCI34614.PubMed CentralView ArticlePubMed
                  6. Dye C: Tuberculosis 2000-2010: control, but not elimination. Int J Tuberc Lung Dis. 2000, 4 (12 Suppl 2): S146-52.PubMed
                  7. Norazmi MN, Sarmiento ME, Acosta A: Recent advances in tuberculosis vaccine development. Current Respiratory Medicine Reviews. 2005, 1 (2): 109-116. 10.2174/1573398054023000.View Article
                  8. Ochoa Azze RF: Inmunoepidemiología y estrategia de vacunación. Ediciones Finlay. 2008, Ciudad de la Habana
                  9. Pearson WR, Lipman DJ: Improved tools for biological sequence comparison. Proceedings of the National Academy of Sciences of the United States of America. 1988, 85 (8): 2444-8. 10.1073/pnas.85.8.2444.PubMed CentralView ArticlePubMed
                  10. Li Pira G, Ivaldi F, Moretti P, Manca F: High throughput T epitope mapping and vaccine development. J Biomed Biotechnol. 2010, Epub 2010: 325720
                  11. Saha S, Raghava GPS: Prediction of Continuous B-cell Epitopes in an Antigen Using Recurrent Neural Network. Proteins. 2006, 65 (1): 40-48. 10.1002/prot.21078.View ArticlePubMed
                  12. Mawuenyega KG, Forst CV, Dobos KM, Belisle JT, Chen J, Bradbury EM: Mycobacterium tuberculosis functional network analysis by global subcellular protein profiling. Molecular Biology of the Cell. 2004, 16: 396-404. 10.1091/mbc.E04-04-0329.View ArticlePubMed
                  13. Yu NY, Wagner JR, Laird MR, Melli G, Rey S, Lo R: PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics. 2010, 26 (13): 1608-15. 10.1093/bioinformatics/btq249.PubMed CentralView ArticlePubMed
                  14. Rashid M, Saha S, Raghava GPS: Support Vector Machine-based method for predicting subcellular localization of mycobacterial proteins using evolutionary information and motifs. BMC Bioinformantics. 2007, 8: 337-10.1186/1471-2105-8-337.View Article
                  15. Romel Calero Ramos: Análisis in silico de proteínas como candidatos para el diagnóstico de la infección humana por MTB. 2011, [Tesis].Universidad de La Habana
                  16. McNamara LA, Yongqun H, Zhenhua Y: Use searching a T epitope predictions to evaluate efficacy and population coverage of the MTB72f vaccine for tuberculosis. BMC Immunology. 2010, 11: 18-10.1186/1471-2172-11-18.PubMed CentralView ArticlePubMed
                  17. Comas Iñaki, Chakravartti Jaidip, Small Peter M: Human T cell epitopes of Mycobacterium tuberculosis are evolutionarily hyperconserved. Nat Genet. 2010, 42 (6): 498-503. 10.1038/ng.590.PubMed CentralView ArticlePubMed
                  18. Longmate J, York J, La Rosa C: Population coverage by HLA class-I restricted cytotoxic T-lymphocyte epitopes. Immunogenetics. 2008, 52 (3-4): 165-73.View Article
                  19. Bhasin M, Raghava GPS: A hybrid approach for predicting promiscuous MHC class I restricted T cell epitopes. J. Biosci. 2006, 32: 31-42.View Article
                  20. Kumar M, Khan FG, Sharma S, Kumar R, Faujdar J, Sharma R: Identification of Mycobacterium tuberculosis genes preferentially expressed during human infection. Microb Pathog. 2011, 50 (1): 31-8. 10.1016/j.micpath.2010.10.003.View ArticlePubMed
                  21. Mark T, John S, Derek P, Srinivas V, Ramaswamy , Graviss EA: Detection of rpoB Mutations Associated with Rifampin Resistance in Mycobacterium tuberculosis Using Denaturing Gradient Gel Electrophoresis. Antimicrob Agents Chemother. 2005, 49 (6): 2200-09. 10.1128/AAC.49.6.2200-2209.2005.View Article
                  22. Banerjee R, Vats P, Sonal D, Sunitha MK, Rajendra J: Comparative Genomics of Cell Envelope Components in Mycobacteria. BMC Microbiol. 2009, 14: 47-
                  23. Talaat AM, Lyons R, Howard ST, Johnston SA: The temporal expression profile of Mycobacterium tuberculosis infection in mice. Proc Natl Acad Sci USA. 2004, 101 (13): 4602-7. 10.1073/pnas.0306023101.PubMed CentralView ArticlePubMed
                  24. Abbas AK, Lichtman AH, Pober JS: Inmunología celular y molecular. Barcelona. Elsevier. 2008, 6
                  25. Karakousis PC, Yoshimatsu T, Lamichhane G, Woolwine SC, Nuermberger EL, Jacques Grosset: Dormancy Phenotype Displayed by Extracellular Mycobacterium tuberculosis within Artificial Granulomas in Mice. J Exp Med. 2004, 200 (5): 647-657. 10.1084/jem.20040646.PubMed CentralView ArticlePubMed

                  Copyright

                  © de la Caridad Addine Ramírez et al; licensee BioMed Central Ltd. 2013

                  This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

                  Advertisement