In silico study to predict promiscuous T cell and B cell-epitopes derived from the vaccine candidate antigens of Plasmodium vivax binding to MHC class-II alleles
Nazam Khan1, Mona N bin-Mwena1, Mashael W Alruways1, Noor Motair M Allehyani1, Maryam Owaid Alanzi1, Shahzad2, Amir Khan3, Rakesh Sehgal4, PK Tripathi4, Umar Farooq5
1 Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Shaqra University, Shaqra, Kingdom of Saudi Arabia
2 Universal Group of Institutions, Lalru, Punjab, India
3 Department of Basic Oral Medicine and Allied Dentistry, College of Dentistry, Taif University, Taif, Kingdom of Saudi Arabia
4 Postgraduate Institute of Medical Education and Research, Chandigarh, India
5 Department of Basic Oral Medicine and Allied Dentistry, College of Dentistry, Taif University, Taif, Kingdom of Saudi Arabia; Faculty of Biotechnology and Applied Sciences, Shoolini University, Himachal Pradesh, India
Correspondence Address:
Umar Farooq
Department of Basic Oral Medicine and Allied Dentistry, College of Dentistry, Taif University, Taif, KSA; Faculty of Biotechnology and Applied Sciences, Shoolini University, Solan, HP, India
Source of Support: None, Conflict of Interest: None
CheckDOI: 10.4103/0972-9062.335726
Malaria is one of the major causes of health and disability globally, even after tremendous efforts to eradicate it. Till date no highly effective vaccine is available for its control. The primary reason for the low efficacy of vaccines is extensive polymorphism in potential vaccine candidate antigen genes and HLA polymorphisms in the human population. This problem can be resolved by developing a vaccine using promiscuous peptides to combine the number of HLA alleles. This study predicted T and B cell epitopes (promiscuous peptides) by targeting PPPK-DHPS and DHFR-TS proteins of Plasmodium vivax, using different in silico tools. Selected peptides were characterized as promiscuous peptides on the basis of their immunogenicity, antigenicity and hydrophobicity. Furthermore, to confirm their immunogenicity, these peptides were utilized for molecular modelling and docking analysis. For determining the requisite affinity with distinct HLA Class-I, and HLA Class-II alleles, only five peptides for DHFR-TS and 3 peptides for PPPK-DHPS were chosen as promiscuous peptides. The D1 peptide has the maximum binding energy with HLA alleles, according to HLA-peptide complex modelling and binding interaction analyses. These findings could lead to the development of epitope-based vaccinations with improved safety and efficacy. These epitopes could be major vaccine targets in P. vivax as they possess a higher number of promiscuous peptides. Also, the B cell epitopes possess maximum affinity towards different alleles as analyzed by docking scores. However, further investigation is warranted in vitro and in vivo.
Keywords: Epitopes; Malaria; HLA; Immunodiagnostic; Promiscuous peptides; Homology modelling; Subunit vaccine; PPPK-DHPS; DHFR-TS
Malaria is a major tropical illness, with an anticipated global incidence of 1.5 billion cases, and 7.6 million deaths expected between 2000 to 2019[1], with roughly 94 percent of cases, and deaths occurring in Africa, where malaria mortality is higher than elsewhere in the world[1]. One of the five species of Plasmodium is Plasmodium vivax (Pv) that causes human infection. Vivax malaria is increasingly recognized as a leading cause of severe morbidity and mortality, with a considerable negative impact on the global health[2],[3]. P. vivax is frequently seen in low frequency but transmissible parasite quantities in the peripheral blood circulation of the host. It can also transform into undetectable dormant liver stages (hypnozoites), to cause recurrent infection and sickness (relapse). It is possibly due to this fact that antimalarial medications are ineffective in completely controlling the condition. Additionally, drug resistance to chloroquine and other anti-malarial medicines[4] is a significant challenge.
Vaccine development can be a successful strategy in controlling malaria due to the parasite’s strong survival in poor environmental circumstances and efficient transmission at lower parasite levels. The selection of a suitable antigen to produce an efficient immune response to control malaria is a primary emphasis in developing a vaccine against the disease[5].
Different proteins involved at specific stages of Plasmodium infection are the primary targets for developing an effective vaccine[6]. Several antigenic proteins have been identified, and their genes have been characterized from different stages of the life cycle of Plasmodium[4],[5],[6],[7]. The presence of extensive antigenic polymorphism in these vaccine candidate proteins is a major hurdle in developing a subunit vaccine. Genetic polymorphism in several vaccine candidate antigen genes i.e., msp-1[4]ama-1[7],[8] and dbp[9],[10], has already been reported from many parts of the world for P. vivax, including India[4],[7]. Extensive polymorphism exhibited by antigens of P. vivax and human leukocyte antigen (HLA) is one of the main reasons for the limited efficiency of the malaria vaccine[4].
Considering all these factors, promiscuous peptides that can bind with different HLA alleles could be utilized to develop peptide-based vaccines[11]. Rather than using a wet lab approach to identify a promiscuous peptide, we could use a computational approach to predict promiscuous peptides by screening and analyzing the target gene sequences; this would then demonstrate optimal binding to HLA alleles prevalent in the inhabitants to be vaccinated.
Checking binding efficiency of the predicted HLA alleles with promiscuous peptides by using in silico tools offers more practical and economic methods[12]. Moreover, computational methods also limit the use of animal models with wet lab experiments which also reduces the cost and further supports the designing of novel vaccine candidates. Furthermore, structure-based research (SAR) utilizing molecular docking is a key method for confirming the binding effectiveness of promiscuous peptides with HLA[13],[14]. Docking can be used to screen many promiscuous peptides to see how they interact with one another and if they can invigorate the immune system. This could be quite useful in determining the peptide immunogenicity.
Keeping these points in view, we planned to discover the T cell and B cell antigenic determinants targeting the potential vaccine candidates DHFR-TS and PPPK-DHPS of P. vivax. The identified promiscuous T cell epitopes were further docked with HLA alleles to know the interaction of these predicted epitopes with HLA alleles.
Material & MethodsRetrieval of the target sequences
The target genes of P. vivax were chosen based on the literature review. These target genes could potentially be used for the development of a vaccine against P. vivax. The amino acid sequence of these vaccine candidates i.e., DHFR-TS, and PPPK-DHPS were fetched from the NCBI (www.ncbi.nlm.nih.gov) using the accession no. DQ517900.1, and EU478871.1, and were compared with PlasmoDB database (www.plasmodb.org/plasmo), and UniProt (www.uniprot.org).
HLA allele selection
The major and frequently occurring HLA class-I alleles i.e., HLA-A*11:01, HLA-A*02:01, HLA-A*02:06, HLA-A*02:11, HLA-A*33:03, HLA-B*01:01, HLA-B*35:01 and HLA-B*40:06 and HLA class-II alleles i.e., DRB1*10:01,DRB1*140:4,DRB1*15:02,DRB1*03:01, DRB1*07:01, DRB1*11:01 and DRB1*15:01[15],[16],[17],[18] among Indian populations were selected for the study. These HLA DR alleles are common in the Indian communities, and our findings are being applied to P. vivax strains, which are common in India. They infect individuals particularly with one of the HLA DR alleles.
Prediction of promiscuous peptides
The T-cell specific epitopes were predicted using the Immune Epitope Database (IEDB) (www.iedb.org)[19] and the NetMHCPan web server (www.cbs.dtu.dk/services/NetMHCpan-3.0)[20]. Artificial neural networks were utilized to predict the binding of peptides to HLA class-I MHC molecules using the NetMHC3.0 server (www.cbs.dtu.dk/services/NetMHC-3.0)[21], NetMHCIIPan 2.4,[22] and Immune Epitope Database and Analysis Resource (IEDB) MHC class-II binding prediction tools NN-align,[23] (C) SMM-align,[24] (D) ARB (average relative binding)[25], and (E) NetMHCII2.2[26] were used to envisage peptides specific for HLA class-II alleles. For each tool, all settings and values were kept as default. Strong binders were defined as epitopes with IC50 values of less than 500 nm, weak binders were defined as peptides with IC50 values between 500–1000 nm, and non-binders were defined as epitopes with IC50 values > 500 nm[27]. Promiscuous peptides were selected using a 50% binding affinity criterion which defined a promiscuous peptide as one with an IC50 value less than 500 nm and 50% of the alleles. Additionally, A peptide with an IC50 value less than 500 nM with all the specified alleles was also classed as a 100% promiscuous peptide binder to find the best binding promiscuous peptides. The promiscuous peptides chosen were then assessed for affinity with various HLA alleles.
Identification of B-cell epitopes
B-cell epitopes were divided into two main categories:linear (continuous) epitopes and conformational (discontinuous) epitopes. However, different studies indicate that about 90% of the residues are conformational, while only 10% are linear. Two different web-based servers BCEpred and ABCpred, were used to predict linear B cell antigenic determinants. We may be able to find the most promising promiscuous peptides suitable for creating a peptides-based subunit vaccination by employing various servers[27]. BCEpred (www.imtech.res.in/raghava/ bcepred/) program analyzes physicochemical parameters such as hydrophilicity, flexibility/mobility, accessibility polarity, and exposed surface to forecast B cell epitopes[28]. The ABCpred ( www.imtech.res.in/raghava/abcpred/ )[29] server also predicted linear B cell epitope regions within an antigen sequence, using an artificial neural network.
Hydrophobicity attribute
Hydrophobicity analysis of selected promiscuous peptides was performed by using the peptide property calculator program[31]. Promiscuous peptides with <50% hydrophobicity are perfectly soluble in aqueous solution and the peptides with hydrophobicity >50% are moderately or insoluble in the aqueous solution. The promiscuous peptides with hydrophobicity >50% were selected further for protein homology modelling analysis.
Antigenicity and allergenicity of promiscuous peptides
The antigenicity of the predicted epitopes was determined by ANTIGENpro (www.scratch.proteomics. ics.uci.edu/) and VaxiJen (www.jenner.ac.uk/VaxiJen)[32]. ANTIGENpro is the first program that uses microarray reactivity data to predict full-length protein antigenicity for five pathogens. It’s an in silico predictor of protein antigenicity that is alignment-free and pathogen-agnostic. It basically gives an idea about the likelihood of a protein to be antigenic. Furthermore, the VaxiJen server automates protein classification based on physicochemical parameters rather than sequence alignment, it can work alone or in conjunction with other servers.
Modelling of promiscuous peptides (epitopes)
The 3D models of promiscuous peptides were generated using the in silico tool PEPFOLD, de novo peptide structure prediction[33],[34]. It’s a de novo peptide structure prediction server used to characterize peptides and protein fragments structurally[35],[36]. PEP-FOLD (www.bioserv.rpbs.univ-paris-diderot.fr/PEP-FOLD) allows for the treatment of 9-36 amino acid cyclic peptides which are linear and with disulfide bond. The best PEPFOLD models differ from the complete NMR structures by an average RMSD of 2.75 Å.
Modelling of HLA-alleles
The 3D models of promiscuous peptides (epitopes) were required to determine their prominent interaction with HLA alleles. All the promiscuous peptides were modeled using UCSF Chimera[36], a flexible tool for interactive visualization and study of molecular structures and related data, including density maps, supramolecular assemblies, sequence alignments, docking findings, trajectories, and conformational ensembles. Usually, the primary protein sequence is utilized for the Peptide Modelling. The 3D models of primary sequence of peptides were generated, and all parameters like Φ and Ψ angles were used with default values. Finally, all the peptide models were minimized by using an Amber field with 100 gradient steps.
Docking interaction of modelled structure of HLA-alleles with epitopes
AutoDock 1.4.6 was used for the 3D Model docking interaction of promiscuous peptides. AutoDock file was prepared by adding polar hydrogens and partial charges, with defining rotatable bonds. Adding polar hydrogens and merging nonpolar hydrogens were also used to make the modeled protein. In addition, the Kollman charge and atom type parameters were included. A grid map was created around the active site as well as the necessary surrounding surface. The Lamarckian genetic search algorithm was used, and the docking run was set at 250 runs. The maximum number of energy evaluations was 25,0000 / run and the maximum number of generations in the genetic algorithm was set at 27,000[36].
Ethical statement: Not applicable
ResultsT-cell epitope prediction
Many peptides were predicted from the two vaccine candidate antigens DHFR-TS and PPPK-DHPS, which possess different binding affinities with different classes of HLA-I and HLA-II alleles. Interestingly, a total of 609 peptides from DHFR-TS and 701 from PPPK-DHPS for HLA-I were predicted, whereas 603 and 695 peptides for HLA class-II allele were obtained from DHFR-TS and PPPK-DHPS, respectively. Promiscuous peptides were categorized on the basis of 100% binding affinity with different HLA class-II alleles. The promiscuous peptides of the DHFR-TS were 5/603 and for PPPK-DHPS 3 of 695. All these peptides have shown 100% binding affinity with all HLA class-II alleles [Table 1]. We have studied the hydrophobicity attribute of promiscuous peptides since it is an important physicochemical property that is used to characterize secondary structures in the proteins. [Table 1] shows that the peptides with hydrophobicity >50% are more soluble in water than those with hydrophobicity <50%. We have observed that D1, D2, and P1 are antigenic against both the bacterial as well as viral pathogens [Table 1]. However, no promiscuous peptides for class HLA-I alleles were found from both the genes DHFR-TS and PPPK-DHPS.
Table 1: Comparative analysis of antigenic properties of promiscuous peptides with bacteria and viruses predicted form DHFR-TS, PPPK-DHPS gene of Plasmodium vivaxLinear B-cell epitope prediction
(i) DHFR-TS, PPPK-DHPS by BCEpred
Antigenic propensity, flexibility, hydrophilicity, accessibility, exposed surface, turns, polarity factors were used to identified B-cell linear epitopes within DHFR-TS and PPPK-DHPS gene sequences. The best-predicted region for B-cell epitope was observed in between 131–135 base pair. This region was present prominently in almost all parameters analyzed to identify B-cell epitopes [Table 2] & [Table 3].
Table 2: Identified B-cell epitopes for DHFR-TS gene using BCEpred serverTable 3: Identified B-cell epitopes for PPPK-DHPS gene using BCEpred server(ii) DHFR-TS, PPPK-DHPS analysis ABCpred
The B-cell epitopes were further identified using AB-Cpred online prediction tool. The threshold was adjusted to 0.8, and the window length was adjusted to 16. The remaining parameters were unchanged. A total of thirty-one B-cell epitopes were found from DHFR-TS region, and the minimum and maximum threshold scores of epitopes ranging from 0.80 to 0.93 [Table 4]. Similarly, a total of twenty-seven B-Cell epitopes were identified from PPPK-DHPS with minimum and maximum threshold score of B cell epitopes ranging from 0.80 to 0.93 [Table 5].
Hydrophobicity of promiscuous peptides
The hydrophobicity analysis revealed that the promiscuous peptides D1, D2, D4, D5, and P1, P2 had confirmed hydrophobicity <50% [Table 1]. Hence, all these peptides D1, D2, D4, D5, and P1, P2 were selected for further structural verification.
Antigenicity and allergenicity
The antigenicity of the promiscuous peptides was predicted using two servers. The VaxiJen 2.0 server (Threshold value for bacterial model: minimum 0.1004 and maximum 0.7127, and threshold value for parasitic model: minimum -0.2915, and maximum 0.6167). The ANTIGENpro server with the minimum and maximum threshold values observed were 0.089440, and 0.190487, respectively. The antigenicity, and allergenicity analysis revealed that all the predicted promiscuous peptides were antigenic in nature with or without adjuvant.
Homology modelling of HLA-alleles
(i) Molecular docking
The molecular docking of the promiscuous peptides with HLA class-II alleles with AutoDock web-based program revealed that the highest affinities were obtained for DHFR-TS promiscuous peptides with docking energies ranging from -8.1 to -9.6 [Table 6]. The promiscuous peptide D1 shows the highest binding energy, i.e., -8.1 with DRB1*01:01 [Figure 2] which reveals that this peptide is best fit in the peptide-binding cleft of HLA, and based on which that peptide could be further utilized for the vaccine development. The docking analysis of other peptides P2, D2, and P1 have shown low binding energy -9.7, -9.6, and -8.6, respectively [Figure 1],[Figure 2],[Figure 3],[Figure 4],[Figure 5]; hence these peptides were not suitable for further study.
Figure 1: Molecular Docking of peptide D1 (LQPAQFIHILGNAHV) with HLA-DRBI*01:01 (3QXA).Figure 2: Molecular Docking of P1 (RLHFLVLNGVPRYRV) peptide with HLA-DRB1*01:01 (3QXA).Figure 3: Molecular Docking of DHFR-TS peptide D2 (QPAQFLIHILGNAHVY) with HLA-DRBI*15:01 (1BX2).Figure 4: Molecular docking of P2 (LHFLVLNGVPRYRVL) peptide with HLA-DRB1*15:01 (1BX2).Figure 5: Pictorial representation of docked structure of P2 peptide (LHFLVLNGVPRYRVL) with HLA-DRB1*01:01 A in surface. and mesh B. form. The peptide is binding in inner groove of HLA molecule. DiscussionMalaria is a serious public health problem in the tropical and sub-tropical regions of developing countries due to its high morbidity and fatality rates in people. Several efforts are currently being made for the development of a malaria vaccine. Several vaccines have been developed and have entered phase 1 and 2 clinical trials. During clinical trials, these vaccines have shown low efficacy, the reason of low efficacy may be due to the presence of antigenic polymorphism in vaccine candidate antigens and HLA DR alleles. Hence predicting promiscuous peptides derived from these vaccine candidate genes that could bind with different HLA alleles that triggers effective immune response against parasites is a challenge.
In silico tools, for identifying vaccine candidate antigens could be an advantage in searching for a novel and potential vaccine candidate that can elicit a better immune response. The vaccine development based on synthetic peptides is already proved to a be successful strategy[39],[40]. Previously, we used in silico tools to predict promiscuous peptides predicted from DBL-3γ, CIDR-1 region of Pfemp-1 antigen, and EBA-175 antigen of P. falciparum. All the predicted peptides have shown immunogenicity in vitro and successfully proliferated T-cells in PBMCs culture[14],[41]. The Flow cytometric analysis of sensitized cultured PBMCs supernatant revealed the presence of anti-inflammatory cytokines in high concentrations[14]. Following this, we also predicted immunogenic promiscuous peptides from the promoter region of the Pfs-25 antigen of P. falciparum by computational analysis[11]. This study aimed to identify B and T cell epitopes from two vaccine candidate antigen genes DHFR-TS and PPPK-DHPS of P. vivax. Considering some of the factors like antigenicity, hydrophobicity, and immunogenicity, only two peptides for each gene were found to be promiscuous in both DHFR-TS and PPPK-DHPS genes. Similarly, promiscuous T-cell epitopes were identified from MSP-1 of P. vivax[37]. Furthermore, these four peptides were selected for structure-based analysis by using docking analysis. The B-cell epitopes were also predicted for the same genes, and the best-predicted region for the B-cell epitope was 131–135, because these regions are present in all parameters. Yoelis et al.[38] also identified T-cell, and B-cell epitopes by using in silico tools. The present study suggests that the PPPK-DHPS and DHFR-TS targeting epitopes could be the potential vaccine targets to develop a subunit vaccine against P. vivax. The study also determined that these two antigens PPPK-DHPS and DHFR-TS possess many promiscuous peptides (T-cell epitopes and B-cell epitopes), and have the highest affinity towards different alleles as analyzed by docking scores. These findings supported the selection of promiscuous peptides as vaccine candidate antigens. However, further investigation using in vitro and in vivo models is required.
ConclusionThere are many vaccines that have been developed for malaria that entered phase I clinical trials, however, only few have reached the phase II clinical trial like Pvs25, and PvCSP based on the sexual stage, and pre-erythrocytic stages, respectively. Recently, another recombinant vaccine VMP001/AS01B based on CSP entered in Phase 1/2a clinical trials[42]. In the prevention of uncomplicated malaria, RTS,S/AS01E was found to be non-inferior to chemoprevention. In comparison to any intervention alone, the combination of these interventions resulted in a significantly lower frequency of fatal malaria[43]. The RTS,S/AS01 malaria vaccine was recently recommended by World Health Organization (WHO) for the prevention of P. falciparum malaria in children residing in WHO-defined moderate to high transmission areas. Only selective vaccines elicit an effective immune response against the challenge. The theory of a peptide-based vaccination for P. vivax, on the other hand, could be investigated. It is based on the prediction of antigenic epitopes that can effectively activate the immune system and protect the host from malarial infection. This computer-aided vaccine design approach is a time-saving and cost-effective technique. In this study, promiscuous peptides were identified from two vaccine candidate antigen genes PPPK-DHPS and DHFR-TS of P. vivax, by using a combined approach of sequence-based and molecular simulation-based. It is an efficient protocol for the rationale design of the vaccine against various infection-causing organisms. Conclusively, more in vitro analysis and in vivo could aid in the development of a peptide-based vaccine.
Conflict of interest: None
AcknowledgementsWe acknowledge the help of Shoolini University, Solan, HP, India for providing the facility and financial support to complete this work.
References
留言 (0)