key: cord-1042858-dwlz3cfl
authors: Sirohi, Preeti Rana; Gupta, Jyoti; Somvanshi, Pallavi; Prajapati, Vijay Kumar; Grover, Abhinav
title: Multiple epitope-based vaccine prediction against SARS-CoV-2 spike glycoprotein
date: 2020-11-27
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1846626
sha: a632cab341b4c38ed7be3a15e076c47726a951c3
doc_id: 1042858
cord_uid: dwlz3cfl

The global emergence of novel coronavirus disease and its rapid global expansion over a short span of time require effective countermeasures to combat it. Development of a specific vaccine can induce an optimal antibody response, thus providing immunity against it. Our study proposes a detailed and comprehensive immunoinformatic approach that can be applied to the currently available coronavirus protein data in the online server for vaccine candidate development. We have identified the receptor binding domain (RBD) of structural spike protein (S1) as a potential target for immunity against COVID- 19 infection. Epitope prediction illustrated cytotoxic T-cell epitopes, helper T-cell epitopes, and B-cell epitopes associated with the target protein. These were joined through specific linkers along with adjuvant beta-defensin located at the N-terminal to create a multi epitope subunit vaccine (MESV). The specificity in the binding of the devised vaccine candidate to the TLR-3 immune cell receptor was evaluated via molecular docking interaction studies. Good docking score combined with robust interactions in the binding cavity certified the stringency of the engineered vaccine. Molecular dynamics simulation data showed minimal variation of the root-mean square deviations (RMSDs) and root-mean-square fluctuations (RMSFs) which confirmed the interaction stability. These results obtained from various in-silico experiments indicate the potency of this vaccine candidate as a probable therapeutic agent against COVID-19. Vaccination strategies targeting conserved epitope-based immune response would be beneficial in providing cross protection across beta-coronaviruses, and such vaccines would be resistant to the ever-evolving viruses. Communicated by Ramaswamy H. Sarma

A novel strain of the coronavirus is suspected to have emerged in Wuhan, China in late 2019 and has resulted in an emergency situation throughout world. The first cases of the novel coronavirus were reported in the Hubei province and later named the disease as COVID-19 outbreak in China (Mackenzie & Smith, 2020) . It was declared a pandemic disease on March 11, 2019 by the World Health Organization (WHO) with more than 21.2 million people being infected and resulting in more than 761,000 deaths worldwide to date (World Health Organization, 2020) due to a severe respiratory infection which induces systemic cytokine storm in the host (Song et al., 2020) . The virus acts mainly by attacking the respiratory system in three key phases: viral replication, hyper-reactive immune response, and subsequent obstruction of pulmonary airflow. According to the reports, perforated holes are created in the lungs giving them a honey-comb like appearance upon viral infection.

The severe acute respiratory syndrome-corona virus-2 (SARS-CoV-2) is a key member from family Coronaviridae, order Nidovirales. It consists of a large genome comprising a positive-sense single-stranded RNA (Fauquet, 2008) . The size of the virus is 65-125 nm with a crown like structure protruding on the surface. The coronaviridae family constitutes four subgroups: alpha, beta, gamma, and delta (Chan et al., 2015) . The novel strain of the coronavirus belongs to beta group which has low pathogenicity Perlman & Netland, 2009) . The virus has specific attributes that allow it to infect interspecies, and it infects a diverse range of hosts-birds to mammals and humans-and causes severe respiratory and enteric diseases (Lee, 2015) . The symptoms are similar to those of the common flu, such as mild to moderate fever, cough, sore throat, and headache in addition to difficulty in breathing (M€ akel€ a et al., 1998; Owusu et al., 2014) .

The virus is comprised of spike glycoprotein, membrane glycoprotein, nucleocapsid, and enveloped non-glycosylated protein (Gorbalenya et al., 2006) . These are essential to synthesize a complete virus and its associated non-structural protein components, which are critical to the process of viral replication and transcription. The spike protein primarily determines the virus infection. Virus interaction is an essential step for initiating infection, and proteolytic cleavage plays a critical role in the peptide release in the host (Belouzard et al., 2012; Ye et al., 2004) . The spike protein of the corona virus is a class I transmembrane that is trimeric, has a molecular weight of 180 kDa, and has a monomer that requires host protease cleavages for its activation, such as CoV which requires two specific cleavages for fusion of spike protein. This last factor is crucial for the interaction, fusion, and internalization of virus to host. The interaction of the spike is associated with cell and tissue tropism as well as pathogenesis. Moreover, tropism is associated with specific cleavage of spike protein, and mutants suggest the alternation in viral pathogenesis. The spike protein consists of two functional domains S1 and S2 (F. Li, 2016) . The RBD of the S1 domain contains the receptor binding motif and plays a key role in host interaction (Tai et al., 2020) . It acts together with the angiotensin converting enzyme-2 (ACE-2) receptor, a known efficient surface receptor of the host. It is the key target of antibodies (W. Li et al., 2003 Li et al., , 2005 .

Currently, there is no specific remedial measure in the form of drug/medicine or vaccine available for the treatment or termination of the spread of SARS-CoV-2 infection. With the daily increasing cases of infection and the high deathtoll, there is an imperative and immediate need to design and develop an efficient vaccine against the deadly virus to order control this global pandemic.

Antigenic epitopes of RBD of spike glycoprotein of novel coronavirus strain have been determined through utilizing immunoinformatic techniques which might be put to use for designing and developing a MESV peptide to suppress the viral infection through eliciting TLR-3 innate immune response. Epitopes, in general, are the short stretches of protein amino acid sequences which have a capability of eliciting a direct, specific and strong immune response, in comparison to the expected response generated by the complete protein itself (Palatnik-de-Sousa et al., 2018) .

In the current study, a strong and prospective vaccine candidate was developed using an immunoinformatic approach, wherein we have identified antigenic epitopes against our target viral protein and combined them. The specificity in binding for the designed vaccine to the TLR-3 receptor was evaluated via molecular docking interaction. Good docking scores combined with robust interactions in the binding cavity certified the stringency of the engineered vaccine. Certain other parameters such as physicochemical properties, antigenicity, allergenicity, and toxicity were also examined as they are the requisites for being a candidate vaccine. Molecular dynamics (MD) simulation data have been generated to investigate the interaction stability. Moreover, strong immunological response was obtained via the T-cell and B-cell mediated humoral response through immune simulation. The results obtained from various in-silico experiments indicate the potency of our vaccine candidate as a probable therapeutic in case of COVID-19.

The spike protein sequence of the SARS coronavirus for subunit vaccine designing was retrieved from NCBI (https:// www.ncbi.nlm.nih.gov) in FASTA format. This downloaded protein was cropped from residue 331-524 which corresponds to its receptor binding domain (RBD). After retrieval, the viral protein sequence was put through pBLAST tool (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE=Proteins) to verify its similarity with human proteome. Next, the antigenic propensity of the protein sequence was done by using an online available antigenic peptides prediction server of the Immunomedicine group (http://imed.med.ucm.es/Tools/antigenic.pl).

(CTL), helper T-lymphocyte (HTL) epitopes with subsequent MHC alleles and B-cell epitopes

The prediction of cytotoxic T-lymphocyte (CTL) epitopes for the preferred SARS-CoV-2 RBD of the protein was done via NetCTL1.2 tool (http://www.cbs.dtu.dk/services/NetCTL/). The helper T-lymphocyte (HTL) epitopes were anticipated with Immune Epitope Databases' (IEDB) MHC-II epitope prediction module (http://tools.iedb.org/mhcii/). The epitopes with low percentile rank were chosen for further consideration. In addition to that, these epitopes were scrutinized for their Th1 type immune response induction propensity following the production of IFN-!. This step was achieved via IFN epitope server (http://crdd.osdd.net/raghava/ifnepitope/). Next, we deployed the ABCPred server (http://crdd.osdd.net/raghava/ abcpred/) for the prediction of B-cell lymphocyte (BCL) epitopes. Finally, all identified epitopes of the CTL, HTL and BCL were examined for their non-toxic/toxic characteristics by deploying the ToxinPred server (http://crdd.osdd.net/ raghava/toxinpred/multi_submit.php).

The top selected epitopes for each cell type were connected sub sequentially to develop the final vaccine construct. Epitope separation was done by utilizing proper linkers/spacer amino acid sequences for the production of an efficient individual function. AAY, GPGPG and KK were the respective linkers employed for joining CTL, HTL and the B-cell epitopes. An immunostimulatory adjuvant (human b-defensin 3) was then added to the sequence towards the N-terminal region. This was followed by another linker sequence (EAAAK) which was added to join the adjuvant with the antigenic epitopes. A single potent vaccine construct was obtained upon subsequent fusion of the adjuvant, epitopes and their specific linkers in a structured way.

The ANTIGENpro tool (http://scratch.proteomics.ics.uci.edu) has been deployed for prediction of antigenicity of the final vaccine candidate. This primordial server can predict the protein antigenicity on the basis of experimentally proven microarray analysis reactivity data. Antigenicity prediction of protein is alignment free and is done based on sequence.

The server is pathogen independent and uses five algorithms of machine learning to make predictions. The non-antigenic/ antigenic character of the protein is decided through SVM (Support Vector Machine) classifier (Magnan et al., 2010) . Next, we made use of two online servers, that is, Algpred (http://crdd.osdd.net/raghava/algpred/) along with AllerTop (https://www.ddg-pharmfac.net/AllerTOP/) to predict the allergenic nature of our newly developed vaccine candidate. The prediction of allergenicity via AlgPred server is based on matching the candidate vaccine sequence with that of the conventional epitope sequences (Saha & Raghava, 2006) while in case of the AllerTop server, machine learning techniques including k nearest neighbors (kNN), Logistic regression (LR), Decision tree (DT), Random forest (RF), Multilayer perceptron (MLP), and naive Bayes (NB) are utilized for prediction. Of all these machine learning techniques, only the kNN method showed five-fold cross validation with 85.3% maximum precision score (Dimitrov et al., 2013) .

Additionally, the designed vaccine candidate was assessed for prediction of its toxicity which was done using ToxinPred server (http://crdd.osdd.net/raghava/toxinpred/). The toxicity prediction through this tool is done based on toxicity of peptides together with the whole protein sequences. It is fed with the input sequence in FASTA format which then works on scoring of position specific basis of the Quantitative Matrix that in return classifies the sequence toxicity (S. Gupta et al., 2015) .

A range of different physicochemical features of this newly designed vaccine candidate were characterized through ProtParam tool (https://web.expasy.org/protparam/) which is present at the ExPASy (Expert Protein Analysis System) bioinformatics resource portal. Various protein parameters together with molecular weight of the peptide, theoretical pI, estimated half-life, extinction coefficient, aliphatic index, instability index and GRAVY index are computed through this server.

We made use of PSIPred 4.0 Protein Sequence Analysis Workbench (http://bioinf.cs.ucl.ac.uk/psipred/) to predict secondary structure features of subunit vaccine candidate. Furthermore, the modeling of tertiary structure was done with Swiss-Model web-server (https://swissmodel.expasy.org/) which has an automated system for generating a 3 D model of a given protein from the amino acid sequence using homology modelling techniques. The only requirement of this server is submission of amino acid sequence of protein under consideration. Selection of a template, its alignment and model building are done automatically by the server. Afterwards, the designed model was evaluated for its quality by using computational approaches like PROCHECK and ProSA (Protein Structural Analysis). The PROCHECK program verifies the structure, figures and the stereochemical behavior of developed model by taking Ramachandran plots as reliable source. The ProSA technique operates by checking energy criteria of the model with respect to a large set of known protein structures of comparable sizes.

Molecular docking of vaccine construct with the human TLR-3 receptor was accomplished through HDOCK web server which was based on intrinsic scoring for determining the binding affinity of our ligand with receptor. The human TLR-3 receptor (PDB ID-3ULV) was downloaded from protein databank.

All-atomistic Molecular dynamics (MD) simulation was accomplished for the vaccine and TLR3 receptor complex using GROMACS 5.0 package which gave details on the dynamics and molecular interactions of the docking complex. Initial optimization of the vaccine-receptor complex was achieved with Schr€ odinger Maestro (Maestro, Schr€ odinger, New York) which was consequently adopted as the starting structure to conduct MD simulation studies. The system was immersed in a cube using simple point charge (SPC) water model and overall charges were neutralized by the addition of counter ions to the system. The energy minimization of this system was carried out by utilizing steepest descent algorithm for 50000 steps until it attained the maximum force < 10.0 kJ/ mol/nm. The electrostatic interactions of a long-range were treated using Particle Mesh Ewald (PME) method. This energy minimized system was further subjected to NVT and NPT equilibration for 100 ps at a constant temperature set to 300 K and 1 bar constant pressure. During equilibration, vrescale, a modified Brendsen thermostat was used for temperature coupling along with Parrinello-Rahman barostat for pressure coupling. The resulting system was subjected to 100 ns MD production run using leap-frog integrator with 2 fs integration time step. The linear constraint solver (LINCS) algorithm was used for restraining the bond lengths. The trajectories were saved every 10 ps and further structural analyses were done using various inbuilt GROMACS scripts. Moreover, the binding energies were evaluated using the molecular mechanics Poisson Boltzmann surface area (MM-PBSA) tool in Gromacs to get detailed insights into molecular interactions of the vaccine candidate and TLR3 complex. It is a widely used method for calculating the binding free energies of the complex from the non-covalent interactions (Kumari et al., 2014) . The g_mmpbsa tool in GROMACS was used after molecular dynamics simulations, the output files obtained were used to post-process binding free energies by the single-trajectory MM-PBSA method. Specifically, for a non-covalent binding interaction in the aqueous phase the binding free energy, DG Binding , is: -

where DG Binding , DG Receptor and DG Ligand are solvation free energies of complex, receptor and ligand, respectively.

A significant way to comprehending the immune system is immune simulation which works by calculating the immune response profile and immunogenicity of probable vaccine construct. An agent-based immune simulator, C-ImmSim (http://kraken.iac.rm.cnr.it/C-IMMSIM/), has been used in our study to carry-out the immune simulation of our designed vaccine construct. The prediction of immune interactions and immune epitopes using C-ImmSim simulator derived from the machine learning techniques and position-specific scoring matrix (PSSM), respectively (Rapin et al., 2010) . Three injections at a gap of four weeks each were given for target product profiling of the newly designed coronavirus vaccine construct. This time period is in accordance with the recommended time interval for most of the currently used vaccines doses that is four weeks. The HLA alleles of host were conserved and time step for injection were computed at 1, 84, and 168 (herein, every time step corresponds to 8 h from real-life, while the first step is injection at time ¼ 0). The volume of simulation was set at 50, while the number of simulation steps was 1000. Remaining simulation specifications were set as default.

The amino acid sequence of target surface spike glycoprotein of SARS-coronavirus was downloaded from National Centre for Biotechonology Information (NCBI). The protein is wellknown to perform key function in interaction, fusion and internalization of virus to host, hence conferring the viral pathogenicity. This sequence was cropped and residues numbered 331 to 524 were retained which corresponds to the RBD of the protein of our interest (Tai et al., 2020) . The viral sequence was checked for homology with the human proteome through the pBLAST server. The results showed no similarity with Homo sapiens. The similarity between pathogen and host should not be more than 40% to avoid chances of cross-reactivity. Thus, the results indicated that the antigenic viral proteins used in designing vaccine candidate will not cross-react upon administration in the host. Next, the protein antigenic propensity average was calculated using the antigenic peptide prediction tool from the Immunomedicine group. It was found to be 1.03, indicating the highly antigenic nature of the peptide which can be further utilized for the prediction of epitopes for developing a multi-epitope vaccine. Figure 1 gives a schematically represents the methodology used for constructing the multi-epitope subunit vaccine.

The Cytotoxic T-lymphocytes (CTLs) have a critical role in inducing the MHC-I based cellular immune responses. These cells neutralize the viral particle infected or damaged host cells by producing the proteins such as perforins, granzymes, etc (Neurath, 2008) . The CTL epitopes of the screened protein were predicted deploying the NetCTL1.2 server (http:// www.cbs.dtu.dk/services/NetCTL/). We considered A2, A3, and B7 supertypes for the CTL epitopes determination as these supertypes majorly cover 88.3% of the total ethnic population (Sidney et al., 2008) . The top ten epitopes obtained to be used for further vaccine designing are shown in Supplementary Table 1 . The Helper T-lymphocytes (HTLs) have a key role in adaptive immunity. These activate the B-cells along with the cytotoxic T cells for production antibody and eventually killing infected/damaged cells (Alberts et al., 2002) . The HTL epitopes for the selected protein were calculated using the prediction tool for MHC-II epitope. The identified epitopes had lower percentile ranks. The HTL epitopes were additionally checked by their potency to induce IFN-! production ability using IFN epitope server. The uppermost immunogenic epitopes of 15-mer were shortlisted based on basis of the percentile rank for each. Along with that, IFN epitope server output gave a positive score to all these epitopes (Supplementary Table 2 ).

Humoral immunity is mainly driven by B-cells which are responsible for the production of antibodies for recognizing and eliminating the antigens (Janeway et al., 2001) . This necessitated the need for predicting the B-cell epitopes before designing our vaccine candidate. The ABCPred server (http://crdd.osdd.net/raghava/abcpred/) was deployed for prediction of these B-cell epitopes (Supplementary Table 3 

Once the epitopes of different immune cells were predicted, we combined them using specific linker sequences in order to create a single vaccine construct. Each epitope of multiepitope subunit vaccine should be should be having the capacity of inducing a potential immune response. The linkers play critical role in locating the epitope and are also key for conferring immunogenicity to the whole vaccine construct. These linker sequences are accountable for the separation of different domains combined in a unique single protein sequence without disrupting the functionality of that protein. Along with this, linkers are instrumental in maintaining stability and flexibility in the protein complex (Chen et al., 2013) .

In the current work, we employed diverse varieties of linker sequences to connect the epitopes. The CTL epitopes were linked using "AAY" (Ojha et al., 2019) derived from the proteasomal cleavage site to eliminated the genesis of newer epitope, HTL epitopes were connected by "GPGPG" linkers (Livingston et al., 2002) used for improving the flexibility of protein and epitopes. The "KK" linkers (Yano et al., 2005) , for the enhancement of the proteasome processing, were used for linking B-cell epitopes. In addition to these, another linker known to form helices, EAAAK, was used for designing the next-generation subunit vaccine.

A human beta-defensin-3 (45 amino acids) adjuvant was attached at the N-terminal location of the designed vaccine to confer immunogenicity to the vaccine construct. The human beta defensin 3 is a TLR-3 agonist capable of eliciting enhancement in the strength of T-cell and B-cell immune responses (Gupta et al., 2020; Mohan et al., 2013) . A total of 10 CTL, 10 HTL, and 10 B-cell epitopes derivative of the receptor binding protein were combined for designing a subunit vaccine (487 aa residues) against SARS-CoV-2. Presently, we drafted the candidate vaccine molecule in the following sequence: "adjuvant-CTL-HTL-B-cell epitopes". Therein, the initial 45 amino acids represent an adjuvant, followed by a 6 residue linker (EAAAK). The region from 51 to 167 amino acids was represented by CTL epitopes while 168-367 and 368-487 amino acid residues were covered by HTL and B-cell epitopes, respectively (Figure 2 ).

An imperative consideration to consider in a vaccine candidate is its immunogenicity, that is, it should have the ability to induce an immune response (cell-mediated/humoral) against the targeted pathogenic organism. Computational findings from our work indicated that the antigenic characteristics of our candidate vaccine construct had a probability score of 0.854 as predicted by the ANTIGENPro server. The preferred average score of antigenic probability of the candidate to be used as vaccine is !0.8 (Pandey et al., 2018) . Thus our obtained score for antigenic probability points towards the potentially antigenic nature of the newly designed subunit vaccine candidate. The antigenic probability for the sequence orientation used in the present study was higher than the cases where the epitopes were differently placed.

The AlgPred server was deployed for allergenicity assessment of the vaccine. The basis for data analysis through this server includes following six approaches:(i) IgE mapping, (ii) MEME/Mast motif, (iii, iv) SVM modules for both amino acid and dipeptide composition, (v) BLAST search on ARPs, and (vi) Fusion technique/hybrid approach involving all the parameters. Of the six approaches, the MEME/Mast motif, BLAST search on ARPs, and IgE mapping approaches indicate that the vaccine construct is non-allergenic in nature, whereas the SVM modules working with both amino acid and dipeptide compositions and hybrid approach recognized allergenicity of the vaccine construct. Therefore, we investigated further and performed the allergenicity analysis using the AllerTop server. The results obtained from this tool validated the non-allergen nature of the vaccine. Thus, analyzing the results from both the servers, we concluded that the vaccine candidate was most likely a non-allergen and it would not produce any allergic reaction during the vaccination protocol.

The newly designed vaccine candidate construct is comprised of 487 amino acids along with an approximate molecular weight of 53.64 kDa. Other physiochemical features of the vaccine candidate were computed to determine the vaccine characteristics. The instability index which should be less than 40 was 17.31, demonstrating the stability of our vaccine construct. The significantly basic nature of the vaccine was determined from the theoretical pI which was 9.84. Next, the protein half-life was computed on the basis of Nterminal residue. The data obtained showed that the vaccines half-life of >30 h in mammalian reticulocytes, suggestive of in vivo stability of the vaccine candidate.

The vaccine was found to be moderately thermostable in nature with the computed aliphatic index score of 57.97. The higher aliphatic index was suggestive of its thermostability. The estimated GRAVY (grand average of hydropathicity) score was À0.415. The negative GRAVY score indicated towards hydrophilic behavior having good interaction with the water molecules. Further evaluated parameters have been tabulated in Table 1 .

The amino acid sequence of our vaccine candidate was used for evaluating the secondary structural features of the model using PSIPred 4 server. Figure 3 illustrates the predicted results for secondary structure. Since the immunoinformatic approach is heavily dependent on the tertiary structural features of the protein, the 3 D-structure prediction is a crucial step in designing a vaccine. The homology modeling of protein structure is a commonly used technique to generate 3 D models for protein in the absence of experimental structures. Thus, the tertiary structure of the vaccine was estimated via employing an online tool SWISS-MODEL which is a fully automated server with a user-friendly web interface that generates reliable models (Figure 4(A) ). Various tools such as the PROCHECK and ProSA were used for validating the 3 D model by checking its structural integrity. Figure 4(B,D) gives information on the quality of our vaccine model generated by homology modelling. Ramachandran plot (Figure 4(D) ) statistics were generated using the PROCHECK server. The output data generated for the vaccine candidate was indicative of a reasonable quality model of the construct with 76.7% (66 aa) of all the residues in the most favored regions and 22.1% (19 aa) were present in additionally permitted regions. Only 1.2% (1 aa) were occurring in the disallowed plot region.

Additionally, the ProSA server was employed to compute the Z-score and generate an energy plot for the 3 D model. The Z-score (dark spot) value was À1.84, which is well within the range defined for known proteins determined through NMR (dark blue) and X-ray (light blue) (Figure 4(E) ), thus predicting a good quality model. The energy plot for the 3 D model has been shown in Figure 4 (F) which illustrates the model quality via plotting energies as a function on the xaxis against the residue sequence position on the y-axis. Low/negative energy values of the residues indicated that the model was in the acceptable range. CoV-2) . The designed vaccine candidate is in the sequential manner like-Adjuvant-PADRE sequence-CTL epitopes-HTL epitopes-B-cell epitopes. In the vaccine sequence, the initial 45 amino acids are occupied by an adjuvant (Black color), followed by next 6 amino acids of EAAAK linker (Red color), then from number 51-167 amino acids represent the CTL epitopes (Blue color) whereas 168-367 and 368-487 are covered by HTL epitopes (Green color) and B-cell epitopes (Brown color), respectively. All the linkers have been represented in red color. 

Docking is a molecular modeling method used for predicting the favored ligand molecule orientation to the target protein for building a stable adduct. Molecular docking studies were conducted to approximately calculate binding energy of our vaccine construct to the TLR-3 receptor ( Figure 5(A,B) ). The binding score of the complex was À217.8 kJ/mol. Furthermore, the vaccine-receptor complex was subjected to molecular dynamics simulation studies. 

An all-atomistic MD simulation study was performed for 100 ns to evaluate the binding mode, molecular dynamics, and structural stability of our vaccine-TLR-3 docked complex ( Figure 5 ). The interactive parameters between vaccine and TLR-3 were determined at atomic-level, and trajectory analysis was performed to examine the variations in root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), and hydrogen bonding over a period of simulation. The convergence of our complex system was studied using RMSD analysis of all the C a atoms present in the simulated complex. The RMSD plot for vaccine-receptor complex showed high energy fluctuations up to 40 ns ( Figure 5(C) ), wherein the resulting RMSD was found to be in the range of 0.2-1.2 nm for the backbone. After 40 ns, the energy was observed to be consistent until the end of the simulation period (100 ns), thus resulting in RMSD values from 1.0-1.2 nm. The magnitude of fluctuations in the simulation trajectories was the basis for our studies. A stable RMSD plot provides for a firm base to perform advanced investigations. The RMSF values for the residual fluctuations were investigated in order to gain insights concerning the alterations in the dynamic behavior of atoms of the docked complex ( Figure 5(D) ).

We further analyzed the values for Rg to further study the comprehensive aspect of protein. The Rg is usually represented as the mass weight root mean square distance for set of atoms distanced based on their common center of mass. Figure 5 (E) shows the Rg plot for the C a atom of protein versus time at a given temperature (300 K). The hydrogen bonds were chiefly accountable for the maintenance of a stable configuration of the vaccine-receptor protein complex. Therefore, NH bond analysis of the vaccine-TLR3 complex was performed to determine the connection between NH bond formation and complex flexibility over a simulation time period of 100 ns. The number of intermolecular hydrogen bonds increased with time ( Figure 5(F) ). The results of H-bond plot were suggestive of reduced flexibility, increased strength of binding, and improved vaccine-receptor stability.

The binding energies were computed using the MM-PBSA tool to get detailed insights into molecular interactions of the vaccine candidate and TLR3 complex. The calculations were done for the equilibrium phase between 90-100 ns from the 100 ns output file of MD simulations. The interactions were found to be the most favorable in terms of the sum of various binding energies. The electrostatic energy distribution in molecular mechanics was derived through the high Vander Waal energy (VDW) _-738.16 ± 38.22 kJ/mol, electrostatic energy À8413.94 ± 346.10 kJ/mol and total binding energy of -6592.13 ± 261.07 kJ/mol, which are considered as most significant contributors in the binding of ligand to the target protein. The polar solvation energy (PBSOL) was calculated to be 2656.29 ± 143.86 kJ/mol and SASA energy was À96.33 ± 4.69 kJ/mol which added a vital contribution to the stable binding of the ligand to target protein.

The designed vaccine construct was finally examined for its ability to generate immune response in the body. This study was performed with an online in-silico server for immune simulation (C-ImmSim). The results obtained were in accordance with a typical immune response generated in the body ( Figure 6 ). Increased levels of IgM with every injection indicated the activation of primary immune response. All the associated immune responses (secondary and tertiary) were identified by looking at the elevated levels of IgM, IgG þ IgM, and IgG1 þ IgG2 antibodies which were followed by subsequent reduction in the concentration of antigen ( Figure  6(A) ). Additionally, long lasting B-cell isotypes were observed, thus indicating a potential to possibly switch the isotypes and memory formation ( Figure 6 (D,E)). Likewise, there were increased levels of the T-helper (T H ) and cytotoxic (T C ) cell populations with corresponding development of memory ( Figure 6 (G,J), respectively). The macrophage activity was found increasing (Figure 6(N) ) whereas the activity levels of dendritic cells were consistent throughout the period of antigenic exposure (Figure 6 (M)). The increased levels of interferon (IFN-c) and interleukin (IL-2) production were indicative of efficient immunoglobulin production, thus implicating towards the generation of a humoral response. Furthermore, various other attributes were analyzed including plasma B lymphocytes (PLB) cell population, CD4þ T regulatory cell population (TR) per state, natural killer (NK) cell population, along with epithelial cells population (EP) per state ( Figure  6 ), wherein a good immune response was observed in the cases of all these cell populations. Overall, the results demonstrated the development of immune memory along with subsequent antigen clearing potential at consecutive exposures. 

The inhibition of viral infections in humans can be best achieved through vaccination. The recent COVID-19 outbreak which has claimed lives worldwide has spurred the need for urgently developing an effective vaccine for treating this deadly infectious disease. Considering the relative ease of application in our investigative experimentations, we made use of the available immunoinformatic techniques for prospective vaccine candidate designing against the SARS-CoV2 antigen. A multi-epitope subunit peptide was constructed with addition of a suitable adjuvant (beta-defensin-3) at the N-terminal to ensure proper delivery along with proficient epitope immuno-processing. This was followed by connecting the CTL epitopes, HTL epitopes, and B-cell epitopes via suitable linker sequences. The effectiveness of the vaccine construct was assessed by analyzing the physicochemical properties through available online servers. Molecular docking interaction studies were performed to evaluate the specificity of binding for the vaccine candidate to the human TLR3-immune cell receptors. A low energy score indicated an efficient interaction that was both strong and specific between the vaccine candidate and interactive cavity of the human TLR-3 receptor. The MD simulations result showed stable and prolonged binding with non-significant deviation in RMS value. Moreover, a strong immunological response was obtained with the T-cell and B-cell mediated humoral response. These encouraging findings through computational analysis are indicative that this newly designed multiple-epitope based subunit vaccine candidate has a high probability of being safe and effective in protecting humans from the SARS-CoV-2 infections. However, this proposed vaccine will need experimental validation in suitable animal models to determine its real potency in inducing immunity against the coronavirus. Prophylactic measures for self-defense against the virus include preventive use of animals or avian in human diet. Also, the maintenance of balance in the ecosystem should be reinvestigated.

In summary, the work done highlights a novel vaccine construct which can elicit a strong immune response against the current COVID-19 outbreak to stop the pandemic. This can serve as a prototype for developing vaccines against the other emerging infectious diseases.

Lymphocytes and the cellular basis of adaptive immunity

Mechanisms of coronavirus cell entry mediated by the viral spike protein

Middle East Respiratory Syndrome Coronavirus: Another Zoonotic Betacoronavirus Causing SARS-Like Disease

Fusion protein linkers: Property, design and functionality

AllerTOP -a server for in silico prediction of allergens

Biochemical aspects of coronavirus replication and virus-host interaction

Taxonomy, Classification and Nomenclature of Viruses. Encyclopedia of Virology

Nidovirales: Evolving the largest RNA virus genome

Receptor-ligand based molecular interaction to discover adjuvant for immune cell TLRs to develop next-generation vaccine

Peptide Toxicity Prediction BT -Computational Peptidology

A, Open Source Drug Discovery Consortium (2014). g_mmpbsa-a GROMACS tool for high-throughput MM-PBSA calculations

Porcine epidemic diarrhea virus: An emerging and reemerging epizootic swine virus

Structure, function, and evolution of coronavirus spike proteins

Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus

Receptor and viral determinants of SARS-coronavirus adaptation to human ACE2

A Rational Strategy to Design Multiepitope Immunogens Based on Multiple Th Lymphocyte Epitopes

COVID-19: A novel zoonotic disease caused by a coronavirus from China: What we know and what we don't

High-throughput prediction of protein antigenicity using protein microarray data

Viruses and bacteria in the etiology of the common cold

Modulation of HIV peptide antigen specific cellular immune response by synthetic aand b-defensin peptides

Immune Response to Viruses: Antibody-Mediated Immunity. Encyclopedia of Virology

Strategic Development of a Next-Generation Multi-Epitope Vaccine To Prevent Nipah Virus Zoonotic Infection

Human coronaviruses associated with upper respiratory tract infections in three rural areas of Ghana

Editorial: Epitope Discovery and Synthetic Vaccine Design

Novel Immunoinformatics Approaches to Design Multi-epitope Subunit Vaccine for Malaria by Investigating Anopheles Salivary Protein

Coronaviruses post-SARS: Update on replication and pathogenesis

Computational immunology meets bioinformatics: The use of prediction tools for molecular binding in the simulation of the immune system

AlgPred: Prediction of allergenic proteins and mapping of IgE epitopes

HLA class I supertypes: A revised and updated classification

Cytokine storm induced by SARS-CoV-2

Characterization of the receptor-binding domain (RBD) of 2019 novel coronavirus: Implication for development of RBD protein as a viral attachment inhibitor and vaccine

Coronavirus disease (COVID-19)

An ingenious design for peptide vaccines

Genetic analysis of determinants for spike glycoprotein assembly into murine coronavirus virions: Distinct roles for charge-rich and cysteine-rich regions of the endodomain

PRS is thankful to Indian Council of Medical Research (ICMR), India for the senior research fellowship. JG is contended to DBT for fellowship. PS is grateful to TERI School of Advanced Studies for research and management. AG is thankful to the University Grant Commission, India for the faculty recharge position.

The authors declare that they have no conflicts of interest.

PRS, JG, PS, VKP and AG conceptualized the study and designed the methods and experimental setup. PRS conducted computational experiments, analyzed the results and drafted the manuscript. JG helped in manuscript compilation. All authors reviewed and approved the manuscript.

Vijay Kumar Prajapati http://orcid.org/0000-0001-6510-0596 Abhinav Grover http://orcid.org/0000-0002-3296-7860