key: cord-034402-64xz1j9i authors: Srivastava, Shivani; Sharma, Suraj Kumar; Srivastava, Vivek; Kumar, Ajay title: Proteomic Exploration of Listeria monocytogenes for the Purpose of Vaccine Designing Using a Reverse Vaccinology Approach date: 2020-10-29 journal: Int J Pept Res Ther DOI: 10.1007/s10989-020-10128-1 sha: doc_id: 34402 cord_uid: 64xz1j9i Listeriosis is a major foodborne infection provoked by a bacterium known as Listeria monocytogenes. It is one of the predominant causes of death in pregnant women, infants, and immunocompromised persons. Despite such fatal effects, until now there is no proper medication or drug available for such a serious foodborne infection. One of the most promising ways to deal with this challenge is vaccination. This present study aims at the prediction of B cell epitopes for subunit vaccine designing against Listeria monocytogenes using a reverse vaccinology approach. Among screened out 299 epitopes of strain F2365 of Listeria monocytogenes, based on the VaxiJen score, the top 20 epitopes were selected. 3D modeling of epitopes and alleles was generated by PEPstrMOD and Swiss Model respectively. Molecular docking reveals 4 epitopes viz., MKFLFPLKL, CEETFGIRL, FLKIDPPIL, and VRHHGGGHK based on binding energy. All 4 epitopes were investigated for non-toxicity, binding affinity, and population coverage. After vigorous investigation, epitope FLKIDPPIL was anticipated as the best vaccine contender. The stability of the FLKIDPPIL-HLA DRB1 _0101 complex was proved by performing the simulation. Here, predicted peptide through the Insilico approach may become a potential remedy against listeriosis, after the wet-lab approach and clinical trials. Changing food habits, advancement in technology regarding the preservation of food products for a longer time, and the ability of microorganisms to grow in adverse conditions are leading to the emergence of the foodborne infection, known as Listeriosis. The genus Listeria consists of seventeen species. Only the three hemolytic species viz., Listeria monocytogenes, Listeria seeligeri, and Listeria ivanovii are considered pathogens. Of these, Listeria monocytogenes is consistently pathogenic and is involved in foodborne outbreaks of listeriosis (Abdelhamed et al. 2019 ). Based on Gram-staining, Listeria monocytogenes comes under the category of Gram-positive. It shows extreme resistance in conditions like very high temperatures or very low temperatures. These bacteria have a rod-like shape and can form small chains (Sallami et al. 2006) . Listeria monocytogenes mainly affects women who are pregnant, infants, elders above 65 years of age, and immunocompromised people (CDC 2019) . Foodborne infection in humans occurs through the consumption of contaminated foods, particularly unpasteurized milk, soft cheeses, vegetables, and prepared meat products. Listeria monocytogenes show completely different behavior in comparison to all other pathogens that cause food contamination. It can multiply at low temperatures in contaminated food. It can be easily transmitted between pregnant women and her newborn either at the time of pregnancy or during delivery (WHO 2019). Pyrexia, cough, cold, headache, and body ache, etc. are the usual symptoms experienced by the patients (Department of Health 2017). Worldwide many countries where food production takes place in absence of proper and better microbiological vigilance and where the percentage of immunocompromised persons are immensely high, Listeria monocytogenes loomed as one of the dominant foodborne pathogens (Thomas et al. 2020 ). Thus, poor surveillance during the production process affects approximately 1600 people every year, and around 260 experience the afterlife (CDC 2019). Listeria monocytogenes consists of two genes viz., chiA and chiB. These two genes play an important role in virulence. A regulatory factor hfq plays a very important role in the formation of biofilm, colony formation, and virulence (Yao et al. 2018) . The Zipper is the name of the mechanism through which Listeria monocytogenes get access to the host cell. In this process, ligands on the surface of bacteria communicate with receptors of the host cells. Internalin A and Internalin B are the ligands on the bacterial surface and E-cadherin and Met are the receptors on host cells with which bacterial ligands interact. This collaboration leads to the rearrangement of actin filaments and invasion of bacteria (Hamon et al. 2006) . When Internalin B-Met interacts together, processes like ubiquitinylation and autophosphorylation takes place (Veiga and Cossart 2005) . In the year 2018, Australia had witnessed around 20 cases of listeriosis between January to April. This minor outbreak had faced around 7 deaths and a single spontaneous abortion (WHO 2018a, b) . National Institute of Communicable diseases (NICD) has proclaimed 978 listeriosis cases between 2017 and 2018 from all provinces of South Africa, Gauteng, Western Cape, and KwaZulu-Natal were mainly hit by this fatal disease known as listeriosis. Around 78% of cases have been reported from the above-mentioned places of South Africa. Out of 674 affected people, 27 have faced death. All these data revealed about the threat of this bacteria and its effect on mankind and society. The percentage of infants that get affected during an outbreak is around 42% (WHO 2018a, b) . Pregnant women can easily get infected with listeriosis through the placenta, still, the establishment of neurolisteriosis is completely occasional. Listeriosis infection in pregnant women is because of the alliance of the quashed immune system and the specificity of bacteria for the placenta (Charlier et al. 2017) . Even after the birth of infants infected with bacteria Listeriosis monocytogenes endurance is possible only with the help of Extracorporeal membrane oxygenation (ECMO) (Lee et al. 2019 ). According to a report of WHO, in India miscarriages and other pregnancy-related disorders is mainly the result of foodborne infection known as listeriosis. Listeriosis is still under-reported in many countries. The ability of Listeria monocytogenes to survive even in harsh conditions is one of the major threats regarding the outbreak of the disease. High fatality rate and frequent outbreaks demand the designing of a vaccine against Listeria monocytogenes, by using the immunoinformatics approach. This study is mainly based on the anticipation of B cell epitopes for the utility of vaccine designing against listeriosis. Previously, a study regarding computational identification and characterization of epitopes has been carried out in the case of the Zika virus, Nipah virus, and bacteria like E. coli (Sharma et al. 2020; Kaushik 2019; Khan et al. 2019) . Considering this approach in this research work, all proteins except hypothetical, putative, and non-structural were retrieved from the UniProtKB database. A potential epitope must not possess any allergic property; therefore, first and foremost allergenicity was checked by using the AlgPred server. NETMHCII 2.3 and VaxiJen server was used to identify B cell epitopes that could bind to MHC II molecules with great stability. Only the top 20 epitopes were selected for further exploration. This selection was done based on the VaxiJen score. 3D modeling of both the epitopes and alleles was performed using PEPstrMOD and Swiss Model. Epitope-allele pair having low binding energy should be selected for the next sequential refining. To do this, molecular docking was performed using AutoDock Vina software. Next to check toxicity, binding affinity, and population coverage Toxin Pred, MHC Pred, and immune epitope database tools were used. The stability of the epitope-allele complex was substantiated by simulation studies. The strategy of the development of subunit vaccines has an upper hand in comparison to traditional vaccines. These next-generation vaccines are extremely specific in eliciting the immune system of the host, can be produced easily in large quantities, and at a comparatively moderate cost. Moreover, peptides consisting of epitopes can be manufactured, purified, and processed easily (Poland et al. 2011 ). For computational identification and characterization of epitopes for the preparation of subunit vaccine designing, complete proteome analysis of Listeria monocytogenes F2365 strain (GenBank accession number AE017262.2) was performed using the UniProtKB database. In comparison to other serotype strains, Listeria monocytogenes strain F2365 belongs to the 4b serotype group and multiplicates more rapidly in monocytes or macrophages (Hasebe et al. 2017 ). Presence of a virulence factor viz. ListeriolysinS (LLS) in the F2365 strain accelerates infection in the intestine and other organs (Quereda et al. 2016) . Listeria monocytogenes F2365 strain is a member of lineage I and comprises a factor known as Internalin B which plays a crucial role in nonpregnant infected animals (Quereda et al. 2018) . All these remarkable features contribute to the pathogenicity of this strain and hence lead to its selection for the study. Excluding hypothetical, putative, and non-structural proteins total of 529 proteins were registered in the UniprotKB database, derived from the different research literature. All these sequences were saved in the FASTA format for further examination. The length of the genome of the F2365 strain is 3,021,822 bp, with GC content of 37.9% approximately (Briers et al. 2011 ). One of the most eminent features in epitope-based vaccine design is that the particular epitope must elicit an immune response only against the target pathogen. Taking this point into consideration, the screened epitope must be non-allergen and thus retrieved proteins were differentiated into allergens and non-allergens by using the AlgPred server (Saha and Raghava 2006) This server segregates non-allergens from allergens and − 0.4 was selected as the cut-off value. Anticipation was done with high accuracy along with sensitivity and specificity of 88.87% and 81.86% respectively. Non-allergens were chosen for another characterization and exploration of antigenic sites for the utility of vaccine designing from the proteome of Listeria monocytogenes. B cell epitopes are typical peptide remnants that bind to the immunoglobulin and thus it becomes immensely important to screen out such epitopes from complete proteome sequence. To accomplish this objective, NETMHCII 2.3 server was used (Jensen et al. 2018) . By making use of artificial neural networks, this server predicts the binding of B cell epitopes with HLA alleles. In this study, three alleles viz. DRB1_0101, DRB1_0701, and DRB1_1301 and locus HLA-DR was chosen. The peptide length was taken at 9 with a threshold set to − 99.9. The potential B cell epitopes were subjected to the Vaxi-Jen server to select those candidates that can strongly bind with MHC II molecules (Doytchinova and Flower 2007) . Only epitopes with a score greater than or equal to 1.1 can bind with MHC II molecules with extreme affinity and be selected. To further proceed with the reverse vaccinology approach, only the top 20 peptides or antigenic sites were chosen. This selection was based on the VaxiJen score. Following allergenicity and prediction of B cell epitopes, modeling of both epitopes and HLA alleles was performed. For the generation of the 3D structure of the selected epitopes, the PEPstrMOD server was used. It offers exclusive advantages to the users to predict the structures of peptides having natural residues, some modified residues, posttranslational modifications, etc. (Singh et al. 2015) . In this research work, filtered epitopes were modeled and saved in the Protein Data Bank (PDB) format for the next sequential investigation. The first fully automated protein homology modeling server known as the Swiss model was used for modeling of HLA alleles (Waterhouse et al. 2018) . The building of models using this server requires four sequential steps. These 4 steps comprise of template selection, its alignment with the target sequence, model building, and its evaluation. In this study 3D structures of three HLA alleles have been performed viz., DRB1_0101, DRB1_0701, and DRB1_1301. To better understand the relationship between anticipated epitopes and their respective alleles, AutoDock Vina software was used to perform molecular docking. It helps us to interpret the synergy between antigenic sites and their corresponding alleles (Trott and Olson 2010) . One of the prerequisites before performing docking is certain modifications both in ligand as well as the receptor, which was performed by AutoDock MGL tools. HLA alleles were selected as receptors viz., DRB1_0101, DRB1_0701, and DRB1_1301. 4AH2, 3C5J, and 6CQL are the crystal structure of these receptors and were retrieved from the Research Collaboratory for Structural Bioinformatics (RCSB) protein data bank. Molecules of water were removed from these receptors and polar hydrogen as well Kollman charges were added to the structure. After modification, the molecule was saved in pdbqt format. Changes were also performed in all 20 ligands and were saved in pdbqt files. All these alterations were performed by AutoDock MGL tools. To perform molecular docking in AutoDock Vina software, 40, 40, 40 were taken as grid box dimensions and energy was calculated at 0 Å. The docking result can be analyzed by a visualization tool called PyMol. 4 epitopes were selected for succeeding rounds of analysis based on negative binding energies where Low binding energy implies good stability. To evaluate the non-toxicity behavior of epitopes Toxin Pred server was used (Gupta et al. 2013) . It is based on machine learning techniques and quantitative matrix scores. Along with toxicity prediction, calculation of physicochemical properties is one of the most notable features of this server. MHC Pred Server was used to vaticinate the binding affinity of epitopes with MHC II molecules. MHC Pred is composed of several models based on structures and its activity, a sturdy multivariate statistical method. Results with articulated by giving IC 50 values (Guan et al. 2003) . IC 50 values less than 500 are considered to be good binders and were chosen for the next and last analysis. Because of the exceptionally heterogeneous behavior of HLA alleles, their frequency of expression varies greatly across the globe, and therefore Population coverage analysis becomes the utmost important step in computational vaccine designing. It was performed using the Immune Epitope Database (IEDB) Population Coverage analysis tool (Bui et al. 2006 ). It is extremely essential to understand the stability of the peptide-allele complex and to analyze that in this research work MD Web server was used (Hospital et al. 2012) . Simulation of 10 ns with an output frequency of 500 steps was set to equilibrate the system. Coarse-grained Brownian dynamics were analyzed for trajectory and output was given in the terms of Root mean square deviation (RMSD) and B-factor values. Both RMSD and B-factor plots corroborate the stability of epitope-allele complex. With time, the world has acknowledged extreme advancement in medicine and technology thus combating some deadly diseases, but still, diseases like listeriosis were left unnoticed. Despite several outbreaks in different parts of the world, there is no legitimate treatment or drug or vaccine available for it. Therefore, it becomes extremely important to predict and characterize some potential vaccine contenders that can evoke a strong immune response and this study is one such step in this direction. Here we have used computational tools to predict B cell epitopes that can elicit an immune response. The first requirement in the reverse vaccinology approach of vaccine designing is to eliminate all non-allergic proteins from a complete proteome set of bacteria, Listeria monocytogenes. The AlgPred server was used to predict allergenicity of retrieved proteins, to get the most capable subunit vaccine candidate. A total of 529 protein sequences of Listeria monocytogenes F2365 strain was retrieved from the UniProtKB database (excluding hypothetical, putative, and non-structural proteins) and were saved in the FASTA format for further analysis. After examination by the AlgPred server, out of 529, a total of 172 proteins were proved to be non-allergens (Table 1) . The result has been summarized in Table 1 . Table 1 consists of protein ID, protein names, and scores of all non-allergens. Non-allergic proteins were analyzed further by using NetMHC II 2.3 server. By selecting peptide lengths 9 and threshold value − 99.9. B cell epitopes were selected. These chosen epitopes were next investigated by the VaxiJen server and the cut-off value was 1.1 Å total of 299 epitopes were found to bind with MHC II molecules ( Table 2 ). All 299 epitopes have a VaxiJen score of ≥ 1.1 and can bind with MHC II molecules with great stability. Among these epitopes, the majority of epitopes were found to bind with DRB1_1301. Based on the high VaxiJen score, among 299 epitopes, only the top 20 epitopes were selected for modeling. The generation of 3D structures of epitopes was performed by PEPstrMOD. 3D modeling of the HLA allele's viz. DRB1_0101, DRB1_0701, and, DRB1_1301 were performed by the Swiss model ( Fig. 1) . For the generation of tertiary structures of DRB1_0101, DRB1_0701 and, DRB1_1301 alleles, proteins having PDB ID 4AH2, 3C5J, and 6CQL were used as templates, respectively. All tertiary structures of HLA alleles were visualized by the PyMOL visualization tool. 3D models have been represented in Fig. 1 . AutoDock Vina software was used to perform molecular docking between 20 nonallergic and antigenic epitopes with their respective alleles. The lowest binding energy was obtained for epitope FLKIDPPIL-DRB1_0101 (− 7.3 kcal/ mol) and the highest binding energy was obtained for epitopes MKGQAGSKK-DRB1_1301 (− 5.1 kcal/mol). As low binding energies imply, high stability of the complex, therefore 4 epitopes based on low binding energy was selected viz., CEETFGIRL, MKFLFPLKL, FLKIDP-PIL, and VRHHGGGHK (Table 3 ). The stable complex of CEETFGIRL-3C5J shows the energy of − 6.7 kcal/mol and 6 hydrogen bonds ( Fig. 2 ) Complexes viz. MKFLFPLKL-4AH2 and FLKIDPPIL-4AH2 shows binding energy of − 6.9 kcal/mol and − 7.3 kcal/mol along with 2 and 6 hydrogen bonds respectively (Figs. 3 and 4). The energy of − 6.7 kcal/mol and 6 hydrogen bonds was shown by epitope VRHHGGGHK along with its receptor 6CQL (Fig. 5) . Most promising vaccine aspirants must not cause any kind of toxicity or vigorous reaction inside the host. So, checking of toxic nature of epitopes is notably important. This prominently important step was performed by Toxin Pred. It was found that all 4 selected epitopes were non-toxic (Table 4 ). All epitopes along with their result of toxicity analysis and physicochemical properties like hydrophobicity, hydrophilicity, and molecular weight were summarized in Table 4. MHC Pred server was used to study the binding affinity of four non-allergic, non-toxic, and antigenic peptides with allele's viz., HLA DRB1_0101, HLA DRB1_0401, and HLA DRB1_0701. Binding affinity was depicted in terms of IC 50 value (Table 5) . Epitopes showing IC 50 value less than 500 nM were considered to be good binders. Epitopes viz., CEETFGIRL and VRHHGGGHK were found to bind with HLA DRB1_0101 and HLA DRB1_0401, respectively. Both FLKIDPPIL and MKFLFPLKL were found to bind with HLA DRB1_0101 and HLA DRB1_0701. Most eligible vaccine contenders must show satisfactorily population coverage in different parts of the world. Both the (Fig. 6) . Epitope CEETFGIRL and VRHHGGGHK shows population coverage of 11.53% and 11.21% worldwide respectively (Figs. 7 and 8) . The final selection of best and most promiscuous vaccine bidders depends on two main factors, one is low binding energy and another one is high population coverage worldwide. Based on these two factors, epitope FLKIDPPIL was refined. To check the stability of complex FLKIDPPIL-4AH2, molecular dynamics simulation was performed by MD Web simulation. RMSD value of FLKIDPPIL-4AH2 was given in between 0.1 and 1.0 Å (Fig. 9) and B factor scores between 1 and 25 Å 2 (Fig. 10 ). Both RMSD values and B factor plot of complex viz., FLKIDPPIL-4AH2 confirm the stability of the epitope. Reverse vaccinology is known by different names like computational biology, immunoinformatics, and many more. It is a combination of immunological research as well as experimental and computational science. It includes computational tools and software to study the immune response of the host against various infectious diseases. Immunoinformatics helps us to understand antigen presentation in host cells, the behavior of the host during the infection cycle, and thus enriches the knowledge about the disease that affects the immune system and its control (Brusic and Petrovsky 2005) . With the help of Insilico tools, antigenic regions can be mapped easily (Davies and Flower 2007) . Previously, finding these antigenic regions are extremely costly and time-consuming methods like Nuclear Magnetic Resonance (NMR) were used. But today, computational vaccinology had made it possible to predict these antigenic regions in a short period and also with extreme accuracy (Potocnakova et al. 2016) . In this exploration and investigation, the prediction of B cell epitopes has been performed by the authors for the designing of the vaccine against listeriosis by using a reverse vaccinology approach. This research work started with the retrieval of a complete proteome sequence of Listeria monocytogenes F2365, from the UniProtKB database. Most promiscuous B cell epitopes must not show allergic properties. Therefore, to remove all allergic proteins from the investigation AlgPred server was used. A total of 529 proteins of the F2365 strain of Listeria monocytogenes have been proclaimed from the UniProtKB database. Out of 529 proteins, 172 have shown non-allergenicity. These 172 non-allergic proteins have been used to find out the best antigenic regions or peptides that can provoke great immune inflammation in the human body, by using NETMHCII 2.3 server. 299 epitopes have been identified by the VaxiJen server that could bind with MHC II molecules with great stability. Based on the VaxiJen score, only the top 20 B cell epitopes were selected for succeeding refining. 3D modeling of all 20 epitopes has been performed by PEPstrMOD and all these tertiary structures have been saved in PDB format. Tertiary structure modeling of alleles was generated with the help of HLA alleles were performed Based on low binding energy, 4 peptides were selected viz., CEETFGIRL, MKFLFPLKL, FLKIDPPIL, and VRH-HGGGHK. CEETFGIRL showed the energy of − 6.7 kcal/ mol and 6 hydrogen bonds. MKFLFPLKL showed the energy of − 6.9 kcal/mol and 2 hydrogen bonds. FLKIDP-PIL showed the energy of − 7.3 kcal/mol and 6 hydrogen bonds. VRHHGGGHK showed the energy of − 6.7 kcal/ mol and 6 hydrogen bonds. These 4 epitopes were selected on low binding energy as low energy means high stability. Most promiscuous B cell epitope which is a nano peptide, must not be toxic and therefore toxicity analysis must be performed. Toxin Pred server is used for this analysis. This server also anticipates various physicochemical properties of the epitopes like molecular weight, hydrophobicity, and coverage, worldwide epitope FLKIDPPIL was selected. To check the binding energy of epitope FLKIDPPIL with its corresponding 4AH2 receptor molecular dynamics simulation study was performed by using MD Web. RMSD and B factor plot was used to interpret the result of the simulation. After all these vigorous steps of the investigation, epitope FLKIDPPIL proved to be the most eligible candidate that should be used for vaccine designing. Reverse vaccinology has been proved as one of the most powerful weapons to combat some deadly bacterial diseases and had shown tremendous results also. First and foremost, a peptide-based vaccine using the reverse vaccinology approach was created against E. coli in the year 1985 (Jacob et al. 1985) . It has been proved effective against tuberculosis (Mustafa 2013) and many more pathogenic diseases. The identification of antigenic peptides by using a reverse vaccinology approach has been found effective against Staphylococcus aureus (Oyama et al. 2019) . From this research work, we found during the identification and characterization of epitopes for the utility of vaccine designing against Listeria monocytogenes, the epitope FLKIDPPIL was non-allergic, non-toxic, highly antigenic, and can provoke a better immune response. Despite major advancements in the field of technology, society and mankind have been plagued by several kinds of lifethreatening diseases. Although vigorous research is going on, on several deadly diseases in various parts of the world. But still, some foodborne diseases are under-reported and Listeriosis is one of them. In such conditions, computational vaccine technology is one of the best alternatives to deal with such diseases. Computational vaccine technology is a boon in research domains as it accelerates the process of epitope screening and vaccine designing and development. It is a branch of vaccinology that is based on the central idea of solving vaccine development by using a computerdriven algorithm. Listeriosis is still under-reported in many countries of the world. Computational vaccine technology is going to create some awareness and will bring out the best treatment and remedy for the disease. In this research work, after performing molecular docking, 4 epitopes were screened out. These 4 epitopes viz., CEETFGIRL, MKFLF-PLKL, FLKIDPPIL, and VRHHGGGHK were screened as the most promiscuous B cell epitopes among 299 antigenic sites identified. Low binding energy and population coverage analysis predicted FLKIDPPIL as the most potent epitope. Epitope FLKIDPPIL can elicit a strong immune response in the host against listeriosis. Further wet lab trials can assure the stability as well as the response of the epitope in vitro and in vivo. Reverse vaccinology can be proved as the most powerful approach to find remedies against diseases like listeriosis. Validation of predicted virulence factors in Listeria monocytogenes identified using comparative genomics Genome sequence of Listeria monocytogenes Scott A, a clinical isolate from a food-borne listeriosis outbreak Immunoinformatics and its relevance to understanding human immune disease Predicting population coverage of T-cell epitope-based diagnostics and vaccines Clinical features and prognostic factors of listeriosis: the MONALISA national prospective cohort study Harnessing bioinformatics to discover new vaccines VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines MHCPred, a server for quantitative prediction of peptide-MHC binding In silico approach for predicting toxicity of peptides and proteins Listeria monocytogenes: a multifaceted model Listeria monocytogenes serotype 4b strains replicate in monocytes/macrophages more than other serotypes MDWeb and MDMoby: an integrated web-based platform for molecular dynamics simulations Priming immunization against cholera toxin and E. coli heat-labile toxin by a cholera toxin short peptide-beta-galactosidase hybrid synthesized in E. coli Improved methods for predicting peptide binding energy to MHC class II molecules Silico identification of epitope-based peptide vaccine for Nipah virus Computational identification and characterization of potential T-cell epitopes for the utility of vaccine design against Enterotoxigenic Escherichia coli Outcomes of neonates with listeriosis supported with extracorporeal membrane oxygenation from 1991 to 2017 Listeria (listeriosis) | Listeria | cdc In silico analysis and experimental validation of Mycobacterium tuberculosis-specific proteins and peptides of Mycobacterium tuberculosis for immunological diagnosis and vaccine development In silico identification of novel peptides with antibacterial activity against multidrug-resistant Staphylococcus aureus Vaccinomics and a new paradigm for the development of preventive vaccines against viral infections An introduction to B-cell epitope mapping and in silico epitope prediction Bacteriocin from epidemic Listeria strains alters the host intestinal microbiota to favor infection Reassessing the role of Internalin B in Listeria monocytogenes virulence using the epidemic strain F2365 AlgPred: prediction of allergenic proteins and mapping of IgE epitopes Heat inactivation of Listeria monocytogenes and Salmonella enterica serovar Typhi in a typical bologna matrix during an industrial cooking-cooling cycle Kaushik V (2020) In-Silico prediction of a peptide-based vaccine against Zika virus PEPstrMOD: structure prediction of peptides containing natural, non-natural and modified residues Outbreak of Listeriosis in South Africa associated with processed meat AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading Listeria hijacks the clathrin-dependent endocytic machinery to invade mammalian cells SWISS-MODEL: homology modeling of protein structures and complexes An essential role for hfq involved in biofilm formation and virulence in serotype 4b Listeria monocytogenes The authors hereby declare they that have no conflict of interest. The authors did not perform any experiments on human or animals.