Protein intrinsic disorder and human papillomaviruses: increased amount of disorder in E6 and E7 oncoproteins from high risk HPVs

It is recognized now that many functional proteins or their long segments are devoid of stable secondary and/or tertiary structure and exist instead as very dynamic ensembles of conformations. They are known by different names including natively unfolded, intrinsically disordered, intrinsically unstructured, rheomorphic, pliable, and different combinations thereof. Many important functions and activities have been associated with these intrinsically disordered proteins (IDPs), including molecular recognition, signaling, and regulation. It is also believed that disorder of these proteins allows function to be readily modified through phosphorylation, acetylation, ubiquitination, hydroxylation, and proteolysis. Bioinformatics analysis revealed that IDPs comprise a large fraction of different proteomes. Furthermore, it is established that the intrinsic disorder is relatively abundant among cancer-related and other disease-related proteins and IDPs play a number of key roles in oncogenesis. There are more than 100 different types of human papillomaviruses (HPVs), which are the causative agents of benign papillomas/warts, and cofactors in the development of carcinomas of the genital tract, head and neck, and epidermis. With respect to their association with cancer, HPVs are grouped into two classes, known as low (e.g., HPV-6 and HPV-11) and high-risk (e.g., HPV-16 and HPV-18) types. The entire proteome of HPV includes six nonstructural proteins [E1, E2, E4, E5, E6, and E7 (the latter two are known to function as oncoproteins in the high-risk HPVs)] and two structural proteins (L1 and L2). To understand whether intrinsic disorder plays a role in the oncogenic potential of different HPV types, we have performed a detailed bioinformatics analysis of proteomes of high-risk and low-risk HPVs with the major focus on E6 and E7 oncoproteins. The results of this analysis are consistent with the conclusion that high-risk HPVs are characterized by the increased amount of intrinsic disorder in transforming proteins E6 and E7.     

A Mathematical Model of Cell Cycle Dysregulation Due to Human Papillomavirus Infection

Human papillomaviruses (HPVs) that infect mucosal epithelium can be classified as high risk or low risk based on their propensity to cause lesions that can undergo malignant progression. HPVs produce the E7 protein that binds to cell cycle regulatory proteins including the retinoblastoma tumor suppressor protein (RB) to modulate cell cycle control. Generally, high-risk HPV E7 proteins bind to RB with a higher affinity than low-risk HPV E7s, but both are able to deactivate RB and trigger S phase progression. In uninfected cells, RB inactivation is a tightly controlled process that must coincide with growth factor stimulation to commit cells to division. High-risk HPV E7 proteins short-circuit this control by decreasing growth factor requirement for cell division. We develop a mathematical model to examine the role that RB binding affinity, growth factor concentration, and E7 concentration have on cell cycle progression. Our model predicts that high RB binding affinity and E7 concentration accelerate the \(\mathrm {G_{1}}\) to S phase transition and weaken the dependence on growth factor. This model thus captures a key step in high-risk HPV oncogenesis.     

The E6 Oncoproteins of High-Risk Papillomaviruses Bind to a Novel Putative GAP Protein, E6TP1, and Target It for Degradation

The high-risk human papillomaviruses (HPVs) are associated with carcinomas of the cervix and other genital tumors. Previous studies have identified two viral oncoproteins, E6 and E7, which are expressed in the majority of HPV-associated carcinomas. The ability of high-risk HPV E6 protein to immortalize human mammary epithelial cells (MECs) has provided a single-gene model to study the mechanisms of E6-induced oncogenic transformation. In this system, the E6 protein targets the p53 tumor suppressor protein for degradation, and mutational analyses have shown that E6-induced degradation of p53 protein is required for MEC immortalization. However, the inability of most dominant-negative p53 mutants to induce efficient immortalization of MECs suggests the existence of additional targets of the HPV E6 oncoprotein. Using the yeast two-hybrid system, we have isolated a novel E6-binding protein. This polypeptide, designated E6TP1 (E6-targeted protein 1), exhibits high homology to GTPase-activating proteins for Rap, including SPA-1, tuberin, and Rap1GAP. The mRNA for E6TP1 is widely expressed in tissues and in vitro-cultured cell lines. The gene for E6TP1 localizes to chromosome 14q23.2-14q24.3 within a locus that has been shown to undergo loss of heterozygosity in malignant meningiomas. Importantly, E6TP1 is targeted for degradation by the high-risk but not the low-risk HPV E6 proteins both in vitro and in vivo. Furthermore, the immortalization-competent but not the immortalization-incompetent HPV16 E6 mutants target the E6TP1 protein for degradation. Our results identify a novel target for the E6 oncoprotein and provide a potential link between HPV E6 oncogenesis and alteration of a small G protein signaling pathway.     

Genetic variations in the DNA replication origins of human papillomavirus family correlate with their oncogenic potential

Human papillomaviruses (HPVs) encompass a large family of viruses that range from benign to highly carcinogenic. The crucial differences between benign and carcinogenic types of HPV remain unknown, except that the two HPV types differ in the frequency of DNA replication. We have systematically analyzed the mechanism of HPV DNA replication initiation in low-risk and high-risk HPVs. Our results demonstrate that HPV-encoded E2 initiator protein and its four binding sites in the replication origin play pivotal roles in determining the destiny of the HPV-infected cell. We have identified strain-specific single nucleotide variations in E2 binding sites found only in the high-risk HPVs. We have demonstrated that these variations result in attenuated formation of the E2-DNA complex. E2 binding to these sites is linked to the activation of the DNA replication origin as well as initiation of DNA replication. Both electrophoretic mobility shift assay and atomic force microscopy studies demonstrated that binding of E2 from either low- or high-risk HPVs with variant binding sequences lacked multimeric E2-DNA complex formation in vitro. These results provided a molecular basis of differential DNA replication in the two types of HPVs and pointed to a correlation with the development of cancer.     

Comparison of the structure and DNA-binding properties of the E2 proteins from an oncogenic and a non-oncogenic human papillomavirus

Human papillomaviruses (HPVS) that infect the genital tract can be divided into two groups: high-risk HPV types, such as HPV 16 and HPV 18, are associated with cancer, low-risk HPV types, such as HPV 6, are associated with benign warts. In both high-risk and low-risk HPV types, the papillomavirus E2 protein binds to four sites within the viral long control region (LCR) and regulates viral gene expression. Here, we present the crystal structure of the minimal DNA-binding domain (DBD) from the HPV 6 E2 protein. We show that the HPV 6 E2 DBD is structurally more similar to the HPV 18 and bovine papillomavirus type 1 (BPV1) E2 proteins than it is to the HPV 16 E2 protein. Using gel retardation assays, we show that the hierarchy of E2 sites within the HPV 16 and HPV 6 LCRs are different. However, despite these differences in structure and site preference, both the HPV 16 and 6 E2 DBDs recognise an extended version of the consensus E2 binding site derived from studies of the BPV1 E2 protein. In both cases, the preferred binding site is 5'AACCGN(4)CGGTT3', where the additional flanking base-pairs are in bold and N(4) represents a four base-pair central spacer. Both of these HPV proteins bind preferentially to E2 sites that contain an A:T-rich central spacer. We show that the preference for an A:T-rich central spacer is due, at least in part, to the need to adopt a DNA conformation that facilitates protein contacts with the flanking base-pairs.     

Conditionally activated E7 proteins of high-risk and low-risk human papillomaviruses induce S phase in postmitotic, differentiated human keratinocytes

The productive program of human papillomaviruses (HPVs) in epithelia is tightly linked to squamous differentiation. The E7 proteins of high-risk HPV genotypes efficiently inactivate the pRB family of proteins that control the cell cycle, triggering S phase in suprabasal keratinocytes. This ability has until now not been demonstrated for the low-risk HPV-6 or HPV-11 E7 proteins. An inducible system in which HPV-16 E7 is fused to the ligand binding domain of the human estrogen receptor (ER) was described by Smith-McCune et al. (K. Smith-McCune, D. Kalman, C. Robbins, S. Shivakumar, L. Yuschenkoff, and J. M. Bishop, Proc. Natl. Acad. Sci. USA 96:6999-7004, 1999). In the absence of hormone, E7ER is cytoplasmic, and upon addition of 17beta-estradiol, it translocates to the nucleus. Using organotypic epithelial raft cultures developed from primary human keratinocytes, we show that 16E7ER promotes either S-phase reentry or p21cip1 accumulation in differentiated keratinocytes in a stochastic manner as early as 6 h postinduction with 17beta-estradiol. A vector expressing the ER moiety alone had no effect. These observations prove unequivocally that the E7 protein drives S-phase reentry in postmitotic, differentiated keratinocytes rather than preventing S-phase exit while the cells ascend through the epithelium. HPV-11 E7ER and, much less efficiently, HPV-6 E7ER also promoted S-phase reentry by differentiated cells upon exposure to 17beta-estradiol. S-phase induction required the consensus pRB binding motif. We propose that the elevated nuclear levels of the low-risk HPV E7 protein afforded by the inducible system account for the positive results. These observations are entirely consistent with the fact that low-risk HPV genotypes replicate in the differentiated strata in patient specimens, as do the high-risk HPVs.     

HPVomics: An integrated resource for the human papillomavirus epitome and therapeutics

Human papillomaviruses (HPVs) belongs to the Papillomaviridae family, which is divided into high-risk (HR), and low-risk (LR) HPVs based on their disease-causing competence. HR-HPVs 16 and 18 are known to cause distinct carcinomas like cervical and head and neck, whereas LR-HPVs are commonly associated with the genital warts. We have developed an integrative platform; HPVomics dedicated to the potential therapeutic regimens targeting all HPV genes including oncoproteins E6, E7 and E5. We primarily focused on eighteen HR-HPVs and eleven LR-HPVs. It mainly deals with therapeutically imperative elements, i.e., vaccine epitopes, siRNAs, sgRNAs, and anti-viral peptides. Simultaneously, it also comprises of genome browser, whole-genome sequences and annotation of HPVs with searching and filtering capabilities. Moreover, we have also developed an integrated support vector machine (SVM) based computational algorithm “HPVepi” for the prediction of HPV epitome. We hope that HPVomics (http://bioinfo.imtech.res.in/manojk/hpvomics/) will assist the scientific community engaged in HPV research.     

Role of the PDZ domain-binding motif of the oncoprotein E6 in the pathogenesis of human papillomavirus type 31

A number of PDZ domain-containing proteins have been identified as binding partners for the oncoprotein E6 of the high-risk type human papillomaviruses (HPVs). These include hDlg, hScrib, MAGI-1, MAGI-2, MAGI-3, and MUPP1. The PDZ domain-binding motif (-X-T-X-V) at the carboxy terminus of E6 is essential for targeting PDZ proteins for proteasomal degradation. The presence of this motif only in the high-risk HPVs suggests its possible role in HPV-induced oncogenesis. To investigate the role of the PDZ domain-binding motif of E6 in the HPV life cycle, two mutant HPV31 genomes were constructed: E6ValDelta, with a deletion of the last amino acid residue of E6 (valine), and E6ETQVDelta, with a deletion of the entire PDZ domain-binding motif of E6 (ETQV). Three human foreskin keratinocyte (HFK) cell lines were established which maintained transfected wild-type HPV31 or either of two mutant genomes. Cells containing either of two mutant genomes were significantly retarded in their growth rates and reduced in their viral copy numbers compared to those transfected with wild-type genomes. Western analysis did not reveal any significant changes in the levels of PDZ proteins following stable transfection of any HPV31 genomes into HFKs. Although the E6ETQVDelta-transfected HFKs exhibited a pattern of morphological differentiation that appeared different from the HPV31 wild-type-transfected HFKs in organotypic raft cultures, immunohistochemical analysis failed to identify substantial changes in the differentiation-dependent membrane localization of hDlg proteins. These results suggest that binding of E6 to PDZ proteins modulates the early viral functions such as proliferation and maintenance of the viral copy number in undifferentiated cells.     

E6AP-dependent degradation of DLG4/PSD95 by high-risk human papillomavirus type 18 E6 protein

In most cervical cancers, DNAs of high-risk mucosotropic human papillomaviruses (HPVs), such as types 16 and 18, are maintained so as to express two viral proteins, E6 and E7, suggesting that they play important roles in carcinogenesis. The carboxy-terminal PDZ domain-binding motif of the E6 proteins is in fact essential for transformation of rodent cells and induction of hyperplasia in E6-transgenic mouse skin. To date, seven PDZ domain-containing proteins, including DLG1/hDLG, which is a human homologue of the Drosophila discs large tumor suppressor (Dlg), have been identified as targets of high-risk HPV E6 proteins. Here, we describe DLG4/PSD95, another human homologue of Dlg, as a novel E6 target. DLG4 was found to be expressed in normal human cells, including cervical keratinocytes, but only to a limited extent in both HPV-positive and HPV-negative cervical cancer cell lines. Expression of HPV18 E6 in HCK1T decreased DLG4 levels more strongly than did HPV16 E6, the carboxy-terminal motif of the proteins being critical for binding and degradation of DLG4 in vitro. DLG4 levels were restored by expression of either E6AP-specific short hairpin RNA or bovine papillomavirus type 1 E2 in HeLa but not CaSki or SiHa cells, reflecting downregulation of DLG4 mRNA as opposed to protein by an HPV-independent mechanism in HPV16-positive cancer lines. The tumorigenicity of CaSki cells was strongly inhibited by forced expression of DLG4, while growth in culture was not inhibited at all. These results suggest that DLG4 may function as a tumor suppressor in the development of HPV-associated cancers.     

Kinetic analysis of the interactions of human papillomavirus E6 oncoproteins with the ubiquitin ligase E6AP using surface plasmon resonance

Cervical cancers evolve from lesions generated by genital human papillomaviruses (HPV). "Low-risk" genital HPVs cause benign proliferations whereas "high-risk" types have the potential to progress into cancer. High-risk HPV E6 oncoproteins interact with the ubiquitin ligase E6AP and target several cellular proteins, including p53 and proteins of the MAGI family, towards ubiquitin-mediated degradation. E6AP, like other E6 binding proteins such as E6BP, IRF-3 and paxillin, interacts with E6 via a consensus leucine-charged motif. Here we have investigated the kinetics of the interactions of a 15-mer peptide containing the LxxvarphiLsh motif of E6AP with E6. For this we have developed a Biacore assay based on antibody-capture on the sensor surface of GST- and/or MBP-E6AP peptide constructs followed by E6 protein injection. Our experiments show that E6 oncoproteins from four major high-risk (16, 18, 33 and 58) HPV types bind to E6AP with equilibrium dissociation constants in the low micromolar range. The kinetic dissociation parameters of these interactions are remarkably similar. On the other hand, low-risk HPV 11 E6 does not interact with E6AP even at relatively high concentrations. We also show that the two zinc-binding domains of E6 are required for E6AP recognition. Finally, we have analysed the binding properties of site-directed mutants of the E6AP-derived peptide. We demonstrate the importance for binding of conserved aliphatic side-chains and the moderate role of the global negative charge of the peptide. This work provides the first quantitative data on an HPV E6-mediated interaction, which support the current models of E6AP-mediated degradation.     

Nonconserved lysine residues attenuate the biological function of the low-risk human papillomavirus E7 protein

The mucosotrophic human papillomaviruses (HPVs) are classified as high-risk (HR) or low-risk (LR) genotypes based on their neoplastic properties. We have demonstrated previously that the E7 protein destabilizes p130, a pRb-related pocket protein, thereby promoting S-phase reentry in postmitotic, differentiated keratinocytes of squamous epithelia, and that HR HPV E7 does so more efficiently than LR HPV E7. The E7 proteins of LR HPV-11 and -6b uniquely possess lysine residues following a casein kinase II phosphorylation motif which is critical for the biological function of E7. We now show that mutations of these lysine residues elevated the efficiency of S-phase reentry, independent of their charge. An 11E7 K39,42R mutation moderately increased the association with and the destabilization of p130. Unexpectedly, polyubiquitination on these lysine residues did not attenuate E7 activity, as their mutation caused elevated proteasomal degradation and decreased protein stability. In this regard, the biologically more potent HR HPV E7 proteins were also less stable than the LR HPV E7 proteins. We infer that these lysine residues impede functional protein-protein interactions. A G22D mutation of 11E7 at the pocket protein binding motif possessed augmented efficiency in promoting S-phase reentry and strongly enhanced association with p130 and pRb. The combined effects of these two classes of 11E7 mutations exhibited an efficiency of S-phase reentry comparable to that of HR HPV E7. Thus, these nonconserved residues are primarily responsible for the differential abilities of LR and HR HPV E7 proteins to promote unscheduled DNA replication in organotypic raft cultures.     

Sterile α Motif Domain Containing 9 Is a Novel Cellular Interacting Partner to Low-Risk Type Human Papillomavirus E6 Proteins

Low-risk type human papillomavirus (HPV) 6 and 11 infection causes recurrent respiratory papillomatosis (RRP) and genital warts. RRP is the most common benign tumor of the larynx in children with frequent relapses. Repeated surgeries are often needed to improve vocal function and prevent life-threatening respiratory obstruction. Currently, there are no effective treatments available to completely eliminate these diseases, largely due to limited knowledge regarding their viral molecular pathogenesis. HPV E6 proteins contribute to cell immortalization by interacting with a variety of cellular proteins, which have been well studied for the high-risk type HPVs related to cancer progression. However, the functions of low-risk HPV E6 proteins are largely unknown. In this study, we report GST-pulldown coupled mass spectrometry analysis with low-risk HPV E6 proteins that identified sterile alpha motif domain containing 9 (SAMD9) as a novel interacting partner. We then confirmed the interaction between HPV-E6 and SAMD9 using co-immunoprecipitation, proximity ligation assay, and confocal immunofluorescence staining. The SAMD9 gene is down-regulated in a variety of neoplasms and deleteriously mutated in normophosphatemic familial tumoral calcinosis. Interestingly, SAMD9 also has antiviral functions against poxvirus. Our study adds to the limited knowledge of the molecular properties of low-risk HPVs and describes new potential functions for the low-risk HPV E6 protein.