key: cord-0774822-kpmwn0c2 authors: Abraham Peele, K.; Srihansa, T.; Krupanidhi, S.; Vijaya Sai, A.; Venkateswarulu, T. C. title: Design of multi-epitope vaccine candidate against SARS-CoV-2: a in-silico study date: 2020-06-01 journal: J Biomol Struct Dyn DOI: 10.1080/07391102.2020.1770127 sha: defec48c714bf474bf4c1b74c6d6128035c5e241 doc_id: 774822 cord_uid: kpmwn0c2 The best therapeutic strategy to find an effective vaccine against SARS-CoV-2 is to explore the target structural protein. In the present study, a novel multi-epitope vaccine is designed using in silico tools that potentially trigger both CD4 and CD8 T-cell immune responses against the novel Coronavirus. The vaccine candidate was designed using B and T-cell epitopes that can act as an immunogen and elicits immune response in the host system. NCBI was used for the retrieval of surface spike glycoprotein, of novel corona virus (SARS-CoV-2) strains. VaxiJen server screens the most important immunogen of all the proteins and IEDB server gives the prediction and analysis of B and T cell epitopes. Final vaccine construct was designed in silico composed of 425 amino acids including the 50S ribosomal protein adjuvant and the construct was computationally validated in terms of antigenicity, allergenicity and stability on considering all critical parameters into consideration. The results subjected to the modeling and docking studies of vaccine were validated. Molecular docking study revealed the protein-protein binding interactions between the vaccine construct and TLR-3 immune receptor. The MD simulations confirmed stability of the binding pose. The immune simulation results showed significant response for immune cells. The findings of the study confirmed that the final vaccine construct of chimeric peptide could able to enhance the immune response against nCoV-19. COVID-19 pandemic is result of the infection caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and it attacks the vital organs of body and targets pneumocytes in lungs which leading to fatal respiratory distress (Elfiky, 2020; Galante et al., 2016; Joshi et al., 2020; Tse et al., 2004; Yi et al., 2020) . SARS, MERS, and SARS-CoV-2 caused diseases are characterized by a lower respiratory ailment like bronchitis, bronchiolitis, and pneumonia (Bogoch et al., 2020; Elfiky & Azzam, 2020) . Coronaviruses are enveloped, single stranded RNA viruses which bears club shaped glycoproteins on their surface Sarma et al., 2020) . Corona virus is contagious and spreads through inhalation, ingestion of viral droplets resulting from coughing and sneezing. The coronavirus genome is comprised of $30000 nucleotides and it encodes with four structural proteins, Nucleocapsid (N) protein, Membrane (M) protein, Spike (S) protein and Envelop (E) protein and several non-structural proteins (nsp) (Boopathi et al., 2020; Gupta et al., 2020; Hasan et al., 2020; . The non-structural protein 16 (nsp16) or 2 0 -OMTase is the crucial protein responsible for viral replication and expression in host cells . The SARS-CoV-2 has a surface spike like glycoprotein (S-Spike), it binds specifically to angiotensinconverting enzyme 2 (ACE2) to access host type-2 pneumocyte cells that could bind to human cells (Hoffmann et al., 2020; Verdecchia et al., 2020) . The enveloped virus covered with spike surface glycoprotein and after binding the host cell is engulfed in to the cell and releases positive sense single stranded RNA in to type-2 pneumocyte. RNA dependent RNA polymerases and proteinases makeup the components of viral proteins like nucleocapsid spike proteins and enzymes necessary for viral replication and damage to host cells. Type-2 pneumocytes that produces a surfactant molecule and decreases surface tension in the alveoli and reduce the collapsing pressure (Belouzard et al., 2012; Weiss & Navas-Martin, 2005) . However, at present researchers are aiming to uncover this spike protein processing proteases for early drug development using various approved inhibitor drugs/other compounds (Aanouz et al., 2020; Elmezayen et al., 2020; Enmozhi et al., 2020; Islam et al., 2020; Muralidharan et al., 2020; Pant et al., 2020; Sinha et al., 2020; Wahedi et al., 2020) . The viral proteins encounters immune cells and release specific cytokines leads to vasodilation, capillary permeability, aleveolar edema and finally increase the collapsing pressure to burst out the pneumocyte (van de Veerdonk et al., 2020) . SARS-CoV-2 surface spike (S) protein contains two subunits S1 and S2, of which S1 is sole responsible for host cell receptor and S2 harbors the membrane fusion machinery. Spike protein from SARS-CoV-2 shares the high structural similarity with SARS-CoV spike (Weiss & Navas-Martin, 2005) . The vaccine development against SARS-CoV-2 spike protein is an important approach and hence, the present study is focused on epitope prediction analysis for construction of vaccine candidate by computational methods. The complete amino acid sequence of spike glycoprotein (NCBI Accession id ¼ "QHD43416.1) of SARS-CoV-2 strain was retrieved from NCBI. Cytotoxic T-cell lymphocyte (CTL) epitopes were predicted for spike glycoprotein using artificial neural network algorithm based online server NetCTL 1.2 which predicts the MHCclass I binding, and then followed by submitting predicted NetCTL generated results to VaxiJen v2.0 (http://www.ddgpharmfac.net/vaxijen/VaxiJen/VaxiJen.html), ToxinPred (http:// crdd.osdd.net/raghava/toxinpred/) servers to predict protective nontoxic antigens (Larsen et al., 2007) . After the screening of epitopes with VaxiJen v2.0 and ToxinPred servers, the resultant epitopes were subjected to immunogenicity prediction using IEDB server (https://www.iedb.org/). IEDB server was used to predict the MHC-II restricted epitopes as it uses special patterns for HLA-DRB1 Ã 01:01, HLA-DPA1 Ã 01/DPB1 Ã 04:01, HLA-DQA1 Ã 03:01/DQB1 Ã 03:02 alleles, further, T-helper 1-type immune response activation and IFNc production was predicted using IFNepitope server (http:// crdd.osdd.net/raghava/ifnepitope/). Toxicity was predicted using ToxinPred server (http://crdd.osdd.net/raghava/toxinpred/) and BCPRED 2.0 online server was used to predict the linear B-cell epitopes of spike protein (http://crdd.osdd. net/raghava/bcepred/). The results of linear B-cell epitopes and HTL epitopes of overlapping regions were assembled and considered as final predicted epitopes. The potential non toxic and probable antigenic vaccine was constructed using selected epitopes. The linear B-cell and HTL epitopes were joined with GPGPG linker peptides and CTL epitopes were joined using AAY linker. The N-terminal position of the vaccine construct was linked with the sequence of 125 amino acid residue 50S ribosomal L7/L12 peptide which acts as an adjuvant and C-terminal portion was linked with HHHHHH (6HIS) linker. The basic property of allergenicity was assessed using an online server, AllerTOP v. 2.0 (https://www. ddg-pharmfac.net/AllerTOP/) and then submitted to VaxiJen server to find out whether the construct vaccine could be a probable antigen to elicit immune response. PSIPRED web tool was employed for finding the secondary structure analysis (http://bioinf.cs.ucl.ac.uk/psipred/). Homology modeling of the vaccine protein tertiary structure was performed using fully automated protein modeling I-TASSER server (https:// zhanglab.ccmb.med.umich.edu/I-TASSER/) and the best model was selected and then optimized using SPDB viewer (https:// spdbv.vital-it.ch/). Loop refinement was done using ModLoop (https://modbase.compbio.ucsf.edu/modloop/). Ramachandran plot and ERRAT server (https://servicesn.mbi.ucla.edu/ERRAT/) analysis were performed for further validation study. The designed vaccine candidate was subjected to docking using GRAMM-X Simulation web server (http://vakser.compbio.ku.edu/ resources/gramm/grammx/) for docking vaccine protein model with TLR-3 (PDB ID: 1ziw) and interaction were visualized using LIGPLUS 1.2 software (https://www.ebi.ac.uk/thornton-srv/ software/LIGPLOT/). The protein & receptor complex was subjected to MD simulations. The MD simulations were done by GROMACS 2018 package to carry out 20 ns simulations using OPLS force field. The TIP3P water model was selected for solvating complexes followed by the addition of ions to neutralize. Periodic boundary conditions were used and Equilibration of the system was done using NVT and NPT ensemble for 100 ps. The trajectory was set to be generated every 2 fs and save every 2 ps. The protein-protein complex result was then analyzed (Enayatkhani et al., 2020) . The immune response profile of vaccine construct was recorded by in silco method C-ImmSim, online simulation server (http://150.146.2.1/C-IMMSIM/index.php). The C-ImmSim model describes both humoral and cellular response of a mammalian immune system against vaccine construct. The target product profile of a prophylactic vaccine, three injections were given at different intervals of four weeks. The time step of simulation corresponds few hours of real life and simulation was performed with default parameters. The sequence of injections is Ag1, Ag2, Ag3 were administered four weeks apart. The simulation volume and simulation steps were set at 1000, (random seed ¼ 12345 with an injection of vaccine containing no LPS. The physiochemical properties of the vaccine construct was assessed using online web tool the ProtParam (https://web. expasy.org/protparam/), where as the solubility, allergenicity and probable antigenic prediction were performed using protsol, AllerTOP v. 2.0 and VaxiJen v2.0 online web servers. J-CAT tool (http://www.jcat.de/) is used for codon optimization of vaccine construct and E.coli (K12) strain is selected as source organism. SnapGene software (https://www.snapgene. com/try-snapgene/) was used for in silico cloning of vaccine construct in to pET-28a vector. The amino acid sequence was used to predict the possible probable antigenic epitopes of linear B-cell, HTL and CTL epitopes for designing the multi-epitope vaccine. The vaccine construct consisted of 425 amino acid residues derived from different peptide sequences. CTL epitopes of 9-mer lengths were predicted using NetCTL1.2 (Table 1) . Based on high binding affinity score, the results were submitted to VaxiJen v2.0 and predicted the 16 protective probable antigens. The non antigenic epitopes were removed and subjected to predict the toxicity using ToxinPred and after removing two toxin epitopes 14 non-allergenic epitopes were selected using toxinpred and, the IEDB immunogenicity server produced the results of seven epitopes and were given in the Table 2 . The predicted probable antigenic HTL epitopes were selected for further screening of toxigenicity prediction using vaxigen 2.0 server and 13 HTL epitopes were selected and further classified as non-toxins using Toxinpred server (Tables 3 and 4 ). The final HTL epitopes were selected as a result of IFN-c inducing epitopes (Table 5 ). The linear B-cell epitopes were used in vaccine construct as overlapping B-cell and T-cell epitopes. The final multi-epitope subunit vaccine model was generated using I-TASSER server. The top finest threading templates for building the protein models were selected (1rqu, 6f0k, 1dd3, 3j4a, 2ftc) and based on high c-score value -0.65 and an estimated TM score of 0.63, model protein was selected. Energy minimization and refinement of the modeled structure was carried out with SPBD viewer and loop regions were identified and refined. Final structure was checked with Ramachandran plot and showed 99% of the residues are in favorable region (Figure 1a) and ERRAT server showed 81% quality score (not shown). The results of ProSA-web obtained for vaccine construct was provided the z-score value of À0.79, as the major parts of the energy plot with N-terminal region and C-terminal region showed highly positive energy values (Figure 1b) . PSIPRED produced the secondary structural information of the vaccine construct (Figure 2a) . The Prosol server provided the solubility prediction calculations and average of all residues produced the value of 0.46 and indicated good solubility of the vaccine construct (Figure 2b) . The physicochemical features of the vaccine were analyzed by ProtParam tool and were given the Molecular weight of 44 kda, therotical pI caluculated value of 5.16. The instability index (II) is computed to be 16.88 and protein classified as stable and the aliphatic index is 94.92%. Grand average of hydropathicity (GRAVY): 0.262. The estimated half-life is: 30 h (mammalian reticulocytes, in vitro). The overall prediction of the vaccine construct is found to be probable antigen with a score of 0.5848 generated by vaxigen 2.0 and AllerTOP 2.0 has classified the construct to be non-allergen and predicted that the nearest protein to be Scarecrow 1 in Oryza sativa ( Figure 3 ). Optimized codons sequence length was found to be 1175 nucleotides and average GC content was 59.1%. The pET28a (þ) vector was used to clone the vaccine construct DNA sequence using SnapGene software (Figure 4 ). The adjuvant which is a 50S ribosomal protein has the ability to stimulate TLR3, Molecular docking using the GRAMXX server (http://vakser.compbio.ku.edu/resources/gramm/grammx/) produced the best protein-protein docking model complex with binding score for molecular docking produced the best structure with global energy -35.98 and interactions and attractive vanderwaals (-32.28) was selected and DIMPLOT of LIGPLUS version 1.2 visualized the interaction between chain A-TLR-3 protein and chain B -Vaccine construct (Figure 5a-c) . The RMSD values of protein-ligand complexes were recorded from 0 to 20 ns. The RMSD values steadily increased from 0 to 5 ns and reached a stable state throughout the simulation. The average RMSD values of the complex were found to be 0.27 nm (Figure 6 ). C-ImmSim considers the successive and successful immune responses of the state of the cell and model the memory of immune cells by a mechanism that increases their half-life. The result of process is that few cells increase their half-life considerably and live longer than other cells. ImmSim server immune simulation results confirmed consistency with actual immune responses. High levels of IgM indicated the primary response. Furthermore, an increase in the B-cell population was characterized by an increase in the expression of immunoglobulins which resulted in a decrease in the concentration of the antigen. Also, there is a consistent rise in Th (helper) cell population with memory development (Figure 7a-c) . It was also observed that the production of IFN-c was stimulated after immunization ( Figure 7d ). The results clearly explained the T cell population was highly responsive as the memory developed and all other immune cell population shown to be consistent. The vaccine candidate against spike viral surface glycoprotein of SARS-CoV-2 was designed by in silico methods. The epitopes predicted with different web servers and adjuvant linkers were used to construct a potent antigenic, nonallergenic vaccine that could elicit strong immune response against SARS-CoV-2. Docking analysis provided the validation in the form of affinity between two molecules (TLR-3 and vaccine) and stability of complex was supported by MD simulations. The in silico immune simulation confirmed immune cell response against antigen clearance rate. The computational cloning by SnapGene confirmed the strong expression of proteins. However, the experimental validation could be essential to ensure to vaccine construct efficacy against COVID-19. Moroccan Medicinal plants as inhibitors of COVID-19: Computational investigations Mechanisms of coronavirus cell entry mediated by the viral spike protein Potential for global spread of a novel coronavirus from China Novel 2019 coronavirus structure, mechanism of action, antiviral drug promises and rule out against its treatment SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: An in silico perspective Novel Guanosine Derivatives against MERS CoV polymerase: An in silico perspective Drug repurposing for coronavirus (COVID-19): In silico screening of known drugs against coronavirus 3CL hydrolase and protease enzymes Reverse vaccinology approach to design a novel multi-epitope vaccine candidate against COVID-19: An in silico study Andrographolide as a potential inhibitor of SARS-CoV-2 main protease:" An in silico approach Coronavirus NL63-induced adult respiratory distress syndrome Figure 7. C-ImmSim server prediction results of immune response after administering vaccine construct; (a) Antigen and immunoglobulins; (b) B-lymphocytes cell population; (c) CD4 þ helper T cells population per state; (d) Induced levels of the cytokine and Simpson index In-silico approaches to detect inhibitors of the human severe acute respiratory syndrome coronavirus envelope protein ion channel A review on the cleavage priming of the spike protein on coronavirus by angiotensin-converting enzyme-2 and furin SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor A molecular modeling approach to identify effective antiviral phytochemicals against the main protease of SARS-CoV-2 Discovery of potential multi-target-directed ligands by targeting host-specific SARS-CoV-2 structurally conserved main protease Targeting SARS-Cov-2: A systematic drug repurposing approach to identify promising inhibitors against 3C-like proteinase and 2'-O-ribose methyltransferase Identification of chymotrypsin-like protease inhibitors of SARS-CoV-2 via integrated computational approach Large-scale validation of methods for cytotoxic T-lymphocyte epitope prediction Computational studies of drug repurposing and synergism of lopinavir, oseltamivir and ritonavir binding with SARS-CoV-2 Protease against COVID-19 Peptide-like and small-molecule inhibitors against Covid-19 Drug targets for corona virus: A systematic review In-silico homology assisted identification of inhibitor of RNA binding against 2019-nCoV N-protein (N terminal domain) An in-silico evaluation of different saikosaponins for their potency against SARS-CoV-2 using NSP15 and fusion spike glycoprotein as targets Pulmonary pathological features in coronavirus associated severe acute respiratory syndrome (SARS) Kinins and cytokines in COVID-19: A comprehensive pathophysiological approach The pivotal link between ACE2 deficiency and SARS-CoV-2 infection Stilbene-based natural compounds as promising drug candidates against COVID-19 Coronavirus pathogenesis and the emerging pathogen severe acute respiratory syndrome coronavirus COVID-19: What has been learned and to be learned about the novel coronavirus disease The authors acknowledge to VFSTR (Deemed to be university) and DST-FIST (LSI-576/2013) networking facility to carry out this work. No potential conflict of interest is reported by the authors.