Next Article in Journal
Neuroprotective Effects of Baicalein, Wogonin, and Oroxylin A on Amyloid Beta-Induced Toxicity via NF-κB/MAPK Pathway Modulation
Next Article in Special Issue
Discovery of Cysteine and Its Derivatives as Novel Antiviral and Antifungal Agents
Previous Article in Journal
Multiple β Forms of Saturated Monoacid Triacylglycerol Crystals
Previous Article in Special Issue
Assessment of Soybean Oil Oxidative Stability from Rapid Analysis of its Minor Component Profile
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 25
Issue 21
10.3390/molecules25215088
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Epitope-Based Immunoinformatics Approach on Nucleocapsid Protein of Severe Acute Respiratory Syndrome-Coronavirus-2

by
Ahmed Rakib
1,
Saad Ahmed Sami
1,
Md. Ashiqul Islam
1,2,
Shahriar Ahmed
1,
Farhana Binta Faiz
1,
Bibi Humayra Khanam
1,
Kay Kay Shain Marma
1,
Maksuda Rahman
1,
Mir Muhammad Nasir Uddin
1,
Firzan Nainu
3,
Talha Bin Emran
4,* and
Jesus Simal-Gandara
5,*
1
Department of Pharmacy, Faculty of Biological Sciences, University of Chittagong, Chittagong 4331, Bangladesh
2
Department of Pharmacy, Mawlana Bhashani Science & Technology University, Santosh, Tangail 1902, Bangladesh
3
Faculty of Pharmacy, Hasanuddin University, Tamalanrea, Kota Makassar, Sulawesi Selatan 90245, Indonesia
4
Department of Pharmacy, BGC Trust University Bangladesh, Chittagong 4381, Bangladesh
5
Nutrition and Bromatology Group, Department of Analytical and Food Chemistry, Faculty of Food Science and Technology, University of Vigo–Ourense Campus, E32004 Ourense, Spain
*
Authors to whom correspondence should be addressed.
Molecules 2020, 25(21), 5088; https://doi.org/10.3390/molecules25215088
Submission received: 19 September 2020 / Revised: 26 October 2020 / Accepted: 28 October 2020 / Published: 2 November 2020
(This article belongs to the Special Issue Food and Drug Analysis Ⅱ)
Browse Figures
Versions Notes

Abstract

:
With an increasing fatality rate, severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) has emerged as a promising threat to human health worldwide. Recently, the World Health Organization (WHO) has announced the infectious disease caused by SARS-CoV-2, which is known as coronavirus disease-2019 (COVID-2019), as a global pandemic. Additionally, the positive cases are still following an upward trend worldwide and as a corollary, there is a need for a potential vaccine to impede the progression of the disease. Lately, it has been documented that the nucleocapsid (N) protein of SARS-CoV-2 is responsible for viral replication and interferes with host immune responses. We comparatively analyzed the sequences of N protein of SARS-CoV-2 for the identification of core attributes and analyzed the ancestry through phylogenetic analysis. Subsequently, we predicted the most immunogenic epitope for the T-cell and B-cell. Importantly, our investigation mainly focused on major histocompatibility complex (MHC) class I potential peptides and NTASWFTAL interacted with most human leukocyte antigen (HLA) that are encoded by MHC class I molecules. Further, molecular docking analysis unveiled that NTASWFTAL possessed a greater affinity towards HLA and also available in a greater range of the population. Our study provides a consolidated base for vaccine design and we hope that this computational analysis will pave the way for designing novel vaccine candidates.

1. Introduction

The present world has witnessed the outbreak of many life-threatening human pathogens including Ebola, Chikungunya, Zika, severe acute respiratory syndrome coronavirus (SARS-CoV), and Middle East respiratory syndrome coronavirus (MERS-CoV) in the 21st century. More recently in late December 2019, a cluster of pneumonia cases was reported in the city of Wuhan, Hubei province, China, which was of unknown cause. Later it was confirmed that these pneumonia cases were due to a novel coronavirus named SARS-CoV-2 (previously named as 2019-nCoV) and the disease condition of this virus is referred to as COVID-19 [1,2,3]. On 11 March, 2020, the World Health Organization (WHO) assessed that COVID-19 can be characterized as a pandemic. The current COVID-19 pandemic is a global concern and is spreading at an alarming rate and as of 26 October, 2020, more than 43.2 million cases along with over 1.16 million deaths have been reported globally [4].
As COVID-19 is mainly a respiratory disease, in most cases it might affect the lungs only. The primary mode of infection is human-to-human transmission through close contact, which occurs via spraying droplets from the infected individual through their cough or sneeze. The symptoms of this coronavirus can be mild to moderate or severe including, fever, cough, and shortness of breath or pneumonia. Respiratory, hepatic and neurological complications can be seen in case of severe cases that can lead to death. It seems that the severity and fatality rate of COVID-19 is milder than that of SARS and MERS. Although diarrhea was presented in about 20–25% of patients with SARS and MERS, intestinal symptoms were rarely reported in patients with COVID-19 [5,6,7]. Multi-organ failure, especially in elderly people and people with underlying health conditions, such as hypertension, cardiovascular disease and diabetes, are exhibiting a higher mortality rate in COVID-19.
Interestingly, SARS-CoV-2 has 82% similarity with the original SARS-CoV virus attributed to the outbreak in 2003 [8]. A mature SARS-CoV-2 virus generally has a polyprotein (the open reading frame 1a and 1b, Orf1ab), four structural proteins such as envelope (E) protein; membrane (M) protein; nucleocapsid (N) protein; spike (S) protein and five accessory proteins (Orf3a, Orf6, Orf7a, Orf8 and Orf10), and, particularly, SARS-CoV-2 encodes an additional glycoprotein having acetyl esterase and hemagglutination (HE) attributes, which identified it distinct to its two predecessors [9]. The functions of accessory proteins may include signal inhibition, apoptosis induction and cell cycle arrest [10]. The S protein on the surface of the viral particle enables the infection of host cells by binding to the host cell receptor angiotensin-converting enzyme 2 (ACE2), utilizing the S-protein’s receptor-binding domain (RBD).
The N protein binds to the RNA genome of the COVID-19 and creates a shell or capsid around the enclosed nucleic acid. The N protein is involved in viral RNA synthesis and folding, which interacts with the viral membrane protein during viral assembly and affects host cell responses including cell cycle and translation. An epitope-based peptide vaccine has been raised in this aspect. The core mechanism of the peptide vaccine is based on the chemical method to synthesize the recognized B-cell and T-cell epitopes that can induce specific immune responses and are immune-dominant. T-cell epitopes are short peptide fragments (8–20 amino acids) while the B-cell epitopes can be proteins [11,12].
Once a mutated virus infects the host cells by escaping the antibodies, it then relies upon the T-cell mediated immunity to fight against the virus. Viral proteins are processed into short peptides inside the infected cells and then loaded onto major histocompatibility complexes (MHC) proteins. After that, the MHC-peptide complexes are presented on the infected cell surface for recognition by specific T-cells. Activated CD8+ T-cells then recognize the infected cells and clear them. T-cell immunity also depends strictly on the MHC-peptide complexes, which are similar to the antigen-antibody association. MHC proteins are encoded by human leukocyte antigen (HLA), which is located among the most genetically variable regions on the human genome. Each HLA allele can only present a certain set of peptides that can be presented on the infected cell surface and recognized by T-cells are called T-cell epitopes. For a vaccine, it is essential to identify T-cell epitopes that originate from conserved regions of the virus T cell responses against the S and N proteins have been reported to be the most dominant and long-lasting [13].
To develop effective diagnostic tests and vaccine, the identification of B-cell and T-cell epitopes for SARS-CoV-2 proteins are critical especially for structural N and S proteins. Both humoral immunity and cellular immunity provided by B-cell antibodies and T-cells respectively are essential for effective vaccines [14,15]. Although humans may mount an antibody response against viruses normally, only neutralizing antibodies can block the entry of viruses into human cells completely [16]. Antibody binding site’s location on a viral protein strongly affects the body’s ability to produce neutralizing antibodies [17]. It is important to understand whether SARS-CoV-2 has potential antibody binding sites (B-cell epitopes) near their interacting surface with its known human entry receptor, ACE2. Besides neutralizing antibodies, human bodies also depend on cytotoxic CD8+ T-cells and helper CD4+ T-cells to clear viruses completely from the body. For antiviral T-cell responses, presentation of viral peptides by human MHC class I and class II is essential [18]. MHC-I analysis includes common alleles for HLA-A, HLA-B and HLA-C. Multiple investigations have indicated that antibodies generated against the N protein of SARS-CoV are a highly immunogenic and abundantly expressed protein during infection [19].
Our group is targeting for immunoinformatics-based vaccine design using bioinformatics and immunoinformatics tools by utilizing different protein sequences of SARS-CoV-2. Recently, we have already established potential B and T-cell epitopes with a greater candidacy profile using the S protein of SARS-CoV-2 [20]. Moreover, other published work also utilized the S protein of SARS-CoV-2 for epitope-based vaccine design [21]. The purpose of our present study is to promote the designing of a vaccine against COVID-19 using in silico methods, considering SARS-CoV-2 N protein. The reason for focusing particularly on the epitopes in the N structural proteins is due to their dominant and long-lasting immune response, which was reported against SARS-CoV previously [22]. Besides, it has been reported that the N protein of many viruses are highly conserved and immunogenic, which expressed extensively in the course of infection [23]. Particularly, it has been reported recently that the N protein and E protein of SARS-CoV-2 are most evolutionarily conserved [24,25]. For the identified T-cell epitopes, we incorporated the information on the associated MHC alleles so that we can provide a list of epitopes that seek to maximize population coverage globally. Therefore, we designed an epitope-based peptide vaccine through utilizing the SARS-CoV-2 N protein (Figure 1) to potentially narrow down the search for potent targets against SARS-CoV-2 using the computational approach with an expectation that the wet laboratory research will validate our result.

2. Results

2.1. Sequence Retrieval and Analysis

We retrieved the SARS-CoV-2 N protein sequence from the NCBI database (Accession No.: QIC53221.1). Then we performed BLASTp using NCBI-BLAST for the N protein of SARS-CoV-2. We searched for a total of 100 homologs with >60% identical sequences. Multiple sequence alignment (MSA) was then performed to find out the conservancy among the target proteins. (Supplementary Data 1), and a phylogenetic tree was constructed to analyze the evolutionary divergence amongst them (Figure S1). From the results of the MSA analysis, it has been confirmed that the protein sequences have a close relationship.

2.2. Antigenic Protein Prediction

The most potent antigenic protein of SARS-CoV-2 N protein was predicted by VaxiJen v2.0, which is based on the auto-cross covariance transformation of protein sequences into uniform vectors of principal amino acid properties. The VaxiJen tool mainly encompasses the physicochemical properties of the protein sequence [26]. The overall antigen prediction score was 0.5002 (probable antigen) at a 0.4 threshold value.

2.3. Toxicity Prediction

Prediction of the toxicity of peptides before considering them, as epitopes are very important for saving both time and to make it cost effective. The toxicity of the selected peptide sequences was assessed using the ToxinPred web server. ToxinPred is a unique tool, which is based on support vector machine (SVM) in predicting toxicity of peptides and several physicochemical properties, including hydrophilicity, hydrophobicity, charge and molecular weight. The results from the ToxinPred tool showed that all of our probable epitopes were found non-toxic (Table 1).

2.4. Protein Structure Prediction and Validation

The secondary structure of the SARS-CoV-2 N protein was predicted using the self-optimized prediction method with alignment (SOPMA), an online server, During prediction, the SOPMA server can be able to locate almost all of the stretches with the regular structure, which investigate the recognition of folding pattern in an efficient way [27]. The secondary structure of a protein describes mainly the α-helix, β-sheets and random coil. SARS-CoV-2 N protein has 419 residues (Figure 2A), of which 89 residues were remained in the α-helix, 70 residues were from the extended strand, 29 residues were observed in the β-sheets, and 219 residues were remained as random coil (Figure 2B,C). For 3D structure, we built a model using the Robetta online server. The Robetta server predicts the tertiary structure of a given protein from the inputted genomic data. The Robetta server utilizes a fully automated implementation of the Rosetta software package for the inference of the structural information of the protein [28]. In the current experiment, the Robetta server predicted five models for the SARS-CoV-2 N protein, which were validated using PROCHECK and PROSA-Z score. From the result of the validation, it has been observed that Model 4 predicted by the Robetta server have possessed 88.4% amino acid residues in the Rama favored region and delineated Z-score of −7.24, which depicted the model as a good quality model (Figure 2D,E). Although the Z-score for model 1 was shown −7.42, it possessed less amino acid residues in the Rama favored region (Figure S2). In addition, we analyzed the Ramachandran plot statistics and Z-score for the crystal structure of SARS-Cov-2 N protein (Resolution: 2.70 Å). The results showed that the Rama favored region for the crystal structure of SARS-CoV-2 N protein was 88.1% and Z-score was −5.06, which was less compared to the model structure (Figure S2). Hence, model 4 could be used for further analysis.

2.5. CD8+ T-Cell Epitope Identification

The NetCTL 1.2 server was utilized for the prediction of T-cell epitopes. The number of T-cell epitopes depended on the length of the sequence. Further the predicted epitopes with strong binding affinities were subjected to several immune filters in order to screen out the best possible epitopes, including conservation among the protein sequences included in the study, should be immunogenic, should be non-allergic and importantly should not overlap with any human proteins. Based on high combinatorial and MHC binding, the top eight epitopes were predicted by the NetCTL server from the selected protein sequence that was selected for further analysis. Using the MHC-I binding prediction tool, which is based on stabilized matrix method (SMM), we selected those MHC-I alleles for which the epitopes showed the highest affinity (half maximal inhibitory concentration, IC50 < 200 nm).
Proteasomes play an important role in cleaving the peptide bond, resulting in the conversion of protein into the peptide. The peptide molecules that are homogeneous to class I MHC molecules and the peptide-MHC molecule after the proteasomal cleavage were presented as T-helper cells after the transportation into the cell membrane. The total score of each epitope–HLA interaction was taken into consideration and higher processing efficiency was meant by obtaining a higher score. The epitope NTASWFTAL interacted with most of the MHC-I alleles including, HLA-A*68:02, HLA-C*16:01, HLA-C*03:03, HLA-C*03:04, HLA-C*12:03, HLA-A*02:06, HLA-C*03:02, HLA-A*26:01 and HLA-C*14:02 (Table 2). Moreover, the MHC-NP prediction tool was used to find the highest probable score of our predicted epitope NTASWFTAL, with a score of 1.11, for HLA-A*68:02. Furthermore, all the predicted epitopes had a maximum identity for conservancy hit and 100% maximum identity was found (Table 2). Additionally, the I-pMHC immunogenicity prediction analysis of the epitope NTASWFTAL was found 0.22775 (Table 2).

2.6. Population Coverage

Population coverage analysis is crucial in determining a peptide sequence as vaccine candidates. Accordingly, epitope-based vaccines can be designed to maximize the population coverage and minimizing the complexity regarding the variability of the population coverage observed in different ethnic groups. In the current study, the cumulative amount of the population coverage was obtained for the predicted epitope NTASWFTAL. Results from the population coverage demonstrated that with 57.16% coverage, East Asia found the highest coverage region. The results of the population coverage were shown in Table 3 and Figures S3–S6.

2.7. Allergenicity Assessment

The AllerTop server was used for the identification of the allergic reaction caused by a vaccine in an individual that might be harmful or life-threatening. The AllerTop server predicts allergenicity based on several factors, including, amino acid descriptors, accounting for residue hydrophobicity, size, abundance, helix- and β-strand forming propensities and a machine learning approach, namely the k nearest neighbors (kNN) method was implemented to classify allergens and non-allergens [29]. The allergenicity of the selected epitope was calculated using the AllerTop tool and predicted as a probable non-allergen.

2.8. Molecular Docking Analysis for HLA and Epitope Interaction

Molecular docking analysis is used for the prediction of a ligand–receptor interaction. The advancement in computational biology techniques in the last few decades have allowed for further development in molecular docking algorithms for determining the flexibility of a protein and currently, molecular docking is considered as widespread tools used in computational biology techniques. In this study, the verification of the interaction between the HLA molecules and our predicted potential epitope was done by molecular docking simulation using AutoDock Vina in PyRx 0.8 software. Among all the MHC class I alleles, only HLA-A*68:02 had a maximum probable score for our most potent epitope NTASWFTAL. Therefore, we carried out the molecular docking study using HLA-A*68:02 (PDB ID: 4I48). The 3D structure of the predicted epitope, NTASWFTAL and HLA-A*68:02 molecules are represented in Figure 3.
We found that our predicted epitope NTASWFTAL interacted with HLA-A*68:02 with strong binding affinities of -9.4 kcal/mol (Table 4). The selected epitope interacted with Arg6, Ser4, Ser2 and Asp30 residues of chain-A and Lys59, Asp60, Ser58 and Gly30 of chain-B through hydrogen bonding (H-bond), whereas Lys7 residue of chain-B form bonds as a result of sharing electrons (which may happen as a result of charge distribution; Figure 4). Further, for the validation of the docking study, we performed molecular docking analysis between HLA-A*68:02 and the 9-mer peptide bound with the crystal structure of HLA-A*68:02, where the peptide was considered as a positive control. Conversely, the molecular docking analysis between the positive control and HLA-A*68:02 showed less binding affinities than the predicted epitope, where the positive control exhibited a docking score of −8.2 kcal/mol (Table 4). Although the positive control formed six hydrogen bonds, the formed hydrogen bond was less than NTASWFTAL (Figure 5). In addition, a salt bridge was formed between the positive control and Asp29 residue from A chain of HLA-A*68:02.

2.9. B-Cell Epitope Prediction

B-cell epitopes play an important role in the development of epitope-based vaccine and allergic research. A dominant linear B-cell epitope can be used in the autoimmune diseases as the target of neutralizing antibody responses [30]. In addition, they are able to induce an antibody that cross reacts with the parent protein. In this study, using the amino acid scale-based method, we predicted the B-cell epitope identification. Different analysis methods were used for the prediction of the continuous B-cell epitope. The results of the B-cell predictions were shown in Table 5, Table 6 and Table 7, Tables S1 and S2 and Figure 6 and Figure 7.
Firstly, BepiPred linear epitope prediction was used, which is regarded as the best single method for predicting linear B-cell epitopes using a Hidden Markov model. The findings from the BepiPred linear epitope prediction showed maximum score of 2.416 and a minimum score of −0.001, where the average scores were displayed as 0.813 (Table S1).
The β-turns were predicted by the Chaus and Fasman β-turn prediction method. The maximum score was found for the amino residues 2–8 (Figure 6) and the minimum score was attributed for amino acid residues 218–224 (Figure 6).
For antigenicity prediction, the Kolaskar and Tongaonkar antigenicity prediction methods were implied. The method evaluates the antigenicity based on the physicochemical properties of amino acids and their abundances in experimentally known epitopes. The average antigenic propensity of our SARS-CoV-2 N protein was 0.988 with a maximum of 1.197 and a minimum of 0.874 (Figure 7). In addition, the average flexibility of 1.035 and a minimum of 0.874 were predicted by the Karplus and Schulz flexibility prediction method. The residues from 238 to 244 were found to be the most flexible with the highest score of 1.161. The Parker hydrophilicity prediction tool predicts the hydrophilicity of the SARS-CoV-2 N protein with an average score of 2.80, a minimum of 0.874 and the region from amino acid residues 77–83 have shown the maximum score, where the maximum value was 7.006 (Figure 7).
For predicting the surface ability, this study included the Emini surface accessibility prediction method. The average surface accessibility was 1.0 and a minimum 0.050 (Figure 6).

3. Discussion

As of yet, it has been reported that the reproduction rate of SARS-CoV-2 is greater than SARS and MERS and the symptoms of the COVID-19 infection include fever with more than 38 °C body temperature along with alveolar edema, leading to difficulty in breathing, whereas mild symptoms perhaps not engender a high fever [31]. Surprisingly, with a high fatality rate, the severity of the infection was found to be more than the infection caused by both SARS and MERS, with multiple organ damage, which was reported not long ago [32].
At present, researchers are examining repurposed compounds from other viral infections to treat SARS-CoV-2. For example, both lopinavir and ritonavir are HIV protease inhibitors but in a lopinavir–ritonavir clinical trial report, the treatment benefit derived was dubious [33]. From recovering patients, several convalescent immunoglobulins are derived, which is currently investigated as a potential treatment for the disease [34]. As there have been no approved treatments for COVID-19 that exists until now, but remdesivir has been used in some emergency cases and evidence also showed that convalescent plasma could be used as treatment without severe adverse effects [34,35]. These treatments are the best hope for striving to keep the mortality rate low before vaccines become widely available.
Despite many potential challenges, vaccine development is a crucial factor in modern biotechnology as vaccines are the most important prerequisites for defending the burden of diseases over the world [36].
With the divulgement of sequence-based technology in genomics and proteomics, enough pieces of information are available regarding different eukaryotic and prokaryotic organisms including viruses. Therefore, utilizing various bioinformatics tools, it is possible to design peptide-based vaccines through comprehensibly studying the epitopes and several studies suggested epitope-based vaccines against different diseases including dengue, chikungunya, Saint Louis encephalitis virus [37,38,39]. Although epitope-dependent vaccine design is quite familiar, little research works are done in the case of SARS-CoV-2. Being an RNA virus, SARS-CoV-2 is different from the DNA virus and the rate of mutation is higher than the DNA viruses and according to various research, it can be assumed that the mutations might occur in the N protein [40]. Recently, N proteins of SARS-CoV-2 are regarded as a primary target for vaccine development as its function includes viral replication and directly associated with the infection process, as a consequence related to the pathogenesis of COVID-19 [41]. Previous research works have already established that N proteins of several viruses and SARS are considered as a potential target for the development of vaccines [42,43,44,45]. Moreover, we already mentioned the detrimental role of SARS-CoV-2 in host–cell responses. This aspect led us to conduct in silico experiments for designing a peptide-based vaccine against the novel SARS-CoV-2.
Earlier, it has been thought that vaccine development primarily relies on B-cell immunity, but recent discovery unveiled that T-cell epitopes are more propitious as a result of a more long-lasting immune response mediated by CD8+ T-cells and due to the antigenic drift, by which an antibody is not able to respond against an antibody [46]. In this study, focusing on MHC class I potential peptide epitopes, we predicted T-cell and B-cell epitopes, which were able to show immune responses in various ways. Many characteristics including antigenicity, toxicity need to take into consideration for identifying a protein sequence-based epitope into a vaccine candidate and the predicted eight epitopes fulfilled the entire criterion. Toxicity analysis is regarded as an important parameter during design of a peptide sequence into a vaccine candidate. For instance, melittin, a major peptide of bee venom, is a promising candidate for cancer therapy, but due to its toxicity, its applicability has met with critical challenges [47]. In the current study, only five potent epitopes have been predicted from the NetCTL 1.2 server and the epitopes were further taken for the progressive analysis. Besides, all peptides except SSPDDQIGY were able to interact with the MHC class I alleles, and NTASWFTAL interacted with the most MHC class I alleles. Amongst them, HLA-A*68:02 possessed the highest probable score. Further, the conservancy of the epitopes, which was predicted by the IEDB conservancy analysis tool delineated that all of our predicted epitopes had the maximum identity of 100%. Apart from this, a computational study unraveled that the targeted epitope NTASWFTAL showed conservancy along with several epitopes from SARS-CoV-2 [48]. Previously, NTASWFTAL has been used in order to determine the ability to elicit the SARS-CoV immune response [49]. Furthermore, a previous study has already demonstrated that NTASWFTAL interacted with most of the HLA supertypes, including, HLA-A*01:01, HLA-A*02:01, HLA-A*03:01, HLA-A*24:02, HLA-A*26:01, HLA-B*07:02, HLA-B*08:01, HLA-B*27:05, HLA-B*39:01, HLA-B*40:01, HLA-B*58:01 and HLA-B*15:01 [50]. The amino sequence GLPNNTASWFTALTQHGK of SARS-COV-2 N protein also demonstrated the characteristics of the B-cell epitope, which includes the targeted epitope NTASWFTAL [51]. Therefore, we took the epitope NTASWFTAL for further analysis due to its maximum interaction with MHC class I alleles and the highest conservancy.
Generally, allergy is considered as an overreaction of the immune system to a previously captured, harmless, normal protein in nature. True allergic reactions to vaccines are rare; however, their identification is crucial because they can be detrimental to the body [52]. Occasionally, the vaccine itself causes hypersensitivity due to the toxoids present in it. Hence, allergenicity is regarded as one of the most noteworthy obstacles in vaccine development. Importantly, T-cells not CD4+ T-cells are involved in an allergic reaction and an allergic reaction is stimulated by type 2 T helper cell along with immunoglobulin E [53]. In this experiment, we assessed the allergenicity using AllerTop 2.0, which is well recognized for its high sensitivity, and able to identify structurally diverse allergens in comparison with the known allergens. AllerTop predicted our selected epitope as non-allergen.
It has been proposed that the T-cell epitopes bind with the MHC molecules and MHC class I molecules generally presented short peptides that are 8–11 amino acid long, whereas MHC class II molecules present longer peptides with 13–17 amino acid residues [54]. In this experiment, we determined the binding (presence of the antigen on the surface) affinity of the predicted epitope using molecular docking analysis and demonstrated that NTASWFTAL interacted with HLA-A*68:02 and found a binding affinity of −9.4 kcal/mol, which depicted a greater interaction with the epitope and the HLA molecule as the more negative energy implied to more binding affinity [55]. In addition, our predicted epitope delineated greater binding affinities to HLA-A*68:02 than its native ligand. Importantly, a study from Zhang reported the highest binding affinity of NTASWFTAL towards the HLA-A2/A0201-restricted T-cell epitopes [56]. The results from the molecular docking studies in the current study also revealed that epitope NTASWFTAL formed H-bond with both chain-A and chain-B of the HLA molecule and attractive charges were also responsible for the binding.
Another factor that is considered as the most prominent one during the process of vaccine development is population coverage, as the distribution of HLA varies according to ethnicity and geographical region. Although after implementation of several clinical studies, genetic variability on a global scale could have an effect on the significant application of the vaccine candidates in humans [57]. Our experiment showed that the epitope NTASWFTAL covered almost all regions of the world, where the highest coverage was observed in East Asia, where COVID-19 was first reported. Interestingly, our findings indicated that our predicted epitope specifically binds with the widespread HLA molecules and the vaccine will be easily employed.
Importantly, the accurate prediction of T-cell epitopes along with B-cell epitopes is a crucial challenge for the immunoinformatics study and it is unlikely that different HLAs are expressed at different frequencies amongst the ethnic groups. However, substantial research in several in silico markers including the matrix-based profile, and regular expressions in the immunoinformatics study provide a cogent way for prediction of several immunobiological phenomena, for instance, the prediction of subcellular localization (SCL) of a protein is identified by several computational tools. Similarly, T-cell epitope identification has been undergone by implying numerous computational methods and various research areas, including cancer therapy and other infections, T-cell epitope identification is presently apparent [58,59,60]. Additionally, experimental methods established for the calculation of the binding interaction between MHC molecules and an antigenic protein is complicated and time-consuming. Hence, several computational tools have been introduced concerning simulation of the experimental methods, and the methods of the MHC binder prediction are based on motifs, quantitative matrices (QMs), ab initio prediction, machine-learning techniques, DiscoTope, etc. [61]. Several algorithms including PePSSI (peptide–MHC prediction of structure through solvated interfaces) and PREDEP (prediction of MHC class I epitopes) are implemented for the structural prediction and side-chain orientation of the binding proteins. In the current study, the prediction of MHC-I binding with T-cell antigenic peptides from the SARS-CoV-2 N protein sequences was done through the SMM algorithm, which incorporates proteasomal cleavage, TAP transport and MHC class I affinity into the final output and recent studies suggested that SMM is more established than other algorithms such as EpiJen and MAPPP [62,63,64].
Recently, other research works have suggested vaccine design from antigenic protein sequences of SARS-CoV-2 through utilizing in silico immunoinformatics-based methodologies. A study from Lee et al. reported a comprehensive list of antigenic peptides for vaccine development against SARS-CoV-2 [65]. However, the findings of the research work represented that the N protein patterns retained from SARS-CoV-2 were unable to interact with HLA alleles. Several other studies also delineated the high binding affinity of predicted epitopes towards HLA-A*24:02 and HLA-A*02:01 alleles respectively [66,67]. Conversely, in the current research work, our predicted epitope NTASWFTAL exhibited greater affinity towards HLA-A*68:02, predicted by NetCTL 1.2 server. Besides, molecular docking simulation unveiled the greater interaction between the predicted epitope and HLA-A*68:02 molecules. Moreover, our current study is in alignment with previous research work, which depicted peptide-based sequence against the S protein of the human coronavirus [36]. However, we cannot rule out the role of MHC class II peptides during the design of the epitope-based vaccine, as it plays a phenomenal role in humoral immunity through helping B-cells.
In addition, the B-cell epitope provides a strong immune response without causing any adverse effects. Generally, B-cell epitopes are either linear (continuous) or conformational (non- continuous) [68]. Importantly, flexible regions are observed in several crucial parts of a protein, including binding sites, catalytic sites, proteolytic cleavage susceptible sites, allosteric sites and most importantly the antigenic part of a protein sequence. Flexibility analysis is one of the major concerns for the identification of the surface residues forming a protein, which is further demonstrated as potential continuous epitopes [69]. For vaccine development, it would be crucial for predicting the antigenic region. In addition, hydrophilic amino acid residues are major determinants of antigenic features of a protein sequence, as the point highest hydrophilicity is located in or adjacent to an antigenic portion of the protein [70]. In this experiment, we also calculated the linear B-cell epitope prediction. It has been documented that peptide vaccines that are able to demonstrate immune responses against foreign particles contain peptides that are comprised of linear B cell epitopes [71]. B cell epitopes carry specific antigens that bind to the B lymphocytes, as a result they are recognized as potential antigenic determinants and are crucial for vaccine design [72]. In addition, B cell epitopes elicited a stronger immune response, but no side effects were observed. Recently, Grifoni et al. predicted B cell epitopes by utilizing the structural proteins of SARS-CoV and SARS-CoV-2 [73]. The Grifoni study predicted the identity of three peptide sequences from 42–62, 153–172 and 355–401 amino acid residues having an identity ≥ 90% [73]. In the current experiment, by using several tools from IEDB database, we predicted several B-cell epitopes from the SARS-CoV-2 N protein. As a consequence, our study predicted several B cell epitopes that were in line with those identified by Grifoni et al. (Table S2). Additionally, one of the predicted B-cell epitope from amino acid residues 154–166 was in agreement with the study from Amrun et al. (Table S2) [74]. Moreover, several studies have reported the characterization of B-cell epitopes from the N protein of many viruses from humans and animals [23,75,76,77].
Recently, immunoinformatics-aided vaccine design has received experimental validation, which targeted multi-epitope protein clusters from Mycobacterium tuberculosis that interacted with HLA class I and II molecules and their prediction was experimentally validated through in vitro studies [78]. On the other hand, our study was more specific than some similar studies, for example, a study from Khan et al. had selected MHC-I alleles for which the epitopes representing higher affinity (IC50 < 500 nm), but in our study, we showed that epitopes for MHC I alleles showing higher affinity (IC50 < 200 nm), as peptides with minimum IC50 values, exhibited greater inhibition [79,80]. In addition, we assessed immunogenicity, allergenicity and toxicity of the selected epitopes. Moreover, B-cell epitopes can pave the way for experimental epitope mapping and also crucial concerning the interpretation of results from several experiments, including ELISA, radioimmunoassay and Western blotting.
Of course, we understood that this research work does not claim to be exhaustive and all-inclusive as it is true that in silico works have its advantages and limitations. However, recently immunoinformatics is regarded as a new branch of computational biology techniques and is effective in the quest of new immunotherapeutics, amalgamating bioinformatics techniques to figure out several unique problems of vaccinology and immunology [81]. Epitope prediction can be regarded as a high parameter in immunoinformatics investigation, and immunoinformatics calculations are considered as the high frontier to develop effective vaccines true of the practical value. However, the experimental validations of the underlying approaches are required to establish a predicted epitope into a vaccine candidate. The accuracy of the predicted computational analysis should be corroborated by more accessible and robust laboratory experiments.

4. Materials and Methods

4.1. Protein Sequence Retrieval

The SARS-CoV-2 N protein sequence was extracted from the NCBI (National Center for Biotechnology Information) (Bethesda, MD, USA) protein database (Accession no.: QIC53221.1, GI: 1811294683) in the FASTA format.

4.2. Sequence Analysis

The understanding of the features, function, structure and evaluation is mainly based on the process of sequence analysis, which depicts the process of subjecting DNA, RNA or peptide sequences to wide ranges of analytical methods. We employed NCBI BLAST (Basic Local Alignment Search Tool) [82] that screens homologous sequences from its database and selects those sequences that are more similar to our SARS-CoV-2 N protein; we also performed multiple sequence alignment (MSA) using the ClustalW (Conway Institute, UCD, Dublin, Ireland) web server with default settings, and a phylogenetic tree was assembled using MEGA6 software [82,83,84].

4.3. Protein Antigenicity and Toxicity Prediction

To determine the potent antigenic protein of the SARS-CoV-2 N protein, we used the online server VaxiJen v2.0, with a default threshold value [85]. All the antigenic proteins of SARS-CoV-2 N protein with their respective scores were obtained then sorted in Notepad++. A single antigenic protein with maximum antigenicity scores was selected for further evaluation. The toxicity of epitopes was analyzed using the ToxinPred web server [86].

4.4. Protein Secondary and Tertiary Structure Prediction

The secondary structure of the SARS-CoV-2 N protein was predicted by using the SOPMA tool (Institute of Biology and Protein Chemistry, Lyon, France), which correctly predicts 69.5% of amino acids for a three-state description of the secondary structure (α-helix, β-sheet and coil) in a whole database [27]. Additionally, we predicted the 3D structure of the protein using Robetta (University of Washington, Seattle, WA, USA) server, which provides automated tools for prediction and analysis of the tertiary structure of the protein [28]. The model was validated using PROCHECK and PROSA web servers [87,88]. In addition, the 3D crystal structure of SARS-CoV-2 N protein (PDB ID: 6M3M) was downloaded from the Protein Data Bank (PDB) database for comparing the modeled 3D structure of the SARS-CoV-2 N protein.

4.5. T-Cell Epitope Prediction

CD8+ T-Cell Epitope Prediction

For the de novo prediction of the T-cell epitope, NetCTL 1.2 server (DTU Health Tech, Kongens Lyngby, Denmark) was used in this experiment, using a 0.95 threshold to maintain the sensitivity and specificity of 0.90 and 0.95, respectively. The tool expands the prediction for 12 MHC-I supertypes and integrates the prediction of peptide MHC-I binding and proteasomal C-terminal cleavage with TAP transport efficiency. These predictions were performed by an artificial neural network, weighted TAP transport efficiency matrix and a combined algorithm for MHC-I binding and proteasomal cleavage efficiency was then used to determine the overall scores and translated into sensitivity/specificity. Based on this overall score, five best peptides (epitopes) were selected for further evaluation.
For the prediction of peptides binding to MHC-I, we used a tool from the Immune Epitope Database (IEDB) (National Institute of Allergy and Infectious Diseases, Bethesda, MD, USA) and calculate IC50 values for peptides binding to specific MHC-I molecules [89]. For the binding analysis, all the frequently used alleles were selected with a word length of nine residues and binding affinity <200 nm for further analysis. Another tool (named as MHC-NP) provided by the IEDB server was used to assess the probability that a given peptide was naturally processed and bound to a given MHC molecule [90].

4.6. Epitope Conservancy and Immunogenicity Prediction

The degree of similarity between the epitope and the target (i.e., given) sequence was elucidated by epitope conservancy. This property of the epitope gave us the promise of its availability in a range of different strains. Hence for the analysis of the epitope conservancy, the web-based tool from IEDB analysis resources was used [91]. Immunogenicity prediction can uncover the degree of influence (or efficiency) of the respective epitope to produce an immunogenic response. The T-cell class I pMHC immunogenicity predictor at IEDB, which uses amino acid properties as well as their position within the peptide to predict the immunogenicity of a class I peptide MHC (pMHC) complex [92].

4.7. Prediction of Population Coverage and Allergenicity Assessment

The population coverage tool from IEDB was applied to determine the population coverage for every single epitope by selecting HLA alleles of the corresponding epitope.
Allergenicity of the predicted epitope was calculated using AllerTop v2.0 (Medical University, Sofia, Bulgaria) [29], which is an alignment-free server, used for in silico based allergenicity prediction of a protein-based on its physiochemical properties.

4.8. HLA and Epitope Interaction Analysis Using Molecular Docking Studies

4.8.1. Epitope Model Generation

The 3D structures of the selected epitopes were predicted by PEP-FOLD, a web-based server [93]. For each sequence, the server predicted five probable structures. The energy of each structure was determined by SWISS-PDB VIEWER and the structure with the lowest energy was chosen for further analysis [94].

4.8.2. Retrieval of the HLA Allele Molecule

The three-dimensional structure of the HLA-A*68:02 (PDB ID: 4I48) was retrieved from Protein Data Bank (RCSB-PDB).

4.8.3. Molecular Docking Analysis

Molecular docking analysis was performed using AutoDock vina (Scripps Research, La Jolla, CA, USA) in PyRx 0.8, by considering the HLA-A*68:02 molecule as the receptor protein and identified epitopes as the ligand molecule [95]. Firstly, we used the protein preparation wizard of UCSF Chimera (Version 1.11.2) to prepare the protein for docking analysis by deleting the attached ligand, adding hydrogens and Gasteiger–Marsili charges [96,97]. The prepared file was then added to the AutoDock wizard of PyRx 0.8 and converted into the pdbqt format. The energy form of the ligand was minimized and converted to the pdbqt format by OpenBabel [98]. The parameters used for the docking simulation were set to the default. The size of the grid box in AutoDock Vina was kept at 50.183 Å × 50.183 Å × 50.183 Å respectively, for X, Y and Z-axis. AutoDock Vina was implemented via the shell script offered by AutoDock Vina developers [99]. Docking results were observed by the negative score in kcal/mol, as the binding affinity of ligands are depicted in negative energies [100,101]. In addition, for validation of the docking approach, we selected 9-mer peptide from the envelope glycoprotein gp160 from human immunodeficiency virus (HIV) type 1 attached with the crystal structure of HLA-A*68:02 as a positive control and performed molecular docking analysis using the aforementioned similar parameters.

4.9. B-Cell Epitope Identification

The prediction of B-cell epitopes was performed to find the potential antigen that assures humoral immunity. To detect the B-cell epitope, various tools from IEDB were used to identify the B-cell antigenicity, together with the Emini surface accessibility prediction, Kolaskar and Tongaonkar antigenicity scale, Karplus and Schulz flexibility prediction and Bepipred linear epitope prediction analysis and since antigenic parts of a protein belonging to the beta-turn regions, the Chou and Fasman beta-turn prediction tool was also used [102,103,104,105,106,107].

5. Conclusions

The advancement in immunoinformatics has now emerged as a potential field for the prediction of epitope-based vaccines. As viruses can delineate both T-cell and humoral immunity, our predicted epitope might suggest enhancing the immunity against SARS-CoV-2. The assumption is based on the basic principles of immunity, which confers the attachment of the virus with the host cell, evoking immune responses and transfers the information to a broad spectrum of T cells and B cells. Our investigated epitopes mimicked the interaction to CD8+ cells antigen presentation using computational approaches. However, our study was an introductory design to predict epitope-based vaccine against SARS-CoV-2 and we hope that this predicted epitope would assist the further laboratory analysis for designing and predicting novel candidates against COVID-19.

Supplementary Materials

The following are available online. Data S1: Multiple sequence alignment of SARS-CoV-2 nucleocapsid protein; Figure S1: Evolutionary divergence analysis of available N proteins of different strains; results are represented in a phylogenetic tree; Figure S2: (A) Ramachandran plot analysis for model 1, (B) Ramachandran plot analysis for model 2, (C) Ramachandran plot analysis for model 3, (D) Ramachandran plot analysis for model 5, (E) Ramachandran plot analysis for the crystal structure of SARS-CoV-2 N protein (PDB ID; 6M3M), (F) Z-score for model 1, (G) Z-score for model 2, (H) Z-score for model 3, (I) Z-score for model 5 and (J) Z-score for the crystal structure of SARS-CoV-2 N protein (PDB ID; 6M3M); Figure S3: Population coverage based on MHC restriction data for (A) Central Africa, (B) East Africa, (C) East Asia and (D) North Africa—using the Immune Epitope Database analysis resource; Figure S4: Population coverage based on MHC restriction data for (A) North Africa, (B) North America, (C) Northeast Asia and (D) Oceania—using the Immune Epitope Database analysis resource; Figure S5: Population coverage based on MHC restriction data for (A) South Africa, (B) South America, (C) South Asia and (D) Southeast Asia—using the Immune Epitope Database analysis resource; Figure S6: Population coverage based on MHC restriction data for (A) Southwest Asia, (B) West Africa and (C) West Indies—using the Immune Epitope Database analysis resource. Table S1: Scores from the combined B-cell linear epitope prediction; Table S2. Predicted B-cell epitopes linear with the previous studies from Grifoni et al. and Amrun et al.

Author Contributions

Conceptualization, A.R., S.A.S., M.A.I., T.B.E. and J.S.-G.; methodology, A.R., S.A.S. and M.A.I.; software, A.R. and S.A.S.; validation, A.R., S.A., M.A.I., T.B.E. and J.S.-G.; formal analysis, A.R., M.A.I., F.B.F., B.H.K., K.K.S.M., M.R., T.B.E., and J.S.-G.; investigation, A.R., S.A.S. and M.A.I.; resources, F.B.F., B.H.K., K.K.S.M. and M.R.; data curation, A.R., S.A.S. and M.A.I.; writing—original draft preparation, A.R., S.A.S. and M.A.I.; writing—review and editing, M.M.N.U., F.N., T.B.E. and J.S.-G.; visualization, M.M.N.U. and T.B.E.; supervision, T.B.E. and J.S.-G.; project administration, M.M.N.U.; funding acquisition, T.B.E. and J.S.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, D.; Hu, B.; Hu, C.; Zhu, F.; Liu, X.; Zhang, J.; Wang, B.; Xiang, H.; Cheng, Z.; Xiong, Y.; et al. Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus--infected pneumonia in Wuhan, China. Jama 2020, 323, 1061–1069. [Google Scholar] [CrossRef] [PubMed]
  2. Li, Q.; Guan, X.; Wu, P.; Wang, X.; Zhou, L.; Tong, Y.; Ren, R.; Leung, K.S.M.; Lau, E.H.Y.; Wong, J.Y.; et al. Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. N. Eng. J. Med. 2020, 382, 1199–1207. [Google Scholar] [CrossRef] [PubMed]
  3. Gralinski, L.E.; Menachery, V.D. Return of the Coronavirus: 2019-nCoV. Viruses 2020, 12, 135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. World Health Organization. Coronavirus Disease 2019 (COVID-19): Situation Report; World Health Organization: Geneva, Switzerland, 2020. [Google Scholar]
  5. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.-M.; Wang, W.; Song, Z.-G.; Hu, Y.; Tao, Z.-W.; Tian, J.-H.; Pei, Y.-Y.; et al. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X.; et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 2020, 395, 497–506. [Google Scholar] [CrossRef] [Green Version]
  7. Chan, J.F.-W.; Yuan, S.; Kok, K.-H.; To, K.K.-W.; Chu, H.; Yang, J.; Xing, F.; Liu, J.; Yip, C.C.-Y.; Poon, R.W.-S.; et al. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: A study of a family cluster. Lancet 2020, 395, 514–523. [Google Scholar] [CrossRef] [Green Version]
  8. Chan, J.F.-W.; Kok, K.-H.; Zhu, Z.; Chu, H.; To, K.K.-W.; Yuan, S.; Yuen, K.-Y. Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan. Emerg. Microbes Infect. 2020, 9, 221–236. [Google Scholar] [CrossRef] [Green Version]
  9. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.-M.; Wang, W.; Hu, Y.; Song, Z.-G.; Tao, Z.-W.; Tian, J.-H.; Pei, Y.-Y.; et al. Complete genome characterisation of a novel coronavirus associated with severe human respiratory disease in Wuhan, China. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  10. Narayanan, K.; Huang, C.; Makino, S. SARS coronavirus accessory proteins. Virus Res. 2008, 133, 113–121. [Google Scholar] [CrossRef]
  11. Dermime, S.; Gilham, D.E.; Shaw, D.M.; Davidson, E.J.; Meziane, E.-K.; Armstrong, A.; Hawkins, R.E.; Stern, P.L. Vaccine and antibody-directed T cell tumour immunotherapy. Biochim. Biophys. Acta BBA-Rev. Cancer 2004, 1704, 11–35. [Google Scholar] [CrossRef]
  12. Meloen, R.H.; Langeveld, J.P.M.; Schaaper, W.M.M.; Slootstra, J.W. Synthetic peptide vaccines: Unexpected fulfillment of discarded hope? Biologicals 2001, 29, 233–236. [Google Scholar] [CrossRef]
  13. Channappanavar, R.; Fett, C.; Zhao, J.; Meyerholz, D.K.; Perlman, S. Virus-specific memory CD8 T cells provide substantial protection from lethal severe acute respiratory syndrome coronavirus infection. J. Virol. 2014, 88, 11034–11044. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rappuoli, R.; Black, S.; Bloom, D.E. Vaccines and global health: In search of a sustainable model for vaccine development and delivery. Sci. Transl. Med. 2019, 11, eaaw2888. [Google Scholar] [CrossRef] [PubMed]
  15. Olsson, S.-E.; Villa, L.L.; Costa, R.L.R.; Petta, C.A.; Andrade, R.P.; Malm, C.; Iversen, O.-E.; Høye, J.; Steinwall, M.; Riis-Johannessen, G.; et al. Induction of immune memory following administration of a prophylactic quadrivalent human papillomavirus (HPV) types 6/11/16/18 L1 virus-like particle (VLP) vaccine. Vaccine 2007, 25, 4931–4939. [Google Scholar] [CrossRef]
  16. Suarez, D.L.; Schultz-Cherry, S. Immunology of avian influenza virus: A review. Dev. Comp. Immunol. 2000, 24, 269–283. [Google Scholar] [CrossRef]
  17. Briney, B.; Sok, D.; Jardine, J.G.; Kulp, D.W.; Skog, P.; Menis, S.; Jacak, R.; Kalyuzhniy, O.; De Val, N.; Sesterhenn, F.; et al. Tailored immunogens direct affinity maturation toward HIV neutralizing antibodies. Cell 2016, 166, 1459–1470. [Google Scholar] [CrossRef] [Green Version]
  18. Pedersen, S.R.; Christensen, J.P.; Buus, S.; Rasmussen, M.; Korsholm, K.S.; Nielsen, M.; Claesson, M.H. Immunogenicity of HLA class I and II double restricted influenza A-derived peptides. PLoS ONE 2016, 11, e0145629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Ying, L.I.N.; Xu, S.; Yang, R.F.; Li, Y.X.; Ji, Y.Y.; He, Y.Y.; De Shi, M.; Wei, L.U.; Shi, T.L.; Jin, W.; et al. Identification of an epitope of SARS-coronavirus nucleocapsid protein. Cell Res. 2003, 13, 141–145. [Google Scholar]
  20. Rakib, A.; Sami, S.A.; Mimi, N.J.; Chowdhury, M.M.; Eva, T.A.; Nainu, F.; Paul, A.; Shahriar, A.; Tareq, A.M.; Emon, N.U.; et al. Immunoinformatics-guided design of an epitope-based vaccine against severe acute respiratory syndrome coronavirus 2 spike glycoprotein. Comput. Biol. Med. 2020, 124, 103967. [Google Scholar] [CrossRef] [PubMed]
  21. Kharisma, V.D.; Ansori, A.N.M. Construction of epitope-based peptide vaccine against SARS-CoV-2: Immunoinformatics study. J. Pure Appl. Microbiol. 2020, 14, 999–1005. [Google Scholar] [CrossRef]
  22. Peng, H.; Yang, L.T.; Wang, L.Y.; Li, J.; Huang, J.; Lu, Z.Q.; Koup, R.A.; Bailer, R.T.; Wu, C. Long-lived memory T lymphocyte responses against SARS coronavirus nucleocapsid protein in SARS-recovered patients. Virology 2006, 351, 466–475. [Google Scholar] [CrossRef] [Green Version]
  23. Seah, J.N.; Yu, L.; Kwang, J. Localization of linear B-cell epitopes on infectious bronchitis virus nucleocapsid protein. Vet. Microbiol. 2000, 75, 11–16. [Google Scholar] [CrossRef]
  24. Zhou, Y.; Hou, Y.; Shen, J.; Huang, Y.; Martin, W.; Cheng, F. Network-based drug repurposing for novel coronavirus 2019-nCoV/SARS-CoV-2. Cell Discov. 2020, 6, 1–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Parashar, N.C.; Poddar, J.; Chakrabarti, S.; Parashar, G. Repurposing of SARS-CoV nucleocapsid protein specific nuclease resistant RNA aptamer for therapeutics against SARS-CoV-2. Infect. Genet. Evol. 2020, 85, 104497. [Google Scholar] [CrossRef] [PubMed]
  26. Dalsass, M.; Brozzi, A.; Medini, D.; Rappuoli, R. Comparison of open-source reverse vaccinology programs for bacterial vaccine antigen discovery. Front. Immunol. 2019, 10, 113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Geourjon, C.; Deleage, G. SOPMA: Significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments. Bioinformatics 1995, 11, 681–684. [Google Scholar] [CrossRef]
  28. Kim, D.E.; Chivian, D.; Baker, D. Protein structure prediction and analysis using the Robetta server. Nucleic Acids Res. 2004, 32, W526–W531. [Google Scholar] [CrossRef] [Green Version]
  29. Dimitrov, I.; Flower, D.R.; Doytchinova, I. AllerTOP-a server for in silico prediction of allergens. BMC Bioinform. 2013, 14, S4. [Google Scholar] [CrossRef] [Green Version]
  30. Saha, S.; Raghava, G.P.S. Prediction of continuous B-cell epitopes in an antigen using recurrent neural network. Proteins Struct. Funct. Bioinform. 2006, 65, 40–48. [Google Scholar] [CrossRef]
  31. Liu, Y.; Gayle, A.A.; Wilder-Smith, A.; Rocklöv, J. The reproductive number of COVID-19 is higher compared to SARS coronavirus. J. Travel Med. 2020, 27, taaa021. [Google Scholar] [CrossRef] [Green Version]
  32. Wang, T.; Du, Z.; Zhu, F.; Cao, Z.; An, Y.; Gao, Y.; Jiang, B. Comorbidities and multi-organ injuries in the treatment of COVID-19. Lancet 2020, 395, e52. [Google Scholar] [CrossRef]
  33. Cao, B.; Wang, Y.; Wen, D.; Liu, W.; Wang, J.; Fan, G.; Ruan, L.; Song, B.; Cai, Y.; Wei, M.; et al. A trial of lopinavir--ritonavir in adults hospitalized with severe Covid-19. N. Engl. J. Med. 2020, 382, 1787–1799. [Google Scholar] [CrossRef]
  34. Chen, L.; Xiong, J.; Bao, L.; Shi, Y. Convalescent plasma as a potential therapy for COVID-19. Lancet Infect. Dis. 2020, 20, 398–400. [Google Scholar] [CrossRef]
  35. Chhabra, H.S.; Bagaraia, V.; Keny, S.; Kalidindi, K.K.V.; Mallepally, A.; Dhillon, M.S.; Malhotra, R.; Rajasekharan, S. COVID-19: Current Knowledge and Best Practices for Orthopaedic Surgeons. Indian J. Orthop. 2020, 54, 917–918. [Google Scholar] [CrossRef] [PubMed]
  36. Oany, A.R.; Emran, A.-A.; Jyoti, T.P. Design of an epitope-based peptide vaccine against spike protein of human coronavirus: An in silico approach. Drug Des. Devel. Ther. 2014, 8, 1139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Hasan, A.; Hossain, M.; Alam, J. A computational assay to design an epitope-based Peptide vaccine against Saint Louis encephalitis virus. Bioinform. Biol. Insights 2013, 7, BBI-S13402. [Google Scholar] [CrossRef]
  38. Chakraborty, S.; Chakravorty, R.; Ahmed, M.; Rahman, A.; Waise, T.M.; Hassan, F.; Rahman, M.; Shamsuzzaman, S. A computational approach for identification of epitopes in dengue virus envelope protein: A step towards designing a universal dengue vaccine targeting endemic regions. Silico Biol. 2010, 10, 235–246. [Google Scholar] [CrossRef]
  39. Islam, R.; Sakib, M.S.; Zaman, A. A computational assay to design an epitope-based peptide vaccine against chikungunya virus. Future Virol. 2012, 7, 1029–1042. [Google Scholar] [CrossRef]
  40. Huang, Y.; Khorchid, A.; Wang, J.; Parniak, M.A.; Darlix, J.-L.; Wainberg, M.A.; Kleiman, L. Effect of mutations in the nucleocapsid protein (NCp7) upon Pr160 (gag-pol) and tRNA (Lys) incorporation into human immunodeficiency virus type 1. J. Virol. 1997, 71, 4378–4384. [Google Scholar] [CrossRef] [Green Version]
  41. Thomas, J.A.; Gorelick, R.J. Nucleocapsid protein function in early infection processes. Virus Res. 2008, 134, 39–63. [Google Scholar] [CrossRef] [Green Version]
  42. Kim, T.W.; Lee, J.H.; Hung, C.-F.; Peng, S.; Roden, R.; Wang, M.-C.; Viscidi, R.; Tsai, Y.-C.; He, L.; Chen, P.-J.; et al. Generation and characterization of DNA vaccines targeting the nucleocapsid protein of severe acute respiratory syndrome coronavirus. J. Virol. 2004, 78, 4638–4645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Sabara, M.; Frenchick, P.J.; Mullin-Ready, K.F. Rotavirus nucleocapsid protein VP6 in vaccine compositions. Biotechnol. Adv. 1995, 13, 803–804. [Google Scholar]
  44. Arthur, L.O.; Bess, J.W., Jr.; Chertova, E.N.; Rossio, J.L.; Esser, M.T.; Benveniste, R.E.; Henderson, L.E.; Lifson, J.D. Chemical inactivation of retroviral infectivity by targeting nucleocapsid protein zinc fingers: A candidate SIV vaccine. AIDS Res. Hum. Retrovir. 1998, 14, S311–S319. [Google Scholar] [PubMed]
  45. Zhao, P.; Cao, J.; Zhao, L.-J.; Qin, Z.-L.; Ke, J.-S.; Pan, W.; Ren, H.; Yu, J.-G.; Qi, Z.-T. Immune responses against SARS-coronavirus nucleocapsid protein induced by DNA vaccine. Virology 2005, 331, 128–135. [Google Scholar] [CrossRef] [Green Version]
  46. Chiou, S.-S.; Fan, Y.-C.; Crill, W.D.; Chang, R.-Y.; Chang, G.-J.J. Mutation analysis of the cross-reactive epitopes of Japanese encephalitis virus envelope glycoprotein. J. Gen. Virol. 2012, 93, 1185–1192. [Google Scholar] [CrossRef] [Green Version]
  47. Yu, Z.; Wu, S.; Zhao, W.; Ding, L.; Shiuan, D.; Chen, F.; Li, J.; Liu, J. Identification and the molecular mechanism of a novel myosin-derived ACE inhibitory peptide. Food Funct. 2018, 9, 364–370. [Google Scholar] [CrossRef] [PubMed]
  48. Mishra, S. Designing of cytotoxic and helper T cell epitope map provides insights into the highly contagious nature of the pandemic novel coronavirus SARS-CoV2. ChemRxiv 2020. [Google Scholar] [CrossRef]
  49. Liniger, M.; Zuniga, A.; Tamin, A.; Azzouz-Morin, T.N.; Knuchel, M.; Marty, R.R.; Wiegand, M.; Weibel, S.; Kelvin, D.; Rota, P.A.; et al. Induction of neutralising antibodies and cellular immune responses against SARS coronavirus by recombinant measles viruses. Vaccine 2008, 26, 2164–2174. [Google Scholar] [CrossRef]
  50. Mishra, S. T cell epitope-based vaccine design for pandemic novel coronavirus 2019-nCoV. ChemRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  51. Noorimotlagh, Z.; Karami, C.; Mirzaee, S.A.; Kaffashian, M.; Mami, S.; Azizi, M. Immune and bioinformatics identification of T cell and B cell epitopes in the protein structure of SARS-CoV-2: A systematic review. Int. Immunopharmacol. 2020, 86, 106738. [Google Scholar] [CrossRef]
  52. Caubet, J.-C.; Ponvert, C. Vaccine allergy. Immunol. Allergy Clin. 2014, 34, 597–613. [Google Scholar] [CrossRef]
  53. Kallinich, T.; Beier, K.C.; Wahn, U.; Stock, P.; Hamelmann, E. T-cell co-stimulatory molecules: Their role in allergic immune reactions. Eur. Respir. J. 2007, 29, 1246–1255. [Google Scholar] [CrossRef] [Green Version]
  54. Alberts, B.; Johnson, A.; Lewis, J.; Raff, M.; Roberts, K.; Walter, P. Molecular Biology of the Cell; Garland Science: New York, NY, USA, 2002. [Google Scholar]
  55. Ahmed, S.; Rakib, A.; Islam, M.A.; Khanam, B.H.; Faiz, F.B.; Paul, A.; Chy, M.N.U.; Bhuiya, N.M.M.A.; Uddin, M.M.N.; Ullah, S.M.A.; et al. In vivo and in vitro pharmacological activities of Tacca integrifolia rhizome and investigation of possible lead compounds against breast cancer through in silico approaches. Clin. Phytosci. 2019, 5, 36. [Google Scholar] [CrossRef]
  56. Zhang, X.W. A combination of epitope prediction and molecular docking allows for good identification of MHC class I restricted T-cell epitopes. Comput. Biol. Chem. 2013, 45, 30–35. [Google Scholar] [CrossRef]
  57. Knudsen, N.P.H.; Nørskov-Lauritsen, S.; Dolganov, G.M.; Schoolnik, G.K.; Lindenstrøm, T.; Andersen, P.; Agger, E.M.; Aagaard, C. Tuberculosis vaccine with high predicted population coverage and compatibility with modern diagnostics. Proc. Natl. Acad. Sci. USA 2014, 111, 1096–1101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Lazoura, E.; Lodding, J.; Farrugia, W.; Ramsland, P.A.; Stevens, J.; Wilson, I.A.; Pietersz, G.A.; Apostolopoulos, V. Enhanced major histocompatibility complex class I binding and immune responses through anchor modification of the non-canonical tumour-associated mucin 1-8 peptide. Immunology 2006, 119, 306–316. [Google Scholar] [CrossRef] [PubMed]
  59. Pietersz, G.A.; Pouniotis, D.S.; Apostolopoulos, V. Design of peptide-based vaccines for cancer. Curr. Med. Chem. 2006, 13, 1591–1607. [Google Scholar] [CrossRef]
  60. Liu, Z.; Xiao, Y.; Chen, Y.-H. Epitope-vaccine strategy against HIV-1: Today and tomorrow. Immunobiology 2003, 208, 423–428. [Google Scholar] [CrossRef]
  61. Vivona, S.; Gardy, J.L.; Ramachandran, S.; Brinkman, F.S.L.; Raghava, G.P.S.; Flower, D.R.; Filippini, F. Computer-aided biotechnology: From immuno-informatics to reverse vaccinology. Trends Biotechnol. 2008, 26, 190–200. [Google Scholar] [CrossRef]
  62. Larsen, M.V.; Lundegaard, C.; Lamberth, K.; Buus, S.; Lund, O.; Nielsen, M. Large-scale validation of methods for cytotoxic T-lymphocyte epitope prediction. BMC Bioinform. 2007, 8, 424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Peters, B.; Bulik, S.; Tampe, R.; Van Endert, P.M.; Holzhütter, H.-G. Identifying MHC class I epitopes by predicting the TAP transport efficiency of epitope precursors. J. Immunol. 2003, 171, 1741–1749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Peters, B.; Sette, A. Generating quantitative models describing the sequence specificity of biological processes with the stabilized matrix method. BMC Bioinform. 2005, 6, 132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Lee, C.H.; Koohy, H. In silico identification of vaccine targets for 2019-nCoV. F1000Research 2020, 9, 145. [Google Scholar] [CrossRef] [PubMed]
  66. Zhou, M.; Xu, D.; Li, X.; Li, H.; Shan, M.; Tang, J.; Wang, M.; Wang, F.-S.; Zhu, X.; Tao, H.; et al. Screening and identification of severe acute respiratory syndrome-associated coronavirus-specific CTL epitopes. J. Immunol. 2006, 177, 2138–2145. [Google Scholar] [CrossRef] [PubMed]
  67. Kiyotani, K.; Toyoshima, Y.; Nemoto, K.; Nakamura, Y. Bioinformatic prediction of potential T cell epitopes for SARS-Cov-2. J. Hum. Genet. 2020, 65, 569–575. [Google Scholar] [CrossRef] [PubMed]
  68. Khalili, S.; Jahangiri, A.; Borna, H.; Ahmadi Zanoos, K.; Amani, J. Computational vaccinology and epitope vaccine design by immunoinformatics. Acta Microbiol. Immunol. Hung. 2014, 61, 285–307. [Google Scholar] [CrossRef]
  69. Vihinen, M.; Torkkila, E.; Riikonen, P. Accuracy of protein flexibility predictions. Proteins Struct. Funct. Bioinform. 1994, 19, 141–149. [Google Scholar] [CrossRef]
  70. Hopp, T.P.; Woods, K.R. Prediction of protein antigenic determinants from amino acid sequences. Proc. Natl. Acad. Sci. USA 1981, 78, 3824–3828. [Google Scholar] [CrossRef] [Green Version]
  71. Nardin, E.H.; Calvo-Calle, J.M.; Oliveira, G.A.; Nussenzweig, R.S.; Schneider, M.; Tiercy, J.-M.; Loutan, L.; Hochstrasser, D.; Rose, K. A totally synthetic polyoxime malaria vaccine containing Plasmodium falciparum B cell and universal T cell epitopes elicits immune responses in volunteers of diverse HLA types. J. Immunol. 2001, 166, 481–489. [Google Scholar] [CrossRef] [Green Version]
  72. Shey, R.A.; Ghogomu, S.M.; Esoh, K.K.; Nebangwa, N.D.; Shintouo, C.M.; Nongley, N.F.; Asa, B.F.; Ngale, F.N.; Vanhamme, L.; Souopgui, J. In-silico design of a multi-epitope vaccine candidate against onchocerciasis and related filarial diseases. Sci. Rep. 2019, 9, 1–18. [Google Scholar] [CrossRef] [Green Version]
  73. Grifoni, A.; Sidney, J.; Zhang, Y.; Scheuermann, R.H.; Peters, B.; Sette, A. A sequence homology and bioinformatic approach can predict candidate targets for immune responses to SARS-CoV-2. Cell Host Microbe 2020, 27, 671–680. [Google Scholar] [CrossRef] [PubMed]
  74. Amrun, S.N.; Lee, C.Y.-P.; Lee, B.; Fong, S.-W.; Young, B.E.; Chee, R.S.-L.; Yeo, N.K.-W.; Torres-Ruesta, A.; Carissimo, G.; Poh, C.M.; et al. Linear B-cell epitopes in the spike and nucleocapsid proteins as markers of SARS-CoV-2 exposure and disease severity. EBioMedicine 2020, 58, 102911. [Google Scholar] [CrossRef]
  75. Vapalahti, O.; Kallio-Kokko, H.; Närvänen, A.; Julkunen, I.; Lundkvist, Å.; Plyusnin, A.; Lehvaslaiho, H.; Brummer-Korvenkontio, M.; Vaheri, A.; Lankinen, H. Human B-cell epitopes of Puumala virus nucleocapsid protein, the major antigen in early serological response. J. Med. Virol. 1995, 46, 293–303. [Google Scholar] [CrossRef] [PubMed]
  76. Lundkvist, Å.K.E.; Kallio-Kokko, H.; Sjölander, K.B.; Lankinen, H.; Niklasson, B.O.; Vaheri, A.; Vapalahti, O. Characterization of Puumala virus nucleocapsid protein: Identification of B-cell epitopes and domains involved in protective immunity. Virology 1996, 216, 397–406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Choi, K.-S.; Nah, J.-J.; Ko, Y.-J.; Kang, S.-Y.; Yoon, K.-J.; Jo, N.-I. Antigenic and immunogenic investigation of B-cell epitopes in the nucleocapsid protein of peste des petits ruminants virus. Clin. Diagn. Lab. Immunol. 2005, 12, 114–121. [Google Scholar] [CrossRef] [Green Version]
  78. Khan, M.K.; Zaman, S.; Chakraborty, S.; Chakravorty, R.; Alam, M.M.; Bhuiyan, T.R.; Rahman, M.J.; Fernández, C.; Qadri, F.; Seraj, Z.I. In silico predicted mycobacterial epitope elicits in vitro T-cell responses. Mol. Immunol. 2014, 61, 16–22. [Google Scholar] [CrossRef] [PubMed]
  79. Khan, F.; Srivastava, V.; Kumar, A. Epitope based peptide prediction from proteome of enterotoxigenic E. coli. Int. J. Pept. Res. Ther. 2018, 24, 323–336. [Google Scholar] [CrossRef]
  80. Lundegaard, C.; Lamberth, K.; Harndahl, M.; Buus, S.; Lund, O.; Nielsen, M. NetMHC-3.0: Accurate web accessible predictions of human, mouse and monkey MHC class I affinities for peptides of length 8–11. Nucleic Acids Res. 2008, 36, W509–W512. [Google Scholar] [CrossRef] [Green Version]
  81. Flower, D.R. Immunoinformatics and the in Silico Prediction of Immunogenicity. In Immunoinformatics; Springer: Berlin/Heidelberg, Germany, 2007; pp. 1–15. [Google Scholar]
  82. Johnson, M.; Zaretskaya, I.; Raytselis, Y.; Merezhuk, Y.; McGinnis, S.; Madden, T.L. NCBI BLAST: A better web interface. Nucleic Acids Res. 2008, 36, W5–W9. [Google Scholar] [CrossRef]
  83. Li, W.; Cowley, A.; Uludag, M.; Gur, T.; McWilliam, H.; Squizzato, S.; Park, Y.M.; Buso, N.; Lopez, R. The EMBL-EBI bioinformatics web and programmatic tools framework. Nucleic Acids Res. 2015, 43, W580–W584. [Google Scholar] [CrossRef] [Green Version]
  84. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Doytchinova, I.A.; Flower, D.R. VaxiJen: A server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinform. 2007, 8, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Gupta, S.; Kapoor, P.; Chaudhary, K.; Gautam, A.; Kumar, R.; Raghava, G.P.S. Peptide toxicity prediction. In Computational Peptidology; Springer: New York, NY, USA, 2015; pp. 143–157. [Google Scholar]
  87. Laskowski, R.A.; MacArthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 1993. [Google Scholar] [CrossRef]
  88. Wiederstein, M.; Sippl, M.J. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007, 35, W407–W410. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Buus, S.; Lauemøller, S.L.; Worning, P.; Kesmir, C.; Frimurer, T.; Corbet, S.; Fomsgaard, A.; Hilden, J.; Holm, A.; Brunak, S. Sensitive quantitative predictions of peptide-MHC binding by a ’Query by Committee’artificial neural network approach. Tissue Antigens 2003, 62, 378–384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Giguère, S.; Drouin, A.; Lacoste, A.; Marchand, M.; Corbeil, J.; Laviolette, F. MHC-NP: Predicting peptides naturally processed by the MHC. J. Immunol. Methods 2013, 400, 30–36. [Google Scholar] [CrossRef]
  91. Bui, H.-H.; Sidney, J.; Li, W.; Fusseder, N.; Sette, A. Development of an epitope conservancy analysis tool to facilitate the design of epitope-based diagnostics and vaccines. BMC Bioinform. 2007, 8, 361. [Google Scholar] [CrossRef] [Green Version]
  92. Moutaftsi, M.; Peters, B.; Pasquetto, V.; Tscharke, D.C.; Sidney, J.; Bui, H.-H.; Grey, H.; Sette, A. A consensus epitope prediction approach identifies the breadth of murine T CD8+-cell responses to vaccinia virus. Nat. Biotechnol. 2006, 24, 817–819. [Google Scholar] [CrossRef]
  93. Maupetit, J.; Derreumaux, P.; Tuffery, P. PEP-FOLD: An online resource for de novo peptide structure prediction. Nucleic Acids Res. 2009, 37, W498–W503. [Google Scholar] [CrossRef] [Green Version]
  94. Guex, N.; Peitsch, M.C. SWISS-MODEL and the Swiss-Pdb Viewer: An environment for comparative protein modeling. Electrophoresis 1997, 18, 2714–2723. [Google Scholar] [CrossRef]
  95. Dallakyan, S. PyRx-python prescription v. 0.8. Scripps Res. Inst. 2008, 2010. [Google Scholar] [CrossRef]
  96. Dunbrack, R.L. Rotamer libraries in the 21st century. Curr. Opin. Struct. Biol. 2002, 12, 431–440. [Google Scholar] [CrossRef]
  97. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  98. O’Boyle, N.M.; Banck, M.; James, C.A.; Morley, C.; Vandermeersch, T.; Hutchison, G.R. Open Babel: An open chemical toolbox. J. Cheminform. 2011, 3, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  99. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2009, 31, 455–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Rakib, A.; Paul, A.; Chy, M.N.U.; Sami, S.A.; Baral, S.K.; Majumder, M.; Tareq, A.M.; Amin, M.N.; Shahriar, A.; Uddin, M.Z.; et al. Biochemical and Computational Approach of Selected Phytocompounds from Tinospora crispa in the Management of COVID-19. Molecules 2020, 25, 3936. [Google Scholar] [CrossRef]
  101. Jahan, I.; Tona, M.R.; Sharmin, S.; Sayeed, M.A.; Tania, F.Z.; Paul, A.; Chy, M.N.U.; Rakib, A.; Emran, T.B.; Simal-Gandara, J. GC-MS Phytochemical Profiling, Pharmacological Properties, and In Silico Studies of Chukrasia velutina Leaves: A Novel Source for Bioactive Agents. Molecules 2020, 25, 3536. [Google Scholar] [CrossRef]
  102. Chou, P.Y.; Fasman, G.D. Empirical predictions of protein conformation. Annu. Rev. Biochem. 1978, 47, 251–276. [Google Scholar]
  103. Emini, E.A.; Hughes, J.V.; Perlow, D.; Boger, J. Induction of hepatitis A virus-neutralizing antibody by a virus-specific synthetic peptide. J. Virol. 1985, 55, 836–839. [Google Scholar] [CrossRef] [Green Version]
  104. Kolaskar, A.S.; Tongaonkar, P.C. A semi-empirical method for prediction of antigenic determinants on protein antigens. FEBS Lett. 1990, 276, 172–174. [Google Scholar] [CrossRef] [Green Version]
  105. Karplus, P.A.; Schulz, G.E. Prediction of chain flexibility in proteins. Naturwissenschaften 1985, 72, 212–213. [Google Scholar] [CrossRef]
  106. Larsen, J.E.P.; Lund, O.; Nielsen, M. Improved method for predicting linear B-cell epitopes. Immunome Res. 2006, 2, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Rini, J.M.; Schulze-Gahmen, U.; Wilson, I.A. Structural evidence for induced fit as a mechanism for antibody-antigen recognition. Science 1992, 255, 959–965. [Google Scholar] [CrossRef]
Molecules 25 05088 g001 550
Figure 1. Workflow of the methodologies used in peptide vaccine design by utilizing SARS-CoV-2 nucleocapsid (N) protein.
Figure 1. Workflow of the methodologies used in peptide vaccine design by utilizing SARS-CoV-2 nucleocapsid (N) protein.
Molecules 25 05088 g001
Molecules 25 05088 g002 550
Figure 2. (A) 3D structure of modeled SARS-CoV-2 N protein; (B) composition of the secondary structure from amino acid residues of SARS-CoV-2 N protein; (C) predicted secondary structure of SARS-CoV-2 N protein; (D) Z-score of the SARS-CoV-2 N protein predicted by PROSA server and (E) Ramachandran plot analysis of the SARS-CoV-2 N protein.
Figure 2. (A) 3D structure of modeled SARS-CoV-2 N protein; (B) composition of the secondary structure from amino acid residues of SARS-CoV-2 N protein; (C) predicted secondary structure of SARS-CoV-2 N protein; (D) Z-score of the SARS-CoV-2 N protein predicted by PROSA server and (E) Ramachandran plot analysis of the SARS-CoV-2 N protein.
Molecules 25 05088 g002
Molecules 25 05088 g003 550
Figure 3. (A) Three-dimensional representation of the predicted epitope, NTASWFTAL, and (B) three-dimensional representation of the HLA-A*68:02 molecule.
Figure 3. (A) Three-dimensional representation of the predicted epitope, NTASWFTAL, and (B) three-dimensional representation of the HLA-A*68:02 molecule.
Molecules 25 05088 g003
Molecules 25 05088 g004 550
Figure 4. 3D representation of molecular docking studies representing the binding affinity of the predicted epitope, NTASWFTAL to the groove of the HLA-A*68:02. The interacting A chain residues are displayed as red ball and stick, interacting B chain residues are displayed as cyan ball and stick, hydrogen bonds are displayed as green dotted lines, alkyl/pi-alkyl bonds are displayed as pink dotted lines and carbon-hydrogen bonds are displayed as white dotted lines.
Figure 4. 3D representation of molecular docking studies representing the binding affinity of the predicted epitope, NTASWFTAL to the groove of the HLA-A*68:02. The interacting A chain residues are displayed as red ball and stick, interacting B chain residues are displayed as cyan ball and stick, hydrogen bonds are displayed as green dotted lines, alkyl/pi-alkyl bonds are displayed as pink dotted lines and carbon-hydrogen bonds are displayed as white dotted lines.
Molecules 25 05088 g004
Molecules 25 05088 g005 550
Figure 5. 3D representation of molecular docking studies representing the binding affinity of the 9-mer peptide from envelope glycoprotein gp160 form HIV type 1 (positive control) to the groove of the HLA-A*68:02. The interacting A chain residues are displayed as red ball and stick, interacting B chain residues are displayed as cyan ball and stick, hydrogen bonds are displayed as green dotted lines, alkyl/pi-alkyl bonds are displayed as pink dotted lines and salt bridges are displayed as gold dotted lines.
Figure 5. 3D representation of molecular docking studies representing the binding affinity of the 9-mer peptide from envelope glycoprotein gp160 form HIV type 1 (positive control) to the groove of the HLA-A*68:02. The interacting A chain residues are displayed as red ball and stick, interacting B chain residues are displayed as cyan ball and stick, hydrogen bonds are displayed as green dotted lines, alkyl/pi-alkyl bonds are displayed as pink dotted lines and salt bridges are displayed as gold dotted lines.
Molecules 25 05088 g005
Molecules 25 05088 g006 550
Figure 6. Combined B-cell linear epitope prediction using (A) Bepipred linear epitope prediction, (B) Chou and Fasman beta-turn prediction and (C) Emini surface accessibility prediction methods.
Figure 6. Combined B-cell linear epitope prediction using (A) Bepipred linear epitope prediction, (B) Chou and Fasman beta-turn prediction and (C) Emini surface accessibility prediction methods.
Molecules 25 05088 g006
Molecules 25 05088 g007 550
Figure 7. Combined B-cell linear epitope prediction using (A) Karplus and Schulz flexibility prediction, (B) Kolaskar and Tongaonkar antigenicity and (C) Parker hydrophilicity prediction methods.
Figure 7. Combined B-cell linear epitope prediction using (A) Karplus and Schulz flexibility prediction, (B) Kolaskar and Tongaonkar antigenicity and (C) Parker hydrophilicity prediction methods.
Molecules 25 05088 g007
Table
Table 1. Toxicity prediction of selected epitopes.
Table 1. Toxicity prediction of selected epitopes.
EpitopesToxicity PredictionSVM ScoreHydrophobicityHydrophilicityMolecular Weight
LSPRWYFYYNon-Toxin−1.08−0.06−1.261294.59
GTTLPKGFYNon-Toxin−1.13−0.01−0.49983.26
DLSPRWYFYNon-Toxin−1.18−0.14−0.671246.5
SSPDDQIGYNon-Toxin−0.33−0.20.3981.1
LLNKHIDAYNon-Toxin−0.81−0.09−0.281086.39
GTDYKHWPQNon-Toxin−0.29−0.29−0.041131.34
SPDDQIGYYNon-Toxin−0.54−0.170.011057.2
NTASWFTALNon-Toxin−10.08-11010.23
Table
Table 2. The potential CD8+ T-cell epitopes along with their interacting MHC class I alleles and total processing score, epitopes conservancy_hits and pMHC-I immunogenicity score.
Table 2. The potential CD8+ T-cell epitopes along with their interacting MHC class I alleles and total processing score, epitopes conservancy_hits and pMHC-I immunogenicity score.
EpitopesNetCTL
Combined Score		Epitope
Conservancy
Hit
(MAX. Identity %)		MHC-I Interaction with
an Affinity of IC50 < 200
and the Total Score (Proteasome
Score, TAP Score, MHC-I Score,
Processing Score)		pMHC-I
Immunogenicity
Score
LSPRWYFYY2.3408100HLA-A*29:02 (1.32),
HLA-A*30:02 (0.8),
HLA-A*01:01 (0.66),
HLA-C*16:01 (0.26)		0.35734
NTASWFTAL0.9521100HLA-A*68:02 (1.11),
HLA-C*16:01 (0.18),
HLA-C*03:03 (0.12),
HLA-C*03:04 (0.12),
HLA-C*12:03 (0.10),
HLA-A*02:06 (0.04),
HLA-C*03:02 (−0.07),
HLA-A*26:01 (−0.13),
HLA-C*14:02 (−0.34)		0.22775
DLSPRWYFY1.4994100HLA-A*29:02 (0.99)0.25933
SPDDQIGYY1.1404100HLA-B*35:01 (0.52)0.06844
SSPDDQIGY0.6895100 0.0634
Table
Table 3. Analysis of the population coverage for the proposed epitope against SARS-CoV-2.
Table 3. Analysis of the population coverage for the proposed epitope against SARS-CoV-2.
PopulationCoverage (%) aAverage Hit bPC90 c
Central Africa35.310.400.15
East Africa39.250.450.16
East Asia57.560.710.24
Europe42.950.500.18
North Africa42.150.490.17
North America45.320.530.18
Northeast Asia48.110.550.19
Oceania31.430.340.15
South Africa33.910.380.15
South America38.660.440.16
South Asia36.530.410.16
Southeast Asia49.450.570.20
Southwest Asia28.530.320.14
West Africa56.220.670.23
West Indies12.890.130.11
Notes: a Projected population coverage. b Average number of epitope hits/HLA combinations recognized by the population. c Minimum number of epitope hits/HLA combinations recognized by 90% of the population.
Table
Table 4. Results of the molecular docking analysis amongst HLA-A*68:02 and the predicted epitope, NTASWFTAL, and 9-mer peptide from envelope glycoprotein gp160 from HIV type 1 (positive control).
Table 4. Results of the molecular docking analysis amongst HLA-A*68:02 and the predicted epitope, NTASWFTAL, and 9-mer peptide from envelope glycoprotein gp160 from HIV type 1 (positive control).
EpitopesDocking Score (kcal/mol)
NTASWFTAL−9.4
Positive Control−8.2
Table
Table 5. List of predicted B cell epitopes from BepiPred linear epitope prediction analysis.
Table 5. List of predicted B cell epitopes from BepiPred linear epitope prediction analysis.
StartEndPeptideLength
361390KTFPPTEPKKDKKKKADETQALPQRQKKQQ30
338347KLDDKDPNFK10
323331EVTPSGTWL9
273287AFGRRGPEQTQGNFG15
232269SKMSGKGQQQQGQTVTKKSAAEASKKPRQKRTATKAYN38
164216GTTLPKGFYAEGSRGGSQASSRSSSRSRNSSRNSTPGSSRGTSPARMAGNGGD53
137154GALNTPKDHIGTRNPANN18
115127TGPEAGLPYGANK13
93104RIRGGDGKMKDL12
5885QHGKEDLKFPRGQGVPINTNSSPDDQIG28
151MSDNGPQNQRNAPRITFGGPSDSTGSNQNGERSGARSKQRRPQGLPNNTAS51
Table
Table 6. List of predicted B-cell epitopes from Kolaskar and Tongaonkar antigenicity prediction method.
Table 6. List of predicted B-cell epitopes from Kolaskar and Tongaonkar antigenicity prediction method.
StartEndPeptideLength
5259WFTALTQH8
6975GQGVPIN7
8389QIGYYRR7
106115PRWYFYYLGT10
130136IIWVATE7
154166NAAIVLQLPQGTT13
217227AALALLLLDRL11
243249GQTVTKK7
267273AYNVTQA7
299315KHWPQIAQFAPSASAFF17
333339YTGAIKL7
347363KDQVILLNKHIDAYKTF17
379385TQALPQR7
389401QQTVTLLPAADLD13
403411FSKQLQQSM9
Table
Table 7. List of predicted B-cell epitopes from Emini surface accessibility prediction method.
Table 7. List of predicted B-cell epitopes from Emini surface accessibility prediction method.
StartEndPeptideLength
411NGPQNQRN8
3642RSKQRRP7
185197RSSSRSRNSSRNS13
254264ASKKPRQKRTA11
277282RGPEQT6
295300GTDYKH6
340346DDKDPNF7
365377PTEPKKDKKKKAD13
384390QRQKKQQ7
Sample Availability: Samples of the compounds are available from the authors.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rakib, A.; Sami, S.A.; Islam, M.A.; Ahmed, S.; Faiz, F.B.; Khanam, B.H.; Marma, K.K.S.; Rahman, M.; Uddin, M.M.N.; Nainu, F.; et al. Epitope-Based Immunoinformatics Approach on Nucleocapsid Protein of Severe Acute Respiratory Syndrome-Coronavirus-2. Molecules 2020, 25, 5088. https://doi.org/10.3390/molecules25215088

AMA Style

Rakib A, Sami SA, Islam MA, Ahmed S, Faiz FB, Khanam BH, Marma KKS, Rahman M, Uddin MMN, Nainu F, et al. Epitope-Based Immunoinformatics Approach on Nucleocapsid Protein of Severe Acute Respiratory Syndrome-Coronavirus-2. Molecules. 2020; 25(21):5088. https://doi.org/10.3390/molecules25215088

Chicago/Turabian Style

Rakib, Ahmed, Saad Ahmed Sami, Md. Ashiqul Islam, Shahriar Ahmed, Farhana Binta Faiz, Bibi Humayra Khanam, Kay Kay Shain Marma, Maksuda Rahman, Mir Muhammad Nasir Uddin, Firzan Nainu, and et al. 2020. "Epitope-Based Immunoinformatics Approach on Nucleocapsid Protein of Severe Acute Respiratory Syndrome-Coronavirus-2" Molecules 25, no. 21: 5088. https://doi.org/10.3390/molecules25215088

APA Style

Rakib, A., Sami, S. A., Islam, M. A., Ahmed, S., Faiz, F. B., Khanam, B. H., Marma, K. K. S., Rahman, M., Uddin, M. M. N., Nainu, F., Emran, T. B., & Simal-Gandara, J. (2020). Epitope-Based Immunoinformatics Approach on Nucleocapsid Protein of Severe Acute Respiratory Syndrome-Coronavirus-2. Molecules, 25(21), 5088. https://doi.org/10.3390/molecules25215088

Article Metrics

Back to TopTop