ABSTRACT
Several serine and threonine residues of the papillomavirus early E2 protein have been found to be phosphorylated. In contrast, only one E2 tyrosine phosphorylation site in BPV-1 (tyrosine 102) and one in HPV-16/31 (tyrosine 138) site have been characterized. Between BPV-1 and HPV-31 E2, 8 of the 11 tyrosines are conserved in the N-terminal domain, suggesting that phosphorylation of tyrosines has an essential role in E2 biology. In this study, we examine the effect of Y102 phosphorylation on HPV-31 E2 biology. Y102 proteins mutated either to the potential phospho-mimetic glutamic acid (Y102E) or to the nonphosphorylated homologue phenylalanine (Y102F) remain nuclear; however, Y102E is more associated with the nuclear matrix fraction. This is consistent with the inability of Y102E to bind TopBP1. Both BPV-1 and HPV-31 Y102E are similar in that neither binds the C terminus of Brd4, but in all other aspects the mutant behaves differently between the two families of papillomaviruses. BPV-1 Y102E was unable to bind E1 and did not replicate in a transient in vitro assay, while HPV-31 Y102E binds E1 and was able to replicate, albeit at lower levels than wild type. To examine the effect of E2 mutations under more native-like infection conditions, a neomycin-selectable marker was inserted into L1/L2 of the HPV-31 genome, creating HPV-31neo. This genome was maintained in every cell line tested for at least 50 days posttransfection/infection. Y102E in both transfection and infection conditions was unable to maintain high episome copy numbers in epithelial cell lines.
IMPORTANCE Posttranslational modifications by phosphorylation can change protein activities, binding partners, or localization. Tyrosine 102 is conserved between delta papillomavirus BPV-1 and alpha papillomavirus HPV-31 E2. We characterized mutations of HPV-31 E2 for interactions with relevant cellular binding partners and replication in the context of the viral genome.
INTRODUCTION
Human papillomavirus 31 (HPV-31) is a small double-stranded DNA virus that can infect mucosal epithelium and, if persistent, can lead to cervical and oropharyngeal cancer. Regulation of both viral transcription as well as replication relies on the viral early protein E2. The N-terminal 200 amino acids of E2 contain the transcriptional activation domain (TAD), which regulates a multitude of functions during the viral life cycle. More than 70 proteins have been reported to interact with E2 TAD (1, 2). The C-terminal 100 amino acids form a high-affinity, sequence-specific DNA-binding domain (DBD). Separating these two regions is the nonconserved hinge domain.
Posttranslational phosphorylations often alter protein stability, induce a conformational change, cause subcellular redistribution, or modify the repertoire of protein-protein interactions. Several phosphorylation sites have been discovered in E2, the first of which were S298 and S301 in BPV-1 in the DBD (3). Preventing phosphorylation at these sites through mutation to alanine increased E2 protein levels and corresponding transient replication and transcription activities (4), suggesting that phosphorylation at these sites is inhibitory.
The absence of any prior identification of tyrosine (Y) phosphorylation in E2 could be attributed to the method of purification of E2 from insect or mouse cells (3, 5, 6). Our lab recently identified the first phosphorylated tyrosine residue, Y102, through mass spectroscopy of E2 purified from human cervical cells (7). Mutation of this tyrosine to the phospho-mimetic glutamic acid (Y102E) inhibited both the transcription and replication activities of E2. The homologous site in HPV-31 E2 is also tyrosine (Y102); however, given that HPV-31 and BPV-1 are only moderately identical, have different tissue tropisms, and are under different transcriptional control (8), it is possible the role of phosphorylation at Y102 differs between HPV-31 and BPV-1. In this study, HPV-31 E2 tyrosine 102 was mutated to glutamic acid (Y102E) as well as to phenylalanine (Y102F) to mimic the phospho-deficient form, and the mutants were examined for effects on replication, transcription, localization, and genome maintenance.
RESULTS
Identification of tyrosine phosphorylation in HPV-31 E2.We discovered the first tyrosine phosphorylation site of the papillomavirus E2 protein at residue 102 through mass spectrometry analysis of the BPV-1 E2 protein (7). Y102 is located in the TAD domain (Fig. 1A) and is conserved between BPV-1, HPV-16, and HPV-31 E2 (Fig. 1B). To determine if HPV-31 E2 can be phosphorylated at Y102, we performed mass spectrometry analysis of HPV-31 E2 overexpressed in HEK 293TT cells. The peptides collected covered 95% of the entire E2 protein (Fig. 1C) and a peptide with phosphorylated Y102 was captured multiple times, along with several additional phosphorylated residues (Fig. 1D). The protein sequence of HPV-31 E2 is only 28% identical to BPV-1, but, despite this low homology, most (8 of 11) of the tyrosines in the N-terminal TAD domain are conserved, including Y102. In contrast, only about 30% of the BPV-1 serine and threonine residues are similarly positioned in the HPV-31 TAD. This preference of tyrosine conservation suggests phosphorylation at these sites is essential for E2 function and regulation. Three computer algorithms also predict that Y102 is among the three most likely phosphorylated residues in HPV-31 (NetPhos 3.1 (9), PhosphoSVM (10), and Viral Phos [http://csb.cse.yzu.edu.tw/ViralPhos/index.html]).
Identification of phosphorylated Y102. (A) Schematic of HPV-31 E2 protein representing the transactivation domain (TAD) and DNA-binding domain (DBD). (B) Conservation of the region surrounding Y102. (C) HPV-31 E2 was purified, digested with trypsin, and submitted for mass spectrometric analysis. Peptide coverage was 95% as highlighted in gray. (D) Phosphorylated amino acids of HPV-31 E2 with high PhosphoRS probabilities (above 75%) are listed on the left. The trypsin fragments carrying the modifications are listed in the middle, with phosphorylation sites in black. HPV-16 and BPV-1 homologue residues are listed on the right, with conserved residues in black and nonconserved in gray. (E) FLAG-HPV-31 E2 (WT or Y102 mutants) was expressed in HEK 293TT cells, immunoprecipitated with custom phospho-specific Y102 antibody, and blotted with FLAG antibody. (F) Y102 mutants localize to the nucleus. N/TERT cells were transfected with FLAG-HPV-31 E2 (WT, Y102E, and Y102F) and HA-HPV-31 E1, and immunostained with FLAG and HA antibodies.
Polyclonal antibodies to phosphorylated HPV-31 E2 Y102 were generated in rabbits as a tool to examine Y102 phosphorylation. We created point mutants in HPV-31 E2 to mimic the phosphorylated state (Y102E) as well as the nonphosphorylated state (Y102F) of E2, and tested their ability to be recognized by the phospho-102 antibody. HEK 293TT cells transfected with the wild-type (WT) as well as mutant FLAG-HPV-31 E2 proteins were immunoprecipitated with the phospho-102 antibody. As expected, Y102E was the predominant form enriched, with minor binding seen for WT and Y102F (Fig. 1E). A synthesized phosphorylated E2 peptide was used as the antigen, but the phosphate group will get cleaved upon injection into the rabbit. As a result, the resultant polyclonal antibody stock will detect nonphosphorylated as well as phosphorylated protein as, indeed, we see in Fig. 1E. Robust recognition of Y102E by the phospho-102 antibody also suggests that this glutamic acid mutation is a close homologue in structure to the true phosphorylated form of Y102. Unfortunately, this antibody was not effective for E2 localization or coimmunoprecipitation experiments and, while useful for immunoblotting, detection of native E2 in infected cells was not reliable.
Nuclear localization of HPV-31 E2 Y102 mutants.E2 resides in the nucleus, which coincides with its role in genome replication and transcription. Two basic-rich nuclear localization signals (NLS) were identified in BPV-1: N-terminal 107 KRCFKKGAR 115 and C-terminal 339 KCYRFRVKKNHRHR 352 (11). As this N-terminal NLS is close to the Y102 phosphorylation site, and since nuclear shuttling often is modulated by adjacent phosphorylation (12, 13), HPV-31 Y102 mutants might have altered localization. Examination of the HPV-31 E2 sequence, however, reveals conservation of only the C-terminal NLS: 306 KCLRYRLSKYK 316. The N-terminal NLS region contains very few basic amino acids, making it unlikely this region confers localization, and, therefore, phosphorylation of Y102 should not interfere with nuclear positioning. To determine if Y102 affects E2 localization, FLAG-tagged WT, Y102E, and Y102F HPV-31 E2 mutants were cotransfected with HA-HPV-31 E1 into N/TERT cells (human foreskin keratinocytes immortalized with H-tert) (14), fixed with paraformaldehyde, and immunostained with M2 (FLAG) and HA antibodies. Localization of WT E2 as well as both mutants was nuclear (Fig. 1F).
The E2-interacting protein topoisomerase II binding protein 1 (TopBP1) is a large scaffolding protein with nine breast cancer gene C-terminal (BRCT) domains (15) that mediate phospho-dependent protein interactions required for the initiation of human DNA replication, initial viral replication, and segregation of the viral genome during mitosis (16–18). Changing the phosphorylation state of E2 could therefore modulate its ability to bind TopBP1. To determine if binding is altered, FLAG-HPV-31 E2 WT and mutants were transfected into HEK 293TT cells, which were then lysed and subjected to immunoprecipitation (IP) with TopBP1 antibody. Immunoblot analysis of the precipitates shows that WT and Y102F bind endogenous TopBP1, while Y102E binding was greatly reduced (Fig. 2A). Performing the reverse IP with FLAG antibody similarly showed that WT and Y102F, but not Y102E, coimmunoprecipitate TopBP1 (Fig. 2B). Previous studies on HPV-16 E2 demonstrated that in the absence of TopBP1, E2 associates with a nuclear matrix fraction (18). With Y102E unable to bind TopBP1, a similar relocalization would be expected to occur. E2 from transfected cells was sequentially extracted according to an established protocol (Fig. 2C), and Western blot analysis showed that, as predicted, the majority of Y102E associated with the nuclear matrix fraction, much more so than WT E2 (Fig. 2D and E). To control for fractionation, lysates were analyzed for tubulin, histone, and p53: cytoplasmic tubulin was present in the cytoplasmic fraction (lane 1); nuclear histones were released upon DNase and moderate salt treatment (lanes 3 and 4); and p53 was present as previously described in both soluble (lane 1) and nuclear matrix fractions (lane 5) (19) (Fig. 2D).
HPV-31 E2 Y102E does not bind TopBP1. (A and B) Lysates of cells transfected with FLAG-HPV-31 E2 (WT, Y102E, and Y102F) were immunoprecipitated either with an antibody to TopBP1 (A) or FLAG (B). (C) Diagram of lysate fractionation performed in parts D and E. (D) HEK 293TT cells transfected with FLAG-HPV-31 E2 WT or Y102E were lysed and fractionated according to the scheme shown in C, run on a gel, and blotted with FLAG, tubulin, histone, and p53 antibodies. (E) Ratio of nuclear matrix bound (lane 4 in D) over soluble (lane 1) HPV-31 E2 from three experiments, normalized to WT. Y102E was significantly different from WT with P < 0.005, and means are expressed ± SEM.
E2 Y102E transcription and DNA replication activities.Transcription of papillomavirus genes is mainly regulated by E2 and its ability to recruit cellular transcription factors to the viral genome. One such factor is the activation domain modulating factor (Amf1) protein, also known as G-protein pathway suppressor (GPS2). Previously, our lab has shown that the beta-sheet region of BPV-1 E2 binds to GPS2, and mutants that cannot bind GPS2 are defective for transcriptional activation (20). To determine if mutation of Y102 effects binding, soluble lysate from HEK 293TT cells transfected with FLAG-HPV-31 E2 constructs and HA-GPS2 was immunoprecipitated with M2 FLAG antibody. GPS2 bound to wild type and Y102F efficiently, while Y102E had significantly decreased binding (Fig. 3A).
Y102E prevents binding to GPS2 and Brd4-CTM. (A) HEK 293TT cells transfected with FLAG-HPV-31 E2 (WT, Y102E, and Y102F) and HA-GPS2 were lysed and immunoprecipitated with HA antibody. (B) Same as A, except GST-CTM was cotransfected instead of GPS2. (C) Ratio of immunoprecipitated GST-CTM normalized to E2 from three experiments. Values are expressed as means ± SEM. *, P < 0.05; **, P < 0.0001.
Brd4 binding to E2 also regulates the transcriptional activities of E2. Since the BPV-1 Y102E mutant was unable to bind the C-terminal motif (CTM) of Brd4, we sought to check if this occurs with HPV-31 E2. GST-tagged CTM was coexpressed with FLAG-HPV-31 E2 WT or Y102 mutants and immunoprecipitated with FLAG antibody. WT and Y102F bound efficiently (Fig. 3B), while Y102E binding was decreased ∼80% compared to WT, as calculated from three replicates (Fig. 3C).
HPV DNA replication requires E2 binding to the viral E1 DNA helicase. E2 transports E1 to binding palindromes at the origin of replication (ori) (nucleotides [nt] 7721 to 100) (21). To determine if Y102 mutants alter E1 binding, lysates from HA-HPV-31 E1 and FLAG-HPV-31 E2 WT/mutant transfected HEK 293TT cells were immunoprecipitated with FLAG antibody. WT, Y102E, and Y102F all bound E1 to similar levels (Fig. 4A). In contrast, BPV-1 E2 Y102E was unable to bind E1, preventing replication of the BPV-1 genome (7). With HPV-31 E2 Y102E still competent for binding to E1, HPV-31 Y102E might retain replication activities. We utilized a luciferase reporter system in which firefly luciferase is encoded on a plasmid in cis with the HPV-31 origin (22). Cotransfection of this reporter together with E1 and E2 plasmids resulted in increased luciferase expression. Since E2 can also transcriptionally activate gene expression, a control plasmid encoding Renilla luciferase was included. Replication activity can then be calculated as firefly expression divided by Renilla expression. Y102F did not affect transient replication (Fig. 4B), while Y102E retained replication, though at 70% of WT.
Y102E supports transient replication. (A) Y102E/F mutants bind E1. HEK 293TT cells transfected with HA-HPV-31 E1 and FLAG-HPV-31 E2 (WT, Y102E, and Y102F) were lysed and immunoprecipitated with FLAG antibody. (B) In vitro luciferase-based transient replication assay. WT/mutant E2 along with E1, pRenilla, and luciferase plasmid with E2-binding sites were transfected into C33a cells. Following 2 days of transfection, luciferase signals analyzed with DuoGlo kit were normalized against Renilla. Replication levels are an average from 6 replicates. Y102E replicates at significantly lower levels than WT (*, P < 0.005), and means are expressed ± SEM. (C) Expression of E2 and actin loading control for part B.
Upon infection of cells with HPV-31, the genome either is maintained as an extrachromosomal DNA entity called an episome that segregates during mitosis, or integrates into the host chromosomes with a breakpoint often in the E2 region, disrupting E2 regulatory functions. While Y102E and Y102F transiently replicate in C33a cells, the abnormally high levels of transfected E1/E2 in this system may not accurately reflect their ability to initiate DNA replication, maintain the entire genome, and segregate. To examine maintenance of genomes in cells, we created HPV-31 quasiviruses and infected HaCaT and N/TERT cells. Isolation of DNA from these cells at 2, 4, and 7 days postinfection revealed that the cells consistently failed to maintain genomes past day 4 (Fig. 5A). Given the success with which addition of a neo cassette within the HPV-18 genome required cells to retain the genome (23), we sought to create a similar selectable HPV-31 genome in which the L1/L2 region was replaced with the SV40-promoter-intron-Neo-β-actin-polyA cassette, creating HPV-31neo (Fig. 5B). N/TERT cells were transfected with HPV-31neo WT, Y102E, and Y102F constructs, maintained for 50 days in Geneticin, and harvested by Hirt extraction, which specifically isolates small episomal DNA. As shown in Fig. 5B, WT and Y102F replicated the genome and were maintained as episomes, but not Y102E. Genomic DNA isolated from these same cell lines show that all three genomes are present. E1̂E4 RNA was produced in all three cell lines (Fig. 5D). Together these data suggest that the E2 Y102E genome was unable to be maintained episomally and integrated early after infection.
HPV-31 E2 Y102E prevents episomal maintenance. (A) HaCaT and N/TERT cells infected with HPV-31 quasivirus were harvested at days 2, 4, and 7 postinfection. Total DNA was isolated and amplified with primers to the long control region (LCR). (B) Schematic of HPV-31neo genome. The inserted neomycin cassette replaced parts of the L1 and L2 genes in the HPV-31 genome. (C) N/TERT cells transfected with WT, Y102E, or Y102F HPV-31neo genomes were stably maintained through Geneticin selection. At 50 days posttransfection, episomal DNA was isolated through Hirt extraction and the copy number was established through qPCR with LCR primers. (D) RNA was extracted from the cells in B and amplified with primers to E1̂E4. (E) Total DNA was isolated from HaCaT cells at 1 and 2 days postinfection with quasivirus, amplified with primers to LCR, and normalized to the WT genome. (F) Hirt extraction of HaCaT cells infected with quasivirus and stably maintaining genome for 50 days. (G) Same as E except NIKS cells were infected for 14 days. (H) Episomal DNA from 14-day NIKS infection isolated through the exonuclease V method. (I) Controls for exonuclease digestion of part H. Linear genomic actin DNA should be completely digested by exonuclease, while circular mitochondrial DNA is more resistant to digestion. Each graph is representative of multiple infections, and means are expressed ± SEM.
In additional experiments, quasivirions made with HPV-31neo genomes packaged into L1/L2 capsid proteins were infected into NIKS and HaCaT cells to similarly examine episomal maintenance, but in a context with more physiologically accurate copy numbers. Initial infections were performed with pseudoviruses packaging a plasmid encoding mCherry to quantify infection rate. Infection of both NIKS and HaCaTs initially only resulted in 25% red-fluorescing cells. More virus could be used to increase infectivity, but we preferred to use a more normal physiological range of infection with less than 50 viral particles/cell. Since several reports suggest that binding of virus first to the extra cellular matrix (ECM) before the cell surface yields efficient infection (24–26), we tried HaCaT ECM for increasing infectivity. First, HaCaT cells were grown to confluence, the cells but not the deposited ECM were removed with ammonium hydroxide/Triton-X, then virus was allowed to bind overnight. The next day, unbound virus was removed, fresh medium was added along with either HaCaT or NIKS cells. Following 2 days of infection, >95% of cells expressed mCherry and fluoresced red, suggesting that HaCaT ECM drastically enhanced infectivity of these cell lines.
With efficient infection conditions established, we infected HaCaTs with quasivirus encoding WT, Y102E, and Y102F HPV-31 genomes. At days one and two, total DNA was isolated through phenol chloroform extraction and quantified through quantitative PCR (qPCR). All three infections showed similar genome content, though Y102F was present at slightly lower levels (Fig. 5E). At 50 days postinfection, when the cells were well established, DNA was extracted from cell pellets using the Hirt method, and taking into consideration that Y102F had slightly lower initial genomes, we found that again Y102E was unable to maintain high episomal copy number compared to WT and Y102F (Fig. 5F). Infection of NIKS cells similarly resulted in low Y102E episomal copy number (Fig. 5G). Primary HFKs infected with Y102E virus did not result in any proliferating colonies.
Hirt separation of DNA into episomal and cellular genomic fractions has been a widely accepted technique, but some of one fraction frequently contaminates the other. A new technique for distinguishing episomal from integrated HPV involves isolating total DNA with a kit, digesting any linear cellular DNA with an exonuclease, and performing qPCR for the HPV genome, a nuclear gene, and a mitochondrial gene (27). Exonuclease III and V were first compared for DNA digestion. Since exonuclease V digestion was more consistent, more efficient, and does not digest nicked DNA (unlike exonuclease III), all experiments were conducted with exonuclease V. The infected NIKS cells in Fig. 5G were analyzed by this exonuclease protocol for comparison to the Hirt method. HPV genomic qPCR of exonuclease-reacted DNA showed episomal DNA levels consistent with Hirt extraction and that Y102E cells did not contain high episomal copy numbers (Fig. 5H). By subtracting uncut from cut DNA, the integrated DNA content can be estimated. As an internal control for exonuclease, actin DNA levels were monitored to confirm that linear genomic DNA is efficiently digested with exonuclease, and mitochondrial DNA levels should remain relatively unchanged, since the circular mitochondrial genome is not a substrate of exonuclease (Fig. 5I).
DISCUSSION
Several serine/threonine phosphorylation sites have been identified in the papillomavirus protein E2 (3, 5, 6, 28, 29), but only two tyrosine phosphorylation sites have been examined: Y102 and Y138 (7, 30). Phosphorylation at Y102 was originally identified in BPV-1 E2 through mass spectrometric analysis (7) and subsequently confirmed in the current work for HPV-31 E2. Mutation of Y102 to the phospho-mimetic residue glutamic acid (Y102E) blocked viral genome replication. Since HPV-31 also has a tyrosine at residue 102, we sought to determine if this inhibitory phenotype is conserved between BPV-1 of the delta papillomavirus family and HPV-31 of the alpha papillomavirus family.
Glutamic acid (E) or phenylalanine (F) mutations were introduced at position Y102 to mimic either the phosphorylated or unphosphorylated state, respectively, of the tyrosine. When introduced into either BPV-1 or HPV-31, these mutations did not alter the normal nuclear localization of E2 and both the Y102E mutants no longer bound the C terminus of Brd4. These are the extent of the similarities. As summarized in Table 1, BPV-1 Y102E mutant was defective while HPV-31 E2 retained binding to E1, directly correlating with the results observed in transient DNA replication assays. Given that the TAD domain of BPV-1 and HPV-31 only share 28% identity, it is possible that the variations in the interface between E2 and its binding partners results in the differences in BPV-1 versus HPV-31 Y102E binding to TopBP1, GPS2, and E1. Other labs have also noted differences between BPV-1 and HPV-31 binding; for example, BPV-1 but not alpha papillomavirus E2 proteins segregate tightly with chromosomes during mitosis (31), and HPV-16 and BPV-1 may have different segregation receptors (18).
Summary of BPV-1 and HPV-31 E2 Y102E mutant phenotypes
Initial quasivirus infections utilized HPV-31 genome without any selection. Regardless of the cell line used, the high genome copy numbers found at 1 or 2 days posttransfection became negligible by days 4 and 7. One method to selectively maintain HPV in cell lines is to cotransfect the HPV-31 genome with another plasmid containing a selection gene, but in our experience the resulting cell lines preferentially retain the selectable plasmid and lose the HPV genome. Only upon substituting the neo cassette for L1/L2 sequences was HPV genome copy number reliably maintained in cells. HPV-31 L2 and L1 genes contain multiple CCCTC DNA sequences which confer CTCF binding (32). Mutation of these sites in L2 restricted episomal maintenance of the HPV-31 genome, as well as genome amplification (32); however, the effect of mutation of the L1 sites has not been examined. HPV-31neo retained the three L2 CCCTC sequences and deleted the two L1 sites. Our observed amplification and episomal maintenance of HPV-31neo genome suggests that these CCCTC sites within L1 are not required for HPV replication.
HPV-31 E2 Y102E no longer bound TopBP1, and, as a result, associated very tightly with the nuclear matrix. Other residues in the same region, N89/E90, have also been characterized as necessary for TopBP1 association (33). While TopBP1 is not essential in transcriptional activation assays, it has been implicated as crucial for either initiation of replication in early infection or segregation during mitosis (16–18). The reduced interaction of Y102E with TopBP1 may prevent episomal maintenance and mitotic segregation in cells infected with HPV-31neo Y102E. This could explain our observation that Y102E always formed at least 10-fold fewer colonies than WT when equal amounts of quasivirus were used for infections. Infections were successfully performed in multiple cell lines, and even though HaCaT cells are notorious for being unable to support long-term HPV episomes (34–36), HaCaT cells can do so with a cis-linked neo selection gene.
Early studies on E2 phosphorylation concluded that at least 20 phosphorylation sites are present in the N terminus (5), none of which were on a tyrosine (3). These experiments were based upon E2 expressed in insect cells that are not capable of fully modifying human proteins; there are only 3 tyrosine kinases in the Spodoptera frugiperda genome, instead of the full complement of 90 human kinases. Also limiting detection in these studies was the absence of tyrosine phosphatase inhibitors in the lysates. Tyrosine phosphorylations are notoriously short-lived, often with a half-life of only a few seconds due to highly active phosphatases (37). Unless the phospho-tyrosine site is bound to SH2/PTP domains, the site goes through hundreds of rounds of phosphorylation/dephosphorylation, which makes detection difficult in the absence of tyrosine phosphatase inhibitors. There may be a bias also to detection of serine/threonine phosphorylations, as 86% of all phosphorylation occurs on serines, 12% on threonine, and only 2% on tyrosines (38). It is therefore possible that additional tyrosines besides Y102 and Y138 can be phosphorylated, and these might regulate replication and genome maintenance.
Through mass spectrometry of HPV-31 E2, we identified tyrosine phosphorylation not only on Y102, but also on Y87 (Fig. 1D). Further studies are required to determine what effect phosphorylation has at this site. Additionally, we found phosphorylation at serine 269, the homologue of which was previously characterized as a minor phosphorylation site in HPV-1 (S281) (39). We have also detected phosphotyrosine at positions Y131 and Y138 in HPV-31 E2, though with lower than 75% probability scores. Our mass spectrometric analysis of BPV-1 E2 has also expanded the list of E2 phosphorylation sites that are fibroblast growth factor receptor 2 (FGFR2) dependent: Y32, Y44, Y131, Y158, Y159, Y169, Y170, and Y262 (40).
Another tyrosine mutant, Y131A, has a very similar phenotype to HPV-31 Y102E. HPV-16 Y131A prevents phosphorylation, increases nuclear matrix association, has no effect on transient replication, and cannot be maintained episomally in HFKs (41). A more direct comparison would require Y131F and Y131E mutations to be generated, but it appears phosphorylation of Y102 mimics nonphosphorylated Y131A. Phosphorylation at one site may restrict phosphorylation at the other, each the target of specific kinases that are activated through different pathways. The hydroxyl of Y102 and tyrosine 159 are only 11 Å apart and face the same direction in the crystal structure of 16E2 (PDB: 2NNU), so it is possible phosphorylation at one site influences subsequent phosphorylation at the other residue.
The kinase which phosphorylates at position 102 has yet to be identified. We previously published that that the FGFRs associate with E2 (40, 42), though none of these phosphorylate at Y102. Of the 90 known tyrosine kinases, a consensus sequence is only recognized for LCK, Abl, Src, ALK, and EGFR (38). HPV-31 E2 Y102 fits the epidermal growth factor receptor (EGFR) consensus sequence, but our previous observations that E2 does not bind EGFR and that EGFR does not inhibit replication (40) imply that EGFR is not the kinase for Y102. The kinase that phosphorylates BPV-1 may differ from HPV-31, as the sequence flanking Y102 differs significantly (Fig. 1B).
MATERIALS AND METHODS
Plasmids.For the replication assay, pFLORI-31, containing the firefly luciferase as well as HPV-31 origin of replication, and pRenilla, containing the Renilla luciferase behind a CMV promoter, were used (43). For immunoprecipitations, the following plasmids were used: 3×FLAG-HPV-31 E1 (codon optimized) in pCMV 3Tag 1a vector (43), pcDNA3-FLAG-HPV-31 E2 (codon optimized) (42), pHA-GPS2 (20), pCN-GST-hBrd4 1224–1362 (CTM) (44), FLAG-TopBP1 (45), p16LL (46), and pBR322-HPV-31, HA-HPV-31 E1.
Cloning.HPV-31neo. First, 31neoA, in which an AgeI restriction site was added to the L1 region (4573 bp) of pBR322-HPV-31 through site directed mutagenesis using primers designed by PrimerX (F: AGCATCTACCGGTACACCAGCA; R: TGCTGGTGTACCGGTAGATGCT). The neo cassette, encompassing the SV40promoter-intron-neomycin-βactin-poly-A-tail, was PCR amplified from the template pBR322-HPV-18neo (23) using primers which add an Sbf1 site (in lower case) to the 5′ end and an AgeI site to the 3′ end (F: GATCcctgcaggGGCCTGAAATAAC; R: GTACaccggtTAAAATACAGCATAGCAA), and inserted into Age/Sbf1-digested 31neoA. The neo cassette contains a HindIII site which was to be removed to permit release of the HPV-31neo genome from the pBR322 backbone through HindIII digestion. Similarly, the neo cassette contains a CsiI site that was removed to allow cut-and-paste insertion of E2 mutants into the HPV-31neo plasmid. Primers for HindIII removal (F: GCTTTTGCAAATCGCTTCTGCCTTC, R: GAAGGCAGAAGCGATTTGCAAAAGC) and CsiI removal (F: TCCACACCTGCTTGCTGACTA, R: TAGTCAGCAAGCAGGTGTGGA) were designed in the PrimerX website. The resulting pBR322-HPV-31neoHC was sequenced to confirm no mutations were incorporated during cloning. This plasmid is referred to as HPV-31neo and is available in the Addgene repository (Plasmid number 153283).
Mutants 31E2 Y102E/F. Site directed mutagenesis was performed following the QuikChange II protocol (Agilent Technologies) to generate HPV-31 E2 Y102E and Y102F mutants in codon-optimized E2 of pcDNA3-FLAG 31 E2, noncodon-optimized E2 of pSB-HA-HPV-31 E2 (8), and endogenous E2 of HPV-31neo. Primers were designed using PrimerX website (Y102E F: CAAACAAGTCTTGAACTGGAGTTAACTGCACCTACAGGG, R: CCCTGTAGGTGCAGTTAACTCCAGTTCAAGACTTGTTTG; Y102F F: CAAACAAGTCTTGAACTGTTCTTAACTGCACCTACAGGG, R: CCCTGTAGGTGCAGTTAAGAACAGTTCAAGACTTGTTTG; codon-optimized Y102E [+XcmI silent mutation for selection] F: CAGCCTCGAGCTGGAGCTGACCGCTCCCA, R: TGGGAGCGGTCAGCTCCAGCTCGAGGCTG; codon-optimized Y102F F: GACCAGCCTCGAGCTGTTCCTGACCGC, R: GCGGTCAGGAACAGCTCGAGGCTGGTC).
Cells and transfections.HEK 293TT, TTF, and C33a cells were cultured in Dulbecco’s modified Eagle medium (DMEM) with 10% fetal bovine serum (FBS; Atlas Biologicals). HaCaT cells were cultured in calcium-free DMEM supplemented with 5% FBS, which was calcium depleted for 1 h in Chelex-100 resin, then supplemented with CaCl2 added to a final concentration of 0.02 mM. N/TERT cells were grown in KSFM, and NIKS cells were grown in F-media. Polyethylenimine (PEI) was used for transfections at a ratio of 2 μg PEI to 1 μg DNA, which is equivalent to an N/P ratio of 15:1. PEI and DNA were diluted in Opti-MEM serum-free medium before placing on cells. At 20 h posttransfection, cells were harvested for lysis or fixed for immunofluorescence. For viral production in HEK 293 TTF cells, lipofectamine 2000 was used at a ratio of 2:1, as PEI tends to precipitate the large DNA. HEK 293TT cells were transfected with FLAG-HPV-31 E2 WT or Y102E and harvested 24 h later. Cell pellets were fractionated as previously described (47).
Transient DNA replication assays were performed with C33a cells in triplicate using Dual-Glo luciferase assay system (Promega) (22), with one modification: 20 μl Dual-Glo reagent was used instead of 50 μl. At 3 days posttransfection, firefly and Renilla luminescence was measured with a PHERAstar plate reader and software.
Sample preparation for mass spectrometry.HEK293 TT cells transfected with FLAG-HPV-31 E2 were grown on 6 × 15-cm dishes, washed in HBSS buffer (10 mM HEPES pH 7.5, 140 mM NaCl, 0.5 mM CaCl2, 0.8 mM MgCl2), and incubated with the irreversible tyrosine phosphatase inhibitor pervanadate (30 μM) for 1 h at 37°C. Cells were lysed in 20 mM HEPES pH 7.5, 150 mM NaCl, 30 μM pervanadate, 7.5 μM trichostatin A (TSA), and 0.525 mg/ml NaF and immunoprecipitated overnight with B201 antibody. Bands were excised from Coomassie-stained polyacrylamide gels for tandem mass spectrometry (MS-MS). The gel bands were subjected to reduction (10 mM DTT) and alkylation (55 mM iodoacetamide) and digested with trypsin (Promega) overnight at 37°C and injected into a C18 column. Peptide spectra were recorded in Orbitrap Velos Pro (Thermo Fisher Scientific) and Dionex UltiMate 3000 RSLC nanoflow system (Thermo Fisher Scientific) and database searches were carried out using Sequest algorithm. From the MS-MS data, individual probability values for each phosphorylation site were calculated using Proteome Discoverer V1.3 equipped with phospho RS 2.0.
Immunoprecipitations and immunoblotting.At 24 h posttransfection, cells were washed in PBS, harvested, and then lysed in lysis buffer (0.5% NP-40, 150 mM NaCl, 20 mM Tris pH 7.5). Insoluble protein was separated by centrifugation, and soluble lysate was quantitated using BCA reagent (Pierce Thermo Scientific). Equal amounts of lysate per immunoprecipitation (IP) were diluted to 500 μl in lysis buffer and incubated for 3 h at 4°C with 20 μl of 50% slurry M2 affinity gel (Sigma), 20 μl of 50% slurry glutathione beads, or 1 μl of TopBP1 antibody + 25 μl of 50% slurry protein A/G Sepharose beads. E1 and Brd4-CTM were also treated with 0.1 μg/μl ethidium bromide as previously described during immunoprecipitation (44). Beads were washed 5× in lysis buffer, boiled in 2× protein sample buffer, run on SDS-PAGE gels, and transferred onto 0.45-μm PVDF membranes (Millipore) in semidry transfer boxes (Bio-Rad). Membranes were blocked in 5% nonfat milk/TBST overnight and incubated with primary antibodies for 1 h at the following dilutions: FLAG: M2 (Sigma; 1:5,000), HA: 12CA5 supernatant 1:1,000, TopBP1 (Bethyl; 1:2,000), or GST (Arbor Vitae; 1:5,000). Light-chain-specific secondary antibodies diluted 1:5,000 were added following three TBST washes. Signals were detected with ECL (Amersham) and ImageQuant LAS 4000 system (GE Healthcare).
For phospho-Y102 antibody production, New Zealand white rabbits were immunized with an HPV-31 E2 phospho-tyrosine peptide (amino acids [aa] 97 to109; TSLELYLTAPTGC) formulated in complete Freund’s adjuvant followed by four injections in incomplete Freund’s adjuvant. Serum was analyzed by enzyme-linked immunosorbent assay (ELISA) to confirm antibody production. The rabbit 10008 antiserum predominantly recognizes phosphorylated Y102, but also detects several cellular phosphotyrosine proteins by immunoblot. This antiserum can be utilized for E2 chromatin immunoprecipitation (ChIP) assays (48).
Cell fractionation.Based upon published methods (47), HEK293 TT cells transfected with FLAG-HPV-31 E2 WT or Y102E were lysed in lysis buffer (15 mM Tris pH 8.0, 150 mM NaCl, 0.5% Triton-x) for 20 min, and then spun at 12,000 rpm for 10 min. The soluble fraction was kept on ice while the pellet was resuspended in lysis buffer and divided into two equal parts, A and B. The A portion was solubilized in an equal volume of 0.8 M NaCl for a final 0.4 M concentration and spun at 12,000 rpm. The soluble fraction was kept on ice. The B portion was digested with Benzonase at room temperature for 30 min, then spun at 7,000 rpm. The soluble fraction was kept on ice, while the pellet was resuspended in lysis buffer with 2 M NaCl and spun at 12,000 rpm. The soluble fraction was kept and the pellet was resuspended in 8 M Urea, 100 mM Tris pH 8.0, 100 mM NaCl.
Immunofluorescence.N/TERT cells were plated onto UV-sterilized glass coverslips at 3 × 105 cells per well of a 12-well plate. Cells were transfected with FLAG-HPV-31 E2 and HA-HPV-31 E1 constructs. After 24 h, these were fixed with 4% paraformaldehyde/phosphate-buffered saline (PBS) and permeabilized in blocking solution (5% goat serum, 0.25% Triton X-100, PBS). Primary M2 FLAG antibody diluted 1:5,000 and rat 3F10 HA antibody diluted 1:500 in blocking solution and incubated for 1 h at room temperature (RT). Secondary antibody was 594 Alexa-Fluor mouse and 488 Alexa-Fluor rat antibody diluted 1:1,000. Slides were mounted using Prolong Gold with DAPI (4′,6-diamidino-2-phenylindole), and visualized using a Nikon Microphot-SA microscope.
Quasivirus preparation and infection.Viral stocks were made according to a modified version of the original method (https://ccrod.cancer.gov/confluence/display/LCOTF/Home). A more detailed description of the method will be reported separately. To generate virus genomes, pBR322-HPV-31neoHC was digested with HindIII (Anza) for 2 h, gel purified to separate from the vector, and religated with T4 DNA ligase for 3 h. HEK293 TTF cells (49) were transfected with 5 μg of p16sheLL (46) and 5 μg of recircularized HPV-31neo. Two days after transfection, cells were resuspended in viral maturation buffer (0.5% Triton-X, 10 mM MgCl2, 5 mM CaCl2, PBS) at 1.5 cell pellet volume equivalents and incubated at 37°C for 2 days. Lysate was spun at 7,000 rpm and the soluble viral preparation was stored at 4°C and used within 1 month.
HaCaT cells were grown to confluence and treated with extracellular matrix (ECM) buffer (170 mM NH4OH, 0.5% Triton, PBS), which removes the cell body but leaves behind the ECM (25). ECM was washed 3× in PBS to remove buffer, then HaCaT medium was added along with quasivirus/pseudovirus (10 multiplicity of infection [MOI] for HaCaT, 40 MOI for NIKS). After a 4-h incubation of virus on ECM, HaCaT or NIKS cells were added. Two days after infection, Geneticin was added at 500 μg/ml for HaCaT or 200 μg/ml for NIKS to select for infected cells. For transfection, 1 μg recircularized HPV-31 genomes as prepared for viral production was transfected into N/TERT cells. Two days posttransfection, cells were cultured in 200 μg/ml Geneticin.
For quantitation of viral titer, 5 μl of virus was diluted 10× in 15 mM Tris pH 7.5, 1 mM MgCl2, 1 μl Benzonase, and digested at 37°C for 2 h. Samples were heated for 5 min to inactivate the DNase and break open the capsids. Proteinase K was added to degrade the proteins at 55°C for 30 min, followed by inactivation with 5 mM phenylmethylsulfonyl fluoride (PMSF). Quantitative PCR was then performed on the viral preparation, along with 500 pg to 0.05 pg genome standards for calibration.
DNA and RNA isolation.Trypsinized cells were counted and washed in PBS, and episomal DNA was isolated using the standard Hirt protocol (50). qPCR primers were used to amplify the HPV genome (F: ACACCATGCATTATACTAAC; and R: TCACTGCTAAATACAGATTC) using Ssofast Evagreen Supermix (Bio-Rad). RNA was isolated using standard TRIzol protocol, converted into cDNA using standard SuperScript III reverse transcriptase (Thermo Fisher), and qPCR primers to E1̂E4 for amplification. For DNA extraction and digestion with exonuclease V, previously described methods were used (27). In brief, DNA was isolated using the Qiagen DNeasy blood and tissue kit according to the standard protocol, and 10 μl of DNA was digested with 1 μl exonuclease V in NEB Cutsmart buffer supplemented with ATP for 3 h at 37°C. Exonuclease was heat-inactivated at 70°C for 30 min, and 1 μl of cut or uncut DNA was amplified using qPCR primers for the HPV genome as described earlier: actin (F: GAGGCACTCTTCCAGCCTTC; R: CGGATGTCCACGTCACACTT) (48), and mitochondria (F: CAGGAGTAGGAGAGAGGGAGGTAAG; R: TACCCATCATAATCGGAGGCTTTGG) (27).
Statistical analysis.For all experiments, means are expressed ± standard error of the mean (SEM). P values of ≤0.05 were considered significant. Student’s t test was performed for E2 fractionation comparing wild type to mutant. Each experiment was performed multiple times.
ACKNOWLEDGMENTS
We appreciate the generosity of the following people for cells and plasmids: Alison McBride (NIAID/NIH) for the codon-optimized HPV-31 E2 and selectable HPV-18 genome, Jacques Archambault (McGill University) for HPV-31 ori luciferase plasmids, John Schiller (NCI) for HEK 293TT cells, Richard Roden for HEK293 TTF cells, Iain Morgan (Virginia Commonwealth University) for N/TERT cells, and Cheng-Ming Chiang (University of Texas Southwestern Medical Center) for the Brd4 constructs. This research was supported by the National Cancer Institute (R01 CA058376 to E.A.) and the National Institute of Allergy and Infectious Diseases (T32 AI007637 and T32 AR062495 to T.G., T32 AI060519 to M.D., F30 AI114284 to S.P.C.). M.D. was also supported by an award from the Ralph W. and Grace M. Showalter Research Trust. The content is solely the responsibility of the authors and does not represent the official views of the NIH.
FOOTNOTES
- Received 1 April 2020.
- Accepted 17 May 2020.
- Accepted manuscript posted online 3 June 2020.
- Copyright © 2020 American Society for Microbiology.
