key: cord-319689-33h22ikl authors: Srivastava, Sukrit; Kamthania, Mohit; Singh, Soni; Saxena, Ajay K; Sharma, Nishi title: Structural basis of development of multi-epitope vaccine against Middle East respiratory syndrome using in silico approach date: 2018-11-21 journal: Infect Drug Resist DOI: 10.2147/idr.s175114 sha: doc_id: 319689 cord_uid: 33h22ikl BACKGROUND: Middle East respiratory syndrome (MERS) is caused by MERS coronavirus (MERS-CoV). Thus far, MERS outbreaks have been reported from Saudi Arabia (2013 and 2014) and South Korea (2015). No specific vaccine has yet been reported against MERS. PURPOSE: To address the urgent need for an MERS vaccine, in the present study, we have designed two multi-epitope vaccines (MEVs) against MERS utilizing several in silico methods and tools. METHODS: The design of both the multi-epitope vaccines (MEVs) are composed of cytotoxic T lymphocyte (CTL) and helper T lymphocyte (HTL) epitopes, screened form thirteen different proteins of MERS-CoV. Both the MEVs also carry potential B-cell linear epitope regions, B-cell discontinuous epitopes as well as interferon-γ-inducing epitopes. Human β-defensin-2 and β-defensin-3 were used as adjuvants to enhance the immune response of MEVs. To design the MEVs, short peptide molecular linkers were utilized to link screened most potential CTL epitopes, HTL epitopes and the adjuvants. Tertiary models for both the MEVs were generated, refined, and further studied for their molecular interaction with toll-like receptor 3. The cDNAs of both MEVs were generated and analyzed in silico for their expression in a mammalian host cell line (human). RESULTS: Screened CTL and HTL epitopes were found to have high propensity for stable molecular interaction with HLA alleles molecules. CTL epitopes were also found to have favorable molecular interaction within the cavity of transporter associated with antigen processing. The selected CTL and HTL epitopes jointly cover upto 94.0% of worldwide human population. Both the CTL and HTL MEVs molecular models have shown to have stable binding and complex formation propensity with toll-like receptor 3. The cDNA analysis of both the MEVs have shown high expression tendency in mammalian host cell line (human). CONCLUSION: After multistage in silico analysis, both the MEVs are predicted to elicit humoral as well as cell mediated immune response. Epitopes of the designed MEVs are predicted to cover large human population worldwide. Hence both the designed MEVs could be tried in vivo as potential vaccine candidates against MERS. Middle East respiratory syndrome (MERS) is a respiratory disease caused by Middle East respiratory syndrome coronavirus (MERS-CoV). MERS involves high fever, cough, difficulty in breathing, chills, chest pain, body aches, sore throat, diarrhea, nausea/vomiting, running nose, renal failure, and pneumonia. The first case of MERS in humans was reported in Saudi Arabia in 2012. After 2012, within a span of only 4 years, MERS-CoV infection was reported from 27 countries. 1 Three MERS outbreaks have already occurred in Saudi Arabia (2013 and 2014) and South Korea (2015). [2] [3] [4] The outbreak reported in South Korea also involved spread of MERS through hospital-to-hospital transit of patients and had the fatality rate as high as 40%. High attack rate and easy spread of MERS indicate an epidemic risk. [5] [6] [7] To date, no specific vaccine is available for MERS. The steep increase in MERS cases and its high mortality rate demand an urgent need for specific and safe MERS vaccine. Thus far, only little information is known about the pathogenesis of MERS-CoV. Hence, an immunoinformatics approach to thoroughly study and screen immunogenic proteins from the available proteome sequence data of MERS-CoV is essential for vaccine design. The proteome of MERS-CoV consists of several important vaccine candidate and drug target proteins. These proteins are involved in infection and pathogenesis of MERS-CoV to human host cells. Spike (S) glycoprotein, in particular, its receptor-binding domain, is involved in virus and host cell interaction. 8 Envelope (E) protein plays an important role in host cell recognition. 9 Nucleocapsid (N) protein is involved in RNA binding during ribonucleocapsid formation by MERS-CoV. 10 Membrane (M) protein has interferon (IFN)-antagonizing properties, thus reducing IFN levels in infected patients. 11 Open reading frame (ORF) proteins also play critical roles in viral infection and pathogenesis. A mutation study of ORFs (ORF3, ORF4a, ORF4b, and ORF5) indicated that ORFs have major implications in viral infection involving disruption of host cell processes, dysregulated IFN pathway activation, and abrupt inflammation. 12 Proteins ORF1a (4P16) and ORF1ab (4WUR) are papain-like proteases (PL(pro)) involved in viral infection and are potential targets for the development of antiviral drugs. 13 PL(pro) facilitate infection by their proteolytic, deubiquitinating, and deISGylating activities suppressing innate immune response from the host cell. Protein ORF1a (4RSP), a protease (3CLpro), facilitates proteolytic activity during viral infection and replication. 14, 15 Protein ORF1ab (5WWP) is a helicase protein of MERS-CoV, and it is one of the most conserved proteins among nidoviruses. Protein ORF8b is a highly conserved protein, and during infection, it induces various immune responses. 16 In the present study, we propose two multi-epitope vaccine (MEV) designs for MERS. The proposed MEVs are composed of cytotoxic T lymphocyte (CTL) and helper T lymphocyte (HTL) epitopes. Both CTL and HTL MEVs contain overlapping regions of linear B-cell epitopes. Both MEVs also contain human β-defensin-2 (hBD-2) and human β-defensin-3 (hBD-3) as adjuvants to enhance the immunogenic response. 17, 18 Methodology To design the MEVs, different in silico tools were used to screen potential CTL, HTL, and B-cell epitopes from 13 MERS-CoV proteins that are involved in virus-host interaction and viral proliferation. All three types of epitopes were studied for overlapping consensus regions. CTL and HTL epitopes showing partial or complete overlapping regions with all or any of the three kinds of epitopes and the epitopes with the highest number of human leukocyte antigen (HLA) allele binders were selected for detailed analysis. The selected CTL and HTL epitopes were further validated for their molecular interaction with their HLA allele binders. Furthermore, the selected CTL epitopes were validated for their molecular interaction with cavity of transporter associated with antigen processing (TAP) to analyze their smooth passage from the cytoplasm to the endoplasmic reticulum (ER) lumen. 19 Three-dimensional (3D) structure model for both MEVs was generated, refined, validated, and analyzed for different physiochemical properties. Both CTL and HTL MEV models were further screened to have several B-cell discontinuous epitopes and IFN-γ-inducing epitopes. Because the immune system broadly recognizes pathogens partially due to the involvement of toll-like receptor 3 (TLR-3) and its signaling cascade, both CTL and HTL MEVs were further studied for their molecular interaction with TLR-3. 20,21 Moreover, the cDNAs of both MEVs were analyzed and were predicted to be highly expressing in a mammalian host cell line (human) ( Figure S1 ). For vaccine development, proteins critically involved in virus-host cell interaction, viral proliferation, and host cell disruption were chosen. Thirteen proteins of MERS-CoV proteome chosen for epitope screening included the following: S glycoprotein; E protein; N protein; M protein; ORFs (ORF3, ORF4a, ORF4b, and ORF5); ORF1a (4P16), a papain-like protease (PL(pro)); ORF1a (4RSP), a protease structural basis of development of multi-epitope vaccine against MeRs (3CLpro); ORF1ab (4WUR), a PL(pro); ORF1ab (5WWP); and ORF8b. For sequence-based epitope screening, fulllength protein sequences of the abovementioned MERS-CoV proteins were retrieved from the NCBI database (National Center for Biotechnology Information, https://www.ncbi. nlm.nih.gov/protein). A total of 2,499 amino acid sequences belonging to different proteins of MERS-CoV of different strains and origins were retrieved. For structure-based epitope screening, available 3D structures of MERS-CoV proteins were retrieved from Protein Data Bank (PDB, http://www.rcsb.org/pdb/home/home.do). For the proteins with no structure available in PDB, homology modeling was performed using the Swiss-model, (http://swissmodel. expasy.org/) (Table S1 ). 22 To screen CTL epitopes, the Immune Epitope Database (IEDB) tool "Proteasomal cleavage/TAP transport/MHC class I combined predictor" (http://tools.iedb.org/processing/) was used. [23] [24] [25] The total score generated by the tool is a combined score of proteasome, TAP (N-terminal interaction), major histocompatibility complex (MHC), and processing analysis scores. The score is generated using the combination of six methods, namely, Consensus, NN-align, SMM-align, Combinatorial Library, Sturniolo, and NetMHCIIpan. The IC50 (nM) value was also obtained using this IEDB tool. Epitopes with high, intermediate, and least affinity of binding to an HLA allele have IC50 values <50 nM, <500 nM, and <5,000 nM, respectively. Immunogenicity of the screened CTL epitopes was also obtained by the "MHC I Immunogenicity" tool of IEDB (http://tools.iedb.org/ immunogenicity/) with all parameters set to default (first, second, and C-terminus amino acids). 26 The tool predicts immunogenicity of a peptide MHC complex on the basis of physiochemical properties of amino acids and their position within the short peptide sequence. To screen HTL epitopes, the IEDB tool "MHC-II Binding Predictions" (http://tools.iedb.org/mhcii/) was used. The percentile rank for each peptide is generated by the combination of three methods (Combinatorial Library, SMM-align, and Sturniolo) by comparing the score of peptide against the scores of other random five million 15-mer peptides from the SwissProt database. [27] [28] [29] [30] The rank for the Consensus method was generated by the median percentile rank of the three methods. The "Population Coverage" tool of IEDB (http://tools. iedb.org/population/) was used to analyze the world human population coverage by the shortlisted 28 CTL and 28 HTL epitopes. 31 Use of multi-epitopes involving both CTL and HTL epitopes would have the higher probability of larger human population coverage worldwide. B-cell epitope prediction sequence-based B-cell epitope prediction Protein sequence-based linear B-cell epitopes were predicted by six different prediction methods available at "B Cell Epitope Prediction Tools" tool of IEDB server (http://tools. iedb.org/bcell/); these tools are based on the propensity scale method and physicochemical properties of the antigenic sequence. These methods include "BepiPred Linear Epitope Prediction" (propensity scale method such as hidden Markov model), "Chou & Fasman Beta-Turn Prediction", "Emini Surface Accessibility Prediction", "Karplus & Schulz Flexibility Prediction", "Kolaskar & Tongaonkar Antigenicity", and "Parker Hydrophilicity Prediction". [32] [33] [34] [35] [36] [37] structure-based B-cell epitope prediction Two structure-based B-cell epitope prediction methods, namely, DiscoTope 2.0 (Structure-based Antibody Prediction tool; http://tools.iedb.org/discotope/) and Ellipro (Antibody Epitope Prediction tool; http://tools.iedb.org/ellipro/), available in IEDB, were used for linear and discontinuous B-cell epitope prediction. 38, 39 The ElliPro method is based on the location of the residue in the protein's 3D structure. The residues lying outside of the ellipsoid covering 90% of the inner core residues of the protein score highest protrusion index (PI) of 0.9 and so on. Discontinuous epitopes predicted by ElliPro are clustered on the basis of distance R in Å between two residues' centers of mass lying outside the largest possible ellipsoid. The larger the value of R, the larger will be the number of discontinuous epitopes clustered. DiscoTope 2.0 is based on the number of contacts of the residues' Cα carbon atom with other Cα carbon atoms in the 3D structure within the 10 Å distance and the propensity of a residue to be a part of an epitope. Toxicity assessment of CTL and HTL epitopes and B-cell epitopes was performed by ToxinPred (http://crdd.osdd.net/ raghava/toxinpred/multi_submit.php) analysis. This analysis allows to identify highly toxic or nontoxic peptides. The analysis was performed by the "support vector machine (SVM) (SwissProt)-based" method with all the parameters set to default. 41 Overlapping residue analysis Multiple Sequence Alignment (MSA) analysis using the Clustal Omega tool (https://www.ebi.ac.uk/Tools/ msa/clustalo/) of European Bioinformatics Institute was performed for all the shortlisted CTL, HTL, and B-cell epitopes from 13 MERS-CoV proteins. 42 MSA by Clustal Omega virtually aligns any number of protein sequences and delivers accurate alignments. epitope selected for molecular interaction study with hla allele and TaP CTL and HTL epitopes were shortlisted for further in silico analysis on the basis of their overlapping sequence regions among all three types of epitopes (CTL, HTL, and B-cell) or complete overlap among any two types of epitopes or the highest number of HLA allele binders. Tertiary structure modeling of hla alleles and selected T-cell epitopes Template-based homology modeling for HLA class I and II allele binders of the shortlisted epitopes was performed by the Swiss-model. 22 Protein sequences of HLA class I and II alleles were retrieved from Immuno Polymorphism Database (IPD-IMGT/HLA) (https://www.ebi.ac.uk/ipd/ imgt/hla/allele.html). 69 Template with high sequence identity was chosen for modeling, and the models thus generated were validated for quality by the Qualitative Model Energy ANalysis (QMEAN) analysis. The QMEAN value is a composite (both global and local [ie, per residue] structure) quality estimate of the generated model. 43 Models with an acceptable QMEAN value with a cutoff of −4.0 were chosen for further studies (Table S2) . "Natural Peptides Module for Beginners" module of the PEPstrMOD tool (http://osddlinux.osdd.net/raghava/ pepstrmod/nat_ss.php) was used to generate tertiary structures of the selected epitopes. 44 Prediction was carried out with a simulation time window of 100 ps, and the peptide environment was set to vacuum. AutoDock Tool 4.2 and AutoDock Vina were used for molecular docking study of the selected epitopes and their respective HLA allele binders. 45, 46 Further, the docked complexes were subjected to MD simulation study by Gromacs 5.1.4 using Optimized Potentials for Liquid Simulations -all atom (OPLS-AA) force field. 47, 48 Molecular interaction analysis (docking) of the selected cTl epitopes with TaP Molecular docking study of the shortlisted CTL epitopes with TAP receptor was performed by AutoDock Vina. 45, 46 For more accurate prediction, instead of homology modeling, the cryo-EM structure of TAP (PDB ID: 5u1d) was used. 49 For docking, the antigen from the TAP cavity of the original structure was removed. From the screened and shortlisted CTL and HTL epitopes, two MEVs were designed using EAAAK and GGGGS as short peptide rigid and flexible linkers, respectively ( Figure 1A and B). To enhance the immune response, hBD-2 (PDB ID: 1FD3, sequence: GIGDPVTCLKSGAICHPVFCPRRYKQIGTCG LPGTKCCKKP) and hBD-3 (PDB ID: 1KJ6, sequence: GII NTLQKYYCRVRGGRCAVLSCLPKEEQIGKCSTRGRK CCRRKK) were used as adjuvants for both MEVs at N and C terminals, respectively. 17, 18, [50] [51] [52] [53] Upon lung infection, the expression of hBD-2 and hBD-3 was found to be increased. β-Defensins are involved in chemotactic activity for memory T cells, monocytes, and immature dendritic cells as well as in degranulation of mast cells. Thus, hBDs enhance innate and adaptive immunity and therefore were chosen here as adjuvants for the design of MEVs. IFn-γ-inducing epitope prediction IFN-γ epitopes with potential to induce the release of IFN-γ from CD4+ T cells from both MEVs were predicted by "IFNepitope" server (http://crdd.osdd.net/raghava/ ifnepitope/scan.php) using the "Motif and SVM hybrid" (MERCI: Motif-EmeRging and with Classes-Identification, and SVM) approach. The tool generates overlapping sequences from the query protein/antigen and uses it for IFN-γ-inducing epitope prediction. The prediction is based on a dataset of IFN-γ-inducing and IFN-γ-noninducing MHC class II binders. 54 For both MEVs, AlgPred was used for allergenicity prediction (http://crdd.osdd.net/raghava/algpred/submission.html). 56 The AlgPred prediction of allergens is based on similarity of the known epitope with any region of the submitted protein. VaxiJen (http://www.ddg-pharmfac.net/vaxijen/ VaxiJen/VaxiJen.html) was used to analyze the probability of antigenicity of both MEVs. 57 The VaxiJen analysis applies an alignment-free approach that is solely based on the physicochemical properties of the sequences of the submitted proteins. To analyze the physiochemical properties of CTL and HTL MEVs, ProtParam (https://web.expasy.org/protparam/) was used. 58 ProtParam analysis performs empirical investigation for a given protein amino acid sequence. Both the CTL and HTL MEVs were further subjected to 3D structure modeling by RaptorX structure prediction server (http://raptorx.uchicago.edu/StructurePrediction/predict/). 59 The quality of the generated model is indicated by its P-value. The P-value is the probability of a predicted model being worse than the best. It indicates relative quality in terms of modeling error by combining global distance test (GDT) and un-normalized GDT, that is, error at each residue. The smaller the P-value, the higher is the quality of the predicted model. The generated MEV models were further refined by ModRefiner and GalaxyRefin. 60, 70 GalaxyRefine performs repeated structure perturbation along with subsequent structural relaxation by dynamics simulation to refine a protein structure. To avoid breaks in model structures, GalaxyRefine uses the triaxial loop closure method. The refined 3D models of both MEVs were then subjected to RAMPAGE analysis for generating a Ramachandran plot (http://mordred.bioc.cam.ac.uk/~rapper/rampage.php). 61 The Ramachandran plot shows the residues that form energetically allowed and disallowed dihedral angles psi (ψ) and phi (ϕ), which are calculated on the basis of van der Waal radius of the amino acid side chain. Both the designed CTL and HTL MEVs were subjected to discontinuous B-cell epitope prediction to analyze structurebased humoral immunogenic potential of both MEVs using the ElliPro tool. 39 Molecular docking and MD simulation study of MeVs and the immunological receptor TlR-3 To study the molecular interaction, the refined models of CTL and HTL MEVs were docked with TLR-3 by PatchDock server (http://bioinfo3d.cs.tau.ac.il/PatchDock/). [62] [63] [64] The 3D structure of human TLR-3 ectodomain (ECD) was retrieved Optimized cDNAs for both MEVs were generated by the Codon Usage Wrangler Tool with the option of a mammalian host cell line (human) as the expression system (http://www. mrc-lmb.cam.ac.uk/ms/methods/codon.html). Further, the GenScript Rare Codon Analysis Tool (https://www.genscript. com/tools/rare-codon-analysis) was used to analyze the cDNAs of both MEVs. The tool provides GC content, codon adaptation index (CAI), and tandem rare codon frequency for cDNA. [65] [66] [67] Results and discussion screening of potential epitopes T-cell epitope prediction screening of cTl epitope On the basis of "total score" and acceptably low IC50 (nM) value of epitope-HLA class I allele pairs, 75 CD8+ CTL epitopes were chosen. Later, 28 epitopes were shortlisted with high scoring and larger number of HLA class I allele binders for further studies. Immunogenicity of the screened CTL epitopes was also determined. The higher the immunogenicity score, the greater is the immunogenicity of the epitope (Tables S3 and S4) . screening of hTl epitopes Screening of HTL epitopes was performed on the basis of "percentile rank". Small percentile rank shows the high affinity of the peptide for its respective HLA allele. Initially, from the 13 MERS-CoV proteins, 70 CD4+ T-cell epitopes with the highest percentile rank were screened, and 28 epitopes with high percentile rank and the highest number of HLA class II allele binders were then shortlisted (Tables S5 and S6 ). Several CTL and HTL epitopes predicted in the present study show overlapping regions with the epitopes detected to induce T-cell responses in previous studies done using peripheral blood mononuclear cells from infected patients. 68 Hence, epitopes screened and reported in the present study could be predicted to originate from the highly immunogenic stretch of MERS-CoV protein sequences. In this study, most geographical regions of the world, and in particular, the countries of Middle East, South Asia, East Asia, and Northeast Asia, were included. This study indicates that the combined use of all shortlisted CTL and HTL epitopes would have an average worldwide population coverage of 94.0%, with an SD of 20.19 (Table S7 ). sequence-based B-cell epitope prediction Initially, a total of 144 B-cell epitopes were screened from 13 MERS-CoV proteins, with epitope length of four or more than four amino acids by the BepiPred Linear Epitope Prediction method. The epitopes predicted by five other methods based on different physicochemical properties showed a significant consensus overlap of amino acid sequences with that of BepiPred Linear Epitope Prediction. From the 144 B-cell epitopes, 12 with the length of 4-19 amino acids were shortlisted (Table S8 and Figure 2 ). Discontinuous and linear epitopes predicted by DiscoTope 2.0 and Ellipro methods showed significant consensus overlap of amino acid sequence with the linear epitopes predicted by the BepiPred linear epitope method (Table S8 and Figure 2 ). This result confirms that the shortlisted BepiPred linear epitopes are highly immunogenic B-cell epitopes. The conservation analysis of shortlisted 28 CTL, 28 HTL, and 12 B-cell epitopes shows that amino acid sequence conservancy of CTL epitopes varied from 72.7% to 100%, that of HTL epitopes varied from (50% for two epitopes) 68.18% to 100%, and that of B-cell epitopes varied from 85.71% to 99.26% (Tables S3, S5 , and S8). This result indicates high conservancy nature of the shortlisted CTL, HTL, and B-cell epitopes. All the shortlisted CTL, HTL, and B-cell epitopes analyzed by the ToxinPred tool were predicted to be nontoxic (Tables S3, S5 , and S8). The analysis was based on the ToxinPred main dataset consisting of 1,805 toxic peptides for its prediction. MSA analysis of all the screened CTL, HTL, and B-cell epitopes from the 13 MERS-CoV protein candidates revealed that several epitopes of CTL, HTL, and B-cell have overlapping amino acid sequence regions. Epitopes with two or more than two residues in the overlapping region are shown in Figure 3 . Overlapping cTl, hTl, and B-cell epitopes. Notes: Overlapping cTl (red), hTl (blue), and B-cell epitopes (green) were sorted out by Msa using clustal Omega at eBI server. The ringed clusters of epitopes involve all three types of epitopes, epitopes with full sequence overlap, and the epitopes with highest number of hla allele binders. Abbreviations: cTl, cytotoxic T lymphocyte; eBI, european Bioinformatics Institute; hla, human leukocyte antigen; hTl, helper T lymphocyte; Msa, multiple sequence alignment. Figure 2 Overlapping regions of the BepiPred linear B-cell epitopes and the epitopes predicted by other methods. Notes: This analysis shows a strong consensus between the BepiPred linear B-cell epitopes and the epitopes predicted by other sequence-based and structure-based epitope prediction methods. Different colors are used to highlight overlapping regions of epitopes predicted by the different methods. epitope selected for molecular interaction study with hla allele and TaP Overlapping CTL and HTL epitopes clustering with all three (CTL, HTL, and B-cell) types of epitopes or with complete epitope sequence overlap or with the highest number of HLA allele binders were identified and selected for their molecular interaction analysis with HLA alleles and TAP ( Figure 3 ). These include CTL epitopes: 47 Molecular interaction analysis of the selected epitopes with hla allele Molecular docking and MD simulation study of the selected epitopes with hla alleles Molecular docking study of all the selected epitopes with their respective HLA allele binders revealed a significantly favorable molecular interaction. Docking complexes thus formed have significantly negative binding energy, and several amino acid residues of epitopes and HLA alleles were involved in hydrogen bond formation ( Figure 4 ). To analyze the stability of binding, docking complexes were further subjected to MD simulation study with an analysis time window of 0.1 ns at the reasonably invariable temperature (~300 K) and pressure (~1 bar). The results of MD simulation for all the epitope-HLA allele complexes showed a very convincing reasonably invariant root mean square deviation (RMSD) value between ~0.1 and 0.2 nm, thus indicating stable complex formation ( Figure 5) . Moreover, the reasonably invariant radius of gyration (Rg) of the complexes ( Figure S2 ) and RMS fluctuation (RMSF) for all the atoms in the complexes ( Figure S3 ) indicate that the epitope-HLA complexes remain very stable in their folded form. B-factor of the epitope and HLA allele complexes is shown in the rainbow color presentation in Figure S4 . Most of the regions of the complexes are stable (blue) with a very small region being acceptably fluctuating (yellow and orange). Molecular docking results show a favorable molecular interaction between the selected CTL epitopes and the TAP cavity. The result reveals several molecular interaction sites for the selected epitopes within the TAP cavity. Two sites of interaction-one close to the cytoplasm and another close to the ER lumen-are shown in Figure 6 . All the interactions showed a significantly negative binding energy with one or more than one hydrogen bond formation at both sites of interaction. From this study, we may predict a smooth passage for selected CTL epitopes from the cytoplasm to the ER lumen throughout the entire passage of the TAP cavity ( Figure 6 ). IFn-γ-inducing epitope prediction IFN-γ is involved in both adaptive and innate immune responses. It stimulates macrophages and natural killer cells. IFN-γ-inducing 15-mer peptide epitopes were predicted from (Table S9 and Figure 7C and G). The AlgPred analysis of CTL and HTL MEVs revealed that both MEVs are nonallergens with the score of −0.95415008 and −0.8647714, respectively, with sensitivity and specificity of 88.87% and 81.86%, respectively, while the default threshold value was −0.4. The VaxiJen analysis indicated both MEVs to be probable ANTIGENS with prediction score of 0.5302 and 0.5097 for CTL MEVs and HTL MEVs, respectively, while the default threshold value was 0.4 for viral proteins. Hence, both MEVs are predicted to be nonallergic and potentially antigenic in nature. ProtParam analysis showed that CTL MEVs have 500 amino acids, 50.4 kDa molecular weight, and 9.72 theoretical pI. The half-life in mammalian reticulocytes, yeast, and Escherichia coli was 30 hours, 20 hours, and 10 hours, respectively; aliphatic index was 61.68, and grand average of hydropathicity (GRAVY) was −0.020, indicating globular and hydrophilic nature of CTL MEVs; the instability index was 50.70, indicating that CTL MEVs are theoretically close to stable in nature. ProtParam study showed HTL MEV has 657 amino acids, 67.583 kDa molecular weight, and 9.29 theoretical pI. The half-life in mammalian reticulocytes, yeast, and E. coli was 30 hours, 20 hours, and 10 hours, respectively; aliphatic index was 94.22, and grand average of hydropathicity (GRAVY) was 0.337, both indicating globular and hydrophilic nature of HTL MEV; the instability index was 37.66, indicating theoretically stable nature of HTL MEV. The model generated by RaptorX for CTL MEV has 4% helix, 37% β-sheet, and 58% coil; the model is 26% exposed, 41% medium, and 31% buried. It has three domains: first domain (1-46 aa, template-1fd3:A); second domain (47-444 aa, template-4dou:A, 4zch:A); and third domain (445-500 aa, template-1kj6:A) ( Figure 7A ). The model generated by RaptorX for HTL MEV has 35% helix, 20% β-sheet, and 44% coil; the model is 26% exposed, 35% medium, and 37% buried. It has four domains: first domain (1-46 aa, template-1fd3:A); second domain (47-546 aa, templates-4ufcA, 2eabA, and 2rdyA); third domain (547-601 aa, template-1tue:B); and fourth domain (602-657 aa, template-1kj6:A) ( Figure 7E ). The P-values of the best template of CTL and HTL MEV models are 3.98e-05 and 5.41e-04, respectively. Good quality alpha proteins and good quality beta proteins should have P-value of less than 10 -3 and 10 -4 respectively. Hence, both the models of CTL and HTL MEVs generated by RaptorX are of good quality. Because the CTL and HTL epitopes used for the design of both MEVs show overlapping with B-cell linear epitopes (Figure 3 ), both the CTL and HTL MEV models also carry overlapping regions of B-cell linear epitopes ( Figure 7B and F) . The structural accuracy of initial model including joining of the gaps was performed by ModRefiner refinment. Galaxy-Refine was used to further refine the 3D models of CTL and HTL MEVs. For both MEVs, refined model 1 was chosen on the basis of best scorings. For CTL MEV refinement, the score of model 1 was as follows: Rama favored was 94.2%, GDT-HA was 0.9500, RMSD was 0.433, MolProbity was 2.211, Clash score was 18.3, and Poor rotamers was 1.2. For HTL MEV refinement, the score of model 1 was as follows: Rama favored was 93.4%, GDT-HA was 0.9502, RMSD was 0.412, MolProbity was 2.262, Clash score was 21.9, and Poor rotamers was 0.6. Here, MolProbity is the log-weighted combination of clash score, percentage Ramachandran not favored, and percentage bad side-chain rotamers. After refinement, all the mentioned parameters improved significantly as compared to those of the initial model (Table S10 ). The RAMPAGE analysis showed that the refined CTL MEV model has 94% residues in the favored region, 3.8% residues in the allowed region, and only 1.4% residues in the outlier region, while the refined HTL MEV model has 94% of residues in the favored region, 5.2% residues in the allowed region, and only 0.8% residues in the outlier region ( Figure 7D and H). Discontinuous B-cell epitope prediction from both MEVs was performed by the ElliPro tool. The screening revealed CTL MEV to have six discontinuous epitopes and HTL MEV to have five discontinuous epitopes. The PI score for CTL MEV epitopes ranged from 0.549 to 0.906 and that for HTL MEV ranged from 0.56 to 0.729 (Table S11 and Figure 7C and G). The higher the score, the greater is the potential of the B-cell discontinuous epitope. The refined models of CTL and HTL MEVs were further analyzed for their interaction with the ECD of human TLR-3 by molecular docking using the PatchDock tool. Docking conformation chosen for the CTL and HTL MEVs showed scores of 18,096 and 23,690, respectively, which were highest among all docked complexes. The highest score indicates the best geometric shape complementarity fitting of ligand and receptor predicted by the tool. The docking complex shows a fitting confirmation of both the MEVs within the ECD of TLR-3 ( Figure 8A and E). Further, the MD simulation analysis of the docked CTL-TLR-3 and HTL-TLR-3 complexes showed a very convincing reasonably stable RMSD value between ~0.4 and 0.5 nm for a given time window of 6 ns at the reasonably invariable temperature (~300 K) and pressure (~1 bar). These results indicate a stable complex formation for both MEVs with TLR-3 ( Figure 8B and F). The reasonably invariant Rg of MEV-TLR-3 complexes ( Figure 8C and G) and RMSF for all the atoms in the complexes ( Figure 8D and H) indicate that both the MEV-TLR-3 complexes are reasonably stable and are properly folded. The B-factor of CTL and HTL MEV complexes with the TLR-3 receptor is shown by rainbow color presentation in Figure 8A and Optimized cDNAs for both CTL and HTL MEVs were generated by the Codon Usage Wrangler Tool with a mammalian host cell line (human) as the choice of the expression system. The GenScript Rare Codon Analysis Tool showed that the GC content of optimized CTL MEV cDNA was 69.74% and CAI score was 1.0 with 0% tandem rare codons. Likewise, the GC content of optimized HTL MEV cDNA was 66.63%, and the CAI score was 1.0 with 0% tandem rare codons. Ideally, the GC content should be 30%-70%, the CAI score indicating the possibility of cDNA expression in the chosen expression system should be between 0.8 and 1.0, and tandem rare codon frequency indicating the percentage of low-frequency codons present in cDNA for the chosen expression system should be <30%. Tandem rare codons may hinder the expression of cDNA or even interrupt the translational machinery. Hence, the optimized cDNAs of both MEVs are predicted to be highly expressing in the mammalian host cell line (human). In this study, we propose the design of two MEVs against MERS-CoV that consisted of CTL and HTL epitopes. The selected CTL and HTL epitopes were validated by in silico methods for their molecular interaction with HLA alleles and TAP (for CTL epitopes). The population coverage by the designed MEVs is as high as 94% of the world population. Both MEVs were found to contain IFN-γ epitopes and B-cell discontinuous epitopes. Moreover, they also showed stable interaction with the immunoreceptor TLR-3. On the basis of the design and in silico validation of both the CTL and HTL MEVs, their joint administration is predicted to induce humoral and cell-mediated immune response. Codon-based cDNAs of both MEVs are predicted to have high expression in the mammalian host cell line (human); hence, they could be cloned and expressed at the laboratory level for in vivo trials on humanized HLA-expressing mice for further study. Available from: www.who.int/emergencies/mers-cov/mers-summary-2016 KSA MERS-CoV Investigation Team. Hospital outbreak of Middle East respiratory syndrome coronavirus MERS-CoV outbreak in Jeddah -a link to health care facilities Hospital Outbreaks of Middle East Respiratory Syndrome Epidemiologic features of the first MERS outbreak in Korea: focus on Pyeongtaek St. Mary's Hospital MERS outbreak in Korea: hospital-to-hospital transmission Reponse: Coronavirus Infections Evaluation of candidate vaccine approaches for MERS-CoV Analysis of the genome sequence and prediction of B-cell epitopes of the envelope protein of Middle East respiratory syndrome-coronavirus Crystallographic analysis of the N-terminal domain of Middle East respiratory syndrome coronavirus nucleocapsid protein Middle East respiratory syndrome coronavirus M protein suppresses type I interferon expression through the inhibition of TBK1-dependent phosphorylation of IRF3 MERS-CoV Accessory ORFs Play Key Role for Infection and Pathogenesis Crystal structure of the papain-like protease of MERS coronavirus reveals unusual, potentially druggable active-site features Ligand-induced Dimerization of Middle East Respiratory Syndrome (MERS) Coronavirus nsp5 Protease (3CLpro): IMPLICATIONS FOR nsp5 REGULATION AND THE DEVELOPMENT OF ANTIVIRALS Structural and mutational analysis of the interaction between the Middle-East respiratory syndrome coronavirus (MERS-CoV) papain-like protease and human ubiquitin Sequencing and Phylogenetic Analyses of Structural and Accessory Proteins of Middle East Respiratory Syndrome Coronavirus from the First Imported Case in China Antiviral mechanisms of human defensins Rhinovirus increases human β-defensin-2 and-3 mRNA expression in cultured bronchial epithelial cells. Pathogens and Disease Assembly and export of MHC class I peptide ligands Innate immunity: structure and function of TLRs Toll-like receptor 3 signaling via TRIF contributes to a protective innate immune response to severe acute respiratory syndrome coronavirus infection The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling Modeling the MHC class I pathway by combining predictions of proteasomal cleavage, TAP transport and MHC class I binding Identifying MHC class I epitopes by predicting the TAP transport efficiency of epitope precursors NetMHCpan, a method for MHC class I binding prediction beyond humans Properties of MHC class I presented peptides that enhance immunogenicity Peptide binding predictions for HLA DR, DP and DQ molecules Quantitative peptide binding motifs for 19 human and mouse MHC class I molecules derived using positional scanning combinatorial peptide libraries Prediction of MHC class II binding affinity using SMM-align, a novel stabilization matrix alignment method Generation of tissue-specific and promiscuous HLA ligand databases using DNA microarrays and virtual HLA class II matrices Predicting population coverage of T-cell epitope-based diagnostics and vaccines Improved method for predicting linear B-cell epitopes Prediction of the secondary structure of proteins from their amino acid sequence Induction of hepatitis A virus-neutralizing antibody by a virus-specific synthetic peptide Prediction of chain flexibility in proteins A semi-empirical method for prediction of antigenic determinants on protein antigens New hydrophilicity scale derived from high-performance liquid chromatography peptide retention data: correlation of predicted surface residues with antigenicity and X-rayderived accessible sites Reliable B cell epitope predictions: impacts of method development and improved benchmarking ElliPro: a new structure-based tool for the prediction of antibody epitopes Development of an epitope conservancy analysis tool to facilitate the design of epitope-based diagnostics and vaccines In silico approach for predicting toxicity of peptides and proteins Fast, scalable generation of highquality protein multiple sequence alignments using Clustal Omega QMEAN: A comprehensive scoring function for model quality assessment PEPstrMOD: structure prediction of peptides containing natural, non-natural and modified residues AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading High performance molecular simulations through multi-level parallelism from laptops to supercomputers Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids Structure of the transporter associated with antigen processing trapped by herpes simplex virus A flexible peptide linker enhances the immunoreactivity of two copies HBsAg preS1 (21-47) fusion protein The structure of human beta-defensin-2 shows evidence of higher order oligomerization Immunoinformatics analysis and in silico designing of a novel multi-epitope peptide vaccine against Staphylococcus aureus Fusion protein linkers: property, design and functionality VaccineDA: Prediction, design and genome-wide screening of oligodeoxynucleotide-based vaccine adjuvants Designing of interferon-gamma inducing MHC class-II binders AlgPred: prediction of allergenic proteins and mapping of IgE epitopes VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines Protein identification and analysis tools on the ExPASy server Template-based protein structure modeling using the RaptorX web server GalaxyWEB server for protein structure prediction and refinement Structure validation by Calpha geometry: phi,psi and Cbeta deviation The molecular structure of the Tolllike receptor 3 ligand-binding domain structural basis of development of multi-epitope vaccine against MeRs Efficient unbound docking of rigid molecules PatchDock and SymmDock: servers for rigid and symmetric docking Designing an efficient multi-epitope oral vaccine against Helicobacter pylori using immunoinformatics and structural vaccinology approaches Synonymous codon usage pattern in glycoprotein gene of rabies virus Computational identification of rare codons of Escherichia coli based on codon pairs preference Recovery from the Middle East respiratory syndrome is associated with antibody and T-cell responses The IPD and IMGT/HLA database: allele variant databases Improving the Physical Realism and Structural Accuracy of Protein Models by a Two-step Atomic-level Energy Minimization The authors acknowledge the Advance Instrumentation Research Facility (AIRF) at Jawaharlal Nehru University, New Delhi, for providing advance computational facility to conduct experiments. The current research was also supported by DST-PURSE grant to Jawaharlal Nehru University. S Srivastava and MK designed the protocol and performed the experiments. S Srivastava, MK, S Singh, AKS, and NS contributed toward data analysis, drafting and revising the article, gave final approval of the version to be published, and agree to be accountable for all aspects of the work. The authors report no conflicts of interest in this work. Infection and Drug Resistance is an international, peer-reviewed openaccess journal that focuses on the optimal treatment of infection (bacterial, fungal and viral) and the development and institution of preventive strategies to minimize the development and spread of resistance. The journal is specifically concerned with the epidemiology of antibiotic resistance and the mechanisms of resistance development and diffusion in both hospitals and the community. The manuscript management system is completely online and includes a very quick and fair peerreview system, which is all easy to use. Visit http://www.dovepress.com/ testimonials.php to read real quotes from published authors.