key: cord-0011421-17mmga74 authors: Dehghani, Behzad; Hasanshahi, Zahra; Hashempour, Tayebeh; Motamedifar, Mohamad title: The possible regions to design Human Papilloma Viruses vaccine in Iranian L1 protein date: 2019-12-24 journal: Biologia (Bratisl) DOI: 10.2478/s11756-019-00386-w sha: 2d0a7ca2d7805d3f24d8848492783f81cc79f752 doc_id: 11421 cord_uid: 17mmga74 Human Papilloma Virus (HPV) genome encodes several proteins, as L1is major capsid protein and L2 is minor capsid protein. Among all HPV types HPV-16 and HPV-18 are the most common high-risk HPV (HR-HPV) types globally and the majority of cases are infected with these types. HPV entry and the initial interaction with the host cell are mainly related to the L1 protein which is the main component of HPV vaccines. The aim of this research was comparison analysis among all Iranian L1 protein sequences submitted in NCBI GenBank to find the major substitutions as well as structural and immune properties of this protein. All sequences HPV L1 protein from Iranian isolates from 2014 to 2016 were selected and obtained from NCBI data bank. “CLC Genomics Workbench” was used to translate alignment. To predict B cell epitopes, we employed several programs. Modification sites such as phosphorylation, glycosylation, and disulfide bonds were determined. Secondary and tertiary structures of all sequences were analyzed. Several mutations were found and major mutations were in amino acid residues 102, 202, 207, 292, 379, and 502. The mentioned mutations showed the minor effect on B cell and physicochemical properties of the L1 protein. Six disulfide bonds were determined in L1 protein and also in several N-link glycosylation and phosphorylation sites. Five L1 loops were determined, which had great potential to be B cell epitopes with high antigenic properties. All in all, this research as the first report from Iran described the tremendous potential of two L1 loops (BC and FG) to induce immune system which can be used as the descent candidate to design a new vaccine against HPV in the Iranian population. In addition, some differences between the reference sequence and Iranian patients’ sequences were determined. It is essential to consider these differences to monitor the effectiveness and efficacy of the vaccine for the Iranian population. Our results provide a vast understanding of L1 protein that can be useful for further studies on HPV infections and new vaccine generations. Globally, 2784 million women aged 15 years and older are at risk of developing cervical cancer; 569,847 women are diagnosed with cervical cancer, and 311,365 die from the disease annually (Bruni et al. 2019) . After breast cancer, it is the second major cancer of women's reproductive organs related to cancer deaths. (Beckmann et al. 2010; Rahimifar et al. 2010) . Studies showed more than 69.4% of this cancer is related to HPV infection as HPV was found in pre-invasive lesions (Bruni et al. 2019) . The rate of Human Papillomavirus (HR HPV) among Iranian women has increased in the recent years; several studies have estimated this rate among healthy women from31% to 34% (Jamdar et al. 2018; Shafaghi et al. 2013; Yousefzadeh et al. 2014 ). This rate was around 38.6% in Iranian patients with cervical infections and 23% with breast cancer (Dadashi et al. 2017; Haghshenas et al. 2016) . HPV is a non-enveloped particle and like other Papillomaviruses (PVs) it is characterized by having an icosahedral capsid of about 55 nm diameter (Doorbar et al. 2015) . HPV genome expresses six nonstructural viral regulatory proteins (E1, E2, E4, E5, E6, and E7) and the virus capsid contains two encoded proteins, the L1 major and L2 minor capsid proteins (Frazer 2009; Zheng and Baker 2006) . L1 protein is highly conserved among different HPV types and only 10% of the L1 sequence is different among HPV types, and it has the ability to self-assemble in virus-like particles (VLPs) (Kim 2016; Xu et al. 2016) ; it seems that a surface cell protein was involved in VLP binding. Based on VLPs, VLP-based vaccines were introduced that offered highly effective protection against HPV infections (Buck et al. 2013 ). L1 has an indispensable role in HPV entry and the initial interaction is remarkably attributable to L1 interactions with proteoglycans (Buck et al. 2013) . This virus is categorized into low-risk or high-risk. Prior studies showed that about 50-60% of human anogenital carcinomas belonged to HPV-16 and around 20% was HPV-18(Vera-Bravo et al. 2003) . Likewise, several studies showed that the majority of HPV infections in Iranian patients were related to HPV-16. Two recommended HPV vaccines, Gardasil and Cervarix contain VLPs assembled from the L1 proteins of different HPV types ( M o ni e e t a l . 2 0 0 8 ; S i d d i q u i a n d P e r r y 2 0 0 6 ) . Consequently, any change in HPV L1 can lead to a decrease in the effectiveness and efficacy of vaccines that are used commonly. The aim of this study, as the first report of the most prevalent substitutions of HPV L1 protein among sequences obtained from Iranian patients, was to find mutations and their effect on the structure, physicochemical properties, and antigenic features. We also attempted to provide a comprehensive view of L1 protein as a major candidate designing vaccines against HPV infections. All 60 L1 sequences from 2014 to 2016 (KM058636-K M 0 5 8 6 6 6 , K P 1 6 0 9 8 8 -K P 1 6 0 9 9 9 , K P 1 6 1 0 0 0 -KP161014, and KX827590) and a reference sequence (K02718.1) were obtained from NCBI databank at (http://www.ncbi.nlm.nih.gov), all sequences were used in this study belonged to HPV-16. Table 1 shows all software which were used in this study. The CLC Genomics Workbench was used for the analysis of the mutations, alignment of the translated peptides and homology among sequences. Based on mutations, the selected sequenced were categorized into 3 groups (Table 2) . To predict signal peptide, we used "Signal-BLAST" (Frank and Sippl 2008) and "predisi" (Hiller et al. 2004 ). Prediction of instability index, aliphatic index, theoretical isoelectric point (pI), and grand average of hydropathy (GRAVY) was done by "Expasy'sProtParam" (Gasteiger et al. 2005 ). "DiANNA" (Ferrè and Clote 2006) and "SCRATCH" (Cheng et al. 2005) were used to predict disulfide bonds prediction. Phosphorylation sites were computed by "NetPhosK" (Blom et al. 2004) , "DISPHOS" (Iakoucheva et al. 2004) , and "NetPhos" (Blom et al. 1999) , and finally, N-glycosylation sites were found by "NetNGlyc" (Gupta et al. 2017) and "GlycoEP" (Chauhan et al. 2013 ). Secondary and tertiary structure "SOPMA" (Geourjon and Deleage 1995) and "Phyre" (Kelley and Sternberg 2009 ) server were reliable programs we employed to predict the secondary structure. "I-TASSER" (Roy et al. 2010) , "Phyre 2 " (Kelley and Sternberg 2009), and "(PS)2-v2" (Chen et al. 2006) were utilized to predict the tertiary structure of the selected sequences. All predicted 3D structures were evaluated for the stereochemistry, reliability, and quality by "Qmean" (Benkert et al. 2008) and "Rammpage". The obtained 3D structures were refined by "GalaxyRefine" (Heo et al. 2013) program. "Discovery Studio" was used to find the loops on HPV L1. "Immuneepitope" (Vita et al. 2018) , "BcePred" (Saha and Raghava 2004) , "ABCpred" (Saha and Raghava 2006b), and "Bepipred" (Jespersen et al. 2017 ) online programs were performed to find B Cell epitopes sites. Allergenic properties were predicted by "AlgPred" (Saha and Raghava 2006a) and "VaxiJen" (Doytchinova and Flower 2007) software computed protective antigens and subunit vaccines prediction. Hydrophilicity, flexibility, and surface accessibility of GP were anticipated using "IEDB" tools. All the selected sequences contained a mutation in amino acid 228, and the most prevalent mutations occurred in positions 102, 202, 207, 292, 379, 502 (Table 2 ). Analysis by two selected software showed that signal-blast and predisi could not find any position for signal peptide. L1 protein with 531 amino acids length and molecular weight of 59,554.02 Da had pI of 8.27 indicated that it is a basic protein. Its' in vivo half-life was estimated to be 30 h in mammalian cells and more than 10 h in E.coli. Its instability index and aliphatic index were computed to be a stable protein (Table 3) . Table 4 illustrates the disulfide bonds results, determining 6 probable disulfide bonds in all sequences; however, because of mutations, different bonds were found in some positions. Prediction of glycosylation sites by several reliable tools showed 4 conserve positions for all sequences (157,242,367, and 421) . The analysis showed that L1 is a highly phosphorylated protein and the software determined a large number of phosphorylation positions in such protein (Table 5 ). In comparison with the reference sequence, 4 new phosphorylation positions (229, 376, 507, and 514) were found in 3 groups of selected sequences. In considering the list of mutations we could not find any related mutation to these new positions. Therefore, other factors may have influenced one software to predict these positions. According to the analysis, the majority of the secondary structures of L1 were random coil (42.5%); extended strand and alpha helix were 28% and 19%, respectively (Table 6 ). Figure 1 shows the tertiary structure of L1 predicted by reliable software. The Qmean and Ramachandran plot results are summarized in Table 7 . The coverage of each predicted structure is shown in Fig. 2 . By using "discovery Studio", 5 loops were found on the L1 structure which was refined by "GalaxyRefine"; the results were summarized in Table 8 and illustrated in Fig. 3 . B cell epitope prediction by using software showed 6 highly potential regions (37-50, 79-86, 194-210, 303-313, 427-442, and 509-530) . "AlgPred" could not find any IgE epitope, and "VexiJen" showed L1 is a Probable antigen. "IEDB" results are displayed in Fig. 3 , with 3 distinct parameters which used to locate regions with the highest antigenic ability and all prediction calculations are based on propensity scales for each of the 20 amino acids. In each diagram, the regions with higher scale than the threshold were favorable regions consisting of higher scored residues. Combination of the results of all three diagrams could define the region with high potential to be a candidate as an antigen. Several studies on Iranian patients with HPV infections illustrated that the majority of cases were infected with HPV-16 and after that HPV-18 is ranked second (Farjadian et al. 2003; Ghaffari et al. 2006; Salehi-Vaziri et al. 2016; Shahramian et al. 2011) . As there were few sequences related to HPV-18, all the sequences used in this study belonged to HPV-16. Mortazavi et al. in 2002 showed that more than 73% of them were infected with HPV-16 and around 12% with HPV-18 and HPV-33 types (Mortazavi et al. 2002) . In 2014, research on 851 Iranian women showed 7.3% of them were infected with HPV-16 (Yousefzadeh et al. 2014) . While numerous investigations aimed to determine HPV and genotyping by using L1 sequences, there are no study reporting mutations in this gene. Therefore, this study was the first report on HPV L1 mutations in Iranian patients and the results are applicable to increase the efficiency of the current vaccine. Antoine Touze et al. in 1998 showed substitution in the amino acid region from residues 83 to 97 seemed to affect the level of expression of the L1 protein that has an important effect on HPV vaccines development and serological tests (Touze et al. 1998 ). In 2014, Fleury et al. found that only a few mutations within the FG loop are sufficient to generate a new serotype escaping vaccination (Fleury et al. 2014) . Shailja Pande in 2008 described that the most frequent variation in amino acid (T353P), threonine, was replaced by nonpolar aliphatic amino acid and may have an effect on the structure or function of the L1 protein (Pande et al. 2008) . It was suggested that it had a significant role in immune recognition and vaccine development strategies and led to conformational changes within epitopes relevant for viral neutralization. In 2003, Yoshiyuki Ishii found that three cysteines (175, 185, and 428) are required for the normal assembly of L1capsids through trimerization and dimerization of L1; also, substitution in the mentioned regions affects the normal function of HPV (Ishii et al. 2003) . The present study showed several mutations in L1 protein; however, we did not find any mutation similar to previous studies in other countries, probably due to different regions or different types of their studies. Lee et al. in 2008 showed the expression rate of HPV L1 capsid protein in the uterine cervical specimens and determined that HPV low-risk group had a higher L1 capsid expression rate than HPV high-risk group; they suggested that L1 capsid expression might be related to favorable disease biology (Lee et al. 2008 ). Expression of L1 in the insect cells Table 5 Phosphorylation sites found by NetPhosK", "DISPHOS" and "NetPhos; several serine, threonine, and tyrosine amino acids being predicted Serine Threonine Tyrosine Reference 49-115-244-308-324-369-375-422-518-519-521 36-65-121-292-320-362-380-439-443-544-516-523-507-517-520-522 53-161-260-268-302- by using recombinant baculoviruses system by Antoine Touze et al. (1998) in 1998 was done, showing the stability of this protein in eukaryotic cells and the yield depend on mutations in a region (83 to 97) that seemed to affect the level of expression. In 2009, Bazan et al. expressed, purified recombinant HPV16 L1 in methylotrophic yeast (Pichia pastoris) that could be useful to produce low-cost vaccines (Bazan et al. 2009 ). Wang in 1999 and Zhang in 1998 produced a fused form in Escherichia coli using an inducible expression system and showed the stability of the L1 protein in this host (Zhang et al. 1998 ). Likewise, our results confirmed the stability of the L1 protein in prokaryotic and eukaryotic hosts and showed it was a thermostable and hydrophilic protein and two selected software did not predict any signal peptide for the L1 protein. Li in 1998 determined the vital role of disulfide bonds in papillomavirus capsid assembly and suggested a conserved position (cysteine (C) 424) and also the mutant one could not assemble in vitro into capsid-like structures (Li et al. 1998 ). In 1998, Martin Sapp et al. described the conserved exclusive disulfide bond of C176 with C427 and confirmed the importance of this bond in the structure of the papillomavirus capsid and DNA packaging (Sapp et al. 1998 ). Conway et al. in 2011 indicated that capsids properly mature and become stabilized over time (10-day to 20-day) (Conway et al. ). In addition, several individual L1 cysteine residues (428, 185, and 175) have an indispensable role in this process. Similar to previous studies, the present results confirmed several cysteine bonds and found some bonds which were omitted as a result of some mutations occurred in such regions in different selected sequences. However, eight cysteine residues (13, 211, 371, 128, 405, 172, 255 and 371) were found in all selected and reference sequences which have the critical role in the structure of the L1 protein. The difference between the positions found in the present study and previous studies may be related to the different geographical regions, the different methods which used, and the different HPV genotypes. It can be suggested that the disulfide bond 13-13,as well as all other disulfide bonds, may play critical roles in constructing L1 pentamer which has interaction with L2 proteins. Although our results showed 4 conserved glycosylation positions in all sequences, based on previous studies we could not find any function for glycosylated regions. Zhou et al. in 1993 showed while the majority of L1 protein localizes in the cell nucleus, glycosylated L1 remains in the endoplasmic reticulum and it is not exported from the cell nor translocated to the cell membrane or the cell nucleus (Zhou et al. 1993 ). Therefore, they concluded that glycosylated L1 was not important in the construction of papillomavirus virion. While it was suggested that N-linked glycosylation played a significant role in VLP binding cells, in 1999 Joyce et al. could not find any concentration dose of tunicamycin to measure the effect of glycosylation (Joyce et al. 1999) . Two conserved amino acid residues (Threonine(T) 340 and T129) were introduced as phosphorylation positions on L1 protein by Buck et al. in 2013 (Buck et al. 2013 ). In addition, Xi et al. showed that L1 was phosphorylated only in the earlier weeks and this seems to be fairly unstable (Xi and Banks 1991) . The present results indicated several phosphorylation sites that were conserved among all selected sequences and the reference sequence. Bioinformatics has provided a great context to define the structure of virus proteins (Behzad et al. 2019; Dehghani et al. 2017 Dehghani et al. , 2019a Moattari et al. 2015; Zahra et al. 2020) . Bioinformatics tools were employed widely in this research to determine the structure of L1 and illustrated that the majority of the secondary structures belonged to the random coil and extended strand. The tertiary structure prediction by using three reliable software showed (PS)2-v2 could be able to construct a more reliable structure with the highest percentage of amino acid residues in the favored region in Ramachandran plot and the highest Qmean score. However, the predicted structure by (PS)2-v2 just covered 85% of the L1 sequence and by considering this fact we can conclude that I-TASSER predicted a more useful structure while it covered the complete sequence. In addition, Phyre 2 showed the lower accuracy in comparison with (PS)2-v2 and its' coverage was around 85%. In our study, by using "discovery studio", the structural analysis showed five distinct loops on L1 protein and the related sequences for each loop. Some mutations were located in these loops, i.e. T202 N and N207 T in EF, T292A in FG, and T379P in H1, which may have a significant effect on the loops structure and the antibodies response. Carter et al. showed multiple neutralizing epitopes on the HPV virion surface, three loops of which (DE, FG, and HI) were most important and vital for binding by neutralizing antibodies (Carter et al. 2006) . The crystal structures of four L1 pentamers were determined by Bishop et al. in 2007 (Bishop et al. 2007 . They determined that the surface loops contain the known epitopes for neutralizing monoclonal antibodies (NmAbs). Bissett et al. found the role of the DE loop and the late region of the FG loop as neutralizing antibodies epitopes while their findings were critical to describe vaccineinduced cross-neutralizing antibodies which can play a vital role in vaccine-induced cross-protection (Bissett et al. 2012 ). Roth et al. by using monoclonal antibodies (mAbs) Fig. 2 Coverage of the predicted tertiary structures by 3 reliable software. The coverage of "I-TASSER" was 100%; it was 90% for "Phyre 2" , and around 85% for "(PS)2-v2" recognized three loops (BC, DE, and FG) as the main epitopes and showed the majority of neutralizing epitopes are located on the tip of the capsomere that were related to these variable loops; this data is crucial for designing intertypic HPV vaccines (Roth et al. 2006) . Comparison between the identified loops, hydrophilicity, flexibility, and surface accessibility revealed that BC, DE, and H1 loops were located in regions with the highest surface accessibility scores. In addition, the highest flexibility scores belonged to DE, EF, and FG loops and 4 loops (BC, DE, FG, and H1) had high hydrophilicity scores. Based on three indices (hydrophilicity, flexibility, and surface accessibility), we found BC, DE, FG, and H1 were the highest antigenic ability among other regions of the L1 protein. The combination of B cell epitope prediction and antigenic properties clarified that BC and FG were the most capable regions of L1 with the highest ability to apply in recombinant vaccines against HPV infections (Table 9 ). By considering the list of mutations, the BC loop was completely conserved and just one codon (292) in FG loop showed substitution among selected sequences. Humoral immunity has always the key role to control virus infections (Ajorloo et al. 2015; Alborzi et al. 2017; Negahdaripour et al. 2017) . B cell epitope prediction showed 6 potential regions, 3 of which were completely consistent with 3 loops: BC, EF, and FG that illustrated the importance of these loops in humoral response against HPV L1 protein. In comparison with previous studies, it can be concluded that FG is the most reported loop which has tremendous potential as an epitope for neutralizing antibodies. It seemed that the mentioned mutations did not have a significant effect on B cell epitope regions and their scores; however, substitution in amino acid 202 had a positive impact on EF loop and increased the score of this region as a potential B cell epitope. To conclude, the present results suggested two regions in HPV L1 protein with high potential to induce immune system which can be used as new candidates for new HPV vaccines. In spite of several substitutions in the studied sequences, they had an insignificant impact on all features of L1 protein; it can be inferred that L1 has a conserved sequence and the existing vaccines are still functional for Iranian patients. However, monitoring mutations is crucial to examine the effectiveness and efficacy of the vaccines. Structural and immune properties of Iranian L1 showed the vital role of 5 loops which lead to a strong immune response against HPV infections. Detection of specific antibodies to HCV-ARF/ CORE+ 1 protein in cirrhotic and non-cirrhotic patients with hepatitis C: a possible association with progressive fibrosis archives of Role of serum level and genetic variation of IL-28B in interferon responsiveness and advanced liver disease in chronic hepatitis C patients Expression and characterization of HPV-16 L1 capsid protein in Pichia pastoris Interaction of human herpesvirus 8 viral interleukin-6 with human interleukin-6 receptor using in silico approach: the potential role in HHV-8 pathogenesis QMEAN: a comprehensive scoring function for model quality assessment Crystal structures of four types of Human Papillomavirus L1 Capsid Proteins understanding the specificity of neutralizing monoclonal antibodies Human papillomavirus antibody reference reagents for use in postvaccination surveillance serology Sequence and structure-based prediction of eukaryotic protein phosphorylation sites Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence CO/IARC information Centre on HPVand Cancer (HPV information Centre). Human Papillomavirus and Related Diseases in the World Summary Report Buck CB, Day PM Identification of human papillomavirus type 16 L1 surface loops required for neutralization by human sera In silico platform for prediction of N-, O-and C-glycosites in eukaryotic protein sequences 2: protein structure prediction server SCRATCH: a protein structure and structural feature prediction server Differentiation-dependent interpentameric disulfide bond stabilizes native human papillomavirus type 16 Epidemiology of human papillomavirus (HPV) infection among Iranian women identified with cervical infections: a systematic review and meta-analysis of national data Functional and structural characterization of Ebola virus glycoprotein (1976-2015)-an in silico study Using immunoinformatics and structural approaches to design a novel HHV8 vaccine Bioinformatics analysis of domain 1 of HCV-Core protein: Iran international journal of peptide research and therapeutics:1-18 Human papillomavirus molecular biology and disease association VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines High risk HPV types in southern Iranian patients with cervical cancer DiANNA 1.1: an extension of the DiANNA web server for ternary cysteine classification Human papillomavirus type 16 pseudovirions with few point mutations in L1 major capsid protein FG loop could escape actual or future vaccination for potential use in gene therapy High-performance signal peptide prediction based on sequence alignment techniques Interaction of human papillomaviruses with the host immune system: a well evolved relationship virology Protein identification and analysis tools on the ExPASy server. The proteomics protocols handbook SOPMA: significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments Prevalence of human papillomavirus genotypes in women with normal and abnormal cervical cytology in Iran Human papillomavirus and breast cancer in Iran: a meta-analysis GalaxyRefine: protein structure refinement driven by side-chain repacking PrediSi: prediction of signal peptides and their cleavage positions The importance of intrinsic disorder for protein phosphorylation Mutational analysis of human papillomavirus type 16 major capsid protein L1: the cysteines affecting the intermolecular bonding and structure of L1-capsids Prevalence of human papillomavirus infection among Iranian women using COBAS HPV DNA testing BepiPred-2.0: improving sequence-based B-cell epitope prediction using conformational epitopes The L1 major capsid protein of human papillomavirus type 11 recombinant virus-like particles interacts with heparin and cell-surface glycosaminoglycans on human keratinocytes Protein structure prediction on the Web: a case study using the Phyre server Human papilloma virus in oral cancer Expression of HPV L1 capsid protein in cervical specimens with HPV infection Intercapsomeric disulfide bonds in papillomavirus assembly and disassembly Silico functional and structural characterization of H1N1 influenza a viruses Hemagglutinin Cervarix™: a vaccine for the prevention of HPV 16, 18-associated cervical cancer The prevalence of human papillomavirus in cervical cancer in Iran A novel HPV prophylactic peptide vaccine, designed by immunoinformatics and structural vaccinology approaches infection Human papillomavirus type 16 variant analysis of E6, E7, and L1 genes and long control region in biopsy samples from cervical cancer patients in North India ctla-4 gene variations may influence cervical cancer susceptibility Characterization of neutralizing epitopes within the major capsid protein of human papillomavirus type 33 I-TASSER: a unified platform for automated protein structure and function prediction BcePred: prediction of continuous B-cell epitopes in antigenic sequences using physico-chemical properties AlgPred: prediction of allergenic proteins and mapping of IgE epitopes Prediction of continuous B-cell epitopes in an antigen using recurrent neural network proteins Distribution of human papillomavirus genotypes in Iranian women according to the severity of the cervical lesion Papillomavirus assembly requires trimerization of the major capsid protein by disulfides between two highly conserved cysteines Human papilloma virus prevalence and types among Iranian women attending regular gynecological visits Reports of Prevalence of HPV infection and high risk HPV genotypes (16, 18), among monogamous and polygamous women Human papillomavirus quadrivalent (types 6, 11, 16, 18) recombinant vaccine (Gardasil®) The L1 major capsid protein of human papillomavirus type 16 variants affects yield of virus-like particles produced in an insect cell expression system Human papillomavirus type 16 and 18 L1 protein peptide binding to VERO and HeLa cells inhibits their VLPs binding The immune epitope database (IEDB): 2018 update Baculovirus expression of the human papillomavirus type 16 capsid proteins: detection of L1-L2 protein complexes Epitomics: IgG-epitome decoding of E6, E7 and L1 proteins from oncogenic human papillomavirus type 58 Human papillomavirus (HPV) prevalence and types among women attending regular gynecological visit in Tehran Antibody development to HCV alternate reading frame protein in liver transplant candidate and its computational analysis Expression of human papillomavirus type 16 L1 protein in Escherichia coli: denaturation, renaturation, and self-assembly of virus-like particlesin vitro Papillomavirus genome structure, expression, and post-transcriptional regulation Frontiers in bioscience Glycosylation of human papillomavirus type 16 L1 protein Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Acknowledgements The authors would like to thank Shiraz University of Medical Sciences, Shiraz, Iran for support (Grant number: 97-01-59-18192), and also Center for Development of Clinical Research of Nemazee Hospital and Dr. Nasrin Shokrpour for editorial assistance and Jean Krugman, retired professor of General Studies, DeVry University, Decatur, Georgia, USA for help in revising the English of the manuscript. Special thanks to Nooshin Zare, Anahid khanoomi Hasanshahi and Rozha Rezaei for their kind assistance and also Dr. Marjan Zare and Dr. Behnam Honarvar for statistical and design consultation, respectively. Conflict of interest Author Behzad Dehghani declares that he has no conflict of interest. Author Zahra Hasanshahi declares that she has no conflict of interest. Author Tayebeh Hashempour declares that she has no conflict of interest. Author Mohamad Motamedifar declares that he has no conflict of interest.Ethical approval This article does not contain any studies with human participants or animals performed by any of the authors.