Abstract
Since 1997, highly pathogenic avian influenza viruses, such as H5N1, have been recognized as a possible pandemic hazard to men and the poultry business. The rapid rate of mutation of H5N1 viruses makes the whole process of designing vaccines extremely challenging. Here, we used an in silico approach to design a multi-epitope vaccine against H5N1 influenza A virus using hemagglutinin (HA) and neuraminidase (NA) antigens. B-cell epitopes, Cytotoxic T lymphocyte (CTL) and Helper T lymphocyte (HTL) were predicted via IEDB, NetMHC-4 and NetMHCII-2.3 respectively. Two adjuvants consisting of Human β-defensin-3 (HβD-3) along with pan HLA DR-binding epitope (PADRE) have been chosen to induce more immune response. Linkers including KK, AAY, HEYGAEALERAG, GPGPGPG and double EAAAK were utilized to link epitopes and adjuvants. This construct encodes a protein having 350 amino acids and 38.46 kDa molecular weight. Antigenicity of ~ 1, the allergenicity of non-allergen, toxicity of negative and solubility of appropriate were confirmed through Vaxigen, AllerTOP, ToxDL and DeepSoluE, respectively. The 3D structure of H5N1 was refined and validated with a Z-Score of - 0.87 and an overall Ramachandran of 99.7%. Docking analysis showed H5N1 could interact with TLR7 (docking score of - 374.08 and by 4 hydrogen bonds) and TLR8 (docking score of - 414.39 and by 3 hydrogen bonds). Molecular dynamics simulations results showed RMSD and RMSF of 0.25 nm and 0.2 for H5N1-TLR7 as well as RMSD and RMSF of 0.45 nm and 0.4 for H5N1-TLR8 complexes, respectively. Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) confirmed stability and continuity of interaction between H5N1-TLR7 with the total binding energy of - 29.97 kJ/mol and H5N1-TLR8 with the total binding energy of - 23.9 kJ/mol. Investigating immune response simulation predicted evidence of the ability to stimulate T and B cells of the immunity system that shows the merits of this H5N1 vaccine proposed candidate for clinical trials.