INTRODUCTION
The early 1B protein E1B-55K has been the focus of many recent investigations into the biology of human adenovirus species C type 5 (HAdV-C5, referred to as C5 in the manuscript). With its multiple functions, E1B-55K is an important regulator of various signaling pathways (
1,
2).
SUMOylation is a reversible posttranslational modification (PTM). An enzymatic cascade covalently binds small ubiquitin-like modifier (SUMO) proteins to a substrate protein with various consequences like regulation of protein-protein interactions, protein stability, functional activation, and subcellular localization (
3). It has previously been shown that HAdV-C5 E1B-55K SUMOylation at a consensus SUMO conjugation motif (SCM) at lysine 104 (K104) regulates numerous functions of this viral protein. E1B-55K SUMOylation regulates its own E3-SUMO-ligase activity, binding to p53 and repression of p53-dependent promoters (
4–8). For example, SUMO-E1B-55K can induce Sp100A SUMOylation to inhibit its function as a tumor suppressor and an activator of p53-dependent transcription (
9). Moreover, SUMOylation regulates E1B-55K-binding to and the subsequent degradation of Daxx, a PML nuclear body (PML-NB)-associated factor that inhibits adenoviral gene expression (
10). SUMO-E1B-55K also interacts with PML IV and V, which are components of PML-NBs (
11). Together with three cellular proteins (cullin 5 [Cul5], elongin B/C and RING-box protein 1 [Rbx1]), E1B-55K and E4orf6 form an E3-ubiquitin-ligase complex that recruits cellular proteins for proteasomal degradation (
12–18). On the other hand, SUMOylation of E1B-55K is reduced by the cellular heterochromatin-associated transcription factor KAP1 through a still unknown mechanism (
19). Finally, HAdV-C5 E1B-55K nuclear localization and intranuclear targeting are SUMOylation-dependent and SUMOylation occurs in close proximity to a nuclear export signal (NES) (
20,
21).
To date, there are only very few reports that describe E1B-55K functions of other HAdV species. For instance, it has been shown that the formation of the E3-ubiquitin-ligase complex with E4orf6 is highly conserved among species. However, heterogeneity exists in the cullin proteins that form part of the complex and in the cellular targets that are ubiquitinated in the presence of different E1B-55K proteins (
22). DNA ligase IV is efficiently degraded in the presence of all tested HAdV E1B-55K proteins and other cellular proteins like Mre11, integrin α-3 and BLM are only degraded by some HAdV species. Regarding p53-binding and -degradation, it has been shown that only HAdV-A12, C5 and F40 E1B-55Ks are able to both bind and degrade this tumor suppressor, while other species bind but do not efficiently degrade p53 (
22,
23). Additionally, a previous report demonstrated that E1B-55K-binding to specific PML isoforms and colocalization with PML-NBs is a conserved function among HAdV species (
24). The localization of E1B-55K in aggresomes (aggregations of misfolded proteins in eukaryotic cells), however, has only been described for HAdV-B16, C5, D9 and E4 (
24,
25). Notably, E1B-55K proteins from different HAdV species differ in size ranging from 472 to 496 amino acid residues (
1,
2). Nevertheless, we decided to keep the conventional nomenclature and use “E1B-55K” for the different proteins inspired by previous literature (
1,
2).
In the present study, we characterize the SUMOylation of the E1B-55K proteins from HAdV species A12, B16 and B34, C5, D9, E4, and F40 as representative members of different HAdV types (
1,
2,
26). We first identify the presence of functional SCMs and further investigate the consequences of SUMOylation on functions of these viral proteins. In addition, we show that a lysine residue at position 101 (K101) of HAdV-C5 E1B-55K affects SUMOylation and nucleo-cytoplasmic shuttling, most probably through regulation of the main SUMOylation-site K104.
DISCUSSION
In this study, we discovered that E1B-55K proteins from almost all HAdV species are highly SUMOylated at conserved SCMs, indicating that this PTM is not only important for HAdV-C5 E1B-55K but also for E1B-55Ks from other HAdV species (
Fig. 1 and
2). Intriguingly, HAdV-F40 E1B-55K is highly SUMOylated at K90 and a consensus SCM is missing (
Fig. 1 and
2). SUMOylation of this protein must occur at a nonconsensus site, as described for several proteins. Some examples are the cellular protein mouse double minute 2 homolog (Mdm2), the cAMP response element-binding protein (CREB) and the Epstein-Barr virus (EBV) transactivator protein BZLF1 (
39–42). Furthermore, this nonconsensus SUMO conjugation site in HAdV-F40 E1B-55K aligns with the SCM in the other HAdV species (
Fig. 1). Slight changes between our alignment and previously published E1B alignments are likely explained by the use of different software (
24,
26). Since accessibility to the substrate plays an important role in SUMOylation, it is conceivable that the position of the lysine residue N-terminal of E1B-55K is critical for a successful modification of the protein (
43). While the SCM is highly conserved among the different HAdV species, a lysine residue K101 and a leucine-rich NES are only found in HAdVs from species C, affecting the SUMOylation and nucleo-cytoplasmic shuttling of E1B-55K (
Fig. 2 and
3). A previous report suggested a model that illustrates the mechanism of E1B-55K shuttling. In this model, HAdV-C5 E1B-55K is SUMOylated at K104 upon entry into the nucleus. There, K104 SUMOylation facilitates targeting of E1B-55K to viral replication centers (RCs) and PML-NBs, simultaneously obstructing CRM1-binding to the NES. At these subnuclear structures, E1B-55K is deSUMOylated via an unknown mechanism, probably involving E4orf6, enabling a CRM1-dependent nuclear export (
21). In our study, we demonstrate that K101 is also involved in SUMOylation and shuttling of HAdV-C5 E1B-55K, but by a so far unknown mechanism. An explanation could be the nuclear localization of HAdV-C5 K101R that further enhances its SUMOylation and could therefore block nuclear export via the NES. Alternatively, SUMOylation could interfere with nuclear-cytoplasmic shuttling of E1B-55K and retain K101R in the nucleus. Despite the use of a CRM1-dependent export, HAdV-C5 E1B-55K is thought to have evolved another, CRM1-independent export pathway (
21). We provide evidence that E1B-55K from most HAdV species exit the nucleus via a CRM1-independent pathway because CRM1 is an export receptor for leucine-rich NES (
44) and the leucine-rich HAdV-C5 E1B-55K NES is only partially conserved and apparently absent in some HAdV species (A12 and F40) (
Fig. 1). Nevertheless, we observed a shift from nuclear to cytoplasmic localizations upon mutational SCM inactivation in E1B-55K from other HAdV species (
Fig. 3 to
5). Most interestingly, we also observed this change in HAdV-F40 E1B-55K, where the complete NES is presumably missing. Thus, it is likely that all tested E1B-55Ks evolved a CRM1-independent export pathway and only HAdV-C5 developed an additional CRM1-dependent mechanism to enable and facilitate nucleo-cytoplasmic shuttling of the protein. This suggests that E1B-55K SUMOylation in these species is regulating nuclear export. Furthermore, these results imply that high SUMOylation of E1B-55Ks from different HAdV species could block nuclear export and enhance nuclear retention, similar to what we observed with the HAdV-C5 K101R-mutant, our prototype mutant for high SUMOylation. It remains elusive whether E1B-55K from other HAdV species are able to actively shuttle between the nucleus and cytoplasm. Yet, our findings indicate that the mechanism is different from HAdV-C5. These remaining questions need to be addressed in follow-up studies with compounds targeting nuclear transport, ideally in the context of HAdV infection.
SUMOylation can have multiple effects on a protein and has been shown to alter protein-protein interactions, protein stability, localization, or activation (
3). In viruses, SUMOylation can additionally serve as a switch between a latent and a lytic infection. The Kaposi’s sarcoma-associated herpesvirus (KSHV) latency-associated nuclear antigen (LANA) is highly SUMOylated and by creating a SUMO-rich milieu required for chromatin condensation, the protein represses the expression of lytic KSHV genes (
45–47). Another example is EBV BZLF1, which also is a target of the SUMO conjugation machinery (
39). SUMOylation and deSUMOylation of BZLF1 are responsible for the onset of latency or the reactivation of the virus, respectively (
48). Notably, both herpesviral proteins are functional homologs of C5 E1B-55K and they all are involved in the repression of gene transcription and inhibition of p53 (
47,
49). We revealed that the lysine residue K101 is specific for HAdV-C and that SUMOylation of E1B-55K from this species is highly affected by K101 (
Fig. 1 and
2). Moreover, wt HAdV-C5 E1B-55K is less SUMOylated compared to E1B-55Ks from most HAdV species (
Fig. 2).
From previous studies, we already know that not all E1B-55K functions are conserved throughout the HAdV species. Although they all form an E3-ubiquitin-ligase complex with E4orf6, they do not degrade the same cellular proteins (
22,
23,
50). Importantly, p53 and Mre11 are only degraded by HAdV-A12, C5 and F40, but not by all HAdV species (
20,
46). Moreover, E1B-55K accumulation in aggresomes near the nucleus, which appears to be inhomogeneous in terms of distribution and morphology (
21,
23,
25,
51), is not conserved among HAdV species (
24). Here, we show that SUMOylation of E1B-55K is another mechanism that is highly conserved among the different HAdV species. Still, the regulatory role of SUMOylation that we know from HAdV-C5 E1B-55K is not entirely conserved. For instance, HAdV-B34 E1B-55Ks evolved a SUMO-independent mechanism to repress p53-mediated transcription (
Fig. 6). Why SUMOylation seems to regulate functions of E1B-55K from some but not all HAdV species could be due to differential phosphorylation (
26) and is an interesting subject for further studies. Likewise, it is critical to verify the observed effects of SUMOylation of different E1B-55K proteins in the context of infection. That would allow to control for the impact of E1A and especially the E1A or E1B-55K binding to the SUMO-conjugating enzyme UBC9 (
26,
52) as well as other viral proteins on the SUMOylation pathway. The use of viruses and virus mutants in infection along with thorough investigations of nuclear import and export pathways that are involved in the SUMO-dependent shuttling will be addressed in follow-up studies. That also applies for studies on the effect of E1B-55K SUMOylation on p53-mediated transcription to confirm the reporter assay findings.
In conclusion, we provide further evidence that the E1B-55K proteins from different HAdV species have conserved PTMs such as SUMOylation; however, E1B-55K SUMOylation in the different HAdV species might have different cellular consequences during infection. While HAdV-C5 still serves as a great molecular model, it is essential to continue investigating the functions of other E1B-55Ks in order to understand tissue tropism and pathogenesis of the different HAdV species.
MATERIALS AND METHODS
Cells.
H1299 cells (ATCC CRL-5803) and HeLa cells (ATCC CCL-2) stably expressing 6×His-SUMO-1 or 6×His-SUMO-2 (
53) were grown in monolayer cultures in Dulbecco’s modified Eagle’s medium (DMEM, Gibco) supplemented with 10% fetal bovine serum (FBS, PAN-Biotech), 100 U penicillin and 100 μg streptomycin (PAN-Biotech) per ml in a 5% CO
2 atmosphere at 37°C. HeLa cells expressing 6×His-SUMO were kept under 2 μM puromycin selection.
Plasmids and transient transfections.
The hemagglutinin (HA)-tagged E1B-55K proteins of the different HAdV species used in this study were expressed from pcDNA-3 as described previously (
22). Mutations in the (putative) SUMO-conjugation motif (SCM) of the respective E1B-55K types (
29,
30) were generated by site directed mutagenesis using primers shown in
Table 1. For transient transfection, subconfluent cells were treated with a mixture of DNA and 25 kDa linear polyethylenimine (PEI; Polysciences Inc., Warrington, PA, USA) as described before (
10).
Reporter gene assays.
For dual-luciferase assays, subconfluent cells were transfected as described above with 0.5 μg pRL-TK renilla under the control of the HSV-TK (herpes simplex virus thymidine kinase) promoter (Promega), 0.5 μg pRE Luc p53 promoter, 0.0015 μg p53 and 0.5 μg of the different E1B-constructs used in this study. Cells were lysed 24 h posttransfection and samples were measured and normalized using the Dual-Glo luciferase assay system (Promega) as described previously (
54).
Antibodies and protein analysis.
Primary antibodies against cellular and ectopically expressed proteins included 6×His mouse MAb (631212; TaKaRa Bio Europe SAS), HA rat MAb 3F10, p53 rabbit pAb FL-393 (sc-6243; Santa Cruz), and β-Actin mouse MAb AC-15 (A5441; Merck).
For protein extracts, cells were resuspended in RIPA lysis buffer (50 mM Tris-HCl/pH 8.0, 150 mM NaCl, 5 mM EDTA/pH 8.0, 0.1% SDS, 1% NP-40, 0.5% sodiumdeoxycholate) freshly supplemented with 1% phenylmethylsulfonyl fluoride, 0.1% aprotinin, 1 μg/ml pepstatin, 1 μg/ml leupeptin, 25 mM iodacetamide, and 25 mM
N-ethylmaleimide. Proteins were boiled 5 min at 95°C in 5× Laemmli. For immunoblotting, equal amounts of protein lysate were separated by SDS-polyacrylamide gel electrophoresis (SDS-PAGE), transferred to nitrocellulose membranes, incubated, and visualized as described previously (
54).
Denaturing, purification, and analysis of SUMO-conjugates.
HeLa cells were transfected with the different HA-tagged E1B-55K plasmids. 48 h posttransfection, cells were lysed in 5 ml GuHCl-buffer (6 M guanidinium-HCl, 0.1 M Na2HPO4, 0.1 M NaH2PO4, 10 mM Tris-HCl/pH 8.0, 20 mM imidazole, 5 mM β-mercaptoethanol). 10% of the cells were lysed with RIPA lysis buffer for total protein analysis as described above. The remaining lysates were incubated overnight at 4°C with 50 μl Ni-nitrilotriacetic acid (Ni-NTA) beads prewashed with GuHCl-buffer. Next, beads were washed once with buffer A (8 M urea, 0.1 M Na2HPO4, 0.1 M NaH2PO4, 10 mM Tris-HCl/pH 8.0, 20 mM imidazole, 5 mM β-mercaptoethanol) and buffer B (8 M urea, 0.1 M Na2HPO4, 0.1 M NaH2PO4, 10 mM Tris-HCl/pH 6.3, 20 mM imidazole, 5 mM β-mercaptoethanol). Afterwards, proteins were eluted and separated by SDS-PAGE as described above.
Indirect immunofluorescence assay.
Cells were grown on glass-coverslips and transfected with the indicated expression vectors. Cells were then fixed with 4% paraformaldehyde (PFA) for 20 min at 4°C, permeabilized in PBS with 0.5% Triton X-100 for 10 min at room temperature (RT) and blocked in Tris-buffered saline-BG (TBS-BG; 5% wt/vol BSA and 5% wt/vol glycine) for 30 min at RT. Coverslips were afterwards incubated for 1 h at RT with the primary antibody diluted in PBS, then washed 3× with TBS-BG. Secondary antibodies were diluted in PBS and coverslips were incubated for 30 min at RT (
Table 2). Coverslips were finally washed with TBS-BG and PBS and mounted in Glow medium. Images were obtained using a confocal spinning disk microscope (Nikon Eclipse Ti-E stand; Yokagawa CSU-W1 spinning disk; 2× Andor888 EM-CCD camera; Nikon 100× NA 1.49 objective) and an inverted fluorescent light microscope (Leica DMI 6000B). Images were analyzed with Fiji/ImageJ (NIH).
Colocalization between HA-E1B-55K and SUMO-2 was quantified with a custom written Matlab app. The app allows to select single cells for which the Manders fractional overlap coefficients M1 and M2 (
55) are calculated as measure of overlap between E1B-55K-foci (green signal) and SUMO-2 (red signal).
The ratio of the nuclear mean fluorescence to the cytoplasmic mean fluorescence was obtained in a three-step process: first, the nuclei were segmented using the U-Net based StarDist convolutional neural network (
56,
57), more specifically the pretrained “versatile fluorescent nuclei” network. Here, the cytoplasmic area was defined by expanding the area of the nuclei and subtracting the original nuclear area – resulting in a small ring-shaped cytoplasmic area around the nuclei. Second, the transfected nuclei (and the corresponding cytoplasmic areas) were determined by thresholding in the red (HA-E1B-55K) channel. Third, the mean nuclear and cytoplasmic fluorescence intensities were measured and their ratios compared.
Statistical analyses.
Statistical analyses were performed using Graph-Pad Prism v9 (GraphPad Software Inc.). Mann-Whitney tests and one-way analysis of variance (ANOVA) with Tukey post hoc tests for multiple comparisons were used to compare E1B-55K intracellular localization, E1B-55K colocalization with SUMO and p53-transactivation. Data were considered significant if P ≤ 0.05.