Human papillomaviruses (HPVs) are double-stranded DNA viruses that are known to be the causative agents of various cellular hyperproliferation-associated abnormalities (Munoz et al., 2003; Tommasino, 2014). More than 200 genotypes have been identified; among these, at least 14 are classified as “high-risk” types based on their pathogenicity in the host, including HPV16, HPV18, HPV31, and HPV45. They are well-defined carcinogens that cause several types of tumors, such as cervical cancer and head and neck squamous cell carcinoma (Kobayashi et al., 2018; Tumban, 2019; Wang et al., 2018). High-risk HPV genomes encode two oncogenic proteins that are responsible for HPV-driven tumorigenesis and malignancy: E6 and E7 (Mittal and Banks, 2017). HPV E6 recognizes p53, a key tumor suppressor protein, and provokes its proteasomal degradation by recruiting the E6AP ubiquitin ligase (Martinez-Zapien et al., 2016; Werness et al., 1990). Another oncogenic protein, HPV E7, contains three conserved regions (CRs). CR1 and CR2 are mainly responsible for targeting and promoting the degradation of pocket protein family members, such as retinoblastoma protein, a well-known cell cycle key regulator (Huh et al., 2007; Lee et al., 1998; Munger et al., 1989). CR3, which covers the C-terminal half of the E7 protein, adopts a homodimeric domain structure (Liu et al., 2006; Ohlenschlager et al., 2006) and regulates multiple viral processes involved in high-risk HPV-mediated tumorigenesis by interacting with several host factors (Avvakumov et al., 2003; Berezutskaya and Bagchi, 1997; Bodily et al., 2011; Brehm et al., 1999; Poirson et al., 2017; Todorovic et al., 2012). We recently determined the crystal structure of HPV18 E7 CR3 bound to the protein tyrosine phosphatase (PTP) domain of PTPN14, and this was the first report of the E7 CR3-mediated complex structure (Yun et al., 2019).

PTPN14 or PTPD2 is a cytosolic protein that contains an N-terminal four-point-one, ezrin-radixin-moesin (FERM) domain, central proline-rich motifs, and a C-terminal catalytic PTP domain. Given that PTPN14 is an important antitumor protein that cooperates with the p53 transcription factor and negatively regulates Hippo signaling-involved cell proliferation (Huang et al., 2013; Liu et al., 2013; Mello et al., 2017; Michaloglou et al., 2013; Wilson et al., 2014), it is not surprising that this protein is targeted by HPV E7 via direct binding, which leads to its proteasomal degradation (Hatterschide et al., 2019; Szalmas et al., 2017; White et al., 2016; Yun et al., 2019). Our previous study also revealed that PTPN21 or PTPD1 can directly interact with CR3 of HPV18 E7 (Yun et al., 2019), consistent with previous systematic proteomic analyses (Rozenblatt-Rosen et al., 2012; White et al., 2012). However, unlike PTPN14, which is defined as a tumor suppressor and particularly controls the Hippo signaling pathway (Huang et al., 2013; Liu et al., 2013; Mello et al., 2017; Michaloglou et al., 2013; Wilson et al., 2014), PTPN21 is known to be upregulated in several types of cancer cells and is associated with EGFR activation and intracellular trafficking, implying that this protein might be potentially oncogenic rather than tumor-suppressive (Carlucci et al., 2010; Plani-Lam et al., 2016; Roda-Navarro and Bastiaens, 2014; Wu et al., 2010). While the association between p53 and PTPN21 remains elusive, PTPN21 was reported to inhibit nuclear factor-κB transcriptional activity (Cho et al., 2017), which is known to play a pivotal role in immune responses and be negatively regulated by p53 (Carra et al., 2020; Jang et al., 2020; Pal et al., 2014; Schneider et al., 2010). Therefore, the biological roles and consequences of the HPV E7‒PTPN21 interaction in HPV-driven tumorigenesis might differ from those of the HPV E7‒PTPN14 interaction and warrant further investigation.

In this study, we elucidated the structural and physiological aspects of the human PTPN21 PTP domain and HPV18 E7 CR3 complex. We found that while the recombinant PTPN21 PTP domain shares high similarity with that of PTPN14 at the molecular level, the two PTP proteins differ from each other in terms of protein expression patterns and consequential cellular effects when expressed in HeLa cervical cancer cells and HaCaT keratinocytes. Therefore, these data expand our understanding of the protein interactions across species that are mediated by the HPV oncoprotein E7.

Preparation of recombinant protein constructs

Each of the DNA fragments encoding the human PTPN21 (residues 873-1174, referred to as PTP_PTPN21) and human PTPN14 (residues 886-1187, referred to as PTP_PTPN14) PTP domains were cloned into the modified pPRoEX HTa plasmid (Invitrogen, USA) with an N-terminal (His)₁₀-linked green fluorescent protein tag and the pET21a plasmid (Novagen, USA), respectively. Each of the E7 DNA fragments encoding CR3s of HPV18 (residues 54-105, referred to as CR3_18E7), HPV16 (residues 45-98), HPV1a (residues 39-93), and HPV11 (residues 45-98) were cloned into the modified pET28a vector (Novagen), with an N-terminal (His)₆-linked maltose-binding protein tag. The PTP_PTPN21 and CR3_18E7 constructs were used as templates for the preparation of a mutant PTPN21 protein containing F1032S·G1043Q·E1083A substitutions, referred to as PTP(SQA)_PTPN21, and a mutant HPV18 E7 CR3 protein containing R84A·L91A substitutions, referred to as CR3(AA)_18E7.

Recombinant protein production

The recombinant proteins used in this study were produced in the Escherichia coli BL21 (DE3) RIL strain (Novagen), which was cultured in Luria-Bertani media. Cells were incubated at 37°C until reaching an optical density of 0.6 at 600 nm, and then cultured for 16 h at 18°C following 0.2 mM isopropyl β-D-thiogalactopyranoside induction of protein expression. Cell lysates were subjected to purification steps separately for single protein purification or after proper mixing for protein complex purification. The purification was carried out using a Ni-NTA affinity chromatography column (Qiagen, Germany), HiTrap Q HP anion exchange chromatography column (GE Healthcare, USA), and HiLoad 26/600 Superdex 75 prep grade size-exclusion chromatography column (GE Healthcare). The (His)₁₀-linked green fluorescent protein, tagged to the PTPN21 proteins, and the (His)₆-linked maltose binding protein, tagged to the E7 proteins, were digested using TEV protease after Ni-NTA affinity chromatography and removed from the samples during the purification process. The protein samples were finally equilibrated with a buffer that consisted of 20 mM Tris-HCl (pH 7.5), 200 mM NaCl, and 1 mM dithiothreitol.

Small angle X-ray scattering (SAXS) experiments

SAXS data were collected on the beamline 4C at the Pohang Accelerator Laboratory at sample-to-detector distances of 1,000 and 4,000 mm (Kim et al., 2017). Measurements were conducted at 4°C with protein samples prepared with 1, 3, and 5 mg/ml, in triplicate. Data were obtained in ten successive frames of 10 s at a flow rate of 0.3 μl/s, using a Microlab 600 advanced syringe pump. An X-ray beam with a wavelength of 0.734 Å was used. Two-dimensional SAXS patterns were averaged and normalized by background subtraction using scattering intensities from the buffer solution. I(q), the scattering intensity data as a function of q, were obtained by radial averaging using the equation q = 4πsinθ/λ (θ, half of the scattering angle; λ, wavelength; 0.01 Å^-1 < q < 0.6 Å^-1). We used the CRYSOL (Svergun et al., 1995) and GNOM (Semenyuk and Svergun, 1991) programs to calculate the SAXS curves and distance distribution function P(r). Molecular envelope reconstructions were performed using the ab-initio shape determination program GASBOR (Svergun et al., 2001). The final refined models were obtained from 30 independently generated models that were compared and selected using the DAMAVER program. The PyMOL program was used for surface rendering. Superimposition of structural diagrams onto reconstructed dummy atoms was archived using the SUPCOMB program (Kozin and Svergun, 2001).

Isothermal titration calorimetry (ITC)

Protein samples containing 100 μM PTP_PTPN21/PTP(SQA)_PTPN21 samples or 5 μM CR3_18E7/CR3(AA)_18E7 were prepared in a solution containing 50 mM Tris-HCl (pH 7.5), 200 mM NaCl, and 10 mM β-mercaptoethanol. We carried out vacuum degasification for 20 min prior to all experiments. The measurements were performed on a VP-ITC microcalorimetry system (MicroCal, UK) at 25°C and analyzed using Origin software (OriginLab, USA) as described previously (Yun et al., 2019).

Size-exclusion chromatography-multi angle light scattering (SEC-MALS)

The SEC-MALS analysis was carried out using a Superdex^TM 75 10/300 GL column (GE Healthcare) with 5 mg/ml protein samples in a buffer that contained 50 mM Tris-HCl (pH 7.5), 200 mM NaCl, and 1 mM dithiothreitol. Optilab T-rEX (Wyatt Technology Corporation, USA) combined with high-performance liquid chromatography (Shimadzu, Japan) and DAWN HELEOS-II (Wyatt Technology Corporation) was used to record the differential refractive index spectra. ASTRA 6 software (Wyatt Technology Corporation) was used to calculate the weight-average molar mass.

Plasmid preparation for the cellular assays

Plasmids p6640 MSCV-P C-FlagHA 16E7-Kozak (#35018), p6641 MSCV-P C-FlagHA 18E7-Kozak (#35019), p6644 MSCV-P C-FlagHA 45E7-Kozak (#35135), MSCV-C 11E7 (#37912), and MSCV-C 6b (#37917) were purchased from Addgene (USA). They were used for transient transfection and as templates for construct preparation. Wild-type (WT) HA-PTPN21 was generously provided by Prof. Sayeon Cho (Chung-Ang University, Korea) and used as a template to generate a binding-defective mutant construct by site-directed mutagenesis (Cho et al., 2017). The constructs were then subcloned into the pcDNA3-Flag plasmid (Addgene).

Cell culture and transfection

HeLa and C33a cervical cancer cells, HaCaT human keratinocytes, and 293T cells were cultured at 37°C and 5% CO₂ in Dulbecco’s modified Eagle’s medium (DMEM; Corning, USA) containing 1% antibiotic-antimycotic (Gibco, USA) and 10% fetal bovine serum (GE Healthcare). Transient plasmid DNA transfection was conducted using Lipofectamine 3000 (Invitrogen) for HeLa and C33a cells and polyethylenimine (Sigma-Aldrich, USA) for 293T cells, according to the manufacturer’s instructions. HaCaT cells stably expressing HPV18 E7 were prepared as previously reported (Yun et al., 2019). For the delivery of siRNAs, HaCaT and HeLa cells were seeded and transfected into a 6-well culture plate with 20 nM SMARTpool ON-TARGET plus PTPN21 (SO-2868072G) and a non-targeting pool (D-001810-10; Dharmacon, USA) using LipoRNAiMAX (Invitrogen). siRNA-treated cells were harvested 48 h post-transfection for reverse transcription-polymerase chain reaction (RT-PCR) analysis.

RT-PCR analysis

Total RNA was extracted using the TRIzol reagent (Ambion, USA), and cDNA was obtained using M-MLV reverse transcriptase (Promega, USA) according to the manufacturer’s protocol. A total of 1 μg RNA was reverse-transcribed using Oligo dT primers for amplification. For quantitative RT-PCR, PTPN21 primers were used in a 2X Real-Time PCR Smart Mix (SolGent, Korea) and EvaGreen (SolGent), which was run on a real-time PCR (Bio-Rad, USA). β-Actin was used for normalization. The following primer sequences were used: PTPN21 (forward), 5′-TGAAACACAAGAATGGAAGGC-3′; PTPN21 (reverse), 5′-TCCTCTTTATTTGCCAGCTCT-3′.

Immunoblotting and immunoprecipitation

For immunoblotting, phosphate-buffered saline-resuspended cell lysates were mixed with 2X SDS sample buffer in a 1:1 ratio and boiled for 10 min. For immunoprecipitation, cells were lysed in a buffer containing 50 mM Tris-HCl (pH 7.5), 150 mM NaCl, 1 mM ethylenediaminetetraacetic acid, 0.5% Triton X-100, and a protease inhibitor cocktail (Roche, Switzerland). Proteins from cell extracts were sequentially immunoprecipitated with anti-Flag M2 agarose (A2220; Sigma-Aldrich) overnight at 4°C and Protein G-Sepharose (GE Healthcare) for 2 h at 4°C. Samples were mixed with the 2X SDS sample buffer in a 1:1 ratio and boiled for 5 min. Immunoprecipitated proteins were analyzed by western blotting using anti-HA (12013819001; Sigma-Aldrich), anti-Flag (F3165; Sigma-Aldrich), and anti-β-actin (A300-491A; Bethyl Laboratories, USA). Proteins were detected using enhanced chemiluminescence reagents (GE Healthcare) according to the manufacturer’s instructions. MG132 was purchased from A.G. Scientific (M-1157).

Cell proliferation assays

For cell growth assays, HeLa and HaCaT cells were transfected with genes coding for PTPN21 constructs for 24 h or with PTPN21-targeting or non-targeting siRNA for 48 h and seeded in a 48-well plate, at 2,000 cells per well. Cells were cultured in 200 μl DMEM supplemented with 1% antibiotic-antimycotic solution and 10% fetal bovine serum, without medium change. The CellTiter-Glo luminescent cell viability assay kit (Promega) was used to analyze cell proliferation by measuring luminescence after a 10-min incubation with a luminometer Wallac Victor^TM X3 (Perkin Elmer, USA).

Migration and invasion assay

Costar Transwell^® chambers containing 8.0 μm polycarbonate membrane (CLS3422; Corning) were used to analyze cell migration and invasion. HeLa and HaCaT cells were collected and suspended in 100 μl of serum-free DMEM at a density of 5 × 10³ cells/ml and 1.5 × 10⁵ cells/ml, respectively, to analyze migration. HeLa cells were also collected and suspended in 100 μl of the same medium, at a density of 5 × 10⁵ cells/ml, to analyze invasion. Medium supplemented with 10% FBS was subsequently added to the wells. Prior to seeding for the invasion assay, a thin layer of Matrigel (354234; Corning) was added to the membrane matrix, and then 100 μl of prepared cells were added on top of the inserts and Matrigel-coated inserts. The chamber plate was incubated at 37°C and 5% CO₂ for 24 h. Invaded cells were fixed with 4% paraformaldehyde (Biosesang, Korea) in phosphate-buffered saline for 30 min, and stained with 0.5% crystal violet (Tokyo Chemical Industry, Japan) for 20 min. Cells were observed and imaged using a microscope Eclipse Ts2 (Nikon, Japan).

HPV18 E7 interacts similarly with PTPN21 and PTPN14

We recently determined the crystal structure of CR3_18E7 bound to PTP_PTPN14 and demonstrated that each homodimeric CR3_18E7 protomer binds a single PTP_PTPN14 molecule and that the two PTP_PTPN14 molecules in a 2:2 complex do not mutually interact with each other (Yun et al., 2019). Our previous study also verified the molecular association of CR3_18E7 with PTP_PTPN21 using biochemical analyses (Yun et al., 2019). PTP_PTPN21 exhibits high amino acid sequence homology (71% identity and 88% similarity) with PTP_PTPN14 (Fig. 1A). Moreover, PTP_PTPN14 residues, which are known to play an important role in binding CR3_18E7, are highly conserved in PTP_PTPN21 (Fig. 1A). Therefore, we hypothesized that CR3_18E7 recognizes PTP_PTPN21 similarly as it recruits PTP_PTPN14 by constituting a 2:2 complex. To verify this hypothesis, we prepared recombinant CR3_18E7 bound to PTP_PTPN21 or PTP_PTPN14 proteins (Fig. 1B) and conducted a SAXS analysis (Table 1). We first confirmed that the experimentally observed scattering intensity curves and distance distribution function P(r) of the two complexes matched well with the theoretical values obtained from a calculation that used the crystal structure of CR3_18E7 bound to PTP_PTPN14 (PDB code 6IWD; Fig. 1C). Moreover, the X-ray crystallography-determined CR3_18E7−PTP_PTPN14 structure could be incorporated into the SAXS analysis-based molecular envelopes of the same protein complex (Fig. 1D, right) and the CR3_18E7−PTP_PTPN21 complex (Fig. 1D, left). These data indicate that the homodimeric CR3_18E7 constitutes the 2:2 complex together with two PTP_PTPN21 molecules similarly as it does with two PTP_PTPN14 molecules.

Modeling and verification of the E7−PTPN21 intermolecular interaction

As PTP_PTPN14 and PTP_PTPN21 share amino acid sequence (Fig. 1A) and CR3_18E7-binding mode (Figs. 1C and 1D) similarities, we modeled intermolecular CR3_18E7−PTP_PTPN21 contacts based on the CR3_18E7−PTP_PTPN14 crystal structure. In this model, the complex formation between the two proteins is mediated by hydrophobic interactions that involve Leu1014, Arg1031, Phe1032, Gly1043, and Trp1058 of PTPN21 and Met61, Phe90, and Leu91 of HPV18 E7 (Fig. 2A). Glu1083 of PTPN21 is expected to form an electrostatic interaction with Arg84 of HPV18 E7, which might also contribute to complex formation (Fig. 2A). Therefore, based on this complex model, we selected Phe1032, Gly1043, and Glu1083 of PTPN21 and Arg84 and Leu91 of HPV18 E7 as potential key residues for intermolecular interaction studies between the two proteins. The ITC measurements confirmed that the high-affinity binding between PTP_PTPN21 and CR3_18E7, with a dissociation constant of 5.43 nM, was completely abrogated by introducing the F1032S∙G1043Q∙E1083A (referred to as SQA in this study) triple mutation in PTPN21 or the R84A∙L91A (referred to as AA in this study) double mutation in HPV18 E7 (Fig. 2B). These findings imply the relevance of the homologous structure-based modeling of CR3_18E7−PTP_PTPN21.

Given the structurally and biochemical similarity between CR3_18E7−PTP_PTPN14 and CR3_18E7−PTP_PTPN21, we hypothesized that PTP_PTPN14 and PTP_PTPN21 are competitive with each other for binding CR3_18E7. To address this hypothesis, recombinant CR3_18E7−PTP_PTPN14 protein complex was incubated with PTP_PTPN21 and then subjected to size-exclusion chromatography. We determined that administration of PTP_PTPN21 into the CR3_18E7−PTP_PTPN14 complex resulted in the isolation of PTP_PTPN14 from the complex (Fig. 2C; marked by a red arrow) and co-elution of CR3_18E7 together with PTP_PTPN21 (Fig. 2C; marked by a green arrow). These results demonstrate that PTP_PTPN21 is able to replace PTP_PTPN14 from the CR3_18E7-bound complex, and that the two PTPs recognize CR3_18E7 in a competitive manner.

E7s from various HPV genotypes interact with PTPN21

Given the molecular similarity between PTP_PTPN14 and PTP_PTPN21, we assumed that PTP_PTPN21 associates with E7s from a wide range of HPV genotypes, as PTP_PTPN14 does (Poirson et al., 2017; Szalmas et al., 2017; White et al., 2016; Yun et al., 2019). To verify this hypothesis, we first carried out SEC-MALS using recombinant PTP_PTPN21 proteins in HPV E7-bound or unbound forms. As shown in Fig. 3A, PTPN21 (~38 kDa) was determined to form a stable complex with E7s (estimated to be 13 kDa as a homodimer) from HPV1a, HPV11, HPV16, and HPV18 with molecular weights of ~86 kDa, which corresponds to a 2:2 complex. Next, we conducted a coimmunoprecipitation assay using PTPN21 and E7 proteins transiently expressed in 293T cells. PTPN21 was coimmunoprecipitated with E7 from both high-risk HPVs, including HPV16, HPV18, and HPV45, and low-risk HPVs, including HPV6b and HPV11 (Fig. 3B). Collectively, these results demonstrate that E7 proteins across a wide range of HPV genotypes are able to interact with PTPN21, as they do with PTPN14.

PTPN21 is not subjected to proteasomal degradation in HeLa cells

In our previous study, the interaction between PTPN14 and HPV18 E7 was verified to induce PTPN14 degradation in HPV18-positive HeLa cervical cancer cells, which was delayed by the introduction of binding-defective mutations or proteasome inhibitor MG132 treatment (Yun et al., 2019). We investigated whether the interaction between PTPN21 and HPV18 E7 results in similar biological consequences. CR3_18E7 and full-length PTPN21 were transiently expressed in HPV18-positive HeLa and HPV-negative C33a cervical cancer cells, in WT and binding-defective mutant forms. The immunoprecipitation results verified that the intermolecular interaction was detected between the WT proteins, but not when binding-defective mutations were introduced into PTPN21 or HPV18 E7, in both HeLa and C33a cells (Fig. 4A). Remarkably, we found that the protein levels of the WT and SQA mutant PTPN21 constructs differed from those of the PTPN14 constructs when they were expressed in HeLa cells. The cellular amount of E7 binding-defective SQA mutant PTPN14 was significantly higher than that of WT, and this protein was sensitive to proteasome inhibitor MG132 treatment (Yun et al., 2019). Contrastively, the WT and binding-defective PTPN21 protein levels were comparable to each other and hardly affected by the treatment of the proteasome inhibitor MG132 when expressed in HeLa cells (Fig. 4B). Consistently, introduction of HPV18 E7 binding-defective mutations did not cause any noticeable difference of the protein level patterns of PTPN21 when expressed in HPV-free 293T cells with increasing transfection doses of full-length or the CR3 domain of HPV18 E7 (Fig. 4C). These results suggest that PTPN21 was not targeted by the E7-induced proteasomal degradation system in HPV18-positive cervical cancer cells, unlike PTPN14, despite the high similarity in their E7-binding mode (Fig. 1).

PTPN21 potentially functions as an oncoprotein rather than a tumor-suppressor

PTPN14 is a well-known tumor suppressor, and our previous study demonstrated that E7 binding-defective mutations remarkably enhance the ability of PTPN14 to delay HeLa cell proliferation, migration, and invasion (Yun et al., 2019). On the other hand, the physiological function of PTPN21 in cervical cancer development is not yet clearly understood, although PTPN21 was reported to be upregulated in several particular cancer types, such as B-cell non-Hodgkin’s gastric lymphoma and bladder cancer (Carlucci et al., 2010; Plani-Lam et al., 2016). Unlike PTPN14, PTPN21 exhibited a proteasomal degradation-evasive character, and E7-binding capability did not notably affect its protein level (Fig. 4). Therefore, we investigated whether PTPN21 could influence cell proliferation, migration, and invasion in cervical cancer cells. First, we examined how the transient expression of the full-length WT and SQA mutant PTPN21 affects the proliferation of HeLa cervical cancer cells and nontumorigenic HaCaT keratinocytes. We found that PTPN21 expression did not restrict cell proliferation in either cell type, irrespective of protein form (Fig. 5A). Consistently, siRNA-mediated suppression of PTPN21 expression hardly affected the growth of HeLa and HaCaT cells (Fig. 5B), which indicates that PTPN21 is not a proliferation regulator of these cell types.

We next analyzed the effects of PTPN21 ablation on HeLa cell migration/invasion and HPV18 E7-expressing or control HaCaT cell migration. Fig. 6A shows that HeLa cell migration and invasion were critically attenuated by the genetic knockdown of PTPN21. Moreover, HaCaT cell migration was significantly repressed with the siPTPN21 treatment, in control cells and in cells stably expressing HPV18 E7 (Fig. 6B). These results collectively imply that PTPN21 positively controls cell migration and invasion. Taken together, this protein, which was bound to HPV18 E7 (Figs. 1-3) and evades proteasomal degradation (Fig. 4), might function as an oncoprotein rather than a tumor-suppressor in HPV18-positive HeLa cervical cancer cells and in keratinocytes that are vulnerable to HPV infection.

The recognition of human proteins by E6- or E7-associated direct interaction and subsequent proteasomal degradation of the recruited host proteins are among the main strategies of HPV for inducing the transformation and immortalization of HPV-infected cells. Based on the determination of the complex crystal structure accompanied by biochemical and cellular analysis, we previously verified that tumor-suppressive PTPN14 is targeted by HPV18 E7 for degradation through a potent direct interaction (Yun et al., 2019). Our previous report also showed that PTPN21 is another binding partner associated with HPV18 E7 (Yun et al., 2019). In this subsequent study, we elucidated the molecular aspects of this novel interaction, which was accomplished by combining biochemical tools, including SAXS (Fig. 1), structure-based modeling (Fig. 2A), and SEC-MALS (Fig. 3). These results collectively demonstrated that PTP_PTPN21 forms a 2:2 complex with CR3_18E7, and defined the binding interface between the two proteins, which was confirmed by binding affinity measurements (Fig. 2B) and immunoprecipitation assays (Fig. 4A) carried out using WT and binding-defective mutant constructs. Despite the high similarity between PTPN14 and PTPN21 at the molecular level, we noticed that the biological consequences of HPV E7-binding differ between the two PTPs. First, while the E7 binding-defective mutant PTPN14 protein level was more than four-times higher than that of WT when expressed in HPV18-positive HeLa cells (Yun et al., 2019), these values were comparable in the case of the PTPN21 constructs (Figs. 4B and 4C). Second, MG132 protease inhibitor treatment led to a more than three-fold increase in WT PTPN14 protein levels in HeLa cells, but it did not induce such a prominent change in PTPN21 expression (Fig. 4B). Third, whereas PTPN14 functioned as a negative regulator of HeLa cell proliferation, migration, and invasion, and E7-binding defective mutations reinforced its tumor-suppressive activities, neither the ectopic expression nor genetic knockdown of PTPN21 relevantly affected HeLa or HaCaT cell proliferation (Fig. 5). Instead, the cellular phenotypes directly associated with tumor development, such as cell migration and invasion, were significantly attenuated by the repression of PTPN21 expression in HeLa and HPV18 E7-expressing HaCaT cells, which strongly suggest that PTPN21 is putatively oncogenic rather than tumor-suppressive, as implied previously (Carlucci et al., 2008; 2010; Plani-Lam et al., 2016; Roda-Navarro and Bastiaens, 2014; Wu et al., 2010). In conclusion, PTPN21 does not appear to be the target of the HPV E7-dependent host protein degradation system, despite its high structural homology with PTPN14 in the E7 binding mode (Fig. 1) and potent binding affinity with HPV18 E7 (Fig. 2).

What might be the biological consequence of the interaction between PTPN21 and HPV18 E7? One possibility is that HPV might control the localization or host protein association of PTPN21 through direct binding to E7, and thus take advantage of its oncogenic properties. Another possibility is that PTPN21 might be a collateral binding partner of HPV E7 due to the high structural similarity between PTP_PTPN14 and PTP_PTPN21. At present, the precise biological function of this interaction during viral infection and E7-mediated tumorigenesis remains unclear. Another unresolved issue is how PTPN21 evades proteasomal degradation despite its potent interaction with HPV18 E7. PTPN14 degradation is known to be mediated by the E3 ubiquitin-protein ligase UBR4/p600, which is recruited through direct interaction with the N-terminal region of HPV E7 (Szalmas et al., 2017; White et al., 2016). Given that PTPN14, but not PTPN21, was subjected to proteasomal degradation (Fig. 4), we suggest that PTPN21 might somehow hinder intermolecular interactions between the N-terminal region of E7 and UBR4/p600. Interestingly, we found that expression of full-length HPV18 E7, but not that of the CR3 domain, stabilizes both WT and CR3 domain-binding defective mutant PTPN21 (Fig. 4C), in which the N-terminal regions of these proteins are presumed to be involved. We thus examined whether PTPN21 has another region that associates with HPV18 E7 and found an unidentified interaction between the PTPN21 N-terminal FERM domain and full-length HPV18 E7 using immunoprecipitation (unpublished data). The relevance and biological importance of this novel interaction are our next objectives to determine.

In this study, we analyzed the molecular basis of how HPV18 E7 associates with human PTPN21 and revealed the physiological aspects of this interaction in HeLa cervical cancer cells and HaCaT keratinocytes. Our study expands our understanding of the diverse host protein recognition mediated by the HPV-encoded oncoprotein E7. We also anticipate that our observations will eventually serve as a rational basis for therapeutic approaches to treat HPV-implicated malignant diseases.

SAXS measurements were performed at the beamline 4C at the Pohang Accelerator Laboratory in Korea. We appreciate Dr. Eunha Hwang (Korea Basic Science Institute, Korea) for help with the SEC-MALS measurements and Dr. Jungwon Hwang and Dr. Myung Hee Kim (Korea Research Institute of Bioscience and Biotechnology, Korea) for help with the ITC measurements. This study was supported by the National Research Foundation of Korea (NRF_2020R1C1C1008451 and NRF_2019M3E5D6063955 to B.K., NRF_2017M3A9G5083321 to S.C.L., and NRF_2019R1C1C1002831 to E.W.L.) and the KRIBB Research Initiative Programs (to B.K.), which were funded by the Ministry of Science and ICT (MSIT) of Republic of Korea.

H.S.L. and M.W.K. performed and analyzed experiments and wrote the manuscript. K.S.J. performed and analyzed experiments. H.C.S. and W.K.K. provided expertise and feedback. S.C.L. provided expertise and feedback and secured funding. S.J.K. provided supervision. E.W.L. conceived and analyzed experiments, secured funding, and provided supervision. B.K. conceived and analyzed experiments, wrote the manuscript, secured funding, and provided supervision.

The authors have no potential conflicts of interest to disclose.

Structural parameters obtained from the SAXS data