- Research note
- Open Access
- Published:
Bioinformatics analysis of calcium-dependent protein kinase 4 (CDPK4) as Toxoplasma gondii vaccine target
BMC Research Notes volume 14, Article number: 50 (2021)
Abstract
Objectives
Toxoplasma gondii (T. gondii), an obligate intracellular apicomplexan parasite, could affect numerous warm-blooded animals, such as humans. Calcium-dependent protein kinases (CDPKs) are essential Ca2+ signaling mediators and participate in parasite host cell egress, outer membrane motility, invasion, and cell division.
Results
Several bioinformatics online servers were employed to analyze and predict the important properties of CDPK4 protein. The findings revealed that CDPK4 peptide has 1158 amino acid residues with average molecular weight (MW) of 126.331 KDa. The aliphatic index and GRAVY for this protein were estimated at 66.82 and – 0.650, respectively. The findings revealed that the CDPK4 protein comprised 30.14% and 34.97% alpha-helix, 59.84% and 53.54% random coils, and 10.02% and 11.49% extended strand with SOPMA and GOR4 tools, respectively. Ramachandran plot output showed 87.87%, 8.40%, and 3.73% of amino acid residues in the favored, allowed, and outlier regions, respectively. Also, several potential B and T-cell epitopes were predicted for CDPK4 protein through different bioinformatics tools. Also, antigenicity and allergenicity evaluation demonstrated that this protein has immunogenic and non-allergenic nature. This paper presents a basis for further studies, thereby provides a fundamental basis for the development of an effective vaccine against T. gondii infection.
Introduction
Toxoplasma gondii (T. gondii), a compulsory intracellular parasite, could affect nearly all warm-blooded animals, such as humans [1, 2]. The current drugs against toxoplasmosis are not effective enough and are associated with serious side effects [3]. Therefore, the finding and development of an effective vaccine have high priority and is critically required to limit T. gondii infection [4,5,6]. Calcium-dependent protein kinases (CDPKs) are a category of serine/threonine kinases found only in plants and protists like ciliates and apicomplexan parasites [7, 8]. Multiple CDPKs (the most important household proteins) have been identified in apicomplexan protists, especially in T. gondii. The CDPKs are essential Ca2+ signaling mediators and participate in parasite host cell egress, outer membrane motility, invasion, and cell division [9,10,11,12]. The CDPK family was regarded as a good choice for anti-Toxoplasma medications and an appropriate option for the design of vaccines [13,14,15,16,17]. There is no report regarding immunization with CDPK4 in experimental animals yet. However, in several studies, vaccination with CDPK1 [13, 18, 19], CDPK2 [20], CDPK3 [15, 21, 22], CDPK5 [14], and CDPK6 [16] induced strong humoral and cellular responses and prolonged the survival time in mouse models.
Predicting epitopes is highly important for the determination of an antigen’s immunogenicity in vaccine development. So, bioinformatics tools and online resources can enable researchers to predict and recognize the potential epitopes of B- and T-cells [23,24,25,26]. Bioinformatics has recently become the preferred interdisciplinary new science for analyzing biological data using defined computer science, mathematics, statistics, physics, biology, medicine technologies and algorithms [23, 26]. Thus, the current study was performed to analyze the several important features of CDPK4 protein using different bioinformatics servers.
Methodology
Retrieval of CDPK4 protein sequence of T. gondii
First, the complete amino acid sequence of T. gondii CDPK4 protein was attained from ToxoDB server (https://toxodb.org/toxo/).
Physicochemical evaluations
The physicochemical characteristics of the CDPK4 protein are of crucial significance in the evaluation of its aliphatic index, half-life, theoretical isoelectric point (pI), a grand average of hydropathicity (GRAVY), and electric charge distribution. The mentioned features were explored by ProtParam server [27, 28].
Prediction of acylation and phosphorylation sites of CDPK4
In order to predict acylation and phosphorylation sites of CDPK4 protein, CSS-Palm and NetPhos online tools were employed, respectively [29, 30].
Transmembrane domains and subcellular location prediction
The transmembrane domains of CDPK4 were examined by the TMHMM server v.2.0. Moreover, the PSORT II server was applied to predict the subcellular position of the CDPK4 protein.
Prediction of secondary structure
Garnier–Osguthorpe–Robson 4 (GOR4) [31], SOPMA [32], and PSIPRED [33] online servers were utilized for prediction of the secondary structure of the CDPK4 protein.
Construction of the 3D model
3D models play a decisive role in the development of vaccines. In this case, the SWISS-MODEL webserver was applied to build the three-dimensional models of the CDPK4 protein through a homology modeling approach [34].
Refining and validating the 3D modeled structure
The proper SWISS-MODEL-constructed model was selected and modified by GalaxyRefine to attain high-quality template-based protein predictions [35]. The Ramachandran plot was utilized to validate the three-dimensional construct of the CDPK4 protein through the use of the SWISS-MODEL software [36]. The quality of the model was checked out using ProSAweb [37, 38].
Prediction of B-cell epitopes
Various databases were employed for the analysis of the linear epitopes of the B-cells, including, Bcepred [39], ABCpred [40], SVMTriP [41], and the immune epitope database (IEDB). Furthermore, discontinuous B-cell epitopes were predicted by ElliPro [42].
Prediction of major histocompatibility complex-I (MHC-I) and MHC-II epitopes
In this study, the IEDB [43] and NetMHCcons 1.1 [44] online services were used to predict the binding affinity of CDPK4 peptides toward the MHC class I. Furthermore, IEDB [45] and NetMHCIIpan 3.2 [46] servers were exploited to examine the 15-mer T- cell epitopes of H-2-IEd, H2IAd, and H2IAb mouse alleles.
Detection of the CTL epitopes
To activate the immune system, an antigen should be first presented on the MHC-I surface. So, the choice of cytotoxic T lymphocyte (CTL) epitopes plays a decisive role in designing a vaccine. To this end, a free web server CTLpred [47] was utilized.
Allergenicity, immunogenicity, and solubility evaluation
The possible allergenicity of the CDPK4 protein was explored utilizing AllerTOP [48] and AlgPred [28] servers. Protein antigenicity of CDPK4 protein was estimated using ANTIGENpro [49] and VaxiJen v. 2.0 [50] servers. Furthermore, we used the SOLpro online site to determine the protein solubility [38].
Results and discussion
Initial overview of the protein CDPK4
The amino acid sequence of CDPK4 protein was determined by the ToxoDB server under the accession ID: TGME49_237890. The CDPK4 protein includes 1158 amino acid residues with an estimated molecular weight of 126.331 KDa (antigens which have MW of < 5–10 KDa are considered as weak immunogens) [51], whereas its theoretical pI was 9.15. The total number of residues with the negative and positive charge was 145 and 167, respectively. Its half-life was predicted as 30 h, > 20 h, and > 10 h in mammalian reticulocytes cells (in vitro), yeast cells, and E. coli, respectively. Based on the instability index results (58.84), the CDPK4 protein is unstable. The relatively good estimated aliphatic index value of 66.82 indicates the thermostability of the protein. According to GRAVY data (− 0.650), CDPK4 protein exhibits hydrophilic features.
Prediction of PTM sites of CDPK4
As it is evident, PTMs have important roles in cellular control processes [52]. Based on the findings, 137 phosphorylation and 21 acylation sites were predicted on CDPK4 sequence, suggesting that these sites may control several functions of the protein and affect protein activity (Fig. 1; Additional file 1: Table S1).
Identification of transmembrane domains and subcellular location
The data obtained from TMHMM server v. 2.0 indicated no transmembrane domain in the CDPK4 sequence (Additional file 2: Figure S1). Furthermore, based on the PSORT II prediction, the subcellular location of CDPK4 is as follows: 82.6% nuclear, 8.7% plasma membrane, 4.3% mitochondrial, and 4.3% cytoskeletal.
Secondary structure assessment
It should be noted that determining the protein secondary structure by introducing special constraints, such as beta-turn or alpha helix, is a key step in the assessment of the tertiary structure. The findings showed the CDPK4 protein comprised 30.14% (349/1158) and 34.97% (405/1158) alpha-helix, 59.84% (693/1158) and 53.54% (620/1158) random coils, and 10.02% (116/1158) and 11.49% (133/1158) extended strand by SOPMA and GOR4 servers, respectively (Additional file 2: Figures S2 and S3). The findings from the PSIPRED server are depicted in Additional file 2: Figure S4. It is apparent that alpha-helix and beta-turn placed in the protein’s internal portions, with high hydrogen bond-energy, will maintain the protein's structure resulting in a better interaction with antibodies [38, 53]. The principal biological behavior of the proteins is focused on their spatial structure. Knowledge of protein structures and awareness the relationships between structures and functions are important [38].
3D model analysis
Following the analysis, five 3D models were established for the CDPK4 sequence among which, the one with the highest identity was chosen. The chosen template exhibited a 34.99% sequence identity. The SWISS-MODEL results are presented in Additional file 2: Figure S5.
Refinement and validation of the tertiary structure
The GalaxyRefine software was employed to refine the tertiary structures. According to the results of the Ramachandran plot and ProSAweb servers, an enhancement was observed in the quality of the three-dimensional structure after the refinement. Prior to the refinement process, validation of the protein using the SWISS-MODEL tool showed that 87.87% of residues were situated in favored regions, while 8.40% and 3.73% of them lied in allowed and outlier regions, respectively, verifying its immunogenic efficiency (Fig. 2c). The post-refinement exploration of Ramachandran plots demonstrated that 95.34% of the residues lied within the favored regions, whereas 3.54% of them were in the allowed regions and only 1.12% of the residues were placed the outlier regions (Fig. 2d). The Z-score is indicative of the model quality; this parameter was − 8.09 in the initial model (based on the ProSA-web), and the majority of residues lied in the favored regions (Fig. 2a). Further improvement in the quality of the 3D structure after refinement can be also inferred from the Z-score value (− 8.15) (Fig. 2b).
Validation of CDPK4 protein 3D structure using Ramachandran plot. a The Z-score plot for 3D structure of predicted protein before refinement with ProSA-web server. b The Z-score plot for 3D structure of predicted protein after refinement with ProSA-web server. c The analysis of Ramachandran plot using SWISS-MODEL server in initial model showed 87.87%, 8.40% and 3.73% of residues were located in favored, allowed and outlier regions, respectively. d The results after refinement were changed as follow: 95.34%, 3.54% and 1.12% of residues were located in favored, allowed and outlier regions, respectively
B-cell epitopes prediction
The epitope prediction could offer invaluable data that can be used to identify the immunogenic peptides. The Bcepred-determined linear B-cell epitopes are presented in Additional file 1: Table S2; while the results obtained from the ABCpred server are listed in Additional file 1: Table S3. The greater peptide score indicates the higher possibility of being an epitope. According to the IEDB online tool, the mean scores of hydrophilicity, antigenicity, beta-turn, bepipred linear epitope prediction, flexibility, and surface accessibility of the CDPK4 protein are 2.381, 1.016, 1.012, 0.350, 1.013, and 1.000, respectively (Additional file 2: Figure S6). The SVMTriP-derived results are also tabulated in Additional file 1: Table S4. The analysis of linear B-cell epitopes demonstrated that the CDPK4 protein contains favorable epitopes and appropriate indices. The Bcepred’s estimation accuracy of models based on various properties differs from 52.92% to 57.53%. This server also assists to forecast epitopes of B-cells using physicochemical features [38, 39]. Another valuable step for the in silico analysis is the identification of the conformational epitopes needed for antibody-antigen interaction [30]. In this case, the application of ElliPro tool resulted in five discontinuous B-cell epitopes (Table 1).
Analysis of MHC-I and MHC-II molecules
The connection of peptides to MHC molecules is an important step in the presentation of antigens to T-cells and also a key step in the choice of potential epitopes. Additional file 1: Tables S5 and S6, respectively indicate the lowest percentile ranks for CDPK4 MHC-I and MHC-II as obtained from the IEDB site. The results derived from NetMHCcons and NetMHCIIpan are also listed in Additional file 1: Tables S7 and S8, respectively. In general, the lower percentile ranks (or IC50 values) indicate the higher level affinity, which represents a better T-cell epitopes and vice versa [28]. Based on the bioinformatics analyses, CDPK4 T-cell epitopes could strongly bind to MHC-I and MHC-II molecules. Since T. gondii is considered as an intracellular protozoa, the cellular immunity mediated by the T-cells have a pivotal role against this microorganism [54]. It is therefore extremely essential for the development of an effective vaccine against T. gondii to explain which type of T-cell-mediated immune response is participated [38, 54].
Prediction of the CTL epitopes
The CTLpred server was utilized to select 10 high-rank and suitable epitopes to analyze the CTL epitope. The details are mentioned in Additional file 1: Table S9.
Allergenicity, immunogenicity, and solubility analysis
CDPK4 protein could exhibit high immunogenicity as its antigenicity score was determined 0.780 and 0.622 (through the use of ANTIGENpro and VaxiJen servers), respectively. AllerTOP and AlgPred servers suggest the non-allergic features of this protein. The ability to determine allergenicity is important to make sure that candidates for vaccines are low in allergenicity [38]. Based on the SOLpro output, the solubility of the CDPK4 protein was determined as 0.7087.
Conclusion
This paper provided a detailed explanation of the fundamental aspects of CDPK4 protein, such as physicochemical characteristics, transmembrane domains, secondary and tertiary structures, B- and T-epitopes, and other features of CDPK4, using bioinformatics servers. Based on the findings, CDPK4 protein revealed an acceptable antigenicity score. Also, this protein contains many good epitopes of B- and T-cells, suggesting that CDPK4 can considered as an appropriate vaccine candidate against T. gondii. This research presented important fundamental and theoretical evidence for further in vivo investigations on the CDPK4 protein to establish an effective vaccine against acute and chronic T. gondii infection.
Limitations
In this paper, only in silico analysis was performed. More studies are recommended for the development of an effective vaccine in vivo using the CDPK4 alone or combined with other antigens in the future. Also, a combination of immunodominant CDPK4 epitopes with various adjuvants and formulations will be useful.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the ToxoDB repository (https://toxodb.org/toxo/app/record/gene/TGME49_237890).
Abbreviations
- 3D:
-
Three-dimensional
- CDPK:
-
Calcium-dependent protein kinase
- CTL:
-
Cytotoxic T lymphocyte
- GOR:
-
Garnier–Osguthorpe–Robson
- GRAVY:
-
Grand average of hydropathicity
- IC50 :
-
Half maximal inhibitory concentration
- IEDB:
-
Immune epitope database
- MHC:
-
Major histocompatibility complex
- MW:
-
Molecular weight
- PDB:
-
Protein data bank
- pI:
-
Isoelectric point
- PTM:
-
Post-translational modification
- T. gondii :
-
Toxoplasma gondii
References
Cenci-Goga BT, Rossitto PV, Sechi P, McCrindle CM, Cullor JS. Toxoplasma in animals, food, and humans: an old parasite of new concern. Foodborne Pathog Dis. 2011;8(7):751–62.
Foroutan M, Fakhri Y, Riahi SM, Ebrahimpour S, Namroodi S, Taghipour A, Spotin A, Gamble HR, Rostami A. The global seroprevalence of Toxoplasma gondii in pigs: a systematic review and meta-analysis. Vet Parasitol. 2019;269:42–52.
Antczak M, Dzitko K, Długońska H. Human toxoplasmosis—searching for novel chemotherapeutics. Biomed Pharmacother. 2016;82:677–84.
Zhang NZ, Wang M, Xu Y, Petersen E, Zhu XQ. Recent advances in developing vaccines against Toxoplasma gondii: an update. Expert Rev Vaccines. 2015;14(12):1609–21.
Foroutan M, Ghaffarifar F, Sharifi Z, Dalimi A, Jorjani O. Rhoptry antigens as Toxoplasma gondii vaccine target. Clin Exp Vaccine Res. 2019;8(1):4–26.
Foroutan M, Zaki L, Ghaffarifar F. Recent progress in microneme-based vaccines development against Toxoplasma gondii. Clin Exp Vaccine Res. 2018;7(2):93–103.
Tzen M, Benarous R, Dupouy-Camet J, Roisin M. A novel Toxoplasma gondii calcium-dependent protein kinase. Parasite. 2007;14(2):141–7.
Billker O, Lourido S, Sibley LD. Calcium-dependent signaling and kinases in apicomplexan parasites. Cell Host Microbe. 2009;5(6):612–22.
Lourido S, Shuman J, Zhang C, Shokat KM, Hui R, Sibley LD. Calcium-dependent protein kinase 1 is an essential regulator of exocytosis in Toxoplasma. Nature. 2010;465(7296):359–62.
Lourido S, Tang K, Sibley LD. Distinct signalling pathways control Toxoplasma egress and host-cell invasion. EMBO J. 2012;31(24):4524–34.
Morlon-Guyot J, Berry L, Chen CT, Gubbels MJ, Lebrun M, Daher W. The T oxoplasma gondii calcium-dependent protein kinase 7 is involved in early steps of parasite division and is crucial for parasite survival. Cell Microbiol. 2014;16(1):95–114.
Uboldi AD, McCoy JM, Blume M, Gerlic M, Ferguson DJ, Dagley LF, Beahan CT, Stapleton DI, Gooley PR, Bacic A. Regulation of starch stores by a Ca2+-dependent protein kinase is essential for viable cyst development in Toxoplasma gondii. Cell Host Microbe. 2015;18(6):670–81.
Chen J, Li Z-Y, Huang S-Y, Petersen E, Song H-Q, Zhou D-H, Zhu X-Q. Protective efficacy of Toxoplasma gondiicalcium-dependent protein kinase 1 (TgCDPK1) adjuvated with recombinant IL-15 and IL-21 against experimental toxoplasmosis in mice. BMC Infect Dis. 2014;14(1):487.
Zhang N-Z, Huang S-Y, Xu Y, Chen J, Wang J-L, Tian W-P, Zhu X-Q. Evaluation of immune responses in mice after DNA immunization with putative Toxoplasma gondii calcium-dependent protein kinase 5. Clin Vaccine Immunol. 2014;21(7):924–9.
Zhang N-Z, Huang S-Y, Zhou D-H, Chen J, Xu Y, Tian W-P, Lu J, Zhu X-Q. Protective immunity against Toxoplasma gondii induced by DNA immunization with the gene encoding a novel vaccine candidate: calcium-dependent protein kinase 3. BMC Infect Dis. 2013;13(1):512.
Zhang N-Z, Xu Y, Wang M, Chen J, Huang S-Y, Gao Q, Zhu X-Q. Vaccination with Toxoplasma gondii calcium-dependent protein kinase 6 and rhoptry protein 18 encapsulated in poly (lactide-co-glycolide) microspheres induces long-term protective immunity in mice. BMC Infect Dis. 2016;16(1):168.
Foroutan M, Ghaffarifar F. Calcium-dependent protein kinases are potential targets for Toxoplasma gondii vaccine. Clin Exp Vaccine Res. 2018;7(1):24–36.
Chen J, Li ZY, Petersen E, Liu WG, Zhu XQ. Co-administration of interleukins 7 and 15 with DNA vaccine improves protective immunity against Toxoplasma gondii. Exp Parasitol. 2016;162:18–23.
Huang SY, Chen K, Wang JL, Yang B, Zhu XQ. Evaluation of protective immunity induced by recombinant calcium-dependent protein kinase 1 (TgCDPK1) protein against acute toxoplasmosis in mice. Microb Pathog. 2019;133:103560.
Chen K, Wang JL, Huang SY, Yang WB, Zhu WN, Zhu XQ. Immune responses and protection after DNA vaccination against Toxoplasma gondii calcium-dependent protein kinase 2 (TgCDPK2). Parasite. 2017;24:41.
Zhang NZ, Gao Q, Wang M, Elsheikha HM, Wang B, Wang JL, Zhang FK, Hu LY, Zhu XQ. Immunization with a DNA vaccine cocktail encoding TgPF, TgROP16, TgROP18, TgMIC6, and TgCDPK3 genes protects mice against chronic toxoplasmosis. Front Immunol. 2018;9:1505.
Wu M, An R, Chen Y, Chen T, Wen H, Yan Q, Shen J, Chen L, Du J. Vaccination with recombinant Toxoplasma gondii CDPK3 induces protective immunity against experimental toxoplasmosis. Acta Trop. 2019;199:105148.
Flower DR, Macdonald IK, Ramakrishnan K, Davies MN, Doytchinova IA. Computer aided selection of candidate vaccine antigens. Immun Res. 2010;6(2):S1.
Wang Y, Wang G, Cai J, Yin H. Review on the identification and role of Toxoplasma gondii antigenic epitopes. Parasitol Res. 2016;115(2):459–68.
Keyvani H, Ahmadi NA, Ranjbar MM, Ataei Kachooei S, Ghorban K, Dadmanesh M. Immunoinformatics study of gp120 of human immunodeficiency virus type 1 subtype CRF35_AD isolated from Iranian patients. Arch Clin Infect Dis. 2016;11(4):e36270.
Ranjbar MM, Nayeb Ali A, Ghorban K, Ghalyanchi Langeroudi A, Dadmanesh M, Amini H-R, Sedighi Moghaddam B. Immnoinformatics: novel view in understanding of immune system function, databases and prediction of immunogenic epitopes. Koomesh. 2015;17(1):18–26.
Gasteiger E, Hoogland C, Gattiker A, Wilkins MR, Appel RD, Bairoch A. Protein identification and analysis tools on the ExPASy server. In: The proteomics protocols handbook. Berlin: Springer; 2005. p. 571–607.
Majidiani H, Soltani S, Ghaffari AD, Sabaghan M, Taghipour A, Foroutan M. In-depth computational analysis of calcium-dependent protein kinase 3 of Toxoplasma gondii provides promising targets for vaccination. Clin Exp Vaccine Res. 2020;9(2):146–58.
Zhou J, Wang L, Zhou A, Lu G, Li Q, Wang Z, Zhu M, Zhou H, Cong H, He S. Bioinformatics analysis and expression of a novel protein ROP48 in Toxoplasma gondii. Acta Parasitol. 2016;61(2):319–28.
Foroutan M, Ghaffarifar F, Sharifi Z, Dalimi A, Pirestani M. Bioinformatics analysis of ROP8 protein to improve vaccine design against Toxoplasma gondii. Infect Genet Evol J Mol Epidemiol Evol Genet Infect Dis. 2018;62:193–204.
Garnier J, Gibrat J-F, Robson B. [32] GOR method for predicting protein secondary structure from amino acid sequence. In: Methods in enzymology. The Netherlands: Elsevier; 1996. p. 540–53.
Deléage G. ALIGNSEC: viewing protein secondary structure predictions within large multiple sequence alignments. Bioinformatics. 2017.
McGuffin LJ, Bryson K, Jones DT. The PSIPRED protein structure prediction server. Bioinformatics. 2000;16(4):404–5.
Guex N, Peitsch MC, Schwede T. Automated comparative protein structure modeling with SWISS-MODEL and Swiss-PdbViewer: a historical perspective. Electrophoresis. 2009;30(S1):S162–73.
Park H, Seok C. Refinement of unreliable local regions in template-based protein models. Proteins Struct Funct Bioinform. 2012;80(8):1974–86.
Bertoni M, Kiefer F, Biasini M, Bordoli L, Schwede T. Modeling protein quaternary structure of homo-and hetero-oligomers beyond binary interactions by homology. Sci Rep. 2017;7(1):1–15.
Wiederstein M, Sippl MJ. ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007;35(Suppl2):W407–10.
Ghaffari AD, Dalimi A, Ghaffarifar F, Pirestani M. Structural predication and antigenic analysis of ROP16 protein utilizing immunoinformatics methods in order to identification of a vaccine against Toxoplasma gondii: An in silico approach. Microb Pathog. 2020;142:104079.
Saha S, Raghava GPS. BcePred: prediction of continuous B-cell epitopes in antigenic sequences using physico-chemical properties. In: International Conference on Artificial Immune Systems. Springer, Berlin, pp. 197–204; 2004.
Saha S, Raghava GPS. Prediction of continuous B-cell epitopes in an antigen using recurrent neural network. Proteins Struct Funct Bioinform. 2006;65(1):40–8.
Yao B, Zhang L, Liang S, Zhang C. SVMTriP: a method to predict antigenic epitopes using support vector machine to integrate tri-peptide similarity and propensity. PLoS ONE. 2012;7(9):e45152.
Ponomarenko J, Bui H-H, Li W, Fusseder N, Bourne PE, Sette A, Peters B. ElliPro: a new structure-based tool for the prediction of antibody epitopes. BMC Bioinform. 2008;9(1):514.
Wang P, Sidney J, Kim Y, Sette A, Lund O, Nielsen M, Peters B. Peptide binding predictions for HLA DR, DP and DQ molecules. BMC Bioinform. 2010;11(1):568.
Karosiene E, Lundegaard C, Lund O, Nielsen M. NetMHCcons: a consensus method for the major histocompatibility complex class I predictions. Immunogenetics. 2012;64(3):177–86.
Wang P, Sidney J, Dow C, Mothé B, Sette A, Peters B. A systematic assessment of MHC class II peptide binding predictions and evaluation of a consensus approach. PLoS Comput Biol. 2008;4(4):e1000048.
Jensen KK, Andreatta M, Marcatili P, Buus S, Greenbaum JA, Yan Z, Sette A, Peters B, Nielsen M. Improved methods for predicting peptide binding affinity to MHC class II molecules. Immunology. 2018;154(3):394–406.
Bhasin M, Raghava G. Prediction of CTL epitopes using QM SVM and ANN techniques. Vaccine. 2004;22(23–24):3195–204.
Dimitrov I, Flower DR, Doytchinova I. AllerTOP—a server for in silico prediction of allergens. BMC Bioinform BioMed Central. 2013;14:S4.
Magnan CN, Zeller M, Kayala MA, Vigil A, Randall A, Felgner PL, Baldi P. High-throughput prediction of protein antigenicity using protein microarray data. Bioinformatics. 2010;26(23):2936–43.
Doytchinova IA, Flower DR. VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinform. 2007;8(1):4.
Berzofsky JA. Immunogenicity and antigen structure. Fundamental immunology. 3rd ed. pp. 235–282; 1993.
Lee T-Y, Hsu JB-K, Chang W-C, Wang T-Y, Hsu P-C, Huang H-D. A comprehensive resource for integrating and displaying protein post-translational modifications. BMC Res Notes. 2009;2(1):111.
Shaddel M, Ebrahimi M, Tabandeh MR. Bioinformatics analysis of single and multi-hybrid epitopes of GRA-1, GRA-4, GRA-6 and GRA-7 proteins to improve DNA vaccine design against Toxoplasma gondii. J Parasit Dis. 2018;42(2):269–76.
El-Kady IM. T-cell immunity in human chronic toxoplasmosis. J Egypt Soc Parasitol. 2011;41(1):17–28.
Acknowledgements
Not applicable.
Funding
This study was financially supported by the Behbahan Faculty of Medical Sciences, Behbahan, Iran (Grant No. 99013).
Author information
Authors and Affiliations
Contributions
MF, SS, ADG, and HM conceived the study and designed the study protocol; MF performed the bioinformatics analysis; ADG rechecked the results; ADG and MF drafted the manuscript; MF, HM, AT, and MS critically revised the manuscript. All authors read and approved the final version of the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study received the approval from the Behbahan Faculty of Medical Sciences Ethical Committee (IR.BHN.REC.1399.034).
Consent for publication
Not applicable.
Competing interests
The authors declare that there is no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1
: Table S1. The acylation sites of CDPK4 sequence. Table S2. Epitopes predicted in CDPK4 protein by different parameters based on Bcepred online server. Table S3. The predicted B-cell epitopes via ABCpred tool. Table S4. Linear B-Cell epitope of the CDPK4 protein by SVMTriP. Table S5. IC50 values for CDPK4 binding to MHC class I molecules obtained using the IEDB. Table S6. IC50 values for CDPK4 binding to MHC class II molecules obtained using the IEDB. Table S7. Details of selected MHC-I T-cell epitope of T. gondii CDPK4 protein sequence using NetMHCcons server. Table S8. Details of selected MHC-II T-cell epitope of T. gondii CDPK4 protein sequence using NetMHCIIpan server. Table S9. Predicted CDPK4 epitopes by CTLpred.
Additional file 2
: Figure S1. Transmembrane domains expected in CDPK4 protein. (A) Some statistics and a list of the location of the predicted transmembrane helices and the predicted location of the intervening loop regions. Length: the length of the protein sequence; Number of predicted TMHs: The number of predicted transmembrane helices; Exp number of AAs in TMHs: The expected number of amino acids in transmembrane helices. If this number is larger than 18 it is very likely to be a transmembrane protein (OR have a signal peptide); Exp number, first 60 AAs: The expected number of amino acids in transmembrane helices in the first 60 amino acids of the protein. If this number more than a few, you should be warned that a predicted transmembrane helix in the N-term could be a signal peptide; Total prob of N-in: The total probability that the N-term is on the cytoplasmic side of the membrane; (B) Analysis of the transmembrane domains of CDPK4. Figure S2. (A) The results of the GOR4 server suggested that CDPK4 contains 34.97% alpha helix (Hh), 11.49% extended strand (Ee) and 53.54% random coils (Cc) in secondary structure; (B) Graphical finding from prediction of secondary structure of CDPK4 using GOR4. Figure S3. (A) The results of the SOPMA server suggested that CDPK4 contains 30.14% alpha helix (Hh), 10.02% extended strand (Ee) and 59.84% random coils (Cc) in secondary structure; (B) Graphical finding from prediction of secondary structure of CDPK4 using SOPMA server. Figure S4. Graphical output from prediction of secondary structure of CDPK4 using PSIPRED tool. Figure S5. Predicted 3D model by the SWISS-MODEL server. Figure S6. Propensity scale plots of CDPK4 protein. (A) Surface accessibility; (B) Antigenicity; (C) Bepipred linear epitope prediction; (D) Beta-turn; (E) Flexibility; (F) Hydrophilicity. On the graphs, the Y-axes indicate the corresponding score for each residue (averaged in the specified window), while the X-axes indicate the residue positions in the sequence. The higher residue score could be interpreted as having a higher likelihood that the residue would be part of the epitope (yellow color on the graphs). Green color (under the threshold) shows the unfavorable regions that are related to the properties of interest.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Foroutan, M., Ghaffari, A.D., Soltani, S. et al. Bioinformatics analysis of calcium-dependent protein kinase 4 (CDPK4) as Toxoplasma gondii vaccine target. BMC Res Notes 14, 50 (2021). https://doi.org/10.1186/s13104-021-05467-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13104-021-05467-1
Keywords
- Toxoplasma gondii
- Calcium-dependent protein kinase 4
- Bioinformatics analysis
- In silico
- Vaccine