Journal of Immunology Research
Volume 2018 (2018), Article ID 6718083, 22 pages
https://doi.org/10.1155/2018/6718083

Research Article

Immunoinformatics Approach for Epitope-Based Peptide Vaccine Design and Active Site Prediction against Polyprotein of Emerging Oropouche Virus

Utpal Kumar Adhikari,1,2 Mourad Tayebi,² and M. Mizanur Rahman¹

¹Department of Biotechnology and Genetic Engineering, Islamic University, Kushtia 7003, Bangladesh
²School of Medicine, Western Sydney University, Campbelltown, NSW, Australia

Correspondence should be addressed to M. Mizanur Rahman; mmrahmanbtg79@hotmail.com

Received 16 June 2018; Accepted 28 August 2018; Published 8 October 2018

Academic Editor: Nejat K. Egilmez

Copyright © 2018 Utpal Kumar Adhikari et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Oropouche virus (OROV) is an emerging pathogen which causes Oropouche fever and meningitis in humans. Several outbreaks of OROV in South America, especially in Brazil, have changed its status as an emerging disease, but no vaccine or specific drug target is available yet. Our approach was to identify the epitope-based vaccine candidates as well as the ligand-binding pockets through the use of immunoinformatics. In this report, we identified both T-cell and B-cell epitopes of the most antigenic OROV polyprotein with the potential to induce both humoral and cell-mediated immunity. Eighteen highly antigenic and immunogenic CD8⁺ T-cell epitopes were identified, including three 100% conserved epitopes (TSSWGCEEY, CSMCGLIHY, and LAIDTGCLY) as the potential vaccine candidates. The selected epitopes showed 95.77% coverage for the mixed Brazilian population. The docking simulation ensured the binding interaction with high affinity. A total of five highly conserved and nontoxic linear B-cell epitopes “NQKIDLSQL,” “HPLSTSQIGDRC,” “SHCNLEFTAITADKIMSL,” “PEKIPAKEGWLTFSKEHTSSW,” and “HHYKPTKNLPHVVPRYH” were selected as potential vaccine candidates. The predicted eight conformational B-cell epitopes represent the accessibility for the entered virus. In the posttherapeutic strategy, ten ligand-binding pockets were identified for effective inhibitor design against emerging OROV infection. Collectively, this research provides novel candidates for epitope-based peptide vaccine design against OROV.

1. Introduction

Oropouche virus (OROV) is the most common Orthobunyavirus of the Bunyaviridae family, an important causative agent of Oropouche fever in human widespread in South America, especially in Brazil. This fever is clinically characterized as an acute febrile urban arboviral disease [1, 2]. Oropouche fever causes a dengue-like illness with sudden onset of high fever (upwards of 105°F/41°C), headache, malaise, myalgia, arthralgia, nausea, vomiting, and dizziness, allied with viremia in the early five to six days. Usually, the patients recover within two to three weeks without any apparent sequel or death [3]. In some patients, OROV can cause viral meningitis with skin rash, neck stiffness, and photophobia [4]. The OROV can enter into the central nervous system (CNS) and can spread into the brain parenchyma tissue through some infective cycles. Lastly, it causes inflammation and severe manifestations of encephalitis [5]. In spite of its clinical relevance to public health, very little is known about the pathogenesis of OROV in humans and in other vertebrate animals. Proenca-Modena and coworkers demonstrated that the MAVS (mitochondrial antiviral-signaling protein), IRF-3 (interferon regulatory transcription factor) and IRF-7, and MAVS-dependent type I IFN signaling pathway (IFNAR) have dominant roles in restricting OROV infection, and this signaling in nonmyeloid cells greatly contributes to the host defense against orthobunya viruses [6].

OROV was first isolated from the blood of a febrile forest worker in Trinidad & Tobago in 1955. The second strain of OROV was isolated and collected from the pool of 177 Coquillettidia venezuelensis mosquitoes in Trinidad in 1960 [1, 7, 8]. Several outbreaks have been recorded in the South Americas, from both urban and rural people of Brazil, Panama, Peru, Ecuador, and Trinidad & Tobago [2]. The virus is deliberated as a public health threat in tropical and subtropical areas of Central and South America, with thirty outbreaks and over half a million (500,000) infected people as of 2005. In 2009, OROV reemerged in the municipalities of Altamira and Santa Barbara, Pará State, and subsequently in Mazagão in the state of Amapá, Brazil [1].

The virus is a single-stranded negative-sense and three-segmented RNA virus. These three segments are large (L), medium (M), and small (S), which encodes the viral RNA-dependent RNA polymerase, the surface glycoproteins (Gn and Gc), and the nucleocapsid (N) protein, respectively [8, 9]. Initially, the partial genome was sequenced for the Brazilian prototype OROV strain BeAn19991, and the complete and corrected genome sequences of the three segments, especially the M segment, were published in 2015 [10, 11]. The M segment comprises a single ORF which encodes a large polyprotein that is cleaved during or after the translation and produces a nonstructural protein (NSm), and two structural surface glycoproteins (Gn and Gc) [2]. These glycoproteins Gn and Gc help in the interaction between the virus particle and the cell receptor. According to Tilston-Lunel et al., the nonstructural protein (NSm) of the M segment is unessential for virus replication in the mosquito and mammalian cell lines [10]. Based on the S segment genetic information, the OROV can consist of four major genotypes (I, II, III, and IV), and Brazil is the main geographical region of these four genotypes [2]. The transmission of OROV is maintained in nature by two distinct cycles such as sylvatic and urban [12]. Bradypus tridactylus (sloth), some wild birds, and nonhuman primates (NHPs) control the sylvatic cycle, whereas the urban cycle is maintained only by humans [3, 8]. In humans, the OROV is transmitted by the midge (Culicoides paraensis) and the mosquito (Culex p. quinquefasciatus), but no evidence for direct transmission from human to human has ever been recorded [12, 13]. Diagnosis of OROV infection is performed by using various types of classic methods such as virus isolation and cell culture and molecular techniques such as serologic assays (CF, HI, and NT tests and in-house enzyme-linked immunosorbent assay) and real-time polymerase chain reaction (RT-PCR) in acute samples for genome detection [14–16].

Before 2016, there were no reports of the occurrence of Oropouche fever in humans in Goiânia city, the State of Goiás, and central Brazil. But the OROV antibody has been found in NHPs in Goiânia city [17]. So, immunity development in NHPs against this OROV should be taken into consideration as a vital warning sign of ongoing virus circulation.

Since effective and specific commercial rapid diagnostic tests are unavailable, there are great concerns for its possible future outbreaks similar to Zika in South American countries, especially in Brazil. Furthermore, there is no specific approved drug or vaccine against OROV available in the market. Although the role of the MAVS-dependent type I IFN signaling pathway has been shown to control OROV replication [6], it is not enough to completely eradicate OROV from the system. In this situation, adaptive immunity development against OROV could play an effective role in the complete elimination of the virus. Adaptive immunity works based on the recognition of specific epitopes by T-cells after reinfection. Consequently, humoral immune response is activated in humans.

In the present study, we have analyzed the complete proteome of the OROV for the identification of the most antigenic protein and its highest immunogenic and antigenic T-cell epitopes along with B-cell epitopes for the development of the epitope-based peptide vaccine using immunoinformatic approaches. We have also predicted the 3D structure of the most antigenic polyprotein as well as inhibitor-binding sites for the docking simulation study using various bioinformatics tools. Ultimately, this study is aimed at assisting the future laboratory efforts in developing effective vaccination for the prevention of OROV infection.

2. Materials and Methods

The flow chart representing the overall procedures of epitope-based vaccine design and ligand-binding pocket prediction for OROV polyprotein and is illustrated in Figure 1.

Figure 1: Schematic representation of epitope-based vaccine design and ligand-binding pocket prediction against the OROV polyprotein.

2.1. Retrieval of the Target Sequence

The FASTA formatted amino acid sequences of OROV were retrieved from the Virus Pathogen Database and Analysis Resource (ViPR) (http://www.viprbrc.org/), an integrated, powerful resource for several virus families and their respective species [18]. This database is integrated with the National Center for Biotechnology Information (NCBI) as well as UniProtKB database and helps to download the FASTA formatted sequence from both GenBank and UniProtKB. In this study, the nonstructural proteins (NSPs) were excluded from the total proteome of OROV. So, a total available 21 polymerase, 6 glycoprotein, 39 nucleocapsids, 20 polyproteins, and 46 nucleoprotein sequences were primarily selected for antigenicity prediction.

2.2. Highest Antigenic Protein Identification

To uncover the highest antigenic protein, the FASTA formatted amino acid sequences of total structural proteins were submitted to the VaxiJen v2.0 server (http://www.ddg-pharmfac.net/vaxijen/VaxiJen/VaxiJen.html) for antigenicity prediction [19]. This approach led to the selction of the highly antigenic protein for further analysis.

2.3. T-Cell Epitope Identification

The T-cell epitopes are typically peptide fragments which are immunodominant and can elicit specific immune responses, important for epitope-based peptide vaccine design. Due to the importance of T-cell epitopes, the NetCTL 1.2 server was used for the prediction of cytotoxic T-lymphocyte (CTL) epitopes derived from the protein sequence. The NetCTL 1.2 server (http://www.cbs.dtu.dk/services/NetCTL/) can predict human CTL epitopes in any given protein. The NetCTL 1.2 server can predict CD8⁺ T-cell epitopes for 12 supertypes such as A1, A2, A3, A24, A26, B7, B8, B27, B39, B44, B58, and B62 taking into account the proteasomal C terminal cleavage, MHC class I binding, and TAP transport efficiency. MHC class I binding, proteasomal C terminal cleavage, and TAP transport efficiency are predicted using artificial neural networks and the weight matrix, respectively [20, 21]. We identified epitopes for all the 12 supertypes. In this study, the threshold value for epitope identification was set at 1.25 which has a sensitivity and specificity of 0.54 and 0.993, respectively. The weight on C terminal cleavage and TAP transport efficiency were used as default parameters. Initially, 128 epitopes were selected based on the highest combined score, but the final selection completed after the antigenicity and immunogenicity prediction of the epitopes through the VaxiJen v2.0 server and IEDB server, respectively.

For the identification of both frequently and nonfrequently occurring MHC-I-binding alleles, the epitopes were analyzed by the stabilized matrix base method (SMM) in the IEDB analysis tool (http://tools.iedb.org/mhci/) [22, 23]. The amino acid length of peptide 9.0 and the IC50 value less than 200 were selected as parameters for the identification of MHC-I-binding alleles. Here, peptides with IC50 are considered as high affinity, <500 nM intermediate affinity, and <5000 nM low affinity. So, the lower IC50 value indicates higher affinity. IEDB is a resourceful server, and this server was also used for the prediction of processing score, proteasomal cleavage, TAP score, and the MHC-I-binding score of the selected epitopes and their respective alleles using the SMM [24, 25]. The SMM-align method was used for the prediction of MHC-II-binding alleles of the best candidate epitope [26].

2.4. Epitope Conservancy and Immunogenicity Prediction

For epitope conservancy and immunogenicity prediction, the conservancy (http://tools.iedb.org/conservancy/) and immunogenicity (http://tools.iedb.org/immunogenicity/) prediction tools of IEDB were used [26, 27]. The conservancy indicates the specific portion of a protein sequence that restrains the epitope and shows availability with a specific level of identity. More immunogenic peptides are superior to the less immunogenic peptides and considered as effective T-cell epitopes. So, the epitope with high immunogenicity was selected for further analysis.

2.5. Prediction of Population Coverage

We used the IEDB population coverage tool (http://tools.iedb.org/population/) for the analysis of the population coverage of the predicted epitopes and their respective MHC HLA-binding alleles [28].

2.6. Allergenicity and Toxicity Extrapolation

Allergenicity was anticipated through AllerTOP v. 2.0 (http://www.ddg-pharmfac.net/AllerTOP/) [28] and AllergenFP 1.0 (http://www.ddg-pharmfac.net/AllergenFP/) [29]. AllergenFP 1.0 has been developed based on a novel alignment-free descriptor-based fingerprint approach. AllergenFP 1.0 can identify both allergens and nonallergens with 87.9% accuracy [29]. On the other hand, AllerTOP v. 2.0 is a robust and strong complimentary approach which has been developed based on the method, -nearest neighbours (kNN), to classify allergens and nonallergens. The comparison among different allergen prediction servers designates that AllerTOP v. 2.0 is the best-performing method with 88.7% accuracy which is followed by AllergenFP 1.0 (87.9%) [28]. The toxicity of the peptides was predicted through the ToxinPred (http://www.imtech.res.in/raghava/toxinpred/) web server [30]. This method was developed based on the machine learning technique and quantitative matrix using different properties of peptides.

2.7. Three-Dimensional Structure Design of the Best Epitope

After analyzing the conservancy of epitopes, combined score, antigenicity, immunogenicity, nontoxicity, and the availability of alleles, we identified the candidate epitopes ¹⁰⁴³TSSWGCEEY¹⁰⁵¹ and ⁴LAIDTGCLY¹² and submitted to the PEP-FOLD peptide structure prediction server (http://bioserv.rpbs.univ-paris-diderot.fr/services/PEP-FOLD/) at the Ressource Parisienne en Bioinformatique Structurale Mobyle Portal [31, 32]. We used the best model for analysis of the interaction with MHC-I- and MHC-II-binding alleles.

2.8. Docking Simulation Study

To know the interaction between binding alleles and predicted epitope, the molecular docking simulation study was performed using AutoDockTools [33] and AutoDock Vina software [34]. For this purpose, the crystal structure of the HLA-B35:01 protein molecule named 3LKR was retrieved from the RCSB Protein Data Bank (PDB) as the PDB format [35]. However, the predicted crystal structure was in a complex form with protein and an epitope. So, the Discovery Studio (version 16.1.0.15350) was used to simplify the complex structure. After the separation, the protein molecule and the peptide (epitope) were then converted to PDBQT format through AutoDockTools. For the identification of binding energy at the binding groove of HLA-B35:01 with an epitope, the space box center was set at 4.496, −12.428, and −35.544 Å in the , , and axes, respectively. The size was set at 30, 24, and 18 Å in the , , and dimensions, respectively, and these analyses were done at a 0.964 Å spacing parameter. Finally, the docking simulation was performed by using AutoDock Vina.

Another docking study was performed as a control for the critical evaluation and scientific acceptance of our docking study. The docking method and parameters were set just like the previously mentioned parameters so that we can easily compare between sample and control. In the control section, the docking simulation was completed between influenza NP418 epitope from the 2009 strain and the HLA-B35:01 protein molecule. The docking parameters for HLA-DR and the selected epitope were set at 20.066, 31.590, and 29.302 Å in the , , and axes, respectively, and the size was set at 24, 36, and 18 Å in the , , and dimensions, respectively. This docking was completed at a 0.769 Å spacing parameter.

2.9. Prediction of B-Cell Epitopes

B-cell epitope identification is the vital step for epitope-based peptide vaccine design. So, the B-cell epitopes were identified from the highest antigenic protein through the online BepiPred-2.0 (http://www.cbs.dtu.dk/services/BepiPred/) [36] and LBtope methods (http://crdd.osdd.net/raghava//lbtope/) [37]. BepiPred-2.0 has been developed based on the random forest algorithm which was trained on epitopes annotated from antigen-antibody protein structures. BepiPred-2.0 can be considered as the best and up-to-date B-cell epitope prediction method because it shows significant results both on epitope data derived from a large collection of linear epitopes taken from the IEDB database and on crystallography-derived structural epitope data [36]. LBtope is another powerful tool for linear B-cell epitope prediction which has been developed based on the experimentally validated non B-cell epitopes obtained from the IEDB database [37]. The Lbtope_confirm dataset, the amino acid length 15-mer, and the scoring threshold of 60% were selected as the parameters for the prediction of linear B-cell epitopes in the LBtope server.

The conformational or discontinuous B-cell epitopes were predicted by using the IEDB conformational B-cell prediction tool ElliPro (http://tools.iedb.org/ellipro/) [38]. The conformational B-cell epitope is the sequence of amino acid or subunits which comprise an antigen and can come in direct interaction with a receptor of the immune system. However, ElliPro is the most comprehensive method because it can predict both the linear and conformational epitopes based on a protein 3D structure and showed the score as a PI (protrusion index) value [38]. The parameters for conformational epitope prediction were set at 0.8 and 7 Å for a minimum score and maximum distance (Angstrom), respectively.

2.10. Ligand-Binding Pocket Prediction

For the assessment of various physiological and physicochemical properties of the selected most antigenic OROV polyprotein, different Bioinformatics tools such as ExPASy ProtParam tools (https://web.expasy.org/protparam/) and the PSIPRED (http://bioinf.cs.ucl.ac.uk/psipred/) server were used. ExPASy ProtParam [39] was used for the prediction of amino acid composition, molecular weight, extinction coefficient, isoelectric point (pI), instability index, aliphatic index, and grand average hydropathicity (GRAVY) value. The secondary structural properties such as alpha helix, extended strands, and random coils were predicted through the PSIPRED server [40]. DISOPRED3, an integrated tool of the PSIPRED server, was used for the enumeration of disordered regions or disease-causing regions [41]. The I-TASSER server (https://zhanglab.ccmb.med.umich.edu/I-TASSER/), a unified platform for automated protein structure prediction [42, 43], was used for the design of the three-dimensional (3D) structure of the OROV polyprotein. The predicted best model was refined two times by using the ModRefiner Server (https://zhanglab.ccmb.med.umich.edu/ModRefiner/) [44] and further refined by the Swiss-PdbViewer version 4.1.0 [45] with a harmonic constraint of 20 kJ/mol/Å². This energy score was applied to all protein atoms, using the steepest descent and conjugate gradient method to reduce the bad contacts between protein atoms and structural water molecules. Computations were done in vacuo with the GROMOS96 43B1 parameter set, without reaction field implementation of Swiss-PdbViewer. To validate the protein model, we used different evaluation tools such as ProSA [46], QMEAN [47], and the ProQ server [48]. The stereochemical property of the predicted protein model was checked using PROCHECK through the PDBsum server [49, 50]. Finally, the refined and evaluated protein structure was submitted to the PMDB database (https://bioinformatics.cineca.it/PMDB/) [51] for a unique accession identity number. The hydrophobicity properties of the protein were calculated through the Discovery Studio (version 16.1.0.15350).

The active site or ligand-binding site analyses make a clear perception of the molecular docking simulation study. The binding pockets of the OROV polyprotein structure were explored by the SiteHound-web (http://scbx.mssm.edu/sitehound/sitehound-web/Input.html) [52] using default parameters. This server identifies ligand-binding sites by computing the interactions between a chemical probe and a protein structure.

3. Results

3.1. Selection of the Highest Antigenic Protein

The query for Oropouche virus structural and nonstructural proteins resulted in a total of 208 hits in the ViPR database. All the structural proteins were evaluated by the VaxiJen v2.0 server, which predicted an overall score for each protein sequence and denoting the protein sequence as antigen or nonantigen. The protein sequence with the GenBank id ALB07207.1 exhibited the highest antigenic score 0.5124 among all the structural proteins. The protein itself is the OROV polyprotein containing 1420 amino acid residues, and this polyprotein was used for further analyses.

3.2. T-Cell Epitope Identification

The NetCTL 1.2 web server was used for the prediction of CD8⁺ T-cell epitopes, and we considered the combined score for the epitope selection. In this study, we unraveled 128 epitopes in total, which achieved the selected threshold value 1.25 (Table S1). The predicted T-cell epitopes were further evaluated by the VaxiJen server and the IEDB immunogenicity prediction tool. It has been repeated twice due to the loss of antigenicity property perceived in some T-cell epitopes. This repetition results facilitated to get a good rate of T-cell epitope prediction. Herein, among 128 primarily selected T-cell epitopes, 59 epitopes were under the threshold value 0.00 on the VaxiJen server, whereas 56 epitopes contained the immunogenicity value less than 0.00 (Table S1). So, the results revealed that all epitopes with the highest combined score could not be the best epitopes. Among 128 epitopes, we selected 37 epitopes which have positive antigenicity and immunogenicity scores, so that these highly antigenic epitopes can interact to the MHC alleles with high affinity and can create an effective immune response (Table S2). Finally, we identified 18 epitopes which contained more than >70% conservancy score, because the conservancy score is an important issue during the development of an effective vaccine. Epitope conservancy among the different strains of the selected protein can provide immunization effectively. So, the higher conservancy of an epitope ensures the better target for improved vaccine design. The selected epitopes, conservancy score, and their position are shown in Table 1. In this study, 18 epitope conservancy analysis results revealed that the epitopes TSSWGCEEY, CSMCGLIHY, and LAIDTGCLY are 100.0% conserved, whereas the rest of the other 15 epitopes are conserved from 71.43% (RQHRGCIRF, YNVAWRTYL, and YLKNHNIDL) to 95.24% (CEEYGCLAI). However, among 18 T-cell epitopes, seven epitopes, including LIHYRPGLK, YLKNHNIDL, CEEYGCLAI, VPRYHSIDV, IPAISGLGV, AISGLGVGY, and RQHRGCIRF, have the VaxiJen score of more than 1 (>1) which indicates the high antigenicity nature of these epitopes (Table 1). The immunogenicity score of the epitope can be a good criterion for the selection of the best epitope. The higher score designates a greater probability of eliciting an effective immune response. Herein, the epitopes RQHRGCIRF (0.18739), LEPIIGDKL (0.18454), and NLPHVVPRY (0.12889) are less immunogenic than the epitopes IFAAIIFAL (0.42861), SLIEICITL (0.3578), YNVAWRTYL (0.34779), IIFALIITK (0.32685), and TSSWGCEEY (0.29444), among 18 epitopes (Table S2).

Table 1: Most potential 18 T-cell epitopes with interacting MHC-I alleles, epitope conservancy score, allergenicity, and toxicity.

Epitopes
Position
Antigenicity
MHC-I alleles
Conservancy
Allergenicity
Toxicity

TSSWGCEEY
1043–1051
0.9867
HLA-B27:20, HLA-C03:03, HLA-A32:07, HLA-A68:23, HLA-C12:03, HLA-B40:13, HLA-A32:15, HLA-C07:01, HLA-B15:02, HLA-B15:17, HLA-B35:01
100.00%
No
Yes

CSMCGLIHY
48–56
0.5400
HLA-B27:20, HLA-A68:23, HLA-C12:03, HLA-A32:07, HLA-B15:17, HLA-B40:13, HLA-B15:02, HLA-B15:01, HLA-C14:02, HLA-A32:15, HLA-B15:03, HLA-A30:02
100.00%
No
Yes

LAIDTGCLY
4–12
0.9551
HLA-C03:03, HLA-C12:03, HLA-B15:17, HLA-A68:23, HLA-B35:01, HLA-A32:07, HLA-A29:02, HLA-B15:01, HLA-B15:03, HLA-B15:02, HLA-B40:13
100.00%
Yes
Yes

SLIEICITL
34–42
0.9867
HLA-A02:11, HLA-A32:07, HLA-A02:12, HLA-A02:16, HLA-A02:01, HLA-A68:23, HLA-A02:02, HLA-A02:50, HLA-B27:20, HLA-C03:03, HLA-A02:03, HLA-A02:06, HLA-C12:03, HLA-A02:17, HLA-B15:02, HLA-C07:02, HLA-C14:02, HLA-A02:19, HLA-A32:15
76.19%
Yes
Yes

IIFALIITK
12–20
0.7267
HLA-C12:03, HLA-A11:01, HLA-A32:07, HLA-B40:13, HLA-A03:01, HLA-B27:20, HLA-A68:23, HLA-A68:01, HLA-A32:15, HLA-C03:03, HLA-A30:01
80.95%
Yes
Yes

LIHYRPGLK
53–61
1.7718
HLA-A32:07, HLA-B27:20, HLA-C12:03, HLA-C03:03, HLA-A68:23, HLA-A30:01, HLA-A03:01, HLA-B40:13, HLA-A32:15
80.95%
Yes
Yes

IFAAIIFAL
8–16
0.6140
HLA-C03:03, HLA-A02:50, HLA-B27:20, HLA-A02:02, HLA-A68:23, HLA-B40:13, HLA-A32:07, HLA-A02:11, HLA-A02:06, HLA-B15:02, HLA-C14:02, HLA-B42:01, HLA-C12:03, HLA-C07:02, HLA-A32:15, HLA-A68:02, HLA-A02:17
80.95%
Yes
Yes

NLPHVVPRY
25–33
0.6878
HLA-B27:20, HLA-A32:07, HLA-A68:23, HLA-C12:03, HLA-A02:17, HLA-B40:13, HLA-C14:02, HLA-C07:02, HLA-A32:15, HLA-B15:02
80.95%
Yes
Yes

NTIPAISGL
54–62
0.7847
HLA-A68:23, HLA-C03:03, HLA-A68:02, HLA-A02:50, HLA-A32:07, HLA-B15:02, HLA-B27:20, HLA-A26:02, HLA-C12:03, HLA-A02:17, HLA-A02:06, HLA-A69:01, HLA-C07:02, HLA-A32:15, HLA-B15:17
80.95%
Yes
Yes

VPRYHSIDV
30–38
1.5628
HLA-A68:23, HLA-B27:20, HLA-C12:03, HLA-A32:07, HLA-B07:02, HLA-B42:01
80.95%
No
Yes

IPAISGLGV
56–64
1.2054
HLA-A02:50, HLA-B27:20, HLA-A32:07, HLA-A68:23, HLA-C12:03, HLA-B42:01, HLA-B07:02, HLA-B40:13, HLA-A32:15
76.19%
Yes
Yes

YLKNHNIDL
26–34
1.7400
HLA-C12:03, HLA-B15:02, HLA-A02:50, HLA-A68:23, HLA-A32:07, HLA-C14:02, HLA-C03:03, HLA-B27:20, HLA-A02:12, HLA-A02:17, HLA-A02:03, HLA-C06:02, HLA-A32:15, HLA-B08:01, HLA-A02:02, HLA-C07:02
71.43%
Yes
Yes

YNVAWRTYL
19–27
0.4500
HLA-C03:03, HLA-A02:50, HLA-B27:20, HLA-B15:02, HLA-A68:23, HLA-A32:15, HLA-C12:03, HLA-C14:02, HLA-A32:07, HLA-A02:17
71.43%
No
Yes

LEPIIGDKL
54–62
0.4388
HLA-A02:17, HLA-A02:50, HLA-C12:03, HLA-A32:07, HLA-A68:23, HLA-B40:13, HLA-B27:20, HLA-B40:01, HLA-A32:15, HLA-B15:02
80.95%
No
Yes

CEEYGCLAI
1048–1056
1.6829
HLA-B27:20, HLA-C12:03, HLA-A32:07, HLA-B40:13, HLA-B40:01, HLA-A02:50, HLA-A68:23, HLA-C05:01, HLA-A32:15, HLA-C03:03
95.24%
Yes
Toxin

NSGPYNVAW
15–23
0.8498
HLA-B27:20, HLA-A32:07, HLA-C12:03, HLA-B40:13, HLA-A32:15, HLA-A68:23, HLA-C07:01, HLA-C03:03, HLA-B58:01, HLA-A02:50, HLA-B57:01
80.95%
No
Yes

RQHRGCIRF
43–51
1.0559
HLA-B27:20, HLA-B15:03, HLA-A32:07, HLA-B40:13, HLA-A68:23, HLA-B15:01, HLA-C12:03
71.43%
Yes
Toxin

AISGLGVGY
58–66
1.1548
HLA-A32:07, HLA-C12:03, HLA-A32:15, HLA-C03:03, HLA-A68:23, HLA-B27:20, HLA-B40:13, HLA-A30:02, HLA-B15:17
76.19%
No
Yes

The selected 18 T-cell epitopes were found to be recognized by the significant MHC class-I molecule such as HLA-A, HLA-B, and HLA-C according to IEDB MHC class 1-binding analysis resource. We selected humans as MHC source species and the SMM method for the prediction of a distinct set of MHC HLA alleles for the humans. This tool gives an output result for HLA-binding affinity of the epitopes in the IC50 nM unit. A lower IC50 value indicates higher binding affinity of the epitopes with the MHC class I molecule. So, in this study, we chose IC50 values less than 200 nM () for ensuring higher affinity (Table 1). For MHC-I processing, the analysis tool of the IEDB generates an overall score for each epitope’s intrinsic potential of being a T-cell epitope based on proteasomal processing, TAP transport, and MHC-I-binding efficiency. Before being presented to the T-cells on the plasma membrane of the cell, the protein is degraded into small peptides in the proteasome by the cytosolic proteases and MHC-1 forms a complex with the peptides as well. Then, the MHC-I peptide complex is transported to the endoplasmic reticulum (ER) through heat shock proteins and transport-associated proteins (TAP). However, the higher the total score of the epitopes with the HLA alleles ensures the presentation to the T-cell and a successful immune response critically depends on it. The epitopes with their respective alleles and total scores are summarized in Table 1 and Table S2. In this in silico analysis, the epitopes “CEEYGCLAI” and “RQHRGCIRF” were found as toxic to humans.

Effective immune response not only depends on the successful recognition of epitopes by HLA molecules with significant affinity but also depends on the antigenicity and immunogenicity score. So, the epitopes that were recognized by the considerable number of HLA alleles and contained the highest immunogenicity, antigenicity value, and nontoxic to human were considered as the potential epitope to induce a strong immune response. So, we selected two 100.0% conserved 9-mer epitopes “TSSWGCEEY” and “LAIDTGCLY” for docking purposes. Here, the epitope “TSSWGCEEY” showed the highest affinity for 11 MHC-1 molecules, including HLA-B27:20 (8.37), HLA-C03:03 (14.05), HLA-A32:07 (14.19), HLA-A68:23 (22.63), HLA-C12:03 (23.03), HLA-B40:13 (65.91), HLA-A32:15 (101.22), HLA-C07:01 (104.75), HLA-B15:02 (126.83), HLA-B15:17 (136.68), and HLA-B35:01 (180.13), whereas the epitope “LAIDTGCLY” exhibited the affinity with HLA-C03:03 (8.35), HLA-C12:03 (9.53), HLA-B15:17 (10.93), HLA-A68:23 (17.90), HLA-B35:01 (25.92), HLA-A32:07 (28.58), HLA-A29:02 (73.92), HLA-B15:01 (122.74), HLA-B15:03 (136.49), HLA-B15:02 (160.04), and HLA-B40:13 (169.80) (Table 1). So, 7 MHC-I alleles HLA-C03:03, HLA-A32:07, HLA-C12:03, HLA-B40:13, HLA-B15:02, HLA-B15:17, and HLA-B35:01 were found as common in both “TSSWGCEEY” and “LAIDTGCLY.” However, we selected HLA-B35:01 for docking purposes due to the availability of this allele in the PDB database. After that, we further used the epitopes “TSSWGCEEY” and “LAIDTGCLY” for the prediction of MHC-II alleles and their respective peptide or CD4⁺ T-cell epitope. The results showed that the epitope “TSSWGCEEY” and respective 15-mer peptides “KEHTSSWGCEEYGCL” and “EHTSSWGCEEYGCLA” have low affinity with the predicted MHC-II alleles (Table 2). On the other hand, the epitope “LAIDTGCLY” was found as the core sequence of 31 15-mer peptides or CD4⁺ T-cell epitopes in the range of the IC50 value 1–3000 nM. Among these 31 peptides, ten peptides exhibited a strong affinity (IC50 value 194 to 611 nM) with HLA-DR alleles, HLA-DRB301:01, HLA-DRB103:01, and HLA-DRB112:01 (Table 2). So, HLA-DR alleles could be the best-binding grooves for the predicted peptides.

Table 2: MHC-II alleles for the desired epitope.

Core sequence
MHC-II alleles
position
Peptide
IC50 value

TSSWGCEEY
HLA-DQA101:01/DQB105:01
1040–1054
KEHTSSWGCEEYGCL
2715

HLA-DQA101:01/DQB105:01
1041–1055
EHTSSWGCEEYGCLA
2801

LAIDTGCLY
HLA-DRB301:01
1049–1063
EEYGCLAIDTGCLYG
194

HLA-DRB301:01
1050–1064
EYGCLAIDTGCLYGS
195

HLA-DRB301:01
1052–1066
GCLAIDTGCLYGSCQ
195

HLA-DRB301:01
1051–1065
YGCLAIDTGCLYGSC
196

HLA-DRB103:01
1049–1063
EEYGCLAIDTGCLYG
401

HLA-DRB103:01
1050–1064
EYGCLAIDTGCLYGS
401

HLA-DRB103:01
1051–1065
YGCLAIDTGCLYGSC
403

HLA-DRB103:01
1052–1066
GCLAIDTGCLYGSCQ
405

HLA-DRB301:01
1054–1068
LAIDTGCLYGSCQDV
601

HLA-DRB301:01
1053–1067
CLAIDTGCLYGSCQD
611

HLA-DRB103:01
1054–1068
LAIDTGCLYGSCQDV
1173

HLA-DRB103:01
1053–1067
CLAIDTGCLYGSCQD
1200

HLA-DRB112:01
1051–1065
YGCLAIDTGCLYGSC
1450

HLA-DRB112:01
1050–1064
EYGCLAIDTGCLYGS
1457

HLA-DPA101:03/DPB102:01
1049–1063
EEYGCLAIDTGCLYG
1545

HLA-DPA101:03/DPB102:01
1051–1065
YGCLAIDTGCLYGSC
1551

HLA-DPA101:03/DPB102:01
1050–1064
EYGCLAIDTGCLYGS
1554

HLA-DRB112:01
1052–1066
GCLAIDTGCLYGSCQ
1561

HLA-DPA101:03/DPB102:01
1052–1066
GCLAIDTGCLYGSCQ
1571

HLA-DRB101:01
1052–1066
GCLAIDTGCLYGSCQ
1597

HLA-DPA102:01/DPB101:01
1050–1064
EYGCLAIDTGCLYGS
2005

HLA-DRB113:02
1050–1064
EYGCLAIDTGCLYGS
2117

HLA-DPA102:01/DPB101:01
1051–1065
YGCLAIDTGCLYGSC
2123

HLA-DRB113:02
1049–1063
EEYGCLAIDTGCLYG
2127

HLA-DRB113:02
1051–1065
YGCLAIDTGCLYGSC
2145

HLA-DPA102:01/DPB101:01
1052–1066
GCLAIDTGCLYGSCQ
2201

HLA-DPA102:01/DPB101:01
1049–1063
EEYGCLAIDTGCLYG
2287

HLA-DRB104:01
1050–1064
EYGCLAIDTGCLYGS
2673

HLA-DRB112:01
1053–1067
CLAIDTGCLYGSCQD
2760

HLA-DRB112:01
1054–1068
LAIDTGCLYGSCQDV
2917

HLA-DRB112:01
1049–1063
EEYGCLAIDTGCLYG
2949

3.3. Analysis of Population Coverage and Allergenicity

MHC HLA allele distribution differs among diverse geographic regions and ethnic groups around the world. Therefore, population coverage must be taken into consideration during the design of an effective vaccine. In this study, for the population coverage, identified MHC-I-binding alleles of 18 epitopes were considered. A significant population coverage was found for the selected 18 epitopes in different geographic regions of the world (Figure 2 and Table S3). These epitopes and their HLA-alleles cover 96.44% of the world population. The highest population coverage was found in the European country Finland (99.69%) which was closely followed by Ireland, Germany, England, Poland, France, Italy, and Sweden with a population coverage of 99.51%, 99.19%, 98.66%, 98.60%, 98.40%, 97.27%, and 96.99%, respectively. The lowest population coverage was found in South Africa (1.53%). The Oropouche fever was first found in South America, especially in Brazil, and several outbreaks have been recorded in this country. So, the population coverage prediction in Brazil is essential for vaccine design. The cumulative population coverage in Brazil is 93.92%, but it is found as 95.77% in mixed Brazilian people (Table S3). However, another South American country Chile has 92.79% population coverage. It varies in the North American country, especially in the United States people, and mostly covers the United States Caucasian race (98.96%) and United States African American (96.18%), United States Hispanic (95.86%), United States Mestizo (95.64%), United States Asian (94.58%), United States Polynesian (93.86%), and United States Amerindian (90.99%) people (Table S3). The population coverage in West Africa is 91.97%, whereas the North African country Morocco and Sudan covers 96.69% and 96.64%, respectively, also followed by the East African country Zambia (97.75%) and Kenya (94.81%). However, the Oceanian country, Australia (Caucasian race), covers about 99.30% populations for the selected 18 epitopes and their respective alleles (Figure 2).

Figure 2: Population coverage of the epitopes and their respective HLA alleles of the OROV polyprotein. The epitopes with different HLA-binding specificities increase the population coverage as the several HLA types are expressed at radically different frequencies in various ethnicities.

Before vaccine design, the allergenicity prediction is a key step, because most vaccines transfer the immune response to the allergic reaction by initiating immunoglobulin E and Type II T-helper cells [53]. In that case, we predicted the allergenicity of the selected epitope through the AllergenFP 1.0 server. The server identified the epitopes, TSSWGCEEY, CSMCGLIHY, VPRYHSIDV, YNVAWRTYL, LEPIIGDKL, NSGPYNVAW, and AISGLGVGY as nonallergen and the rest of the other epitopes as a potential allergen to humans (Table 1). On the contrary, AllerTOP v. 2.0 showed the allergenicity nature for the 100.0% conserved epitope TSSWGCEEY and 80.95% conserved epitope VPRYHSIDV. However, the results were found as the same for the rest of the other five epitopes (CSMCGLIHY, YNVAWRTYL, LEPIIGDKL, NSGPYNVAW, and AISGLGVGY) in both AllergenFP 1.0 and AllerTOP v. 2.0 servers (data not shown). The predicted three-dimensional structures of the T-cell epitopes were deposited to the Protein Model Database (PMDB). As a result, the database provided a unique accession number for each of the submissions of the epitope. The structure of the epitopes could be identified by using the accession number of PM0081248 and PM0081251 for T-cell epitopes TSSWGCEEY and LAIDTGCLY, respectively, in the PMDB database.

3.4. Docking Simulation Study

The molecular docking simulation study anticipated nine (9) potential binding affinities between the HLA-B35:01 and epitope TSSWGCEEY. At the same time, it was also found in the interaction between HLA-B35:01 and epitope LAIDTGCLY. The finest model was selected by higher binding affinity. The best interrelating model of epitope TSSWGCEEY and the HLA molecule contained a binding affinity of −8.0 kcal/mol. On the other hand, the epitope LAIDTGCLY and HLA molecule contained a binding affinity of −8.4 kcal/mol. The epitope position in the cartoon and surface structure and the binding interactions are shown in Figure 3. The epitope and the HLA molecule interacted by the hydrogen bonds, and we found nine hydrogen bonds after the docking simulation between HLA-B35:01 and epitope TSSWGCEEY. In this study, the amino acid residues, namely, Tyr-7, Tyr-9, Tyr-99, Tyr-159, Trp-147, Thr-69, and Thr-73, acted as the hydrogen bond donors or acceptors with the epitope residues (Figure 3). Herein, eight hydrogen bonds, 89 nonbonded contacts, and the interacting residues, Thr-143, Ser-77, Gln-155, Tyr-9, Asn-70, and Tyr-59, were found for HLA-B35:01 and epitope LAIDTGCLY.

Figure 3: Docking simulation study (sample). (a) The docking results as the surface structure of the MHC-I HLA allele “HLA-B35:01” (tint light pink color) and the sticky form of epitope “TSSWGCEEY” (orange olive color). The black color box represents the position of epitope in the surface structure. (b) Epitope in the cartoon structure. (c) Sticky form of epitope “TSSWGCEEY” (orange olive color) with residue’s position. (d) Binding interaction of the epitope residues and the HLA-B35:01 residues through the hydrogen bonds. Hydrogen bonds and the distances are shown in green color. (e) The docking results as the surface structure of the MHC-I HLA allele “HLA-B35:01” (tint light blue color) and the sticky form of epitope “LAIDTGCLY” (orange olive color). The black color box represents the position of epitope in the surface structure. (f) Epitope LAIDTGCLY in the cartoon structure. (g) Sticky form of epitope “LAIDTGCLY” (orange olive color) with residue’s position. (h) Binding interaction of the epitope residues and the HLA-B35:01 residues through the hydrogen bonds. Hydrogen bonds and the distances are shown in green color.

The cartoon and surface structures were prepared and visualized by the PyMOL molecular graphics system.

In the control section, the docking between HLA-B35:01 and the influenza NP418 epitope (LPFERATVM) from the 2009 strain was performed. The best affinity or binding energy between these two molecules was found −8.3 kcal/mol, which is very close to the sample. The amino acid residues, namely, Tyr-7, Tyr-9, Tyr-99, Tyr-159, Tyr-171, Ser-77, and Gln-155, are involved in the interaction with the epitope residues (Figure S1). So, the results revealed that Tyr-7, Tyr-9, Tyr-99, and Tyr-159 are conserved residues in the docking simulation of the sample (TSSWGCEEY) and control and ensuring the accuracy of the binding affinity. Besides, Tyr-9, Ser-77, and Gln-155 residues were found as common residues in the docking simulation of the sample LAIDTGCLY (binding affinity of −8.4 kcal/mol) and control which also ensures the accuracy and the acceptance of the docking simulation.

The epitope (peptide) binding to the MHC class II HLA-DR molecule is important for the proper functioning of a vaccine candidate. So, we predicted the alleles for our selected epitope and performed a docking study between those molecules. Herein, we found the best binding affinity −7.9 kcal/mol for the docking study among the predicted nine binding affinity results. This result indicates a good affinity between HLA-DR and the epitope LAIDTGCLY (Figure 4). During the analysis of the interaction between HLA-DR and LAIDTGCLY, we found five hydrogen bonds and 125 nonbonded contacts in the combined structure. However, the amino acid residues, namely, Ser-53, Gln-9, Asn-62 (2 times), and His-13, are involved in the interaction with the epitope residues (Figure 4).

Figure 4: Docking simulation study. (a) The docking results as the surface structure of the MHC-II (HLA-DR) (lime green color) and the sticky form of epitope “LAIDTGCLY” (red color). Black color box represents the position of epitope in the surface structure. (b) Epitope in the cartoon structure. (c) Sticky form of epitope “LAIDTGCLY” (red color) with residue position. (d) Binding interaction of the epitope residues and the HLA-DR residues through the hydrogen bonds. Hydrogen bonds and the distances are shown in green color.

3.5. Analysis of Linear and Conformational B-Cell Epitopes

A B-cell epitope is a precise portion of an antigen recognized by either a specific B-cell receptor or the elicited antibody in a humoral response [23, 54]. The B-cell epitopes are divided into two main categories such as continuous or linear B-cell epitope and discontinuous or conformational B-cell epitope. It has been found that most of the B-cell epitopes are conformational epitopes, and the amount of this epitope is more than 90% [55].

The prediction of B-cell epitopes from the antigenic protein is the major step of epitope-based vaccine design. Thus, to get essential B-cell epitope candidates in the polyprotein of OROV, we performed an in silico analysis through the web server BepiPred-2.0 and LBtope.

Based on the BepiPred-2.0 method, a total of 31 linear B-cell epitopes with varied length were identified from the polyprotein of OROV. Among all these predicted epitopes, 8 B-cell epitopes ⁸¹VKDWSECNPVLDKL⁹⁴, ³³⁴AGERMYKLEELADD³⁴⁷, ⁵⁰⁴YIMGDTCVSKSELYSDLLSKNLITEYDKKYFEPDTVN⁵⁴⁰, ⁸⁷⁵CSEQRFPIHKNWFKKCNWDKTHKEFTTMRQINYNDITSYRKAIES⁹¹⁹, ⁹²²GTDLMTHHYKPTKNLPHVVPRYH⁹⁴⁴, ¹⁰⁹⁵ETYCNDMDILEPIIGD¹¹¹⁰, ¹¹⁸⁶FNDNYQSCT¹¹⁹⁴, and ¹²³¹DINYKIFTKSV¹²⁴¹ were nonantigenic as predicted by the VaxiJen v2.0 server. These epitopes had 5.56% to 55.56% conservancy level, except FNDNYQSCT (100% conservancy score), among the specified polyprotein sequences (Table 3). So, these epitopes could not be the vaccine candidates. On the other hand, the epitopes ¹⁷HPLSTSQIGDRC²⁸ (94.44%), ³⁵FKEMNLSVGL⁴⁴ (100%), ²⁸²SCTEALKVHRMGKDCLGYKS³⁰¹ (100%), ⁵⁶⁶ECEMLSLKKNSGPYN⁵⁸⁰ (100%), ⁵⁹⁴LCSRHNY⁶⁰⁰ (100%), ¹⁰²⁶PEKIPAKEGWLTFSKEHTSSW¹⁰⁴⁶ (100%), and ¹¹¹⁹TQTNQLP¹¹²⁵ (100%) can be considered as vaccine candidates due to their antigenicity and high conservancy score (Table 3). After the prediction of allergenicity (AllerTOP 2.0 and AllergenFP v. 1.0) and toxicity, it can be said that the epitopes ¹⁷HPLSTSQIGDRC²⁸, ¹⁰²⁶PEKIPAKEGWLTFSKEHTSSW¹⁰⁴⁶, and ⁷⁸⁸SHCNLEFTAITADKIMSL⁸⁰⁵ could be the potential linear B-cell epitopes for vaccine design (Table 3).

Table 3: Linear B-cell epitopes of the polyprotein of OROV predicted by BepiPred-2.0.

B-cell epitopes
Position
Antigenicity score
Conservancy (%)
Toxicity
Allergenicity

AllerTOP 2.0
AllergenFP 1.0

HPLSTSQIGDRC
17–28
1.2943
94.44
No
No
No

FKEMNLSVGL
35–44
1.9084
100.00
No
Yes
Yes

FSKNTLALEA
62–71
0.4645
5.56
No
Yes
No

VKDWSECNPVLDKL
81–94
−0.0129
5.56
No
No
No

SCTEALKVHRMGKDCLGYKS
282–301
0.4510
100.00
No
Yes
No

AGERMYKLEELADD
334–347
0.2927
38.89
No
Yes
No

IHYRPGLKVDLTK
404–416
0.4322
33.33
No
Yes
No

FDEQQSSGFEYEVFLKDMHIQRESCKFGPR
427–456
0.6740
5.56
No
Yes
Yes

DCLSKNIKITYQELH
484–498
1.2448
11.11
No
Yes
No

YIMGDTCVSKSELYSDLLSKNLITEYDKKYFEPDTVN
504–540
0.1143
5.56
No
No
No

ECEMLSLKKNSGPYN
566–580
0.6930
100.00
No
Yes
No

LCSRHNY
594–600
1.2777
100.00
No
Yes
Yes

HSMCKNTDIDFNNEIETYYKANAAAY
610–635
0.7190
5.56
No
Yes
No

NISLGRSTRSAQHSIMTNEIPKSNPFTDYSYSSVNIKECMSP
699–740
0.6399
5.56
No
Yes
Yes

KKRDSVPHT
747–755
0.9721
11.11
No
No
Yes

KYKAFDWPEIETVQKGQKLCL
765–785
1.3196
33.33
No
Yes
No

SHCNLEFTAITADKIMSL
788–805
1.0518
83.33
No
No
No

NCYKESFTAQPADMQTGI
807–824
0.8661
61.11
No
No
No

CITLEDKA
835–842
2.5070
5.56
No
No
Yes

KYYYSDGKEHAKDGSIN
851–867
0.6483
33.33
No
No
No

CSEQRFPIHKNWFKKCNWDKTHKEFTTMRQINYNDITSYRKAIES
875–919
0.1740
27.78
No
No
No

GTDLMTHHYKPTKNLPHVVPRYH
922–944
0.3102
72.22
No
No
No

AQYQKVYTTGPSISINIEHNEK
999–1020
0.6311
33.33
No
Yes
No

PEKIPAKEGWLTFSKEHTSSW
1026–1046
0.5512
100.00
No
No
No

ETYCNDMDILEPIIGD
1095–1110
0.1167
100.00
No
No
Yes

TQTNQLP
1119–1125
0.9980
100.00
No
Yes
Yes

TGQINDVGNTALQCGSIQVVNGSTIGSGNPK
1137–1167
0.8145
33.33
No
Yes
Yes

FNDNYQSCT
1186–1194
0.0858
100.00
No
No
Yes

DARNDLIPYRKG
1197–1208
1.4056
5.56
No
No
No

DINYKIFTKSV
1231–1241
0.3093
55.56
No
No
No

TVKPYNQKIDLSQLDESNYIKEEDLQCG
1322–1349
0.8134
38.89
No
No
No

During the prediction by LBtope, a total of 19 linear B-cell epitopes were found from the OROV polyprotein. The epitopes ¹³IAYGHPLSTSQI²⁴, ¹¹⁰PRMYTCRAACDI¹²¹, ⁴³²SSGFEYEVFLKDMH⁴⁴⁵, ⁵¹⁹DLLSKNLITEYDKKYFEPDTVNDQFNKIEFAQD⁵⁵¹, ⁶¹⁹DFNNEIETYYKANAAAYRSDFYTI⁶⁴², ⁷⁰⁶TRSAQHSIMTN⁷¹⁶, ⁸⁴⁸CGDKYYYSDGKEHAKDGS⁸⁶⁵, and ¹⁰⁸⁹CITLPHETYCNDMDILEPII¹¹⁰⁸ are nonantigenic among the 19 epitopes due to their low antigenicity score in the VaxiJen server (Table 4). On the contrary, the epitopes ⁸⁸⁹KCNWDKT⁸⁹⁵, ⁹²⁸HHYKPTKNLPHVVPRYH⁹⁴⁴, ¹²⁷³CHCKTNCED¹²⁸¹, and ¹³²⁷NQKIDLSQL¹³³⁵ were found as antigenic and 100% conserved among all the predicted epitopes. However, the epitopes ⁸⁸⁹KCNWDKT⁸⁹⁵ and ¹²⁷³CHCKTNCED¹²⁸¹ could not be the potential vaccine candidates because these two B-cell epitopes were identified as toxic to the human (Table 4).

Table 4: Linear B-cell epitopes of the polyprotein of OROV predicted by LBtope.

B-cell epitopes
Position
Antigenicity scores
Conservancy (%)
Toxicity
Allergenicity

AllerTOP 2.0
AllergenFP 1.0

IAYGHPLSTSQI
13–24
0.3474
50.00
No
Yes
No

TTAFSKNTLALEAT
59–72
0.6595
5.56
No
Yes
No

PRMYTCRAACDI
110–121
−0.5422
5.56
No
Yes
No

GLKVDLTKTN
409–418
0.5124
33.33
No
Yes
Yes

SSGFEYEVFLKDMH
432–445
0.1548
38.89
No
No
No

QRESCKFGPRVNHFRNTK
447–464
1.1279
38.89
No
No
No

DLLSKNLITEYDKKYFEPDTVNDQFNKIEFAQD
519–551
0.2055
33.33
No
No
No

DFNNEIETYYKANAAAYRSDFYTI
619–642
0.1701
5.56
No
Yes
No

TRSAQHSIMTN
706–716
0.1229
38.89
No
Yes
No

HTNHLLCKMDNKYKAFDWPEIETVQK
754–779
1.0858
16.67
No
Yes
No

CITLEDKAWP
835–844
2.0270
5.56
No
No
Yes

CGDKYYYSDGKEHAKDGS
848–865
0.3358
33.33
No
Yes
No

KCNWDKT
889–895
0.4951
100.00
Yes
Yes
No

HHYKPTKNLPHVVPRYH
928–944
0.6022
100.00
No
No
No

VYTTGPSIS
1004–1013
0.4965
38.89
No
Yes
Yes

CITLPHETYCNDMDILEPII
1089–1108
0.2420
100.00
No
Yes
Yes

DNYQSCTRLDARNDLIPYRKG
1188–1208
1.0163
5.56
No
Yes
No

CHCKTNCED
1273–1281
1.1907
100.00
Yes
No
Yes

NQKIDLSQL
1327–1335
0.5667
100.00
No
Yes
Yes

After the similarity search between BepiPred-2.0 and LBtope predicted B-cell epitopes, we identified 18 sequences as most conserved B-cell epitopes. However, 7 B-cell epitopes were found as nonantigenic and less conserved (5.56% to 38.89%), except the epitopes ¹⁰⁹⁵ETYCNDMDILEPII¹¹⁰⁸ (100%) and ¹¹⁸⁸DNYQSCT¹¹⁹⁴ (100%) (Table 5). The predicted epitope ⁹²⁸HHYKPTKNLPHVVPRYH⁹⁴⁴ is 100% conserved, which is also followed by ¹³²⁷NQKIDLSQL¹³³⁵ (100%), ⁸⁸⁹KCNWDKT⁸⁹⁵ (100%), and ¹⁷HPLSTSQI²⁴ (94.44%). The analysis of antigenicity, conservancy, toxicity, and allergenicity of B-cell epitopes revealed that the epitope ⁹²⁸HHYKPTKNLPHVVPRYH⁹⁴⁴could be the most potential B-cell epitope candidate for peptide-based vaccine design because this epitope is antigenic, 100% conserved, nontoxic, and nonallergenic to the human. The exact position and length of the predicted B-cell epitopes are shown in Table 5.

Table 5: Most conserved B-cell epitope sequences predicted by both BepiPred-2.0 and LBtope.

B-cell epitopes
Position
Length
Antigenicity score
Conservancy (%)
Toxicity
Allergenicity

AllerTOP 2.0
AllergenFP 1.0

HPLSTSQI
17–24
8
0.7973
94.44
No
Yes
Yes

FSKNTLALEA
62–71
10
0.4645
5.56
No
Yes
No

GLKVDLTK
409–416
8
−0.8608
33.33
No
Yes
Yes

SSGFEYEVFLKDMH
432–445
14
0.1548
38.89
No
No
No

QRESCKFGPR
447–456
10
2.2285
38.89
No
No
Yes

DLLSKNLITEYDKKYFEPDTVN
519–540
22
0.2587
33.33
No
No
No

DFNNEIETYYKANAAAY
619–635
17
0.1778
5.56
No
No
No

TRSAQHSIMTN
706–716
11
0.1229
38.89
No
Yes
No

KYKAFDWPEIETVQK
765–779
15
1.1837
33.33
No
Yes
No

CITLEDKA
835–842
8
2.5070
5.56
No
No
Yes

KYYYSDGKEHAKDGS
851–865
15
0.5485
33.33
No
No
No

KCNWDKT
889–895
7
0.4951
100.00
Yes
Yes
No

HHYKPTKNLPHVVPRYH
928–944
17
0.6022
100.00
No
No
No

VYTTGPSIS
1004–1013
9
0.4965
38.89
No
Yes
Yes

ETYCNDMDILEPII
1095–1108
14
0.1553
100.00
No
No
No

DNYQSCT
1188–1194
7
−0.1560
100.00
No
No
Yes

DARNDLIPYRKG
1197–1208
12
1.4056
5.56
No
No
No

NQKIDLSQL
1327–1335
9
0.5667
100.00
No
Yes
Yes

In this study, a total of 8 discontinuous or conformational B-cell epitopes were predicted using the ElliPro tool of IEDB. The epitopes having a protrusion index (PI) score above 0.8 were taken into consideration. The PI value of the 8 predicted epitopes ranged from 0.837 to 0.906. The residues with higher scores indicate better solvent accessibility. Herein, the highest probability of a discontinuous epitope was calculated at 90.6% (PI score: 0.906) and is shown in Figure 5(a). Conformational epitopes and their individual residues, residue position, length, and the scores are shown in Table 6, whereas the positions of epitopes on 3D structures are exhibited in Figures 5(a)–5(h).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 5: Three-dimensional representation of conformational or discontinuous epitopes (a–h) of the highest antigenic polyprotein of OROV. The epitopes are represented by yellow surface, and the bulk of the polyprotein is represented in grey sticks.

Table 6: Predicted conformational B-cell epitopes of the polyprotein of OROV.

Epitope no.
Residues and their positions
No. of residues
Score
3D structure

1.
V992, K993, K994, A995, V996, Y997, K998, A999, Q1000, Y1001, Q1002, K1003, V1004, T1006, T1007, G1008, P1009, S1010, I1011, S1012, I1013, N1014, I1015, E1016, H1017, N1018, E1019, K1020, C1021, T1022, G1023, H1024, C1025, P1026, E1027, K1028, I1029, P1030, A1031, K1032, E1033, G1034, W1035, L1036, T1037, K1040, E1041, S1044, S1045, W1046, G1047, C1048, E1049, E1050, Y1051, G1052, C1053, L1054, A1055, I1056, D1057, T1058, G1059, C1060, L1061, Y1062, G1063, S1064, C1065, Q1066, D1067, V1068, I1069, R1070, P1071, E1072, L1073, D1074, I1075, K1077, K1078, I1079, G1080, S1081, E1082, D1110, K1111, I1140
88
0.906
Figure 5(a)

2.
Y76, R77, I80, V81, D83, W84, S85, E86, C87, N88, P89, V90, L91, D92, K93, L94, G95, N96, M98, N105, G106, H107, L108, V109, P110, R111, M112, Y113, T114, C115, R116, A117, A118, C119, D120, I121, R122, L123, N124, K125, D126, N127, A128, E129, I130, I131, L132, S133, S134, T135, K136, L137, N138, H139, F140, E141, I142, V143, G144, T145, T146, S147, T148, S149, G150, W151, F152, K153, N154
69
0.895
Figure 5(b)

3.
M738, S739, P740, S742, C760, K761, M762, D763, N764, K765, Y766, K767, A768, F769, D770, W771, P772, E773, I774, E775, T776, V777, Q778, K779, G780, Q781, K782, L783, C784, L785, G786, D787, S788, H789, C790, K810, E811, S812, F813, T814, A815, Q816, P817
43
0.875
Figure 5(c)

4.
Y563, T565, E566, C567, E568, M569, L570, S571, L572, K573, K574, N575, S576, G577, P578
15
0.864
Figure 5(d)

5.
G43, L44, G45, E46, I47, C48, V49, K50, D51, D52, I53, S54, I55
13
0.845
Figure 5(e)

6.
I1385, I1388, V1389, L1390, V1394, L1395, L1396, P1397, I1398, G1399, R1400, L1401, L1402, K1403, A1404, F1405, L1406, I1407, K1408, E1412
20
0.843
Figure 5(f)

7.
M292, D295, C296, L297, G298
5
0.838
Figure 5(g)

8.
P332, I333, A334, G335, E336, R337, M338, Y339, K340, E342, E343, A345, D346, F387, N388, Y389, L390, S397, C398, S399, M400, C401, G402, L403, I404, H405, Y406, R407, P408, G409, L410, K411, V412, D413, L414, T415, K416, T417, N418, K419, C420, G421, S422
43
0.837
Figure 5(h)

3.6. Ligand-Binding Pocket Analysis

The physicochemical properties of the OROV polyprotein were found as 8.01 isoelectric points, 37.22 instability index value, 88.99 aliphatic index values, and −0.138 GRAVY value. The other features are shown in Table 7. These results revealed that the protein is naturally mildly basic, stable, and hydrophilic. The secondary structural features showed that the random coils with 45.99% is dominated by alpha helix (31.62%) and extended strands (22.39%) (Figure S2). The unstructured or disordered regions can play an important role in the catalytic activity of the protein that was predicted by the DISOPRED3 server. The disordered regions were found in the ⁵²⁸EY⁵²⁹, ⁵³⁵E, ⁷⁰⁵STRSAQHS⁷¹², ⁷²³PF⁷²⁴, and ⁷²⁷YS⁷²⁸ positions. The disordered regions and their confidence score are shown in Figure 6.

Table 7: Different physicochemical properties of the Oropouche virus (OROV) polyprotein.

Features
Value

Number of amino acids
1420

Molecular weight
161275.63

Theoretical pI
8.07

Total number of negatively charged residues (Asp + Glu)
151

Total number of positively charged residues (Arg + Lys)
161

Total number of atoms
22,605

Extinction coefficients
177,070

Instability index
37.22

Aliphatic index
88.99

Grand average of hydropathicity (GRAVY)
−0.138

Figure 6: Disordered regions of the OROV polyprotein. Here, amino acids in the input sequence are considered disordered when the blue line is above the grey dashed line, that is, the confidence score is higher than 0.5. The orange line shows the confidence of disordered protein-binding residue predictions.

The 3D structure of the OROV polyprotein is not available in the PDB database. So, a 3D structure of the OROV polyprotein was predicted by the I-TASSER server and visualized by the PyMOL tools (Figure 7). Different evaluation scores of the I-TASSER server, such as C-score (−1.02), estimated TM score (0.59 ± 0.14), estimated RMSD (12.3 ± 4.4 Å), cluster density (0.1329), and number of decoys (213), ensured the good quality of the predicted model. Energy minimization through the ModRefiner server and Swiss-PdbViewer software tools removes the bad contacts of the predicted protein model. The total force field energy of the protein structure was −36896.734 kJ/mol after energy minimization. Here, PROCHECK was used to check the stereochemical quality of the structure by considering residue geometry and overall structural geometry. The results showed 72.2% in most favored regions, 21.9% in the additional allowed region, 3.6% in the generously allowed region, and 2.3% in the disallowed region in the Ramachandran plot statistics (Figure S3). To check the potential errors of the protein 3D model, ProSA was used, and it predicts the negative -score of −6.67, which ensures the good quality of the model (Figure S4). The ProQ server showed an LG score of 3.651 (LG indicates very good model) and MaxSub score of 0.500 ( indicates very good model), ensuring the “very good” quality of the predicted protein model. Another evaluation tool QMEAN gave different scores with an environmental profile graph for the OROV polyprotein structure (Figure S5). The refined and evaluated protein model was submitted to the PMDB database, and the server provided a unique accession id PM0081042 for the submitted model. The hydrophobicity nature greatly impacts the drug design and delivery system. The hydrophobicity plot of the OROV polyprotein is shown in Figure S6.

Figure 7: Three-dimensional structure of the OROV polyprotein. This protein model was predicted by the I-TASSER server and visualized by the PyMOL molecular graphics system.

The ligand-binding pockets of OROV polyprotein were not reported yet. So, the SiteHound server was used to predict the top 10 ligand-binding pockets of the OROV polyprotein which ranked by energy (Table 8). The surface volume of the binding pockets was also explored in the center -axis, -axis, and -axis. The residues which are involved in the formation of protein cavities were explored for copious binding of new ligands such as epitope. The energy was ranged from −1677.27 to −940.23 which ensures the efficiency of the predicted top 10 ligand-binding pockets. These ligand-binding pockets could be the target for molecular docking simulation in the future.

Table 8: Top 10 binding pockets with binding residues, energy, and volume of OROV polyprotein.

Rank
Energy (kcal/mol)
Energy Range
Volume (Å)
Certer X
Certer Y
Certer Z
Residues

1
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Gly-323, Leu-324, Leu-326, Gln-430, Gln-431, Ser-432, Ser-433, Gly-434, Phe-435, Glu-436, Phe-453, Lys-647, Thr-648, Ala-649, Phe-650, Cys-835, Ile-836, Leu-838

2
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Asp-750, Ser-751, Val-752, Ser-828, Asp-830, Ile-832, Gly-833, Ile-836, Thr-837, Cys-870, Leu-871, Ala-872, Asn-873, Lys-874, Glu-877, Phe-880, Pro-881, His-928, His-929, Tyr-930, Lys-931, Pro-932

3
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Lys-744, Cys-745, Phe-746, Lys-747, Lys-748, Leu-758, Leu-759, Cys-760, Thr-776, Gln-778, Lys-779, Gly-780, Gln-781, Lys-782, Leu-783, His-789, Phe-813

4
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Leu-34, Phe-35, Lys-36, Glu-37, Asn-39, Phe-178, His-179, Cys-181, Ala-264, Ser-265, Pro-267, Phe-268, Phe-281, Cys-309, Ser-311, Lys-312, Trp-314

5
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Phe-79, Asp-83, Glu-86, Met-98, Val-99, Leu-100, Asn-101, Val-102, Asp-103, Asn-105, Gly-106, Val-109, Tyr-113, Arg-116, Asp-120, Asp-126, Asn-127, Ala-128, Glu-129, Ile-131, Leu-132, Asn-138, Ile-142

6
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Ile-472, Ser-475, Phe-476, Val-581, Arg-584, Thr-585, Leu-587, Lys-588, Ala-630, Ala-632, Ala-633, Arg-636, Ser-637, Asp-638, Phe-639, Tyr-640, Thr-641

7
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Ser-912, Tyr-913, Ile-958, Asn-959, Gly-960, Phe-961, Ile-962, Gln-963, Cys-1301, Lys-1302, Ser-1303, Lys-1304, Lys-1324, Leu-1332, Ser-1333, Gln-1334, Phe-1364, Lys-1365, Gly-1366

8
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Ile-962, Gln-963, Asn-964, Thr-965, Ile-966, Lys-1304, Asn-1311, Gln-1320, Leu-1321, Thr-1322, Lys-1324, Leu-1367, Phe-1368, Ser-1369, Tyr-1413

9
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Tyr-661, Val-690, Leu-691, Lys-692, Phe-693, Ala-694, Leu-697, Leu-698, Ser-701, Arg-914, Lys-915, Glu-918, Ser-919, Gly-922, Thr-923

10
−1677.27
(−14.79, −8.90)
156.00
123.724
99.245
115.842
Ile-255, Asn-256, Ser-257, Ser-259, Asn-260, Cys-261, Leu-262, Phe-268, Cys-271, Pro-272, Cys-283, Thr-284, Ala-286, Leu-287, His-290, Lys-310, Ser-311, Lys-312

4. Discussion

In recent years, many diseases have emerged due to the occurrence of several outbreaks through the different types of newer viruses. So, vaccine development against these emerging diseases within a short time is very crucial to protecting the people from the rising viral attacks. Vaccines are the pharmacological products which can provide the finest cost-benefit ratio in the prevention or treatment of diseases. However, an effective vaccine progression and production are costly and can take years to be completed. So, the researchers have tried for many years to minimize the cost and time for the development of vaccines. At this time, there are different strategies available for the design and development of effective and safe new-generation vaccines based on the Bioinformatics approaches [56, 57]. The next-generation sequencing and progressive genomics and proteomics technologies have brought about a great change in computational immunology. However, the advancement of newer immunoinformatics tools has made a broader way in developing the vaccine or vaccine candidates through the satisfactory understanding of the immune response of the human body against an organism within a short time [58–60].

After recovery from an acute viral or microbial infection, a person generally achieves long-term or even lifelong immunity to the same viral disease or pathogen [61, 62]. Numerous acute viral infections exhort a protective humoral response which is characterized by the presence of a preexisting antibody. However, the duration of this preexisting antibody response differs among viral infections and can exist for a few months to many years in the host [63]. Slifka and Ahmed reported that the antibody persistence of acute viral infection in humans, such as Chikungunya, Rift Valley fever virus, dengue, yellow fever, and measles, are about 30 years, 12 years, 32 years, 75 years, and 65 years, respectively [63]. According to Tesh and Vasconcelos, OROV could not be isolated from the patient’s serum during the recurrent illness, but the detectable humoral antibodies are usually found [64]. However, no significant fatalities have been stated for OROV, and the lifelong immunity ensures the recovery from this emerging virus [64].

B-cell memory consists of humoral memory (long-lived antibody-secreting plasma cells) and reactive memory (long-lived memory B-cells which can react quickly to a recurrent antigenic challenge) [61]. Mainly memory B-cells afford the potential to respond fast and with more adapted antibodies against reinfection [61, 65]. Welsh et al. concluded that memory B-cell responses are long-lived even in the absence of antigen, but these responses are controlled and influenced by other infections in the host. This memory B-cells are responsible for providing long-term protective immunity or contribute to immunopathology on infection with homologous or heterologous viruses [66]. B-cell memory is virus-specific, but the exact cellular mechanisms of virus-specific antibody response to acute viral infection remain unclear [67]. Hebeis et al. demonstrated that the activation of virus-specific memory B-cells to secrete immunoglobulin G (Ig G) is autonomous of cognate or bystander T-cell help [68]. The role of virus-specific B-cell memory has been found in mumps virus [69], human influenza virus [70], rotavirus [71, 72], and dengue virus [73].

In recent trends, the surface glycoprotein (GP) is considered as the initial focus for the development of epitope-based peptide vaccine, as they are involved in the interaction between cell receptor and virus particle, thus playing a significant role in the pathogenesis of the disease [74]. To our knowledge, the subunit or epitope/peptide-based vaccine approaches have not been used yet for the prediction of a potential vaccine against this neglected emerging OROV. However, the virus-like particle assay has been developed for this OROV that could be used as a tool to develop a potential recombinant vaccine [11, 74]. The population coverage information of the epitopes and their respective HLA alleles are crucial to being the potential vaccine candidates. The computational prediction of this HLA allele is much easier and cost-effective than the experimental one. Therefore, in this study, different types of immunoinformatics and molecular simulation tools were used to find the potential vaccine candidates present in the OROV polyprotein.

In the present study, an immunoinformatics-driven approach was incorporated to screen emergent immunogen against OROV proteome. The results revealed that the OROV polyprotein was the most antigenic protein with the highest antigenic scores (0.5124) compared with other viral proteins of OROV. To our knowledge, no immunoinformatics study for the OROV polyprotein has been performed yet for the identification of a potential vaccine target. Although we found the polyprotein as a potential antigenic protein, another OROV protein may be the source of vaccine candidate identification. However, we identified the potential T-cell epitopes from the most antigenic polyprotein, as they play a key role in the creation of a defensive immune response against different pathogenic infections [75]. Various successful studies have been performed for the epitope-based peptide vaccine design against West Nile virus [76], Zika virus [77], dengue virus [78], Chikungunya virus [79], Rift valley fever virus [80], shigellosis [81], and so on. Primarily, we identified 146 epitopes in total using 12 supertypes (A1, A2, A3, A24, A26, B7, B8, B27, B39, B44, B58, and B62) and selected 128 epitopes based on the highest combined score. Finally, we identified 18 epitopes by maintaining some critical criteria such as antigenicity, immunogenicity, and most important conservancy scores of the epitopes (Table 1). We gave priority to the higher conservancy score of the epitopes because higher conservancy ensures the extensive protection across numerous strains or even species. Our selected 18 T-cell epitopes were 70.0% to 100.0% conserved which could ensure the acceptance of these predicted epitopes and indicate them as a potential vaccine candidate. Hossain et al. showed 65.66% conserved epitope “TLKNLNDNY” as a potential vaccine candidate for universal vaccine design against rotavirus [82]. On the contrary, Hasan et al. revealed 69.49% conserved epitope “YYYELYPTM” as a possible vaccine lead for vaccine design against Chikungunya virus [79]. Since our predicted 18 epitopes are highly conserved (3 epitopes—100% conserved, 1 epitope—95.24% conserved, 8 epitopes—80.95% conserved, 3 epitopes—76.19% conserved, and 3 epitopes—71.43% conserved), these epitopes could be considered as potential vaccine candidates for vaccine design against OROV. Predicted 100.0% conserved epitopes TSSWGCEEY, CSMCGLIHY, and LAIDTGCLY were found as antigenic, but the antigenicity was lower than the epitopes, CEEYGCLAI (95.24%), LIHYRPGLK (80.95%), VPRYHSIDV (80.95%), IPAISGLGV (76.19%), and YLKNHNIDL (71.43%). However, the epitope CSMCGLIHY could be the most important epitope in comparison with the rest of the other epitopes, because it is not allergenic (as predicted by both AllergenFP 1.0 and AllerTOP v. 2.0 servers) and nontoxic to the host. At present, the researchers are using different immunoinformatics approaches for the prediction of high potential CD4⁺ T-cell epitopes and CD8⁺ T-cell epitopes which can bind to different MHC-I and MHC-II molecules. de Freitas e Silva et al. identified 10 top-ranked epitopes from the proteome of Leishmania braziliensis using immunoinformatics tools and also synthesized and assessed in vitro to stimulate peripheral blood mononuclear cells (PBMC) from posttreated cutaneous leishmaniasis patients and found that 50% of these epitopes are immunogenic [83]. Khan et al. identified the T-cell epitopes from the Mycobacterium tuberculosis proteome using the in silico approaches for in vitro analyses and found T-cell responses during experimental evaluation of the identified epitopes [84]. So, the analysis clearly revealed that the immunoinformatics or in silico approaches could make a great contribution in vaccine lead identification for experimental study and/or the design of an effective vaccine.

The HLA alleles maintain the response to T-cell epitopes, and these alleles are extremely polymorphic in numerous ethnic populations. Therefore, MHC-I and MHC-II HLA allele specificity of T-cell epitopes should be considered as the initial criteria to induce extensive immune responses in diverse ethnic human populations [85, 86]. However, to get more population coverage, the T-cell epitope should bind with more MHC HLA alleles. So, we took the selected 18 epitopes and their respective MHC-I-binding HLA alleles for significant population coverage. The results revealed that the selected epitopes and their alleles have ideal population coverage in different geographic regions around the world. The highest population coverage was recorded at 99.69% in Finland, and those epitopes and their respective HLA-alleles cover 96.44% of the world population. The OROV outbreak occurred in the South American countries, especially in Brazil in the highest proportions. So, the vaccine candidates for this geographic region are essential to protecting individuals from OROV infection. The population coverage was found at 95.77% among the mixed Brazilian population, which indicates the acceptability of the epitopes for the OROV epidemic region. It was observed that the population coverage (1.53%) was found lower than the other countries (Figure 2). South Africa consists of different ethnic groups, but the HLA distribution in the allele frequency net database (AFND) is comparatively low which reflects the underrepresentation of this region. Tshabalala et al. reviewed that HLA-C alleles are mostly found in Caucasian South Africans, Black South Africans, Shona Zimbabweans, Tamil South Africans, and Zambians while HLA-B allele frequencies are found in Shonas and Black South Africans [87]. Cao et al. reported 9 different HLA alleles, HLA-A 02:02, HLA-B15:03, HLA-A34:02, HLA-A36:01, HLA-B42:01, HLA-B53:01, HLA-B57:03, HLA-B58:02, and HLA-A74:01, as unique African alleles [88]. On the other hand, a very low amount of HLA class-II alleles is available in South African populations. Among all HLA class-II alleles, HLA-DRB1 alleles are available in most South Africans, while HLA-DP data were found in only Zimbabwean Shonas and Black South Africans, but no alleles were reported for DQA1 in AFND [87, 89]. So, the unavailability of our predicted HLA-alleles in the South African population could be the main cause of the lowest population coverage in South Africa. The 3D structure of the epitopes TSSWGCEEY and LAIDTGCLY were designed for docking purposes and compared with a control (HLA-B35:01 and influenza NP418 epitope) to measure the docking efficiency with a specific HLA allele HLA-B35:01. The docking results showed almost similar binding affinity or energy for both the control (−8.3 kcal/mol) and the sample (−8.4 kcal/mol and −8.0 kcal/mol for the epitopes “LAIDTGCLY” and “TSSWGCEEY,” respectively) which designates the satisfactory accuracy of the predicted epitopes to form an interaction with the MHC-I allele. However, another docking was performed between the MHC-II allele and the selected epitope “LAIDTGCLY” which also exhibits the satisfactory binding affinity. So, the results indicate the epitope “LAIDTGCLY” as a novel vaccine candidate of OROV.

In this study, five linear B-cell epitopes ¹⁷HPLSTSQIGDRC²⁸ (94.44%), ⁷⁸⁸SHCNLEFTAITADKIMSL⁸⁰⁵ (83.33%), ⁹²⁸HHYKPTKNLPHVVPRYH⁹⁴⁴ (100%), ¹⁰²⁶PEKIPAKEGWLTFSKEHTSSW¹⁰⁴⁶ (100%), and ¹³²⁷NQKIDLSQL¹³³⁵ (100%) showed the highest conservancy among diverse OROV isolates compared with other epitopes, proposing that these epitopes might be the ideal vaccine candidates based on the fact that an epitope should be as conservative as possible to provide extensive protection among different virus isolates. These epitopes were also identified as nontoxic to humans. However, the cross processing of the BepiPred-2.0 and LBtope predicted B-cell epitope data revealed that the 17-mer epitope ⁹²⁸HHYKPTKNLPHVVPRYH⁹⁴⁴ is the most potent B-cell epitope. Besides, the discontinuous B-cell epitopes were drawn on the polyprotein of OROV (Figure 5). All the predicted discontinuous or conformational epitopes were placed on the surface of the OROV polyprotein, representing the accessibility for the entered virus.

As recommended by the antigenicity analysis through the VaxiJen server, the polyprotein of OROV seemed to be a suitable vaccine target. The important information was found from the physicochemical properties and secondary structural features of the protein by the ExPASy ProtParam and PSIPRED servers. The protein was found to be basic, stable, and hydrophilic in nature and could sustain high temperature due to a higher aliphatic index value. The abundance of random coils in the OROV polyprotein structure ensured the higher conservancy and stability of the protein [39, 90]. This protein contained disordered regions in some specific positions. The most disordered region was positioned in the 705 to 712 (STRSAQHS) location in the protein. The predicted 3D structure of the OROV polyprotein by the I-TASSER server was acceptable. The energy minimization and evaluation results ensured the good quality of the predicted model. The evaluated model was submitted to the PMDB database for further research in the future. To support the future drug discovery and drug design, ligand-binding pockets and hydrophobicity analyses were accomplished. However, ten possible ligand-binding pockets and their respective amino acid residues were found. The results suggested that the binding residues of these pockets could be the future research prospect for effective inhibitor design against emerging OROV infection.

Peptide’s half-life prediction is one of the key challenges in successful peptide vaccine design as its bioavailability depended on the extent of the half-life. So it is vital to identify the optimal half-life of the peptide for ensuring their optimal function [91] post-vaccination. Recently, a database called “PEPlife” has been developed as a repository of the half-life of peptides [92]; however, it is customary to test and measure for the peptide half-life to ensure optimal protection post-vaccination. Cytokines play an important role in elucidating humoral immune response against antiviral activity, and methods have been developed to predict the different interleukins such as IFN-γ [93], IL-4 [94], and IL-10 [95]. In general, protein-based vaccines have limited immunogenicity and reactogenicity compared to live-attenuated and inactivated vaccines, although protein-based vaccines have greater advantages in relation to safety and cost-effectiveness [96, 97]. An overarching strategy to achieve maximum protection against viral infections is to design a successful peptide-based vaccine following the identification of crucial epitopes via the immunoinformatic approach combined with effective adjuvant choice. Adjuvants are substances added to vaccine formulations that influence the onset, strength, and longevity of specific immune responses to antigens and induce protection against infection [96, 98]. Moreover, adjuvants receive much attention due to their ability to selectively modulate the immune response to elicit humoral and/or cell-mediated immune response [97]. A number of approved vaccine adjuvants such as aluminium salts, oil-in-water (o/w) emulsions, virosomes, liposomes, monophosphoryl lipid A, and TLR agonist are available, but at this time in silico adjuvant prediction tools are not available due to the unavailability of sufficient data and exact molecular mechanism of vaccine adjuvants.

5. Conclusion

This study provided valuable insights for the identification of vaccine candidates and ligand-binding pockets against the neglected and emerging OROV infection through the immunoinformatics as well as Bioinformatics approaches. Based on this strategy, the T-cell and B-cell epitopes in the polyprotein of OROV were mapped and selected as putative OROV vaccine targets. The predicted epitopes exhibited T-cell and B-cell selectivity, higher conservancy, nonallergenicity, nontoxicity, higher population coverage, and significant interaction with MHC class I allele with good affinity. Nonetheless, these findings will provide primary data for molecular docking and epitope-based peptide vaccine research in the future.

Abbreviations

OROV: | Oropouche virus
CNS: | Central nervous system
CTL: | Cytotoxic T-lymphocytes
FAO: | Food and Agriculture Organization
GRAVY: | Grand average hydropathicity
HMM: | Hidden Markov model
IRF-3: | Interferon regulatory transcription factor
MAVS: | Mitochondrial antiviral signaling protein
MHC-I: | Major histocompatibility complex class I
MHC-II: | Major histocompatibility complex class II
PI: | Protrusion index
PMDB: | Protein model database
PDB: | Protein Data Bank
SMM: | Stabilized matrix base method
SVM: | Support vector machine
TAP: | Transport-associated protein
WHO: | World Health Organization.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no competing interests.

Acknowledgments

The authors are thankful to Professor A.K.M Akhtarul Islam, Department of English, Islamic University, Kushtia, for editing and checking the spelling and grammatical error of this manuscript.

Supplementary Materials

Supplementary 1. Figure S1: docking simulation study (control). (A) The docking results as the surface structure of the MHC-I HLA allele “HLA-B35:01” (cyan color) and the sticky form of epitope “LPFERATVM” (red color). The black color box represents the position of epitope in the surface structure. (B) Epitope in the cartoon structure. (C) Sticky form of epitope “LPFERATVM” (red color) with residue position. (D) Binding interaction of the epitope residues and the HLA-B35:01 residues through the hydrogen bonds. Hydrogen bonds and the distances are shown in green color.

Supplementary 2. Figure S2: secondary structure of the OROV polyprotein, ALB07207.

Supplementary 3. Figure S3: Ramachandran plot analysis of OROV polyprotein. Here, the red region indicates the favored region, the yellow region for the allowed region, light yellow for the generously allowed region, and white for the disallowed region. Phi and Psi angles determine torsion angles.

Supplementary 4. Figure S4: ProSA—web results of the OROV polyprotein. (A) ProSA—web z-scores of all protein chains in PDB determined by X-ray crystallography (light blue) or NMR spectroscopy (dark blue) with respect to their length. The z-score of the OROV polyprotein is highlighted as large dots. (B) Energy plot of the OROV polyprotein. Residue energies averaged over a sliding window are plotted as a function of the central residue in the window. (C) Jmol Cα trace of the OROV polyprotein. Residues are colored from blue to red in the order of increasing residue energy.

Supplementary 5. Figure S5: QMEAN score and QMEANDisCo score of the OROV polyprotein. (A) QMEAN4 score with linear combination of four statistical potential terms. (B) QMEAN6 score with two additional agreement terms evaluating the consistency of structural features with sequence-based predictions. (C) Predicted local similarity to target structure. (D) QMEANDisCo scoring plot. It compares local QMEAN scores with local DisCo scores and the resulting QMEANDisCo scores. Depending on the situation of finding homologous, DisCo is not necessarily defined everywhere.

Supplementary 6. Figure S6: hydrophobicity plot analysis of OROV polyprotein, ALB07207.

Supplementary 7. Table S1: total 128 epitopes predicted for 12 super types. (√) indicating the acceptance of the epitope which has both positive antigenicity and immunogenicity scores.

Supplementary 8. Table S2: 37 CD8⁺ T-cell epitopes and their combined score, antigenicity, immunogenicity, conservancy, and MHC-I interaction.

Supplementary 9. Table S3: population coverage for the selected 18 epitopes and their HLA alleles.

References

H. B. Vasconcelos, M. R. T. Nunes, L. M. N. Casseb et al., “Molecular epidemiology of Oropouche virus, Brazil,” Emerging Infectious Diseases, vol. 17, no. 5, pp. 800–806, 2011.

J. F. T. da Rosa, W. M. de Souza, F. de Paula Pinheiro et al., “Oropouche virus: clinical, epidemiological, and molecular aspects of a neglected Orthobunyavirus,” The American Journal of Tropical Medicine and Hygiene, vol. 96, no. 5, pp. 1019–1030, 2017.

F. P. Pinheiro, A. P. A. Travassos da Rosa, J. F. S. Travassos da Rosa et al., “Oropouche virus. I. A review of clinical, epidemiological, and ecological findings,” The American Journal of Tropical Medicine and Hygiene, vol. 30, no. 1, pp. 149–160, 1981.

F. P. Pinheiro, A. G. Rocha, R. B. Freitas et al., “Meningitis associated with Oropouche virus infections,” Revista do Instituto de Medicina Tropical de São Paulo, vol. 24, no. 4, pp. 246–251, 1982.

R. Santos, L. Bueno-Júnior, R. Ruggiero et al., “Spread of Oropouche virus into the central nervous system in mouse,” Viruses, vol. 6, no. 10, pp. 3827–3836, 2014.

J. L. Proenca-Modena, R. Sesti-Costa, A. K. Pinto et al., “Oropouche virus infection and pathogenesis are restricted by MAVS, IRF-3, IRF-7, and type I interferon signaling pathways in nonmyeloid cells,” Journal of Virology, vol. 89, no. 9, pp. 4720–4737, 2015.

C. R. Anderson, L. Spence, W. G. Downs, and T. H. Aitken, “Oropouche virus: a new human disease agent from Trinidad, West Indies,” The American Journal of Tropical Medicine and Hygiene, vol. 10, no. 4, pp. 574–578, 1961.

M. R. T. Nunes, L. C. Martins, S. G. Rodrigues et al., “Oropouche virus isolation, southeast Brazil,” Emerging Infectious Diseases, vol. 11, no. 10, pp. 1610–1613, 2005.

R. M. Elliott, “Orthobunyaviruses: recent genetic and structural insights,” Nature Reviews Microbiology, vol. 12, no. 10, pp. 673–685, 2014.

N. L. Tilston-Lunel, J. Hughes, G. O. Acrani et al., “Genetic analysis of members of the species Oropouche virus and identification of a novel M segment sequence,” Journal of General Virology, vol. 96, no. 7, pp. 1636–1650, 2015.

G. O. Acrani, N. L. Tilston-Lunel, M. Spiegel et al., “Establishment of a minigenome system for Oropouche virus reveals the S genome segment to be significantly longer than reported previously,” Journal of General Virology, vol. 96, pp. 513–523, 2015.

F. Pinheiro, A. Travassos da Rosa, M. Gomes, J. LeDuc, and A. Hoch, “Transmission of Oropouche virus from man to hamster by the midge Culicoides paraensis,” Science, vol. 215, no. 4537, pp. 1251–1253, 1982.

A. L. Hoch, F. P. Pinheiro, D. R. Roberts, and M. L. Gomes, “Laboratory transmission of Oropouche virus by Culex quinquefasciatus Say,” Bulletin of the Pan American Health Organization, vol. 21, no. 1, pp. 55–61, 1987.

M. L. Moreli, V. H. Aquino, A. C. R. Cruz, and L. T. M. Figueiredo, “Diagnosis of Oropouche virus infection by RT-nested-PCR,” Journal of Medical Virology, vol. 66, no. 1, pp. 139–142, 2002.

M. F. Saeed, M. Nunes, P. F. Vasconcelos et al., “Diagnosis of Oropouche virus infection using a recombinant nucleocapsid protein-based enzyme immunoassay,” Journal of Clinical Microbiology, vol. 39, no. 7, pp. 2445–2452, 2001.

M. Weidmann, V. Rudaz, M. R. T. Nunes, P. F. C. Vasconcelos, and F. T. Hufert, “Rapid detection of human pathogenic orthobunyaviruses,” Journal of Clinical Microbiology, vol. 41, no. 7, pp. 3299–3305, 2003.

M. M. Gibrail, F. S. Fiaccadori, M. Souza et al., “Detection of antibodies to Oropouche virus in non-human primates in Goiânia City, Goiás,” Revista da Sociedade Brasileira de Medicina Tropical, vol. 49, no. 3, pp. 357–360, 2016.

B. E. Pickett, E. L. Sadat, Y. Zhang et al., “ViPR: an open bioinformatics database and analysis resource for virology research,” Nucleic Acids Research, vol. 40, no. D1, pp. D593–D598, 2012.

I. A. Doytchinova and D. R. Flower, “VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines,” BMC Bioinformatics, vol. 8, no. 1, p. 4, 2007.

M. V. Larsen, C. Lundegaard, K. Lamberth, S. Buus, O. Lund, and M. Nielsen, “Large-scale validation of methods for cytotoxic T-lymphocyte epitope prediction,” BMC Bioinformatics, vol. 8, no. 1, p. 424, 2007.

M. V. Larsen, C. Lundegaard, K. Lamberth et al., “An integrative approach to CTL epitope prediction: a combined algorithm integrating MHC class I binding, TAP transport efficiency, and proteasomal cleavage predictions,” European Journal of Immunology, vol. 35, no. 8, pp. 2295–2303, 2005.

C. Lundegaard, K. Lamberth, M. Harndahl, S. Buus, O. Lund, and M. Nielsen, “NetMHC-3.0: accurate web accessible predictions of human, mouse and monkey MHC class I affinities for peptides of length 8–11,” Nucleic Acids Research, vol. 36, Supplement 2, pp. W509–W512, 2008.

B. Peters, J. Sidney, P. Bourne et al., “The design and implementation of the immune epitope database and analysis resource,” Immunogenetics, vol. 57, no. 5, pp. 326–336, 2005.

B. Peters, S. Bulik, R. Tampe, P. M. van Endert, and H.-G. Holzhütter, “Identifying MHC class I epitopes by predicting the TAP transport efficiency of epitope precursors,” Journal of Immunology, vol. 171, no. 4, pp. 1741–1749, 2003.

S. Tenzer, B. Peters, S. Bulik et al., “Modeling the MHC class I pathway by combining predictions of proteasomal cleavage, TAP transport and MHC class I binding,” Cellular and Molecular Life Sciences CMLS, vol. 62, no. 9, pp. 1025–1037, 2005.

M. Nielsen, C. Lundegaard, and O. Lund, “Prediction of MHC class II binding affinity using SMM-align, a novel stabilization matrix alignment method,” BMC Bioinformatics, vol. 8, no. 1, p. 238, 2007.

J. J. A. Calis, M. Maybeno, J. A. Greenbaum et al., “Properties of MHC class I presented peptides that enhance immunogenicity,” PLoS Computational Biology, vol. 9, no. 10, article e1003266, 2013.

I. Dimitrov, I. Bangov, D. R. Flower, and I. Doytchinova, “AllerTOP v.2—a server for in silico prediction of allergens,” Journal of Molecular Modeling, vol. 20, no. 6, p. 2278, 2014.

I. Dimitrov, L. Naneva, I. Doytchinova, and I. Bangov, “AllergenFP: allergenicity prediction by descriptor fingerprints,” Bioinformatics, vol. 30, no. 6, pp. 846–851, 2014.

S. Gupta, P. Kapoor, K. Chaudhary et al., “In silico approach for predicting toxicity of peptides and proteins,” PLoS One, vol. 8, no. 9, article e73957, 2013.

P. Thévenet, Y. Shen, J. Maupetit, F. Guyon, P. Derreumaux, and P. Tufféry, “PEP-FOLD: an updated de novo structure prediction server for both linear and disulfide bonded cyclic peptides,” Nucleic Acids Research, vol. 40, no. W1, pp. W288–W293, 2012.

Y. Shen, J. Maupetit, P. Derreumaux, and P. Tufféry, “Improved PEP-FOLD approach for peptide and miniprotein structure prediction,” Journal of Chemical Theory and Computation, vol. 10, no. 10, pp. 4745–4758, 2014.

G. M. Morris, R. Huey, W. Lindstrom et al., “AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility,” Journal of Computational Chemistry, vol. 30, no. 16, pp. 2785–2791, 2009.

O. Trott and A. J. Olson, “AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading,” Journal of Computational Chemistry, vol. 31, no. 2, pp. 455–461, 2010.

S. Gras, L. Kedzierski, S. A. Valkenburg et al., “Cross-reactive CD8⁺ T-cell immunity between the pandemic H1N1-2009 and H1N1-1918 influenza A viruses,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 28, pp. 12599–12604, 2010.

M. C. Jespersen, B. Peters, M. Nielsen, and P. Marcatili, “BepiPred-2.0: improving sequence-based B-cell epitope prediction using conformational epitopes,” Nucleic Acids Research, vol. 45, no. W1, pp. W24–W29, 2017.

H. Singh, H. R. Ansari, and G. P. S. Raghava, “Improved method for linear B-cell epitope prediction using antigen’s primary sequence,” PLoS One, vol. 8, no. 5, article e62216, 2013.

J. Ponomarenko, H.-H. Bui, W. Li et al., “ElliPro: a new structure-based tool for the prediction of antibody epitopes,” BMC Bioinformatics, vol. 9, no. 1, p. 514, 2008.

E. Gasteiger, C. Hoogland, A. Gattiker et al., “Protein identification and analysis tools on the ExPASy server,” in The Proteomics Protocols Handbook, J. M. Walker, Ed., pp. 571–607, Humana Press, Totowa, NJ, USA, 2005.

L. J. McGuffin, K. Bryson, and D. T. Jones, “The PSIPRED protein structure prediction server,” Bioinformatics, vol. 16, no. 4, pp. 404-405, 2000.

D. T. Jones and D. Cozzetto, “DISOPRED3: precise disordered region predictions with annotated protein-binding activity,” Bioinformatics, vol. 31, no. 6, pp. 857–863, 2015.

A. Roy, A. Kucukural, and Y. Zhang, “I-TASSER: a unified platform for automated protein structure and function prediction,” Nature Protocols, vol. 5, no. 4, pp. 725–738, 2010.

Y. Zhang, “I-TASSER server for protein 3D structure prediction,” BMC Bioinformatics, vol. 9, no. 1, p. 40, 2008.

D. Xu and Y. Zhang, “Improving the physical realism and structural accuracy of protein models by a two-step atomic-level energy minimization,” Biophysical Journal, vol. 101, no. 10, pp. 2525–2534, 2011.

N. Guex, M. C. Peitsch, and T. Schwede, “Automated comparative protein structure modeling with SWISS-MODEL and Swiss-PdbViewer: a historical perspective,” Electrophoresis, vol. 30, no. S1, pp. S162–S173, 2009.

M. Wiederstein and M. J. Sippl, “ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins,” Nucleic Acids Research, vol. 35, Supplement 2, pp. W407–W410, 2007.

P. Benkert, M. Künzli, and T. Schwede, “QMEAN server for protein model quality estimation,” Nucleic Acids Research, vol. 37, Supplement 2, pp. W510–W514, 2009.

B. Wallner and A. Elofsson, “Can correct protein models be identified?” Protein Science, vol. 12, no. 5, pp. 1073–1086, 2003.

T. A. P. de Beer, K. Berka, J. M. Thornton, and R. A. Laskowski, “PDBsum additions,” Nucleic Acids Research, vol. 42, no. D1, pp. D292–D296, 2014.

U. K. Adhikari and M. M. Rahman, “Comparative analysis of amino acid composition in the active site of nirk gene encoding copper-containing nitrite reductase (CuNiR) in bacterial spp.,” Computational Biology and Chemistry, vol. 67, pp. 102–113, 2017.

T. Castrignanò, P. D. De Meo, D. Cozzetto, I. G. Talamo, and A. Tramontano, “The PMDB protein model database,” Nucleic Acids Research, vol. 34, Supplement 1, no. 90001, pp. D306–D309, 2006.

M. Hernandez, D. Ghersi, and R. Sanchez, “SITEHOUND-web: a server for ligand binding site identification in protein structures,” Nucleic Acids Research, vol. 37, Supplement 2, pp. W413–W416, 2009.

T. M. McKeever, S. A. Lewis, C. Smith, and R. Hubbard, “Vaccination and allergic disease: a birth cohort study,” American Journal of Public Health, vol. 94, no. 6, pp. 985–989, 2004.

P. Sun, H. Ju, Z. Liu et al., “Bioinformatics resources and tools for conformational B-cell epitope prediction,” Computational and Mathematical Methods in Medicine, vol. 2013, Article ID 943636, 11 pages, 2013.

M. H. V. Van Regenmortel, “Antigenicity and immunogenicity of synthetic peptides,” Biologicals, vol. 29, no. 3-4, pp. 209–213, 2001.

R. A. R. María, C. V. J. Arturo, J. A. Alicia, M. L. G. Paulina, and A. O. Gerardo, “The impact of bioinformatics on vaccine design and development,” in Vaccines, F. Afrin, H. Hemeg, and H. Ozbak, Eds., InTech, Rijeka, Croatia, 2017.

K. L. Seib, X. Zhao, and R. Rappuoli, “Developing vaccines in the era of genomics: a decade of reverse vaccinology,” Clinical Microbiology and Infection, vol. 18, pp. 109–116, 2012.

A. S. De Groot and R. Rappuoli, “Genome-derived vaccines,” Expert Review of Vaccines, vol. 3, no. 1, pp. 59–76, 2004.

B. Korber, M. LaBute, and K. Yusim, “Immunoinformatics comes of age,” PLoS Computational Biology, vol. 2, no. 6, article e71, 2006.

A. W. Purcell, J. McCluskey, and J. Rossjohn, “More than one reason to rethink the use of peptides in vaccine design,” Nature Reviews Drug Discovery, vol. 6, no. 5, pp. 404–414, 2007.

R. Ahmed and D. Gray, “Immunological memory and protective immunity: understanding their relation,” Science, vol. 272, no. 5258, pp. 54–60, 1996.

I. J. Amanna, N. E. Carlson, and M. K. Slifka, “Duration of humoral immunity to common viral and vaccine antigens,” The New England Journal of Medicine, vol. 357, no. 19, pp. 1903–1915, 2007.

M. Slifka and R. Ahmed, “Long-term humoral immunity against viruses: revisiting the issue of plasma cell longevity,” Trends in Microbiology, vol. 4, no. 10, pp. 394–400, 1996.

R. B. Tesh and P. F. C. Vasconcelos, “Sandfly fever, Oropouche fever, and other bunyavirus infections,” in Tropical Infectious Diseases: Principles, Pathogens and Practice, pp. 481-482, Elsevier, Philadelphia, PA, USA, Third edition, 2011.

T. Dorner and A. Radbruch, “Antibodies and B cell memory in viral immunity,” Immunity, vol. 27, no. 3, pp. 384–392, 2007.

R. M. Welsh, L. K. Selin, and E. Szomolanyi-Tsuda, “Immunological memory to viral infections,” Annual Review of Immunology, vol. 22, no. 1, pp. 711–743, 2004.

K.-Y. A. Huang, C. K. F. Li, E. Clutterbuck et al., “Virus-specific antibody secreting cell, memory B-cell, and sero-antibody responses in the human influenza challenge model,” The Journal of Infectious Diseases, vol. 209, no. 9, pp. 1354–1361, 2014.

B. J. Hebeis, K. Klenovsek, P. Rohwer et al., “Activation of virus-specific memory B cells in the absence of T cell help,” Journal of Experimental Medicine, vol. 199, no. 4, pp. 593–602, 2004.

C. Vandermeulen, L. Verhoye, S. Vaidya et al., “Detection of mumps virus-specific memory B cells by transfer of peripheral blood mononuclear cells into immune-deficient mice,” Immunology, vol. 131, no. 1, pp. 33–39, 2010.

H. M. Joo, Y. He, A. Sundararajan, L. Huan, and M. Y. Sangster, “Quantitative analysis of influenza virus-specific B cell memory generated by different routes of inactivated virus vaccination,” Vaccine, vol. 28, no. 10, pp. 2186–2194, 2010.

C. F. Narvaez, N. Feng, C. Vasquez et al., “Human rotavirus-specific IgM memory B cells have differential cloning efficiencies and switch capacities and play a role in antiviral immunity in vivo,” Journal of Virology, vol. 86, no. 19, pp. 10829–10840, 2012.

D. Herrera, O. L. Rojas, C. Duarte-Rey, R. D. Mantilla, J. Ángel, and M. A. Franco, “Simultaneous assessment of rotavirus-specific memory B cells and serological memory after B cell depletion therapy with rituximab,” PLoS One, vol. 9, no. 5, article e97087, 2014.

S. A. Smith, A. R. de Alwis, N. Kose, R. S. Jadi, A. M. de Silva, and J. E. Crowe, “Isolation of dengue virus-specific memory B cells with live virus antigen from human subjects following natural infection reveals the presence of diverse novel functional groups of antibody clones,” Journal of Virology, vol. 88, no. 21, pp. 12233–12241, 2014.

N. L. Tilston-Lunel, G. O. Acrani, R. E. Randall, and R. M. Elliott, “Generation of recombinant Oropouche viruses lacking the nonstructural protein NSm or NSs,” Journal of Virology, vol. 90, no. 5, pp. 2616–2627, 2016.

M. T. Esser, R. D. Marchese, L. S. Kierstead et al., “Memory T cells and vaccines,” Vaccine, vol. 21, no. 5-6, pp. 419–430, 2003.

M. V. Larsen, A. Lelic, R. Parsons et al., “Identification of CD8⁺ T cell epitopes in the West Nile virus polyprotein by reverse-immunology using NetCTL,” PLoS One, vol. 5, no. 9, article e12697, 2010.

A. Alam, S. Ali, S. Ahamad, M. Z. Malik, and R. Ishrat, “From ZikV genome to vaccine: in silico approach for the epitope-based peptide vaccine against Zika virus envelope glycoprotein,” Immunology, vol. 149, no. 4, pp. 386–399, 2016.

S. Chakraborty, R. Chakravorty, M. Ahmed et al., “A computational approach for identification of epitopes in dengue virus envelope protein: a step towards designing a universal dengue vaccine targeting endemic regions,” In Silico Biology, vol. 10, no. 5-6, pp. 235–246, 2010.

M. A. Hasan, M. A. Khan, A. Datta, M. H. H. Mazumder, and M. U. Hossain, “A comprehensive immunoinformatics and target site study revealed the corner-stone toward Chikungunya virus treatment,” Molecular Immunology, vol. 65, no. 1, pp. 189–204, 2015.

U. K. Adhikari and M. M. Rahman, “Overlapping CD8 + and CD4 + T-cell epitopes identification for the progression of epitope-based peptide vaccine from nucleocapsid and glycoprotein of emerging Rift Valley fever virus using immunoinformatics approach,” Infection, Genetics and Evolution, vol. 56, pp. 75–91, 2017.

A. R. Oany, T. Pervin, M. Mia et al., “Vaccinomics approach for designing potential peptide vaccine by targeting Shigella spp. serine protease autotransporter subfamily protein SigA,” Journal of Immunology Research, vol. 2017, Article ID 6412353, 14 pages, 2017.

M. U. Hossain, A. Hashem, C. A. Keya, and M. Salimullah, “Therapeutics insight with inclusive immunopharmacology explication of human rotavirus a for the treatment of diarrhea,” Frontiers in Pharmacology, vol. 7, p. 153, 2016.

R. de Freitas e Silva, L. F. G. R. Ferreira, M. Z. Hernandes et al., “Combination of in silico methods in the search for potential CD4⁺ and CD8⁺ T cell epitopes in the proteome of Leishmania braziliensis,” Frontiers in Immunology, vol. 7, p. 327, 2016.

M. K. Khan, S. Zaman, S. Chakraborty et al., “In silico predicted mycobacterial epitope elicits in vitro T-cell responses,” Molecular Immunology, vol. 61, no. 1, pp. 16–22, 2014.

K. Maenaka and E. Y. Jones, “MHC superfamily structure and the immune system,” Current Opinion in Structural Biology, vol. 9, no. 6, pp. 745–753, 1999.

L. J. Stern and D. C. Wiley, “Antigenic peptide binding by class I and class II histocompatibility proteins,” Structure, vol. 2, no. 4, pp. 245–251, 1994.

M. Tshabalala, J. Mellet, and M. S. Pepper, “Human leukocyte antigen diversity: a southern African perspective,” Journal of Immunology Research, vol. 2015, Article ID 746151, 11 pages, 2015.

K. Cao, A. M. Moormann, K. E. Lyke et al., “Differentiation between African populations is evidenced by the diversity of alleles and haplotypes of HLA class I loci,” Tissue Antigens, vol. 63, no. 4, pp. 293–325, 2004.

M. Paximadis, T. Y. Mathebula, N. L. Gentle et al., “Human leukocyte antigen class I (A, B, C) and II (DRB1) diversity in the black and Caucasian South African population,” Human Immunology, vol. 73, no. 1, pp. 80–92, 2012.

K. Guruprasad, B. V. B. Reddy, and M. W. Pandit, “Correlation between stability of a protein and its dipeptide composition: a novel approach for predicting in vivo stability of a protein from its primary sequence,” Protein Engineering, Design and Selection, vol. 4, no. 2, pp. 155–161, 1990.

H. Jenssen and S. I. Aspmo, “Serum stability of peptides,” in Peptide-Based Drug Design, L. Otvos, Ed., pp. 177–186, Humana Press, Totowa, NJ, USA, 2008.

D. Mathur, S. Prakash, P. Anand et al., “PEPlife: a repository of the half-life of peptides,” Scientific Reports, vol. 6, no. 1, article 36617, 2016.

S. K. Dhanda, P. Vir, and G. P. S. Raghava, “Designing of interferon-gamma inducing MHC class-II binders,” Biology Direct, vol. 8, no. 1, p. 30, 2013.

S. K. Dhanda, S. Gupta, P. Vir, and G. P. S. Raghava, “Prediction of IL4 inducing peptides,” Clinical and Developmental Immunology, vol. 2013, article 263952, 9 pages, 2013.

G. Nagpal, S. S. Usmani, S. K. Dhanda et al., “Computer-aided designing of immunosuppressive peptides based on IL-10 inducing potential,” Scientific Reports, vol. 7, no. 1, article 42851, 2017.

S. Lee and M. T. Nguyen, “Recent advances of vaccine adjuvants for infectious diseases,” Immune Network, vol. 15, no. 2, pp. 51–57, 2015.

S. R. Bonam, C. D. Partidos, S. K. M. Halmuthur, and S. Muller, “An overview of novel adjuvants designed for improving vaccine efficacy,” Trends in Pharmacological Sciences, vol. 38, no. 9, pp. 771–793, 2017.

V. E. J. C. Schijns and E. C. Lavelle, “Trends in vaccine adjuvants,” Expert Review of Vaccines, vol. 10, no. 4, pp. 539–550, 2011.