Identification of immunodominant epitopes in allelic variants VK210 and VK247 of Plasmodium Vivax Circumsporozoite immunogen

• Published in: Infection, genetics and evolution genetics and evolution, 2021-12, Vol.96, p.105120-105120, Article 105120
• Main Authors: Naz, Shumaila, Ahmad, Sajjad, Abbasi, Sumra Wajid, Ismail, Saba, Waseem, Shahid, Tahir ul Qamar, Muhammad, Almatroudi, Ahmad, Ali, Zain
• Format: Article
• Language: English
• Subjects: Binding energy calculations; Circumsporozoite protein; Immunodominant Epitopes - genetics; Immunodominant Epitopes - immunology; Immunoinformatics; Molecular docking; Molecular dynamic simulation; Plasmodium vivax; Plasmodium vivax - physiology; Protozoan Proteins - genetics; Protozoan Proteins - immunology
• ISSN: 1567-1348; 1567-7257
• DOI: 10.1016/j.meegid.2021.105120

Plasmodium vivax-induced malaria is among the leading causes of morbidity and mortality in sub-tropical and tropical regions and infect 2.85 billion people globally. The continual rise and propagation of resistance against anti-malarial drugs is a prerequisite to develop a potent vaccine candidate for Plasmodium vivax (P. vivax). Circumsporozoite protein (CSP) is an important immunogen of malaria parasite that has the conserved CSP structure as an immune dominant B-cell epitope. In current study, we focused on designing multi-epitope vaccines (MEVs) using various immunoinformatics tools against Pakistani based allelic variants VK210 and VK247 of P. vivax CSP (PvCSP) gene. Antigenicity, allergic potential and physicochemical parameters of both PvCSP variants were assessed for the designed MEVs and they were within acceptable range suitable for post experimental investigations. The three-dimensional structures of both MEVs have been predicted ab initio, optimized, and validated by using different online servers. The both MEVs candidates were stable and free from aggregation-prone regions. The stability of both MEVs had been improved by a disulfide engineering approach. To estimate the binding energy and stability of the MEVs, molecular docking simulation and binding free energy calculations with TLR-4 immune receptor have been conducted. The docking score of PvCSP210 and PvCSP247 for TLR-4 was −6.34 kJ/mol and − 2.3 kJ/mol, respectively. For PvCSP210-TLR4 system, mean RMSD was 4.96 Å while PvCSP247-TLR4 system, average RMSD was 4.49 Å. The binding free energy of PvCSP210-TLR4 complex and PvCSP247-TLR4 complex was −50.49/−117.15 kcal/mol (MMGBSA/MMPSA) and −52.94/−96.26 kcal/mol (MMGBSA/MMPSA), respectively. The expression of both MEVs produced in Escherichia coli K12 expression system by in silico cloning was significant. Immune simulation revealed that the proposed MEVs induce strong humoral and cellular immunological responses, in addition to significant production of interleukins and cytokines. In conclusions, we believed that the MEVs proposed in current research, using combine approach of immunoinformatics, structural biology and biophysical approaches, could induce protective and effective immune responses against P. vivax and the experimental validation of our findings could contribute to the development of potential malaria vaccine. [Display omitted] •Immunoinformatics study of allelic variants of Plasmodium vivax CSP (PvCSP) gene.•Prediction of anti