key: cord-0928211-m4ev5opy authors: Haq, Azaz Ul; Khan, Abbas; Khan, Jafar; Irum, Shamaila; Waheed, Yasir; Ahmad, Sajjad; Nizam-Uddin, N.; Albutti, Aqel; Zaman, Nasib; Hussain, Zahid; Ali, Syed Shujait; Waseem, Muhammad; Kanwal, Fariha; Wei, Dong-Qing; Wang, Qian title: Annotation of Potential Vaccine Targets and Design of a Multi-Epitope Subunit Vaccine against Yersinia pestis through Reverse Vaccinology and Validation through an Agent-Based Modeling Approach date: 2021-11-15 journal: Vaccines (Basel) DOI: 10.3390/vaccines9111327 sha: 9740657b8aa0b3a01d048f7be64dd94b0b47b21a doc_id: 928211 cord_uid: m4ev5opy Yersinia pestis is responsible for plague and major pandemics in Asia and Europe. This bacterium has shown resistance to an array of drugs commonly used for the treatment of plague. Therefore, effective therapeutics measurements, such as designing a vaccine that can effectively and safely prevent Y. pestis infection, are of high interest. To fast-track vaccine development against Yersinia pestis, herein, proteome-wide vaccine target annotation was performed, and structural vaccinology-assisted epitopes were predicted. Among the total 3909 proteins, only 5 (rstB, YPO2385, hmuR, flaA1a, and psaB) were shortlisted as essential vaccine targets. These targets were then subjected to multi-epitope vaccine design using different linkers. EAAK, AAY, and GPGPG as linkers were used to link CTL, HTL, and B-cell epitopes, and an adjuvant (beta defensin) was also added at the N-terminal of the MEVC. Physiochemical characterization, such as determination of the instability index, theoretical pI, half-life, aliphatic index, stability profiling, antigenicity, allergenicity, and hydropathy of the ensemble, showed that the vaccine is highly stable, antigenic, and non-allergenic and produces multiple interactions with immune receptors upon docking. In addition, molecular dynamics simulation confirmed the stable binding and good dynamic properties of the vaccine–TLR complex. Furthermore, in silico and immune simulation of the developed MEVC for Y. pestis showed that the vaccine triggered strong immune response after several doses at different intervals. Neutralization of the antigen was observed at the third day of injection. Conclusively, the vaccine designed here for Y. pestis produces an immune response; however, further immunological testing is needed to unveil its real efficacy. Yersinia pestis is an etiologic agent of plague, which is infamous for three major pandemics killing millions of people [1] . Of these, the major devastating pandemic occurred in Europe, affecting a significant portion of the population [2] . The pathogen also has excellent potential for use in bioterrorism [3] . Y. pestis is a Gram-negative bacterium that transfers from fleas to humans and rodents [1] . Experimental studies performed on mice showed that after incubation of Y. pestis in mice, high growth of the pathogen was evident, leading to death of the mice in a few days [4] . It was noticed that Y. pestis caused inflammation and affected the function of the immune system [5] . This pathogenicity is the outcome of contributions from lipopolysaccharides (LPSs) that hold vital importance in tackling the host immune system [6] . The pathogen has two different plasmids that encode a number of virulent proteins. One is the pPCP1 plasmid, which encodes an activator of virulence essential plasminogen (Pla). The second pMT1 plasmid codes for F1 capsular protein and murine toxin Ymt, which are significant in plague infection transmission. Ymt encodes murine toxin and is lethal but required for its survival in fleas [6] . Furthermore, Y. pestis has machinery that inhibits MAPK pathway activation and the release of the proinflammatory cytokine TNF-α, according to recent findings. The same study also found that two Y. pestis strains, the fully virulent Kimberley53 and the pgm attenuated EV76 strain, are less capable of inducing macrophage apoptosis, impairing NF-B activation, and activating caspase pathways than the virulent Y. enterocolitica O:8 strain under varied infection circumstances. These limitations correlate with inefficient translocation of YopJ by Y. pestis, which leads to slow accumulation of the effector in the target cells. In addition, they also found that Y. enterocolitica-derived YopP can provide effective apoptotic potential to Y. pestis, and they suggested that this can be attributed to the fact that Y. enterocolitica O:8 YopP is better adapted to translocation by the Yersinia TTSS than YopJ [7] . Y. pestis interacts with Toll-like receptor (TLR) immune cells. TLRs are receptors present on human innate immune cells. TLRs can detect a large number of pathogenic molecules called PAMPs (pathogen-associated molecular patterns) on the bacterial cell wall components, LPS (lipopolysaccharide), and nucleic acid [8] . TLRs activate the hosts' initial line of protection against the pathogen. Besides LcrV, Y. pestis has LPS as a major component in its outer membrane. This works as a ligand for Toll-like receptor 4 (TLR-4) and plays a basic role in virulence; this is why we used TLR-4 for our study [9] . In traditional vaccinology, microbes are cultured, followed by pathogenic antigen identification, characterization, isolation, inactivation, and reinjection into the host to provoke immune responses [10] . These traditional methods are expensive and time-consuming, and also not appropriate for pathogens like Y. pestis. Advancements in system biology, DNA sequencing, genomics, and proteomics have facilitated a better understanding of organism pathogenicity, contributing multifold to vaccine design. Such techniques have good accuracy and cost-effectiveness and have been successfully applied to many pathogens in the recent past with promising results [11] [12] [13] [14] [15] [16] [17] [18] . Several types of vaccines based on the strategy of epitope predictions have been developed against human pathogens. This includes the development of a potential epitope (AMA-1) based vaccine candidate against Plasmodium vivax [19] . Recently, in 2019, a computationally predicted multi-epitope candidate vaccine was evaluated against Acinetobacter baumannii with potent IgG antibody-specific immune response in mice [20] . Similarly, the utility of computationally predicted B-cell epitopes was also confirmed in diagnostics against Trypsonoma vivax [21] . For instance, such approaches Vaccines 2021, 9, 1327 3 of 19 have been deployed against several human pathogenic viruses, including Ebola virus [22] , Marburg virus [23] , Crimean-Congo hemorrhagic fever virus [24] , and Mokola Rabies virus [25] . Moreover, Ziwei et al. predicted T-cell and B-cell epitopes against SARS-CoV-2 and tested them experimentally, which resulted in a highly immunogenic response in experimental mice [26] . This reflects the potential implications of reverse-vaccinology-based approaches in vaccine design against human pathogenic viruses. Keeping in view the importance of subtractive proteomics in vaccine target mining and reverse-vaccinology-based vaccine design, in the present study, we also aimed to design a multi-epitope subunit vaccine for Y. pestis by screening the whole proteome. Using subtractive proteomics followed by a reverse vaccinology pipeline, potential virulent and highly antigenic proteins from the Y. pestis proteome to design B-cell, HTL, and T-cell epitopes were identified. These epitopes were used in the design of a vaccine construct, which was then modeled and docked with TLR-4 immune receptor and subjected to biophysical analysis. Overall, the results presented in this study are interesting and may be useful to accelerating the development of an effective vaccine to combat this deadly pathogen. The study was initiated by retrieving Yersinia pestis' proteome from UniProt, available under accession number UP000000815 (Consortium, 2014) [27] . Yersinia pestis' proteins that are homologous to the human proteome in structure and function generate autoimmunity and cross-reactions; thus, they are not desired in vaccine development. Therefore, all the proteins were used in NCBI's BLASTp against the human reference genome by setting the E-value threshold to 0.000001 [28] . The non-homologous proteins were then used in redundancy checks to discard duplicated sequences in the proteome. For this, a CD-HIT server was used, keeping the cut-off value at 0.8 (80%). The output file containing only non-redundant proteins was analyzed further in the downward analysis [29] . Essential proteins are vital for the survival of an organism, and their disruption leads to cell death. The Database of Essential Genes (DEG) contains all essential genes of viruses, bacteria, and fungi [30] . Therefore, the non-redundant set of proteins was checked via BLASTp against the DEG, and only proteins with hits (E-value 0.000001) were considered for further analysis. Proteins of small size (<100 amino acids) were manually removed from the essential proteome as they do not harbor potential epitopes. Protein subcellular localization was performed to explore the locations of the essential proteins, as surface proteins of bacterial cells are exposed and easily recognized by the host immune system. We used the CELLO online server, keeping the class bacteria and the type Gram-negative [31] . We prioritized the cell membrane, peripheral membrane, extracellular, and outer membrane proteins for subunit-based vaccine design. On the contrary, cytoplasmic proteins contain hydrophobic pockets for binding to small drug molecules; thus, we removed the cytoplasmic proteins from the vaccine design perspective [17, 18, 32] . Virulent proteins have pathogenic characteristics like adherence to the host cells, invasion, toxin action, and actin-based motility, and they are involved actively in secretion systems [33] . The Virulence Factor Database (VFD) was used to predict the virulent propensity of the proteins filtered in the previous step using PSI/PHI BLAST, and those with >30% identity and >100% bit score were allowed to enter the next phase [34] . Antigenicity is the capacity of proteins (epitopes) to bind with immune cells and generate productive adaptive immune responses [35] . The capacity of a foreign particle (antigen) to attach to or interact with the products of the final cell-mediated response, such as B-cell or T-cell receptors, is referred to as "antigenicity". Antigenicity is the term used to define whether and how well the substance attaches to immune cells, beginning the immune response process. Antigenicity may also be described as the ability of an antigen or a hapten to bind to a B-cell's or T-cell's receptors. To trigger an immune response, an antigen binds to the receptor of an immune cell, such as a B-cell or T-cell. The surface of every antigen or hapten has a unique biological key called an epitope. It is this epitope that binds with a corresponding receptor or antibody. When T-cell receptors interact with a major histocompatibility complex (MHC) molecule, they identify linear amino acid sequences inside a protein antigen, commonly known as epitopes. Hence, it is important to estimate the antigenic potential of each protein and their potential epitopes. The antigenicity of protein was checked via the VaxiJen server at 0.4 threshold [36] . Five proteins (rstB, YPO2385, hmuR, flaA1, and psaB) were prioritized. These proteins have a length of more than 100 amino acids, are virulent and antigenic in nature, and have a molecular weight of more than 100 kbs. The five selected proteins, rstB, YPO2385, hmuR, flaA1, and psaB, were used for the prediction of epitopes. Cytotoxic T lymphocyte (CTL) and helper T lymphocyte (HTL) epitopes were predicted by using two different servers. For CTL epitope prediction, the NetCTL 1.2 server was used, and HTL epitopes were predicted by the IEDB MHC II server [37] . Epitopes from each protein were selected on the basis of the highest combined score among selected epitopes. The E-value was set to 0.75. Seven alleles were selected as an HLA reference set. Epitopes of 15mer in size were selected, three from each protein. The AlgPred server determined the allergenicity of each HTL epitope. A final vaccine construct was designed by linking the epitopes at a ratio of 2:3 CTL and HTL, comprising 8 CTL and 12 HTL. Overlapping HTL epitopes were not used in the construct design. Adjuvant was joined to the CTL epitope peptide at the N-terminal using the EAAAK linker, whereas CTL epitopes were fused by the AAY linker. The last CTL epitope was joined to HTL by a GPGPG linker. Ratio: 2:3 2 × 4 = 8 CTL epitopes 3 × 4 = 12 HTL epitopes The vaccine construct was further scanned for allergenic sequences by analyzing it in the AlgPred 2.0 server (https://webs.iiitd.edu.in/raghava/algpred2/, accessed on 30 December 2020) [38] . The AlgPred server uses different algorithms such as SVM + MAST + ARPs BLAST + IgE epitope to properly evaluate allergenic sequences with up to 85% accuracy. The antigenic nature of the vaccine is very important and is required to boost the immune system. The antigenicity of the vaccine construct was evaluated by using a freely accessible server, VaxiJen v 2.0. The accuracy of this server varies from 70% to 80%. The final vaccine construct was analyzed further in the ProtParam server to estimate its instability index, theoretical pI, half-life, aliphatic index, stability profiling, and hydropathy; the Grand Average and other properties were calculated [39] . The three-dimensional (3D) model of the constructed vaccine was modeled using the Robetta (http://robetta.edu/StructurePrediction/predict/ (accessed on 05 January 2021)) webserver. Robetta is good at predicting query protein 3D structures without close similarity (<30%) of the input sequence [40] . A total of five models were generated, and the one with the best prediction score was chosen. The vaccine construct secondary structure was predicted by the PSIPRED 3.0 server [41] . The tertiary structure of the vaccine was refined using GalaxyRefine to improve its local and global structural quality [42] . The 3D structure of the vaccine construct was then validated by using three freely available online servers: RAMPAGE (http://mordred.bioc.cam.ac.uk//~rapper/rampage. php, accessed on 05 January 2021) [43] , ProSA-web, accessed on https://prosa./services. came.sbg.ac.at/prosa.php, accessed on 05 January 2021) [44] , and ERRAT (http://servicesn. mbi.ucla.edu//ERRAT/, accessed on 05 January 2021) Lengths and Angles [45] . B-cell epitopes trigger humoral immunity to produce antibodies specific to the antigen, and memory cells are significant in memorizing the pathogen and responding on reencounter [46] . Linear and conformational B-cell epitopes were mapped for the subunit vaccine by using the B-cell epitope prediction server of BCPREDS [47] . Conformational B-cell epitopes were predicted by IEDB MHC Ellipro [48] . Molecular docking was used to predict the binding conformation and interactions of the vaccine construct with the TLR-4 immune receptor. Docking was performed using a blind docking approach through the HAWKDOCK server [49] . The server functions by using a geometry-based molecular docking algorithm that gives the best docked intermolecular conformation [49] . A molecular dynamics simulation of the vaccine-TLR complex was performed to check the stability of the docked complex. The AMBER20 simulation package was used [50] . Two steps of energy minimization followed by heating, equilibration, and production were performed. A total 100 ns simulation was performed. Default parameters were used as given in the previous study [11, 13, 51] . The RMSD and RMSF were calculated to check the binding stability and flexibility. Vaccines 2021, 9, 1327 6 of 19 Here, N is the number of atoms, m i is the mass of atom i, x i is the coordinate vector for target atom i, Y i is the coordinate vector for reference atom i, and M is the total mass. If the RMSD is not mass-weighted, all m i = 1 and M = N. The RMSF is a measure of the deviation between the position of particle i and some reference position: where T is the time over which one wants to average and ri ref is the reference position of particle i. This reference position will be the time-averaged position of the same particle i. The codon usage of the vaccine construct was adjusted according to the Escherichia coli system to obtain better expression of the vaccine sequence. For this, Jcat software was used first to generate a reverse transcript of the input vaccine sequence, with calculation of the GC content. An appropriate CAI value (needs to be 0.9+) and GC content (45-70%) are indications of a well-optimized sequence [52] . To carry the vaccine sequence into the expression system, suitable restriction enzyme sites were created to ease the insertion of the vaccine into the vector. SnapGene software was used for cloning purposes [53] . In order to understand the dynamics of the human immune system in response to foreign particles, a server that uses agent-based modeling, C-ImmSim, was used to predict the relationships between the human immune system and the foreign particle [54] . The production of cytokines and other substances like interferon and antibodies was estimated by applying the PSSM method. Moreover, the response for T helper cell 1 and T helper cell 2 (Th1 and Th2) was also predicted with the server's default parameter measure of diversity or Simps Index [55] . The immune simulation uses the following mathematical model to estimate the interaction of the vaccine and therapeutics. The above linear mathematical model was used to model the immune simulation of the designed vaccine [56] . The pathogen proteome comprises a total of 3909 proteins. These proteins were subsequently subjected to the designed vaccine framework to prioritize potential vaccine targets and design a vaccine ensemble. The whole methodological flow is shown in Figure 1 . Y. pestis proteins that were homologous to the human proteome were discarded from the complete proteome; thus, 1021 proteins were filtered as non-homologous. This set of proteins was considered for additional analysis as they are not able to elicit any non-specific responses. In total, 704 proteins were found to be non-redundant and were therefore allowed into further steps of the vaccine design framework. These proteins are singly represented in the proteome and are critical to the pathogen's core functional pathways. Figure 1 . Workflow of the complete methodology, shown in eight steps. First, the selected proteins' antigenicity was checked. Secondly, virulence proteins were screened, giving five in number. The third step was B-and T-cell epitope prediction; the fourth step was the construction of the vaccine; the fifth step was determining the different physicochemical parameters of the vaccine construct, the sixth step was molecular docking; the seventh step was in silico cloning; and the eighth and last step was an immune simulation. Y. pestis proteins that were homologous to the human proteome were discarded from the complete proteome; thus, 1021 proteins were filtered as non-homologous. This set of proteins was considered for additional analysis as they are not able to elicit any non-specific responses. In total, 704 proteins were found to be non-redundant and were therefore allowed into further steps of the vaccine design framework. These proteins are singly represented in the proteome and are critical to the pathogen's core functional pathways. The essential proteins play a central role in pathogen survival and are good candidates for vaccine design. A total of 368 essential proteins were found after excluding nonessential and small-sized proteins. Subcellular localization analysis revealed 73 periplasmic, 46 outer membrane, 152 inner membrane, and 11 extracellular proteins. Such proteins come in close contact with the host during interactions and contain antigenic determinants in vaccine design. Figure 1 . Workflow of the complete methodology, shown in eight steps. First, the selected proteins' antigenicity was checked. Secondly, virulence proteins were screened, giving five in number. The third step was B-and T-cell epitope prediction; the fourth step was the construction of the vaccine; the fifth step was determining the different physicochemical parameters of the vaccine construct, the sixth step was molecular docking; the seventh step was in silico cloning; and the eighth and last step was an immune simulation. The essential proteins play a central role in pathogen survival and are good candidates for vaccine design. A total of 368 essential proteins were found after excluding non-essential and small-sized proteins. Subcellular localization analysis revealed 73 periplasmic, 46 outer membrane, 152 inner membrane, and 11 extracellular proteins. Such proteins come in close contact with the host during interactions and contain antigenic determinants in vaccine design. Virulent proteins stimulate infectious pathways of the host, allowing an efficient response from the immune system to the antigenic region of the proteins. The five selected proteins' (rstB, YPO2385, hmuR, flaA1a, and psaB) amino acid sequences were retrieved from the gene bank by using their gene ids (1175141, 1175217, 1173129, 1173578, and 1174146) for the multi-epitope subunit vaccine design ( Figure 2 ). Virulent proteins stimulate infectious pathways of the host, allowing an efficient response from the immune system to the antigenic region of the proteins. The five selected proteins' (rstB, YPO2385, hmuR, flaA1a, and psaB) amino acid sequences were retrieved from the gene bank by using their gene ids (1175141, 1175217, 1173129, 1173578, and 1174146) for the multi-epitope subunit vaccine design ( Figure 2 ). Antigenic proteins that have the capacity to bind with the immune cells' receptors efficiently were predicted from among the virulent proteins. Among the total, only five highly antigenic proteins were selected, shown in Table 1 . The physiochemical properties of the five selected vaccine targets (rstB, YPO2385, hmuR, flaA1a, and psaB) were then determined. These physiochemical properties, such as virulence, antigenicity, mass, and total number of amino acids, are given in Table 1 . Antigenic proteins that have the capacity to bind with the immune cells' receptors efficiently were predicted from among the virulent proteins. Among the total, only five highly antigenic proteins were selected, shown in Table 1 . The physiochemical properties of the five selected vaccine targets (rstB, YPO2385, hmuR, flaA1a, and psaB) were then determined. These physiochemical properties, such as virulence, antigenicity, mass, and total number of amino acids, are given in Table 1 . Epitope prediction is an essential step in vaccine design as the immune system recognizes and binds to epitopes to activate specific immunity. In total, 72 CTL epitopes were predicted for the five selected proteins; these can be split into 14 epitopes for sensor kinase protein, 10 epitopes for putative exported protein, 26 epitopes for hemin receptor, 9 epitopes for flagellin, and 13 epitopes for chaperone protein PsaB. Only 8 CTL epitopes were selected because of their high MHC binding score and non-allergenic nature. The selected set of epitopes comprises two epitopes from sensor kinase protein, one from Vaccines 2021, 9, 1327 9 of 19 putative exported protein, two from hemin receptor, one from flagellin, and two from chaperone protein PsaB, as tabulated in Table 2 . A total of 12 HTL epitopes were selected, giving a 2:3 ratio; three of the HTL epitopes were selected from sensor kinase protein, two from putative exported protein, three from hemin receptor, and two each from flagellin and chaperone protein PsaB, as shown in Table 3 . B-cells, upon maturation, produce antibodies and are essential for long-lasting immunity against a particular antigen. In total, 16 linear B-cell epitopes with score >0.6 were selected. Likewise, the ElliPro suite for the prediction of conformational B-cell epitopes with a score of 0.7 was used. B-cell epitopes were predicted from the selected five proteins. A total of 17 epitopes were predicted from Chaperone protein, out of which 4 epitopes were chosen on the basis of a score above 0.9; a total of 24 epitopes were predicted from Histidine kinase, out of which 2 epitopes were selected. Twenty-eight epitopes were predicted from the Hemin receptor, out of which two were selected with a high score. From 24 epitopes predicted from Flagellin protein, 4 epitopes were chosen, and of 21 epitopes predicted from Putative exported protein, 5 epitopes were selected on the basis of a score above 0.9. The predicted linear B epitopes are given in Table 4 . The vaccine ensemble was generated by joining epitopes at a ratio of 2:3; in total, 8 CTL and 12 HTL epitopes were added using linkers. The overlapping HTL epitopes were excluded during vaccine design. The adjuvant was joined to CTL epitopes by using EAAAK linkers at the N terminal, whereas CTL epitopes were joined via AAY linker. The last CTL epitope was joined with HTL by GPGPG linkers, and GPGPG linkers were also used among HTL epitopes, as given in Figure 3 . The vaccine ensemble was generated by joining epitopes at a ratio of 2:3; in total, 8 CTL and 12 HTL epitopes were added using linkers. The overlapping HTL epitopes were excluded during vaccine design. The adjuvant was joined to CTL epitopes by using EAAAK linkers at the N terminal, whereas CTL epitopes were joined via AAY linker. The last CTL epitope was joined with HTL by GPGPG linkers, and GPGPG linkers were also used among HTL epitopes, as given in Figure 3 . By using the Robetta server, 3D structures of the vaccine were constructed. The modulated structure is shown in Figure 4 . A comparative modeling approach was used, and multiple templates were selected for the 3D structure model design. The confidence levels of the templates for 3D structure design were all above 97.9%. The 3D structure predicted By using the Robetta server, 3D structures of the vaccine were constructed. The modulated structure is shown in Figure 4 . A comparative modeling approach was used, and multiple templates were selected for the 3D structure model design. The confidence levels of the templates for 3D structure design were all above 97.9%. The 3D structure predicted by Robetta, as shown, possesses alpha-helices, beta-sheets, and loops. The visual structure analysis shows that it has proper folding. Figure 3 . Arrangement of the final vaccine construct design comprising CTL and HTL epitopes using three different linkers: AAY, GPGPG, and EAAAK. By using the Robetta server, 3D structures of the vaccine were constructed. The modulated structure is shown in Figure 4 . A comparative modeling approach was used, and multiple templates were selected for the 3D structure model design. The confidence levels of the templates for 3D structure design were all above 97.9%. The 3D structure predicted by Robetta, as shown, possesses alpha-helices, beta-sheets, and loops. The visual structure analysis shows that it has proper folding. The results of refinement analysis of the vaccine are tabulated in Table 5 . It was concluded that Model 2 was the best model for docking. Model 2 ( Figure 4 ) had a high GDT-HA score of 0.9583, and its root mean square deviation value was 0.402. The MolProbity score was 1.400, the clash score was improved, poor rotamers were 0.6, and the Rama favored level was calculated to be 98.5%. The results of refinement analysis of the vaccine are tabulated in Table 5 . It was concluded that Model 2 was the best model for docking. Model 2 ( Figure 4 ) had a high GDT-HA score of 0.9583, and its root mean square deviation value was 0.402. The MolProbity score was 1.400, the clash score was improved, poor rotamers were 0.6, and the Rama favored level was calculated to be 98.5%. Table 5 . The tertiary structure of the vaccine was refined using the GalaxyRefine server, and Model 2 was selected due to its good Rama favored score of 98.54 and low mol probity of 1.40. The 3D structure was validated using three available online servers. The results of analysis using the Ramachandran and PROSA-web servers are given in Figure 5 . Panel (A) shows the Z-score (−4.78) predicted by PROSA-web, while in panel (B), the Ramachandran plot shows the amino acid distribution in the favored region as 91.8%, with 7.9% allowed and 0.3% disallowed. These results suggest that our predicted structure has a properly folded topology and could be used for further processes. The 3D structure was validated using three available online servers. The results of analysis using the Ramachandran and PROSA-web servers are given in Figure 5 . Panel (A) shows the Z-score (−4.78) predicted by PROSA-web, while in panel (B), the Ramachandran plot shows the amino acid distribution in the favored region as 91.8%, with 7.9% allowed and 0.3% disallowed. These results suggest that our predicted structure has a properly folded topology and could be used for further processes. The vaccine construct is non-allergenic with a score of −1.52; isoelectric point of 9.34; molecular weight of 49.70 kDa; half-life of >30 h (mammalian reticulocytes, in vitro), >20 h (yeast, in vivo), or >10 h (E. coli, in vivo); instability index of 35.03; GRAVY score of −0.225; and aliphatic index of 73.14. All these parameters strongly support the construct as a good vaccine candidate. From an antigenicity point of view, the vaccine is a good antigen (score 0.7548). The vaccine molecule is also non-allergic (score 1.523). The PSIPRED webserver was used to predict the secondary structure of the vaccine construct shown in Figure 6 . The secondary structure of the vaccine revealed 35.7% alpha helices (including 4% transmembrane helices), 14.3% beta sheets, and 50% coils. The vaccine construct is non-allergenic with a score of −1.52; isoelectric point of 9.34; molecular weight of 49.70 kDa; half-life of >30 h (mammalian reticulocytes, in vitro), >20 h (yeast, in vivo), or >10 h (E. coli, in vivo); instability index of 35.03; GRAVY score of −0.225; and aliphatic index of 73.14. All these parameters strongly support the construct as a good vaccine candidate. From an antigenicity point of view, the vaccine is a good antigen (score 0.7548). The vaccine molecule is also non-allergic (score 1.523). The PSIPRED webserver was used to predict the secondary structure of the vaccine construct shown in Figure 6 . The secondary structure of the vaccine revealed 35.7% alpha helices (including 4% transmembrane helices), 14.3% beta sheets, and 50% coils. Vaccines 2021, 9, x FOR PEER REVIEW 13 Figure 6 . Secondary structure of the vaccine, which shows 35.7% helices, 14.3% B sheets, and 50% coils. Solution 2 of vaccine-TLR-4 was selected as the best-docked complex. The select was based on FireDock results with global energy of −15.47, attractive van der Waals fo of −11.27, repulsive van der Waals force of −4.39, and atomic contact energy of −3.98. interactions given in Figure 7 show that the vaccine candidate is significantly involved salt bridges, hydrogen bonding, and disulfide bridges. Solution 2 of vaccine-TLR-4 was selected as the best-docked complex. The selection was based on FireDock results with global energy of −15.47, attractive van der Waals force of −11.27, repulsive van der Waals force of −4.39, and atomic contact energy of −3.98. The interactions given in Figure 7 show that the vaccine candidate is significantly involved in salt bridges, hydrogen bonding, and disulfide bridges. Solution 2 of vaccine-TLR-4 was selected as the best-docked complex. The selection was based on FireDock results with global energy of −15.47, attractive van der Waals force of −11.27, repulsive van der Waals force of −4.39, and atomic contact energy of −3.98. The interactions given in Figure 7 show that the vaccine candidate is significantly involved in salt bridges, hydrogen bonding, and disulfide bridges. A molecular dynamics simulation of the vaccine-TLR complex revealed that the complex is stable with no significant convergence, while the RMSF results show that the vaccine structure possesses higher flexibility than the TLR. Our results suggest that our complex remained stable during the 100 ns simulation, and this signifies our vaccine's efficiency. The RMSD and RMSF values of the complex are given in Figure 8 . A molecular dynamics simulation of the vaccine-TLR complex revealed that the complex is stable with no significant convergence, while the RMSF results show that the vaccine structure possesses higher flexibility than the TLR. Our results suggest that our complex remained stable during the 100 ns simulation, and this signifies our vaccine's efficiency. The RMSD and RMSF values of the complex are given in Figure 8 . JCat software was utilized to quantify the expression of vaccine in Escherichia coli O6: k15:H31 (strain 536/UPEC). The vaccine is 1400 nucleotides long; its CAI is 1.0 with a GC content of 55.19%. This indicates a better expression level of the vaccine in E. coli sO6: JCat software was utilized to quantify the expression of vaccine in Escherichia coli O6: k15:H31 (strain 536/UPEC). The vaccine is 1400 nucleotides long; its CAI is 1.0 with a GC content of 55.19%. This indicates a better expression level of the vaccine in E. coli sO6: k15:H31 (strain 536/UPEC). A GC content range of 35% to 70% is reported to be appropriate for better expression. Restriction sites like EcoR1 and Pas1 were added to the 5 and 3 ends of the optimized vaccine sequence and cloned in pET-53-Dast vector as given in Figure 9 . JCat software was utilized to quantify the expression of vaccine in Escherichia coli O6: k15:H31 (strain 536/UPEC). The vaccine is 1400 nucleotides long; its CAI is 1.0 with a GC content of 55.19%. This indicates a better expression level of the vaccine in E. coli sO6: k15:H31 (strain 536/UPEC). A GC content range of 35% to 70% is reported to be appropriate for better expression. Restriction sites like EcoR1 and Pas1 were added to the 5′ and 3′ ends of the optimized vaccine sequence and cloned in pET-53-Dast vector as given in Figure 9 . An immune simulation study was performed on web server C-ImmSim 1.0, which predicts possible dynamics of the host immune system in response to vaccine administration. The host immune system dynamics are summarized in Figure 10 . Panel A shows the humoral immunity response to the antigen and reported immunoglobulin response: IgG, IgG1, IgG2, and IgM produced against the vaccine. According to the literature, most plague patients have been reported with igG1 and IgM. In parts B, C, and D, production of B-and T-cells in cells per cubic millimeter is shown to be very high in the presence of vaccine; the B-cell concentration is high and remains so because they produce antibiotics, staying active for a longer time. The graph for the T-cells shows a downward trend after a few days. Panel E shows macrophages, and panel F shows the interleukin and cytokine concentration, which was very high for the vaccine inside immune cells for the first few days. plague patients have been reported with igG1 and IgM. In parts B, C, and D, production of B-and T-cells in cells per cubic millimeter is shown to be very high in the presence of vaccine; the B-cell concentration is high and remains so because they produce antibiotics, staying active for a longer time. The graph for the T-cells shows a downward trend after a few days. Panel E shows macrophages, and panel F shows the interleukin and cytokine concentration, which was very high for the vaccine inside immune cells for the first few days. For infectious disease, antibiotic controls are traditionally considered; however, due to their decreasing effectiveness, the scientific community has begun the search for a new approach, including vaccines that could replace antibiotics. Vaccinations were initially used for infectious disease control and played a significant role in this, but recently, due to our greater dependency on them, scientists have sought more effective, cheap, and safe techniques for designing vaccines. Because traditional vaccine technology harbors the risk of virulence reoccurrence and is not very cost-effective, a new technology, called immuneinformatics, has been developed, requiring knowledge of genomics and proteomics to bring new vaccines to the market. This uses only a particular part of the pathogen (epitope) which can only boost the immune system. Immune informatics technology uses computational software, online servers, and databases. It gives good-quality results and is low cost, with good safety and greater effectiveness. Several types of vaccines based on For infectious disease, antibiotic controls are traditionally considered; however, due to their decreasing effectiveness, the scientific community has begun the search for a new approach, including vaccines that could replace antibiotics. Vaccinations were initially used for infectious disease control and played a significant role in this, but recently, due to our greater dependency on them, scientists have sought more effective, cheap, and safe techniques for designing vaccines. Because traditional vaccine technology harbors the risk of virulence reoccurrence and is not very cost-effective, a new technology, called immuneinformatics, has been developed, requiring knowledge of genomics and proteomics to bring new vaccines to the market. This uses only a particular part of the pathogen (epitope) which can only boost the immune system. Immune informatics technology uses computational software, online servers, and databases. It gives good-quality results and is low cost, with good safety and greater effectiveness. Several types of vaccines based on the strategy of epitope predictions have been developed against human pathogens. This includes the development of a potential epitope (AMA-1) based vaccine candidate against Plasmodium vivax [19] . Recently, in 2019, a computationally predicted multi-epitope candidate vaccine was evaluated against Acinetobacter baumannii with potent IgG antibody-specific immune response in mice [20] . Similarly, the utility of computationally predicted B-cell epitopes was also confirmed in diagnostics against Trypsonoma vivax [21] . For instance, such approaches have been deployed against several human pathogenic viruses, including Ebola virus [22] , Marburg virus [23] , Crimean-Congo hemorrhagic fever virus [24] , and Mokola Rabies virus [25] . Moreover, Ziwei et al. predicted T-cell and B-cell epitopes against SARS-CoV-2 and tested them experimentally, which resulted in a highly immunogenic response in experimental mice [26] . This reflects the potential applications of in silico reversevaccinology-based approaches in vaccine design against human pathogenic viruses. The proteome is the protein factory of an organism and contains the proteins of traits present in an organism. This is why we exploited the proteome of Yersinia pestis to design a vaccine. Highly antigenic and immunogenic proteins rstB, YPO2385, hmuR, flaA1, and psaB were selected using online servers, and the epitopes for adaptive immune cells, B lymphocytes, and T-cell lymphocytes were predicted. A foreign organism first deals with innate immune cells, so we added adjuvants that interact with the TLR of innate immune cells. Epitopes were predicted through online servers, which showed the binding potency of the epitopes with MHC-I and MHC-II of T-cells, and predicted B-cell linear epitopes that interacted with immune B-cells. The designed vaccine's physiochemical properties were evaluated. Although the designed vaccine candidate exhibits non-allergenic features, the "cytokine storm" phenomenon cannot be avoided in infectious diseases, but this is rare. Based on recent developments in the field, new therapies will likely move away from targeting single soluble mediators and focus more generally on inflammatory cascades. However, what signaling cascades are targeted will depend both on the infection and on the targeted patient population. Current experimental anti-inflammatory strategies being tested in animal models include COX inhibitors, PAR2 agonists, sphingosine 1 phosphate modulators, modulation of production/signaling of oxidized phospholipids, Slit-Robo4 pathway modulators, and TLR-4 antagonists. These anti-inflammatory agents can be supplied as supplementary agents to overcome the problem of cytokine storm [57] . Furthermore, the 2D structure was designed using the PSI PRED server, and the 3D structure was predicted by the Robetta server and refined by the GALAXY refinement server. Further corrections were made using Chimera software. Then, more of the interaction with the TLR of innate immune cells was checked by docking analysis using PatchDock. Stability and fluctuations were achieved by MD simulation. High levels of expression were needed. Thus, the Jcat server was used as it automatically performs reverse transcription, CAI calculations, and calculations of GC contents for codon optimization. E. coli was used as an expression vector. Vaccine stability is required for better expression in E. coli; therefore, a bisulfate bond was added with the vaccine construct. Immuno-informatics approaches were used in the present study to give a stable, very effective, and highly expressible vaccine in the E. coli expression vector. Clinical trials are required to check its effectiveness in human beings. The designed MEVC construct can be further used in an experimental setup to evaluate the in vitro and in vivo potency, safety, and efficacy of the vaccine construct. The plague is one of the deadliest diseases remaining in the present day, and the agent of this disease is Yersinia pestis. In the study presented herein, we used immune informatics techniques to design a vaccine from the epitopes of Yersinia pestis for better control. Such a vaccine is capable of providing immunity for both types of immune systems. Different techniques were used, such as molecular docking, simulation for thermal stability, and in silico cloning, to verify its thermal stability inside the human body, its better expression in vectors for mass production, and its immune response boost inability. The current study has the only limitation of lacking experimental testing, which can be performed by experimental scientists to confirm the validity of the designed candidate and its use for therapeutic purposes. Data Availability Statement: All the data will be provided on reasonable request. Pneumonic plague: The darker side of Yersinia pestis Plague, camels, and lice Plague: Recognition, treatment, and prevention Human Neutrophil Isolation and Degranulation Responses to Yersinia pestis Infection Yersinia pestis Pla protein thwarts T cell defense against plague Genome sequence of Yersinia pestis KIM Interaction of Yersinia pestis with macrophages: Limitations in YopJ-dependent apoptosis Insight into bacterial virulence mechanisms against host immune response via the Yersinia pestis-human protein-protein interaction network Interaction between Yersinia pestis and the host immune system Reverse vaccinology and vaccines for serogroup B Neisseria meningitidis. In Hot Topics in Infection and Immunity in Children II Computational vaccinology based development of multi-epitope subunit vaccine for protection against the Norovirus infections Immunoinformatic and systems biology approaches to predict and validate peptide vaccines against Epstein-Barr virus (EBV) Combined drug repurposing and virtual screening strategies with molecular dynamics simulation identified potent inhibitors for SARS-CoV-2 main protease (3CLpro) Computational identification, characterization and validation of potential antigenic peptide vaccines from hrHPVs E6 proteins using immunoinformatics and computational systems biology approaches Phylogenetic analysis and structural perspectives of RNA-dependent RNA-polymerase inhibition from SARs-CoV-2 with natural products Structural Insights into the mechanism of RNA recognition by the N-terminal RNA-binding domain of the SARS-CoV-2 nucleocapsid phosphoprotein Immunoinformatics approaches to explore Helicobacter Pylori proteome (Virulence Factors) to design B and T cell multi-epitope subunit vaccine Immunoinformatics and structural vaccinology driven prediction of multi-epitope vaccine against Mayaro virus and validation through in-silico expression Identification of a highly antigenic linear B cell epitope within Plasmodium vivax apical membrane antigen 1 (AMA-1) Design and evaluation of a multi-epitope assembly peptide vaccine against Acinetobacter baumannii infection in mice A comparative in silico linear B-cell epitope prediction and characterization for South American and African Trypanosoma vivax strains In silico designed ebola virus T-cell multi-epitope DNA vaccine constructions are immunogenic in mice Vaccinomics strategy for developing a unique multi-epitope monovalent vaccine against Marburg marburgvirus Towards the first multi-epitope recombinant vaccine against Crimean-Congo hemorrhagic fever virus: A computer-aided vaccine design approach Epitope-based peptide vaccine design against Mokola rabies virus glycoprotein G utilizing in silico approaches T and B cell Epitope analysis of SARS-CoV-2 S protein based on immunoinformatics and experimental research UniProt Knowledgebase: A hub of integrated protein data Fast and accurate NCBI BLASTP: Acceleration with multiphase FPGA-based prefiltering Accelerated for clustering the next-generation sequencing data DEG: A database of essential genes CELLO2GO: A web server for protein subCELlular LOcalization prediction with functional gene ontology annotation Immunogenomics Guided Design of Immunomodulatory Multi-Epitope Subunit Vaccine against the SARS-CoV-2 new Variants, and its Validation through in Silico Cloning and Immune Simulation Yersinia pestis and plague: An updated view on evolution, virulence determinants, immune subversion, vaccination, and diagnostics VFDB: A reference database for bacterial virulence factors Genetic diversity and antigenicity analysis of Streptococcus pneumoniae pneumolysin isolated from children with pneumococcal infection VaxiJen Dataset of Bacterial Immunogens: An Update Large-scale validation of methods for cytotoxic T-lymphocyte epitope prediction AlgPred: Prediction of allergenic proteins and mapping of IgE epitopes Structure prediction and functional characterization of secondary metabolite proteins of Ocimum Prediction of CASP6 structures using automated Robetta protocols The PSIPRED protein structure prediction server GalaxyWEB server for protein structure prediction and refinement Data set for phylogenetic tree and RAMPAGE Ramachandran plot analysis of SODs in Gossypium raimondii and G. arboreum. Data Brief ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins Limitations of structure evaluation tools errat Predicting linear B-cell epitopes using string kernels ElliPro: A new structure-based tool for the prediction of antibody epitopes PatchDock and SymmDock: Servers for rigid and symmetric docking An overview of the Amber biomolecular simulation package Structural and free energy landscape of novel mutations in ribosomal protein S1 (rpsA) associated with pyrazinamide resistance JCat: A novel tool to adapt codon usage of a target gene to its potential expression host TSPO PET for glioma imaging using the novel ligand 18 F-GE-180: First results in patients with glioblastoma Immune system simulation online Comparing In Silico Results to In Vivo and Ex Vivo of Influenza-Specific Immune Responses after Vaccination or Infection in Humans An Der Heiden, U. A basic mathematical model of the immune response Cytokine storms in infectious diseases Acknowledgments: The computations were partially performed at the Center for High-Performance Computing, Shanghai Jiao Tong University. We acknowledge their help. The authors declare no conflict of interest.