CORD-19:ffe133ed880d6c77ae340c5374817232e21f8315 JSONTXT 9 Projects

Rational Design of Peptide Vaccines Against Multiple Types of Human Papillomavirus Supplementary Issue: Novel Data Analytics in Biomedical Informatics and Translational Medicine Abstract Human papillomavirus (HPV) occurs in many types, some of which cause cervical, genital, and other cancers. While vaccination is available against the major cancer-causing HPV types, many others are not covered by these preventive measures. Herein, we present a bioinformatics study for the designing of multivalent peptide vaccines against multiple HPV types as an alternative strategy to the virus-like particle vaccines being used now. Our technique of rational design of peptide vaccines is expected to ensure stability of the vaccine against many cycles of mutational changes, elicit immune response, and negate autoimmune possibilities. Using the L1 capsid protein sequences, we identified several peptides for potential vaccine design for HPV 16, 18, 33, 35, 45, and 11 types. Although there are concerns about the epitope-binding affinities for the peptides identified in this process, the technique indicates possi bilities of multivalent, adjuvanted, peptide vaccines against a wider range of HPV types, and tailor-made different combinations of the peptides to address frequency variations of types over different population groups as required for prophylaxis and at lower cost than are in use at the present time. Certain types of cancers such as liver cancer, cervical cancer, Hodgkin's lymphoma, and others are caused by viruses, for some of which vaccines are well established. Vaccines operate by triggering immune response in an individual. It has been shown that under some circumstances immunocompetent individuals can protect themselves against primary tumor development. 1 It is also established that while a healthy immune system can act as a tumor suppressor, in other cases it might facilitate tumor progression. 2 Nevertheless, given the wide variety of viral cancers and progression of the disease, using the immune system to prevent cancer through antiviral vaccines could be a reasonable procedure. Some viral cancers such as those caused by human papillomavirus (HPV) and hepatitis B viruses have been effectively prevented through such techniques. Papillomaviruses affect mammalian, avian, and reptilian hosts. Over 170 types of HPVs have been identified. 3 The large variety of HPV types is believed to have evolved from the large-scale migration of human populations and appears to follow their major migration patterns. 4 The virus is generally spread by sexual activities and can lie dormant for years before symptoms develop. In the majority of cases, the virus causes warts on the body, but a small fraction of the infections develops into cancers. It is estimated that over 27,000 people develop HPV cancers every year in the USA, usually detected at a fairly late stage 5 ; an analysis of all cancers worldwide in 2002 revealed that of about 1.9 million cases of cancer attributable to infections, around 561,200 were due to HPVs, 6 with an incidence ratio of 2:7 between numbers in the developed and developing countries. The wide varieties of HPVs are causative agents for various types of cancers and other diseases. Of the 40 types of the HPVs that affect genital areas, 7 HPV types 16, 18, 33, 35 , and 45 are believed to be the primary agents, with HPV16 and HPV18 identified as being responsible for about 70% of such cancers. One study conducted among a large sample of patients with histologically confirmed squamous cell cervical cancers 8 determined 15 HPV types (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 68, 73 , and 82) classified as high risk, three (26, 53 , and 66) as probable high risk, and several others (6, 11, 40, 42, 43, 44, 54, 61, 70, 72, 81) as low-risk types. Some HPV types also cause oropharyngeal cancers: 90% of HPV-positive oropharyngeal cancers are attributable to HPV16, 8, 9 and HPV6 and HPV11 have also been determined to cause such cancers 10 ; these findings are important since HPV-positive oropharyngeal squamous cell cancers (OPSCCs) are reported to have doubled in the period 1984-2004, whereas HPV-negative OPSCCs decreased by 50% while overall 18% of oral and oropharyngeal cancers worldwide are HPV associated. 10 Other types of HPVs are known to cause various ailments: HPV6 and 11 cause recurrent respiratory papillomatosis, 11 HPV types 2, 27, and 57 can cause skin warts, while HPV 1, 2, 27, and 63 were found on clinically normal skin. 12 Moreover, the proportion of various types of HPV varies across countries and continents, and thus affects the outcome of standard preventive medications. 8 HPV is a DNA virus with a genome length of about 8000 base pairs comprising nine segments that code for eight proteins and one long control region. There are basically six early proteins and two late proteins, namely, L1 and L2, which are outer capsid proteins, at approximately 95:5 ratio. The L1 proteins form a pentameric assembly on the virion shell, the individual proteins being held together by disulfide bonds; these proteins form the basis for vaccines against papillomaviruses. Unlike the RNA viruses that lack proofreading mechanism, DNA viruses are capable of correcting replication errors and are much more stable to mutational changes than their RNA counterparts. Thus, DNA viruses are more conserved, and vaccines are more effective over a longer period. The direct linkage of some types of HPV to cervical cancers has led to much research on prophylactic vaccination against the most common HPV types that are causative agents for cervical cancers. 13 Antiviral vaccines have been developed against the most common of the HPV types that cause these cancers, viz., HPV16 and HPV18; Gardasil (Merck & Co., Inc.) protects against high-risk HPV types 16 and 18 and lowrisk types 6 and 11; Cervarix (GSK) protects against HPV types 16 and 18 through an AS04 adjuvanted recombinant L1 protein vaccine, which was found to also protect against HPV types 31 and 45. 13 The efficacy of these vaccines is very high but limited to only the selected HPVs due to type specificity, inhibiting their preventive action against a wide range of other HPV-related cancers 3 ; the two vaccines leave out a coverage of 20%-30% of the HPV-associated cancer-causing types. This situation has been mitigated to some extent by the launch of a 9-valent Gardasil vaccine (9vHPV) that protects against HPV types covered by the older quadrivalent HPV vaccine (HPV6, 11, 16, and 18) and five additional oncogenic types (HPV31, 33, 45, 52, and 58). 14 A controlled test determined that the 9vHPV was noninferior to the previous quadrivalent HPV in preventing infections and generating antibodies against the particular types of HPVs it was designed against, but was not effective against other types. 14 The HPV vaccines belong to the class of recombinant virus-like particle (VLP) vaccines, which are self-assembling bionanoparticles that mimic the structure of the original virus particle exposing multiple epitopes on the surface. They do not carry any genetic material, thus providing a safer alternative to live-attenuated or inactivated vaccines. The first VLP vaccine license was granted in 1986 for the hepatitis B virus by the US Food and Drug Administration, 15 and VLP vaccines for two more viruses, such as HPV and hepatitis E, have secured licenses. The fact that over the last 30 years only three such vaccines have been marketed while many new and re-emerging viruses have been recognized indicates how difficult and expensive it must be to develop VLP vaccines. Although the VLP vaccines have proved very effective in practice, they are generally beyond the reach of patients in developing countries, and also the nature of the VLP vaccine can inhibit attempts at personalized medication, which may be required in certain population groups. 16 The availability of genomic data, understanding of immune responses and immunogenetic variations, new developments in bioinformatics and computer applications, and systems biology approach are fashioning new directions in vaccine development, away from the traditional techniques of live-attenuated or inactivated virus or recombinant VLP vaccines. This is a move away from the one-size-fits-all concept of mass-applicable vaccines to one of personalized vaccines where the vaccines, ideally, can be readily designed to be individual, gender-specific, or race/community/population-specific and keep costs in control. 16, 17 This strategy is admittedly geared at the moment to hypervariable viruses such as human immunodeficiency virus (HIV), human coronavirus (HCV), influenza, and others, to which group the HPV does not belong, but the principles certainly can apply. Rational design of vaccines based on genomic and immunogenomic information and the science of reverse vaccinology, as yet an evolving technology, could be a pointer to the future. At the immediate level, peptide vaccines may be such an alternative strategy to HPV VLP vaccines that could be explored. The idea is to scan the viral genome for the protein antigens that can elicit immune response and then synthesize them into a peptide vaccine. A more focused approach is to precisely locate the epitope regions within these antigens and utilize them to elicit the immune response. 18 The recent advancements in the technological and bioinformatics fields enable computer-based approaches for this purpose. Synthetic peptide vaccines, which are just beginning to draw wide interest at this time, are easier to construct, the bioinformatics approaches lead to shorter development times, such vaccines can be tailor-made to suit individual applications and avoid autoimmune threats, and overall are cheaper to design and market. 18 It is interesting to note that the US National Institute of Health projects portal, Clini-calTrials.gov (https://clinicaltrials.gov/), reports, as of mid April 2016, 562 projects under trial for peptide vaccines, 106 projects for VLP vaccines, and 74 projects for synthetic peptide vaccines. We have carried out computer-based studies to identify target regions in surface proteins of influenza and rotavirus. [19] [20] [21] Our interest in applying these techniques to HPV is to see whether it is feasible to design peptide vaccines against the disease that afflicts many more people in the developing nations where the economic burden is much more acute than in the developed nations. There is also the advantage with peptide vaccines that, in theory, multivalent vaccines can be more easily designed and marketed at more affordable prices than multivalent VLP vaccines against a wide range of HPV types. This makes it possible to tailor-make different combinations of the peptides to address frequency variations of types over different population groups and therefore reach a wider percentage of affected people to suit individual communities in those countries where the populations are higher and much less homogeneous than in the Western nations. database. In this work, L1 being the major component of the capsid assembly and the focus of HPV vaccine development, we restrict our approach to consideration of the L1 protein only and postpone consideration of the other capsid assembly component, the L2 protein, to a future work. We have downloaded 222 sequences of the L1 capsid gene of various HPV types from GenBank, updated till March 22, 2016 . Table 1 lists the number of sequences of the different HPV types forming the dataset we have used for the purposes of the present exercise. Annotated HPV data are also available through the specialized Papillomavirus Episteme database maintained by the National Institute for Allergic and Infectious Diseases. Given the high incidence of papillomavirus infections, one could expect more data to accumulate over time. Method of analysis. The method of analysis has been described in detail in our previous papers, especially in the study by Ghosh et al. 20 Initially, we scan a library of sequences of the designated protein to determine segments that are unchanged or least changed among the various strains, and then, couple these with average solvent accessibility (ASA) profiles of the sequences to select those segments that are most conserved and having highest solvent accessibility profile. We then determine epitope regions within these segments with acceptable binding affinity to human leukocyte antigens and finally select those peptides that poses no autoimmune threat; in the process, the 3D crystal structure of the surface protein is also utilized to ensure that the selected peptide regions are not covered against solvent accessibility due to neighboring proteins that together constitute the quaternary structure. As an optional step, the process may start with observations from graphical representations of DNA/RNA sequences in a 2D Cartesian coordinate system to determine similarities between the selected gene sequences or conserved segments within the selected sequences. In this system, a sequence is represented by a point for each nucleotide in succession plotted in a series of steps with one step in the negative x-direction for an adenine, positive y-direction for a cytosine, in the positive x-direction for a guanine, and in the negative y-direction for a thymine. 22 This generates a plot of the sequence that reflects the base distribution in the sequence. For a quantitative measure, we take the weighted center of mass (µ x , µ y ) of the points and the distance of the centre of mass (c.m.) from the center as a graph radius g R , which is considered as the characteristic of the sequence. 23 For protein sequences, this concept is extended to a rectangular coordinate system of 20 dimensions, where each amino acid (aa) is associated with a particular direction and the protein sequence is plotted by taking successive steps in the directions dictated by the aa sequence. 24 The plot is in abstract space, but we can again define a weighted center of mass and a protein graph radius, p R , just as for g R . The special property of g R and p R is that in each case when two sequences result in the same graph radius, albeit with some adjustments for p R as specified in the study by Ghosh et al. 20 , the sequences are found to be identical. 25 This property has been used in several instances previously 19, 20, 26 and is used in the current work to determine protein sequence segment varieties by determining the number of instances where the p R 's are identical (see more details below). To determine the target peptides that can be considered to elicit immune response, the steps in the process are, briefly: Global view of L1 gene sequences to determine similar and dissimilar sequences. 2D graphical representation plots are generated for the L1 gene sequences of different HPV types and visually inspected to determine those that appear similar and can be clubbed together for analysis. Identification of protein conserved segments. Viruses generally have high mutation rates, more for RNA viruses than for DNA viruses. The mutation rates for the dsDNA HPVs are very low, but nevertheless, mutations do occur leading to the various strains that populate the papillomavirus database -higher differences between types and smaller differences within each type; a phylogenetic tree of the protein sequences of the primary HPV types we are interested in this work, with five strains from each type, brings out this feature explicitly (Fig. 1) . Now, given the strict type specificity of the neutralizing antibodies and the mutations in aa composition, deletions, and insertions that may take place, 27 this implies that in specific cases neutralizing antibodies raised against one strain may not be viable against another strain, although the success of the current HPV vaccines show that this may happen seldom. Our interest is to identify those peptide segments that change the least between strains to reduce the possibility of such instances and examine whether one or more of such segments could form the basis for peptide vaccines. The methodology we follow for the purpose is to map the protein sequence using a window of 12-14 aas and sliding in increments of one aa at a time; ie, compute p R for the segment with aa nos. 1-12, say, for a chosen window size of 12 aa, then another p R for the stretch between aa nos. 2-13 and so on, until the end of the aa sequence is reached. At each stage, the peptide graph radius, p R , is computed and stored (see Refs. 19 and 20 for details). This is then repeated for each protein sequence of the selected group of sequences and then the entire lot is scanned at each mapped aa position to determine the number of different p R -values, among all the sequences. Since each p R value represents a particular peptide sequence, the number of values of the p R at each aa position map out a protein variety profile of the aas covered by the sliding window. The minima of the profile indicate regions where variety, and by inference changes in aas, are relatively least, indicating most conserved regions. The size of the window covering a number of aas needs careful consideration. In our previous analysis of the neuraminidase protein, 19 we had considered peptide lengths of 8, 10, and 12 aas; for rotavirus where mutational changes occur very frequently, it was considered prudent to consider window sizes of 12-14 aas 20 ; 14 residue peptides were found to be more potent antigenically compared with smaller peptides, 28 and we had used that size for the analysis of influenza hemagglutinin. 21 It may be noted that while peptide lengths of 6-20 aas are used in scanning for B-cell epitopes, peptides of 10+ residues may contain overlapping linear B-cell epitopes 29 ; cytotoxic T-cell epitopes have a limited length (8-11 residues), whereas for helper T-cells, longer peptide lengths are used. 30 Considering these issues and the fact that the HPV viruses are dsDNA type, which therefore are comparatively more stable, we fixed window length for this exercise at 12 aas. Solvent accessibility profile. The protein sequences are next mapped for solvent accessibility by subjecting them to an ASA server such as ITASSER (http://zhanglab.ccmb.med.umich. edu/I-TASSER/) or SABLE (http://sable.cchmc.org/) or any other suitable web-based server. Taking an average over 12 aas of solvent accessibility index at each aa position, and averaging over several sequences for each HPV type, this is mapped out to reveal the most solvent accessible regions by the maxima of the graph. First cut peptide segments. Mapping both the graphs, namely, peptide variability profile and solvent accessibility profile, together, we scan for those regions where the protein variability is least and solvent accessibility is among the highest. These are the first list of conserved solvent accessible segments of the protein. Confirmation from 3D structures and second cut choice. To confirm that the regions identified in the previous step are indeed surface situated, we examine the identified regions in a protein 3D structure, eg, as shown through a software such as Cn3D4.3 available at NIH website. 31 Marking out the identified segments on one protein of the 3D structure will reveal how much of the segments lie on the surface. If composite structures of this protein with the other proteins are available, then we can further fine-tune the selection process to identify those segments that really lie on an exposed surface of the capsid protein and not get covered by a neighboring protein of the virion. Determination of epitope segments. Next, we consider which of the segments from the second-stage choice contains epitope regions to elicit necessary immune response. For this, we use two epitope-prediction tools, viz., the Immune Epitope Database (IEDB) and Analysis Resource server [32] [33] [34] [35] and ABCpred server, 36, 37 which, as mentioned earlier, 21 have given good results. These tools indicate binding affinities of T-cell epitopes and B-cell epitopes (with 66% accuracy), major histocompatibility complex (MHC) I, and also of discontinuous epitopes through Ellipro. Autoimmune test. Finally, we test the peptides that have remained after the third stage above to ensure that they do not possess any autoimmune threats and are unique peptides by themselves. Each peptide segment is subjected to a proteinprotein BLAST to ensure that the selected peptides are not homologous to human proteins to cause autoimmune threats. These segments then become our recommendations for peptide vaccines subject to final verification in wet laboratories. That autoimmune threat can be a potential hazard has been explored in one of our earlier reports. 38 Our initial attempt at determining a common motif of conserved peptide sequences in the L1 gene of HPV, of all types available in our database as given in Table 1 , met with little success. While signals for low variability were observed, these were heavily influenced by the preponderance of HPV16 data in our database; in fact, detailed inspection of the peptide stretches corresponding to the signal regions showed little homology between HPV16 and HPV18 sequences. This is observable from the nucleotide plots also (Fig. 2) , which shows close similarity between HPV16, HPV33, and HPV35, and between HPV18 and HPV45, but the two groups have distinctly different plots, indicating significant differences in base distribution. In our previous analyses for vaccine targets, we had considered large numbers of sequences of only one primary type, viz., the neuraminidase 19 and hemagglutinin 21 of influenza A genomes, or VP7 proteins of rotavirus. 20 In the present instance, we have to consider different types of the HPV where the similarities within sequences of each type are greater than similarities between sequences of different types; we note that, as remarked by Chen et al. 39 , if the difference of the L1 gene of a cloned sequence from any other characterized type becomes 10% or more, a new type is claimed to have been established, else these are all strains of the same type group. The phylogenetic tree (Fig. 1) shows explicitly that the strains of each type belong to individual groups and the different types form separate clades: the HPV18 and 45 belong to one clade and HPV16, 33, and 35 to another. The 2D graphs (Fig. 2 ) also show at a glance the close similarities between HPV16, 33, and 35, and also between HPV18 and 45, and indications of where the two groups are different. Thus, our earlier experience could not be repeated here, which has to be considered as intertype differences. However, there can be large differences within a type too: eg, HPV type 16 from Iran (KM058666) and African type 2 variant (AF472509) are quite different, but because these are dsDNA types of viruses, we expect they will be regionally, ie, within the population groups, stable. Considering the preponderance of data, only the former variety of HPV16 and related sequences were used for our analyses. The analyses for the three high-risk types, HPV16, 33, and 35, were done together using the p R technique mentioned in the "Materials and methods" section with a 12 aa window, and the p R -values at each aa position were scanned to determine the number of different values. These numbers were then smoothed over a range of 12 positions, and a running average computed over the complete sequence. Next, ASA profiles were determined for several samples of the selected HPV types using the SABLE server, which has been found to be quite reliable. Taking an average of the solvent accessibility numbers determined residuewise for each sequence from SABLE, and then taking the running average of 12 aa, we get a smoothed out solvent accessibility profile whose peaks indicate highest solvent accessibility around that residue position. By our logic, the regions where the lowest variety numbers for aa stretches match with the highest values of the ASA profiles represent potentially the best situated highly conserved surface accessible peptide segments in the protein sequence. Figure 3 shows the plots of these two profiles of the high-risk HPV16, 33, and 35 types, where the five regions so identified are marked by short horizontal lines for further analysis. To identify in terms of aas the regions we had determined in the varieties vs ASA profiles, we scanned the aa sequences of the HPV16, 33, and 35 of those regions. The three HPV types were found to have very close similarities in the peptide stretches, extending, in fact, to several aas beyond the 12 aa windows. aa variety vs ASA profile curves were also drawn similarly for HPV18 and HPV45 and appropriate peptide segments were identified. Table 2 lists the peptides for all identified segments for the five high-risk HPV types, and also of the low-risk HPV11 determined by a similar method. Since the HPV virion capsomers are composed of five L1 proteins, and one L2 each, it is possible that a surface accessible antigen in one protein may get covered by the folds of a neighboring protein. Step three of our search for vaccine target regions is to verify the surface accessibility of the identified regions through protein crystal structure. For this exercise, we used 2R5H for HPV16, 2R5I for HPV18, 2R5J for HPV35, and 2R5K for HPV11 available from the Protein Data Bank and individually checked for each of the regions in space-fill models; since no structure was available for HPV33 and HPV45, detailed comparison had to be omitted here. Figures 4-7 (HPV16, HPV18, HPV35 , and HPV11, respectively) give snapshots of the identified peptides highlighted in yellow in one of the proteins of the pentameric structure of the capsid assembly. In the case of region 1 of HPV16, which was surface accessible for a solitary HPV protein, this peptide stretch is almost entirely covered by the adjacent protein and, therefore, has to be rejected as a vaccine target candidate. Region 2 was found to be entirely viable. In the case of regions 3 and 4, the original 12 aa stretches were found to be surface situated, but 5 aas at the 5′-end were found to be covered by the neighboring residues, ie, even if these extended segments belong to a group of conserved residues, the parts beyond the original 12aa stretches are not viable as candidate vaccine targets. Interestingly, the two regions, 3 and 4, are on two sides of the protein, but both are surface accessible. These segments that are not available for consideration as potential vaccine targets are marked in red in Table 2 . Region 5 listed in the table could not be tested in a similar manner due to the nonavailability of these stretches in the protein crystals. Notes: these are listed regionwise to show peptide segment similarities. the segments marked in red are covered by neighboring proteins/sequence overlaps. since HPV33 crystal structure was not available, overlap regions could not be ascertained, but can be expected to be similar to those for HPV35; similar considerations applied to HPV45 compared with HPV18. The next stage is to determine if the identified regions have epitope potential to enable immune response. For this, we used the IEDB server to test for T-cell MHC II epitope regions using two alleles, namely, HLA-DRB1-01-01 and HLA-DRB1-04-01. Table 3 shows that all the regions we had identified passed the MHC II test, although in some instances one or the other allele gave results that were just outside the permissible range for generating immune response. However, none of the regions of any of the three HPV types failed in both alleles. Likewise, we tested for B-cell humoral response through the ABCpred server, which has been known to predict epitope regions for B-cell response with about 66% accuracy. As can be seen from the results extracted in Table 4 , we find that in this instance all our regions passed this test with high acceptability ranks, except for region 4 for HPV35 whose score was outside the critical range. Both, linear and conformational epitopes, are found on the L1 protein surface, but it is well established that most of the neutralizing antibody production are through conformational epitopes. 27, 40 Ludmerer et al. 41 showed that a distinct surface-exposed linear epitope overlaps the neighboring conformational epitope in both HPV11 and HPV16 VLPs and that neutralizing epitopes also bind to a linear epitope. Combita et al. 42 experimenting with mice model confirmed the existence of linear epitopes inducing cross-neutralization, but also inferred that such cross-neutralization will not exceed 1% of the effects of the dominant conformational epitopes. Type-specific and cross-reactive linear epitopes have been determined on the L2 protein surface also, 43 but for the present purposes, we restrict our analyses to the L1 protein epitopes only. Our tests above identified linear epitopes but since the papillomavirus L1 proteins have mostly conformational epitopes, we tested for discontinuous epitopes on the aa sequences also through Ellipro available in the IEDB web server. We found that for all the HPV types we investigated, almost all the identified peptides had overlaps with the conformational epitopes. The figure 5. Display in space-fill rendering of HPV18 L1 protein pentameric structure (2R5I). The conserved surface-exposed segments identified by comparison of protein variability and ASA profile analysis as given in Table 2 peptides that had the overlaps are marked in Table 5 . Table 5 also shows the percentile rank of binding affinity for HLA for MHC Class I predicted by the IEDB server; a number of peptide regions show low ranks, which indicate high binding affinity. The remaining test was envisaged for autoimmune possibilities. In the protein-protein BLAST analysis using the nonredundant protein sequences from GenBank, PDB, SWISSPROT, etc., the results were found to match the super papillomavirus family with no matches with any non-HPV protein. The final list of identified peptide stretches that met all our criteria and could serve as potential peptide vaccine targets are given in Table 5 . We have ended up with four possible surface accessible highly conserved peptide regions in each of the HPV16, 18, 33, 35, 45 , and HPV11, with some being common to several types. Some show less binding affinities in T-cell or B-cell activity, other regions appear to have purely linear epitopes with no overlaps with any discontinuous epitopes, but by and large, every peptide segment identified in this exercise show prospects of eliciting immune response; in some instances, cross-neutralization between types seem possible. What would be required, however, are wet laboratory experiments and adjuvant combinations to validate this model of alternate strategies and improve and expand upon the possible multivalent peptide vaccine approach to the prevention of HPV-induced cancers. cervical, genital, and other cancers and the sufferings these cause, and the large variety of the virus, such preparations are to be strongly advocated. Notes: "Y" represents acceptable binding affinities. The discontinuous epitopes for HPV18 could not be determined due to mismatch with homologs.

Annnotations TAB TSV DIC JSON TextAE-old TextAE

  • Denotations: 7
  • Blocks: 0
  • Relations: 0