key: cord-0992835-cu3dqdw7
authors: Joshi, Amit; Ray, Nillohit Mitra; Singh, Joginder; Upadhyay, Atul Kumar; Kaushik, Vikas
title: T-cell epitope-based vaccine designing against Orthohantavirus: a causative agent of deadly cardio-pulmonary disease
date: 2021-12-07
journal: Netw Model Anal Health Inform Bioinform
DOI: 10.1007/s13721-021-00339-x
sha: d9a9da5604d80074787b9f0ceceea318d65f8805
doc_id: 992835
cord_uid: cu3dqdw7

Orthohantavirus, a zoonotic virus responsible for causing human cardio-pulmonary disease, is proven to be a fatal disease. Due to the paucity of regimens to cure the disease and efficient management to eradicate this deadly virus, there is a constant need to expand in-silico approaches belonging to immunology domain to formulate best feasible peptide-based vaccine against it. In lieu of that, we have predicted and validated an epitope of nine-residue-long sequence “MIGLLSSRI”. The predicted epitope has shown best interactions with HLA alleles of MHC Class II proteins, namely HLA DRB1_0101, DRB1_0401, DRB1_0405, DRB1_0701, DRB1_0901, DRB1_1302, and DRB1_1501. The structure of the epitope was modeled by deploying PEPFOLD 3.5 and verified by Ramachandran plot analysis. Molecular docking and simulation studies reveal that this epitope has satisfactory binding scores, ACE value and global energies for docked complexes along with selectable range of RMSD and RMSF values. Also, the predicted epitope “MIGLLSSRI” exhibits population coverage of more than 62% in world population and maximum of 70% in the United States of America. In this intensive study, we have used many tools like AllergenFP, NETMHCII 3.2, VaxiJen, ToxinPred, PEPFOLD 3.5, DINC, IEDB-Population coverage, MHCPred and JCat server. Most of these tools are based on modern innovative statistical algorithms like HMM, ANN, ML, etc. that help in better predictions of putative candidates for vaccine crafting. This innovative methodology is facile, cost-effective and time-efficient, which could facilitate designing of a vaccine against this virus.

Hantavirus belongs to the family of Bunyaviridae with key features including enveloped single-stranded negative sense RNA; they are mostly associated with zoonosis (Guo et al. 2013 ). Hantavirus pulmonary syndrome and hemorrhagic fever renal syndrome are two common ailments associated with interference of human systems with Hantavirus. The causative agents of two human illnesses are rodent-borne hantaviruses (genus: Orthohantavirus): hemorrhagic fever with kidney disorders (HFRS) in Europe and hantavirus pulmonary syndrome (HPS) in North and South America. HPS mortality rates in South America can range from 35 to 50% (Maes et al. 2009; Ferro et al. 2020 ). An MHC class II molecule binds to epitopes that are produced after proteolytic processes on antigens coming from pathogens within the antigen-presenting cells (APCs), including B cells, macrophages and dendritic cells. Aftermath, MHC-II holding pathogenic epitopes binds to CD4 + T cells to trigger an adaptive immune response against similar kind of pathogens. Hantavax, a vaccine derived from rodent brain formulations with adjuvant like aluminum hydroxide, was used in Korea, but its production was halted and discouraged by many other countries due to use of rodent brain (Schmaljohn 2009 ). The proteasome is largely responsible for the processing of cytosolic proteins. TAP transports the short peptides into the ER, where they are assembled with MHC-I molecules, CD4 + T cell epitopes are moreover less variant in comparison to B cells and CD8 + T cell epitopes, also availability of data supports the multiple MHC Class II HLA variants makes this choice preferable (Blum et al. 2013; Bibi et al. 2021; Sanchez-Trincado et al. 2017) . Exogenous proteins can also be supplied into this pathway via retro-translocation from phagosomes in some antigen-presenting cells, notably dendritic cells, a process known as cross-presentation. The retro-translocation channels may be recruited from the ER, where they are utilized to degrade misfolded membranespanning or cytoplasmic proteins via ER-associated degradation. MHC-II molecules, on the other hand, are largely responsible for presenting exogenous proteins. Antigens are absorbed through phagocytosis, macro-pinocytosis, and endocytosis, and then transported to a mature or late endosomal compartment, in which they are handled and placed into MHC-II molecules. Ribavirin (1-β-d-ribofuranosyl-1,2,4-triazole-3-carboxamide), another potential drug, was also not approved by Food and Drug Administration (FDA) because of poor experimental results (Mahmud-Al-Rafat et al. 2015) . So, there is a great challenge and scope for developing T cell epitopes that can possibly interact directly with MHC molecules and elicit an immunogenic response against Hantavirus. This study is primarily focused in developing epitope-based vaccine against Hantavirus by deploying modern in-silico immunoinformatic approaches. Adaptive immunity plays a crucial role in annihilating viral agents and it is more established with epitope-based peptide vaccines (Zhang 2018) . The spread of Hantaviruses has also been reported to occur from contaminated aerosols inhalation (Jiang et al. 2014) .

In present day times, immunotherapy is an amazing and proficient procedure to control harmful ailments. The prediction of peptide-based epitope vaccine is always a good approach for the treatment of tumors or viral diseases. There are numerous multi-epitope vaccines targeting various diseases including diseases caused by HIV, dengue Sun et al. 2017) , hepatitis B and C, Ebola, Chikungunya, Avian flu-A (H7N9), Zikavirus (Sharma et al. 2020) , classical pig fever, Nipah virus (Kaushik 2020; Liang et al. 2018; Qu et al. 2018) , and Norovirus. A perfect multiepitope vaccine candidate can be used to incorporate a progressive immunological change and it can also evoke either a cellular or a humoral resistant reaction against the infection. A well-designed epitope-based peptide vaccine ought to turn into incredible prophylactic agents against the viral infections. In this investigation, a progression of immunoinformatic approaches has been applied against the proteome of Hantavirus for crafting an epitope-based vaccine. We have vigorously screened the repertoire of epitope by computational methods and validated by molecular simulation of the selected epitope. In addition, the codon usage and in-silico cloning, and the molecular simulations of the best antibody receptor docked complex were successively considered in our study. In present study, we crafted epitope-based vaccine strategy against this harmful pathogen by deploying best modern in-silico tools based on artificial intelligence and neural networks. The aim of study was to predict epitopebased vaccine peptide that holds therapeutic properties against Hantavirus disease. It is a fast and effective way for developing vaccine over hit-and-trial approaches of wet-lab analysis. It provides ease of doing for wet-lab validations and testing for novel vaccination programs.

Primarily, the IEDB (Kim et al. 2012) and ViPR databases were explored to determine essential proteins related to Orthohantavirus. Coding sequence for surface glycoprotein (Accession no. KF537002.1) of Orthohantavirus, and NCBI-genbank accession number (AIA08876.1) for selected protein (SG). Molecular weight of SG Protein was found as 126,191 Da, and Theoretical pI as 7.73 using ExPASY ProtParam tool. The coding DNA sequence and the amino acid sequences were retrieved from NCBI-GenBank database. After that, RSCU values and codon usage analysis of selected coding sequence was done using DnaSP6 software (Rozas et al. 2017) to analyze more validation in selection. For the coding sequence under consideration, ENC-effective number of codons (its value should be in the range of 20-61; 20 indicates high biasness, while 61 shows no biasness) (Wright 1990 ) was computed.

AllergenFP server was deployed to analyze allergenicity of selected protein. NETMHC-II Pan 3.2 server (Jensen et al. 2018 ) was used to screen out single best interacting epitope to maximum number of HLA-DRB alleles, it is based on artificial neural network (ANN) algorithms. Predicted epitope was subjected to VaxiJen server (Doytchinova and Flower 2007) for deriving best antigenicity criterion (> 1.1).

After intensive analysis of epitope sequence, its structure was constructed using PEPFOLD 3.5 server (Lamiable et al. 2016 ). The obtained modeled structure was validated and observed for biochemical and toxicity analysis by deploying MolProbity (Chen et al. 2010) , Ramachandran plot analysis and ToxinPred (Gupta et al. 2013) , QM and SVM scores for toxicity analysis.

RCSB-PDB database (Berman et al. 2000) was deployed to get MHC Class II HLA alleles' structure that exhibits strong interaction to screened and analyzed epitopes. Then molecular docking analysis was performed by DINC server (Antunes et al. 2017 ) and PatchDock-Firedock server (Schneidman-Duhovny et al. 2005 ) to obtain atomic contact energies and binding scores for docked complexes. Molecular simulation studies were performed using Gromacs software (Van Der Spoel et al. 2005 ) used for coarse-grained Brownian dynamics in GROMOS 9643a1 force field environment. GROMOS 9643a1 force field was utilized to evaluate thermodynamic properties throughout the MD simulation procedure. To eliminate very large energy configuration (steric conflicts or bonded components far off from the equilibrium), energy reduction was used. Some big protein structures can result in poor connections and unphysical fluctuations in total energy in later simulations; therefore, a quick reduction can eliminate these issues. After that, the subsystems were subjected through a 100-picosecond NVT ensemble MD_Simulation at 300 K temperature. Consequently, at 300 K and 1 bar pressure, NPT equilibrium was reached. In the binding procedure, a 10 ns MD_Simulation was performed on the selected conformational geometry of the HLA receptor to epitope vaccine complex. Trajectory analysis for 10 ns was used to evaluate RMSD values and RMSF analysis of docked complex (Joshi et al. 2020b ).

Finally, IEDB population coverage analysis (> 50%) based on restriction database (Bui et al. 2006 ) and MHCPred 2.0 server (Guan et al. 2003 ) logIC 50 (1-10) quantitative prediction analysis was deployed to get more validation of selected putative T cell epitope for vaccine crafting. JCat (Java Codon Adaptation Tool) (Grote et al. 2005 ) server was deployed to obtain optimized codon sequence for epitope in E. coli k12 strain to get CAI near to best score 1.0 and GCcontent codon percentage in 30-70% range to ensured high expression rate of vaccine candidate in E. coli k12 strain. 

NETMHC-II Pan 3.2 server (Jensen et al. 2018 ) was used to screen out single best interacting Epitope to maximum number of HLA-DRB alleles, it is based on artificial neural network (ANN) algorithms. Predicted epitope subjected to VaxiJen server (Doytchinova and Flower 2007) for deriving best antigenicity criterion (> 1.1). NETMHC-II PAN 3.2 server and VaxiJen results in Table 1 depict the 1-log50k value, affinity, and antigenicity criterion of selection for epitope and HLA alleles. Here, we found MIGLLSSRI as a putative epitope that exhibits interactions with all selected HLA-DRB alleles of MHC-II. As the NETMHC server for screening epitopes is based on ANN, makes it efficient in prediction of epitopes with considered HLA alleles. Figure 2 is a graphical representation of 1-log50k (affinity) values for better understanding of epitope selection. The structure for epitope MIGLLSSRI is obtained from the fast structure prediction De Novo method by deploying PEP-FOLD3.5, it was found to be helical in shape. Figure 3 shows its predicted structure viewed in PyMOL and probability plot obtained from PEP-FOLD server.

Toxicity analysis and Ramachandran plot analysis provide essential information about peptide toxicity and secondary structure. ToxinPred server designing peptide module allows us to generate all possible scores and help in predicting whether the considered peptide is toxic or not, it is based on QM (− 33.50) and SVM (− 0.56) scores (Table 2) . Molprobity was deployed to generate Ramachandran plot for selected epitope, and it was found to have good acceptance (more than 85% in allowed region) of structure (Fig. 4) . Also, half-life prediction was estimated to 30 h in mammalian reticulocytes by deploying ExPASy tool (ProtParam) (Gasteiger et al. 2005 ).

Structures of HLA-DRB alleles were retrieved from RCSB-PDB database enlisted in Table 3 . The selected HLA alleles are DRB*01:01, DRB*04:01, DRB*04:05, DRB*07:01, DRB*09:01, DRB*13:01, DRB*15:01.

Molecular docking studies assist in examining the best binding scores, atomic contact energies and global energies for docked complexes (Table 4 ). For this, DINC server was used as it is based on Autodock software and also Patch-Dock along with Fire-Dock screening to generate more effective aspects of interaction. A comparison of energy values against docking complexes is graphically represented in Fig. 5 .

The best docked interaction complex 1FV1-MIGLLSSRI is represented in Fig. 6 , as it shows a binding score of − 7.4 kcal/mol, ACE value − 10.87, and global energy − 52.59. This epitope MIGLLSSRI exhibits better interactions with other HLA alleles too. This makes the epitope the most suitable candidate for crafting peptide-based vaccine. 

Molecular simulation investigations assisted to reveal the interaction of MIGLLSSRI epitope on the basis of RMSD value per residue. RMSD values (0-1 nm) and RMSF values (1-1.5 nm) for docked complexes are represented by results of Gromacs tool in Fig. 7 . Trajectory analysis for docked complex dynamic simulation was conducted for 10 ns. The RMSD and RMSF plots clearly represent less fluctuation under the simulation time span that makes the interaction more stable between HLA allele and selected epitope MIGLLSSRI.

IEDB Population conservancy analysis outcomes were 62.96%, 70.16%, 66.48%, 59.61% population coverage for World, USA, Europe and India, respectively. These are satisfactory results for considering MIGLLSSRI epitope for crafting vaccine (Fig. 8) . MHCPred server has three HLA-DRB1 alleles (01*01, 04*01, 07*01) and results on LogIC 50 value (Table 5) were also found satisfactory and provide more confidence in epitope selection.

JCat server results show GC content of about 48.148% and CAI value 1.0 with improved or optimized codon sequence in E. coli K12, i.e., "ATG ATC GGT CTG CTG TCT TCT CGT ATC" for epitope MIGLLSSRI. This clearly indicates a high expression of vaccine candidate in E. coli K12 strain (Ali et al. 2017) . So, this epitope can be easily used for vaccine formulations against Hantavirus after wet-lab validation.

Orthohantavirus group of virus is a major etiological threat as associated with Hantavirus cardio-pulmonary syndrome in North American (Richardson et al. 2013) , Korean and Chinese populations. Immunologically HLA haplotypes of hosts are major determinants to develop B-cell and T-cell immune response (Hooper et al. 2013 ).

The current study primarily focuses on screening peptide from glycoprotein present in Orthohantavirus, which can exhibit interaction with MHC Class II HLA alleles, and elicit adaptive immune response in the human host domain. For this purpose, NETMHC-II pan 3.2 server, molecular docking, molecular simulations, and population coverage were applied. The surface or envelope glycoproteins create heterodimers, which then form tetrameric spikes, which are the lattice's building components. Viral entry is achieved by receptor-mediated endocytosis and endosomal membrane fusion, and glycoproteins are the exclusive targets of neutralizing antibodies (Serris et al. 2020), so we targeted SG protein to develop (Sun et al. 2017; Liang et al. 2018; Qu et al. 2018) . To determine the above issues, we propose to build an all-inclusive hereditary designed novel subunit protein antibody against Hantavirus by consolidating bioinformatics strategies, i.e., viral surface protein structure science information. The most elevated antigenic proteins were exposed to the NetMHC server and an aggregate of 3769 epitopes with a joined score (1-log50k ≥ 0.6) were anticipated. The anticipated epitopes were considered to explore their antigenicity, immunogenicity, conservancy, harmfulness, and their separate MHC HLA alleles. From that point, the anticipated epitopes were subjected to docking-simulation parameters to screen whether the epitopes effectively fulfill all the rules. Afterward, these epitopes were tried for biochemical parameters testing and population conservancy analysis to determine the susceptible nature. Lastly, the best epitope (i.e., MIGLLSSRI) from considered surface glycoprotein was chosen as an immunization applicant that has attributes superior to other epitopes of that protein. In recent studies, it was found that SG protein was found to be main factor for immunogenic peptide development by deploying proteomic analysis study (Abdulla et al. 2021) . Similar methodology was applied for successful prediction of SARS-CoV-2 epitopes in recent studies (Akhtar et al. 2020; Joshi et al. 2020a ). Many other recent studies successfully demonstrated that epitope-based vaccine designing was successful against Dengue virus (Krishnan et al. 2020 , Candida fungus (Akhtar et al. 2021a ), Human cytomegalovirus (Akhtar et al. 2021b) , Canine circovirus (Jain et al. 2021 ) and bacterium Tropheryma whipplei . This study will open new prospects in developing vaccine agents against Hantavirus as it is economically very beneficial and is a viable alternative to cumbersome, long hit-and-trial wet-lab experiments.

Immunoinformatic methodologies are very useful in the modern era, and can be easily implemented, but extensive studies in revealing further new epitopes and wet-lab screening is a must for any immunotherapy. Our predicted epitope "MIGLLSSRI" passes all in-silico parameters and makes a way directing modern vaccine designers to go for specificity rather than getting involved in long hit-andtrial experiments. The molecular dynamics and simulation studies along with population coverage and expressive behavior of this epitope mark it as a possible choice in crafting a vaccine against Hantavirus. 

A comprehensive screening of the whole proteome of hantavirus and designing a multi-epitope subunit vaccine for cross-protection against hantavirus: structural vaccinology and immunoinformatics study

Immuno-Informatics quest against COVID-19/SARS-COV-2: determining putative T-cell epitopes for vaccine prediction

Insilico design of a multivalent epitope-based vaccine against Candida auris

Design of a novel and potent multivalent epitope based human cytomegalovirus peptide vaccine: an immunoinformatics approach

Exploring dengue genome to construct a multi-epitopebased subunit vaccine by utilizing immunoinformatics approach to battle against dengue infection

DINC 2.0: a new protein-peptide docking webserver using an incremental approach

RCSB Protein Data Bank: structural biology views for basic and applied research

In silico analysis of epitope-based vaccine candidate against tuberculosis using reverse vaccinology

Pathways of antigen processing

Predicting population coverage of T-cell epitope-based diagnostics and vaccines

MolProbity: all-atom structure validation for macromolecular crystallography

VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines

Hantavirus pulmonary syndrome outbreaks associated with climate variability in Northwestern Argentina

Protein Identification and Analysis Tools on the ExPASy Server

JCat: a novel tool to adapt codon usage of a target gene to its potential expression host

MHCPred: a server for quantitative prediction of peptide-MHC binding

Phylogeny and origins of hantaviruses harbored by bats, insectivores, and rodents

In silico approach for predicting toxicity of peptides and proteins

A novel Sin Nombre virus DNA vaccine and its inclusion in a candidate pan-hantavirus vaccine against hantavirus pulmonary syndrome (HPS) and hemorrhagic fever with renal syndrome (HFRS)

An immunoinformatics study: designing multivalent T-cell epitope vaccine against canine circovirus

Improved methods for predicting peptide binding affinity to MHC class II molecules

Development of a SYBR Green I based one-step real-time PCR assay for the detection of Hantaan virus

In-silico proteomic exploratory quest: crafting T-cell epitope vaccine against Whipple's disease

Epitope based vaccine prediction for SARS-COV-2 by deploying immunoinformatics approach

Molecular docking and simulation investigation: effect of beta-sesquiphellandrene with ionic integration on SARS-CoV2 and SFTS viruses

In silico identification of epitope-based peptide vaccine for Nipah virus

Immune epitope database analysis resource

2020) T cell epitope designing for dengue peptide vaccine using docking and molecular simulation studies

Immunoinformatics designed T cell multi epitope dengue peptide vaccine derived from non structural proteome

PEP-FOLD3: faster de novo structure prediction for linear peptides in solution and in complex

Defining of MAbs-neutralizing sites on the surface glycoproteins Gn and Gc of a hantavirus using vesicular stomatitis virus pseudotypes and site-directed mutagenesis

Recombinant Zika virus envelope protein elicited protective immunity against Zika virus in immunocompetent mice

Recent approaches in hantavirus vaccine development

Effects of anthropogenic events and viral persistence on rodent reservoirs of hantavirus infection understanding host-pathogen interactions facilitates novel approaches to intervention strategies

Insect cell-produced recombinant protein subunit vaccines protect against Zika virus infection

Human exposure to particulate matter potentially contaminated with Sin Nombre virus

DNA sequence polymorphism analysis of large data sets

Fundamentals and Methods for T-and B-Cell Epitope Prediction

Vaccines for hantaviruses

PatchDock and SymmDock: servers for rigid and symmetric docking

The hantavirus surface glycoprotein lattice and its fusion control mechanism

Kaushik V (2020) In-silico prediction of peptide based vaccine against Zika virus

Elaboration of tetravalent antibody responses against dengue viruses using a subunit vaccine comprised of a single consensus dengue envelope sequence

GROMACS: fast, flexible, and free

The 'effective number of codons' used in a gene

Multi-epitope vaccines: a promising strategy against tumors and viral infections

The authors acknowledge computational services