https://doi.org/10.5812/archcid.69241

In Silico Prediction of T and B Cell Epitopes of SAG1-Related Sequence 3 (SRS3) Gene for Developing Toxoplasma gondii Vaccine

authors:

avatar Abolfazl Mirzadeh 1 , avatar Geita Saadatnia 1 , * , avatar Majid Golkar 2 , avatar Jalal Babaie 2 , 3 , avatar Samira Amiri 2 , avatar Asiyeh Yoosefy 2

Department of Biotechnology, Iranian Research Organization for Science and Technology (IROST), Tehran, Iran
Molecular Parasitology Laboratory, Department of Parasitology, Pasteur Institute of Iran, Tehran, Iran
Marand Health High Education Complex, Tabriz University of Medical Sciences, Marand, Iran
Archives of Clinical Infectious Diseases: Vol.15, issue 6; e69241
published online: December 26, 2020
article type: Research Article
received: April 10, 2018
revised: September 3, 2020
accepted: September 5, 2020

how to cite: Mirzadeh A, Saadatnia G, Golkar M, Babaie J, Amiri S, et al. In Silico Prediction of T and B Cell Epitopes of SAG1-Related Sequence 3 (SRS3) Gene for Developing Toxoplasma gondii Vaccine. Arch Clin Infect Dis. 2020;15(6):e69241. https://doi.org/10.5812/archcid.69241.

Abstract

Toxoplasmosis is a worldwide infection that can lead to serious problems in immune-compromised individuals and fetuses. A DNA vaccine strategy would be an ideal tool against Toxoplasma gondii. One of the necessary measures to provide an effective vaccine is the selection of proteins with high antigenicity. The SAG1-related sequence 3 (SRS3) protein is a major surface antigen in T. gondii that can be used as a vaccine candidate. In the present study, bioinformatics and computational methods were utilized to predict protein characteristics, as well as secondary and tertiary structures. The in silico approach is highly suited to analyze, design, and evaluate DNA vaccine strategies. Hence, in silico prediction was used to identify B and T cell epitopes and compare the antigenicity of SRS3 and other candidate genes of Toxoplasma previously applied in the production of vaccines. The results of the analysis theoretically showed that SRS3 has multiple epitopes with high antigenicity, proposing that SRS3 is a promising immunogenic candidate for the development of DNA vaccines against toxoplasmosis.

1. Background

Toxoplasma gondii is an obligate intracellular protozoan parasite of humans and other warm-blooded animals (1). This parasite is responsible for congenital infections and can lead to miscarriage, hydrocephalus, blindness, mental retardation, stillbirth, and spontaneous abortions (2-4). Tissue cysts remain largely inactive for the life of the host but can reactivate and cause severe disease in immunocompromised patients, such as those with HIV/AIDS (5, 6). None of the current drug treatments can affect tissue cysts. Due to the lack of an effective drug for toxoplasmosis, researchers strive to develop a vaccine against this disease (7).

Recent research has shown that DNA-based vaccines can draw out strong responses of Th1 type CD4+ T cells and CD8+ cytotoxic T cells. In addition, since all naturally occurring T and B cell epitopes are encoded for any protein by the DNA construct, it can induce both humoral and cellular immune responses. Hence, DNA vaccines could be a suitable option for effective protection against T. gondii infection compared to traditional vaccines (8, 9). Many reports have shown that some single antigens were successful as DNA vaccines such as SAG1, ROP8, ROP18, ROP54, GRA1, MIC3, TgCyP, and GRA6 (9-14). The most common method for the design of DNA vaccines is the selection of parasite proteins that are involved in the host cell invasion process and can stimulate an immunological response against toxoplasmosis (8).

Surface antigens like SAG1, SAG2, and SAG3 have been investigated widely in vaccine experiments (9, 15). Immunization with SAG1, SAG2A, and SAG3 has been tested in BALB/c mice and indicated to elicit Th1 cellular immune responses against all three antigens (11). One of the main surface antigens is SAG1-related sequence 3 (SRS3). It is a member of an SRS superfamily (glycosyl-phosphatidylinositol-linked proteins), which is structurally related to highly immunogenic surface antigen SAG1 with the highest overall amino acid identity among surface antigens of Toxoplasma (16).

A DNA vaccine based on SAG1 can protect against fatal infection and decrease the number of cysts. Besides, SAG1 and SRS3 are members of the same gene family and share high homology. In addition, mature SRS3 and SAG1 contain 12 cysteine residues and some identical short peptides (16-18). These findings suggest that they have similar folding, common overall topology, and the same function in stimulating the host system (16, 18, 19).

Since molecular experiments are highly expensive, bioinformatics approaches are needed for further justifying before conducting any laboratory work. Nowadays, advances in the in silico computational methods have been impressive in predicting the possible antigenic determinants. The in silico approach has numerous advantages to designing a safe, effective, and widespread vaccine against infections (20-25).

2. Objectives

In the present work, we applied bioinformatics tools with the in silico approach to predict physical and chemical properties, signal peptide, a transmembrane domain, secondary and tertiary structures, and B cell and T cell epitopes. Finally, the antigenicity of the SRS3 protein was evaluated by the Vaxigen server to prove protective immunity against T. gondii in humans.

3. Methods

3.1. Physical and Chemical Properties and Posttranslational Patterns and Motifs

The functions of all proteins depend on their structures that rely on physical and chemical parameters. Therefore, to create new drugs and vaccines, we need to study the chemical, classical, biological, physical, mathematical, and informatics properties of these molecules and work together in a new area known as bioinformatics. Physical and chemical properties, molecular weight, theoretical isoelectric point (pI), instability index, grand average of hydropathicity (GRAVY), and aliphatic index of the protein were analyzed by ProtParam. Motif analyses were done by PROSITE. The trans-membrane structure and subcellular location of SRS3 within eukaryotic cells were analyzed using servers listed in Table 1. Identifying the areas of signal peptides was made by signalP version 4.0 (26). Phylogenetic analysis was performed by multiple-alignment with ClustalW using GENEIOUS version 4.8.5 software.

Table 1.

Subcellular Localization Predictors and Their Results on SRS3 Sequence

ServerURLPredicted Localization
CELLOhttp://cello.life.nctu.edu.tw/Extracellular, cytoplasmic
WoLFPSORT http://wolfpsort.org/Extracellular
PSORTII http://psort.hgc.jp/form2.htmlCytoplasmic, nuclear
SLPFAhttp://sunflower.kuicr.kyoto-u.ac.jp/~tamura/slpfa.htmlExtracellular
SLP-Localhttp://sunflower.kuicr.kyoto-u.ac.jp/~smatsuda/slplocal.htmlExtracellular
SherLoc2http://abi.inf.uni-tuebingen.de/Services/SherLoc2Extracellular

3.2. Prediction of RNA Secondary Structure

The mRNA stability of the SRS3 construct and RNA secondary structure in accordance with the optimized gene were analyzed by Mfold web server at http://mfold.rna.albany.edu/?q=mfold/RNA-Folding-Form.

3.3. Prediction of Secondary and Tertiary Structures

The secondary structure of the SRS3 protein was predicted by SOPMA (27). The prediction of the tertiary structure, function, and folding of the SRS3 protein was carried out by the protein homology/analogy recognition engine (Phyre2) (28). Ligand-binding site prediction in SRS3 was done using 3DLigandSite (29). The validation of tertiary structure prediction and modeling was performed using ProSA-web (30).

3.4. Prediction of T and B Cell Epitopes for SRS3 Protein and Comparison of Antigenicity with Applied Genes in Vaccines

The T cell epitopes were predicted using several prediction online servers, such as Immune Epitope Database, SYFPEITHI, ProPred, BIMAS, SVMHC, NetMHC, EpiJen, Rankpep, NetCTL1.2.

The protein sequence of SRS3 was analyzed for the prediction of B cell epitopes using the Bioinformatics Predicted Antigenic Peptides (BPAP) system, ABCPred, Bcepred 1.0 server, and BepiPred at http://www.cbs.dtu.dk/services/BepiPred/.

The final B and T cell predicted epitopes were evaluated using the VaxiJen 2.0 server, which is used to evaluate the final B and T cell predicted epitopes (31, 32). Moreover, Vaxijen was utilized for the antigenicity prediction of SRS3 and comparison with applied genes in vaccines.

4. Results

4.1. Physical and Chemical Characteristics Study of tSRS3 Using Protparam Program

The SRS3 protein contains 345 amino acid residues. The molecular weight of SRS3 was 35 kDa and the acidity of the protein was indicated by the pI value. Theoretical PI was 5.76. The instability index was obtained as 43.17. The Grand Average of Hydropathicity (GRAVY) and aliphatic index of the SRS3 protein were 0.140 and 83.36, respectively.

4.2. Prediction of Modification Sites

We used PROSITE to analyze the motif. We found 12 hits on the SRS3 sequence, suggesting that there will be 12 potential Posttranslational Modification (PTM) sites on the sequence. They included four N-myristoylation sites, four protein kinase C phosphorylation sites, three casein kinase II phosphorylation sites, and one N-glycosylation site.

4.3. Prediction of Trans-Membrane Structure and Sub-Cellular Localization

Predicting that a protein has trans-membrane segments gives more information about the function of the protein. In this study, five out of six advanced servers predicted extracellular localization for the protein sequence expressed in the eukaryotic cell program. The CELLO and PSORT II servers predicted two localizations for the hypothetical protein. The results showed that SRS3 is an outside protein (Table 1). Signal peptides and N-terminal trans-membrane helices are hydrophobic, but trans-membrane helices typically have longer hydrophobic regions. Furthermore, trans-membrane helices do not have cleavage sites; however, the cleavage-site pattern is not sufficient to distinguish the two types of sequence. Therefore, signalP (version 4.0) was used to discriminate between signal peptides and transmembrane regions. The results showed that the score of predicted signal peptide sequence was 0.329, and the cutoff was predicted to be 0.5, indicating that there was no signal peptide sequence and signal peptidase cleavage site on the amino acid sequence of SRS3 (Figure 1).

Signal peptide predicted within the SRS3 sequence. The C-score is high at the position immediately after the cleavage site. The S-score distinguishes between positions within signal peptides from those of the mature part of the proteins. The Y-score is a geometric average of the C-score and the slope of the S-score.
Signal peptide predicted within the SRS3 sequence. The C-score is high at the position immediately after the cleavage site. The S-score distinguishes between positions within signal peptides from those of the mature part of the proteins. The Y-score is a geometric average of the C-score and the slope of the S-score.
Figure 1.

Signal peptide predicted within the SRS3 sequence. The C-score is high at the position immediately after the cleavage site. The S-score distinguishes between positions within signal peptides from those of the mature part of the proteins. The Y-score is a geometric average of the C-score and the slope of the S-score.

4.4. Phylogenetic Analysis of SRS3

Phylogenetic analyses indicated that the SRS3 gene expressed in stage-specific tachyzoite is the closest paralog to SAG1 and BSR4, as shown in Figure 2. In mRNA stability analyses, 28 secondary structures were anticipated for optimized SRS3 mRNA, among which, a model with ΔG = -248.60 kcal/mol possessed the best-predicted stability. The value of the total ΔG of the two base pairs was -0.70 kcal/mol. Among the 28 models, just two models formed a single-stranded structure at the first 20 nucleotides of the 5’ end. The calculated ΔG of these two models was -247.50 and -247.70 kcal/mol. The best structure predicted on native mRNA had ΔG of -240.30 kcal/mol.

Evolutionary relationship between representative SRS family members from the life cycle
Evolutionary relationship between representative SRS family members from the life cycle
Figure 2.

Evolutionary relationship between representative SRS family members from the life cycle

4.5. Prediction of Secondary and Tertiary Structures of tSRS3

The predicted secondary structure of the tSRS3 protein is shown in Figure 3. The results revealed that tSRS3 is a coil structure-based protein. The rates of the α-helix, β-Turn, extended strand, and random coil were 13.74%, 7.22%, 31.30%, and 47.71%. The functions of proteins depend mainly on their structures. We predicted the protein fold recognition and 3D-structure related to tSRS3 using the Phyre2 server based on the protein primary structure. The results of the tertiary structure of tSRS3 showed a protein with two main domains. The confidence score (C-score) for estimating the quality of predicted models by Phyre2 was 100%, as shown in Figure 4A. The prediction of ligand-binding sites in tSRS3 was obtained using 3DLigandSite. The tertiary model of the predicted tSRS3 protein was examined using the more sensitive structure alignment program “Matching Molecular Models Obtained from Theory (MAMMOTH)”. The result identified 10 ligand clusters; among them, cluster 10 was the most significant predictor of one ligand called ZINC ION (ZN) with an average mammoth score of 30.8, as shown in Figure 4B. The residues predicted to form part of the binding site are shown in shape B (blue glob). Amino acids and their residue numbers include VAL (96), THR (97), LEU (98), PRO (99), GLY (170), CYS (171), and THR (172). Besides, ProSA was used to check 3D models of the SRS3 structure for potential errors. The z score indicates the overall model quality. The z score of -5.7 was calculated for the input structure that is within the range of scores generally found for native proteins of similar size, as shown in Figure 4C.

The SRS3 secondary structure prediction of alpha helix (Hh), extended strand (Ee), beta turn (Tt), and random coil (Cc)
The SRS3 secondary structure prediction of alpha helix (Hh), extended strand (Ee), beta turn (Tt), and random coil (Cc)
Figure 3.

The SRS3 secondary structure prediction of alpha helix (Hh), extended strand (Ee), beta turn (Tt), and random coil (Cc)

A, The tertiary structure prediction of SRS3 by Phyre2; 324 residues (98.7% of SRS3 sequences) were modeled with 100.0% confidence by the single highest scoring template (BSR4), a, α-helix; b, random coil; c, β-folded; B, ligand-binding site results for the SRS3 protein. Blue colors list amino acid residues observed in cluster 10 of the predicted binding site; C, ProSA-web z-score SRS3 protein plot. The z-score indicates overall model quality. The Z-scores of SRS3 protein chain A in PDB were determined by X-ray crystallography (Pale blue) or NMR spectroscopy (bold blue) concerning their lengths. The plot shows results with a z-score of ≤ 10. The z-score of SRS3 is -5.7 and is shown as a large dot. The value of -5.7 is in the range of native conformations.
A, The tertiary structure prediction of SRS3 by Phyre2; 324 residues (98.7% of SRS3 sequences) were modeled with 100.0% confidence by the single highest scoring template (BSR4), a, α-helix; b, random coil; c, β-folded; B, ligand-binding site results for the SRS3 protein. Blue colors list amino acid residues observed in cluster 10 of the predicted binding site; C, ProSA-web z-score SRS3 protein plot. The z-score indicates overall model quality. The Z-scores of SRS3 protein chain A in PDB were determined by X-ray crystallography (Pale blue) or NMR spectroscopy (bold blue) concerning their lengths. The plot shows results with a z-score of ≤ 10. The z-score of SRS3 is -5.7 and is shown as a large dot. The value of -5.7 is in the range of native conformations.
Figure 4.

A, The tertiary structure prediction of SRS3 by Phyre2; 324 residues (98.7% of SRS3 sequences) were modeled with 100.0% confidence by the single highest scoring template (BSR4), a, α-helix; b, random coil; c, β-folded; B, ligand-binding site results for the SRS3 protein. Blue colors list amino acid residues observed in cluster 10 of the predicted binding site; C, ProSA-web z-score SRS3 protein plot. The z-score indicates overall model quality. The Z-scores of SRS3 protein chain A in PDB were determined by X-ray crystallography (Pale blue) or NMR spectroscopy (bold blue) concerning their lengths. The plot shows results with a z-score of ≤ 10. The z-score of SRS3 is -5.7 and is shown as a large dot. The value of -5.7 is in the range of native conformations.

4.6. In Silico Prediction of T Cell Epitopes for SRS3 Antigen

Several online prediction software was utilized to identify HLA-A*1101 and HLA A*0201 restricted epitopes from T. gondii SRS3 with the highest numbers (Table 2). The results of antigenicity among the outstanding epitopes of Toxoplasma antigens by Vaxijen showed that SRS3 epitopes were highly antigenic (Table 3).

Table 2.

In Silico Prediction of CD8+ Final T Cell Epitopes Related to SRS3 Protein, Based on Their Predicted Binding Affinity to HLA-A*02 and HLA-A*03 Supertype Moleculesa

MHC-I AllelePositionPeptide SequencesServerVaxijen Scores
HLA-A*1101167 - 175IVGCTQEGKBIMAS1.7863
HLA-A*1101167 - 175IVGCTQEGKProPred-I
HLA-A*1101167 - 175IVGCTQEGKEpiJen
HLA-A*1101247 - 256LQENCDEKYRANKPEP0.985
HLA-A*1101269 - 277VSQAGTSNKConsensus (ann/smm)0.5653
HLA-A*1101246 - 256PLQENCDEKSVM1.1607
HLA-A*1101144 - 153DSQTVTCAYNetCTL1.20.864
HLA-A*020115 - 24STLFFVALLAConsensus1.6329
HLA-A*020116 - 23LFFVALLASVMHC2.0760
HLA-A*0201112 - 121GMTDPQSCVBIMAS1.2069
HLA-A*0201179 - 188CLLHVTLSAEpiJen0.8877

aAntigenicity of full-length protein epitopes was calculated using VaxiJen.

Table 3.

Comparison of Antigenicity of Peptide Candidates Utilized for Screening of CD8+ T Cells with SRS3a

HLA-AAntigenPeptide SequencesPositionVaxijen Scores
HLA-A*1101SAG1KSFKDILPK*224 - 2320.4842
HLA-A*1101GRA6AMLTAFFLR*164 - 1721.9220
HLA-A*1101GRA7RSFKDLLKK*134 - 142-0.2503
HLA-A*1101SPASSAYVFSVK*250 - 2580.8154
HLA-A*1101SAG2CSTFWPCLLR*13 - 210.8148
HLA-A*1101SAG2DALLEHVFIEK63 - 720.9153
HLA-A*1101BSR4TTRFVEIFPK284 - 2931.0332
HLA-A*1101SAG3VVGHVTLNK80 - 88-0.4365
HLA-A*1101SRS3IVGCTQEGK**167 - 1751.7863
HLA-A*0201GRA6FMGVLVNSL*29 - 371.9980
HLA-A*0201SPAGLAAAVVAV*82 - 900.5290
HLA-A*0201GRA3FLVPFVVFL*25 - 332.2250
HLA-A*0201SAG2CFLSLSLLVI*38 - 460.7125
HLA-A*0201SAG2DFMIAFISCFA*180 - 1891.1592
HLA-A*0201SAG2xFMIVSISLV*351 - 3590.1095
HLA-A*0201SAG3FLTDYIPGA*136 - 1441.4166
HLA-A*0201SRS3LFFVALLA**16 - 232.0760

aCandidate epitopes were used in the vaccine * and the elected epitope of SRS3 protein **. Antigenicity of full-length protein epitopes was calculated using VaxiJen.

4.7. In Silico Prediction of B Cell Epitopes for SRS3 Antigen and Evaluation of Antigenicity

The B cell linear epitopes were predicted using BcePred, BAP, and ABCpred servers. The best epitopes were selected based on physicochemical properties of amino acids (accessibility, hydrophilicity, flexibility, polarity, exposed surface, and turns and strong antigenicity) (Table 4). In general, the final epitopes of SRS3 had BCpreds, ABCpred, BepiPred, and VaxiJen cutoff values of > 0.8, > 0.8, and > 0.5, respectively. The results of the evaluation of SRS3 protein antigenicity and comparison with Toxoplasma antigens used in DNA vaccines are shown in Table 5.

Table 4.

Final B Cell Linear Epitopes Based on Ability of Antigenicity. Antigenicity of Full-Length Protein Epitopes Was Calculated Using VaxiJen

B Cell EpitopePrediction ServerAmino Acid PositionsScoresVaxijen Scores
HRVTVSPKDPTFTLDCSSDGBCPred2100.9981.2979
VCPAGMTDPQSCVSPSGTIPBCPred1080.9750.9927
GSCTSQGTSSPMQVCTCIVPBCPred490.9550.8694
MQPASFTTSYCPABCPred2310.8790.7269
TGPLQENCDEKYASILABCpred2440.721.3467
TLSARASSKDSQTVTCABCpred1840.741.3257
SLLSGDGASAVAWVGDABCpred1230.691.1767
PEAQKKVLVGCKASTAABCpred2880.721.0681
GETSNHAVVVSAEKNQABCpred710.870.9510
ELTITSSSPSSAWFLSABCpred3130.730.9448
GTIPFASLLSGDGASAABCpred1240.770.8582
TLDCSSDBAP222High1.7363
WWQVSQABAP266High1.6518
NHAVVVSAEBAP75High1.2996
QKKVLVGCKASTABAP291High1.1761
ASAVAWVGDEBAP130High1.1352
PQSCVSPSGTIPFASLLSGBAP116High0.9466
HFAEAAGHQGSCTSQGTSSPMBepiPred40High1.084
SLLSGDGASAVAWVGDEVNBepiPred130High1.046
Table 5.

Evaluation of Antigenicity and Comparision With Toxoplasma gondii Antigenes by Vaxijena

AntigenEvaluation of Antigenicity by Vaxijen
SAG10.7846
SRS30.6919
SAG2A0.5861
SAG30.6270
SAG5A0.6728
ROP10.5595
ROP20.5905
ROP80.6753
GRA1b0.4926
GRA20.6902
GRA6c0.8949
GRA70.5553
Tgcyp0.6872
TGME490.5526

aThe server assessment threshold is 0.5.

bNon-antigen.

cHigh score = High antigenicity.

5. Discussion

Vaccination with DNA vaccines against toxoplasmosis is more effective than other vaccines in inducing the activation and proliferation of CD4+ and CD8+ T cells and eliciting specific antibodies essential to control chronic infection (9). The identification and verification of immunodominant recombinant antigens is the best way to achieve effective vaccines against toxoplasmosis (33). The development of DNA vaccines is a complicated process that needs molecular biology and immunology (34). In recent years, computational methods along with in silico studies compared to wet-lab methods have shown strong potential for the development of effective vaccines (35, 36). As a crucial step in the preparation of vaccines and antibodies, we require the necessary data for the identification of B and T cell epitopes. Therefore, in the present study, SRS3 was selected as a candidate vaccine and examined by the comprehensive analysis of bioinformatics with the in silico approach for the first time.

As known, SRS3 is one of the main surface antigens of tachyzoites that plays a key role in the process of invasion to the host cell (16, 37). Expression analysis using quantitative real-time PCR showed that the SRS3 gene was up-regulated in vivo relative to their expression levels in vitro (37). In addition, this antigen creates a very strong reactivity with IgG and IgM antibodies against Toxoplasma infection (37). The objectives of this study were to assess the antigenicity of SRS3 as a vaccine candidate and predict the B and T cell epitopes of tSRS3. Hence, to better understand the structure and function of SRS3, the basic sequence features were analyzed and the secondary and tertiary structures of SRS3 were studied. Protein post-translational modification (PTM) plays a basic role in cellular control mechanisms (38). We found 12 PTM sites in SRS3 including four N-myristoylation sites, four protein kinase C phosphorylation, three casein kinase II phosphorylation sites, and one N-glycosylation site. The subcellular localization prediction indicated that SRS3 is an extracellular protein. In addition, analysis of the signal peptide cleavage site by signal P (version 4.0) showed that SRS3 lacks signal peptides, contrary to the prediction made by Manger et al. (16). However, the nascent polypeptide related to a primary SRS3 mRNA transcript contains hydrophobic regions in N- and C-terminals that are not present in mature SRS3 proteins on the surface of T. gondii. The results of phylogenetic analyses displayed that this gene has the most conservation in the sequence of proteins to SAG1 and BSR4 among surface antigens of Toxoplasma, as shown in Figure 2.

The secondary structure of SRS3 represented that the rate of the random coil was above 47%. Thus, SRS3 is a coil structure-based protein that leads to an increased likelihood of forming an antigenic epitope. The third structure of SRS3 was carried out by homology modeling (Phyre 2), which has been widely used in many areas of structure-based analysis and showed high similarity with BSR4. Thus, SRS3 has a folding similar to BSR4 and SAG1. Therefore, SRS3 plays a functional role like SAG1 and BSR4 in invading the host cell. Ligand-binding site prediction identified one ligand called ZINC ION (ZN) with an average mammoth score of 30.8. The ideal DNA vaccine to protect against toxoplasmosis in humans should be able to elicit a protective T helper cell type 1 immune response and generate long-lived IFN-γ-producing CD8+ T cells. Furthermore, the vaccine should be able to induce a strong B cell and T cell responses. Hence, the identification of the precise location of T cell and B cell epitopes in the protein is important for vaccine design and development. In a previous study, T. gondii proteins were analyzed to recognize CD8+ T cell epitopes that were based on their predicted binding affinity to HLA-A*02 and HLA-A*03 supertype molecules (39, 40). Therefore, in the present study, we used the most important computational prediction tools for the identification of T. gondii SRS3 protein HLA-A*0201 and HLA-A*1101-restricted epitopes. Several potential T cell epitopes on SRS3 were determined with high-affinity binding to HLA-A*02 and HLA-A*03 supertype molecules (Table 2). In addition, from an antigenic point of view, the comparison with the identified peptides of Toxoplasma antigens showed that the peptide of SRS3 is extremely antigenic (Table 3). The results of linear B cell epitope prediction showed 17 B-cells for the SRS3 protein (Table 4).

Vaxijen was used to select the final B and T cell predicted epitopes and evaluate SRS3 antigenicity (Table 5). The results revealed that SRS3 antigenicity could stimulate the immune system against toxoplasmosis more appropriately than did other T. gondii antigens that were used to construct a DNA vaccine (9, 10, 41-43). The results of the analysis theoretically determined that SRS3 is a potent antigen of strong antigenicity and has multiple epitopes with high antigenicity. Thus, SRS3 may be a promising immunogenic candidate for the development of DNA vaccines against toxoplasmosis.

5.1. Conclusions

In this study, we used bioinformatics methods with an in silico approach to predict protein characteristics, secondary and tertiary structures, and functions and finally confirmed the high antigenicity of SRS3 by the Vaxigen server for the first time. This study provides an appropriate basis for further studies on SRS3 in vaccine research.

References

  • 1.

    Dubey J. Louis M, Kami K, editors. History and Life Cycle of Toxoplasma gondii. Toxoplasma gondii: the model apicomplexan: perspectives and methods. 2006/11/30 ed. Elsevier; 2007. https://doi.org/10.1016/B978-012369542-0/50003-9.

  • 2.

    Olariu TR, Remington JS, McLeod R, Alam A, Montoya JG. Severe congenital toxoplasmosis in the United States: clinical and serologic findings in untreated infants. Pediatr Infect Dis J. 2011;30(12):1056-61. [PubMed ID: 21956696]. https://doi.org/10.1097/INF.0b013e3182343096.

  • 3.

    Kravetz JD, Federman DG. Toxoplasmosis in pregnancy. Am J Med. 2005;118(3):212-6. [PubMed ID: 15745715]. https://doi.org/10.1016/j.amjmed.2004.08.023.

  • 4.

    Wong SY, Remington JS. Toxoplasmosis in pregnancy. Clin Infect Dis. 1994;18(6):853-61. quiz 862. [PubMed ID: 8086543]. https://doi.org/10.1093/clinids/18.6.853.

  • 5.

    Melzer TC, Cranston HJ, Weiss LM, Halonen SK. Host Cell Preference of Toxoplasma gondii Cysts in Murine Brain: A Confocal Study. J Neuroparasitology. 2010;1. [PubMed ID: 21625284]. [PubMed Central ID: PMC3103221]. https://doi.org/10.4303/jnp/N100505.

  • 6.

    Suzuki Y, Sa Q, Gehman M, Ochiai E. Interferon-gamma- and perforin-mediated immune responses for resistance against Toxoplasma gondii in the brain. Expert Rev Mol Med. 2011;13. e31. [PubMed ID: 22005272]. [PubMed Central ID: PMC3372998]. https://doi.org/10.1017/S1462399411002018.

  • 7.

    Suzuki Y, Wang X, Jortner BS, Payne L, Ni Y, Michie SA, et al. Removal of Toxoplasma gondii cysts from the brain by perforin-mediated activity of CD8+ T cells. Am J Pathol. 2010;176(4):1607-13. [PubMed ID: 20167872]. [PubMed Central ID: PMC2843452]. https://doi.org/10.2353/ajpath.2010.090825.

  • 8.

    Liu Q, Singla LD, Zhou H. Vaccines against Toxoplasma gondii: status, challenges and future directions. Hum Vaccin Immunother. 2012;8(9):1305-8. [PubMed ID: 22906945]. [PubMed Central ID: PMC3579912]. https://doi.org/10.4161/hv.21006.

  • 9.

    Lim SS, Othman RY. Recent advances in Toxoplasma gondii immunotherapeutics. Korean J Parasitol. 2014;52(6):581-93. [PubMed ID: 25548409]. [PubMed Central ID: PMC4277020]. https://doi.org/10.3347/kjp.2014.52.6.581.

  • 10.

    Sun XM, Zou J, A A, Yan WC, Liu XY, Suo X, et al. DNA vaccination with a gene encoding Toxoplasma gondii GRA6 induces partial protection against toxoplasmosis in BALB/c mice. Parasit Vectors. 2011;4:213. [PubMed ID: 22070984]. [PubMed Central ID: PMC3229464]. https://doi.org/10.1186/1756-3305-4-213.

  • 11.

    Jongert E, Roberts CW, Gargano N, Forster-Waldl E, Petersen E. Vaccines against Toxoplasma gondii: challenges and opportunities. Mem Inst Oswaldo Cruz. 2009;104(2):252-66. [PubMed ID: 19430651]. https://doi.org/10.1590/s0074-02762009000200019.

  • 12.

    Yuan ZG, Zhang XX, Lin RQ, Petersen E, He S, Yu M, et al. Protective effect against toxoplasmosis in mice induced by DNA immunization with gene encoding Toxoplasma gondii ROP18. Vaccine. 2011;29(38):6614-9. [PubMed ID: 21762755]. https://doi.org/10.1016/j.vaccine.2011.06.110.

  • 13.

    Gong P, Cao L, Guo Y, Dong H, Yuan S, Yao X, et al. Toxoplasma gondii: Protective immunity induced by a DNA vaccine expressing GRA1 and MIC3 against toxoplasmosis in BALB/c mice. Exp Parasitol. 2016;166:131-6. [PubMed ID: 27059254]. https://doi.org/10.1016/j.exppara.2016.04.003.

  • 14.

    Yang WB, Zhou DH, Zou Y, Chen K, Liu Q, Wang JL, et al. Vaccination with a DNA vaccine encoding Toxoplasma gondii ROP54 induces protective immunity against toxoplasmosis in mice. Acta Trop. 2017;176:427-32. [PubMed ID: 28935555]. https://doi.org/10.1016/j.actatropica.2017.09.007.

  • 15.

    Chuang SC, Ko JC, Chen CP, Du JT, Yang CD. Induction of long-lasting protective immunity against Toxoplasma gondii in BALB/c mice by recombinant surface antigen 1 protein encapsulated in poly (lactide-co-glycolide) microparticles. Parasit Vectors. 2013;6:34. [PubMed ID: 23398973]. [PubMed Central ID: PMC3584932]. https://doi.org/10.1186/1756-3305-6-34.

  • 16.

    Manger ID, Hehl AB, Boothroyd JC. The surface of Toxoplasma tachyzoites is dominated by a family of glycosylphosphatidylinositol-anchored antigens related to SAG1. Infect Immun. 1998;66(5):2237-44. [PubMed ID: 9573113]. [PubMed Central ID: PMC108187]. https://doi.org/10.1128/IAI.66.5.2237-2244.1998.

  • 17.

    Radke JR, Gubbels MJ, Jerome ME, Radke JB, Striepen B, White MW. Identification of a sporozoite-specific member of the Toxoplasma SAG superfamily via genetic complementation. Mol Microbiol. 2004;52(1):93-105. [PubMed ID: 15049813]. https://doi.org/10.1111/j.1365-2958.2003.03967.x.

  • 18.

    Mirzadeh A, Saadatnia G, Golkar M, Babaie J, Noordin R. Production of refolded Toxoplasma gondii recombinant SAG1-related sequence 3 (SRS3) and its use for serodiagnosis of human toxoplasmosis. Protein Expr Purif. 2017;133:66-74. [PubMed ID: 28263855]. https://doi.org/10.1016/j.pep.2017.03.001.

  • 19.

    Mirzadeh A, Valadkhani Z, Yoosefy A, Babaie J, Golkar M, Esmaeili Rastaghi AR, et al. Expression, purification and in vitro refolding of the recombinant truncated Saposin-like protein 2 antigen for development of diagnosis of human fascioliasis. Acta Trop. 2017;171:163-71. [PubMed ID: 28300559]. https://doi.org/10.1016/j.actatropica.2017.03.002.

  • 20.

    Moradi J, Tabrizi M, Izad M, Mosavari N, Feizabadi MM. Designing a Novel Multi-epitope DNA- Based Vaccine Against Tuberculosis: In Silico Approach. Jundishapur J Microbiol. 2017;10(3). https://doi.org/10.5812/jjm.43950.

  • 21.

    Farhadi T, Ranjbar MM. Designing and modeling of complex DNA vaccine based on MOMP of Chlamydia trachomatis: an in silico approach. Network Model Anal Health Inf Bioinf. 2016;6(1). https://doi.org/10.1007/s13721-016-0142-5.

  • 22.

    Jahangiri A, Rasooli I, Gargari SL, Owlia P, Rahbar MR, Amani J, et al. An in silico DNA vaccine against Listeria monocytogenes. Vaccine. 2011;29(40):6948-58. [PubMed ID: 21791233]. https://doi.org/10.1016/j.vaccine.2011.07.040.

  • 23.

    Jahangiri A, Amani J, Halabian R, Imani fooladi AA. In Silico Analyses of Staphylococcal Enterotoxin B as a DNA Vaccine for Cancer Therapy. Int J Peptide Res Ther. 2017;24(1):131-42. https://doi.org/10.1007/s10989-017-9595-3.

  • 24.

    Jahangiri A, Rasooli I, Owlia P, Fooladi AA, Salimian J. In silico design of an immunogen against Acinetobacter baumannii based on a novel model for native structure of Outer membrane protein A. Microb Pathog. 2017;105:201-10. [PubMed ID: 28242426]. https://doi.org/10.1016/j.micpath.2017.02.028.

  • 25.

    Toussaint NC, Kohlbacher O. Towards in silico design of epitope-based vaccines. Expert Opin Drug Discov. 2009;4(10):1047-60. [PubMed ID: 23480396]. https://doi.org/10.1517/17460440903242283.

  • 26.

    Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785-6. [PubMed ID: 21959131]. https://doi.org/10.1038/nmeth.1701.

  • 27.

    Geourjon C, Deleage G. SOPMA: significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments. Comput Appl Biosci. 1995;11(6):681-4. [PubMed ID: 8808585]. https://doi.org/10.1093/bioinformatics/11.6.681.

  • 28.

    Kelley LA, Sternberg MJ. Protein structure prediction on the Web: a case study using the Phyre server. Nat Protoc. 2009;4(3):363-71. [PubMed ID: 19247286]. https://doi.org/10.1038/nprot.2009.2.

  • 29.

    Wass MN, Kelley LA, Sternberg MJ. 3DLigandSite: predicting ligand-binding sites using similar structures. Nucleic Acids Res. 2010;38(Web Server issue):W469-73. [PubMed ID: 20513649]. [PubMed Central ID: PMC2896164]. https://doi.org/10.1093/nar/gkq406.

  • 30.

    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(Web Server issue):W407-10. [PubMed ID: 17517781]. [PubMed Central ID: PMC1933241]. https://doi.org/10.1093/nar/gkm290.

  • 31.

    Soria-Guerra RE, Nieto-Gomez R, Govea-Alonso DO, Rosales-Mendoza S. An overview of bioinformatics tools for epitope prediction: implications on vaccine development. J Biomed Inform. 2015;53:405-14. [PubMed ID: 25464113]. https://doi.org/10.1016/j.jbi.2014.11.003.

  • 32.

    Ahmad TA, Eweida AE, Sheweita SA. B-cell epitope mapping for the design of vaccines and effective diagnostics. Trials Vaccinol. 2016;5:71-83. https://doi.org/10.1016/j.trivac.2016.04.003.

  • 33.

    Wang HL, Zhang TE, Yin LT, Pang M, Guan L, Liu HL, et al. Partial protective effect of intranasal immunization with recombinant Toxoplasma gondii rhoptry protein 17 against toxoplasmosis in mice. PLoS One. 2014;9(9). e108377. [PubMed ID: 25255141]. [PubMed Central ID: PMC4177930]. https://doi.org/10.1371/journal.pone.0108377.

  • 34.

    Nascimento IP, Leite LC. Recombinant vaccines and the development of new vaccine strategies. Braz J Med Biol Res. 2012;45(12):1102-11. [PubMed ID: 22948379]. [PubMed Central ID: PMC3854212]. https://doi.org/10.1590/s0100-879x2012007500142.

  • 35.

    Flower DR, Davies MN, Doytchinova IA. Identification of Candidate Vaccine Antigens In Silico. Immunomic Discovery of Adjuvants and Candidate Subunit Vaccines. 2013. p. 39-71. https://doi.org/10.1007/978-1-4614-5070-2_3.

  • 36.

    Brusic V, Bajic VB, Petrovsky N. Computational methods for prediction of T-cell epitopes--a framework for modelling, testing, and applications. Methods. 2004;34(4):436-43. [PubMed ID: 15542369]. https://doi.org/10.1016/j.ymeth.2004.06.006.

  • 37.

    Amerizadeh A, Idris ZM, Khoo BY, Kotresha D, Yunus MH, Karim IZ, et al. Identification of Toxoplasma gondii in-vivo induced antigens by cDNA library immunoscreening with chronic toxoplasmosis sera. Microb Pathog. 2013;54:60-6. [PubMed ID: 23044055]. https://doi.org/10.1016/j.micpath.2012.09.006.

  • 38.

    Zhao Y, Jensen ON. Modification-specific proteomics: strategies for characterization of post-translational modifications using enrichment techniques. Proteomics. 2009;9(20):4632-41. [PubMed ID: 19743430]. [PubMed Central ID: PMC2892724]. https://doi.org/10.1002/pmic.200900398.

  • 39.

    Cong H, Mui EJ, Witola WH, Sidney J, Alexander J, Sette A, et al. Towards an immunosense vaccine to prevent toxoplasmosis: protective Toxoplasma gondii epitopes restricted by HLA-A*0201. Vaccine. 2011;29(4):754-62. [PubMed ID: 21095258]. [PubMed Central ID: PMC3014376]. https://doi.org/10.1016/j.vaccine.2010.11.015.

  • 40.

    Cong H, Mui EJ, Witola WH, Sidney J, Alexander J, Sette A, et al. Human immunome, bioinformatic analyses using HLA supermotifs and the parasite genome, binding assays, studies of human T cell responses, and immunization of HLA-A*1101 transgenic mice including novel adjuvants provide a foundation for HLA-A03 restricted CD8+T cell epitope based, adjuvanted vaccine protective against Toxoplasma gondii. Immunome Res. 2010;6:12. [PubMed ID: 21129215]. [PubMed Central ID: PMC3009956]. https://doi.org/10.1186/1745-7580-6-12.

  • 41.

    Parthasarathy S, Fong MY, Ramaswamy K, Lau YL. Protective immune response in BALB/c mice induced by DNA vaccine of the ROP8 gene of Toxoplasma gondii. Am J Trop Med Hyg. 2013;88(5):883-7. [PubMed ID: 23509124]. [PubMed Central ID: PMC3752752]. https://doi.org/10.4269/ajtmh.12-0727.

  • 42.

    Lu G, Wang L, Zhou A, Han Y, Guo J, Song P, et al. Epitope analysis, expression and protection of SAG5A vaccine against Toxoplasma gondii. Acta Trop. 2015;146:66-72. [PubMed ID: 25792417]. https://doi.org/10.1016/j.actatropica.2015.03.013.

  • 43.

    Bai Y, He S, Zhao G, Chen L, Shi N, Zhou H, et al. Toxoplasma gondii: bioinformatics analysis, cloning and expression of a novel protein TgIMP1. Exp Parasitol. 2012;132(4):458-64. [PubMed ID: 23026454]. https://doi.org/10.1016/j.exppara.2012.09.015.