key: cord-0071697-l0ex9quc
authors: Lari, Alireza; Lari, Niloofar; Biabangard, Atefeh
title: Immunoinformatics Approach to Design a Novel Subunit Vaccine Against Visceral Leishmaniasis
date: 2021-12-16
journal: Int J Pept Res Ther
DOI: 10.1007/s10989-021-10344-3
sha: d59bbfd872419a91f89abbadffaea6c4446f3a96
doc_id: 71697
cord_uid: l0ex9quc

Visceral leishmaniasis (VL) infection is mostly caused by Leishmania donovani and affects countries worldwide. Despite the need for a safe and effective vaccine against leishmaniasis due to the increased drug resistance, however, no vaccine has yet been licensed for clinical use. This study revolves around the immunoinformatics approach to design a multi-epitope vaccine against VL infection. In this case, the proteome of L. donovani has been investigated, and three host non-homologous and antigenic extracellular secretory proteins have been identified as potential vaccine candidates with low transmembrane helices (≤ 1). The multi-epitope subunit vaccine construct consists of T-cell (cytotoxic T-lymphocyte (CTL) and helper T-lymphocyte (HTL)) epitopes accompanied by appropriate adjuvant and linkers. A 372-amino acid vaccine construct has been established with specific characteristics, such as soluble, stable, antigenic, non-allergenic, non-toxic, and non-host homologous. Besides, the tertiary structure of the designed vaccine was modeled and validated. Also, the stability and affinity of the vaccine- TLR4 complex were confirmed by using molecular docking and molecular dynamics (MD) simulation. In addition, in silico immunization assay showed the efficiency of this candidate vaccine to stimulate an effective immune response. Furthermore, the refined vaccine was optimized and cloned in the pET28a (+) vector, and its successful expression was confirmed virtually. However, the experimental validation is required to verify the multi-epitope vaccine efficacy against VL infection.

Leishmaniasis is a vector-borne disease caused by infection with Leishmania parasites that predominantly affects the world's most vulnerable and poorest populations. Leishmania parasites are transmitted to humans by the bite of infected female sandflies (Alvar et al. 2006) . The clinical manifestations of leishmaniasis depend on the causative Leishmania species. There are mainly three forms of leishmaniasis, including cutaneous, mucocutaneous, and visceral. The latest estimate of the World Health Organization (WHO) in 2019 demonstrated that more than 90% of the new cases of visceral leishmaniasis (VL) were reported in ten countries, including Brazil, Ethiopia, Eritrea, India, Iraq, Kenya, Nepal, Somalia, South Sudan, and Sudan (https:// www. who. int/ news-room/ fact-sheets/ detail/ leish mania sis). The VL, caused by Leishmania donovani and Leishmania infantum, is the most severe form of the disease. Also, it can lead to death if not treated. The mortality rate in VL is higher than in other types of disease. This issue is related to the emergence of resistance in the Leishmania strains against most antibiotics (Moore et al. 2010) .

Various VL vaccines have been tested in animal models, such as the fucose mannose ligand (FML)-saponin (Santos et al. 2003) , Leishmania excreted/secreted antigens (Lemesre et al. 2005) , and liposome-based vaccines (Afrin et al. 1997; Sharma et al. 2006) . Although vaccines are the most effective way to prevent diseases, there is still no approved vaccine against human VL (Kedzierski et al. 2006 ). This situation is due to the pathogen antigenic variation and gene polymorphism (Ghorbani et al. 2018; Kumar et al. 2014) . One of the computational methods for detecting protective antigens by screening the proteome of a pathogen is reverse vaccination (RV), which is highly recommended for vaccine development (Rappuoli 2000) . Therefore, researchers' have focused on selecting of the appropriate pathogenic antigen in designing vaccines. Immunoinformatics approaches are cost-effective methods used for designing a peptide-based vaccine. This vaccine contains the antigen part of pathogens that can induce the immunogenicity and affect the vaccine efficacy (Bertholet et al. 2009 ). Thus, developing a novel multi-epitope subunit vaccine by predicting T-cell epitopes (Patronov et al. 2013) and employing computational tools are efficient ways to terminate the virulence of this pathogen.

Some proteins in the proteome of Leishmania have been identified as unknown/hypothetical proteins (HPs). These proteins were potentially utilized as drug targets or vaccine candidates (Sinha et al. 2017) . Since some HPs substantially affect the pathogenesis of leishmaniasis, some of these proteins are used in a recombinant vaccine against infection with Leishmania parasites (Duarte et al. 2017; Martins et al. 2016; Ribeiro et al. 2018) . Besides, the membrane and secretory proteins play a crucial role in stimulating a host-specific immune response (Khatoon et al. 2019) . The extracellular secretory proteins of Leishmania contain several antigenic components and induce an immune response (John et al. 2012; Khan et al. 2020) . Therefore, these proteins are very helpful for novel vaccine development.

Several recent studies have investigated the secretory or membrane proteins of Leishmania using computational tools (Khan et al. 2020; Khatoon et al. 2017 Khatoon et al. , 2019 Vakili et al. 2018) . In these studies, a multi-epitope vaccine was designed against leishmaniasis. The aim of this study was to pinpoint the potential T-cell epitopes from three secretory HPs of the genome of L. donovani that can be further joined through proper linkers and adjuvant to design a multi-epitope subunit vaccine. Then, the designed vaccine's physicochemical properties, antigenicity, and allergenicity were predicted by in silico approaches. The secondary and tertiary structures of the designed vaccine were also predicted and validated. Molecular docking and molecular dynamics simulation have been performed to confirm the binding interaction and stability of the docked complex. Finally, in silico cloning into the pET28a (+) vector has been performed to confirm the expression efficiency of the final vaccine construct. The proposed vaccine is also expected to be effective in stimulating humoral and cellular immunity against VL infection.

In this study, protein sequences of L. donovani were derived from the TriTrypDB (http:// tritr ypdb. org/ tritr ypdb/) to identify new HP candidates for vaccine design. The amino acid sequences of three HPs of L. donovani (LDBPK_120390, LDBPK_300860, and LDBPK_312610) have been chosen as secretory proteins. This procedure has been accomplished by predicting the score derived from the SignalP (Armenteros et al. 2019 ) server and evaluating the combined prediction of SecretomeP (Bendtsen et al. 2004 ) and TMHMM 2.0 (Krogh et al. 2001) . Besides, CELLO2GO (Yu et al. 2014) and PSORTb v3.0 (Yu et al. 2010) have been employed to identify the subcellular localization of hypothetical L. donovani proteins. Also, the antigenic property has been measured for these HP sequences using VaxiJen 2.0 (Doytchinova et al. 2007 ).

Cytotoxic T lymphocytes (CTLs) play a crucial role in host defense against foreign antigens. Thus, it is essential to predict CTL epitopes for the vaccine design. The mentioned HP sequences of L. donovani were submitted to the NetCTL 1.2 (http:// www. cbs. dtu. dk/ servi ces/ NetCTL/) server to predict CTLs. In this case, the threshold value was set at 0.75. This server allows prediction of CTL epitopes based on MHC-I binding peptides, proteasomal C-terminal cleavage, and Transporter Associated with Antigen Processing (TAP)associated efficiency with antigen processing (Larsen et al. 2007 ).

Since the helper T cells had significant influences on adaptive immunity, the IEDB database (http:// tools. iedb. org/ mhcii/) has been utilized to predict helper T lymphocyte (HTL) epitopes based on their receptor binding affinity (Vita et al. 2019) . In this case, epitopes with the lowest percentile rank represent a high affinity of MHC-II. Also, the 15-mer epitopes in mouse alleles were retrieved for all HP sequences in FASTA format.

In this study, the predicted epitopes have been combined to construct a multi-epitope subunit vaccine. This procedure has been performed using suitable linkers and adjuvants. Besides, the selected CTLs and HTLs epitopes were linked via AAY and GPGPG linkers to connect various components of a multi-epitope vaccine construct. Also, the EAAAK linker has been considered to add the Mycobacterium tuberculosis 50 S ribosomal L7/L12 (Uniprot ID: P9WHE3.1) to the N-terminal sequence as an adjuvant (Lee et al. 2014 ). Thus, the immunogenicity of the vaccine has been enhanced.

Antigenicity is defined as the ability of a chemical structure to produce an immune response. In general, it is recognized as an antigen. Besides, allergenicity is characterized as the ability of small antigens in creating immune-mediated and allergic reactions. In this case, these two properties were identified as significant factors affecting vaccine development. Also, ANTIGENpro (Magnan et al. 2010 ) and Aller-TOP 2.0 (Dimitrov et al. 2014 ) servers have been measured the antigenic and allergenic properties of the designed vaccine protein. In addition, the toxicity of the designed vaccine construct has been checked by the ToxinPred (Gupta et al. 2013) server.

ExPASy ProtParam tool (https:// web. expasy. org/ protp aram/) has been used to calculate a set of physicochemical parameters such as theoretical isoelectric point (pI), amino acid composition, in vitro and in vivo half-life, instability, aliphatic index, molecular weight (MW), and grand average of hydropathicity (GRAVY) for the subunit vaccine protein (Gasteiger et al. 2005) . Also, the Protein-Sol (https:// prote in-sol. manch ester. ac. uk/) server has been employed to determine the solubility of the designed vaccine protein. This server predicts the vaccine solubility based on properties like isoelectric pH and amino acid compositions. In this case, protein scored ≥ 0.45 is expected to be soluble compared to the average solubility of E. coli proteins (Hebditch et al. 2017) .

Furthermore, the NCBI BlastP has been considered to compare the whole multi-epitope vaccine sequence with Homo sapiens proteome for preventing autoimmunity. In this case, non-host homologous proteins were screened with default parameters (Mahram et al. 2015) .

PSIPRED 4.0 (Buchan and Jones 2019) online web server has been considered to predict the secondary structure of the multi-epitope subunit vaccine construct. Besides, the homology modeling tool in the RaptorX (Kallberg et al. 2012) server has been utilized to prepare the tertiary model of the subunit vaccine protein. Also, the Galaxy Refine (Heo et al. 2013 ) server has been employed to refine the 3D model of the multi-epitope vaccine construct and improve the structure quality. In this study, the quality of the refined vaccine construct has been assessed by the ProSA-web server (https:// prosa. servi ces. came. sbg. ac. at/ prosa. php). It predicts the overall quality of the protein model in the form of a Z-score (Wiederstein et al. 2007) . Also, the Ramachandran plot was generated using the PROCHECK (https:// saves. mbi. ucla. edu/) server to illustrate the Ramachandran plot with allowed and disallowed regions of amino acids.

The epitopes of the humoral immune response, particularly B-cell epitopes, are antigenic determinants and considerably affect vaccine designing. In the present study, these linear B-cell epitopes (20-mer peptides) of the final vaccine construct were predicted using the BCPREDS server at default parameters. This server utilizes a support vector machine (SVM) based on kernel methods for predicting linear B-cell epitopes (El-Manzalawy et al. 2008) . Also, the Ellipro server (https:// tools. iedb. org/ ellip ro) has been employed to identify discontinuous B-cell epitopes of the refined final vaccine construct. This server implements three algorithms based on their protrusion index (PI) values and allocates a PI value to the predicted epitopes (Ponomarenko et al. 2008 ). Besides, interferon-gamma (IFN-γ) plays a vital role in innate and adaptive immunity. In this case, the IFNepitope (https:// webs. iiitd. edu. in/ ragha va/ ifnep itope/ index. php) server has been considered for predicting the IFN-γ epitope of the final vaccine. This server predicts epitopes via the hybrid algorithms of motif and SVM (Dhanda et al. 2013 ).

The Computed Altas of Surface Topography of proteins (CASTp) server V3.0 has been employed to predict the protein binding site and hydrophobic interaction on the protein surface (Dundas et al. 2006) . This situation initiates an interaction between the final vaccine construct and the TLR4 immune receptor. Then, the PatchDock web server (http:// bioin fo3d. cs. tau. ac. il/ Patch Dock/) has been considered for the molecular docking of the final vaccine construct with TLR4 (PDB ID: 4G8A) receptor. This procedure has been carried out to demonstrate an immune response development and analyze the interaction between the ligand and receptor (Schneidman-Duhovny et al. 2005) . Afterward, the Fast Interaction Refinement in Molecular Docking (FireDock) server has been utilized to refine and re-score the rigid body molecular docking complexes (Andrusier et al. 2007 ). This situation has been led to achieving the best structure. This server provides the ten best solutions for the procedure. Besides, the refined models are configured based on several factors such as atomic contact energy, partial electrostatics, Van Der Waals interaction (vdW), and binding energy estimations. Also, the TLR4 receptor-vaccine interaction was visualized using the UCSF Chimera (Pettersen et al. 2004) and LigPlot+ (Laskowski et al. 2011) tools.

In addition, the stability of the interaction between TLR-4 receptor and the vaccine construct was determined by GROMACS 5.1.1 package (Van Der Spoel et al. 2005) . This process has been performed using the GROMOS96 43a1 force field parameter set to implement molecular dynamics simulations. Before simulation, residues 980 to 1200 of TLR4 receptor that were located considerably far from the vaccine-TLR4 complex (based on docking analysis) were manually removed to reduce the time of simulation and further improve MD analysis. In this case, the structure was initially solvated with transferable intermolecular potential with 3 points (TIP3P) water model. Sodium and chloride ions were also added to neutralize the system. Besides, the energy has been minimized to provide a proper molecular arrangement in space. It was conducted by eliminating high-energy configurations. Also, the simulation process has been run under canonical (NVT) and isothermal-isobaric (NPT) conditions at 300 K and 1 bar pressure with the duration of 50 ps. Then, the stabilized structure has been subjected to the MD simulation of 50 ns. Finally, the root mean square deviation (RMSD), the root mean square fluctuation (RMSF), the radius of gyration (Rg) (Lobanov et al. 2008) , and hydrogen bonds (Joshi et al. 2021) were calculated for all the systems. The MD complex analysis were determined by using graphs generated by Xmgrace software.

The immunological responses of the designed vaccine construct have been evaluated by in silico method C-ImmSim (https:// www. iac. cnr. it/ ~filip po/c-immsim/ index. html) online dynamic immune simulation server (Rapin et al. 2010) .

Reverse translation and codon optimization of the final vaccine construct were conducted by the Java Codon Adaptation Tool (JCat) server (http:// www. jcat. de/) to attain an improved nucleotide sequence and make it compatible with widely used prokaryotic expression system; E. coli K12 (Grote et al. 2005) . Besides, three options of avoiding rho-independent transcription terminators, prokaryote ribosome binding sites, and cleavage sites restriction enzymes were selected. The JCat output includes the codon adaptation index (CAI) value and the percentage of GC content. The desired CAI value is ranged from 0.8 to 1.0, and the optimal percentage of GC content ranges from 30 to 70%. Also, SnapGene 5.1.5 software cloned the optimized codon sequence into E. coli pET-28a (+) vector. Afterward, the restriction sites of XhoI and NdeI were incorporated in the N-and C-terminal sequences, respectively.

In this case, the amino acid sequences of three secretory proteins of L. donovani (LDBPK_120390, LDBPK_300860, and LDBPK_312610) were retrieved from the NCBI protein database to design an immunogenic subunit vaccine. The CELLO2GO (Yu et al. 2014 ) and PSORTb v3.0 (Yu et al. 2010 ) servers predicted that these secreted proteins were localized in extracellular location. These HPs were found with an antigenic probability of greater than 0.5. This value was in agreement with VaxiJen 2.0 (Doytchinova et al. 2007 ).

A total of five CTL epitopes with high-ranked binding affinity scores were predicted for the three nominated proteins (Table 1) . This process has been completed by the NetCTL 1.2 server (Larsen et al. 2007 ). Likewise, nine HTL epitopes based on their lowest percentile rank with MHC-II were predicted from the IEDB server (Vita et al. 2019) for mouse alleles (H2-1Ad, H2-1Ed, and H2-1Ab) to design the vaccine (Table 2) . 

The vaccine sequence has been constructed by combining five CTL epitopes and nine HTL epitopes with the AAY and GPGPG linkers, respectively. Also, the 50 S ribosomal L7/ L12 (P9WHE3) was added to the N-terminal of the construct using the EAAAK linker, which is ideal for achieving maximum immune response (Fig. 1 ).

In the designed vaccine construct, the immunological properties such as antigenicity, allergenicity, and toxicity were evaluated by the three servers of ANTIGENpro ( demonstrated that the designed vaccine protein was antigenic, non-allergenic, and non-toxic. Also, the scaled solubility value of the vaccine construct was 0.496. This value was higher than the average soluble E. coli protein derived from the experimental solubility dataset (Niwa et al. 2009 ). Besides, the NCBI BlastP has been utilized to search for similarities between the multi-epitope vaccine and the Homo sapiens proteome. The results showed no similarities between them. This situation reveals that the designed vaccine was non-homologous with the host proteome. Figure 2 depicts the secondary structure prediction of the subunit vaccine construct using the PSIPRED (Buchan and Jones 2019) server. It contains alpha-helix (38%), beta strands (17%), and coils (45%). Besides, the RaptorX (Kallberg et al. 2012 ) server has been utilized to construct a 3D structure of the subunit vaccine. The best template for Fig. 2 Schematic illustration of secondary structure of the subunit vaccine construct including alpha-helix, beta strands, and coils modeling the designed vaccine was the crystal structure of Thermotoga maritime (PDB ID: 1dd3A) with the p-value of 6.19e−05. Then, the GalaxyRefine (Heo et al. 2013) server refined the vaccine structure and generated five models. Among these generated models, model 1 (Fig. 3a) has been chosen due to its properties such as GDT-HA (0.9348), RMSD (0.458), MolProbity (1.734), Clash score (7.2), Poor rotamers (0.4), and Rama favored (95.1).

In this case, the final vaccine construct has been validated by the ProSA-web and Ramachandran plots after the model refinement. As shown in Fig. 3b , the ProSA-web server provided a Z-score of − 5.67 for the input vaccine protein. Figure 3c illustrates the Ramachandran plot of the vaccine construct, and only 1.7% of the residues are in the disallowed regions. These issues confirm the quality of the designed vaccine.

Two servers, BCPred (El-Manzalawy et al. 2008) and ElliPro (Ponomarenko et al. 2008) , have been employed to predict the linear and conformational B-cell binding epitopes among the primary input sequence of the final vaccine construct. Figure 4 depicts information about the eight epitopes. These epitopes had a BCPreds cut-off score of more than 0.9. A total of 66 c Ramachandran plot showing 98.3% of the residues are present in the favored and allowed regions, and 1.7% of residues are in disallowed region residues were predicted as the conformational B-cell epitopes with a probability score greater than 0.8 and residue numbers ranging from 303 to 370. These scores start with Tyrosine residue (score 0.85) and end with Tryptophan residue (score 0.98). The findings showed that the designed vaccine could enhance humoral immunity (Khatoon et al. 2017; Shi et al. 2015) .

The IFNepitope (Dhanda et al. 2013 ) server has been utilized to predict 364 epitopes (15-mer) based on SVM. This case included both negative and positive predictions. In this regard, two IFN-γ inducing epitopes were selected. These epitopes had a prediction score of one (Fig. 4) .

The CASTp server analyzed the active binding sites of the protein. It could serve as a possible site for interacting with the TLR4 receptor. Besides, this server predicted the potential binding site with 2073.8 Å 2 surface area and 2980.2 Å 3 molecular surface volume. Also, the PatchDock (Schneidman-DuhovnyInbar et al. 2005) server has been considered to perform the molecular docking of the vaccine construct against TLR4 at the predicted binding pocket. In this case, the top ten models were refined via their binding energy using the FireDock (Andrusier et al. 2007 ) server. Figure 5 displays the selected final model (model No. 2) for molecular docking based on the binding energy.

In order to validate structural stability of the subunit vaccine-TLR4 complex, molecular dynamics simulation was performed using a 50 ns MD simulation study. Obtained potential energy for the complex was − 7.72982e + 6 KJ/ mol. The complex attained equilibrium during the simulation trajectory after about 20 ns of simulation and acquired a value between 1.9 and 2.23 nm. The RMSD plot shows the complex stabilized in the last 30 ns of the simulation (Fig. 6a) . The RMSF plot of the TLR4-vaccine complex indicated that most of the peaks are belong to residues, which are located in the loop regions and the terminal residues of strand (Fig. 6b) . Because the loop structure is much more flexible, the loop region consisting of Val288-Ala372 residues showed larger RMSF values compared to α-helix structures. The overall predicted RMSF graph suggested, the ligand have a little fluctuated (less than 2 nm) throughout the simulation time. The compactness of vaccine and TLR4 complex evaluated and the average Rg was found to be 3.34 Å. As shown in Fig. 6c , we could assume that this complex did not induce noticeable fluctuation after 50 ns of simulation and remain compact after 25 ns that is in agreement with the RMSD values. The number of hydrogen bonds formed between the vaccine and the TLR4 through 50 ns simulation period is shown in Fig. 6d . The observed bonding graph, demonstrated that hydrogen bond interaction patterns increased slightly and then remained constant during the last 30 ns of simulation time. All evidences indicate that receptor tends to bind to vaccine and remain stable as previous RMSD and Rg represented.

The refined vaccine protein was subjected to in silico immune simulations to check the proficiency of the vaccine construct to the adaptive immune system. The simulation results showed that responses (secondary and tertiary) produced high levels of IgG and IgM antibodies (IgM + IgG and IgM), while the antigen level decreased (Fig. 7a) . Besides, Fig. 7b exhibits an effective antiviral response by the vaccine in the elevated levels of IFN-γ and interleukin-2. Also, Fig. 7c displays the estimation of the B-cell population in the elevated levels of B memory cells and B-cell isotype IgM. In addition, a high response was observed in the T cells (both T h and T c ) populations ( Fig. 7d and e) . The obtained results suggest that the proposed vaccine construct can effectively trigger the immune response and provide immunity against VL infections.

The Jcat (Grote et al. 2005 ) server has been utilized to adapt the codon for multi-epitope vaccine as per E. coli K12 strain. In this case, the Codon Adaptation Index (CAI) of the improved sequence was 0.968 with GC contents of 56.45%. Also, the length of the optimized codon sequence was 1,116 nucleotides. The results showed the possibility of good expression of the vaccine candidate in the E. coli host. Afterward, the SnapGene software generated the sequence of the recombinant plasmid to insert the optimized sequence into the pET28a (+) vector (Fig. 8) . This procedure has been performed for in silico cloning.

The VL infection is a serious public health problem, which annually affects millions of people across the globe. Leishmaniasis treatment is confined to several anti-leishmaniasis drugs, but these drugs showed increased drug resistance and toxicity. Therefore, there is no effective controlling approach for VL (Nagill et al. 2011) . Despite the efforts of scientists to develop a vaccine against leishmaniasis, there is still no effective vaccine against this disease. Several antigens (e.g., cysteine protease C and kinetoplastid membrane protein-11) have been evaluated as targets for vaccines against Leishmania (Basu et al. 2005; Das et al. 2014) . Also, recombinant vaccines (e.g., Leish111f and Leishmune) have been developed as vaccine candidates against leishmaniasis (Jain et al. 2015) . Nevertheless, no vaccine has been approved against leishmaniasis for use in humans (Kedzierski et al. 2006) .

Advances in bioinformatics tools have been led to the development of subunit vaccines containing specific immunogenic components of pathogens that stimulated an immune response (Adhikari et al. 2018; Pandey et al. 2018; Rahman et al. 2020; Shey et al. 2019 ). Since T-cells play a crucial role in controlling VL infection, most researchers have designed the vaccine via identifying L. donovani antigenic epitopes and their interaction with host MHC alleles and immune cells (Stager et al. 2012) . In this regard, most recent studies have focused on the L. donovani genome sequences to identify specific immunogenic components for designing novel vaccines using immunoinformatics methods (Atapour et al. 2021; Khan et al. 2020; Singh et al. 2020; Vakili et al. 2018) . Study on the pathogenic secretory proteins revealed that these proteins were antigenic in nature and had diverse functions during infection (Khatoon et al. 2017) . Identification of these proteins and understanding of their antigenic epitopes with host MHC alleles can be used to produce an effective vaccine against leishmaniasis (Kaye et al. 2011) . For example, the literature survey showed that amastigote-associated potential secretory proteins were an effective vaccine design against VL (Khan et al. 2020) . Also, the KMP-11 protein sequence (as a target antigen) and B-type flagellin (as adjuvant) have been used for designing vaccine construct (Atapour et al. 2021) .

This study emphasized the computational approaches for predicting the binding of peptides to MHC molecules. This procedure is performed to select the three hypothetical secretory proteins in L. donovani. These are used as effective multi-epitope vaccine candidates. In the present study, three L. donovani hypothetical secretory proteins have been scanned to predict the presence of a signal peptide, subcellular localization, and antigenicity. The potential epitopes of these proteins were selected against T-cell (Cytotoxic and Helper T-lymphocyte) epitopes using the available online tools. In this case, five CTL and nine HTL epitopes were chosen as part of the developed multi-epitope vaccine construct. This procedure was conducted by joining suitable linkers and an adjuvant together. None of these predicted epitopes have been used earlier that shows the novelty of our study.

The designed vaccine protein has been subjected to BLASTP search against human proteins to prevent the undesired pathogen cross-reactivity with the host. The analysis results demonstrated that the designed vaccine did not have any homology with the host proteome. Besides, the secondary structure of the vaccine construct including alpha-helices, beta-sheets, and coils was predicted. Also, the RaptorX and GalaxyRefine servers have been utilized to generate the 3D structure and model refinement, respectively. The Ramachandran plot showed that the refined tertiary structure of the vaccine improved to the desired level. Meanwhile, the final vaccine construct had physicochemical and antigenic properties with no probable allergenicity or toxicity.

Besides, molecular docking and MD simulation were performed to understand the binding affinity and stability of the TLR4-vaccine complex. The vaccine-TLR4 complex showed stability and compactness during MD simulation. Also, an immune simulation was carried out for the designed vaccine. This procedure has been accomplished to characterize its ability to produce immunological responses against the VL. The simulation results revealed that the final multi-epitope vaccine construct elicited unique immune responses. In addition, the in silico optimization of codon has been implemented for the final vaccine construct. The values of CAI and GC contents demonstrated that the vaccine construct could be expressed in the E. coli host. The in silico cloning of the optimized vaccine showed its better integration for the expression in a bacterial system. This study focuses on the design of a multi-epitope subunit vaccine using immunoinformatics tools, which suggests that the final vaccine construct may have a good protective effect against VL infection.

One of the major causes of mortality in tropical and subtropical regions across the world is VL. In this study, we utilized immunoinformatics tools to design a novel multi-epitope vaccine against VL. Based on the results, our multi-epitope subunit vaccine is immunogenic and non-allergic, which has the potential to induce cellular and humoral immune responses. Although in silico results showed the effectiveness of this vaccine, immunogenic potency of the proposed vaccine needs to be experimentally validated.

The author(s) declare no competing interests in this manuscript.

Immunoinformatics approach for epitope-based peptide vaccine design and active site prediction against polyprotein of emerging oropouche virus

Adjuvanticity and protective immunity elicited by Leishmania donovani antigens encapsulated in positively charged liposomes

Leishmaniasis and poverty

FireDock: fast interaction refinement in molecular docking

SignalP 5.0 improves signal peptide predictions using deep neural networks

Designing of a novel fusion protein vaccine candidate against human visceral leishmaniasis (vl) using immunoinformatics and structural approaches

Kinetoplastid membrane protein-11 DNA vaccination induces complete protection against both pentavalent antimonial-sensitive and -resistant strains of Leishmania donovani that correlates with inducible nitric oxide synthase activity and IL-4 generation: evidence for mixed Th1-and Th2-like responses in visceral leishmaniasis

Feature-based prediction of non-classical and leaderless protein secretion

Optimized subunit vaccine protects against experimental leishmaniasis

The PSIPRED protein analysis workbench: 20 years on

Combining cationic liposomal delivery with MPL-TDM for cysteine protease cocktail vaccination against Leishmania donovani: evidence for antigen synergy and protection

Designing of interferon-gamma inducing MHC class-II binders

AllerTOP v. 2-a server for in silico prediction of allergens

VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines

A vaccine composed of a hypothetical protein and the eukaryotic initiation factor 5a from Leishmania braziliensis cross-protection against Leishmania amazonensis infection

CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues

Predicting linear B-cell epitopes using string kernels

Protein identification and analysis tools on the ExPASy server

Leishmaniasis in humans: drug or vaccine therapy?

JCat: a novel tool to adapt codon usage of a target gene to its potential expression host

In silico approach for predicting toxicity of peptides and proteins

Protein-Sol: a web tool for predicting protein solubility from sequence

GalaxyRefine: Protein structure refinement driven by side-chain repacking

Vaccines for visceral leishmaniasis: a review

A reverse vaccinology approach for the identification of potential vaccine candidates from Leishmania spp

Molecular docking and molecular dynamics simulation approach to screen natural compounds for inhibition of Xanthomonas oryzae pv. Oryzae by targeting peptide deformylase

Template-based protein structure modeling using the RaptorX web server

Visceral leishmaniasis: immunology and prospects for a vaccine

Leishmania vaccines: progress and problems

Abd El Wahed A, Mondal D (2020) An immunoinformatic approach driven by experimental proteomics: in silico design of a subunit candidate vaccine targeting secretory proteins of Leishmania donovani amastigotes

Exploratory algorithm to devise multi-epitope subunit vaccine by investigating Leishmania donovani membrane proteins

Exploring Leishmania secretory proteins to design B and T cell multi-epitope subunit vaccine using immunoinformatics approach

Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes

Vaccines to prevent leishmaniasis

Large-scale validation of methods for cytotoxic T-lymphocyte epitope prediction

LigPlot+: multiple ligand-protein interaction diagrams for drug discovery

A potential protein adjuvant derived from Mycobacterium tuberculosis Rv0652 enhances dendritic cells-based tumor immunotherapy

Protection against experimental visceral leishmaniasis infection in dogs immunized with purified excreted secreted antigens of Leishmania infantum promastigotes

Radius of gyration as an indicator of protein structure compactness

High-throughput prediction of protein antigenicity using protein microarray data

Fast and accurate NCBI BLASTP: acceleration with multiphase FPGA-based prefiltering

A new Leishmania-specific hypothetical protein, LiHyT, used as a vaccine antigen against visceral leishmaniasis

Treatment of visceral leishmaniasis

Vaccine candidates for leishmaniasis: a review

Bimodal protein solubility distribution revealed by an aggregation analysis of the entire ensemble of Escherichia coli proteins

Novel immunoinformatics approaches to design multi-epitope subunit vaccine for malaria by investigating anopheles salivary protein

T-cell epitope vaccine design by immunoinformatics

UCSF Chimera-a visualization system for exploratory research and analysis

ElliPro: a new structure-based tool for the prediction of antibody epitopes

Vaccine design from the ensemble of surface glycoprotein epitopes of SARS-COV-2: an immunoinformatics approach. Vaccines (Basel)

Computational immunology meets bioinformatics: the use of prediction tools for molecular binding in the simulation of the immune system

Reverse vaccinology

Evaluation of a Leishmania hypothetical protein administered as DNA vaccine or recombinant protein against Leishmania infantum infection and its immunogenicity in humans

Immunotherapy against murine experimental visceral leishmaniasis with the FML-vaccine

PatchDock and SymmDock: servers for rigid and symmetric docking

Mohammad O (2006) Non PC liposome entrapped promastigote antigens elicit parasite specific CD8+ and CD4+ T-cell immune response and protect hamsters against visceral leishmaniasis

In-silico design of a multi-epitope vaccine candidate against onchocerciasis and related filarial diseases

Epitope-based vaccine target screening against highly pathogenic mers-cov: an in silico approach applied to emerging infectious diseases

Putative drug and vaccine target identification in Leishmania donovani membrane proteins using naive bayes probabilistic classifier

CD8(+) T cells in leishmania infections: friends or foes?

Immunoinformatics-aided design of a potential multi-epitope peptide vaccine against Leishmania infantum

GROMACS: fast, flexible, and free

The immune epitope database (IEDB): 2018 update

ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins

CELLO2GO: a web server for protein subcellular localization prediction with functional gene ontology annotation

PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes