The human papillomavirus (HPV) 58 is considered to be the second most predominant genotype in cervical cancer incidents in China. HPV type-restriction, non-targeted delivery, and the highcost of existing vaccines necessitate continuing research on the HPV vaccine. We aimed to explore the papillomaviral proteome in order to identify potential candidates for a chimeric vaccine against cervix papilloma using computational immunology and structural vaccinology approaches. Two overlapped epitope segments (23–36) and (29–42) from the N-terminal region of the HPV58 minor capsid protein L2 are selected as capable of inducing both cellular and humoral immunity. In total, 318 amino acid lengths of the vaccine construct SGD58 contain adjuvants (Flagellin and RS09), two Th epitopes, and linkers. SGD58 is a stable protein that is soluble, antigenic, and non-allergenic. Homology modeling and the structural refinement of the best models of SGD58 and TLR5 found 96.8% and 93.9% favored regions in Rampage, respectively. The docking results demonstrated a HADDOCK score of −62.5 ± 7.6, the binding energy (−30 kcal/mol) and 44 interacting amino acid residues between SGD58-TLR5 complex. The docked complex are stable in 100 ns of simulation. The coding sequences of SGD58 also show elevated gene expression in Escherichia coli with 1.0 codon adaptation index and 59.92% glycine-cysteine content. We conclude that SGD58 may prompt the creation a vaccine against cervix papilloma.

Immunomics tools and databases play an important role in the designing of prophylactic or therapeutic vaccines against pathogenic bacteria and viruses. Therefore, we aimed to illustrate the different immunological databases and web servers used to design a chimeric vaccine candidate against human cervix papilloma. Initially, cellular immunity inducing major histocompatibility complex class I and II epitopes from L2 protein of papilloma 58 strain were predicted using the IEDB, NetMHC, and Tepi tools. Then, the overlapped segments from the above analysis were used to calculate efficiency on interferon-gamma and humoral immunity production. In addition, the allergenicity, antigenicity, cross-reactivity with human proteomes, and epitope conservancy of elite segments were determined. The chimeric vaccine candidate (SGD58) was constructed with two different overlapped peptide segments (23–36) and (29–42), adjuvants (flagellin and RS09), two Th epitopes, and amino acid linkers. The results of homology modeling demonstrated that SGD58 have 88.6% of favored regions based on Ramachandran plot. Protein–protein docking with Swarm Dock reveals SGD58 with receptor complex have −54.74 kcal/mol of binding energy with more than 20 interacting residues. Docked complex are stable in 100ns of molecular dynamic simulation. Further, coding sequences of SGD58 also show elevated gene expression in E. coli. In conclusion, SGD58 may prompt vaccine against cervix papilloma. This study provides insight of vaccine design against different pathogenic microbes as well.

Flagellin activates innate, cellular and humoral immunity through Toll-like receptor 5 (TLR5), and this interaction can be exploited to improve vaccines against oncogenic human papillomavirus (HPV). The binding mechanism between flagellin-fabricated peptide vaccines (FFPVs) and TLR5, however, remains unclear. Here, three FFPV candidates (C58, C52 and C16) were designed from overlapped epitope segments in the conserved N-terminal region of the L2 minor capsid protein of HPV58, HPV52 and HPV16. These segments were predicted to induce cytotoxic T-lymphocyte, T-helper, IFNγ and B-cell responses. Each construct also carried two TLR agonist adjuvants (flagellin and RS09), two universal T-helper epitopes and defined peptide linkers. Physicochemical and immunological properties were evaluated, and three-dimensional models of the FFPVs and TLR5 were generated, refined and validated. Protein–protein docking with HADDOCK and PRODIGY ranked the TLR5-C58 complex highest, ahead of TLR5-C16 and TLR5-C52. To test the docking result, the TLR5-C58 complex was simulated for [100 ns] in explicit solvent, and the root mean square deviation, total energy, hydrogen bonds and key-residue distances were analyzed. The simulation identified a binding mode with ten key residues in the leucine-rich repeat (LRR) region of TLR5, and the binding free energies agreed with experimental reports. These findings offer a structural basis for the rational design of L2-based vaccines against cervical cancer.