The rationale for targeting dense granule antigens, or GRAs, rests on their central role in the parasite’s survival strategy. Once T. gondii invades a host cell, it seals itself inside a parasitophorous vacuole and secretes more than 70 GRA proteins that decorate the vacuole membrane, the intravacuolar network and even the host cell nucleus. These effectors allow the parasite to manipulate host cell function and evade destruction. The research team focused on six GRAs selected for their recent discovery and documented roles in virulence or immunogenicity: GRA15, a well-characterised immunogen that activates the NF-κB pathway; GRA60, a virulence factor acting alongside ROP18; GRA76, identified in 2024 as a growth factor; GRA83, which provokes innate immune responses; and two novel effectors, GRA-α and GRA-β, shown to suppress parasite replication and cyst formation.

No approved vaccine exists for human toxoplasmosis. The only licensed product, the live-attenuated Toxovax used in sheep, cannot be given to people because of its residual pathogenic potential. Existing drug therapies for the acute tachyzoite stage often carry side effects and do not eliminate the dormant tissue cysts that cause lifelong infection. This therapeutic gap has pushed researchers toward reverse vaccinology, the computational discipline that screens pathogen genomes for the most promising antigenic fragments before any experiment begins. By filtering candidates in silico first, laboratories can save years of work, reduce costs and avoid animal testing of constructs destined to fail.

The epitope-mapping stage was deliberately stringent. The team retrieved the amino acid sequences of the six GRAs from the ToxoDB database and predicted linear B-cell epitopes using three independent web servers, retaining only fragments identified by at least two tools. Each surviving epitope was then scored for antigenicity with VaxiJen, screened for allergenicity with AllergenFP and checked for water solubility. In parallel, the researchers predicted cytotoxic T-lymphocyte epitopes of nine or ten residues bound by MHC class I molecules, and helper T-lymphocyte epitopes of fifteen residues bound by MHC class II, selecting the strongest allele-epitope combinations by percentile rank. CTL candidates underwent additional testing for immunogenicity and interferon-gamma induction potential, while HTL candidates were screened for toxicity as well.

From this filtering emerged eighteen final epitopes: six B-cell, six CTL and six HTL segments drawn from the six GRAs. The team then assembled four vaccine constructs, each combining the same epitope payload with a different immunostimulatory adjuvant. RS09-Toxo carried the synthetic TLR-4 agonist peptide RS-09; RSFN-Toxo was double-adjuvanted with RS-09 and human interferon-gamma; 50S-Toxo used a Mycobacterium tuberculosis ribosomal protein; and Rpf-Toxo employed the M. tuberculosis resuscitation-promoting factor E, RpfE. Flexible glycine-serine linkers, proteasomal cleavage sites and the rigid EAAAK spacer held the domains apart, and a C-terminal histidine tag enabled purification. The adjuvant choices were not arbitrary: all favour the Th1-type immunity that is critical for clearing intracellular pathogens such as T. gondii.

Biochemical profiling separated the field quickly. All four constructs scored highly for antigenicity, with values ranging from 0.9148 for Rpf-Toxo to 2.1188 for RS09-Toxo, and none was predicted to be allergenic. However, the ProtParam analysis revealed a critical weakness in the two most antigenic designs: RS09-Toxo and RSFN-Toxo carried instability indices of 60.23 and 50.51, well above the threshold of 40 that marks a protein as unstable in laboratory conditions. Only 50S-Toxo, with an instability index of 24.83, and Rpf-Toxo, at 35.41, passed the stability test, and both showed respectable aliphatic indices suggesting moderate thermal tolerance. All four constructs were hydrophilic, as indicated by negative GRAVY scores, and all exceeded the solubility threshold of the Protein-Sol predictor.

With Rpf-Toxo selected as the lead candidate, the team turned to structural biology. Homology modelling with I-TASSER produced a three-dimensional model that was then refined with the GalaxyRefine server, yielding improved geometry: 86.1 percent of residues fell into favoured Ramachandran regions, the MolProbity score dropped to 2.145 and the clash score fell to 9.6. The researchers candidly note that this figure remains below the 90 percent expected of high-resolution crystal structures, a limitation they attribute to the flexible spacers inherent in multi-epitope constructs. Molecular docking against the human TLR-4 receptor, the innate immune sensor central to Th1-polarised responses against the parasite, produced a stable predicted complex with a binding energy of −20.6 kcal/mol, although the authors caution that the corresponding dissociation constant of 7.3 × 10⁻¹⁶ M likely reflects overestimation by the PRODIGY algorithm rather than a physically realistic affinity.

Immune simulation with the C-ImmSim server, modelling a three-dose schedule at days 1, 84 and 186, painted an encouraging picture. Combined IgG and IgM titres peaked above 170,000 at day 65, memory B-cells reached 500 cells per cubic millimetre and settled at a stable set-point of around 200, and the helper T-cell population included roughly 1,800 memory cells that persisted through day 350. Interferon-gamma levels climbed to approximately 380,000 ng/mL, a signature consistent with the Th1-skewed response needed to eliminate intracellular tachyzoites, and natural killer cell activity peaked near 380 cells per cubic millimetre. Six conformational B-cell epitopes were also predicted on the folded surface of the construct, with ElliPro scores between 0.515 and 0.859, suggesting multiple antibody-accessible sites.

Safety and manufacturability checks completed the computational dossier. A BLASTp search against the human proteome found no significant homology, reducing the risk of autoimmune cross-reactivity, and the sequence was reverse-translated and codon-optimised for expression in E. coli K12, raising the codon adaptation index toward 1.00 and adjusting GC content to 65.63 percent. Restriction sites for Eco53kI and EcoRV, a Shine-Dalgarno ribosome-binding sequence and a stop codon were added to enable directional cloning into the pET-28a(+) vector for recombinant production.

The authors are emphatic that every result remains theoretical until confirmed experimentally. They list a series of limitations, including algorithmic bias in epitope prediction, the simplified nature of immune simulations, the absence of molecular dynamics analysis of the vaccine-receptor complex, and the lack of conservation analysis across genetically diverse T. gondii strains, which will be essential for assessing broad cross-protection. They also flag the inherent risks of multi-epitope design, such as epitope competition or immune masking, which flexible linkers can mitigate but not guarantee against. The proposed next steps include molecular dynamics simulations of the TLR-4 complex, advanced delivery platforms such as liposomes and nanoparticles, and challenge studies in murine models. If those experiments succeed, Rpf-Toxo could represent a meaningful step toward a long-sought human vaccine against one of medicine’s most persistent parasites.

Subject of Research: Immunoinformatic design of a multi-epitope vaccine candidate against Toxoplasma gondii based on dense granule antigens

Article Title: Rpf‐Toxo: A Preliminary Computationally Designed Dense Granule Antigen‐Based Multi‐Epitope Vaccine Against Toxoplasma gondii

Article References: Safari, M. H., Hosseini, S. A., Ghiabi, S., Ghiabi, S., Siamian, D., Majidiani, H., Irannejad, H., & Asghari, A. (2026). Rpf‐Toxo: A Preliminary Computationally Designed Dense Granule Antigen‐Based Multi‐Epitope Vaccine Against Toxoplasma gondii. Veterinary Medicine and Science, 12(5), Article e71189. https://doi.org/10.1002/vms3.71189

Image Credits: AI Generated

DOI: 10.1002/vms3.71189

Keywords: Toxoplasma gondii, toxoplasmosis, multi-epitope vaccine, dense granule antigens, reverse vaccinology, immunoinformatics, TLR-4, RpfE adjuvant, molecular docking, immune simulation, epitope prediction, vaccine design

The burden of sleeping sickness is considerable. The disease progresses in stages, beginning with fever, headaches and lymphadenopathy as parasites multiply in the blood and lymph, and advancing to neurological involvement once the parasites cross the blood-brain barrier, producing sleep disturbances, cognitive decline and, if untreated, death. The parasite’s most formidable weapon is antigenic variation: a dense coat of variant surface glycoproteins, or VSGs, that is continually reshuffled through gene conversion, allowing the parasite to stay one step ahead of the host antibody response. This immune evasion strategy, well documented in the literature, has long frustrated conventional vaccine development, which is one reason the researchers turned to conserved, functionally essential proteins as alternative targets.

Specifically, the team selected three T. brucei proteins as the basis for their construct: the variant surface glycoprotein itself, heat shock protein 70, and the vacuolar transporter chaperone complex. Heat shock protein 70 is a highly conserved molecular chaperone central to protein folding and stress responses, while the vacuolar transporter chaperone complex is involved in polyphosphate synthesis and acidocalcisome function, processes essential to parasite survival. By mining these proteins for immunogenic peptides, the researchers aimed to build a construct that combines surface exposure with conservation across strains, increasing the likelihood that an immune response raised against the vaccine would recognize the parasite before it can establish infection.

The design workflow followed the now-standard logic of reverse vaccinology. Protein sequences were retrieved from the UniProt database and aligned with tools such as Clustal Omega to assess conservation. Cytotoxic T lymphocyte epitopes were predicted for MHC class I presentation, and helper T lymphocyte epitopes for MHC class II, using the Immune Epitope Database analysis resource and related servers. B-cell epitopes were predicted with linear and discontinuous methods, including ElliPro for structure-based antibody epitope mapping. Each candidate epitope was then filtered for allergenicity with AllerTOP, for toxicity with dedicated peptide toxicity predictors, and for antigenicity with VaxiJen, ensuring that only immunogenic, non-allergenic and non-toxic peptides entered the final construct.

Population coverage analysis, which estimates how many people worldwide carry HLA alleles capable of presenting the chosen epitopes, returned a striking result: the vaccine candidate achieved 100 percent global population coverage. This metric matters because a vaccine that only fits a narrow slice of HLA diversity would leave large populations unprotected. Biophysical characterization of the final multi-epitope protein showed an aliphatic index of 71.23, indicating good thermal stability, and a GRAVY score of minus 0.719, indicating a hydrophilic, soluble protein likely to fold and express well. Solubility and instability assessments supported the view that the construct should behave as a stable, expressible protein in a cellular context.

Structural modeling came next. The tertiary structure of the vaccine construct was predicted and evaluated with a TM-score of 0.65 plus or minus 0.13 and a C-score of minus 0.50, values consistent with a reliable fold. The model was then refined, and validation metrics confirmed its quality: a Ramachandran score of 86.8 percent, meaning the vast majority of residues occupy favored or allowed backbone conformations, and a ProSA Z-score of minus 5.26, within the range expected for proteins of comparable size. Disulfide engineering was considered to further stabilize the fold, and secondary structure predictions from PSIPRED and SOPMA were used to cross-check the modeled architecture.

To test whether the vaccine could actually engage the innate immune sensors that trigger adaptive responses, the team docked the construct against Toll-like receptors 2 and 4, key pattern-recognition receptors on antigen-presenting cells. Molecular docking predicted strong binding, with energy scores of minus 1013.5 kJ/mol for TLR-2 and minus 1002.8 kJ/mol for TLR-4. These docked complexes were then subjected to molecular dynamics simulation, principal component analysis, dynamic cross-correlation matrix analysis and MM-GBSA binding free energy calculations, all of which supported the stability and favorable energetics of the receptor-vaccine interactions. In practical terms, the simulations suggest the vaccine construct should bind robustly to the very receptors that initiate the innate immune cascade.

Immune simulation provided the most direct readout of the construct’s potential immunogenicity. Using computational immune system modeling, the researchers predicted robust humoral and cell-mediated responses, including elevated B lymphocyte and T lymphocyte populations and rising titers of IgM and IgG antibodies over the simulated immunization course. Cytokine profiles indicated activation of both Th1-type and Th2-type pathways, the dual signature generally desired in a prophylactic vaccine. While such simulations are approximations of a vastly more complex biological reality, they serve as a critical screening step, allowing weak candidates to be discarded before any laboratory resource is spent.

Because the platform is mRNA, the team also optimized the nucleic acid sequence itself. Codon optimization for expression in Escherichia coli strain K12, conducted for in-silico cloning into the pET-28a(+) vector, yielded a codon adaptation index of 0.9688 and a GC content of 44.70 percent, both indicative of high expression potential. In-silico cloning confirmed that the construct could be inserted into the vector without disrupting restriction sites. Finally, minimum free energy analysis of the mRNA sequence was used to evaluate the structural integrity and stability of the transcript, an important consideration since mRNA secondary structure influences translation efficiency and vaccine performance.

The authors are careful to frame the work as a computational proof of concept rather than a finished vaccine. As they conclude, the in-silico designed candidate demonstrated strong structural stability, favorable receptor interactions and promising immunogenic potential against T. brucei, but experimental validation and in-vivo studies are required to verify its safety and efficacy. That caveat applies to the entire field of computational vaccinology: docking scores and immune simulations can prioritize candidates and dramatically shorten development timelines, but only animal studies and clinical trials can establish whether a designed construct protects against real infection. Still, for a disease with no licensed vaccine, in which drug therapy is costly, logistically difficult and increasingly challenged by resistance, a rationally designed mRNA candidate that clears every computational hurdle represents a meaningful step forward, and a template for applying the same pipeline to other neglected tropical parasites.

Subject of Research: Computational immunoinformatics design of an mRNA vaccine candidate against the sleeping sickness parasite Trypanosoma brucei

Article Title: Immunoinformatics Approach for the Designing of a Novel mRNA Vaccine Candidate Against Trypanosoma brucei

Article References: Nil, M. S., Khandker, S., Sadaf, S., Ahmed, N., Reja, S., Sayfullah, M., Ferdous, J., Saha, S., Alamri, A., Khan, M. S., Ahmed, S., Wajed, S., Mahdeen, A. A., & Siddiquee, N. H. (2026). Immunoinformatics Approach for the Designing of a Novel mRNA Vaccine Candidate Against Trypanosoma brucei. Acta Parasitologica, 71(5), Article 205. https://doi.org/10.1007/s11686-026-01388-w

Image Credits: AI Generated

DOI: 10.1007/s11686-026-01388-w

Keywords: Trypanosoma brucei, human African trypanosomiasis, mRNA vaccine, immunoinformatics, reverse vaccinology, epitope prediction, molecular docking, molecular dynamics simulation, TLR-2, TLR-4, codon optimization, neglected tropical disease