key: cord-0969754-j0runrkf
authors: Kaliamurthi, Satyavani; Selvaraj, Gurudeeban; Kaushik, Aman Chandra; Gu, Ke-Ren; Wei, Dong-Qing
title: Designing of CD8(+) and CD8(+)-overlapped CD4(+) epitope vaccine by targeting late and early proteins of human papillomavirus
date: 2018-10-02
journal: Biologics
DOI: 10.2147/btt.s177901
sha: 4aa64b68fec8e053f6abcc1594eee14507853618
doc_id: 969754
cord_uid: j0runrkf

BACKGROUND AND AIM: Human papillomavirus (HPV) is an oncogenic agent that causes over 90% of cases of cervical cancer in the world. Currently available prophylactic vaccines are type specific and have less therapeutic efficiency. Therefore, we aimed to predict the potential species-specific and therapeutic epitopes from the protein sequences of HPV45 by using different immunoinformatics tools. METHODS: Initially, we determined the antigenic potential of late (L1 and L2) and early (E1, E2, E4, E5, E6, and E7) proteins. Then, major histocompatibility complex class I-restricted CD8(+) T-cell epitopes were selected based on their immunogenicity. In addition, epitope conservancy, population coverage (PC), and target receptor-binding affinity of the immunogenic epitopes were determined. Moreover, we predicted the possible CD8(+), nested interferon gamma (IFN-γ)-producing CD4(+), and linear B-cell epitopes. Further, antigenicity, allergenicity, immunogenicity, and system biology-based virtual pathway associated with cervical cancer were predicted to confirm the therapeutic efficiency of overlapped epitopes. RESULTS: Twenty-seven immunogenic epitopes were found to exhibit cross-protection (≥55%) against the 15 high-risk HPV strains (16, 18, 31, 33, 35, 39, 51, 52, 56, 58, 59, 68, 69, 73, and 82). The highest PC was observed in Europe (96.30%), North America (93.98%), West Indies (90.34%), North Africa (90.14%), and East Asia (89.47%). Binding affinities of 79 docked complexes observed as global energy ranged from −10.80 to −86.71 kcal/mol. In addition, CD8(+) epitope-overlapped segments in CD4(+) and B-cell epitopes demonstrated that immunogenicity and IFN-γ-producing efficiency ranged from 0.0483 to 0.5941 and 0.046 to 18, respectively. Further, time core simulation revealed the overlapped epitopes involved in pRb, p53, COX-2, NF-X1, and HPV45 infection signaling pathways. CONCLUSION: Even though the results of this study need to be confirmed by further experimental peptide sensitization studies, the findings on immunogenic and IFN-γ-producing CD8(+) and overlapped epitopes provide new insights into HPV vaccine development.

Human papillomavirus (HPV) is an oncogenic, mucosa-infecting, and etiologic agent that causes over 90% of cases of cervical cancer in the world. 1, 2 The genome of HPV consists of three regions, viz, noncoding region (extended control), and nonstructural (E1-E7) and structural regions (L1, L2) of the open reading frames. 3, 4 Briefly, E1 protein is the ATP-dependent DNA helicase which is involved in replication. E2 protein acts as a regulator of viral gene transcription and genome replication. E4 and E5 proteins are involved in the early and late phases of the virus life cycle and immune modulation.

E6 and E7 proteins are involved in cell differentiation and targeting numerous negative regulators of the cell cycle (tumor suppressor genes p105Rb and p53). The assemblage of L1 and L2 proteins in capsomers results in the formation of icosahedral capsids around the viral genetic material during formation of progeny. 5, 6 Identification of CD8 + antigenic determinants is crucial in vaccine development. 7 Cytotoxic T lymphocytes (CTLs) are also known as killer cells or CD8 cells matured from the thymus of a specific primary lymphoid organ of the human immune system. CD8 + T cells play a vital role in cell-mediated immune responses against invading intracellular foreign substances like viruses. 8 In addition to CD8 + T cells, the peptides of CD4 + T cells are presented on major histocompatibility complex (MHC) class II, which are located in the antigen-presenting cells (APCs), to trigger immune response. Therefore, CD4 + cells are also known as T helper (Th) cells. Initially, APC identifies the pathogen-associated molecular pattern by toll-like receptors (TLRs) which results in the cleavage of antigen into small fragments. Then, MHC class II peptides complex activates the naive CD4 + Th cells and triggers CD8 + -mediated cellular immune response or B cell-mediated humoral immune response. Interestingly, when the antigen complexes with MHC class I directly, it results in CD8 + -mediated cellular immune response. 9, 10 The detection of CTL epitopes and the stimulation of immune response play a critical role in an individual's immune system.

Currently available prophylactic vaccines such as Gardasil (HPV6, 11, 16, and 18) and Gardasil 9 (HPV31, 33, 45, 52, and 58) could provide protection to women by inducing significant immunity, but these vaccines are limited to those aged between 9 and 26 years. 11 However, these prophylactic vaccines demonstrated no therapeutic benefits in patients who already had HPV infections. 12, 13 The US Food and Drug Administration (FDA) and the Centers for Disease Control and Prevention (CDCP) have recommended that HPV vaccines are safe but can cause more pain, uneasiness, soreness, and Guillain-Barre syndrome in patients. 14, 15 However, to control the prevalence rate of high-grade HPV infections, species-specific screening and vaccination will be needed. 16 Out of 150 types of HPV, the 13 high-risk HPV (hrHPV) types (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68) which induce uncontrolled cell growth in the cervix uteri region have been identified. 17, 18 Among these, HPV16 and 18 were reported as the major strains to cause cervical cancers (about 67.4%) around the world. [19] [20] [21] In fact, the contributions of other hrHPV types are limited due to their geographical distribution. HPV45 is a member of HPV18 alpha papillomaviridae family and a species of α7. [22] [23] [24] Moreover, the prevalence rate of HPV45 in cases of cervical cancer was found to be high varying from 1.6% to 37.4% which was recorded in different geographical regions such as Ghana, Portugal, Serbia, Spain, Pakistan, USA, Japan, Burkina Faso, Northern parts of China, Southeast regions of Brazil, and Northeast of India. 18, [25] [26] [27] [28] [29] [30] [31] [32] Difference in the host genetic framework and oncogenicity of the circulating variants might play a critical role in geographical distribution of HPV45. Together with HPV16 and 18, HPV45 is highly recorded in cases of squamous cell carcinoma (75%) and adenocarcinoma (94%), and in cases of invasive cervical cancer and cervical intraepithelial neoplasia grade 3. 2, 16, 33 Based on the level of enrichment in cervical cancer cases compared to women with normal cytological features, HPV45 has been suggested as the third highest oncogenic agent next to HPV16 and 18. 1, 30, 34, 35 In comparison to other hrHPV types, limited reports exist about the peptide research in HPV45.

Immunogenic peptides are one of the excellent resources to develop therapeutic vaccines against hrHPV strains. Instant accessibility of open databases and informatic tools suggest efficient methods to develop vaccines with reduced time, labor, and experimental cost. 36, 37 Several studies reported that virtual peptides of HPV16 and 18 oncoproteins (E6 and E7) efficiently induce immune response in cervical cancer patients as revealed by their peripheral blood mononuclear cell samples. [38] [39] [40] [41] [42] However, cross-protection against all the hrHPV strains could not be ignored. 40 The cross-reactivity of CTL peptide is suspected to play an essential role in producing immune response against multiple strains. 43 Selection of immunogenic epitopes from the whole genome of hrHPV strains may enable designing of novel CTL vaccines. Therefore, we identified the potential immunogenic CD8 + , CD8 + -overlapped CD4 + , and B-cell epitopes from the antigenic proteins of HPV45 and anticipated an epitopes pool with the capability to demonstrate immunogenic responses against all 15 hrHPV strains. The workflow of the study is illustrated in Figure 1 .

The whole amino acid sequence of HPV45 genome was retrieved from UniProt Knowledgebase (UniProtKB) database. The amino acid sequences of complete list of nonstructural and structural proteins of HPV45, which included sequences of replication initiation protein E1, regulatory proteins E2, E4, E5, E6, and E7, major capsid protein L1, and minor capsid protein L2, were obtained in October 2017. 44 Biologics: Targets and Therapy 2018:12 submit your manuscript | www.dovepress.com

CD8 + and its overlapped CD4 + epitopes of HPV In total, sequences of 139 proteins were extracted which included E1 (11), E2 (12), E4 (04), E5 (03), E6 (27) , E7 (26) , L1 (32) , and L2 (24) . Alignment-independent VaxiJen 2.0 tool was employed to determine the antigenicity of the complete list of selected sequences. The protective antigenic score was calculated based on the auto cross-covariance (ACC), which is the transformation of protein sequences into normal vectors with essential amino acid properties. 45 (1)

where A jj (l) denotes ACC, n is the number of amino acids in a sequence, i is the amino acid position, and l is the lag; index j (1-3) was used for the Z scales. Antigenic score of each protein was predicted with the respective threshold value (>0.4). In addition, the same method was followed for the overlapped epitopes.

Prediction of cD8 + cTl epitopes NetCTL 1.2 server was employed to predict CTL epitopes from the input sequence. The integrated prediction (PI) scores of the CTL epitopes with the presence of MHC ligand were obtained as a result. The threshold value 0.75 with a sensitivity of 0.80 and a specificity of 0.970 was set for CTL epitope prediction. 46 The Immune Epitope Database (IEDB) combined tool was used to predict the processing of antigen by MHC class I antigen processing and presentation pathway. NetMHC Artificial Neural Network (ANN) 3.4 was used to find out the natural potential (proteasomal cleavage, transporter associated with antigen processing (TAP) transport, and MHC class I binding) of being a CTL epitope. 47 Fifty-one frequently occurring MHC class I alleles, including human leukocyte antigen (HLA)-A (16), HLA-B (23), HLA-C (10), and HLA-E (2), were selected as input. In the TAP transport predictions, the peptide length (9mer) and 0.2 alpha factors were defined as threshold. The MHC class I alleles showing a binding affinity of IC 50 <200 nm were selected and used for further study.

selection and evaluation of cD8 + T-cell epitopes

The immunogenicity of the optimized MHC class I peptides was determined using T-cell class I pMHC immunogenicity prediction tool. Based on the properties and position of amino acids within the peptide, the immunogenicity of an MHC peptide was predicted using the following formula. 48

where S denotes the immunogenicity score of a peptide ligand (L) presented on an HLA molecule (H). For every position (p) in the ligand (L), the log enrichment score (E) for the amino acid at that position A (L, p) weighted by the importance of that position I p was summed. M is the masking of anchor positions on the HLA. In the present analysis, the first, second, and C-terminal amino acids were masked. The same method was followed for the overlapped epitopes.

Epitopes conservancy analysis tools were used to find out the degree of conservancy of an epitope within a set of given protein sequences at a given identity level. 49 A set of sequences of epitopes (e) and proteins {p} were given as input. The degree of conservation percentage of e within p was estimated as the fraction of {p} that matched the aligned e above a chosen identity level. The conservancy analysis was performed for structural and nonstructural protein sequences of 15 hrHPV strains (HPV16, 18, 31, 33, 35, 39, 51, 52, 56, 58, 59, 68, 69, 73, and 82) .

The IEDB PC tool was used to calculate the PC of predicted CTL epitopes. It calculates the fraction of individuals predicted to respond to an epitope set based on the HLA genotype frequencies and MHC binding and T-cell restriction data. 50 It provided allele frequencies for 21 different ethnicities in 115 countries that assembled into 16 different geographical areas. The number of epitope sets, epitopes with MHC-restricted alleles, MHC class selection, and the population's area were selected as per the user's choice during the analysis.

Prediction of peptide structure and homology modeling PEP-FOLD3 was used to predict the chemical structure of the selected antigenic proteins. 51 Hidden Markov model suboptimal conformation sampling approach was applied to predict the structure from the input peptide sequence. For confirmation, 200 simulation runs were selected and the different cluster models were sorted by the sOPEP energy. 52 Phyre2 and Iterative Threading ASSEmbly Refinement (I-TASSER) were used to predict the three-dimensional (3D) structure of target receptors. The Phyre2 is a protein modeling and analysis engine which uses highly innovative remote homology detection methods to accurately construct the 3D structure of targets within a duration of 30 minutes to 2 hours. 53 I-TASSER is another vital protein structure prediction tool that identifies the respective structural templates from the Protein Data Bank (PDB) through local meta-threading server (LOMETS) with full-length atomic models created by the repetitive template fragment 

PatchDock tool was employed for shape representation, surface patch matching, and filtering and scoring. 55 Based on the ranking, the elite peptide models were selected for docking with targets. Further, the prepared HLA molecules were used for docking with the predicted epitopes. Then, the results were refined by FireDock server according to the energy function. In the final refinement, a full interface side chain was optimized with an atomic radius of 0.85 to facilitate decreasing the number of conflicts. 56 The binding interactions of the docked structures were evaluated using UCSF Chimera 1.11.2, a highly extensible program for interactive visualization (data not shown). 57 Prediction and identification of CD8 + epitopes overlapping cD4 + and B-cell epitopes Prediction of interferon-gamma (iFn-γ)-producing cD4 + T-cell epitopes IEDB consensus method was used to predict CD4 + T-cell epitopes by combining the following: neural network alignment (Net MHCII 2.2), stabilization matrix method alignment, and CombLib or Sturniolo method. IEDB is a comprehensive dataset consisting of over 10,000 unpublished MHC peptide binding affinities, 29 crystal structures, and 664 experimentally proven CD4 + peptides. The major population in the tool is HLA-DRB locus, which consists of 51 alleles. HLA-DRB1, HLA-DRB3, HLA-DRB4, and HLA-DRB5 subtypes were selected for the analysis. As a result, the lowest percentile rank peptides exhibited higher binding affinity. 58 The antigen was presented by the APCs to the CD4 + receptor of Th cells, and then Th cells producing IFN-γ were activated. Moreover, the IFN-γ-producing CD4 + MHC class II-restricted epitopes were analyzed using IFN epitope tool. The prediction program was used with the motif-based model or support vector machine or hybrid algorithms. It showed an accuracy of 81.39% for the prediction of IFN-γ-producing epitopes. 59 The dataset consists of 3,705 IFN-γ-inducing and 6,728 non-IFN-γ-inducing MHC class II binders.

The B-cell epitopes of the selected antigenic sequences were predicted using ABCpred tool. The results were predicted based on the experimentally trained B-cell epitope datasets. ABCpred is the first server created based on recurrent neural network (ie, Jordan network, a machine-based technique) using fixed-length patterns with an accuracy of 65.93%. 60 evaluation of selected overlapping epitopes

AllergenFP v.1.0 is an essential tool to predict allergenicity of peptide sequences based on novel descriptor fingerprint approach with an accuracy of 87%, a specificity of 89%, and a sensitivity of 86%. It contains 2,427 allergenic and 2,427 non-allergenic amino acid segments. 61 analysis of cross-reactivity with human proteomes

The similarity analysis was performed to find out the epitopes that match with the human proteome by using the integrated open database, Protein Information Resource, with peptidematching program. 62, 63 systems biology-based virtual cancer pathway prediction and kinetics A literature survey was performed to collect the required information about the selected protein sequences to design and predict the peptides involved in biochemical pathways. [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] A mathematical and computational model of the entire HPV45 infection signaling pathway interacting with specific genes in the presence of peptides was developed and was visualized using CellDesigner v4.4 modeling tool. 75 Nodes in the pathway represented the entities, and edges represented connectivity of one node to another, which was closely related to each other. 76 The inhibitory concentration (1 µg) of proposed peptides was collected from the published experimental data and pharmacokinetic studies. 64 Mass action kinetics (V=KΠi S) was applied for pharmacokinetic study. Time course and kinetics simulation was performed to concern transition time of cancer and complete Table S1 .

cD8 + cTl epitopes In total, 80 CTL epitopes were predicted from the early proteins E1 (19) , E2 (17), E4 (01), E6 (05), and E7 (14), and late proteins L1 (07) and L2 (17) of HPV45 based on MHC class I supertype by NetCTL 1.2 (threshold value >0.75). Among the 80 epitopes, only 50 CTL epitopes showed a higher binding affinity (IC 50 <200 nm) for MHC class I alleles as found using NetMHC ANN 3.4. Table 1 depicts the proteasome score, TAP score, MHC score, processing score, total score, and binding affinity of predicted epitopes for different alleles as determined using NetMHC ANN 3.4. According to the results, 17 out of 19 epitopes of E1 protein, such as ATDTGSDMV 40 

In the present findings, 27 out of 50 epitopes were found to be immunogenic. A list of immunogenic epitopes from structural and nonstructural proteins is illustrated in Table S2 .

While designing epitope-based vaccines, the use of conserved epitopes may be helpful to provide broader protection across multiple strains and to develop cross-reactive vaccines. Table  S3 shows the detailed results of conservancy (≥50%) analysis of the 27 immunogenic epitopes and their positions, and protein sub-sequences and their respective cross-protection against hrHPV strains. The conservation of immunogenic epitopes across the 15 hrHPV strains is shown in Figure 2 . 

CD8 + and its overlapped CD4 + epitopes of HPV 

CD8 + and its overlapped CD4 + epitopes of HPV Notes: a cTl prediction by netcTl 1.2 (threshold >0.75). b Prediction by netMhc ann 3.4; proteasome cleavage score is interpreted as logarithm of the total amount of cleavage site usage liberating the peptide c-terminus; TaP score estimates an effective log(ic 50 ) value for binding to TaP of a peptide or its n-terminal prolonged precursors; Mhc binding prediction score is identical to the class i -log(ic 50 ) values; processing score combines the proteasomal cleavage and TaP transport predictions, and predicts a quantity proportional to the amount of peptide present in the endoplasmic reticulum, where a peptide can bind to multiple Mhc molecules; total score combines the proteasomal cleavage, TaP transport, and Mhc-binding predictions and predicts a quantity proportional to the amount of peptide presented by Mhc molecules on the cell surface. Abbreviations: cTl, cytotoxic T lymphocyte; hla, human leukocyte antigen; hPV, human papillomavirus; Mhc, major histocompatibility complex; TaP, transporter associated with antigen processing.

Biologics: Targets and Therapy 2018:12 submit your manuscript | www.dovepress.com

Kaliamurthi et al America (80.95%). About 93.68% of cumulative PC was achieved around the world (Figure 3 ). 

CD8 + and its overlapped CD4 + epitopes of HPV Mhc class ii-restricted cD4 + epitopes and iFn-γ production

The high-ranked MHC class II CD4 + epitopes with lower percentile value <1% and their IFN-γ-producing capacities are illustrated in Table S4 . Among them, the high-ranked CD4 + epitopes and the overlapped potential immunogenic CD8 + T-cell epitope regions are highlighted in bold.

In total, 197 B-cell epitopes (16mer amino acid length) were predicted from N-and C-terminal regions of the E1 (61) Notes: a FireDock is a rigid body scoring energy for best-ranked solution number, global energy, attractive and repulsive van der Waals energy, atomic contact energy, and hydrogen bond energy of the peptide and Mhc receptor protein complex. b amino acid residues of peptide and Mhc receptor protein docked complex and the initial distance between hydrogen bond donor and acceptor analyzed by chimera. Abbreviations: ac, atomic contact; aVDW, attractive van der Waals; hB, hydrogen bond; hla, human leukocyte antigen; Mhc, major histocompatibility complex; rVDW, repulsive van der Waals.

Biologics: Targets and Therapy 2018:12 submit your manuscript | www.dovepress.com

CD8 + and its overlapped CD4 + epitopes of HPV V ( F T F P H A F P F ) D K N G N P 5 7 8 -5 9 3 ( 0 . 7 6 ) a n d Overlapping epitopes Table 3 illustrates the epitope sequence, position in the selected protein, allergenicity, antigenicity, immunogenicity, and IFNγ-producing efficiency of 23 CD8 + T-cell epitopes overlapping CD4 + and B-cell epitopes. All the selected epitopes were found to be immunogenic. Among the 23 epitopes, 19 were identified as non-allergenic and 15 were antigenic with a threshold value of 0.4. In addition, 19 epitopes had the efficiency to induce IFN-γ production. The similarity level zero indicated no unique matches of query epitopes found against the Homo sapiens (Taxonomy ID: 9606). Moreover, these peptides would eliminate the risk of inducing autoimmune disorders.

The predicted peptides of E1, E2, E6, and L1 inhibit or trigger the P13K/AKT signaling pathway, MAPK pathway, overexpression of COX-2, downregulation of tumor suppressor p53, transcriptional repressor NF-X1 gene, multifunctional regulator gene (MYC), tumor suppressor pRb, higher expression of minichromosome maintenance proteins, proliferating cell nuclear antigen, P16 protein, antigen Ki-67, alternate reading frame protein (P14ARF), and MDM gene which are involved in uncontrolled cell division, virus immortalization, MAPK pathway, apoptosis, notch pathway, transcription, and telomerase reverse transcription ( Figure 4A ). These genes and pathways are actively involved in the development of HPV causing cervical cancers including cervical intraepithelial neoplasia and invasive cervical carcinoma. The results Figure  4B , which indicate the concentration of entities and time in different colors. After 30 seconds of time course simula-tion, all entities were constant. The results indicated that the proposed peptides of HPV45 could work as best anticancer peptides through the inhibition of various pathways involved in the cervical cancer development. The kinetics simulation was also performed using different concentrations of HPV45 

CD8 + and its overlapped CD4 + epitopes of HPV peptides for cross-validation. The results obtained at the concentration of 0.0300 µm support the literature information. Time course and kinetics simulation prediction demonstrated the biological behavior of HPV45 in the presence of peptides where HPV45 interacts with different molecules, which play an important role in cervical cancer inhibition.

HPV infection is solely responsible for 99.7% of cervical cancer cases. Fitzmaurice et al 78 reported the agestandardized incidence and death rates of cervical cancer in the developing and developed countries from 1990 to 2013. Exploitation of bioinformatic techniques has led to an emerging concept in advanced vaccine design against different diseases. 79 Immunoinformatics approaches have played a vital role in the qualitative screening of multiple genomes to predict target epitopes, which provides a rational design to enhance the T-cell immune responses in a cost-effective manner. 42, 80, 81 The possibility of antigen identification at the molecular level has enabled sensible designing of peptide vaccines. An advantage of the epitope-based vaccines is that they have the potential to induce high dosage of immunogenicity at a lower cost. Therefore, the present study aimed at the prediction of potential immunogenic epitopes from the whole genome of hrHPV45 strain for designing epitopebased vaccines.

VaxiJen is the first server developed for the prediction of antigenicity of proteins of viruses, bacteria, fungi, parasites, and tumor cells. 45 The potential antigen classification is completely based on the physicochemical properties of proteins and not related to sequence alignment. Therefore, we utilized VaxiJen tool in the present study and predicted potential antigenic proteins with 87% accuracy. The reliable predictions of CTL epitopes are vital for consistent drug development. 82, 83 The CD8 + receptors consist of single CD8α and CD8β chains. Furthermore, the CD8 + CTLs identify the epitopes which are presented by the MHC class I established on all the nucleated cells. This process is facilitated by the binding of CD8 + to α3 conserved region of MHC class I. Recently, elution studies by Bassani-Sternberg et al 84 proved that the 9mer peptides are naturally optimized by the antigen processing process. This result supports the present study.

MHC class I-presenting peptides are synthesized from endogenous processing pathway-derived proteins. 85 We predicted 50 potential CD8 + MHC class I-restricted 9mer epitopes based on endogenous processing pathway (proteasomal cleavage, TAP transport, and MHC class I binding algorithm) by using combined methods in the IEDB. The α-chain of MHC class I contains the peptidebinding domain, immunoglobulin-like domain, and a transmembrane region. The heavy chain of the MHC class I molecule is coded by the genes at HLA (A, B, and C) loci. The MHC genes are highly polymorphic in nature; it means the residues in the N-and C-terminal regions of MHC proteins could generate clefts for specific peptide binding. MHC class I alleles exhibit broader binding specificity for peptides between 1,000 and 10,000 amino acids with high promiscuity. 86 Neural networks based on computational approaches have played a vital role in the prediction of the binding affinity of a peptide for specific MHC molecules. 87 In the present study, the FTFPHAFPF, CVSYWGVYY, LTAEVMSYI, and FSDPSIIEV antigenic epitopes were found to interact with more number of HLA alleles. This specific high-binding affinity is absolutely desired because the efficiency of an epitope in vaccine development greatly depends on the extent of its interaction with HLA alleles. Moreover, activated CD8 + cells kill the targets through three crucial mechanisms including secretion of antitumor and antiviral (IFN-γ) substances, 88 membrane disruption by granulysin, and Fas-mediated killing. 89 Recent studies reported the prevalence rate of multiple hrHPV genotypes in cervical cancer 90 and the poor survival rate associated with invasive cervical carcinomas. Moreover, the hrHPV strains co-infection has been shown to increase the death rate of HIV patients. 91 Conserved epitopes provide powerful immunization for prolonged period. Sabah et al 92 found that the single 9mer CTL epitope ETSVHEIEL was conserved (100%) in all E6 proteins of HPV58 isolates. While designing an epitope-based vaccine, the use of conserved epitopes may be helpful to provide broader cross-protection across multiple strains. The present study demonstrated that 14 immunogenic CD8 + epitopes showed cross-protection against more than 12 hrHPVs.

Clifford et al 93 reported that the HPV45-positive invasive cervical cancer rate was higher in Africa (15.5%) and Europe (6.3%). The highest grade of cervical cancer lesions was also recorded in Asia (2.8%) and Latin America (8.1%). Moreover, cervical cancer is the second fatal cancer contributing to 14% of cancers in Indian women. 94 The American Cancer Society estimated that there were 12,820 cases of invasive uterine cervical cancer in 2017. 95 The agestandardized incidence rate of cervical cancer in the overall geographical region of China is around 98.9%. 96 

Kaliamurthi et al findings showed that the highest PC of the pooled epitopes was recorded in Europe, North America, West Indies, West Indies, North Africa, East Asia, Northeast Asia, West Africa, South Asia, Southwest Asia, Southeast Asia, Oceania, East Africa, Central Africa, South America, Central America, and South Africa. These PC findings showed that the peptide vaccine from HPV45 would be effective for immunizing a huge population over an extensive geographical region.

There is a lack of scientific reports on binding cleft interactions of CD8 + epitopes of HPV45 with HLA molecules. Considering this fact, the binding affinity of the potential 27 immunogenic CD4 + epitopes for receptors such as HLA-A*(02:06, 0201, 11:01, 68), HLA-B*(35:01, 58:01), and HLA-C*05:01 was analyzed using docking tools. It is notable that the prophylactic vaccines from HPV L1 protein are used to treat HPV infections in more than 100 countries. 97 In the present study, we found strong binding affinity between the predicted L1 epitope LTAEVMSYI and three different MHC class I alleles HLA-A*(02:06, 02:01, 68, 58:01). Moreover, studies reported that multi-epitope HPV16 E6/E7 vaccine enhanced the clearance of HPVpositive cancer in the mouse model. 98, 99 TA-CIN vaccine was synthesized from a single fusion protein of HPV16 (L2, E6, and E7), and the clinical Phase I study in healthy volunteers showed significant production of IFN-γ which confirmed T-cell immune response in patients. 100 In the present study, minor capsid protein FSDPSIIEV epitope showed affinity for a peptide-binding groove with allele-specific pockets in MHC class I HLA-A*02:01, HLA-A*01:01, HLA-A*68:02, and HLA-C*15:02 alleles. In the case of E6 CD8 + epitopes, RTEVYQFAF showed potential binding interactions with HLA-A*32:01 and HLA-B*58:01 with the highest global binding energy of -61.05 and -58.52 kcal/ mol, respectively. Moreover, E2 is the violent regulator of extra-chromosomal replication and E1 protein expression of the viruses. In addition, the E2 proteins from hrHPV cause apoptosis in HPV-transformed cells compared to E2 from low-risk HPV strains. From the present findings, the E1 epitope FTFPHAFPF exhibited excellent affinity for a peptide-binding groove with allele-specific pockets in MHC class I HLA-A*(32:01, 26:01, 68:02), HLA-B*15:01, and HLA-C*(07:02, 14:02) alleles. The E2 CD8 + epitope CVSYWGVYY exhibited excellent affinity for MHC class I HLA-A*(29:02, 30:02, 11:01, 68:01) alleles.

Both MHC class I and II molecules share some superficial similarities and bind to the 9mer peptides. 101 However, vital differences exist between the MHC class I and II molecules. The capped nature of the MHC class I peptide-binding groove with allele-specific pockets does not allow variation in the length of peptides (9mer), which is called register shifting. This is because the peptide-binding groove of MHC class I molecules is closed at each end. In contrast, the peptidebinding groove of MHC class II molecules is open at both ends which allows for binding of more extended peptides (ranging from >15 to 24 AA in length). Sercarz and Maverakis 102 reported that longer peptides have peptide-flanking residues that lie outside of the peptide-binding groove in MHC class II molecules and might interact with peptides in another distal location. MHC class II proteins primarily present peptides derived from endocytosis of extracellular proteins (exogenous processing pathway). In the present study, promising 15mer CD4 + epitopes were predicted by consensus approaches from the late and early proteins of HPV45. The predicted CD4 + epitopes possessed good binding affinity (<1% lowest percentile rank) for MHC class II alleles. CD4 + immune responses are associated with the production of IFN-γ or IL-2. Furthermore, the IFN-γ-producing CD4 + epitopes were found by using machine learning hybrid method. B cell-mediated humoral immunity involves the recognition of antigens circulating in the body fluid. Recently, Paolini et al 103 reported the production of CD8 + and CD4 + lymphocytes induced by anticancer pVAX-E5CP and pVAX-E5MultiCP E5 vaccine (carrying the whole E5 gene or multi-epitope) in the preclinical cancer model. Bristo et al 104 reported that immunization with a single peptide that contained both the CD4 + and CD8 + CTL epitopes (ie, 9mer CD8 + MHC class I-restricted peptides nested within the 13mer CD4 + epitopes) elicited the response and production of T cells in an animal model. We also predicted the CD8 + -overlapped (9mer) epitopes within the ideal B-cell (16mer) and CD4 + cell (15mer) epitope regions. The results suggested that the CD8 + , nested CD4 + , or B-cell epitopes may elicit CTL-induced cellular immunity or T cell-induced adaptive immunity or B cellmediated humoral immunity. These identified epitopes from HPV45 may have substantial implications in the peptide vaccine-based immunotherapy.

Each HPV type is immunologically distinct, which is necessary to find suitable target for the development of therapeutic vaccines. We conclude that the predicted CD8 + , nested IFN-γproducing CD4 + , and linear B-cell epitopes can possibly induce both cell-mediated and humoral immunity by CTLs, Th1 cells, and B cells in the virus-infected host system. Moreover, the CD8 + epitopes from HPV45 proteins have the potential to be effectively used for the development of unique cross-protective Biologics: Targets and Therapy 2018:12 submit your manuscript | www.dovepress.com

CD8 + and its overlapped CD4 + epitopes of HPV peptide vaccines against hrHPV-attributed cervical cancer. However, further experimental peptide sensitization studies are needed to confirm our findings on humoral and cell-mediated immune response of the overlapped epitopes. This may enable designing of novel peptide-based therapeutic vaccines with improved efficiency for better protection.

Solution structure of the partially folded high-risk human papilloma virus 45 oncoprotein E7

Retrospective International Survey and HPV Time Trends Study Group. Human papillomavirus genotype attribution in invasive cervical cancer: a retrospective crosssectional worldwide study

The Papillomavirus E2 proteins

The role of human papillomavirus vaccines in cervical cancer: Prevention and treatment

The role of human papillomaviruses in oncogenesis

Human papillomavirus: gene expression, regulation and prospects for novel diagnostic methods and antiviral therapies

CD8+ T cell effector mechanisms in resistance to infection

Molecular interactions of coreceptor CD8 and MHC class I: the molecular basis for functional coordination with the T-cell receptor

The ins and outs of MHC class II-mediated antigen processing and presentation

Peptide-based synthetic vaccines

Use of 9-valent human papillomavirus (HPV) vaccine: updated HPV vaccination recommendations of the advisory committee on immunization practices

The human papillomavirus (HPV) vaccine, HPV related diseases and cervical cancer in the post-reproductive years

Human papilloma virus vaccines: Current scenario

Pharyngeal-cervical-brachial variant of Guillain-Barre syndrome

Knowledge of HPV and HPV Vaccine among Women Ages 19 to 26. Women's Health Issues

Don't forget HPV-45 in cervical cancer screening

International Agency for Research on Cancer Multicenter Cervical Cancer Study Group. Epidemiologic classification of human papillomavirus types associated with cervical cancer

Population-based study for human papillomavirus (HPV) infection in young women in Japan: A multicenter study by the Japanese human papillomavirus disease education research survey group (J-HERS)

Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical lesions: a meta-analysis update

HPV distribution in cervical cancer in Portugal. A retrospective study from 1928 to

The causal relation between human papillomavirus and cervical cancer

Classification of papillomaviruses

Characterization of human papillomavirus type 45, a new type 18-related virus of the genital tract

Evolution and taxonomic classification of alphapapillomavirus 7 complete genomes: HPV18, HPV39, HPV45, HPV59, HPV68 and HPV70

Prevalence of human papillomavirus genotypes among women with cervical cancer in Ghana

High-risk human papilloma virus genotypes in cervical carcinoma of Serbian women: Distribution and association with pathohistological findings

Differential presence of Papillomavirus variants in cervical cancer: an analysis for HPV33, HPV45 and HPV58

Human papillomavirus infection in females with normal cervical cytology: Genotyping and phylogenetic analysis among women in Punjab, Pakistan

Prevalence and genotype characterization of Human Papillomaviruses among HIV-seropositive in

Prevalence of human papillomavirus genotypes and relative risk of cervical cancer in China: a systematic review and meta-analysis

Genetic diversity of human papillomavirus types 35, 45 and 58 in cervical cancer in Brazil

Nested multiplex PCR based detection of human papillomavirus in cervical carcinoma patients of North-East India

HPV DNA methylation at the early promoter and E1/E2 integrity: A comparison between HPV16, HPV18 and HPV45 in cervical cancer

Unique LCR variations among lineages of HPV16, 18 and 45 isolates from women with normal cervical cytology in Ghana

Human papillomavirus types in 115,789 HPV-positive women: a meta-analysis from cervical infection to cancer

Towards peptide vaccines against Zika virus: Immunoinformatics combined with molecular dynamics simulations to predict antigenic epitopes of Zika viral proteins

T-cell epitope vaccine design by immunoinformatics

Identification of a naturally processed HLA-A*0201 HPV18 E7 T cell epitope by tumor cell mediated in vitro vaccination

Identification of novel epitopes from human papillomavirus type 18 E7 that can sensitize PBMCs of multiple HLA class I against human cervical cancer

Sequence-based approach for rapid identification of cross-clade CD8+ T-cell vaccine candidates from all high-risk HPV strains

Identification of the impact on T-and B-cell epitopes of human papillomavirus type-16 E6 and E7 variant in Southwest China

Computational identification, characterization and validation of potential antigenic peptide vaccines from hrHPVs E6 proteins using immunoinformatics and computational systems biology approaches

Use of 9-valent human papillomavirus (HPV) vaccine: updated HPV vaccination recommendations of the advisory committee on immunization practices

UniProt: the universal protein knowledgebase

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

Large-scale validation of methods for cytotoxic T-lymphocyte epitope prediction

Modeling the MHC class I pathway by combining predictions of proteasomal cleavage, TAP transport and MHC class I binding

Properties of MHC class I presented peptides that enhance immunogenicity

Development of an epitope conservancy analysis tool to facilitate the design of epitopebased diagnostics and vaccines

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

PEPstr: a de novo method for tertiary structure prediction of small bioactive peptides

Improved PEP-FOLD approach for peptide and miniprotein structure prediction

The Phyre2 web portal for protein modeling, prediction and analysis

TASSER-Lite: an automated tool for protein comparative modeling

PatchDock and SymmDock: servers for rigid and symmetric docking

FireDock: a web server for fast interaction refinement in molecular docking

Enhancing UCSF Chimera through web services

A systematic assessment of MHC class II peptide binding predictions and evaluation of a consensus approach

Designing of interferon-gamma inducing MHC class-II binders

Prediction of continuous B-cell epitopes in an antigen using recurrent neural network

AllergenFP: allergenicity prediction by descriptor fingerprints

Proteome-guided search for influenza A B-cell epitopes

UniProt Consortium. A fast Peptide Match service for UniProt Knowledgebase

Human papillomavirus type 16 E7 peptide-directed CD8+ T cells from patients with cervical cancer are cross-reactive with the coronavirus NS2 protein

Disruption of either the E1 or the E2 regulatory gene of human papillomavirus type 16 increases viral immortalization capacity

The role of globular heads of the C1q receptor in HPV 16 E2-induced human cervical squamous carcinoma cell apoptosis is associated with p38 MAPK/JNK activation

Roles of PI3K/ Akt and c-Jun signaling pathways in human papillomavirus type 16 oncoprotein-induced HIF-1α, VEGF, and IL-8 expression and in vitro angiogenesis in non-small cell lung cancer cells

Cyclooxygenase-2 transcription is regulated by human papillomavirus 16 E6 and E7 oncoproteins: evidence of a corepressor/coactivator exchange

GRIM-19 disrupts E6/E6AP complex to rescue p53 and induce apoptosis in cervical cancers

Multiplex Identification of Human Papillomavirus 16 DNA Integration Sites in Cervical Carcinomas

The Biology and Life-Cycle of Human Papillomaviruses

12 submit your manuscript | www.dovepress.com Dovepress Dovepress Biologics: Targets and Therapy CD8 + and its overlapped CD4 + epitopes of HPV

P53 and PCNA is Positively Correlated with HPV Infection in Laryngeal Epitheliopapillomatous Lesions in Patiets with Different Ethnic Backgrounds in Xinjiang

Regulation of p14ARF expression by HPV-18 E6 variants

The prognostic signifi cance of p16, Ki-67, p63, and CK17 expression determined by immunohistochemical staining in cervical intraepithelial neoplasia 1

CellDesigner: a process diagram editor for gene-regulatory and biochemical networks

Deciphering the Biochemical Pathway and Pharmacokinetic Study of Amyloid βeta-42 with Superparamagnetic Iron Oxide Nanoparticles (SPIONs) Using Systems Biology Approach

COPASI--a COmplex PAthway SImulator

The Global Burden of Cancer

Rare diseases: Drug discovery and informatics resource

Screening and structure-based modeling of T-cell epitopes of Nipah virus proteome: an immunoinformatic approach for designing peptide-based vaccine

Computer-aided analysis of structural properties and epitopes of Iranian HPV-16 E7 oncoprotein

A phase I trial of a human papillomavirus (HPV) peptide vaccine for women with highgrade cervical and vulvar intraepithelial neoplasia who are HPV 16 positive

Rates of CTL killing in persistent viral infection in vivo

Mass spectrometry of human leukocyte antigen class I peptidomes reveals strong effects of protein abundance and turnover on antigen presentation

Allele-dependent processing pathways generate the endogenous human leukocyte antigen (HLA) class I peptide repertoire in transporters associated with antigen processing (TAP)-deficient cells

Computational methods for prediction of T-cell epitopes -a framework for modelling, testing, and applications

Prediction of peptides binding to MHC class I and II alleles by temporal motif mining

Antigen-Specific Regulation of T Cell-Mediated Cytokine Production

Gamma Interferon Expression in CD8+ T Cells Is a Marker for Circulating Cytotoxic T Lymphocytes That Recognize an HLA A2-Restricted Epitope of Human Cytomegalovirus Phosphoprotein pp65

Lillsunde-Larsson G. Findings of multiple HPV genotypes in cervical carcinoma are associated with poor cancer-specific survival in a Swedish cohort of cervical cancer primarily treated with radiotherapy

Detection and genotyping of human papillomavirus (HPV) in HIV-infected women and its relationship with HPV/HIV co-infection

Designing of Epitope-Focused Vaccine by Targeting E6 and E7 Conserved Protein Sequences: An Immuno-Informatics Approach in Human Papillomavirus 58 Isolates

Carcinogenicity of human papillomavirus (HPV) types in HIV-positive women: A meta-analysis from HPV infection to cervical cancer

HPV and Cervical Cancer Epidemiology -Current Status of HPV Vaccination in India

Cancer statistics

Cancer statistics in China

Cervarix: a vaccine for the prevention of HPV 16, 18-associated cervical cancer

Towards a consensus on datasets and evaluation metrics for developing B-cell epitope prediction tools

Design, immune responses and anti-tumor potential of an HPV16 E6E7 multi-epitope vaccine

Enhancement of human papillomavirus (HPV) type 16 E6 and E7-specific T-cell immunity in healthy volunteers through vaccination with TA-CIN, an HPV16 L2E7E6 fusion protein vaccine

MHC class I and class II structures

Mhc-guided processing: binding of large antigen fragments

HPV 16 E5 oncoprotein is expressed in early stage carcinogenesis and can be a target of immunotherapy

Identification of a ras oncogene peptide that contains both CD4(+) and CD8(+) T cell epitopes in a nested configuration and elicits both T cell subset responses by peptide or DNA immunization

The 

All authors made substantial contributions to conception and design, acquisition of data, or analysis and interpretation of data; took part in drafting the article or revising it critically for important intellectual content; gave final approval of the version to be published; and agree to be accountable for all aspects of the work.

The authors report no conflicts of interest in this work.

Biologics: Targets and Therapy 2018:12 submit your manuscript | www.dovepress.com

Submit your manuscript here: https://www.dovepress.com/biologics-targets-and-therapy-journal Biologics: Targets and Therapy is an international, peer-reviewed journal focusing on the patho-physiological rationale for and clinical application of Biologic agents in the management of autoimmune diseases, cancers or other pathologies where a molecular target can be identified. This journal is indexed on PubMed Central, EMBase, and Scopus. The manuscript management system is completely online and includes a very quick and fair peer-review system, which is all easy to use. Visit http://www.dovepress.com/testimonials.php to read real quotes from published authors.