Next Article in Journal
Secondary Metabolism Rearrangements in Linum usitatissimum L. after Biostimulation of Roots with COS Oligosaccharides from Fungal Cell Wall
Next Article in Special Issue
Computer-Aided Screening for Potential Coronavirus 3-Chymotrypsin-like Protease (3CLpro) Inhibitory Peptides from Putative Hemp Seed Trypsinized Peptidome
Previous Article in Journal
Functional Protein Composition in Femoral Glands of Sand Lizards (Lacerta agilis)
Previous Article in Special Issue
A Comprehensive Computational Investigation into the Conserved Virulent Proteins of Shigella species Unveils Potential Small-Interfering RNA Candidates as a New Therapeutic Strategy against Shigellosis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 27
Issue 7
10.3390/molecules27072375
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

Potential Immunogenic Activity of Computationally Designed mRNA- and Peptide-Based Prophylactic Vaccines against MERS, SARS-CoV, and SARS-CoV-2: A Reverse Vaccinology Approach

by
Taimoor Khan
1,
Abbas Khan
1,
Jawad Khaliq Ansari
2,
Muzammil Hasan Najmi
2,
Dong-Qing Wei
1,3,4,
Khalid Muhammad
5,* and
Yasir Waheed
2,*
1
Department of Bioinformatics and Biological Statistics, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
2
Foundation University Medical College, Foundation University Islamabad, Islamabad 46000, Pakistan
3
Peng Cheng Laboratory, Vanke Cloud City Phase I Building 8, Xili Street, Nashan District, Shenzhen 518055, China
4
State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
5
Department of Biology, College of Sciences, United Arab Emirates University, Al Ain 15551, United Arab Emirates
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(7), 2375; https://doi.org/10.3390/molecules27072375
Submission received: 8 February 2022 / Revised: 26 March 2022 / Accepted: 28 March 2022 / Published: 6 April 2022
(This article belongs to the Special Issue Structural Biology Approach in Drug Discovery)
Browse Figures
Versions Notes

Abstract

:
The continued emergence of human coronaviruses (hCoVs) in the last few decades has posed an alarming situation and requires advanced cross-protective strategies against these pandemic viruses. Among these, Middle East Respiratory Syndrome coronavirus (MERS-CoV), Severe Acute Respiratory Syndrome coronavirus (SARS-CoV), and Severe Acute Respiratory Syndrome coronavirus-2 (SARS-CoV-2) have been highly associated with lethality in humans. Despite the challenges posed by these viruses, it is imperative to develop effective antiviral therapeutics and vaccines for these human-infecting viruses. The proteomic similarity between the receptor-binding domains (RBDs) among the three viral species offers a potential target for advanced cross-protective vaccine designs. In this study, putative immunogenic epitopes including Cytotoxic T Lymphocytes (CTLs), Helper T Lymphocytes (HTLs), and Beta-cells (B-cells) were predicted for each RBD-containing region of the three highly pathogenic hCoVs. This was followed by the structural organization of peptide- and mRNA-based prophylactic vaccine designs. The validated 3D structures of these epitope-based vaccine designs were subjected to molecular docking with human TLR4. Furthermore, the CTL and HTL epitopes were processed for binding with respective human Lymphocytes Antigens (HLAs). In silico cloning designs were obtained for the prophylactic vaccine designs and may be useful in further experimental designs. Additionally, the epitope-based vaccine designs were evaluated for immunogenic activity through immune simulation. Further studies may clarify the safety and efficacy of these prophylactic vaccine designs through experimental testing against these human-pathogenic coronaviruses.

1. Introduction

Human coronaviruses (hCoVs) form a diverse group of single-stranded RNA viruses [1]. Among the four genera of CoVs (alpha, beta, gamma, and delta), α-coronavirus and β-coronavirus are highly pathogenic to humans [2]. This is reflected by the severe disease outcomes caused in humans due to highly pathogenic hCoVs that mainly involve MERS-CoV, SARS-CoV, and SARS-CoV-2 [3,4]. These three spherical-shaped hCoVs affect human lungs by targeting the lower respiratory tract and may ultimately lead to organ dysfunction and fatality [5]. Previously, the emergence and transmission of MERS-CoV and SARS-CoV caused catastrophic effects around the globe [6,7]. Despite the associated threats of these pandemic-causing hCoVs, the recent global outbreak of SARS-CoV-2 further aggravated the situation [8]. These human-infecting CoVs have a genome size spanning from 26 to 32 kb [9,10]. This includes MERS-CoV, which has a genome size of approximately 30 kb, while SARS-CoV and SARS-CoV-2 have a smaller size of approximately 29 kb [11]. The genomic similarity of the recently emerged SARS-CoV-2 with SARS-CoV and MERS-CoV is approximately 79.5% and 50%, respectively [12]. Moreover, these hCoVs have a typical genomic organization that consists of 5′-methylated caps and 3′-polyadenylated tails [8,13]. The 3′-terminal genomic region is mainly responsible for encoding the viral life cycle regulatory proteins (spike protein, nucleocapsid protein, membrane protein, and envelope protein). On the other hand, the 5′-terminal genomic region encodes for the nonstructural proteins that are involved in viral replication [13,14]. The structural and functional profiling of related proteins of hCoVs provides a strong basis to understand viral pathogenesis. On the basis of this knowledge, advanced cross-protective strategies may be designed to combat these human-pathogenic CoVs.
The proteomic organization of hCoVs includes the spike protein, which has an important structural and functional role. The spike protein mainly plays a key role in mediating the CoV’s interaction with the host cell [15]. Spike proteins contain two main subunits called S1 and S2. The receptor-binding domain (RBD) in the S1 subunit mainly regulates the binding of hCoVs with host receptors. Moreover, the S2 subunit mediates the fusion of hCoVs with host membranes and regulates the penetration of the viral genome into the cytoplasm of the host cells [16]. The host cell entry of the three highly pathogenic hCoVs is mainly mediated by the binding of the RBD to the host’s cell surface functional receptors [17]. Despite the important function of binding receptors, spike proteins also regulate the process of hCoVs’ host cell entry [18]. Additionally, for spike protein priming, these viruses deploy the cellular serine protease TMPRSS2 and the endosomal cysteine proteases cathepsin B and L to mediate the entry of the virus into host cells [17,19].
The proteomic similarity of the three highly pathogenic hCoVs may be utilized to set up precise strategies to develop cross-protective therapeutics against these deadly pathogens [20,21,22]. Although remarkable advances have been made in therapeutics development against lethal hCoVs [23,24,25], advanced efforts are needed to develop prophylactic vaccines against them. Previously, there were no available cross-protective vaccines against pandemic-causing hCoVs. Herein, immunoinformatics-based approaches were adapted to design a prophylactic vaccine against three highly pathogenic hCoVs. In this work, the RBD-containing region of the three species was targeted for potential vaccine designs (peptides and mRNA). This was achieved through domain-specific potential epitope screening followed by the structural design of two different prophylactic vaccine designs. Herein, several linkers (EAAK, AAY, GPGPG, and KK) were used to join the adjuvant and different immune epitopes together [24]. The peptide-based prophylactic vaccine design was evaluated for potential immunogenic activity through molecular docking, in silico cloning, and immune simulation approaches. Additionally, the predicted epitopes were evaluated through molecular docking to explore the interaction patterns of different epitopes with HLA (Human Leukocyte Antigen) molecules. Our study provided in silico, strategy-based prophylactic vaccine designs against three hCoVs. The aim of the study was to design a prophylactic vaccine that included epitopes from the different SARS species. Due to the cross-protective properties of the designed vaccine, it is suitable for being used as an advanced therapeutic against any of these viruses. Moreover, such a vaccine design provides protection from coinfections of SARS from different strains and may be able to produce immunity against a wider range of viral species. Further studies with experimental validations may clarify the clinical manifestations associated with the vaccine designs against these highly pathogenic human coronaviruses.

2. Results

2.1. Sequences Retrieval

The spike glycoprotein containing the RBD sequences of MERS-CoV, SARS-CoV, and SARS-CoV-2 was downloaded from UniProt. The accession details of each domain of the three different coronaviruses that were evaluated for further screening for putative T-cell, B-cell, HTL, and IFN epitopes are given in Table 1.

2.2. Prioritization of Putative Vaccine Epitopes

Different epitopes recognized by the immune system were characterized as a prerequisite for peptide-based vaccine designs. Herein, each RBD protein sequence of MERS-CoV, SARS-CoV, and SARS CoV-2 was subjected to putative vaccine epitope screening. This was followed by the selection of T-cell, HTL, and B-cell epitopes for the final prophylactic vaccine designs.

2.2.1. T-Cell Epitopes

An immune response mainly requires T-cell-specific recognition of antigenic epitopes for robust induction. These T-cell-recognized epitopes are vital to combat invading pathogens from developing infections [26]. Moreover, the potential of T-cell epitope prediction based on binding specificity through in silico approaches offers wide utilities [27]. T-cell epitopes were predicted for each RBD sequence of the three different coronaviruses. The total number of identified T-cell epitopes for MERS-RBD was 11. In addition, nine epitopes were predicted for SARS-RBD and for SARS-CoV-2-RBD. All the identified T-cell epitopes were further screened to select putative vaccine epitopes. After evaluating the potential antigenicity and allergenicity status of each predicted T-cell epitope, only the one with high immunogenic potential (Table 2) was included in further vaccine-designing procedures.

2.2.2. HTL Epitopes

HTL epitopes play a pivotal role in the induction of a protective immune response against invading pathogens [28]. They are involved in B-cell activation and maturation. The inclusion of HTL epitopes in peptide-based vaccine constructs helps in the development of safe and effective therapies against human pathogens [29]. Herein, the HTL epitope screening for each RBD of the different CoV species was performed. The numbers of identified HTL epitopes for MERS, SARS-CoV, and SARS-CoV-2 were 1351, 1456, and 1463, respectively. All the identified HTL epitopes were further screened to select putative vaccine epitopes. After mapping HTL epitopes on the basis of antigenic potential, selected epitopes (Table 3) were included in further vaccine design procedures.

2.2.3. B-Cell Epitopes

B-cell epitope mapping plays a key role in robust immune response induction. This is due to the important role of B-cells in the regulation of an immune response through the production of antibodies against the invading pathogen [30,31]. This B-cell response is triggered by the identification of epitope sequences by the antigen. Moreover, the identification of highly immunogenic B-cell epitopes is a prerequisite in peptide-based vaccine designs [32,33,34]. The total number of B-cell epitopes was predicted for each CoV species. The numbers of identified B-cell epitopes for MERS, SARS-CoV, and SARS-CoV-2 were 18, 24, and 23, respectively. All the identified B-cell epitopes were further screened to select putative vaccine epitopes. After mapping, selected T-cell epitopes with high immunogenic potential (Table 4) were included in further vaccine design procedures.

2.3. Epitopes–HLAs Molecular Docking

Peptides are presented to T-cells by HLA molecules, which results in the T-cell immune response. HLA molecules are highly polymorphic in nature; hence, peptides capable of interacting with HLA molecules are potent vaccine candidates [35]. In silico molecular docking was carried out using the HawkDock server to analyze the binding behavior of the CTL and HTL epitopes with their respective HLAs. The list of species-specific CTLs and HTLs with their respective HLAs is provided in Table 5. The 3D models of nine HLAs were retrieved from the protein databank, which included HLA-B*57:01 (PDB ID: 3X11), HLA-B*35:01 (PDB ID: 4PRA), HLA-A*01:01 (PDB ID: 6MTM), HLA-B*58:01 (5VWJ), HLA-B*15:01 (6UZP), HLA-A*24:02 (7MJA), HLA-A*02:06 (3OXR), HLA-B*51:01 (4MJI), and HLA-A*68:01 (6PBH), and prepared for docking. Higher binding energies obtained after the docking analysis of the studied CTL and HTL epitopes with their corresponding HLAs indicated their potential ability to be utilized in the vaccine design. The details, including peptide sequences, corresponding HLAs, and binding energies for each of the docking complexes, are given in Table 5. We characterized CTLs as MHC Class I (MHC-I) and HTLs as MHC Class II (MHC-II) binding epitopes with their corresponding HLAs.
For the docking complexes of the CTL epitopes with their corresponding HLAs, the total binding free energies were −27.86 kcal/mol for HLA-B*57:01-MTEQLQMGF, −20.65 kcal/mol for HLA-B*35:01-NATKFPSVY, and −21.87 kcal/mol for HLA-A*01:01-VGGNYNYLY. The CTL–HLA docking complexes are shown in Figure 1; the binding free energy calculations and interaction patterns are given in Table S1 (see Supplementary Material). For the docking complexes of the HTL epitopes with their corresponding HLAs, the total binding free energies were −37.1 kcal/mol for HLA-B*58:01-LSMKSDLSVSSAGPI, −31.49 kcal/mol for HLA-B*15:01-QFNYKQSFSNPTCLI, −33.44 kcal/mol for HLA-A*24:02-NYNYKYRYLRHGKLR, −35.7 kcal/mol for HLA-A*02:06-SGDVVRFPNITNLCP, −29.69 kcal/mol for HLA-B*51:01-TESIVRFPNITNLCP, and −46.27 kcal/mol for HLA-A*68:01-RFASVYAWNRKRISN. The HTL–HLA docking complexes are shown in Figure S1 (see Supplementary Material); the binding patterns of HTL–HLAs are given in Figure S2 (see Supplementary Material), and the binding free energy calculations are given in Table S2 (see Supplementary Material).

2.4. Prophylactic MEVC Designs

The use of computational vaccinology approaches accelerates the design of protective vaccines with increased safety and efficacy. Advanced trends including the development of predisposed biological resources have been adapted to accelerate vaccine designing against human-pathogenic species [22,23,24]. The identification of highly specific vaccine constructs on the basis of the inclusion of immunogenic epitopes [21] is vital for the induction of robust immune responses against pathogens. These vaccine constructs are required to be highly immunogenic and nonallergenic. Moreover, the use of proper adjuvants and linkers is required for the proper designing of the final vaccine construct. After screening, all the putative epitopes were subjected to the construction of multiepitope-based and mRNA-based constructs. Two different prophylactic vaccines named MEVC-mRNA-CoV and MEVC-CoV were designed against the three target CoV species. This was achieved through joining the different putative epitopes (Table 2, Table 3 and Table 4) with the help of different linkers, i.e., EAAK, AAY, GPGPG, and KK, for epitope-based vaccine designs. The use of these linkers increases structural stability by avoiding self-folding and enhances the immunogenicity of protein sequences [36]. All the constructed vaccine sequences were further evaluated for immunogenicity potential. This was performed through the analysis of each vaccine sequence on the basis of antigenicity and allergenicity status. The details of all the designed multiepitope-based vaccine constructs along with immunogenicity status are given in Table 5.

2.4.1. mRNA-Based Vaccine Construction

An mRNA-based vaccine was designed by targeting the RBD of the selected three hCoV species. This comprised joining together the epitopes predicted for the RBDs of all the studied species in the final vaccine construct (MEVC-mRNA-CoV). The final vaccine construct contained four CTL epitopes, four HTL epitopes, and six linear B-cell epitopes. The selected CTL, HTL, and linear B-cell epitopes were placed together using AAY, PMGLP, and GGGGS linkers, respectively. The 5′ m7G cap was placed right before the NCA-7d (5′ UTR) followed by the Kozak sequence and signal peptide (tissue Plasminogen Activator). The mRNA vaccine sequence spanned from the tPA to the last LBL epitope, where a stop codon was added right before the S27a+R3U (3′ UTR) and poly (A) tail (120 nucleotides long). The selected CTL, HTL, and B-cell epitopes are provided in Table 2, Table 3 and Table 4, respectively. The arrangements of the mRNA vaccine components are graphically illustrated in Figure 2.

2.4.2. Multiepitope-Based Vaccine Construct

The final vaccine construct containing the putative CTL, HTL, and B-cell epitopes was joined using different linkers, as shown in Figure 3A. Next, the linear vaccine constructs were subjected to 3D modeling for the depiction of the MEVC structures. By utilizing the structural modeling package Robetta, all the MEVC structures were modeled on the basis of the amino acid sequences. This was performed by submitting these sequences into the online server using the Robetta Fold option, and the obtained structures were refined using PyMOL. The different 3D models generated for each MEVC are given in Figure 3B. Moreover, all the modeled structures for each MEVC were validated through Ramachandran and ProSA-Web analyses. These analyses revealed a Z score of −5.24 (Figure 3C) and a large number (>80%) of residues in the favored region (Figure 3D), reflecting high confidence for the structure. Moreover, the calculated ERRAT quality factor was >85, and approximately 80 percent of the residues averaged a 3D–1D score ≥0.2. All these positive scores indicated a stable nature and validated the vaccine structure.
The total length of MEVC-Prophylactic-CoV was 329 amino acids, and it had a high antigenicity score of 0.73 and showed a nonallergenic nature. The final prophylactic vaccine was evaluated for numerous physiochemical properties to ensure the safety and efficacy of the vaccine design. During this process, JCat (http://www.jcat.de, accessed on 26 November 2021) was utilized; after providing the input protein sequence, we selected E. coli strain-K12 as the target organism. This was followed by further analyses and choosing the “Avoid rho-independent transcription terminators” and “Avoid prokaryotic ribosome binding sites” options in the online server. The GC content (51.67%), CAI value (0.96), molecular weight (35kd), theoretical pI (9.82), and several other analyzed parameters of MEVC-Prophylactic-CoV are shown in Table S3 (see Supplementary Material). These results reflected the hydrophilic and stable nature of the designed MEVC.

2.5. In Silico Cloning of the MEVC

Cloning is vital to acquire the optimized DNA sequence and insert it properly in the E. coli expression vector before proceeding into further experimental procedures. Using in silico modeling approaches, we designed the corresponding clone for the MEVC-CoV construct. This was achieved through codon optimization, followed by confirmation of putative expression in the E. coli (K-12) strain as represented by the CAI value of the construct. The optimized sequence was added to a Pet28a + vector with the utility of SnapGene software (from Insightful Science; available at snapgene.com). The specific sites of restriction enzymes selected for insertion of the optimized sequence into the vector were EcoR1 and Xho1. The final plasmid map obtained for MEVC-CoV is shown in Figure 4.

2.6. Immune Simulation of the MEVC

Using an immune simulation approach, the potential immune response induction with predicted antibody titers produced against the MEVC-prophylactic-CoV construct was evaluated. As was predicted (Figure 5), the injected antigen, after achieving the highest antigen count at day 2, was slowly neutralized until day 5. Afterwards, a strong antibody (IgM and IgG) response was observed, and the highest achieved titer of >6000 counts/mL (arbitrary scale) occurred between 10 and 15 days. A high of approximately >6000 counts/mL for combined IgM and IgG, followed by an IgM-specific antibody titer of approximately 5000 counts/mL, was obtained for the MEVC. However, further demonstrations may clarify the depicted higher immunogenic potential of the designed MEVC and the development of protective immunity against hCoVs. The final immune simulation graph obtained for the prophylactic MEVC is shown in Figure 5. Moreover, the prophylactic vaccine design was checked for the potential induction of responses by interleukins (ILs), interferon (IFN), transforming growth factors (TGFs), and transforming necrosis factor (TNF). Among these factors, the highest response produced against the antigenic vaccine was by IFN-γ with a maximum peak (>400,000) observed between days 10 and 15.

2.7. Molecular Docking of MEVC with Human TLR4

The potential interactions of the designed prophylactic vaccine with human immune receptor TLR4 through molecular docking analysis were also explored. Toll-like receptors (TLRs) play a crucial role in the activation of innate immunity by detecting conserved pathogen-associated molecular patterns (PAMPs) and the production of an adaptive immune response against the invading pathogens. Specifically, TLR4 was implicated in the recognition of hCoV structural proteins related to immune response induction. The docking complexes for each of the designed prophylactic MEVCs and TLR4 are shown in Figure 6. A docking score of −42.36 kcal/mol was observed for the MEVC–TLR4 complex. Furthermore, the PDB sum [25] was used to determine the interaction patterns between the prophylactic MEVC and human TLR4. Additionally, one hydrogen bond and two salt bridges were found in the MEVC-CoV + TLR4 docking complex. The details of the binding scores and interaction patterns, including the number of hydrogen bonds and salt bridges for the docking complexes, are given in Figure 6.

3. Materials and Methods

3.1. RBD Sequence Retrieval

The amino acid sequences of the receptor-binding domains of the three target hCoVs (MERS, SARS-CoV, and SARS-CoV-2) were downloaded from the Universal Protein Knowledge Base (https://www.uniprot.org/, accessed on 2 June 2021) [37].

3.2. Putative Vaccine Epitope Screening

Putative epitopes specific to RBD sequence of each studied species were characterized on the basis of the utility of various online available servers. This was followed by immunogenic potential (antigenicity and allergenicity nature) evaluation of the individually characterized immune epitopes (CTL, HTL, B-cell) for each RBD of the studied hCoV species.

3.2.1. Prediction of CTL Epitopes

The prediction of CTL epitopes for each RBD of the three hCoV species was carried out with NetCTL 1.2 server (http://www.cbs.dtu.dk/services/NetCTL/, accessed on 3 June 2021) [38] at default threshold of 0.75. Epitope selection was performed on the basis of the obtained descending order of the combined score. The server predicts MHC-I receptor-specific binding affinity of the input epitope sequence. It utilizes an artificial neural network (ANN) approach to quantify the binding affinity of epitope toward MHC-I. The server also offers wide sensitivity and follows a standard calculation method based on 216 experimentally identified epitopes covering all 12 confirmed HLA supertypes [38].

3.2.2. Prediction of HTL Epitopes

The predictions of HTL epitopes for each RBD of the selected viral species were conducted using IEDB MHC-II module (http://www.iedb.org/, accessed on 8 June 2021) [39], corresponding to a reference set of seven HLAs. Herein, default parameters were used to allow the prediction of epitopes covering the maximum population possible. IEDB MHC-II module assigns an IC50 value for each HTL epitope that is inversely related to the binding affinity for receptor MHC-II. An IC50 score less than 50 nM means high binding affinity; less than 500 nM means intermediate binding affinity, whereas an IC50 value less than 5000 nM is related to the low binding affinity of an epitope to the MHC-II receptor [39]. The binding affinity of HTL epitopes predicted by IEDB is inversely related to the percentile rank.

3.2.3. B-Cell Epitope Prediction

The prediction of B-cell epitopes for each RBD of the studied hCoV species was performed using ABCpred (https://webs.iiitd.edu.in/raghava/abcpred/ABC_submission.html, accessed on 15 June 2021) server [40]. This server employs a new method (the kernel method) for the prediction of linear B-cell epitopes. SVM (support vector machine) is a renowned component of kernel techniques that is composed of several algorithms used for pattern analysis. The output performance of ABCpred server (AUC = 0.75) is based on an SVM in combination with amino acid pair antigenicity (AAP) (AUC = 0.7), which is utilized for the prediction of linear B-cell epitopes.

3.3. Immunogenic and MHC-I Binding Potential Evaluation

The epitopes, before inclusion in the final vaccine constructs, were evaluated for immunogenic potential (antigenicity and allergenicity) through using two different servers. Similarly, MHC-I binding potentials of the CTL and HTL epitopes were also evaluated.

3.3.1. Antigenicity Profiling

VaxiJen server (http://www.ddg-pharmfac.net/vaxijen/VaxiJen/VaxiJen.html/, accessed on 18 June 2021) [41] was employed to predict the antigenic profile of putative epitopes and final peptide-based vaccine designs. The server’s antigenic capacity prediction is based on the physicochemical features of a specific protein or peptide, and it can be used instead of alignment techniques. It also offers a quite reliable accuracy of 70–89% for antigenicity profiling. Identification of antigenic sequences is based on higher scores (>0.4) than the threshold. On the basis of these selection criteria, only highly antigenic epitopes were short-listed in final vaccine designs.

3.3.2. Allergenicity Profiling

Putative vaccine epitopes and final vaccine constructs were subjected to allergenic profiling using AllerTOP 2.0 server (https://www.ddg-pharmfac.net/AllerTOP/index.html, accessed on 21 June 2021) [42]. AllerTOP 2.0 is an improved version; numerous features were incorporated to improve the performance of the method. The server calculates the allergenic potential of a given amino acid sequence with high specificity by utilizing the k-nearest neighbors algorithm (kNN, k = 1) based on a training set containing 2427 known allergens and nonallergens. The resultant nonallergenic epitopes were further considered in vaccine designs.

3.3.3. MHC-I Binders Profiling

To predict putative binding MHC class I alleles corresponding to short-listed CTL and HTL epitopes in the final vaccine design, IEDB MHC-I binding tool (http://tools.iedb.org/mhci/, accessed on 29 June 2021) [43] was utilized. The tool uses the IEDB server’s specific consensus algorithm for prediction of binding HLA for specific epitopes.

3.4. Molecular Docking of CTL and HTL Epitopes and HLAs

The HawkDock server was used for in silico molecular docking of selected T-cell and HTL epitopes with respective HLAs. The server was provided with 3D structures (pdb format) of the HLAs as receptors and the epitopes as ligands. Herein, the 3D models of Human Leukocyte Antigens corresponding to selected T-cell and HTL epitopes were downloaded from the Protein Databank (PDB) [44]. Furthermore, the peptide structures were modeled with an online, available (https://webs.iiitd.edu.in/raghava/pepstrmod/, accessed on 9 July 2021) tool. Additionally, the HawkDock server (http://cadd.zju.edu.cn/hawkdock/, accessed on 14 July 2021) [45] was utilized to predict the binding patterns and calculate the MM/GBSA scores for the best docking complexes.

3.5. Vaccine Construction

After short-listing several putative epitopes with qualifying parameters, the process of vaccine construction took place. This included the construction of an mRNA- and a multiepitope-based prophylactic vaccine against the highly pathogenic hCoVs.

3.5.1. mRNA Vaccine Sequence Construction

The formulation of the mRNA vaccine required the addition of several critical components in the open reading frame of the construct. This included 5′ and 3′ UTRs, Kozak sequence, putative epitopes, suitable linkers, and a stop codon. The start codon (AUG) is included in the Kozak sequence [46], whereas the stop codon can be added at the end [47]. Moreover, it is important to select rigid, cleavable, and flexible linkers to allow independent operation of each subunit of the mRNA vaccine and avoid potential interference [48]. In accordance with previous studies [49,50], suitable linkers were chosen on the basis of optimal length and other physical properties. The epitopes (CTL, HTL and B-cell) in the mRNA construct were joined through linkers (AAY, PMGLP, and GGGGS, respectively) [51,52,53]. Additionally, cellular facilitation of antigen presentation of epitopes was ensured with the addition of the tissue plasminogen activator (tPA) secretory signal sequence. This was achieved after acquiring the added tPA protein sequence from UniProtKB by using accession ID: P00750 [54]. Furthermore, to ensure stability of the mRNA-based vaccine design, the addition of elements specific to eukaryotic mRNAs [55,56] was performed. This included the addition of a 5′ cap, poly (A) tail, and 5′ and 3′ untranslated regions. It is vital to choose the optimal length of the poly (A) tail required in efficient translation [57]. Herein, the length of the poly (A) tail was kept between the suggested range of 115–150 nucleotides, in line with previous reports [58]. Moreover, as poly (A) tails operate in tandem with 5′ m7G cap sequences [59], the NCA-7d at the 5′ UTR and S27a + R3U at the 3′ UTR regions were added in the construct to stabilize the proposed mRNA vaccine [60,61].

3.5.2. Multiepitope Vaccine Construction, Structural Modeling, and Validation

All the short-listed epitopes for vaccine designs were subjected to the process of multiepitope-based vaccine construction. A total of 3 CTL and 6 HTL and 6 B-cell epitopes were included in the final prophylactic multiepitope vaccine construct against hCoVs. This was achieved through the use of different linkers (AAY, GPGPG, and KK) to join the epitope sequences in the final constructs. Moreover, adjuvant human beta defensin 2 (HbD-2) was added at the N terminal to improve the stability of the vaccine construct. The finalized MEVC was evaluated for antigenicity and allergenicity status. Furthermore, the 3D structures for the designed MEVC (multiepitope vaccine construct) were modeled through using Robetta server. The structures were also validated through the use of ProSA-web (https://prosa.Servi-ces.came.sbg.acat/prosa.php, accessed on 21 July 2021) [62], PROCHECK (https://servicesn.mbi.ucla.edu/PROCHECK/, accessed on 23 July 2021), and VERIFY3D (https://servicesn.mbi.ucla.edu/Verify3D/, accessed on 29 July 2021). Additionally, the physicochemical properties were evaluated for each of the designed vaccines through Expasy [63].

3.6. Molecular Docking Analysis

The designed prophylactic-MEVC-CoV was processed for molecular docking analysis with human TLR4 using the HawkDock server. The server is a multifunctional program that integrates the HawkRank scoring function, the ATTRACT algorithm for docking, and the MM/GBSA free energy decomposition analysis. The molecular mechanical energies and scale solvent methods are combined in the MM/GBSA method [64]. Moreover, the docking analysis of selected T-cell epitopes with corresponding HLAs was also performed. Herein, each HLA 3D structure was retrieved from the Protein Databank (PDB), and HawkDock server (http://cadd.zju.edu.cn/hawkdock/, accessed on 10 August 2021) was utilized to calculate the MM/GBSA score for the best docking poses.

3.7. In Silico Cloning and Codon Optimization

The optimized DNA sequence for the final prophylactic-MEVC-CoV design was obtained through the utility of JCat server. The CAI values and GC content for the final vaccine construct sequence were also obtained. The DNA sequence was then inserted into a protein expression plasmid vector (Pet28a+) through using SnapGene. Restriction enzyme (Xho1 and EcoR1) sites were selected to perform in silico cloning and acquire the final plasmid map.

3.8. Immune Simulation

In silico immune simulation was carried out to predict the real-life immune system response to the prophylactic-MEVC-CoV design using C-ImmSim server [65]. This simulator employs machine learning (ML) and a position-specific scoring matrix (PSSM) for the prediction of immune system and epitope interactions. The smallest suggested gap between the first and second dose of most vaccines currently in use is 4 weeks [66]. Three injections, each containing 1000 units of the vaccine, were given four weeks apart for our immune simulation. For calculating simulation durations, the C-ImmSim server employs a time-step scale. Each time step on this scale corresponds to 8 h in real life. The total number of time steps for simulation was customized to 1050, and the injection points were set at time steps 1, 84, and 168. The remaining parameters were left at their default values.

4. Conclusions

In conclusion, an immunoinformatics approach was utilized to short-list putative immune epitopes for RBDs that were utilized in the final construction of an MEVC for three different hCoVs. Different epitopes, including T-cell, B-cell, and HTL epitopes, were mapped and included in the final vaccine design against the target RBDs. Two different prophylactic vaccines were designed including an mRNA-based vaccine and a proteome-wide multiepitope-based vaccine. Additionally, the short-listed T-cell and HTL epitopes were analyzed through molecular docking with their corresponding HLAs, which demonstrated high binding energies. This step was followed by structural modeling and the evaluation of the potential immunization ability of the constructed MEVC against hCoVs. This study provides new insights into the development of future vaccines against pandemic coronaviruses. The characterization of epitopes and MEVC designs with their evaluated immunogenic potential suggested that experimental processing of these vaccines against hCoVs is needed. However, further experiments may validate the immune reinforcement potential of these putative vaccine candidates.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27072375/s1. Figure S1—Showing docking complexes of HTL epitopes for each CoV specie with respective HLAs; Figure S2—Showing interaction patterns of individual HTL epitopes with respective HLAs; Table S1—Showing details of binding scores and identified interaction patterns between the different immune epitopes and respective HLAs; Table S2—Showing details free energy calculations for the different CTL and HTL epitopes and respective HLAs; Table S3—Showing physiochemical properties including of the proposed MEVC.

Author Contributions

Conceptualization, T.K., A.K. and Y.W.; methodology, T.K., A.K., D.-Q.W. and Y.W.; software, T.K., A.K. and D.-Q.W.; validation, T.K., M.H.N. and K.M.; formal analysis, T.K., A.K., J.K.A., M.H.N., D.-Q.W., K.M. and Y.W.; investigation, T.K., A.K., K.M. and Y.W.; resources, D.-Q.W., J.K.A., K.M. and Y.W.; data curation, A.K. and T.K.; writing—original draft preparation, T.K., A.K. and Y.W.; writing—review and editing, T.K., A.K., K.M. and Y.W.; visualization, J.K.A., M.H.N., K.M. and Y.W.; supervision, Y.W. and K.M.; project administration, Y.W.; funding acquisition, J.K.A. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

J.K.A. and Y.W. acknowledge funding from Foundation University Islamabad internal funding grants from FUI ORIC. K.M.’s work is supported by United Arab Emirates University Start-up Grant #G00003347 and UAEU-UPAR Grant #G00003458.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are thankful to the administrative staff of Foundation University, Islamabad.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Weiss, S.R. Forty years with coronaviruses. J. Exp. Med. 2020, 217, e20200537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Coleman, C.M.; Frieman, M.B. Coronaviruses: Important emerging human pathogens. J. Virol. 2014, 88, 5209–5212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Meo, S.; Alhowikan, A.; Al-Khlaiwi, T.; Meo, I.; Halepoto, D.; Iqbal, M.; Usmani, A.; Hajjar, W.; Ahmed, N. Novel coronavirus 2019-nCoV: Prevalence, biological and clinical characteristics comparison with SARS-CoV and MERS-CoV. Eur. Rev. Med. Pharm. Sci. 2020, 24, 2012–2019. [Google Scholar]
  4. Khan, A.; Waris, H.; Rafique, M.; Suleman, M.; Mohammad, A.; Ali, S.S.; Khan, T.; Waheed, Y.; Liao, C.; Wei, D.-Q. The Omicron (B. 1.1. 529) variant of SARS-CoV-2 binds to the hACE2 receptor more strongly and escapes the antibody response: Insights from structural and simulation data. Int. J. Biol. Macromol. 2022, 200, 438–448. [Google Scholar] [CrossRef]
  5. Vos, L.; Bruyndonckx, R.; Zuithoff, N.; Little, P.; Oosterheert, J.; Broekhuizen, B.; Lammens, C.; Loens, K.; Viveen, M.; Butler, C. Lower respiratory tract infection in the community: Associations between viral aetiology and illness course. Clin. Microbiol. Infect. 2021, 27, 96–104. [Google Scholar] [CrossRef]
  6. Zhong, N.; Zheng, B.; Li, Y.; Poon, L.; Xie, Z.; Chan, K.; Li, P.; Tan, S.; Chang, Q.; Xie, J. Epidemiology and cause of severe acute respiratory syndrome (SARS) in Guangdong, People’s Republic of China, in February, 2003. Lancet 2003, 362, 1353–1358. [Google Scholar] [CrossRef] [Green Version]
  7. Azhar, E.I.; Hui, D.S.; Memish, Z.A.; Drosten, C.; Zumla, A. The middle east respiratory syndrome (MERS). Infect. Dis. Clin. 2019, 33, 891–905. [Google Scholar] [CrossRef]
  8. Zhou, P.; Yang, X.-L.; Wang, X.-G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.-R.; Zhu, Y.; Li, B.; Huang, C.-L.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [Green Version]
  9. Saberi, A.; Gulyaeva, A.A.; Brubacher, J.L.; Newmark, P.A.; Gorbalenya, A.E. A planarian nidovirus expands the limits of RNA genome size. PLoS Pathog. 2018, 14, e1007314. [Google Scholar] [CrossRef] [Green Version]
  10. Su, S.; Wong, G.; Shi, W.; Liu, J.; Lai, A.C.; Zhou, J.; Liu, W.; Bi, Y.; Gao, G.F. Epidemiology, genetic recombination, and pathogenesis of coronaviruses. Trends Microbiol. 2016, 24, 490–502. [Google Scholar] [CrossRef] [Green Version]
  11. Al-Qahtani, A.A. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2): Emergence, history, basic and clinical aspects. Saudi J. Biol. Sci. 2020, 27, 2531–2538. [Google Scholar] [CrossRef] [PubMed]
  12. Lu, R.; Zhao, X.; Li, J.; Niu, P.; Yang, B.; Wu, H.; Wang, W.; Song, H.; Huang, B.; Zhu, N. Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding. Lancet 2020, 395, 565–574. [Google Scholar] [CrossRef] [Green Version]
  13. Yang, D.; Leibowitz, J.L. The structure and functions of coronavirus genomic 3′ and 5′ ends. Virus Res. 2015, 206, 120–133. [Google Scholar] [CrossRef]
  14. Wu, A.; Peng, Y.; Huang, B.; Ding, X.; Wang, X.; Niu, P.; Meng, J.; Zhu, Z.; Zhang, Z.; Wang, J. Genome composition and divergence of the novel coronavirus (2019-nCoV) originating in China. Cell Host Microbe 2020, 27, 325–328. [Google Scholar] [CrossRef] [Green Version]
  15. Schoeman, D.; Fielding, B.C. Coronavirus envelope protein: Current knowledge. Virol. J. 2019, 16, 69. [Google Scholar] [CrossRef] [Green Version]
  16. Ou, X.; Liu, Y.; Lei, X.; Li, P.; Mi, D.; Ren, L.; Guo, L.; Guo, R.; Chen, T.; Hu, J. Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV. Nat. Commun. 2020, 11, 1620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.-H.; Nitsche, A. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 2020, 181, 271–280.e8. [Google Scholar] [CrossRef]
  18. Cui, J.; Li, F.; Shi, Z.-L. Origin and evolution of pathogenic coronaviruses. Nat. Rev. Microbiol. 2019, 17, 181–192. [Google Scholar] [CrossRef] [Green Version]
  19. Hoffmann, M.; Kleine-Weber, H.; Pöhlmann, S. A multibasic cleavage site in the spike protein of SARS-CoV-2 is essential for infection of human lung cells. Mol. Cell 2020, 78, 779–784.e5. [Google Scholar] [CrossRef]
  20. Llanes, A.; Restrepo, C.M.; Caballero, Z.; Rajeev, S.; Kennedy, M.A.; Lleonart, R. Betacoronavirus genomes: How genomic information has been used to deal with past outbreaks and the COVID-19 pandemic. Int. J. Mol. Sci. 2020, 21, 4546. [Google Scholar] [CrossRef]
  21. Borbone, N.; Piccialli, G.; Roviello, G.N.; Oliviero, G. Nucleoside analogs and nucleoside precursors as drugs in the fight against SARS-CoV-2 and other coronaviruses. Molecules 2021, 26, 986. [Google Scholar] [CrossRef] [PubMed]
  22. George, P.J.; Tai, W.; Du, L.; Lustigman, S. The potency of an anti-MERS coronavirus subunit vaccine depends on a unique combinatorial adjuvant formulation. Vaccines 2020, 8, 251. [Google Scholar] [CrossRef] [PubMed]
  23. Peeri, N.C.; Shrestha, N.; Rahman, M.S.; Zaki, R.; Tan, Z.; Bibi, S.; Baghbanzadeh, M.; Aghamohammadi, N.; Zhang, W.; Haque, U. The SARS, MERS and novel coronavirus (COVID-19) epidemics, the newest and biggest global health threats: What lessons have we learned? Int. J. Epidemiol. 2020, 49, 717–726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Khan, A.; Khan, S.; Saleem, S.; Nizam-Uddin, N.; Mohammad, A.; Khan, T.; Ahmad, S.; Arshad, M.; Ali, S.S.; Suleman, M. Immunogenomics guided design of immunomodulatory multi-epitope subunit vaccine against the SARS-CoV-2 new variants, and its validation through in silico cloning and immune simulation. Comput. Biol. Med. 2021, 133, 104420. [Google Scholar] [CrossRef]
  25. Chowdhury, K.H.; Chowdhury, M.; Mahmud, S.; Tareq, A.M.; Hanif, N.B.; Banu, N.; Reza, A.; Emran, T.B.; Simal-Gandara, J. Drug repurposing approach against novel coronavirus disease (COVID-19) through virtual screening targeting SARS-CoV-2 main protease. Biology 2021, 10, 2. [Google Scholar] [CrossRef]
  26. Khan, T.; Khan, A.; Wei, D.-Q. MMV-db: Vaccinomics and RNA-based therapeutics database for infectious hemorrhagic fever-causing mammarenaviruses. Database 2021, 2021, baab063. [Google Scholar] [CrossRef]
  27. Khan, T.; Khan, A.; Nasir, S.N.; Ahmad, S.; Ali, S.S.; Wei, D.-Q. CytomegaloVirusDb: Multi-Omics knowledge database for Cytomegaloviruses. Comput. Biol. Med. 2021, 135, 104563. [Google Scholar] [CrossRef]
  28. Alexander, J.; Fikes, J.; Hoffman, S.; Franke, E.; Sacci, J.; Appella, E.; Chisari, F.V.; Guidotti, L.G.; Chesnut, R.W.; Livingston, B. The optimization of helper T lymphocyte (HTL) function in vaccine development. Immunol. Res. 1998, 18, 79–92. [Google Scholar] [CrossRef]
  29. Khan, A.; Khan, S.; Ahmad, S.; Anwar, Z.; Hussain, Z.; Safdar, M.; Rizwan, M.; Waseem, M.; Hussain, A.; Akhlaq, M. HantavirusesDB: Vaccinomics and RNA-based therapeutics database for the potentially emerging human respiratory pandemic agents. Microb. Pathog. 2021, 160, 105161. [Google Scholar] [CrossRef]
  30. Humayun, F.; Cai, Y.; Khan, A.; Farhan, S.A.; Khan, F.; Rana, U.I.; binte Qamar, A.; Fawad, N.; Shamas, S. Structure-guided design of multi-epitopes vaccine against variants of concern (VOCs) of SARS-CoV-2 and validation through In silico cloning and immune simulations. Comput. Biol. Med. 2022, 140, 105122. [Google Scholar] [CrossRef]
  31. Gomes, L.R.; Durans, A.M.; Napoleão-Pêgo, P.; Waterman, J.A.; Freitas, M.S.; De Sá, N.B.; Pereira, L.V.; Furtado, J.S.; Aquino, R.G.; Machado, M.C. Multiepitope Proteins for the Differential Detection of IgG Antibodies against RBD of the Spike Protein and Non-RBD Regions of SARS-CoV-2. Vaccines 2021, 9, 986. [Google Scholar] [CrossRef] [PubMed]
  32. Das, B.K.; Chakraborty, D. Epitope-Based Potential Vaccine Candidate for Humoral and Cell-Mediated Immunity to Combat Severe Acute Respiratory Syndrome Coronavirus 2 Pandemic. J. Phys. Chem. Lett. 2020, 11, 9920–9930. [Google Scholar] [CrossRef] [PubMed]
  33. Fu, D.; Zhang, G.; Wang, Y.; Zhang, Z.; Hu, H.; Shen, S.; Wu, J.; Li, B.; Li, X.; Fang, Y. Structural basis for SARS-CoV-2 neutralizing antibodies with novel binding epitopes. PLoS Biol. 2021, 19, e3001209. [Google Scholar] [CrossRef] [PubMed]
  34. Hussain, A.; Hasan, A.; Babadaei, M.M.N.; Bloukh, S.H.; Chowdhury, M.E.; Sharifi, M.; Haghighat, S.; Falahati, M. Targeting SARS-CoV2 spike protein receptor binding domain by therapeutic antibodies. Biomed. Pharmacother. 2020, 130, 110559. [Google Scholar] [CrossRef] [PubMed]
  35. Poland, G.; Ovsyannikova, I.; Jacobson, R.; Smith, D. Heterogeneity in vaccine immune response: The role of immunogenetics and the emerging field of vaccinomics. Clin. Pharmacol. Ther. 2007, 82, 653–664. [Google Scholar] [CrossRef]
  36. Shamriz, S.; Ofoghi, H.; Moazami, N. Effect of linker length and residues on the structure and stability of a fusion protein with malaria vaccine application. Comput. Biol. Med. 2016, 76, 24–29. [Google Scholar] [CrossRef]
  37. Consortium, U. UniProt: A worldwide hub of protein knowledge. Nucleic Acids Res. 2019, 47, D506–D515. [Google Scholar] [CrossRef] [Green Version]
  38. Stranzl, T.; Larsen, M.V.; Lundegaard, C.; Nielsen, M. NetCTLpan: Pan-specific MHC class I pathway epitope predictions. Immunogenetics 2010, 62, 357–368. [Google Scholar] [CrossRef] [Green Version]
  39. Singh, S.P.; Mishra, B.N. Major histocompatibility complex linked databases and prediction tools for designing vaccines. Hum. Immunol. 2016, 77, 295–306. [Google Scholar] [CrossRef]
  40. 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]
  41. 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]
  42. Dimitrov, I.; Bangov, I.; Flower, D.R.; Doytchinova, I. AllerTOP v. 2—A server for in silico prediction of allergens. J. Mol. Modeling 2014, 20, 2278. [Google Scholar] [CrossRef] [PubMed]
  43. Calis, J.J.; Maybeno, M.; Greenbaum, J.A.; Weiskopf, D.; De Silva, A.D.; Sette, A.; Keşmir, C.; Peters, B. Properties of MHC class I presented peptides that enhance immunogenicity. PLoS Comput. Biol. 2013, 9, e1003266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The protein data bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Weng, G.; Wang, E.; Wang, Z.; Liu, H.; Zhu, F.; Li, D.; Hou, T. HawkDock: A web server to predict and analyze the protein–protein complex based on computational docking and MM/GBSA. Nucleic Acids Res. 2019, 47, W322–W330. [Google Scholar] [CrossRef]
  46. Kozak, M. Point mutations define a sequence flanking the AUG initiator codon that modulates translation by eukaryotic ribosomes. Cell 1986, 44, 283–292. [Google Scholar] [CrossRef]
  47. Liu, Q. Comparative analysis of base biases around the stop codons in six eukaryotes. Biosystems 2005, 81, 281–289. [Google Scholar] [CrossRef]
  48. Chen, X.; Zaro, J.L.; Shen, W.-C. Fusion protein linkers: Property, design and functionality. Adv. Drug Deliv. Rev. 2013, 65, 1357–1369. [Google Scholar] [CrossRef] [Green Version]
  49. Pandey, R.K.; Bhatt, T.K.; Prajapati, V.K. Novel immunoinformatics approaches to design multi-epitope subunit vaccine for malaria by investigating anopheles salivary protein. Sci. Rep. 2018, 8, 1125. [Google Scholar] [CrossRef] [Green Version]
  50. 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, 4409. [Google Scholar] [CrossRef] [Green Version]
  51. Schönefuß, A.; Wendt, W.; Schattling, B.; Schulten, R.; Hoffmann, K.; Stuecker, M.; Tigges, C.; Lübbert, H.; Stichel, C. Upregulation of cathepsin S in psoriatic keratinocytes. Exp. Dermatol. 2010, 19, e80–e88. [Google Scholar] [CrossRef] [PubMed]
  52. Riese, R.J.; Mitchell, R.N.; Villadangos, J.A.; Shi, G.-P.; Palmer, J.T.; Karp, E.R.; De Sanctis, G.T.; Ploegh, H.L.; Chapman, H.A. Cathepsin S activity regulates antigen presentation and immunity. J. Clin. Investig. 1998, 101, 2351–2363. [Google Scholar] [CrossRef] [PubMed]
  53. Michel-Todó, L.; Reche, P.A.; Bigey, P.; Pinazo, M.-J.; Gascón, J.; Alonso-Padilla, J. In silico Design of an Epitope-Based Vaccine Ensemble for Chagas Disease. Front. Immunol. 2019, 10, 2698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Ahammad, I.; Lira, S.S. Designing a novel mRNA vaccine against SARS-CoV-2: An immunoinformatics approach. Int. J. Biol. Macromol. 2020, 162, 820–837. [Google Scholar] [CrossRef] [PubMed]
  55. Gallie, D.R. The cap and poly (A) tail function synergistically to regulate mRNA translational efficiency. Genes Dev. 1991, 5, 2108–2116. [Google Scholar] [CrossRef] [Green Version]
  56. Munroe, D.; Jacobson, A. mRNA poly (A) tail, a 3′enhancer of translational initiation. Mol. Cell. Biol. 1990, 10, 3441–3455. [Google Scholar]
  57. Zhao, Y.; Moon, E.; Carpenito, C.; Paulos, C.M.; Liu, X.; Brennan, A.L.; Chew, A.; Carroll, R.G.; Scholler, J.; Levine, B.L. Multiple injections of electroporated autologous T cells expressing a chimeric antigen receptor mediate regression of human disseminated tumor. Cancer Res. 2010, 70, 9053–9061. [Google Scholar] [CrossRef] [Green Version]
  58. Holtkamp, S.; Kreiter, S.; Selmi, A.; Simon, P.; Koslowski, M.; Huber, C.; Türeci, O.Z.; Sahin, U. Modification of antigen-encoding RNA increases stability, translational efficacy, and T-cell stimulatory capacity of dendritic cells. Blood 2006, 108, 4009–4017. [Google Scholar] [CrossRef]
  59. Bernstein, P.; Ross, J. Poly (A), poly (A) binding protein and the regulation of mRNA stability. Trends Biochem. Sci. 1989, 14, 373–377. [Google Scholar] [CrossRef]
  60. Wang, Z.; Day, N.; Trifillis, P.; Kiledjian, M. An mRNA stability complex functions with poly (A)-binding protein to stabilize mRNA in vitro. Mol. Cell. Biol. 1999, 19, 4552–4560. [Google Scholar] [CrossRef] [Green Version]
  61. Pourseif, M.M.; Parvizpour, S.; Jafari, B.; Dehghani, J.; Naghili, B.; Omidi, Y. A domain-based vaccine construct against SARS-CoV-2, the causative agent of COVID-19 pandemic: Development of self-amplifying mRNA and peptide vaccines. BioImpacts 2021, 11, 65. [Google Scholar] [CrossRef]
  62. 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]
  63. Gasteiger, E.; Hoogland, C.; Gattiker, A.; Wilkins, M.R.; Appel, R.D.; Bairoch, A. Protein identification and analysis tools on the ExPASy server. In The Proteomics Protocols Handbook; Humana Press: Totowa, NJ, USA, 2005; pp. 571–607. [Google Scholar]
  64. Sun, H.; Li, Y.; Tian, S.; Xu, L.; Hou, T. Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. Phys. Chem. Chem. Phys. 2014, 16, 16719–16729. [Google Scholar] [CrossRef]
  65. Rapin, N.; Lund, O.; Bernaschi, M.; Castiglione, F. Computational immunology meets bioinformatics: The use of prediction tools for molecular binding in the simulation of the immune system. PLoS ONE 2010, 5, e9862. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Castiglione, F.; Mantile, F.; De Berardinis, P.; Prisco, A. How the interval between prime and boost injection affects the immune response in a computational model of the immune system. Comput. Math. Methods Med. 2012, 2012, 842329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Molecules 27 02375 g001 550
Figure 1. Docking complexes of selected CTL epitopes for each CoV species with respective HLAs. Whereas (A–C) shows the different-colored (red, cyan, and orange) complexes, green, purple, and marine blue represent each class of HLA molecule and the small peptide (epitope) used for docking.
Figure 1. Docking complexes of selected CTL epitopes for each CoV species with respective HLAs. Whereas (A–C) shows the different-colored (red, cyan, and orange) complexes, green, purple, and marine blue represent each class of HLA molecule and the small peptide (epitope) used for docking.
Molecules 27 02375 g001
Molecules 27 02375 g002 550
Figure 2. The final mRNA vaccine construct with organized order of elements from left (5′) to right (3′). Each structural element including 5′ cap, 5′ UTR, Kozak sequence, signal protein, CTL epitopes, HTL epitopes, B-cell epitopes, 3′ UTR, and a poly (A) tail are represented by different color scheme.
Figure 2. The final mRNA vaccine construct with organized order of elements from left (5′) to right (3′). Each structural element including 5′ cap, 5′ UTR, Kozak sequence, signal protein, CTL epitopes, HTL epitopes, B-cell epitopes, 3′ UTR, and a poly (A) tail are represented by different color scheme.
Molecules 27 02375 g002
Molecules 27 02375 g003 550
Figure 3. Topological organization, 3D structures, and validation of the MEVC-CoV. (A) represents the topological organization of the vaccine design with subunits and linkers colored differently; (B) shows the 3D modeled structures of the vaccine construct; and (C,D) show the ProSA-Web and Ramachandran plot validation, respectively.
Figure 3. Topological organization, 3D structures, and validation of the MEVC-CoV. (A) represents the topological organization of the vaccine design with subunits and linkers colored differently; (B) shows the 3D modeled structures of the vaccine construct; and (C,D) show the ProSA-Web and Ramachandran plot validation, respectively.
Molecules 27 02375 g003
Molecules 27 02375 g004 550
Figure 4. Constructed plasmid map Pet28a+ with blue-colored insert for the designed MEVC.
Figure 4. Constructed plasmid map Pet28a+ with blue-colored insert for the designed MEVC.
Molecules 27 02375 g004
Molecules 27 02375 g005 550
Figure 5. A graphical image of the immune response against the studied pathogenic viruses. Figure (A) represents the antibody titers produced, while (B) represents cytokine and interleukin production upon the assessment of immune response induction potential of the designed prophylactic-MEVC-CoV.
Figure 5. A graphical image of the immune response against the studied pathogenic viruses. Figure (A) represents the antibody titers produced, while (B) represents cytokine and interleukin production upon the assessment of immune response induction potential of the designed prophylactic-MEVC-CoV.
Molecules 27 02375 g005
Molecules 27 02375 g006 550
Figure 6. Docking complex and interaction patterns of MEVC-TLR4. Figure (A) represents the docking complex for prophylactic-MEVC-CoV. Figure (B) represents the binding interactions of MEVC with human TLR4.
Figure 6. Docking complex and interaction patterns of MEVC-TLR4. Figure (A) represents the docking complex for prophylactic-MEVC-CoV. Figure (B) represents the binding interactions of MEVC with human TLR4.
Molecules 27 02375 g006
Table
Table 1. The accession ID, domain name, protein name, and position of the selected sequences in the spike RBD of each studied hCoV species.
Table 1. The accession ID, domain name, protein name, and position of the selected sequences in the spike RBD of each studied hCoV species.
Accession IDDomain NameProtein NameOrganismPosition of RBD in the Spike (GP)
K9N5Q8BetaCoV S1-CTDSpike glycoproteinMERS381–587
P59594RBDSpike glycoproteinSARS-CoV306–527
P0DTC2RBDSpike glycoproteinSARS-CoV-2319–541
Table
Table 2. The details of total identified T-cell epitopes and short-listed ones for the vaccine design against each RBD of MERS, SARS, and SARS-CoV-2.
Table 2. The details of total identified T-cell epitopes and short-listed ones for the vaccine design against each RBD of MERS, SARS, and SARS-CoV-2.
Target SpeciesNumber of Predicted T-Cell EpitopesImmunogenic T-Cell EpitopesPosition in the RBD
MERS11MTEQLQMGF183–191
SARS-CoV9NATKFPSVY25–33
SARS-CoV-29VGGNYNYLY127–135
Table
Table 3. The details of total identified HTL epitopes and short-listed epitopes for prophylactic vaccine design against each RBD of MERS, SARS, and SARS-CoV-2.
Table 3. The details of total identified HTL epitopes and short-listed epitopes for prophylactic vaccine design against each RBD of MERS, SARS, and SARS-CoV-2.
Target SpeciesNumber of Predicted
HTL Epitopes		Immunogenic
HTL Epitopes		Position in the RBD
MERS1351LSMKSDLSVSSAGPI
QFNYKQSFSNPTCLI		70–84
86–100
SARS-CoV1456NYNYKYRYLRHGKLR
SGDVVRFPNITNLCP		130–144
5–19
SARS-CoV-21463TESIVRFPNITNLCP
RFASVYAWNRKRISN		5–19
28–42
Table
Table 4. The details of total identified B-cell epitopes and short-listed epitopes for the vaccine design against each target RBD sequence of MERS, SARS-CoV, and SARS-CoV-2.
Table 4. The details of total identified B-cell epitopes and short-listed epitopes for the vaccine design against each target RBD sequence of MERS, SARS-CoV, and SARS-CoV-2.
Target SpeciesNumber of Predicted B-Cell EpitopesImmunogenic B-Cell EpitopesPosition in the RBD
MERS18TCLILATVPHNLTTIT
TVAMTEQLQMGFGITV		97–112
180–195
SARS-CoV24MGCVLAWNTRNIDATS
TTTGIGYQPYRVVVLS		112–127
180–195
SARS-CoV-223FSTFKCYGVSPTKLND
TGCVIAWNSNNLDSKV		56–71
112–127
Table
Table 5. Characterization of CTLs as MHC Class I (MHC-I) and HTLs as MHC Class II (MHC-II) binding epitopes with corresponding HLAs.
Table 5. Characterization of CTLs as MHC Class I (MHC-I) and HTLs as MHC Class II (MHC-II) binding epitopes with corresponding HLAs.
SpeciesType of PeptidePeptide SequencePosition in the RBDHLABinding Energies
MERSCTLMTEQLQMGF183–191HLA-B*57:01−27.86 kcal/mol
SARS-CoVCTLNATKFPSVY25–33HLA-B*35:01−20.65 kcal/mol
SARS-CoV-2CTLVGGNYNYLY127–135HLA-A*01:01−21.87 kcal/mol
MERSHTLLSMKSDLSVSSAGPI70–84HLA-B*58:01−37.1 kcal/mol
MERSHTLQFNYKQSFSNPTCLI86–100HLA-B*15:01−31.49 kcal/mol
SARS-CoVHTLNYNYKYRYLRHGKLR130–144HLA-A*24:02−33.44 kcal/mol
SARS-CoVHTLSGDVVRFPNITNLCP5–19HLA-A*02:06−35.7 kcal/mol
SARS-CoV-2HTLTESIVRFPNITNLCP5–19HLA-B*51:01−29.69 kcal/mol
SARS-CoV-2HTLRFASVYAWNRKRISN28–42HLA-A*68:01−46.27 kcal/mol
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khan, T.; Khan, A.; Ansari, J.K.; Najmi, M.H.; Wei, D.-Q.; Muhammad, K.; Waheed, Y. Potential Immunogenic Activity of Computationally Designed mRNA- and Peptide-Based Prophylactic Vaccines against MERS, SARS-CoV, and SARS-CoV-2: A Reverse Vaccinology Approach. Molecules 2022, 27, 2375. https://doi.org/10.3390/molecules27072375

AMA Style

Khan T, Khan A, Ansari JK, Najmi MH, Wei D-Q, Muhammad K, Waheed Y. Potential Immunogenic Activity of Computationally Designed mRNA- and Peptide-Based Prophylactic Vaccines against MERS, SARS-CoV, and SARS-CoV-2: A Reverse Vaccinology Approach. Molecules. 2022; 27(7):2375. https://doi.org/10.3390/molecules27072375

Chicago/Turabian Style

Khan, Taimoor, Abbas Khan, Jawad Khaliq Ansari, Muzammil Hasan Najmi, Dong-Qing Wei, Khalid Muhammad, and Yasir Waheed. 2022. "Potential Immunogenic Activity of Computationally Designed mRNA- and Peptide-Based Prophylactic Vaccines against MERS, SARS-CoV, and SARS-CoV-2: A Reverse Vaccinology Approach" Molecules 27, no. 7: 2375. https://doi.org/10.3390/molecules27072375

APA Style

Khan, T., Khan, A., Ansari, J. K., Najmi, M. H., Wei, D. -Q., Muhammad, K., & Waheed, Y. (2022). Potential Immunogenic Activity of Computationally Designed mRNA- and Peptide-Based Prophylactic Vaccines against MERS, SARS-CoV, and SARS-CoV-2: A Reverse Vaccinology Approach. Molecules, 27(7), 2375. https://doi.org/10.3390/molecules27072375

Article Metrics

Back to TopTop