key: cord-0821676-kd0xyk7d authors: Waqas, Muhammad; Haider, Ali; Rehman, Abdur; Qasim, Muhammad; Umar, Ahitsham; Sufyan, Muhammad; Akram, Hafiza Nisha; Mir, Asif; Razzaq, Roha; Rasool, Danish; Tahir, Rana Adnan; Sehgal, Sheikh Arslan title: Immunoinformatics and Molecular Docking Studies Predicted Potential Multiepitope-Based Peptide Vaccine and Novel Compounds against Novel SARS-CoV-2 through Virtual Screening date: 2021-02-26 journal: Biomed Res Int DOI: 10.1155/2021/1596834 sha: f01eb317e90a2ea49e635c5c30bcbe81fc190749 doc_id: 821676 cord_uid: kd0xyk7d BACKGROUND: Coronaviruses (CoVs) are enveloped positive-strand RNA viruses which have club-like spikes at the surface with a unique replication process. Coronaviruses are categorized as major pathogenic viruses causing a variety of diseases in birds and mammals including humans (lethal respiratory dysfunctions). Nowadays, a new strain of coronaviruses is identified and named as SARS-CoV-2. Multiple cases of SARS-CoV-2 attacks are being reported all over the world. SARS-CoV-2 showed high death rate; however, no specific treatment is available against SARS-CoV-2. METHODS: In the current study, immunoinformatics approaches were employed to predict the antigenic epitopes against SARS-CoV-2 for the development of the coronavirus vaccine. Cytotoxic T-lymphocyte and B-cell epitopes were predicted for SARS-CoV-2 coronavirus protein. Multiple sequence alignment of three genomes (SARS-CoV, MERS-CoV, and SARS-CoV-2) was used to conserved binding domain analysis. RESULTS: The docking complexes of 4 CTL epitopes with antigenic sites were analyzed followed by binding affinity and binding interaction analyses of top-ranked predicted peptides with MHC-I HLA molecule. The molecular docking (Food and Drug Regulatory Authority library) was performed, and four compounds exhibiting least binding energy were identified. The designed epitopes lead to the molecular docking against MHC-I, and interactional analyses of the selected docked complexes were investigated. In conclusion, four CTL epitopes (GTDLEGNFY, TVNVLAWLY, GSVGFNIDY, and QTFSVLACY) and four FDA-scrutinized compounds exhibited potential targets as peptide vaccines and potential biomolecules against deadly SARS-CoV-2, respectively. A multiepitope vaccine was also designed from different epitopes of coronavirus proteins joined by linkers and led by an adjuvant. CONCLUSION: Our investigations predicted epitopes and the reported molecules that may have the potential to inhibit the SARS-CoV-2 virus. These findings can be a step towards the development of a peptide-based vaccine or natural compound drug target against SARS-CoV-2. There are a variety of human diseases with unknown etiology. A viral parentage has been purposed for numerous diseases and also has significance to search new viruses [1] . Various difficulties have been faced which scrutinize new viruses, such as some viruses do not replicate in vitro and have cytopathic effects (CPE). The viruses that are unable to replicate in vitro leads to the failure of virus discovery. The DNA-amplified restriction fragment length polymorphism (cDNA-AFLP 4) technique helps to identify the new viruses including the discovery of new coronavirus [1] . Coronaviruses, a genus of the Coronaviridae family, are enveloped viruses recognized as of large plus RNA strand genome. The size of RNA is 27-32 kb and polyadenylated. There are three groups of coronaviruses that are serologically distinct. Viruses are characterized within each group by their genomic sequence and host range [2] . Coronaviruses have been discovered in mice, turkeys, cats, horse, and humans and cause many diseases including respiratory tract and gastroenteritis [2] . Two human viruses (HCoV-229E, HCoV-OC43) were identified in the mid-1960s and are known to cause the common cold. The recently identified SARS-CoV can cause a lifethreatening pneumonia and is the most pathogenic human coronaviruses identified thus far [3] . SARS-CoV is probable to occupy in animal source and recently initiated the epidemic in humans through zoonotic transmission [4] . SARS-CoV is the first membrane of a fourth group of coronaviruses [5] . In Wuhan (Hubei province, China), multiple patients associated to Hunan south China seafood market diagnosed with third zoonotic human coronavirus (CoV) of the century emerged in 31st of December 2019. CoV is similar to severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) infections including fever, lung infiltration, and difficulty breathing [6] . After an extensive speculation about the causative agent of CoV, the identification of novel CoV was announced by the Chinese Center for Disease Control (CDS) on 19th of January 2020 [7] . The novel CoV, SARS-CoV-2, was insulated from a single patient and later corroborated from 16 more patients [8] . The viral pneumonia of SARS-CoV-2 was quickly predicted as the likely causative agent, while not yet confirmed. The first sequence of SARS-CoV-2 has been submitted after its conformation [9] . Later, five more sequences of SARS-CoV-2 were deposited to the GSAID database on 11th of January from Chinese institutes [10] (Supplementary 1); multiple sequence alignment of SARS-CoV, MERS-CoV, and SARS-CoV-2 carried out and conserved part in DNA, as well as protein sequence, was observed. Hundreds of human deaths were linked with infection having significant morbidities with the age>50. Various clinical symptoms have been highlighted such as dry cough, leukopenia, fever, and shortness of breath. The extracorporeal membrane oxygenation of the patients considered severe cases and need supportive care. The infection of SARS-CoV-2 in elderly patients are less virulent as compared to SARS-CoV (10% mortality) and MERS-CoV (35% mortality) [11] . 1.1. Origin. The source of the SARS-CoV-2 is still unclear, although the initial cases have been associated with the Huanan South China Seafood Market. The early patients present in the Market got the virus through either human-to-human transmission or a more widespread animal source [11] . The samples from the infected market showed positive results for the novel coronavirus while no specific animal association has been identified [12] . Through codon analyses, it is suggested that the snakes might be the possible source of the viral infection [13] , although the assertion has been disputed by others [14] including possible animal vectors, and the researchers are trying to discover the source of SARS-CoV-2. Coronavirus was thought to infect humans and bats more effectively as both are more related to Coronavirus lifecycle [15] . It has been evidenced that several bats are capable of infecting human cells without intermediate adaptation [16] . The human serology data shows the association of bat CoV proteins leads to zoonotic transmission of SARS-like bat coronavirus for deadliest out breaks [17] . MERS-CoV is also a zoonotic virus and have the origin from the bats [18] . The zoonotic contacts of camel has been evidenced in primary cases of MERS-CoV [19] . These lessons from SARS and MERS highlight the importance of rapidly finding source for SARS-CoV-2 in order to stem the ongoing outbreak [19] . With low patient data, who may be most sensitive to SARS-CoV-2 is difficult to make robust resolution. Disease severity such as SARS-CoV and MERS-CoV equated strongly to host the condition including biological sex, age, and the overall health [20] , and similar findings have been observed in early patients of SARS-CoV-2. The SARS-and MERS-CoV infection leads to increase the severity and death rate in people over the age of 50 years [21] . The observed patients having novel CoV had poor health conditions including diabetes, kidney or heart function issues, and hypertension that make them more susceptible for MERS-CoV outbreak, while diabetes, smoking, cardiovascular disease, hypertension, and other chronic illness have also been observed. In the majority of deaths and corresponding to findings in animal models [22] , the results indicate that vigilance is essential for these weak patients following SARS-CoV-2 infection [22] . 1.3. Insights from the Sequence. Dr. Zhang's group at Fudan University and many other groups in China instance the dedication and increased the capacity of the scientific infrastructures in China by rapid sequencing of nearly 30,000 nucleotide of the (COVID) genome [23] . The whole genome analyses of SARS-CoV-2 showed~80% nucleotide identity to the original SARS epidemic virus. The two different bat SARS-like CoVs (ZC45 and ZXC21) shared~89% identity with the genome of SARS-CoV-2 [24] . It has been observed that the novel CoV showed recombination with previously identified bat coronaviruses through phylogenetic analyses [25] . A CoV sequence of bat (RaTG3) having 92% sequence identity with the novel virus supports the bat origins for the SARS-CoV-2 [14] . The SARS-CoV-2 spike protein has roughly 75% amino acid identity with SARS-CoV [26] while the SARS-CoV-2 receptor-binding domain (RBD) is 73% conserved with spike RBD of SARS-CoV by narrowing analysis relative to the epidemic RBD [27] . The receptor-binding domain of SARS-CoV-2 was capable of binding with ACE2 in the context of the SARS-CoV spike protein [28] . Features and Lifecycle of the Coronavirus. Coronaviruses have unique club-like spikes, and the RNA genome is larger than other virus which leads to a unique mode of replication. Coronaviruses contain~30 kb of positive-strand RNA genome [29] . The significant features of coronavirus genomes include a 5 ′ caped end which plays an important role in the replication of RNA, as 5 ′ end has a leader sequence along with a UTR region, possessing essential loops. The 3 ′ poly-A tail end has essential structures for RNA genome synthesis and replication [30] . These two modifications allow RNA viruses for translation of replication (replicase) proteins [23] . A coronavirus genome has significant parts and helps for the synthesis and replications of whole genome (Figure 1 ) [31] . The conformed cases of virus have been confirmed by 25 countries [32] [33] [34] Tables 1 and 2 (Supplementary 1) . Our current study is aimed at exploring and identifying potential B-and T-cell epitopes through immunoinformatics approaches which help to design effective vaccine against deadly SARS-CoV-2. In addition, the study is aimed at pointing out specific peptides from coronaviral proteome, which have ability to bind with major histocompatibility complex (MHC), one of the most crucial step in vaccine designing. Different bioinformatics tools are applied to follow immunoinformatics approach. Retrieval. The primary amino acid sequence of coronavirus protein was extracted from the crystal structure of SARS-CoV-2 main protease in complex with an inhibitor N3 from Protein Data Bank (PDB ID: 6LU7) [35] . The individual sequence length of corona viral protein was 306 amino acids from the genome polyprotein, and a three-dimensional (3D) structure was determined by X-ray diffraction having 2.16 Å resolution. The physiochemical properties of the selected protein were evaluated by using ProtParam [36] . MSA is performed on all three full-length genomes (SARS-CoV = NC_004718, MERS-CoV = NC_019843.3, and SARS-CoV-2 = NC_ 045512.2), all genomic sequences taken by GenBank [37, 38] and multiple sequence alignment carried out by Clustal Omega [39, 40] . The conserved parts were labeled by using WebLogo3 [41] . The interaction of the antigen B-cell epitope with Blymphocyte classifies the B-lymphocytes to differentiate into the two types of cells as memory cells and antibody-secreting plasma [42] . The accessibility and hydrophilic nature were considered the key features of the B-cell [43] by accessing the immune epitope database and analysis resource (IEDB) (http://www.iedb.org/) as stated by flexibility prediction of Karplus and Schulz [44] , hydrophilicity prediction of Parker et al. [43] , antigenicity scale of Kolaskar and Tongaonkar [45] , and Emini et al. surface accessibility prediction [46] . The conformational B-cell epitopes were predicted by employing ElliPro (http://tools.immuneepitope.org/ toolsElliPro/) [46] from the IEDB analysis resource. This analysis resource incorporates three diverse algorithms comprising protein shape approximation [47] , residues protrusion index (pI) [48] , and the adjacent residue clustering based on pI. Prediction. CTL epitopes were predicted by employing the NetCTL.1.2 server [49] . MHC molecules act as an antigen and utilize their surface to activate the CTLs. The NetCTL.1.2 server was employed to integrate the proteasomal C-terminal cleavage, MHC class I binding prediction, and transporter associated with antigen processing (TAP) transport efficiency. The sequences of the organism in FASTA format were submitted to the server, and afterwards, peptide lengths and human leukocyte antigen (HLA) alleles were selected and observed. Additionally, the T-cell epitope prediction and weight matrix algorithm were used for the TAP transport efficiency prediction, and artificial neural network was implemented to predict the proteasomal C-terminal cleavage and MHC class-I binding. 2.5. World Population Coverage Analysis. The world population coverage analysis was performed by utilizing IEDB server by utilizing the selected CTL epitopes which were searched against respective allele sets, and major world populations were covered by this analysis. The key purpose for this coverage analyses were to analyze whether the selected candidates were suitable for major populations or not. The analyses were performed against China, Iran, Japan, Korea, and some other countries which were being affected by the coronavirus in 2020 viral outbreak [50] . Studies. The predicted CTL epitope peptides of SARS-CoV-2 with antigenic residues were selected for the molecular docking analyses. The PEP-FOLD3 server [51] was employed to model the 3D structures of the selected peptides with 200 simulation runs to sample the conformations. The conformational models clustered by PEP-FOLD3 server were evaluated on the basis of sOPEP energy scores [52] . Afterwards, the peptides with higher scores were selected for molecular docking experiments with MHC class I binding molecule comprising HLA-B (PDB ID: 3VCL) through the PatchDock docking server [53] . All the docked complexes which showed the undesirable penetrations of the receptor's atoms into the ligand were rejected, and the geometric shape complementarity score was applied to classify the other complexes. Subsequently, the FireDock server [54, 55] was utilized to refine the docked complexes and also predict the score of the docking outputs. The FireDock server supports to rectify the scoring and flexibility issues generated during the docking calculations by fast rigid-body docking tools [56] . The molecular visualization programs PyMOL [45] (Schrodinger, Inc.) and UCSF Chimera 1.11 [46] were employed to analyze and identify the hydrogen-bonding interactions of the docked complexes. The observed results suggested that the followed strategy 3 BioMed Research International has the capability to identify the effective epitope-based vaccines against coronavirus SARS-CoV-2 [42, 57, 58] . The FDA-approved library was selected for virtual screening and molecular docking analyses. The selected library has 1615 FDA-approved compounds, and all the compounds were minimized through UCSF Chimera and Chemdraw to obtain the stable configurations; all these drugs were previously derived from the ZINC database. The selected library was docked against nonstructural corona virus protein (PDB: 6LU7) involved in the replication of SARS-CoV-2 genome. The molecular docking analyses were carried out through Molecular Operating Environment (MOE) [59] , AutoDock tools, and AutoDock Vina [60] . Molecular docking analyses were performed having parameters as rescoring function 1, rescoring function 2, London dG = 10, placement: triangle matcher, retain: 2, and refinement: force field = 10 for MOE. The best hits were selected based on S-score and root-mean-square deviation (RMSD) values. The admetSAR server [61] , Molinspiration [62] , and Osiris explorer [63] were used to calculate the chemical and physical properties of drug-like hits. The interacting residues were analyzed and visualized through the UCSF Chimera and Ligplot tool [64] . Replicase protein, NSP1, spikes, membrane, nucleocapsid and envelope proteins were retrieved by utilizing UniProt KB [65, 66] . HTL and CTL epitopes from the selected Figure 1 : The organization of Coronavirus genome, which contains a 5′ end, a leader sequence, replicase protein (important for replication of whole genome), spikes, envelope, membrane, nucleocapsid, and a 3′UTR poly-A-tail end. proteins were predicted by using the NETCTL server and ABCpred server [67] . Their physiochemical properties, antigenicity, toxicity, and immunogenicity were predicted by using the ProtParam, Vaxijen, Toxinpred, and IEDB servers, respectively [68, 69] . An adjuvant-based MEV construct was designed manually by using the selected 28 epitopes, and 3D structures were predicted by using RaptorX [70] . Structure validation was carried out by the SAVES server, and the refined structures were docked with TLR3 and TLR8 by using the HADDOCK server [71, 72] . The viral pneumonia with unknown etiology had an outbreak recently in Wuhan, China [13] . Severe acute respiratory syndrome (SARS), Middle East respiratory syndrome (MERS), influenza virus, and adenovirus were not involved in the outbreak of viral pneumonia [73] . The virological.org sequenced the viral RNA genome, and World Health Organization (WHO) [74] reported the designation on 10th of January 2020. Based on genetic properties, the Coronavirinae family consists four genera including alpha-coronavirus, genus beta-coronavirus, genus gamma-coronavirus, and genus delta-coronavirus (Supplementary 1) [75] . CoVs have considered as minimal responsible pathogens causing "colds" in humans. Two extremely pathogenic CoVs named as SARS-CoV and MERS-CoV were emerged from the livestock reservoirs and caused deadly outbreaks in the 21st century. A new strain of CoV was identified named as SARS-CoV-2 in Wuhan city on December 31st, 2019. Due to the rapid changing situation, the final dimension and impact of this outbreak are currently uncertain [76] . The novel virus infects the host cells rapidly, proven through recombination of various genome practices. For this infection, no reliable mediation is currently available. The preventative measures are urgently needed due to the significant global disease burden resultant of SARS-CoV-2 [77] . A variety of tools and servers have resulted through recent advancement in immunological bioinformatics, which lessens the time and cost of traditional vaccine advancement. The development of an effective multiple-epitope vaccine remains difficult, due to problems in the selection of suitable antigen candidates and immune-dominant epitopes. Thus, it is important to predict the appropriate antigen epitopes of a targeted protein by immune-informatics approaches for designing a multiple-epitope vaccine [48] . The main target is to use immune-informatics approaches and the prediction of peptide vaccine through recognizing CTL epitopes. The discovery of novel vaccines is possible through [78] . To analyze the complete spectrum of the potential antigen, immune-informatics approaches help, and furthermore, complications regarding in vitro expression of antigen and pathogen culturing can also be evaded. By means of computational methods, the immune research groups have reported various vaccine candidates, having promising preclinical outputs [79] . In current efforts, CTL epitopes have been identified to design the peptide vaccine against HLA-B protein [80] . The development of epitope-based vaccines targets the structural proteins of SARS-CoV-2, and CTL epitopes of the target proteins were predicted to support the host's immune response. One nonstructural protein (PDB: 6LU7) stands with the reason to use this nonstructural protein due to involvement in the replication of the virus [81] [82] [83] [84] [85] [86] [87] . The antigenicity and allergenicity of CTL epitopes were observed through Vaxijen and Allergen F.P 1.0 [88] . The population coverage estimation of predicted epitopes was calculated, and 0.5639 coverage with average hits of 4.0 for MHC class I and 0.2462 coverage with average hits of 0.91 for MHC class II (Table 1) were observed in China. The peptides were designed against eight epitopes by utilizing PEP-FOLD3. The molecular docking analyses of the selected eight peptides were performed through PatchDock and further refined through FireDock [53] [54] [55] to identify the effective binding sites. 3.1. Surface Accessibility Analysis for SARS-CoV-2. A peptide with surface accessibility probability of >1.0 reflects more probable chances for a peptide to be found on the surface [43] . Numerous peptides were predicted, and the topranked predicted peptides of SARS-CoV-2 on the basis of surface probability (y-axis) and sequence position (x-axis) were selected for further analyses (Figure 2(a) ). The maximum surface probability score of 8.254 was observed that ranges from 97 to 102 amino acids with the hexapeptide sequence of KTPKYK, while the lowest score was 0.285 from 246 to 251 residues with the hexapeptide sequence of HVDILG (Supplementary 2). Flexibility for Protein SARS-CoV-2. The Karplus and Schulz flexibility method was utilized to calculate and analyze the atomic vibrational motions in the protein structure designated through B-factor and temperature. The stability and organization of the structure depend upon the Bfactor values. The quality of the predicted models depends upon the B-factor values as a lower B-factor value is [44] . The surface flexibility outputs for SARS-CoV-2 were critically analyzed (Figure 2(b) ), and it was observed that the minimum and maximum flexibility scores were 0.983 and 1.082 with the heptapeptide sequences of 129 AMRPNFT 135 and 106 IQPGQTF 112, respectively (Supplementary 2). The hydrophilicity scale process of Parker was carried out to observe the peptides hydrophilicity based on the peptide retention times through HPLC on reversed phase column. Immunological analyses have revealed the association of antigenic sites with the hydrophilic regions [43] . Parker's hydrophilicity of SARS-CoV-2-predicted peptides in graphical form was analyzed (Figure 2(c) ), where hydrophilicity is plotted along the y-axis and residues position is plotted along the x-axis. It was observed that the Parker hydrophilicity prediction has a maximum hydrophilicity score of 5.329 which ranges from 92 to 98 with the sequence of heptapeptide 92 DTANPKT 98 while the minimum hydrophilicity score was -4.257 which ranges from 204 to 210 with the peptide sequence 204 VLAWLYA 210 (Supplementary 2). The correlation among the protein structure antigenicity, epitope prediction, accessibility, and flexibility within 3D structure was determined through ElliPro [89] . The significant properties including protein-antibody interactions were analyzed to differentiate the predicted epitopes. The five top-ranked conformational epitopes for SARS-CoV-2 having ≥0.6 score were observed and selected for further analyses. The pI (isoelectric point value) [89] score was observed to analyze the percentage of the atoms which extends over the molecular bulk and also liable for the antibody binding. The pI value 5.95 was observed for 6LU7. The six top-ranked conformational predicted epitopes along with residues name, length, and locations were critically analyzed (Table 2) , and the score was observed between 0.51 and 0.78. The comparative molecular docking analyses were executed for 8 top-ranked selected CTL epitopes of SARS-CoV-2 out of 87 designed peptides with MHC class I HLB. The strong binding affinities have been observed for all the selected CTL epitopes having Van der Waals (VdW) energy values ranging from -23.45 to -32.62 kcal/mol, and the observed global energy was -29.63 to -50.38 kcal/mol ( Table 3 ). The molecular docking analyses of the 8 selected CTL predicted epitopes (GTDLEGNFY, TVNVLAWLY, GSVGFNIDY, QTFSVLACY, DYDCVSFCY, TANPKTPKY, SEDMLNPNY, and LLEDEFTPF) were carried out, and effective binding affinities with HLA-B were observed . TYR7, TYR9, TYR59, ARG62, ILE66, GLN70, THR73, SER77, TYR99, TYR116, THR143, TRP147, GLU152, GLN155, ARG156, TYR159, GLU163, TYR171 TYR9, ARG62, ILE66, GLN70, THR73, ASP74, SER97, TYR99, TYR116, GLU152, GLN155, ARG156, TYR159, GLU163 TYR7, TYR9, ARG62, ILE66, ALA69, GLN70, THR73, TYR99, TYR116, TRP147, GLU152, GLN155, ARG156, ALA158, TYR159 TYR9, ARG62, GLN65, ILE66, ALA69, GLN70, THR73, SER77, TYR99, ASP114, TYR116, TRP147 BioMed The top-ranked four docked complexes were visualized (Figure 3) , and similar binding pocket has been observed in all the selected peptides. It was observed that Tyr9, Ile66, Gln70, Tyr99, Tyr116, and Arg156 residues were conserved in all the selected peptides. Figure 4 ) from the selected library were common from each selected docking tool and docking approach having least binding energies (Table 4 ). Almost all the docked compounds from the FDA library bound on similar binding site. The four top-ranked complexes were elucidated (Figure 4) , and similar binding pocket was revealed in comparison with molecular docking analyses. The selected compounds may have the potential to inhibit the replication of SARS-CoV-2. It was elucidated that all the compounds bound at the domain II of SARS-CoV-2. It was observed that Asp153, Phe294, Ile152, Asn151, Val104, Arg105, Gln107, Gln110, and Ile106 residues showed effective binding interactions with all the docked compounds of the FDA library. In an effort to understand the insights of the binding interactions between the docked compounds and amino acid residues of SARS-CoV-2, a plot of interactional analyses was generated by utilizing Ligplot and UCSF Chimera ( Figure 5 ). The FDA library has all the compounds approved by the FDA and utilized for different diseases. The FDA library's aim was to select the available compounds to inhibit the replication of SARS-CoV-2 in minimal time frame. Molinspiration, admetSAR online server, and Osiris explorer were utilized for absorption, distribution, metabolism, excretion, and toxicity (ADMET) analyses of the selected compounds ( Table 4 ). The aqueous solubility prediction (defined water at 25°C) of the selected library revealed that the scrutinized molecules can be soluble in water. It was observed that the compounds have the ability to follow Lipinski's rule of five and also have less values of LogP involved in effective oral bioavailability. All the selected nine compounds showed similar binding site and highest binding affinity (Supplementary 5). The amino acid sequences of SARS-CoV-2 vaccine-target proteins (replicase protein, NSp1, envelope, membrane, nucleocapsid, and spike protein) were retrieved and saved in 11 BioMed Research International FASTA format. The VaxiJen server was used to analyze the antigenicity of the selected proteins. Spike protein was observed as the most antigenic protein, followed by E, M, NSp1, N, and replicase proteins with antigenic values of 0.7185, 0.6502, 0.6441, 0.6131, 0.6025, and 0.5102, respectively. The 3D models of the selected proteins were predicted in order to select the suitable quality models, and the predicted structures were further refined by galaxy refine server followed by the Ramachandran plot validations. Therefore, good-quality models were selected for further analyses. There was no suitable structure predicted for spike protein because of the small number of residues. 3.11. HLA-B7 Allele and Epitope Interaction Analyses. To construct a subunit vaccine, the selected epitopes should be 100% conserved, overlapping, and antigenic [91, 92] . Therefore, a total of 50 conserved/antigenic epitopes from the selected proteins overlapping in all 3 categories (B-cell, T-cell, and IFN-Γ) were selected for further validation of their interactions with a common human allele. The 3D structures of the selected epitopes were predicted by using PEP-FOLD. The binding patterns of the selected epitopes with a common conserved allele HLA-B7 were analyzed through molecular docking, and it was found that only 28 epitopes bound deep inside in the HLA-B7 binding pocket. Each bound epitope to HLA-B7 depicts stronger than -10.00 kcal/mol docking affinity. All the 28 selected epitopes showed their binding efficiency as well as their suitability to be used in multiplepitope-based vaccine construct (Table 5) . 3.12. Construction of Multiepitope-Based Vaccine. All 28 selected epitopes (replicase 3, NSp1 3, envelope 2, membrane 5, nucleocapsid 6, and spikes 9) were analyzed for interinteractions and further used to develop an MEV construct. An adjuvant (45 amino acid long ß defensin) was linked with the help of EAAAK linker at the start (to the N-terminal of BioMed Research International the MEV). The EAAAK linker reduces the interaction with other protein regions with efficient separation and increases the stability. The immunogenicity of the vaccine may increase with an adjuvant. Epitopes were merged together based on their interactional compatibility in sequential manner with AAY and GPGPG linkers, respectively. AAY and GPGPG prevent the generation of junctional epitopes, which is a major concern in the design of multiepitope vaccines. Contrarily, multiepitope vaccines facilitate the immunization and presentation of the epitopes. The final vaccine construct comprises of 479 amino acids ( Figure 6 ). 3.13. Evaluation of Multiepitope Vaccine. BlastP was performed for the proteome of Homo sapiens, and it was observed that MEV is nonhomologous. Proteins having less than 37% identity was generally considered nonhomologous [93, 94] . However, MEV showed no similarity (higher or equal to 37%) with the proteins of human. The allergenicity, antigenicity, and toxicity of the vaccine construct were evaluated. It was observed that MEV is highly antigenic (0.6741 at 0.5% threshold), nonallergenic, and nontoxic. Furthermore, the physiochemical properties of the SARS-CoV-2 MEV construct were determined by using ProtParam. To determine the tertiary structure of the vaccine, RaptorX was used and the structure was refined by Galaxy ( Figure 7) . The selected structure showed that 96.3% amino acids were in allowed region, 3.7% of residues in permitted region, and 0.0% in outer region according to the Ramachandran plot analyses. Further analyses revealed that qRMSD was 0.428, poor rotamers were 0%, MolProbity was 1.889, clash score was 13.6, and Z score was -2.25. In addition, the refined structure showed 0 errors with PROCHECK validation. The refined structure showed 85.7143% of the overall quality factor through ERRAT. The results showed the reliability of the selected structure. The Ramachandran plot analyses of the predicted MEV structure showed that 96.3% of residues were present in favorable region. Vaccine against TLR3 and TLR8. An appropriate association between immune receptor molecules and the antigen molecule is essential to activate an immune responsiveness [95] . HADDOCK has been used to perform the molecular docking analyses of the MEV with human immune receptors TLR3 and TLR8. TLR3 and TLR8 can efficiently induce the immune response after virus recognition [33, 34] . The molecular docking analyses showed effective binding interactions between MEV and TLR3/TLR8. The binding scores of MEV-TLR3 and MEV-TLR8 were observed as -293.90 kcal/mol and -283.20 kcal/mol, respectively. It was observed that MEV generated 11 hydrogen bonds within the range of 3.00 Å with TLR3. MEV-interacting amino acids with hydrogen bonding to TLR3 are shown in green-colored stick representation, while similarly, TLR3 amino acids interacting through hydrogen bonding with MEV are shown in redcolored stick representation (Figure 8 ). It was observed that MEV made 9 hydrogen bond interactions within the range of 3.00 Å with TLR8. Similar to TLR3, MEV-interacting amino acids with hydrogen bonding to TLR8 are shown in green-colored stick representation, while TLR8 amino acids interacting through hydrogen bonding with MEV are shown in red-colored stick representation (Figure 9 ). The need of dealing with coronaviruses has been increased since its recent breakout affecting millions of human lives. This SARS-CoV-2 viral outbreak became an emergency in different regions of the world [96] . As an immediate response, numerous efforts have been made to design the peptide-based vaccine against SARS-CoV-2. Peptide inhibitors are of great interest to develop vaccines [97, 98] . The peptide targets are more superior than traditional ligandbased drugs including less toxicity, fewer side-effects, and their ultrafast action. Immunoinformatics methodologies are helping researchers by reducing the workload of laboratory trials; additionally, these approaches are less timeconsuming and cost-efficient than traditional approaches [99] [100] [101] . Since the last decade, there has been much progress in in silico drug designing [102] . Numerous biological complications are being solved by the implementation of different bioinformatics approaches [80, 102, 103] . The potential CTL epitopes have been predicted for nonstructural protein (PDB: 6LU7) of SARS-CoV-2. The molecular docking tools are applied to analyze MHC-1 and ligandbinding affinities for the selected peptides [104] . Other evidences like C-terminal cleavage affinities also validate the binding affinity of peptide-MHC-I complexes. In this study, eight peptides were reported as the potential targets with effective MHC-I protein (HLA-B) interactions. Based on global energy scores, four peptides were selected having maximum binding affinities and antigenicity, increasing the probability of the potential vaccine targets for the observed residues to be a promising target. Surface accessibility and surface flexibility, as well as hydrophobicity and antigenicity , GIINTLQKYYCRVRGGRCAVLSCLPKEEQIGKCSTRG RKCCRRKKEAAAKGSVGFNIDYAAYLLEDEFTPFAAY HVGEIPVAYAAYLSEARQHLKAAYLVKPSFYVYAAYLV GLMWLSYAAYAGDSGFAAYAAYLSPRWYFYYAAYSSP DDQIGYAAYWTAGAAAYYAAYCNDPFLGVYAAYITDA VDCALAAYSTQDLFLPFAAYQLTPTWRVYAAYVLPFN DGVYAAY GPGPGTLNGLWLDDVVYCPR GPGPG GPGPGVLLFLAFVVFLL VTLGPGPG GPGPGCLLQFAYAN RNRFLYGPGPG GPGPGQIGYYRR ATRRIRGGGPGPG GPGPGATK AYNVTQAFGRRGGPGPG GPGP GQSLLIVNNATNVVIKGPGPGINITRFQTLLALHRS 14 BioMed Research International GLU78, ILE79, PRO80, ALA84, ALA85, TYR86, SER88, GLU89, PHE104, VAL112, LEU114, TRP116, TRP139, TYR140, PHE141, TYR158, TRP159, ALA161, GLY162, PHE222, TRP250, PHE291, PRO321, LEU476, HIS477, ARG478, SER479 VAL168, SER188, ALA195, ASN196, SER198, LYS200, PRO214, GLY215, HIS218, ALA219, ARG222, ALA295, TRP296, PRO298, GLN299, GLU301, TYR302, HIS319, PHE322, ASN323, LYS355, CYS356, ILE379, HIS406, HIS410, HIS432, GLU434, GLU456, ASN457 MEV vaccine MEV-interacting residues TLR3 TLR3-interacting residues Interactions Figure 8 : All interacting residues from MEV are shown in green color, and the rest of all red residues are TLR3-interacting residues. PRO80, ALA84 , TYR86, LEU114, TRP116, TRP139, TYR140, PHE141, TYR142, TYR143, TYR158, TRP159, THR160, ALA161, LEU199, LEU201, PHE222, TRP250, ASP252, PHE291, TRP318, PRO321, PRO323, CYS325, GLY342, PRO343, GLY344, ALA345, TYR347, PRO363, GLY364, GLN365, TYR368, TYR369, ARG370, THR406, LYS407, ARG417, PRO421, GLY422, PRO423, ASP426, ALA427, ALA428, ASN452, ALA453, ARG469, GLN471, THR472, LEU473, LEU474, ALA475, LEU476, HIS477, ARG478, SER479 PHE261, ASN262, PRO264, PHE346, TYR353, ARG375, ILE403, TYR424, SER426, GLU427, ARG429, PHE470, LEU490, ASN491, SER492, PHE494, SER513, ALA514, SER516, ALA518, ASN539, ARG541, TYR563, SER565, HIS566, TYR567, PHE568, ARG569, ALA571, HIS593, ASN595, TYR597, THR598, GLU612, VAL614, ARG619, ILE622, ASN625, ARG643, ASP645, SER647, LEU648, ARG650, LYS652, HIS653, HIS670, ASN672, ASP673, ASN674, MET675, LYS677, GLY697, ASN698, LYS699, LEU701, HIS721, ASN722, ARG723 MEV vaccine MEV-interacting residues TLR8 TLR8-interacting residues Interactions Figure 9 : All interacting residues from MEV shown in green color and residues of TLR8 interacting residues in red color. for SARS-CoV-2 nonstructural protein were calculated and cross-verified using the IEDB server [105] . Based on an extensive literature review, it was observed that the selected peptides were not reported against SARS-CoV-2. The predicted peptides were modeled through PEP-FOLD3 server and docked to MHC-1 using PatchDock and further refined with FireDock. PyMOL and UCSF Chimera 1.11 were used to analyze the interactions of the docked complexes [46] . The S-value is a scoring function based upon the affinity of the ligand with the receptor [59] . The compounds having higher S-value with lower values of RMSD can be developed as potential inhibitors for a target protein [106] . For further evaluation, the binding energy of these selected hits were identified. The binding affinity showed the polar interaction of the hits with the binding site of receptor, and the value observed between 5 and 15 kcal/mol is considered a strong interaction among the ligands and the receptor [107, 108] . The molecular docking was also carried out using AutoDock and AutoDock Vina [109, 110] . Multiepitope vaccine construct revealed effective binding affinities against TLR3 and TLR8. The construct contains multiple epitopes from replicase, NSp1, N, E, M, and S coronavirus proteins. Various studies have been conducted by using immunoinformatics approach leading to efficient results [111] [112] [113] [114] [115] . The aim of our work was to identify the effective peptidebased inhibitors against SARS-CoV-2 nonstructural protein (PDB: 6LU7), which plays an important role in viral genome replication. Epitopes were designed, and then molecular docking was performed against MHC-I; interactional analyses of the selected docked complexes were carried out. In conclusion, four CTL epitopes (GTDLEGNFY, TVNVLAWLY, GSVGFNIDY, and QTFSVLACY) and four FDA-scrutinized compounds indicated potential targets as a peptide vaccine and potential biomolecule against deadly SARS-CoV-2, respectively. On the other hand, a multiepitope vaccine was also designed using different epitopes of coronavirus proteins joined by linkers and led by an adjuvant, which can be a possible potential MEV against coronavirus. Our findings can be a step towards the development of a peptide-based vaccine or natural compound drug target against SARS-CoV-2 which is one of the trending issues nowadays due to the exponentially increasing death rate all over the world. Cytopathic effects SARS-CoV: Severe acute respiratory syndrome coronavirus MERS-CoV: Middle East respiratory syndrome coronavirus RBD: Receptor-binding domain MHC: Major histocompatibility complex HLA: Human leukocyte antigen MOE: Molecular Operating Environment CTL: Cytotoxic T-lymphocyte pI: Isoelectric point ADMET: Absorption distribution metabolism elimination toxicity. Authors have no conflicts of interest form anyone. MW, AH, AR, MQA, AU, MS, HNA, AM, RR, and DR performed the computational analyses, RAT analyzed the data, and SAS conceived the project, analyzed the results, and drafted the manuscript. Viral induced demyelination Characterization of a coronavirus isolated from a diarrheic foal Clinical progression and viral load in a community outbreak of coronavirusassociated SARS pneumonia: a prospective study Virology: SARS virus infection of cats and ferrets Unique and conserved features of genome and proteome of SARS-coronavirus, an early split-off from the coronavirus group 2 lineage Return of the Coronavirus: 2019-nCoV Real-time tentative assessment of the epidemiological characteristics of novel 16 Note from the editors: World Health Organization declares novel coronavirus (2019-nCoV) sixth public health emergency of international concern Discovery of SARS-CoV-2 antiviral drugs through large-scale compound repurposing Estimating the potential total number of novel coronavirus (2019-nCoV) cases in Wuhan City Molecular diagnosis of a novel coronavirus (2019-nCoV) causing an outbreak of pneumonia Homologous recombination within the spike glycoprotein of the newly identi?ed coronavirus may boost cross-species transmission from snake to human nCoV's relationship to bat coronaviruses and recombination signals no snakesJanuary 2020 Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats SARS-like WIV1-CoV poised for human emergence Serological evidence of bat SARS-related coronavirus infection in humans, China Identification of a severe acute respiratory syndrome coronavirus-like virus in a leafnosed bat in Nigeria Evidence for camel-to-human transmission of MERS coronavirus Middle East respiratory syndrome: emergence of a pathogenic human coronavirus Epidemiological, demographic, and clinical characteristics of 47 cases of Middle East respiratory syndrome coronavirus disease from Saudi Arabia: A descriptive study Risk factors for fatal Middle East respiratory syndrome coronavirus infections in Saudi Arabia: Analysis of the WHO Line List Architecture of the SARS coronavirus prefusion spike Determine the potential epitope based peptide vaccine against novel SARS-CoV-2 targeting structural proteins using immunoinformatics approaches Viral shedding patterns of coronavirus in patients with probable severe acute respiratory syndrome Discovery of a novel coronavirus associated with the recent pneumonia outbreak in humans and its potential bat origin Jumping species-a mechanism for coronavirus persistence and survival Synthetic recombinant bat SARS-like coronavirus is infectious in cultured cells and in mice Antagonism of the interferon-induced OAS-RNase L pathway by murine coronavirus ns2 protein is required for virus replication and liver pathology Cryo-electron tomography of mouse hepatitis virus: insights into the structure of the coronavirion Modular organization of SARS coronavirus nucleocapsid protein World Health Organization The application of toll like receptors for cancer therapy TLR, NLR agonists, and other immune modulators as infectious disease vaccine adjuvants RCSB Protein Data Bank: biological macromolecular structures enabling research and education in fundamental biology, biomedicine, biotechnology and energy Protein identification and analysis tools in the ExPASy server Clustal omega Clustal omega for making accurate alignments of many protein sequences WebLogo: a sequence logo generator Epitope recognition by diverse antibodies suggests conformational convergence in an antibody response New hydrophilicity scale derived from high-performance liquidchromatography peptide retention data -correlation of predicted surface residues with antigenicity and X-ray-derived accessible sites Prediction of chain flexibility in proteins bcl::Cluster : A method for clustering biological molecules coupled with visualization in the Pymol Molecular Graphics System UCSF chimera -a visualization system for exploratory research and analysis Induction of hepatitis A virus-neutralizing antibody by a virusspecific synthetic peptide Proteome-wide screening for designing a multi-epitope vaccine against emerging pathogenElizabethkingia anophelisusing immunoinformatic approaches The Immune Epitope Database (IEDB): 2018 update PEP-FOLD3: faster de novo structure prediction for linear peptides in solution and in complex A coarsegrained protein force field for folding and structure prediction PPDock-Portal Patch Dock: a web server for drug virtual screen and visualizing the docking structure by GP and X-score FireDock: a web server for fast interaction refinement in molecular docking FireDock: fast interaction refinement in molecular docking Solving and analyzing side-chain positioning problems using linear and integer programming Editorial: epitope discovery and synthetic vaccine design Epitope-based peptide vaccine design and target site depiction against Middle East respiratory syndrome coronavirus: an immune-informatics study Medicinal chemistry and the Molecular Operating Environment (MOE): application of QSAR and molecular docking to drug discovery Small-molecule library screening by docking with PyRx Estimation of ADME properties with substructure pattern recognition Molecular docking, PASS analysis, bioactivity score prediction, synthesis, characterization and biological activity evaluation of a functionalized 2-butanone thiosemicarbazone ligand and its complexes Synthesis, in vitro antifungal evaluation and in silico study of 3-azolyl-4-chromanone phenylhydrazones Ligplot -a program to generate schematic diagrams of protein ligand interactions UniProt archive The Universal Protein Resource (UniProt): an expanding universe of protein information Identification of CD8+ T cell epitopes in the West Nile virus polyprotein by reverse-immunology using NetCTL The IDB and IEDB: intron sequence and evolution databases IEDB-3D: structural data within the immune epitope database RaptorX: exploiting structure information for protein alignment by statistical inference HAD-DOCK(2P2I): a biophysical model for predicting the binding 18 BioMed Research International affinity of protein-protein interaction inhibitors Solvated protein-DNA docking using HADDOCK Outbreak of pneumonia of unknown etiology in Wuhan, China: the mystery and the miracle Fenner's veterinary virology Host factors in coronavirus replication Adaptive evolution influences the infectious dose of MERS-CoV necessary to achieve severe respiratory disease Reverse vaccinology and subtractive genomics reveal new therapeutic targets againstMycoplasma pneumoniae: a causative agent of pneumonia Harnessing bioinformatics to discover new vaccines Immunoinformatics and molecular docking studies reveal potential epitope-based peptide vaccine against DENV-NS3 protein A conserved virulence region within alphacoronavirus nsp1 The hepatitis C viral nonstructural protein 5A stabilizes growth-regulatory human transcripts A hypervariable region within the 3' cis-acting element of the murine coronavirus genome is nonessential for RNA synthesis but affects pathogenesis Kaposi's sarcoma-associated herpesvirus nonstructural membrane protein pK15 recruits the class II phosphatidylinositol 3-kinase PI3K-C2α to activate productive viral replication Nonstructural protein Pns4 of rice dwarf virus is essential for viral infection in its insect vector Nonstructural protein Pns12 of rice dwarf virus is a principal regulator for viral replication and infection in its insect vector Viperin inhibits classical swine fever virus replication by interacting with viral nonstructural 5A protein Aller-genFP: allergenicity prediction by descriptor fingerprints ElliPro: a new structure-based tool for the prediction of antibody epitopes ZINC-a free database of commercially available compounds for virtual screening Use of defined TLR ligands as adjuvants within human vaccines TLR ligand-peptide conjugate vaccines: toward clinical application Recognition of analogous and homologous protein folds: analysis of sequence and structure conservation 1 TM-align: a protein structure alignment algorithm based on the TM-score Efficient immunization and cross-priming by vaccine adjuvants containing TLR3 or TLR9 agonists complexed to cationic liposomes A new look at an old disease -smallpox and biotechnology Peptides as therapeutic agents for dengue virus Towards peptide vaccines against Zika virus: immunoinformatics combined with molecular dynamics simulations to predict antigenic epitopes of Zika viral proteins Computational design of peptide ligands Computational design of peptide ligands for ochratoxin A Computational design of peptide ligands to target the intermolecular interaction between viral envelope protein and pediatric receptor Pharmacoinformatics, adaptive evolution, and elucidation of six novel compounds for schizophrenia treatment by targeting DAOA (G72) isoforms Structural, phylogenetic and docking studies of D-amino acid oxidase activator (DAOA), a candidate schizophrenia gene From ZikV genome to vaccine: in silico approach for the epitopebased peptide vaccine against Zika virus envelope glycoprotein Predicting affinity and specificity of antigenic peptide binding to major histocompatibility class I molecules Discovery of novel dengue NS2B/NS3 protease inhibitors using pharmacophore modeling and molecular docking based virtual screening of the ZINC database Applications of in silico methods for design and development of drugs targeting protein-protein interactions A comparative study of the efficiency of HCV NS3/4A protease drugs against different HCV genotypes using in silico approaches AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading Comparing AutoDock and Vina in ligand/decoy discrimination for virtual screening Design of a multiepitope-based peptide vaccine against the E protein of human COVID-19: an immunoinformatics approach Design of a peptide-based subunit vaccine against novel coronavirus SARS-CoV-2 Design of an epitopebased peptide vaccine against spike protein of human coronavirus: an in silico approach Development of epitope-based peptide vaccine against novel coronavirus 2019 (SARS-COV-2): immunoinformatics approach Reverse vaccinology approach to design a novel multi-epitope vaccine candidate against COVID-19: anin silicostudy Authors are thankful to Mr. Jonathan Javid for the help in molecular docking analyses.